Evidence for hybridization-driven heteroplasmy maintained across generations in a ricefish endemic to a Wallacean ancient lake.

Heteroplasmy, the presence of multiple mitochondrial DNA (mtDNA) haplotypes within cells of an individual, is caused by mutation or paternal leakage. However, heteroplasmy is usually resolved to homoplasmy within a few generations because of germ-line bottlenecks; therefore, instances of heteroplasmy are limited in nature. Here, we report heteroplasmy in the ricefish species Oryzias matanensis, endemic to Lake Matano, an ancient lake in Sulawesi Island, in which one individual was known to have many heterozygous sites in the mitochondrial NADH dehydrogenase subunit 2 (ND2) gene. In this study, we cloned the ND2 gene for some additional individuals with heterozygous sites and demonstrated that they are truly heteroplasmic. Phylogenetic analysis revealed that the extra haplotype within the heteroplasmic O. matanensis individuals clustered with haplotypes of O. marmoratus, a congeneric species inhabiting adjacent lakes. This indicated that the heteroplasmy originated from paternal leakage due to interspecific hybridization. The extra haplotype was unique and contained two non-synonymous substitutions. These findings demonstrate that this hybridization-driven heteroplasmy was maintained across generations for a long time to the extent that the extra mitochondria evolved within the new host.

Ancient Admixture in Freshwater Halfbeaks of the Genus Nomorhamphus in Southeast Sulawesi.

Freshwater halfbeaks of the genus Nomorhamphus (Zenarchopteridae) uniquely diversified on Sulawesi Island, where tectonic movements have been very active since the Pliocene. Most species of this genus have quite limited distributions, which indicates that geographic isolations have contributed to their diversification. In this study, we demonstrated that secondary contacts and resultant admixtures between long-isolated species/populations may have also been important. We found that the mitochondrial phylogeny of a group of Nomorhamphus in Southeast Sulawesi was discordant with the nuclear phylogeny. Most notably, individuals in the upper and lower streams of the Moramo River, a small river in this region, clustered with each other in the mitochondrial phylogeny but not in the nuclear phylogeny; in the latter, the lower-stream individuals formed a clade with individuals in the Anduna River, a different river with no present water connection to the Moramo River. Phylogenetic network and population structure analyses using genomic data obtained from RNA-seq revealed that the lower-stream Moramo population admixed with the upper-stream Moramo lineage in ancient times. These findings indicate that the observed mito-nuclear discordance was caused by mitochondrial introgression and not incomplete lineage sorting. The phylogenetic network also revealed several other admixtures between ancient lineages. Repeated admixtures were also evidenced by topological incongruence in population trees estimated using the RNA-seq data. We propose that activities of many fault systems dissecting Southeast Sulawesi caused repeated secondary contact.

Species divergence and repeated ancient hybridization in a Sulawesian lake system.

An increasing volume of empirical studies demonstrated that hybridization between distant lineages may have promoted speciation in various taxa. However, the timing, extent and direction of introgressive hybridization remain unknown in many cases. Here, we report a possible case in which repeated hybridization promoted divergence of Oryzias ricefishes (Adrianichthyidae) on Sulawesi, an island of Wallacea. Four Oryzias species are endemic to the Malili Lake system in central Sulawesi, which is composed of five tectonic lakes; of these, one lake is inhabited by two species. Morphological and population genomic analyses of genome-wide single-nucleotide polymorphisms revealed that these two sympatric species are phylogenetically sister to but substantially reproductively isolated from each other. Analyses of admixture and comparison of demographic models revealed that the two sympatric species experienced several substantial introgressions from outgroup populations that probably occurred soon after they had secondary contact with each other in the lake. However, the ratio of migrants from the outgroups was estimated to be different between the two species, which is consistent with the hypothesis that these introgressions aided their divergence or prevented them from forming a hybrid swarm. Repeated lake fragmentations and fusions may have promoted diversification of this freshwater fish species complex that is endemic to this ancient lake system.

