Museum specimens tell the history of rhinoceroses.

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.

Drivers and dynamics of a massive adaptive radiation in cichlid fishes.

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.

Towards reliable detection of introgression in the presence of among-species rate variation.

The role of interspecific hybridization has recently seen increasing attention, especially in the context of diversification dynamics. Genomic research has now made it abundantly clear that both hybridization and introgression - the exchange of genetic material through hybridization and backcrossing - are far more common than previously thought. Besides cases of ongoing or recent genetic exchange between taxa, an increasing number of studies report "ancient introgression" - referring to results of hybridization that took place in the distant past. However, it is not clear whether commonly used methods for the detection of introgression are applicable to such old systems, given that most of these methods were originally developed for analyses at the level of populations and recently diverged species, affected by recent or ongoing genetic exchange. 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 old systems, 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 showed that some commonly applied statistical methods, including the D-statistic and certain tests based on sets of local phylogenetic trees, can produce false-positive signals of introgression between divergent taxa that have different rates of evolution. These misleading signals are caused by the presence of homoplasies occurring at different rates in different lineages. 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 along the genome, and implemented this test in the program Dsuite.

Genomic stability through time despite decades of exploitation in cod on both sides of the Atlantic.

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.

Evolution of male pregnancy associated with remodeling of canonical vertebrate immunity in seahorses and pipefishes.

A fundamental problem for the evolution of pregnancy, the most specialized form of parental investment among vertebrates, is the rejection of the nonself-embryo. Mammals achieve immunological tolerance by down-regulating both major histocompatibility complex pathways (MHC I and II). Although pregnancy has evolved multiple times independently among vertebrates, knowledge of associated immune system adjustments is restricted to mammals. All of them (except monotremata) display full internal pregnancy, making evolutionary reconstructions within the class mammalia meaningless. Here, we study the seahorse and pipefish family (syngnathids) that have evolved male pregnancy across a gradient from external oviparity to internal gestation. We assess how immunological tolerance is achieved by reconstruction of the immune gene repertoire in a comprehensive sample of 12 seahorse and pipefish genomes along the "male pregnancy" gradient together with expression patterns of key immune and pregnancy genes in reproductive tissues. We found that the evolution of pregnancy coincided with a modification of the adaptive immune system. Divergent genomic rearrangements of the MHC II pathway among fully pregnant species were identified in both genera of the syngnathids: The pipefishes (Syngnathus) displayed loss of several genes of the MHC II pathway while seahorses (Hippocampus) featured a highly divergent invariant chain (CD74). Our findings suggest that a trade-off between immunological tolerance and embryo rejection accompanied the evolution of unique male pregnancy. That pipefishes survive in an ocean of microbes without one arm of the adaptive immune defense suggests a high degree of immunological flexibility among vertebrates, which may advance our understanding of immune-deficiency diseases.

Divergence and hybridization in sea turtles: Inferences from genome data show evidence of ancient gene flow between species.

Reconstructing past events of hybridization and population size changes are required to understand speciation mechanisms and current patterns of genetic diversity, and ultimately contribute to species' conservation. Sea turtles are ancient species currently facing anthropogenic threats including climate change, fisheries, and illegal hunting. Five of the seven extant sea turtle species are known to currently hybridize, especially along the Brazilian coast where some populations can have ~32%-42% of hybrids. Although frequently observed today, it is not clear what role hybridization plays in the evolutionary diversification of this group of reptiles. In this study, we generated whole genome resequencing data of the five globally distributed sea turtle species to estimate a calibrated phylogeny and the population size dynamics, and to understand the role of hybridization in shaping the genomes of these ancient species. Our results reveal discordant species divergence dates between mitochondrial and nuclear genomes, with a high frequency of conflicting trees throughout the nuclear genome suggesting that some sea turtle species frequently hybridized in the past. The reconstruction of the species' demography showed a general decline in effective population sizes with no signs of recovery, except for the leatherback sea turtle. Furthermore, we discuss the influence of reference bias in our estimates. We show long-lasting ancestral gene flow events within Chelonioidea that continued for millions of years after initial divergence. Speciation with gene flow is a common pattern in marine species, and it raises questions whether current hybridization events should be considered as a part of these species' evolutionary history or a conservation issue.

Estimating uncertainty in divergence times among three-spined stickleback clades using the multispecies coalescent.

Incomplete lineage sorting (ILS) can lead to biased divergence time estimates. To explore if and how ILS has influenced the results of a recent 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, since 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.

BEAST 2.5: An advanced software platform for Bayesian evolutionary analysis.

Elaboration of Bayesian phylogenetic inference methods has continued at pace in recent years with major new advances in nearly all aspects of the joint modelling of evolutionary data. It is increasingly appreciated that some evolutionary questions can only be adequately answered by combining evidence from multiple independent sources of data, including genome sequences, sampling dates, phenotypic data, radiocarbon dates, fossil occurrences, and biogeographic range information among others. Including all relevant data into a single joint model is very challenging both conceptually and computationally. Advanced computational software packages that allow robust development of compatible (sub-)models which can be composed into a full model hierarchy have played a key role in these developments. Developing such software frameworks is increasingly a major scientific activity in its own right, and comes with specific challenges, from practical software design, development and engineering challenges to statistical and conceptual modelling challenges. BEAST 2 is one such computational software platform, and was first announced over 4 years ago. Here we describe a series of major new developments in the BEAST 2 core platform and model hierarchy that have occurred since the first release of the software, culminating in the recent 2.5 release.

