Giant lungfish genome elucidates the conquest of land by vertebrates.

Lungfishes belong to lobe-fined fish (Sarcopterygii) that, in the Devonian period, 'conquered' the land and ultimately gave rise to all land vertebrates, including humans1-3. Here we determine the chromosome-quality genome of the Australian lungfish (Neoceratodus forsteri), which is known to have the largest genome of any animal. The vast size of this genome, which is about 14× larger than that of humans, is attributable mostly to huge intergenic regions and introns with high repeat content (around 90%), the components of which resemble those of tetrapods (comprising mainly long interspersed nuclear elements) more than they do those of ray-finned fish. The lungfish genome continues to expand independently (its transposable elements are still active), through mechanisms different to those of the enormous genomes of salamanders. The 17 fully assembled lungfish macrochromosomes maintain synteny to other vertebrate chromosomes, and all microchromosomes maintain conserved ancient homology with the ancestral vertebrate karyotype. Our phylogenomic analyses confirm previous reports that lungfish occupy a key evolutionary position as the closest living relatives to tetrapods4,5, underscoring the importance of lungfish for understanding innovations associated with terrestrialization. Lungfish preadaptations to living on land include the gain of limb-like expression in developmental genes such as hoxc13 and sall1 in their lobed fins. Increased rates of evolution and the duplication of genes associated with obligate air-breathing, such as lung surfactants and the expansion of odorant receptor gene families (which encode proteins involved in detecting airborne odours), contribute to the tetrapod-like biology of lungfishes. These findings advance our understanding of this major transition during vertebrate evolution.

Chromatin organization and behavior in HRAS-transformed mouse fibroblasts.

In higher eukaryotic cells, a string of nucleosomes, where long genomic DNA is wrapped around core histones, are rather irregularly folded into a number of condensed chromatin domains, which have been revealed by super-resolution imaging and Hi-C technologies. Inside these domains, nucleosomes fluctuate and locally behave like a liquid. The behavior of chromatin may be highly related to DNA transaction activities such as transcription and repair, which are often upregulated in cancer cells. To investigate chromatin behavior in cancer cells and compare those of cancer and non-cancer cells, we focused on oncogenic-HRAS (Gly12Val)-transformed mouse fibroblasts CIRAS-3 cells and their parental 10T1/2 cells. CIRAS-3 cells are tumorigenic and highly metastatic. First, we found that HRAS-induced transformation altered not only chromosome structure, but also nuclear morphology in the cell. Using single-nucleosome imaging/tracking in live cells, we demonstrated that nucleosomes are locally more constrained in CIRAS-3 cells than in 10T1/2 cells. Consistently, heterochromatin marked with H3K27me3 was upregulated in CIRAS-3 cells. Finally, Hi-C analysis showed enriched interactions of the B-B compartment in CIRAS-3 cells, which likely represents transcriptionally inactive chromatin. Increased heterochromatin may play an important role in cell migration, as they have been reported to increase during metastasis. Our study also suggests that single-nucleosome imaging provides new insights into how local chromatin is structured in living cells.

The shh limb enhancer is activated in patterned limb regeneration but not in hypomorphic limb regeneration in Xenopus laevis.

Xenopus young tadpoles regenerate a limb with the anteroposterior (AP) pattern, but metamorphosed froglets regenerate a hypomorphic limb after amputation. The key gene for AP patterning, shh, is expressed in a regenerating limb of the tadpole but not in that of the froglet. Genomic DNA in the shh limb-specific enhancer, MFCS1 (ZRS), is hypermethylated in froglets but hypomethylated in tadpoles: shh expression may be controlled by epigenetic regulation of MFCS1. Is MFCS1 specifically activated for regenerating the AP-patterned limb? We generated transgenic Xenopus laevis lines that visualize the MFCS1 enhancer activity with a GFP reporter. The transgenic tadpoles showed GFP expression in hoxd13-and shh-expressing domains of developing and regenerating limbs, whereas the froglets showed no GFP expression in the regenerating limbs despite having hoxd13 expression. Genome sequence analysis and co-transfection assays using cultured cells revealed that Hoxd13 can activate Xenopus MFCS1. These results suggest that MFCS1 activation correlates with regeneration of AP-patterned limbs and that re-activation of epigenetically inactivated MFCS1 would be crucial to confer the ability to non-regenerative animals for regenerating a properly patterned limb.

Development of a heat-stable alkaline phosphatase reporter system for cis-regulatory analysis and its application to 3D digital imaging of Xenopus embryonic tissues.

