Search

In the Japanese hagfish,Eptatretus burgeri, approximately 21% of the genomic DNA in germ cells (2n=52) consists of 16 chromosomes (eliminated [E]-chromosomes) that are eliminated from presumptive somatic cells (2n=36). To uncover the eliminated genome (E-genome), we have identified 16 eliminated repetitive DNA families from eight hagfish species, with 11 of these repeats being selectively amplified in the germline genome ofE. burgeri. Furthermore, we have demonstrated that six of these sequences, namely EEEb1–6, are exclusively localized on all 16 E-chromosomes. This has led to the hypothesis that the eight pairs of E-chromosomes are derived from one pair of ancestral chromosomes via multiple duplication events over a prolonged evolutionary period. NGS analysis has recently facilitated the re-assembly of two distinct draft genomes ofE. burgeri, derived from the testis and liver. This advancement allows for the prediction of not only nonrepetitive eliminated sequences but also over 100 repetitive and eliminated sequences, accomplished through K-mer-based analysis. In this study, we report four novel eliminated repetitive DNA sequences (designated as EEEb7–10) and confirm the relative chromosomal localization of all eliminated repeats (EEEb1–10) by fluorescencein situhybridization (FISH). With the exception of EEEb10, all sequences were exclusively detected on EEEb1-positive chromosomes. Surprisingly, EEEb10 was detected as an intense signal on EEEb1-positive chromosomes and as a scattered signal on other chromosomes in germ cells. The study further divided the eight pairs of E-chromosomes into six groups based on the signal distribution of each DNA family, and fiber-FISH experiments showed that the EEEb2–10 family was dispersed in the EEEb1-positive extended chromatin fiber. These findings provide new insights into the mechanisms underlying chromosome elimination and the evolution of E-chromosomes, supporting our previous hypothesis.

Genome sizes are known to vary within and among closely related species, but the knowledge about genomic factors contributing to the variation and their impacts on gene functions is limited to only a small number of species. This study identified a more than twofold heritable genome size variation among the unicellular Zygnematophycean alga,Closterium peracerosum-strigosum-littorale(C. psl.) complex, based on short-read sequencing analysis of 22 natural strains and F1segregation analysis. Sixde novoassembled genomes revealed that genome size variation is largely attributable to genome-wide copy number variation (CNV) among strains rather than mating type-linked genomic regions or specific repeat sequences such as rDNA. Notably, about 30% of genes showed CNV even between strains that can mate with each other. Transcriptome and gene ontology analysis demonstrated that CNV is distributed nonrandomly in terms of gene functions, such that CNV was more often observed in the gene set with stage-specific expression. Furthermore, in about 30% of these genes with CNV, the expression level does not increase proportionally with the gene copy number, suggesting presence of dosage compensation, which was overrepresented in genes involved in basic biological functions, such as translation. Nonrandom patterns in gene duplications and corresponding expression changes in terms of gene functions may contribute to maintaining the high level of CNV associated with extensive genome size variation in theC. psl. complex, despite its possible detrimental effects.

SummaryThe formation of RAD51/DMC1 filaments on single-stranded (ss)DNAs, which is essential for homology search and strand exchange in DNA double-strand break (DSB) repair and for the protection of stalled DNA replication forks, is tightly regulated in time and space. FIGNL1 AAA+++ ATPase plays positive and negative roles in the RAD51-mediated recombination in human cells. However, the role of FIGNL1 in gametogenesis remains unsolved. Here, we characterized a germ-line-specific conditional knockout (cKO) mouse ofFIGNL1. TheFignl1cKO male mice showed defective chromosome synapsis and impaired meiotic DSB repair with the accumulation of RAD51/DMC1 on chromosomes in mid-meiotic prophase I, supporting a role of FIGNL1 in a post-assembly stage of RAD51/DMC1 filaments.Fignl1cKO spermatocytes accumulate RAD51 and DMC1 ensembles on chromosomes not only in early meiotic prophase I but also in meiotic S-phase. These RAD51/DMC1 assemblies are independent of meiotic DSB formation. Finally, we showed that purified FIGNL1 dismantles RAD51 filament on double-stranded (ds)DNA as well as ssDNA. These results suggest a critical role of FIGNL1 to limit the uncontrolled assembly of RAD51 and DMC1 on native dsDNAs during the meiotic S-phase and meiotic prophase I.

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.

Transposable elements (TEs) are among the most dynamic parts of genomes. Since TEs are potentially deleterious, eukaryotes silence them through epigenetic mechanisms such as DNA methylation and RNAi. We have previously reported that Arabidopsis TEs, calledVANDALs, counteract epigenetic silencing through a group of sequence-specific anti-silencing proteins, VANCs. VANC proteins bind to noncoding regions of specificVANDALcopies and induce a loss of silent chromatin marks. Sequence-specific anti-silencing allows these TEs to proliferate with minimum host damage. Here, we show that RNAi efficiently targets noncoding regions ofVANDALTEs to silence themde novo. Target motifs of VANC, in turn, evolved to escape RNAi. Escaping RNAi could be the primary event leading to the differentiation of sequence-specific anti-silencing systems. We propose that this selfish behaviour of TEs paradoxically could make them less harmful to the host.

