CRISPR-Cas9-Mediated Gene Knockout in a Non-Model Sea Urchin, Heliocidaris crassispina.

Sea urchins have been used as model organisms in developmental biology research and the genomes of several sea urchin species have been sequenced. Recently, genome editing technologies have become available for sea urchins, and methods for gene knockout using the CRISPRCas9 system have been established. Heliocidaris crassispina is an important marine fishery resource with edible gonads. Although H. crassispina has been used as a biological research material, its genome has not yet been published, and it is a non-model sea urchin for molecular biology research. However, as recent advances in genome editing technology have facilitated genome modification in non-model organisms, we applied genome editing using the CRISPR-Cas9 system to H. crassispina. In this study, we targeted genes encoding ETS transcription factor (HcEts) and pigmentation-related polyketide synthase (HcPks1). Gene fragments were isolated using primers designed by inter-specific sequence comparisons within Echinoidea. When Ets gene was targeted using two sgRNAs, one successfully introduced mutations and impaired skeletogenesis. In the Pks1 gene knockout, when two sgRNAs targeting the close vicinity of the site corresponding to the target site that showed 100% mutagenesis efficiency of the Pks1 gene in Hemicentrotus pulcherrimus, mutagenesis was not observed. However, two other sgRNAs targeting distant sites efficiently introduced mutations. In addition, Pks1 knockout H. crassispina exhibited an albino phenotype in the pluteus larvae and adult sea urchins after metamorphosis. This indicates that the CRISPRCas9 system can be used to modify the genome of the non-model sea urchin H. crassispina.

Assembly of continuous high-resolution draft genome sequence of Hemicentrotus pulcherrimus using long-read sequencing.

The update of the draft genome assembly of sea urchin, Hemicentrotus pulcherrimus, which is widely studied in East Asia as a model organism of early development, was performed using Oxford nanopore long-read sequencing. The updated assembly provided ~600-Mb genome sequences divided into 2,163 contigs with N50 = 516 kb. BUSCO completeness score and transcriptome model mapping ratio (TMMR) of the present assembly were obtained as 96.5% and 77.8%, respectively. These results were more continuous with higher resolution than those by the previous version of H. pulcherrimus draft genome, HpulGenome_v1, where the number of scaffolds = 16,251 with a total of ~100 Mb, N50 = 143 kb, BUSCO completeness score = 86.1%, and TMMR = 55.4%. The obtained genome contained 36,055 gene models that were consistent with those in other echinoderms. Additionally, two tandem repeat sequences of early histone gene locus containing 47 copies and 34 copies of all histone genes, and 185 of the homologous sequences of the interspecifically conserved region of the Ars insulator, ArsInsC, were obtained. These results provide further advance for genome-wide research of development, gene regulation, and intranuclear structural dynamics of multicellular organisms using H. pulcherrimus.

Possibilities of skin coat color-dependent risks and risk factors of squamous cell carcinoma and deafness of domestic cats inferred via RNA-seq data.

The transcriptome data of skin cells from domestic cats with brown, orange, and white coats were analyzed using a public database to investigate the possible relationship between coat color-related gene expression and squamous cell carcinoma risk, as well as the mechanism of deafness in white cats. We found that the ratio of the expression level of genes suppressing squamous cell carcinoma to that of genes promoting squamous cell carcinoma might be considerably lower than the theoretical estimation in skin cells with orange and white coats in white-spotted cat. We also found the possibility of the frequent production of KIT lacking the first exon (d1KIT) in skin cells with white coats, and d1KIT production exhibited a substantial negative correlation with the expression of SOX10, which is essential for melanocyte formation and adjustment of hearing function. Additionally, the production of d1KIT was expected to be due to the insulating activity of the feline endogenous retrovirus 1 (FERV1) LTR in the first intron of KIT by its CTCF binding sequence repeat. These results contribute to basic veterinary research to understand the relationship between cat skin coat and disease risk, as well as the underlying mechanism.

The crucial role of CTCF in mitotic progression during early development of sea urchin.