Discrepancy in the degree of population differentiation between color-morph frequencies and neutral genetic loci in the damselfly Ischnura senegalensis in Okinawa Island, Japan.

Evaluation of relative contribution of natural selection and stochastic processes to population differentiation has been of great interest in evolutionary biology. In a damselfly, Ischnura senegalensis, females show color dimorphism (gynochrome vs. androchrome), and color-morph frequencies are known to greatly vary among local populations within Okinawa Island, a small island of Ryukyu Archipelago, Japan. In this study, to examine the effects of natural selection and stochastic processes on the within-island variation in color-morph frequencies, we compared the degree of population differentiation at the color-morph locus with that at a mitochondrial DNA region and ten nuclear microsatellite loci. F ST values at the neutral loci were close to zero, indicating presence of sufficient gene flow (dispersal of adult individuals) between the local populations. In contrast, F ST values at the color-morph locus were significantly different from zero. These results suggest that variation in female color-morph frequencies observed among local populations in Okinawa Island has been caused by divergent selection acting on the phenotype and/or genes tightly linked with the color locus.

A comparative study of natural variation in hemolymph glucose levels under different dietary sugar conditions in Drosophila melanogaster and D. simulans.

Physiological responses to environmental changes play important roles in adaptive evolution. In particular, homeostatic regulatory systems that maintain constant circulating glucose levels are crucial in animals. However, variation in circulating glucose levels and the genetic effects on phenotypic variation in natural populations remain to be clarified. Here, we investigated the hemolymph glucose levels in natural populations of Drosophila melanogaster and its sibling species, D. simulans, in Japan. We quantified hemolymph glucose concentrations in third instar larvae of 27 lines for each species, which were reared on either glucose-free or glucose-rich food. In both species, genetic variation was not a major component of phenotypic variation on either glucose-free or glucose-rich food. The hemolymph glucose concentrations were much higher in D. simulans than in D. melanogaster. Genetic variance was larger in D. simulans than in D. melanogaster. The observed differences between the two species may be associated with the much more recent colonization history of D. simulans populations in Japan and/or the tolerance to environmental stresses. Our findings suggest that natural selection acting on hemolymph glucose levels in D. melanogaster is different from that in D. simulans.

Speciation history of three closely related oak gall wasps, Andricus mukaigawae, A. kashiwaphilus, and A. pseudoflos (Hymenoptera: Cynipidae) inferred from nuclear and mitochondrial DNA sequences.

The Andricus mukaigawae complex of oak gall wasps is composed of cyclically parthenogenetic species: A. mukaigawae and Andricus kashiwaphilus, and a parthenogenetic species, Andricus pseudoflos. The component species differ in life history, host plant, karyotype, and asexual gall shape, although little difference is found in the external morphology of asexual adults. To understand the speciation history of this species complex, DNA sequences of one mitochondrial region and nine nuclear gene regions were investigated. The genetic relationship among the species suggested that a loss of sex occurred after host shift. Unexpectedly, two or three distinct groups in the parthenogenetic species, A. pseudoflos, were revealed by both mitochondrial and nuclear DNA data. Gene flow in nuclear genes from the species not infected by Wolbachia (A. kashiwaphilus) to the species infected by it (A. mukaigawae) was suggested by a method based on coalescent simulations. On the other hand, gene flow in mitochondrial genes was suggested to be in the opposite direction. These findings indicate possible involvement of Wolbachia infection in the speciation process of the A. mukaigawae complex.

Seasonal changes in the long-distance linkage disequilibrium in Drosophila melanogaster.

