Congenital anaemia associated with loss-of-function variants in DNA polymerase epsilon 1.

DNA polymerase epsilon (Pol ε), a component of the core replisome, is involved in DNA replication. Although genetic defects of Pol ε have been reported to cause immunodeficiency syndromes, its role in haematopoiesis remains unknown. Here, we identified compound heterozygous variants (p.[Asp1131fs];[Thr1891del]) in POLE, encoding Pol ε catalytic subunit A (POLE1), in siblings with a syndromic form of severe congenital transfusion-dependent anaemia. In contrast to Diamond-Blackfan anaemia, marked reticulocytopenia or marked erythroid hypoplasia was not found. Their bone marrow aspirates during infancy revealed erythroid dysplasia with strongly positive TP53 in immunostaining. Repetitive examinations demonstrated trilineage myelodysplasia within 2 years from birth. They had short stature and facial dysmorphism. HEK293 cell-based expression experiments and analyses of patient-derived induced pluripotent stem cells (iPSCs) disclosed a reduced mRNA level of Asp1131fs-POLE1 and defective nuclear translocation of Thr1891del-POLE1. Analysis of iPSCs showed compensatory mRNA upregulation of the other replisome components and increase of the TP53 protein, both suggesting dysfunction of the replisome. We created Pole-knockout medaka fish and found that heterozygous fishes were viable, but with decreased RBCs. Our observations expand the phenotypic spectrum of the Pol ε defect in humans, additionally providing unique evidence linking Pol ε to haematopoiesis.

A resource of single-cell gene expression profiles in a planarian Dugesia japonica.

The freshwater planarian Dugesia japonica maintains an abundant heterogeneous cell population called neoblasts, which include adult pluripotent stem cells. Thus, it is an excellent model organism for stem cell and regeneration research. Recently, many single-cell RNA sequencing (scRNA-seq) databases of several model organisms, including other planarian species, have become publicly available; these are powerful and useful resources to search for gene expression in various tissues and cells. However, the only scRNA-seq dataset for D. japonica has been limited by the number of genes detected. Herein, we collected D. japonica cells, and conducted an scRNA-seq analysis. A novel, automatic, iterative cell clustering strategy produced a dataset of 3,404 cells, which could be classified into 63 cell types based on gene expression profiles. We introduced two examples for utilizing the scRNA-seq dataset in this study using D. japonica. First, the dataset provided results consistent with previous studies as well as novel functionally relevant insights, that is, the expression of DjMTA and DjP2X-A genes in neoblasts that give rise to differentiated cells. Second, we conducted an integrative analysis of the scRNA-seq dataset and time-course bulk RNA-seq of irradiated animals, demonstrating that the dataset can help interpret differentially expressed genes captured via bulk RNA-seq. Using the R package "Seurat" and GSE223927, researchers can easily access and utilize this dataset.

Advances in developmental biology using next-generation sequencing (NGS): Workshop for next-generation researchers.

This is a meeting report of "The workshop of research and techniques using next-generation sequencing (NGS) for developmental biology 2022." About 20 researchers attended the workshop. We discuss their NGS studies and techniques.

Fillable and unfillable gaps in plant transcriptome under field and controlled environments.

The differences between plants grown in field and in controlled environments have long been recognized. However, few studies have addressed the underlying molecular mechanisms. To evaluate plant responses to fluctuating environments using laboratory equipment, we developed SmartGC, a high-performance growth chamber that reproduces the fluctuating irradiance, temperature and humidity of field environments. We analysed massive transcriptome data of rice plants grown under field and SmartGC conditions to clarify the differences in plant responses to field and controlled environments. Rice transcriptome dynamics in SmartGC mimicked those in the field, particularly during the morning and evening but those in conventional growth chamber conditions did not. Further analysis revealed that fluctuation of irradiance affects transcriptome dynamics in the morning and evening, while fluctuation of temperature affects transcriptome dynamics only in the morning. We found upregulation of genes related to biotic and abiotic stress, and their expression was affected by environmental factors that cannot be mimicked by SmartGC. Our results reveal fillable and unfillable gaps in the transcriptomes of rice grown in field and controlled environments and can accelerate the understanding of plant responses to field environments for both basic biology and agricultural applications.