CCCTC-binding factor (CTCF), an insulator protein with 11 zinc fingers, is enriched at the boundaries of topologically associated domains (TADs) in eukaryotic genomes. In this study, we isolated and analyzed the cDNAs encoding HpCTCF, the CTCF homolog in the sea urchin Hemicentrotus pulcherrimus, to investigate its expression patterns and functions during the early development of sea urchin. HpCTCF contains nine zinc fingers corresponding to fingers 2-10 of the vertebrate CTCF. Expression pattern analysis revealed that HpCTCF mRNA was detected at all developmental stages and in the entire embryo. Upon expressing the HpCTCF-GFP fusion protein in early embryos, we observed its uniform distribution within interphase nuclei. However, during mitosis, it disappeared from the chromosomes and subsequently reassembled on the chromosome during telophase. Moreover, the morpholino-mediated knockdown of HpCTCF resulted in mitotic arrest during the morula to blastula stage. Most of the arrested chromosomes were not phospholylated at serine 10 of histone H3, indicating that mitosis was arrested at the telophase by HpCTCF depletion. Furthermore, impaired sister chromatid segregation was observed using time-lapse imaging of HpCTCF-knockdown embryos. Thus, HpCTCF is essential for mitotic progression during the early development of sea urchins, especially during the telophase-to-interphase transition. However, the normal development of pluteus larvae in CRISPR-mediated HpCTCF-knockout embryos suggests that disruption of zygotic HpCTCF expression has little effect on embryonic and larval development.

Mathematical model of structural changes in nuclear speckle.

Nuclear speckles are nuclear bodies consisting of populations of small and irregularly shaped droplet-like molecular condensates that contain various splicing factors. Recent experiments have revealed the following structural features of nuclear speckles: (I) Each molecular condensate contains SON and SRRM2 proteins, and MALAT1 non-coding RNA surrounds these condensates; (II) During normal interphase of the cell cycle in multicellular organisms, these condensates are broadly distributed throughout the nucleus. In contrast, when cell transcription is suppressed, the condensates fuse and form strongly condensed spherical droplets; (III) SON is dispersed spatially in MALAT1 knocked-down cells and MALAT1 is dispersed in SON knocked-down cells because of the collapse of the nuclear speckles. However, the detailed interactions among the molecules that are mechanistically responsible for the structural variation remain unknown. In this study, a coarse-grained molecular dynamics model of the nuclear speckle was developed by considering the dynamics of SON, SRRM2, MALAT1, and pre-mRNA as representative components of the condensates. The simulations reproduced the structural changes, which were used to predict the interaction network among the representative components of the condensates.

Molecular dynamics analysis of biomolecular systems including nucleic acids.

With the recent progress in structural biology and genome biology, structural dynamics of molecular systems that include nucleic acids has attracted attention in the context of gene regulation. The structure-function relationship is an important topic that highlights the importance of the physicochemical properties of nucleotides, as well as that of amino acids in proteins. Simulations are a useful tool for the detailed analysis of molecular dynamics that complement experiments in molecular biology; however, molecular simulation of nucleic acids is less well developed than that of proteins partly due to the physical nature of nucleic acids. In this review, we briefly describe the current status and future directions of the field as a guide to promote collaboration between experimentalists and computational biologists.

Methylome data derived from maternal-zygotic DNA methyltransferase 3aa-/- zebrafish.

Genomic DNA methylation is an epigenetic marker mediated by DNA methyltransferases (Dnmts); in vertebrates, it comprises of a maintenance DNA methyltransferase, Dnmt1, and two de novo DNA methyltransferases (Dnmt3a and Dnmt3b). In zebrafish, there are two homologs of the mammalian Dnmt3a: Dnmt3aa and Dnmt3ab. A knockout (KO) mutant of zebrafish dnmt3aa was generated using the CRISPR/Cas9 genome-editing system as a new model for DNA methylation research. Since zebrafish dnmt3aa KO mutants were viable and fertile, a maternal-zygotic dnmt3aa deficient mutant (MZdnmt3aa) was generated. We performed whole-genome bisulfite sequencing (WGBS) to reveal the DNA methylation profile using this mutant and identified genomic regions with altered CpG methylation as differentially methylated regions (DMRs) in this mutant compared to those in the wild-type fish. We provided novel raw and processed datasets using the MZdnmt3aa KO mutant, and the raw data of WGBS are available through the Gene Expression Omnibus (GEO), accession number GSE178690.