The Iberian ribbed newt (Pleurodeles waltl) constitutes a central model for probing the basis of regeneration. Here, we present the sequencing and chromosome-scale assembly of the 20.3GbP.waltlgenome, which exhibits the highest level of contiguity and completeness among giant genomes. We uncover that DNA transposable elements are the major contributors to its expansion, with hAT transposons comprising a large portion of repeats. Several hATs are actively transcribed and differentially expressed during adultP.waltllimb regeneration, along with domesticated hAT transposons of the ZBED transcription factor family. Despite its size, syntenic relationships are conserved across the genome. As an example we show the high degree of conservation of the regeneration-associated Tig1 locus with several neighbouring genes. Together, theP.waltlgenome provides a fundamental resource for the study of regenerative, developmental and evolutionary principles.

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, for insertion of long transgenes such as fluorescent reporters in human cells, the optimal type of HDR donors remains unclear. In this study, we established a simple and efficient CRISPR-mediated knock-in method for long transgenes using linear dsDNA and ssDNA donors, and systematically compared the performance of these two donors for endogenous gene tagging in human non-transformed diploid cells. Quantification using flow cytometry revealed higher efficiency of fluorescent tagging with dsDNA donors than with ssDNA. By analyzing knock-in outcomes using long-read amplicon sequencing and a classification framework, a variety of mis-integration events were detected regardless of the donor type. Importantly, the ratio of precise insertion was higher with dsDNA donors than with ssDNA. Moreover, in off-target integration analyses, dsDNA and ssDNA were comparably prone to non-homologous integration. These results indicate that ssDNA is not superior to dsDNA as long HDR donors for gene knock-in in human cells.

De novo mutations accumulate with zygotic cell divisions. However, the occurrence of these mutations and the way they are inherited by somatic cells and germ cells remain unclear. Here, we present a novel method to reconstruct cell lineages. We identified mosaic mutations in mice using deep whole-genome sequencing and reconstructed embryonic cell lineages based on the variant allele frequencies of the mutations. The reconstructed trees were confirmed using nuclear transfer experiments and the genotyping of approximately 50 offspring of each tree. The most detailed tree had 32 terminal nodes and showed cell divisions from the fertilized egg to germ cell– and somatic cell–specific lineages, indicating at least five independent cell lineages that would be selected as founders of the primordial germ cells. The contributions of each lineage to germ cells and offspring varied widely. At the emergence of the germ cell–specific lineages, 10–15 embryonic mutations had accumulated, suggesting that the pregastrulation mutation rate is 1.0 mutation per mitosis. Subsequent mutation rates were 0.7 for germ cells and 13.2 for tail fibroblasts. Our results show a new framework to assess embryonic lineages; further, we suggest an evolutionary strategy for preserving heterogeneity owing to postzygotic mutations in offspring.

Gain of alternative splicing gives rise to functional diversity in proteins and underlies the complexity and diversity of biological aspects. However, it is still not fully understood how alternatively spliced genes develop the functional novelty. To this end, we infer the evolutionary history of the doublesex gene, the key transcriptional factor in the sexual differentiation of arthropods. doublesex is controlled by sex-specific splicing and promotes both male and female differentiation in some holometabolan insects. In contrast, doublesex promotes only male differentiation in some hemimetabolan insects. Here, we investigate ancestral states of doublesex using Thermobia domestica belonging to Zygentoma, the sister group of winged insects. We find that doublesex of T. domestica expresses sex-specific isoforms but is only necessary for male differentiation of sexual morphology. This result ensures the hypothesis that doublesex was initially only used to promote male differentiation during insect evolution. However, T. domestica doublesex has a short female-specific region and upregulates the expression of vitellogenin homologs in females, suggesting that doublesex may have already controlled some aspects of feminization in the common ancestor of winged insects. Reconstruction of the ancestral sequence and prediction of the protein structure show that the female-specific isoform of doublesex has a long C-terminal disordered region in holometabolan insects, but not in non-holometabolan species. We propose that doublesex acquired a female-specific isoform and then underwent a change in the protein motif structure, which became essential for female differentiation in sexual dimorphisms.