Seasonal environmental changes have the potential to influence the genetic structure of species with a short generation time, such as Drosophila. We previously found the seasonal change in linkage disequilibrium (LD) between the chemoreceptor (Cr) genes in a local Japanese population (Kyoto [KY]). This could be caused by fluctuation in the population size or selection in temporally heterogeneous environments or both. Here, we analyzed the scale of LD between 51 X-linked polymorphisms (10 Cr and 41 non-Cr gene markers) in the 2 seasonal samples from the KY population and an autumn sample from 106 localities in and around Japan (Ja03au). Many of the non-Cr genes have receptor function but fewer functional connections to each other. The magnitude of LD in Ja03au did not significantly differ from that in the KY autumn sample. The lack of local differentiation was confirmed in an autumn sample from another local Japanese population. On the other hand, the magnitude of LD was significantly larger in spring than in autumn in the 2 independent KY samples. This suggests that reduction in the population size during winter increased the magnitude of LD in spring in the mainland population in Japan. Long-distance LD could be a useful measure for assessing seasonal fluctuation in effective population size.

Levels and patterns of DNA variation in two sympatric mangrove species, Rhizophora apiculata and R. mucronata from Thailand.

In mangrove species the past geomorphic changes in coastal regions and reproductive systems are important factors of their distribution and genetic structure of populations. However, very little is known about the levels of genetic variation of Rhiozophora species in Southeast Asia. In this study, we surveyed levels and patterns of genetic variation as well as population structure of two sympatric mangrove species, Rhizophora apiculata and R. mucronata in Thailand, using five nuclear genes and two cpDNA regions. In all investigated DNA regions, nucleotide variation within species was low, while nucleotide divergence between the two species was considerable. The nuclear genes evolved 10 times faster than the cpDNA regions. In both R. apiculata and R. mucronata, significant positive F(IS) values were found, indicating deviation from Hardy-Weinberg proportions and a deficiency of heterozygotes. In both species, we found significant genetic differentiation between populations. However, the pattern of population differentiation (F(ST)) of R. apiculata differed from that of R. mucronata. Our results suggest that the two investigated species have different demographic history, even though they are sympatric and have similar reproductive systems.

Genetic variation in populations of Allothrombium pulvinum (Acari: Trombidiidae) from Northern Iran revealed by mitochondrial coxI and nuclear rDNA ITS2 sequences.

Allothrombium pulvinum Ewing is a common natural enemy of aphids and some other arthropods. So far, there are no studies that have addressed genetic variation of this predatory mite. We investigated genetic variation of A. pulvinum across its whole known range in Iran. A 410 bp portion of the mitochondrial cytochrome c oxidase subunit I gene (coxI) and 797-802 bp portion of the internal transcribed spacer 2 of rDNA (ITS2) were sequenced for 55 individuals from 11 populations, resulting in 12 and 26 haplotypes, respectively. In the coxI region, haplotype and nucleotide diversities varied among populations from 0.00 to 0.90 and from 0.0000 to 0.0110, respectively. In the ITS2 region they varied from 0.20 to 0.91 and from 0.0006 to 0.0023, respectively. For both gene regions the highest haplotype and nucleotide diversities were detected in population Mahmoud Abad from northern Iran. Statistically significant population differentiation (F (ST)) was detected in most pair-wise population comparisons. The results of population differentiation for both gene regions were generally congruent indicating that A. pulvinum from Iran consists of genetically different populations. This suggests that A. pulvinum comprises at least two geographically distinct populations or even more than one species. This study is an initial step towards understanding genetic variation of A. pulvinum, a taxon for which little molecular information is available. More intensive sampling and analysis of additional DNA regions are necessary for more detailed classification of this taxon.

Association between duplicated maltase genes and the transcriptional regulation for the carbohydrate changes in Drosophila melanogaster.

Gene duplication could promote phenotypic and genetic adaptation to various environments. To understand the effects of gene duplication on transcriptional regulation associated with environmental changes, we focused on the starch hydrolysis pathway, in which amylase enzymes together with maltase enzymes hydrolyze starch into glucose. Drosophila genomes involve ten duplicated Maltase genes. We examined the levels of transcription of the nine of these genes in 36 lines of Drosophila melanogaster collected from a natural population. In the investigated population, the levels of transcription were different between the two dietary carbohydrate sources, glucose and starch. At the transcriptional level, a single Maltase gene, which transcribes the specific Maltase transcripts, worked together with an Amylase gene in the pathway. The three of nine genes responded to carbohydrate changes, and the degree of the response was similar to Amylase gene. Our results suggest that gene duplication could increase capacity of the transcriptional regulation associated with environmental changes.

