Reticulate evolution and marine organisms: the final frontier?

Topology Testing and Demographic Modeling Illuminate a Novel Speciation Pathway in the Greater Caribbean Sea Following the Formation of the Isthmus of Panama

Recent genomic analyses have highlighted the prevalence of speciation with gene flow in many taxa and have underscored the importance of accounting for these reticulate evolutionary processes when constructing species trees and generating parameter estimates. This is especially important for deepening our understanding of speciation in the sea where fast-moving ocean currents, expanses of deep water, and periodic episodes of sea level rise and fall act as soft and temporary allopatric barriers that facilitate both divergence and secondary contact. Under these conditions, gene flow is not expected to cease completely while contemporary distributions are expected to differ from historical ones. Here, we conduct range-wide sampling for Pederson's cleaner shrimp (Ancylomenes pedersoni), a species complex from the Greater Caribbean that contains three clearly delimited mitochondrial lineages with both allopatric and sympatric distributions. Using mtDNA barcodes and a genomic ddRADseq approach, we combine classic phylogenetic analyses with extensive topology testing and demographic modeling (10 site frequency replicates × 45 evolutionary models × 50 model simulations/replicate = 22,500 simulations) to test species boundaries and reconstruct the evolutionary history of what was expected to be a simple case study. Instead, our results indicate a history of allopatric divergence, secondary contact, introgression, and endemic hybrid speciation that we hypothesize was driven by the final closure of the Isthmus of Panama and the strengthening of the Gulf Stream Current ~3.5 Ma. The history of this species complex recovered by model-based methods that allow reticulation differs from that recovered by standard phylogenetic analyses and is unexpected given contemporary distributions. The geologically and biologically meaningful insights gained by our model selection analyses illuminate what is likely a novel pathway of species formation not previously documented that resulted from one of the most biogeographically significant events in Earth's history.

Phylogenomic analyses of all species of swordtail fishes (genus Xiphophorus) show that hybridization preceded speciation

Hybridization has been recognized to play important roles in evolution, however studies of the genetic consequence are still lagging behind in vertebrates due to the lack of appropriate experimental systems. Fish of the genus Xiphophorus are proposed to have evolved with multiple ancient and ongoing hybridization events. They have served as an informative research model in evolutionary biology and in biomedical research on human disease for more than a century. Here, we provide the complete genomic resource including annotations for all described 26 Xiphophorus species and three undescribed taxa and resolve all uncertain phylogenetic relationships. We investigate the molecular evolution of genes related to cancers such as melanoma and for the genetic control of puberty timing, focusing on genes that are predicted to be involved in pre-and postzygotic isolation and thus affect hybridization. We discovered dramatic size-variation of some gene families. These persisted despite reticulate evolution, rapid speciation and short divergence time. Finally, we clarify the hybridization history in the entire genus settling disputed hybridization history of two Southern swordtails. Our comparative genomic analyses revealed hybridization ancestries that are manifested in the mosaic fused genomes and show that hybridization often preceded speciation.

Contrasting effect of hybridization on genetic differentiation in three rockfish species with similar life history

Hybridization can provide evolutionary benefits (e.g., population resilience to climate change) through the introduction of adaptive alleles and increase of genetic diversity. Nevertheless, management strategies may be designed based only on the parental species within a hybrid zone, without considering the hybrids. This can lead to ineffective spatial management of species, which can directly harm population diversity and negatively impact food webs. Three species of rockfish (Brown Rockfish (Sebastes caurinus), Copper Rockfish (S. auriculatus), and Quillback Rockfish (S. maliger)) are known to hybridize within Puget Sound, Washington, but genetic data from these species are used to infer population structure in the entire genus, including in species that do not hybridize. The goal of this project was to estimate the hybridization rates within the region and determine the effect of hybridization on geographic patterns of genetic structure. We sequenced 290 Brown, Copper, and Quillback rockfish using restriction-site associated DNA sequencing (RADseq) from four regions within and outside Puget Sound, Washington. We show that (i) hybridization within Puget Sound was asymmetrical, not recent, widespread among individuals, and relatively low level within the genome, (ii) hybridization affected population structure in Copper and Brown rockfish, but not in Quillback Rockfish and (iii) after taking hybridization into account we found limited directional dispersal in Brown and Copper rockfish, and evidence for two isolated populations in Quillback Rockfish. Our results suggest that rockfish population structure is species-specific, dependent on the extent of hybridization, and cannot be inferred from one species to another despite similar life history.

Phylogenomics analyses of all species of Swordtails (Genus Xiphophorus ) highlights hybridization precedes speciation

