Enigmatic Nodal and Lefty gene repertoire discrepancy: Latent evolutionary history revealed by vertebrate-wide phylogeny

Homology in vertebrate body plans is traditionally ascribed to the high-level conservation of regulatory components within the genetic programs governing them, particularly during the "phylotypic stage." However, advancements in embryology and molecular phylogeny have unveiled the dynamic nature of gene repertoires responsible for early development. Notably, the Nodal and Lefty genes, members of the transforming growth factor-beta superfamily producing intercellular signaling molecules and crucial for left-right (L-R) symmetry breaking, exhibit distinctive features within their gene repertoires. These features encompass among-species gene repertoire variations resulting from gene gain and loss, as well as gene conversion. Despite their significance, these features have been largely unexplored in a phylogenetic context, but accumulating genome-wide sequence information is allowing the scrutiny of these features. It has exposed hidden paralogy between Nodal1 and Nodal2 genes resulting from differential gene loss in amniotes. In parallel, the tandem cluster of Lefty1 and Lefty2 genes, which was thought to be confined to mammals, is observed in sharks and rays, with an unexpected phylogenetic pattern. This article provides a comprehensive review of the current understanding of the origins of these vertebrate gene repertoires and proposes a revised nomenclature based on the elucidated history of vertebrate genome evolution.

Evolutionary origins of bitter taste receptors in jawed vertebrates

Taste is a sense that detects information about nutrients and toxins in foods. Of the five basic taste qualities, bitterness is associated with the detection of potentially harmful substances like plant alkaloids. In bony vertebrates, type 2 taste receptors (T2Rs), which are G-protein-coupled receptors (GPCRs), act as bitter taste receptors1,2. In vertebrates, six GPCR gene families are described as chemosensory receptor genes, encoding taste receptor families (T1Rs and T2Rs) and olfactory receptor families (ORs, V1Rs, V2Rs, and TAARs). These families of receptors have been found in all major jawed vertebrate lineages, except for the T2Rs, which are confined to bony vertebrates3. Therefore, T2Rs are believed to have emerged later than the other chemosensory receptor genes in the bony vertebrate lineage. So far, only the genomes of two cartilaginous fish species have been mined for TAS2R genes, which encode T2Rs4. Here, we identified novel T2Rs in elasmobranchs, namely selachimorphs (sharks) and batoids (rays, skates, and their close relatives) by an exhaustive search covering diverse cartilaginous fishes. Using functional and mRNA expression analyses, we demonstrate that their T2Rs are expressed in the oral taste buds and contribute to the detection of bitter compounds. This finding indicates the early origin of T2Rs in the common ancestor of jawed vertebrates.

Chromosomal DNA sequences of the Pacific saury genome: versatile resources for fishery science and comparative biology

Pacific saury (Cololabis saira) is a commercially important small pelagic fish species in Asia. In this study, we conducted the first-ever whole genome sequencing of this species, with single molecule, real-time (SMRT) sequencing technology. The obtained high-fidelity (HiFi) long-read sequence data, which amount to ~30-folds of its haploid genome size that was measured with quantitative PCR (1.17 Gb), were assembled into contigs. Scaffolding with Hi-C reads yielded a whole genome assembly containing 24 chromosome-scale sequences, with a scaffold N50 length of 47.7 Mb. Screening of repetitive elements including telomeric repeats was performed to characterize possible factors that need to be resolved towards 'telomere-to-telomere' sequencing. The larger genome size than in medaka, a close relative in Beloniformes, is at least partly explained by larger repetitive element quantity, which is reflected in more abundant tRNAs, in the Pacific saury genome. Protein-coding regions were predicted using transcriptome data, which resulted in 22,274 components. Retrieval of Pacific saury homologs of aquaporin (AQP) genes known from other teleost fishes validated high completeness and continuity of the genome assembly. These resources are available at https://treethinkers.nig.ac.jp/saira/ and will assist various molecular-level studies in fishery science and comparative biology.

Transcriptomic dysregulation and autistic-like behaviors in Kmt2c haploinsufficient mice rescued by an LSD1 inhibitor

Recent studies have consistently demonstrated that the regulation of chromatin and gene transcription plays a pivotal role in the pathogenesis of neurodevelopmental disorders. Among many genes involved in these pathways, KMT2C, encoding one of the six known histone H3 lysine 4 (H3K4) methyltransferases in humans and rodents, was identified as a gene whose heterozygous loss-of-function variants are causally associated with autism spectrum disorder (ASD) and the Kleefstra syndrome phenotypic spectrum. However, little is known about how KMT2C haploinsufficiency causes neurodevelopmental deficits and how these conditions can be treated. To address this, we developed and analyzed genetically engineered mice with a heterozygous frameshift mutation of Kmt2c (Kmt2c+/fs mice) as a disease model with high etiological validity. In a series of behavioral analyses, the mutant mice exhibit autistic-like behaviors such as impairments in sociality, flexibility, and working memory, demonstrating their face validity as an ASD model. To investigate the molecular basis of the observed abnormalities, we performed a transcriptomic analysis of their bulk adult brains and found that ASD risk genes were specifically enriched in the upregulated differentially expressed genes (DEGs), whereas KMT2C peaks detected by ChIP-seq were significantly co-localized with the downregulated genes, suggesting an important role of putative indirect effects of Kmt2c haploinsufficiency. We further performed single-cell RNA sequencing of newborn mouse brains to obtain cell type-resolved insights at an earlier stage. By integrating findings from ASD exome sequencing, genome-wide association, and postmortem brain studies to characterize DEGs in each cell cluster, we found strong ASD-associated transcriptomic changes in radial glia and immature neurons with no obvious bias toward upregulated or downregulated DEGs. On the other hand, there was no significant gross change in the cellular composition. Lastly, we explored potential therapeutic agents and demonstrate that vafidemstat, a lysine-specific histone demethylase 1 (LSD1) inhibitor that was effective in other models of neuropsychiatric/neurodevelopmental disorders, ameliorates impairments in sociality but not working memory in adult Kmt2c+/fs mice. Intriguingly, the administration of vafidemstat was shown to alter the vast majority of DEGs in the direction to normalize the transcriptomic abnormalities in the mutant mice (94.3 and 82.5% of the significant upregulated and downregulated DEGs, respectively, P < 2.2 × 10-16, binomial test), which could be the molecular mechanism underlying the behavioral rescuing. In summary, our study expands the repertoire of ASD models with high etiological and face validity, elucidates the cell-type resolved molecular alterations due to Kmt2c haploinsufficiency, and demonstrates the efficacy of an LSD1 inhibitor that might be generalizable to multiple categories of psychiatric disorders along with a better understanding of its presumed mechanisms of action.

A vertebrate-wide catalogue of T1R receptors reveals diversity in taste perception

Taste is a vital chemical sense for feeding behaviour. In mammals, the umami and sweet taste receptors comprise three members of the taste receptor type 1 (T1R/TAS1R) family: T1R1, T1R2 and T1R3. Because their functional homologues exist in teleosts, only three TAS1R genes generated by gene duplication are believed to have been inherited from the common ancestor of bony vertebrates. Here, we report five previously uncharacterized TAS1R members in vertebrates, TAS1R4, TAS1R5, TAS1R6, TAS1R7 and TAS1R8, based on genome-wide survey of diverse taxa. We show that mammalian and teleost fish TAS1R2 and TAS1R3 genes are paralogues. Our phylogenetic analysis suggests that the bony vertebrate ancestor had nine TAS1Rs resulting from multiple gene duplications. Some TAS1Rs were lost independently in descendent lineages resulting in retention of only three TAS1Rs in mammals and teleosts. Combining functional assays and expression analysis of non-teleost fishes we show that the novel T1Rs form heterodimers in taste-receptor cells and recognize a broad range of ligands such as essential amino acids, including branched-chain amino acids, which have not been previously considered as T1R ligands. This study reveals diversity of taste sensations in both modern vertebrates and their ancestors, which might have enabled vertebrates to adapt to diverse habitats on Earth.

