Amphioxus functional genomics and the origins of vertebrate gene regulation

The Catalan initiative for the Earth BioGenome Project: contributing local data to global biodiversity genomics

The Catalan Initiative for the Earth BioGenome Project (CBP) is an EBP-affiliated project network aimed at sequencing the genome of the >40 000 eukaryotic species estimated to live in the Catalan-speaking territories (Catalan Linguistic Area, CLA). These territories represent a biodiversity hotspot. While covering less than 1% of Europe, they are home to about one fourth of all known European eukaryotic species. These include a high proportion of endemisms, many of which are threatened. This trend is likely to get worse as the effects of global change are expected to be particularly severe across the Mediterranean Basin, particularly in freshwater ecosystems and mountain areas. Following the EBP model, the CBP is a networked organization that has been able to engage many scientific and non-scientific partners. In the pilot phase, the genomes of 52 species are being sequenced. As a case study in biodiversity conservation, we highlight the genome of the Balearic shearwater Puffinus mauretanicus, sequenced under the CBP umbrella.

Annelid methylomes reveal ancestral developmental and aging-associated epigenetic erosion across Bilateria

BACKGROUND

DNA methylation in the form of 5-methylcytosine (5mC) is the most abundant base modification in animals. However, 5mC levels vary widely across taxa. While vertebrate genomes are hypermethylated, in most invertebrates, 5mC concentrates on constantly and highly transcribed genes (gene body methylation; GbM) and, in some species, on transposable elements (TEs), a pattern known as "mosaic". Yet, the role and developmental dynamics of 5mC and how these explain interspecies differences in DNA methylation patterns remain poorly understood, especially in Spiralia, a large clade of invertebrates comprising nearly half of the animal phyla.

RESULTS

Here, we generate base-resolution methylomes for three species with distinct genomic features and phylogenetic positions in Annelida, a major spiralian phylum. All possible 5mC patterns occur in annelids, from typical invertebrate intermediate levels in a mosaic distribution to hypermethylation and methylation loss. GbM is common to annelids with 5mC, and methylation differences across species are explained by taxon-specific transcriptional dynamics or the presence of intronic TEs. Notably, the link between GbM and transcription decays during development, alongside a gradual and global, age-dependent demethylation in adult stages. Additionally, reducing 5mC levels with cytidine analogs during early development impairs normal embryogenesis and reactivates TEs in the annelid Owenia fusiformis.

CONCLUSIONS

Our study indicates that global epigenetic erosion during development and aging is an ancestral feature of bilateral animals. However, the tight link between transcription and gene body methylation is likely more important in early embryonic stages, and 5mC-mediated TE silencing probably emerged convergently across animal lineages.

The developmental and evolutionary characteristics of transcription factor binding site clustered regions based on an explainable machine learning model

Gene expression is temporally and spatially regulated by the interaction of transcription factors (TFs) and cis-regulatory elements (CREs). The uneven distribution of TF binding sites across the genome poses challenges in understanding how this distribution evolves to regulate spatio-temporal gene expression and consequent heritable phenotypic variation. In this study, chromatin accessibility profiles and gene expression profiles were collected from several species including mammals (human, mouse, bovine), fish (zebrafish and medaka), and chicken. Transcription factor binding sites clustered regions (TFCRs) at different embryonic stages were characterized to investigate regulatory evolution. The study revealed dynamic changes in TFCR distribution during embryonic development and species evolution. The synchronization between TFCR complexity and gene expression was assessed across species using RegulatoryScore. Additionally, an explainable machine learning model highlighted the importance of the distance between TFCR and promoter in the coordinated regulation of TFCRs on gene expression. Our results revealed the developmental and evolutionary dynamics of TFCRs during embryonic development from fish, chicken to mammals. These data provide valuable resources for exploring the relationship between transcriptional regulation and phenotypic differences during embryonic development.

Molecular Dating of the Teleost Whole Genome Duplication (3R) Is Compatible With the Expectations of Delayed Rediploidization

Vertebrate evolution has been punctuated by three whole genome duplication events that have been implicated causally in phenotypic evolution, from the origin of phenotypic novelties to explosive diversification. Arguably, the most dramatic of these is the 3R whole genome duplication event associated with the origin of teleost fishes which comprise more than half of all living vertebrate species. However, tests of a causal relationship between whole genome duplication and teleost diversification have proven difficult due to the challenge of establishing the timing of these phenomena. Here we show, based on molecular clock dating of concatenated gene alignments, that the 3R whole genome duplication event occurred in the early-middle Permian (286.18 to 267.20 million years ago; Ma), 52.02 to 12.84 million years (Myr) before the divergence of crown-teleosts in the latest Permian-earliest Late Triassic (254.36 to 234.16 Ma) and long before the major pulses of teleost diversification in Ostariophysi and Percomorpha (56.37 to 100.17 Myr and at least 139.24 to 183.29 Myr later, respectively). The extent of this temporal gap between putative cause and effect precludes 3R as a deterministic driver of teleost diversification. However, these age constraints remain compatible with the expectations of a prolonged rediploidization process following whole genome duplication which, through the effects of chromosome rearrangement and gene loss, remains a viable mechanism to explain the evolution of teleost novelties and diversification.

Conservation of cis-Regulatory Syntax Underlying Deuterostome Gastrulation

Throughout embryonic development, the shaping of the functional and morphological characteristics of embryos is orchestrated by an intricate interaction between transcription factors and cis-regulatory elements. In this study, we conducted a comprehensive analysis of deuterostome cis-regulatory landscapes during gastrulation, focusing on four paradigmatic species: the echinoderm Strongylocentrotus purpuratus, the cephalochordate Branchiostoma lanceolatum, the urochordate Ciona intestinalis, and the vertebrate Danio rerio. Our approach involved comparative computational analysis of ATAC-seq datasets to explore the genome-wide blueprint of conserved transcription factor binding motifs underlying gastrulation. We identified a core set of conserved DNA binding motifs associated with 62 known transcription factors, indicating the remarkable conservation of the gastrulation regulatory landscape across deuterostomes. Our findings offer valuable insights into the evolutionary molecular dynamics of embryonic development, shedding light on conserved regulatory subprograms and providing a comprehensive perspective on the conservation and divergence of gene regulation underlying the gastrulation process.

Evolution of diapause in the African turquoise killifish by remodeling the ancient gene regulatory landscape

Suspended animation states allow organisms to survive extreme environments. The African turquoise killifish has evolved diapause as a form of suspended development to survive a complete drought. However, the mechanisms underlying the evolution of extreme survival states are unknown. To understand diapause evolution, we performed integrative multi-omics (gene expression, chromatin accessibility, and lipidomics) in the embryos of multiple killifish species. We find that diapause evolved by a recent remodeling of regulatory elements at very ancient gene duplicates (paralogs) present in all vertebrates. CRISPR-Cas9-based perturbations identify the transcription factors REST/NRSF and FOXOs as critical for the diapause gene expression program, including genes involved in lipid metabolism. Indeed, diapause shows a distinct lipid profile, with an increase in triglycerides with very-long-chain fatty acids. Our work suggests a mechanism for the evolution of complex adaptations and offers strategies to promote long-term survival by activating suspended animation programs in other species.

Amphioxus (Branchiostoma lanceolatum) in the North Adriatic Sea: ecological observations and spawning behavior