Hybridization has been recognized as an important driving force for evolution, however studies of the genetic consequence and its cause are still lagging behind in vertebrates due to the lack of appropriate experimental systems. Fish of the central American genus Xiphophorus were proposed to have evolved with multiple ancient and ongoing hybridization events, and served as a valuable research model in evolutionary biology and in biomedical research on human disease for more than a century. Here, we provide the complete genome resource and its annotation of all 26 Xiphophorus species. On this dataset we resolved the so far conflicting phylogeny. Through comparative genomic analyses we investigated the molecular evolution of genes related to melanoma, for a main sexually selected trait and for the genetic control of puberty timing, which are predicted to be involved in pre-and postzygotic isolation and thus to influence the probability of interspecific hybridization in Xiphophorus . We demonstrate dramatic size-variation of some gene families across species, despite the reticulate evolution and short divergence time. Finally, we clarify the hybridization history in the genus Xiphophorus genus, settle the long dispute on the hybridization origin of two Southern swordtails, highlight hybridizations precedes speciation, and reveal the distribution of hybridization ancestry remaining in the fused genome.

Genetic divergences and hybridization within the Sebastes inermis complex

The Sebastes inermis complex includes three sympatric species (Sebastes cheni, viz Sebastes inermis, and Sebastes ventricosus) with clear ecomorphological differences, albeit incomplete reproductive isolation. The presence of putative morphological hybrids (PMH) with plausibly higher fitness than the parent species indicates the need to confirm whether hybridization occurs within the complex. In this sense, we assessed the dynamics of genetic divergence and hybridization within the species complex using a panel of 10 microsatellite loci, and sequences of the mitochondrial control region (D-loop) and the intron-free rhodopsin (RH1) gene. The analyses revealed the presence of three distinct genetic clusters, large genetic distances using D-loop sequences, and distinctive mutations within the RH1 gene. These results are consistent with the descriptions of the three species. Two microsatellite loci had signatures of divergent selection, indicating that they are linked to genomic regions that are crucial for speciation. Furthermore, nonsynonymous mutations within the RH1 gene detected in S. cheni and "Kumano" (a PMH) suggest dissimilar adaptations related to visual perception in dim-light environments. The presence of individuals with admixed ancestry between two species confirmed hybridization. The presence of nonsynonymous mutations within the RH1 gene and the admixed ancestry of the "Kumano" morphotype highlight the potential role of hybridization in generating novelties within the species complex. We discuss possible outcomes of hybridization within the species complex, considering hybrid fitness and assortative mating. Overall, our findings indicate that the genetic divergence of each species is maintained in the presence of hybridization, as expected in a scenario of speciation-with-gene-flow.

Genetic features of bivalve transmissible neoplasia in blue mussels from the Kola Bay (Barents Sea) suggest a recent trans-Arctic migration of the cancer lineages

Ecology and biogeography of bivalve transmissible neoplasia (BTN) are underexplored due to its recent discovery and a challenging diagnostics. Blue mussels harbour two evolutionary lineages of BTN, MtrBTN1 and MtrBTN2, both derived from Mytilus trossulus. MtrBTN1 has been found only in M. trossulus from North Pacific. MtrBTN2 parasitizes different Mytilus spp. worldwide. BTN in M. trossulus in the Atlantic sector has never been studied. We looked for BTN in mussels from the Barents Sea using flow cytometry of cells, qPCR with primers specific to cancer-associated alleles and sequencing of mtDNA and nuclear loci. Both MtrBTN1 and MtrBTN2 were present in our material, though their prevalence was low (~0.4%). All cancers parasitized M. trossulus except one, MtrBTN1, which was found in a hybrid between M. trossulus and M. edulis. The mtDNA haplotypes found in both lineages were nearly identical to those known from the Northwest Pacific but not from elsewhere. Our results suggest that these two lineages may have arrived in the Barents Sea in recent decades with the maritime transport along the Northern Sea Route. A young evolutionary age of MtrBTN1 seems to indicate that it is an emerging disease in the process of niche expansion. Comparing the new and the published sequence data on tumour suppressor p53, we proved that the prevalence of BTN in mussels can reach epizootic levels. The finding of diverse recombinants between paternally and maternally inherited mtDNAs in somatic tissues of M. trossulus was an unexpected result of our study.

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.

Some Examples of the Use of Molecular Markers for Needs of Basic Biology and Modern Society

