Search

The salamander limb regenerates only the missing portion. Each limb segment can only form segments equivalent to- or more distal to their own identity, relying on a property termed “positional information”. How positional information is encoded in limb cells has been unknown. By cell-type-specific chromatin profiling of upper arm, lower arm, and hand, we found segment-specific levels of histone H3K27me3 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. This means that, in the hand segment, the permissive chromatin state of the hand homeoprotein geneHoxA13engages with regeneration regulatory elements, bypassing the upper limb program.

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.

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 13(2),17(3)-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.

Significance The axolotl is an important model organism because it is a tetrapod with a similar body plan to humans. Unlike humans, the axolotl regenerates limbs and other complex tissues. Therefore, the axolotl contributes to understanding evolution, development, and regeneration. With sophisticated tools for gene modification and tissue labeling, a fully assembled genome sequence was a sorely missing resource. Assembly was difficult because the genome size is 10× that of humans. Here, we use a cross-linking strategy called Hi-C to link together fragmented genome sequences to chromosome scale. We show that gene regulation occurs over very large genomic distances and that mitotic chromosomes are packaged efficiently., 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.

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., A chromosome-quality genome of the lungfish Neoceratodus fosteri sheds light on the development of obligate air-breathing and the gain of limb-like gene expression in lobed fins, providing insights into the water-to-land transition in vertebrate evolution.

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.

Summary 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., Highlights • Distinct migratory and proliferative dynamics in axolotl connective tissue subtypes • Timing of dermal migration into the blastema biases contribution to skeleton or dermis • Participating source zone of 50–500 μm, depending on cell type and amputation site • PDGF-BB is a connective tissue pro-migratory signal necessary for blastema formation, Currie et al. use live imaging of brainbow transgenic axolotls to identify which connective tissues build the blastema for regeneration and the source zone from which they migrate. Each connective tissue cell type has highly characteristic migration kinetics that choreograph their fate and tissue contribution during regeneration.

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. Electronic supplementary material The online version of this article (10.1186/s40851-018-0102-2) contains supplementary material, which is available to authorized users.

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. 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. 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 .

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.

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.

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.

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.

Recent accumulation of genome sequences indicates that vertebrate genomes retained many duplicate genes, called paralogues, after whole-genome duplications (WGDs) occurring in the early stages of chordate evolution. Although paralogous gene pairs probably had the same expression patterns immediately after WGDs, they often showed distinct expression patterns in modern vertebrates. Nucleotide sequence changes in noncoding genomic regions, including cis-regulatory elements (CRE) such as enhancers and silencers, are the primary cause of the divergent paralogue expression. Such CRE evolution can be classified as follows: (1) enhancer innovation, (2) silencer innovation, and (3) ancestral CRE degeneration. In this chapter, we first show examples of the degeneration of ancestral CREs, with reference to the duplication–degeneration–complementation model, and subsequently introduce recently discovered evidence for silencer innovation and conservation of ancestral pleiotropic enhancers. Finally, we discuss the possible involvement of enhancer innovation in the divergent paralogue expression.

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.

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.

Recent studies have revealed that epigenetic modifications, such as histone methylation, histone acetylation and DNA methylation, play fundamental roles in gene expression from the genome. Among them, the trimethylation/demethylation of histone H3 lysine 27 (H3K27) is associated with silencing/activation of some developmental regulatory genes, such as Hox clusters. In mammals, the trimethylation of H3K27 is mediated by the Polycomb repressive complex 2 (PRC2) that consists of three core subunits, Ezh2, Eed and Suz12 proteins, whereas the demethylation is mediated by Jmjd3 and Utx proteins. We isolated Xenopus orthologs of Ezh2, Eed, Jmjd3 and Utx, and examined their expression and function to explore roles for the H3K27 trimethylation/demethylation in frog development. Their expression all initially localizes in the animal hemisphere at the blastula stages, and subsequently remains faintly in the ectoderm at the gastrula stages. After the neural tube closure, expression of Ezh2, Eed, Jmjd3 and Utx becomes remarkable in the brain, eye, spinal cord, and migrating head neural crest cells. Overexpression of either Ezh2 or Jmjd3 perturbed expression of neural crest genes, Slug and Twist, but did not apparently affect gene expression in the central nervous system, such as Otx2 and Six3 in the forebrain and midbrain, Krox20 in the hindbrain, and HoxB9 in the spinal cord, at the neurula stages. Alcian blue staining revealed that these over-expressed embryos later developed abnormal cranial cartilage structures at the tadpole stages. These results showed that both the genes for trimethylation and demethylation of H3K27 are co-expressed in the neural tissues, and suggests that the regulated balance between these two activities is critical for the neural crest development.

