Seasonal diet partition among top predators of a small island, Iriomote Island in the Ryukyu Archipelago, Japan.

Small islands tend to lack predators because species at higher trophic levels often cannot survive. However, two exceptional top predators-the Iriomote cat Prionailurus bengalensis iriomotensis, and the Crested Serpent Eagle Spilornis cheela perplexus-live on the small Iriomote Island in the Ryukyu Archipelago. To understand how these predators coexist with limited resources, we focused on their seasonal diets between which conflicts are considered to occur. To compare the diets, we used DNA metabarcoding analysis of faecal samples. In the summer, we identified 16 unique prey items from Iriomote cat faecal samples and 15 unique prey items from Crested Serpent Eagle faecal samples. In the winter, we identified 37 and 14, respectively. Using a non-metric multidimensional scaling and a permutational multivariate analysis of variance, our study reveals significant differences in the diet composition at the order level between the predators during both seasons. Furthermore, although some prey items at the species-to-order level overlapped between them, the frequency of occurrence of most prey items differed in both seasons. These results suggest that this difference in diets is one of the reasons why the Iriomote cat and the Crested Serpent Eagle are able to coexist on such a small island.

MitoFish, MitoAnnotator, and MiFish Pipeline: Updates in 10 Years.

MitoFish, MitoAnnotator, and MiFish Pipeline are comprehensive databases of fish mitochondrial genomes (mitogenomes), accurate annotation software of fish mitogenomes, and a web platform for metabarcoding analysis of fish mitochondrial environmental DNA (eDNA), respectively. The MitoFish Suite currently receives over 48,000 visits worldwide every year; however, the performance and usefulness of the online platforms can still be improved. Here, we present essential updates on these platforms, including an enrichment of the reference data sets, an enhanced searching function, substantially faster genome annotation and eDNA analysis with the denoising of sequencing errors, and a multisample comparative analysis function. These updates have made our platform more intuitive, effective, and reliable. These updated platforms are freely available at http://mitofish.aori.u-tokyo.ac.jp/.

Population admixtures in medaka inferred by multiple arbitrary amplicon sequencing.

Cost-effective genotyping can be achieved by sequencing PCR amplicons. Short 3-10 base primers can arbitrarily amplify thousands of loci using only a few primers. To improve the sequencing efficiency of the multiple arbitrary amplicon sequencing (MAAS) approach, we designed new primers and examined their efficiency in sequencing and genotyping. To demonstrate the effectiveness of our method, we applied it to examining the population structure of the small freshwater fish, medaka (Oryzias latipes). We obtained 2987 informative SNVs with no missing genotype calls for 67 individuals from 15 wild populations and three artificial strains. The estimated phylogenic and population genetic structures of the wild populations were consistent with previous studies, corroborating the accuracy of our genotyping method. We also attempted to reconstruct the genetic backgrounds of a commercial orange mutant strain, Himedaka, which has caused a genetic disturbance in wild populations. Our admixture analysis focusing on Himedaka showed that at least two wild populations had genetically been contributed to the nuclear genome of this mutant strain. Our genotyping methods and results will be useful in quantitative assessments of genetic disturbance by this commercially available strain.

Omics profiles of fecal and oral microbiota change in irritable bowel syndrome patients with diarrhea and symptom exacerbation.