Application of molecular genetic markers appeared to be very fruitful in achieving many goals, including (i) proving the theoretic basements of general biology and (ii) assessment of worldwide biodiversity. Both are provided in the present meta-analysis and a review as the main signal. One of the basic current challenges in modern biology in the face of new demands in the 21st century is the validation of its paradigms such as the synthetic theory of evolution (STE) and biological species concept (BSC). Another of most valuable goals is the biodiversity assessment for a variety of social needs including free web-based information resources about any living being, renovation of museum collections, nature conservation that recognized as a global project, iBOL, as well as resolving global trading problems such as false labeling of species specimens used as food, drug components, entertainment, etc. The main issues of the review are focused on animals and combine four items. (1) A combination of nDNA and mtDNA markers best suits the identification of hybrids and estimation of genetic introgression. (2) The available facts on nDNA and mtDNA diversity seemingly make introgression among many taxa obvious, although it is evident, that introgression may be quite restricted or asymmetric, thus, leaving at least the "source" taxon (taxa) intact. (3) If we consider sexually reproducing species in marine and terrestrial realms introgressed, as it is still evident in many cases, then we should recognize that the BSC, in view of the complete lack of gene flow among species, is inadequate because many zoological species are not biological ones yet. However, vast modern molecular data have proven that sooner or later they definitely become biological species. (4) An investigation into the fish taxa divergence using the BOLD database shows that most gene trees are basically monophyletic and interspecies reticulations are quite rare.

Nuclear microsatellites reveal population genetic structuring and fine-scale pattern of hybridization in the Japanese mantis shrimp Oratosquilla oratoria

The interplay between historical and contemporary processes can produce complex patterns of genetic differentiation in the marine realm. Recent mitochondrial and nuclear sequence analyses revealed cryptic speciation in the Japanese mantis shrimp Oratosquilla oratoria. Herein, we applied nuclear microsatellite markers to examine patterns and causes of genetic differentiation in this morphotaxon. Population structure analyses revealed two genetically divergent and geographically structured clades in O. oratoria, one dominating the temperate zone of the Northwestern (NW) Pacific and the other occurring in the subtropical and tropical waters where are influenced by the Kuroshio Current. Two sympatric zones, one around the Changjiang Estuary in China coast and the other in the northern Japan Sea, were demonstrated to be hybrid zones where introgressive hybridization occurred asymmetrically. The interaction between historical climate shifts and contemporary factors (e.g., freshwater discharge, temperature gradient and isolation by distance) may contribute to the present-day genetic architecture in the Japanese mantis shrimp. Range shift induced by climate changes and oceanographic factors may promote hybridization and gene flow between the O. oratoria complex. Our results provide insights into the interacting mechanisms that give rise to diversification and speciation of coastal species in the NW Pacific.

Genotyping by sequencing can reveal the complex mosaic genomes in gene pools resulting from reticulate evolution: a case study in diploid and polyploid citrus

BACKGROUND AND AIMS

Reticulate evolution, coupled with reproductive features limiting further interspecific recombinations, results in admixed mosaics of large genomic fragments from the ancestral taxa. Whole-genome sequencing (WGS) data are powerful tools to decipher such complex genomes but still too costly to be used for large populations. The aim of this work was to develop an approach to infer phylogenomic structures in diploid, triploid and tetraploid individuals from sequencing data in reduced genome complexity libraries. The approach was applied to the cultivated Citrus gene pool resulting from reticulate evolution involving four ancestral taxa, C. maxima, C. medica, C. micrantha and C. reticulata.

METHODS

A genotyping by sequencing library was established with the restriction enzyme ApeKI applying one base (A) selection. Diagnostic single nucleotide polymorphisms (DSNPs) for the four ancestral taxa were mined in 29 representative varieties. A generic pipeline based on a maximum likelihood analysis of the number of read data was established to infer ancestral contributions along the genome of diploid, triploid and tetraploid individuals. The pipeline was applied to 48 diploid, four triploid and one tetraploid citrus accessions.

KEY RESULTS

Among 43 598 mined SNPs, we identified a set of 15 946 DSNPs covering the whole genome with a distribution similar to that of gene sequences. The set efficiently inferred the phylogenomic karyotype of the 53 analysed accessions, providing patterns for common accessions very close to that previously established using WGS data. The complex phylogenomic karyotypes of 21 cultivated citrus, including bergamot, triploid and tetraploid limes, were revealed for the first time.

CONCLUSIONS

The pipeline, available online, efficiently inferred the phylogenomic structures of diploid, triploid and tetraploid citrus. It will be useful for any species whose reproductive behaviour resulted in an interspecific mosaic of large genomic fragments. It can also be used for the first generations of interspecific breeding schemes.

Species delimitation and evolutionary reconstruction within an integrative taxonomic framework: A case study on Rhinolophus macrotis complex (Chiroptera: Rhinolophidae)

