Introgressive hybridization in the west Pacific pen shells (genus Atrina): Restricted interspecies gene flow within the genome.

A compelling interest in marine biology is to elucidate how species boundaries between sympatric free-spawning marine invertebrates such as bivalve molluscs are maintained in the face of potential hybridization. Hybrid zones provide the natural resources for us to study the underlying genetic mechanisms of reproductive isolation between hybridizing species. Against this backdrop, we examined the occurrence of introgressive hybridization (introgression) between two bivalves distributed in the western Pacific margin, Atrina japonica and Atrina lischkeana, based on single-nucleotide polymorphisms (SNPs) derived from restriction site-associated DNA sequencing. Using 1066 ancestry-informative SNP sites, we also investigated the extent of introgression within the genome to search for SNP sites with reduced interspecies gene flow. A series of our individual-level clustering analyses including the principal component analysis, Bayesian model-based clustering, and triangle plotting based on ancestry-heterozygosity relationships for an admixed population sample from the Seto Inland Sea (Japan) consistently suggested the presence of specimens with varying degrees of genomic admixture, thereby implying that the two species are not completely isolated. The Bayesian genomic cline analysis identified 10 SNP sites with reduced introgression, each of which was located within a genic region or an intergenic region physically close to a functional gene. No, or very few, heterozygotes were observed at these sites in the hybrid zone, suggesting that selection acts against heterozygotes. Accordingly, we raised the possibility that the SNP sites are within genomic regions that are incompatible between the two species. Our finding of restricted interspecies gene flow at certain genomic regions gives new insight into the maintenance of species boundaries in hybridizing broadcast-spawning molluscs.

Accumulation of gene copy number variations during the early phase of free-spawning abalone speciation.

The genetic basis of speciation in free-spawning marine invertebrates is poorly understood. Although gene copy number variations (GCNVs) and nucleotide variations possibly trigger the speciation of these organisms, empirical evidence for such a hypothesis is limited. In this study, we searched for genomic signatures of GCNVs that may contribute to the speciation of Western Pacific abalone species. Whole-genome sequencing data suggested the existence of significant amounts of GCNVs in closely related abalones, Haliotis discus and H. madaka, in the early phase of speciation. In addition, the degree of interspecies genetic differentiation in the genes where GCNVs were estimated was higher than that in other genes, suggesting that nucleotide divergence also accumulates in the genes with GCNVs. GCNVs in some genes were also detected in other related abalone species, suggesting that these GCNVs are derived from both ancestral and de novo mutations. Our findings suggest that GCNVs have been accumulated in the early phase of free-spawning abalone speciation.

Genomic Evidence for Speciation with Gene Flow in Broadcast Spawning Marine Invertebrates.

How early stages of speciation in free-spawning marine invertebrates proceed is poorly understood. The Western Pacific abalones, Haliotis discus, H. madaka, and H. gigantea, occur in sympatry with shared breeding season and are capable of producing viable F1 hybrids in spite of being ecologically differentiated. Population genomic analyses revealed that although the three species are genetically distinct, there is evidence for historical and ongoing gene flow among these species. Evidence from demographic modeling suggests that reproductive isolation among the three species started to build in allopatry and has proceeded with gene flow, possibly driven by ecological selection. We identified 27 differentiation islands between the closely related H. discus and H. madaka characterized by high FST and dA, but not high dXY values, as well as high genetic diversity in one H. madaka population. These genomic signatures suggest differentiation driven by recent ecological divergent selection in presence of gene flow outside of the genomic islands of differentiation. The differentiation islands showed low polymorphism in H. gigantea, and both high FST, dXY, and dA values between H. discus and H. gigantea, as well as between H. madaka and H. gigantea. Collectively, the Western Pacific abalones appear to occupy the early stages speciation continuum, and the differentiation islands associated with ecological divergence among the abalones do not appear to have acted as barrier loci to gene flow in the younger divergences but appear to do so in older divergences.

Estimation of tuna population by the improved analytical pipeline of unique molecular identifier-assisted HaCeD-Seq (haplotype count from eDNA).

