Interspecific hybridization causes long-term phylogenetic discordance between nuclear and mitochondrial genomes in freshwater fishes

Comparative mitogenomics of marine angelfishes (F: Pomacanthidae)

The targeted capture of ultraconserved elements (UCEs) has substantially increased the amount of genetic data available for phylogenomic reconstructions. These capture datasets frequently contain mitochondrial DNA as a by-product, often in the form of complete mitogenomes. These can be efficiently harvested to expand existing datasets without additional costs. Here, we present new mitochondrial genomes for six marine angelfish species (F: Pomacanthidae), assembled and annotated from off-target UCE reads. We provide the first comparative analysis of all mitochondrial genomes available for the Pomacanthidae. Results showed that the average length of pomacanthid mitogenomes is 16.8 kbp. Total GC and AT content varied between 44.5% and 46.3%, and 53.7% and 55.5%, respectively. The architecture of angelfish mitogenomes was comparable to that seen in other fish species with 13 protein-coding genes (PCGs), 22 transfer RNA genes, two ribosomal RNA genes and the control region. All 13 PCGs evolved under purifying selection, highlighting a high level of selection pressure and gene expression to preserve genetic integrity. The ND6 and ATP8 genes had the highest ratio of non-synonymous (dN) to synonymous (dS) substitutions, indicating a relaxation of purifying selection constraints. Finally, these newly assembled mitogenomes will allow further investigations of the population genetics, systematics and evolutionary biology of one of the most prominent reef fish family in the aquarium trade.

Toward telomere-to-telomere cat genomes for precision medicine and conservation biology

Genomic data from species of the cat family Felidae promise to stimulate veterinary and human medical advances, and clarify the coherence of genome organization. We describe how interspecies hybrids have been instrumental in the genetic analysis of cats, from the first genetic maps to propelling cat genomes toward the T2T standard set by the human genome project. Genotype-to-phenotype mapping in cat models has revealed dozens of health-related genetic variants, the molecular basis for mammalian pigmentation and patterning, and species-specific adaptations. Improved genomic surveillance of natural and captive populations across the cat family tree will increase our understanding of the genetic architecture of traits, population dynamics, and guide a future of genome-enabled biodiversity conservation.

A comprehensive phylogenomic study unveils evolutionary patterns and challenges in the mitochondrial genomes of Carcharhiniformes: A focus on Triakidae

The complex evolutionary patterns in the mitochondrial genome (mitogenome) of the most species-rich shark order, the Carcharhiniformes (ground sharks) has led to challenges in the phylogenomic reconstruction of the families and genera belonging to the order, particularly the family Triakidae (houndsharks). The current state of Triakidae phylogeny remains controversial, with arguments for both monophyly and paraphyly within the family. We hypothesize that this variability is triggered by the selection of different a priori partitioning schemes to account for site and gene heterogeneity within the mitogenome. Here we used an extensive statistical framework to select the a priori partitioning scheme for inference of the mitochondrial phylogenomic relationships within Carcharhiniformes, tested site heterogeneous CAT + GTR + G4 models and incorporated the multi-species coalescent model (MSCM) into our analyses to account for the influence of gene tree discordance on species tree inference. We included five newly assembled houndshark mitogenomes to increase resolution of Triakidae. During the assembly procedure, we uncovered a 714 bp-duplication in the mitogenome of Galeorhinus galeus. Phylogenetic reconstruction confirmed monophyly within Triakidae and the existence of two distinct clades of the expanded Mustelus genus. The latter alludes to potential evolutionary reversal of reproductive mode from placental to aplacental, suggesting that reproductive mode has played a role in the trajectory of adaptive divergence. These new sequences have the potential to contribute to population genomic investigations, species phylogeography delineation, environmental DNA metabarcoding databases and, ultimately, improved conservation strategies for these ecologically and economically important species.

Fish biogeography and hybridization: do contemporary distributions predict introgression history?

Freshwater ecosystems frequently house diverse assemblages of closely related fish taxa, which can be particularly prone to hybridization and introgression. While extensive introgression may be expected among biogeographically proximate lineages, recent analyses imply that contemporary distributions do not always accurately predict hybridization history. Here, we use the ABBA-BABA approach to test biogeographic hypotheses regarding the extent of hybridization in the recent evolution of New Zealand's species-rich freshwater Galaxias vulgaris fish complex. Genome-wide comparisons reveal significant increases in introgression associated with increasing geographic overlap of taxa. The estimator DP, which assesses the net proportion of a genome originating from introgression, shows a particularly strong relationship with biogeographic overlap (R2 = .43; p = .005). Our analyses nevertheless reveal surprisingly substantial signatures of introgression among taxa that currently have disjunct distributions within drainages (e.g., separate subcatchments). These "anomalies" imply that current biogeography is not always an accurate predictor of introgression history. Our study suggests that both modern and ancient biogeographic shifts, including recent anthropogenic range fragmentation and tectonically driven riven capture events, have influenced introgression histories in this dynamic freshwater fish radiation.

A community genomics approach to natural hybridization

Hybridization is a complicated, oft-misunderstood process. Once deemed unnatural and uncommon, hybridization is now recognized as ubiquitous among species. But hybridization rates within and among communities are poorly understood despite the relevance to ecology, evolution and conservation. To clarify, we examined hybridization across 75 freshwater fish communities within the Ozarks of the North American Interior Highlands (USA) by single nucleotide polymorphism (SNP) genotyping 33 species (N = 2865 individuals; double-digest restriction site-associated DNA sequencing (ddRAD)). We found evidence of hybridization (70 putative hybrids; 2.4% of individuals) among 18 species-pairs involving 73% (24/33) of study species, with the majority being concentrated within one family (Leuciscidae/minnows; 15 species; 66 hybrids). Interspecific genetic exchange-or introgression-was evident from 24 backcrossed individuals (10/18 species-pairs). Hybrids occurred within 42 of 75 communities (56%). Four selected environmental variables (species richness, protected area extent, precipitation (May and annually)) exhibited 73-78% accuracy in predicting hybrid occurrence via random forest classification. Our community-level assessment identified hybridization as spatially widespread and environmentally dependent (albeit predominantly within one diverse, omnipresent family). Our approach provides a more holistic survey of natural hybridization by testing a wide range of species-pairs, thus contrasting with more conventional evaluations.