Epigenetic-structural changes in X chromosomes promote Xic pairing during early differentiation of mouse embryonic stem cells.

X chromosome inactivation center (Xic) pairing occurs during the differentiation of embryonic stem (ES) cells from female mouse embryos, and is related to X chromosome inactivation, the circadian clock, intra-nucleus architecture, and metabolism. However, the mechanisms underlying the identification and approach of X chromosome pairs in the crowded nucleus are unclear. To elucidate the driving force of Xic pairing, we developed a coarse-grained molecular dynamics model of intranuclear chromosomes in ES cells and in cells 2 days after the onset of differentiation (2-day cells) by considering intrachromosomal epigenetic-structural feature-dependent mechanics. The analysis of the experimental data showed that X-chromosomes exhibit the rearrangement of their distributions of open/closed chromatin regions on their surfaces during cell differentiation. By simulating models where the excluded volume effects of closed chromatin regions are stronger than those of open chromatin regions, such rearrangement of open/closed chromatin regions on X-chromosome surfaces promoted the mutual approach of the Xic pair. These findings suggested that local intrachromosomal epigenetic features may contribute to the regulation of cell species-dependent differences in intranuclear architecture.

Mathematical model of chromosomal dynamics during DNA double strand break repair in budding yeast.

During the repair of double-strand breaks (DSBs) in DNA, active mobilizations for conformational changes in chromosomes have been widely observed in eukaryotes, from yeast to animal and plant cells. DSB-damaged loci in the yeast genome showed increased mobility and relocation to the nuclear periphery. However, the driving forces behind DSB-induced chromatin dynamics remain unclear. In this study, mathematical models of normal and DSB-damaged yeast chromosomes were developed to simulate their structural dynamics. The effects of histone degradation in the whole nucleus and the change in the physical properties of damaged loci due to the binding of SUMOylated repair proteins were considered in the model of DSB-induced chromosomes based on recent experimental results. The simulation results reproduced DSB-induced changes to structural and dynamical features by which the combination of whole nuclear histone degradation and the rigid structure formation of repair protein accumulations on damaged loci were suggested to be primary contributors to the process by which damaged loci are relocated to the nuclear periphery.

Identification of aberrant transcription termination at specific gene loci with DNA hypomethylated transcription termination sites caused by DNA methyltransferase deficiency.

CpG methylation of genomic DNA is a well-known repressive epigenetic marker in eukaryotic transcription, and DNA methylation of promoter regions is correlated with gene silencing. In contrast to the promoter regions, the function of DNA methylation during transcription termination remains to be elucidated. A recent study revealed that mouse DNA methyltransferase 3a (Dnmt3a) mainly functions in de novo methylation in the promoter and gene body regions, including transcription termination sites (TTSs), during development. To investigate the relationship between DNA methylation overlapping the TTSs and transcription termination, we performed bioinformatics analysis using six pre-existing Dnmt-/- mouse cell datasets: four types of neurons (three Dnmt3a-/- and one Dnmt1-/- mutants) and two types of embryonic fibroblasts (MEFs) (Dnmt3a-/- and Dnmt3b-/- mutants). Combined analyses using methylome and transcriptome data revealed that read counts downstream of hypomethylated TTSs were increased in three types of neurons (two Dnmt3a-/- and one Dnmt1-/- mutants). Among these, an increase in chimeric transcripts downstream of the TTSs was observed in Dnmt3a-/- mature olfactory sensory neurons and Dnmt3a-/- agouti-related peptide (protein)-producing neurons, thereby indicating that read-through occurs in hypomethylated TTSs at specific gene loci in these two mutants. Conversely, in Dnmt3a-/- MEFs, we detected reductions in read counts downstream of hypomethylated TTSs. These results indicate that the hypomethylation of TTSs can both positively and negatively regulate transcription termination, dependent on Dnmt and cell types. This study is the first to identify the aberrant termination of transcription at specific gene loci with DNA hypomethylated TTSs attributable to Dnmt deficiency.