Many studies have investigated the ability to identify species from environmental DNA (eDNA). However, even when individual species are identified, the accurate estimation of their abundances by traditional eDNA analyses has been still difficult. We previously developed a novel analytical method called HaCeD-Seq (Haplotype Count from eDNA), which focuses on the mitochondrial D-loop sequence. The D-loop is a rapidly evolving sequence and has been used to estimate the abundance of eel species in breeding water. In the current study, we have further improved this method by applying unique molecular identifier (UMI) tags, which eliminate the PCR and sequencing errors and extend the detection range by an order of magnitude. Based on this improved HaCeD-Seq pipeline, we computed the abundance of Pacific bluefin tuna (Thunnus orientalis) in aquarium tanks at the Tokyo Sea Life Park (Kasai, Tokyo, Japan). This tuna species is commercially important but is at high risk of resource depletion. With the developed UMI tag method, 90 out of 96 haplotypes (94%) were successfully detected from Pacific bluefin tuna eDNA. By contrast, only 29 out of 96 haplotypes (30%) were detected when UMI tags were not used. Our findings indicate the potential for conducting non-invasive fish stock surveys by sampling eDNA.

Microsatellite DNA markers applicable to paternity inference in the androdioecious gooseneck barnacle Octolasmis warwickii (Lepadiformes: Poecilasmatidae).

The gooseneck barnacle Octolasmis warwickii has a rare sexual system called androdioecy, in which hermaphrodites and dwarf males co-occur. It has been hypothesized that dwarf males can coexist with conspecific hermaphrodites when dwarf males are capable of leaving more offspring than hermaphrodites via male reproduction. This hypothesis of reproductive superiority of dwarf males can be validated by comparing the reproductive success between dwarf males and hermaphrodites through DNA marker-based parentage testing. In the present study, we developed microsatellite DNA markers for O. warwickii, and evaluated the power of these markers to infer parentage based on simulation analysis. Using next generation sequencing, we obtained 344 microsatellite sequences suitable for designing primer sets for amplification in polymerase chain reaction (PCR). Of these, we examined the PCR amplification efficiency of 54 primer sets, of which 11 passed our primer screening in a population sample (n = 35). The developed markers exhibited moderate to high levels of polymorphisms, and met Hardy-Weinberg equilibrium with little evidence of significant allelic association to each other. Our simulated paternity inference suggested that the combinational use of the markers allows a high resolution of parentage (success rate of > 99.9%) if all candidate fathers are available.

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.

Development and characterization of 19 novel microsatellite markers in the Pacific seaweed pipefish Syngnathus schlegeli using next-generation sequencing.

Syngnathids (pipefishes, seahorses and seadragons) are vulnerable to human-mediated habitat perturbation. The Pacific seaweed pipefish Syngnathus schlegeli has a large distribution in the northwestern Pacific, where deterioration, loss and fragmentation of its seagrass habitat are occurring through coastal development. So far, few studies have been conducted to access the genetic structure and conservation status of S. schlegeli because of the low number of genetic markers currently described. Nineteen polymorphic microsatellite markers were developed for S. schlegeli using next-generation sequencing, and characterized in 32 individuals. The mean number of alleles was 14, with 2-28 alleles per locus. The estimates of observed heterozygosity (HO) and expected heterozygosity (HE) varied depending on the locus, ranging from 0.063 to 1.000, and from 0.062 to 0.969, respectively. Seventeen of the 19 microsatellites conformed to Hardy-Weinberg equilibrium. These new microsatellite markers should provide a wealth of information for studies on conservation genetics and the behavioral ecology of S. schlegeli.

Novel full-spectral flow cytometry with multiple spectrally-adjacent fluorescent proteins and fluorochromes and visualization of in vivo cellular movement.

Flow cytometric analysis with multicolor fluoroprobes is an essential method for detecting biological signatures of cells. Here, we present a new full-spectral flow cytometer (spectral-FCM). Unlike conventional flow cytometer, this spectral-FCM acquires the emitted fluorescence for all probes across the full-spectrum from each cell with 32 channels sequential PMT unit after dispersion with prism, and extracts the signals of each fluoroprobe based on the spectral shape of each fluoroprobe using unique algorithm in high speed, high sensitive, accurate, automatic and real-time. The spectral-FCM detects the continuous changes in emission spectra from green to red of the photoconvertible protein, KikGR with high-spectral resolution and separates spectrally-adjacent fluoroprobes, such as FITC (Emission peak (Em) 519 nm) and EGFP (Em 507 nm). Moreover, the spectral-FCM can measure and subtract autofluorescence of each cell providing increased signal-to-noise ratios and improved resolution of dim samples, which leads to a transformative technology for investigation of single cell state and function. These advances make it possible to perform 11-color fluorescence analysis to visualize movement of multilinage immune cells by using KikGR-expressing mice. Thus, the novel spectral flow cytometry improves the combinational use of spectrally-adjacent various FPs and multicolor fluorochromes in metabolically active cell for the investigation of not only the immune system but also other research and clinical fields of use.