Genomic Basis of Transcriptome Dynamics in Rice under Field Conditions.

How genetic variations affect gene expression dynamics of field-grown plants remains unclear. Expression quantitative trait loci (eQTL) analysis is frequently used to find genomic regions underlying gene expression polymorphisms. This approach requires transcriptome data for the complete set of the QTL mapping population under the given conditions. Therefore, only a limited range of environmental conditions is covered by a conventional eQTL analysis. We sampled sparse time series of field-grown rice from chromosome segment substitution lines (CSSLs) and conducted RNA sequencing (RNA-Seq). Then, by using statistical analysis integrating meteorological data and the RNA-Seq data, we identified 1,675 eQTLs leading to polymorphisms in expression dynamics under field conditions. A genomic region on chromosome 11 influences the expression of several defense-related genes in a time-of-day- and scaled-age-dependent manner. This includes the eQTLs that possibly influence the time-of-day- and scaled-age-dependent differences in the innate immunity between Koshihikari and Takanari. Based on the eQTL and meteorological data, we successfully predicted gene expression under environments different from training environments and in rice cultivars with more complex genotypes than the CSSLs. Our novel approach of eQTL identification facilitated the understanding of the genetic architecture of expression dynamics under field conditions, which is difficult to assess by conventional eQTL studies. The prediction of expression based on eQTLs and environmental information could contribute to the understanding of plant traits under diverse field conditions.

What is the role of PIWI family proteins in adult pluripotent stem cells? Insights from asexually reproducing animals, planarians.

Planarians have a remarkable regenerative ability owing to their adult pluripotent stem cells (aPSCs), which are called "neoblasts." Planarians maintain a considerable number of neoblasts throughout their adulthood to supply differentiated cells for the maintenance of tissue homeostasis and asexual reproduction (fission followed by regeneration). Thus, planarians serve as a good model to study the regulatory mechanisms of in vivo aPSCs. In asexually reproducing invertebrates, such as sponge, Hydra, and planaria, piwi family genes are the markers most commonly expressed in aPSCs. While piwi family genes are known as guardians against transposable elements in the germline cells of animals that only sexually propagate, their functions in the aPSC system have remained elusive. In this review, we introduce recent knowledge on the PIWI family proteins in the aPSC system in planarians and other organisms and discuss how PIWI family proteins contribute to the regulation of the aPSC system.

Identification and characterization of a fibroblast growth factor gene in the planarian Dugesia japonica.

Planarians belong to the phylum Platyhelminthes and can regenerate their missing body parts after injury via activation of somatic pluripotent stem cells called neoblasts. Previous studies suggested that fibroblast growth factor (FGF) signaling plays a crucial role in the regulation of head tissue differentiation during planarian regeneration. To date, however, no FGF homologues in the Platyhelminthes have been reported. Here, we used a planarian Dugesia japonica model and identified an fgf gene termed Djfgf, which encodes a putative secreted protein with a core FGF domain characteristic of the FGF8/17/18 subfamily in bilaterians. Using Xenopus embryos, we found that DjFGF has FGF activity as assayed by Xbra induction. We next examined Djfgf expression in non-regenerating intact and regenerating planarians. In intact planarians, Djfgf was expressed in the auricles in the head and the pharynx. In the early process of regeneration, Djfgf was transiently expressed in a subset of differentiated cells around wounds. Notably, Djfgf expression was highly induced in the process of head regeneration when compared to that in the tail regeneration. Furthermore, assays of head regeneration from tail fragments revealed that combinatorial actions of the anterior extracellular signal-regulated kinase (ERK) and posterior Wnt/ß-catenin signaling restricted Djfgf expression to a certain anterior body part. This is the region where neoblasts undergo active proliferation to give rise to their differentiating progeny in response to wounding. The data suggest the possibility that DjFGF may act as an anterior counterpart of posteriorly localized Wnt molecules and trigger neoblast responses involved in planarian head regeneration.

