Deeply conserved synteny resolves early events in vertebrate evolution

Comparative genomics illuminates karyotype and sex chromosome evolution of sharks

Chondrichthyes is an important lineage to reconstruct the evolutionary history of vertebrates. Here, we analyzed genome synteny for six chondrichthyan chromosome-level genomes. Our comparative analysis reveals a slow evolutionary rate of chromosomal changes, with infrequent but independent fusions observed in sharks, skates, and chimaeras. The chondrichthyan common ancestor had a proto-vertebrate-like karyotype, including the presence of 18 microchromosome pairs. The X chromosome is a conversed microchromosome shared by all sharks, suggesting a likely common origin of the sex chromosome at least 181 million years ago. We characterized the Y chromosomes of two sharks that are highly differentiated from the X except for a small young evolutionary stratum and a small pseudoautosomal region. We found that shark sex chromosomes lack global dosage compensation but that dosage-sensitive genes are locally compensated. Our study on shark chromosome evolution enhances our understanding of shark sex chromosomes and vertebrate chromosome evolution.

Allopolyploid subgenome identification and implications for evolutionary analysis

Whole-genome duplications (WGDs) are widespread genomic events in eukaryotes that are hypothesized to contribute to the evolutionary success of many lineages, including flowering plants, Saccharomyces yeast, and vertebrates. WGDs generally can be classified into autopolyploids (ploidy increase descended from one species) or allopolyploids (ploidy increase descended from multiple species). Assignment of allopolyploid progenitor species (called subgenomes in the polyploid) is important to understanding the biology and evolution of polyploids, including the asymmetric subgenome evolution following hybridization (biased fractionation). Here, I review the different methodologies used to identify the ancestors of allopolyploid subgenomes, discuss the advantages and disadvantages of these methods, and outline the implications of how these methods affect the subsequent evolutionary analysis of these genomes.

A plurality of morphological characters need not equate with phylogenetic accuracy: A rare genomic change refutes the placement of Solifugae and Pseudoscorpiones in Haplocnemata

Recent advances in higher-level invertebrate phylogeny have leveraged shared features of genomic architecture to resolve contentious nodes across the tree of life. Yet, the interordinal relationships within Chelicerata have remained recalcitrant given competing topologies in recent molecular analyses. As such, relationships between topologically unstable orders remain supported primarily by morphological cladistic analyses. Solifugae, one such unstable chelicerate order, has long been thought to be the sister group of Pseudoscorpiones, forming the clade Haplocnemata, on the basis of eight putative morphological synapomorphies. The discovery, however, of a shared whole genome duplication placing Pseudoscorpiones in Arachnopulmonata provides the opportunity for a simple litmus test evaluating the validity of Haplocnemata. Here, we present the first developmental transcriptome of a solifuge (Titanopuga salinarum) and survey copy numbers of the homeobox genes for evidence of systemic duplication. We find that over 70% of the identified homeobox genes in T. salinarum are retained in a single copy, while representatives of the arachnopulmonates retain orthologs of those genes as two or more copies. Our results refute the placement of Solifugae in Haplocnemata. Subsequent reevaluation of putative interordinal morphological synapomorphies among chelicerates reveals a high incidence of homoplasy, reversals, and inaccurate coding within Haplocnemata and other small clades, as well as Arachnida more broadly, suggesting existing morphological character matrices are insufficient to resolve chelicerate phylogeny.

Neuropeptides and receptors in the cephalochordate: A crucial model for understanding the origin and evolution of vertebrate neuropeptide systems

Genomes and transcriptomes from diverse organisms are providing a wealth of data to explore the evolution and origin of neuropeptides and their receptors in metazoans. While most neuropeptide-receptor systems have been extensively studied in vertebrates, there is still a considerable lack of understanding regarding their functions in invertebrates, an extraordinarily diverse group that account for the majority of animal species on Earth. Cephalochordates, commonly known as amphioxus or lancelets, serve as the evolutionary proxy of the chordate ancestor. Their key evolutionary position, bridging the invertebrate to vertebrate transition, has been explored to uncover the origin, evolution, and function of vertebrate neuropeptide systems. Amphioxus genomes exhibit a high degree of sequence and structural conservation with vertebrates, and sequence and functional homologues of several vertebrate neuropeptide families are present in cephalochordates. This review aims to provide a comprehensively overview of the recent findings on neuropeptides and their receptors in cephalochordates, highlighting their significance as a model for understanding the complex evolution of neuropeptide signaling in vertebrates.

A genome resource for the marine annelid Platynereis dumerilii

The marine annelid Platynereis dumerilii is a model organism used in many research areas including evolution and development, neurobiology, ecology and regeneration. Here we present the genomes of P. dumerilii and of the closely related P. massiliensis and P. megalops, to facilitate comparative genomic approaches and help explore Platynereis biology. We used long-read sequencing technology and chromosomal-conformation capture along with extensive transcriptomic resources to obtain and annotate a draft genome assembly of ~1.47 Gbp for P. dumerilii, of which more than half represent repeat elements. We predict around 29,000 protein-coding genes, with relatively large intron sizes, over 38,000 non-coding genes, and 580 miRNA loci. We further explore the high genetic variation (~3% heterozygosity) within the Platynereis species complex. Gene ontology reveals the most variable loci to be associated with pigmentation, development and immunity. The current work sets the stage for further development of Platynereis genomic resources.

Genome organization: Raison d'être of ancestral linkage groups

The genomes of extant organisms contain conserved blocks of regions that can be traced back to ancient ancestors, yet the evolutionary pressures that maintained such genomic segments remain unclear. New research on a curious organism with two different genomes sheds light on why our genomes are organized as they are.

Chromosome-level genome assemblies of 2 hemichordates provide new insights into deuterostome origin and chromosome evolution