Partial exogastrulation due to apical-basal polarity of F-actin distribution disruption in sea urchin embryo by omeprazole.

Gastrulation is a universal process in the morphogenesis of many animal embryos. Although morphological and molecular events in gastrulation have been well studied, the mechanical driving forces and underlying regulatory mechanisms are not fully understood. Here, we investigated the gastrulation of embryos of a sea urchin, Hemicentrotus pulcherrimus, which involves the invagination of a single-layered vegetal plate into the blastocoel. We observed that omeprazole, a proton pump inhibitor capable of perturbing the left-right asymmetry of sea urchin embryo, induced "partial exogastrulation" where the secondary invagination proceeds outward. During early gastrulation, intracellular apical-basal polarity of F-actin distribution in vegetal half was higher than those in animal half, while omeprazole treatment disturbed the apical-basal polarity of F-actin distribution in vegetal half. Furthermore, gastrulation stopped and even partial exogastrulation did not occur when F-actin polymerization or degradation in whole embryo was partially inhibited via RhoA or YAP1 knockout. A mathematical model of the early gastrulation reproduced the shapes of both normal and exogastrulating embryos using cell-dependent cytoskeletal features based on F-actin. Additionally, such cell position-dependent intracellular F-actin distributions might be regulated by intracellular pH distributions. Therefore, apical-basal polarity of F-actin distribution disrupted by omeprazole may induce the partial exogastrulation via anomalous secondary invagination.

Ribosomal protein L5 facilitates rDNA-bundled condensate and nucleolar assembly.

The nucleolus is the site of ribosome assembly and formed through liquid-liquid phase separation. Multiple ribosomal DNA (rDNA) arrays are bundled in the nucleolus, but the underlying mechanism and significance are unknown. In the present study, we performed high-content screening followed by image profiling with the wndchrm machine learning algorithm. We revealed that cells lacking a specific 60S ribosomal protein set exhibited common nucleolar disintegration. The depletion of RPL5 (also known as uL18), the liquid-liquid phase separation facilitator, was most effective, and resulted in an enlarged and un-separated sub-nucleolar compartment. Single-molecule tracking analysis revealed less-constrained mobility of its components. rDNA arrays were also unbundled. These results were recapitulated by a coarse-grained molecular dynamics model. Transcription and processing of ribosomal RNA were repressed in these aberrant nucleoli. Consistently, the nucleoli were disordered in peripheral blood cells from a Diamond-Blackfan anemia patient harboring a heterozygous, large deletion in RPL5 Our combinatorial analyses newly define the role of RPL5 in rDNA array bundling and the biophysical properties of the nucleolus, which may contribute to the etiology of ribosomopathy.

Structural dynamics of DNA depending on methylation pattern.

DNA methylation is associated with a number of biological phenomena, and plays crucial roles in epigenetic regulation of eukaryotic gene expression. It is also suggested that DNA methylation alters the mechanical properties of DNA molecules, which is likely to affect epigenetic regulation. However, it has not been systematically investigated how methylation changes the structural and dynamic features of DNA. In this research, to elucidate the effects of methylation on DNA mechanics, a fully atomic molecular dynamics simulation of double-stranded DNA with several methylation patterns was performed. Through the analysis of the relative positioning of the nucleotides (base-step variables), characteristic changes in terms of local flexibility were observed, which further affected the overall DNA geometry and stiffness. These findings may serve as a basis for a discussion on methylation-dependent DNA dynamics in physiological conditions.

Ultrasensitive Change in Nucleosome Binding by Multiple Phosphorylations to the Intrinsically Disordered Region of the Histone Chaperone FACT.

