Author Correction: Elephant shark genome provides unique insights into gnathostome evolution.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

The seahorse genome and the evolution of its specialized morphology.

Seahorses have a specialized morphology that includes a toothless tubular mouth, a body covered with bony plates, a male brood pouch, and the absence of caudal and pelvic fins. Here we report the sequencing and de novo assembly of the genome of the tiger tail seahorse, Hippocampus comes. Comparative genomic analysis identifies higher protein and nucleotide evolutionary rates in H. comes compared with other teleost fish genomes. We identified an astacin metalloprotease gene family that has undergone expansion and is highly expressed in the male brood pouch. We also find that the H. comes genome lacks enamel matrix protein-coding proline/glutamine-rich secretory calcium-binding phosphoprotein genes, which might have led to the loss of mineralized teeth. tbx4, a regulator of hindlimb development, is also not found in H. comes genome. Knockout of tbx4 in zebrafish showed a 'pelvic fin-loss' phenotype similar to that of seahorses.

Conservation as well as divergence in Mcidas function underlies the differentiation of multiciliated cells in vertebrates.

Multiciliated cells (MCCs) differentiate hundreds of motile cilia that beat to drive fluid movement over various kinds of epithelia. In Xenopus, mice and human, the coiled-coil containing protein Mcidas (Mci) has been shown to be a key transcriptional regulator of MCC differentiation. We have examined Mci function in the zebrafish, another model organism that is widely used to study ciliary biology. We show that zebrafish mci is expressed specifically in the developing MCCs of the kidney tubules, but surprisingly, not in those of the nasal placodes. Mci proteins lack a DNA binding domain and associate with the cell-cycle transcription factors E2f4/5 for regulating MCC-specific gene expression. We found that while the zebrafish Mci protein can complex with the E2f family members, its sequence as well as the requirement and sufficiency for MCC differentiation has diverged significantly from Mci homologues of the tetrapods. We also provide evidence that compared to Gmnc, another related coiled-coil protein that has recently been shown to regulate MCC development upstream of Mci, the Mci protein originated later within the vertebrate lineage. Based on these data, we argue that in contrast to Gmnc, which has a vital role in the genetic circuitry that drives MCC formation, the requirement of Mci, at least in the zebrafish, is not obligatory.

Elephant shark genome provides unique insights into gnathostome evolution.

The emergence of jawed vertebrates (gnathostomes) from jawless vertebrates was accompanied by major morphological and physiological innovations, such as hinged jaws, paired fins and immunoglobulin-based adaptive immunity. Gnathostomes subsequently diverged into two groups, the cartilaginous fishes and the bony vertebrates. Here we report the whole-genome analysis of a cartilaginous fish, the elephant shark (Callorhinchus milii). We find that the C. milii genome is the slowest evolving of all known vertebrates, including the 'living fossil' coelacanth, and features extensive synteny conservation with tetrapod genomes, making it a good model for comparative analyses of gnathostome genomes. Our functional studies suggest that the lack of genes encoding secreted calcium-binding phosphoproteins in cartilaginous fishes explains the absence of bone in their endoskeleton. Furthermore, the adaptive immune system of cartilaginous fishes is unusual: it lacks the canonical CD4 co-receptor and most transcription factors, cytokines and cytokine receptors related to the CD4 lineage, despite the presence of polymorphic major histocompatibility complex class II molecules. It thus presents a new model for understanding the origin of adaptive immunity.

Lampreys, the jawless vertebrates, contain only two ParaHox gene clusters.

ParaHox genes (Gsx, Pdx, and Cdx) are an ancient family of developmental genes closely related to the Hox genes. They play critical roles in the patterning of brain and gut. The basal chordate, amphioxus, contains a single ParaHox cluster comprising one member of each family, whereas nonteleost jawed vertebrates contain four ParaHox genomic loci with six or seven ParaHox genes. Teleosts, which have experienced an additional whole-genome duplication, contain six ParaHox genomic loci with six ParaHox genes. Jawless vertebrates, represented by lampreys and hagfish, are the most ancient group of vertebrates and are crucial for understanding the origin and evolution of vertebrate gene families. We have previously shown that lampreys contain six Hox gene loci. Here we report that lampreys contain only two ParaHox gene clusters (designated as α- and ß-clusters) bearing five ParaHox genes (Gsxα, Pdxα, Cdxα, Gsxß, and Cdxß). The order and orientation of the three genes in the α-cluster are identical to that of the single cluster in amphioxus. However, the orientation of Gsxß in the ß-cluster is inverted. Interestingly, Gsxß is expressed in the eye, unlike its homologs in jawed vertebrates, which are expressed mainly in the brain. The lamprey Pdxα is expressed in the pancreas similar to jawed vertebrate Pdx genes, indicating that the pancreatic expression of Pdx was acquired before the divergence of jawless and jawed vertebrate lineages. It is likely that the lamprey Pdxα plays a crucial role in pancreas specification and insulin production similar to the Pdx of jawed vertebrates.

