Inferring the demographic history of Japanese cedar, Cryptomeria japonica, using amplicon sequencing.

The evolution of a species depends on multiple forces, such as demography and natural selection. To understand the trajectory and driving forces of evolution of a target species, it is first necessary to uncover that species' population history, such as past and present population sizes, subdivision and gene flow, by using appropriate genetic markers. Cryptomeria japonica is a long-lived monoecious conifer species that is distributed in Japan. There are two main lines (omote-sugi and ura-sugi), which are distinguished by apparent differences in morphological traits that may have contributed to their local adaptation. The evolution of these morphological traits seems to be related to past climatic changes in East Asia, but no precise estimate is available for the divergence time of these two lines and the subsequent population dynamics in this species. Here, we analyzed the nucleotide variations at 120 nuclear genes in 94 individuals by using amplicon sequencing in combination with high-throughput sequencing technologies. Our analysis indicated that the population on Yakushima Island, the southern distribution limit of C. japonica in Japan, diverged from the other populations 0.85 million years ago (MYA). The divergence time of the other populations on mainland Japan was estimated to be 0.32 MYA suggesting that the divergence of omote-sugi and ura-sugi might have occurred before the last glacial maximum. Although we found modest levels of gene flow between the present populations, the long-term isolation and environmental heterogeneity caused by climatic changes might have contributed to the differentiation of the lines and their local adaptation.

Inferences of population structure and demographic history for Taxodium distichum, a coniferous tree in North America, based on amplicon sequencing analysis.

PREMISE OF THE STUDY: Studies of natural genetic variation can elucidate the genetic basis of phenotypic variation and the past population structure of species. Our study species, Taxodium distichum, is a unique conifer that inhabits the flood plains and swamps of North America. Morphological and ecological differences in two varieties, T. distichum var. distichum (bald cypress) and T. distichum var. imbricarium (pond cypress), are well known, but little is known about the level of genetic differentiation between the varieties and the demographic history of local populations. METHODS: We analyzed nucleotide polymorphisms at 47 nuclear loci from 96 individuals collected from the Mississippi River Alluvial Valley (MRAV), and Gulf Coastal populations in Texas, Louisiana, and Florida using high-throughput DNA sequencing. Standard population genetic statistics were calculated, and demographic parameters were estimated using a composite-likelihood approach. KEY RESULTS: Taxodium distichum in North America can be divided into at least three genetic groups, bald cypress in the MRAV and Texas, bald cypress in Florida, and pond cypress in Florida. The levels of genetic differentiation among the groups were low but significant. Several loci showed the signatures of positive selection, which might be responsible for local adaptation or varietal differentiation. CONCLUSIONS: Bald cypress was genetically differentiated into two geographical groups, and the boundary was located between the MRAV and Florida. This differentiation could be explained by population expansion from east to west. Despite the overlap of the two varieties' ranges, they were genetically differentiated in Florida. The estimated demographic parameters suggested that pond cypress split from bald cypress during the late Miocene.

Difference in flowering time can initiate speciation of nocturnally flowering species.

Isolation mechanisms that prevent gene flow between populations prezygotically play important roles in achieving speciation. In flowering plants, the nighttime flowering system provides a mechanism for isolation from diurnally flowering species. Although this system has long been of interest in evolutionary biology, the evolutionary process leading to this system has yet to be elucidated because of the lack of good model species. However, the genetic mechanisms underlying the differences in flowering times and the traits that attract pollinators between a pair of diurnally and nocturnally flowering species have recently been identified in a few cases. This identification enables us to build a realistic model for theoretically studying the evolution of a nocturnally flowering species. In this study, based on previous experimental data, we assumed a model in which two loci control the flowering time and one locus determines a trait that attracts pollinators. Using this model, we evaluated the possibility of the evolution of a nocturnally flowering species from a diurnally flowering ancestor through simulations. We found that a newly emerging nighttime flowering flower exhibited a sufficiently high fitness, and the evolution of a nocturnally flowering species from a diurnally flowering species could be achieved when hybrid viability was intermediate to low, even in a completely sympatric situation. Our results suggest that the difference in flowering time can act as a magic trait that induces both natural selection and assortative mating and would play an important role in speciation between diurnally and nocturnally flowering species pairs.

Speciation through sensory drive in cichlid fish.