High-yielding rice Takanari has superior photosynthetic response to a commercial rice Koshihikari under fluctuating light.

Leaves within crop canopies experience variable light over the course of a day, which greatly affects photosynthesis and crop productivity. Little is known about the mechanisms of the photosynthetic response to fluctuating light and their genetic control. Here, we examined gas exchange, metabolite levels, and chlorophyll fluorescence during the photosynthetic induction response in an Oryza sativa indica cultivar with high yield (Takanari) and a japonica cultivar with lower yield (Koshihikari). Takanari had a faster induction response to sudden increases in light intensity than Koshihikari, as demonstrated by faster increases in net CO2 assimilation rate, stomatal conductance, and electron transport rate. In a simulated light regime that mimicked a typical summer day, the faster induction response in Takanari increased daily CO2 assimilation by 10%. The faster response of Takanari was explained in part by its maintenance of a larger pool of Calvin-Benson cycle metabolites. Together, the rapid responses of electron transport rate, metabolic flux, and stomatal conductance in Takanari contributed to the greater daily carbon gain under fluctuating light typical of natural environments.

Searching for non-transposable targets of planarian nuclear PIWI in pluripotent stem cells and differentiated cells.

Nuclear PIWIs together with their guide RNAs (piRNAs) epigenetically silence various genes including transposons in many organisms. In planarians, the nuclear piwi family gene, DjpiwiB is specifically transcribed in adult pluripotent stem cells (adult PSC, neoblast), but not in differentiated cells. However, the protein accumulates in the nuclei of both neoblasts and their descendant differentiated cells. Interestingly, PIWI(DjPiwiB)-piRNA complexes are indispensable for the repression of transposable genes at the onset of differentiation from neoblasts. Here, we conducted a comparative transcriptome analysis between control and DjpiwiB(RNAi) animals to identify non-transposable target genes of the DjPiwiB-piRNA complexes. Using bioinformatic analyses and RNAi we demonstrate that DjPiwiB-piRNA complexes are required for the proper expression of Djmcm2 and Djhistone h4 in neoblasts and that DjPiwiB-piRNA complexes regulate the transient expression of Djcalu during neoblast differentiation. Thus, DjPiwiB-piRNA complexes regulate the correct expression patterns during neoblast self-renewal and differentiation.

Heterogeneity of chromatoid bodies in adult pluripotent stem cells of planarian Dugesia japonica.

The robust regenerative ability of planarians is known to be dependent on adult pluripotent stem cells called neoblasts. One of the morphological features of neoblasts is cytoplasmic ribonucleoprotein granules (chromatoid bodies: CBs), which resemble germ granules present in germline cells in other animals. Previously, we showed by immuno-electron microscopic analysis that DjCBC-1, a planarian Me31B/Dhh1/DDX6 homologue, which is a component of ribonucleoprotein granules, was localized in CBs in the planarian Dugesia japonica. Also, recently it was reported using another planarian species that Y12 antibody recognizing symmetrical dimethylarginine (sDMA) specifically binds to CBs in which histone mRNA is co-localized. Here, we showed by double immunostaining and RNA interference (RNAi) that DjCBC-1-containing CBs and Y12-immunoreactive CBs are distinct structures, suggesting that CBs are composed of heterogeneous populations. We also found that the Y12-immunoreactive CBs specifically contained a cytoplasmic type of planarian PIWI protein (DjPiwiC). We revealed by RNAi experiments that Y12-immunoreactive CBs may have anti-transposable element activity involving the DjPiwiC protein in the neoblasts.

Molecular markers for X-ray-insensitive differentiated cells in the Inner and outer regions of the mesenchymal space in planarian Dugesia japonica.