Facilitates chromatin transcription (FACT) is a histone chaperone that functions as a nucleosome remodeler and a chaperone. The two subunits of FACT, Spt16 and SSRP1, mediate multiple interactions between the subunits and components of the nucleosome. Among the interactions, the role of the DNA-binding domain in SSRP1 has not been characterized. We reported previously that the DNA-binding domain in Drosophila SSRP1 (dSSRP1) has multiple casein kinase II phosphorylation sites, and the DNA binding affinity of the domain changes sigmoidally in response to the degree of phosphorylation ("ultrasensitive response"). In this report, we explored the molecular mechanisms for the ultrasensitive response of the DNA-binding domain in dSSRP1 using the shortest fragment (AB-HMG, residues 434-624) responsible for nucleosome binding. AB-HMG contains two intrinsically disordered (ID) regions: the N-terminal part rich in acidic residues (AID) and the C-terminal part rich in basic residues (BID) followed by the HMG box. NMR and coarse-grained molecular dynamics simulations revealed a phosphorylation-dependent change in intramolecular contacts between the AID and BID-HMG, which is mediated by a hinge bending motion of AB-HMG to enable the ultrasensitive response. Ultrasensitivity generates two distinct forms of dSSRP1, which are high- and low-affinity nucleosome-binding forms. Drosophila FACT (dFACT) switches function according to the degree of phosphorylation of the AID in dSSRP1. We propose that dFACT in various phosphorylation states functions cooperatively to facilitate gene regulation in the context of the chromatin.

Affinity of rhodopsin to raft enables the aligned oligomer formation from dimers: Coarse-grained molecular dynamics simulation of disk membranes.

The visual photopigment protein rhodopsin (Rh) is a typical G protein-coupled receptor (GPCR) that initiates the phototransduction cascade in retinal disk membrane of rod-photoreceptor cells. Rh molecule has a tendency to form dimer, and the dimer tends to form rows, which is suggested to heighten phototransduction efficiency in single-photon regime. In addition, the dimerization confers Rh an affinity for lipid raft, i.e. raftophilicity. However, the mechanism by which Rh-dimer raftophilicity contributes to the organization of the higher order structure remains unknown. In this study, we performed coarse-grained molecular dynamics simulations of a disk membrane model containing unsaturated lipids, saturated lipids with cholesterol, and Rh-dimers. We described the Rh-dimers by two-dimensional particle populations where the palmitoyl moieties of each Rh exhibits raftophilicity. We simulated the structuring of Rh in a disk for two types of Rh-dimer, i.e., the most and second most stable Rh dimers, which exposes the raftophilic regions at the dimerization-interface (H1/H8 dimer) and two edges away from the interface (H4/H5 dimer), respectively. Our simulations revealed that only the H1/H8 dimer could form a row structure. A small number of raftophilic lipids recruited to and intercalated in a narrow space between H1/H8 dimers stabilize the side-by-side interaction between dimers in a row. Our results implicate that the nano-sized lipid raft domains act as a "glue" to organize the long row structures of Rh-dimers.

Histone Tail Dynamics in Partially Disassembled Nucleosomes During Chromatin Remodeling.

Nucleosomes are structural units of the chromosome consisting of DNA wrapped around histone proteins, and play important roles in compaction and regulation of the chromatin structure. While the structure and dynamics of canonical nucleosomes have been studied extensively, those of nucleosomes in intermediate states, that occur when their structure or positioning is modulated, have been less understood. In particular, the dynamic features of partially disassembled nucleosomes have not been discussed in previous studies. Using all-atom molecular dynamics simulations, in this study, we investigated the dynamics and stability of nucleosome structures lacking a histone-dimer. DNA in nucleosomes lacking a histone H2A/H2B dimer was drastically deformed due to loss of local interactions between DNA and histones. In contrast, conformation of DNA in nucleosomes lacking H3/H4 was similar to the canonical nucleosome, as the H2A C-terminal domain infiltrated the space originally occupied by the dissociated H3/H4 histones and restricted DNA dynamics in close proximity. Our results suggest that, besides histone chaperones, the intrinsic dynamics of nucleosomes support the exchange of H2A/H2B, which is significantly more frequent than that of H3/H4.

Establishment of knockout adult sea urchins by using a CRISPR-Cas9 system.