Different waves of postglacial recolonisation and genomic structure of bank vole populations in NE Poland

Previous studies indicated that in some species phylogeographic patterns obtained in the analysis of nuclear and mitochondrial DNA (mtDNA) markers can be different. Such mitonuclear discordance can have important evolutionary and ecological consequences. In the present study, we aimed to check whether there was any discordance between mtDNA and nuclear DNA in the bank vole population in the contact zone of its two mtDNA lineages. We analysed the population genetic structure of bank voles using genome-wide genetic data (SNPs) and diversity of sequenced heart transcriptomes obtained from selected individuals from three populations inhabiting areas outside the contact zone. The SNP genetic structure of the populations confirmed the presence of at least two genetic clusters, and such division was concordant with the patterns obtained in the analysis of other genetic markers and functional genes. However, genome-wide SNP analyses revealed the more detailed structure of the studied population, consistent with more than two bank vole recolonisation waves, as recognised previously in the study area. We did not find any significant differences between individuals representing two separate mtDNA lineages of the species in functional genes coding for protein-forming complexes, which are involved in the process of cell respiration in mitochondria. We concluded that the contemporary genetic structure of the populations and the width of the contact zone were shaped by climatic and environmental factors rather than by genetic barriers. The studied populations were likely isolated in separate Last Glacial Maximum refugia for insufficient amount of time to develop significant genetic differentiation.

Evaluating a species phylogeny using ddRAD SNPs: Cyto-nuclear discordance and introgression in the salmonid genus Thymallus (Salmonidae)

Hybridization and introgression are very common among freshwater fishes due to the dynamic nature of hydrological landscapes. Cyclic patterns of allopatry and secondary contact provide numerous opportunities for interspecific gene flow, which can lead to discordant paths of evolution for mitochondrial and nuclear genomes. Here, we used double digest restriction-site associated DNA sequencing (ddRADseq) to obtain a genome-wide single nucleotide polymorphism (SNP) dataset comprehensive for allThymallus (Salmonidae)species to infer phylogenetic relationships and evaluate potential recent and historical gene flow among species. The newly obtained nuclear phylogeny was largely concordant with a previously published mitogenome-based topology but revealed a few cyto-nuclear discordances. These incongruencies primarily involved the placement of internal nodes rather than the resolution of species, except for one European species where anthropogenic stock transfers are thought to be responsible for the observed pattern. The analysis of four contact zones where multiple species are found revealed a few cases of mitochondrial capture and limited signals of nuclear introgression. Interestingly, the mechanisms restricting interspecific gene flow might be distinct; while in zones of secondary contact, small-scale physical habitat separation appeared as a limiting factor, biologically based reinforcement mechanisms are presumed to be operative in areas where species presumably evolved in sympatry. Signals of historical introgression were largely congruent with the routes of species dispersal previously inferred from mitogenome data. Overall, the ddRADseq dataset provided a robust phylogenetic reconstruction of the genus Thymallus including new insights into historical hybridization and introgression, opening up new questions concerning their evolutionary history.

Intergeneric hybridization of two stickleback species leads to introgression of membrane-associated genes and invasive TE expansion

Interspecific hybridization has occurred relatively frequently during the evolution of vertebrates. This process usually abolishes reproductive isolation between the parental species. Moreover, it results in the exchange of genetic material and can lead to hybridogenic speciation. Hybridization between species has predominately been observed at the interspecific level, whereas intergeneric hybridization is rarer. Here, using whole-genome sequencing analysis, we describe clear and reliable signals of intergeneric introgression between the three-spined stickleback (Gasterosteus aculeatus) and its distant mostly freshwater relative the nine-spined stickleback (Pungitius pungitius) that inhabit northwestern Russia. Through comparative analysis, we demonstrate that such introgression phenomena apparently take place in the moderate-salinity White Sea basin, although it is not detected in Japanese sea stickleback populations. Bioinformatical analysis of the sites influenced by introgression showed that they are located near transposable elements, whereas those in protein-coding sequences are mostly found in membrane-associated and alternative splicing-related genes.

Endangered Schizothoracin Fish in the Tarim River Basin Are Threatened by Introgressive Hybridization

Simple Summary

Interspecific hybridization and introgression may threaten the survival of endangered species. Big-head Schizothoracin (Aspiorhynchus laticeps) and Tarim Schizothoracin (Schizothorax biddulphi) are two recognized types of fish of Schizothoracinae which belong to the family Cyprinidae. Big-head Schizothoracin and Tarim Schizothoracin are sympatrically distributed in the Tarim River basin in Xinjiang, China, and have been listed as first-class and second-class national key protected animals, respectively. Based on morphological characteristics, Schizothorax esocinus (another Schizothoracin that occurs in the Tarim River basin) is speculated to be hybrid offspring of Big-head Schizothoracin and Tarim Schizothoracin, but there is no direct genetic evidence for this point. In this study, the hybridization status of Schizothorax esocinus was confirmed by comparing genetic constitutions of mitochondrial DNA (mtDNA) and inter transcribed spacer (ITS2) encoded by the nuclear genome. Extensive hybridization and introgression between Big-head Schizothoracin and Tarim Schizothoracin was detected. Since both Big-head Schizothoracin and Tarim Schizothoracin are critically endangered fishes, risk management due to hybridization is recommended to be considered as a part of current conservation programs.

Abstract

Big-head Schizothoracin (Aspiorhynchus laticeps) and Tarim Schizothoracin (Schizothorax biddulphi) are locally sympatric in the Tarim River Basin. Although another Schizothoracin (Schizothorax esocinus) in Tarim River basin has been speculated to be hybrid offspring of Big-head Schizothoracin and Tarim Schizothoracin, there was no genetic evidence. Previous studies on the genetics and evolution of Schizothoracins in Xinjiang Province were mostly based on mitochondrial DNA (mtDNA), whose characteristics of maternal inheritance made it hard to answer the question of whether there was hybridization and introgression between Big-head Schizothoracin and Tarim Schizothoracin. In this study, cytochrome b (cytb) gene of mtDNA and internal transcribed spacer 2 (ITS2) that is encoded by the nuclear genome were genotyped within the entire samples at the same time. Our results confirmed that Schizothorax esocinus was the hybrid offspring of Big-head Schizothoracin and Tarim Schizothoracin. The heterozygous ITS2 genotypes and/or Aspiorhynchus laticeps-like mtDNA were also detected in a subset of samples that should have been identified as pure Schizothorax biddulphi based on morphology. The ITS2 is characterized by multi-copy, concert evolution, and biparental inheritance. Thus, by comparing with mtDNA data, broad-scale bidirectional hybridization and introgression between Big-head Schizothoracin and Tarim Schizothoracin were revealed. Although interspecific hybridization may play a positive role in ecology and evolution, interspecific hybrids could threaten their parental species by the swamping of genetics and demography. As both parents of hybridization are critically endangered fishes, in this case, it is urgently necessary to strengthen the scientific assessment of the risks of the hybrids and the control of the hybridization and introgression between Aspiorhynchus laticeps and Schizothorax biddulphi in the Tarim River Basin.