Planarian's strong regenerative ability is dependent on stem cells (called neoblasts) that are X-ray-sensitive and proliferative stem cells. In addition to neoblasts, another type of X-ray-sensitive cells was newly identified by recent research. Thus, planarian's X-ray-sensitive cells can be divided into at least two populations, Type 1 and Type 2, the latter corresponding to planarian's classically defined "neoblasts". Here, we show that Type 1 cells were distributed in the outer region (OR) immediately underneath the muscle layer at all axial levels from head to tail, while the Type 2 cells were distributed in a more internal region (IR) of the mesenchymal space at the axial levels from neck to tail. To elucidate the biological significance of these two regions, we searched for genes expressed in differentiated cells that were locate close to these X-ray-sensitive cell populations in the mesenchymal space, and identified six genes mainly expressed in the OR or IR, named OR1, OR2, OR3, IR1, IR2 and IR3. The predicted amino acid sequences of these genes suggested that differentiated cells expressing OR1, OR3, IR1, or IR2 provide Type 1 and Type 2 cells with specific extracellular matrix (ECM) environments.

Multiple Neuropeptide-Coding Genes Involved in Planarian Pharynx Extension.

Planarian feeding behavior involves three steps: moving toward food, extending the pharynx from their planarian's ventral side after arriving at the food, and ingesting the food through the pharynx. Although pharynx extension is a remarkable behavior, it remains unknown what neuronal cell types are involved in its regulation. To identify neurons involved in regulating pharynx extension, we quantitatively analyzed pharynx extension and sought to identify these neurons by RNA interference (RNAi) and in situ hybridization. This assay, when performed using planarians with amputation of various body parts, clearly showed that the head portion is indispensable for inducing pharynx extension. We thus tested the effects of knockdown of brain neurons such as serotonergic, GABAergic, and dopaminergic neurons by RNAi, but did not observe any effects on pharynx extension behavior. However, animals with RNAi of the Prohormone Convertase 2 (PC2, a neuropeptide processing enzyme) gene did not perform the pharynx extension behavior, suggesting the possible involvement of neuropeptide(s in the regulation of pharynx extension. We screened 24 neuropeptide-coding genes, analyzed their functions by RNAi using the pharynx extension assay system, and identified at least five neuropeptide genes involved in pharynx extension. These was expressed in different cells or neurons, and some of them were expressed in the brain, suggesting complex regulation of planarian feeding behavior by the nervous system.

Extremely low-frequency electromagnetic fields facilitate both osteoblast and osteoclast activity through Wnt/ß-catenin signaling in the zebrafish scale.

Electromagnetic fields (EMFs) have received widespread attention as effective, noninvasive, and safe therapies across a range of clinical applications for bone disorders. However, due to the various frequencies of devices, their effects on tissues/cells are vary, which has been a bottleneck in understanding the effects of EMFs on bone tissue. Here, we developed an in vivo model system using zebrafish scales to investigate the effects of extremely low-frequency EMFs (ELF-EMFs) on fracture healing. Exposure to 10 millitesla (mT) of ELF-EMFs at 60 Hz increased the number of both osteoblasts and osteoclasts in the fractured scale, whereas 3 or 30 mT did not. Gene expression analysis revealed that exposure to 10 mT ELF-EMFs upregulated wnt10b and Wnt target genes in the fractured scale. Moreover, ß-catenin expression was enhanced by ELF-EMFs predominantly at the fracture site of the zebrafish scale. Inhibition of Wnt/ß-catenin signaling by IWR-1-endo treatment reduced both osteoblasts and osteoclasts in the fractured scale exposed to ELF-EMFs. These results suggest that ELF-EMFs promote both osteoblast and osteoclast activity through activation of Wnt/ß-catenin signaling in fracture healing. Our data provide in vivo evidence that ELF-EMFs generated with a widely used commercial AC power supply have a facilitative effect on fracture healing.

The RNA-binding protein Msi2 regulates autophagy during myogenic differentiation.

