Massive horizontal gene transfer and the evolution of nematomorph-driven behavioral manipulation of mantids.

To complete their life cycle, a wide range of parasites must manipulate the behavior of their hosts.1 This manipulation is a well-known example of the "extended phenotype,2" where genes in one organism have phenotypic effects on another organism. Recent studies have explored the parasite genes responsible for such manipulation of host behavior, including the potential molecular mechanisms.3,4 However, little is known about how parasites have acquired the genes involved in manipulating phylogenetically distinct hosts.4 In a fascinating example of the extended phenotype, nematomorph parasites have evolved the ability to induce their terrestrial insect hosts to enter bodies of water, where the parasite then reproduces. Here, we comprehensively analyzed nematomorphs and their mantid hosts, focusing on the transcriptomic changes associated with host manipulations and sequence similarity between host and parasite genes to test molecular mimicry. The nematomorph's transcriptome changed during host manipulation, whereas no distinct changes were found in mantids. We then discovered numerous possible host-derived genes in nematomorphs, and these genes were frequently up-regulated during host manipulation. Our findings suggest a possible general role of horizontal gene transfer (HGT) in the molecular mechanisms of host manipulation, as well as in the genome evolution of manipulative parasites. The evidence of HGT between multicellular eukaryotes remains scarce but is increasing and, therefore, elucidating its mechanisms will advance our understanding of the enduring influence of HGT on the evolution of the web of life.

Interploidy gene flow involving the sexual-asexual cycle facilitates the diversification of gynogenetic triploid Carassius fish.

Asexual vertebrates are rare and at risk of extinction due to their restricted adaptability through the loss of genetic recombination. We explore the mechanisms behind the generation and maintenance of genetic diversity in triploid asexual (gynogenetic) Carassius auratus fish, which is widespread in East Asian fresh waters and exhibits one of the most extensive distribution among asexual vertebrates despite its dependence on host sperm. Our analyses of genetic composition using dozens of genetic markers and genome-wide transcriptome sequencing uncover admixed genetic composition of Japanese asexual triploid Carassius consisting of both the diverged Japanese and Eurasian alleles, suggesting the involvement of Eurasian lineages in its origin. However, coexisting sexual diploid relatives and asexual triploids in Japan show regional genetic similarity in both mitochondrial and nuclear markers. These results are attributed to a unique unidirectional gene flow from diploids to sympatric triploids, with the involvement of occasional sexual reproduction. Additionally, the asexual triploid shows a weaker population structure than the sexual diploid, and multiple triploid lineages coexist in most Japanese rivers. The generated diversity via repeated interploidy gene flow as well as an increased establishment of immigrants is assumed to offset the cost of asexual reproduction and might contribute to the successful broad distribution of this asexual vertebrate.

Ecosystem size predicts the probability of speciation in migratory freshwater fish.

Predicting speciation is a fundamental goal of research in evolutionary ecology. The probability of speciation is often positively correlated with ecosystem size. Although the mechanisms driving this correlation are generally difficult to identify, a shared geographical and ecological context provides a suitable condition to study the mechanisms that promote speciation in large ecosystems by reducing the number of factors to be considered. Here, we determined the correlation between speciation and ecosystem size, and discuss the underlying mechanisms of this relationship, using a probable parallel ecotype formation for freshwater fish. Our population genetic analysis revealed that speciation of the landlocked goby, Rhinogobius sp. YB, of the Ryukyu Archipelago, Japan, from its migratory ancestor, R. brunneus, occurred in parallel across five islands. Logistic regression analysis showed that speciation probability could be predicted using island size. The results suggest that ecosystem size predicts the occurrence of adaptation and reproductive isolation, probably through its association with three possible factors: divergent selection strength, population persistence, and occurrence probability of habitat separation.

A new targeted capture method using bacterial artificial chromosome (BAC) libraries as baits for sequencing relatively large genes.

To analyze a specific genome region using next-generation sequencing technologies, the enrichment of DNA libraries with targeted capture methods has been standardized. For enrichment of mitochondrial genome, a previous study developed an original targeted capture method that use baits constructed from long-range polymerase chain reaction (PCR) amplicons, common laboratory reagents, and equipment. In this study, a new targeted capture method is presented, that of bacterial artificial chromosome (BAC) double capture (BDC), modifying the previous method, but using BAC libraries as baits for sequencing a relatively large gene. We applied the BDC approach for the 214 kb autosomal region, ring finger protein 213, which is the susceptibility gene of moyamoya disease (MMD). To evaluate the reliability of BDC, cost and data quality were compared with those of a commercial kit. While the ratio of duplicate reads was higher, the cost was less than that of the commercial kit. The data quality was sufficiently the same as that of the kit. Thus, BDC can be an easy, low-cost, and useful method for analyzing individual genome regions with substantial length.