The European amphioxus (Branchiostoma lanceolatum) is a member of the chordate subphylum Cephalochordata, and, as such, a key model organism for providing insights into the origin and evolution of vertebrates. Despite its significance and global distribution, detailed characterizations of natural populations of cephalochordates are still very limited. This study investigates the abundance, habitat, and spawning behavior of amphioxus in the North Adriatic Sea. Across 32 sampled sites, adult amphioxus were consistently present, reaching densities exceeding 300 individuals m- 2. DNA barcoding confirmed the species as B. lanceolatum, and environmental analyses revealed an amphioxus preference for slightly gravelly sand with low silt content and a correlation between amphioxus density and the presence of specific macroinvertebrate taxa. Remarkably, the amphioxus population was breeding in early spring and possibly late fall, in contrast to the typical late spring/early summer spawning season described for other populations of European amphioxus. Amphioxus adults kept in captivity maintained the spawning seasonality of their place of origin, suggesting the possibility of extending the overall spawning season of European amphioxus in laboratory settings by exploiting populations from diverse geographic origins. This study thus expands our understanding of B. lanceolatum ecology and reproduction in the Mediterranean Sea, emphasizing the role of the North Adriatic Sea as a substantial reservoir.

Single-cell phylotranscriptomics of developmental and cell type evolution

Single-cell phylotranscriptomics is an emerging tool to reveal the molecular and cellular mechanisms of evolution. We summarize its utility in studying the hourglass pattern of ontogenetic evolution and for understanding the evolutionary history of cell types. The developmental hourglass model suggests that the mid-embryonic stage is the most conserved period of development across species, which is supported by morphological and molecular studies. Single-cell phylotranscriptomic analysis has revealed previously underappreciated heterogeneity in transcriptome ages among lineages and cell types throughout development, and has identified the lineages and tissues that drive the whole-organism hourglass pattern. Single-cell transcriptome age analyses also provide important insights into the origin of germ layers, the different selective forces on tissues during adaptation, and the evolutionary relationships between cell types.

Evolution of tissue-specific expression of ancestral genes across vertebrates and insects

Regulation of gene expression is arguably the main mechanism underlying the phenotypic diversity of tissues within and between species. Here we assembled an extensive transcriptomic dataset covering 8 tissues across 20 bilaterian species and performed analyses using a symmetric phylogeny that allowed the combined and parallel investigation of gene expression evolution between vertebrates and insects. We specifically focused on widely conserved ancestral genes, identifying strong cores of pan-bilaterian tissue-specific genes and even larger groups that diverged to define vertebrate and insect tissues. Systematic inferences of tissue-specificity gains and losses show that nearly half of all ancestral genes have been recruited into tissue-specific transcriptomes. This occurred during both ancient and, especially, recent bilaterian evolution, with several gains being associated with the emergence of unique phenotypes (for example, novel cell types). Such pervasive evolution of tissue specificity was linked to gene duplication coupled with expression specialization of one of the copies, revealing an unappreciated prolonged effect of whole-genome duplications on recent vertebrate evolution.

An amphioxus neurula stage cell atlas supports a complex scenario for the emergence of vertebrate head mesoderm

The emergence of new structures can often be linked to the evolution of novel cell types that follows the rewiring of developmental gene regulatory subnetworks. Vertebrates are characterized by a complex body plan compared to the other chordate clades and the question remains of whether and how the emergence of vertebrate morphological innovations can be related to the appearance of new embryonic cell populations. We previously proposed, by studying mesoderm development in the cephalochordate amphioxus, a scenario for the evolution of the vertebrate head mesoderm. To further test this scenario at the cell population level, we used scRNA-seq to construct a cell atlas of the amphioxus neurula, stage at which the main mesodermal compartments are specified. Our data allowed us to validate the presence of a prechordal-plate like territory in amphioxus. Additionally, the transcriptomic profile of somite cell populations supports the homology between specific territories of amphioxus somites and vertebrate cranial/pharyngeal and lateral plate mesoderm. Finally, our work provides evidence that the appearance of the specific mesodermal structures of the vertebrate head was associated to both segregation of pre-existing cell populations, and co-option of new genes for the control of myogenesis.

Rewiring of the epigenome and chromatin architecture by exogenously induced retinoic acid signaling during zebrafish embryonic development

Retinoic acid (RA) is the ligand of RA receptors (RARs), transcription factors that bind to RA response elements. RA signaling is required for multiple processes during embryonic development, including body axis extension, hindbrain antero-posterior patterning and forelimb bud initiation. Although some RA target genes have been identified, little is known about the genome-wide effects of RA signaling during in vivo embryonic development. Here, we stimulate the RA pathway by treating zebrafish embryos with all-trans-RA (atRA) and use a combination of RNA-seq, ATAC-seq, ChIP-seq and HiChIP to gain insight into the molecular mechanisms by which exogenously induced RA signaling controls gene expression. We find that RA signaling is involved in anterior/posterior patterning, central nervous system development, and the transition from pluripotency to differentiation. AtRA treatment also alters chromatin accessibility during early development and promotes chromatin binding of RARαa and the RA targets Hoxb1b, Meis2b and Sox3, which cooperate in central nervous system development. Finally, we show that exogenous RA induces a rewiring of chromatin architecture, with alterations in chromatin 3D interactions involving target genes. Altogether, our findings identify genome-wide targets of RA signaling and provide a molecular mechanism by which developmental signaling pathways regulate target gene expression by altering chromatin topology.

Evolutionary analysis of gene ages across TADs associates chromatin topology with whole-genome duplications

Topologically associated domains (TADs) are interaction subnetworks of chromosomal regions in 3D genomes. TAD boundaries frequently coincide with genome breaks while boundary deletion is under negative selection, suggesting that TADs may facilitate genome rearrangements and evolution. We show that genes co-localize by evolutionary age in humans and mice, resulting in TADs having different proportions of younger and older genes. We observe a major transition in the age co-localization patterns between the genes born during vertebrate whole-genome duplications (WGDs) or before and those born afterward. We also find that genes recently duplicated in primates and rodents are more frequently essential when they are located in old-enriched TADs and interact with genes that last duplicated during the WGD. Therefore, the evolutionary relevance of recent genes may increase when located in TADs with established regulatory networks. Our data suggest that TADs could play a role in organizing ancestral functions and evolutionary novelty.

Quantitative proteome dynamics across embryogenesis in a model chordate

The evolution of gene expression programs underlying the development of vertebrates remains poorly characterized. Here, we present a comprehensive proteome atlas of the model chordate Ciona, covering eight developmental stages and ∼7,000 translated genes, accompanied by a multi-omics analysis of co-evolution with the vertebrate Xenopus. Quantitative proteome comparisons argue against the widely held hourglass model, based solely on transcriptomic profiles, whereby peak conservation is observed during mid-developmental stages. Our analysis reveals maximal divergence at these stages, particularly gastrulation and neurulation. Together, our work provides a valuable resource for evaluating conservation and divergence of multi-omics profiles underlying the diversification of vertebrates.

Gene networks and the evolution of olfactory organs, eyes, hair cells and motoneurons: a view encompassing lancelets, tunicates and vertebrates