BACKGROUND: Irritable bowel syndrome (IBS) is a disorder of gut-brain interaction, including dysregulation of the hypothalamic-pituitary-adrenal axis with salivary cortisol changes. However, the role of gastrointestinal microbiota during IBS symptom exacerbation remains unclear. We tested the hypothesis that the microbial species, gene transcripts, and chemical composition of fecal and oral samples are altered during the exacerbation of IBS symptoms. METHODS: Fecal, salivary, and dental plaque samples were collected at baseline from 43 men with IBS with diarrhea (IBS-D) and 40 healthy control (HC) men. Samples in the IBS-D patients were also collected during symptom exacerbation. The composition of the fecal microbiota was determined by analyzing the 16S rRNA gene, RNA-based metatranscriptome, and metabolites in samples from HC and IBS patients with and without symptom exacerbation. Oral samples were also analyzed using omics approaches. RESULTS: The fecal microbiota during IBS symptom exacerbation exhibited significant differences in the phylogenic pattern and short-chain fatty acid compared with fecal samples during defecation when symptoms were not exacerbated. Although there were no significant differences in the phylogenic pattern of fecal microbiota abundance between HCs and IBS-D patients, significant differences were detected in the expression patterns of bacterial transcriptomes related to butyrate production and neuroendocrine hormones, including tryptophan-serotonin-melatonin synthesis and glutamine/GABA. The composition of plaque microbiota was different between HC and IBS-D patients during normal defecation. CONCLUSIONS: Our findings suggest that colonic host-microbial interactions are altered in IBS-D patients during exacerbation of symptoms. There were no overlaps between feces and oral microbiomes.

Analysis of human clinical and environmental Leptospira to elucidate the eco-epidemiology of leptospirosis in Yaeyama, subtropical Japan.

BACKGROUND: Leptospirosis, a zoonosis caused by species in the spirochete genus Leptospira, is endemic to the Yaeyama region in Okinawa, subtropical Japan. Species of the P1 subclade "virulent" group, within the genus Leptospira, are the main etiological agents of leptospirosis in Okinawa. However, their environmental persistence is poorly understood. This study used a combination of bacterial isolation and environmental DNA (eDNA) metabarcoding methods to understand the eco-epidemiology of leptospirosis in this endemic region. FINDINGS: Polymerase chain reaction (PCR) characterized twelve human clinical L. interrogans isolates belonging to the P1 subclade "virulent" subgroup and 11 environmental soil isolates of the P1subclade "low virulent" subgroup (genetically related to L. kmetyi, n = 1; L. alstonii, n = 4; L. barantonii, n = 6) from the Yaeyama region targeting four virulence-related genes (lipL32, ligA, ligB and lpxD1). Clinical isolates were PCR positive for at least three targeted genes, while all environmental isolates were positive only for lipL32. Analysis of infected renal epithelial cells with selected clinical and environmental strains, revealed the disassembly of cell-cell junctions for the Hebdomadis clinical strain serogroup. Comparison of leptospiral eDNA during winter and summer identified operational taxonomic units corresponding to the species isolated from soil samples (L. kmetyi and L. barantonii) and additional P2 subclade species (L. licerasiae, L. wolffii-related, among others) that were not detected by soil cultivation. Total Leptospira read counts were higher in summer than in winter and the analysis of leptospiral/animal eDNA relationship suggested Rattus spp. as a potential reservoir animal. CONCLUSION: Our study demonstrated high environmental Leptospira diversity in the Yaeyama region, particularly during summer, when most of the leptospirosis cases are reported. In addition, several Leptospira species with pathogenic potential were identified that have not yet been reported in Yaeyama; however, the environmental persistence of P1 subclade species previously isolated from human clinical cases in this region was absent, suggesting the need of further methodology development and surveillance.

Development and evaluation of a loop-mediated isothermal amplification (LAMP) assay for quick identification of the Japanese salamander Hynobius tokyoensis.

Species identification using molecular techniques has recently become common for various taxa. Loop-mediated isothermal amplification (LAMP) is one of the easiest and least expensive molecular identification methods. Although few studies have developed LAMP assays for amphibians, we believe that LAMP is also useful for identifying endangered amphibians. Hynobius tokyoensis and H. lichenatus occur in Honshu, Japan, and have parapatric distributions. They are similar morphologically, especially at early developmental stages, including eggs and larvae. Hynobius tokyoensis has been listed as a national endangered species in Japan since 2020, and unambiguous identification of these species is therefore important for their conservation and management. In this study, we developed a LAMP primer set for the mitochondrial cytochrome b region to detect H. tokyoensis, and we evaluated the LAMP assay using total genomic DNA from four H. tokyoensis and three H. lichenatus individuals from across most of their ranges. Our LAMP primer set could distinguish these two species. This study should help to establish LAMP assays for other endangered species and morphologically similar species.