Xenopus is one of the essential model systems for studying vertebrate development. However, one drawback of this system is that, because of the opacity of Xenopus embryos, 3D imaging analysis is limited to surface structures, explant cultures, and post-embryonic tadpoles. To develop a technique for 3D tissue/organ imaging in whole Xenopus embryos, we identified optimal conditions for using placental alkaline phosphatase (PLAP) as a transgenic reporter and applied it to the correlative light microscopy and block-face imaging (CoMBI) method for visualization of PLAP-expressing tissues/organs. In embryos whose endogenous alkaline phosphatase activities were heat-inactivated, PLAP staining visualized various tissue-specific enhancer/promoter activities in a manner consistent with green fluorescent protein (GFP) fluorescence. Furthermore, PLAP staining appeared to be more sensitive than GFP fluorescence as a reporter, and the resulting expression patterns were not mosaic, in striking contrast to the mosaic staining pattern of ß-galactosidase expressed from the lacZ gene that was introduced by the same transgenesis method. Owing to efficient penetration of alkaline phosphatase substrates, PLAP activity was detected in deep tissues, such as the developing brain, spinal cord, heart, and somites, by whole-mount staining. The stained embryos were analyzed by the CoMBI method, resulting in the digital reconstruction of 3D images of the PLAP-expressing tissues. These results demonstrate the efficacy of the PLAP reporter system for detecting enhancer/promoter activities driving deep tissue expression and its combination with the CoMBI method as a powerful approach for 3D digital imaging analysis of specific tissue/organ structures in Xenopus embryos.

The giant axolotl genome uncovers the evolution, scaling, and transcriptional control of complex gene loci.

Vertebrates harbor recognizably orthologous gene complements but vary 100-fold in genome size. How chromosomal organization scales with genome expansion is unclear, and how acute changes in gene regulation, as during axolotl limb regeneration, occur in the context of a vast genome has remained a riddle. Here, we describe the chromosome-scale assembly of the giant, 32 Gb axolotl genome. Hi-C contact data revealed the scaling properties of interphase and mitotic chromosome organization. Analysis of the assembly yielded understanding of the evolution of large, syntenic multigene clusters, including the Major Histocompatibility Complex (MHC) and the functional regulatory landscape of the Fibroblast Growth Factor 8 (Axfgf8) region. The axolotl serves as a primary model for studying successful regeneration.

Asymmetrically reduced expression of hand1 homeologs involving a single nucleotide substitution in a cis-regulatory element.

During vertebrate evolution, whole genome duplications resulted in a number of duplicated genes, some of which eventually changed their expression patterns and/or levels via alteration of cis-regulatory sequences. However, the initial process involved in such cis-regulatory changes remains unclear. Therefore, we investigated this process by analyzing the duplicated hand1 genes of Xenopus laevis (hand1.L and hand1.S), which were generated by allotetraploidization 17-18 million years ago, and compared these with their single ortholog in the ancestral-type diploid species X. tropicalis. A dN/dS analysis indicated that hand1.L and hand1.S are still under purifying selection, and thus, their products appear to retain ancestral functional properties. RNA-seq and in situ hybridization analyses revealed that hand1.L and hand1.S have similar expression patterns to each other and to X. tropicalis hand1, but the hand1.S expression level was much lower than the hand1.L expression level in the primordial heart. A comparative sequence analysis, luciferase reporter analysis, ChIP-PCR analysis, and transgenic reporter analysis showed that a single nucleotide substitution in the hand1.S promoter was responsible for the reduced expression in the heart. These findings demonstrated that a small change in the promoter sequence can trigger diversification of duplicated gene expression prior to diversification of their encoded protein functions in a young duplicated genome.

Co-accumulation of cis-regulatory and coding mutations during the pseudogenization of the Xenopus laevis homoeologs six6.L and six6.S.

Common models for the evolution of duplicated genes after genome duplication are subfunctionalization, neofunctionalization, and pseudogenization. Although the crucial roles of cis-regulatory mutations in subfunctionalization are well-documented, their involvement in pseudogenization and/or neofunctionalization remains unclear. We addressed this issue by investigating the evolution of duplicated homeobox genes, six6.L and six6.S, in the allotetraploid frog Xenopus laevis. Based on a comparative expression analysis, we observed similar eye-specific expression patterns for the two loci and their single ortholog in the ancestral-type diploid species Xenopus tropicalis. However, we detected lower levels of six6.S expression than six6.L expression. The six6.S enhancer sequence was more highly diverged from the orthologous enhancer of X. tropicalis than the six6.L enhancer, and showed weaker activity in a transgenic reporter assay. Based on a phylogenetic analysis of the protein sequences, we observed greater divergence between X. tropicalis Six6 and Six6.S than between X. tropicalis Six6 and Six6.L, and the observed mutations were reminiscent of a microphthalmia mutation in human SIX6. Misexpression experiments showed that six6.S has weaker eye-enlarging activity than six6.L, and targeted disruption of six6.L reduced the eye size more significantly than that of six6.S. These results suggest that enhancer attenuation stimulates the accumulation of hypomorphic coding mutations, or vice versa, in one duplicated gene copy and facilitates pseudogenization. We also underscore the value of the allotetraploid genome of X. laevis as a resource for studying latent pathogenic mutations.

