Development of a time-series shotgun metagenomics database for monitoring microbial communities at the Pacific coast of Japan.

Although numerous metagenome, amplicon sequencing-based studies have been conducted to date to characterize marine microbial communities, relatively few have employed full metagenome shotgun sequencing to obtain a broader picture of the functional features of these marine microbial communities. Moreover, most of these studies only performed sporadic sampling, which is insufficient to understand an ecosystem comprehensively. In this study, we regularly conducted seawater sampling along the northeastern Pacific coast of Japan between March 2012 and May 2016. We collected 213 seawater samples and prepared size-based fractions to generate 454 subsets of samples for shotgun metagenome sequencing and analysis. We also determined the sequences of 16S rRNA (n = 111) and 18S rRNA (n = 47) gene amplicons from smaller sample subsets. We thereafter developed the Ocean Monitoring Database for time-series metagenomic data ( http://marine-meta.healthscience.sci.waseda.ac.jp/omd/ ), which provides a three-dimensional bird's-eye view of the data. This database includes results of digital DNA chip analysis, a novel method for estimating ocean characteristics such as water temperature from metagenomic data. Furthermore, we developed a novel classification method that includes more information about viruses than that acquired using BLAST. We further report the discovery of a large number of previously overlooked (TAG)n repeat sequences in the genomes of marine microbes. We predict that the availability of this time-series database will lead to major discoveries in marine microbiome research.

Whole-Genome Sequencing of 84 Japanese Eels Reveals Evidence against Panmixia and Support for Sympatric Speciation.

The Japanese eel (Anguilla japonica), European eel (Anguilla anguilla), and American eel (Anguilla rostrata) are migratory, catadromous, temperate zone fish sharing several common life cycle features. The population genetics of panmixia in these eel species has already been investigated. Our extensive population genetics analysis was based on 1400 Gb of whole-genome sequence (WGS) data from 84 eels. It demonstrated that a Japanese eel group from the Kuma River differed from other populations of the same species. Even after removing the potential adapted/selected single nucleotide polymorphism (SNP) data, and with very small differences (fixation index [Fst] = 0.01), we obtained results consistently indicating that panmixia does not occur in Japanese eels. The life cycle of the Japanese eel is well-established and the Kuma River is in the center of its habitat. Nevertheless, simple reproductive isolation is not the probable cause of non-panmixia in this species. We propose that the combination of spawning area subdivision, philopatry, and habitat preference/avoidance accounts for the non-panmixia in the Japanese eel population. We named this hypothesis the "reproductive isolation like subset mapping" (RISM) model. This finding may be indicative of the initial stages of sympatric speciation in these eels.

Changes in free amino acid concentrations and associated gene expression profiles in the abdominal muscle of kuruma shrimp (Marsupenaeus japonicus) acclimated at different salinities.

Shrimps inhabiting coastal waters can survive in a wide range of salinity. However, the molecular mechanisms involved in their acclimation to different environmental salinities have remained largely unknown. In the present study, we acclimated kuruma shrimp (Marsupenaeus japonicus) at 1.7%, 3.4% and 4.0% salinities. After acclimating for 6, 12, 24 and 72 h, we determined free amino acid concentrations in their abdominal muscle, and performed RNA sequencing analysis on this muscle. The concentrations of free amino acids were clearly altered depending on salinity after 24 h of acclimation. Glutamine and alanine concentrations were markedly increased following the increase of salinity. In association with such changes, many genes related to amino acid metabolism changed their expression levels. In particular, the increase of the expression level of the gene encoding glutamate-ammonia ligase, which functions in glutamine metabolism, appeared to be associated with the increased glutamine concentration at high salinity. Furthermore, the increased alanine concentration at high salinity was likely associated with the decrease in the expression levels of the the gene encoding alanine-glyoxylate transaminase. Thus, there is a possibility that changes in the concentration of free amino acids for osmoregulation in kuruma shrimp are regulated by changes in the expression levels of genes related to amino acid metabolism.

Ultrahigh-Density Linkage Map Construction Using Low-Coverage Whole-Genome Sequencing of a Doubled Haploid Population: Case Study of Torafugu (Takifugu rubripes).

Next-generation sequencing enables genome-wide genotyping of a large population and further facilitates the construction of a genetic linkage map. Low-coverage whole-genome sequencing has been employed for genetic linkage map construction in several species. However, this strategy generally requires available high-quality reference genomes and/or designed inbred pedigree lines, which restrict the scope of application for non-model and unsequenced species. Here, using torafugu (Takifugu rubripes) as a test model, we propose a new strategy for ultrahigh-density genetic linkage map construction using low-coverage whole-genome sequencing of a haploid/doubled haploid (H/DH) population without above requirements. Low-coverage (≈1×) whole-genome sequencing data of 165 DH individuals were used for de novo assembly and further performed single nucleotide polymorphisms (SNPs) calling, resulting in the identification of 1,070,601 SNPs. Based on SNP genotypes and de novo assembly, genotypes were associated with short DNA segments and an ultrahigh-density linkage map was constructed containing information of 802,277 SNPs in 3090 unique positions. Comparative analyses showed near-perfect concordance between the present linkage map and the latest published torafugu genome (FUGU5). This strategy would facilitate ultrahigh-density linkage map construction in various sexually reproducing organisms for which H/DH populations can be generated.

High-Throughput Sequencing of the Expressed Torafugu (Takifugu rubripes) Antibody Sequences Distinguishes IgM and IgT Repertoires and Reveals Evidence of Convergent Evolution.

B-cell antigen receptor (BCR) or antibody diversity arises from somatic recombination of immunoglobulin (Ig) gene segments and is concentrated within the Ig heavy (H) chain complementarity-determining region 3 (CDR-H3). We performed high-throughput sequencing of the expressed antibody heavy-chain repertoire from adult torafugu. We found that torafugu use between 70 and 82% of all possible V (variable), D (diversity), and J (joining) gene segment combinations and that they share a similar frequency distribution of these VDJ combinations. The CDR-H3 sequence repertoire observed in individuals is biased with the preferential use of a small number of VDJ, dominated by sequences containing inserted nucleotides. We uncovered the common CDR-H3 amino-acid (aa) sequences shared by individuals. Common CDR-H3 sequences feature highly convergent nucleic-acid recombination compared with private ones. Finally, we observed differences in repertoires between IgM and IgT, including the unequal usage frequencies of V gene segment and the biased number of nucleotide insertion/deletion at VDJ junction regions that leads to distinct distributions of CDR-H3 lengths.

Global gene expression analysis of the muscle tissues of medaka acclimated to low and high environmental temperatures.

Medaka (Oryzias latipes) is a temperate eurythermal fish that is able to survive over a wide range of water temperatures ranging from near zero to over 30°C throughout the year; it maintains its normal physiological and biochemical processes through temperature acclimation. To determine the mechanisms involved in temperature acclimation of fish, the fast skeletal muscle tissues of medaka underwent global gene expression analysis using next-generation sequencing. Ten individuals were placed into two aquariums at 24°C. While the water temperature of one aquarium was decreased to 10°C, that of the other aquarium was increased to 30°C; these temperatures were subsequently maintained for 5weeks. RNA sequencing (RNA-Seq) analyses revealed that 11 genes involved in energy metabolism and muscle atrophy were significantly highly expressed in the 10°C-acclimated fish. Meanwhile, significantly higher expression levels were observed for 20 genes encoding myofibrillar proteins and heat shock proteins in the 30°C-acclimated fish. Moreover, 1103 genes had at least fourfold differential expression between the acclimation groups. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses provided important information: although the expression of genes related to metabolic processes were activated, muscle atrophy occurred in the 10°C-acclimated fish, and muscle cells divided actively in the 30°C-acclimated fish and avoided thermal stress by expressing heat shock proteins. Therefore, RNA-Seq analyses with the available genome database will be useful for better understanding the molecular mechanisms involved in the temperature acclimation of fish.

Local-Level Genetic Diversity and Structure of Matsutake Mushroom (Tricholoma matsutake) Populations in Nagano Prefecture, Japan, Revealed by 15 Microsatellite Markers.

The annual yield of matsutake mushrooms (Tricholoma matsutake) has consistently decreased in Japan over the past few decades. We used 15 polymorphic and codominant simple sequence repeat (SSR) markers, developed using next-generation sequencing, to carry out genetic analyses of 10 populations in Nagano, Japan. Using the SSRs, we identified 223 genotypes, none of which was observed in more than one population. The mean expected heterozygosity and standardized allelic richness values were 0.67 and 4.05, respectively. Many alleles appeared in only one of the 10 populations; 34 of these private alleles were detected with a mean number per population of 3.4. The fixation index (FST) and standardized genetic differentiation (G'ST) values were 0.019 and 0.028, respectively. Analysis of molecular variance (AMOVA) showed that the contribution of among population, among genets within a population, and within genets variation to the total variation was 2.91%, 11.62%, and 85.47%, respectively, with genetic differentiation being detected for all sources. Twenty-eight of 45 pairwise FST values were significantly larger than zero, and no pattern of isolation by distance was detected among the 10 populations. Bayesian-based clustering did not show clear differences among populations. These results suggest that reestablishment of a colony would be best accomplished by transplantation within a field; if this is not possible, then transplantation from within several dozen kilometers will cause little damage to the original population genetic structure.

Different gene expression profiles between normal and thermally selected strains of rainbow trout, Oncorhynchus mykiss, as revealed by comprehensive transcriptome analysis.

A high-temperature selected (HT) strain of rainbow trout was established from the Donaldson (DS) strain by traditional selective breeding in Japan. The aim of this study is to investigate genes related to upper temperature tolerance in this strain utilizing next generation sequencer (NGS), and to establish comprehensive and comparable datasets in brain, liver, muscle, heart and gill tissues between the HT and DS strains. After assembling, clustering and filtering, 242,530 contigs were obtained. Among them, 7624 transcripts had at least 10 counts in expression analysis in all tissues and used as references. BLASTX homology search showed that 7329 transcripts matched with known genes. Compared to the DS strain, the HT strain expressed 90, 775, 349, 188 and 194 genes 2 folds or more in brain, liver, muscle, heart, and gill, respectively in the case of fish before heat-exposure treatment. Meanwhile, the HT strain expressed 292, 363, 433, 322 and 211 genes 2 folds or more in brain, liver, muscle, heart, and gill, respectively in case of fish after heat-exposure treatment. Many of heat shock protein family genes and transcription factor AP-1 related genes were highly expressed in all tissues of the HT strain compared with the DS strain. These results suggest that these genes play key roles in upper temperature tolerance. These comprehensive and comparable datasets will offer broad visions for upper temperature tolerance in fish species.

Characterization of the torafugu (Takifugu rubripes) immunoglobulin heavy chain gene locus.

In this study, we investigated the immunoglobulin heavy (IGH) gene locus of torafugu (Takifugu rubripes) from publicly available assembly sequences and presented an annotated locus map, including the IGHV genes, pseudogenes, and IGHC genes. Three new IGHV gene families (IGHV3-IGHV5) were discovered. We observed the interspersion of IGHV1 and IGHV2 family members and that they often intermingled with each other, while other family members were further interspersed. Conservation of the promoter and recombination signal sequences (RSS) was observed in a family-specific manner. In addition to known variable region genes present on chromosome 5 (current torafugu genome assembly), we found 34 additional IGHV genes on scaffold 287 and three novel potentially functional IGHD genes on scaffold 483. In total, the variable region of the torafugu IGH locus consists of at least 48 IGHV genes, seven IGHD genes, and six IGHJ genes. IGHC genes have also been mapped in this study, with three genes encoding immunoglobulin classes: IgT, IgM, and IgD. We confirmed the expression of newly identified IGHV3 family sequences in the spleen and kidney of adult torafugu and found a favorable IGHV segment usage by IgM and IgT. Possible structural variation in the IGHδ locus was observed based on the current torafugu assembly. The complete characterization of the torafugu IGH locus will facilitate detailed studies of large-scale mechanisms associated with the recombination of the variable region genes and will offer insights into the genetic basis of the potential diversity in the antibody response observed in torafugu.

Draft Genome Sequence of Thermoanaerobacterium saccharolyticum Strain NTOU1, a Thermophilic Bacterium Isolated from Marine Shallow Hydrothermal Vents.

Thermoanaerobacterium saccharolyticum strain NTOU1 has the ability to utilize several kinds of sugars in lignocellulosic biomass to produce ethanol more efficiently than other bacteria. Here, we report the draft genome sequence and annotation of this strain, which may provide insights into the possible genes and metabolic pathways related to ethanol production.

Dramatic improvement in genome assembly achieved using doubled-haploid genomes.

Improvement in de novo assembly of large genomes is still to be desired. Here, we improved draft genome sequence quality by employing doubled-haploid individuals. We sequenced wildtype and doubled-haploid Takifugu rubripes genomes, under the same conditions, using the Illumina platform and assembled contigs with SOAPdenovo2. We observed 5.4-fold and 2.6-fold improvement in the sizes of the N50 contig and scaffold of doubled-haploid individuals, respectively, compared to the wildtype, indicating that the use of a doubled-haploid genome aids in accurate genome analysis.

Draft Genome Sequence of Aeromonas hydrophila Strain Ae34, Isolated from a Septicemic and Moribund Koi Carp (Cyprinus carpio koi), a Freshwater Aquarium Fish.

Aeromonas hydrophila is an important opportunistic pathogen that infects a variety of aquatic and terrestrial animals, including humans. We report here the draft genome sequence of A. hydrophila Ae34, a multidrug-resistant isolate from the kidney of a moribund koi carp (Ciprinus carpio koi) with signs of hemorrhagic septicemia.

Transcriptional responses of a bicarbonate-tolerant monocot, Puccinellia tenuiflora, and a related bicarbonate-sensitive species, Poa annua, to NaHCO3 stress.

Puccinellia tenuiflora is an alkaline salt-tolerant monocot found in saline-alkali soil in China. To identify the genes which are determining the higher tolerance of P. tenuiflora compared to bicarbonate sensitive species, we examined the responses of P. tenuiflora and a related bicarbonate-sensitive Poeae plant, Poa annua, to two days of 20 mM NaHCO3 stress by RNA-seq analysis. We obtained 28 and 38 million reads for P. tenuiflora and P. annua, respectively. For each species, the reads of both unstressed and stressed samples were combined for de novo assembly of contigs. We obtained 77,329 contigs for P. tenuiflora and 115,335 contigs for P. annua. NaHCO3 stress resulted in greater than two-fold absolute expression value changes in 157 of the P. tenuiflora contigs and 1090 of P. annua contigs. Homologs of the genes involved in Fe acquisition, which are important for the survival of plants under alkaline stress, were up-regulated in P. tenuiflora and down-regulated in P. annua. The smaller number of the genes differentially regulated in P. tenuiflora suggests that the genes regulating bicarbonate tolerance are constitutively expressed in P. tenuiflora.