Tetra-repeat microsatellite markers for the masu salmon (Oncorhynchus masou masou) and its application in cross-subspecies amplification.

We developed tetranucleotide-repeat microsatellite markers for the masu salmon (Oncorhynchus masou) complex. 454 pyrosequencing was used to discover repeat motifs, and seven polymorphic microsatellite-primer sets were identified. The number of alleles detected at each locus ranged from four to 24 and the expected heterozygosity varied from 0.57 to 0.92. Cross-subspecies amplification for O. m. masou, O. m. ishikawae and O. m. subsp. was successful. These microsatellites can be utilized in studies of genetic structure, genetic diversity, and intra- and inter-subspecific hybridization, making a contribution to conservation and management of the Oncorhynchus masou complex.

Complete mitochondrial DNA sequence of the ark shell Scapharca broughtonii: an ultra-large metazoan mitochondrial genome.

The complete mitochondrial (mt) genome of the ark shell Scapharca broughtonii was determined using long PCR and a genome walking sequencing strategy with genus-specific primers. The S. broughtonii mt genome (GenBank accession number AB729113) contained 12 protein-coding genes (the atp8 gene is missing, as in most bivalves), 2 ribosomal RNA genes, and 42 transfer tRNA genes, in a length of 46,985 nucleotides for the size of mtDNA with only one copy of the heteroplasmic tandem repeat (HTR) unit. Moreover the S. broughtonii mt genome shows size variation; these genomes ranged in size from about 47 kb to about 50 kb because of variation in the number of repeat sequences in the non-coding region. The mt-genome of S. broughtonii is, to date, the longest reported metazoan mtDNA sequence. Sequence duplication in non-coding region and the formation of HTR arrays were two of the factors responsible for the ultra-large size of this mt genome. All the tRNA genes were found within the S. broughtonii mt genome, unlike the other bivalves usually lacking one or more tRNA genes. Twelve additional specimens were used to analyze the patterns of tandem repeat arrays by PCR amplification and agarose electrophoresis. Each of the 12 specimens displayed extensive heteroplasmy and had 8-10 length variants. The motifs of the HTR arrays are about 353-362 bp and the number of repeats ranges from 1 to 11.

Coevolution of interacting fertilization proteins.

Reproductive proteins are among the fastest evolving in the proteome, often due to the consequences of positive selection, and their rapid evolution is frequently attributed to a coevolutionary process between interacting female and male proteins. Such a process could leave characteristic signatures at coevolving genes. One signature of coevolution, predicted by sexual selection theory, is an association of alleles between the two genes. Another predicted signature is a correlation of evolutionary rates during divergence due to compensatory evolution. We studied female-male coevolution in the abalone by resequencing sperm lysin and its interacting egg coat protein, VERL, in populations of two species. As predicted, we found intergenic linkage disequilibrium between lysin and VERL, despite our demonstration that they are not physically linked. This finding supports a central prediction of sexual selection using actual genotypes, that of an association between a male trait and its female preference locus. We also created a novel likelihood method to show that lysin and VERL have experienced correlated rates of evolution. These two signatures of coevolution can provide statistical rigor to hypotheses of coevolution and could be exploited for identifying coevolving proteins a priori. We also present polymorphism-based evidence for positive selection and implicate recent selective events at the specific structural regions of lysin and VERL responsible for their species-specific interaction. Finally, we observed deep subdivision between VERL alleles in one species, which matches a theoretical prediction of sexual conflict. Thus, abalone fertilization proteins illustrate how coevolution can lead to reproductive barriers and potentially drive speciation.

Genetic differences between hatchery stocks and natural populations in Pacific abalone (Haliotis discus) estimated using microsatellite dNA markers.

