Three amphioxus reference genomes reveal gene and chromosome evolution of chordates

The genome sequence of the amphioxus, Branchiostoma lanceolatum (Pallas, 1774)

We present a genome assembly from a specimen of Branchiostoma lanceolatum (Amphioxus; Chordata; Leptocardii; Amphioxiformes; Branchiostomatidae). The assembly contains two haplotypes with total lengths of 468.40 megabases and 465.81 megabases, respectively. Most of haplotype 1 (99.34%) is scaffolded into 19 chromosomal pseudomolecules. Haplotype 2 is a scaffold level assembly. The mitochondrial genome has also been assembled and is 15.14 kilobases in length.

Evolution and regulation of animal sex chromosomes

Animal sex chromosomes typically carry the upstream sex-determining gene that triggers testis or ovary development and, in some species, are regulated by global dosage compensation in response to functional decay of the Y chromosome. Despite the importance of these pathways, they exhibit striking differences across species, raising fundamental questions regarding the mechanisms underlying their evolutionary turnover. Recent studies of non-model organisms, including insects, reptiles and teleosts, have yielded a broad view of the diversity of sex chromosomes that challenges established theories. Moreover, continued studies in model organisms with recently developed technologies have characterized the dynamics of sex determination and dosage compensation in three-dimensional nuclear space and at single-cell resolution. Here, we synthesize recent insights into sex chromosomes from a variety of species to review their evolutionary dynamics with respect to the canonical model, as well as their diverse mechanisms of regulation.

Single-cell profiling of the amphioxus digestive tract reveals conservation of endocrine cells in chordates

Despite their pivotal role, the evolutionary origins of vertebrate digestive systems remain enigmatic. We explored the cellular characteristics of the amphioxus (Branchiostoma floridae) digestive tract, a model for the presumed primitive chordate digestive system, using bulk tissue companioned with single-cell RNA sequencing. Our findings reveal segmentation and a rich diversity of cell clusters, and we highlight the presence of epithelial-like, ciliated cells in the amphioxus midgut and describe three types of endocrine-like cells that secrete insulin-like, glucagon-like, and somatostatin-like peptides. Furthermore, Pdx, Ilp1, Ilp2, and Ilpr knockout amphioxus lines revealed that, in amphioxus, Pdx does not influence Ilp expression. We also unravel similarity between amphioxus Ilp1 and vertebrate insulin-like growth factor 1 (Igf1) in terms of predicted structure, effects on body growth and amino acid metabolism, and interactions with Igf-binding proteins. These findings indicate that the evolutionary alterations involving the regulatory influence of Pdx over insulin gene expression could have been instrumental in the development of the vertebrate digestive system.

Single-cell Transcriptomic Studies Unveil Potential Nodes of the Notochord Gene Regulatory Network

Transcription factors (TFs) are DNA-binding proteins able to modulate the timing, location, and levels of gene expression by binding to regulatory DNA regions. Therefore, the repertoire of TFs present in the genome of a multicellular organism and the expression of variable constellations of TFs in different cellular cohorts determine the distinctive characteristics of developing tissues and organs. The information on tissue-specific assortments of TFs, their cross-regulatory interactions, and the genes/regulatory regions targeted by each TF is summarized in gene regulatory networks (GRNs), which provide genetic blueprints for the specification, development, and differentiation of multicellular structures. In this study, we review recent transcriptomic studies focused on the complement of TFs expressed in the notochord, a distinctive feature of all chordates. We analyzed notochord-specific datasets available from organisms representative of the three chordate subphyla, and highlighted lineage-specific variations in the suite of TFs expressed in their notochord. We framed the resulting findings within a provisional evolutionary scenario, which allows the formulation of hypotheses on the genetic/genomic changes that sculpted the structure and function of the notochord on an evolutionary scale.

The genomes of all lungfish inform on genome expansion and tetrapod evolution