Deuterostomes are a monophyletic group of animals that includes Hemichordata, Echinodermata (together called Ambulacraria), and Chordata. The diversity of deuterostome body plans has made it challenging to reconstruct their ancestral condition and to decipher the genetic changes that drove the diversification of deuterostome lineages. Here, we generate chromosome-level genome assemblies of 2 hemichordate species, Ptychodera flava and Schizocardium californicum, and use comparative genomic approaches to infer the chromosomal architecture of the deuterostome common ancestor and delineate lineage-specific chromosomal modifications. We show that hemichordate chromosomes (1N = 23) exhibit remarkable chromosome-scale macrosynteny when compared to other deuterostomes and can be derived from 24 deuterostome ancestral linkage groups (ALGs). These deuterostome ALGs in turn match previously inferred bilaterian ALGs, consistent with a relatively short transition from the last common bilaterian ancestor to the origin of deuterostomes. Based on this deuterostome ALG complement, we deduced chromosomal rearrangement events that occurred in different lineages. For example, a fusion-with-mixing event produced an Ambulacraria-specific ALG that subsequently split into 2 chromosomes in extant hemichordates, while this homologous ALG further fused with another chromosome in sea urchins. Orthologous genes distributed in these rearranged chromosomes are enriched for functions in various developmental processes. We found that the deeply conserved Hox clusters are located in highly rearranged chromosomes and that maintenance of the clusters are likely due to lower densities of transposable elements within the clusters. We also provide evidence that the deuterostome-specific pharyngeal gene cluster was established via the combination of 3 pre-assembled microsyntenic blocks. We suggest that since chromosomal rearrangement events and formation of new gene clusters may change the regulatory controls of developmental genes, these events may have contributed to the evolution of diverse body plans among deuterostomes.

Pan-evolutionary and regulatory genome architecture delineated by an integrated macro- and microsynteny approach

