Genomic screening of fish-specific genes in gnathostomes and their functions in fin development.

In this study, we comprehensively searched for fish-specific genes in gnathostomes that contribute to development of the fin, a fish-specific trait. Many previous reports suggested that animal group-specific genes are often important for group-specific traits. Clarifying the roles of fish-specific genes in fin development of gnathostomes, for example, can help elucidate the mechanisms underlying the formation of this trait. We first identified 91 fish-specific genes in gnathostomes by comparing the gene repertoire in 16 fish and 35 tetrapod species. RNA-seq analysis narrowed down the 91 candidates to 33 genes that were expressed in the developing pectoral fin. We analyzed the functions of approximately half of the candidate genes by loss-of-function analysis in zebrafish. We found that some of the fish-specific and fin development-related genes, including fgf24 and and1/and2, play roles in fin development. In particular, the newly identified fish-specific gene qkia is expressed in the developing fin muscle and contributes to muscle morphogenesis in the pectoral fin as well as body trunk. These results indicate that the strategy of identifying animal group-specific genes is functional and useful. The methods applied here could be used in future studies to identify trait-associated genes in other animal groups.

Hagfish genome elucidates vertebrate whole-genome duplication events and their evolutionary consequences.

Polyploidy or whole-genome duplication (WGD) is a major event that drastically reshapes genome architecture and is often assumed to be causally associated with organismal innovations and radiations. The 2R hypothesis suggests that two WGD events (1R and 2R) occurred during early vertebrate evolution. However, the timing of the 2R event relative to the divergence of gnathostomes (jawed vertebrates) and cyclostomes (jawless hagfishes and lampreys) is unresolved and whether these WGD events underlie vertebrate phenotypic diversification remains elusive. Here we present the genome of the inshore hagfish, Eptatretus burgeri. Through comparative analysis with lamprey and gnathostome genomes, we reconstruct the early events in cyclostome genome evolution, leveraging insights into the ancestral vertebrate genome. Genome-wide synteny and phylogenetic analyses support a scenario in which 1R occurred in the vertebrate stem-lineage during the early Cambrian, and 2R occurred in the gnathostome stem-lineage, maximally in the late Cambrian-earliest Ordovician, after its divergence from cyclostomes. We find that the genome of stem-cyclostomes experienced an additional independent genome triplication. Functional genomic and morphospace analyses demonstrate that WGD events generally contribute to developmental evolution with similar changes in the regulatory genome of both vertebrate groups. However, appreciable morphological diversification occurred only in the gnathostome but not in the cyclostome lineage, calling into question the general expectation that WGDs lead to leaps of bodyplan complexity.

Distinguishing Evolutionary Conservation from Derivedness.

While the concept of "evolutionary conservation" has enabled biologists to explain many ancestral features and traits, it has also frequently been misused to evaluate the degree of changes from a common ancestor, or "derivedness". We propose that the distinction of these two concepts allows us to properly understand phenotypic and organismal evolution. From a methodological aspect, "conservation" mainly considers genes or traits which species have in common, while "derivedness" additionally covers those that are not commonly shared, such as novel or lost traits and genes to evaluate changes from the time of divergence from a common ancestor. Due to these differences, while conservation-oriented methods are effective in identifying ancestral features, they may be prone to underestimating the overall changes accumulated during the evolution of certain lineages. Herein, we describe our recently developed method, "transcriptomic derivedness index", for estimating the phenotypic derivedness of embryos with a molecular approach using the whole-embryonic transcriptome as a phenotype. Although echinoderms are often considered as highly derived species, our analyses with this method showed that their embryos, at least at the transcriptomic level, may not be much more derived than those of chordates. We anticipate that the future development of derivedness-oriented methods could provide quantitative indicators for finding highly/lowly evolvable traits.

How can recapitulation be reconciled with modern concepts of evolution?

To understand Haeckel's idea of recapitulation with modern evolutionary biology, one has to realize how evolutionarily conserved embryonic stages appear sequentially in developmental processes as chains of causality. Whether the idea of evolution was accepted or not, Haeckel and von Baer commonly saw an importance of a particularly conserved mid-embryonic stage in biphasic development of metazoans, the phylotype, that defines an animal phylum as the developmental source of a basic body plan. In an evolutionary context, the phylotypic stage was once understood by Haeckel to reflect the common ancestor of animal phyla, which went through hypermorphosis independently into various phyla. Recent comprehensive molecular studies, however, accumulated data to refute this idea. The conserved embryonic pattern does not reflect an ancestral adult morphology but appears to have arisen primarily as an embodiment of developmental constraints established through evolutionary processes. How the developmental burden results in a nested series of constraints will solve the recapitulative tendency of developmental programs.