Noncanonical Hox, Etv4, and Gli3 gene activities give insight into unique limb patterning in salamanders.

Limb development in salamanders is unique among tetrapods in significant ways. Not only can salamanders regenerate lost limbs repeatedly and throughout their lives, but also the preaxial zeugopodial element and digits form before the postaxial ones and, hence, with a reversed polarity compared to all other tetrapods. Moreover, in salamanders with free-swimming larval stages, as exemplified by the axolotl (Ambystoma mexicanum), each digit buds independently, instead of undergoing a paddle stage. Here, we report gene expression patterns of Hoxa and d clusters, and other crucial transcription factors during axolotl limb development. During early phases of limb development, expression patterns are mostly similar to those reported for amniotes and frogs. Likewise, Hoxd and Shh regulatory landscapes are largely conserved. However, during late digit-budding phases, remarkable differences are present: (i) the Hoxd13 expression domain excludes developing digits I and IV, (ii) we expand upon previous observation that Hoxa11 expression, which traditionally marks the zeugopodium, extends distally into the developing digits, and (iii) Gli3 and Etv4 show prolonged expression in developing digits. Our findings identify derived patterns in the expression of key transcription factors during late phases of salamander limb development, and provide the basis for a better understanding of the unique patterning of salamander limbs.

Epigenetic modification maintains intrinsic limb-cell identity in Xenopus limb bud regeneration.

Many amphibians can regenerate limbs, even in adulthood. If a limb is amputated, the stump generates a blastema that makes a complete, new limb in a process similar to developmental morphogenesis. The blastema is thought to inherit its limb-patterning properties from cells in the stump, and it retains the information despite changes in morphology, gene expression, and differentiation states required by limb regeneration. We hypothesized that these cellular properties are maintained as epigenetic memory through histone modifications. To test this hypothesis, we analyzed genome-wide histone modifications in Xenopus limb bud regeneration. The trimethylation of histone H3 at lysine 4 (H3K4me3) is closely related to an open chromatin structure that allows transcription factors access to genes, whereas the trimethylation of histone H3 at lysine 27 (H3K27me3) is related to a closed chromatin state that blocks the access of transcription factors. We compared these two modification profiles by high-throughput sequencing of samples prepared from the intact limb bud and the regenerative blastema by chromatin immunoprecipitation. For many developmental genes, histone modifications at the transcription start site were the same in the limb bud and the blastema, were stable during regeneration, and corresponded well to limb properties. These results support our hypothesis that histone modifications function as a heritable cellular memory to maintain limb cell properties, despite dynamic changes in gene expression during limb bud regeneration in Xenopus.

Application of the cis-regulatory region of a heat-shock protein 70 gene to heat-inducible gene expression in the ascidian Ciona intestinalis.

Temporally controlled induction of gene expression is a useful technique for analyzing gene function. To make such a technique possible in Ciona intestinalis embryos, we employed the cis-regulatory region of the heat-shock protein 70 (HSP70) gene Ci-HSPA1/6/7-like for heat-inducible gene expression in C. intestinalis embryos. We showed that Ci-HSPA1/6/7-like becomes heat shock-inducible by the 32-cell stage during embryogenesis. The 5'-upstream region of Ci-HSPA1/6/7-like, which contains heat-shock elements indispensable for heat-inducible gene expression, induces the heat shock-dependent expression of a reporter gene in the whole embryo from the 32-cell to the middle gastrula stages and in progressively restricted areas of embryos in subsequent stages. We assessed the effects of heat-shock treatments in different conditions on the normality of embryos and induction of transgene expression. We evaluated the usefulness of this technique through overexpression experiments on the well-characterized, developmentally relevant gene, Ci-Bra, and showed that this technique is applicable for inferring the gene function in C. intestinalis.

A chromatin code for limb segment identity in axolotl limb regeneration.

