Search

The availability of public genomic resources can greatly assist biodiversity assessment, conservation, and restoration efforts by providing evidence for scientifically informed management decisions. Here we survey the main approaches and applications in biodiversity and conservation genomics, considering practical factors, such as cost, time, prerequisite skills, and current shortcomings of applications. Most approaches perform best in combination with reference genomes from the target species or closely related species. We review case studies to illustrate how reference genomes can facilitate biodiversity research and conservation across the tree of life. We conclude that the time is ripe to view reference genomes as fundamental resources and to integrate their use as a best practice in conservation genomics., Trends in Genetics, 39 (7), ISSN:0168-9525

The role of interspecific hybridization in the context of diversification dynamics has recently seen increasing attention. Genomic research has now made it abundantly clear that both hybridization and introgression – the exchange of genetic material through hybridization – are far more common than previously thought. Moreover, even highly divergent species were found to hybridize and backcross. These findings raise the question whether commonly used methods for the detection of introgression are applicable to such divergent hybridizing species, given that most of these methods were originally developed for analyses at the level of populations and recently diverged species. In particular, the assumption of constant evolutionary rates, which is implicit in many commonly used approaches, is more likely to be violated as evolutionary divergence increases. To test the limitations of introgression detection methods when being applied to divergent species, we simulated thousands of genomic datasets under a wide range of settings, with varying degrees of among-species rate variation and introgression. Using these simulated datasets, we were able to show that commonly applied statistical methods, including the D-statistic and tests based on sets of phylogenetic trees, produce false-positive signals of introgression between highly divergent taxa when these have different rates of evolution. These misleading signals are caused by the presence of homoplasies that occur at different rates when rate variation is present. To distinguish between the patterns caused by rate variation and genuine introgression, we developed a new test that is based on the expected clustering of introgressed sites and implemented this test in the program Dsuite.

Repetitive DNA make up a considerable fraction of most eukaryotic genomes. In fish, transposable element (TE) activity have coincided with rapid species diversification. Here, we annotated the repetitive content in 100 genome assemblies, covering the major branches of the diverse lineage of teleost fish. We investigated if TE content correlates with family level net diversification rates and found support for a weak negative correlation. Further, we found that TE content, the degree of parental care and short tandem repeat (STR) content contributed to genome size variability. In contrast to TEs, STR content showed a negative relationship with genome size. STR content did not correlate with TE content, which implies independent evolutionary paths. Last, marine and freshwater fish have large differences in STR content. The most extreme propagation was found in the genomes of codfish species and Atlantic herring. Such a high density of STRs is likely to increase the mutational load, which we propose could be counterbalanced by high fecundity as seen in codfishes and herring.

Currents are unique drivers of oceanic phylogeography and so determine the distribution of marine coastal species, along with past glaciations and sea level changes. Here, we reconstruct the worldwide colonization history of eelgrass (Zostera marinaL.), the most widely distributed marine flowering plant or seagrass from its origin in the Northwest Pacific, based on nuclear and chloroplast genomes. We identified two divergent Pacific clades with evidence for admixture along the East Pacific coast. Multiple west to east (trans-Pacific) colonization events support the key role of the North Pacific Current. Time-calibrated nuclear and chloroplast phylogenies yielded concordant estimates of the arrival ofZ. marinain the Atlantic through the Canadian Arctic, suggesting that eelgrass-based ecosystems, hotspots of biodiversity and carbon sequestration, have only been present since ∼208 Kya (thousand years ago). Mediterranean populations were founded ∼53 Kya while extant distributions along western and eastern Atlantic shores coincide with the end of the Last Glacial Maximum (∼20 Kya). The recent colonization and 5-to 7-fold lower genomic diversity of Atlantic compared to the Pacific populations raises concern and opportunity about how Atlantic eelgrass might respond to rapidly warming coastal oceans.