The developmental hourglass model and recapitulation: An attempt to integrate the two models.

Recapitulation is a hypothetical concept that assumes embryogenesis of an animal parallels its own phylogenetic history, sequentially developing from more ancestral features to more derived ones. This concept predicts that the earliest developmental stage of various animals should represent the most evolutionarily conserved patterns. Recent transcriptome-based studies, on the other hand, have reported that mid-embryonic, organogenetic periods show the highest level of conservation (the developmental hourglass model). This, however, does not rule out the possibility that recapitulation would still be detected after the mid-embryonic period. In accordance with this, recapitulation-like morphological features are enriched in late developmental stages. Moreover, our recent chromatin accessibility-based study provided molecular evidence for recapitulation in the mid-to-late embryogenesis of vertebrates, as newly evolved gene regulatory elements tended to be activated at late embryonic stages. In this review, we revisit the recapitulation hypothesis, together with recent molecular-based studies that support the developmental hourglass model. We contend that the recapitulation hypothesis does not entirely contradict the developmental hourglass model and that these two may even coexist in later embryonic stages of vertebrates. Finally, we review possible mechanisms underlying the recapitulation pattern of chromatin accessibility together with the hourglass-like evolutionary conservation in vertebrate embryogenesis.

Derivedness Index for Estimating Degree of Phenotypic Evolution of Embryos: A Study of Comparative Transcriptomic Analyses of Chordates and Echinoderms.

Species retaining ancestral features, such as species called living fossils, are often regarded as less derived than their sister groups, but such discussions are usually based on qualitative enumeration of conserved traits. This approach creates a major barrier, especially when quantifying the degree of phenotypic evolution or degree of derivedness, since it focuses only on commonly shared traits, and newly acquired or lost traits are often overlooked. To provide a potential solution to this problem, especially for inter-species comparison of gene expression profiles, we propose a new method named "derivedness index" to quantify the degree of derivedness. In contrast to the conservation-based approach, which deals with expressions of commonly shared genes among species being compared, the derivedness index also considers those that were potentially lost or duplicated during evolution. By applying our method, we found that the gene expression profiles of penta-radial phases in echinoderm tended to be more highly derived than those of the bilateral phase. However, our results suggest that echinoderms may not have experienced much larger modifications to their developmental systems than chordates, at least at the transcriptomic level. In vertebrates, we found that the mid-embryonic and organogenesis stages were generally less derived than the earlier or later stages, indicating that the conserved phylotypic period is also less derived. We also found genes that potentially explain less derivedness, such as Hox genes. Finally, we highlight technical concerns that may influence the measured transcriptomic derivedness, such as read depth and library preparation protocols, for further improvement of our method through future studies. We anticipate that this index will serve as a quantitative guide in the search for constrained developmental phases or processes.

MeCP2 controls neural stem cell fate specification through miR-199a-mediated inhibition of BMP-Smad signaling.

Rett syndrome (RTT) is a severe neurological disorder, with impaired brain development caused by mutations in MECP2; however, the underlying mechanism remains elusive. We know from previous work that MeCP2 facilitates the processing of a specific microRNA, miR-199a, by associating with the Drosha complex to regulate neuronal functions. Here, we show that the MeCP2/miR-199a axis regulates neural stem/precursor cell (NS/PC) differentiation. A shift occurs from neuronal to astrocytic differentiation of MeCP2- and miR-199a-deficient NS/PCs due to the upregulation of a miR-199a target, Smad1, a downstream transcription factor of bone morphogenetic protein (BMP) signaling. Moreover, miR-199a expression and treatment with BMP inhibitors rectify the differentiation of RTT patient-derived NS/PCs and development of brain organoids, respectively, suggesting that facilitation of BMP signaling accounts for the impaired RTT brain development. Our study illuminates the molecular pathology of RTT and reveals the MeCP2/miR-199a/Smad1 axis as a potential therapeutic target for RTT.

Genomic insights of body plan transitions from bilateral to pentameral symmetry in Echinoderms.

