CTCF looping is established during gastrulation in medaka embryos.

Chromatin looping plays an important role in genome regulation. However, because ChIP-seq and loop-resolution Hi-C (DNA-DNA proximity ligation) are extremely challenging in mammalian early embryos, the developmental stage at which cohesin-mediated loops form remains unknown. Here, we study early development in medaka (the Japanese killifish, Oryzias latipes) at 12 time points before, during, and after gastrulation (the onset of cell differentiation) and characterize transcription, protein binding, and genome architecture. We find that gastrulation is associated with drastic changes in genome architecture, including the formation of the first loops between sites bound by the insulator protein CTCF and a large increase in the size of contact domains. In contrast, the binding of the CTCF is fixed throughout embryogenesis. Loops form long after genome-wide transcriptional activation, and long after domain formation seen in mouse embryos. These results suggest that, although loops may play a role in differentiation, they are not required for zygotic transcription. When we repeated our experiments in zebrafish, loops did not emerge until gastrulation, that is, well after zygotic genome activation. We observe that loop positions are highly conserved in synteny blocks of medaka and zebrafish, indicating that the 3D genome architecture has been maintained for >110-200 million years of evolution.

Bioarchaeological study of ancient Teotihuacans based on complete mitochondrial genome sequences and diet isotopes.

BACKGROUND: The Teotihuacan civilisation was the largest one in ancient Mesoamerica. The Teotihuacan city was born in the north-eastern Basin of Mexico around the second century BC, reached its peak in the fourth century AD, and had cultural influence throughout Mesoamerica. At its peak, the size of the city reached more than 20 km2, and the total population is estimated to have increased from 100,000 to 200,000. However, knowledge of the genetic background of the Teotihuacan people is still limited. AIM: We aimed to determine the mitogenome sequences of the Teotihuacan human remains and compare the ancient and present Mesoamericans. In addition, we aimed to identify the food habits of ancient Teotihuacans. SUBJECTS AND METHODS: We determined the mitogenome sequences of human remains dated to 250-636 cal AD using target enrichment-coupled next generation sequencing. We also performed stable isotope analysis. RESULTS: We successfully obtained nearly full-length sequences newly unearthed from a civilian dwelling in the Teotihuacan site. Teotihuacan mitochondrial DNA was classified into the haplogroups in present and ancient Mesoamericans. In addition, Teotihuacan individuals had a diet dependent on C4 plants such as maize. CONCLUSION: Genetic diversity varied among the Teotihuacans.

Prevotella lacticifex sp. nov., isolated from the rumen of cows.

The genus Prevotella plays an important role in polysaccharide degradation and fermentation in the rumen. To further understand the function of the phylogenetically diverse genus Prevotella, it is necessary to explore the individual characteristics at the species level. In this study, Gram-negative anaerobic bacterial strains isolated from the rumen of Holstein cows were identified. Strain R5019T was classified within the genus Prevotella based on 16S rRNA gene sequence-based phylogenetic analysis. The values of 16S rRNA gene sequence similarity, average nucleotide identity and digital DNA-DNA hybridization between strain R5019T and its phylogenetically nearest species Prevotella multisaccharivorax PPPA20T were 89.8, 82.6, and 29.3 %, respectively. The genome size of R5019T was estimated to be ca. 4.19 Mb with a genomic G+C content of 49.5 mol%. The major cellular fatty acids and menaquinones were C15 : 0 anteiso and C17 : 0 anteiso and MK-11 and MK-12, respectively. Succinate, lactate, malate, acetate and formate were produced as the fermentation end products using glucose. Based on phylogenetic, physiological, biochemical and genomic differences between 11 strains and other phylogenetically related Prevotella species, a novel species, Prevotella lacticifex sp. nov., is proposed within the genus Prevotella. The type strain is R5019T (=JCM 34664T=DSM 112675T).

Investigation of the Genetic Diversity of a Rice Core Collection of Japanese Landraces using Whole-Genome Sequencing.