Using dense locality sampling resolves the subtle genetic population structure of the dispersive fish species Plecoglossus altivelis.

In dispersive species with continuous distributions, genetic differentiation between local populations is often absent or subtle and thus difficult to detect. To incorporate such subtle differentiation into management plans, it may be essential to analyse many samples from many localities using adequate numbers of high-resolution genetic markers. Here, we evaluated the usefulness of dense locality sampling in resolving genetic population structure in the ayu (Plecoglossus altivelis), a dispersive fish important in Japanese inland fisheries. Genetic variability in, and differentiation between, ayu populations around the Japan-Ryukyu Archipelago were investigated in 4746 individuals collected from 120 localities by genotyping 12 microsatellite markers. These individuals represented the two subspecies of ayu, namely the Ryukyuan subspecies (Plecoglossus altivelis ryukyuensis) and both amphidromous and landlocked forms of the nominotypical subspecies (P. a. altivelis) along the archipelago. We successfully detected an absence of genetic differentiation within the landlocked form and subtle but significant differentiation and clear geographic patterns of genetic variation among populations of the amphidromous form, which had been considered genetically homogeneous. This suggests that dense locality sampling effectively resolves subtle differences in genetic population structure, reducing stochastic deviation in the detection of genetic differentiation and geographic patterns in local populations of this dispersive species. Resampling analyses based on empirical data sets clearly demonstrate the effectiveness of increasing the number of locality samples for stable and reliable estimations of genetic fixation indices. The genetic population structure observed within the amphidromous form provides useful information for identifying management or conservation units in ayu.

Parallel Mitogenome Sequencing Alleviates Random Rooting Effect in Phylogeography.

Reliably rooted phylogenetic trees play irreplaceable roles in clarifying diversification in the patterns of species and populations. However, such trees are often unavailable in phylogeographic studies, particularly when the focus is on rapidly expanded populations that exhibit star-like trees. A fundamental bottleneck is known as the random rooting effect, where a distant outgroup tends to root an unrooted tree "randomly." We investigated whether parallel mitochondrial genome (mitogenome) sequencing alleviates this effect in phylogeography using a case study on the Sea of Japan lineage of the intertidal goby Chaenogobius annularis Eighty-three C. annularis individuals were collected and their mitogenomes were determined by high-throughput and low-cost parallel sequencing. Phylogenetic analysis of these mitogenome sequences was conducted to root the Sea of Japan lineage, which has a star-like phylogeny and had not been reliably rooted. The topologies of the bootstrap trees were investigated to determine whether the use of mitogenomes alleviated the random rooting effect. The mitogenome data successfully rooted the Sea of Japan lineage by alleviating the effect, which hindered phylogenetic analysis that used specific gene sequences. The reliable rooting of the lineage led to the discovery of a novel, northern lineage that expanded during an interglacial period with high bootstrap support. Furthermore, the finding of this lineage suggested the existence of additional glacial refugia and provided a new recent calibration point that revised the divergence time estimation between the Sea of Japan and Pacific Ocean lineages. This study illustrates the effectiveness of parallel mitogenome sequencing for solving the random rooting problem in phylogeographic studies.

Evolutionary origin of the Scombridae (tunas and mackerels): members of a paleogene adaptive radiation with 14 other pelagic fish families.

Uncertainties surrounding the evolutionary origin of the epipelagic fish family Scombridae (tunas and mackerels) are symptomatic of the difficulties in resolving suprafamilial relationships within Percomorpha, a hyperdiverse teleost radiation that contains approximately 17,000 species placed in 13 ill-defined orders and 269 families. Here we find that scombrids share a common ancestry with 14 families based on (i) bioinformatic analyses using partial mitochondrial and nuclear gene sequences from all percomorphs deposited in GenBank (10,733 sequences) and (ii) subsequent mitogenomic analysis based on 57 species from those targeted 15 families and 67 outgroup taxa. Morphological heterogeneity among these 15 families is so extraordinary that they have been placed in six different perciform suborders. However, members of the 15 families are either coastal or oceanic pelagic in their ecology with diverse modes of life, suggesting that they represent a previously undetected adaptive radiation in the pelagic realm. Time-calibrated phylogenies imply that scombrids originated from a deep-ocean ancestor and began to radiate after the end-Cretaceous when large predatory epipelagic fishes were selective victims of the Cretaceous-Paleogene mass extinction. We name this clade of open-ocean fishes containing Scombridae "Pelagia" in reference to the common habitat preference that links the 15 families.

