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/.

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.

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.

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.

Estimating copy numbers of alleles from population-scale high-throughput sequencing data.

BACKGROUND: With the recent development of microarray and high-throughput sequencing (HTS) technologies, a number of studies have revealed catalogs of copy number variants (CNVs) and their association with phenotypes and complex traits. In parallel, a number of approaches to predict CNV regions and genotypes are proposed for both microarray and HTS data. However, only a few approaches focus on haplotyping of CNV loci. RESULTS: We propose a novel approach to infer copy unit alleles and their numbers in each sample simultaneously from population-scale HTS data by variational Bayesian inference on a generative probabilistic model inspired by latent Dirichlet allocation, which is a well studied model for document classification problems. In simulation studies, we evaluated concordance between inferred and true copy unit alleles for lower-, middle-, and higher-copy number dataset, in which precision and recall were ≥ 0.9 for data with mean coverage ≥ 10× per copy unit. We also applied the approach to HTS data of 1123 samples at highly variable salivary amylase gene locus and a pseudogene locus, and confirmed consistency of the estimated alleles within samples belonging to a trio of CEPH/Utah pedigree 1463 with 11 offspring. CONCLUSIONS: Our proposed approach enables detailed analysis of copy number variations, such as association study between copy unit alleles and phenotypes or biological features including human diseases.

HLA-VBSeq: accurate HLA typing at full resolution from whole-genome sequencing data.

BACKGROUND: Human leucocyte antigen (HLA) genes play an important role in determining the outcome of organ transplantation and are linked to many human diseases. Because of the diversity and polymorphisms of HLA loci, HLA typing at high resolution is challenging even with whole-genome sequencing data. RESULTS: We have developed a computational tool, HLA-VBSeq, to estimate the most probable HLA alleles at full (8-digit) resolution from whole-genome sequence data. HLA-VBSeq simultaneously optimizes read alignments to HLA allele sequences and abundance of reads on HLA alleles by variational Bayesian inference. We show the effectiveness of the proposed method over other methods through the analysis of predicting HLA types for HLA class I (HLA-A, -B and -C) and class II (HLA-DQA1,-DQB1 and -DRB1) loci from the simulation data of various depth of coverage, and real sequencing data of human trio samples. CONCLUSIONS: HLA-VBSeq is an efficient and accurate HLA typing method using high-throughput sequencing data without the need of primer design for HLA loci. Moreover, it does not assume any prior knowledge about HLA allele frequencies, and hence HLA-VBSeq is broadly applicable to human samples obtained from a genetically diverse population.

SUGAR: graphical user interface-based data refiner for high-throughput DNA sequencing.

BACKGROUND: Next-generation sequencers (NGSs) have become one of the main tools for current biology. To obtain useful insights from the NGS data, it is essential to control low-quality portions of the data affected by technical errors such as air bubbles in sequencing fluidics. RESULTS: We develop a software SUGAR (subtile-based GUI-assisted refiner) which can handle ultra-high-throughput data with user-friendly graphical user interface (GUI) and interactive analysis capability. The SUGAR generates high-resolution quality heatmaps of the flowcell, enabling users to find possible signals of technical errors during the sequencing. The sequencing data generated from the error-affected regions of a flowcell can be selectively removed by automated analysis or GUI-assisted operations implemented in the SUGAR. The automated data-cleaning function based on sequence read quality (Phred) scores was applied to a public whole human genome sequencing data and we proved the overall mapping quality was improved. CONCLUSION: The detailed data evaluation and cleaning enabled by SUGAR would reduce technical problems in sequence read mapping, improving subsequent variant analysis that require high-quality sequence data and mapping results. Therefore, the software will be especially useful to control the quality of variant calls to the low population cells, e.g., cancers, in a sample with technical errors of sequencing procedures.

TIGAR2: sensitive and accurate estimation of transcript isoform expression with longer RNA-Seq reads.

BACKGROUND: High-throughput RNA sequencing (RNA-Seq) enables quantification and identification of transcripts at single-base resolution. Recently, longer sequence reads become available thanks to the development of new types of sequencing technologies as well as improvements in chemical reagents for the Next Generation Sequencers. Although several computational methods have been proposed for quantifying gene expression levels from RNA-Seq data, they are not sufficiently optimized for longer reads (e.g. >250 bp). RESULTS: We propose TIGAR2, a statistical method for quantifying transcript isoforms from fixed and variable length RNA-Seq data. Our method models substitution, deletion, and insertion errors of sequencers based on gapped-alignments of reads to the reference cDNA sequences so that sensitive read-aligners such as Bowtie2 and BWA-MEM are effectively incorporated in our pipeline. Also, a heuristic algorithm is implemented in variational Bayesian inference for faster computation. We apply TIGAR2 to both simulation data and real data of human samples and evaluate performance of transcript quantification with TIGAR2 in comparison to existing methods. CONCLUSIONS: TIGAR2 is a sensitive and accurate tool for quantifying transcript isoform abundances from RNA-Seq data. Our method performs better than existing methods for the fixed-length reads (100 bp, 250 bp, 500 bp, and 1000 bp of both single-end and paired-end) and variable-length reads, especially for reads longer than 250 bp.