The Rice Core Collection of Japanese Landraces (JRC) consisting of 50 accessions was developed by the genebank at the National Agriculture and Food Research Organization (NARO) in 2008. As a Japanese landrace core collection, the JRC has been used for many research projects, including screening for different phenotypes and allele mining for target genes. To understand the genetic diversity of Japanese Landraces, we performed whole-genome resequencing of these 50 accessions and obtained a total of 2,145,095 single nucleotide polymorphism (SNPs) and 317,832 insertion-deletions (indels) by mapping against the Oryza sativa ssp. japonica Nipponbare genome. A JRC phylogenetic tree based on 1,394 representative SNPs showed that JRC accessions were divided into two major groups and one small group. We used the multiple genome browser, TASUKE+, to examine the haplotypes of flowering genes and detected new mutations in these genes. Finally, we performed genome-wide association studies (GWAS) for agronomical traits using the JRC and another core collection, the World Rice Core Collection (WRC), comprising 69 accessions also provided by the NARO genebank. In leaf blade width, a strong peak close to NAL1, a key gene for the regulation of leaf width, and, in heading date, a peak near HESO1 involved in flowering regulation were observed in GWAS using the JRC. They were also detected in GWAS using the combined JRC + WRC. Thus, JRC and JRC + WRC are suitable populations for GWAS of particular traits.

Development of an SSR marker-based genetic linkage map and identification of a QTL associated with flowering time in Eustoma.

Lisianthus (Eustoma grandiflorum) is an important floricultural crop cultivated worldwide. Despite its commercial importance, few DNA markers are available for molecular genetic research. In this study, we constructed a genetic linkage map and to detect quantitative trait loci (QTLs) for important agronomic traits of lisianthus. To develop simple sequence repeat (SSR) markers, we used 454-pyrosequencing technology to obtain genomic shotgun sequences and subsequently identified 8263 putative SSRs. A total of 3990 primer pairs were designed in silico and 1189 unique primer pairs were extracted through a BLAST search. Amplification was successful for more than 1000 primer pairs, and ultimately 278 SSR markers exhibited polymorphism between the two lisianthus accessions evaluated. Based on these markers, a genetic linkage map was constructed using a breeding population derived from crosses between the two accessions, for which flowering time differed (>140 days when grown under 20°C). We detected one QTL associated with flowering time (phenotypic variance, 27%; LOD value, 3.7). The SSR marker located at this QTL may account for variation in flowering time among accessions (i.e., three accessions whose nodes of the first flower were over 30 had late-flowering alleles of this QTL).

Unlinking the methylome pattern from nucleotide sequence, revealed by large-scale in vivo genome engineering and methylome editing in medaka fish.

The heavily methylated vertebrate genomes are punctuated by stretches of poorly methylated DNA sequences that usually mark gene regulatory regions. It is known that the methylation state of these regions confers transcriptional control over their associated genes. Given its governance on the transcriptome, cellular functions and identity, genome-wide DNA methylation pattern is tightly regulated and evidently predefined. However, how is the methylation pattern determined in vivo remains enigmatic. Based on in silico and in vitro evidence, recent studies proposed that the regional hypomethylated state is primarily determined by local DNA sequence, e.g., high CpG density and presence of specific transcription factor binding sites. Nonetheless, the dependency of DNA methylation on nucleotide sequence has not been carefully validated in vertebrates in vivo. Herein, with the use of medaka (Oryzias latipes) as a model, the sequence dependency of DNA methylation was intensively tested in vivo. Our statistical modeling confirmed the strong statistical association between nucleotide sequence pattern and methylation state in the medaka genome. However, by manipulating the methylation state of a number of genomic sequences and reintegrating them into medaka embryos, we demonstrated that artificially conferred DNA methylation states were predominantly and robustly maintained in vivo, regardless of their sequences and endogenous states. This feature was also observed in the medaka transgene that had passed across generations. Thus, despite the observed statistical association, nucleotide sequence was unable to autonomously determine its own methylation state in medaka in vivo. Our results apparently argue against the notion of the governance on the DNA methylation by nucleotide sequence, but instead suggest the involvement of other epigenetic factors in defining and maintaining the DNA methylation landscape. Further investigation in other vertebrate models in vivo will be needed for the generalization of our observations made in medaka.