Cyclostomes Lack Clustered Protocadherins.

The brain, comprising billions of neurons and intricate neural networks, is arguably the most complex organ in vertebrates. The diversity of individual neurons is fundamental to the neuronal network complexity and the overall function of the vertebrate brain. In jawed vertebrates, clustered protocadherins provide the molecular basis for this neuronal diversity, through stochastic and combinatorial expression of their various isoforms in individual neurons. Based on analyses of transcriptomes from the Japanese lamprey brain and sea lamprey embryos, genome assemblies of the two lampreys, and brain expressed sequence tags of the inshore hagfish, we show that extant jawless vertebrates (cyclostomes) lack the clustered protocadherins. Our findings indicate that the clustered protocadherins originated from a nonclustered protocadherin in the jawed vertebrate ancestor, after the two rounds of whole-genome duplication. In the absence of clustered protocadherins, cyclostomes might have evolved novel molecules or mechanisms for generating neuronal diversity which remains to be discovered.

Identification of three somatostatin genes in lampreys.

Somatostatins (SSs) are a structurally diverse family of neuropeptides that play important roles in the regulation of growth, development and metabolism in vertebrates. It has been recently proposed that the common ancestor of gnathostomes possessed three SS genes, namely SS1, SS2 and SS5. SS1 and SS2 are still present in most extant gnathostome species investigated so far while SS5 primarily occurs in chondrichthyes, actinopterygians and actinistia but not in tetrapods. Very little is known about the repertoire of SSs in cyclostomes, which are extant jawless vertebrates. In the present study, we report the cloning of the cDNAs encoding three distinct lamprey SS variants that we call SSa, SSb and SSc. SSa and SSb correspond to the two SS variants previously characterized in lamprey, while SSc appears to be a totally novel one. SSa exhibits the same sequence as gnathostome SS1. SSb differs from SSa by only one substitution (Thr12→Ser). SSc exhibits a totally unique structure (ANCRMFYWKTMAAC) that shares only 50% identity with SSa and SSb. SSa, SSb and SSc precursors do not exhibit any appreciable sequence similarity outside the C-terminal region containing the SS sequence. Phylogenetic analyses failed to clearly assign orthology relationships between lamprey and gnathostome SS genes. Synteny analysis suggests that the SSc gene arose before the split of the three gnathostome genes SS1, SS2 and SS5.

Sequencing of Pax6 loci from the elephant shark reveals a family of Pax6 genes in vertebrate genomes, forged by ancient duplications and divergences.

Pax6 is a developmental control gene essential for eye development throughout the animal kingdom. In addition, Pax6 plays key roles in other parts of the CNS, olfactory system, and pancreas. In mammals a single Pax6 gene encoding multiple isoforms delivers these pleiotropic functions. Here we provide evidence that the genomes of many other vertebrate species contain multiple Pax6 loci. We sequenced Pax6-containing BACs from the cartilaginous elephant shark (Callorhinchus milii) and found two distinct Pax6 loci. Pax6.1 is highly similar to mammalian Pax6, while Pax6.2 encodes a paired-less Pax6. Using synteny relationships, we identify homologs of this novel paired-less Pax6.2 gene in lizard and in frog, as well as in zebrafish and in other teleosts. In zebrafish two full-length Pax6 duplicates were known previously, originating from the fish-specific genome duplication (FSGD) and expressed in divergent patterns due to paralog-specific loss of cis-elements. We show that teleosts other than zebrafish also maintain duplicate full-length Pax6 loci, but differences in gene and regulatory domain structure suggest that these Pax6 paralogs originate from a more ancient duplication event and are hence renamed as Pax6.3. Sequence comparisons between mammalian and elephant shark Pax6.1 loci highlight the presence of short- and long-range conserved noncoding elements (CNEs). Functional analysis demonstrates the ancient role of long-range enhancers for Pax6 transcription. We show that the paired-less Pax6.2 ortholog in zebrafish is expressed specifically in the developing retina. Transgenic analysis of elephant shark and zebrafish Pax6.2 CNEs with homology to the mouse NRE/Pα internal promoter revealed highly specific retinal expression. Finally, morpholino depletion of zebrafish Pax6.2 resulted in a "small eye" phenotype, supporting a role in retinal development. In summary, our study reveals that the pleiotropic functions of Pax6 in vertebrates are served by a divergent family of Pax6 genes, forged by ancient duplication events and by independent, lineage-specific gene losses.

