Search

In most cases of sporadic colorectal cancers, tumorigenesis is a multistep process, involving genomic alterations in parallel with morphologic changes. In addition, accumulating evidence suggests that the human gut microbiome is linked to the development of colorectal cancer. Here we performed fecal metagenomic and metabolomic studies on samples from a large cohort of 616 participants who underwent colonoscopy to assess taxonomic and functional characteristics of gut microbiota and metabolites. Microbiome and metabolome shifts were apparent in cases of multiple polypoid adenomas and intramucosal carcinomas, in addition to more advanced lesions. We found two distinct patterns of microbiome elevations. First, the relative abundance of Fusobacterium nucleatum spp. was significantly (P

Colonization of new ecological niches has triggered large adaptive radiations. Although some lineages have made use of such opportunities, not all do so. The factors causing this variation among lineages are largely unknown. Here, we show that deficiency in docosahexaenoic acid (DHA), an essential w-3 fatty acid, can constrain freshwater colonization by marine fishes. Our genomic analyses revealed multiple independent duplications of the fatty acid desaturase gene Fads2 in stickleback lineages that subsequently colonized and radiated in freshwater habitats, but not in close relatives that failed to colonize. Transgenic manipulation of Fads2 in marine stickleback increased their ability to synthesize DHA and survive on DHA-deficient diets. Multiple freshwater ray-finned fishes also show a convergent increase in Fads2 copies, indicating its key role in freshwater colonization., This project was supported by JSPS KAKENHI grants 15H02418 and 23113007 (to J.K.), grant 16H06279 (to A.T.), and grants 26870824 and 16K07469 (to A.I.); the Asahi Glass Foundation (to J.K.); the Sumitomo Foundation (to A.I.); and JSPS postdoctoral fellowships 11J04816 to A.I. and 16J06812 to N.K. The work by N.L. and C.J.S. was funded by internal Eawag funds. B.M. was supported by SNF grant 31003A 175614. O.S. was supported by SNSF grant 31003A_163338. D.M. was supported by SNSF grant PDAMP3_123135 (to O.S).

In 2008, the genome assembly and gene models for the domestic silkworm, Bombyx mori, were published by a Japanese and Chinese collaboration group. However, the genome assembly contains a non-negligible number of misassembled and gap regions due to the presence of many repetitive sequences within the silkworm genome. The erroneous genome assembly occasionally causes incorrect gene prediction. Here we performed hybrid assembly based on 140 × deep sequencing of long (PacBio) and short (Illumina) reads. The remaining gaps in the initial genome assembly were closed using BAC and Fosmid sequences, giving a new total length of 460.3 Mb, with 30 gap regions and an N50 comprising 16.8 Mb in scaffolds and 12.2 Mb in contigs. More RNA-seq and piRNA-seq reads were mapped on the new genome assembly compared with the previous version, indicating that the new genome assembly covers more transcribed regions, including repetitive elements. We performed gene prediction based on the new genome assembly using available mRNA and protein sequence data. The number of gene models was 16,880 with an N50 of 2154 bp. The new gene models reflected more accurate coding sequences and gene sets than old ones. The proportion of repetitive elements was also reestimated using the new genome assembly, and was calculated to be 46.8% in the silkworm genome. The new genome assembly and gene models are provided in SilkBase ( http://silkbase.ab.a.u-tokyo.ac.jp ).

Identification of early biomarkers of stress is important for preventing mood and anxiety disorders. Saliva is an easy-to-collect and non-invasive diagnostic target. The aim of this study was to characterize the changes in salivary whole microRNAs (miRNAs) and metabolites in mice subjected to subchronic and mild social defeat stress (sCSDS). In this study, we identified seven upregulated and one downregulated miRNAs/PIWI-interacting RNA (piRNA) in the saliva of sCSDS mice. One of them, miR-208b-3p, which is reported as a reliable marker for myocardial infarction, was upregulated in the saliva of sCSDS mice. Histological analysis showed frequent myocardial interstitial fibrosis in the heart of such mice. In addition, gene ontology and pathway analyses suggested that the pathways related to energy metabolism, such as the oxidative phosphorylation and the pentose phosphate pathway, were significantly related to the miRNAs affected by sCSDS in saliva. In contrast, salivary metabolites were not significantly changed in the sCSDS mice, which is consistent with our previous metabolomic study on the plasma of sCSDS mice. Taken in the light of previous studies, the present study provides novel potential stress biomarkers for future diagnosis using saliva.

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., 論文

The evolution of the functionality of genes and genetic systems is a major source of animal diversity. Its best example is insect sex differentiation systems: promoting male and female differentiation (dual-functionality) or only male differentiation (single-functionality). However, the evolutionary origin of such functional diversity is largely unknown. Here, we investigate the ancestral functions ofdoublesex, a key factor of insect sex differentiation system, using the apterygote insect,Thermobia domestica, and reveal that itsdoublesexis essential for only males at the phenotypic level, but contributes to promoting female-specificvitellogeninexpression in females. This functional discordance between the phenotypic and transcription-regulatory levels inT. domesticashows a new type of functionality of animal sex differentiation systems. Then, we examine how the sex differentiation system transited from the single-functionality to the dual-functionality in phenotypes and uncover that a conserved female-specific motif ofdoublesexis detected in taxa with the dual-functionaldoublesex. It is estimated that the role of the sex differentiation system for female phenotypes may have evolved through accumulating mutations in the protein motif structures that led to the enhancement of its transcription-regulatory function.

The crown-of-thorns starfish (COTS) is a coral predator that is widely distributed in Indo-Pacific Oceans. A previous phylogenetic study using partial mitochondrial sequences suggested that COTS had diverged into four distinct species, but a nuclear genome-based analysis to confirm this was not conducted. To address this, COTS species nuclear genome sequences were analysed here, sequencing Northern Indian Ocean (NIO) and Red Sea (RS) species genomes for the first time, followed by a comparative analysis with the Pacific Ocean (PO) species. Phylogenetic analysis and ADMIXTURE analysis revealed clear divergences between the three COTS species. Furthermore, within the PO species, the phylogenetic position of the Hawaiian sample was further away from the other Pacific-derived samples than expected based on the mitochondrial data, suggesting that it may be a PO subspecies. The pairwise sequentially Markovian coalescent model showed that the trajectories of the population size diverged by region during the Mid-Pleistocene transition when the sea-level was dramatically decreased, strongly suggesting that the three COTS species experienced allopatric speciation. Analysis of the orthologues indicated that there were remarkable genes with species-specific positive selection in the genomes of the PO and RS species, which suggested that there may be local adaptations in the COTS species.

SummaryTermites are model social organisms characterized by a polyphenic caste system. Subterranean termites (Rhinotermitidae) are ecologically and economically important species, including acting as destructive pests. Rhinotermitidae occupies an important evolutionary position within the clade representing an intermediate taxon between the higher (Termitidae) and lower (other families) termites. Here, we report the genome, transcriptome and methylome of the Japanese subterranean termite Reticulitermes speratus. The analyses highlight the significance of gene duplication in social evolution in this termite. Gene duplication associated with caste-biased gene expression is prevalent in the R. speratus genome. Such duplicated genes encompass diverse categories related to social functions, including lipocalins (chemical communication), cellulases (wood digestion and social interaction), lysozymes (social immunity), geranylgeranyl diphosphate synthase (social defense) and a novel class of termite lineage-specific genes with unknown functions. Paralogous genes were often observed in tandem in the genome, but the expression patterns were highly variable, exhibiting caste biases. Some duplicated genes assayed were expressed in caste-specific organs, such as the accessory glands of the queen ovary and frontal glands in soldier heads. We propose that gene duplication facilitates social evolution through regulatory diversification leading to caste-biased expression and subfunctionalization and/or neofunctionalization that confers caste-specialized functions.Significance StatementTermites are model social organisms characterized by a sophisticated caste system, where distinct castes arise from the same genome. Our genomics data of Japanese subterranean termite provides insights into the evolution of the social system, highlighting the significance of gene duplication. Gene duplication associated with caste-biased gene expression is prevalent in the termite genome. Many of the duplicated genes were related to social functions, such as chemical communication, social immunity and defense, and they often expressed in caste-specific organs. We propose that gene duplication facilitates social evolution through regulatory diversification leading to caste-biased expression and functional specialization. In addition, since subterranean termites are ecologically and economically important species including destructive pests in the world, our genomics data serves as a foundation for these studies.

The level of γδ T cells varies between species because of differences in cellular dynamics and signals from thymic niche., αβ and γδ T cells are two distinct sublineages that develop in the vertebrate thymus. Thus far, their differentiation from a common progenitor is mostly understood to be regulated by intrinsic mechanisms. However, the proportion of αβ/γδ T cells varies in different vertebrate taxa. How this process is regulated in species that tend to produce a high frequency of γδ T cells is unstudied. Using an in vivo teleost model, the medaka, we report that progenitors first enter a thymic niche where their development into γδ T cells is favored. Translocation from this niche, mediated by chemokine receptor Ccr9b, is a prerequisite for their differentiation into αβ T cells. On the other hand, the thymic niche also generates opposing gradients of the cytokine interleukin-7 and chemokine Ccl25a, and, together, they influence the lineage outcome. We propose a previously unknown mechanism that determines the proportion of αβ/γδ lineages within species.

Laboratory mouse strains have mosaic genomes derived from at least three major subspecies that are distributed in Eurasia. Here, we describe genomic variations in ten inbred strains: Mus musculus musculus-derived BLG2/Ms, NJL/Ms, CHD/Ms, SWN/Ms, and KJR/Ms; M. m. domesticus-derived PGN2/Ms and BFM/Ms; M. m. castaneus-derived HMI/Ms; and JF1/Ms and MSM/Ms, which were derived from a hybrid between M. m. musculus and M. m. castaneus. These strains were established by Prof. Moriwaki in the 1980s and are collectively named the “Mishima Battery”. These strains show large phenotypic variations in body size and in many physiological traits. We resequenced the genomes of the Mishima Battery strains and performed a comparative genomic analysis with dbSNP data. More than 81 million nucleotide coordinates were identified as variant sites due to the large genetic distances among the mouse subspecies; 8,062,070 new SNP sites were detected in this study, and these may underlie the large phenotypic diversity observed in the Mishima Battery. The new information was collected in a reconstructed genome database, termed MoG+ that includes new application software and viewers. MoG+ intuitively visualizes nucleotide variants in genes and intergenic regions, and amino acid substitutions across the three mouse subspecies. We report statistical data from the resequencing and comparative genomic analyses and newly collected phenotype data of the Mishima Battery, and provide a brief description of the functions of MoG+, which provides a searchable and unique data resource of the numerous genomic variations across the three mouse subspecies. The data in MoG+ will be invaluable for research into phenotype-genotype links in diverse mouse strains.

Significance Evolutionary transitions between species with separate sexes and species in which individuals have both sex functions have wide-ranging biological implications. It is largely unknown how such transitions occur in systems with haploid male- and female-determining chromosomes in algae and bryophytes. We investigated such a transition in the algal genus Volvox by making whole-genome sequences of two closely related species, one of which is heterothallic (with distinct males and females) and the other homothallic (with only bisexual, self-compatible individuals). The heterothallic species harbors a sex-determining region (SDR), while the homothallic species retains a nearly intact female-derived SDR-like region and separate regions containing key male genes. Thus, an ancestral female has probably become homothallic by acquiring genes that confer male functions., Transitions between separate sexes (dioecy) and other mating systems are common across eukaryotes. Here, we study a change in a haploid dioecious green algal species with male- and female-determining chromosomes (U and V). The genus Volvox is an oogamous (with large, immotile female gametes and small, motile male gametes) and includes both heterothallic species (with distinct male and female genotypes, associated with a mating-type system that prevents fusion of gametes of the same sex) and homothallic species (bisexual, with the ability to self-fertilize). We date the origin of an expanded sex-determining region (SDR) in Volvox to at least 75 Mya, suggesting that homothallism represents a breakdown of dioecy (heterothallism). We investigated the involvement of the SDR of the U and V chromosomes in this transition. Using de novo whole-genome sequences, we identified a heteromorphic SDR of ca 1 Mbp in male and female genotypes of the heterothallic species Volvox reticuliferus and a homologous region (SDLR) in the closely related homothallic species Volvox africanus, which retained several different hallmark features of an SDR. The V. africanus SDLR includes a large region resembling the female SDR of the presumptive heterothallic ancestor, whereas most genes from the male SDR are absent. However, we found a multicopy array of the male-determining gene, MID, in a different genomic location from the SDLR. Thus, in V. africanus, an ancestrally female genotype may have acquired MID and thereby gained male traits.

Serpentinization of ultramafic rocks provides molecular hydrogen (H

Primordial microorganisms are postulated to have emerged in H2-rich alkaline Hadean serpentinite-hosted environments with homoacetogenesis as a core metabolism. However, investigation of two modern serpentinization-active analogues of early Earth reveals that conventional H2-/CO2-dependent homoacetogenesis is thermodynamically unfavorable in situ due to picomolar CO2 levels. Through metagenomics and thermodynamics, we discover unique taxa capable of metabolism adapted to the habitat. This included a novel deep-branching phylum, “Ca. Lithoacetigenota”, that exclusively inhabits Hadean analogues and harbors genes encoding alternative modes of H2-utilizing lithotrophy. Rather than CO2, these metabolisms utilize reduced carbon compounds detected in situ presumably serpentinization-derived: formate and glycine. The former employs a partial homoacetogenesis pathway and the latter a distinct pathway mediated by a rare selenoprotein – the glycine reductase. A survey of serpentinite-hosted system microbiomes shows that glycine reductases are diverse and nearly ubiquitous in Hadean analogues. “Ca. Lithoacetigenota” glycine reductases represent a basal lineage, suggesting that catabolic glycine reduction is an ancient bacterial innovation for gaining energy from geogenic H2 even under serpentinization-associated hyperalkaline, CO2-poor conditions. This draws remarkable parallels with ancestral archaeal H2-driven methyl-reducing methanogenesis recently proposed. Unique non-CO2-reducing metabolic strategies presented here may provide a new view into metabolisms that supported primordial life and the diversification of LUCA towards Archaea and Bacteria.

An obligately anaerobic, Gram-stain-negative, rod-shaped bacterium, designated strain 2CBH44T , was isolated from the fecal sample of a healthy Japanese man. This strain was initially assigned as a novel species of the genus Coprobacter based on the 16S rRNA gene sequence similarities compared with other Coprobacter species. The 16S rRNA gene sequence analysis revealed strain 2CBH44T had relatively low 16S rRNA gene sequence similarity (97.5%) to Coprobacter secundus 177T . However, strain 2CBH44T showed 96.9% average nucleotide identity value with C. secundus 177T , indicating that strain 2CBH44T and C. secundus 177T belong to the same species. On the other hand, the digital DNA-DNA hybridization value between strain 2CBH44T and C. secundus 177T was 73.5%, indicating that strain 2CBH44T is a subspecies of C. secundus. Other anaerobic, Gram-stain-variable, rod-shaped bacterium, designated strain 12CBH8T , was also isolated from human feces. Strain 12CBH8T had significantly low 16S rRNA gene sequence similarities (

Some sea slugs sequester chloroplasts from algal food in their intestinal cells and photosynthesize for months. This phenomenon, kleptoplasty, poses a question of how the chloroplast retains its activity without the algal nucleus. There have been debates on the horizontal transfer of algal genes to the animal nucleus. To settle the arguments, this study reported the genome of a kleptoplastic sea slug, Plakobranchus ocellatus, and found no evidence of photosynthetic genes encoded on the nucleus. Nevertheless, it was confirmed that light illumination prolongs the life of mollusk under starvation. These data presented a paradigm that a complex adaptive trait, as typified by photosynthesis, can be transferred between eukaryotic kingdoms by a unique organelle transmission without nuclear gene transfer. Our phylogenomic analysis showed that genes for proteolysis and immunity undergo gene expansion and are up-regulated in chloroplast-enriched tissue, suggesting that these molluskan genes are involved in the phenotype acquisition without horizontal gene transfer.

Clostridioides difficile causes nosocomial antibiotic-associated diarrhea on a global scale. Susceptibility to C. difficile infection (CDI) is influenced by the composition and metabolism of gut microbiota, which in turn are affected by diet. However, the mechanism underlying the interplay between diet and gut microbiota that modulates susceptibility to CDI remains unclear. Here, we show that a soy protein diet increases the mortality of antibiotic-treated, C. difficile-infected mice while also enhancing the intestinal levels of amino acids (aas) and relative abundance of Lactobacillus genus. Indeed, Ligilactobacillus murinus-mediated fermentation of soy protein results in the generation of aas, thereby promoting C. difficile growth, and the process involves the anchored cell wall proteinase PrtP. Thus, mutual interaction between dietary protein and the gut microbiota is a critical factor affecting host susceptibility to CDI, suggesting that dietary protein sources can be an important determinant in controlling the disease.

Asobara japonica is an endoparasitic wasp that parasitizes Drosophila flies. It synthesizes various toxic components in the venom gland and injects them into host larvae during oviposition. To identify and characterize these toxic components for enabling parasitism, we performed the whole-genome sequencing (WGS) and devised a protocol for RNA interference (RNAi) with A. japonica. Because it has a parthenogenetic lineage due to Wolbachia infection, we generated a clonal strain from a single wasp to obtain highly homogenous genomic DNA. The WGS analysis revealed that the estimated genome size was 322 Mb with a heterozygosity of 0.132%. We also performed RNA-seq analyses for gene annotation. Based on the qualified WGS platform, we cloned ebony-Aj, which encodes the enzyme N-β-alanyl dopamine synthetase, which is involved in melanin production. The microinjection of double-stranded RNA (dsRNA) targeting ebony-Aj led to body colour changes in adult wasps, phenocopying ebony-Dm mutants. Furthermore, we identified putative venom genes as a target of RNAi, confirming that dsRNA injection-based RNAi specifically suppressed the expression of the target gene in wasp adults. Taken together, our results provide a powerful genetic toolkit for studying the molecular mechanisms of parasitism.

Genome-wide association studies (GWASs) can reveal genetic variations associated with a phenotype in the absence of any hypothesis of candidate genes. The problem of false-positive sites linked with the responsible site might be bypassed in bacteria with a high homologous recombination rate, such as Helicobacter pylori, which causes gastric cancer (GC). We conducted a GWAS followed by regression-based prediction of GC and duodenal ulcer H. pylori strains. We identified 14 single nucleotide polymorphisms (11 amino acid changes) that, combined, allowed effective disease discrimination. They were often informative of the underlying molecular mechanisms, such as electric charge alteration at the ligand-binding pocket, alteration in subunit interaction, and mode-switching of DNA methylation. We also identified three novel virulence factors/oncoprotein candidates. These results provide both defined targets for further informatic and experimental analyses to gain insights into GC pathogenesis and a basis for identifying a set of biomarkers for application in clinical settings.

Dosage compensation (DC) on the X chromosome is a mechanism to counteract the deleterious effects by gene loss from the Y chromosome. However, DC cannot work efficiently if the X chromosome also degenerates. This indeed occurs in the neo-sex chromosomes inDrosophila miranda, where neo-X as well as neo-Y chromosomes are under accelerated pseudogenization. To examine the generality of this pattern, we investigated the evolution of two additional neo-sex chromosomes that independently emerged inD. albomicansandD. americanaand compared their evolutionary processes with that inD. miranda. Comparative genomic and transcriptomic analyses revealed that the pseudogenization rate on neo-X is also accelerated in the two species (though lesser extent inD. americana). We also found that neo-X-linked genes whose neo-Y homologs are pseudogenized tend to be upregulated more stringently than those whose neo-Y homologs remain functional. Moreover, the genes under strong functional constraints and highly expressed in the testis tended to remain functional on neo-X and neo-Y, respectively. Focusing on theD. mirandaandD. albomicansneo-sex chromosomes that independently emerged from the same autosome, we further found that the same genes tend to have been pseudogenized in parallel on neo-Y. Those genes includeIdgf6andJhI-26whose functions seem to be unnecessary or could be even harmful for males. These results indicate that neo-sex chromosomes inDrosophilashare a common evolutionary trajectory after their emergence, which may be applicable to other sex chromosomes in a variety of organisms to avoid being an evolutionary dead-end.

Dosage compensation (DC) on the X Chromosome counteracts the deleterious effects of gene loss on the Y Chromosome. However, DC is not efficient if the X Chromosome also degenerates. This indeed occurs in Drosophila miranda, in which both the neo-Y and the neo-X are under accelerated pseudogenization. To examine the generality of this pattern, we investigated the evolution of two additional neo-sex chromosomes that emerged independently in D. albomicans and D. americana and reanalyzed neo-sex chromosome evolution in D. miranda. Comparative genomic and transcriptomic analyses revealed that the pseudogenization rate on the neo-X is also accelerated in D. albomicans and D. americana although to a lesser extent than in D. miranda. In males, neo-X-linked genes whose neo-Y-linked homologs are pseudogenized tended to be up-regulated more than those whose neo-Y-linked homologs remain functional. Moreover, genes under strong functional constraint and genes highly expressed in the testis tended to remain functional on the neo-X and neo-Y, respectively. Focusing on the D. miranda and D. albomicans neo-sex chromosomes that emerged independently from the same autosome, we further found that the same genes tend to become pseudogenized in parallel on the neo-Y. These genes include Idgf6 and JhI-26, which may be unnecessary or even harmful in males. Our results indicate that neo-sex chromosomes in Drosophila share a common evolutionary trajectory after their emergence, which may prevent sex chromosomes from being an evolutionary dead end.

Genome sequencing of a unique frog (Glandirana rugosa) having XY/ZW systems within the species and comparative genomics with other six frogs were performed using a batch-learning self-organizing map, which is unsupervised AI for oligonucleotide compositions, to clarify its genome characteristics., The Japanese wrinkled frog (Glandirana rugosa) is unique in having both XX-XY and ZZ-ZW types of sex chromosomes within the species. The genome sequencing and comparative genomics with other frogs should be important to understand mechanisms of turnover of sex chromosomes within one species or during a short period. In this study, we analyzed the newly sequenced genome of G. rugosa using a batch-learning self-organizing map which is unsupervised artificial intelligence for oligonucleotide compositions. To clarify genome characteristics of G. rugosa, we compared its short oligonucleotide compositions in all 1-Mb genomic fragments with those of other six frog species (Pyxicephalus adspersus, Rhinella marina, Spea multiplicata, Leptobrachium leishanense, Xenopus laevis, and Xenopus tropicalis). In G. rugosa, we found an Mb-level large size of repeat sequences having a high identity with the W chromosome of the African bullfrog (P. adspersus). Our study concluded that G. rugosa has unique genome characteristics with a high CG frequency, and its genome is assumed to heterochromatinize a large size of genome via methylataion of CG.

Abyssinian black-and-white colobus (Colobus guereza) inhabits west, central, and east Africa and lives in social groups in the rain forests. Colobus guereza has a unique foregut digestive system similar to that of ruminants and absorbs organic acids, including Short-Chain Fatty Acids (SCFAs), as energy sources derived from microbial fermentation of plant materials in the gastrointestinal tract. In this study, the gastrointestinal metabolic and microbial features of a male Colobus guereza singly housed in a zoo were characterized, and each digesta from each segment of the digestive tract was collected and subjected to biochemical, microbiome, and metagenome analyses. In this case report, high levels of acetate and propionate were observed in the foregut, while a relatively high level of lactate was detected in the small intestine. Moreover, in the hindgut and the feces, acetate was dominant compared to the other SCFAs. SCFAs analysis indicated that Colobus guereza obtains energy via SCFAs, especially acetate, fermentation in the foregut and hindgut. A metagenome analysis revealed that each part of the digestive tract of Colobus guereza has a unique microbiota. Similar to ruminants, Prevotella and Selenomonas were the dominant genera in the foregut, which may indicate microbial fermentation of plant materials in the foregut of Colobus guereza.

Altered brain energy metabolism associated with increase in lactate levels and the resultant decrease in pH have been increasingly implicated in multiple neuropsychiatric disorders, such as schizophrenia, bipolar disorder, autism spectrum disorder and neurodegenerative disorders. Although it is controversial, change of pH/ lactate level as a primary feature of these diseases, rather than a result of confounding factors such as medication and agonal state, has been evidenced. Animal models that can be studied without such confounding factors inherent to humans are a suitable alternative to understand the controversy. However, the knowledge in animal models regarding brain pH and lactate and their relation to behavioral outcomes is limited in the context of neuropsychiatric disease conditions. In this study, we investigated the common occurrence of changes in the pH and lactate levels in the brain in animal models by analyzing 65 animal models related to neuropsychiatric and neurodegenerative diseases with 1,239 animals. Additionally, we evaluated the behavioral phenotypes relative to the chemical changes in the brain. Among the models, 27 and 24 had significant changes in brain pH and lactate levels, respectively, including Shank2 KO mice, Clock mutant mice, serotonin transporter KO mice, mice with a paternal duplication of human chromosome 15q11-13, Fmr1 KO mice, BTBR mice, APP-J20 Tg mice, social defeat stress-exposed mice, corticosterone-treated mice, and streptozotocin-induced diabetic mice. Meta-analysis of the data revealed a highly significant negative correlation between brain pH and lactate levels, suggestive of increased lactate levels causing decreased brain pH. Statistical learning algorithm based on the comprehensive data has revealed that the increased brain lactate levels can be predominantly predicted by the indices for the percentage of correct response in working memory test, with a significant simple, negative correlation. Our results suggest that brain energy metabolism is commonly altered in many animal models of neuropsychiatric and neurodegenerative diseases, which may be associated with working memory performance. We consider our study to be an essential step suggesting that the brain endophenotypes serve as a basis for the transdiagnostic characterization of the biologically heterogeneous and debilitating cognitive illnesses. Based on these results, we are openly accepting collaborations to extend these findings and to test the hypotheses generated in this study using more animal models. We welcome any mice/rat models of diseases with or without any behavioral phenotypes.

In animal gonads, transposable elements (TEs) are actively repressed to preserve genome integrity through the Piwi-interacting RNA (piRNA) pathway. In mice, piRNAs are most abundantly expressed in male germ cells, and form effector complexes with three distinct PIWI proteins. The depletion of individual Piwi genes causes male-specific sterility owing to severe defects in spermatogenesis with no discernible phenotype in female mice. Unlike mice, most other mammals have four PIWI genes, some of which are expressed in the ovary. Here, purification of PIWI complexes from oocytes of the golden hamster revealed that the size of the piRNAs loaded onto PIWIL1 changed during oocyte maturation. In contrast, PIWIL3, an ovary-specific PIWI in most mammals, associates with short piRNAs only in metaphase II oocytes, which coincides with intense phosphorylation of the protein. An improved high-quality genome assembly and annotation revealed that PIWIL1- and PIWIL3-associated piRNAs appear to share the 5′- ends of common piRNA precursors and are mostly derived from unannotated sequences with a diminished contribution from TE-derived sequences, most of which correspond to endogenous retroviruses (ERVs). Although binding sites for the transcription factor A-Myb are identified in the transcription start site regions of the testis piRNA clusters, the piRNA clusters in the ovary show no well-defined binding motifs in their upstream regions. These results show that hamster piRNA clusters are transcribed by different transcriptional factors in the ovary and testis, resulting in the generation of sex-specific piRNAs. Our findings show the complex and dynamic nature of biogenesis of piRNAs in hamster oocytes, and together with the new genome sequence generated, serve as the foundation for developing useful models to study the piRNA pathway in mammalian oocytes.Highlights-The size of PIWIL1-associated piRNAs changes during oocyte maturation-Phosphorylation of PIWIL3 in MII oocytes coincides with its association with small 19-nt piRNAs-Improved high-quality genome assembly and annotation identifies young endogenous retroviruses as major targets of piRNAs in hamster oocytes-PIWIL1- and PIWIL3-associated piRNAs share the 5′-ends of the common piRNA precursors in oocytes

The ancestral tomato species are known to possess genes that are valuable for improving traits in breeding. Here, we aimed to construct high-quality de novo genome assemblies of Solanum pimpinellifolium ‘LA1670’ and S. lycopersicum var. cerasiforme ‘LA1673’, originating from Peru. The Pacific Biosciences (PacBio) long-read sequences with 110× and 104× coverages were assembled and polished to generate 244 and 202 contigs spanning 808.8 Mbp for ‘LA1670’ and 804.5 Mbp for ‘LA1673’, respectively. After chromosome-level scaffolding with reference guiding, 14 scaffold sequences corresponding to 12 tomato chromosomes and 2 unassigned sequences were constructed. High-quality genome assemblies were confirmed using the Benchmarking Universal Single-Copy Orthologs and long terminal repeat assembly index. The protein-coding sequences were then predicted, and their transcriptomes were confirmed. The de novo assembled genomes of S. pimpinellifolium and S. lycopersicum var. cerasiforme were predicted to have 71,945 and 75,230 protein-coding genes, including 29,629 and 29,185 non-redundant genes, respectively, as supported by the transcriptome analysis results. The chromosome-level genome assemblies coupled with transcriptome data sets of the two accessions would be valuable for gaining insights into tomato domestication and understanding genome-scale breeding.

Epstein–Barr virus (EBV) establishes lifelong latent infection in the majority of healthy individuals, while it is a causative agent for various diseases, including some malignancies. Recent high-throughput sequencing results indicate that there are substantial levels of viral genome heterogeneity among different EBV strains. However, the extent of EBV strain variation among asymptomatically infected individuals remains elusive. Here, we present a streamlined experimental strategy to clone and sequence EBV genomes derived from human tonsillar tissues, which are the reservoirs of asymptomatic EBV infection. Complete EBV genome sequences, including those of repetitive regions, were determined for seven tonsil-derived EBV strains. Phylogenetic analyses based on the whole viral genome sequences of worldwide non-tumour-derived EBV strains revealed that Asian EBV strains could be divided into several distinct subgroups. EBV strains derived from nasopharyngeal carcinoma-endemic areas constitute different subgroups from a subgroup of EBV strains from non-endemic areas, including Japan. The results could be consistent with biased regional distribution of EBV-associated diseases depending on the different EBV strains colonizing different regions in Asian countries.

Tameness is a major element of animal domestication and involves two components: motivation to approach humans (active tameness) and reluctance to avoid humans (passive tameness). To understand the behavioral and genetic mechanisms of active tameness in mice, we had previously conducted selective breeding for long durations of contact and heading toward human hands in an active tameness test using a wild‐derived heterogeneous stock. Although the study showed a significant increase in contacting and heading with the 12th generation of breeding, the effect on other behavioral indices related to tameness and change of gene expression levels underlying selective breeding was unclear. Here, we analyzed nine tameness‐related traits at a later stage of selective breeding and analyzed how gene expression levels were changed by the selective breeding. We found that five traits, including contacting and heading, showed behavioral change in the selective groups comparing to the control through the generations. Furthermore, we conducted cluster analyses to evaluate the relationships among the nine traits and found that contacting and heading combined in an independent cluster in the selected groups, but not in the control groups. RNA‐Seq of hippocampal tissue revealed differential expression of 136 genes between the selection and control groups, while the pathway analysis identified the networks associated with these genes. These results suggest that active tameness was hidden in the control groups but became apparent in the selected populations by selective breeding, potentially driven by changes in gene expression networks., Tameness can be divided into two components: active tameness and passive tameness. We conducted selective breeding for higher scores of active tameness using a wild‐derived heterogeneous stock. In the analyses of nine behavioral traits related to tameness, we found that five traits showed changes in the selective groups compared to the control groups through the generations. We conducted cluster analyses to evaluate the relationship among the nine traits and RNA‐Seq analysis to characterize the molecular network related to tameness. The results suggested that active tameness was hidden in the control groups but became apparent in the selected populations by selective breeding, potentially driven by changes in gene expression networks.

A new disease in rice that is characterized by leaf bleaching was recently identified in some fields in the Mekong Delta region of Vietnam. The present study was the first to isolate and identify the pathogen of this disease. We confirmed that leaf bleaching symptoms were due to infection with Methylobacterium indicum bacteria using molecular biology approaches. A full-length genome analysis of pathogenic Methylobacterium strain VL1 revealed that it comprises a single chromosome and six plasmids, with a total size of 7.05‍ ‍Mbp and GC content of 70.5%. The genomic features of VL1 were similar to those of the non-pathogenic M. indicum strain SE2.11T; however, VL1 possessed additional unique genes, including those related to homoserine lactone biosynthesis. We established a loop-mediated isothermal amplification (LAMP) assay using the unique sequences of VL1 as target sequences for the rapid and simple detection of pathogenic M. indicum strains. Our initial evaluation demonstrated that the LAMP assay successfully distinguished between pathogenic and non-pathogenic strains infecting rice plants in a rapid and sensitive manner. The present results provide insights into the pathogenesis and development of control measures for novel rice diseases.

We present here the second complete genome of anaerobic ammonium oxidation (anammox) bacterium, Candidatus (Ca.) Brocadia pituitae, along with those of a nitrite oxidizer and two incomplete denitrifiers from the anammox bacterial community (ABC) metagenome. Although NO2− reduction to NO is considered to be the first step in anammox, Ca. B. pituitae lacks nitrite reductase genes (nirK and nirS) responsible for this reaction. Comparative genomics of Ca. B. pituitae with Ca. Kuenenia stuttgartiensis and six other anammox bacteria with nearly complete genomes revealed that their core genome structure contains 1,152 syntenic orthologues. But nitrite reductase genes were absent from the core, whereas two other Brocadia species possess nirK and these genes were horizontally acquired from multiple lineages. In contrast, at least five paralogous hydroxylamine oxidoreductase genes containing candidate ones (hao2 and hao3) encoding another nitrite reductase were observed in the core. Indeed, these two genes were also significantly expressed in Ca. B. pituitae as in other anammox bacteria. Because many nirS and nirK genes have been detected in the ABC metagenome, Ca. B. pituitae presumably utilises not only NO supplied by the ABC members but also NO and/or NH2OH by self-production for anammox metabolism.

Here, we report the complete genome sequence of Faecalibacillus intestinalis JCM 34082, a member of the family Erysipelotrichaceae that was isolated from feces from a healthy Japanese woman. The genome assembly comprised 2,869,982 bp, with a G+C content of 29.8%.

Here, we report the complete genome sequence of Faecalibacillus intestinalis JCM 34082, a member of the family Erysipelotrichaceae that was isolated from feces from a healthy Japanese woman. The genome assembly comprised 2,869,982 bp, with a G+C content of 29.8%.

Prevention of stress‐induced adverse effects is important for animals and humans to maintain their quality of life (QOL). Stress decreases the productivity of farm animals and induces abnormal behaviors, which is one of the major problems in animal welfare. In humans, stress increases the risk of mental illness which adversely impacts QOL. Stress is, thus, a common health problem for both animals and humans, and stress prevention and promotion of stress resilience could improve animal and human health and QOL. Among various stresses, psychosocial stress experienced by individuals is particularly difficult to prevent and it could, thus, prove beneficial to attempt to increase resilience to psychosocial stress. There exist a few critical interventions for promoting such resilience, environmental enrichment being one. However, this review describes recent progress in nutritional interventions that could confer resilience to psychosocial stress. The efficacy of this intervention is studied in the social defeat model mouse, which is a standard model for studying psychosocial stress. Several nutrients were found to rescue stress vulnerability using the models. Furthermore, probiotics and prebiotics became crucial dietary interventions for combating psychosocial stress. Collectively, dietary intake of appropriate nutrients will be more important for maintaining QOL in animals and humans.

Forward genetic analysis in a parasitic plant revealed that ethylene signaling is essential for host invasion., Parasitic plants form a specialized organ, a haustorium, to invade host tissues and acquire water and nutrients. To understand the molecular mechanism of haustorium development, we performed a forward genetics screening to isolate mutants exhibiting haustorial defects in the model parasitic plant Phtheirospermum japonicum. We isolated two mutants that show prolonged and sometimes aberrant meristematic activity in the haustorium apex, resulting in severe defects on host invasion. Whole-genome sequencing revealed that the two mutants respectively have point mutations in homologs of ETHYLENE RESPONSE 1 (ETR1) and ETHYLENE INSENSITIVE 2 (EIN2), signaling components in response to the gaseous phytohormone ethylene. Application of the ethylene signaling inhibitors also caused similar haustorial defects, indicating that ethylene signaling regulates cell proliferation and differentiation of parasite cells. Genetic disruption of host ethylene production also perturbs parasite invasion. We propose that parasitic plants use ethylene as a signal to invade host roots.

Genome characterization of California poppy (Eschscholzia californica cv. “Hitoezaki”), which produces pharmaceutically important benzylisoquinoline alkaloids (BIAs), was carried out using the draft genome sequence. The numbers of tRNA and rRNA genes were close to those of the other plant species tested, whereas the frequency of repetitive sequences was distinct from those species. Comparison of the predicted genes with those of Amborella trichopoda, Nelumbo nucifera, Solanum lycopersicum, and Arabidopsis thaliana, and analyses of gene ontology and Kyoto Encyclopedia of Genes and Genomes pathway indicated that the enzyme genes involved in BIA biosynthesis were highly enriched in the California poppy genome. Further comparative analysis using the genome information of Papaver somniferum and Aquilegia coerulea, both BIA-producing plants, revealed that many genes encoding BIA biosynthetic enzymes, transcription factors, transporters, and candidate proteins, possibly related to BIA biosynthesis, were specifically distributed in these plant species.

The cryptic parasite Sparganum proliferum proliferates in humans and invades tissues and organs. Only scattered cases have been reported, but S. proliferum infection is always fatal. However, S. proliferum’s phylogeny and life cycle remain enigmatic. To investigate the phylogenetic relationships between S. proliferum and other cestode species, and to examine the mechanisms underlying pathogenicity, we sequenced the entire genomes of S. proliferum and a closely related non–life-threatening tapeworm Spirometra erinaceieuropaei. Additionally, we performed larvae transcriptome analyses of S. proliferum plerocercoid to identify genes involved in asexual reproduction in the host. The genome sequences confirmed that the S. proliferum has experienced a clearly distinct evolutionary history from S. erinaceieuropaei. Moreover, we found that nonordinal extracellular matrix coordination allows asexual reproduction in the host, and loss of sexual maturity in S. proliferum are responsible for its fatal pathogenicity to humans. Our high-quality reference genome sequences should be valuable for future studies of pseudophyllidean tapeworm biology and parasitism.

An increasing volume of empirical studies demonstrated that hybridization between distant lineages may have promoted speciation in various taxa. However, the timing, extent and direction of introgressive hybridization remain unknown in many cases. Here, we report a possible case in which repeated hybridization promoted divergence of Oryzias ricefishes (Adrianichthyidae) on Sulawesi, an island of Wallacea. Four Oryzias species are endemic to the Malili Lake system in central Sulawesi, which is composed of five tectonic lakes; of these, one lake is inhabited by two species. Morphological and population genomic analyses of genome-wide single-nucleotide polymorphisms revealed that these two sympatric species are phylogenetically sister to but substantially reproductively isolated from each other. Analyses of admixture and comparison of demographic models revealed that the two sympatric species experienced several substantial introgressions from outgroup populations that probably occurred soon after they had secondary contact with each other in the lake. However, the ratio of migrants from the outgroups was estimated to be different between the two species, which is consistent with the hypothesis that these introgressions aided their divergence or prevented them from forming a hybrid swarm. Repeated lake fragmentations and fusions may have promoted diversification of this freshwater fish species complex that is endemic to this ancient lake system.

The cryptic parasite Sparganum proliferum proliferates in humans and invades tissues and organs. Only scattered cases have been reported, but S. proliferum infection is always fatal. However, S. proliferum’s phylogeny and life cycle remain enigmatic. To investigate the phylogenetic relationships between S. proliferum and other cestode species, and to examine the mechanisms underlying pathogenicity, we sequenced the entire genomes of S. proliferum and a closely related non–life-threatening tapeworm Spirometra erinaceieuropaei. Additionally, we performed larvae transcriptome analyses of S. proliferum plerocercoid to identify genes involved in asexual reproduction in the host. The genome sequences confirmed that the S. proliferum has experienced a clearly distinct evolutionary history from S. erinaceieuropaei. Moreover, we found that nonordinal extracellular matrix coordination allows asexual reproduction in the host, and loss of sexual maturity in S. proliferum are responsible for its fatal pathogenicity to humans. Our high-quality reference genome sequences should be valuable for future studies of pseudophyllidean tapeworm biology and parasitism., Kikuchi et al. sequence the genome of the fatal tapeworm Sparganum proliferum and a closely related non–life-threatening tapeworm Spirometra erinaceieuropaei, and describe its genomic features suggesting the natural history and molecular mechanisms underlying pathogenicity. Their findings indicate that nonordinal extracellular matrix coordination is important for its asexual reproduction, and suggest that loss of sexual maturity contributes to the fatal pathogenicity of S. proliferum to humans.

Funder: The excavation of the Ikawazu Jomon individual was supported by Grant-in-Aid for Scientific Research (B) (25284157) to YY. The Ikawazu Jomon genome project was organized by HI, and TH & HO who were supported by MEXT KAKENHI Grant Numbers 16H06408 and 17H05132, by Grant-in-Aid for Scientific Research on Innovative Areas (Cultural History of Paleoasia), and by Grant-in-Aid for Challenging Exploratory Research (23657167) and for Scientific Research (B) (17H03738). The Ikawazu Jomon genome sequencing was supported by JSPS KAKENHI Grant Number 16H06279 to ATo, and partly supported in the CHOZEN project in Kanazawa University, and in the Cooperative Research Project Program of the Medical Institute of Bioregulation, Kyushu University. Computations for the Ikawazu Jomon genome were partially performed on the NIG supercomputer at ROIS National Institute of Genetics., Anatomically modern humans reached East Asia more than 40,000 years ago. However, key questions still remain unanswered with regard to the route(s) and the number of wave(s) in the dispersal into East Eurasia. Ancient genomes at the edge of the region may elucidate a more detailed picture of the peopling of East Eurasia. Here, we analyze the whole-genome sequence of a 2,500-year-old individual (IK002) from the main-island of Japan that is characterized with a typical Jomon culture. The phylogenetic analyses support multiple waves of migration, with IK002 forming a basal lineage to the East and Northeast Asian genomes examined, likely representing some of the earliest-wave migrants who went north from Southeast Asia to East Asia. Furthermore, IK002 shows strong genetic affinity with the indigenous Taiwan aborigines, which may support a coastal route of the Jomon-ancestry migration. This study highlights the power of ancient genomics to provide new insights into the complex history of human migration into East Eurasia.

Sexual selection drives rapid phenotypic diversification of mating traits. However, we know little about the causative genes underlying divergence in sexually selected traits. Here, we investigate the genetic basis of male mating trait diversification in the medaka fishes (genus Oryzias) from Sulawesi, Indonesia. Using linkage mapping, transcriptome analysis, and genome editing, we identify csf1 as a causative gene for red pectoral fins that are unique to male Oryzias woworae. A cis-regulatory mutation enables androgen-induced expression of csf1 in male fins. csf1-knockout males have reduced red coloration and require longer for mating, suggesting that coloration can contribute to male reproductive success. Contrary to expectations, non-red males are more attractive to a predatory fish than are red males. Our results demonstrate that integrating genomics with genome editing enables us to identify causative genes underlying sexually selected traits and provides a new avenue for testing theories of sexual selection., The mutations underlying sexually selected traits like the red fins on a male medaka fish can be hard to pinpoint. Using a new genome, transcriptomics and gene editing, Ansai et al. find that the gene csf1 causes male fins to be red, which attracts females and, surprisingly, is less attractive to predators.