Rickettsiae are obligate intracellular bacteria that have small genomes as a result of reductive evolution. Many Rickettsia species of the spotted fever group (SFG) cause tick-borne diseases known as "spotted fevers". The life cycle of SFG rickettsiae is closely associated with that of the tick, which is generally thought to act as a bacterial vector and reservoir that maintains the bacterium through transstadial and transovarial transmission. Each SFG member is thought to have adapted to a specific tick species, thus restricting the bacterial distribution to a relatively limited geographic region. These unique features of SFG rickettsiae allow investigation of how the genomes of such biologically and ecologically specialized bacteria evolve after genome reduction and the types of population structures that are generated. Here, we performed a nationwide, high-resolution phylogenetic analysis of Rickettsia japonica, an etiological agent of Japanese spotted fever that is distributed in Japan and Korea. The comparison of complete or nearly complete sequences obtained from 31 R. japonica strains isolated from various sources in Japan over the past 30 years demonstrated an extremely low level of genomic diversity. In particular, only 34 single nucleotide polymorphisms were identified among the 27 strains of the major lineage containing all clinical isolates and tick isolates from the three tick species. Our data provide novel insights into the biology and genome evolution of R. japonica, including the possibilities of recent clonal expansion and a long generation time in nature due to the long dormant phase associated with tick life cycles.

Recent studies on paleogenomics have reported some Paleolithic and Neolithic genomes that have provided new insights into the human population history in East and Northeast Asia. However, there remain some cases where more recent migration events need to be examined to elucidate the detailed formation process of local populations. Although the area around northern Japan is one of the regions archaeologically suggested to have been affected by migration waves after the Neolithic period, the genetic source of these migrations are still unclear. Thus, genomic data from such past migrant populations would be highly informative to clarify the detailed formation process of local populations in this region. Here, we report the genome sequence of a 900-year-old adult female (NAT002) belonging to the prehistoric Okhotsk people, who have been considered to be the past migrants to northern Japan after the Neolithic period. We found a close relationship between NAT002 and modern Lower Amur populations and past admixture events between the Amur, Jomon, and Kamchatka ancestries. The admixture dating suggested migration of Amur-related ancestry at approximately 1,600 BP, which is compatible with the archaeological evidence regarding the settlement of the Okhotsk people. Our results also imply migration of Kamchatka-related ancestry at approximately 2,000 BP. In addition, human leukocyte antigen (HLA) typing detected the HLA-B*40 allele, which is reported to increase the risk of arthritis, suggesting the genetic vulnerability of NAT002 to hyperostosis, which was observed around her chest clavicle.

We report the complete genome sequence of Staphylococcus aureus strain JP025, which was isolated from a furunculosis sample from a Japanese patient. The strain carried two hyaluronate lyase genes, JP025hysA and JP025hysB, on the chromosome and was classified as sequence type 30.

Heterothallic strains of the Closterium peracerosum-strigosum-littorale (C. psl.) complex have two sexes, mating-type plus (mt+) and mating-type minus (mt-). Conjugation between these two sexes is regulated by two sex pheromones, protoplast-release-inducing protein (PR-IP) and PR-IP Inducer, which are produced by mt+ and mt- cells, respectively. PR-IP mediates the release of protoplasts from mt- cells during mating. In this study, we examined the mechanism of action of CpRLP1 (receptor-like protein 1), which was previously identified in a cDNA microarray analysis as one of the PR-IP-inducible genes. Using CRISPR/Cas9 technology, we generated CpRLP1 knockout mutants in mt- cells of the C. psl. complex. When the knockout mt- cells were mixed with wild-type mt+ cells, conjugation was severely reduced. Many cells released protoplasts without pairing, suggesting a loss of synchronization between the two mating partners. Furthermore, the knockout mutants were hypersensitive to PR-IP. We conclude that CpRLP1 is a negative regulator of PR-IP that regulates the timing of protoplast release in conjugating C. psl. cells. As the first report of successful gene knockout in the class Charophyceae, this study provides a basis for research aimed at understanding the ancestral roles of genes that are indispensable for the development of land plants.

