Sea anemone genomes reveal ancestral metazoan chromosomal macrosynteny

Cnidaria (sea anemones, jellyfish, corals and hydra) form a close sister group to Bilateria. Within this clade, the sea anemone Nematostella vectensis has emerged as a slow evolving model for investigating characteristics of the cnidarian-bilaterian common ancestor, which diverged near the Cambrian explosion. Here, using long read sequencing and high throughput chromosome conformation capture, we generate high quality chromosome-level genome assemblies for N. vectensis and the closely related edwardsiid sea anemone, Scolanthus callimorphus. In both cases we find a robust set of 15 chromosomes comprising a stable linkage group detectable within all major clades of sequenced cnidarian genomes. Further, both genomes show remarkable chromosomal conservation with chordates. In contrast with Bilateria, we report that extended Hox and NK gene clusters are chromosomally linked but do not retain a tight spatial conservation. Accordingly, there is a lack of evidence for topologically associated domains, which have been implicated in the evolutionary pressure to retain tight microsyntenic gene clusters. We also uncover ultra-conserved noncoding elements at levels previously undetected in non-chordate lineages. Both genomes are accessible through an actively updated genome browser and database at https://simrbase.stowers.org

Developmental constraint shaped genome evolution and erythrocyte loss in Antarctic fishes following paleoclimate change

In the frigid, oxygen-rich Southern Ocean (SO), Antarctic icefishes (Channichthyidae; Notothenioidei) evolved the ability to survive without producing erythrocytes and hemoglobin, the oxygen-transport system of virtually all vertebrates. Here, we integrate paleoclimate records with an extensive phylogenomic dataset of notothenioid fishes to understand the evolution of trait loss associated with climate change. In contrast to buoyancy adaptations in this clade, we find relaxed selection on the genetic regions controlling erythropoiesis evolved only after sustained cooling in the SO. This pattern is seen not only within icefishes but also occurred independently in other high-latitude notothenioids. We show that one species of the red-blooded dragonfish clade evolved a spherocytic anemia that phenocopies human patients with this disease via orthologous mutations. The genomic imprint of SO climate change is biased toward erythrocyte-associated conserved noncoding elements (CNEs) rather than to coding regions, which are largely preserved through pleiotropy. The drift in CNEs is specifically enriched near genes that are preferentially expressed late in erythropoiesis. Furthermore, we find that the hematopoietic marrow of icefish species retained proerythroblasts, which indicates that early erythroid development remains intact. Our results provide a framework for understanding the interactions between development and the genome in shaping the response of species to climate change.

De novo assembly of the goldfish (Carassius auratus) genome and the evolution of genes after whole-genome duplication

For over a thousand years, the common goldfish (Carassius auratus) was raised throughout Asia for food and as an ornamental pet. As a very close relative of the common carp (Cyprinus carpio), goldfish share the recent genome duplication that occurred approximately 14 million years ago in their common ancestor. The combination of centuries of breeding and a wide array of interesting body morphologies provides an exciting opportunity to link genotype to phenotype and to understand the dynamics of genome evolution and speciation. We generated a high-quality draft sequence and gene annotations of a “Wakin” goldfish using 71X PacBio long reads. The two subgenomes in goldfish retained extensive synteny and collinearity between goldfish and zebrafish. However, genes were lost quickly after the carp whole-genome duplication, and the expression of 30% of the retained duplicated gene diverged substantially across seven tissues sampled. Loss of sequence identity and/or exons determined the divergence of the expression levels across all tissues, while loss of conserved noncoding elements determined expression variance between different tissues. This assembly provides an important resource for comparative genomics and understanding the causes of goldfish variants.

Evidence for purifying selection on conserved noncoding elements in the genome of Drosophila melanogaster

Identification of regulatory regions within genomes is a key challenge for understanding the influence of functional traits on species evolution. Development of whole-genome sequencing programs in the early 2000s has been associated with the rise of comparative genomics as a powerful tool for predicting functional regions in genomes, using structural characteristics of DNA sequences. Conserved Noncoding Elements (CNEs, i.e. untranslated sequences highly similar across divergent species) were early identified as candidate regulatory regions in metazoans, plants, and fungi. CNEs have been repeatedly discriminated from mutational cold-spots, but only few studies have assessed their functional relevance at a genome-wide scale. In the present study, we used a whole-genome alignment of 27 insect species to build a coherent mapping of CNEs in the genome of Drosophila melanogaster. We then exploited polymorphism data from 48 European populations of D. melanogaster to estimate levels of nucleotide diversity and adaptive evolution in CNEs. We show here that about 35% of D. melanogaster autosomes fall into conserved noncoding regions that exhibit reduced genetic diversity and undergo purifying selection. In addition, we report six insertions of Transposable Elements (TEs) in the genome of D. melanogaster showing high levels of conservation across Drosophila species. Half of these insertions are located in untranslated transcribed regions (UTRs) of genes involved in developmental pathways and thus represent potential relics of ancient TE domestication events.

Convergent regulatory evolution and loss of flight in paleognathous birds

A core question in evolutionary biology is whether convergent phenotypic evolution is driven by convergent molecular changes in proteins or regulatory regions. We combined phylogenomic, developmental, and epigenomic analysis of 11 new genomes of paleognathous birds, including an extinct moa, to show that convergent evolution of regulatory regions, more so than protein-coding genes, is prevalent among developmental pathways associated with independent losses of flight. A Bayesian analysis of 284,001 conserved noncoding elements, 60,665 of which are corroborated as enhancers by open chromatin states during development, identified 2355 independent accelerations along lineages of flightless paleognaths, with functional consequences for driving gene expression in the developing forelimb. Our results suggest that the genomic landscape associated with morphological convergence in ratites has a substantial shared regulatory component.

CNEr: a toolkit for exploring extreme noncoding conservation

Conserved Noncoding Elements (CNEs) are elements exhibiting extreme noncoding conservation in Metazoan genomes. They cluster around developmental genes and act as long-range enhancers, yet nothing that we know about their function explains the observed conservation levels. Clusters of CNEs coincide with topologically associating domains (TADs), indicating ancient origins and stability of TAD locations. This has suggested further hypotheses about the still elusive origin of CNEs, and has provided a comparative genomics-based method of estimating the position of TADs around developmentally regulated genes in genomes where chromatin conformation capture data is missing. To enable researchers in gene regulation and chromatin biology to start deciphering this phenomenon, we developedCNEr, a R/Bioconductor toolkit for large-scale identification of CNEs and for studying their genomic properties. We applyCNErto two novel genome comparisons - fruit fly vs tsetse fly, and two sea urchin genomes - and report novel insights gained from their analysis. We also show how to reveal interesting characteristics of CNEs by coupling CNEr with existing Bioconductor packages.CNEris available at Bioconductor (https://bioconductor.org/packages/CNEr/) and maintained at github (https://github.com/ge11232002/CNEr).