While the inference of species trees from molecular sequences has become a common type of analysis in studies of species diversification, few programs so far allow for the use of single-nucleotide polymorphisms (SNPs) for the same purpose. In this book chapter, I discuss the use of the Bayesian program SNAPP, which infers the species tree by mathematically integrating over all possible genealogies at each SNP. In particular, I focus on a molecular clock model developed for SNAPP, allowing the inference of divergence times together with the species tree topology and the population size, directly from SNP datasets in variant call format. With the growing availability of SNP datasets for multiple closely related species, this approach is becoming increasingly relevant for the reconstruction of the temporal framework of recent species diversification.

In the NHMO zoological collections, specimens from the Arctic include about 9,000 mammals and 7,100 birds, whereas the Insect Collection holds about 105,000 specimens plus more than hundred jars with unsorted material. The Fish Collection contains approximately 1,400 specimens, while the Herptile Collection (amphibians & reptiles) holds only thirty-one specimens of three taxa. Many of these specimens originate from expeditions to E Greenland, N Canada, Svalbard, Novaya Zemlya, Finnmark, and NE Siberia in the period 1898 to 1966. Furthermore, the DNA Bank has about 5,600 tissue and extracted DNA samples, mostly sampled from wild animals during the last decades but also from specimens in the voucher collections. Most of the Arctic specimens have been digitized and are available in online data portals like GBIF, except for the Insect Collection, where only the type material and about 30 percent of the total specimens are digitized.

Antarctic notothenioid fishes are recognised as one of the rare examples of adaptive radiation in the marine system. Withstanding the freezing temperatures of Antarctic waters, these fishes have diversified into over 100 species within no more than 10-20 million years. However, the exact species richness of the radiation remains contested. In the genus

Adaptive radiation is the likely source of much of the ecological and morphological diversity of life1–4. How adaptive radiations proceed and what determines their extent remains unclear in most cases1,4. Here we report the in-depth examination of the spectacular adaptive radiation of cichlid fishes in Lake Tanganyika. On the basis of whole-genome phylogenetic analyses, multivariate morphological measurements of three ecologically relevant trait complexes (body shape, upper oral jaw morphology and lower pharyngeal jaw shape), scoring of pigmentation patterns and approximations of the ecology of nearly all of the approximately 240 cichlid species endemic to Lake Tanganyika, we show that the radiation occurred within the confines of the lake and that morphological diversification proceeded in consecutive trait-specific pulses of rapid morphospace expansion. We provide empirical support for two theoretical predictions of how adaptive radiations proceed, the ‘early-burst’ scenario1,5 (for body shape) and the stages model1,6,7 (for all traits investigated). Through the analysis of two genomes per species and by taking advantage of the uneven distribution of species in subclades of the radiation, we further show that species richness scales positively with per-individual heterozygosity, but is not correlated with transposable element content, number of gene duplications or genome-wide levels of selection in coding sequences. Analyses of molecular, anatomical, pigmentation and ecological characteristics of nearly all of the approximately 240 species of cichlid fishes in Lake Tanganyika show that the massive adaptive radiation occurred within the confines of the lake through trait-specific pulses of accelerated evolution.

Stochastic simulation is a key tool in population genetics, since the models involved are often analytically intractable and simulation is usually the only way of obtaining ground-truth data to evaluate inferences. Because of this, a large number of specialized simulation programs have been developed, each filling a particular niche, but with largely overlapping functionality and a substantial duplication of effort. Here, we introduce msprime version 1.0, which efficiently implements ancestry and mutation simulations based on the succinct tree sequence data structure and the tskit library. We summarize msprime’s many features, and show that its performance is excellent, often many times faster and more memory efficient than specialized alternatives. These high-performance features have been thoroughly tested and validated, and built using a collaborative, open source development model, which reduces duplication of effort and promotes software quality via community engagement.

Progress in genome sequencing now enables the large-scale generation of reference genomes. Various international initiatives aim to generate reference genomes representing global biodiversity. These genomes provide unique insights into genomic diversity and architecture, thereby enabling comprehensive analyses of population and functional genomics, and are expected to revolutionize conservation genomics., Trends in Ecology & Evolution, 37 (3), ISSN:0169-5347, ISSN:1872-8383