Theoretically, divergent selection on sensory systems can cause speciation through sensory drive. However, empirical evidence is rare and incomplete. Here we demonstrate sensory drive speciation within island populations of cichlid fish. We identify the ecological and molecular basis of divergent evolution in the cichlid visual system, demonstrate associated divergence in male colouration and female preferences, and show subsequent differentiation at neutral loci, indicating reproductive isolation. Evidence is replicated in several pairs of sympatric populations and species. Variation in the slope of the environmental gradients explains variation in the progress towards speciation: speciation occurs on all but the steepest gradients. This is the most complete demonstration so far of speciation through sensory drive without geographical isolation. Our results also provide a mechanistic explanation for the collapse of cichlid fish species diversity during the anthropogenic eutrophication of Lake Victoria.

Correlation between nuptial colors and visual sensitivities tuned by opsins leads to species richness in sympatric Lake Victoria cichlid fishes.

Reproductive isolation that prevents interspecific hybridization between closely related coexisting species maintains sympatric species diversity. One of the reproductive isolations is mate choice based on color signals (breeding color perceived by color vision). This is well known in several animal taxa, yet little is known about its genetic and molecular mechanism. Lake Victoria cichlid fishes are thought to be an example of sympatric species diversity. In the species inhabiting different light environments in rocky shore, speciation by sensory drive through color signals has been proposed by analyses of the long wavelength-sensitive (LWS) opsin gene and the male nuptial coloration. However, the genetic and molecular mechanism of how diversity of sympatric species occurring in the same habitat is maintained remains unknown. To address this issue, we determined nucleotide sequences of eight opsins of six sympatric species collected from a sandy-muddy shore--an ideal model system for studying sympatric species. Among eight opsins, the LWS and RH1 alleles were diversified and one particular allele is dominant or fixed in each species, and we propose that this is due to natural selection. The functions of their LWS alleles were also diversified as shown by absorption measurements of reconstituted visual pigments. To analyze the relationship between nuptial coloration and the absorption of LWS pigments, we systematically evaluated and defined nuptial coloration. We showed that the coloration was species specific with respect to hue and significantly differentiated by the index values of hue (dominant wavelength: λ(d)). The λ(d) value of the male nuptial coloration correlated with the absorption of LWS pigments from all the species, suggesting that reproductive isolation through mate choice using color signals may prevent sympatric interspecific hybridization, thereby maintaining the species diversity in sympatric species in Lake Victoria.

Symbiosis between hydra and chlorella: molecular phylogenetic analysis and experimental study provide insight into its origin and evolution.

Although many physiological studies have been reported on the symbiosis between hydra and green algae, very little information from a molecular phylogenetic aspect of symbiosis is available. In order to understand the origin and evolution of symbiosis between the two organisms, we compared the phylogenetic relationships among symbiotic green algae with the phylogenetic relationships among host hydra strains. To do so, we reconstructed molecular phylogenetic trees of several strains of symbiotic chlorella harbored in the endodermal epithelial cells of viridissima group hydra strains and investigated their congruence with the molecular phylogenetic trees of the host hydra strains. To examine the species specificity between the host and the symbiont with respect to the genetic distance, we also tried to introduce chlorella strains into two aposymbiotic strains of viridissima group hydra in which symbiotic chlorella had been eliminated in advance. We discussed the origin and history of symbiosis between hydra and green algae based on the analysis.

Difference in flowering time as an isolating barrier.

Although many theoretical studies have reported strong effects of different flowering times on reproductive isolation, such studies have all focused on the different flowering time within a season, and the subsequently developed models are difficult to apply to the cases of diurnal- and nocturnal-flowering species pairs. The different flowering times within a day differ from those within a season because of the simultaneous opening and closing of the flowers for each species and the carry-over of the pollen from early to later times. In this study, we consider pollinator-mediated, diurnal- and nocturnal-flowering plants and build a new model to study the effects of the different flowering times within a day on reproductive isolation. We assume two loci, each with two alleles, which determine the opening and closing times of flowers, respectively. We numerically calculate the changes in the frequencies of the gametes in a model incorporating the reductions in hybrid viability, flowering costs, recombination rate and degree of dominance at each locus. We found that the early-opening flowers had a much higher fitness than the late-opening flowers and that the maintenance of the two species was difficult even if their flowering times were not overlapping. Therefore, some other mechanisms, such as pollinator preference, may be required to explain the coexistence of closely related diurnal and nocturnal flowers.

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.