The evolutionarily conserved RNA helicase DDX6 is a central player of post-transcriptional regulation, but its role during embryogenesis remains elusive. We here demonstrated that DDX6 enables proper cell lineage specification from pluripotent cells by analyzing Ddx6 KO mouse embryos and in vitro epiblast-like cell (EpiLC) induction system. Our study unveiled a great impact of DDX6-mediated RNA regulation on signaling pathways. Deletion of Ddx6 caused the aberrant transcriptional upregulation of the negative regulators of BMP signaling, which accompanied with enhanced Nodal signaling. Ddx6△/△ pluripotent cells acquired higher pluripotency with a strong inclination toward neural lineage commitment. During gastrulation, abnormally expanded Nodal expression in the primitive streak likely promoted endoderm cell fate specification while inhibiting mesoderm development. We further clarified the mechanism how DDX6 regulates cell fate determination of pluripotent cells by genetically dissecting major DDX6 pathways: processing body (P-body) formation, translational repression, mRNA decay, and miRNA-mediated silencing. P-body-related functions were dispensable, but the miRNA pathway was essential for the DDX6 function. DDX6 may prevent aberrant transcriptional upregulation of the negative regulators of BMP signaling by repressing translation of certain transcription factors through the interaction with miRNA-induced silencing complexes (miRISCs). Overall, this delineates how DDX6 affects development of the three primary germ layers during early mouse embryogenesis and the underlying mechanism of DDX6 function.Author summaryGene expression occurs through the two steps: transcription (DNA to RNA) and translation (RNA to protein). Cells have very sophisticated regulatory processes working on various levels for the accurate gene expression. Post-transcriptional regulation, which includes all RNA-related controls, is crucial because it enables fine-tuning and rapid alteration of gene expression. RNA- binding proteins and non-coding RNAs are the two main players of post-transcriptional regulation. DDX6, the subject of our study, is an RNA-binding protein, more specifically an RNA helicase, which can unwind or rearrange RNA secondary structures. Its diverse molecular and cellular functions have been reported, but its embryogenic role is unknown. Here, we describe DDX6 function during early mouse embryogenesis and the underlying mechanism using genetic methodology. DDX6 enables proper cell lineage specification of pluripotent stem cells by mainly regulating BMP signaling through miRNA-mediated gene silencing. As DDX6- mediated RNA regulation affected signaling pathways, the loss of Ddx6 had a wide impact on developmental processes from pluripotency to embryo patterning. In addition, we identified which DDX6 molecular function is essential during early embryogenesis by genetically dissecting its main pathways.

Recent studies have revealed a surprising diversity of sex chromosomes in vertebrates. However, the detailed mechanism of their turnover is still elusive. To understand this process, it is necessary to compare closely related species in terms of sex-determining genes and the chromosomes harboring them. Here, we explored the genusTakifugu, in which one strong candidate sex-determining gene,Amhr2, has been identified. To trace the processes involved in transitions in the sex determination system in this genus, we studied 12 species and found that while theAmhr2locus likely determines sex in the majority ofTakifuguspecies, three species have acquired sex-determining loci at different chromosomal locations. Nevertheless, the generation of genome assemblies for the three species revealed that they share a portion of the male-specific supergene that contains a candidate sex-determining gene,GsdfY, along with genes that potentially play a role in male fitness. The shared supergene span approximately 100 kb and are flanked by two duplicated regions characterized by CACTA transposable elements. These results suggest that the shared supergene has taken over the role of sex-determining locus fromAmhr2in lineages leading to the three species, and repeated translocations of the supergene underlie the turnover of sex chromosomes in these lineages. These findings highlight the underestimated role of a mobile supergene in the turnover of sex chromosomes in vertebrates.SignificanceAlthough turnover of sex chromosomes is very common in many vertebrate lineages, the transition process is still elusive. We studied the sex-determining region (SDR) of 12 congeneric fish species. We found that while nine species retained their ancestral SDR, three species had acquired derived SDRs. Although the derived SDRs resided in three different chromosomes, they harbored a shared supergene flanked by two putative transposable elements. The results highlight the underestimated role of a mobile supergene in turnover of sex chromosomes in vertebrates.

The paper nautilus,Argonauta argo, also known as the greater argonaut, is a species of octopods distinctly characterized by its pelagic lifestyle and by the presence of a spiral-shaped shell-like eggcase in females. The eggcase functions by protecting the eggs laid inside it, and by building and keeping air intakes for buoyancy. To reveal the genomic background of the species’ adaptation to pelagic lifestyle and the acquisition of its shell-like eggcase, we sequenced the draft genome sequence of the species. The genome size was 1.1 Gb, which is the smallest among the cephalopods known to date, with the top 215 scaffolds (average length 5,064,479 bp) covering 81% (1.09 Gb) of the total assembly. A total of 26,433 protein-coding genes were predicted from 16,802 assembled scaffolds. From these, we identified nearly intact HOX, Parahox, Wnt clusters and some gene clusters probably related to the pelagic lifestyle, such asreflectin, tyrosinase, andopsin. For example,opsinmight have undergone an extensive duplication in order to adapt to the pelagic lifestyle, as opposed to other octopuses, which are mostly the benthic. Our gene models also discovered several genes homologous to those related to calcified shell formation in Conchiferan Mollusks, such as Pif-like, SOD, and TRX. Interestingly, comparative genomics analysis revealed that the homologous genes for such genes were also found in the genome of the octopus, which does not have a shell, as well as the basal cephalopodsNautilus. Therefore, the draft genome sequence ofA. argowe presented here had not only helped us to gain further insights into the genetic background of the dynamic recruitment and dismissal of genes for the formation of an important, converging extended phenotypic structure such as the shell and the shell-like eggcase, but also the evolution of lifestyles in Cephalopods and the octopods, from benthic to pelagic.