Key developmental pathways and gene networks underlie the formation of sensory cell types and structures involved in chemosensation, vision and mechanosensation, and of the efferents these sensory inputs can activate. We describe similarities and differences in these pathways and gene networks in selected species of the three main chordate groups, lancelets, tunicates, and vertebrates, leading to divergent development of olfactory receptors, eyes, hair cells and motoneurons. The lack of appropriately posited expression of certain transcription factors in lancelets and tunicates prevents them from developing vertebrate-like olfactory receptors and eyes, although they generate alternative structures for chemosensation and vision. Lancelets and tunicates lack mechanosensory cells associated with the sensation of acoustic stimuli, but have gravisensitive organs and ciliated epidermal sensory cells that may (and in some cases clearly do) provide mechanosensation and thus the capacity to respond to movement relative to surrounding water. Although functionally analogous to the vertebrate vestibular apparatus and lateral line, homology is questionable due to differences in the expression of the key transcription factors Neurog and Atoh1/7, on which development of vertebrate hair cells depends. The vertebrate hair cell-bearing inner ear and lateral line thus likely represent major evolutionary advances specific to vertebrates. Motoneurons develop in vertebrates under the control of the ventral signaling molecule hedgehog/sonic hedgehog (Hh,Shh), against an opposing inhibitory effect mediated by dorsal signaling molecules. Many elements of Shh-signaling and downstream genes involved in specifying and differentiating motoneurons are also exhibited by lancelets and tunicates, but the repertoire of MNs in vertebrates is broader, indicating greater diversity in motoneuron differentiation programs.

Extensive DNA methylome rearrangement during early lamprey embryogenesis

DNA methylation (5mC) is a repressive gene regulatory mark widespread in vertebrate genomes, yet the developmental dynamics in which 5mC patterns are established vary across species. While mammals undergo two rounds of global 5mC erasure, teleosts, for example, exhibit localized maternal-to-paternal 5mC remodeling. Here, we studied 5mC dynamics during the embryonic development of sea lamprey, a jawless vertebrate which occupies a critical phylogenetic position as the sister group of the jawed vertebrates. We employed 5mC quantification in lamprey embryos and tissues, and discovered large-scale maternal-to-paternal epigenome remodeling that affects ~30% of the embryonic genome and is predominantly associated with partially methylated domains. We further demonstrate that sequences eliminated during programmed genome rearrangement (PGR), are hypermethylated in sperm prior to the onset of PGR. Our study thus unveils important insights into the evolutionary origins of vertebrate 5mC reprogramming, and how this process might participate in diverse developmental strategies.

Hagfish genome elucidates vertebrate whole-genome duplication events and their evolutionary consequences

Polyploidy or whole-genome duplication (WGD) is a major event that drastically reshapes genome architecture and is often assumed to be causally associated with organismal innovations and radiations. The 2R hypothesis suggests that two WGD events (1R and 2R) occurred during early vertebrate evolution. However, the timing of the 2R event relative to the divergence of gnathostomes (jawed vertebrates) and cyclostomes (jawless hagfishes and lampreys) is unresolved and whether these WGD events underlie vertebrate phenotypic diversification remains elusive. Here we present the genome of the inshore hagfish, Eptatretus burgeri. Through comparative analysis with lamprey and gnathostome genomes, we reconstruct the early events in cyclostome genome evolution, leveraging insights into the ancestral vertebrate genome. Genome-wide synteny and phylogenetic analyses support a scenario in which 1R occurred in the vertebrate stem-lineage during the early Cambrian, and 2R occurred in the gnathostome stem-lineage, maximally in the late Cambrian-earliest Ordovician, after its divergence from cyclostomes. We find that the genome of stem-cyclostomes experienced an additional independent genome triplication. Functional genomic and morphospace analyses demonstrate that WGD events generally contribute to developmental evolution with similar changes in the regulatory genome of both vertebrate groups. However, appreciable morphological diversification occurred only in the gnathostome but not in the cyclostome lineage, calling into question the general expectation that WGDs lead to leaps of bodyplan complexity.

The hagfish genome and the evolution of vertebrates

As the only surviving lineages of jawless fishes, hagfishes and lampreys provide a crucial window into early vertebrate evolution1-3. Here we investigate the complex history, timing and functional role of genome-wide duplications4-7 and programmed DNA elimination8,9 in vertebrates in the light of a chromosome-scale genome sequence for the brown hagfish Eptatretus atami. Combining evidence from syntenic and phylogenetic analyses, we establish a comprehensive picture of vertebrate genome evolution, including an auto-tetraploidization (1RV) that predates the early Cambrian cyclostome-gnathostome split, followed by a mid-late Cambrian allo-tetraploidization (2RJV) in gnathostomes and a prolonged Cambrian-Ordovician hexaploidization (2RCY) in cyclostomes. Subsequently, hagfishes underwent extensive genomic changes, with chromosomal fusions accompanied by the loss of genes that are essential for organ systems (for example, genes involved in the development of eyes and in the proliferation of osteoclasts); these changes account, in part, for the simplification of the hagfish body plan1,2. Finally, we characterize programmed DNA elimination in hagfish, identifying protein-coding genes and repetitive elements that are deleted from somatic cell lineages during early development. The elimination of these germline-specific genes provides a mechanism for resolving genetic conflict between soma and germline by repressing germline and pluripotency functions, paralleling findings in lampreys10,11. Reconstruction of the early genomic history of vertebrates provides a framework for further investigations of the evolution of cyclostomes and jawed vertebrates.

Sea lamprey enlightens the origin of the coupling of retinoic acid signaling to vertebrate hindbrain segmentation

Retinoic acid (RA) is involved in antero-posterior patterning of the chordate body axis and, in jawed vertebrates, has been shown to play a major role at multiple levels of the gene regulatory network (GRN) regulating hindbrain segmentation. Knowing when and how RA became coupled to the core hindbrain GRN is important for understanding how ancient signaling pathways and patterning genes can evolve and generate diversity. Hence, we investigated the link between RA signaling and hindbrain segmentation in the sea lamprey Petromyzon marinus, an important jawless vertebrate model providing clues to decipher ancestral vertebrate features. Combining genomics, gene expression, and functional analyses of major components involved in RA synthesis (Aldh1as) and degradation (Cyp26s), we demonstrate that RA signaling is coupled to hindbrain segmentation in lamprey. Thus, the link between RA signaling and hindbrain segmentation is a pan vertebrate feature of the hindbrain and likely evolved at the base of vertebrates.

Developmental hourglass: Verification by numerical evolution and elucidation by dynamical-systems theory

Determining the general laws between evolution and development is a fundamental biological challenge. Developmental hourglasses have attracted increased attention as candidates for such laws, but the necessity of their emergence remains elusive. We conducted evolutionary simulations of developmental processes to confirm the emergence of the developmental hourglass and unveiled its establishment. We considered organisms consisting of cells containing identical gene networks that control morphogenesis and evolved them under selection pressure to induce more cell types. By computing the similarity between the spatial patterns of gene expression of two species that evolved from a common ancestor, a developmental hourglass was observed, that is, there was a correlation peak in the intermediate stage of development. The fraction of pleiotropic genes increased, whereas the variance in individuals decreased, consistent with previous experimental reports. Reduction of the unavoidable variance by initial or developmental noise, essential for survival, was achieved up to the hourglass bottleneck stage, followed by diversification in developmental processes, whose timing is controlled by the slow expression dynamics conserved among organisms sharing the hourglass. This study suggests why developmental hourglasses are observed within a certain phylogenetic range of species.

Evolution of the gene regulatory network of body axis by enhancer hijacking in amphioxus