Supergenes are sets of genes that are inherited as a single marker and encode complex phenotypes through their joint action. They are identified in an increasing number of organisms, yet their origins and evolution remain enigmatic. In Atlantic cod, four megabase-scale supergenes have been identified and linked to migratory lifestyle and environmental adaptations. Here we investigate the origin and maintenance of these four supergenes through analysis of whole-genome-sequencing data, including a new long-read-based genome assembly for a non-migratory Atlantic cod individual. We corroborate the finding that chromosomal inversions underlie all four supergenes, and we show that they originated at different times between 0.40 and 1.66 million years ago. We reveal gene flux between supergene haplotypes where migratory and stationary Atlantic cod co-occur and conclude that this gene flux is driven by gene conversion, on the basis of an increase in GC content in exchanged sites. Additionally, we find evidence for double crossover between supergene haplotypes, leading to the exchange of an ~275 kilobase fragment with genes potentially involved in adaptation to low salinity in the Baltic Sea. Our results suggest that supergenes can be maintained over long timescales in the same way as hybridizing species, through the selective purging of introduced genetic variation.

Stochastic simulation is a key tool in population genetics, since the models involved are often analytically intractable and simulation is usually the only way of obtaining ground-truth data to evaluate inferences. Because of this necessity, a large number of specialised simulation programs have been developed, each filling a particular niche, but with largely overlapping functionality and a substantial duplication of effort. Here, we introducemsprimeversion 1.0, which efficiently implements ancestry and mutation simulations based on the succinct tree sequence data structure andtskitlibrary. We summarisemsprime’s many features, and show that its performance is excellent, often many times faster and more memory efficient than specialised alternatives. These high-performance features have been thoroughly tested and validated, and built using a collaborative, open source development model, which reduces duplication of effort and promotes software quality via community engagement.

Supergenes are sets of genes that are inherited as a single marker and encode complex phenotypes through their joint action. They are identified in an increasing number of organisms, yet their origins and evolution remain enigmatic. In Atlantic cod, four large supergenes have been identified and linked to migratory lifestyle and environmental adaptations. Here, we investigate the origin and fate of these four supergenes through analysis of whole-genome-sequencing data, including a new long-read-based genome assembly for a non-migratory Atlantic cod individual. We corroborate that chromosomal inversions underlie all four supergenes, and show that they originated separately, between 0.40 and 1.66 million years ago. While introgression was not involved in the origin of the four supergenes, we reveal gene flow between inverted and noninverted supergene haplotypes, occurring both through gene conversion and double crossover. Moreover, the presence of genes linked to salinity adaptations in a sequence transferred through double crossover indicates that these sequences exchanged between the haplotypes are subject to selection. Our results suggest that the fate of supergenes is comparable to that of hybridizing species, by depending on the degree to which separation is maintained through purging of introduced genetic variation.

Significance Both theory and experiments suggest that fishing can drive the evolution of an earlier maturation age. However, determining whether changes in the wild are the result of fisheries-induced evolution has been difficult. Temporal, genome-wide datasets can directly reveal responses to selection. Here, we investigate the genomes of two wild Atlantic cod populations from samples that pre- and postdate periods of intensive fishing. Although phenotypic changes suggest fisheries-induced evolution, we do not find evidence for any strong genomic change or loss of genetic diversity. While evolution could have occurred through undetectable frequency changes at many loci, the irreversible loss of late-maturing genotypes appears unlikely. Instead, we suggest that recovery of former phenotypes is possible with reduced fishing pressure., The mode and extent of rapid evolution and genomic change in response to human harvesting are key conservation issues. Although experiments and models have shown a high potential for both genetic and phenotypic change in response to fishing, empirical examples of genetic responses in wild populations are rare. Here, we compare whole-genome sequence data of Atlantic cod (Gadus morhua) that were collected before (early 20th century) and after (early 21st century) periods of intensive exploitation and rapid decline in the age of maturation from two geographically distinct populations in Newfoundland, Canada, and the northeast Arctic, Norway. Our temporal, genome-wide analyses of 346,290 loci show no substantial loss of genetic diversity and high effective population sizes. Moreover, we do not find distinct signals of strong selective sweeps anywhere in the genome, although we cannot rule out the possibility of highly polygenic evolution. Our observations suggest that phenotypic change in these populations is not constrained by irreversible loss of genomic variation and thus imply that former traits could be reestablished with demographic recovery.