Genetic variations within and between nine hatchery stocks and seven natural populations of abalone including Ezo-abalone (Haliotis discus hannai) and Kuro-abalone (H. d. discus) were assayed with nine microsatellite markers. Marked reductions of genetic variability in the hatchery stocks were recognized in the allelic diversity and mean heterozygosity compared with the natural populations. Thirteen of 16 significant HWE deviations in hatchery stocks revealed heterozygotes excess, while all natural populations did not show such a tendency. Highly significant F (ST) values were observed for all cases between the hatchery stocks, and between the hatchery stocks and natural populations. Genetic distance (D (A)) between each hatchery stock and the geographically proximal population (mean +/- SD, 0.108 +/- 0.035) were similar to those estimated for between the natural Ezo-abalone and Kuro-abalone (0.101 +/- 0.021). The self-assignment test, which allocated individuals to their own stock with a high success rate, provided evidence of solid genetic differences among the nine hatchery stocks. These results suggests that the allelic composition and diversity in the natural populations was not necessarily reflected in the hatchery stocks owing to population bottleneck and genetic drift through seedling process, and thus the seedling and stocking practice of these hatchery stocks should take much notice of the results to conserve the genetic diversity of natural populations.

Linkage maps for the Pacific abalone (genus Haliotis) based on microsatellite DNA markers.

This study presents linkage maps for the Pacific abalone (Haliotis discus hannai) based on 180 microsatellite DNA markers. Linkage mapping was performed using three F1 outbred families, and a composite linkage map for each sex was generated by incorporating map information from the multiple families. A total of 160 markers are placed on the consolidated female map and 167 markers on the male map. The numbers of linkage groups in the composite female and male maps are 19 and 18, respectively; however, by aligning the two maps, 18 linkage groups are formed, which are consistent with the haploid chromosome number of H. discus hannai. The female map spans 888.1 cM (Kosambi) with an average spacing of 6.3 cM; the male map spans 702.4 cM with an average spacing of 4.7 cM. However, we encountered several linkage groups that show a high level of heterogeneity in recombination rate between families even within the same sex, which reduces the precision of the consolidated maps. Nevertheless, we suggest that the composite maps are of significant potential use as a scaffold to further extend the coverage of the H. discus hannai genome with additional markers.

Segregation and linkage analysis of 75 novel microsatellite DNA markers in pair crosses of Japanese abalone (Haliotis discus hannai) using the 5'-tailed primer method.

We present novel microsatellite markers of the Japanese abalone (Haliotis discus hannai) for general mapping studies in this species. A total of 75 microsatellite markers were developed, and the allele-transmission patterns of these markers were studied in three families generated by pair crosses. For allele scoring, we employed the 5'-tailed primer polymerase chain reaction (PCR) technique, which substantially reduces the cost for fluorescent labeling of primers. Of the 225 possible marker-family combinations (75 markers x 3 families), 18 cases of informative null-allele segregation were inferred. When such null-allele segregations were allowed, more than 70% of the 75 markers in the families turned out to be markers with an expected segregation ratio of 1:1:1:1, allowing maximal exploitation of the codominant nature of microsatellite markers. There were 16 instances of segregation distortion at the 5% significance level. The test for independence of segregation assigned the 75 markers into 17 linkage groups, which is in close agreement with the haploid chromosome number of H. discus hannai (n = 18). Six markers could not be placed into any linkage group. We suggest that these markers could help construct a H. discus hannai linkage map.

Development of novel microsatellite DNA markers from the Pacific oyster Crassostrea gigas.

We document the potential of novel microsatellites as a genetic tool in furthering our understanding of the Crassostrea gigas genetic structure. From the microsatellite-enriched libraries we constructed, 123 repeat regions that had sufficient sequence information to design polymerase chain reaction primer sets were isolated. From these, 9 primer pairs were screened in a C. gigas population of 67 individuals to evaluate the genetic variability. All but 1 of the 9 loci showed allelic variation (number of alleles, 2-20; observed heterozygosity, 0.119-0.925; unbiased expected heterozygosity, 0.139-0.914). Considerable discrepancy of genotypic proportions from the Hardy-Weinberg equilibrium was observed at 1 locus with an apparent heterozygote deficiency. Several loci were successfully amplified in 3 other related species with the appropriate allele size: 6 loci in C. sikamea, 4 loci in C. ariakensis, and 5 loci in C. nippona.