BackgroundThe intestinal microbiome is closely related to host health, and metatranscriptomic analysis can assess the functional activity of microbiomes by quantifying the bacterial gene expression level, which helps to elucidate the interaction between the microbiome and the environment. However, functional changes in the microbiome along the host intestinal tract remain unknown, and previous analytical methods have limitations, such as potentially overlooking unknown genes due to dependence on existing databases and being unable to take full advantage of metatranscriptome to reveal the functional change among multiple environments.ResultTo close these gaps, we develop a novel method that integrates metagenome and metatranscriptome to analyze the functional activity of microbiomes between intestinal sites. This method reconstructs a reference metagenomic sequence across multiple intestinal sites, allowing the gene expression levels of microbiome including unknown bacterial genes to be compared between multiple sites. As a result of applying this method to metatranscriptomic analysis in the intestinal tract of common marmoset, the reconstructed metagenome covered most of the expressed genes and it revealed that the changes in bacterial gene expressions among the caecum, transverse colon, and faeces were more dynamic and sensitive to environmental shifts than its abundance. In typical, the coenzyme synthesis gene and antibacterial resistance gene were more highly expressed in the caecum and transverse colon than in faeces, while there was no significant change in abundance of these genes.ConclusionOur findings demonstrate that an analytical method that integrates metagenome and metatranscriptome in multiple intestinal sites captures functional changes in the microbiomes at the gene resolution level.

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.

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.

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.

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.

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.

Restriction–modification (RM) systems are typically regarded as “primitive immune systems” in bacteria. The roles of methylation in gene regulation, segregation, and mismatch repair are increasingly recognized. To analyze methyltransferase (MTase) diversity in Streptococcus pyogenes, we compared the RM system distribution in eight new complete genome sequences obtained here and in the database-deposited complete genome sequences of 51 strains. The MTase gene distribution showed that type I MTases often change DNA sequence specificity via switching target recognition domains between strains. The type II MTases in the included strains fell into two groups: a prophage-dominant one and a CRISPR-dominant one. Some highly variable type II MTases were found in the prophage region, suggesting that MTases acquired from phage DNA can generate methylome diversity. Additionally, to investigate the possible contribution of DNA methylation to phenotype, we compared the methylomes and transcriptomes from the four most closely related strains, the results of which suggest that phage-derived methylases possibly regulate the methylome, and, hence, regulate expression levels in S. pyogenes. Our findings will benefit further experimental work on the relationship between virulence genes and pathogenicity in S. pyogenes.

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 slugPlakobranchus ocellatusand 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 DNA-independent phenotype acquisition.

Like animals, plants accommodate a rich diversity of microbes, typically without discernible disease symptoms. How their pathogenesis is prevented in the host remains obscure. Here, we show that the root-infecting fungus Colletotrichum fructicola of the C. gloeosporioides clade (CgE), isolated from field-grown healthy Brassicaceae plants, inhibits growth of pathogenic fungi in Arabidopsis thaliana, in a phosphate status-dependent manner. Loss of host ethylene signaling or phytoalexins, camalexin or indole glucosinolates, however, allows CgE to display pathogenesis, suggesting host contributions to endophytic CgE colonization and benefit. Compared to a closely-related C. gloeosporioides pathogen (CgP), CgE is characterized by genome expansion and >700 fungal genes (4.34%) specifically induced in the host roots when co-inoculated with CgP, including genes related to fungal secondary metabolism. This may underlie antimicrobial tolerance of CgE and its dominance over pathogenic fungi within the host, pointing to a role for fungus-fungus competition in asymptomatic fungal colonization in plants.

BackgroundThe 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, the S. proliferum phylogeny and lifecycle are still an enigma.ResultsTo investigate the phylogenetic relationships between S. proliferum and other cestode species, and to examine the underlying mechanisms of pathogenicity, we sequenced the entire S. proliferum genome. Additionally, S. proliferum plerocercoid larvae transcriptome analyses were performed to identify genes involved in asexual reproduction in the host. The genome sequences confirmed that the S. proliferum genetic sequence is distinct from that of the closely related Spirometra erinaceieuropaei. Moreover, nonordinal extracellular matrix coordination allows for asexual reproduction in the host and loss of sexual maturity in S. proliferum is related to its fatal pathogenicity in humans.ConclusionsThe high-quality reference genome sequences generated should prove valuable for future studies of pseudophyllidean tapeworm biology and parasitism.

