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.

Single-cell transcriptome atlas of Drosophila gastrula 2.0.

During development, positional information directs cells to specific fates, leading them to differentiate with their own transcriptomes and express specific behaviors and functions. However, the mechanisms underlying these processes in a genome-wide view remain ambiguous, partly because the single-cell transcriptomic data of early developing embryos containing accurate spatial and lineage information are still lacking. Here, we report a single-cell transcriptome atlas of Drosophila gastrulae, divided into 77 transcriptomically distinct clusters. We find that the expression profiles of plasma-membrane-related genes, but not those of transcription-factor genes, represent each germ layer, supporting the nonequivalent contribution of each transcription-factor mRNA level to effector gene expression profiles at the transcriptome level. We also reconstruct the spatial expression patterns of all genes at the single-cell stripe level as the smallest unit. This atlas is an important resource for the genome-wide understanding of the mechanisms by which genes cooperatively orchestrate Drosophila gastrulation.

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.

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.

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.

Hair follicle epidermal stem cells define a niche for tactile sensation.

The heterogeneity and compartmentalization of stem cells is a common principle in many epithelia, and is known to function in epithelial maintenance, but its other physiological roles remain elusive. Here we show transcriptional and anatomical contributions of compartmentalized epidermal stem cells in tactile sensory unit formation in the mouse hair follicle. Epidermal stem cells in the follicle upper-bulge, where mechanosensory lanceolate complexes innervate, express a unique set of extracellular matrix (ECM) and neurogenesis-related genes. These epidermal stem cells deposit an ECM protein called EGFL6 into the collar matrix, a novel ECM that tightly ensheathes lanceolate complexes. EGFL6 is required for the proper patterning, touch responses, and αv integrin-enrichment of lanceolate complexes. By maintaining a quiescent original epidermal stem cell niche, the old bulge, epidermal stem cells provide anatomically stable follicle-lanceolate complex interfaces, irrespective of the stage of follicle regeneration cycle. Thus, compartmentalized epidermal stem cells provide a niche linking the hair follicle and the nervous system throughout the hair cycle.

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.

Shark genomes provide insights into elasmobranch evolution and the origin of vertebrates.

Modern cartilaginous fishes are divided into elasmobranchs (sharks, rays and skates) and chimaeras, and the lack of established whole-genome sequences for the former has prevented our understanding of early vertebrate evolution and the unique phenotypes of elasmobranchs. Here we present de novo whole-genome assemblies of brownbanded bamboo shark and cloudy catshark and an improved assembly of the whale shark genome. These relatively large genomes (3.8-6.7 Gbp) contain sparse distributions of coding genes and regulatory elements and exhibit reduced molecular evolutionary rates. Our thorough genome annotation revealed Hox C genes previously hypothesized to have been lost, as well as distinct gene repertories of opsins and olfactory receptors that would be associated with adaptation to unique underwater niches. We also show the early establishment of the genetic machinery governing mammalian homoeostasis and reproduction at the jawed vertebrate ancestor. This study, supported by genomic, transcriptomic and epigenomic resources, provides a foundation for the comprehensive, molecular exploration of phenotypes unique to sharks and insights into the evolutionary origins of vertebrates.

CTCF binding landscape in jawless fish with reference to Hox cluster evolution.

The nuclear protein CCCTC-binding factor (CTCF) contributes as an insulator to chromatin organization in animal genomes. Currently, our knowledge of its binding property is confined mainly to mammals. In this study, we identified CTCF homologs in extant jawless fishes and performed ChIP-seq for the CTCF protein in the Arctic lamprey. Our phylogenetic analysis suggests that the lamprey lineage experienced gene duplication that gave rise to its unique paralog, designated CTCF2, which is independent from the previously recognized duplication between CTCF and CTCFL. The ChIP-seq analysis detected comparable numbers of CTCF binding sites between lamprey, chicken, and human, and revealed that the lamprey CTCF protein binds to the two-part motif, consisting of core and upstream motifs previously reported for mammals. These findings suggest that this mode of CTCF binding was established in the last common ancestor of extant vertebrates (more than 500 million years ago). We analyzed CTCF binding inside Hox clusters, which revealed a reinforcement of CTCF binding in the region spanning Hox1-4 genes that is unique to lamprey. Our study provides not only biological insights into the antiquity of CTCF-based epigenomic regulation known in mammals but also a technical basis for comparative epigenomic studies encompassing the whole taxon Vertebrata.