MitoFish and MitoAnnotator: a mitochondrial genome database of fish with an accurate and automatic annotation pipeline.

Mitofish is a database of fish mitochondrial genomes (mitogenomes) that includes powerful and precise de novo annotations for mitogenome sequences. Fish occupy an important position in the evolution of vertebrates and the ecology of the hydrosphere, and mitogenomic sequence data have served as a rich source of information for resolving fish phylogenies and identifying new fish species. The importance of a mitogenomic database continues to grow at a rapid pace as massive amounts of mitogenomic data are generated with the advent of new sequencing technologies. A severe bottleneck seems likely to occur with regard to mitogenome annotation because of the overwhelming pace of data accumulation and the intrinsic difficulties in annotating sequences with degenerating transfer RNA structures, divergent start/stop codons of the coding elements, and the overlapping of adjacent elements. To ease this data backlog, we developed an annotation pipeline named MitoAnnotator. MitoAnnotator automatically annotates a fish mitogenome with a high degree of accuracy in approximately 5 min; thus, it is readily applicable to data sets of dozens of sequences. MitoFish also contains re-annotations of previously sequenced fish mitogenomes, enabling researchers to refer to them when they find annotations that are likely to be erroneous or while conducting comparative mitogenomic analyses. For users who need more information on the taxonomy, habitats, phenotypes, or life cycles of fish, MitoFish provides links to related databases. MitoFish and MitoAnnotator are freely available at http://mitofish.aori.u-tokyo.ac.jp/ (last accessed August 28, 2013); all of the data can be batch downloaded, and the annotation pipeline can be used via a web interface.

Mitogenomic sequences and evidence from unique gene rearrangements corroborate evolutionary relationships of myctophiformes (Neoteleostei).

BACKGROUND: A skewed assemblage of two epi-, meso- and bathypelagic fish families makes up the order Myctophiformes - the blackchins Neoscopelidae and the lanternfishes Myctophidae. The six rare neoscopelids show few morphological specializations whereas the divergent myctophids have evolved into about 250 species, of which many show massive abundances and wide distributions. In fact, Myctophidae is by far the most abundant fish family in the world, with plausible estimates of more than half of the oceans combined fish biomass. Myctophids possess a unique communication system of species-specific photophore patterns and traditional intrafamilial classification has been established to reflect arrangements of photophores. Myctophids present the most diverse array of larval body forms found in fishes although this attribute has both corroborated and confounded phylogenetic hypotheses based on adult morphology. No molecular phylogeny is available for Myctophiformes, despite their importance within all ocean trophic cycles, open-ocean speciation and as an important part of neoteleost divergence. This study attempts to resolve major myctophiform phylogenies from both mitogenomic sequences and corroborating evidence in the form of unique mitochondrial gene order rearrangements. RESULTS: Mitogenomic evidence from DNA sequences and unique gene orders are highly congruent concerning phylogenetic resolution on several myctophiform classification levels, corroborating evidence from osteology, larval ontogeny and photophore patterns, although the lack of larval morphological characters within the subfamily Lampanyctinae stands out. Neoscopelidae is resolved as the sister family to myctophids with Solivomer arenidens positioned as a sister taxon to the remaining neoscopelids. The enigmatic Notolychnus valdiviae is placed as a sister taxon to all other myctophids and exhibits an unusual second copy of the tRNA-Met gene - a gene order rearrangement reminiscent of that found in the tribe Diaphini although our analyses show it to be independently derived. Most tribes are resolved in accordance with adult morphology although Gonichthyini is found within a subclade of the tribe Myctophini consisting of ctenoid scaled species. Mitogenomic sequence data from this study recognize 10 reciprocally monophyletic lineages within Myctophidae, with five of these clades delimited from additional rearranged gene orders or intergenic non-coding sequences. CONCLUSIONS: Mitogenomic results from DNA sequences and unique gene orders corroborate morphology in phylogeny reconstruction and provide a likely scenario for the phylogenetic history of Myctophiformes. The extent of gene order rearrangements found within the mitochondrial genomes of myctophids is unique for phylogenetic purposes.