Sex is a fundamental trait that is determined, depending on the species, by different environmental and/or genetic factors, including various types of sex chromosomes. While the functioning and emergence of sex chromosomes have been explored in species scattered across the eukaryotic tree of life, little is known about tempo and mode of sex chromosome evolution in closely related species. Here, we examine the dynamics of sex chromosome evolution in an archetypical example of adaptive radiation, the cichlid fishes of African Lake Tanganyika. Through inspection of male and female genomes from 244 cichlid taxa and the analysis of transcriptomes from 66 taxa, we identify signatures of sex chromosomes in 79 taxa, involving 12 different linkage groups. We estimate that Tanganyikan cichlids have the highest rates of sex chromosome turnover and heterogamety transitions known to date. We further show that the recruitment of chromosomes as sex chromosomes is not at random and that some chromosomes have convergently emerged as sex chromosomes in cichlids, which provides empirical support to the “limited options” hypothesis of sex chromosome evolution.

Rhinoceroses are facing an uncertain future due to human-induced population-size reductions. In this issue of Cell, Liu et al. analyze rhinoceros genomes to study their past and find support for an early divergence between African and Asian species, signatures of past hybridization, and evidence of loss of genetic diversity.

Adaptive radiation is the likely source of much of the ecological and morphological diversity of life

Cichlid fishes are celebrated for their vast taxonomic, phenotypic, and ecological diversity; however, a central aspect of their evolution — the timeline of their diversification — remains contentious. Here, we generate draft genome assemblies of 14 species representing the global cichlid diversity and integrate these into a new phylogenomic hypothesis of cichlid and teleost evolution that we time-calibrate with 58 re-evaluated fossil constraints and a new Bayesian model accounting for fossil-assignment uncertainty. Our results support cichlid diversification long after the breakup of the supercontinent Gondwana and lay the foundation for precise temporal reconstructions of the exceptional continental cichlid adaptive radiations., Cichlids are a model for adaptive radiation, but the timing of their diversification is debated. Here the authors assemble 14 cichlid genomes, introduce a Bayesian approach to account for fossil-assignment uncertainty, and present a dated phylogenomic hypothesis of cichlid and teleost evolution.

The closure of the Isthmus of Panama has long been considered to be one of the best defined biogeographic calibration points for molecular divergence-time estimation. However, geological and biological evidence has recently cast doubt on the presumed timing of the initial isthmus closure around 3 Ma but has instead suggested the existence of temporary land bridges as early as the Middle or Late Miocene. The biological evidence supporting these earlier land bridges was based either on only few molecular markers or on concatenation of genome-wide sequence data, an approach that is known to result in potentially misleading branch lengths and divergence times, which could compromise the reliability of this evidence. To allow divergence-time estimation with genomic data using the more appropriate multispecies coalescent (MSC) model, we here develop a new method combining the single-nucleotide polymorphism-based Bayesian species-tree inference of the software SNAPP with a molecular clock model that can be calibrated with fossil or biogeographic constraints. We validate our approach with simulations and use our method to reanalyze genomic data of Neotropical army ants (Dorylinae) that previously supported divergence times of Central and South American populations before the isthmus closure around 3 Ma. Our reanalysis with the MSC model shifts all of these divergence times to ages younger than 3 Ma, suggesting that the older estimates supporting the earlier existence of temporary land bridges were artifacts resulting at least partially from the use of concatenation. We then apply our method to a new restriction-site associated DNA-sequencing data set of Neotropical sea catfishes (Ariidae) and calibrate their species tree with extensive information from the fossil record. We identify a series of divergences between groups of Caribbean and Pacific sea catfishes around 10 Ma, indicating that processes related to the emergence of the isthmus led to vicariant speciation already in the Late Miocene, millions of years before the final isthmus closure.

How new species form in the ocean, and thus what determines the diversity of fish in the sea, is not well understood. A study in Caribbean coral-reef fishes sheds light on the genomic underpinnings of diversification in the marine realm.