Integrating ecological and molecular data to investigate species maintenance and interspecific hybridization between the redbreast sunfish (Lepomis auritus) and the pumpkinseed sunfish (L. gibbosus)

True sunfishes (genus Lepomis) of North America are renowned for hybridization, although hybrids often comprise of only a fraction of the sunfish population in natural systems. Species boundaries may be maintained through premating mechanisms (e.g., ethological and habitat isolation), which may be facilitated by the promiscuous sunfish mating system; primarily, singular males guard nests which multiple females visit. Natural hybridization occurs between redbreast (L. auritus (Linnaeus, 1758)) and pumpkinseed (L. gibbosus (Linnaeus, 1758)) sunfishes, but there is no published research into its extent and causes. In this study, we assessed the relationship of ecological variables to the prevalence of hybridization in a sunfish population in a freshwater lake. Nests of both species were evaluated for differences in substrate, date, water depth, temperature, and dissolved oxygen. Sequencing of one fragment of nuclear DNA (calmodulin intron 4) and one fragment of mitochondrial DNA (12S) was employed on adult males and eight offspring from each nest to determine genetic identity and hybridization direction. Nonparametric tests found a significant difference in substrate between species’ nests (p < 0.01), warranting further investigation into nest substrate as a mechanism for species maintenance. Furthermore, we observed asymmetrical hybridization with the less abundant species (redbreast sunfish) disproportionately affected.

Genomic data and multi-species demographic modelling uncover past hybridization between currently allopatric freshwater species

Evidence for ancient interspecific gene flow through hybridization has been reported in many animal and plant taxa based on genetic markers. The study of genomic patterns of closely related species with allopatric distributions allows the assessment of the relative importance of vicariant isolating events and past gene flow. Here, we investigated the role of gene flow in the evolutionary history of four closely related freshwater fish species with currently allopatric distributions in western Iberian rivers-Squalius carolitertii, S. pyrenaicus, S. torgalensis and S. aradensis-using a population genomics dataset of 23,562 SNPs from 48 individuals, obtained through genotyping by sequencing (GBS). We uncovered a species tree with two well-differentiated clades: (i) S. carolitertii and S. pyrenaicus; and (ii) S. torgalensis and S. aradensis. By using D-statistics and demographic modelling based on the site frequency spectrum, comparing alternative demographic scenarios of hybrid origin, secondary contact and isolation, we found that the S. pyrenaicus North lineage is likely the result of an ancient hybridization event between S. carolitertii (contributing ~84%) and S. pyrenaicus South lineage (contributing ~16%), consistent with a hybrid speciation scenario. Furthermore, in the hybrid lineage, we identify outlier loci potentially affected by selection favouring genes from each parental lineage at different genomic regions. Our results suggest that ancient hybridization can affect speciation and that freshwater fish species currently in allopatry are useful to study these processes.

The roles of aridification and sea level changes in the diversification and persistence of freshwater fish lineages

While the influence of Pleistocene climatic changes on divergence and speciation has been well-documented across the globe, complex spatial interactions between hydrology and eustatics over longer timeframes may also determine species evolutionary trajectories. Within the Australian continent, glacial cycles were not associated with changes in ice cover and instead largely resulted in fluctuations from moist to arid conditions across the landscape. We investigated the role of hydrological and coastal topographic changes brought about by Plio-Pleistocene climatic changes on the biogeographic history of a small Australian freshwater fish, the southern pygmy perch Nannoperca australis. Using 7958 ddRAD-seq (double digest restriction-site associated DNA) loci and 45,104 filtered SNPs, we combined phylogenetic, coalescent and species distribution analyses to assess the various roles of aridification, sea level and tectonics and associated biogeographic changes across southeast Australia. Sea-level changes since the Pliocene and reduction or disappearance of large waterbodies throughout the Pleistocene were determining factors in strong divergence across the clade, including the initial formation and maintenance of a cryptic species, N. 'flindersi'. Isolated climatic refugia and fragmentation due to lack of connected waterways maintained the identity and divergence of inter- and intraspecific lineages. Our historical findings suggest that predicted increases in aridification and sea level due to anthropogenic climate change might result in markedly different demographic impacts, both spatially and across different landscape types.

Testing the utility of alternative metrics of branch support to address the ancient evolutionary radiation of tunas, stromateoids, and allies (Teleostei: Pelagiaria)

The use of high-throughput sequencing technologies to produce genome-scale datasets was expected to settle some long-standing controversies across the Tree of Life, particularly in areas where short branches occur at deep timescales. Instead, these datasets have often yielded many well-supported but conflicting topologies, and highly variable gene-tree distributions. A variety of branch-support metrics beyond the nonparametric bootstrap are now available to assess how robust a phylogenetic hypothesis may be, as well as new methods to quantify gene-tree discordance. We applied multiple branch support metrics to an ancient group of marine fishes (Teleostei: Pelagiaria) whose interfamilial relationships have proven difficult to resolve due to a rapid accumulation of lineages very early in its history. We analyzed hundreds of loci including published UCE data and newly generated exonic data along with their flanking regions to represent all 16 extant families for more than 150 out of 284 valid species in the group. Branch support was lower for interfamilial relationships (except the SH-like aLRT and aBayes methods) regardless of the type of marker used. Several nodes that were highly supported with bootstrap had very low site and gene-tree concordance, revealing underlying conflict. Despite this conflict, we were able to identify four consistent interfamilial clades, each comprised of two or three families. Combining exons with their flanking regions also produced increased branch lengths in the deep branches of the pelagiarian tree. Our results demonstrate the limitations of employing current metrics of branch support and species-tree estimation when assessing the confidence of ancient evolutionary radiations and emphasize the necessity to embrace alternative measurements to explore phylogenetic uncertainty and discordance in phylogenomic datasets.