Rice Varieties in Archaic East Asia: Reduction of Its Diversity from Past to Present Times.

The Asian cultivated rice, Oryza sativa, is one of the most important crops feeding more than a third of global population. In spite of the studies for several decades, the origin and domestication history of rice varietal groups, japonica and indica, have not been fully unveiled. Genetic information of ancient rice remains is essential for direct and exclusive insight into the domestication history of rice. We performed ancient DNA analysis of 950- to 2,800-year-old rice remains excavated from Japan and Korea. We found the presence of both japonica- and indica-type varieties in the Yayoi period and the middle ages of Japan and the middle part of Korea Peninsula 2,000 years ago. It is popularly considered that japonica has been exclusively cultivated in northern part of East Asia including Japan and Korea. Our result disclosed unexpectedly wide diversity of rice varieties in archaic East Asia. The present results from ancient rice DNA reveal an exclusive insight for the domestication history of rice which is not provided as far as contemporary rice.

Imputation approach for deducing a complete mitogenome sequence from low-depth-coverage next-generation sequencing data: application to ancient remains from the Moon Pyramid, Mexico.

It is considered that more than 15 depths of coverage are necessary for next-generation sequencing (NGS) data to obtain reliable complete nucleotide sequences of the mitogenome. However, it is difficult to satisfy this requirement for all nucleotide positions because of problems obtaining a uniform depth of coverage for poorly preserved materials. Thus, we propose an imputation approach that allows a complete mitogenome sequence to be deduced from low-depth-coverage NGS data. We used different types of mitogenome data files as panels for imputation: a worldwide panel comprising all the major haplogroups, a worldwide panel comprising sequences belonging to the estimated haplogroup alone, a panel comprising sequences from the population most closely related to an individual under investigation, and a panel comprising sequences belonging to the estimated haplogroup from the population most closely related to an individual under investigation. The number of missing nucleotides was drastically reduced in all the panels, but the contents obtained by imputation were quite different among the panels. The efficiency of the imputation method differed according to the panels used. The missing nucleotides were most credibly imputed using sequences of the estimated haplogroup from the population most closely related to the individual under investigation as a panel.

Whole-genome analysis of herbicide-tolerant mutant rice generated by Agrobacterium-mediated gene targeting.

Gene targeting (GT) is a technique used to modify endogenous genes in target genomes precisely via homologous recombination (HR). Although GT plants are produced using genetic transformation techniques, if the difference between the endogenous and the modified gene is limited to point mutations, GT crops can be considered equivalent to non-genetically modified mutant crops generated by conventional mutagenesis techniques. However, it is difficult to guarantee the non-incorporation of DNA fragments from Agrobacterium in GT plants created by Agrobacterium-mediated GT despite screening with conventional Southern blot and/or PCR techniques. Here, we report a comprehensive analysis of herbicide-tolerant rice plants generated by inducing point mutations in the rice ALS gene via Agrobacterium-mediated GT. We performed genome comparative genomic hybridization (CGH) array analysis and whole-genome sequencing to evaluate the molecular composition of GT rice plants. Thus far, no integration of Agrobacterium-derived DNA fragments has been detected in GT rice plants. However, >1,000 single nucleotide polymorphisms (SNPs) and insertion/deletion (InDels) were found in GT plants. Among these mutations, 20-100 variants might have some effect on expression levels and/or protein function. Information about additive mutations should be useful in clearing out unwanted mutations by backcrossing.

TASUKE: a web-based visualization program for large-scale resequencing data.