Convergent evolution of body color between sympatric freshwater fishes via different visual sensory evolution.

Although there are many examples of color evolution potentially driven by sensory drive, only few studies have examined whether distinct species inhabiting the same environments evolve similar body colors via shared sensory mechanisms. In this study, we tested whether two sympatric freshwater fish taxa, halfbeaks of the genus Nomorhamphus and ricefishes of the genus Oryzias in Sulawesi Island, converge in both body color and visual sensitivity. After reconstructing the phylogeny separately for Nomorhamphus and Oryzias using transcriptome-wide sequences, we demonstrated positive correlations of body redness between these two taxa across environments, even after phylogenetic corrections, which support convergent evolution. However, substantial differences were observed in the expression profiles of opsin genes in the eyes between Nomorhamphus and Oryzias. Particularly, the expression levels of the long wavelength-sensitive genes were negatively correlated between the taxa, indicating that they have different visual sensitivities despite living in similar light environments. Thus, the convergence of body colorations between these two freshwater fish taxa was not accompanied by convergence in opsin sensitivities. This system presents a case in which body color convergence can occur between sympatric species via different mechanisms.

Short 5'-flanking regions of the Amy gene of Drosophila kikkawai affect amylase gene expression and respond to food environments.

Evolution of the duplicated genes and regulation in gene expression is of great interest, especially in terms of adaptation. Molecular population genetic and evolutionary studies on the duplicated amylase genes of Drosophila species have suggested that their 5'-flanking (cis-regulatory) regions play an important role in evolution of these genes. For better understanding of evolution of the duplicated amylase genes and gene expression, we studied functional significance of the Amy1 gene of Drosophila kikkawai using in vitro deletion mutagenesis followed by P-element-mediated germline transformation. We found that a 1.6-kb of the 5'-flanking region can produce strikingly higher level of larval amylase activity on starch food compared with that on glucose food. We found two cis-regulatory elements, which increase larval amylase activity on starch food. We also found a larval cis-regulatory element, which responds to the food difference. This food-response element is necessary for the function of the element increasing larval activity on starch food. A 5-bp deletion in a putative GRE caused high amylase activity, indicating a cis-regulatory element decreasing amylase activity. These cis-regulatory elements identified in the 5'-flanking region could be the targets of natural selection.

Phylogeography of Eurasian Larix species inferred from nucleotide variation in two nuclear genes.

Larch (Larix Mill.) is one of the most widely distributed tree genera in Eurasia. To determine population structure and to verify classification of five species and three varieties of the Eurasian Larix species, we investigated levels and patterns of nucleotide variation of two nuclear gene regions: the 4-coumarate coenzyme A ligase (4CL) and the coumarate 3-hydroxylase (C3H). In the 4CL region nucleotide diversity at silent sites (pi(sil)) varied between 0.0020 in L. gmelinii to 0.0116 in L. gmelinii var. japonica and in the C3H region between 0.0019 in L. kaempferi to 0.0066 in L. gmelinii var. japonica. In both gene regions statistically significant population differentiation (F(ST)) was detected among adjacent refugial populations of some species suggesting limited gene flow and/or long time isolation of some refugial populations. On the other hand, populations of L. sukaczewii from northwestern Russia, which was glaciated 20,000 years ago showed no differentiation. This result is consistent with recent postglacial origin of these populations. Haplotype composition of some of the investigated Eurasian Larix species suggested that they are considerably diverged. Some haplotypes were unique to individual species. Our results indicate that more intensive sampling especially from known refugial regions is necessary for inferring correct classification of Eurasian Larix species and inferring their postglacial migration.

The P element invaded rapidly and caused hybrid dysgenesis in natural populations of Drosophila simulans in Japan.