Evidence for at least six Hox clusters in the Japanese lamprey (Lethenteron japonicum).

Cyclostomes, comprising jawless vertebrates such as lampreys and hagfishes, are the sister group of living jawed vertebrates (gnathostomes) and hence an important group for understanding the origin and diversity of vertebrates. In vertebrates and other metazoans, Hox genes determine cell fate along the anteroposterior axis of embryos and are implicated in driving morphological diversity. Invertebrates contain a single Hox cluster (either intact or fragmented), whereas elephant shark, coelacanth, and tetrapods contain four Hox clusters owing to two rounds of whole-genome duplication ("1R" and "2R") during early vertebrate evolution. By contrast, most teleost fishes contain up to eight Hox clusters because of an additional "teleost-specific" genome duplication event. By sequencing bacterial artificial chromosome (BAC) clones and the whole genome, here we provide evidence for at least six Hox clusters in the Japanese lamprey (Lethenteron japonicum). This suggests that the lamprey lineage has experienced an additional genome duplication after 1R and 2R. The relative age of lamprey and human paralogs supports this hypothesis. Compared with gnathostome Hox clusters, lamprey Hox clusters are unusually large. Several conserved noncoding elements (CNEs) were predicted in the Hox clusters of lamprey, elephant shark, and human. Transgenic zebrafish assay indicated the potential of CNEs to function as enhancers. Interestingly, CNEs in individual lamprey Hox clusters are frequently conserved in multiple Hox clusters in elephant shark and human, implying a many-to-many orthology relationship between lamprey and gnathostome Hox clusters. Such a relationship suggests that the first two rounds of genome duplication may have occurred independently in the lamprey and gnathostome lineages.

A survey of ancient conserved non-coding elements in the PAX6 locus reveals a landscape of interdigitated cis-regulatory archipelagos.

Biological differences between cell types and developmental processes are characterised by differences in gene expression profiles. Gene-distal enhancers are key components of the regulatory networks that specify the tissue-specific expression patterns driving embryonic development and cell fate decisions, and variations in their sequences are a major contributor to genetic disease and disease susceptibility. Despite advances in the methods for discovery of putative cis-regulatory sequences, characterisation of their spatio-temporal enhancer activities in a mammalian model system remains a major bottle-neck. We employed a strategy that combines gnathostome sequence conservation with transgenic mouse and zebrafish reporter assays to survey the genomic locus of the developmental control gene PAX6 for the presence of novel cis-regulatory elements. Sequence comparison between human and the cartilaginous elephant shark (Callorhinchus milii) revealed several ancient gnathostome conserved non-coding elements (agCNEs) dispersed widely throughout the PAX6 locus, extending the range of the known PAX6 cis-regulatory landscape to contain the full upstream PAX6-RCN1 intergenic region. Our data indicates that ancient conserved regulatory sequences can be tested effectively in transgenic zebrafish even when not conserved in zebrafish themselves. The strategy also allows efficient dissection of compound regulatory regions previously assessed in transgenic mice. Remarkable overlap in expression patterns driven by sets of agCNEs indicates that PAX6 resides in a landscape of multiple tissue-specific regulatory archipelagos.

An ancient genomic regulatory block conserved across bilaterians and its dismantling in tetrapods by retrogene replacement.

Developmental genes are regulated by complex, distantly located cis-regulatory modules (CRMs), often forming genomic regulatory blocks (GRBs) that are conserved among vertebrates and among insects. We have investigated GRBs associated with Iroquois homeobox genes in 39 metazoans. Despite 600 million years of independent evolution, Iroquois genes are linked to ankyrin-repeat-containing Sowah genes in nearly all studied bilaterians. We show that Iroquois-specific CRMs populate the Sowah locus, suggesting that regulatory constraints underlie the maintenance of the Iroquois-Sowah syntenic block. Surprisingly, tetrapod Sowah orthologs are intronless and not associated with Iroquois; however, teleost and elephant shark data demonstrate that this is a derived feature, and that many Iroquois-CRMs were ancestrally located within Sowah introns. Retroposition, gene, and genome duplication have allowed selective elimination of Sowah exons from the Iroquois regulatory landscape while keeping associated CRMs, resulting in large associated gene deserts. These results highlight the importance of CRMs in imposing constraints to genome architecture, even across large phylogenetic distances, and of gene duplication-mediated genetic redundancy to disentangle these constraints, increasing genomic plasticity.

