Michael Matschiner

TANDEM: integrating automated allele binning into genetics and genomics workflows

SUMMARY Computer programs for the statistical analysis of microsatellite data use allele length variation to infer, e.g. population genetic parameters, to detect quantitative trait loci or selective sweeps. However, observed allele lengths are usually inaccurate and may deviate from the expected periodicity of repeats. The common practice of rounding to the nearest whole number frequently results in miscalls and underestimations of allelic richness. Manual sorting of allele lengths into discrete classes, a process called binning, is tedious and error-prone. Here, we present a new program for the automated binning of microsatellite allele lengths to overcome these problems and to facilitate high-throughput allele binning. AVAILABILITY www.evolution.unibas.ch/salzburger/software.htm. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.

On the origin and trigger of the notothenioid adaptive radiation

Adaptive radiation is usually triggered by ecological opportunity, arising through (i) the colonization of a new habitat by its progenitor; (ii) the extinction of competitors; or (iii) the emergence of an evolutionary key innovation in the ancestral lineage. Support for the key innovation hypothesis is scarce, however, even in textbook examples of adaptive radiation. Antifreeze glycoproteins (AFGPs) have been proposed as putative key innovation for the adaptive radiation of notothenioid fishes in the ice-cold waters of Antarctica. A crucial prerequisite for this assumption is the concurrence of the notothenioid radiation with the onset of Antarctic sea ice conditions. Here, we use a fossil-calibrated multi-marker phylogeny of nothothenioid and related acanthomorph fishes to date AFGP emergence and the notothenioid radiation. All time-constraints are cross-validated to assess their reliability resulting in six powerful calibration points. We find that the notothenioid radiation began near the Oligocene-Miocene transition, which coincides with the increasing presence of Antarctic sea ice. Divergence dates of notothenioids are thus consistent with the key innovation hypothesis of AFGP. Early notothenioid divergences are furthermore congruent with vicariant speciation and the breakup of Gondwana.

Gene flow by larval dispersal in the Antarctic notothenioid fish Gobionotothen gibberifrons

The diversification of the teleost suborder Notothenioidei (Perciformes) in Antarctic waters provides one of the most striking examples of a marine adaptive radiation. Along with a number of adaptations to the cold environment, such as the evolution of antifreeze glycoproteins, notothenioids diversified into eight families and at least 130 species. Here, we investigate the genetic population structure of the humped rockcod (Gobionotothen gibberifrons), a benthic notothenioid fish. Six populations were sampled at different locations around the Scotia Sea, comprising a large part of the species’ distribution range (N = 165). Our analyses based on mitochondrial DNA sequence data (352 bp) and eight microsatellite markers reveal a lack of genetic structuring over large geographic distances (ΦST ≤ 0.058, FST ≤ 0.005, P values nonsignificant). In order to test whether this was due to passive larval dispersal, we used GPS‐tracked drifter trajectories, which approximate movement of passive surface particles with ocean currents. The drifter data indicate that the Antarctic Circumpolar Current (ACC) connects the sampling locations in one direction only (west–east), and that passive transport is possible within the 4‐month larval period of G. gibberifrons. Indeed, when applying the isolation‐with‐migration model in IMA, strong unidirectional west‐east migration rates are detected in the humped rockcod. This leads us to conclude that, in G. gibberifrons, genetic differentiation is prevented by gene flow via larval dispersal with the ACC.

Parallel ecological diversification in Antarctic notothenioid fishes as evidence for adaptive radiation

Antarctic notothenioid fishes represent a rare example of a marine species flock. They evolved special adaptations to the extreme environment of the Southern Ocean including antifreeze glycoproteins. Although lacking a swim bladder, notothenioids have diversified from their benthic ancestor into a wide array of water column niches, such as epibenthic, semipelagic, cryopelagic and pelagic habitats. Applying stable carbon (C) and nitrogen (N) isotope analyses to gain information on feeding ecology and foraging habitats, we tested whether ecological diversification along the benthic–pelagic axis followed a single directional trend in notothenioids, or whether it evolved independently in several lineages. Population samples of 25 different notothenioid species were collected around the Antarctic Peninsula, the South Orkneys and the South Sandwich Islands. The C and N stable isotope signatures span a broad range (mean δ13C and δ15N values between −25.4‰ and −21.9‰ and between 8.5‰ and 13.8‰, respectively), and pairwise niche overlap between four notothenioid families was highly significant. Analysis of isotopic disparity‐through‐time on the basis of Bayesian inference and maximum‐likelihood phylogenies, performed on a concatenated mitochondrial (cyt b) and nuclear gene (myh6, Ptr and tbr1) data set (3148 bp), showed that ecological diversification into overlapping feeding niches has occurred multiple times in parallel in different notothenioid families. This convergent diversification in habitat and trophic ecology is a sign of interspecific competition and characteristic for adaptive radiations.