BackgroundPopulation outbreaks of the crown-of-thorns starfish (Acanthaster plancisensu lato; COTS), a primary predator of reef-building corals in the Indo-Pacific Ocean, are major concerns in coral reef management. While biological and ecological knowledge of COTS has been accumulating since the 1960s, little is known about its associated bacteria. The aim of this study was to provide fundamental information on dominant COTS-associated bacteria through a multifaceted molecular approach.MethodsA total of 205 COTS individuals from 17 locations throughout the Indo-Pacific Ocean were examined for the presence of COTS-associated bacteria. We conducted 16S rRNA metabarcoding of COTS to determine the bacterial profiles of different parts of the body, and generated a full-length 16S rRNA gene sequence from a single dominant bacterium, which we designated COTS27. We performed phylogenetic analysis to determine the taxonomy, screening of COTS27 across the Indo-Pacific, FISH to visualize it within the COTS tissues, and reconstruction of the chromosome from the hologenome sequence data.ResultsWe discovered that a single bacterium exists at high densities in the subcuticular space in COTS forming a biofilm-like structure between the cuticle and the epidermis. COTS27 belongs to a clade that presumably represents a distinct order (so-called marine spirochetes) in the phylumSpirochaetesand is universally present in COTS throughout the Indo-Pacific Ocean. The reconstructed genome of COTS27 includes some genetic traits that are probably linked to adaptation to marine environments and evolution as an extracellular endosymbiont in subcuticular spaces.ConclusionsCOTS27 can be found in three allopatrically speciated COTS species, ranging from northern Red Sea to the Pacific, implying that symbiotic relationship arose before the speciation (approximately 2 million years ago). The universal association of COTS27 with COTS and nearly mono-specific association at least with the Indo-Pacific COTS potentially provides a useful model system for studying symbiont-host interactions in marine invertebrates.

How evolutionary novelties have arisen is one of the central questions in evolutionary biology. Pre-existing gene regulatory networks or signaling pathways have been shown to be co-opted for building novel traits in several organisms. However, the structure of entire gene regulatory networks and evolutionary events of gene co-option for emergence of a novel trait are poorly understood. In this study, we used a novel wing pigmentation pattern of the polka-dotted fruit fly, and identified the complete set of genes for pigmentation pattern formation byde novogenome sequencing and transcriptome analyses. In pigmentation areas of wings, 151 genes were positively or negatively regulated bywingless, a master regulator of wing pigmentation. Genes for neural development, Wnt signaling, Dpp signaling, Zinc finger transcription factors, and effectors (such as enzymes) for melanin pigmentation were included among these 151 genes. None of the known regulatory genes that regulate pigmentation pattern formation in other fruit fly species were included. Our results suggest that the novel pigmentation pattern of the polka-dotted fruit fly emerged through multi-step co-options of multiple gene regulatory networks, signaling pathways, and effector genes, rather than recruitment of one large gene circuit.

Summary: Targeted knockin mediated by double-stranded DNA cleavage is accompanied by unwanted insertions and deletions (indels) at on-target and off-target sites. A nick-mediated approach scarcely generates indels but exhibits reduced efficiency of targeted knockin. Here, we demonstrate that tandem paired nicking, a method for targeted knockin involving two Cas9 nickases that create nicks at the homologous regions of the donor DNA and the genome in the same strand, scarcely creates indels at the edited genomic loci, while permitting the efficiency of targeted knockin largely equivalent to that of the Cas9-nuclease-based approach. Tandem paired nicking seems to accomplish targeted knockin by DNA recombination analogous to Holliday’s model and creates intended genomic changes without introducing additional nucleotide changes, such as silent mutations. Targeted knockin through tandem paired nicking neither triggers significant p53 activation nor occurs preferentially in p53-suppressed cells. These properties of tandem paired nicking demonstrate its utility in precision genome engineering. : Hyodo et al. describe a method of precise genome editing using Cas9 nickases, which barely introduces unwanted insertions and deletions at the edited genomic loci. This method neither activates nor is inhibited by p53 signalling, allowing safe genome editing and suggesting promise for this method in clinical genome engineering. Keywords: CRISPR/Cas9, nickase, tandem paired nicking, knockin, Holliday’s model, crossover, non-crossover, p53, p21, tandem nicking

Chloroplasts in photosynthetic eukaryotes originated from a cyanobacterial endosymbiosis far more than 1 billion years ago1-3. Due to this ancientness, it remains unclear how this evolutionary process proceeded. To unveil this mystery, we analysed the whole genome sequence of a photosynthetic rhizarian amoeba4, Paulinella micropora5,6, which has a chloroplast-like organelle that originated from another cyanobacterial endosymbiosis7-10 about 0.1 billion years ago11. Here we show that the predacious amoeba that engulfed cyanobacteria evolved into a photosynthetic organism very quickly in the evolutionary time scale, probably aided by the drastic genome reorganization activated by large DNA virus. In the endosymbiotic evolution of eukaryotic cells, gene transfer from the endosymbiont genome to the host nucleus is essential for the evolving host cell to control the endosymbiont-derived organelle12. In P. micropora, we found that the gene transfer from the free-living and endosymbiotic bacteria to the amoeba nucleus was rapidly activated but both simultaneously ceased within the initiation period of the endosymbiotic evolution, suggesting that the genome reorganization drastically proceeded and completed. During this period, large DNA virus appeared to have infected the amoeba, followed by the rapid amplification and diversification of virus-related genes. These findings led us to re-examine the conventional endosymbiotic evolutionary scenario that exclusively deals with the host and the symbiont, and to extend it by incorporating a third critical player, large DNA virus, which activates the drastic gene transfer and genome reorganization between them. This Paulinella version of the evolutionary hypothesis deserves further testing of its generality in evolutionary systems and could shed light on the unknown roles of large DNA viruses13 in the evolution of terrestrial life.