Cryptomeria japonica D. Don is one of the most commercially important conifers in Japan. However, the allergic disease caused by its pollen is a severe public health problem in Japan. Since large-scale analysis of expressed sequence tags (ESTs) in the male strobili of C. japonica should help us to clarify the overall expression of genes during the process of pollen development, we constructed a full-length enriched cDNA library that was derived from male strobili at various developmental stages.

Sex chromosomes turn over rapidly in some taxonomic groups, where closely related species have different sex chromosomes. Although there are many examples of sex chromosome turnover, we know little about the functional roles of sex chromosome turnover in phenotypic diversification and genomic evolution. The sympatric pair of Japanese threespine stickleback (Gasterosteus aculeatus) provides an excellent system to address these questions: the Japan Sea species has a neo-sex chromosome system resulting from a fusion between an ancestral Y chromosome and an autosome, while the sympatric Pacific Ocean species has a simple XY sex chromosome system. Furthermore, previous quantitative trait locus (QTL) mapping demonstrated that the Japan Sea neo-X chromosome contributes to phenotypic divergence and reproductive isolation between these sympatric species. To investigate the genomic basis for the accumulation of genes important for speciation on the neo-X chromosome, we conducted whole genome sequencing of males and females of both the Japan Sea and the Pacific Ocean species. No substantial degeneration has yet occurred on the neo-Y chromosome, but the nucleotide sequence of the neo-X and the neo-Y has started to diverge, particularly at regions near the fusion. The neo-sex chromosomes also harbor an excess of genes with sex-biased expression. Furthermore, genes on the neo-X chromosome showed higher non-synonymous substitution rates than autosomal genes in the Japan Sea lineage. Genomic regions of higher sequence divergence between species, genes with divergent expression between species, and QTL for inter-species phenotypic differences were found not only at the regions near the fusion site, but also at other regions along the neo-X chromosome. Neo-sex chromosomes can therefore accumulate substitutions causing species differences even in the absence of substantial neo-Y degeneration.

Caenorhabditis elegans mechanoreceptors located in ASG sensory neurons have been found to sense ambient temperature, which is a key trait for animal survival. Here, we show that experimental loss of xanthine dehydrogenase (XDH-1) function in AIN and AVJ interneurons results in reduced cold tolerance and atypical neuronal response to changes in temperature. These interneurons connect with upstream neurons such as the mechanoreceptor-expressing ASG. Ca2+ imaging revealed that ASG neurons respond to warm temperature via the mechanoreceptor DEG-1, a degenerin/epithelial Na+ channel (DEG/ENaC), which in turn affects downstream AIN and AVJ circuits. Ectopic expression of DEG-1 in the ASE gustatory neuron results in the acquisition of warm sensitivity, while electrophysiological analysis revealed that DEG-1 and human MDEG1 were involved in warm sensation. Taken together, these results suggest that cold tolerance is regulated by mechanoreceptor-mediated circuit calculation.

Invasive species pose a major threat to biological diversity. Although introduced populations often experience population bottlenecks, some invasive species are thought to be originated from hybridization between multiple populations or species, which can contribute to the maintenance of high genetic diversity. Recent advances in genome sequencing enable us to trace the evolutionary history of invasive species even at whole-genome level and may help to identify the history of past hybridization that may be overlooked by traditional marker-based analysis. Here, we conducted whole-genome sequencing of eight threespine stickleback (Gasterosteus aculeatus) individuals, four from a recently introduced crater lake population and four of the putative source population. We found that both populations have several small genomic regions with high genetic diversity, which resulted from introgression from a closely related species (Gasterosteus nipponicus). The sizes of the regions were too small to be detected with traditional marker-based analysis or even some reduced-representation sequencing methods. Further amplicon sequencing revealed linkage disequilibrium around an introgression site, which suggests the possibility of selective sweep at the introgression site. Thus, interspecies introgression might predate introduction and increase genetic variation in the source population. Whole-genome sequencing of even a small number of individuals can therefore provide higher resolution inference of history of introduced populations.