A tribal level phylogeny of Lake Tanganyika cichlid fishes based on a genomic multi-marker approach

Graphical abstract

Evolution of the immune system influences speciation rates in teleost fishes

Teleost fishes constitute the most species-rich vertebrate clade and exhibit extensive genetic and phenotypic variation, including diverse immune defense strategies. The genomic basis of a particularly aberrant strategy is exemplified by Atlantic cod, in which a loss of major histocompatibility complex (MHC) II functionality coincides with a marked expansion of MHC I genes. Through low-coverage genome sequencing (9–39×), assembly and comparative analyses for 66 teleost species, we show here that MHC II is missing in the entire Gadiformes lineage and thus was lost once in their common ancestor. In contrast, we find that MHC I gene expansions have occurred multiple times, both inside and outside this clade. Moreover, we identify an association between high MHC I copy number and elevated speciation rates using trait-dependent diversification models. Our results extend current understanding of the plasticity of the adaptive immune system and suggest an important role for immune-related genes in animal diversification.

Comparative population genetics of seven notothenioid fish species reveals high levels of gene flow along ocean currents in the southern Scotia Arc, Antarctica

The Antarctic fish fauna is characterized by high endemism and low species diversity with one perciform suborder, the Notothenioidei, dominating the whole species assemblage on the shelves and slopes. Notothenioids diversified in situ through adaptive radiation and show a variety of life history strategies as adults ranging from benthic to pelagic modes. Their larval development is unusually long, lasting from a few months to more than a year, and generally includes a pelagic larval stage. Therefore, the advection of eggs and larvae with ocean currents is a key factor modulating population connectivity. Here, we compare the genetic population structures and gene flow of seven ecologically distinct notothenioid species of the southern Scotia Arc based on nuclear microsatellites and mitochondrial DNA sequences (D-loop/cytochrome b). The seven species belong to the families Nototheniidae (Gobionotothen gibberifrons, Lepidonotothen squamifrons, Trematomus eulepidotus, T. newnesi) and Channichthyidae (Chaenocephalus aceratus, Champsocephalus gunnari, Chionodraco rastrospinosus). Our results show low-population differentiation and high gene flow for all investigated species independent of their adult life history strategies. In addition, gene flow is primarily in congruence with the prevailing ocean current system, highlighting the role of larval dispersal in population structuring of notothenioids.

Disentangling incomplete lineage sorting and introgression to refine species-tree estimates for Lake Tanganyika cichlid fishes

Abstract Adaptive radiation is thought to be responsible for the evolution of a great portion of the past and present diversity of life. Instances of adaptive radiation, characterized by the rapid emergence of an array of species as a consequence to their adaptation to distinct ecological niches, are important study systems in evolutionary biology. However, because of the rapid lineage formation in these groups, and occasional gene flow between the participating species, it is often difficult to reconstruct the phylogenetic history of species that underwent an adaptive radiation. In this study, we present a novel approach for species‐tree estimation in rapidly diversifying lineages, where introgression is known to occur, and apply it to a multimarker data set containing up to 16 specimens per species for a set of 45 species of East African cichlid fishes (522 individuals in total), with a main focus on the cichlid species flock of Lake Tanganyika. We first identified, using age distributions of most recent common ancestors in individual gene trees, those lineages in our data set that show strong signatures of past introgression. This led us to formulate three hypotheses of introgression between different lineages of Tanganyika cichlids: the ancestor of Boulengerochromini (or of Boulengerochromini and Bathybatini) received genomic material from the derived H‐lineage; the common ancestor of Cyprichromini and Perissodini experienced, in turn, introgression from Boulengerochromini and/or Bathybatini; and the Lake Tanganyika Haplochromini and closely related riverine lineages received genetic material from Cyphotilapiini. We then applied the multispecies coalescent model to estimate the species tree of Lake Tanganyika cichlids, but excluded the lineages involved in these introgression events, as the multispecies coalescent model does not incorporate introgression. This resulted in a robust species tree, in which the Lamprologini were placed as sister lineage to the H‐lineage (including the Eretmodini), and we identify a series of rapid splitting events at the base of the H‐lineage. Divergence ages estimated with the multispecies coalescent model were substantially younger than age estimates based on concatenation, and agree with the geological history of the Great Lakes of East Africa. Finally, we formally tested the three hypotheses of introgression using a likelihood framework, and find strong support for introgression between some of the cichlid tribes of Lake Tanganyika.