Gene expression profiling of granule cells and Purkinje cells in the zebrafish cerebellum.

The structure of the neural circuitry of the cerebellum, which functions in some types of motor learning and coordination, is generally conserved among vertebrates. However, some cerebellar features are species specific. It is not clear which genes are involved in forming these conserved and species-specific structures and functions. This study uses zebrafish transgenic larvae expressing fluorescent proteins in granule cells, Purkinje cells, or other cerebellar neurons and glial cells to isolate each type of cerebellar cells by fluorescence-activated cell sorting and to profile their gene expressions by RNA sequencing and in situ hybridization. We identify genes that are upregulated in granule cells or Purkinje cells, including many genes that are also expressed in mammalian cerebella. Comparison of the transcriptomes in granule cells and Purkinje cells in zebrafish larvae reveals that more developmental genes are expressed in granule cells, whereas more neuronal-function genes are expressed in Purkinje cells. We show that some genes that are upregulated in granule cells or Purkinje cells are also expressed in the cerebellum-like structures. Our data provide a platform for understanding the development and function of the cerebellar neural circuits in zebrafish and the evolution of cerebellar circuits in vertebrates. J. Comp. Neurol. 525:1558-1585, 2017. © 2016 Wiley Periodicals, Inc.

Structural, super-resolution microscopy analysis of paraspeckle nuclear body organization.

Paraspeckles are nuclear bodies built on the long noncoding RNA Neat1, which regulates a variety of physiological processes including cancer progression and corpus luteum formation. To obtain further insight into the molecular basis of the function of paraspeckles, we performed fine structural analyses of these nuclear bodies using structural illumination microscopy. Notably, paraspeckle proteins are found within different layers along the radially arranged bundles of Neat1 transcripts, forming a characteristic core-shell spheroidal structure. In cells lacking the RNA binding protein Fus, paraspeckle spheroids are disassembled into smaller particles containing Neat1, which are diffusely distributed in the nucleoplasm. Sequencing analysis of RNAs purified from paraspeckles revealed that AG-rich transcripts associate with Neat1, which are distributed along the shell of the paraspeckle spheroids. We propose that paraspeckles sequester core components inside the spheroids, whereas the outer surface associates with other components in the nucleoplasm to fulfill their function.

Optimization and cost-saving in tagmentation-based mate-pair library preparation and sequencing.

In de novo genome sequencing, mate-pair reads are crucial for scaffolding assembled contigs. However, preparation of mate-pair libraries is not a trivial task, even when using one of the latest approaches, the Nextera Mate Pair Sample Prep Kit from Illumina. To reduce cost and enhance library yield and fidelity when using this kit, we have modified the manufacturer's protocol based on (i) variable tagmentation conditions, (ii) intensive DNA shearing to decrease library insert length, and (iii) sequencing on an Illumina HiSeq with >150 cycles. Finally, we provide additional suggestions for further improvement in the application of this kit.

A LIM-homeobox gene is required for differentiation of Wnt-expressing cells at the posterior end of the planarian body.

Planarians have high regenerative ability, which is dependent on pluripotent adult somatic stem cells called neoblasts. Recently, canonical Wnt/ß-catenin signaling was shown to be required for posterior specification, and Hedgehog signaling was shown to control anterior-posterior polarity via activation of the Djwnt1/P-1 gene at the posterior end of planarians. Thus, various signaling molecules play an important role in planarian stem cell regulation. However, the molecular mechanisms directly involved in stem cell differentiation have remained unclear. Here, we demonstrate that one of the planarian LIM-homeobox genes, Djislet, is required for the differentiation of Djwnt1/P-1-expressing cells from stem cells at the posterior end. RNA interference (RNAi)-treated planarians of Djislet [Djislet(RNAi)] show a tail-less phenotype. Thus, we speculated that Djislet might be involved in activation of the Wnt signaling pathway in the posterior blastema. When we carefully examined the expression pattern of Djwnt1/P-1 by quantitative real-time PCR during posterior regeneration, we found two phases of Djwnt1/P-1 expression: the first phase was detected in the differentiated cells in the old tissue in the early stage of regeneration and then a second phase was observed in the cells derived from stem cells in the posterior blastema. Interestingly, Djislet is expressed in stem cell-derived DjPiwiA- and Djwnt1/P-1-expressing cells, and Djislet(RNAi) only perturbed the second phase. Thus, we propose that Djislet might act to trigger the differentiation of cells expressing Djwnt1/P-1 from stem cells.