Surprising Pseudogobius: Molecular systematics of benthic gobies reveals new insights into estuarine biodiversity (Teleostei: Gobiiformes)

Snubnose gobies (genus Pseudogobius: Gobionellinae) are ubiquitous to, and important components of, estuarine ecosystems of the Indo-west Pacific. These small benthic fishes occur in freshwater, brackish and marine habitats such as mangroves, sheltered tide pools and lowland streams, and represent a model group for understanding the biodiversity and biogeography of estuarine fauna. To develop the species-level framework required for a concurrent morphological taxonomic appraisal, we undertook thorough sampling around the extensive Australian coastline, referenced to international locations, as part of a molecular systematic review using both nuclear and mitochondrial markers. The results indicate that while there are currently eight recognised species, the true diversity is close to double this, with a hotspot of endemism located in Australia. Complicated patterns were observed in southern Australia owing to two differing zones of introgression/admixture. Key drivers of diversity in the group appear to include plate tectonics, latitude, and historic barriers under glacial maxima, where an interplay between ready dispersal and habitat specialisation has led to regional panmixia but frequent geographic compartmentalisation within past and present landscapes. The findings have significant implications for biodiversity conservation, coastal and estuarine development, the basic foundations of field ecology, and for applied use such as in biomonitoring.

Genomic evidence supports the introgression between two sympatric stickleback species inhabiting the White Sea basin

Interspecies hybridization is driven by a complex interplay of factors where introgression plays an important role. In the present study, the transfer of genetic material, between two quite distant fish species from different genera, through spontaneous hybridization was documented with dedicated molecular and bioinformatics tools. We investigate the genomic landscape of putative stickleback-relative introgression by carefully analyzing the tractable transposable elements (TE) on the admixed genome of some individuals of two sympatric stickleback species inhabiting northwestern Russia, namely the three-spined (Gasterosteus aculeatus) and the nine-spined (Pungitius pungitius) sticklebacks. Our data revealed that unique TE amplification types exist, supporting our proposed hypothesis that infers on the interspecific introgression. By running a restriction site-associated DNA sequencing (RAD-Seq) with eight samples of G. aculeatus and P. pungitius and subjecting further the results to a contrasting analysis by variated bioinformatic tools, we identified the related introgression-linked markers. The admixture nature observed in a single sample of the nine-spined stickleback demonstrated the possible traces of remote introgression between these two species. Our work reveals the potential that introgression has on providing particular variants at a high-frequency speed while linking blocks of sequence with multiple functional mutations. However, even though our results are of significant interest, an increased number of samples displaying the introgression are required to further ascertain our conclusions.

Phylotranscriptomic evidence for pervasive ancient hybridization among Old World salamanders

Hybridization can leave genealogical signatures in an organism's genome, originating from the parental lineages and persisting over time. This potentially confounds phylogenetic inference methods that aim to represent evolution as a strictly bifurcating tree. We apply a phylotranscriptomic approach to study the evolutionary history of, and test for inter-lineage introgression in the Salamandridae, a Holarctic salamanders group of interest in studies of toxicity and aposematism, courtship behavior, and molecular evolution. Although the relationships between the 21 currently recognized salamandrid genera have been the subject of numerous molecular phylogenetic studies, some branches have remained controversial and sometimes affected by discordances between mitochondrial vs. nuclear trees. To resolve the phylogeny of this family, and understand the source of mito-nuclear discordance, we generated new transcriptomic (RNAseq) data for 20 salamandrids and used these along with published data, including 28 mitochondrial genomes, to obtain a comprehensive nuclear and mitochondrial perspective on salamandrid evolution. Our final phylotranscriptomic data set included 5455 gene alignments for 40 species representing 17 of the 21 salamandrid genera. Using concatenation and species-tree phylogenetic methods, we find (1) Salamandrina sister to the clade of the "True Salamanders" (consisting of Chioglossa, Mertensiella, Lyciasalamandra, and Salamandra), (2) Ichthyosaura sister to the Near Eastern genera Neurergus and Ommatotriton, (3) Triturus sister to Lissotriton, and (4) Cynops paraphyletic with respect to Paramesotriton and Pachytriton. Combining introgression tests and phylogenetic networks, we find evidence for introgression among taxa within the clades of "Modern Asian Newts" and "Modern European Newts". However, we could not unambiguously identify the number, position, and direction of introgressive events. Combining evidence from nuclear gene analysis with the observed mito-nuclear phylogenetic discordances, we hypothesize a scenario with hybridization and mitochondrial capture among ancestral lineages of (1) Lissotriton into Ichthyosaura and (2) Triturus into Calotriton, plus introgression of nuclear genes from Triturus into Lissotriton. Furthermore, both mitochondrial capture and nuclear introgression may have occurred among lineages assigned to Cynops. More comprehensive genomic data will, in the future, allow testing this against alternative scenarios involving hybridization with other, extinct lineages of newts.

Genetic Evidence for a Mixed Composition of the Genus Myoxocephalus (Cottoidei: Cottidae) Necessitates Generic Realignment

Sculpin fishes belonging to the family Cottidae represent a large and complex group, inhabiting a wide range of freshwater, brackish-water, and marine environments. Numerous studies based on analysis of their morphology and genetic makeup frequently provided controversial results. In the present work, we sequenced complete mitochondrial (mt) genomes and fragments of nuclear ribosomal DNA (rDNA) of the fourhorn sculpin Myoxocephalus quadricornis and some related cottids to increase the power of phylogenetic and taxonomic analyses of this complex fish group. A comparison of the My. quadricornis mt genomes obtained by us with other complete mt genomes available in GenBank has revealed a surprisingly low divergence (3.06 ± 0.12%) with Megalocottus platycephalus and, at the same time, a significantly higher divergence (7.89 ± 0.16%) with the species of the genus Myoxocephalus. Correspondingly, phylogenetic analyses have shown that My. quadricornis is clustered with Me. platycephalus but not with the Myoxocephalus species. Completely consistent patterns of divergence and tree topologies have been obtained based on nuclear rDNA. Thus, the multi-gene data in the present work indicates obvious contradictions in the relationships between the Myoxocephalus and Megalocottus species studied. An extensive phylogenetic analysis has provided evidence for a closer affinity of My. quadricornis with the species of the genus Megalocottus than with the species of the genus Myoxocephalus. A recombination analysis, along with the additional GenBank data, excludes introgression and/or incorrect taxonomic identification as the possible causative factors responsible for the observed closer affinity between the two species from different genera. The above facts necessitate realignment of the genera Myoxocephalus and Megalocottus. The genetic data supports the two recognized genera, Myoxocephalus and Megalocottus, but suggests changing their compositions through transferring My. quadricornis to the genus Megalocottus. The results of the present study resolve the relationships within a complex group of sculpin fishes and show a promising approach to phylogenetic systematics (as a key organizing principle in biodiversity research) for a better understanding of the taxonomy and evolution of fishes and for supplying relevant information to address various fish biodiversity conservation and management issues.