Truncated radial glia as a common precursor in the late corticogenesis of gyrencephalic mammals

The diversity of neural stem cells is a hallmark of the cerebral cortex development in gyrencephalic mammals, such as Primates and Carnivora. Among them, ferrets are a good model for mechanistic studies. However, information on their neural progenitor cells (NPC), termed radial glia (RG), is limited. Here, we surveyed the temporal series of single-cell transcriptomes of progenitors regarding ferret corticogenesis and found a conserved diversity and temporal trajectory between human and ferret NPC, despite the large timescale difference. We found truncated RG (tRG) in ferret cortical development, a progenitor subtype previously described in humans. The combination of in silico and in vivo analyses identified that tRG differentiate into both ependymal and astrogenic cells. Via transcriptomic comparison, we predict that this is also the case in humans. Our findings suggest that tRG plays a role in the formation of adult ventricles, thereby providing the architectural bases for brain expansion.

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.

Control of Synaptic Levels of Nicotinic Acetylcholine Receptor by the Sequestering Subunit Dα5 and Secreted Scaffold Protein Hig

The presentation of nicotinic acetylcholine receptors (nAChRs) on synaptic membranes is crucial for generating cholinergic circuits, some of which are associated with memory function and neurodegenerative disorders. Although the physiology and structure of nAChR, a cation channel comprising five subunits, have been extensively studied, little is known about how the receptor levels in interneuronal synapses are determined and which nAChR subunits participate in the regulatory process in cooperation with synaptic cleft matrices and intracellular proteins. By a genetic screen of Drosophila, we identified mutations in the nAChR subunit Dα5 gene as suppressors that restored the mutant phenotypes of hig, which encodes a secretory matrix protein localized to cholinergic synaptic clefts in the brain. Only the loss of function of Dα5 among the 10 nAChR subunits suppressed hig mutant phenotypes in both male and female flies. Dα5 behaved as a lethal factor when Hig was defective; loss of Dα5 in hig mutants rescued lethality, upregulating Dα6 synaptic levels. By contrast, levels of Dα5, Dα6, and Dα7 subunits were all reduced in hig mutants. These three subunits have distinct properties for interaction with Hig or trafficking, as confirmed by chimeric subunit experiments. Notably, the chimeric Dα5 protein, which has the extracellular sequences that display no positive interaction with Hig, exhibited abnormal distribution and lethality even in the presence of Hig. We propose that the sequestering subunit Dα5 functions by reducing synaptic levels of nAChR through internalization, and this process is blocked by Hig, which tethers Dα5 to the synaptic cleft matrix.SIGNIFICANCE STATEMENT Because the cholinergic synapse is one of the major synapses that generate various brain functions, numerous studies have sought to reveal the physiology and structure of the nicotinic acetylcholine receptor (nAChR). However, little is known about how synaptic levels of nAChR are controlled and which nAChR subunits participate in the regulatory process in cooperation with synaptic cleft matrices. By a genetic screen of Drosophila, we identified mutations in the nAChR subunit Dα5 gene as suppressors that restored the mutant phenotypes of hig, which encodes a secretory matrix protein localized to cholinergic synaptic clefts. Our data indicate that Dα5 functions in reducing synaptic levels of nAChR, and this process is blocked by Hig, which tethers Dα5 to the synaptic cleft matrix.

The impact of local genomic properties on the evolutionary fate of genes

Functionally indispensable genes are likely to be retained and otherwise to be lost during evolution. This evolutionary fate of a gene can also be affected by factors independent of gene dispensability, including the mutability of genomic positions, but such features have not been examined well. To uncover the genomic features associated with gene loss, we investigated the characteristics of genomic regions where genes have been independently lost in multiple lineages. With a comprehensive scan of gene phylogenies of vertebrates with a careful inspection of evolutionary gene losses, we identified 813 human genes whose orthologs were lost in multiple mammalian lineages: designated 'elusive genes.' These elusive genes were located in genomic regions with rapid nucleotide substitution, high GC content, and high gene density. A comparison of the orthologous regions of such elusive genes across vertebrates revealed that these features had been established before the radiation of the extant vertebrates approximately 500 million years ago. The association of human elusive genes with transcriptomic and epigenomic characteristics illuminated that the genomic regions containing such genes were subject to repressive transcriptional regulation. Thus, the heterogeneous genomic features driving gene fates toward loss have been in place and may sometimes have relaxed the functional indispensability of such genes. This study sheds light on the complex interplay between gene function and local genomic properties in shaping gene evolution that has persisted since the vertebrate ancestor.

Phylogenetic and functional properties of hagfish neurohypophysial hormone receptors distinct from their jawed vertebrate counterparts

Vertebrate neurohypophysial hormones, i.e., vasopressin- and oxytocin-family peptides, exert versatile physiological actions via distinct G protein-coupled receptors. The neurohypophysial hormone receptor (NHR) family was classically categorized into four subtypes (V1aR, V1bR, V2R and OTR), while recent studies have identified seven subtypes (V1aR, V1bR, V2aR, V2bR, V2cR, V2dR and OTR; V2aR corresponds to the conventional V2R). The vertebrate NHR family were diversified via multiple gene duplication events at different scales. Despite intensive research effort in non-osteichthyes vertebrates such as cartilaginous fish and lamprey, the molecular phylogeny of the NHR family has not been fully understood. In the present study, we focused on the inshore hagfish (Eptatretus burgeri), another group of cyclostomes, and Arctic lamprey (Lethenteron camtschaticum) for comparison. Two putative NHR homologs, which were previously identified only in silico, were cloned from the hagfish and designated as ebV1R and ebV2R. In vitro, ebV1R, as well as two out of five Arctic lamprey NHRs, increased intracellular Ca2+ in response to exogenous neurohypophysial hormones. None of the examined cyclostome NHRs altered intracellular cAMP levels. Transcripts of ebV1R were detected in multiple tissues including the brain and gill, with intense hybridization signals in the hypothalamus and adenohypophysis, while ebV2R was predominantly expressed in the systemic heart. Similarly, Arctic lamprey NHRs showed distinct expression patterns, underscoring the multifunctionality of VT in the cyclostomes as in the gnathostomes. These results and exhaustive gene synteny comparisons provide new insights into the molecular and functional evolution of the neurohypophysial hormone system in vertebrates.

Gse1, a component of the CoREST complex, is required for placenta development in the mouse