Ficus erecta, a wild relative of common fig (F. carica), is a donor of Ceratocystis canker resistance in fig breeding programs. Interspecific hybridization followed by recurrent backcrossing is an effective method to transfer the resistance trait from wild to cultivated fig; however, this is time consuming and labor-intensive for trees, especially for gynodioecious plants such as fig. In this study, genome resources were developed for F. erecta to facilitate fig breeding programs. The genome sequence of F. erecta was determined using single-molecule real-time sequencing technology. The resultant assembly spanned 331.6 Mb with 538 contigs and an N50 length of 1.9 Mb, from which 51,806 high-confidence genes were predicted. Pseudomolecule sequences corresponding to the chromosomes of F. erecta were established with a genetic map based on single nucleotide polymorphisms from double-digest restriction-site associated DNA sequencing. Subsequent linkage analysis and whole genome resequencing identified a candidate gene for the Ceratocystis canker resistance trait. Genome-wide genotyping analysis enabled selection of female lines that possessed resistance and effective elimination of donor genome from progeny. The genome resources provided in this study will accelerate and enhance disease resistance breeding programs in fig.

C. elegans mechanoreceptors located in the ASG sensoryneuron have been found to sense temperature — a key trait for animal survival. Experimental loss of xanthine dehydrogenase (XDH-1) function in the AIN and AVJ interneurons resulted in reduced cold tolerance and atypical neuronal response to changes in temperature. These interneurons are synapse with upstream neurons such as the mechanoreceptor-expressing ASG. Ca2+ imaging revealed that ASG responsiveness to temperature change via mechanoreceptor DEG-1, a Degenerin/Epithelial Sodium Channel (DEG/ENaC), affects downstream AIN and AVJ circuits. Ectopic expression of DEG-1 in the ASE gustatory neuron resulted in acquisition of thermosensitivity, while electrophysiological analysis revealed that DEG-1 was involved in temperature sensation. Together, these results suggest that cold tolerance is regulated by mechanoreceptor-mediated circuit calculation.

Anatomical modern humans reached East Asia by >40,000 years ago (kya). However, key questions still remain elusive with regard to the route(s) and the number of wave(s) in the dispersal into East Eurasia. Ancient genomes at the edge of East Eurasia may shed light on the detail picture of peopling to East Eurasia. Here, we analyze the whole-genome sequence of a 2.5 kya individual (IK002) characterized with a typical Jomon culture that started in the Japanese archipelago >16 kya. The phylogenetic analyses support multiple waves of migration, with IK002 forming a lineage basal to the rest of the ancient/present-day East Eurasians examined, likely to represent some of the earliest-wave migrants who went north toward East Asia from Southeast Asia. Furthermore, IK002 has the extra genetic affinity with the indigenous Taiwan aborigines, which may support a coastal route of the Jomon-ancestry migration from Southeast Asia to the Japanese archipelago. This study highlight the power of ancient genomics with the isolated population to provide new insights into complex history in East Eurasia.

SummaryIn animals, cellularization of a coenocyte is a specialized form of cytokinesis that results in the formation of a polarized epithelium during early embryonic development. It is characterized by coordinated assembly of an actomyosin network, which drives inward membrane invaginations. However, whether coordinated cellularization driven by membrane invagination exists outside animals is not known. To that end, we investigate cellularization in the ichthyosporean Sphaeroforma arctica, a close unicellular relative of animals. We show that the process of cellularization involves coordinated inward plasma membrane invaginations dependent on an actomyosin network, and reveal the temporal order of its assembly. This leads to the formation of a polarized layer of cells resembling an epithelium. We show that this epithelium-like stage is associated with tightly regulated transcriptional activation of genes involved in cell adhesion. Hereby we demonstrate the presence of a selforganized, clonally-generated, polarized layer of cells in a unicellular relative of animals.