The genomes of living lungfishes can inform on the molecular-developmental basis of the Devonian sarcopterygian fish-tetrapod transition. We de novo sequenced the genomes of the African (Protopterus annectens) and South American lungfishes (Lepidosiren paradoxa). The Lepidosiren genome (about 91 Gb, roughly 30 times the human genome) is the largest animal genome sequenced so far and more than twice the size of the Australian (Neoceratodus forsteri)1 and African2 lungfishes owing to enlarged intergenic regions and introns with high repeat content (about 90%). All lungfish genomes continue to expand as some transposable elements (TEs) are still active today. In particular, Lepidosiren's genome grew extremely fast during the past 100 million years (Myr), adding the equivalent of one human genome every 10 Myr. This massive genome expansion seems to be related to a reduction of PIWI-interacting RNAs and C2H2 zinc-finger and Krüppel-associated box (KRAB)-domain protein genes that suppress TE expansions. Although TE abundance facilitates chromosomal rearrangements, lungfish chromosomes still conservatively reflect the ur-tetrapod karyotype. Neoceratodus' limb-like fins still resemble those of their extinct relatives and remained phenotypically static for about 100 Myr. We show that the secondary loss of limb-like appendages in the Lepidosiren-Protopterus ancestor was probably due to loss of sonic hedgehog limb-specific enhancers.

Annelid Comparative Genomics and the Evolution of Massive Lineage-Specific Genome Rearrangement in Bilaterians

The organization of genomes into chromosomes is critical for processes such as genetic recombination, environmental adaptation, and speciation. All animals with bilateral symmetry inherited a genome structure from their last common ancestor that has been highly conserved in some taxa but seemingly unconstrained in others. However, the evolutionary forces driving these differences and the processes by which they emerge have remained largely uncharacterized. Here, we analyze genome organization across the phylum Annelida using 23 chromosome-level annelid genomes. We find that while many annelid lineages have maintained the conserved bilaterian genome structure, the Clitellata, a group containing leeches and earthworms, possesses completely scrambled genomes. We develop a rearrangement index to quantify the extent of genome structure evolution and show that, compared to the last common ancestor of bilaterians, leeches and earthworms have among the most highly rearranged genomes of any currently sampled species. We further show that bilaterian genomes can be classified into two distinct categories-high and low rearrangement-largely influenced by the presence or absence, respectively, of chromosome fission events. Our findings demonstrate that animal genome structure can be highly variable within a phylum and reveal that genome rearrangement can occur both in a gradual, stepwise fashion, or rapid, all-encompassing changes over short evolutionary timescales.

Characterization of Myxovirus resistance (Mx) gene from Chinese seabass Lateolabrax maculatus: Insights into the evolution and function of Mx genes

Chinese seabass (Lateolabrax maculatus) stands out as one of the most sought-after and economically significant species in aquaculture within China. Diseases of L. maculatus occur frequently due to the degradation of the germplasm, the aggravation of environmental pollution of water, and the reproduction of pathogenic microorganisms, inflicting considerable economic losses on the Chinese seabass industry. The Myxovirus resistance (Mx) gene plays pivotal roles in the antiviral immune response ranging from mammals to fish. However, the function of the Mx gene in L. maculatus is still unknown. Firstly, the origin and evolutionary history of Mx proteins was elucidated in this study. Subsequently, an Mx gene from L. maculatus (designed as LmMxA gene) was identified, and its functions in combating antiviral and antibacterial threats were investigated. Remarkably, our findings suggested that while Mx group genes were present in chordates, DYN group genes were present in everything from single-celled animals to humans. Furthermore, our investigation revealed that the LmMxA mRNA level increased in the kidney, spleen and liver subsequent to Vibrio anguillarum and poly(I:C) challenged. Immunofluorescence analysis indicated that LmMxA is predominantly localization in the nucleus and the cytoplasm. Notably, the expression of MAVS, IFN1 and Mx1 increased when LmMxA was overexpression within the EPC cells. Moreover, through assessment via cytopathic effect (CPE), virus titer, and antibacterial activity, it becomes evident that LmMxA exerts a dual role in bolstering both antiviral and antibacterial immune responses. These compelling findings laid the foundation for further exploring the mechanism of LmMxA in response to innate immunity of L. maculatus.