Gse1 is a component of the CoREST complex that acts as an H3K4 and H3K9 demethylase and regulates gene expression. Here, we examined the expression and role of Gse1 in mouse development. Gse1 is expressed in male and female germ cells and plays both maternal and zygotic roles. Thus, maternal deletion of Gse1 results in a high incidence of prenatal death, and zygotic deletion leads to embryonic lethality from embryonic day 12.5 (E12.5) and perinatal death. Gse1 is expressed in the junctional zone and the labyrinth of the developing placenta. Cultured trophoblast stem cells lacking Gse1 showed impaired in vitro cell differentiation into spongiotrophoblasts. Gse1 mutant (Gse1^Δex3/Δex3) placenta begins to exhibit histological defects from E14.5, being deficient in MCT4⁺ syncytiotrophoblast II. The number of various cell types was largely maintained in the mutant placenta at E10.5, but several genes were upregulated in giant trophoblasts at E10.5. Placenta-specific deletion of Gse1 with Tat-Cre suggested that defects in Gse1^Δex3/Δex3 embryos are due to placental function deficiency. These results suggest that Gse1 is required for placental development in mice, and in turn, is essential for embryonic development.

Whale shark rhodopsin adapted to deep-sea lifestyle by a substitution associated with human disease

Spectral tuning of visual pigments often facilitates adaptation to new environments, and it is intriguing to study the visual ecology of pelagic sharks with secondarily expanded habitats. The whale shark, which dives into the deep sea of nearly 2,000 meters besides near-surface filter feeding, was previously shown to possess the 'blue-shifted' rhodopsin (RHO), which is a signature of deep-sea adaptation. In this study, our spectroscopy of recombinant whale shark RHO mutants revealed that this blue shift is caused dominantly by an unprecedented spectral tuning site 94. In humans, the mutation at the site causes congenital stationary night blindness (CSNB) by reducing the thermal stability of RHO. Similarly, the RHO of deep-diving whale shark has reduced thermal stability, which was experimentally shown to be achieved by site 178 and 94. RHOs having the natural substitution at site 94 are also found in some Antarctic fishes, suggesting that the blue shift by the substitution at the CSNB site associated with the reduction in thermal stability might be allowed in cold-water deep-sea habitats.

Egg Yolk Protein Homologs Identified in Live-Bearing Sharks: Co-Opted in the Lecithotrophy-to-Matrotrophy Shift?

Reproductive modes of vertebrates are classified into two major embryonic nutritional types: yolk deposits (i.e., lecithotrophy) and maternal investment (i.e., matrotrophy). Vitellogenin (VTG), a major egg yolk protein synthesized in the female liver, is one of the molecules relevant to the lecithotrophy-to-matrotrophy shift in bony vertebrates. In mammals, all VTG genes are lost following the lecithotrophy-to-matrotrophy shift, and it remains to be elucidated whether the lecithotrophy-to-matrotrophy shift in nonmammalians is also associated with VTG repertoire modification. In this study, we focused on chondrichthyans (cartilaginous fishes)-a vertebrate clade that underwent multiple lecithotrophy-to-matrotrophy shifts. For an exhaustive search of homologs, we performed tissue-by-tissue transcriptome sequencing for two viviparous chondrichthyans, the frilled shark Chlamydoselachus anguineus and the spotless smooth-hound Mustelus griseus, and inferred the molecular phylogeny of VTG and its receptor very low-density lipoprotein receptor (VLDLR), across diverse vertebrates. As a result, we identified either three or four VTG orthologs in chondrichthyans including viviparous species. We also showed that chondrichthyans had two additional VLDLR orthologs previously unrecognized in their unique lineage (designated as VLDLRc2 and VLDLRc3). Notably, VTG gene expression patterns differed in the species studied depending on their reproductive mode; VTGs are broadly expressed in multiple tissues, including the uterus, in the two viviparous sharks, and in addition to the liver. This finding suggests that the chondrichthyans VTGs do not only function as the yolk nutrient but also as the matrotrophic factor. Altogether, our study indicates that the lecithotrophy-to-matrotrophy shift in chondrichthyans was achieved through a distinct evolutionary process from mammals.

Transcriptomic exploration of the Coleopteran wings reveals insight into the evolution of novel structures associated with the beetle elytron

The acquisition of novel traits is central to organismal evolution, yet the molecular mechanisms underlying this process are elusive. The beetle forewings (elytra) are evolutionarily modified to serve as a protective shield, providing a unique opportunity to study these mechanisms. In the past, the orthologs of genes within the wing gene network from Drosophila studies served as the starting point when studying the evolution of elytra (candidate genes). Although effective, candidate gene lists are finite and only explore genes conserved across species. To go beyond candidate genes, we used RNA sequencing and explored the wing transcriptomes of two Coleopteran species, the red flour beetle (Tribolium castaneum) and the Japanese stag beetle (Dorcus hopei). Our analysis revealed sets of genes enriched in Tribolium elytra (57 genes) and genes unique to the hindwings, which possess more "typical" insect wing morphologies (29 genes). Over a third of the hindwing-enriched genes were "candidate genes" whose functions were previously analyzed in Tribolium, demonstrating the robustness of our sequencing. Although the overlap was limited, transcriptomic comparison between the beetle species found a common set of genes, including key wing genes, enriched in either elytra or hindwings. Our RNA interference analysis for elytron-enriched genes in Tribolium uncovered novel genes with roles in forming various aspects of morphology that are unique to elytra, such as pigmentation, hardening, sensory development, and vein formation. Our analyses deepen our understanding of how gene network evolution facilitated the emergence of the elytron, a unique structure critical to the evolutionary success of beetles.

Evolution of Vertebrate Hormones and Their Receptors: Insights from Non-Osteichthyan Genomes

Homeostatic control and reproductive functions of humans are regulated at the molecular levels largely by peptide hormones secreted from endocrine and/or neuroendocrine cells in the central nervous system and peripheral organs. Homologs of those hormones and their receptors function similarly in many vertebrate species distantly related to humans, but the evolutionary history of the endocrine system involving those factors has been obscured by the scarcity of genome DNA sequence information of some taxa that potentially contain their orthologs. Focusing on non-osteichthyan vertebrates, namely jawless and cartilaginous fishes, this article illustrates how investigating genome sequence information assists our understanding of the diversification of vertebrate gene repertoires in four broad themes: (a) the presence or absence of genes, (b) multiplication and maintenance of paralogs, (c) differential fates of duplicated paralogs, and (d) the evolutionary timing of gene origins.

Squalomix: shark and ray genome analysis consortium and its data sharing platform

The taxon Elasmobranchii (sharks and rays) contains one of the long-established evolutionary lineages of vertebrates with a tantalizing collection of species occupying critical aquatic habitats. To overcome the current limitation in molecular resources, we launched the Squalomix Consortium in 2020 to promote a genome-wide array of molecular approaches, specifically targeting shark and ray species. Among the various bottlenecks in working with elasmobranchs are their elusiveness and low fecundity as well as the large and highly repetitive genomes. Their peculiar body fluid composition has also hindered the establishment of methods to perform routine cell culturing required for their karyotyping. In the Squalomix consortium, these obstacles are expected to be solved through a combination of in-house cytological techniques including karyotyping of cultured cells, chromatin preparation for Hi-C data acquisition, and high fidelity long-read sequencing. The resources and products obtained in this consortium, including genome and transcriptome sequences, a genome browser powered by JBrowse2 to visualize sequence alignments, and comprehensive matrices of gene expression profiles for selected species are accessible through https://github.com/Squalomix/info.

Molecular and morphological investigations on the renal mechanisms enabling euryhalinity of red stingray Hemitrygon akajei