Echinoderms are an exceptional group of bilaterians that develop pentameral adult symmetry from a bilaterally symmetric larva. However, the genetic basis in evolution and development of this unique transformation remains to be clarified. Here we report newly sequenced genomes, developmental transcriptomes, and proteomes of diverse echinoderms including the green sea urchin (L. variegatus), a sea cucumber (A. japonicus), and with particular emphasis on a sister group of the earliest-diverged echinoderms, the feather star (A. japonica). We learned that the last common ancestor of echinoderms retained a well-organized Hox cluster reminiscent of the hemichordate, and had gene sets involved in endoskeleton development. Further, unlike in other animal groups, the most conserved developmental stages were not at the body plan establishing phase, and genes normally involved in bilaterality appear to function in pentameric axis development. These results enhance our understanding of the divergence of protostomes and deuterostomes almost 500 Mya.

Recapitulation-like developmental transitions of chromatin accessibility in vertebrates.

The relationship between development and evolution has been a central theme in evolutionary developmental biology. Across the vertebrates, the most highly conserved gene expression profiles are found at mid-embryonic, organogenesis stages, whereas those at earlier and later stages are more diverged. This hourglass-like pattern of divergence does not necessarily rule out the possibility that gene expression profiles that are more evolutionarily derived appear at later stages of development; however, no molecular-level evidence of such a phenomenon has been reported. To address this issue, we compared putative gene regulatory elements among different species within a phylum. We made a genome-wide assessment of accessible chromatin regions throughout embryogenesis in three vertebrate species (mouse, chicken, and medaka) and estimated the evolutionary ages of these regions to define their evolutionary origins on the phylogenetic tree. In all the three species, we found that genomic regions tend to become accessible in an order that parallels their phylogenetic history, with evolutionarily newer gene regulations activated at later developmental stages. This tendency was restricted only after the mid-embryonic, phylotypic periods. Our results imply a phylogenetic hierarchy of putative regulatory regions, in which their activation parallels the phylogenetic order of their appearance. One evolutionary mechanism that may explain this phenomenon is that newly introduced regulatory elements are more likely to survive if activated at later stages of embryogenesis. Possible relationships between this phenomenon and the so-called recapitulation are discussed.

Embryonic lethality is not sufficient to explain hourglass-like conservation of vertebrate embryos.

BACKGROUND: Understanding the general trends in developmental changes during animal evolution, which are often associated with morphological diversification, has long been a central issue in evolutionary developmental biology. Recent comparative transcriptomic studies revealed that gene expression profiles of mid-embryonic period tend to be more evolutionarily conserved than those in earlier or later periods. While the hourglass-like divergence of developmental processes has been demonstrated in a variety of animal groups such as vertebrates, arthropods, and nematodes, the exact mechanism leading to this mid-embryonic conservation remains to be clarified. One possibility is that the mid-embryonic period (pharyngula period in vertebrates) is highly prone to embryonic lethality, and the resulting negative selections lead to evolutionary conservation of this phase. Here, we tested this "mid-embryonic lethality hypothesis" by measuring the rate of lethal phenotypes of three different species of vertebrate embryos subjected to two kinds of perturbations: transient perturbations and genetic mutations. RESULTS: By subjecting zebrafish (Danio rerio), African clawed frog (Xenopus laevis), and chicken (Gallus gallus) embryos to transient perturbations, namely heat shock and inhibitor treatments during three developmental periods [early (represented by blastula and gastrula), pharyngula, and late], we found that the early stages showed the highest rate of lethal phenotypes in all three species. This result was corroborated by perturbation with genetic mutations. By tracking the survival rate of wild-type embryos and embryos with genetic mutations induced by UV irradiation in zebrafish and African clawed frogs, we found that the highest decrease in survival rate was at the early stages particularly around gastrulation in both these species. CONCLUSION: In opposition to the "mid-embryonic lethality hypothesis," our results consistently showed that the stage with the highest lethality was not around the conserved pharyngula period, but rather around the early period in all the vertebrate species tested. These results suggest that negative selection by embryonic lethality could not explain hourglass-like conservation of animal embryos. This highlights the potential contribution of alternative mechanisms such as the diversifying effect of positive selections against earlier and later stages, and developmental constraints which lead to conservation of mid-embryonic stages.