Hidden diversity in a reef-dwelling sea slug, Pteraeolidia ianthina (Nudibranchia, Aeolidina), in the Northwestern Pacific.

The population genetic structure of a reef-dwelling aeolid nudibranch, Pteraeolidia ianthina ( Angas, 1864 ) (Gastropoda, Mollusca), was investigated by analyzing the nuclear ribosomal RNA gene region, including 18S and 5.8S rRNA genes and the ITS1 region (545 bp). Among 235 individuals from 10 localities in the northwestern Pacific, two genetically distinct groups were detected: Groups A and B, which are separated by a minimum sequence difference of 5.3% (p-distance). The marked genetic differentiation of these two groups, even in localities where both co-occur, suggests that they are distinct species, although intergroup hybrids were observed in very low frequency. Within-population genetic variability was low in Group A, but high in Group B. Geographic analysis of between-population genetic distance and "allele" sharing data indicated three subgroups within Group A (temperate Japan, Ryukyus + Philippines, and Guam), whose distribution ranges appear to have some relationship to climate zones (temperate and subtropics + tropics) and ocean currents (Kuroshio and North Equatorial Currents). Within Group B detected only from three subtropic and tropic localities, however, genetic distances between populations were much smaller than those in Group A, and thus no significant geographic pattern was found. These differences in genetic variability and population structure between Groups A and B may reflect differences in biological features between them, such as the duration of the pelagic larval stage. Furthermore, some morphological differences were observed between Groups A and B. These observations provide further support for the suggestion that the two groups are separate species.

Low genetic differentiation between two geographically separated populations of demersal gadiform fishes in the Southern Hemisphere.

The distribution patterns of many fishes between the three continents (Africa, Australia, and South America) in the Southern Hemisphere have been uncovered to be influenced by mostly vicariance or historical dispersal. Although some demersal fishes with intercontinental distribution are suggested to be more influenced by current/recent dispersal, few genetic studies have been made for demersal fishes so far. To provide more information for such fishes, genetic divergence was analyzed for two pairs of gadiform species and subspecies distributed around Australasia and South America: the blue grenadier, Macruronus novaezelandiae (from New Zealand) and the Patagonian grenadier, M. magellanicus (from South America) as well as two subspecies of the southern blue whiting, Micromesistius australis pallidus (from New Zealand) and M. a. australis (from South America). The sequence analyses of two mitochondrial DNA regions showed no divergence between Australasian and South American populations of the grenadiers and the southern blue whiting. The microsatellite DNA analysis also indicated significant but very minimal genetic differentiation between the two geographic populations of each pair. These results imply rather recent separation of the two geographic populations. Current/recent dispersal may be an important common factor for determining the distribution of demersal fishes in the Southern Hemisphere. Nonetheless, low but significant genetic differentiation observed requires treating the two populations of the economically important grenadiers and southern blue whiting, respectively, as different stocks for proper resource management.

Ayu (Plecoglossus altivelis) in a contact zone between amphidromous and landlocked forms: genetic analyses of populations in the Yodo River system.

The migratory fish ayu, Plecoglossus altivelis, has two forms: an amphidromous form found throughout the Japanese islands and a landlocked form found in Lake Biwa. Because Lake Biwa is one of the sources of the Yodo River, which flows into Osaka Bay, it is likely that there has been contact between the two forms within the Yodo River system. To examine the effects of such contact, 573 ayu were collected from the Yodo River system and several other rivers and compared through an extensive investigation of microsatellite genotyping. Multivariate and individual admixture analyses of microsatellite data showed that samples of upstream migrants collected from the Yodo River in May 2003 included many individuals assigned to the amphidromous form, although fish collections in June and July 2003 were primarily composed of individuals of the landlocked form. Samples of drifting larvae, which we considered to be derived from eggs spawned in the Yodo River system, were mostly assigned to the amphidromous form. These results imply that an indigenous population of the amphidromous form exists and reproduces in the Yodo River system. it is likely that upstream migrants of the landlocked Lake Biwa form are supplied annually as well-developed larvae drifting out from Lake Biwa, but contribute very little to the next generation.