Species delimitation and evolutionary reconstruction remain challenging for non-model species that have experienced reticulate evolution and exhibit conflicting patterns of differentiation among multiple lines of evidence, such as mitochondrial and nuclear data and phenotypes. Here, we applied an integrative taxonomic approach to a case study of Rhinolophus macrotis complex, whose taxonomic status remains controversial, to provide insight into the systematics and evolutionary history of these species. By integrating traditional genetic markers with different modes of inheritance, genome-wide SNPs as well as phenotypic characteristics, we clarified the presence of three closely related species, R. episcopus, R. siamensis, and R. osgoodi, within this complex, and proposed a new taxonomic treatment for R. osgoodi. Our results suggested that hybridization and introgression are the main causes of low mtDNA divergence in these species. Combined with the demographic inference, we deduced that glacial-interglacial cycles drove geographic isolation and secondary contacts of these species, then promoted hybridization and lineage fusion among them, finally resulting in a reticulate evolutionary pattern. Overall, our study highlights the importance of combining multiple types of data to delimit species, especially those with conserved morphology, and to reveal the sophisticated processes of speciation.

Reticulate evolution in eukaryotes: Origin and evolution of the nitrate assimilation pathway

Genes and genomes can evolve through interchanging genetic material, this leading to reticular evolutionary patterns. However, the importance of reticulate evolution in eukaryotes, and in particular of horizontal gene transfer (HGT), remains controversial. Given that metabolic pathways with taxonomically-patchy distributions can be indicative of HGT events, the eukaryotic nitrate assimilation pathway is an ideal object of investigation, as previous results revealed a patchy distribution and suggested that the nitrate assimilation cluster of dikaryotic fungi (Opisthokonta) could have been originated and transferred from a lineage leading to Oomycota (Stramenopiles). We studied the origin and evolution of this pathway through both multi-scale bioinformatic and experimental approaches. Our taxon-rich genomic screening shows that nitrate assimilation is present in more lineages than previously reported, although being restricted to autotrophs and osmotrophs. The phylogenies indicate a pervasive role of HGT, with three bacterial transfers contributing to the pathway origin, and at least seven well-supported transfers between eukaryotes. In particular, we propose a distinct and more complex HGT path between Opisthokonta and Stramenopiles than the one previously suggested, involving at least two transfers of a nitrate assimilation gene cluster. We also found that gene fusion played an essential role in this evolutionary history, underlying the origin of the canonical eukaryotic nitrate reductase, and of a chimeric nitrate reductase in Ichthyosporea (Opisthokonta). We show that the ichthyosporean pathway, including this novel nitrate reductase, is physiologically active and transcriptionally co-regulated, responding to different nitrogen sources; similarly to distant eukaryotes with independent HGT-acquisitions of the pathway. This indicates that this pattern of transcriptional control evolved convergently in eukaryotes, favoring the proper integration of the pathway in the metabolic landscape. Our results highlight the importance of reticulate evolution in eukaryotes, by showing the crucial contribution of HGT and gene fusion in the evolutionary history of the nitrate assimilation pathway.

One of the most relevant findings in evolution was that lineages, either genes or genomes, can evolve through interchanging genetic material. For example, exon shuffling can lead to genes with complete novel functions, and genomes can acquire novel functionalities by means of horizontal gene transfer (HGT). Whereas HGT is known to be an important driver of metabolic remodelling and ecological adaptations in Bacteria, its importance and prevalence in eukaryotes remains controversial. We show that HGT played a major role in the origin and evolution of the eukaryotic nitrate assimilation pathway, with several bacteria-to-eukaryote and eukaryote-to-eukaryote transfers promoting the acquisition of this ecologically-relevant pathway to autotrophs and to distinct groups of osmotrophs. Moreover, we also show that gene fusion was important in this evolutionary history, underlying the origin of the canonical eukaryotic nitrate reductase, but also of a non-canonical nitrate reductase that we describe in Ichthyosporea, a poorly-characterized eukaryotic group that includes many parasitic species. In conclusion, our results highlight the importance of reticulate evolution in eukaryotes, by showing the contribution of HGT and gene fusion in the evolutionary history of the nitrate assimilation pathway.

Natural re-colonization and admixture of wolves (Canis lupus) in the US Pacific Northwest: challenges for the protection and management of rare and endangered taxa