Population genetic structure and evolution of Batesian mimicry in Papilio polytes from the Ryukyu Islands, Japan, analyzed by genotyping-by-sequencing.

Batesian mimicry is a striking example of Darwinian evolution, in which a mimetic species resembles toxic or unpalatable model species, thereby receiving protection from predators. In some species exhibiting Batesian mimicry, nonmimetic individuals coexist as polymorphism in the same population despite the benefits of mimicry. In a previous study, we proposed that the abundance of mimics is limited by that of the models, leading to polymorphic Batesian mimicry in the swallowtail butterfly, Papilio polytes, on the Ryukyu Islands in Japan. We found that their mimic ratios (MRs), which varied among the Islands, were explained by the model abundance of each habitat, rather than isolation by distance or phylogenetic constraint based on the mitochondrial DNA (mtDNA) analysis. In the present study, this possibility was reexamined based on hundreds of nuclear single nucleotide polymorphisms (SNPs) of 93 P. polytes individuals from five Islands of the Ryukyus. We found that the population genetic and phylogenetic structures of P. polytes largely corresponded to the geographic arrangement of the habitat Islands, and the genetic distances among island populations show significant correlation with the geographic distances, which was not evident by the mtDNA-based analysis. A partial Mantel test controlling for the present SNP-based genetic distances revealed that the MRs of P. polytes were strongly correlated with the model abundance of each island, implying that negative frequency-dependent selection interacting with model species shaped and maintained the mimetic polymorphism. Taken together, our results support the possibility that predation pressure, not isolation by distance or other neutral factors, is a major driving force of evolution of the Batesian mimicry in P. polytes from the Ryukyus.

Understanding leptospirosis eco-epidemiology by environmental DNA metabarcoding of irrigation water from two agro-ecological regions of Sri Lanka.

BACKGROUND: Leptospirosis is one of the most significant zoonoses across the world not only because of its impact on human and animal health but also because of the economic and social impact on agrarian communities. Leptospirosis is endemic in Sri Lanka where paddy farming activities, the use of draught animals in agriculture, and peridomestic animals in urban and rural areas play important roles in maintaining the infection cycle of pathogenic Leptospira, especially concerning animals as a potential reservoir. In this study, an environmental DNA (eDNA) metabarcoding methodology was applied in two different agro-ecological regions of Sri Lanka to understand the eco-epidemiology of leptospirosis. METHODOLOGY/PRINCIPAL FINDINGS: Irrigation water samples were collected in Kandy District (wet zone mid-country region 2) and Girandurukotte, Badulla District (intermediate zone low-country region 2); and analysed for the presence of pathogenic Leptospira, associated microbiome and the potential reservoir animals. Briefly, we generated PCR products for high-throughput sequencing of multiple amplicons through next-generation sequencing. The analysis of eDNA showed different environmental microbiomes in both regions and a higher diversity of Leptospira species circulating in Kandy than in Girandurukotte. Moreover, the number of sequence reads of pathogenic Leptospira species associated with clinical cases such as L. interrogans was higher in Kandy than in Girandurukotte. Kandy also showed more animal species associated with pathogenic bacterial species than Girandurukotte. Finally, several pathogenic bacterial species including Arcobacter cryaerophilus, responsible for abortion in animals, was shown to be associated with pathogenic Leptospira. CONCLUSIONS/SIGNIFICANCE: Leptospirosis has been considered to be endemic in wet regions, consistently, leptospiral sequences were detected strongly in Kandy. The great Leptospira species diversity in Kandy observed in this study shows that the etiological agents of leptospirosis in Sri Lanka might be underestimated. Furthermore, our eDNA metabarcoding can be used to discriminate bacterial and animal species diversity in different regions and to explore environmental microbiomes to identify other associated bacterial pathogens in the environment.

Evidence for frequency-dependent selection maintaining polymorphism in the Batesian mimic Papilio polytes in multiple islands in the Ryukyus, Japan.