Phylogenomics of an extra-Antarctic notothenioid radiation reveals a previously unrecognized lineage and diffuse species boundaries.

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.

Ticket to spawn: Combining economic and genetic data to evaluate the effect of climate and demographic structure on spawning distribution in Atlantic cod.

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.

Bayesian Divergence-Time Estimation with Genome-Wide Single-Nucleotide Polymorphism Data of Sea Catfishes (Ariidae) Supports Miocene Closure of the Panamanian Isthmus.

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.

Millions of Years Behind: Slow Adaptation of Ruminants to Grasslands.

The Late Cretaceous appearance of grasses, followed by the Cenozoic advancement of grasslands as dominant biomes, has contributed to the evolution of a range of specialized herbivores adapted to new diets, as well as to increasingly open and arid habitats. Many mammals including ruminants, the most diversified ungulate suborder, evolved high-crowned (hypsodont) teeth as an adaptation to tooth-wearing diets and habitats. The impact of different causes of tooth wear is still a matter of debate, and the temporal pattern of hypsodonty evolution in relation to the evolution of grasslands remains unclear. We present an improved time-calibrated molecular phylogeny of Cetartiodactyla, with phylogenetic reconstruction of ancestral ruminant diets and habitats, based on characteristics of extant taxa. Using this timeline, as well as the fossil record of grasslands, we conduct phylogenetic comparative analyses showing that hypsodonty in ruminants evolved as an adaptation to both diet and habitat. Our results demonstrate a slow, perhaps constrained, evolution of hypsodonty toward estimated optimal states, excluding the possibility of immediate adaptation. This augments recent findings that slow adaptation is not uncommon on million-year time scales.

Disentangling Incomplete Lineage Sorting and Introgression to Refine Species-Tree Estimates for Lake Tanganyika Cichlid Fishes.

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.

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

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.

Fitchi: haplotype genealogy graphs based on the Fitch algorithm.

UNLABELLED: : In population genetics and phylogeography, haplotype genealogy graphs are important tools for the visualization of population structure based on sequence data. In this type of graph, node sizes are often drawn in proportion to haplotype frequencies and edge lengths represent the minimum number of mutations separating adjacent nodes. I here present Fitchi, a new program that produces publication-ready haplotype genealogy graphs based on the Fitch algorithm. AVAILABILITY AND IMPLEMENTATION: http://www.evoinformatics.eu/fitchi.htm CONTACT: : michaelmatschiner@mac.com SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Genomics of speciation and introgression in Princess cichlid fishes from Lake Tanganyika.

How variation in the genome translates into biological diversity and new species originate has endured as the mystery of mysteries in evolutionary biology. African cichlid fishes are prime model systems to address speciation-related questions for their remarkable taxonomic and phenotypic diversity, and the possible role of gene flow in this process. Here, we capitalize on genome sequencing and phylogenomic analyses to address the relative impacts of incomplete lineage sorting, introgression and hybrid speciation in the Neolamprologus savoryi-complex (the 'Princess cichlids') from Lake Tanganyika. We present a time-calibrated species tree based on whole-genome sequences and provide strong evidence for incomplete lineage sorting in the early phases of diversification and multiple introgression events affecting different stages. Importantly, we find that the Neolamprologus chromosomes show centre-to-periphery biases in nucleotide diversity, sequence divergence, GC content, incomplete lineage sorting and rates of introgression, which are likely modulated by recombination density and linked selection. The detection of heterogeneous genomic landscapes has strong implications on the genomic mechanisms involved in speciation. Collinear chromosomal regions can be protected from gene flow and harbour incompatibility genes if they reside in lowly recombining regions, and coupling can evolve between nonphysically linked genomic regions (chromosome centres in particular). Simultaneously, higher recombination towards chromosome peripheries makes these more dynamic, evolvable regions where adaptation polymorphisms have a fertile ground. Hence, differences in genome architecture could explain the levels of taxonomic and phenotypic diversity seen in taxa with collinear genomes and might have contributed to the spectacular cichlid diversity observed today.

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