Basal vertebrates clarify the evolutionary history of ciliopathy-associated genes Tmem138 and Tmem216.

Recently, Lee et al. (Lee JH, Silhavy JL, Lee JE, et al. (30 co-authors). 2012. Evolutionarily assembled cis-regulatory module at a human ciliopathy locus. Science (335:966-969.) demonstrated that mutation in either of the transmembrane protein encoding genes, TMEM138 or TMEM216, causes phenotypically indistinguishable ciliopathy. Furthermore, on the basis of the observation that their orthologs are linked in a head-to-tail configuration in other mammals and Anolis, but present on different scaffolds or chromosomes in Xenopus tropicalis and zebrafish, the authors concluded that the two genes were joined by chromosomal rearrangement at the evolutionary amphibian-to-reptile transition to form a functional module. We have sequenced these gene loci in a cartilaginous fish, the elephant shark, and found that the two genes together with a related gene (Tmem80) constitute a tandem cluster. This suggests that the two genes were already linked in the vertebrate ancestor and then rearranged independently in Xenopus and zebrafish. Analyses of the coelacanth and lamprey genomes support this hypothesis. Our study highlights the importance of basal vertebrates as critical reference genomes.

Elephant shark (Callorhinchus milii) provides insights into the evolution of Hox gene clusters in gnathostomes.

We have sequenced and analyzed Hox gene clusters from elephant shark, a holocephalian cartilaginous fish. Elephant shark possesses 4 Hox clusters with 45 Hox genes that include orthologs for a higher number of ancient gnathostome Hox genes than the 4 clusters in tetrapods and the supernumerary clusters in teleost fishes. Phylogenetic analysis of elephant shark Hox genes from 7 paralogous groups that contain all of the 4 members indicated an ((AB)(CD)) topology for the order of Hox cluster duplication, providing support for the 2R hypothesis (i.e., 2 rounds of whole-genome duplication during the early evolution of vertebrates). Comparisons of noncoding sequences of the elephant shark and human Hox clusters have identified a large number of conserved noncoding elements (CNEs), which represent putative cis-regulatory elements that may be involved in the regulation of Hox genes. Interestingly, in fugu more than 50% of these ancient CNEs have diverged beyond recognition in the duplicated (HoxA, HoxB, and HoxD) as well as the singleton (HoxC) Hox clusters. Furthermore, the b-paralogs of the duplicated fugu Hox clusters are virtually devoid of unique ancient CNEs. In contrast to fugu Hox clusters, elephant shark and human Hox clusters have lost fewer ancient CNEs. If these ancient CNEs are indeed enhancers directing tissue-specific expression of Hox genes, divergence of their sequences in vertebrate lineages might have led to altered expression patterns and presumably the functions of their associated Hox genes.

Discovery of a genetic module essential for assigning left-right asymmetry in humans and ancestral vertebrates.

The vertebrate left-right axis is specified during embryogenesis by a transient organ: the left-right organizer (LRO). Species including fish, amphibians, rodents and humans deploy motile cilia in the LRO to break bilateral symmetry, while reptiles, birds, even-toed mammals and cetaceans are believed to have LROs without motile cilia. We searched for genes whose loss during vertebrate evolution follows this pattern and identified five genes encoding extracellular proteins, including a putative protease with hitherto unknown functions that we named ciliated left-right organizer metallopeptide (CIROP). Here, we show that CIROP is specifically expressed in ciliated LROs. In zebrafish and Xenopus, CIROP is required solely on the left side, downstream of the leftward flow, but upstream of DAND5, the first asymmetrically expressed gene. We further ascertained 21 human patients with loss-of-function CIROP mutations presenting with recessive situs anomalies. Our findings posit the existence of an ancestral genetic module that has twice disappeared during vertebrate evolution but remains essential for distinguishing left from right in humans.

Whole genome sequence of the rifamycin B-producing strain Amycolatopsis mediterranei S699.

Amycolatopsis mediterranei S699 is an actinomycete that produces an important antibiotic, rifamycin B. Semisynthetic derivatives of rifamycin B are used for the treatment of tuberculosis, leprosy, and AIDS-related mycobacterial infections. Here, we report the complete genome sequence (10.2 Mb) of A. mediterranei S699, with 9,575 predicted coding sequences.

Reconstruction of proto-vertebrate, proto-cyclostome and proto-gnathostome genomes provides new insights into early vertebrate evolution.

Ancient polyploidization events have had a lasting impact on vertebrate genome structure, organization and function. Some key questions regarding the number of ancient polyploidization events and their timing in relation to the cyclostome-gnathostome divergence have remained contentious. Here we generate de novo long-read-based chromosome-scale genome assemblies for the Japanese lamprey and elephant shark. Using these and other representative genomes and developing algorithms for the probabilistic macrosynteny model, we reconstruct high-resolution proto-vertebrate, proto-cyclostome and proto-gnathostome genomes. Our reconstructions resolve key questions regarding the early evolutionary history of vertebrates. First, cyclostomes diverged from the lineage leading to gnathostomes after a shared tetraploidization (1R) but before a gnathostome-specific tetraploidization (2R). Second, the cyclostome lineage experienced an additional hexaploidization. Third, 2R in the gnathostome lineage was an allotetraploidization event, and biased gene loss from one of the subgenomes shaped the gnathostome genome by giving rise to remarkably conserved microchromosomes. Thus, our reconstructions reveal the major evolutionary events and offer new insights into the origin and evolution of vertebrate genomes.

Antigen receptor repertoires of one of the smallest known vertebrates.

The rules underlying the structure of antigen receptor repertoires are not yet fully defined, despite their enormous importance for the understanding of adaptive immunity. With current technology, the large antigen receptor repertoires of mice and humans cannot be comprehensively studied. To circumvent the problems associated with incomplete sampling, we have studied the immunogenetic features of one of the smallest known vertebrates, the cyprinid fish Paedocypris sp. "Singkep" ("minifish"). Despite its small size, minifish has the key genetic facilities characterizing the principal vertebrate lymphocyte lineages. As described for mammals, the frequency distributions of immunoglobulin and T cell receptor clonotypes exhibit the features of fractal systems, demonstrating that self-similarity is a fundamental property of antigen receptor repertoires of vertebrates, irrespective of body size. Hence, minifish achieve immunocompetence via a few thousand lymphocytes organized in robust scale-free networks, thereby ensuring immune reactivity even when cells are lost or clone sizes fluctuate during immune responses.

Seadragon genome analysis provides insights into its phenotype and sex determination locus.

The iconic phenotype of seadragons includes leaf-like appendages, a toothless tubular mouth, and male pregnancy involving incubation of fertilized eggs on an open "brood patch." We de novo-sequenced male and female genomes of the common seadragon (Phyllopteryx taeniolatus) and its closely related species, the alligator pipefish (Syngnathoides biaculeatus). Transcription profiles from an evolutionary novelty, the leaf-like appendages, show that a set of genes typically involved in fin development have been co-opted as well as an enrichment of transcripts for potential tissue repair and immune defense genes. The zebrafish mutants for scpp5, which is lost in all syngnathids, were found to lack or have deformed pharyngeal teeth, supporting the hypothesis that the loss of scpp5 has contributed to the loss of teeth in syngnathids. A putative sex-determining locus encoding a male-specific amhr2y gene shared by common seadragon and alligator pipefish was identified.

Chromosome-level genome assembly of the coastal horseshoe crab (Tachypleus gigas).

Horseshoe crabs, represented by only four extant species, have existed for around 500 million years. However, their existence is now under threat because of anthropogenic activities. The availability of genomic resources for these species will be valuable in planning appropriate conservation measures. Whole-genome sequences are currently available for three species. In this study, we have generated a chromosome-level genome assembly of the fourth species, the Asian coastal horseshoe crab (Tachypleus gigas; genome size 2.0 Gb). The genome assembly has a scaffold N50 value of 140 Mb with ~97% of the assembly mapped to 14 scaffolds representing 14 chromosomes of T. gigas. In addition, we have generated the complete mitochondrial genome sequence and deep-coverage transcriptome assemblies for four tissues. A total of 26,159 protein-coding genes were predicted in the genome. The T. gigas genome contains five Hox clusters similar to the mangrove horseshoe crab (Carcinoscorpius rotundicauda), suggesting that the common ancestor of horseshoe crabs already possessed five Hox clusters. Phylogenomic and divergence time analysis suggested that the American and Asian horseshoe crab lineages shared a common ancestor around the Silurian period (~436 Ma). Comparison of the T. gigas genome with those of other horseshoe crab species with chromosome-level assemblies provided insights into the chromosomal rearrangement events that occurred during the emergence of these species. The genomic resources of T. gigas will be useful for understanding their genetic diversity and population structure and would help in designing strategies for managing and conserving their stocks across Asia.