Supplementing food with functional small molecules has been shown to prevent diseases and improve the quality of life, especially in elderly people. Citrus fruits and citrus fruit-products are popular food supplements across the world. In this study, we focused on a Japanese citrus fruit, Citrus tumida hort. ex Tanaka (C. tumida), and elucidated the effects of supplementation of the peels of immature C. tumida (PIC) on food intake, body and fat tissue weights, and metabolic profiles of plasma and liver in mice. Supplementation with 5% (w/w) PIC for 4 weeks significantly suppressed body weight gain and decreased adipose tissue weight, including that of the epididymal, perirenal, and subcutaneous fats. Metabolome analyses using capillary electrophoresis time-of-flight mass spectrometry showed that the level of 2-hydroxyvaleric acid was reduced in the blood plasma of mice fed with PIC. Supplementation with PIC significantly elevated the levels of dipeptides (Thr-Asp, Ser-Glu, and Ala-Ala), glucuronic acid (and/or galacturonic acid-2), and S-methylglutathione, and significantly reduced the levels of betaine aldehyde in the liver. In conclusion, PIC supplementation affects the metabolism of fatty acids, pectin, glutathione, and choline. Our study demonstrates the potential beneficial effects of PIC, especially in metabolic syndrome and obesity. PIC may be developed as a functional food and used in the treatment of these diseases. Nutritional and metabolome studies are effective in studying the effects of specific dietary supplements and will contribute to the development of functional foods.

SummaryFor over a thousand years throughout Asia, the common goldfish (Carassius auratus) was raised for both food and as an ornamental pet. Selective breeding over more than 500 years has created a wide array of body and pigmentation variation particularly valued by ornamental fish enthusiasts. As a very close relative of the common carp (Cyprinus carpio), goldfish shares the recent genome duplication that occurred approximately 14-16 million years ago (mya) in their common ancestor. The combination of centuries of breeding and a wide array of interesting body morphologies is an exciting opportunity to link genotype to phenotype as well as understanding the dynamics of genome evolution and speciation. Here we generated a high-quality draft sequence of a “Wakin” goldfish using 71X PacBio long-reads. We identified 70,324 coding genes and more than 11,000 non-coding transcripts. We found that the two sub-genomes in goldfish retained extensive synteny and collinearity between goldfish and zebrafish. However, “ohnologous” genes were lost quickly after the carp whole-genome duplication, and the expression of 30% of the retained duplicated gene diverged significantly across seven tissues sampled. Loss of sequence identity and/or exons determined the divergence of the expression across all tissues, while loss of conserved, non-coding elements determined expression variance between different tissues. This draft assembly also provides an important resource for comparative genomics with the very commonly used zebrafish model (Danio rerio), and for understanding the underlying genetic causes of goldfish variants.

How genetic information is modified to generate phenotypic variation within a species is one of the central questions in evolutionary biology. Here we focus on the striking intraspecific diversity of more than 200 aposematic elytral (forewing) colour patterns of the multicoloured Asian ladybird beetle,Harmonia axyridis, which is regulated by a tightly linked genetic locush. Our loss-of-function analyses, genetic association studies,de novogenome assemblies, and gene expression data reveal that the GATA transcription factor genepannieris the major regulatory gene located at thehlocus, and suggest that repeated inversions and cis-regulatory modifications atpannierled to the expansion of colour pattern variation inH. axyridis. Moreover, we show that the colour patterning function ofpannieris conserved in the seven spotted ladybird beetle,Coccinella septempunctata, suggesting thatH. axyridis’extraordinary intra-specific variation may have arisen from ancient modifications in a conserved elytral colour patterning mechanisms in ladybird beetles.

The germline mutation rate is an important parameter that affects the amount of genetic variation and the rate of evolution. However, neither the rate of germline mutations in laboratory mice nor the biological significance of the mutation rate in mammalian populations is clear. Here we studied genome-wide mutation rates and the long-term effects of mutation accumulation on phenotype in more than 20 generations of wild-type C57BL/6 mice and mutator mice, which have high DNA replication error rates. We estimated the base-substitution mutation rate to be 5.4 × 10−9 (95% confidence interval = 4.6 × 10−9–6.5 × 10−9) per nucleotide per generation in C57BL/6 laboratory mice, about half the rate reported in humans. The mutation rate in mutator mice was 17 times that in wild-type mice. Abnormal phenotypes were 4.1-fold more frequent in the mutator lines than in the wild-type lines. After several generations, the mutator mice reproduced at substantially lower rates than the controls, exhibiting low pregnancy rates, lower survival rates, and smaller litter sizes, and many of the breeding lines died out. These results provide fundamental information about mouse genetics and reveal the impact of germline mutation rates on phenotypes in a mammalian population.

In germ cells, early embryos, and stem cells of animals, PIWI-interacting RNAs (piRNAs) have an important role in silencing retrotransposons, which are vicious genomic parasites, through transcriptional and post-transcriptional mechanisms. To examine whether the piRNA pathway can be used to silence genes of interest in germ cells, we have generated knock-in mice in which a foreign DNA fragment was inserted into a region generating pachytene piRNAs. The knock-in sequence was transcribed, and the resulting RNA was processed to yield piRNAs in postnatal testes. When reporter genes possessing a sequence complementary to portions of the knock-in sequence were introduced, they were greatly repressed after the time of pachytene piRNA generation. This repression mainly occurred at the post-transcriptional level, as degradation of the reporter RNAs was accelerated. Our results show that the piRNA pathway can be used as a tool for sequence-specific gene silencing in germ cells and support the idea that the piRNA generating regions serve as traps for retrotransposons, enabling the host cell to generate piRNAs against active retrotransposons.