Admixture resulting from natural dispersal processes can potentially generate novel phenotypic variation that may facilitate persistence in changing environments or result in the loss of population-specific adaptations. Yet, under the US Endangered Species Act, policy is limited for management of individuals whose ancestry includes a protected taxon; therefore, they are generally not protected under the Act. This issue is exemplified by the recently re-established grey wolves of the Pacific Northwest states of Washington and Oregon, USA. This population was likely founded by two phenotypically and genetically distinct wolf ecotypes: Northern Rocky Mountain (NRM) forest and coastal rainforest. The latter is considered potentially threatened in southeast Alaska and thus the source of migrants may affect plans for their protection. To assess the genetic source of the re-established population, we sequenced a ~ 300 bp portion of the mitochondrial control region and ~ 5 Mbp of the nuclear genome. Genetic analysis revealed that the Washington wolves share ancestry with both wolf ecotypes, whereas the Oregon population shares ancestry with NRM forest wolves only. Using ecological niche modelling, we found that the Pacific Northwest states contain environments suitable for each ecotype, with wolf packs established in both environmental types. Continued migration from coastal rainforest and NRM forest source populations may increase the genetic diversity of the Pacific Northwest population. However, this admixed population challenges traditional management regimes given that admixture occurs between an adaptively distinct ecotype and a more abundant reintroduced interior form. Our results emphasize the need for a more precise US policy to address the general problem of admixture in the management of endangered species, subspecies, and distinct population segments.

Analysis of sequence divergence in redfin (Cypriniformes, Cyprinidae, Tribolodon) based on mtDNA and nDNA markers with inferences in systematics and genetics of speciation

To clarify relationship of species of the genus Tribolodon in the Russian part of their distribution ranges, two mitochondrial markers (Co-1 and Cyt-b), a nuclear marker (Rho), and a gene marker of rDNA internal transcribed spacer (ITS-1,2) were used. Depending on the marker, different numbers of species groups were detected by the ABGD method, but in combination with the analysis of phylograms, these data generally support the known species clusters and regional intraspecies groups. A complex analysis of sequences from three redfin species within the area of the study, based on four marker genes and using the methods of molecular phylogenetics, ordination of genetic distances, recombinant analysis, and population genetic approaches, has revealed clusters of three commonly recognized species, regional intraspecific groups or individuals of local populations, and few hybrid individuals. DNA barcoding technique proved to be efficient with the use of two mtDNA markers: Co-1 and Cyt-b. It has been found that analysis of insertions and substitutions within the ITS-1,2 gene marker is also suitable for identification of Tribolodon species. Results of the studies of local groups do not confirm a sufficient level of differences for defining any new taxa of a species rank in the genus Tribolodon.

Barcode index number, taxonomic rank and modes of speciation: examples from fish

Species delimitation by DNA sequence data or DNA barcoding is successful, as confirmed by the vast BOLD data base. However, the theory that would explain this fact has not been developed yet. An approach based on Barcoding Index Number (BIN), suggested in the assignment, allows delimiting of taxa of three ranks (species, genera, and families) and statistical validation with a high precision of delimiting (over 80%), as well as shows for majority of Co-1-based single gene trees good correspondence between their topology and conventional taxa content for analyzed fish species (R² ≈ 0.84-0.98). Knowledge of deviations from these data can help to find out new taxa and improve biodiversity description. It is concluded that delimiting is successful for bulk of cases because the geographic mode of speciation prevails in nature. It takes a long time for new taxa to form in isolation, which allows accumulation of random mutations and many different nucleotide substitutions between them that can be detected by molecular markers and give unique DNA barcodes. The use of BIN approach, described here, can aid greatly in making this important question clearer especially under wider examination of other organisms.

Cryptic diversity, geographical endemism and allopolyploidy in NE Pacific seaweeds

Background

Molecular markers are revealing a much more diverse and evolutionarily complex picture of marine biodiversity than previously anticipated. Cryptic and/or endemic marine species are continually being found throughout the world oceans, predominantly in inconspicuous tropical groups but also in larger, canopy-forming taxa from well studied temperate regions. Interspecific hybridization has also been found to be prevalent in many marine groups, for instance within dense congeneric assemblages, with introgressive gene-flow being the most common outcome. Here, using a congeneric phylogeographic approach, we investigated two monotypic and geographically complementary sister genera of north-east Pacific intertidal seaweeds (Hesperophycus and Pelvetiopsis), for which preliminary molecular tests revealed unexpected conflicts consistent with unrecognized cryptic diversity and hybridization.

Results

The three recovered mtDNA clades did not match a priori species delimitations. H. californicus was congruent, whereas widespread P. limitata encompassed two additional narrow-endemic species from California - P. arborescens (here genetically confirmed) and P. hybrida sp. nov. The congruence between the genotypic clusters and the mtDNA clades was absolute. Fixed heterozygosity was apparent in a high proportion of loci in P. limitata and P. hybrida, with genetic analyses showing that the latter was composed of both H. californicus and P. arborescens genomes. All four inferred species could be distinguished based on their general morphology.

Conclusions