SUMMARY: Because an enormous amount of sequence data is being collected, a method to effectively display sequence variation information is urgently needed. tasuke is a web application that visualizes large-scale resequencing data generated by next-generation sequencing technologies and is suitable for rapid data release to the public on the web. The variation and read depths of multiple genomes, as well as annotations, can be shown simultaneously at various scales. We demonstrate the use of TASUKE by applying it to 50 rice and 100 human genome resequencing datasets. AVAILABILITY AND IMPLEMENTATION: The tasuke program package and user manual are available from http://tasuke.dna.affrc.go.jp/. CONTACT: taitoh@affrc.go.jp.

Domain 3a mutation of VPS33A suppresses larval arrest phenotype in the loss of VPS45 in Caenorhabditis elegans.

The Sec1/Munc18 (SM) protein VPS45 is a key regulator of SNARE-mediated membrane fusion in endosomal trafficking, but its precise role remains unknown. To understand the function of VPS45 in vivo , we performed a genetic suppressor screen in Caenorhabditis elegans . We found that the temperature-sensitive lethality caused by the loss of VPS-45 can be suppressed by a mutation in another SM protein, VPS33A. The VPS33A M376I mutation is located in domain 3a, which is predicted to be essential for SNARE complex assembly. These results highlight the functional importance of domain 3a in endosomal SM proteins and its role in specific membrane fusion.

Plastiome: Plastisphere-enriched mobile resistome in aquatic environments.

Aquatic microplastics (MPs) act as reservoirs for microbial communities, fostering the formation of a mobile resistome encompassing diverse antibiotic (ARGs) and biocide/metal resistance genes (BMRGs), and mobile genetic elements (MGEs). This collective genetic repertoire, referred to as the "plastiome," can potentially perpetuate environmental antimicrobial resistance (AMR). Our study examining two Japanese rivers near Tokyo revealed that waterborne MPs are primarily composed of polyethylene and polypropylene fibers and sheets of diverse origin. Clinically important genera like Exiguobacterium and Eubacterium were notably enriched on MPs. Metagenomic analysis uncovered a 3.46-fold higher enrichment of ARGs on MPs than those in water, with multidrug resistance genes (MDRGs) and BMRGs prevailing, particularly within MPs. Specific ARG and BMRG subtypes linked to resistance to vancomycin, beta-lactams, biocides, arsenic, and mercury showed selective enrichment on MPs. Network analysis revealed intense associations between host genera with ARGs, BMRGs, and MGEs on MPs, emphasizing their role in coselection. In contrast, river water exhibited weaker associations. This study underscores the complex interactions shaping the mobile plastiome in aquatic environments and emphasizes the global imperative for research to comprehend and effectively control AMR within the One Health framework.

Metagenomic analyses of 7000 to 5500 years old coprolites excavated from the Torihama shell-mound site in the Japanese archipelago.

Coprolites contain various kinds of ancient DNAs derived from gut micro-organisms, viruses, and foods, which can help to determine the gut environment of ancient peoples. Their genomic information should be helpful in elucidating the interaction between hosts and microbes for thousands of years, as well as characterizing the dietary behaviors of ancient people. We performed shotgun metagenomic sequencing on four coprolites excavated from the Torihama shell-mound site in the Japanese archipelago. The coprolites were found in the layers of the Early Jomon period, corresponding stratigraphically to 7000 to 5500 years ago. After shotgun sequencing, we found that a significant number of reads showed homology with known gut microbe, viruses, and food genomes typically found in the feces of modern humans. We detected reads derived from several types of phages and their host bacteria simultaneously, suggesting the coexistence of viruses and their hosts. The food genomes provide biological evidence for the dietary behavior of the Jomon people, consistent with previous archaeological findings. These results indicate that ancient genomic analysis of coprolites is useful for understanding the gut environment and lifestyle of ancient peoples.

A calf with hind limb paralysis and dysstasia and a genome sequence analysis of an isolated Clostridium perfringens toxinotype E strain.