Evolution of the number of LRRs in plant disease resistance genes.

The largest group of plant resistance (R) genes contain the regions that encode the nucleotide-binding site (NBS) and leucine-rich repeat (LRR) domains (NBS-LRR genes). To gain new resistance, amino acid substitutions and changes in number of the LRRs that recognize the presence of pathogens are considered important. In this study, we focus on the evolution of the number of LRRs and analyze the genome data of five plant species, Arabidopsis thaliana, Oryza sativa, Medicago truncatula, Lotus japonicus and Populus trichocarpa. We first categorized the NBS-LRR genes in each species into groups and subgroups based on the phylogenetic relationships of their NBS domain sequences. Then we estimated the evolutionary rate of the number of LRRs relative to the synonymous divergence in the NBS domain sequences by a maximum likelihood method assuming the single stepwise mutation model. The estimates ranged from 4.5 to 600 and differed between groups in the same species or between species. This indicated different roles played by different groups of the NBS-LRR genes within a species or the effects of various life history characteristics, such as generation time, of the species. We also tested the fit of the model to the data using the variance of number of LRRs in each subgroup. In some subgroups in some plants (16 out of 174 subgroups), the results of simulation using the estimated rates significantly deviated from the observed data. Those subgroups may have undergone different modes of selection from the other subgroups.

Evaluating the genetic diversity in two tropical leguminous trees, Dalbergia cochinchinensis and D. nigrescens, in lowland forests in Cambodia and Thailand using MIG-seq.

It is vital to measure the levels of genetic diversity and differentiation between populations in a species to understand the current genetic structure and evolution of the species. Here, MIG-seq (multiplexed inter-simple sequence repeat genotyping by sequencing) was employed to assess the genetic variation in two tropical leguminous tree species, Dalbergia cochinchinensis and D. nigrescens, in Cambodia and Thailand. Sequence data for 255-618 loci, each with an approximate length of 100 bp, were obtained, and the nucleotide diversity, Tajima's D and FST were computed. The estimates calculated from the data obtained by MIG-seq were compared to those obtained by Sanger sequencing of nine nuclear coding genes in D. cochinchinensis in our previous study. The nucleotide diversity at the MIG-seq loci was slightly higher than that at silent sites in the coding loci, whereas the FST values at the MIG-seq loci were generally lower than those at the coding loci, although the differences were not significant. Moreover, nucleotide diversities within populations of the two species were similar to each other, at approximately 0.005. Three and four population clusters were genetically recognized in D. cochinchinensis and D. nigrescens, respectively. Although the populations were differentiated from each other, the levels of differentiation among them, as measured by FST, were higher in D. cochinchinensis than in D. nigrescens. This indicates higher levels of gene flow between the populations in the latter species. We recommend using MIG-seq for quick surveys of genetic variation because it is cost-effective and results in smaller variance in the estimates of population genetic parameters.

Genetic structure and putative selective sweep in the pioneer tree, Zanthoxylum ailanthoides.

Zanthoxylum ailanthoides Siebold & Zucc. is one of the most frequently encountered pioneer trees in Japanese warm-temperate evergreen oak forests. Our previous study in one region of Japan suggested high levels of population differentiation and putative natural selection acting on one of the nuclear loci analyzed. Here, we extend our analysis to study the genetic structure of 10 populations of Z. ailanthoides across Japan using 9 simple sequence repeat (SSR) loci for a better understanding of its genetic structure. First, the southernmost population (Kagoshima) in the samples was found to have the highest genetic diversity, suggesting there was a glacial refugium at or near the location of the population. Second, relatively strong genetic differentiation was found among populations, and there was a positive correlation between genetic distances and geographic distances (Mantel test; P < 0.001). Based on this information, we analyzed nucleotide variation at the putatively selected locus homologous to the gene encoding the ADP-glucose pyrophosphorylase large subunit (agpL). Despite the strong genetic differentiation among populations suggested by the SSR loci, the agpL locus was monomorphic in almost all populations analyzed. The results of this study strongly supported the possibility of a selective sweep at or near the agpL locus.

Divergent selection on opsins drives incipient speciation in Lake Victoria cichlids.