This study confirmed additional diversity and reticulation within NE Pacific Hesperophycus/Pelvetiopsis, including the validity of the much endangered, modern climatic relict P. arborescens, and the identification of a new, stable allopolyploid species (P. hybrida) with clearly discernable ancestry (♀ H. californicus x ♂ P. arborescens), morphology, and geographical distribution. Allopolyploid speciation is otherwise completely unknown in brown seaweeds, and its unique occurrence within this genus (P. limitata possibly representing a second example) remains enigmatic. The taxonomic separation of Hesperophycus and Pelvetiopsis is not supported and the genera should be synonymized; we retain only the latter. The transitional coastline between Point Conception and Monterey Bay represented a diversity hotspot for the genus and the likely sites of extraordinary evolutionary events of allopolyploid speciation at sympatric range contact zones. This study pinpoints how much diversity (and evolutionary processes) potentially remains undiscovered even on a conspicuous seaweed genus from the well-studied Californian intertidal shores let alone in other, less studied marine groups and regions/depths.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0878-2) contains supplementary material, which is available to authorized users.

Introgression and selection shaped the evolutionary history of sympatric sister-species of coral reef fishes (genus: Haemulon)

Closely related marine species with large overlapping ranges provide opportunities to study mechanisms of speciation, particularly when there is evidence of gene flow between such lineages. Here, we focus on a case of hybridization between the sympatric sister-species Haemulon maculicauda and H. flaviguttatum, using Sanger sequencing of mitochondrial and nuclear loci, as well as 2422 single nucleotide polymorphisms (SNPs) obtained via restriction site-associated DNA sequencing (RADSeq). Mitochondrial markers revealed a shared haplotype for COI and low divergence for CytB and CR between the sister-species. On the other hand, complete lineage sorting was observed at the nuclear loci and most of the SNPs. Under neutral expectations, the smaller effective population size of mtDNA should lead to fixation of mutations faster than nDNA. Thus, these results suggest that hybridization in the recent past (0.174-0.263 Ma) led to introgression of the mtDNA, with little effect on the nuclear genome. Analyses of the SNP data revealed 28 loci potentially under divergent selection between the two species. The combination of mtDNA introgression and limited nuclear DNA introgression provides a mechanism for the evolution of independent lineages despite recurrent hybridization events. This study adds to the growing body of research that exemplifies how genetic divergence can be maintained in the presence of gene flow between closely related species.

Studying the genetic basis of speciation in high gene flow marine invertebrates

A growing number of genes responsible for reproductive incompatibilities between species (barrier loci) exhibit the signals of positive selection. However, the possibility that genes experiencing positive selection diverge early in speciation and commonly cause reproductive incompatibilities has not been systematically investigated on a genome-wide scale. Here, I outline a research program for studying the genetic basis of speciation in broadcast spawning marine invertebrates that uses a priori genome-wide information on a large, unbiased sample of genes tested for positive selection. A targeted sequence capture approach is proposed that scores single-nucleotide polymorphisms (SNPs) in widely separated species populations at an early stage of allopatric divergence. The targeted capture of both coding and non-coding sequences enables SNPs to be characterized at known locations across the genome and at genes with known selective or neutral histories. The neutral coding and non-coding SNPs provide robust background distributions for identifying FST-outliers within genes that can, in principle, identify specific mutations experiencing diversifying selection. If natural hybridization occurs between species, the neutral coding and non-coding SNPs can provide a neutral admixture model for genomic clines analyses aimed at finding genes exhibiting strong blocks to introgression. Strongylocentrotid sea urchins are used as a model system to outline the approach but it can be used for any group that has a complete reference genome available.

Molecular systematics and DNA barcoding of Altai osmans, oreoleuciscus (pisces, cyprinidae, and leuciscinae), and their nearest relatives, inferred from sequences of cytochrome b (Cyt-b), cytochrome oxidase c (Co-1), and complete mitochondrial genome

Mitochondrial DNA (mtDNA) at the protein-coding Cyt-b gene along with data retrieved from GenBank for Co-1 gene fragments and complete mitochondrial genome (mitogenome) of Altai osmans and the nearest relatives of Leuciscinae fish species were compared for the estimation of variability and phylogenetic tree building. Phylogenetic trees were built by four techniques: Bayesian (BA), maximum likelihood (ML), maximum parsimony (MP), and neighbor-joining (NJ). Resolution of Cyt-b trees for species of two genera (Oreoleuciscus and Phoxinus) was quite distinct at all the approaches. For Tribolodon, the single gene trees were not well resolved; however, the mitogenome tree was resolved. Species identification on per individual basis (DNA barcoding) was high for both Cyt-b and Co-1 genes. The trees built using the data for 13 protein mitochondrial genes revealed a complicated phylogenetic pattern within the subfamily Leuciscinae. Scores of the average p-distances at three taxonomic levels were considerably different: (1) 1.16 ± 0.96, (2) 8.21 ± 1.01, and (3) 16.41 ± 0.85 for Cyt-b and (1) 1.04 ± 0.78, (2) 8.30 ± 0.92, and (3) 10.74 ± 0.79 for 13 protein genes of mitogenome, where (1) is intraspecies, (2) is intragenus, and (3) is intrasubfamily levels. Data on mitogenome distances were summarized for the taxonomic hierarchy for the first time. A concordant increase in distance score with growth of the rank of taxa (having the minimum score at the intraspecies level), both for a single gene and the whole mitogenome, substantiates the concept that speciation in the subfamily Leuciscinae in most cases follows the geographic mode. The distinct clustering of Altai osmans, Oreoleuciscus potanini and O. humilis, in the Cyt-b and Co-1 gene trees with small overall genetic distances, obtained for both genes, allows us to consider these taxa as separate but genetically sister species.