Batesian mimicry is a well-studied adaptation for predation avoidance, in which a mimetic species resembles an unpalatable model species. Batesian mimicry can be under positive selection because of the protection gained against predators, due to resemblance to unpalatable model species. However, in some mimetic species, nonmimetic individuals are present in populations, despite the benefits of mimicry. The mechanism for evolution of such mimetic polymorphism remains an open question. Here, we address the hypothesis that the abundance of mimics is limited by that of the models, leading to mimetic polymorphism. In addition, other forces such as the effects of common ancestry and/or isolation by distance may explain this phenomenon. To investigate this question, we focused on the butterfly, Papilio polytes, that exhibits mimetic polymorphism on multiple islands of the Ryukyus, Japan, and performed field surveys and genetic analysis. We found that the mimic ratio of P. polytes was strongly correlated with the model abundance observed on each of the five islands, suggesting negative frequency-dependent selection is driving the evolution of polymorphism in P. polytes populations. Molecular phylogenetic analysis indicated that the southern island populations are the major source of genetic diversity, and the middle and northern island populations arose by relatively recent migration. This view was also supported by mismatch distribution and Tajima's D analyses, suggesting a recent population expansion on the middle and northern islands, and stable population persistence on the southern islands. The frequency of the mimetic forms within P. polytes populations is thus explained by variations in the model abundance rather than by population structure. Thus, we propose that predation pressure, rather than neutral forces, have shaped the Batesian mimicry polymorphism in P. polytes observed in the Ryukyus.

Environmental DNA metabarcoding to detect pathogenic Leptospira and associated organisms in leptospirosis-endemic areas of Japan.

Leptospires, which cause the zoonotic disease leptospirosis, persist in soil and aqueous environments. Several factors, including rainfall, the presence of reservoir animals, and various abiotic and biotic components interact to influence leptospiral survival, persistence, and pathogenicity in the environment. However, how these factors modulate the risk of infection is poorly understood. Here we developed an approach using environmental DNA (eDNA) metabarcoding for detecting the microbiome, vertebrates, and pathogenic Leptospira in aquatic samples. Specifically, we combined 4 sets of primers to generate PCR products for high-throughput sequencing of multiple amplicons through next-generation sequencing. Using our method to analyze the eDNA of leptospirosis-endemic areas in northern Okinawa, Japan, we found that the microbiota in each river shifted over time. Operating taxonomic units corresponding to pathogenic L. alstonii, L. kmetyi, and L. interrogans were detected in association with 12 nonpathogenic bacterial species. In addition, the frequencies of 11 of these species correlated with the amount of rainfall. Furthermore, 10 vertebrate species, including Sus scrofa, Pteropus dasymallus, and Cynops ensicauda, showed high correlation with leptospiral eDNA detection. Our eDNA metabarcoding method is a powerful tool for understanding the environmental phase of Leptospira and predicting human infection risk.

MitoFish and MiFish Pipeline: A Mitochondrial Genome Database of Fish with an Analysis Pipeline for Environmental DNA Metabarcoding.

Fish mitochondrial genome (mitogenome) data form a fundamental basis for revealing vertebrate evolution and hydrosphere ecology. Here, we report recent functional updates of MitoFish, which is a database of fish mitogenomes with a precise annotation pipeline MitoAnnotator. Most importantly, we describe implementation of MiFish pipeline for metabarcoding analysis of fish mitochondrial environmental DNA, which is a fast-emerging and powerful technology in fish studies. MitoFish, MitoAnnotator, and MiFish pipeline constitute a key platform for studies of fish evolution, ecology, and conservation, and are freely available at http://mitofish.aori.u-tokyo.ac.jp/ (last accessed April 7th, 2018).

Collaborative environmental DNA sampling from petal surfaces of flowering cherry Cerasus × yedoensis 'Somei-yoshino' across the Japanese archipelago.