Depression is associated with taste disorders; however, the mechanisms by which mental stress affects taste perception are not well understood. This study aimed to elucidate the effects of psychosocial stress on peripheral taste-sensing systems using a mouse depression model. Male mice were subjected to subchronic and mild social defeat stress (sCSDS). Results showed that sCSDS significantly increased body weight, food and water intake, and social avoidance behavior and that sCSDS did not change reward-seeking behavior on sucrose preference but tended to decrease pheromonal preference for female urine. Furthermore, sCSDS downregulated the mRNA levels of sweet and umami taste receptor subunits, i.e., sweet taste receptor type 1 members 2 and 3 (T1R2 and T1R3), but not the umami taste receptor subunit, i.e., taste receptor type 1 member 1 (T1R1), in the circumvallate papillae of mice. It is known that sucrose preference is mediated by the gut-brain axis without taste perception; thus, it was considered that sCSDS affected the peripheral taste-sensing systems, rather than the central reward systems, which mediate sucrose preference. This is the first study to report that psychosocial stress affects peripheral sweet and umami taste-sensing systems.

Contact domains of chromatin serve as a fundamental unit to regulate action of enhancers for target genes. Looping between a pair of CCCTC-binding factor (CTCF)-binding sites in convergent orientations underlies the formation of contact domains, while those in divergent orientations establish domain boundaries. However, every CTCF site is not necessarily engaged in loop or boundary structures, leaving functions of CTCF in varied genomic contexts still elusive. The locus containing Tfap2c and Bmp7 encompasses two contact domains separated by a region between the two genes, termed transition zone (TZ), characterized by two arrays of CTCF sites in divergent configuration. In this study, we created deletion and inversion alleles of these and other regions across the locus and investigated how they impinge on the conformation.

Recent advances in CRISPR technology have enabled us to perform gene knock-in in various species and cell lines. CRISPR-mediated knock-in requires donor DNA which serves as a template for homology-directed repair (HDR). For knock-in of short sequences or base substitutions, ssDNA donors are frequently used among various other forms of HDR donors, such as linear dsDNA. However, partly due to the complexity of long ssDNA preparation, it remains unclear whether ssDNA is the optimal type of HDR donors for insertion of long transgenes such as fluorescent reporters in human cells.