Interspecific plastidial recombination in the diatom genus Pseudo-nitzschia

Plastids are usually uni-parentally inherited and genetic recombination between these organelles is seldom observed. The genus Pseudo-nitzschia, a globally relevant marine diatom, features bi-parental plastid inheritance in the course of sexual reproduction. This observation inspired the recombination detection we pursued in this paper over a ~1,400-nucleotide-long region of the plastidial rbcL, a marker used in both molecular taxonomy and phylogenetic studies in diatoms. Among all the rbcL-sequences available in web-databases for Pseudo-nitzschia, 42 haplotypes were identified and grouped in five clusters by Bayesian phylogeny. Signs of hybridization were evident in four of five clusters, at both intra- and interspecific levels, suggesting that, in diatoms, (i) plastidial recombination is not absent and (ii) hybridization can play a role in speciation of Pseudo-nitzschia spp.

Possible natural hybridization of two morphologically distinct species of Acropora (Cnidaria, Scleractinia) in the Pacific: fertilization and larval survival rates

Natural hybridization of corals in the Indo-Pacific has been considered rather rare. However, field studies have observed many corals with intermediate interspecific or unusual morphologies. Given that the existence of F1 hybrids with intermediate interspecific morphologies has been proven in the Caribbean, hybrids may also inhabit the Indo-Pacific and occur more frequently than expected. In this study, we focused on two morphologically different species, Acropora florida and A. intermedia, and performed crossing experiments at Akajima Island, Japan. Results showed that these species could hybridize in both directions via eggs and sperm, but that fertilization rates significantly differed according to which species provided eggs. These results are similar to those reported from the Caribbean. Although all embryos developed normally to the planular larval stage, the developmental processes of some hybrid embryos were delayed by approximately 1 h compared with conspecific embryos, suggesting that fertilization occurred 1 h later in interspecific crosses than in intraspecific crosses. More successful hybridization could occur under conditions with low numbers of conspecific colonies. Additionally, a comparison of survival rates between hybrid and intraspecific larvae revealed that intra- and interspecific larvae produced from eggs of A. florida survived for significantly longer than those produced from eggs of A. intermedia. Considering these data, under specific conditions, hybrids can be expected to be produced and survive in nature in the Pacific. Furthermore, we identified one colony with intermediate morphology between A. florida and A. intermedia in the field. This colony was fertilized only by eggs of A. florida, with high fertilization rates, suggesting that this colony would be a hybrid of these two species and might be backcrossed.

DNA variation in the phenotypically-diverse brown alga Saccharina japonica

BACKGROUND

Saccharina japonica (Areschoug) Lane, Mayes, Druehl et Saunders is an economically important and highly morphologically variable brown alga inhabiting the northwest Pacific marine waters. On the basis of nuclear (ITS), plastid (rbcLS) and mitochondrial (COI) DNA sequence data, we have analyzed the genetic composition of typical Saccharina japonica (TYP) and its two common morphological varieties, known as the "longipes" (LON) and "shallow-water" (SHA) forms seeking to clarify their taxonomical status and to evaluate the possibility of cryptic species within S. japonica.

RESULTS

The data show that the TYP and LON forms are very similar genetically in spite of drastic differences in morphology, life history traits, and ecological preferences. Both, however, are genetically quite different from the SHA form. The two Saccharina lineages are distinguished by 109 fixed single nucleotide differences as well as by seven fixed length polymorphisms (based on a 4,286 bp concatenated dataset that includes three gene regions). The GenBank database reveals a close affinity of the TYP and LON forms to S. japonica and the SHA form to S. cichorioides. The three gene markers used in the present work have different sensitivity for the algal species identification. COI gene was the most discriminant gene marker. However, we have detected instances of interspecific COI recombination reflecting putative historical hybridization events between distantly related algal lineages. The recombinant sequences show highly contrasted level of divergence in the 5'- and 3'- regions of the gene, leading to significantly different tree topologies depending on the gene segment (5'- or 3'-) used for tree reconstruction. Consequently, the 5'-COI "barcoding" region (~ 650 bp) can be misleading for identification purposes, at least in the case of algal species that might have experienced historical hybridization events.