Sea urchins are used as a model organism for research on developmental biology and gene regulatory networks during early development. Gene knockdown by microinjection of morpholino antisense oligonucleotide (MASO) has been used to analyze gene function in early sea urchin embryos. However, as the effect of MASO is not long lasting, it is impossible to perturb genes expressed during late development by MASO. Recent advances in genome editing technologies have enabled gene modification in various organisms. We previously reported genome editing in the sea urchin Hemicentrotus pulcherrimus using zinc-finger nuclease (ZFN) and transcription activator-like effector nuclease (TALEN); however, the efficiencies of these technologies were not satisfactory. Here, we applied clustered regularly interspaced short palindromic repeat (CRISPR)-CRISPR-associated nuclease 9 (Cas9) technology to knock out the Pks1 gene in H. pulcherrimus. When sgRNAs targeting Pks1, which is required for the biosynthesis of larval pigment, were microinjected into fertilized eggs with SpCas9 mRNA, high-efficiency mutagenesis was achieved within 24 hr post fertilization and SpCas9/sgRNA-injected pluteus larvae had an albino phenotype. One of the sgRNAs yielded 100% mutagenesis efficiency, and no off-target effect was detected. In addition, the albino phenotype was maintained in juvenile sea urchins after metamorphosis, and the knockout sea urchins survived for at least one year and grew to albino adult sea urchins. These findings suggest that knockout adult sea urchins were successfully established and the CRISPR-Cas9 system is a feasible method for analyzing gene functions from late developmental to adult stage.

Verification of mathematical models of response threshold through statistical characterisation of the foraging activity in ant societies.

The concept of response threshold (RT) has been developed to explain task allocation in social insect colonies, wherein individual workers engage in tasks depending on their responsiveness to the task-related stimulus. Moreover, a mathematical model of RT has been proposed to explain data obtained from task allocation experiments; however, its applicability range warrants clarification through adequate quantitative analysis. Hence, we used an automatic measuring system to count passage events between a nest chamber and a foraging arena in five colonies of ants, Camponotus japonicus. The events were measured using radio-frequency identification tags attached to all workers of each colony. Here, we examined the detailed forms of i) labour distribution during foraging among workers in each colony and ii) the persistence of rank-order of foraging among workers. We found that labour distribution was characterized by a generalized gamma-distribution, indicating that only few workers carried out a large part of the workload. The rank-order of foraging activity among workers in each colony was maintained for a month and collapsed within a few months. We compared the obtained data with testable predictions of the RT model. The comparison indicated that proper evaluation of the mathematical model is required based on the obtained data.

Insulator Activities of Nucleosome-Excluding DNA Sequences without Bound Chromatin Looping Proteins.

Chromosomes consist of various domains with different transcriptional activities separated by chromatin boundary sequences such as insulator sequences. Recent studies suggested that CTCF or other chromatin loop-forming protein binding sequences represented typical insulators. Alternatively, some long nucleosome-excluding DNA sequences were also reported to exhibit insulator activities in yeast and sea urchin chromosomes, although specific binding of loop-forming proteins was not expected for them. However, the mechanism of the insulator activities of these sequences and the possibilities of similar insulators existing in other organisms remained unclear. In this study, we first constructed and performed simulations of a coarse-grained chromatin model containing nucleosome-rich and nucleosome-excluding DNA regions. We found that a long nucleosome-excluding region between two nucleosome-rich regions could markedly hinder the associations of two neighboring chromatin regions owing to the stronger long-term-averaged rigidity of the nucleosome-excluding region compared to that of nucleosome-rich regions. Subsequent analysis of the genome-wide nucleosome positioning, protein binding, and DNA rigidity in human cells revealed that some nucleosome-excluding rigid DNA sequences without bound chromatin looping proteins could exhibit insulator activities, functioning as chromatin boundaries in various regions of human chromosomes.

Diverse stochasticity leads a colony of ants to optimal foraging.

A mathematical model of garden ants (Lasius japonicus) is introduced herein to investigate the relationship between the distribution of the degree of stochasticity in following pheromone trails and the group foraging efficiency. Numerical simulations of the model indicate that depending on the systematic change of the feeding environment, the optimal distribution of stochasticity shifts from a mixture of almost deterministic and mildly stochastic ants to a contrasted mixture of almost deterministic ants and highly stochastic ants. In addition, the interaction between the stochasticity and the pheromone path regulates the dynamics of the foraging efficiency optimization. Stochasticity could strengthen the collective efficiency when the variance in the sensitivity to pheromone for ants is introduced in the model.