Evolutionary origin of the chordate nervous system revealed by amphioxus developmental trajectories

The common ancestor of all vertebrates had a highly sophisticated nervous system, but questions remain about the evolution of vertebrate neural cell types. The amphioxus, a chordate that diverged before the origin of vertebrates, can inform vertebrate evolution. Here we develop and analyse a single-cell RNA-sequencing dataset from seven amphioxus embryo stages to understand chordate cell type evolution and to study vertebrate neural cell type origins. We identified many new amphioxus cell types, including homologues to the vertebrate hypothalamus and neurohypophysis, rooting the evolutionary origin of these structures. On the basis of ancestor-descendant reconstruction of cell trajectories of the amphioxus and other species, we inferred expression dynamics of transcription factor genes throughout embryogenesis and identified three ancient developmental routes forming chordate neurons. We characterized cell specification at the mechanistic level and generated mutant lines to examine the function of five key transcription factors involved in neural specification. Our results show three developmental origins for the vertebrate nervous system: an anterior FoxQ2-dependent mechanism that is deeply conserved in invertebrates, a less-conserved route leading to more posterior neurons in the vertebrate spinal cord and a mechanism for specifying neuromesoderm progenitors that is restricted to chordates. The evolution of neuromesoderm progenitors may have led to a dramatic shift in posterior neural and mesodermal cell fate decisions and the body elongation process in a stem chordate.

40 years of the homeobox: mechanisms of Hox spatial-temporal collinearity in vertebrates

Animal body plans are established during embryonic development by the Hox genes. This patterning process relies on the differential expression of Hox genes along the head-to-tail axis. Hox spatial collinearity refers to the relationship between the organization of Hox genes in clusters and the differential Hox expression, whereby the relative order of the Hox genes within a cluster mirrors the spatial sequence of expression in the developing embryo. In vertebrates, the cluster organization is also associated with the timing of Hox activation, which harmonizes Hox expression with the progressive emergence of axial tissues. Thereby, in vertebrates, Hox temporal collinearity is intimately linked to Hox spatial collinearity. Understanding the mechanisms contributing to Hox temporal and spatial collinearity is thus key to the comprehension of vertebrate patterning. Here, we provide an overview of the main discoveries pertaining to the mechanisms of Hox spatial-temporal collinearity.

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.

Gonadal transcriptome sequencing reveals sexual dimorphism in expression profiling of sex-related genes in Asian arowana (Scleropages formosus)

Asia arowana (Scleropages formosus) is an ornamental fish with high economic value, while its sex determination mechanism is still poorly understood. By far, no morphological evidence or molecular marker has been developed for effective distinguishment of genders, which poses a critical challenge to our captive breeding efforts. In this study, we sequenced gonadal transcriptomes of adult Asian arowanas and revealed differential expression profiling of sex-related genes. Based on the comparative transcriptomics analysis of testes (n = 3) and ovaries (n = 3), we identified a total of 8,872 differentially expressed genes (DEGs) and 18,490 differentially expressed transposable elements (TEs) between male and female individuals. Interestingly, the expression of TEs usually has been more significantly testis-biased than related coding genes. As expected, several genes related to females (such as foxl2 and cyp19a1a) are significantly transcribed in the ovary, and some genes related to male gonad development (such as dmrt1, gsdf and amh) are highly expressed in the testis. This sexual dimorphism is valuable for ascertaining the differential expression patterns of sex-related genes and enriching the genetic resources of this economically important species. These valuable genetic materials thereby provide instructive references for gender identification and one-to-one breeding practices so as to expand fish numbers for a rapid elevation of economic value.

Identification, Expression and Evolutional Analysis of Two cyp19-like Genes in Amphioxus