Most cartilaginous fishes live in seawater (SW), but a few exceptional elasmobranchs (sharks and rays) are euryhaline and can acclimate to freshwater (FW) environments. The plasma of elasmobranchs is high in NaCl and urea concentrations, which constrains osmotic water loss. However, these euryhaline elasmobranchs maintain high levels of plasma NaCl and urea even when acclimating to low salinity, resulting in a strong osmotic gradient from external environment to body fluid. The kidney consequently produces a large volume of dilute urine to cope with the water influx. In the present study, we investigated the molecular mechanisms of dilute urine production in the kidney of Japanese red stingray, Hemitrygon akajei, transferred from SW to low-salinity environments. We showed that red stingray maintained high plasma NaCl and urea levels by reabsorbing more osmolytes in the kidney when transferred to low salinity. RNA-seq and qPCR analyses were conducted to identify genes involved in NaCl and urea reabsorption under the low-salinity conditions, and the upregulated gene expressions of Na⁺-K⁺-Cl^- cotransporter 2 (nkcc2) and Na⁺/K⁺-ATPase (nka) were found in the FW-acclimated individuals. These upregulations occurred in the early distal tubule (EDT) in the bundle zone of the kidney, which coils around the proximal and collecting tubules to form the highly convoluted structure of batoid nephron. Considering the previously proposed model for urea reabsorption, the upregulation of nkcc2 and nka not only causes the reabsorption of NaCl in the EDT, but potentially also supports enhanced urea reabsorption and eventually the production of dilute urine in FW-acclimated individuals. We propose advantageous characteristics of the batoid-type nephron that facilitate acclimation to a wide range of salinities, which might have allowed the batoids to expand their habitats.

Why sequence all eukaryotes?

Life on Earth has evolved from initial simplicity to the astounding complexity we experience today. Bacteria and archaea have largely excelled in metabolic diversification, but eukaryotes additionally display abundant morphological innovation. How have these innovations come about and what constraints are there on the origins of novelty and the continuing maintenance of biodiversity on Earth? The history of life and the code for the working parts of cells and systems are written in the genome. The Earth BioGenome Project has proposed that the genomes of all extant, named eukaryotes-about 2 million species-should be sequenced to high quality to produce a digital library of life on Earth, beginning with strategic phylogenetic, ecological, and high-impact priorities. Here we discuss why we should sequence all eukaryotic species, not just a representative few scattered across the many branches of the tree of life. We suggest that many questions of evolutionary and ecological significance will only be addressable when whole-genome data representing divergences at all of the branchings in the tree of life or all species in natural ecosystems are available. We envisage that a genomic tree of life will foster understanding of the ongoing processes of speciation, adaptation, and organismal dependencies within entire ecosystems. These explorations will resolve long-standing problems in phylogenetics, evolution, ecology, conservation, agriculture, bioindustry, and medicine.

Whole-Genome Duplications and the Diversification of the Globin-X Genes of Vertebrates

Globin-X (GbX) is an enigmatic member of the vertebrate globin gene family with a wide phyletic distribution that spans protostomes and deuterostomes. Unlike canonical globins such as hemoglobins and myoglobins, functional data suggest that GbX does not have a primary respiratory function. Instead, evidence suggests that the monomeric, membrane-bound GbX may play a role in cellular signaling or protection against the oxidation of membrane lipids. Recently released genomes from key vertebrates provide an excellent opportunity to address questions about the early stages of the evolution of GbX in vertebrates. We integrate bioinformatics, synteny, and phylogenetic analyses to characterize the diversity of GbX genes in nonteleost ray-finned fishes, resolve relationships between the GbX genes of cartilaginous fish and bony vertebrates, and demonstrate that the GbX genes of cyclostomes and gnathostomes derive from independent duplications. Our study highlights the role that whole-genome duplications (WGDs) have played in expanding the repertoire of genes in vertebrate genomes. Our results indicate that GbX paralogs have a remarkably high rate of retention following WGDs relative to other globin genes and provide an evolutionary framework for interpreting results of experiments that examine functional properties of GbX and patterns of tissue-specific expression. By identifying GbX paralogs that are products of different WGDs, our results can guide the design of experimental work to explore whether gene duplicates that originate via WGDs have evolved novel functional properties or expression profiles relative to singleton or tandemly duplicated copies of GbX.

Genetic Mechanism for the Cyclostome Cerebellar Neurons Reveals Early Evolution of the Vertebrate Cerebellum

The vertebrate cerebellum arises at the dorsal part of rhombomere 1, induced by signals from the isthmic organizer. Two major cerebellar neuronal subtypes, granule cells (excitatory) and Purkinje cells (inhibitory), are generated from the anterior rhombic lip and the ventricular zone, respectively. This regionalization and the way it develops are shared in all extant jawed vertebrates (gnathostomes). However, very little is known about early evolution of the cerebellum. The lamprey, an extant jawless vertebrate lineage or cyclostome, possesses an undifferentiated, plate-like cerebellum, whereas the hagfish, another cyclostome lineage, is thought to lack a cerebellum proper. In this study, we found that hagfish Atoh1 and Wnt1 genes are co-expressed in the rhombic lip, and Ptf1a is expressed ventrally to them, confirming the existence of r1's rhombic lip and the ventricular zone in cyclostomes. In later stages, lamprey Atoh1 is downregulated in the posterior r1, in which the NeuroD increases, similar to the differentiation process of cerebellar granule cells in gnathostomes. Also, a continuous Atoh1-positive domain in the rostral r1 is reminiscent of the primordium of valvula cerebelli of ray-finned fishes. Lastly, we detected a GAD-positive domain adjacent to the Ptf1a-positive ventricular zone in lampreys, suggesting that the Ptf1a-positive cells differentiate into some GABAergic inhibitory neurons such as Purkinje and other inhibitory neurons like in gnathostomes. Altogether, we conclude that the ancestral genetic programs for the formation of a distinct cerebellum were established in the last common ancestor of vertebrates.

Technical considerations in Hi-C scaffolding and evaluation of chromosome-scale genome assemblies

The recent development of ecological studies has been fueled by the introduction of massive information based on chromosome‐scale genome sequences, even for species for which genetic linkage is not accessible. This was enabled mainly by the application of Hi‐C, a method for genome‐wide chromosome conformation capture that was originally developed for investigating the long‐range interaction of chromatins. Performing genomic scaffolding using Hi‐C data is highly resource‐demanding and employs elaborate laboratory steps for sample preparation. It starts with building a primary genome sequence assembly as an input, which is followed by computation for genome scaffolding using Hi‐C data, requiring careful validation. This article presents technical considerations for obtaining optimal Hi‐C scaffolding results and provides a test case of its application to a reptile species, the Madagascar ground gecko (Paroedura picta). Among the metrics that are frequently used for evaluating scaffolding results, we investigate the validity of the completeness assessment of chromosome‐scale genome assemblies using single‐copy reference orthologues.

The whale shark genome reveals patterns of vertebrate gene family evolution

Chondrichthyes (cartilaginous fishes) are fundamental for understanding vertebrate evolution, yet their genomes are understudied. We report long-read sequencing of the whale shark genome to generate the best gapless chondrichthyan genome assembly yet with higher contig contiguity than all other cartilaginous fish genomes, and studied vertebrate genomic evolution of ancestral gene families, immunity, and gigantism. We found a major increase in gene families at the origin of gnathostomes (jawed vertebrates) independent of their genome duplication. We studied vertebrate pathogen recognition receptors (PRRs), which are key in initiating innate immune defense, and found diverse patterns of gene family evolution, demonstrating that adaptive immunity in gnathostomes did not fully displace germline-encoded PRR innovation. We also discovered a new toll-like receptor (TLR29) and three NOD1 copies in the whale shark. We found chondrichthyan and giant vertebrate genomes had decreased substitution rates compared to other vertebrates, but gene family expansion rates varied among vertebrate giants, suggesting substitution and expansion rates of gene families are decoupled in vertebrate genomes. Finally, we found gene families that shifted in expansion rate in vertebrate giants were enriched for human cancer-related genes, consistent with gigantism requiring adaptations to suppress cancer.