Diversity and disparity through time in the adaptive radiation of Antarctic notothenioid fishes.

According to theory, adaptive radiation is triggered by ecological opportunity that can arise through the colonization of new habitats, the extinction of antagonists or the origin of key innovations. In the course of an adaptive radiation, diversification and morphological evolution are expected to slow down after an initial phase of rapid adaptation to vacant ecological niches, followed by speciation. Such ‘early bursts’ of diversification are thought to occur because niche space becomes increasingly filled over time. The diversification of Antarctic notothenioid fishes into over 120 species has become one of the prime examples of adaptive radiation in the marine realm and has likely been triggered by an evolutionary key innovation in the form of the emergence of antifreeze glycoproteins. Here, we test, using a novel time‐calibrated phylogeny of 49 species and five traits that characterize notothenioid body size and shape as well as buoyancy adaptations and habitat preferences, whether the notothenioid adaptive radiation is compatible with an early burst scenario. Extensive Bayesian model comparison shows that phylogenetic age estimates are highly dependent on model choice and that models with unlinked gene trees are generally better supported and result in younger age estimates. We find strong evidence for elevated diversification rates in Antarctic notothenioids compared to outgroups, yet no sign of rate heterogeneity in the course of the radiation, except that the notothenioid family Artedidraconidae appears to show secondarily elevated diversification rates. We further observe an early burst in trophic morphology, suggesting that the notothenioid radiation proceeds in stages similar to other prominent examples of adaptive radiation.

Bayesian Phylogenetic Estimation of Clade Ages Supports Trans-Atlantic Dispersal of Cichlid Fishes.

Abstract Divergence‐time estimation based on molecular phylogenies and the fossil record has provided insights into fundamental questions of evolutionary biology. In Bayesian node dating, phylogenies are commonly time calibrated through the specification of calibration densities on nodes representing clades with known fossil occurrences. Unfortunately, the optimal shape of these calibration densities is usually unknown and they are therefore often chosen arbitrarily, which directly impacts the reliability of the resulting age estimates. As possible solutions to this problem, two nonexclusive alternative approaches have recently been developed, the “fossilized birth‐death” (FBD) model and “total‐evidence dating.” While these approaches have been shown to perform well under certain conditions, they require including all (or a random subset) of the fossils of each clade in the analysis, rather than just relying on the oldest fossils of clades. In addition, both approaches assume that fossil records of different clades in the phylogeny are all the product of the same underlying fossil sampling rate, even though this rate has been shown to differ strongly between higher level taxa. We here develop a flexible new approach to Bayesian age estimation that combines advantages of node dating and the FBD model. In our new approach, calibration densities are defined on the basis of first fossil occurrences and sampling rate estimates that can be specified separately for all clades. We verify our approach with a large number of simulated data sets, and compare its performance to that of the FBD model. We find that our approach produces reliable age estimates that are robust to model violation, on par with the FBD model. By applying our approach to a large data set including sequence data from over 1000 species of teleost fishes as well as 147 carefully selected fossil constraints, we recover a timeline of teleost diversification that is incompatible with previously assumed vicariant divergences of freshwater fishes. Our results instead provide strong evidence for transoceanic dispersal of cichlids and other groups of teleost fishes.