The mechanism of sex determination and differentiation in animals remains a central focus of reproductive and developmental biology research, and the regulation of sex differentiation in amphioxus remains poorly understood. Cytochrome P450 Family 19 Subfamily A member 1 (CYP19A1) is a crucial sex differentiation gene that catalyzes the conversion of androgens into estrogens. In this study, we identified two aromatase-like genes in amphioxus: cyp19-like1 and cyp19-like2. The cyp19-like1 is more primitive and may represent the ancestral form of cyp19 in zebrafish and other vertebrates, while the cyp19-like2 is likely the result of gene duplication within amphioxus. To gain further insights into the expression level of these two aromatase-like, we examined their expression in different tissues and during different stages of gonad development. While the expression level of the two genes differs in tissues, both are highly expressed in the gonad primordium and are primarily localized to microsomal membrane systems. However, as development proceeds, their expression level decreases significantly. This study enhances our understanding of sex differentiation mechanisms in amphioxus and provides valuable insights into the formation and evolution of sex determination mechanisms in vertebrates.

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.

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.

TransIntegrator: capture nearly full protein-coding transcript variants via integrating Illumina and PacBio transcriptomes

Genes have the ability to produce transcript variants that perform specific cellular functions. However, accurately detecting all transcript variants remains a long-standing challenge, especially when working with poorly annotated genomes or without a known genome. To address this issue, we have developed a new computational method, TransIntegrator, which enables transcriptome-wide detection of novel transcript variants. For this, we determined 10 Illumina sequencing transcriptomes and a PacBio full-length transcriptome for consecutive embryo development stages of amphioxus, a species of great evolutionary importance. Based on the transcriptomes, we employed TransIntegrator to create a comprehensive transcript variant library, namely iTranscriptome. The resulting iTrancriptome contained 91 915 distinct transcript variants, with an average of 2.4 variants per gene. This substantially improved current amphioxus genome annotation by expanding the number of genes from 21 954 to 38 777. Further analysis manifested that the gene expansion was largely ascribed to integration of multiple Illumina datasets instead of involving the PacBio data. Moreover, we demonstrated an example application of TransIntegrator, via generating iTrancriptome, in aiding accurate transcriptome assembly, which significantly outperformed other hybrid methods such as IDP-denovo and Trinity. For user convenience, we have deposited the source codes of TransIntegrator on GitHub as well as a conda package in Anaconda. In summary, this study proposes an affordable but efficient method for reliable transcriptomic research in most species.

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.

Shifting landscapes: the role of 3D genomic organizations in gene regulatory strategies

3D genome folding enables the physical storage of chromosomes into the compact volume of a cell's nucleus, allows for the accurate segregation of chromatin to daughter cells, and has been shown to be tightly coupled to the way in which genetic information is converted into transcriptional programs [1-3]. Importantly, this link between chromatin architecture and gene regulation is a selectable feature in which modifications to chromatin organization accompany, or perhaps even drive the establishment of new regulatory strategies with enduring impacts on animal body plan complexity. Here, we discuss the nature of different 3D genome folding systems found across the tree of life, with particular emphasis on metazoans, and the relative influence of these systems on gene regulation. We suggest how the properties of these folding systems have influenced regulatory strategies employed by different lineages and may have catalyzed the partitioning and specialization of genetic programs that enabled multicellularity and organ-grade body plan complexity.

Emergence of the Synucleins

This study establishes the origin and evolutionary history of the synuclein genes. A combination of phylogenetic analyses of the synucleins from twenty-two model species, characterization of local synteny similarities among humans, sharks and lampreys, and statistical comparisons among lamprey and human chromosomes, provides conclusive evidence for the current diversity of synuclein genes arising from the whole-genome duplications (WGDs) that occurred in vertebrates. An ancestral synuclein gene was duplicated in a first WGD, predating the diversification of all living vertebrates. The two resulting genes are still present in agnathan vertebrates. The second WGD, specific to the gnathostome lineage, led to the emergence of the three classical synuclein genes, SNCA, SNCB and SNCG, which are present in all jawed vertebrate lineages. Additional WGDs have added new genes in both agnathans and gnathostomes, while some gene losses have occurred in particular species. The emergence of synucleins through WGDs prevented these genes from experiencing dosage effects, thus avoiding the potential detrimental effects associated with individual duplications of genes that encode proteins prone to aggregation. Additional insights into the structural and functional features of synucleins are gained through the analysis of the highly divergent synuclein proteins present in chondrichthyans and agnathans.