Population genetics of wild Macaca fascicularis with low-coverage shotgun sequencing of museum specimens

OBJECTIVES

Long-tailed macaques (Macaca fascicularis) are widely distributed throughout the mainland and islands of Southeast Asia, making them a useful model for understanding the complex biogeographical history resulting from drastic changes in sea levels throughout the Pleistocene. Past studies based on mitochondrial genomes (mitogenomes) of long-tailed macaque museum specimens have traced their colonization patterns throughout the archipelago, but mitogenomes trace only the maternal history. Here, our objectives were to trace phylogeographic patterns of long-tailed macaques using low-coverage nuclear DNA (nDNA) data from museum specimens.

METHODS

We performed population genetic analyses and phylogenetic reconstruction on nuclear single nucleotide polymorphisms (SNPs) from shotgun sequencing of 75 long-tailed macaque museum specimens from localities throughout Southeast Asia.

RESULTS

We show that shotgun sequencing of museum specimens yields sufficient genome coverage (average ~1.7%) for reconstructing population relationships using SNP data. Contrary to expectations of divergent results between nuclear and mitochondrial genomes for a female philopatric species, phylogeographical patterns based on nuclear SNPs proved to be closely similar to those found using mitogenomes. In particular, population genetic analyses and phylogenetic reconstruction from the nDNA identify two major clades within M. fascicularis: Clade A includes all individuals from the mainland along with individuals from northern Sumatra, while Clade B consists of the remaining island-living individuals, including those from southern Sumatra.

CONCLUSIONS

Overall, we demonstrate that low-coverage sequencing of nDNA from museum specimens provides enough data for examining broad phylogeographic patterns, although greater genome coverage and sequencing depth would be needed to distinguish between very closely related populations, such as those throughout the Philippines.

Gene flow and species delimitation in fishes of Western North America: Flannelmouth (Catostomus latipinnis) and Bluehead sucker (C. Pantosteus discobolus)

The delimitation of species boundaries, particularly those obscured by reticulation, is a critical step in contemporary biodiversity assessment. It is especially relevant for conservation and management of indigenous fishes in western North America, represented herein by two species with dissimilar life histories codistributed in the highly modified Colorado River (i.e., flannelmouth sucker, Catostomus latipinnis; bluehead sucker, C. (Pantosteus) discobolus). To quantify phylogenomic patterns and examine proposed taxonomic revisions, we first employed double-digest restriction site-associated DNA sequencing (ddRAD), yielding 39,755 unlinked SNPs across 139 samples. These were subsequently evaluated with multiple analytical approaches and by contrasting life history data. Three phylogenetic methods and a Bayesian assignment test highlighted similar phylogenomic patterns in each, but with considerable difference in presumed times of divergence. Three lineages were detected in bluehead sucker, supporting elevation of C. (P.) virescens to species status and recognizing C. (P.) discobolus yarrowi (Zuni bluehead sucker) as a discrete entity. Admixture in the latter necessitated a reevaluation of its contemporary and historic distributions, underscoring how biodiversity identification can be confounded by complex evolutionary histories. In addition, we defined three separate flannelmouth sucker lineages as ESUs (evolutionarily significant units), given limited phenotypic and genetic differentiation, contemporary isolation, and lack of concordance (per the genealogical concordance component of the phylogenetic species concept). Introgression was diagnosed in both species, with the Little Colorado and Virgin rivers in particular. Our diagnostic methods, and the agreement of our SNPs with previous morphological, enzymatic, and mitochondrial work, allowed us to partition complex evolutionary histories into requisite components, such as isolation versus secondary contact.

Patterns of mtDNA introgression suggest population replacement in Palaearctic whiskered bat species

Secondary contacts can play a major role in the evolutionary histories of species. Various taxa diverge in allopatry and later on come into secondary contact during range expansions. When they meet, their interactions and the extent of gene flow depend on the level of their ecological differentiation and the strength of their reproductive isolation. In this study, we present the multilocus phylogeography of two cryptic whiskered bat species, Myotis mystacinus and M. davidii, with a particular focus on their putative sympatric zone. Our findings suggest that M. mystacinus and M. davidii evolved in allopatry and came into secondary contact during range expansions. Individuals in the area of secondary contact, in Anatolia and the Balkans, have discordant population assignments based on the mitochondrial and the nuclear datasets. These observed patterns suggest that the local M. mystacinus populations hybridized with expanding M. davidii populations, which resulted in mitochondrial introgression from the former. In the introgression area, M. mystacinus individuals with concordant nuclear and mitochondrial genotypes were identified in relatively few locations, suggesting that the indigenous populations might have been largely replaced by invading M. davidii. Changing environmental conditions coupled with ecological competition is the likely reason for this replacement. Our study presents one possible example of a historical population replacement that was captured in phylogeographic patterns.

Rapid and Predictable Evolution of Admixed Populations Between Two Drosophila Species Pairs

The consequences of hybridization are varied, ranging from the origin of new lineages, introgression of some genes between species, to the extinction of one of the hybridizing species. We generated replicate admixed populations between two pairs of sister species of Drosophila: D. simulans and D. mauritiana; and D. yakuba and D. santomea Each pair consisted of a continental species and an island endemic. The admixed populations were maintained by random mating in discrete generations for over 20 generations. We assessed morphological, behavioral, and fitness-related traits from each replicate population periodically, and sequenced genomic DNA from the populations at generation 20. For both pairs of species, species-specific traits and their genomes regressed to those of the continental species. A few alleles from the island species persisted, but they tended to be proportionally rare among all sites in the genome and were rarely fixed within the populations. This paucity of alleles from the island species was particularly pronounced on the X-chromosome. These results indicate that nearly all foreign genes were quickly eliminated after hybridization and that selection against the minor species genome might be similar across experimental replicates.