Chromosomal-scale de novo genome assemblies of Cynomolgus Macaque and Common Marmoset

Cynomolgus macaque (Macaca fascicularis) and common marmoset (Callithrix jacchus) have been widely used in human biomedical research. Long-standing primate genome assemblies used the human genome as a reference for ordering and orienting the assembled fragments into chromosomes. Here we performed de novo genome assembly of these two species without any human genome-based bias observed in the genome assemblies released earlier. We assembled PacBio long reads, and the resultant contigs were scaffolded with Hi-C data, which were further refined based on Hi-C contact maps and alternate de novo assemblies. The assemblies achieved scaffold N50 lengths of 149 Mb and 137 Mb for cynomolgus macaque and common marmoset, respectively. The high fidelity of our assembly is also ascertained by BAC-end concordance in common marmoset. Our assembly of cynomolgus macaque outperformed all the available assemblies of this species in terms of contiguity. The chromosome-scale genome assemblies produced in this study are valuable resources for non-human primate models and provide an important baseline in human biomedical research.

Cooperation of LIM domain-binding 2 (LDB2) with EGR in the pathogenesis of schizophrenia

Genomic defects with large effect size can help elucidate unknown pathologic architecture of mental disorders. We previously reported on a patient with schizophrenia and a balanced translocation between chromosomes 4 and 13 and found that the breakpoint within chromosome 4 is located near the LDB2 gene. We show here that Ldb2 knockout (KO) mice displayed multiple deficits relevant to mental disorders. In particular, Ldb2 KO mice exhibited deficits in the fear-conditioning paradigm. Analysis of the amygdala suggested that dysregulation of synaptic activities controlled by the immediate early gene Arc is involved in the phenotypes. We show that LDB2 forms protein complexes with known transcription factors. Consistently, ChIP-seq analyses indicated that LDB2 binds to > 10,000 genomic sites in human neurospheres. We found that many of those sites, including the promoter region of ARC, are occupied by EGR transcription factors. Our previous study showed an association of the EGR family genes with schizophrenia. Collectively, the findings suggest that dysregulation in the gene expression controlled by the LDB2-EGR axis underlies a pathogenesis of subset of mental disorders.

Measuring potential effects of the developmental burden associated with the vertebrate notochord

The notochord functions primarily as a supporting tissue to maintain the anteroposterior axis of primitive chordates, a function that is replaced entirely by the vertebral column in many vertebrates. The notochord still appears during vertebrate embryogenesis and plays a crucial role in the developmental pattern formation of surrounding structures, such as the somites and neural tube, providing the basis for the vertebrate body plan. The indispensable role of the notochord has often been referred to as the developmental burden and used to explain the evolutionary conservation of notochord; however, the existence of this burden has not been successfully exemplified so far. Since the adaptive value of target tissues appears to result in the evolutionary conservation of upstream structures through the developmental burden, we performed comparative gene expression profiling of the notochord, somites, and neural tube during the mid‐embryonic stages in turtles and chicken to measure their evolutionary conservation. When compared with the somites and neural tube, overall gene expression profiles in the notochord showed significantly lower or merely comparable levels of conservation. However, genes involved in inductive signalings, such as the sonic hedgehog (Shh) cascade and the formation of functional primary cilia, showed relatively higher levels of conservation in all the three structures analyzed. Collectively, these results suggest that shh signals are critical as the inductive source and receiving structures, possibly constituting the inter‐dependencies of developmental burden.

Potential evolutionary effects toward notochord by developmental burden was evaluated by Laser Micro Dissection RNAseq (LMDseq). Notochord was less conserved than neural tube and somites; however, genes in sonic hedgehog (shh) signaling cascade was found to be evolutionarily conserved (not only in notochord but also in somites and neural tube). These results suggest that Shh signals are critical as the inductive source and receiving structures, possibly constituting the inter‐dependencies of developmental burden. Further studies that directly measure the burden required to verify the hypothesis are awaited.

Developmental hourglass and heterochronic shifts in fin and limb development

How genetic changes are linked to morphological novelties and developmental constraints remains elusive. Here, we investigate genetic apparatuses that distinguish fish fins from tetrapod limbs by analyzing transcriptomes and open-chromatin regions (OCRs). Specifically, we compared mouse forelimb buds with the pectoral fin buds of an elasmobranch, the brown-banded bamboo shark (Chiloscyllium punctatum). A transcriptomic comparison with an accurate orthology map revealed both a mass heterochrony and hourglass-shaped conservation of gene expression between fins and limbs. Furthermore, open-chromatin analysis suggested that access to conserved regulatory sequences is transiently increased during mid-stage limb development. During this stage, stage-specific and tissue-specific OCRs were also enriched. Together, early and late stages of fin/limb development are more permissive to mutations than middle stages, which may have contributed to major morphological changes during the fin-to-limb evolution. We hypothesize that the middle stages are constrained by regulatory complexity that results from dynamic and tissue-specific transcriptional controls.

Animals come in all shapes and sizes. This diversity arose through genetic mutations during evolution, but it is unclear exactly how these variations led to the formation of new shapes. There is increasing evidence to suggest that not all shapes are possible and that variability between animals is limited by a phenomenon known as “developmental constraint”. These limitations direct parts of the body towards a specific shape as they develop in the embryo. Therefore, understanding the mechanisms underlying these developmental constraints could help explain how different body shapes evolved.

The limbs of humans and other mammals evolved from the fins of fish, and this transition is often used to study the role developmental constraints play in evolution. This is an ideal model as there is already a detailed fossil record mapping this evolutionary event, and data pinpointing some of the genes involved in the development of limbs and fins. But this data is incomplete, and a full comparison between the genes activated in the fin and the limb during embryonic development had not been achieved. This is because most fish used for research have undergone recent genetic changes, making it hard to spot which genetic differences are linked to the evolution of the limb.

To overcome this barrier, Onimaru et al. compared genetic data from the developing limbs of mice to the developing fins of the brown-banded bamboo shark, which evolves much slower than other fish. This revealed that although many genes commonly played a role in the development of the fin and the limb in the embryo, the activity of these shared genes was not the same. For example, genes that switched on in the late stages of limb development, switched off in the late stages of fin development. But in the middle of development, those differences were relatively small and both species activated very similar sets of genes. Many of these genes were pleiotropic, which means they have important roles in other tissues and therefore mutate less often. This suggests that the mid-stage of limb development is under the strongest level of constraint.

Darwin’s theory of natural selection explains that mutations drive evolution. But the theory cannot predict what kinds of new body shapes new mutations will produce. Understanding how the activity levels of different genes affect development could help to fill this knowledge gap. This has potential medical applications, for example, understanding why some genetic changes cause more serious problems than others. This work suggests that mutations in genes that are active during the mid-stage of limb development may have the most serious impact.

The comparative genomic landscape of adaptive radiation in crater lake cichlid fishes