The Gasterosteidae fish family hosts several species that are important models for eco-evolutionary, genetic and genomic research. In particular, a wealth of genetic and genomic data has been generated for the three-spined stickleback (Gasterosteus aculeatus), the ‘ecology’s supermodel’, while the genomic resources for the nine-spined stickleback (Pungitius pungitius) have remained relatively scarce. Here, we report a high-quality chromosome-level genome assembly of P. pungitius consisting of 5,303 contigs (N50 = 1.2 Mbp) with a total size of 521 Mbp. These contigs were mapped to 21 linkage groups using a high-density linkage map, yielding a final assembly with 98.5% BUSCO completeness. A total of 25,062 protein-coding genes were annotated, and ca. 23% of the assembly was found to consist of repetitive elements. A comprehensive analysis of repetitive elements uncovered centromeric-specific tandem repeats and provided insights into the evolution of retrotransposons. A multigene phylogenetic analysis inferred a divergence time of about 26 million years (MYA) between nine- and three-spined sticklebacks, which is far older than the commonly assumed estimate of 13 MYA. Compared to the three-spined stickleback, we identified an additional duplication of several genes in the hemoglobin cluster. Sequencing data from populations adapted to different environments indicated potential copy number variations in hemoglobin genes. Furthermore, genome-wide synteny comparisons between three- and nine-spined sticklebacks identified chromosomal rearrangements underlying the karyotypic differences between the two species. The high-quality chromosome-scale assembly of the nine-spined stickleback genome obtained with long-read sequencing technology provides a crucial resource for comparative and population genomic investigations of stickleback fishes and teleosts.

Cichlid fishes are one of the most important model systems for evolutionary biology. Unfortunately, however, the timeline of cichlid diversification is still insufficiently known and limits our understanding of the mechanisms that generated their spectacular diversity. The uncertainty regarding this timeline stems from a decades-old controversy surrounding the phylogeographic history of cichlid fishes. Did cichlid subfamilies diversify as the result of Gondwanan vicariance, as supported by their distribution on former Gondwanan landmasses? Or did they diverge much more recently through oceanic dispersal, as suggested by the fossil record? While a large number of studies have already addressed this question with molecular-clock analyses, no single conclusion has emerged from these investigations. Here, I review the molecular evidence for Gondwanan vicariance or trans-Atlantic dispersal resulting from these studies. I discuss the weaknesses and strengths of each study, aiming to promote the formation of consensus on the matter and to prevent the repetition of previously made mistakes. I find that after accounting for inappropriate calibration strategies and saturation in mitochondrial datasets, the molecular evidence points to trans-Atlantic dispersal long after continental separation, probably around 75–60 Ma.

Patterson’s D, also known as the ABBA-BABA statistic, and related statistics such as the f4-ratio, are commonly used to assess evidence of gene flow between populations or closely related species. Currently available implementations require custom file formats and are impractical to evaluate all gene flow hypotheses across datasets with many populations or species.Dsuite is a fast C++ implementation, allowing genome scale calculations of the D and f4-ratio statistics across all combinations of tens or hundreds of populations or species directly from a variant call format (VCF) file. Furthermore, the program can provide evidence of whether introgression is confined to specific loci and aid in interpretation of a system of f4-ratio results by implementing the ‘f-branch’ method.Dsuite is available at https://github.com/millanek/Dsuite, is straightforward to use, substantially more computationally efficient than other comparable programs, and presents a novel suite of tools and statistics, including some not previously available in any software package.Thus, Dsuite facilitates assessment of evidence for gene flow, especially across large genomic datasets.

Fish catch color with rods Vertebrates are typically thought to have a consistent system for processing light, in which multiple cone opsins permit color vision during the day, but a single rod opsin provides only monochrome vision in the dark. Musilova et al. analyzed more than 100 deep-sea fish genomes and found a previously unknown proliferation of rod opsin genes that generate rod opsin photopigments that are tuned to different wavelengths of light. These receptors may allow the fish to perceive bioluminescent signals that pervade their deep-sea environment. These results reveal a previously undescribed visual system that allows for color vision in the dark. Science , this issue p. 588