Constrained vertebrate evolution by pleiotropic genes.

Despite morphological diversification of chordates over 550 million years of evolution, their shared basic anatomical pattern (or 'bodyplan') remains conserved by unknown mechanisms. The developmental hourglass model attributes this to phylum-wide conserved, constrained organogenesis stages that pattern the bodyplan (the phylotype hypothesis); however, there has been no quantitative testing of this idea with a phylum-wide comparison of species. Here, based on data from early-to-late embryonic transcriptomes collected from eight chordates, we suggest that the phylotype hypothesis would be better applied to vertebrates than chordates. Furthermore, we found that vertebrates' conserved mid-embryonic developmental programmes are intensively recruited to other developmental processes, and the degree of the recruitment positively correlates with their evolutionary conservation and essentiality for normal development. Thus, we propose that the intensively recruited genetic system during vertebrates' organogenesis period imposed constraints on its diversification through pleiotropic constraints, which ultimately led to the common anatomical pattern observed in vertebrates.

Evolutionary acquisition of promoter-associated non-coding RNA (pancRNA) repertoires diversifies species-dependent gene activation mechanisms in mammals.

BACKGROUND: Recent transcriptome analyses have shown that long non-coding RNAs (ncRNAs) play extensive roles in transcriptional regulation. In particular, we have reported that promoter-associated ncRNAs (pancRNAs) activate the partner gene expression via local epigenetic changes. RESULTS: Here, we identify thousands of genes under pancRNA-mediated transcriptional activation in five mammalian species in common. In the mouse, 1) pancRNA-partnered genes confined their expression pattern to certain tissues compared to pancRNA-lacking genes, 2) expression of pancRNAs was significantly correlated with the enrichment of active chromatin marks, H3K4 trimethylation and H3K27 acetylation, at the promoter regions of the partner genes, 3) H3K4me1 marked the pancRNA-partnered genes regardless of their expression level, and 4) C- or G-skewed motifs were exclusively overrepresented between-200 and-1 bp relative to the transcription start sites of the pancRNA-partnered genes. More importantly, the comparative transcriptome analysis among five different mammalian species using a total of 25 counterpart tissues showed that the overall pancRNA expression profile exhibited extremely high species-specificity compared to that of total mRNA, suggesting that interspecies difference in pancRNA repertoires might lead to the diversification of mRNA expression profiles. CONCLUSIONS: The present study raises the interesting possibility that the gain and/or loss of gene-activation-associated pancRNA repertoires, caused by formation or disruption of the genomic GC-skewed structure in the course of evolution, finely shape the tissue-specific pattern of gene expression according to a given species.

miR-199a Links MeCP2 with mTOR Signaling and Its Dysregulation Leads to Rett Syndrome Phenotypes.

Rett syndrome (RTT) is a neurodevelopmental disorder caused by MECP2 mutations. Although emerging evidence suggests that MeCP2 deficiency is associated with dysregulation of mechanistic target of rapamycin (mTOR), which functions as a hub for various signaling pathways, the mechanism underlying this association and the molecular pathophysiology of RTT remain elusive. We show here that MeCP2 promotes the posttranscriptional processing of particular microRNAs (miRNAs) as a component of the microprocessor Drosha complex. Among the MeCP2-regulated miRNAs, we found that miR-199a positively controls mTOR signaling by targeting inhibitors for mTOR signaling. miR-199a and its targets have opposite effects on mTOR activity, ameliorating and inducing RTT neuronal phenotypes, respectively. Furthermore, genetic deletion of miR-199a-2 led to a reduction of mTOR activity in the brain and recapitulated numerous RTT phenotypes in mice. Together, these findings establish miR-199a as a critical downstream target of MeCP2 in RTT pathogenesis by linking MeCP2 with mTOR signaling.

Gene activation-associated long noncoding RNAs function in mouse preimplantation development.

In mice, zygotic activation occurs for a wide variety of genes, mainly at the 2-cell stage. Long noncoding RNAs (lncRNAs) are increasingly being recognized as modulators of gene expression. In this study, directional RNA-seq of MII oocytes and 2-cell embryos identified more than 1000 divergently transcribed lncRNA/mRNA gene pairs. Expression of these bidirectional promoter-associated noncoding RNAs (pancRNAs) was strongly associated with the upregulation of their cognate genes. Conversely, knockdown of three abundant pancRNAs led to reduced mRNA expression, accompanied by sustained DNA methylation even in the presence of enzymes responsible for DNA demethylation. In particular, microinjection of siRNA against the abundant pancRNA partner of interleukin 17d (Il17d) mRNA at the 1-cell stage caused embryonic lethality, which was rescued by supplying IL17D protein in vitro at the 4-cell stage. Thus, this novel class of lncRNAs can modulate the transcription machinery in cis to activate zygotic genes and is important for preimplantation development.