A central goal of evolutionary developmental biology is to decipher the evolutionary pattern of gene regulatory networks (GRNs) that control embryonic development, and the mechanism underlying GRNs evolution. The Nodal signaling that governs the body axes of deuterostomes exhibits a conserved GRN orchestrated principally by Nodal, Gdf1/3, and Lefty. Here we show that this GRN has been rewired in cephalochordate amphioxus. We found that while the amphioxus Gdf1/3 ortholog exhibited nearly no embryonic expression, its duplicate Gdf1/3-like, linked to Lefty, was zygotically expressed in a similar pattern as Lefty. Consistent with this, while Gdf1/3-like mutants showed defects in axial development, Gdf1/3 mutants did not. Further transgenic analyses showed that the intergenic region between Gdf1/3-like and Lefty could drive reporter gene expression as that of the two genes. These results indicated that Gdf1/3-like has taken over the axial development role of Gdf1/3 in amphioxus, possibly through hijacking Lefty enhancers. We finally demonstrated that, to compensate for the loss of maternal Gdf1/3 expression, Nodal has become an indispensable maternal factor in amphioxus and its maternal mutants caused axial defects as Gdf1/3-like mutants. We therefore demonstrated a case that the evolution of GRNs could be triggered by enhancer hijacking events. This pivotal event has allowed the emergence of a new GRN in extant amphioxus, presumably through a stepwise process. In addition, the co-expression of Gdf1/3-like and Lefty achieved by a shared regulatory region may have provided robustness during body axis formation, which provides a selection-based hypothesis for the phenomena called developmental system drift.

Topological structures and syntenic conservation in sea anemone genomes

There is currently little information about the evolution of gene clusters, genome architectures and karyotypes in early branching animals. Slowly evolving anthozoan cnidarians can be particularly informative about the evolution of these genome features. Here we report chromosome-level genome assemblies of two related anthozoans, the sea anemones Nematostella vectensis and Scolanthus callimorphus. We find a robust set of 15 chromosomes with a clear one-to-one correspondence between the two species. Both genomes show chromosomal conservation, allowing us to reconstruct ancestral cnidarian and metazoan chromosomal blocks, consisting of at least 19 and 16 ancestral linkage groups, respectively. We show that, in contrast to Bilateria, the Hox and NK clusters of investigated cnidarians are largely disintegrated, despite the presence of staggered hox/gbx expression in Nematostella. This loss of microsynteny conservation may be facilitated by shorter distances between cis-regulatory sequences and their cognate transcriptional start sites. We find no clear evidence for topologically associated domains, suggesting fundamental differences in long-range gene regulation compared to vertebrates. These data suggest that large sets of ancestral metazoan genes have been retained in ancestral linkage groups of some extant lineages; yet, higher order gene regulation with associated 3D architecture may have evolved only after the cnidarian-bilaterian split.

Evolution of the Spider Homeobox Gene Repertoire by Tandem and Whole Genome Duplication

Gene duplication generates new genetic material that can contribute to the evolution of gene regulatory networks and phenotypes. Duplicated genes can undergo subfunctionalization to partition ancestral functions and/or neofunctionalization to assume a new function. We previously found there had been a whole genome duplication (WGD) in an ancestor of arachnopulmonates, the lineage including spiders and scorpions but excluding other arachnids like mites, ticks, and harvestmen. This WGD was evidenced by many duplicated homeobox genes, including two Hox clusters, in spiders. However, it was unclear which homeobox paralogues originated by WGD versus smaller-scale events such as tandem duplications. Understanding this is a key to determining the contribution of the WGD to arachnopulmonate genome evolution. Here we characterized the distribution of duplicated homeobox genes across eight chromosome-level spider genomes. We found that most duplicated homeobox genes in spiders are consistent with an origin by WGD. We also found two copies of conserved homeobox gene clusters, including the Hox, NK, HRO, Irx, and SINE clusters, in all eight species. Consistently, we observed one copy of each cluster was degenerated in terms of gene content and organization while the other remained more intact. Focussing on the NK cluster, we found evidence for regulatory subfunctionalization between the duplicated NK genes in the spider Parasteatoda tepidariorum compared to their single-copy orthologues in the harvestman Phalangium opilio. Our study provides new insights into the relative contributions of multiple modes of duplication to the homeobox gene repertoire during the evolution of spiders and the function of NK genes.

Integration of fluorescence in situ hybridization and chromosome-length genome assemblies revealed synteny map for guinea pig, naked mole-rat, and human

Descriptions of karyotypes of many animal species are currently available. In addition, there has been a significant increase in the number of sequenced genomes and an ever-improving quality of genome assembly. To close the gap between genomic and cytogenetic data we applied fluorescent in situ hybridization (FISH) and Hi-C technology to make the first full chromosome-level genome comparison of the guinea pig (Cavia porcellus), naked mole-rat (Heterocephalus glaber), and human. Comparative chromosome maps obtained by FISH with chromosome-specific probes link genomic scaffolds to individual chromosomes and orient them relative to centromeres and heterochromatic blocks. Hi-C assembly made it possible to close all gaps on the comparative maps and to reveal additional rearrangements that distinguish the karyotypes of the three species. As a result, we integrated the bioinformatic and cytogenetic data and adjusted the previous comparative maps and genome assemblies of the guinea pig, naked mole-rat, and human. Syntenic associations in the two hystricomorphs indicate features of their putative ancestral karyotype. We postulate that the two approaches applied in this study complement one another and provide complete information about the organization of these genomes at the chromosome level.

Spiralian genomics and the evolution of animal genome architecture

Recent developments in sequencing technologies have greatly improved our knowledge of phylogenetic relationships and genomic architectures throughout the tree of life. Spiralia, a diverse clade within Protostomia, is essential for understanding the evolutionary history of parasitism, gene conversion, nervous systems and animal body plans. In this review, we focus on the current hypotheses of spiralian phylogeny and investigate the impact of long-read sequencing on the quality of genome assemblies. We examine chromosome-level assemblies to highlight key genomic features that have driven spiralian evolution, including karyotype, synteny and the Hox gene organization. In addition, we show how chromosome rearrangement has influenced spiralian genomic structures. Although spiralian genomes have undergone substantial changes, they exhibit both conserved and lineage-specific features. We recommend increasing sequencing efforts and expanding functional genomics research to deepen insights into spiralian biology.

Evolutionary conservation of embryonic DNA methylome remodelling in distantly related teleost species

Methylation of cytosines in the CG context (mCG) is the most abundant DNA modification in vertebrates that plays crucial roles in cellular differentiation and identity. After fertilization, DNA methylation patterns inherited from parental gametes are remodelled into a state compatible with embryogenesis. In mammals, this is achieved through the global erasure and re-establishment of DNA methylation patterns. However, in non-mammalian vertebrates like zebrafish, no global erasure has been observed. To investigate the evolutionary conservation and divergence of DNA methylation remodelling in teleosts, we generated base resolution DNA methylome datasets of developing medaka and medaka-zebrafish hybrid embryos. In contrast to previous reports, we show that medaka display comparable DNA methylome dynamics to zebrafish with high gametic mCG levels (sperm: ∼90%; egg: ∼75%), and adoption of a paternal-like methylome during early embryogenesis, with no signs of prior DNA methylation erasure. We also demonstrate that non-canonical DNA methylation (mCH) reprogramming at TGCT tandem repeats is a conserved feature of teleost embryogenesis. Lastly, we find remarkable evolutionary conservation of DNA methylation remodelling patterns in medaka-zebrafish hybrids, indicative of compatible DNA methylation maintenance machinery in far-related teleost species. Overall, these results suggest strong evolutionary conservation of DNA methylation remodelling pathways in teleosts, which is distinct from the global DNA methylome erasure and reestablishment observed in mammals.