With the rapid advances in next-generation sequencing (NGS), datasets for DNA polymorphisms among various species and strains have been produced, stored, and distributed. However, reliability varies among these datasets because the experimental and analytical conditions used differ among assays. Furthermore, such datasets have been frequently distributed from the websites of individual sequencing projects. It is desirable to integrate DNA polymorphism data into one database featuring uniform quality control that is distributed from a single platform at a single place. DNA polymorphism annotation database (DNApod; http://tga.nig.ac.jp/dnapod/) is an integrated database that stores genome-wide DNA polymorphism datasets acquired under uniform analytical conditions, and this includes uniformity in the quality of the raw data, the reference genome version, and evaluation algorithms. DNApod genotypic data are re-analyzed whole-genome shotgun datasets extracted from sequence read archives, and DNApod distributes genome-wide DNA polymorphism datasets and known-gene annotations for each DNA polymorphism. This new database was developed for storing genome-wide DNA polymorphism datasets of plants, with crops being the first priority. Here, we describe our analyzed data for 679, 404, and 66 strains of rice, maize, and sorghum, respectively. The analytical methods are available as a DNApod workflow in an NGS annotation system of the DNA Data Bank of Japan and a virtual machine image. Furthermore, DNApod provides tables of links of identifiers between DNApod genotypic data and public phenotypic data. To advance the sharing of organism knowledge, DNApod offers basic and ubiquitous functions for multiple alignment and phylogenetic tree construction by using orthologous gene information.

While regulation of gene-enhancer interaction is intensively studied, its application remains limited. Here, we reconstituted arrays of CTCF-binding sites and devised a synthetic topological insulator with tetO for chromatin-engineering (STITCH). By coupling STITCH with tetR linked to the KRAB domain to induce heterochromatin and disable the insulation, we developed a drug-inducible system to control gene activation by enhancers. In human induced pluripotent stem cells, STITCH inserted between MYC and the enhancer down-regulated MYC. Progressive mutagenesis of STITCH led to a preferential escalation of the gene-enhancer interaction, corroborating the strong insulation ability of STITCH. STITCH also altered epigenetic states around MYC. Time-course analysis by drug induction uncovered deposition and removal of H3K27me3 repressive marks follows and reflects, but does not precede and determine, the expression change. Finally, STITCH inserted near NEUROG2 impaired the gene activation in differentiating neural progenitor cells. Thus, STITCH should be broadly useful for functional genetic studies.

Depressive disorders are partially attributed to chronic inflammation associated with the tryptophan (Trp)-kynurenine (Kyn) pathway. Recent evidence suggests that anti-inflammatory agents may reduce the risk of depression. The present study aimed to elucidate the potential of the citrus flavonoid hesperidin, which exhibits anti-inflammatory activity, in suppressing the Trp-Kyn pathway in the brain, using a lipopolysaccharide (LPS)-induced inflammation mouse model. Dietary hesperidin was found to suppress activation of the Trp-Kyn pathway in the prefrontal cortex. In addition, it reduced systemic LPS-induced signs of illness, such as low skin temperature and enhanced leukocyte count in the blood. However, dietary supplementation with hesperidin did not improve body weight loss, food intake, water intake, or splenic increases in leukocyte numbers in the LPS model. Collectively, the results suggest that dietary hesperidin can partially regulate central and peripheral events linked to inflammation in LPS mouse models.

Three coenzyme A (CoA) molecular species, i.e., acetyl-CoA, malonyl-CoA, and nonesterified CoA (CoASH), in 13 types of fasted rat tissue were analyzed. A relatively larger pool size of total CoA, consisting of acetyl-CoA, malonyl-CoA, and CoASH, was observed in the medulla oblongata, liver, heart, and brown adipose tissue. Focusing on changes in the CoA pool size in response to the nutrient composition of the diet given, total CoA pools in rats continuously fed a high-fat diet for 4 weeks were significantly higher in the hypothalamus, cerebellum, and kidney, and significantly lower in the liver and skeletal muscle than those of rats fed a high-carbohydrate or high-protein diet. In particular, reductions in the liver were remarkable and were caused by decreased CoASH levels. Consequently, the total CoA pool size was reduced by approximately one-fifth of the hepatic contents of rats fed the other diets. In the hypothalamus, which monitors energy balance, all three CoA molecular species measured were at higher levels when rats were fed the high-fat diet. Thus, it was of interest that feeding rats a high-fat diet affected the behaviors of CoA pools in the hypothalamus, liver, and skeletal muscle, suggesting a significant relationship between CoA pools, especially malonyl-CoA and/or CoASH pools, and lipid metabolism in vivo.

Speciation is a continuous process and analysis of species pairs at different stages of divergence provides insight into how it unfolds. Previous genomic studies on young species pairs have revealed peaks of divergence and heterogeneous genomic differentiation. Yet less known is how localised peaks of differentiation progress to genome-wide divergence during the later stages of speciation in the presence of persistent gene flow. Spanning the speciation continuum, stickleback species pairs are ideal for investigating how genomic divergence builds up during speciation. However, attention has largely focused on young postglacial species pairs, with little knowledge of the genomic signatures of divergence and introgression in older stickleback systems. The Japanese stickleback species pair, composed of the Pacific Ocean three-spined stickleback (Gasterosteus aculeatus) and the Japan Sea stickleback (G. nipponicus), which co-occur in the Japanese islands, is at a late stage of speciation. Divergence likely started well before the end of the last glacial period and crosses between Japan Sea females and Pacific Ocean males result in hybrid male sterility. Here we use coalescent analyses and Approximate Bayesian Computation to show that the two species split approximately 0.68-1 million years ago but that they have continued to exchange genes at a low rate throughout divergence. Population genomic data revealed that, despite gene flow, a high level of genomic differentiation is maintained across the majority of the genome. However, we identified multiple, small regions of introgression, occurring mainly in areas of low recombination rate. Our results demonstrate that a high level of genome-wide divergence can establish in the face of persistent introgression and that gene flow can be localized to small genomic regions at the later stages of speciation with gene flow.

ATP-dependent chromatin remodeling SWI/SNF complexes exist in three subcomplexes: canonical BAF (cBAF), polybromo BAF (PBAF), and a newly described non-canonical BAF (ncBAF). While cBAF and PBAF regulate fates of multiple cell types, roles for ncBAF in hematopoietic stem cells (HSCs) have not been investigated. Motivated by recent discovery of disrupted expression of BRD9, an essential component of ncBAF, in multiple cancers, including clonal hematopoietic disorders, we evaluate here the role of BRD9 in normal and malignant HSCs. BRD9 loss enhances chromatin accessibility, promoting myeloid lineage skewing while impairing B cell development. BRD9 significantly colocalizes with CTCF, whose chromatin recruitment is augmented by BRD9 loss, leading to altered chromatin state and expression of myeloid-related genes within intact topologically associating domains. These data uncover ncBAF as critical for cell fate specification in HSCs via three-dimensional regulation of gene expression and illuminate roles for ncBAF in normal and malignant hematopoiesis.

Herein, we present the first high-quality long-read-based chromosome-level genome assemblies and gene annotations of the genomes of three endangered Tokudaia species: Tokudaia osimensis, Tokudaia tokunoshimensis, and Tokudaia muenninki. These species, which are endemic to different islands of the Ryukyu Islands, Japan, exhibited unique karyotypes and sex chromosomal characteristics. The genome assemblies generated using PacBio, Illumina, and Hi-C sequence data consisted of 13 (corresponded to 12 autosomes and one X chromosome), 23 (corresponded to 22 autosomes and one X chromosome), and 23 (corresponded to 21 autosomes and the neo- and ancestral X regions) chromosome-level scaffolds that contained 2,445, 2,477, and 2,661 Mbp of sequence data, respectively. Annotations of protein-coding genes were performed using RNA-Seq-based, homology-based, and Ab initio methods. BUSCO completeness values for every species exceeded 96% for genomes and 98% for genes. These data can be an important resource for contributing to our understanding of species genomes resulting from allopatric speciation and provide insights into mammalian sex-determination mechanisms and sex chromosome evolution.

Zic3 encodes a zinc finger protein essential for the development of meso-ectodermal tissues. In mammals, Zic3 has important roles in the development of neural tube, axial skeletons, left-right body axis, and in maintaining pluripotency of ES cells. Here we characterized cis-regulatory elements required for Zic3 expression. Enhancer activities of human-chicken-conserved noncoding sequences around Zic1 and Zic3 were screened using chick whole-embryo electroporation. We identified enhancers for meso-ectodermal tissues. Among them, a mesodermal enhancer (Zic3-ME) in distant 3' flanking showed robust enhancement of reporter gene expression in the mesodermal tissue of chicken and mouse embryos, and was required for mesodermal Zic3 expression in mice. Zic3-ME minimal core region is included in the DNase hypersensitive region of ES cells, mesoderm, and neural progenitors, and was bound by T (Brachyury), Eomes, Lef1, Nanog, Oct4, and Zic2. Zic3-ME is derived from an ancestral sequence shared with a sequence encoding a mitochondrial enzyme. These results indicate that Zic3-ME is an integrated cis-regulatory element essential for the proper expression of Zic3 in vertebrates, serving as a hub for a gene regulatory network including Zic3.

With the accumulation of knowledge on the relation between psychological stress and gut microbiota, there is growing interest in controlling stress and/or mood disorders via probiotic supplementation. We aimed to examine the effect of probiotic Bifidobacterium bifidum TMC3115 (TMC3115) supplementation using a sub-chronic and mild social defeat stress murine model in this study. TM3115 supplementation maintained body weight gain and alleviated a polydipsia-like symptom induced by the stress. In the analyses of fecal and cecal bacterial profiles, expansions of Proteobacteria in stressed mice and increases in Actinobacteria and Bifidobacterium in mice supplemented with TMC3115 were observed. There was no marked difference in the diversity of cecal bacteria between the tested mice. Elevated serum levels of inflammatory markers such as tumor necrosis factor (TNF)-α and interleukin (IL)-6 were observed in the stressed mice, while TMC3115 only reduced the IL-6 level. These findings suggest that TMC3115 supplementation confers tolerance to psychosocial stress in the host through modulation of the gut microbiota and alleviation of stress-induced inflammatory responses. Furthermore, it may be expected to exert prevention and treatment of disorders related to peripheral IL-6, including depression.