Recent studies have shown that environmental DNA is found almost everywhere. Flower petal surfaces are an attractive tissue to use for investigation of the dispersal of environmental DNA in nature as they are isolated from the external environment until the bud opens and only then can the petal surface accumulate environmental DNA. Here, we performed a crowdsourced experiment, the "Ohanami Project", to obtain environmental DNA samples from petal surfaces of Cerasus × yedoensis 'Somei-yoshino' across the Japanese archipelago during spring 2015. C. × yedoensis is the most popular garden cherry species in Japan and clones of this cultivar bloom simultaneously every spring. Data collection spanned almost every prefecture and totaled 577 DNA samples from 149 collaborators. Preliminary amplicon-sequencing analysis showed the rapid attachment of environmental DNA onto the petal surfaces. Notably, we found DNA of other common plant species in samples obtained from a wide distribution; this DNA likely originated from the pollen of the Japanese cedar. Our analysis supports our belief that petal surfaces after blossoming are a promising target to reveal the dynamics of environmental DNA in nature. The success of our experiment also shows that crowdsourced environmental DNA analyses have considerable value in ecological studies.

Identification of somatic mutations in postmortem human brains by whole genome sequencing and their implications for psychiatric disorders.

AIM: Somatic mutations in the human brain are hypothesized to contribute to the functional diversity of brain cells as well as the pathophysiology of neuropsychiatric diseases. However, there are still few reports on somatic mutations in non-neoplastic human brain tissues. This study attempted to unveil the landscape of somatic mutations in the human brain. METHODS: We explored the landscape of somatic mutations in human brain tissues derived from three individuals with no neuropsychiatric diseases by whole-genome deep sequencing at a depth of around 100. The candidate mutations underwent multi-layered filtering, and were validated by ultra-deep target amplicon sequencing at a depth of around 200 000. RESULTS: Thirty-one somatic mutations were identified in the human brain, demonstrating the utility of whole-genome sequencing of bulk brain tissue. The mutations were enriched in neuron-expressed genes, and two-thirds of the identified somatic single nucleotide variants in the brain tissues were cytosine-to-thymine transitions, half of which were in CpG dinucleotides. CONCLUSION: Our developed filtering and validation approaches will be useful to identify somatic mutations in the human brain. The vulnerability of neuron-expressed genes to mutational events suggests their potential relevance to neuropsychiatric diseases.

Environmental DNA enables detection of terrestrial mammals from forest pond water.

Terrestrial animals must have frequent contact with water to survive, implying that environmental DNA (eDNA) originating from those animals should be detectable from places containing water in terrestrial ecosystems. Aiming to detect the presence of terrestrial mammals using forest water samples, we applied a set of universal PCR primers (MiMammal, a modified version of fish universal primers) for metabarcoding mammalian eDNA. The versatility of MiMammal primers was tested in silico and by amplifying DNAs extracted from tissues. The results suggested that MiMammal primers are capable of amplifying and distinguishing a diverse group of mammalian species. In addition, analyses of water samples from zoo cages of mammals with known species composition suggested that MiMammal primers could successfully detect mammalian species from water samples in the field. Then, we performed an experiment to detect mammals from natural ecosystems by collecting five 500-ml water samples from ponds in two cool-temperate forests in Hokkaido, northern Japan. MiMammal amplicon libraries were constructed using eDNA extracted from water samples, and sequences generated by Illumina MiSeq were subjected to data processing and taxonomic assignment. We thereby detected multiple species of mammals common to the sampling areas, including deer (Cervus nippon), mouse (Mus musculus), vole (Myodes rufocanus), raccoon (Procyon lotor), rat (Rattus norvegicus) and shrew (Sorex unguiculatus). Many previous applications of the eDNA metabarcoding approach have been limited to aquatic/semiaquatic systems, but the results presented here show that the approach is also promising even for forest mammal biodiversity surveys.

Environmental DNA metabarcoding reveals local fish communities in a species-rich coastal sea.