Amphioxus as a model to study the evolution of development in chordates

Cephalochordates and tunicates represent the only two groups of invertebrate chordates, and extant cephalochordates - commonly known as amphioxus or lancelets - are considered the best proxy for the chordate ancestor, from which they split around 520 million years ago. Amphioxus has been an important organism in the fields of zoology and embryology since the 18th century, and the morphological and genomic simplicity of cephalochordates (compared to vertebrates) makes amphioxus an attractive model for studying chordate biology at the cellular and molecular levels. Here we describe the life cycle of amphioxus, and discuss the natural histories and habitats of the different species of amphioxus. We also describe their use as laboratory animal models, and discuss the techniques that have been developed to study different aspects of amphioxus.

Linking Vertebrate Gene Duplications to the New Head Hypothesis

Vertebrates have diverse morphologies and various anatomical novelties that set them apart from their closest invertebrate relatives. A conspicuous head housing a large brain, paired sense organs, and protected by a skeleton of cartilage and bone is unique to vertebrates and is a defining feature of this taxon. Gans and Northcutt (1980s) proposed that the evolution of this "new head" was dependent on two key developmental innovations: neural crest cells (NCCs) and ectodermal placodes. NCCs are migratory embryonic cells that form bone, cartilage, and neurons in the new head. Based on genome size, Ohno (1970s) proposed a separate hypothesis, stating that vertebrate genome content was quadrupled via two rounds (2R) of whole genome duplications (WGDs), and the surplus of genetic material potentiated vertebrate morphological diversification. While both hypotheses offer explanations for vertebrate success, it is unclear if, and how, the "new head" and "2R" hypotheses are linked. Here, we consider both hypotheses and evaluate the experimental evidence connecting the two. Overall, evidence suggests that while the origin of the NC GRN predates the vertebrate WGDs, these genomic events may have potentiated the evolution of distinct genetic subnetworks in different neural crest subpopulations. We describe the general composition of the NC GRN and posit that its increased developmental modularity facilitated the independent evolution of NC derivatives and the diversification of the vertebrate head skeleton. Lastly, we discuss experimental strategies needed to test whether gene duplications drove the diversification of neural crest derivatives and the "new head".

Chromatin accessibility, not 5mC methylation covaries with partial dosage compensation in crows