The salamander limb correctly regenerates missing limb segments because connective tissue cells have segment-specific identities, termed "positional information". How positional information is molecularly encoded at the chromatin level has been unknown. Here, we performed genome-wide chromatin profiling in mature and regenerating axolotl limb connective tissue cells. We find segment-specific levels of histone H3K27me3 as the major positional mark, especially at limb homeoprotein gene loci but not their upstream regulators, constituting an intrinsic segment information code. During regeneration, regeneration-specific regulatory elements became active prior to the re-appearance of developmental regulatory elements. In the hand, the permissive chromatin state of the homeoprotein gene HoxA13 engages with the regeneration program bypassing the upper limb program. Comparison of regeneration regulatory elements with those found in other regenerative animals identified a core shared set of transcription factors, supporting an ancient, conserved regeneration program.

Chromosome Conformation Capture for Large Genomes.

The gigantic 32Gb Axolotl genome inspires fascinating questions such as: how such a big genome is organized and packed in nuclei and how regulation of gene transcription can happen over such large genomic distances. Currently, there are many technical challenges when we investigate chromatin architecture in axolotl. For example, probing promoter-enhancer interactions in such a large genome. Chromatin capture methods (e.g., Chromatin Conformation Capture) have been used in a variety of species. The large size of the axolotl nuclei and its genome requires the adaptation of such methods. Here, we describe a detailed protocol for high-throughput genome-wide conformation capture (Hi-C) using axolotl limb cells. This Hi-C library preparation protocol can also be used to prepare libraries from other nonmodel organisms such as Lungfish and Cephalopods. We believe that our protocol could be useful for a variety of animal systems including other salamanders.

Emergence of novel cephalopod gene regulation and expression through large-scale genome reorganization.

Coleoid cephalopods (squid, cuttlefish, octopus) have the largest nervous system among invertebrates that together with many lineage-specific morphological traits enables complex behaviors. The genomic basis underlying these innovations remains unknown. Using comparative and functional genomics in the model squid Euprymna scolopes, we reveal the unique genomic, topological, and regulatory organization of cephalopod genomes. We show that coleoid cephalopod genomes have been extensively restructured compared to other animals, leading to the emergence of hundreds of tightly linked and evolutionary unique gene clusters (microsyntenies). Such novel microsyntenies correspond to topological compartments with a distinct regulatory structure and contribute to complex expression patterns. In particular, we identify a set of microsyntenies associated with cephalopod innovations (MACIs) broadly enriched in cephalopod nervous system expression. We posit that the emergence of MACIs was instrumental to cephalopod nervous system evolution and propose that microsyntenic profiling will be central to understanding cephalopod innovations.

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.

Taming chlorophylls by early eukaryotes underpinned algal interactions and the diversification of the eukaryotes on the oxygenated Earth.

Extant eukaryote ecology is primarily sustained by oxygenic photosynthesis, in which chlorophylls play essential roles. The exceptional photosensitivity of chlorophylls allows them to harvest solar energy for photosynthesis, but on the other hand, they also generate cytotoxic reactive oxygen species. A risk of such phototoxicity of the chlorophyll must become particularly prominent upon dynamic cellular interactions that potentially disrupt the mechanisms that are designed to quench photoexcited chlorophylls in the phototrophic cells. Extensive examination of a wide variety of phagotrophic, parasitic, and phototrophic microeukaryotes demonstrates that a catabolic process that converts chlorophylls into nonphotosensitive 132,173-cyclopheophorbide enols (CPEs) is phylogenetically ubiquitous among extant eukaryotes. The accumulation of CPEs is identified in phagotrophic algivores belonging to virtually all major eukaryotic assemblages with the exception of Archaeplastida, in which no algivorous species have been reported. In addition, accumulation of CPEs is revealed to be common among phototrophic microeukaryotes (i.e., microalgae) along with dismantling of their secondary chloroplasts. Thus, we infer that CPE-accumulating chlorophyll catabolism (CACC) primarily evolved among algivorous microeukaryotes to detoxify chlorophylls in an early stage of their evolution. Subsequently, it also underpinned photosynthetic endosymbiosis by securing close interactions with photosynthetic machinery containing abundant chlorophylls, which led to the acquisition of secondary chloroplasts. Our results strongly suggest that CACC, which allowed the consumption of oxygenic primary producers, ultimately permitted the successful radiation of the eukaryotes throughout and after the late Proterozoic global oxygenation.

Draft genome of Dugesia japonica provides insights into conserved regulatory elements of the brain restriction gene nou-darake in planarians.