Divergent natural selection acting on ecological traits, which also affect mate choice, is a key element of ecological speciation theory, but has not previously been demonstrated at the molecular gene level to our knowledge. Here we demonstrate parallel evolution in two cichlid genera under strong divergent selection in a gene that affects both. Strong divergent natural selection fixed opsin proteins with different predicted light absorbance properties at opposite ends of an environmental gradient. By expressing them and measuring absorbance, we show that the reciprocal fixation adapts populations to divergent light environments. The divergent evolution of the visual system coincides with divergence in male breeding coloration, consistent with incipient ecological by-product speciation.

Effects of cyclic changes in population size on neutral genetic diversity.

Recurrent changes in population size are often observed in nature, influencing the efficiency of selection and consequently affecting organismal evolution. Thus, it is important to know whether such changes occurred in the past history of a focal population of evolutionary interests. Here, we focused on cyclic changes in population size and investigated the distributional properties of Tajima's D and its power to distinguish a cyclic change model compared with the standard neutral model, changing the frequency and magnitude of the cyclic change. With very low or very high frequencies of the cycle, the distribution of Tajima's D was similar to that in a constant size population, as demonstrated by previous theoretical works. Otherwise, its mean was negative or positive, and its variance was smaller or larger depending on the time of sampling. The detection rate of the cyclic change against the constancy in size by Tajima's D depended on the sample size, the number of loci, and the time of sampling in addition to the frequency and amplitude of the cycle. Using sequence data of several tens of loci, the detection rate was fairly high if the frequency was intermediate and the sampling was made when population size was large; otherwise, the detection rate was not high. We also found that cyclic change could be discriminated from simple expansion or shrinkage of a population by Tajima's D only if the frequency was in a limited range and the sampling was made when the population was large.

Horizontal Transmission of Symbiotic Green Algae Between Hydra Strains.

Some hydra strains belonging to the vulgaris group show a symbiotic relationship with green algae Chlorococcum sp. The symbiotic green algae can escape from the host polyps and can form swimming zoospores (which have two flagella) in culture solution. We observed that co-culture with the symbiotic polyps caused horizontal transmission of the symbionts into some non-symbiotic hydra strains that have no symbionts in nature and that belong not only to the vulgaris group but also to other hydra species groups. Although most of the horizontal transmission has ended in transient symbioses, a newly formed symbiosis between the symbiotic Chlorococcum sp. and strain 105 of Hydra vulgaris (Hydra magnipapillata) has been sustained for more than five years and has caused morphological and behavioral changes in the host polyps. We named this strain 105G. The asexual proliferation rate by budding increased under light conditions, although the feeding activity decreased and the polyp size was reduced in strain 105G. This new symbiosis between Chlorococcum sp. and strain 105G of H. vulgaris provides us with an intriguing research system for investigating the origin of symbiosis.

Understanding and monitoring the consequences of human impacts on intraspecific variation.

Intraspecific variation is a major component of biodiversity, yet it has received relatively little attention from governmental and nongovernmental organizations, especially with regard to conservation plans and the management of wild species. This omission is ill-advised because phenotypic and genetic variations within and among populations can have dramatic effects on ecological and evolutionary processes, including responses to environmental change, the maintenance of species diversity, and ecological stability and resilience. At the same time, environmental changes associated with many human activities, such as land use and climate change, have dramatic and often negative impacts on intraspecific variation. We argue for the need for local, regional, and global programs to monitor intraspecific genetic variation. We suggest that such monitoring should include two main strategies: (i) intensive monitoring of multiple types of genetic variation in selected species and (ii) broad-brush modeling for representative species for predicting changes in variation as a function of changes in population size and range extent. Overall, we call for collaborative efforts to initiate the urgently needed monitoring of intraspecific variation.

Gene genealogy and properties of test statistics of neutrality under population growth.

We consider the Wright-Fisher model with exponential population growth and investigate effects of population growth on the shape of genealogy and the distributions of several test statistics of neutrality. In the limiting case as the population grows rapidly, the rapid-growth-limit genealogy is characterized. We obtained approximate expressions for expectations and variances of test statistics in the rapid-growth-limit genealogy and star genealogy. The distributions in the star genealogy are narrower than those in the cases of the simulated and rapid-growth-limit genealogies. The expectations and variances of the test statistics are monotone decreasing functions of the time length of the expansion, and the higher power of R(2) against population growth is suggested to be due to their smaller variances rather than to change of the expectations. We also investigated by simulation how quickly the distributions of test statistics approach those of the rapid-growth-limit genealogy.