The species-flocks of cichlid fishes in the East African Great Lakes Victoria, Malawi and Tanganyika constitute the most diverse extant adaptive radiations in vertebrates. Lake Tanganyika, the oldest of the lakes, harbors the morphologically and genetically most diverse assemblage of cichlids and contains the highest number of endemic cichlid genera of all African lakes. Based on morphological grounds, the Tanganyikan cichlid species have been grouped into 12-16 distinct lineages, so-called tribes. While the monophyly of most of the tribes is well established, the phylogenetic relationships among the tribes remain largely elusive. Here, we present a new tribal level phylogenetic hypothesis for the cichlid fishes of Lake Tanganyika that is based on the so far largest set of nuclear markers and a total alignment length of close to 18kb. Using next-generation amplicon sequencing with the 454 pyrosequencing technology, we compiled a dataset consisting of 42 nuclear loci in 45 East African cichlid species, which we subjected to maximum likelihood and Bayesian inference phylogenetic analyses. We analyzed the entire concatenated dataset and each marker individually, and performed a Bayesian concordance analysis and gene tree discordance tests. Overall, we find strong support for a position of the Oreochromini, Boulengerochromini, Bathybatini and Trematocarini outside of a clade combining the substrate spawning Lamprologini and the mouthbrooding tribes of the 'H-lineage', which are both strongly supported to be monophyletic. The Eretmodini are firmly placed within the 'H-lineage', as sister-group to the most species-rich tribe of cichlids, the Haplochromini. The phylogenetic relationships at the base of the 'H-lineage' received less support, which is likely due to high speciation rates in the early phase of the radiation. Discordance among gene trees and marker sets further suggests the occurrence of past hybridization and/or incomplete lineage sorting in the cichlid fishes of Lake Tanganyika.

Hapsolutely: a user-friendly tool integrating haplotype phasing, network construction, and haploweb calculation.

Motivation: Haplotype networks are a routine approach to visualize relationships among alleles. Such visual analysis of single-locus data is still of importance, especially in species diagnosis and delimitation, where a limited amount of sequence data usually are available and sufficient, along with other datasets in the framework of integrative taxonomy. In diploid organisms, this often requires separating (phasing) sequences with heterozygotic positions, and typically separate programs are required for phasing, reformatting of input files, and haplotype network construction. We therefore developed Hapsolutely, a user-friendly program with an ergonomic graphical user interface that integrates haplotype phasing from single-locus sequences with five approaches for network/genealogy reconstruction. Results: Among the novel options implemented, Hapsolutely integrates phasing and graphical reconstruction steps of haplotype networks, supports input of species partition data in the common SPART and SPART-XML formats, and calculates and visualizes haplowebs and fields for recombination, thus allowing graphical comparison of allele distribution and allele sharing among subsets for the purpose of species delimitation. The new tool has been specifically developed with a focus on the workflow in alpha-taxonomy, where exploring fields for recombination across alternative species partitions may help species delimitation. Availability and implementation: Hapsolutely is written in Python, and integrates code from Phase, SeqPHASE, and PopART in C++ and Haxe. Compiled stand-alone executables for MS Windows and Mac OS along with a detailed manual can be downloaded from https://www.itaxotools.org; the source code is openly available on GitHub (https://github.com/iTaxoTools/Hapsolutely).

First chromosome-level genome assembly of a ribbon worm from the Hoplonemertea clade, Emplectonema gracile, and its structural annotation.

Genome-wide information has so far been unavailable for ribbon worms of the clade Hoplonemertea, the most species-rich class within the phylum Nemertea. While species within Pilidiophora, the sister clade of Hoplonemertea, possess a pilidium larval stage and lack stylets on their proboscis, Hoplonemertea species have a planuliform larva and are armed with stylets employed for the injection of toxins into their prey. To further compare these developmental, physiological, and behavioral differences from a genomic perspective, the availability of a reference genome for a Hoplonemertea species is crucial. Such data will be highly useful for future investigations towards a better understanding of molecular ecology, venom evolution, and regeneration not only in Nemertea, but also in other marine invertebrate phyla. To this end, we herein present the annotated chromosome-level genome assembly for Emplectonema gracile (Nemertea; Hoplonemertea; Monostilifera; Emplectonematidae), an easily collected nemertean well-suited for laboratory experimentation. The genome has an assembly size of 157.9 Mbp. Hi-C scaffolding yielded chromosome level scaffolds, with a scaffold N50 of 10.0 Mbp and a score of 95.1% for complete BUSCO genes found as a single copy. Annotation predicted 20,684 protein-coding genes. The high-quality reference genome reaches an Earth BioGenome standard level of 7.C.Q50.

The genetics of niche-specific behavioral tendencies in an adaptive radiation of cichlid fishes.

Behavior is critical for animal survival and reproduction, and possibly for diversification and evolutionary radiation. However, the genetics behind adaptive variation in behavior are poorly understood. In this work, we examined a fundamental and widespread behavioral trait, exploratory behavior, in one of the largest adaptive radiations on Earth, the cichlid fishes of Lake Tanganyika. By integrating quantitative behavioral data from 57 cichlid species (702 wild-caught individuals) with high-resolution ecomorphological and genomic information, we show that exploratory behavior is linked to macrohabitat niche adaptations in Tanganyikan cichlids. Furthermore, we uncovered a correlation between the genotypes at a single-nucleotide polymorphism upstream of the AMPA glutamate-receptor regulatory gene cacng5b and variation in exploratory tendency. We validated this association using behavioral predictions with a neural network approach and CRISPR-Cas9 genome editing.