BACKGROUND: Planarians are non-parasitic Platyhelminthes (flatworms) famous for their regeneration ability and for having a well-organized brain. Dugesia japonica is a typical planarian species that is widely distributed in the East Asia. Extensive cellular and molecular experimental methods have been developed to identify the functions of thousands of genes in this species, making this planarian a good experimental model for regeneration biology and neurobiology. However, no genome-level information is available for D. japonica, and few gene regulatory networks have been identified thus far. RESULTS: To obtain whole-genome information on this species and to study its gene regulatory networks, we extracted genomic DNA from 200 planarians derived from a laboratory-bred asexual clonal strain, and sequenced 476 Gb of data by second-generation sequencing. Kmer frequency graphing and fosmid sequence analysis indicated a complex genome that would be difficult to assemble using second-generation sequencing short reads. To address this challenge, we developed a new assembly strategy and improved the de novo genome assembly, producing a 1.56 Gb genome sequence (DjGenome ver1.0, including 202,925 scaffolds and N50 length 27,741 bp) that covers 99.4% of all 19,543 genes in the assembled transcriptome, although the genome is fragmented as 80% of the genome consists of repeated sequences (genomic frequency ≥ 2). By genome comparison between two planarian genera, we identified conserved non-coding elements (CNEs), which are indicative of gene regulatory elements. Transgenic experiments using Xenopus laevis indicated that one of the CNEs in the Djndk gene may be a regulatory element, suggesting that the regulation of the ndk gene and the brain formation mechanism may be conserved between vertebrates and invertebrates. CONCLUSION: This draft genome and CNE analysis will contribute to resolving gene regulatory networks in planarians. The genome database is available at: http://www.planarian.jp.

Live Imaging of Axolotl Digit Regeneration Reveals Spatiotemporal Choreography of Diverse Connective Tissue Progenitor Pools.

Connective tissues-skeleton, dermis, pericytes, fascia-are a key cell source for regenerating the patterned skeleton during axolotl appendage regeneration. This complexity has made it difficult to identify the cells that regenerate skeletal tissue. Inability to identify these cells has impeded a mechanistic understanding of blastema formation. By tracing cells during digit tip regeneration using brainbow transgenic axolotls, we show that cells from each connective tissue compartment have distinct spatial and temporal profiles of proliferation, migration, and differentiation. Chondrocytes proliferate but do not migrate into the regenerate. In contrast, pericytes proliferate, then migrate into the blastema and give rise solely to pericytes. Periskeletal cells and fibroblasts contribute the bulk of digit blastema cells and acquire diverse fates according to successive waves of migration that choreograph their proximal-distal and tissue contributions. We further show that platelet-derived growth factor signaling is a potent inducer of fibroblast migration, which is required to form the blastema.

Comparative expression analysis of the H3K27 demethylases, JMJD3 and UTX, with the H3K27 methylase, EZH2, in Xenopus.

The regulated removal of the gene-silencing epigenetic mark, trimethylation of lysine 27 of histone H3 (H3K27me3), has been shown to be critical for tissue-specific activation of developmental genes; however, the extent of embryonic expression of its demethylases, JMJD3 and UTX, has remained unclear. In this study, we investigated the expression of jmjd3 and utx genes in Xenopus embryos in parallel with that of the H3K27 methylase gene, ezh2. At the blastula stage, jmjd3, utx and ezh2 showed similar expression patterns in the animal cap and marginal zone that give rise to the ectoderm and mesoderm, respectively. The three genes maintained similar expression patterns in the neural plate, preplacodal ectoderm and axial mesoderm during the gastrula and neurula stages. Later, expression was maintained in the developing brain and cranial sensory tissues, such as the eye and ear, of tailbud embryos. These findings suggest that the H3K27 demethylases and methylase may function continuously for progressive switching of genetic programs during neural development, a model involving the simultaneous action of both of the demethylases for the de-repression of silent genes and the methylase for the silencing of active genes.

Evolution of a tissue-specific silencer underlies divergence in the expression of pax2 and pax8 paralogues.

Recent studies underscore a role for the differential degeneration of enhancers in the evolutionary diversification of paralogue expression. However, no one has reported evidence for the involvement of innovative cis-regulatory changes. Here we show that silencer innovation diversified expression of the vertebrate paralogues, pax2 and pax8. pax2 shows multi-tissue expression, as does the ancestral amphioxus orthologue, pax2/5/8, whereas pax8 expression localizes to a subset of pax2-expressing tissues. We reveal that both pax2 and pax8 retain ancestral enhancers capable of directing pax2-like, multi-tissue expression. However, a silencer within the pax8 proximal promoter suppresses pleiotropic enhancer activity outside the pax8-expressing tissues. In contrast, the combination of the pax2 proximal promoter with either the pax8 or pax2 enhancer recapitulates pax2-like expression, as in the amphioxus pax2/5/8 promoter. We propose that silencer innovation, rather than enhancer degeneration, was crucial for the divergent expression of paralogues with pleiotropic enhancers inherited from their common progenitor.