Transposable elements not only can change genomic positions and disperse across the gene pool, but also can jump to another species through horizontal transmission. Of late, the P element, a DNA transposon in insects, was shown to cross the genetic boundary from Drosophila melanogaster to D. simulans in Europe around 2006. To understand the dynamics of transposable elements, especially in the early stages of invasion, we examined 63 lines of D. simulans from 11 natural populations in Japan established in 1976-2015. Based on PCR analyses, P elements were demonstrated to exist in Japan in 2008 and later. One copy of the full-length P element was identified and mapped to a site on chromosome 3 L in a genome. All of 18 copies of P elements examined shared "A" at the nucleotide position 2040, which is representative of the direct descendants of the original P element that invaded in D. simulans. We also found that some lines having P elements can induce intensive gonadal dysgenesis in D. simulans at 29°C. Our present results imply that P elements in D. simulans arrived at the east end of Asia just a few years later than or almost simultaneously to the initial invasion in Europe, Africa, and North America, suggesting a more astonishingly rapid spread than previously assumed.

Comprehensive molecular phylogeny of the sub-family Dipterocarpoideae (Dipterocarpaceae) based on chloroplast DNA sequences.

Dipterocarpoideae, the largest sub-family of well-known plant family Dipterocarpaceae, dominates in South Asian rain forests. Although several previous studies addressed the phylogeny of the Dipterocarpaceae family, relationships among many of its genera from the Dipterocarpoideae sub-family are still not well understood. In particular, little is known about the relationships of the genera Vateriopsis, Stemonoporus, Vateria and inconsistence remains between phylogenetic results and taxonomic classifications of Shorea and Hopea species. We studied molecular phylogeny of the sub-family Dipterocarpoideae using the trnL-trnF spacer, trnL intron and the matK gene sequences of chloroplast DNA (cpDNA). This study is the first comprehensive phylogeny reconstruction for the sub-family Dipterocarpoideae based on cpDNA, as it includes most genera (14) and a large number of species (79) with most species endemic to Sri Lanka, as well as one species from Seychelles and one species from the genus Monotes from Madagascar. Phylogenetic trees were constructed using the Neighbor Joining (NJ) and Maximum Likelihood (ML) methods using combined set of sequences including all three cpDNA regions. The topologies of the NJ and ML trees were to a certain extent, consistent with the current taxonomy of Dipterocarpoideae based on morphology and with previous molecular phylogenies based on cpDNA. Furthermore, our results provided new evidence regarding the relationships of the following genera: Vateriopsis and Stemonoporus and about the validity of the previous morphology based classifications of Shorea species. In addition, the topology of our trees was consistent with the classification of Shorea species proposed by Maury (1978), Maury-Lechon (1979) and Symington (1943). Finally, our results provided evidence for the affinity of the genus Monotes to Asian Dipterocarpoideae rather than to Tiliaceae and indicated that it is a good candidate for outgroup species for future studies of the former sub-family.

A single-amino-acid change of the gustatory receptor gene, Gr5a, has a major effect on trehalose sensitivity in a natural population of Drosophila melanogaster.

Variation in trehalose sensitivity and nucleotide sequence polymorphism of the Gr5a gene encoding the gustatory receptor to sugar trehalose were investigated in 152 male lines of Drosophila melanogaster collected from a natural population. Among the observed 59 segregating sites, some pairs of sites showed significant linkage disequilibrium. A single SNP, which results in the Ala218Thr amino acid change, was significantly associated with trehalose sensitivity, as previously suggested. Threonine at amino acid position 218 was found to be the ancestral form in D. melanogaster, suggesting that low trehalose sensitivity was an ancestral form with respect to the receptor function. There was large genetic variation in trehalose sensitivity. It was continuously distributed, indicating that trehalose sensitivity measured by the behavioral assay is a quantitative trait. These results suggest that apart from the Gr5a gene, other genetic factors contribute to variation in trehalose sensitivity. Nucleotide diversity (pi) and nucleotide variation (theta) per site were 0.00874 and 0.00590, respectively. Fu and Li's test and the MK test showed no significant departure from the expectation of selective neutrality in the Gr5a gene. However, we rejected selective neutrality by Tajima's test and Fay and Wu's test with the observed level of recombination. We discuss possible causes of the observed pattern of nucleotide variation in the gustatory receptor Gr5a gene.