Comparative phylogeography of trans-Andean freshwater fishes based on genome-wide nuclear and mitochondrial markers

The Neotropical region represents one of the greatest biodiversity hot spots on earth. Despite its unparalleled biodiversity, regional comparative phylogeographic studies are still scarce, with most focusing on model clades (e.g. birds) and typically examining a handful of loci. Here, we apply a genome-wide comparative phylogeographic approach to test hypotheses of codiversification of freshwater fishes in the trans-Andean region. Using target capture methods, we examined exon data for over 1,000 loci combined with complete mitochondrial genomes to study the phylogeographic history of five primary fish species (>150 individuals) collected from eight major river basins in Northwestern South America and Lower Central America. To assess their patterns of genetic structure, we inferred genealogical concordance taking into account all major aspects of phylogeography (within loci, across multiple genes, across species and among biogeographic provinces). Based on phylogeographic concordance factors, we tested four a priori biogeographic hypotheses, finding support for three of them and uncovering a novel, unexpected pattern of codiversification. The four emerging inter-riverine patterns are as follows: (a) Tuira + Atrato, (b) Ranchería + Catatumbo, (c) Magdalena system and (d) Sinú + Atrato. These patterns are interpreted as shared responses to the complex uplifting and orogenic processes that modified or sundered watersheds, allowing codiversification and speciation over geological time. We also find evidence of cryptic speciation in one of the species examined and instances of mitochondrial introgression in others. These results help further our knowledge of the historical geographic factors shaping the outstanding biodiversity of the Neotropics.

Deconstructing an infamous extinction crisis: Survival of Partula species on Moorea and Tahiti

Eleven of eighteen Society Island Partula species endemic to the Windward Island subgroup (Moorea and Tahiti) have been extirpated by an ill-advised biological control program. The conservation status of this critically endangered tree snail radiation is of considerable import, but is clouded by taxonomic uncertainty due to the extensive lack of congruence among species designations, diagnostic morphologies, and molecular markers. Using a combination of museum, captive, and remnant wild snails, we obtained the first high-resolution nuclear genomic perspective of the evolutionary relationships and survival of fourteen Windward Island Partula species, totaling 93 specimens. We analyzed ~1,607-28,194 nuclear genomic loci collected with the double digest restriction-site associated sequencing method. Results from phylogenomic trees, species estimation, and population assignment tests yielded monophyly of the Windward Island subgroup. Within this group, two well-supported clades encompassing five species complexes were recovered. Clade 1 was restricted to Tahiti and contained two species complexes: "P. affinis" (three species) and "P. otaheitana" (five species). Clade 2 occurred on Moorea and on Tahiti and consisted of three species complexes: one Tahitian, "P. clara/P. hyalina"; the other two, "P. taeniata" (three species) and "P. suturalis" (six species), Moorean. Our genomic results largely corroborated previous mitochondrial DNA survival estimates for Moorea and Tahiti, with all five species complexes having members surviving in captivity and/or as remnant wild populations, although the details vary in each case. Continued, proactive conservation and management may yet ensure a phylogenetically representative survival of the fabled Partula species of Moorea and Tahiti.

Integrative approach detects natural hybridization of sympatric lambaris species and emergence of infertile hybrids

Despite its relevance for ecology, evolution and conservation of species, natural hybridization and hybrids biology are still poorly studied in freshwater fish. Here, we tested the hypothesis that sympatric species Astyanax paranae and A. fasciatus are able to interbreed in the natural environment and presented evidence for the first record of hybridization between these species. We analyzed anatomical traits, gametogenesis, reproductive biology, and genetic variations of the COI and S7 genes of both species and putative hybrids. Intermediate morphometric and meristic features were observed in hybrids when compared to A. paranae and A. fasciatus. Overlap in reproductive season was showed for these species, with greater reproductive activity from August to January, but hybrids did not present any sign of gonadal maturation. Oogonia and perinucleolar follicles as well as spermatogonia and primary spermatocytes were found in hybrids, but previtellogenic and vitellogenic follicles, spermatids, and spermatozoa were absent. Moreover, several alterations in gametogenesis were detected, such as interrupted meiosis in both males and females, vacuolated and degenerated germ cells, increased interstitial tissue, and presence of immune cells. Molecular analyses supported the hypothesis of hybridization between A. paranae and A. fasciatus. Overall, our multidisciplinary approach also provides strong evidence that hybrids are infertile.

Bidirectional mitochondrial introgression between Korean cobitid fish mediated by hybridogenetic hybrids

Genomic introgression through interspecific hybridization has been observed in some species of the freshwater fish family Cobitidae. Within this family, a Cobitis hankugensis-Iksookimia longicorpa diploid-triploid hybrid species complex on the Korean peninsula is unique in displaying hybridogenesis, a unisexual reproduction mode that allows hybrids to mediate the transfer of mitochondrial DNA (but not nuclear DNA) between the two parent species. However, populations of the parental species in the wild have never been examined for the potential effect of introgression on their genomes. To address the genetic consequences of unisexual hybridization on the parental species, we examined genetic structure of the two parental species, C. hankugensis and I. longicorpa, in three independent natural habitats where they coexist with their hybrid complex using DNA sequence data of one mitochondrial gene and three nuclear genes. We found that mitochondrial introgression between the two species was extensive in all the examined localities, while there was no evidence of nuclear introgression across the species boundary. This result indicates that the hybridogenetic individuals mediate mitochondrial introgression from one species to the other, producing mito-nuclear mosaic genomes such as C. hankugensis nuclear genomes associated with I. longicorpa mitochondrial DNA and the reverse. The direction and degree of introgression varied among the three localities, but the underlying mechanisms for this observation proved elusive. Introgression might depend on which species serves as the predominant sperm or ovum donor or the environmental conditions of the localities. The present study suggests that introgressive hybridization between pure C. hankugensis and I. longicorpa species is highly likely where the two species co-occur with hybridogenetic individuals, but the consequence of introgression could be variable due to the history and environmental characteristics of particular populations across the parental species' ranges.

Convergent evolution misled taxonomy in schizothoracine fishes (Cypriniformes: Cyprinidae)