Factors ranging from ecological opportunity to genome composition might explain why only some lineages form adaptive radiations. While being rare, particular systems can provide natural experiments within an identical ecological setting where species numbers and phenotypic divergence in two closely related lineages are notably different. We investigated one such natural experiment using two de novo assembled and 40 resequenced genomes and asked why two closely related Neotropical cichlid fish lineages, the Amphilophus citrinellus species complex (Midas cichlids; radiating) and Archocentrus centrarchus (Flyer cichlid; nonradiating), have resulted in such disparate evolutionary outcomes. Although both lineages inhabit many of the same Nicaraguan lakes, whole-genome inferred demography suggests that priority effects are not likely to be the cause of the dissimilarities. Also, genome-wide levels of selection, transposable element dynamics, gene family expansion, major chromosomal rearrangements and the number of genes under positive selection were not markedly different between the two lineages. To more finely investigate particular subsets of the genome that have undergone adaptive divergence in Midas cichlids, we also examined if there was evidence for 'molecular pre-adaptation' in regions identified by QTL mapping of repeatedly diverging adaptive traits. Although most of our analyses failed to pinpoint substantial genomic differences, we did identify functional categories containing many genes under positive selection that provide candidates for future studies on the propensity of Midas cichlids to radiate. Our results point to a disproportionate role of local, rather than genome-wide factors underlying the propensity for these cichlid fishes to adaptively radiate.

Cell culture-based karyotyping of orectolobiform sharks for chromosome-scale genome analysis

Karyotyping, traditionally performed using cytogenetic techniques, is indispensable for validating genome assemblies whose sequence lengths can be scaled up to chromosome sizes using modern methods. Karyotype reports of chondrichthyans are scarce because of the difficulty in cell culture. Here, we focused on carpet shark species and the culture conditions for fibroblasts and lymphocytes. The utility of the cultured cells enabled the high-fidelity characterization of their karyotypes, namely 2n = 102 for the whale shark (Rhincodon typus) and zebra shark (Stegostoma fasciatum), and 2n = 106 for the brownbanded bamboo shark (Chiloscyllium punctatum) and whitespotted bamboo shark (C. plagiosum). We identified heteromorphic XX/XY sex chromosomes for the two latter species and demonstrated the first-ever fluorescence in situ hybridization of shark chromosomes prepared from cultured cells. Our protocols are applicable to diverse chondrichthyan species and will deepen the understanding of early vertebrate evolution at the molecular level.

Multifaceted Hi-C benchmarking: what makes a difference in chromosome-scale genome scaffolding?

Background

Hi-C is derived from chromosome conformation capture (3C) and targets chromatin contacts on a genomic scale. This method has also been used frequently in scaffolding nucleotide sequences obtained by de novo genome sequencing and assembly, in which the number of resultant sequences rarely converges to the chromosome number. Despite its prevalent use, the sample preparation methods for Hi-C have not been intensively discussed, especially from the standpoint of genome scaffolding.

Results

To gain insight into the best practice of Hi-C scaffolding, we performed a multifaceted methodological comparison using vertebrate samples and optimized various factors during sample preparation, sequencing, and computation. As a result, we identified several key factors that helped improve Hi-C scaffolding, including the choice and preparation of tissues, library preparation conditions, the choice of restriction enzyme(s), and the choice of scaffolding program and its usage.

Conclusions

This study provides the first comparison of multiple sample preparation kits/protocols and computational programs for Hi-C scaffolding by an academic third party. We introduce a customized protocol designated “inexpensive and controllable Hi-C (iconHi-C) protocol,” which incorporates the optimal conditions identified in this study, and demonstrate this technique on chromosome-scale genome sequences of the Chinese softshell turtle Pelodiscus sinensis.

Facilitated NaCl Uptake in the Highly Developed Bundle of the Nephron in Japanese Red Stingray Hemitrygon akajei Revealed by Comparative Anatomy and Molecular Mapping

Batoidea (rays and skates) is a monophyletic subgroup of elasmobranchs that diverged from the common ancestor with Selachii (sharks) about 270 Mya. A larger number of batoids can adapt to low-salinity environments, in contrast to sharks, which are mostly stenohaline marine species. Among osmoregulatory organs of elasmobranchs, the kidney is known to be dedicated to urea retention in ureosmotic cartilaginous fishes. However, we know little regarding urea reabsorbing mechanisms in the kidney of batoids. Here, we performed physiological and histological investigations on the nephrons in the red stingray (Hemitrygon akajei) and two shark species. We found that the urine/plasma ratios of salt and urea concentrations in the stingray are significantly lower than those in cloudy catshark (Scyliorhinus torazame) under natural seawater, indicating that the kidney of stingray more strongly reabsorbs these osmolytes. By comparing the three-dimensional images of nephrons between stingray and banded houndshark (Triakis scyllium), we showed that the tubular bundle of stingray has a more compact configuration. In the compact tubular bundle of stingray kidney, the distal diluting tubule was highly developed and frequently coiled around the proximal and collecting tubules. Furthermore, co-expression of NKAα1 (Na⁺/K ⁺-ATPase) and NKCC2 (Na⁺- K⁺-2Cl^- cotransporter 2) mRNAs was prominent in the coiled diluting segment. These findings imply that NaCl reabsorption is greatly facilitated in the stingray kidney, resulting in a higher reabsorption rate of urea. Lowering the loss of osmolytes in the glomerular filtrate is likely favorable to the adaptability of batoids to a wide range of environmental salinity.

Predicting gene regulatory regions with a convolutional neural network for processing double-strand genome sequence information