Skeletal muscle development is a highly ordered process orchestrated transcriptionally by the myogenic regulatory factors. However, the downstream molecular mechanisms of myogenic regulatory factor functions in myogenesis are not fully understood. Here, we identified the RNA-binding protein Musashi2 (Msi2) as a myogenin target gene and a post-transcriptional regulator of myoblast differentiation. Msi2 knockdown in murine myoblasts blocked differentiation without affecting the expression of MyoD or myogenin. Msi2 overexpression was also sufficient to promote myoblast differentiation and myocyte fusion. Msi2 loss attenuated autophagosome formation via down-regulation of the autophagic protein MAPL1LC3/ATG8 (LC3) at the early phase of myoblast differentiation. Moreover, forced activation of autophagy effectively suppressed the differentiation defects incurred by Msi2 loss. Consistent with its functions in myoblasts in vitro, mice deficient for Msi2 exhibited smaller limb skeletal muscles, poorer exercise performance, and muscle fiber-type switching in vivo. Collectively, our study demonstrates that Msi2 is a novel regulator of mammalian myogenesis and establishes a new functional link between muscular development and autophagy regulation.

Establishment of a zebrafish inbred strain, M-AB, capable of regular breeding and genetic manipulation.

Inbred strains of organisms are genetically highly uniform and thus useful for life science research. We have previously reported the ongoing generation of the zebrafish IM strain from the India (IND) strain through full sib-pair mating for 16 generations. However, the IM fish laid a small number of offspring and had a short lifespan, implying the need for discreet care in breeding. Here, we report the subsequent establishment of IM strain as well as the generation of a new inbred zebrafish strain, Mishima-AB (M-AB). M-AB was derived from the *AB strain by full sib-pair mating for over 20 generations, which fulfills the general criterion for the establishment of an inbred strain. In contrast to the IM case, maintenance of the M-AB strain by sib-pair mating required almost no special handling. Genome sequencing of IM individuals from the 47th generation and M-AB individuals from the 27th generation revealed that SNP-based genomic heterogeneity across whole-genome nucleotides was 0.008% and 0.011%, respectively. These percentages were much lower than those of the parental IND (0.197%) and *AB (0.086%) strains. These results indicate that the genomes of these inbred strains were highly homogenous. We also demonstrated the successful microinjection of antisense morpholinos, CRISPR/Cas9, and foreign genes into M-AB embryos at the 1-cell stage. Overall, we report the establishment of a zebrafish inbred strain, M-AB, which is capable of regular breeding and genetic manipulation. This strain will be useful for the analysis of genetically susceptible phenotypes such as behaviors, microbiome features and drug susceptibility.

Activation of XBP1 but not ATF6α rescues heart failure induced by persistent ER stress in medaka fish.

The unfolded protein response is triggered in vertebrates by ubiquitously expressed IRE1α/ß (although IRE1ß is gut-specific in mice), PERK, and ATF6α/ß, transmembrane-type sensor proteins in the ER, to cope with ER stress, the accumulation of unfolded and misfolded proteins in the ER. Here, we burdened medaka fish, a vertebrate model organism, with ER stress persistently from fertilization by knocking out the AXER gene encoding an ATP/ADP exchanger in the ER membrane, leading to decreased ATP concentration-mediated impairment of the activity of Hsp70- and Hsp90-type molecular chaperones in the ER lumen. ER stress and apoptosis were evoked from 4 and 6 dpf, respectively, leading to the death of all AXER-KO medaka by 12 dpf because of heart failure (medaka hatch at 7 dpf). Importantly, constitutive activation of IRE1α signaling-but not ATF6α signaling-rescued this heart failure and allowed AXER-KO medaka to survive 3 d longer, likely because of XBP1-mediated transcriptional induction of ER-associated degradation components. Thus, activation of a specific pathway of the unfolded protein response can cure defects in a particular organ.

Planarian PIWI-piRNA Interaction Analysis Using Immunoprecipitation and piRNA Sequencing.

The freshwater planarian Dugesia japonica is a good in vivo model for studying the function of piwi genes in adult pluripotent stem cell (aPSC) due to their abundant aPSCs. Generally, PIWI family proteins encoded by piwi genes bind to small noncoding RNAs called piRNAs (PIWI-interacting piRNAs). The analysis of PIWI-piRNA complexes in the planarian is useful for revealing the functions of piwi genes in the aPSC system. In this chapter, we present an immunoprecipitation protocol for PIWI-piRNA complexes from whole planarians.

DeLTa-Seq: direct-lysate targeted RNA-Seq from crude tissue lysate.