Epigenetic setting and reprogramming for neural cell fate determination and differentiation.

In the mammalian brain, epigenetic mechanisms are clearly involved in the regulation of self-renewal of neural stem cells and the derivation of their descendants, i.e. neurons, astrocytes and oligodendrocytes, according to the developmental timing and the microenvironment, the 'niche'. Interestingly, local epigenetic changes occur, concomitantly with genome-wide level changes, at a set of gene promoter regions for either down- or upregulation of the gene. In addition, intergenic regions also sensitize the availability of epigenetic modifiers, which affects gene expression through a relatively long-range chromatinic interaction with the transcription regulatory machineries including non-coding RNA (ncRNA) such as promoter-associated ncRNA and enhancer ncRNA. We show that such an epigenetic landscape in a neural cell is statically but flexibly formed together with a variable combination of generally and locally acting nuclear molecules including master transcription factors and cell-cycle regulators. We also discuss the possibility that revealing the epigenetic regulation by the local DNA-RNA-protein assemblies would promote methodological innovations, e.g. neural cell reprogramming, engineering and transplantation, to manipulate neuronal and glial cell fates for the purpose of medical use of these cells.

Bidirectional promoters are the major source of gene activation-associated non-coding RNAs in mammals.

BACKGROUND: The majority of non-coding RNAs (ncRNAs) involved in mRNA metabolism in mammals have been believed to downregulate the corresponding mRNA expression level in a pre- or post-transcriptional manner by forming short or long ncRNA-mRNA duplex structures. Information on non-duplex-forming long ncRNAs is now also rapidly accumulating. To examine the directional properties of transcription at the whole-genome level, we performed directional RNA-seq analysis of mouse and chimpanzee tissue samples. RESULTS: We found that there is only about 1% of the genome where both the top and bottom strands are utilized for transcription, suggesting that RNA-RNA duplexes are not abundantly formed. Focusing on transcription start sites (TSSs) of protein-coding genes revealed that a significant fraction of them contain switching-points that separate antisense- and sense-biased transcription, suggesting that head-to-head transcription is more prevalent than previously thought. More than 90% of head-to-head type promoters contain CpG islands. Moreover, CCG and CGG repeats are significantly enriched in the upstream regions and downstream regions, respectively, of TSSs located in head-to-head type promoters. Genes with tissue-specific promoter-associated ncRNAs (pancRNAs) show a positive correlation between the expression of their pancRNA and mRNA, which is in accord with the proposed role of pancRNA in facultative gene activation, whereas genes with constitutive expression generally lack pancRNAs. CONCLUSIONS: We propose that single-stranded ncRNA resulting from head-to-head transcription at GC-rich sequences regulates tissue-specific gene expression.

Single-stranded noncoding RNAs mediate local epigenetic alterations at gene promoters in rat cell lines.

A growing number of noncoding RNAs (ncRNAs) are thought to be involved in sequence-specific alterations of epigenetic processes, mostly causing gene repression. In this study, promoter-associated ncRNAs (pancRNAs >200 nucleotides in size) that were endogenously generated from the sense strand at Map2b, antisense strand at Nefl, and both strands at Vim were investigated regarding their epigenetic potential as positive or negative regulators in rat pheochromocytoma (PC12) and fibroblast (normal rat kidney) cell lines. The respective antisense pancRNAs were associated with several active chromatin marks at the Nefl and Vim promoters. Forced expression of fragments expressing the antisense pancRNAs caused sequence-specific DNA demethylation, whereas a decrease of expression induced methylation of the same sequences. In contrast, perturbing the expression of the two sense pancRNAs did not change the DNA methylation status. These results suggest that a fraction of naturally occurring ncRNAs acts in cis as a single-stranded form and that the transcriptional orientation of pancRNA is important for the establishment of sequence-specific epigenetic modifications consistent with open chromatin structure.