About 500 million years ago, a new type of adaptive immune defense emerged in basal jawed vertebrates, accompanied by morphological innovations, including the thymus. Did these evolutionary novelties arise de novo or from elaboration of ancient genetic networks? We reconstructed the genetic changes underlying thymopoiesis by comparative genome and expression analyses in chordates and basal vertebrates. The derived models of genetic networks were experimentally verified in bony fishes. Ancestral networks defining circumscribed regions of the pharyngeal epithelium of jawless vertebrates expanded in cartilaginous fishes to incorporate novel genes, notably those encoding chemokines. Correspondingly, novel networks evolved in lymphocytes of jawed vertebrates to control the expression of additional chemokine receptors. These complementary changes enabled unprecedented Delta/Notch signaling between pharyngeal epithelium and lymphoid cells that was exploited for specification to the T cell lineage. Our results provide a framework elucidating the evolution of key features of the adaptive immune system in jawed vertebrates.

Significance Lampreys and hagfishes (cyclostomes) are the only living group of jawless vertebrates and therefore are important for the study of vertebrate evolution. We have characterized Hox clusters in the Japanese lamprey ( Lethenteron japonicum ), and shown that it contains at least six Hox clusters as compared with four Hox clusters in tetrapods. This suggests that the lamprey lineage has undergone an additional round of genome duplication compared with tetrapods. Several conserved noncoding elements (CNEs) were predicted in the Hox clusters of lamprey, elephant shark, and human. Transgenic assay of CNEs demonstrated their potential to function as cis -regulatory elements. Thus, these CNEs may represent part of the core set of cis -regulatory elements that were present in the common ancestor of vertebrates.

With our increasing ability for generating whole-genome sequences, comparative analysis of whole genomes has become a powerful tool for understanding the structure, function, and evolutionary history of human and other vertebrate genomes. By virtue of their position basal to bony vertebrates, cartilaginous fishes (class Chondrichthyes) are a valuable outgroup in comparative studies of vertebrates. Recently, a holocephalan cartilaginous fish, the elephant shark, Callorhinchus milii (Subclass Holocephali: Order Chimaeriformes), has been proposed as a model genome, and low-coverage sequence of its genome has been generated. Despite such an increasing interest, the evolutionary history of the modern holocephalans-a previously successful and diverse group but represented by only 39 extant species-and their relationship with elasmobranchs and other jawed vertebrates has been poorly documented largely owing to a lack of well-preserved fossil materials after the end-Permian about 250 Ma. In this study, we assembled the whole mitogenome sequences for eight representatives from all the three families of the modern holocephalans and investigated their phylogenetic relationships and evolutionary history. Unambiguously aligned sequences from these holocephalans together with 17 other vertebrates (9,409 nt positions excluding entire third codon positions) were subjected to partitioned maximum likelihood analysis. The resulting tree strongly supported a single origin of the modern holocephalans and their sister-group relationship with elasmobranchs. The mitogenomic tree recovered the most basal callorhinchids within the chimaeriforms, which is sister to a clade comprising the remaining two families (rhinochimaerids and chimaerids). The timetree derived from a relaxed molecular clock Bayesian method suggests that the holocephalans originated in the Silurian about 420 Ma, having survived from the end-Permian (250 Ma) mass extinction and undergoing familial diversifications during the late Jurassic to early Cretaceous (170-120 Ma). This postulated evolutionary scenario agrees well with that based on the paleontological observations.

The Asian arowana (Scleropages formosus), one of the world’s most expensive cultivated ornamental fishes, is an endangered species. It represents an ancient lineage of teleosts: the Osteoglossomorpha. Here, we provide a high-quality chromosome-level reference genome of a female golden-variety arowana using a combination of deep shotgun sequencing and high-resolution linkage mapping. In addition, we have also generated two draft genome assemblies for the red and green varieties. Phylogenomic analysis supports a sister group relationship between Osteoglossomorpha (bonytongues) and Elopomorpha (eels and relatives), with the two clades together forming a sister group of Clupeocephala which includes all the remaining teleosts. The arowana genome retains the full complement of eight Hox clusters unlike the African butterfly fish (Pantodon buchholzi), another bonytongue fish, which possess only five Hox clusters. Differential gene expression among three varieties provides insights into the genetic basis of colour variation. A potential heterogametic sex chromosome is identified in the female arowana karyotype, suggesting that the sex is determined by a ZW/ZZ sex chromosomal system. The high-quality reference genome of the golden arowana and the draft assemblies of the red and green varieties are valuable resources for understanding the biology, adaptation and behaviour of Asian arowanas.