BACKGROUND: Quantification of gene expression such as RNA-Seq is a popular approach to study various biological phenomena. Despite the development of RNA-Seq library preparation methods and sequencing platforms in the last decade, RNA extraction remains the most laborious and costly step in RNA-Seq of tissue samples of various organisms. Thus, it is still difficult to examine gene expression in thousands of samples. RESULTS: Here, we developed Direct-RT buffer in which homogenization of tissue samples and direct-lysate reverse transcription can be conducted without RNA purification. The DTT concentration in Direct-RT buffer prevented RNA degradation but not RT in the lysates of several plant tissues, yeast, and zebrafish larvae. Direct reverse transcription on these lysates in Direct-RT buffer produced comparable amounts of cDNA to those synthesized from purified RNA. To maximize the advantage of the Direct-RT buffer, we integrated Direct-RT and targeted RNA-Seq to develop a cost-effective, high-throughput quantification method for the expressions of hundreds of genes: DeLTa-Seq (Direct-Lysate reverse transcription and Targeted RNA-Seq). The DeLTa-Seq method could drastically improve the efficiency and accuracy of gene expression analysis. DeLTa-Seq analysis of 1056 samples revealed the temperature-dependent effects of jasmonic acid and salicylic acid in Arabidopsis thaliana. CONCLUSIONS: The DeLTa-Seq method can realize large-scale studies using thousands of animal, plant, and microorganism samples, such as chemical screening, field experiments, and studies focusing on individual variability. In addition, Direct-RT is also beneficial for gene expression analysis in small tissues from which it is difficult to purify enough RNA for the experiments.

Effect of differences in light source environment on transcriptome of leaf lettuce (Lactuca sativa L.) to optimize cultivation conditions.

When used in closed-type plant factories, light-emitting diode (LED) illumination systems have the particular advantages of low heat emission and high luminous efficiency. The effects of illumination quality and intensity on the growth and morphogenesis of many plant species have been examined, but improvements are needed to optimize the illumination systems for better plant products with lower resource investments. In particular, new strategies are needed to reduce the wastage of plant products related to leaf senescence, and to better control the ingredients and appearance of leafy vegetables. Although the quality of light is often altered to change the characteristics of plant products, the transcriptional status underlying the physiological responses of plants to light has not been established. Herein, we performed a comprehensive gene expression analysis using RNA-sequencing to determine how red, blue, and red/blue LEDs and fluorescent light sources affect transcriptome involved in the leaf aging of leaf lettuce. The RNA-sequencing profiling revealed clear differences in the transcriptome between young and old leaves. Red LED light caused large variation between the two age classes, while a pure or mixed blue LED light spectrum induced fewer transcriptome differences between young and old leaves. Collectively, the expression levels of genes that showed homology with those of other model organisms provide a detailed physiological overview, incorporating such characteristics as the senescence, nutrient deficiency, and anthocyanin synthesis of the leaf lettuce plants. Our findings suggest that transcriptome profiles of leaf lettuce grown under different light sources provide helpful information to achieve better growth conditions for marketable and efficient green-vegetable production, with improved wastage control and efficient nutrient inputs.

Direct-TRI: High-throughput RNA-extracting Method for All Stages of Zebrafish Development.

Recent popularization of next-generation sequencing enables conducting easy transcriptome analysis. Nevertheless, substantial RNA isolation work prior to RNA sequencing, as well as the high cost involved, still makes the routine use of large-scale transcriptome analysis difficult. For example, conventional phenol-chloroform RNA extraction cannot be easily applied to hundreds of samples. Therefore, we developed Direct-TRI, a new cost-effective and high throughput RNA-extraction method that uses a commercial guanidine-phenol-based RNA extraction reagent and a 96-well silica column plate. We applied Direct-TRI to zebrafish whole larvae and juvenile samples and obtained comparable RNA qualities by several different homogenization methods such as vortexing, manual homogenizing, and freezing/crushing. Direct-TRI enabled the extraction of 192 RNA samples in an hour with a cost of less than a dollar per sample. Direct-TRI is useful for large-scale transcriptome studies, manipulating hundreds of zebrafish individuals, and may be used with other animal samples.