Incomplete lineage sorting (ILS) can lead to biased divergence time estimates. To explore if and how ILS has influenced the results of a recent study of worldwide phylogeny of three-spined sticklebacks (Gasterosteus aculeatus), we estimated divergence times among major clades by applying both a concatenation approach and the multispecies coalescent (MSC) model to single-nucleotide polymorphisms. To further test the influence of different calibration strategies, we applied different calibrations to the root and to younger nodes in addition to the ones used in the previous study. Both the updated calibrations and the application of the MSC model influenced divergence time estimates, sometimes significantly. The new divergence time estimates were more ancient than in the previous study for older nodes, whereas the estimates of younger nodes were not strongly affected by the re-analyses. However, given the applied MSC method employs a simple substitution model and cannot account for changes in population size, we suggest that different analytical approaches and calibration strategies should be used in order to explore uncertainty in divergence time estimates. This study provides a valuable reference timeline for the ages of worldwide three-spined stickleback populations and emphasizes the need to embrace, rather than obscure, uncertainties around divergence time estimates.

Background: The impressive adaptive radiation of notothenioid fishes in Antarctic waters is generally thought to have been facilitated by an evolutionary key innovation, antifreeze glycoproteins, permitting the rapid evolution of more than 120 species subsequent to the Antarctic glaciation. By way of contrast, the second-most species-rich notothenioid genus, Patagonotothen, which is nested within the Antarctic clade of Notothenioidei, is almost exclusively found in the non-Antarctic waters of Patagonia. While the drivers of the diversification of Patagonotothen are currently unknown, they are unlikely to be related to antifreeze glycoproteins, given that water temperatures in Patagonia are well above freezing point. Here we performed a phylogenetic analysis based on genome-wide single nucleotide polymorphisms (SNPs) derived from restriction site-associated DNA sequencing (RADseq) in a total of twelve Patagonotothen species. Results: We present a well-supported, time-calibrated phylogenetic hypothesis including closely and distantly related outgroups, confirming the monophyly of the genus Patagonotothen with an origin approximately 3 million years ago and the paraphyly of both the sister genus Lepidonotothen and the family Notothenidae. Our phylogenomic and population genetic analyses highlight a previously unrecognized linage and provide evidence for shared genetic variation between some closely related species. We also provide a mitochondrial phylogeny showing mitonuclear discordance. Conclusions: Based on a combination of phylogenomic and population genomic approaches, we provide evidence for the existence of a new, potentially cryptic, Patagonotothen species, and demonstrate that genetic boundaries between some closely related species are diffuse, likely due to recent introgression and/or incomplete linage sorting. The detected mitonuclear discordance highlights the limitations of relying on a single locus for species barcoding. In addition, our time-calibrated phylogenetic hypothesis shows that the early burst of diversification roughly coincides with the onset of the intensification of Quaternary glacial cycles and that the rate of species accumulation may have been stepwise rather than constant. Our phylogenetic framework not only advances our understanding of the origin of a high-latitude marine radiation, but also provides the basis for the study of the ecology and life history of the genus Patagonotothen, as well as for their conservation and commercial management. Fil: Ceballos, Santiago Guillermo. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Austral de Investigaciones Científicas; Argentina. Universidad Nacional de Tierra del Fuego; Argentina Fil: Roesti, Marius. University of Basel,; Suiza. University of British Columbia; Canadá Fil: Matschiner, Michael. University of Basel,; Suiza. University of Oslo; Noruega Fil: Fernandez, Daniel Alfredo. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Austral de Investigaciones Científicas; Argentina. Universidad Nacional de Tierra del Fuego; Argentina Fil: Damerau, Malte. Johann Heinrich von Thünen-Institut ; Alemania Fil: Hanel, Reinhold. Johann Heinrich von Thünen-Institut ; Alemania Fil: Salzburger, Walter. University of Basel,; Suiza