Highly specialized grade (HSG; genera Gymnocypris, Oxygymnocypris, Schizopygopsis, Platypharodon and Chuanchia) of the Schizothoracinae (Cypriniformes: Cyprinidae) are endemic to the Qinghai-Tibet Plateau (QTP). Previously, two distinct ecomorphs were recognized according to trophic traits. One was a limnetic omnivore with normal lower jaw morphology, terminal mouth, and moderate or dense gill rakers, mostly inhabiting in open water of lakes, including Gymnocypris and Oxygymnocypris. Another was a benthic feeder with inferior mouth, sparse gill rakers and sharp horny sheath on the lower jaw for scraping of attached prey off hard substrates, including Schizopygopsis, Platypharodon and Chuanchia. However, traditional taxonomy of HSG based on these trophic traits presented extensive conflicts with the molecular studies in recent years. The possible cause could be convergent evolution in morphology, retention of ancestral polymorphisms or mitochondrial introgression, but these hypotheses could not be assessed due to incomplete taxon sampling and only mitochondrial data employed in previous works. Here, we conducted the most comprehensive molecular analysis on HSG fishes to date, using four mitochondrial loci and 152,464 genome-wide SNPs, and including 21 of 24 putative species and one undescribed Schizopygopsis species. Both SNP and mtDNA trees confirmed extensive paraphyly of genera Gymnocypris and Schizopygopsis, where species often were clustered together by watershed instead of by genus. Basal split into the north clade B and the south clade C (ca. 3.03 Ma) approximately by the Tanggula-Tanitawen Mountains in SLAF tree coincided with a violent uplift of the QTP during the phase A of 'Qingzang movement' (ca. 3.6 Ma). Ancestral state reconstruction of the trophic ecomorph showed that the limnetic omnivore ecomorph had evolved repeatedly in clade B and C. Furthermore, we presented a striking case of convergent evolution between two 'subspecies' Gymnocypris chui chui and G. chui longimandibularis, which had diverged as early as two million years ago (ca. 2.42 Ma). Ecological analyses revealed that similar food utilization, particularly in zooplankton, was the main underlying driving force. This work showed an example of taxonomy with the most extensive errors at the genus/species levels due to convergent evolution and suggested that trophic traits could be misleading in fish taxonomy. Therefore, we propose a major generic revision for HSG species.

Genetic analysis reveals Finnish Formica fennica populations do not form a separate genetic entity from F. exsecta

Coptoformica Müller, 1923 is a subgenus of Formica Linnaeus, 1758 that consists of c. a dozen species of ants that typically inhabit open grassy habitats and build small nest mounds. The most recent addition to the group is Formica fennica Seifert, 2000. The description was based on morphological characters, but the species status has not been confirmed by molecular methods. In this study, we use thirteen DNA microsatellite markers and a partial mitochondrial COI gene sequence to assess the species status of F. fennica, by comparing the genetic variation among samples identified as F. fennica and six other boreal Formica (Coptoformica) species. Most of the species studied form separate, discontinuous clusters in phylogenetic and spatial analyses with only little intraspecific genetic variation. However, both nuclear and mitochondrial markers fail to separate the species pair F. exsecta Nylander, 1846 and F. fennica despite established morphological differences. The genetic variation within the F. exsecta/fennica group is extensive, but reflects spatial rather than morphological differences. Finnish F. fennica populations studied so far should not be considered a separate species, but merely a morph of F. exsecta.

Use of molecular data in zoo and aquarium collection management: Benefits, challenges, and best practices

The global zoo and aquarium community widely recognizes that its animal collections and cooperative breeding programs are facing a sustainability crisis. It has become commonly accepted that numerous priority species cannot be maintained unless new management strategies are adopted. While molecular data have the potential to greatly improve management across a range of scenarios, they have been generally underutilized by the zoo and aquarium community. This failure to effectively apply molecular data to collection management has been due, in part, to a paucity of resources within the community on which to base informed decisions about when the use of such data is appropriate and what steps are necessary to successfully integrate data into management. Here, we identify three broad areas of inquiry where molecular data can inform management: 1) taxonomic identification; 2) incomplete or unknown pedigrees; and 3) hereditary disease. Across these topics, we offer a discussion of the advantages, limitations, and considerations for applying molecular data to ex situ animal populations in a style accessible to zoo and aquarium professionals. Ultimately, we intend for this compiled information to serve as a resource for the community to help ensure that molecular projects directly and effectively benefit the long-term persistence of ex situ populations.

Systematics of Clupeiformes and testing for ecological limits on species richness in a trans-marine/freshwater clade

ABSTRACT Clupeiformes (herring, sardines, shad, anchovies and allies) are a globally distributed clade with nearly 400 marine, freshwater, and diadromous species. Although best known as filter feeding fishes that form large schools, this group occupies a diverse array of trophic guilds and habitats. Theory suggests that species richness in clades is modulated by ecological limits, which results in diversity-dependent clade growth, a pattern that most clades exhibit. As a trans-marine/freshwater clade that has undergone repeated transitions between marine and freshwaters, Clupeiformes are an excellent system for investigating the interplay between ecological diversity and macroevolutionary dynamics. In this study we review the systematics of Clupeiformes and explore discordance in phylogenetic relationships and divergence times between mitochondrial and nuclear loci. We then use comparative methods to test whether ecological limits regulate diversity in Clupeiformes. We find discordance in phylogenetic relationships at various taxonomic scales, but also considerable agreement between genomes. Our results suggest that trans-marine/freshwater clades are able to circumvent ecological limits on clade growth at regional, but not on local scales. Our study demonstrates that phylogenies are a critical link between ecology and macroevolutionary dynamics, and suggests habitat transitions can play a key role in shaping diversity patterns, particularly in the neotropics.

Complete paternally inherited mitogenomes of two freshwater mussels Unio pictorum and Sinanodonta woodiana (Bivalvia: Unionidae)