The cartilaginous fishes reside at the base of the gnathostome lineage as the oldest extant group of jawed vertebrates. Recently, the genome of the elephant shark, Callorhinchus milii, a chimaerid holocephalan, has been sequenced and therefore becomes the first cartilaginous fish to be analyzed in this way. The chimaeras have been largely neglected and very little is known about the visual systems of these fishes. By searching the elephant shark genome, we have identified gene fragments encoding a rod visual pigment, Rh1, and three cone visual pigments, the middle wavelength-sensitive or Rh2 pigment, and two isoforms of the long wavelength-sensitive or LWS pigment, LWS1 and LWS2, but no evidence for the two short wavelength-sensitive cone classes, SWS1 and SWS2. Expression of these genes in the retina was confirmed by RT-PCR. Full-length coding sequences were used for in vitro expression and gave the following peak absorbances: Rh1 496 nm, Rh2 442 nm, LWS1 499 nm, and LWS2 548 nm. Unusually, therefore, for a deep-sea fish, the elephant shark possesses cone pigments and the potential for trichromacy. Compared with other vertebrates, the elephant shark Rh2 and LWS1 pigments are the shortest wavelength-shifted pigments of their respective classes known to date. The mechanisms for this are discussed and we provide experimental evidence that the elephant shark LWS1 pigment uses a novel tuning mechanism to achieve the short wavelength shift to 499 nm, which inactivates the chloride-binding site. Our findings have important implications for the present knowledge of color vision evolution in early vertebrates.

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.

Background: Vasopressin and oxytocin are mammalian neurohypophysial hormones with distinct functions. Vasopressin is involved mainly in osmoregulation and oxytocin is involved primarily in parturition and lactation. Jawed vertebrates contain at least one homolog each of vasopressin and oxytocin, whereas only a vasopressin-family hormone, vasotocin, has been identified in jawless vertebrates. The genes encoding vasopressin and oxytocin are closely linked tail-to-tail in eutherian mammals whereas their homologs in chicken, Xenopus and coelacanth (vasotocin and mesotocin) are linked tail-to-head. In contrast, their pufferfish homologs, vasotocin and isotocin, are located on the same strand of DNA with isotocin located upstream of vasotocin and separated by five genes. These differences in the arrangement of the two genes in different bony vertebrate lineages raise questions about their origin and ancestral arrangement. To trace the origin of these genes, we have sequenced BAC clones from the neurohypophysial gene loci in a cartilaginous fish, the elephant shark (Callorhinchus milii), and in a jawless vertebrate, the Japanese lamprey (Lethenteron japonicum). We have also analyzed the neurohypophysial hormone gene locus in an invertebrate chordate, the amphioxus (Branchiostoma floridae).

The evolutionary history of horseshoe crabs, spanning approximately 500 million years, is characterized by remarkable morphological stasis and a low species diversity with only four extant species. Here we report a chromosome-level genome assembly for the mangrove horseshoe crab (Carcinoscorpius rotundicauda) using PacBio reads and Hi-C data. The assembly spans 1.67 Gb with contig N50 of 7.8 Mb and 98% of the genome assigned to 16 chromosomes. The genome contains five Hox clusters with 34 Hox genes, the highest number reported in any invertebrate. Detailed analysis of the genome provides evidence that suggests three rounds of whole-genome duplication (WGD), raising questions about the relationship between WGD and species radiation. Several gene families, particularly those involved in innate immunity, have undergone extensive tandem duplication. These expanded gene families may be important components of the innate immune system of horseshoe crabs, whose amebocyte lysate is a sensitive agent for detecting endotoxin contamination.