Clostridium perfringens toxinotype E infections are rare in calves, and the development of intestinal lesions were commonly observed. In 2012, a 6-day-old calf in Japan exhibited swelling with emphysema on the right gluteal region, sudden paralysis of the hind limb and dysstasia. A pathological examination revealed myositis of the gluteal muscle and neuritis of the ischiatic nerve. C. perfringens type E strain CP118 was isolated from the affected muscle. However, the intestinal symptoms and lesions that commonly develop in type E infections in calves were not detected in the present case. Genome analyses revealed that CP118 possessed 16 virulence-related genes, including enterotoxin, and was closely related to other type E and F strains. Particularly, CP118 was more closely related to type E strains from humans, including a food poisoning case, than calf isolates, suggesting its potential to cause food poisoning in humans and, thus, its importance as a potential risk to public health. Since CP118 did not possess the reported toxin genes associated with neuropathy, pyogenic inflammation caused by CP118 and/or other bacteria may have damaged the ischiatic nerve, resulting in neuropathy. Alternatively, unidentified CP118 toxins may have caused the neuropathy. This is the first study to report C. perfringens type E infection with peripheral neuropathy. The distribution of all the reported virulence-related genes in the C. perfringens population as well as the details of this rare case will provide further insights into C. perfringens type E infections.

Detection of non-reference porcine endogenous retrovirus loci in the Vietnamese native pig genome.

The Vietnamese native pig (VnP)-a porcine breed with a small body-has proven suitable as a biomedical animal model. Here, we demonstrate that, compared to other breeds, VnPs have fewer copies of porcine endogenous retroviruses (PERVs), which pose a risk for xenotransplantation of pig organs to humans. More specifically, we sought to characterize non-reference PERVs (nrPERVs) that were previously unidentified in the reference genome. To this end, we used whole-genome sequencing data to identify nrPERV loci with long terminal repeat (LTR) sequences in VnPs. RetroSeq was used to estimate nrPERV loci based on the most current porcine reference genome (Sscrofa11.1). LTRs were detected using de novo sequencing read assembly near the loci containing the target site duplication sequences in the inferred regions. A total of 21 non-reference LTR loci were identified and separated into two subtypes based on phylogenetic analysis. Moreover, PERVs within the detected LTR loci were identified, the presence of which was confirmed using conventional PCR and Sanger sequencing. These novel loci represent previously unknown PERVs as they have not been identified in the porcine reference genome. Thus, our RetroSeq method accurately detects novel PERV loci, and can be applied for development of a useful biomedical model.

Antimicrobial Resistance Genes in Bacteria Isolated From Japanese Honey, and Their Potential for Conferring Macrolide and Lincosamide Resistance in the American Foulbrood Pathogen Paenibacillus larvae.

American foulbrood (AFB) is the most serious bacterial disease of honey bee brood. Spores of the causative agent Paenibacillus larvae are ingested by bee larvae via brood foods and germinated cells proliferate in the larval midgut. In Japan, a macrolide antibiotic, tylosin, is used as the approved prophylactic for AFB. Although tylosin-resistant P. larvae has yet to be found in Japan, it may emerge in the future through the acquisition of macrolide resistance genes from other bacteria, and bacteria latent in brood foods, such as honey, may serve as a source of resistance genes. In this study, to investigate macrolide resistance genes in honey, we attempted to isolate tylosin-resistant bacteria from 53 Japanese honey samples and obtained 209 isolates from 48 samples in the presence of 1 µg/ml of tylosin. All isolates were Gram-positive spore-forming bacteria mainly belonging to genera Bacillus and Paenibacillus, and 94.3% exhibited lower susceptibility to tylosin than Japanese P. larvae isolates. Genome analysis of 50 representative isolates revealed the presence of putative macrolide resistance genes in the isolates, and some of them were located on mobile genetic elements (MGEs). Among the genes on MGEs, ermC on the putative mobilizable plasmid pJ18TS1mac of Oceanobacillus strain J18TS1 conferred tylosin and lincomycin resistance to P. larvae after introducing the cloned gene using the expression vector. Moreover, pJ18TS1mac was retained in the P. larvae population for a long period even under non-selective conditions. This suggests that bacteria in honey is a source of genes for conferring tylosin resistance to P. larvae; therefore, monitoring of bacteria in honey may be helpful to predict the emergence of tylosin-resistant P. larvae and prevent the selection of resistant strains.