The forthcoming massive genome data generated by the Earth BioGenome Project will open up a new era of comparative genomics, for which genome synteny analysis provides an important framework. Profiling genome synteny represents an essential step in elucidating genome architecture, regulatory blocks/elements and their evolutionary history. Here we describe PanSyn, ( https://github.com/yhw320/PanSyn ), the most comprehensive and up-to-date genome synteny pipeline, providing step-by-step instructions and application examples to demonstrate its usage. PanSyn inherits both basic and advanced functions from existing popular tools, offering a user-friendly, highly customized approach for genome macrosynteny analysis and integrated pan-evolutionary and regulatory analysis of genome architecture, which are not yet available in public synteny software or tools. The advantages of PanSyn include: (i) advanced microsynteny analysis by functional profiling of microsynteny genes and associated regulatory elements; (ii) comprehensive macrosynteny analysis, including the inference of karyotype evolution from ancestors to extant species; and (iii) functional integration of microsynteny and macrosynteny for pan-evolutionary profiling of genome architecture and regulatory blocks, as well as integration with external functional genomics datasets from three- or four-dimensional genome and ENCODE projects. PanSyn requires basic knowledge of the Linux environment and Perl programming language and the ability to access a computer cluster, especially for large-scale genomic comparisons. Our protocol can be easily implemented by a competent graduate student or postdoc and takes several days to weeks to execute for dozens to hundreds of genomes. PanSyn provides yet the most comprehensive and powerful tool for integrated evolutionary and functional genomics.

The Evolution of Transglutaminases Underlies the Origin and Loss of Cornified Skin Appendages in Vertebrates

Transglutaminases (TGMs) cross-link proteins by introducing covalent bonds between glutamine and lysine residues. These cross-links are essential for epithelial cornification which enables tetrapods to live on land. Here, we investigated which evolutionary adaptations of vertebrates were associated with specific changes in the family of TGM genes. We determined the catalog of TGMs in the main clades of vertebrates, performed a comprehensive phylogenetic analysis of TGMs, and localized the distribution of selected TGMs in tissues. Our data suggest that TGM1 is the phylogenetically oldest epithelial TGM, with orthologs being expressed in the cornified teeth of the lamprey, a basal vertebrate. Gene duplications led to the origin of TGM10 in stem vertebrates, the origin of TGM2 in jawed vertebrates, and an increasing number of epithelium-associated TGM genes in the lineage leading to terrestrial vertebrates. TGM9 is expressed in the epithelial egg tooth, and its evolutionary origin in stem amniotes coincided with the evolution of embryonic development in eggs that are surrounded by a protective shell. Conversely, viviparous mammals have lost both the epithelial egg tooth and TGM9. TGM3 and TGM6 evolved as regulators of cornification in hair follicles and underwent pseudogenization upon the evolutionary loss of hair in cetaceans. Taken together, this study reveals the gain and loss of vertebrate TGM genes in association with the evolution of cornified skin appendages and suggests an important role of TGM9 in the evolution of amniotes.

Plant genome evolution in the genus Eucalyptus is driven by structural rearrangements that promote sequence divergence

Genomes have a highly organized architecture (nonrandom organization of functional and nonfunctional genetic elements within chromosomes) that is essential for many biological functions, particularly gene expression and reproduction. Despite the need to conserve genome architecture, a high level of structural variation has been observed within species. As species separate and diverge, genome architecture also diverges, becoming increasingly poorly conserved as divergence time increases. However, within plant genomes, the processes of genome architecture divergence are not well described. Here we use long-read sequencing and de novo assembly of 33 phylogenetically diverse, wild and naturally evolving Eucalyptus species, covering 1-50 million years of diverging genome evolution to measure genome architectural conservation and describe architectural divergence. The investigation of these genomes revealed that following lineage divergence, genome architecture is highly fragmented by rearrangements. As genomes continue to diverge, the accumulation of mutations and the subsequent divergence beyond recognition of rearrangements become the primary driver of genome divergence. The loss of syntenic regions also contribute to genome divergence but at a slower pace than that of rearrangements. We hypothesize that duplications and translocations are potentially the greatest contributors to Eucalyptus genome divergence.

Acceleration of genome rearrangement in clitellate annelids

Comparisons of multiple metazoan genomes have revealed the existence of ancestral linkage groups (ALGs), genomic scaffolds sharing sets of orthologous genes that have been inherited from ancestral animals for hundreds of millions of years (Simakov et al. 2022; Schultz et al. 2023) These ALGs have persisted across major animal taxa including Cnidaria, Deuterostomia, Ecdysozoa and Spiralia. Notwithstanding this general trend of chromosome-scale conservation, ALGs have been obliterated by extensive genome rearrangements in certain groups, most notably including Clitellata (oligochaetes and leeches), a group of easily overlooked invertebrates that is of tremendous ecological, agricultural and economic importance (Charles 2019; Barrett 2016). To further investigate these rearrangements, we have undertaken a comparison of 12 clitellate genomes (including four newly sequenced species) and 11 outgroup representatives. We show that these rearrangements began at the base of the Clitellata (rather than progressing gradually throughout polychaete annelids), that the inter-chromosomal rearrangements continue in several clitellate lineages and that these events have substantially shaped the evolution of the otherwise highly conserved Hox cluster.

Whole-genome Sequencing Reveals Autooctoploidy in Chinese Sturgeon and Its Evolutionary Trajectories

The order Acipenseriformes, which includes sturgeons and paddlefishes, represents "living fossils" with complex genomes that are good models for understanding whole-genome duplication (WGD) and ploidy evolution in fishes. Here, we sequenced and assembled the first high-quality chromosome-level genome for the complex octoploid Acipenser sinensis (Chinese sturgeon), a critically endangered species that also represents a poorly understood ploidy group in Acipenseriformes. Our results show that A. sinensis is a complex autooctoploid species containing four kinds of octovalents (8n), a hexavalent (6n), two tetravalents (4n), and a divalent (2n). An analysis taking into account delayed rediploidization reveals that the octoploid genome composition of Chinese sturgeon results from two rounds of homologous WGDs, and further provides insights into the timing of its ploidy evolution. This study provides the first octoploid genome resource of Acipenseriformes for understanding ploidy compositions and evolutionary trajectories of polyploid fishes.

Chromosomal-level genome assembly of golden birdwing Troides aeacus (Felder & Felder, 1860)

The golden birdwing Troides aeacus (Lepidoptera, Papilionidae), a significant species in Asia, faces habitat loss due to urbanization and human activities, necessitating its protection. However, the lack of genomic resources hinders our understanding of their biology and diversity, and impedes our conservation efforts based on genetic information or markers. Here, we present the first chromosomal-level genome assembly of T. aeacus using PacBio SMRT and Omni-C scaffolding technologies. The assembled genome (351 Mb) contains 98.94% of the sequences anchored to 30 pseudo-molecules. The genome assembly has high sequence continuity with contig length N50 = 11.67 Mb and L50 = 14, and scaffold length N50 = 12.2 Mb and L50 = 13. A total of 24,946 protein-coding genes were predicted, with high BUSCO score completeness (98.8% and 94.7% of genome and proteome BUSCO, respectively. This genome offers a significant resource for understanding the swallowtail butterfly biology and carrying out its conservation.

An evolutionary timeline of the oxytocin signaling pathway

Oxytocin is a neuropeptide associated with both psychological and somatic processes like parturition and social bonding. Although oxytocin homologs have been identified in many species, the evolutionary timeline of the entire oxytocin signaling gene pathway has yet to be described. Using protein sequence similarity searches, microsynteny, and phylostratigraphy, we assigned the genes supporting the oxytocin pathway to different phylostrata based on when we found they likely arose in evolution. We show that the majority (64%) of genes in the pathway are 'modern'. Most of the modern genes evolved around the emergence of vertebrates or jawed vertebrates (540 - 530 million years ago, 'mya'), including OXTR, OXT and CD38. Of those, 45% were under positive selection at some point during vertebrate evolution. We also found that 18% of the genes in the oxytocin pathway are 'ancient', meaning their emergence dates back to cellular organisms and opisthokonta (3500-1100 mya). The remaining genes (18%) that evolved after ancient and before modern genes were classified as 'medium-aged'. Functional analyses revealed that, in humans, medium-aged oxytocin pathway genes are highly expressed in contractile organs, while modern genes in the oxytocin pathway are primarily expressed in the brain and muscle tissue.

Exploring the patterns of evolution: Core thoughts and focus on the saltational model

The Modern Synthesis, a pillar in biological thought, united Darwin's species origin concepts with Mendel's laws of character heredity, providing a comprehensive understanding of evolution within species. Highlighting phenotypic variation and natural selection, it elucidated the environment's role as a selective force, shaping populations over time. This framework integrated additional mechanisms, including genetic drift, random mutations, and gene flow, predicting their cumulative effects on microevolution and the emergence of new species. Beyond the Modern Synthesis, the Extended Evolutionary Synthesis expands perspectives by recognizing the role of developmental plasticity, non-genetic inheritance, and epigenetics. We suggest that these aspects coexist in the plant evolutionary process; in this context, we focus on the saltational model, emphasizing how saltation events, such as dichotomous saltation, chromosomal mutations, epigenetic phenomena, and polyploidy, contribute to rapid evolutionary changes. The saltational model proposes that certain evolutionary changes, such as the rise of new species, may result suddenly from single macromutations rather than from gradual changes in DNA sequences and allele frequencies within a species over time. These events, observed in domesticated and wild higher plants, provide well-defined mechanistic bases, revealing their profound impact on plant diversity and rapid evolutionary events. Notably, next-generation sequencing exposes the likely crucial role of allopolyploidy and autopolyploidy (saltational events) in generating new plant species, each characterized by distinct chromosomal complements. In conclusion, through this review, we offer a thorough exploration of the ongoing dissertation on the saltational model, elucidating its implications for our understanding of plant evolutionary processes and paving the way for continued research in this intriguing field.

scRNA-seq analysis of cells comprising the amphioxus notochord

Cephalochordates occupy a key phylogenetic position for deciphering the origin and evolution of chordates, since they diverged earlier than urochordates and vertebrates. The notochord is the most prominent feature of chordates. The amphioxus notochord features coin-shaped cells bearing myofibrils. Notochord-derived hedgehog signaling contributes to patterning of the dorsal nerve cord, as in vertebrates. However, properties of constituent notochord cells remain unknown at the single-cell level. We examined these properties using Iso-seq analysis, single-cell RNA-seq analysis, and in situ hybridization (ISH). Gene expression profiles broadly categorize notochordal cells into myofibrillar cells and non-myofibrillar cells. Myofibrillar cells occupy most of the central portion of the notochord, and some cells extend the notochordal horn to both sides of the ventral nerve cord. Some notochord myofibrillar genes are not expressed in myotomes, suggesting an occurrence of myofibrillar genes that are preferentially expressed in notochord. On the other hand, non-myofibrillar cells contain dorsal, lateral, and ventral Müller cells, and all three express both hedgehog and Brachyury. This was confirmed by ISH, although expression of hedgehog in ventral Müller cells was minimal. In addition, dorsal Müller cells express neural transmission-related genes, suggesting an interaction with nerve cord. Lateral Müller cells express hedgehog and other signaling-related genes, suggesting an interaction with myotomes positioned lateral to the notochord. Ventral Müller cells also expressed genes for FGF- and EGF-related signaling, which may be associated with development of endoderm, ventral to the notochord. Lateral Müller cells were intermediate between dorsal/ventral Müller cells. Since vertebrate notochord contributes to patterning and differentiation of ectoderm (nerve cord), mesoderm (somite), and endoderm, this investigation provides evidence that an ancestral or original form of vertebrate notochord is present in extant cephalochordates.

Convergent and divergent evolution of microRNA-mediated regulation in metazoans

The evolution of microRNAs (miRNAs) has been studied extensively to understand their roles in gene regulation and evolutionary processes. This review focuses on how miRNA-mediated regulation has evolved in bilaterian animals, highlighting both convergent and divergent evolution. Since animals and plants display significant differences in miRNA biogenesis and target recognition, the 'independent origin' hypothesis proposes that miRNA pathways in these groups independently evolved from the RNA interference (RNAi) pathway, leading to modern miRNA repertoires through convergent evolution. However, recent evidence raises the alternative possibility that the miRNA pathway might have already existed in the last common ancestor of eukaryotes, and that the differences in miRNA pathway and miRNA repertoires among animal and plant lineages arise from lineage-specific innovations and losses of miRNA pathways, miRNA acquisition, and loss of miRNAs after eukaryotic divergence. The repertoire of miRNAs has considerably expanded during bilaterian evolution, primarily through de novo creation and duplication processes, generating new miRNAs. Although ancient functionally established miRNAs are rarely lost, many newly emerged miRNAs are transient and lineage specific, following a birth-death evolutionary pattern aligning with the 'out-of-the-testis' and 'transcriptional control' hypotheses. Our focus then shifts to the convergent molecular evolution of miRNAs. We summarize how miRNA clustering and seed mimicry contribute to this phenomenon, and we review how miRNAs from different sources converge to degrade maternal messenger RNAs (mRNAs) during animal development. Additionally, we describe how miRNAs evolve across species due to changes in sequence, seed shifting, arm switching, and spatiotemporal expression patterns, which can result in variations in target sites among orthologous miRNAs across distant strains or species. We also provide a summary of the current understanding regarding how the target sites of orthologous miRNAs can vary across strains or distantly related species. Although many paralogous miRNAs retain their seed or mature sequences after duplication, alterations can occur in the seed or mature sequences or expression patterns of paralogous miRNAs, leading to functional diversification. We discuss our current understanding of the functional divergence between duplicated miRNAs, and illustrate how the functional diversification of duplicated miRNAs impacts target site evolution. By investigating these topics, we aim to enhance our current understanding of the functions and evolutionary dynamics of miRNAs. Additionally, we shed light on the existing challenges in miRNA evolutionary studies, particularly the complexity of deciphering the role of miRNA-mediated regulatory network evolution in shaping gene expression divergence and phenotypic differences among species.

Comparative genomics reveals the dynamics of chromosome evolution in Lepidoptera

Chromosomes are a central unit of genome organization. One-tenth of all described species on Earth are butterflies and moths, the Lepidoptera, which generally possess 31 chromosomes. However, some species display dramatic variation in chromosome number. Here we analyse 210 chromosomally complete lepidopteran genomes and show that the chromosomes of extant lepidopterans are derived from 32 ancestral linkage groups, which we term Merian elements. Merian elements have remained largely intact through 250 million years of evolution and diversification. Against this stable background, eight lineages have undergone extensive reorganization either through numerous fissions or a combination of fusion and fission events. Outside these lineages, fusions are rare and fissions are rarer still. Fusions often involve small, repeat-rich Merian elements and the sex-linked element. Our results reveal the constraints on genome architecture in Lepidoptera and provide a deeper understanding of chromosomal rearrangements in eukaryotic genome evolution.

Highly Resolved Genomes of Two Closely Related Lineages of the Rodent Louse Polyplax serrata with Different Host Specificities

Sucking lice of the parvorder Anoplura are permanent ectoparasites with specific lifestyle and highly derived features. Currently, genomic data are only available for a single species, the human louse Pediculus humanus. Here, we present genomes of two distinct lineages, with different host spectra, of a rodent louse Polyplax serrata. Genomes of these ecologically different lineages are closely similar in gene content and display a conserved order of genes, with the exception of a single translocation. Compared with P. humanus, the P. serrata genomes are noticeably larger (139 vs. 111 Mbp) and encode a higher number of genes. Similar to P. humanus, they are reduced in sensory-related categories such as vision and olfaction. Utilizing genome-wide data, we perform phylogenetic reconstruction and evolutionary dating of the P. serrata lineages. Obtained estimates reveal their relatively deep divergence (∼6.5 Mya), comparable with the split between the human and chimpanzee lice P. humanus and Pediculus schaeffi. This supports the view that the P. serrata lineages are likely to represent two cryptic species with different host spectra. Historical demographies show glaciation-related population size (Ne) reduction, but recent restoration of Ne was seen only in the less host-specific lineage. Together with the louse genomes, we analyze genomes of their bacterial symbiont Legionella polyplacis and evaluate their potential complementarity in synthesis of amino acids and B vitamins. We show that both systems, Polyplax/Legionella and Pediculus/Riesia, display almost identical patterns, with symbionts involved in synthesis of B vitamins but not amino acids.

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.

The origin, evolution, and molecular diversity of the chemokine system

Chemokine signalling performs key functions in cell migration via chemoattraction, such as attracting leukocytes to the site of infection during host defence. The system consists of a ligand, the chemokine, usually secreted outside the cell, and a chemokine receptor on the surface of a target cell that recognises the ligand. Several noncanonical components interact with the system. These include a variety of molecules that usually share some degree of sequence similarity with canonical components and, in some cases, are known to bind to canonical components and/or to modulate cell migration. Whereas canonical components have been described in vertebrate lineages, the distribution of the noncanonical components is less clear. Uncertainty over the relationships between canonical and noncanonical components hampers our understanding of the evolution of the system. We used phylogenetic methods, including gene-tree to species-tree reconciliation, to untangle the relationships between canonical and noncanonical components, identify gene duplication events, and clarify the origin of the system. We found that unrelated ligand groups independently evolved chemokine-like functions. We found noncanonical ligands outside vertebrates, such as TAFA "chemokines" found in urochordates. In contrast, all receptor groups are vertebrate-specific and all-except ACKR1-originated from a common ancestor in early vertebrates. Both ligand and receptor copy numbers expanded through gene duplication events at the base of jawed vertebrates, with subsequent waves of innovation occurring in bony fish and mammals.

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.

Hox genes control homocercal caudal fin development and evolution

Ancient bony fishes had heterocercal tails, like modern sharks and sturgeons, with asymmetric caudal fins and a vertebral column extending into an elongated upper lobe. Teleost fishes, in contrast, developed a homocercal tail characterized by two separate equal-sized fin lobes and the body axis not extending into the caudal fin. A similar heterocercal-to-homocercal transition occurs during teleost ontogeny, although the underlying genetic and developmental mechanisms for either transition remain unresolved. Here, we investigated the role of hox13 genes in caudal fin formation as these genes control posterior identity in animals. Analysis of expression profiles of zebrafish hox13 paralogs and phenotypes of CRISPR/Cas9-induced mutants showed that double hoxb13a and hoxc13a mutants fail to form a caudal fin. Furthermore, single mutants display heterocercal-like morphologies not seen since Mesozoic fossil teleosteomorphs. Relaxation of functional constraints after the teleost genome duplication may have allowed hox13 duplicates to neo- or subfunctionalize, ultimately contributing to the evolution of a homocercal tail in teleost fishes.

Gene modelling and annotation for the Hawaiian bobtail squid, Euprymna scolopes

Coleoid cephalopods possess numerous complex, species-specific morphological and behavioural adaptations, e.g., a uniquely structured nervous system that is the largest among the invertebrates. The Hawaiian bobtail squid (Euprymna scolopes) is one of the most established cephalopod species. With its recent publication of the chromosomal-scale genome assembly and regulatory genomic data, it also emerges as a key model for cephalopod gene regulation and evolution. However, the latest genome assembly has been lacking a native gene model set. Our manuscript describes the generation of new long-read transcriptomic data and, made using this combined with a plethora of publicly available transcriptomic and protein sequence data, a new reference annotation for E. scolopes.

Comparative Genome Microsynteny Illuminates the Fast Evolution of Nuclear Mitochondrial Segments (NUMTs) in Mammals

The escape of DNA from mitochondria into the nuclear genome (nuclear mitochondrial DNA, NUMT) is an ongoing process. Although pervasively observed in eukaryotic genomes, their evolutionary trajectories in a mammal-wide context are poorly understood. The main challenge lies in the orthology assignment of NUMTs across species due to their fast evolution and chromosomal rearrangements over the past 200 million years. To address this issue, we systematically investigated the characteristics of NUMT insertions in 45 mammalian genomes and established a novel, synteny-based method to accurately predict orthologous NUMTs and ascertain their evolution across mammals. With a series of comparative analyses across taxa, we revealed that NUMTs may originate from nonrandom regions in mtDNA, are likely found in transposon-rich and intergenic regions, and unlikely code for functional proteins. Using our synteny-based approach, we leveraged 630 pairwise comparisons of genome-wide microsynteny and predicted the NUMT orthology relationships across 36 mammals. With the phylogenetic patterns of NUMT presence-and-absence across taxa, we constructed the ancestral state of NUMTs given the mammal tree using a coalescent method. We found support on the ancestral node of Fereuungulata within Laurasiatheria, whose subordinal relationships are still controversial. This study broadens our knowledge on NUMT insertion and evolution in mammalian genomes and highlights the merit of NUMTs as alternative genetic markers in phylogenetic inference.

Three foxg1 paralogues in lampreys and gnathostomes-brothers or cousins?

Foxg1 is a key regulator of the early development of the vertebrate forebrain and sensory organs. In this study, we describe for the first time three foxg1 paralogues in lamprey, representative of one of two basally diverged lineages of vertebrates-the agnathans. We also first describe three foxg1 genes in sterlet-representative of one of the evolutionarily ancient clades of gnathostomes. According to the analysis of local genomic synteny, three foxg1 genes of agnathans and gnathostomes have a common origin as a result of two rounds of genomic duplications in the early evolution of vertebrates. At the same time, it is difficult to reliably establish pairwise orthology between foxg1 genes of agnathans and gnathostomes based on the analysis of phylogeny and local genomic synteny, as well as our studies of the spatiotemporal expression of foxg1 genes in the river lamprey Lampetra fluviatilis and the sterlet Acipenser ruthenus. Thus, the appearance of three foxg1 paralogues in agnathans and gnathostomes could have occurred either as a result of two rounds of duplication of the vertebrate common ancestor genome (2R hypothesis) or as a result of the first common round followed by subsequent independent polyploidizations in two evolutionary lineages (1R hypothesis).

The genome of Litomosoides sigmodontis illuminates the origins of Y chromosomes in filarial nematodes

Heteromorphic sex chromosomes are usually thought to have originated from a pair of autosomes that acquired a sex-determining locus and subsequently stopped recombining, leading to degeneration of the sex-limited chromosome. The majority of nematode species lack heteromorphic sex chromosomes and determine sex using an X-chromosome counting mechanism, with males being hemizygous for one or more X chromosomes (XX/X0). Some filarial nematode species, including important parasites of humans, have heteromorphic XX/XY karyotypes. It has been assumed that sex is determined by a Y-linked locus in these species. However, karyotypic analyses suggested that filarial Y chromosomes are derived from the unfused homologue of an autosome involved in an X-autosome fusion event. Here, we generated a chromosome-level reference genome for Litomosoides sigmodontis, a filarial nematode with the ancestral filarial karyotype and sex determination mechanism (XX/X0). By mapping the assembled chromosomes to the rhabditid nematode ancestral linkage (or Nigon) elements, we infer that the ancestral filarial X chromosome was the product of a fusion between NigonX (the ancestrally X-linked element) and NigonD (ancestrally autosomal). In the two filarial lineages with XY systems, there have been two independent X-autosome chromosome fusion events involving different autosomal Nigon elements. In both lineages, the region shared by the neo-X and neo-Y chromosomes is within the ancestrally autosomal portion of the X, confirming that the filarial Y chromosomes are derived from the unfused homologue of the autosome. Sex determination in XY filarial nematodes therefore likely continues to operate via the ancestral X-chromosome counting mechanism, rather than via a Y-linked sex-determining locus.

Divergence time estimates for the hypoxia-inducible factor-1 alpha (HIF1α) reveal an ancient emergence of animals in low-oxygen environments

Unveiling the tempo and mode of animal evolution is necessary to understand the links between environmental changes and biological innovation. Although the earliest unambiguous metazoan fossils date to the late Ediacaran period, molecular clock estimates agree that the last common ancestor (LCA) of all extant animals emerged ~850 Ma, in the Tonian period, before the oldest evidence for widespread ocean oxygenation at ~635-560 Ma in the Ediacaran period. Metazoans are aerobic organisms, that is, they are dependent on oxygen to survive. In low-oxygen conditions, most animals have an evolutionarily conserved pathway for maintaining oxygen homeostasis that triggers physiological changes in gene expression via the hypoxia-inducible factor (HIFa). However, here we confirm the absence of the characteristic HIFa protein domain responsible for the oxygen sensing of HIFa in sponges and ctenophores, indicating the LCA of metazoans lacked the functional protein domain as well, and so could have maintained their transcription levels unaltered under the very low-oxygen concentrations of their environments. Using Bayesian relaxed molecular clock dating, we inferred that the ancestral gene lineage responsible for HIFa arose in the Mesoproterozoic Era, ~1273 Ma (Credibility Interval 957-1621 Ma), consistent with the idea that important genetic machinery associated with animals evolved much earlier than the LCA of animals. Our data suggest at least two duplication events in the evolutionary history of HIFa, which generated three vertebrate paralogs, products of the two successive whole-genome duplications that occurred in the vertebrate LCA. Overall, our results support the hypothesis of a pre-Tonian emergence of metazoans under low-oxygen conditions, and an increase in oxygen response elements during animal evolution.

hsd17b1 is a key gene for ovarian differentiation of the Siberian sturgeon

This is the first work using gonads from undifferentiated, genetically-sexed Siberian sturgeon describing expression changes in genes related to steroid synthesis and female and male sex differentiation. One factor identified as relevant for ovarian differentiation was the gene coding for the enzyme Hsd17b1, which converts estrone into estradiol-17β. hsd17b1 was highly activated in female gonads at 2.5 months of age, around the onset of sex differentiation, preceding activation of two other genes involved in estrogen production (cyp19a1 and foxl2). hsd17b1 was also strongly repressed in males. Two known foxl2 paralogs are found in Siberian sturgeon-foxl2 and foxl2l-but only foxl2 appeared to be associated with ovarian differentiation. With regard to the male pathway, neither 11-oxygenated androgens nor classic male genes (amh, dmrt1, sox9, and dhh) were found to be involved in male sex differentiation, leaving open the question of which genes participate in early male gonad development in this ancient fish. Taken together, these results indicate an estrogen-dependence of female sex differentiation and 11-oxygenated androgen-independence of male sex differentiation.

Maternal dominance contributes to subgenome differentiation in allopolyploid fishes

Teleost fishes, which are the largest and most diverse group of living vertebrates, have a rich history of ancient and recent polyploidy. Previous studies of allotetraploid common carp and goldfish (cyprinids) reported a dominant subgenome, which is more expressed and exhibits biased gene retention. However, the underlying mechanisms contributing to observed 'subgenome dominance' remains poorly understood. Here we report high-quality genomes of twenty-one cyprinids to investigate the origin and subsequent subgenome evolution patterns following three independent allopolyploidy events. We identify the closest extant relatives of the diploid progenitor species, investigate genetic and epigenetic differences among subgenomes, and conclude that observed subgenome dominance patterns are likely due to a combination of maternal dominance and transposable element densities in each polyploid. These findings provide an important foundation to understanding subgenome dominance patterns observed in teleost fishes, and ultimately the role of polyploidy in contributing to evolutionary innovations.

Evolution of lysine-specific demethylase 1 and REST corepressor gene families and their molecular interaction

Lysine-specific demethylase 1A (LSD1) binds to the REST corepressor (RCOR) protein family of corepressors to erase transcriptionally active marks on histones. Functional diversity in these complexes depends on the type of RCOR included, which modulates the catalytic activity of the complex. Here, we studied the duplicative history of the RCOR and LSD gene families and analyzed the evolution of their interaction. We found that RCOR genes are the product of the two rounds of whole-genome duplications that occurred early in vertebrate evolution. In contrast, the origin of the LSD genes traces back before to the divergence of animals and plants. Using bioinformatics tools, we show that the RCOR and LSD1 interaction precedes the RCOR repertoire expansion that occurred in the last common ancestor of jawed vertebrates. Overall, we trace LSD1-RCOR complex evolution and propose that animal non-model species offer advantages in addressing questions about the molecular biology of this epigenetic complex.

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 expression and regulation of the nuclear hormone receptor ERR gene family in the chordate lineage

The Estrogen Related Receptor (ERR) nuclear hormone receptor genes have a wide diversity of roles in vertebrate development. In embryos, ERR genes are expressed in several tissues, including the central and peripheral nervous systems. Here we seek to establish the evolutionary history of chordate ERR genes, their expression and their regulation. We examine ERR expression in mollusc, amphioxus and sea squirt embryos, finding the single ERR orthologue is expressed in the nervous system in all three, with muscle expression also found in the two chordates. We show that most jawed vertebrates and lampreys have four ERR paralogues, and that vertebrate ERR genes were ancestrally linked to Estrogen Receptor genes. One of the lamprey paralogues shares conserved expression domains with jawed vertebrate ERRγ in the embryonic vestibuloacoustic ganglion, eye, brain and spinal cord. Hypothesising that conserved expression derives from conserved regulation, we identify a suite of pan-vertebrate conserved non-coding sequences in ERR introns. We use transgenesis in lamprey and chicken embryos to show that these sequences are regulatory and drive reporter gene expression in the nervous system. Our data suggest an ancient association between ERR and the nervous system, including expression in cells associated with photosensation and mechanosensation. This includes the origin in the vertebrate common ancestor of a suite of regulatory elements in the 3' introns that drove nervous system expression and have been conserved from this point onwards.

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.

Widespread Deviant Patterns of Heterozygosity in Whole-Genome Sequencing Due to Autopolyploidy, Repeated Elements, and Duplication

Most population genomic tools rely on accurate single nucleotide polymorphism (SNP) calling and filtering to meet their underlying assumptions. However, genomic complexity, resulting from structural variants, paralogous sequences, and repetitive elements, presents significant challenges in assembling contiguous reference genomes. Consequently, short-read resequencing studies can encounter mismapping issues, leading to SNPs that deviate from Mendelian expected patterns of heterozygosity and allelic ratio. In this study, we employed the ngsParalog software to identify such deviant SNPs in whole-genome sequencing (WGS) data with low (1.5×) to intermediate (4.8×) coverage for four species: Arctic Char (Salvelinus alpinus), Lake Whitefish (Coregonus clupeaformis), Atlantic Salmon (Salmo salar), and the American Eel (Anguilla rostrata). The analyses revealed that deviant SNPs accounted for 22% to 62% of all SNPs in salmonid datasets and approximately 11% in the American Eel dataset. These deviant SNPs were particularly concentrated within repetitive elements and genomic regions that had recently undergone rediploidization in salmonids. Additionally, narrow peaks of elevated coverage were ubiquitous along all four reference genomes, encompassed most deviant SNPs, and could be partially associated with transposons and tandem repeats. Including these deviant SNPs in genomic analyses led to highly distorted site frequency spectra, underestimated pairwise FST values, and overestimated nucleotide diversity. Considering the widespread occurrence of deviant SNPs arising from a variety of sources, their important impact in estimating population parameters, and the availability of effective tools to identify them, we propose that excluding deviant SNPs from WGS datasets is required to improve genomic inferences for a wide range of taxa and sequencing depths.

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.

Emerging questions on the mechanisms and dynamics of 3D genome evolution in spiralians

Information on how 3D genome topology emerged in animal evolution, how stable it is during development, its role in the evolution of phenotypic novelties and how exactly it affects gene expression is highly debated. So far, data to address these questions are lacking with the exception of a few key model species. Several gene regulatory mechanisms have been proposed, including scenarios where genome topology has little to no impact on gene expression, and vice versa. The ancient and diverse clade of spiralians may provide a crucial testing ground for such mechanisms. Sprialians have followed distinct evolutionary trajectories, with some clades experiencing genome expansions and/or large-scale genome rearrangements, and others undergoing genome contraction, substantially impacting their size and organisation. These changes have been associated with many phenotypic innovations in this clade. In this review, we describe how emerging genome topology data, along with functional tools, allow for testing these scenarios and discuss their predicted outcomes.

Symmetric expression of ohnologs encoding conserved antiviral responses in tetraploid common carp suggest absence of subgenome dominance after whole genome duplication

Allopolyploids often experience subgenome dominance, with one subgenome showing higher levels of gene expression and greater gene retention. Here, we address the functionality of both subgenomes of allotetraploid common carp (Cyprinus carpio) by analysing a functional network of interferon-stimulated genes (ISGs) crucial in anti-viral immune defence. As an indicator of subgenome dominance we investigated retainment of a core set of ohnologous ISGs. To facilitate our functional genomic analysis a high quality genome was assembled (WagV4.0). Transcriptome data from an in vitro experiment mimicking a viral infection was used to infer ISG expression. Transcriptome analysis confirmed induction of 88 ISG ohnologs on both subgenomes. In both control and infected states, average expression of ISG ohnologs was comparable between the two subgenomes. Also, the highest expressing and most inducible gene copies of an ohnolog pair could be derived from either subgenome. We found no strong evidence of subgenome dominance for common carp.

Evolution and phylogenetic distribution of endo-α-mannosidase

While glycans underlie many biological processes, such as protein folding, cell adhesion, and cell-cell recognition, deep evolution of glycosylation machinery remains an understudied topic. N-linked glycosylation is a conserved process in which mannosidases are key trimming enzymes. One of them is the glycoprotein endo-α-1,2-mannosidase which participates in the initial trimming of mannose moieties from an N-linked glycan inside the cis-Golgi. It is unique as the only endo-acting mannosidase found in this organelle. Relatively little is known about its origins and evolutionary history; so far it was reported to occur only in vertebrates. In this work, a taxon-rich bioinformatic survey to unravel the evolutionary history of this enzyme, including all major eukaryotic clades and a wide representation of animals, is presented. The endomannosidase was found to be more widely distributed in animals and other eukaryotes. The protein motif changes in context of the canonical animal enzyme were tracked. Additionally, the data show the two canonical vertebrate endomannosidase genes, MANEA and MANEAL, arose at the second round of the two vertebrate genome duplications and one more vertebrate paralog, CMANEAL, is uncovered. Finally, a framework where N-glycosylation co-evolved with complex multicellularity is described. A better understanding of the evolution of core glycosylation pathways is pivotal to understanding biology of eukaryotes in general, and the Golgi apparatus in particular. This systematic analysis of the endomannosidase evolution is one step toward this goal.

The genome of the deep-sea anemone Actinernus sp. contains a mega-array of ANTP-class homeobox genes

Members of the phylum Cnidaria include sea anemones, corals and jellyfish, and have successfully colonized both marine and freshwater habitats throughout the world. The understanding of how cnidarians adapt to extreme environments such as the dark, high-pressure deep-sea habitat has been hindered by the lack of genomic information. Here, we report the first chromosome-level deep-sea cnidarian genome, of the anemone Actinernus sp., which was 1.39 Gbp in length and contained 44 970 gene models including 14 806 tRNA genes and 30 164 protein-coding genes. Analyses of homeobox genes revealed the longest chromosome hosts a mega-array of Hox cluster, HoxL, NK cluster and NKL homeobox genes; until now, such an array has only been hypothesized to have existed in ancient ancestral genomes. In addition to this striking arrangement of homeobox genes, analyses of microRNAs revealed cnidarian-specific complements that are distinctive for nested clades of these animals, presumably reflecting the progressive evolution of the gene regulatory networks in which they are embedded. Also, compared with other sea anemones, circadian rhythm genes were lost in Actinernus sp., which likely reflects adaptation to living in the dark. This high-quality genome of a deep-sea cnidarian thus reveals some of the likely molecular adaptations of this ecologically important group of metazoans to the extreme deep-sea environment. It also deepens our understanding of the evolution of genome content and organization of animals in general and cnidarians in particular, specifically from the viewpoint of key developmental control genes like the homeobox-encoding genes, where we find an array of genes that until now has only been hypothesized to have existed in the ancient ancestor that pre-dated both the cnidarians and bilaterians.

Ohnologs and SSD Paralogs Differ in Genomic and Expression Features Related to Dosage Constraints

Gene duplication is recognized as a critical process in genome evolution; however, many questions about this process remain unanswered. Although gene duplicability has been observed to differ by duplication mechanism and evolutionary rate, there is so far no broad characterization of its determinants. Many features correlate with this difference in duplicability; however, our ability to exploit these observations to advance our understanding of the role of duplication in evolution is hampered by limitations within existing work. In particular, the existence of methodological differences across studies impedes meaningful comparison. Here, we use consistent definitions of duplicability in the human lineage to explore these associations, allow resolution of the impact of confounding factors, and define the overall relevance of individual features. Using a classifier approach and controlling for the confounding effect of duplicate longevity, we find a subset of gene features important in differentiating genes duplicable by small-scale duplication from those duplicable by whole-genome duplication, revealing critical roles for gene dosage and expression costs in duplicability. We further delve into patterns of functional enrichment and find a lack of constraint on duplicate retention in any context for genes duplicable by small-scale duplication.

Molecular Evolution of Aryl Hydrocarbon Receptor Signaling Pathway Genes

The Aryl hydrocarbon receptor is an ancient transcriptional factor originally discovered as a sensor of dioxin. In addition to its function as a receptor of environmental toxicants, it plays an important role in development. Although a significant amount of research has been carried out to understand the AHR signal transduction pathway and its involvement in species' susceptibility to environmental toxicants, none of them to date has comprehensively studied its evolutionary origins. Studying the evolutionary origins of molecules can inform ancestral relationships of genes. The vertebrate genome has been shaped by two rounds of whole-genome duplications (WGD) at the base of vertebrate evolution approximately 600 million years ago, followed by lineage-specific gene losses, which often complicate the assignment of orthology. It is crucial to understand the evolutionary origins of this transcription factor and its partners, to distinguish orthologs from ancient non-orthologous homologs. In this study, we have investigated the evolutionary origins of proteins involved in the AHR pathway. Our results provide evidence of gene loss and duplications, crucial for understanding the functional connectivity of humans and model species. Multiple studies have shown that 2R-ohnologs (genes and proteins that have survived from the 2R-WGD) are enriched in signaling components relevant to developmental disorders and cancer. Our findings provide a link between the AHR pathway's evolutionary trajectory and its potential mechanistic involvement in pathogenesis.

The genome and sex-dependent responses to temperature in the common yellow butterfly, Eurema hecabe

BACKGROUND

Lepidoptera (butterflies and moths) is one of the most geographically widespread insect orders in the world, and its species play important and diverse ecological and applied roles. Climate change is one of the biggest challenges to biodiversity this century, and lepidopterans are vulnerable to climate change. Temperature-dependent gene expression differences are of relevance under the ongoing climate crisis. However, little is known about how climate affects gene expression in lepidopterans and the ecological consequences of this, particularly with respect to genes with biased expression in one of the sexes. The common yellow butterfly, Eurema hecabe (Family Pieridae), is one of the most geographically widespread lepidopterans that can be found in Asia, Africa, and Australia. Nevertheless, what temperature-dependent effects there may be and whether the effects differ between the sexes remain largely unexplored.

RESULTS

Here, we generated high-quality genomic resources for E. hecabe along with transcriptomes from eight developmental stages. Male and female butterflies were subjected to varying temperatures to assess sex-specific gene expression responses through mRNA and microRNA transcriptomics. We find that there are more temperature-dependent sex-biased genes in females than males, including genes that are involved in a range of biologically important functions, highlighting potential ecological impacts of increased temperatures. Further, by considering available butterfly data on sex-biased gene expression in a comparative genomic framework, we find that the pattern of sex-biased gene expression identified in E. hecabe is highly species-specific, rather than conserved across butterfly species, suggesting that sex-biased gene expression responses to climate change are complex in butterflies.

CONCLUSIONS

Our study lays the foundation for further understanding of differential responses to environmental stress in a widespread lepidopteran model and demonstrates the potential complexity of sex-specific responses of lepidopterans to climate change.

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.

The Impact of Whole Genome Duplication on the Evolution of the Arachnids

The proliferation of genomic resources for Chelicerata in the past 10 years has revealed that the evolution of chelicerate genomes is more dynamic than previously thought, with multiple waves of ancient whole genome duplications affecting separate lineages. Such duplication events are fascinating from the perspective of evolutionary history because the burst of new gene copies associated with genome duplications facilitates the acquisition of new gene functions (neofunctionalization), which may in turn lead to morphological novelties and spur net diversification. While neofunctionalization has been invoked in several contexts with respect to the success and diversity of spiders, the overall impact of whole genome duplications on chelicerate evolution and development remains imperfectly understood. The purpose of this review is to examine critically the role of whole genome duplication on the diversification of the extant arachnid orders, as well as assess functional datasets for evidence of subfunctionalization or neofunctionalization in chelicerates. This examination focuses on functional data from two focal model taxa: the spider Parasteatoda tepidariorum, which exhibits evidence for an ancient duplication, and the harvestman Phalangium opilio, which exhibits an unduplicated genome. I show that there is no evidence that taxa with genome duplications are more successful than taxa with unduplicated genomes. I contend that evidence for sub- or neofunctionalization of duplicated developmental patterning genes in spiders is indirect or fragmentary at present, despite the appeal of this postulate for explaining the success of groups like spiders. Available expression data suggest that the condition of duplicated Hox modules may have played a role in promoting body plan disparity in the posterior tagma of some orders, such as spiders and scorpions, but functional data substantiating this postulate are critically missing. Spatiotemporal dynamics of duplicated transcription factors in spiders may represent cases of developmental system drift, rather than neofunctionalization. Developmental system drift may represent an important, but overlooked, null hypothesis for studies of paralogs in chelicerate developmental biology. To distinguish between subfunctionalization, neofunctionalization, and developmental system drift, concomitant establishment of comparative functional datasets from taxa exhibiting the genome duplication, as well as those that lack the paralogy, is sorely needed.

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".

Elasmobranch genome sequencing reveals evolutionary trends of vertebrate karyotype organization

Genomic studies of vertebrate chromosome evolution have long been hindered by the scarcity of chromosome-scale DNA sequences of some key taxa. One of those limiting taxa has been the elasmobranchs (sharks and rays), which harbor species often with numerous chromosomes and enlarged genomes. Here, we report the chromosome-scale genome assembly for the zebra shark Stegostoma tigrinum, an endangered species that has a relatively small genome among sharks (3.71 Gb), as well as for the whale shark Rhincodon typus Our analysis using a male-female comparison identified an X Chromosome, the first genomically characterized shark sex chromosome. The X Chromosome harbors the Hox C cluster whose intact linkage has not been shown for an elasmobranch fish. The sequenced shark genomes show a gradualism of chromosome length with remarkable length-dependent characteristics-shorter chromosomes tend to have higher GC content, gene density, synonymous substitution rate, and simple tandem repeat content as well as smaller gene length and lower interspersed repeat content. We challenge the traditional binary classification of karyotypes as with and without so-called microchromosomes. Even without microchromosomes, the length-dependent characteristics persist widely in nonmammalian vertebrates. Our investigation of elasmobranch karyotypes underpins their unique characteristics and provides clues for understanding how vertebrate karyotypes accommodate intragenomic heterogeneity to realize a complex readout. It also paves the way to dissecting more genomes with variable sizes to be sequenced at high quality.

Advances in the Xenopus immunome: Diversification, expansion, and contraction

Xenopus is a genus of African clawed frogs including two species, X. tropicalis and X. laevis that are extensively used in experimental biology, immunology, and biomedical studies. The availability of fully sequenced and annotated Xenopus genomes is strengthening genome-wide analyses of gene families and transgenesis to model human diseases. However, inaccuracies in genome annotation for genes involved in the immune system (i.e., immunome) hamper immunogenetic studies. Furthermore, advanced genome technologies (e.g., single-cell and RNA-Seq) rely on well-annotated genomes. The annotation problems of Xenopus immunome include a lack of established orthology across taxa, merged gene models, poor representation in gene pages on Xenbase, misannotated genes and missing gene IDs. The Xenopus Research Resource for Immunobiology in collaboration with Xenbase and a group of investigators are working to resolve these issues in the latest versions of genome browsers. In this review, we summarize the current problems of previously misannotated gene families that we have recently resolved. We also highlight the expansion, contraction, and diversification of previously misannotated gene families.