Codon bias differentiates between the duplicated amylase loci following gene duplication in Drosophila.

We examined the pattern of synonymous substitutions in the duplicated Amylase (Amy) genes (called the Amy1- and Amy3-type genes, respectively) in the Drosophila montium species subgroup. The GC content at the third synonymous codon sites of the Amy1-type genes was higher than that of the Amy3-type genes, while the GC content in the 5'-flanking region was the same in both genes. This suggests that the difference in the GC content at third synonymous sites between the duplicated genes is not due to the temporal or regional changes in mutation bias. We inferred the direction of synonymous substitutions along branches of a phylogeny. In most lineages, there were more synonymous substitutions from G/C (G or C) to A/T (A or T) than from A/T to G/C. However, in one lineage leading to the Amy1-type genes, which is immediately after gene duplication but before speciation of the montium species, synonymous substitutions from A/T to G/C were predominant. According to a simple model of synonymous DNA evolution in which major codons are selectively advantageous within each codon family, we estimated the selection intensity for specific lineages in a phylogeny on the basis of inferred patterns of synonymous substitutions. Our result suggested that the difference in GC content at synonymous sites between the two Amy-type genes was due to the change of selection intensity immediately after gene duplication but before speciation of the montium species.

Nucleotide variation of seven genes in Drosophila kikkawai.

We examined levels and patterns of nucleotide variation in 21 strains of Drosophila kikkawai from Miyako island, Japan for the partial regions of the following seven nuclear genes: Adh, Ddc, esc, ksr, Pgi, su(f), and Tpi. The nucleotide variation at total sites (pi(t)) ranged from 0.0013 in the ksr, to 0.0173 in the Adh. The nucleotide divergence at total sites (K(t)) between D. kikkawai and D. lini ranged from 0.0286 in the Tpi to 0.0687 in the su(f). The levels of nucleotide polymorphism and divergence were heterogeneous among the investigated gene regions. The HKA test, which tests imbalance between the intra and interspecific nucleotide variation, showed that the intraspecific nucleotide variation in the Pgi region was much lower than the interspecific variation, while intraspecific variation in the Tpi region was only slightly lower than interspecific variation. The MK test showed an excess of low frequency replacement polymorphic changes in the Adh region, suggesting that most replacement mutations are deleterious. Fay and Wu's test detected an excess of newly arisen variants in the Ddc region. In total, four of the seven gene regions showed significant deviation from the neutrality.

High prevalence of non-synonymous substitutions in mtDNA of cichlid fishes from Lake Victoria.

When a population size is reduced, genetic drift may fix slightly deleterious mutations, and an increase in nonsynonymous substitution is expected. It has been suggested that past aridity has seriously affected and decreased the populations of cichlid fishes in Lake Victoria, while geographical studies have shown that the water levels in Lake Tanganyika and Lake Malawi have remained fairly constant. The comparably stable environments in the latter two lakes might have kept the populations of cichlid fishes large enough to remove slightly deleterious mutations. The difference in the stability of cichlid fish population sizes between Lake Victoria and the Lakes Tanganyika and Malawi is expected to have caused differences in the nonsynonymous/synonymous ratio, ω (=dN/dS), of the evolutionary rate. Here, we estimated ω and compared it between the cichlids of the three lakes for 13 mitochondrial protein-coding genes using maximum likelihood methods. We found that the lineages of the cichlids in Lake Victoria had a significantly higher ω for several mitochondrial loci. Moreover, positive selection was indicated for several codons in the mtDNA of the Lake Victoria cichlid lineage. Our results indicate that both adaptive and slightly deleterious molecular evolution has taken place in the Lake Victoria cichlids' mtDNA genes, whose nonsynonymous sites are generally conserved.