With advances in sequencing technology, a vast amount of genomic sequence information has become available. However, annotating biological functions particularly of non-protein-coding regions in genome sequences without experiments is still a challenging task. Recently deep learning-based methods were shown to have the ability to predict gene regulatory regions from genome sequences, promising to aid the interpretation of genomic sequence data. Here, we report an improvement of the prediction accuracy for gene regulatory regions by using the design of convolution layers that efficiently process genomic sequence information, and developed a software, DeepGMAP, to train and compare different deep learning-based models (https://github.com/koonimaru/DeepGMAP). First, we demonstrate that our convolution layers, termed forward- and reverse-sequence scan (FRSS) layers, integrate both forward and reverse strand information, and enhance the power to predict gene regulatory regions. Second, we assessed previous studies and identified problems associated with data structures that caused overfitting. Finally, we introduce visualization methods to examine what the program learned. Together, our FRSS layers improve the prediction accuracy for gene regulatory regions.

Nodal paralogues underlie distinct mechanisms for visceral left-right asymmetry in reptiles and mammals

Unidirectional fluid flow generated by motile cilia at the left-right organizer (LRO) breaks left-right (L-R) symmetry during early embryogenesis in mouse, frog and zebrafish. The chick embryo, however, does not require motile cilia for L-R symmetry breaking. The diversity of mechanisms for L-R symmetry breaking among vertebrates and the trigger for such symmetry breaking in non-mammalian amniotes have remained unknown. Here we examined how L-R asymmetry is established in two reptiles, Madagascar ground gecko and Chinese softshell turtle. Both of these reptiles appear to lack motile cilia at the LRO. The expression of the Nodal gene at the LRO in the reptilian embryos was found to be asymmetric, in contrast to that in vertebrates such as mouse that are dependent on cilia for L-R patterning. Two paralogues of the Nodal gene derived from an ancient gene duplication are retained and expressed differentially in cilia-dependent and cilia-independent vertebrates. The expression of these two Nodal paralogues is similarly controlled in the lateral plate mesoderm but regulated differently at the LRO. Our in-depth analysis of reptilian embryos thus suggests that mammals and non-mammalian amniotes deploy distinct strategies dependent on different Nodal paralogues for rendering Nodal activity asymmetric at the LRO.

Loss of olfaction in sea snakes provides new perspectives on the aquatic adaptation of amniotes

Marine amniotes, a polyphyletic group, provide an excellent opportunity for studying convergent evolution. Their sense of smell tends to degenerate, but this process has not been explored by comparing fully aquatic species with their amphibious relatives in an evolutionary context. Here, we sequenced the genomes of fully aquatic and amphibious sea snakes and identified repertoires of chemosensory receptor genes involved in olfaction. Snakes possess large numbers of the olfactory receptor (OR) genes and the type-2 vomeronasal receptor (V2R) genes, and expression profiling in the olfactory tissues suggests that snakes use the ORs in the main olfactory system (MOS) and the V2Rs in the vomeronasal system (VNS). The number of OR genes has decreased in sea snakes, and fully aquatic species lost MOS which is responsible for detecting airborne odours. By contrast, sea snakes including fully aquatic species retain a number of V2R genes and a well-developed VNS for smelling underwater. This study suggests that the sense of smell also degenerated in sea snakes, particularly in fully aquatic species, but their residual olfactory capability is distinct from that of other fully aquatic amniotes. Amphibious species show an intermediate status between terrestrial and fully aquatic snakes, implying their importance in understanding the process of aquatic adaptation.

Comprehensive analysis of genes contributing to euryhalinity in the bull shark, Carcharhinus leucas; Na+-Cl- co-transporter is one of the key renal factors upregulated in acclimation to low-salinity environment

Most cartilaginous fishes live principally in seawater (SW) environments, but a limited number of species including the bull shark, Carcharhinus leucas, inhabit both SW and freshwater (FW) environments during their life cycle. Euryhaline elasmobranchs maintain high internal urea and ion levels even in FW environments, but little is known about the osmoregulatory mechanisms that enable them to maintain internal homeostasis in hypoosmotic environments. In the present study, we focused on the kidney because this is the only organ that can excrete excess water from the body in a hypoosmotic environment. We conducted a transfer experiment of bull sharks from SW to FW and performed differential gene expression analysis between the two conditions using RNA-sequencing. A search for genes upregulated in the FW-acclimated bull shark kidney indicated that the expression of the Na⁺-Cl^- cotransporter (NCC; Slc12a3) was 10 times higher in the FW-acclimated sharks compared with that in SW sharks. In the kidney, apically located NCC was observed in the late distal tubule and in the anterior half of the collecting tubule, where basolateral Na⁺/K⁺-ATPase was also expressed, implying that these segments contribute to NaCl reabsorption from the filtrate for diluting the urine. This expression pattern was not observed in the houndshark, Triakis scyllium, which had been transferred to 30% SW; this species cannot survive in FW environments. The salinity transfer experiment combined with a comprehensive gene screening approach demonstrates that NCC is a key renal protein that contributes to the remarkable euryhaline ability of the bull shark.

Evolution of nodal and nodal-related genes and the putative composition of the heterodimers that trigger the nodal pathway in vertebrates

Nodal is a signaling molecule that belongs to the transforming growth factor-β superfamily that plays key roles during the early stages of development of animals. In vertebrates Nodal forms an heterodimer with a GDF1/3 protein to activate the Nodal pathway. Vertebrates have a paralog of nodal in their genomes labeled Nodal-related, but the evolutionary history of these genes is a matter of debate, mainly because of the presence of a variable numbers of genes in the vertebrate genomes sequenced so far. Thus, the goal of this study was to investigate the evolutionary history of the Nodal and Nodal-related genes with an emphasis in tracking changes in the number of genes among vertebrates. Our results show the presence of two gene lineages (Nodal and Nodal-related) that can be traced back to the ancestor of jawed vertebrates. These lineages have undergone processes of differential retention and lineage-specific expansions. Our results imply that Nodal and Nodal-related duplicated at the latest in the ancestor of gnathostomes, and they still retain a significant level of functional redundancy. By comparing the evolution of the Nodal/Nodal-related with GDF1/3 gene family, it is possible to infer that there are several types of heterodimers that can trigger the Nodal pathway among vertebrates.

Inference of the ancestral vertebrate phenotype through vestiges of the whole-genome duplications

Inferring the phenotype of the last common ancestor of living vertebrates is a challenging problem because of several unresolvable factors. They include the lack of reliable out-groups of living vertebrates, poor information about less fossilizable organs and specialized traits of phylogenetically important species, such as lampreys and hagfishes (e.g. secondary loss of vertebrae in adult hagfishes). These factors undermine the reliability of ancestral reconstruction by traditional character mapping approaches based on maximum parsimony. In this article, we formulate an approach to hypothesizing ancestral vertebrate phenotypes using information from the phylogenetic and functional properties of genes duplicated by genome expansions in early vertebrate evolution. We named the conjecture as ‘chronological reconstruction of ohnolog functions (CHROF)’. This CHROF conjecture raises the possibility that the last common ancestor of living vertebrates may have had more complex traits than currently thought.

Co-option of the PRDM14–CBFA2T complex from motor neurons to pluripotent cells during vertebrate evolution

Gene regulatory networks underlying cellular pluripotency are controlled by a core circuitry of transcription factors in mammals, including POU5F1. However, the evolutionary origin and transformation of pluripotency-related transcriptional networks have not been elucidated in deuterostomes. PR domain-containing protein 14 (PRDM14) is specifically expressed in pluripotent cells and germ cells, and required for establishing embryonic stem cells (ESCs) and primordial germ cells in mice. Here, we compared the functions and expression patterns of PRDM14 orthologues within deuterostomes. Amphioxus PRDM14 and zebrafish PRDM14, but not sea urchin PRDM14, compensated for mouse PRDM14 function in maintaining mouse ESC pluripotency. Interestingly, sea urchin PRDM14 together with sea urchin CBFA2T, an essential partner of PRDM14 in mouse ESCs, complemented the self-renewal defect in mouse Prdm14 KO ESCs. Contrary to the Prdm14-expression pattern in mouse embryos, Prdm14 was expressed in motor neurons of amphioxus embryos as observed in zebrafish embryos. Thus, Prdm14 expression in motor neurons was conserved in non-tetrapod deuterostomes and the co-option of the PRDM14-CBFA2T complex from motor neurons into pluripotent cells may have maintained the transcriptional network for pluripotency during vertebrate evolution.

Embryonic transcriptome sequencing of the ocellate spot skate Okamejei kenojei

Chondrichthyans (cartilaginous fishes) exhibit highly variable reproductive styles, categorized as viviparity and oviparity. Among these, species with oviparity provide an enormous potential of molecular experimentation with stable sample supply which does not demand the sacrifices of live mothers. Cartilaginous fishes are divided into two subclasses, chimaeras (Holocephali) and elasmobranchs (Elasmobranchii), and the latter consists of two monophyletic groups, Batoidea (rays, skates and torpedoes) and Selachimorpha (sharks). Here we report transcriptome assemblies of the ocellate spot skate Okamejei kenojei, produced by strand-specific RNA-seq of its embryonic tissues. We obtained a total of 325 million illumina short reads from libraries prepared using four different tissue domains and assembled them all together. Our assembly result confirmed the species authenticity and high continuity of contig sequences. Also, assessment of its coverage of pre-selected one-to-one orthologs supported high diversity of transcripts in the assemblies. Our products are expected to provide a basis of comparative molecular studies encompassing other chondrichthyan species with emerging genomic and transcriptomic sequence information.

A de novo transcriptome assembly of the zebra bullhead shark, Heterodontus zebra

Although cartilaginous fishes have played crucial roles in various fields, including evolutionary biology, marine ecology, bioresources, and aquarium exhibitions, molecular information for these species is poorly available. The present study reports a transcriptome assembly from an embryo of the zebra bullhead shark (Heterodontus zebra), produced by paired-end RNA sequencing. Transcriptome data is generated with a de novo transcriptome assembler, Trinity. Amino acid sequences are predicted from the assemblies, using TransDecoder. Because cartilaginous fishes serve as the outgroup of bony vertebrates, the data would contribute to comparative analyses of a various biological fields. In addition, this study would be useful for conservation biology, such as transcriptome-based population genetics.

gVolante for standardizing completeness assessment of genome and transcriptome assemblies

Motivation

Along with the increasing accessibility to comprehensive sequence information, such as whole genomes and transcriptomes, the demand for assessing their quality has been multiplied. To this end, metrics based on sequence lengths, such as N50, have become a standard, but they only evaluate one aspect of assembly quality. Conversely, analyzing the coverage of pre-selected reference protein-coding genes provides essential content-based quality assessment, but the currently available pipelines for this purpose, CEGMA and BUSCO, do not have a user-friendly interface to serve as a uniform environment for assembly completeness assessment.

Results

Here, we introduce a brand-new web server, gVolante, which provides an online tool for (i) on-demand completeness assessment of sequence sets by means of the previously developed pipelines CEGMA and BUSCO and (ii) browsing pre-computed completeness scores for publicly available data in its database section. Completeness assessments performed on gVolante report scores based on not just the coverage of reference genes but also on sequence lengths (e.g. N50 scaffold length), allowing quality control in multiple aspects. Using gVolante, one can compare the quality of original assemblies between their multiple versions (obtained through program choice and parameter tweaking, for example) and evaluate them in comparison to the scores of public resources found in the database section.

Availability and implementation

gVoalte is freely available at https://gvolante.riken.jp/.

Loss of Fam60a, a Sin3a subunit, results in embryonic lethality and is associated with aberrant methylation at a subset of gene promoters

We have examined the role of Fam60a, a gene highly expressed in embryonic stem cells, in mouse development. Fam60a interacts with components of the Sin3a-Hdac transcriptional corepressor complex, and most Fam60a^-/- embryos manifest hypoplasia of visceral organs and die in utero. Fam60a is recruited to the promoter regions of a subset of genes, with the expression of these genes being either up- or down-regulated in Fam60a^-/- embryos. The DNA methylation level of the Fam60a target gene Adhfe1 is maintained at embryonic day (E) 7.5 but markedly reduced at E9.5 in Fam60a^-/- embryos, suggesting that DNA demethylation is enhanced in the mutant. Examination of genome-wide DNA methylation identified several differentially methylated regions, which were preferentially hypomethylated, in Fam60a^-/- embryos. Our data suggest that Fam60a is required for proper embryogenesis, at least in part as a result of its regulation of DNA methylation at specific gene promoters.

A staging table for the embryonic development of the brownbanded bamboo shark (Chiloscyllium punctatum)

BACKGROUND

Studying cartilaginous fishes (chondrichthyans) has helped us understand vertebrate evolution and diversity. However, resources such as genome sequences, embryos, and detailed staging tables are limited for species within this clade. To overcome these limitations, we have focused on a species, the brownbanded bamboo shark (Chiloscyllium punctatum), which is a relatively common aquarium species that lays eggs continuously throughout the year. In addition, because of its relatively small genome size, this species is promising for molecular studies.

RESULTS

To enhance biological studies of cartilaginous fishes, we establish a normal staging table for the embryonic development of the brownbanded bamboo shark. Bamboo shark embryos take around 118 days to reach the hatching period at 25°C, which is approximately 1.5 times as fast as the small-spotted catshark (Scyliorhinus canicula) takes. Our staging table divides the embryonic period into 38 stages. Furthermore, we found culture conditions that allow early embryos to grow in partially opened egg cases.

CONCLUSIONS

In addition to the embryonic staging table, we show that bamboo shark embryos exhibit relatively fast embryonic growth and are amenable to culture, key characteristics that enhance their experimental utility. Therefore, the present study is a foundation for cartilaginous fish research. Developmental Dynamics 247:712-723, 2018. © 2017 Wiley Periodicals, Inc.

Molecular cytogenetic characterization of chromosome site-specific repetitive sequences in the Arctic lamprey (Lethenteron camtschaticum, Petromyzontidae)

All extant lamprey karyotypes are characterized by almost all dot-shaped microchromosomes. To understand the molecular basis of chromosome structure in lampreys, we performed chromosome C-banding and silver staining and chromosome mapping of the 18S–28S and 5S ribosomal RNA (rRNA) genes and telomeric TTAGGG repeats in the Arctic lamprey (Lethenteron camtschaticum). In addition, we cloned chromosome site-specific repetitive DNA sequences and characterized them by nucleotide sequencing, chromosome in situ hybridization, and filter hybridization. Three types of repetitive sequences were detected; a 200-bp AT-rich repetitive sequence, LCA-EcoRIa that co-localized with the 18S–28S rRNA gene clusters of 3 chromosomal pairs; a 364-bp AT-rich LCA-EcoRIb sequence that showed homology to the EcoRI sequence family from the sea lamprey (Petromyzon marinus), which contains short repeats as centromeric motifs; and a GC-rich 702-bp LCA-ApaI sequence that was distributed on nearly all chromosomes and showed significant homology with the integrase-coding region of a Ty3/Gypsy family long terminal repeat (LTR) retrotransposon. All three repetitive sequences are highly conserved within the Petromyzontidae or within Petromyzontidae and Mordaciidae. Molecular cytogenetic characterization of these site-specific repeats showed that they may be correlated with programed genome rearrangement (LCA-EcoRIa), centromere structure and function (LCA-EcoRIb), and site-specific amplification of LTR retroelements through homogenization between non-homologous chromosomes (LCA-ApaI).

Correction to 'Evolutionary diversification of retinoic acid receptor ligand-binding pocket structure by molecular tinkering'

[This corrects the article DOI: 10.1098/rsos.150484.].

The evolutionally-conserved function of group B1 Sox family members confers the unique role of Sox2 in mouse ES cells

BACKGROUND

In mouse ES cells, the function of Sox2 is essential for the maintenance of pluripotency. Since the Sox-family of transcription factors are well conserved in the animal kingdom, addressing the evolutionary origin of Sox2 function in pluripotent stem cells is intriguing from the perspective of understanding the origin of pluripotency.

RESULTS

Here we approach this question using a functional complementation assay in inducible Sox2-null ES cells. Assaying mouse Sox proteins from different Groups, we found that only Group B1 and Group G proteins were able to support pluripotency. Interestingly, invertebrate homologs of mammalian Group B1 Sox proteins were able to replace the pluripotency-associated function of mouse Sox2. Moreover, the mouse ES cells rescued by the Drosophila SoxNeuro protein are able to contribute to chimeric embryos.

CONCLUSIONS

These data indicate that the function of mouse Sox2 supporting pluripotency is based on an evolutionally conserved activity of the Group B1 Sox family. Since pluripotent stem cell population in developmental process could be regarded as the evolutional novelty in vertebrates, it could be regarded as a co-optional use of their evolutionally conserved function.