Isolation and characterization of 12 microsatellite loci for cross-species amplification in the cyprinid gudgeons Squalidus chankaensis and S. japonicus.

Twelve dinucleotide markers were successfully isolated and characterized from a microsatellite-enriched genomic library obtained for the gudgeon Squalidus chankaensis biwae. These markers were also available for the congeners S. c. tsuchigae and S. japonicus from Japan, which had five to 46 alleles and an expected heterozygosity ranging from zero to 0.946. Linkage equilibrium was observed at all loci, and most loci did not show significant deviation from Hardy-Weinberg equilibrium. The isolated microsatellite markers will be useful for genetic diversity studies of Squalidus populations.

Isolation and characterization of 39 microsatellite loci in the endangered Japanese loach Leptobotia curta.

A microsatellite-enriched genomic library was obtained for the endangered Japanese loach Leptobotia curta, and 39 dinucleotide markers were successfully isolated and characterized. These markers had between one and nine alleles, with expected heterozygosity ranging from 0 to 0.839, in a population from the Lake Biwa-Yodo River system of Japan. Linkage equilibrium was observed in most loci, and only one locus showed significant deviation from Hardy-Weinberg equilibrium. These microsatellite markers will be useful for genetic diversity studies of wild and captive L. curta populations.

Unexpected ceiling of genetic differentiation in the control region of the mitochondrial DNA between different subspecies of the ayu Plecoglossus altivelis.

Sequence analyses of the non-coding, control region (CR) and coding region of the ND4-tRNA(Ser) genes in the mitochondrial DNA (mtDNA) were conducted for populations of the ayu Plecoglossus altivelis altivelis and the Ryukyu-ayu P. a. ryukyuensis. The level of genetic differentiation between the two subspecies evaluated from the CR data was substantially low, when comparing with that estimated from ND4-tRNA(Ser) gene region data, as well as those from nuclear genome data sets. By contrast, the differentiation between subspecies in the ND4-tRNA(Ser) gene region was substantial, being consistent with the results from the previous nuclear genome analyses. Results of UPGMA and minimum spanning network analyses also implied the unexpected ceiling of genetic differentiation in the CR. These results suggest that the CR does not reflect accurately the level of overall genetic differentiation between the populations of the ayu, but other coding regions of the mtDNA do reflect it so that the mtDNA on the whole may function as a rich source of useful markers for genetic assessment of populations of this species.

Major patterns of higher teleostean phylogenies: a new perspective based on 100 complete mitochondrial DNA sequences.

A recent preliminary study using complete mitochondrial DNA sequences from 48 species of teleosts has suggested that higher teleostean phylogenies should be reinvestigated on the basis of more intensive taxonomic sampling. As a second step towards the resolution of higher teleostean phylogenies, which have been described as the "(unresolved) bush at the top of the tree," we reanalyzed their relationships using mitogenomic data from 100 purposefully chosen species that fully represented all of the higher teleostean orders, except for the Batrachoidiformes. Unweighted and weighted maximum parsimony analyses were conducted with the data set that comprised concatenated nucleotide sequences from 12 protein-coding genes (excluding 3rd codon positions) and 21 transfer RNA (tRNA) genes (stem regions only) from each species. The resultant trees were well resolved and largely congruent, with most internal branches being supported by high statistical values. All major, comprehensive groups above ordinal level as currently defined in higher teleosts (with the exception of the Neoteleostei and several monotypic groups), such as the Eurypterygii, Ctenosquamata, Acanthomorpha, Paracanthopterygii, Acanthopterygii, and Percomorpha, appeared to be nonmonophyletic in the present tree. Such incongruities largely resulted from differences in the placement and/or limits of the orders Ateleopodiformes, Lampridiformes, Polymixiiformes, Ophidiiformes, Lophiiformes, Beryciformes, Stephanoberyciformes, and Zeiformes, long-standing problematic taxa in systematic ichthyology. Of these, the resulting phylogenetic positions of the Ophidiiformes and Lophiiformes were totally unexpected, because, although they have consistently been considered relatively primitive groups within higher teleosts (Paracanthopterygii), they were confidently placed within a crown group of teleosts, herein called the Percomorpha. It should be noted that many unexpected, but highly supported relationships were found within the Percomorpha, being highly promising for the next investigative step towards resolution of this remarkably diversified group of teleosts.