Population dynamics in the Japanese Archipelago since the Pleistocene revealed by the complete mitochondrial genome sequences.

The Japanese Archipelago is widely covered with acidic soil made of volcanic ash, an environment which is detrimental to the preservation of ancient biomolecules. More than 10,000 Palaeolithic and Neolithic sites have been discovered nationwide, but few skeletal remains exist and preservation of DNA is poor. Despite these challenging circumstances, we succeeded in obtaining a complete mitogenome (mitochondrial genome) sequence from Palaeolithic human remains. We also obtained those of Neolithic (the hunting-gathering Jomon and the farming Yayoi cultures) remains, and over 2,000 present-day Japanese. The Palaeolithic mitogenome sequence was not found to be a direct ancestor of any of Jomon, Yayoi, and present-day Japanese people. However, it was an ancestral type of haplogroup M, a basal group of the haplogroup M. Therefore, our results indicate continuity in the maternal gene pool from the Palaeolithic to present-day Japanese. We also found that a vast increase of population size happened and has continued since the Yayoi period, characterized with paddy rice farming. It means that the cultural transition, i.e. rice agriculture, had significant impact on the demographic history of Japanese population.

Draft Genome Sequence of Brucella ceti Isolated in the Western Pacific Ocean.

In 2018, Brucella ceti was isolated from a bottlenose dolphin from the western Pacific Ocean. Here, we report a draft genome sequence of the isolate BD1442 of sequence type 27, which is the only sequence type known to have been isolated from human clinical cases.

Same Locus for Non-shattering Seed Pod in Two Independently Domesticated Legumes, Vigna angularis and Vigna unguiculata.

Loss of pod shattering is one of the most important domestication-related traits in legume crops. The non-shattering phenotypes have been achieved either by disturbed formation of abscission layer between the valves, or by loss of helical tension in sclerenchyma of endocarp, that split open the pods to disperse the seeds. During domestication, azuki bean (Vigna angularis) and yard-long bean (Vigna unguiculata cv-gr. Sesquipedalis) have reduced or lost the sclerenchyma and thus the shattering behavior of seed pods. Here we performed fine-mapping with backcrossed populations and narrowed the candidate genomic region down to 4 kbp in azuki bean and 13 kbp in yard-long bean. Among the genes located in these regions, we found MYB26 genes encoded truncated proteins in azuki bean, yard-long bean, and even cowpea. As such, our findings indicate that independent domestication on the two legumes has selected the same locus for the same traits. We also argue that MYB26 could be a target gene for improving shattering phenotype in other legumes, such as soybean.

TASUKE+: a web-based platform for exploring genome-wide association studies results and large-scale resequencing data.

Recent revolutionary advancements in sequencing technologies have made it possible to obtain mass quantities of genome-scale sequence data in a cost-effective manner and have drastically altered molecular biological studies. To utilize these sequence data, genome-wide association studies (GWASs) have become increasingly important. Hence, there is an urgent need to develop a visualization tool that enables efficient data retrieval, integration of GWAS results with diverse information and rapid public release of such large-scale genotypic and phenotypic data. We developed a web-based genome browser TASUKE+ (https://tasuke.dna.affrc.go.jp/), which is equipped with the following functions: (i) interactive GWAS results visualization with genome resequencing data and annotation information, (ii) PCR primer design, (iii) phylogenetic tree reconstruction and (iv) data sharing via the web. GWAS results can be displayed in parallel with polymorphism data, read depths and annotation information in an interactive and scalable manner. Users can design PCR primers for polymorphic sites of interest. In addition, a molecular phylogenetic tree of any region can be reconstructed so that the overall relationship among the examined genomes can be understood intuitively at a glance. All functions are implemented through user-friendly web-based interfaces so that researchers can easily share data with collaborators in remote places without extensive bioinformatics knowledge.