Validation of multiple single nucleotide variation calls by additional exome analysis with a semiconductor sequencer to supplement data of whole-genome sequencing of a human population.

BACKGROUND: Validation of single nucleotide variations in whole-genome sequencing is critical for studying disease-related variations in large populations. A combination of different types of next-generation sequencers for analyzing individual genomes may be an efficient means of validating multiple single nucleotide variations calls simultaneously. RESULTS: Here, we analyzed 12 independent Japanese genomes using two next-generation sequencing platforms: the Illumina HiSeq 2500 platform for whole-genome sequencing (average depth 32.4×), and the Ion Proton semiconductor sequencer for whole exome sequencing (average depth 109×). Single nucleotide polymorphism (SNP) calls based on the Illumina Human Omni 2.5-8 SNP chip data were used as the reference. We compared the variant calls for the 12 samples, and found that the concordance between the two next-generation sequencing platforms varied between 83% and 97%. CONCLUSIONS: Our results show the versatility and usefulness of the combination of exome sequencing with whole-genome sequencing in studies of human population genetics and demonstrate that combining data from multiple sequencing platforms is an efficient approach to validate and supplement SNP calls.

A statistical variant calling approach from pedigree information and local haplotyping with phase informative reads.

MOTIVATION: Variant calling from genome-wide sequencing data is essential for the analysis of disease-causing mutations and elucidation of disease mechanisms. However, variant calling in low coverage regions is difficult due to sequence read errors and mapping errors. Hence, variant calling approaches that are robust to low coverage data are demanded. RESULTS: We propose a new variant calling approach that considers pedigree information and haplotyping based on sequence reads spanning two or more heterozygous positions termed phase informative reads. In our approach, genotyping and haplotyping by the assignment of each read to a haplotype based on phase informative reads are simultaneously performed. Therefore, positions with low evidence for heterozygosity are rescued by phase informative reads, and such rescued positions contribute to haplotyping in a synergistic way. In addition, pedigree information supports more accurate haplotyping as well as genotyping, especially in low coverage regions. Although heterozygous positions are useful for haplotyping, homozygous positions are not informative and weaken the information from heterozygous positions, as majority of positions are homozygous. Thus, we introduce latent variables that determine zygosity at each position to filter out homozygous positions for haplotyping. In performance evaluation with a parent-offspring trio sequencing data, our approach outperforms existing approaches in accuracy on the agreement with single nucleotide polymorphism array genotyping results. Also, performance analysis considering distance between variants showed that the use of phase informative reads is effective for accurate variant calling, and further performance improvement is expected with longer sequencing data. CONTACT: kojima@megabank.tohoku.ac.jp .

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.

Complete genome sequence of a metallo-ß-lactamase-producing Aeromonas dhakensis strain, RYU-Ah62, isolated from a patient with an abdominal abscess.

Aeromonas dhakensis is highly virulent but often misidentified in clinical settings. The entire genome sequence of a metallo-ß-lactamase-producing A. dhakensis strain from a clinical specimen has been presented in this study. The genome comprised a single chromosome of 4.89 Mbp with 61.6% G + C content.

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 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.

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.

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.

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.

Transfer RNA gene arrangement and codon usage in vertebrate mitochondrial genomes: a new insight into gene order conservation.

BACKGROUND: Mitochondrial (mt) gene arrangement has been highly conserved among vertebrates from jawless fishes to mammals for more than 500 million years. It remains unclear, however, whether such long-term persistence is a consequence of some constraints on the gene order. RESULTS: Based on the analysis of codon usage and tRNA gene positions, we suggest that tRNA gene order of the typical vertebrate mt-genomes may be important for their translational efficiency. The vertebrate mt-genome encodes 2 rRNA, 22 tRNA, and 13 transmembrane proteins consisting mainly of hydrophobic domains. We found that the tRNA genes specifying the hydrophobic residues were positioned close to the control region (CR), where the transcription efficiency is estimated to be relatively high. Using 47 vertebrate mt-genome sequences representing jawless fishes to mammals, we further found a correlation between codon usage and tRNA gene positions, implying that highly-used tRNA genes are located close to the CR. In addition, an analysis considering the asymmetric nature of mtDNA replication suggested that the tRNA loci that remain in single-strand for a longer time tend to have more guanine and thymine not suffering deamination mutations in their anticodon sites. CONCLUSIONS: Our analyses imply the existence of translational constraint acting on the vertebrate mt-gene arrangement. Such translational constraint, together with the deamination-related constraint, may have contributed to long-term maintenance of gene order.