CONCLUSION

Taking into account the potential roles of phenotypic plasticity in evolution, we conclude that the TYP and LON forms represent examples of algae phenotypic diversification that enables successful adaptation to contrasting shallow- and deep-water marine environments, while the SHA form is very similar to S. cichorioides and should be considered a different species. Practical applications for algal management and conservation are briefly considered.

Weak prezygotic isolating mechanisms in threatened Caribbean Acropora corals

The Caribbean corals, Acropora palmata and A. cervicornis, recently have undergone drastic declines primarily as a result of disease. Previous molecular studies have demonstrated that these species form a hybrid (A. prolifera) that varies in abundance throughout the range of the parental distribution. There is variable unidirectional introgression across loci and sites of A. palmata genes flowing into A. cervicornis. Here we examine the efficacy of prezygotic reproductive isolating mechanisms within these corals including spawning times and choice and no-choice fertilization crosses. We show that these species have subtly different mean but overlapping spawning times, suggesting that temporal isolation is likely not an effective barrier to hybridization. We found species-specific differences in gametic incompatibilities. Acropora palmata eggs were relatively resistant to hybridization, especially when conspecific sperm are available to outcompete heterospecific sperm. Acropora cervicornis eggs demonstrated no evidence for gametic incompatibility and no evidence of reduced viability after aging four hours. This asymmetry in compatibility matches previous genetic data on unidirectional introgression.

The deep-sea microfossil record of macroevolutionary change in plankton and its study

The deep-sea planktonic microfossil record (foraminifera, coccolithophores, diatoms, radiolaria and dinoflagellates) provides a unique resource for palaeobiology. Despite some geographical gaps due to poor regional preservation, and intermittant time intervals lost to erosion, most time periods for each Cenozoic planktonic biogeographical province are preserved. Vast numbers of specimens and numerous deep-sea cores provide abundant material and the opportunity to tightly integrate macroevolutionary and palaeoenvironmental data. Current documentation of this record is mixed. Catalogues for foraminifera and coccolithophores offer nearly complete species-level clade histories, but taxonomy for siliceous microfossils is incomplete. Published occurrence data is primarily stratigraphic and covers only a fraction of the total preserved diversity. Age models for some sections are excellent (accuracy c. 100 kya) but for many other sections are still poor. Taxonomic errors, age model errors and reworking displace fossil occurrences in time, complicating palaeobiological analysis. With additional taxonomic work, careful collection of whole fauna/floral assemblage occurrence data, improved age models, and the development of better data filtering and analysis tools to deal with data outliers the deep-sea microfossil record can deliver its promise of providing the most complete, detailed record of macroevolutionary change available to science.

Supplementary material:

Supplementary Appendix is available at http://www.geolsoc.org.uk/SUP18485

Reticulate evolution: frequent introgressive hybridization among Chinese hares (genus lepus) revealed by analyses of multiple mitochondrial and nuclear DNA loci

BACKGROUND

Interspecific hybridization may lead to the introgression of genes and genomes across species barriers and contribute to a reticulate evolutionary pattern and thus taxonomic uncertainties. Since several previous studies have demonstrated that introgressive hybridization has occurred among some species within Lepus, therefore it is possible that introgressive hybridization events also occur among Chinese Lepus species and contribute to the current taxonomic confusion.

RESULTS

Data from four mtDNA genes, from 116 individuals, and one nuclear gene, from 119 individuals, provides the first evidence of frequent introgression events via historical and recent interspecific hybridizations among six Chinese Lepus species. Remarkably, the mtDNA of L. mandshuricus was completely replaced by mtDNA from L. timidus and L. sinensis. Analysis of the nuclear DNA sequence revealed a high proportion of heterozygous genotypes containing alleles from two divergent clades and that several haplotypes were shared among species, suggesting repeated and recent introgression. Furthermore, results from the present analyses suggest that Chinese hares belong to eight species.

CONCLUSION

This study provides a framework for understanding the patterns of speciation and the taxonomy of this clade. The existence of morphological intermediates and atypical mitochondrial gene genealogies resulting from frequent hybridization events likely contribute to the current taxonomic confusion of Chinese hares. The present study also demonstrated that nuclear gene sequence could offer a powerful complementary data set with mtDNA in tracing a complete evolutionary history of recently diverged species.