Freshwater bivalves from the family Unionidae usually have two very divergent mitogenomes, inherited according to the doubly uniparental model. The early divergence of these two mitogenomic lineages gives a unique opportunity to use two mitogenomic data sets in a single phylogenetic context. However, the number of complete sequences of the maternally inherited mitogenomes of these animals available in GenBank greatly exceeds that of the paternally inherited mitogenomes. This is a problem for phylogenetic reconstruction because it limits the use of both mitogenomic data sets. Moreover, since long branch attraction phenomenon can bias reconstructions if only a few but highly divergent taxa are considered, the shortage of the faster evolving paternally inherited mitogenome sequences is a real problem. Here we provide, for the first time, complete sequences of the M mitogenomes sampled from Polish populations of two species: native Unio pictorum and invasive Sinanodonta woodiana. It increases the available set of mitogenomic pairs to 18 species per family, and allows unambiguous reconstruction of phylogenetic relationships among them. The reconstructions based on M and F mitogenomes which were separated for many millions of years, and subject to differing evolutionary dynamics, are fully congruent.

Hidden diversity in the genus Tethya: comparing molecular and morphological techniques for species identification

Correctly determining species' identity is critical for estimating biodiversity and effectively managing marine populations, but is difficult for species that have few morphological traits or are highly plastic. Sponges are considered a taxonomically difficult group because they lack multiple consistent diagnostic features, which coupled with their common phenotypic plasticity, makes the presence of species complexes likely, but difficult to detect. Here, we investigated the evolutionary relationship of Tethya spp. in central New Zealand using both molecular and morphological techniques to highlight the potential for cryptic speciation in sponges. Phylogenetic reconstructions based on two mitochondrial markers (rnl, COI-ext) and one nuclear marker (18S) revealed three genetic clades, with one clade representing Tethya bergquistae and two clades belonging to what was a priori thought to be a single species, Tethya burtoni. Eleven microsatellite markers were also used to further resolve the T. burtoni group, revealing a division consistent with the 18S and rnl data. Morphological analysis based on spicule characteristics allowed T. bergquistae to be distinguished from T. burtoni, but revealed no apparent differences between the T. burtoni clades. Here, we highlight hidden genetic diversity within T. burtoni, likely representing a group consisting of incipient species that have undergone speciation but have yet to express clear morphological differences. Our study supports the notion that cryptic speciation in sponges may go undetected and diversity underestimated when using only morphology-based taxonomy, which has broad scale implications for conservation and management of marine systems.

The puzzling phylogeography of the haplochromine cichlid fish Astatotilapia burtoni

Astatotilapia burtoni is a member of the "modern haplochromines," the most species-rich lineage within the family of cichlid fishes. Although the species has been in use as research model in various fields of research since almost seven decades, including developmental biology, neurobiology, genetics and genomics, and behavioral biology, little is known about its spatial distribution and phylogeography. Here, we examine the population structure and phylogeographic history of A. burtoni throughout its entire distribution range in the Lake Tanganyika basin. In addition, we include several A. burtoni laboratory strains to trace back their origin from wild populations. To this end, we reconstruct phylogenetic relationships based on sequences of the mitochondrial DNA (mtDNA) control region (d-loop) as well as thousands of genomewide single nucleotide polymorphisms (SNPs) derived from restriction-associated DNA sequencing. Our analyses reveal high population structure and deep divergence among several lineages, however, with discordant nuclear and mtDNA phylogenetic inferences. Whereas the SNP-based phylogenetic hypothesis uncovers an unexpectedly deep split in A. burtoni, separating the populations in the southern part of the Lake Tanganyika basin from those in the northern part, analyses of the mtDNA control region suggest deep divergence between populations from the southwestern shoreline and populations from the northern and southeastern shorelines of Lake Tanganyika. This phylogeographic pattern and mitochondrial haplotype sharing between populations from the very North and the very South of Lake Tanganyika can only partly be explained by introgression linked to lake-level fluctuations leading to past contact zones between otherwise isolated populations and large-scale migration events.

Coupling molecular data and experimental crosses sheds light about species delineation: a case study with the genus Ciona

Molecular studies sometimes reveal evolutionary divergence within accepted species. Such findings can initiate taxonomic revision, as exemplified in the formerly recognized species Ciona intestinalis. While an increasing number of studies have examined the ecology, reproductive barriers and genetics of C. intestinalis and C. robusta, there are still much uncertainties regarding other species of this genus. Using experimental crosses and mitochondrial data, we investigated the evolutionary relationships among four native and introduced Ciona spp., found in sympatry in the Mediterranean Sea or English Channel. Outcome of 62 bi-parental reciprocal crosses between C. intestinalis, C. robusta, C. roulei and C. edwardsi showed that C. edwardsi is reproductively isolated from the other taxa, which is in agreement with its distinct location in the phylogenetic tree. Conversely, hybrids are easily obtained in both direction when crossing C. intestinalis and C. roulei, reinforcing the hypothesis of two genetically differentiated lineages but likely being from a same species. Altogether, this study sheds light on the evolutionary relationship in this complex genus. It also calls for further investigation notably based on genome-wide investigation to better describe the evolutionary history within the genus Ciona, a challenging task in a changing world where biological introductions are shuffling species distribution.

Divergent clades or cryptic species? Mito-nuclear discordance in a Daphnia species complex

BACKGROUND

Genetically divergent cryptic species are frequently detected by molecular methods. These discoveries are often a byproduct of molecular barcoding studies in which fragments of a selected marker are used for species identification. Highly divergent mitochondrial lineages and putative cryptic species are even detected in intensively studied animal taxa, such as the crustacean genus Daphnia. Recently, eleven such lineages, exhibiting genetic distances comparable to levels observed among well-defined species, were recorded in the D. longispina species complex, a group that contains several key taxa of freshwater ecosystems. We tested if three of those lineages represent indeed distinct species, by analyzing patterns of variation of ten nuclear microsatellite markers in six populations.

RESULTS

We observed a discordant pattern between mitochondrial and nuclear DNA, as all individuals carrying one of the divergent mitochondrial lineages grouped at the nuclear level with widespread, well-recognized species coexisting at the same localities (Daphnia galeata, D. longispina, and D. cucullata).

CONCLUSIONS

A likely explanation for this pattern is the introgression of the mitochondrial genome of undescribed taxa into the common species, either in the distant past or after long-distance dispersal. The occurrence of highly divergent but rare mtDNA lineages in the gene pool of widespread species would suggest that hybridization and introgression in the D. longispina species complex is frequent even across large phylogenetic distances, and that discoveries of such distinct clades must be interpreted with caution. However, maintenance of ancient polymorphisms through selection is another plausible alternative that may cause the observed discordance and cannot be entirely excluded.