Vertebrate vision is accomplished through a set of light-sensitive photopigments, which are located in the photoreceptors of the retina and consist of a visual opsin protein bound to a chromophore. In dim-light, vertebrates generally rely upon a single rod opsin (RH1) for obtaining visual information. By inspecting 101 fish genomes, we found that three deep-sea teleost lineages have independently expanded theirRH1gene repertoires. Amongst these, the silver spinyfin (Diretmus argenteusJohnson 1863) stands out as having the highest number of visual opsins known for animals to date (2 cone and 38 rod opsins). Spinyfins simultaneously express up to 14RH1sencoding for photopigments with different peak spectral sensitivities (λmax=448-513 nm) that cover the range of the residual daylight, as well as the bioluminescence spectrum present in the deep-sea. Our findings present novel molecular and functional evidence for the recurrent evolution of multiple rod opsin-based vision in vertebrates.SHORT ABSTRACTContrary to the single rod opsin used by most vertebrates, some fishes use multiple rod opsins for vision in the dimly lit deep-sea.

Climate warming and harvesting affect the dynamics of species across the globe through a multitude of mechanisms, including distribution changes. In fish, migrations to and distribution on spawning grounds are likely influenced by both climate warming and harvesting. The Northeast Arctic (NEA) cod (Gadus morhua) performs seasonal migrations from its feeding grounds in the Barents Sea to spawning grounds along the Norwegian coast. The distribution of cod between the spawning grounds has historically changed at decadal scales, mainly due to variable use of the northern and southern margins of the spawning area. Based on historical landing records, two major hypotheses have been put forward to explain these changes: climate and harvesting. Climate could affect the distribution through, for example, spatial habitat shifts. Harvesting could affect the distribution through impacting the demographic structure. If demographic structure is important, theory predicts increasing spawner size with migration distance. Here, we evaluate these hypotheses with modern data from a period (2000–2016) of increasing temperature and recovering stock structure. We first analyze economic data from the Norwegian fisheries to investigate geographical differences in size of spawning fish among spawning grounds, as well as interannual differences in mean latitude of spawning in relation to changes in temperature and demographic parameters. Second, we analyze genetically determined fish sampled at the spawning grounds to unambiguously separate between migratory NEA cod and potentially smaller sized coastal cod of local origin. Our results indicate smaller spawners farther away from the feeding grounds, hence not supporting the hypothesis that harvesting is a main driver for the contemporary spawning ground distribution. We find a positive correlation between annual mean spawning latitude and temperature. In conclusion, based on contemporary data, there is more support for climate compared to harvesting in shaping spawning ground distribution in this major fish stock in the North Atlantic Ocean., Climate warming and harvesting affect the dynamics of species across the globe through a multitude of mechanisms, including distribution changes of spawning fish. For the Northeast Arctic cod, two major hypotheses have been put forward to explain these changes: climate and harvesting. Here, we evaluate these hypotheses by combining two sources of modern data: economic data from the Norwegian fisheries and genetically determined individual fish sampled at the spawning grounds. Our results indicate that demographic size truncation due to harvesting is currently not the dominating factor in shaping spawning ground distribution, but our results provide support for the climate hypothesis.

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 Contact: michael.matschiner@unibas.ch Supplementary information: Supplementary data are available at Bioinformatics online.

The world’s smallest fishes belong to the genus Paedocypris. These miniature fishes are endemic to an extreme habitat: the peat swamp forests in Southeast Asia, characterized by highly acidic blackwater. This threatened habitat is home to a large array of fishes, including a number of miniaturized but also developmentally truncated species. Especially the genus Paedocypris is characterized by profound, organism-wide developmental truncation, resulting in sexually mature individuals of

Recent developments in the field of genomics have provided new and powerful insights into population structure and dynamics that are essential for the conservation of biological diversity. As a commercially highly valuable species, the yellowfin tuna (Thunnus albacares) is intensely exploited throughout its distribution in tropical oceans around the world, and is currently classified as near threatened. However, conservation efforts for this species have so far been hampered by limited knowledge of its population structure, due to incongruent results of previous investigations. Here, we use whole-genome sequencing in concert with a draft genome assembly to decipher the global population structure of the yellowfin tuna, and to investigate its demographic history. We detect significant differentiation of Atlantic and Indo-Pacific yellowfin tuna populations as well as the possibility of a third diverged yellowfin tuna group in the Arabian Sea. We further observe evidence for past population expansion as well as asymmetric gene flow from the Indo-Pacific to the Atlantic.