During the chordate evolution, genome duplications have produced a number of paralogous gene pairs. While genetic studies often discover functional compensation between paralog genes, little is known about regulatory mechanisms underlying their interactions. In this study we addressed this issue by studying regulation of vertebrate Pax2, Pax5 and Pax8 genes that have evolved from a single ancestral gene of early cephalochordates and exhibit overlapping but different expression during eye, ear, brain and kidney development. Comparison of human, chicken and Xenopus genomic sequences in the Pax2/5/8 loci and high-throughput transgenic analysis in Xenopus identified 16 enhancers for Pax2, 7 for Pax5, and 3 for Pax8. Interestingly, Pax5 and Pax8 are associated with enhancers that are active in the tissues expressing Pax2, but not Pax5 or Pax8. One of such Pax5 enhancers, which is active in the pronephros, shares sequence features with a Pax2 pronephric enhancer. These two appear to have evolved from the common ancestral enhancer through the genome duplication. This Pax5 enhancer implied a hidden potential of this gene to be activated in the pronephros, and we actually found that Pax2-knockdowned embryos express Pax5 in the pronephros to rescue expression of WT1, a downstream target of Pax2 in normal embryos. These results revealed a dynamic balancing mechanism between the paralogs, which involves ectopic gene activation and may reflect regulatory mechanisms of their ancestral gene.

Ground based observations of NO_2 and O_3 using visible spectrometers have been made at Syowa Station (69°S) since March 1990. The NO_2 slant column amounts observed at the solar zenith angle of 90°(sunrise and sunset) exhibit a large seasonal variation, decreasing with shortening daylight hours in fall, reaching a minimum of about 1×(10)^(cm)^ in midwinter and increasing to a maximum of about 16×(10)^(cm)^ in midsummer. The recovery of NO_2 in spring is 2-3 times slower than the fall decay, because the heterogeneous chemistry on PSCs, the conversion of NO_x into HNO_3,is effective from midwinter to early spring, and NO_y is removed from the stratosphere through gravitational sedimentation of PSCs particles. In late spring the PSCs decrease and NO_2 increases because of the transport of air from lower latitudes. Following the eruption of Mt. Pinatubo (15°N) in June 1991,the NO_2 amounts in midsummer of 1991 were lower by 20-30% than those in 1990. The NO_2 amount in midsummer increased year after year and recovered to the 1990 level by 1994. The NO_2 amount in fall also recovered since 1992; however, the rate of recovery is smaller than in summer. The NO_2 level in fall of 1993 was as low as those in 1992. On the other hand, the rate of recovery in the box model using the observed aerosol surface areas by SAGE (S. SOLOMON et al., J. Geophys. Res., 99,3509,1994) is larger than that observed at McMurdo and Syowa Stations. This may suggest that another factor delays the NO_2 recovery in fall. NO_2 levels in winter have a positive correlation with temperature at 20 and 70 hPa (about 25 and 18km, respectively). In colder winters such as 1990 and 1993,the NO_2 levels were lower than those in other years. This may be due to the difference in the location of the polar vortex boundary relative to Syowa Station. The NO_2 amount and the temperature are considered to be lower deep inside the vortex as compared to near the boundary. In spring, the NO_2 and O_3 amounts and the temperature largely change when the vortex boundary crosses over Syowa Station. The vortex boundary can be clearly seen in the total O_3 maps obtained by TOMS. The monthly mean O_3 values suggest that Syowa Station was located deep inside of the vortex in October 1992 when the NO_2 level was lower than in other years. In contrast, Syowa Station was located near the edge of the vortex in October 1991 when the NO_2 level was much higher.