We report here the isolation and sequencing of 10 Y-specific tammar wallaby (Macropus eugenii) BAC clones, revealing five hitherto undescribed tammar wallaby Y genes (in addition to the five genes already described) and several pseudogenes. Some genes on the wallaby Y display testis-specific expression, but most have low widespread expression. All have partners on the tammar X, along with homologs on the human X. Nonsynonymous and synonymous substitution ratios for nine of the tammar XY gene pairs indicate that they are each under purifying selection. All 10 were also identified as being on the Y in Tasmanian devil (Sarcophilus harrisii; a distantly related Australian marsupial); however, seven have been lost from the human Y. Maximum likelihood phylogenetic analyses of the wallaby YX genes, with respective homologs from other vertebrate representatives, revealed that three marsupial Y genes (HCFC1X/Y, MECP2X/Y, and HUWE1X/Y) were members of the ancestral therian pseudoautosomal region (PAR) at the time of the marsupial/eutherian split; three XY pairs (SOX3/SRY, RBMX/Y, and ATRX/Y) were isolated from each other before the marsupial/eutherian split, and the remaining three (RPL10X/Y, PHF6X/Y, and UBA1/UBE1Y) have a more complex evolutionary history. Thus, the small marsupial Y chromosome is surprisingly rich in ancient genes that are retained in at least Australian marsupials and evolved from testis–brain expressed genes on the X.

The centromere is essential for faithful chromosome segregation by providing the site for kinetochore assembly. Although the role of the centromere is conserved throughout evolution, the DNA sequences associated with centromere regions are highly divergent among species and it remains to be determined how centromere DNA directs kinetochore formation. Despite the active use of chicken DT40 cells in studies of chromosome segregation, the sequence of the chicken centromere was unclear. Here, we performed a comprehensive analysis of chicken centromere DNA which revealed unique features of chicken centromeres compared with previously studied vertebrates. Centromere DNA sequences from the chicken macrochromosomes, with the exception of chromosome 5, contain chromosome-specific homogenous tandem repetitive arrays that span several hundred kilobases. In contrast, the centromeres of chromosomes 5, 27, and Z do not contain tandem repetitive sequences and span non-tandem-repetitive sequences of only ∼30 kb. To test the function of these centromere sequences, we conditionally removed the centromere from the Z chromosome using genetic engineering and have shown that that the non-tandem-repeat sequence of chromosome Z is a functional centromere.

Previous studies revealed thatIgf2andMpr/Igf2rare imprinted in eutherian mammals and marsupials but not in monotremes or birds.Igf2lies in a large imprinted cluster in eutherians, and its imprinting is regulated by long-range mechanisms. As a step to understand how the imprinted cluster evolved, we have determined a 490-kb chicken sequence containing the orthologs of mammalianAscl2/Mash2, Ins2andIgf2. We found that most of the genes in this region are conserved between chickens and mammals, maintaining the same transcriptional polarities and exon–intron structures. However,H19, an imprinted noncoding transcript, was absent from the chicken sequence. ChickenASCL2/CASH4andINS, the orthologs of the imprinted mammalian genes, showed biallelic expression, further supporting the notion that imprinting evolved after the divergence of mammals and birds. TheH19imprinting center and many of the local regulatory elements identified in mammals were not found in chickens. Also, a large segment of tandem repeats and retroelements identified between the two imprinted subdomains in mice was not found in chickens. Our findings show that the imprinted genes were clustered before the emergence of imprinting and that the elements associated with imprinting probably evolved after the divergence of mammals and birds.

MSM/Ms is an inbred strain derived from the Japanese wild mouse, Mus musculus molossinus. It is believed that subspecies molossinus has contributed substantially to the genome constitution of common laboratory strains of mice, although the majority of their genome is derived from the west European M. m. domesticus. Information on the molossinus genome is thus essential not only for genetic studies involving molossinus but also for characterization of common laboratory strains. Here, we report the construction of an arrayed bacterial artificial chromosome (BAC) library from male MSM/Ms genomic DNA, covering ∼1× genome equivalent. Both ends of 176,256 BAC clone inserts were sequenced, and 62,988 BAC-end sequence (BES) pairs were mapped onto the C57BL/6J genome (NCBI mouse Build 30), covering 2,228,164 kbp or 89% of the total genome. Taking advantage of the BES map data, we established a computer-based clone screening system. Comparison of the MSM/Ms and C57BL/6J sequences revealed 489,200 candidate single nucleotide polymorphisms (SNPs) in 51,137,941 bp sequenced. The overall nucleotide substitution rate was as high as 0.0096. The distribution of SNPs along the C57BL/6J genome was not uniform: The majority of the genome showed a high SNP rate, and only 5.2% of the genome showed an extremely low SNP rate (percentage identity = 0.9997); these sequences are likely derived from the molossinus genome.