Environmental DNA (eDNA) metabarcoding has emerged as a potentially powerful tool to assess aquatic community structures. However, the method has hitherto lacked field tests that evaluate its effectiveness and practical properties as a biodiversity monitoring tool. Here, we evaluated the ability of eDNA metabarcoding to reveal fish community structures in species-rich coastal waters. High-performance fish-universal primers and systematic spatial water sampling at 47 stations covering ~11 km2 revealed the fish community structure at a species resolution. The eDNA metabarcoding based on a 6-h collection of water samples detected 128 fish species, of which 62.5% (40 species) were also observed by underwater visual censuses conducted over a 14-year period. This method also detected other local fishes (≥23 species) that were not observed by the visual censuses. These eDNA metabarcoding features will enhance marine ecosystem-related research, and the method will potentially become a standard tool for surveying fish communities.

Comparison of Boiling and Robotics Automation Method in DNA Extraction for Metagenomic Sequencing of Human Oral Microbes.

The rapid improvement of next-generation sequencing performance now enables us to analyze huge sample sets with more than ten thousand specimens. However, DNA extraction can still be a limiting step in such metagenomic approaches. In this study, we analyzed human oral microbes to compare the performance of three DNA extraction methods: PowerSoil (a method widely used in this field), QIAsymphony (a robotics method), and a simple boiling method. Dental plaque was initially collected from three volunteers in the pilot study and then expanded to 12 volunteers in the follow-up study. Bacterial flora was estimated by sequencing the V4 region of 16S rRNA following species-level profiling. Our results indicate that the efficiency of PowerSoil and QIAsymphony was comparable to the boiling method. Therefore, the boiling method may be a promising alternative because of its simplicity, cost effectiveness, and short handling time. Moreover, this method was reliable for estimating bacterial species and could be used in the future to examine the correlation between oral flora and health status. Despite this, differences in the efficiency of DNA extraction for various bacterial species were observed among the three methods. Based on these findings, there is no "gold standard" for DNA extraction. In future, we suggest that the DNA extraction method should be selected on a case-by-case basis considering the aims and specimens of the study.

Rapid genome reshaping by multiple-gene loss after whole-genome duplication in teleost fish suggested by mathematical modeling.

Whole-genome duplication (WGD) is believed to be a significant source of major evolutionary innovation. Redundant genes resulting from WGD are thought to be lost or acquire new functions. However, the rates of gene loss and thus temporal process of genome reshaping after WGD remain unclear. The WGD shared by all teleost fish, one-half of all jawed vertebrates, was more recent than the two ancient WGDs that occurred before the origin of jawed vertebrates, and thus lends itself to analysis of gene loss and genome reshaping. Using a newly developed orthology identification pipeline, we inferred the post-teleost-specific WGD evolutionary histories of 6,892 protein-coding genes from nine phylogenetically representative teleost genomes on a time-calibrated tree. We found that rapid gene loss did occur in the first 60 My, with a loss of more than 70-80% of duplicated genes, and produced similar genomic gene arrangements within teleosts in that relatively short time. Mathematical modeling suggests that rapid gene loss occurred mainly by events involving simultaneous loss of multiple genes. We found that the subsequent 250 My were characterized by slow and steady loss of individual genes. Our pipeline also identified about 1,100 shared single-copy genes that are inferred to have become singletons before the divergence of clupeocephalan teleosts. Therefore, our comparative genome analysis suggests that rapid gene loss just after the WGD reshaped teleost genomes before the major divergence, and provides a useful set of marker genes for future phylogenetic analysis.

Rare variant discovery by deep whole-genome sequencing of 1,070 Japanese individuals.

The Tohoku Medical Megabank Organization reports the whole-genome sequences of 1,070 healthy Japanese individuals and construction of a Japanese population reference panel (1KJPN). Here we identify through this high-coverage sequencing (32.4 × on average), 21.2 million, including 12 million novel, single-nucleotide variants (SNVs) at an estimated false discovery rate of <1.0%. This detailed analysis detected signatures for purifying selection on regulatory elements as well as coding regions. We also catalogue structural variants, including 3.4 million insertions and deletions, and 25,923 genic copy-number variants. The 1KJPN was effective for imputing genotypes of the Japanese population genome wide. These data demonstrate the value of high-coverage sequencing for constructing population-specific variant panels, which covers 99.0% SNVs of minor allele frequency ≥0.1%, and its value for identifying causal rare variants of complex human disease phenotypes in genetic association studies.