Contrasting patterns of DNA variation in natural populations of two related conifers, Cryptomeria japonica and Taxodium distichum (Cupressaceae sensu lato).

We investigated DNA variation within and between two closely related conifers, Cryptomeria japonica and Taxodium distichum, at nuclear loci encoding ferredoxin, glutamyl-tRNA reductase, lycopene beta cyclase, and phosphoribosylanthranilate transferase. Average nucleotide diversity at silent sites was estimated to be 0.0035 (SE 0.0012) in C. japonica and 0.0058 (SE 0.0006) in T. distichum. One population in C. japonica was differentiated from the others but generally there was not much differentiation among populations or varieties within the two species. However, the two species seemed to differ in frequency spectra of DNA polymorphisms. Excesses of intermediate-frequency variants were found in C. japonica, whereas excesses of both rare and high-frequency variants were found in T. distichum, which suggested different histories of population structures in the two species. Deviations from the standard neutral expectations in DNA polymorphisms were found by applications of neutrality tests. The results show that actions of selection to respective loci seem to differ between the two species, indicating differences of interaction among evolutionary factors.

Comparative phylogeography of two bitterlings, Tanakia lanceolata and T. limbata (Teleostei, Cyprinidae), in Kyushu and adjacent districts of western Japan, based on mitochondrial DNA analysis.

To understand the geographical patterns of genetic variation in freshwater fishes in western Japan, the genetic structures of populations of Tanakia lanceolata and T. limbata (Teleostei, Cyprinidae) in this area were investigated using mitochondrial DNA (mtDNA) cytb sequences. Neighbor-joining trees of mtDNA haplotypes revealed four and three genetically divergent groups in T. lanceolata and T. limbata, respectively. Each group was restricted to one or the other of the geographical regions in the area studied. The patterns of geographical divergence in the two species showed some similarities, which seem to reflect common historical events experienced by freshwater fishes distributed in western Japan. On the other hand, dissimilarities were also found in the patterns, indicating that species-specific historical processes also occurred. Within one region, T. lanceolata was less differentiated than T. limbata, suggesting a difference in the dispersal abilities of the two bitterlings. In addition, several individuals in Kunichika River on Shikoku island were morphologically identified as T. lanceolata but had mtDNA haplotypes of T. limbata. We tentatively suggest that these individuals are hybrids of the two species, but further studies employing nuclear markers are necessary to validate this hypothesis.

A neutral model with fluctuating population size and its effective size.

We consider a diffusion model with neutral alleles whose population size is fluctuating randomly. For this model, the effects of fluctuation of population size on the effective size are investigated. The effective size defined by the equilibrium average heterozygosity is larger than the harmonic mean of population size but smaller than the arithmetic mean of population size. To see explicitly the effects of fluctuation of population size on the effective size, we investigate a special case where population size fluctuates between two distinct states. In some cases, the effective size is very different from the harmonic mean. For this concrete model, we also obtain the stationary distribution of the average heterozygosity. Asymptotic behavior of the effective size is obtained when the population size is large and/or autocorrelation of the fluctuation is weak or strong.

DNA variation in a conifer, Cryptomeria japonica (Cupressaceae sensu lato).

We investigated the nucleotide variation of a conifer, Cryptomeria japonica, and the divergence between this species and its closest relative, Taxodium distichum, at seven nuclear loci (Acl5, Chi1, Ferr, GapC, HemA, Lcyb, and Pat). Samples of C. japonica were collected from three areas, Kantou-Toukai, Hokuriku, and Iwate. No apparent geographic differentiation was found among these samples. However, the frequency spectrum of the nucleotide polymorphism revealed excesses of intermediate-frequency variants, which suggests that the population was not panmictic and a constant size in the past. The average nucleotide diversity, pi, for silent sites was 0.00383. However, values of pi for silent sites vary among loci. Comparisons of polymorphism to divergence among loci (the HKA test) showed that the polymorphism at the Acl5 locus was significantly lower. We also observed a nearly significant excess of replacement polymorphisms at the Lcyb locus. These results suggested possibilities of natural selection acting at some of the loci. Intragenic recombination was detected only once at the Chi1 locus and was not detected at the other loci. The low level of population recombination rate, 4Nr, seemed to be due to both low level of recombination, r, and small population size, N.