SummaryRecent genomic investigations have revealed hybridization to be an important source of variation, the working material of natural selection1,2. Hybridization can spur adaptive radiations3, transfer adaptive variation across species boundaries4, and generate species with novel niches5. Yet, the limits to viable hybrid genome formation are poorly understood. Here we investigated to what extent hybrid genomes are free to evolve or whether they are restricted to a specific combination of parental alleles by sequencing the genomes of four isolated island populations of the homoploid hybrid Italian sparrowPasser italiae6,7. Based on 61 Italian sparrow genomes from Crete, Corsica, Sicily and Malta, and 10 genomes of each of the parent speciesP. domesticusandP. hispaniolensis,we report that a variety of novel and fully functional hybrid genomic combinations have arisen on the different islands, with differentiation in candidate genes for beak shape and plumage colour. There are limits to successful genome fusion, however, as certain genomic regions are invariably inherited from the same parent species. These regions are overrepresented on the Z-chromosome and harbour candidate incompatibility loci, including DNA-repair and mito-nuclear genes; loci that may drive the general reduction of introgression on sex chromosomes8. Our findings demonstrate that hybridization is a potent process for generating novel variation, but variation is limited by DNA-repair and mito-nuclear genes, which play an important role in reproductive isolation and thus contribute to speciation.

The closure of the Isthmus of Panama has long been considered to be one of the best defined biogeographic calibration points for molecular divergence-time estimation. However, geological and biological evidence has recently cast doubt on the presumed timing of the initial isthmus closure around 3 Ma but has instead suggested the existence of temporary land bridges as early as the Middle or Late Miocene. The biological evidence supporting these earlier land bridges was based either on only few molecular markers or on concatenation of genome-wide sequence data, an approach that is known to result in potentially misleading branch lengths and divergence times, which could compromise the reliability of this evidence. To allow divergence-time estimation with genomic data using the more appropriate multi-species coalescent model, we here develop a new method combining the SNP-based Bayesian species-tree inference of the software SNAPP with a molecular clock model that can be calibrated with fossil or biogeographic constraints. We validate our approach with simulations and use our method to reanalyze genomic data of Neotropical army ants (Dorylinae) that previously supported divergence times of Central and South American populations before the isthmus closure around 3 Ma. Our reanalysis with the multi-species coalescent model shifts all of these divergence times to ages younger than 3 Ma, suggesting that the older estimates supporting the earlier existence of temporary land bridges were artifacts resulting at least partially from the use of concatenation. We then apply our method to a new RAD-sequencing data set of Neotropical sea catfishes (Ariidae) and calibrate their species tree with extensive information from the fossil record. We identify a series of divergences between groups of Caribbean and Pacific sea catfishes around 10 Ma, indicating that processes related to the emergence of the isthmus led to vicariant speciation already in the Late Miocene, millions of years before the final isthmus closure.

Teleost fishes comprise more than half of all vertebrate species, yet genomic data are only available for 0.2% of their diversity. Here, we present whole genome sequencing data for 66 new species of teleosts, vastly expanding the availability of genomic data for this important vertebrate group. We report on de novo assemblies based on low-coverage (9–39×) sequencing and present detailed methodology for all analyses. To facilitate further utilization of this data set, we present statistical analyses of the gene space completeness and verify the expected phylogenetic position of the sequenced genomes in a large mitogenomic context. We further present a nuclear marker set used for phylogenetic inference and evaluate each gene tree in relation to the species tree to test for homogeneity in the phylogenetic signal. Collectively, these analyses illustrate the robustness of this highly diverse data set and enable extensive reuse of the selected phylogenetic markers and the genomic data in general. This data set covers all major teleost lineages and provides unprecedented opportunities for comparative studies of teleosts. Machine-accessible metadata file describing the reported data (ISA-Tab format)