The evolution of genetic sex determination is often accompanied by degradation of the sex-limited chromosome. Male heterogametic systems have evolved convergent, epigenetic mechanisms restoring the resulting imbalance in gene dosage between diploid autosomes (AA) and the hemizygous sex chromosome (X). Female heterogametic systems (AAf Zf, AAm ZZm) tend to only show partial dosage compensation (0.5 < Zf:AAf < 1) and dosage balance (0.5 Collapse

Key Words Collapse

MESH Headings

• Animals
• Female
• Male
• Chromatin/genetics
• Crows/genetics
• Epigenesis, Genetic
• Methylation
• Dosage Compensation, Genetic
• Sex Chromosomes

Collapse

Grants

• HORIZON EUROPE European Research Council
• Swedish Research Council
• Knut och Alice Wallenbergs Stiftelse

Collapse

Affiliation(s)

Ana Catalán Department of Evolutionary Biology, Evolutionary Biology Centre (EBC), Uppsala University, Uppsala, Sweden Division of Evolutionary Biology, LMU Munich, Planegg-Martinsried, Germany
Justin Merondun Division of Evolutionary Biology, LMU Munich, Planegg-Martinsried, Germany
Ulrich Knief Division of Evolutionary Biology, LMU Munich, Planegg-Martinsried, Germany Evolutionary Biology & Ecology,Faculty of Biology, University of Freiburg, Freiburg, Germany
Jochen B. W. Wolf Department of Evolutionary Biology, Evolutionary Biology Centre (EBC), Uppsala University, Uppsala, Sweden Division of Evolutionary Biology, LMU Munich, Planegg-Martinsried, Germany

Collapse

Transposon waves at the water-to-land transition

The major transitions in vertebrate evolution are associated with significant genomic reorganizations. In contrast to the evolutionary processes that occurred at the origin of vertebrates or prior to the radiation of teleost fishes, no whole-genome duplication events occurred during the water-to-land transition, and it remains an open question how did genome dynamics contribute to this prominent evolutionary event. Indeed, the recent sequencing of sarcopterygian and amphibian genomes has revealed that the extant lineages immediately preceding and succeeding this transition harbor an exceptional number of transposable elements and it is tempting to speculate that these sequences might have catalyzed the adaptations that enabled vertebrates to venture into land. Here, we review the genome dynamics associated with the major transitions in vertebrate evolution and discuss how the highly repetitive genomic landscapes revealed by recent efforts to characterize the genomes of amphibians and sarcopterygians argue for turbulent genome dynamics occurring before the water-to-land transition and possibly enabling it.

Sea lamprey enlightens the origin of the coupling of retinoic acid signaling to vertebrate hindbrain segmentation

Retinoic acid (RA) is involved in antero-posterior patterning of the chordate body axis and, in jawed vertebrates, has been shown to play a major role at multiple levels of the gene regulatory network (GRN) regulating hindbrain segmentation. Knowing when and how RA became coupled to the core hindbrain GRN is important for understanding how ancient signaling pathways and patterning genes can evolve and generate diversity. Hence, we investigated the link between RA signaling and hindbrain segmentation in the sea lamprey Petromyzon marinus, an important jawless vertebrate model providing clues to decipher ancestral vertebrate features. Combining genomics, gene expression, and functional analyses of major components involved in RA synthesis (Aldh1as) and degradation (Cyp26s), we demonstrate that RA signaling is coupled to hindbrain segmentation in lamprey. Thus, the link between RA signaling and hindbrain segmentation is a pan vertebrate feature of the hindbrain and likely evolved at the base of vertebrates.

Identification and Characterization of the Very-Low-Density Lipoprotein Receptor Gene from Branchiostoma belcheri: Insights into the Origin and Evolution of the Low-Density Lipoprotein Receptor Gene Family

Low-density lipoprotein receptors (LDLRs) are a class of cell-surface endocytosis receptors that are mainly involved in cholesterol homeostasis and cellular signal transduction. Very-low-density lipoprotein receptors (VLDLRs), which are members of the LDLR family, have been regarded as multi-function receptors that fulfill diverse physiological functions. However, no VLDLR gene has been identified in protochordates to date. As a representative protochordate species, amphioxi are the best available example of vertebrate ancestors. Identifying and characterizing the VLDLR gene in amphioxi has high importance for exploring the evolutionary process of the LDLR family. With this study, a new amphioxus VLDLR gene (designated AmphiVLDLR) was cloned and characterized using RACE-PCR. The 3217 nt transcript of the AmphiVLDLR had a 2577 nt ORF, and the deduced 858 amino acids were highly conserved within vertebrate VLDLRs according to their primary structure and three-dimensional structure, both of which contained five characteristic domains. In contrast to other vertebrate VLDLRs, which had a conserved genomic structure organization with 19 exons and 18 introns, the AmphiVLDLR had 13 exons and 12 introns. The results of a selective pressure analysis showed that the AmphiVLDLR had numerous positive selection sites. Furthermore, the tissue expression of AmphiVLDLR using RT-qPCR showed that AmphiVLDLR RNA expression levels were highest in the gills and muscles, moderate in the hepatic cecum and gonads, and lowest in the intestines. The results of the evolutionary analysis demonstrated that the AmphiVLDLR gene is a new member of the VLDLR family whose family members have experienced duplications and deletions over the evolutionary process. These results imply that the functions of LDLR family members have also undergone differentiation. In summary, we found a new VLDLR gene homolog (AmphiVLDLR) in amphioxi. Our results provide insight into the function and evolution of the LDLR gene family.

Genomic and Transcriptional Profiles of Kelch-like (klhl) Gene Family in Polyploid Carassius Complex

Genome duplication supplies raw genetic materials and has been thought to be essential for evolutionary innovation and ecological adaptation. Here, we select Kelch-like (klhl) genes to study the evolution of the duplicated genes in the polyploid Carassius complex, including amphidiploid C. auratus and amphitriploid C. gibelio. Phylogenetic, chromosomal location and read coverage analyses indicate that most of Carassius klhl genes exhibit a 2:1 relationship with zebrafish orthologs and confirm two rounds of polyploidy, an allotetraploidy followed by an autotriploidy, occurred during Carassius evolution. The lineage-specific expansion and biased retention/loss of klhl genes are also found in Carassius. Transcriptome analyses across eight adult tissues and seven embryogenesis stages reveal varied expression dominance and divergence between the two species. The expression of klhls in response to Carassius herpesvirus 2 infection shows different expression changes corresponding to distinct herpesvirus resistances in three C. gibelio gynogenetic clones. Finally, we find that most C. gibelio klhl genes possess three alleles except eight genes that have lost one or two alleles due to genome rearrangement. The allele expression bias is prosperous for Cgklhl genes and varies during embryogenesis owning to the sequential expression manner of the alleles. The current study provides global insights into the genomic and transcriptional evolution of duplicated genes in a given superfamily resulting from multiple rounds of polyploidization.

The hagfish genome and the evolution of vertebrates

As the only surviving lineages of jawless fishes, hagfishes and lampreys provide a critical window into early vertebrate evolution. Here, we investigate the complex history, timing, and functional role of genome-wide duplications in vertebrates in the light of a chromosome-scale genome of the brown hagfish Eptatretus atami. Using robust chromosome-scale (paralogon-based) phylogenetic methods, we confirm the monophyly of cyclostomes, document an auto-tetraploidization (1RV) that predated the origin of crown group vertebrates ~517 Mya, and establish the timing of subsequent independent duplications in the gnathostome and cyclostome lineages. Some 1RV gene duplications can be linked to key vertebrate innovations, suggesting that this early genomewide event contributed to the emergence of pan-vertebrate features such as neural crest. The hagfish karyotype is derived by numerous fusions relative to the ancestral cyclostome arrangement preserved by lampreys. These genomic changes were accompanied by the loss of genes essential for organ systems (eyes, osteoclast) that are absent in hagfish, accounting in part for the simplification of the hagfish body plan; other gene family expansions account for hagfishes' capacity to produce slime. Finally, we characterise programmed DNA elimination in somatic cells of hagfish, identifying protein-coding and repetitive elements that are deleted during development. As in lampreys, the elimination of these genes provides a mechanism for resolving genetic conflict between soma and germline by repressing germline/pluripotency functions. Reconstruction of the early genomic history of vertebrates provides a framework for further exploration of vertebrate novelties.

Bacterial origin of a key innovation in the evolution of the vertebrate eye

The vertebrate eye was described by Charles Darwin as one of the greatest potential challenges to a theory of natural selection by stepwise evolutionary processes. While numerous evolutionary transitions that led to the vertebrate eye have been explained, some aspects appear to be vertebrate specific with no obvious metazoan precursor. One critical difference between vertebrate and invertebrate vision hinges on interphotoreceptor retinoid-binding protein (IRBP, also known as retinol-binding protein, RBP3), which enables the physical separation and specialization of cells in the vertebrate visual cycle by promoting retinoid shuttling between cell types. While IRBP has been functionally described, its evolutionary origin has remained elusive. Here, we show that IRBP arose via acquisition of novel genetic material from bacteria by interdomain horizontal gene transfer (iHGT). We demonstrate that a gene encoding a bacterial peptidase was acquired prior to the radiation of extant vertebrates >500 Mya and underwent subsequent domain duplication and neofunctionalization to give rise to vertebrate IRBP. Our phylogenomic analyses on >900 high-quality genomes across the tree of life provided the resolution to distinguish contamination in genome assemblies from true instances of horizontal acquisition of IRBP and led us to discover additional independent transfers of the same bacterial peptidase gene family into distinct eukaryotic lineages. Importantly, this work illustrates the evolutionary basis of a key transition that led to the vertebrate visual cycle and highlights the striking impact that acquisition of bacterial genes has had on vertebrate evolution.

Analysis of the P. lividus sea urchin genome highlights contrasting trends of genomic and regulatory evolution in deuterostomes

Sea urchins are emblematic models in developmental biology and display several characteristics that set them apart from other deuterostomes. To uncover the genomic cues that may underlie these specificities, we generated a chromosome-scale genome assembly for the sea urchin Paracentrotus lividus and an extensive gene expression and epigenetic profiles of its embryonic development. We found that, unlike vertebrates, sea urchins retained ancestral chromosomal linkages but underwent very fast intrachromosomal gene order mixing. We identified a burst of gene duplication in the echinoid lineage and showed that some of these expanded genes have been recruited in novel structures (water vascular system, Aristotle's lantern, and skeletogenic micromere lineage). Finally, we identified gene-regulatory modules conserved between sea urchins and chordates. Our results suggest that gene-regulatory networks controlling development can be conserved despite extensive gene order rearrangement.

The little skate genome and the evolutionary emergence of wing-like fins

Skates are cartilaginous fish whose body plan features enlarged wing-like pectoral fins, enabling them to thrive in benthic environments1,2. However, the molecular underpinnings of this unique trait remain unclear. Here we investigate the origin of this phenotypic innovation by developing the little skate Leucoraja erinacea as a genomically enabled model. Analysis of a high-quality chromosome-scale genome sequence for the little skate shows that it preserves many ancestral jawed vertebrate features compared with other sequenced genomes, including numerous ancient microchromosomes. Combining genome comparisons with extensive regulatory datasets in developing fins-including gene expression, chromatin occupancy and three-dimensional conformation-we find skate-specific genomic rearrangements that alter the three-dimensional regulatory landscape of genes that are involved in the planar cell polarity pathway. Functional inhibition of planar cell polarity signalling resulted in a reduction in anterior fin size, confirming that this pathway is a major contributor to batoid fin morphology. We also identified a fin-specific enhancer that interacts with several hoxa genes, consistent with the redeployment of hox gene expression in anterior pectoral fins, and confirmed its potential to activate transcription in the anterior fin using zebrafish reporter assays. Our findings underscore the central role of genome reorganization and regulatory variation in the evolution of phenotypes, shedding light on the molecular origin of an enigmatic trait.

Three amphioxus reference genomes reveal gene and chromosome evolution of chordates

The slow-evolving invertebrate amphioxus has an irreplaceable role in advancing our understanding of the vertebrate origin and innovations. Here we resolve the nearly complete chromosomal genomes of three amphioxus species, one of which best recapitulates the 17 chordate ancestor linkage groups. We reconstruct the fusions, retention, or rearrangements between descendants of whole-genome duplications, which gave rise to the extant microchromosomes likely existed in the vertebrate ancestor. Similar to vertebrates, the amphioxus genome gradually establishes its three-dimensional chromatin architecture at the onset of zygotic activation and forms two topologically associated domains at the Hox gene cluster. We find that all three amphioxus species have ZW sex chromosomes with little sequence differentiation, and their putative sex-determining regions are nonhomologous to each other. Our results illuminate the unappreciated interspecific diversity and developmental dynamics of amphioxus genomes and provide high-quality references for understanding the mechanisms of chordate functional genome evolution.

Annelid functional genomics reveal the origins of bilaterian life cycles

Indirect development with an intermediate larva exists in all major animal lineages1, which makes larvae central to most scenarios of animal evolution2-11. Yet how larvae evolved remains disputed. Here we show that temporal shifts (that is, heterochronies) in trunk formation underpin the diversification of larvae and bilaterian life cycles. We performed chromosome-scale genome sequencing in the annelid Owenia fusiformis with transcriptomic and epigenomic profiling during the life cycles of this and two other annelids. We found that trunk development is deferred to pre-metamorphic stages in the feeding larva of O. fusiformis but starts after gastrulation in the non-feeding larva with gradual metamorphosis of Capitella teleta and the direct developing embryo of Dimorphilus gyrociliatus. Accordingly, the embryos of O. fusiformis develop first into an enlarged anterior domain that forms larval tissues and the adult head12. Notably, this also occurs in the so-called 'head larvae' of other bilaterians13-17, with which the O. fusiformis larva shows extensive transcriptomic similarities. Together, our findings suggest that the temporal decoupling of head and trunk formation, as maximally observed in head larvae, facilitated larval evolution in Bilateria. This diverges from prevailing scenarios that propose either co-option9,10 or innovation11 of gene regulatory programmes to explain larva and adult origins.

Computational approaches to understand transcription regulation in development

Gene regulatory networks (GRNs) serve as useful abstractions to understand transcriptional dynamics in developmental systems. Computational prediction of GRNs has been successfully applied to genome-wide gene expression measurements with the advent of microarrays and RNA-sequencing. However, these inferred networks are inaccurate and mostly based on correlative rather than causative interactions. In this review, we highlight three approaches that significantly impact GRN inference: (1) moving from one genome-wide functional modality, gene expression, to multi-omics, (2) single cell sequencing, to measure cell type-specific signals and predict context-specific GRNs, and (3) neural networks as flexible models. Together, these experimental and computational developments have the potential to significantly impact the quality of inferred GRNs. Ultimately, accurately modeling the regulatory interactions between transcription factors and their target genes will be essential to understand the role of transcription factors in driving developmental gene expression programs and to derive testable hypotheses for validation.

Impacts of microplastics and the associated plastisphere on physiological, biochemical, genetic expression and gut microbiota of the filter-feeder amphioxus

Oceanic plastic pollution is of major concern to marine organisms, especially filter feeders. However, limited is known about the toxic effects of the weathered microplastics instead of the pristine ones. This study evaluates the effects of weathered polystyrene microplastic on a filter-feeder amphioxus under starvation conditions via its exposure to the microplastics previously deployed in the natural seawater allowing for the development of a mature biofilm (so-called plastisphere). The study focused on the integration of physiological, histological, biochemical, molecular, and microbiota impacts on amphioxus. Overall, specific alterations in gene expression of marker genes were observed to be associated with oxidative stresses and immune systems. Negligible impacts were observed on antioxidant biochemical activities and gut microbiota of amphioxus, while we highlighted the potential transfer of 12 bacterial taxa from the plastisphere to the amphioxus gut microbiota. Moreover, the classical perturbation of body shape detected in control animals under starvation conditions (a slim and curved body) but not for amphioxus exposed to microplastic, indicates that the microorganisms colonizing plastics could serve as a nutrient source for this filter-feeder, commitment with the elevated proportions of goblet cell-like structures after the microplastic exposure. The multidisciplinary approach developed in this study underlined the trait of microplastics that acted as vectors for transporting microorganisms from the plastisphere to amphioxus.

Comparative analysis of genome-scale, base-resolution DNA methylation profiles across 580 animal species

Methylation of cytosines is a prototypic epigenetic modification of the DNA. It has been implicated in various regulatory mechanisms across the animal kingdom and particularly in vertebrates. We mapped DNA methylation in 580 animal species (535 vertebrates, 45 invertebrates), resulting in 2443 genome-scale DNA methylation profiles of multiple organs. Bioinformatic analysis of this large dataset quantified the association of DNA methylation with the underlying genomic DNA sequence throughout vertebrate evolution. We observed a broadly conserved link with two major transitions-once in the first vertebrates and again with the emergence of reptiles. Cross-species comparisons focusing on individual organs supported a deeply conserved association of DNA methylation with tissue type, and cross-mapping analysis of DNA methylation at gene promoters revealed evolutionary changes for orthologous genes. In summary, this study establishes a large resource of vertebrate and invertebrate DNA methylomes, it showcases the power of reference-free epigenome analysis in species for which no reference genomes are available, and it contributes an epigenetic perspective to the study of vertebrate evolution.

Active DNA demethylation of developmental cis-regulatory regions predates vertebrate origins

DNA methylation [5-methylcytosine (5mC)] is a repressive gene-regulatory mark required for vertebrate embryogenesis. Genomic 5mC is tightly regulated through the action of DNA methyltransferases, which deposit 5mC, and ten-eleven translocation (TET) enzymes, which participate in its active removal through the formation of 5-hydroxymethylcytosine (5hmC). TET enzymes are essential for mammalian gastrulation and activation of vertebrate developmental enhancers; however, to date, a clear picture of 5hmC function, abundance, and genomic distribution in nonvertebrate lineages is lacking. By using base-resolution 5mC and 5hmC quantification during sea urchin and lancelet embryogenesis, we shed light on the roles of nonvertebrate 5hmC and TET enzymes. We find that these invertebrate deuterostomes use TET enzymes for targeted demethylation of regulatory regions associated with developmental genes and show that the complement of identified 5hmC-regulated genes is conserved to vertebrates. This work demonstrates that active 5mC removal from regulatory regions is a common feature of deuterostome embryogenesis suggestive of an unexpected deep conservation of a major gene-regulatory module.

Identification of Nodal-dependent enhancer of amphioxus Chordin sufficient to drive gene expression into the chordate dorsal organizer

The core molecular mechanisms of dorsal organizer formation during gastrulation are highly conserved within the chordate lineage. One of the key characteristics is that Nodal signaling is required for the organizer-specific gene expression. This feature appears to be ancestral, as evidenced by the presence in the most basally divergent chordate amphioxus. To provide a better understanding of the evolution of organizer-specific gene regulation in chordates, we analyzed the cis-regulatory sequence of amphioxus Chordin in the context of the vertebrate embryo. First, we generated stable zebrafish transgenic lines, and by using light-sheet fluorescent microscopy, characterized in detail the expression pattern of GFP driven by the cis-regulatory sequences of amphioxus Chordin. Next, we performed a 5'deletion analysis and identified an enhancer sufficient to drive the expression of the reporter gene into a chordate dorsal organizer. Finally, we found that the identified enhancer element strongly depends on Nodal signaling, which is consistent with the well-established role of this pathway in the regulation of the expression of dorsal organizer-specific genes across chordates. The enhancer identified in our study may represent a suitable simple system to study the interplay of the evolutionarily conserved regulatory mechanisms operating during early chordate development.

Parallel evolution of amphioxus and vertebrate small-scale gene duplications

BACKGROUND

Amphioxus are non-vertebrate chordates characterized by a slow morphological and molecular evolution. They share the basic chordate body-plan and genome organization with vertebrates but lack their 2R whole-genome duplications and their developmental complexity. For these reasons, amphioxus are frequently used as an outgroup to study vertebrate genome evolution and Evo-Devo. Aside from whole-genome duplications, genes continuously duplicate on a smaller scale. Small-scale duplicated genes can be found in both amphioxus and vertebrate genomes, while only the vertebrate genomes have duplicated genes product of their 2R whole-genome duplications. Here, we explore the history of small-scale gene duplications in the amphioxus lineage and compare it to small- and large-scale gene duplication history in vertebrates.

RESULTS

We present a study of the European amphioxus (Branchiostoma lanceolatum) gene duplications thanks to a new, high-quality genome reference. We find that, despite its overall slow molecular evolution, the amphioxus lineage has had a history of small-scale duplications similar to the one observed in vertebrates. We find parallel gene duplication profiles between amphioxus and vertebrates and conserved functional constraints in gene duplication. Moreover, amphioxus gene duplicates show levels of expression and patterns of functional specialization similar to the ones observed in vertebrate duplicated genes. We also find strong conservation of gene synteny between two distant amphioxus species, B. lanceolatum and B. floridae, with two major chromosomal rearrangements.

CONCLUSIONS

In contrast to their slower molecular and morphological evolution, amphioxus' small-scale gene duplication history resembles that of the vertebrate lineage both in quantitative and in functional terms.

A novel nuclear receptor subfamily enlightens the origin of heterodimerization

Background

Nuclear receptors are transcription factors of central importance in human biology and associated diseases. Much of the knowledge related to their major functions, such as ligand and DNA binding or dimerization, derives from functional studies undertaken in classical model animals. It has become evident, however, that a deeper understanding of these molecular functions requires uncovering how these characteristics originated and diversified during evolution, by looking at more species. In particular, the comprehension of how dimerization evolved from ancestral homodimers to a more sophisticated state of heterodimers has been missing, due to a too narrow phylogenetic sampling. Here, we experimentally and phylogenetically define the evolutionary trajectory of nuclear receptor dimerization by analyzing a novel NR7 subgroup, present in various metazoan groups, including cnidarians, annelids, mollusks, sea urchins, and amphioxus, but lost in vertebrates, arthropods, and nematodes.

Results

We focused on NR7 of the cephalochordate amphioxus B. lanceolatum. We present a complementary set of functional, structural, and evolutionary analyses that establish that NR7 lies at a pivotal point in the evolutionary trajectory from homodimerizing to heterodimerizing nuclear receptors. The crystal structure of the NR7 ligand-binding domain suggests that the isolated domain is not capable of dimerizing with the ubiquitous dimerization partner RXR. In contrast, the full-length NR7 dimerizes with RXR in a DNA-dependent manner and acts as a constitutively active receptor. The phylogenetic and sequence analyses position NR7 at a pivotal point, just between the basal class I nuclear receptors that form monomers or homodimers on DNA and the derived class II nuclear receptors that exhibit the classical DNA-independent RXR heterodimers.

Conclusions

Our data suggest that NR7 represents the “missing link” in the transition between class I and class II nuclear receptors and that the DNA independency of heterodimer formation is a feature that was acquired during evolution. Our studies define a novel paradigm of nuclear receptor dimerization that evolved from DNA-dependent to DNA-independent requirements. This new concept emphasizes the importance of DNA in the dimerization of nuclear receptors, such as the glucocorticoid receptor and other members of this pharmacologically important oxosteroid receptor subfamily. Our studies further underline the importance of studying emerging model organisms for supporting cutting-edge research.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01413-0.

Evolution of stickleback spines through independent cis-regulatory changes at HOXDB

Understanding the mechanisms leading to new traits or additional features in organisms is a fundamental goal of evolutionary biology. We show that HOXDB regulatory changes have been used repeatedly in different fish genera to alter the length and number of the prominent dorsal spines used to classify stickleback species. In Gasterosteus aculeatus (typically 'three-spine sticklebacks'), a variant HOXDB allele is genetically linked to shortening an existing spine and adding an additional spine. In Apeltes quadracus (typically 'four-spine sticklebacks'), a variant HOXDB allele is associated with lengthening a spine and adding an additional spine in natural populations. The variant alleles alter the same non-coding enhancer region in the HOXDB locus but do so by diverse mechanisms, including single-nucleotide polymorphisms, deletions and transposable element insertions. The independent regulatory changes are linked to anterior expansion or contraction of HOXDB expression. We propose that associated changes in spine lengths and numbers are partial identity transformations in a repeating skeletal series that forms major defensive structures in fish. Our findings support the long-standing hypothesis that natural Hox gene variation underlies key patterning changes in wild populations and illustrate how different mutational mechanisms affecting the same region may produce opposite gene expression changes with similar phenotypic outcomes.

Drosophila as a Model Organism to Study Basic Mechanisms of Longevity

The spatio-temporal regulation of gene expression determines the fate and function of various cells and tissues and, as a consequence, the correct development and functioning of complex organisms. Certain mechanisms of gene activity regulation provide adequate cell responses to changes in environmental factors. Aside from gene expression disorders that lead to various pathologies, alterations of expression of particular genes were shown to significantly decrease or increase the lifespan in a wide range of organisms from yeast to human. Drosophila fruit fly is an ideal model system to explore mechanisms of longevity and aging due to low cost, easy handling and maintenance, large number of progeny per adult, short life cycle and lifespan, relatively low number of paralogous genes, high evolutionary conservation of epigenetic mechanisms and signalling pathways, and availability of a wide range of tools to modulate gene expression in vivo. Here, we focus on the organization of the evolutionarily conserved signaling pathways whose components significantly influence the aging process and on the interconnections of these pathways with gene expression regulation.

Evolutionary analysis of swimming speed in early vertebrates challenges the 'New Head Hypothesis'

The ecological context of early vertebrate evolution is envisaged as a long-term trend towards increasingly active food acquisition and enhanced locomotory capabilities culminating in the emergence of jawed vertebrates. However, support for this hypothesis has been anecdotal and drawn almost exclusively from the ecology of living taxa, despite knowledge of extinct phylogenetic intermediates that can inform our understanding of this formative episode. Here we analyse the evolution of swimming speed in early vertebrates based on caudal fin morphology using ancestral state reconstruction and evolutionary model fitting. We predict the lowest and highest ancestral swimming speeds in jawed vertebrates and microsquamous jawless vertebrates, respectively, and find complex patterns of swimming speed evolution with no support for a trend towards more active lifestyles in the lineage leading to jawed groups. Our results challenge the hypothesis of an escalation of Palaeozoic marine ecosystems and shed light into the factors that determined the disparate palaeobiogeographic patterns of microsquamous versus macrosquamous armoured Palaeozoic jawless vertebrates. Ultimately, our results offer a new enriched perspective on the ecological context that underpinned the assembly of vertebrate and gnathostome body plans, supporting a more complex scenario characterized by diverse evolutionary locomotory capabilities reflecting their equally diverse ecologies.