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.

Ghost introgression in ricefishes of the genus Adrianichthys in an ancient Wallacean lake.

Because speciation might have been promoted by ancient introgression from an extinct lineage, it is important to detect the existence of 'ghost introgression' in focal taxa and examine its contribution to their diversification. In this study, we examined possible ghost introgression and its contributions to the diversification of ricefishes of the genus Adrianichthys in Lake Poso, an ancient lake on Sulawesi Island, in which some extinctions are known to have occurred. Population-genomic analysis revealed that two extant Adrianichthys species, A. oophorus and A. poptae are reproductively isolated from each other. Comparisons of demographic models demonstrated that introgression from a ghost population, which diverged from the common ancestor of A. oophorus and A. poptae, is essential for reconstructing the demographic history of Adrianichthys. The best model estimated that the divergence of the ghost population greatly predated the divergence between A. oophorus and A. poptae, and that the ghost population secondarily contacted the two extant species within Lake Poso more recently. Genome scans and simulations detected a greatly divergent locus, which cannot be explained without ghost introgression. This locus was also completely segregated between A. oophorus and A. poptae. These findings suggest that variants that came from a ghost population have contributed to the divergence between A. oophorus and A. poptae, but the large time-lag between their divergence and ghost introgression indicates that the contribution of introgression may be restricted.

Seasonal gene expression signatures of delayed fertilization in Fagaceae.

In the family Fagaceae, fertilization is delayed by several weeks to 1 year after pollination, leading to 1- or 2-year fruiting species depending on whether fruiting occurs in the same or the next year after flowering. To investigate physiological responses underlying the regulation of delayed fertilization, we monitored seasonal changes in genome-wide gene expression in tissues including leaves and buds over 2 years under natural conditions in one- (Quercus glauca) and 2-year fruiting species (Lithocarpus edulis). Genes associated with metabolic changes in response to winter cold, photosynthesis and cell proliferation, which are essential for survival and growth, showed highly conserved seasonal expression profiles between species. However, seasonal expression profiles diverged between species in genes associated with pollination, an important process contributing to the origin and maintenance of the reproductive barrier between plant species. By comparing seasonal progression of ovule development and gene expression in pistillate flowers, we revealed that ovules started developing after winter in the 2-year fruiting species, which could be linked to the activation of genes involved in fertilization and female gametophyte development after winter. These findings suggest that the 2-year fruiting species may have evolved a requirement of winter cold to prevent fertilization before winter and facilitate fertilization and embryo development in the following spring when temperature rises. This study offers new possibilities to explore the evolution of reproductive strategies in Fagaceae.

Multiple colonizations and hybridization of a freshwater fish group on a satellite island of Sulawesi.

Repeated colonizations and resultant hybridization may increase lineage diversity on an island if introgression occurs only in a portion of the indigenous island lineage. Therefore, to precisely understand how island biodiversity was shaped, it is essential to reconstruct the history of secondary colonization and resultant hybridization both in time and space. In this study, we reconstructed the history of multiple colonizations of the Oryzias woworae species group, a freshwater fish group of the family Adrianichthyidae, from Sulawesi Island to its southeast satellite island, Muna Island. Phylogenetic and species tree analyses using genome-wide single-nucleotide polymorphisms revealed that all local populations on Muna Island were monophyletic, but that there were several genetically distinct lineages within the island. Population structure and phylogenetic network analyses demonstrated that colonization of this island occurred more than once, and that secondary colonization and resultant introgressive hybridization occurred only in one local population on the island. The spatially heterogeneous introgression induced by the multiple colonizations were also supported by differential admixture analyses. In addition, the differential admixture analyses detected reverse colonization from Muna Island to the Sulawesi mainland. Coalescence-based demographic inference estimated that these mutual colonizations occurred during the middle to late Quaternary period, during which sea level repeatedly declined; this indicates that the colonizations occurred via land bridges. We conclude that these mutual colonizations between Muna Island and the Sulawesi mainland, and the resultant spatially heterogeneous introgression shaped the current biodiversity of this species group in this area.

Symbiont specificity differs among green hydra strains.

The symbiotic hydra Hydra viridissima has a stable symbiotic relationship with the green alga Chlorella. This hydra appears to cospeciate with the symbiotic alga, and some strains are known to have strain-specific host/symbiont combinations. To investigate the mechanism of the specificity between host and symbiont, we explored the effect of the removal or exchange of symbionts in two distantly related H. viridissima strains (K10 and M9). In the K10 strain, severe morphological and behavioural changes were found in symbiont-removed and symbiont-exchanged polyps. Interestingly, both polyps showed a similar gene expression pattern. The gene ontology (GO) enrichment analysis revealed that the removal or exchange of symbionts caused the downregulation of genes involved in the electron transport chain and the upregulation of genes involved in translation in the K10 strain. On the other hand, symbiont-removed and symbiont-exchanged M9 polyps showed modest changes in their morphology and behaviour compared with the K10 strain. Furthermore, the patterns of the gene expression changes in the M9 strain were quite different between the symbiont-removed and symbiont-exchanged polyps. Our results suggested that the regulation of energy balance is one of the crucial mechanisms for maintaining symbiotic relationships in green hydra, and this mechanism differs between the strains.

Deeply divergent freshwater fish species within a single river system in central Sulawesi.

Sulawesi is a biodiversity hotspot for ricefishes (Adrianichthyidae), with many species endemic to the central part of this island in single ancient lakes or lake systems. Frequent vicariance by lake fragmentation since the Pliocene may be largely responsible for diversification in this family. In this study, we demonstrate that not only lacustrine species but also riverine species in this area are also deeply divergent even within a single river system. A mitochondrial phylogeny revealed that a ricefish population newly discovered from Cerekang River is sister to Oryzias dopingdopingensis Mandagi, Mokodongan, Tanaka, & Yamahira, another riverine species endemic to Doping-doping River. However, the Cerekang Oryzias was genetically isolated from O. dopingdopingensis, despite that Cerekang River and Doping-doping River share a connection across estuarine waters. This separation was supported by phylogenomic trees and population genetic structure analyses based on genome-wide single nucleotide polymorphisms. Coalescent-based demographic inference demonstrated that the ancestral population of these two riverine ricefishes had experienced a substantial population decrease and subsequently diverged into two sub-populations. Because the Cerekang Oryzias was also morphologically distinguished from O. dopingdopingensis, we described it as a new species, O. landangiensis. We infer that O. landangiensis and O. dopingdopingensis are of lake-origin and are relic species which were left in these rivers after the lake disappeared, and that they have lost their dispersal ability when inhabiting the ancient lake. The lost dispersal ability possibly contributed to the formation of the biodiversity hotspot for this fish group on this island.

Pollinator sharing, copollination, and speciation by host shifting among six closely related dioecious fig species.

The obligate pollination mutualism between figs (Ficus, Moraceae) and pollinator wasps (Agaonidae, Hymenoptera) is a classic example of cospeciation. However, examples of phylogenetic incongruencies between figs and their pollinators suggest that pollinators may speciate by host shifting. To investigate the mechanism of speciation by host shifting, we examined the phylogenetic relationships and population genetic structures of six closely related fig species and their pollinators from southern China and Taiwan-Ryukyu islands using various molecular markers. The results revealed 1) an extraordinary case of pollinator sharing, in which five distinct fig species share a single pollinator species in southern China; 2) two types of copollination, namely, sympatric copollination by pollinator duplication or pollinator migration, and allopatric copollination by host migration and new pollinator acquisition; 3) fig species from southern China have colonized Taiwan repeatedly and one of these events has been followed by host shifting, reestablishment of host specificity, and pollinator speciation, in order. Based on our results, we propose a model for pollinator speciation by host shifting in which the reestablishment of host-specificity plays a central role in the speciation process. These findings provide important insights into understanding the mechanisms underlying pollinator speciation and host specificity in obligate pollination mutualism.

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.

Copy number analyses of DNA repair genes reveal the role of poly(ADP-ribose) polymerase (PARP) in tree longevity.

Long-lived organisms are exposed to the risk of accumulating mutations due to DNA damage. Previous studies in animals have revealed the positive relationship between the copy number of DNA repair genes and longevity. However, the role of DNA repair in the lifespan of plants remains poorly understood. Using the recent accumulation of the complete genome sequences of diverse plant species, we performed systematic comparative analyses of the copy number variations of DNA repair genes in 61 plant species with different lifespans. Among 121 DNA repair gene families, PARP gene family was identified as a unique gene that exhibits significant expansion in trees compared to annual and perennial herbs. Among three paralogs of plant PARPs, PARP 1 showed a close association with growth rate. PARPs catalyze poly(ADP-ribosyl)ation and play pivotal roles in DNA repair and antipathogen defense. Our study suggests the conserved role of PARPs in longevity between plants and animals.

Mitochondrial introgression by ancient admixture between two distant lacustrine fishes in Sulawesi Island.

Sulawesi, an island located in a biogeographical transition zone between Indomalaya and Australasia, is famous for its high levels of endemism. Ricefishes (family Adrianichthyidae) are an example of taxa that have uniquely diversified on this island. It was demonstrated that habitat fragmentation due to the Pliocene juxtaposition among tectonic subdivisions of this island was the primary factor that promoted their divergence; however, it is also equally probable that habitat fusions and resultant admixtures between phylogenetically distant species may have frequently occurred. Previous studies revealed that some individuals of Oryzias sarasinorum endemic to a tectonic lake in central Sulawesi have mitochondrial haplotypes that are similar to the haplotypes of O. eversi, which is a phylogenetically related but geologically distant (ca. 190 km apart) adrianichthyid endemic to a small fountain. In this study, we tested if this reflects ancient admixture of O. eversi and O. sarasinorum. Population genomic analyses of genome-wide single-nucleotide polymorphisms revealed that O. eversi and O. sarasinorum are substantially reproductively isolated from each other. Comparison of demographic models revealed that the models assuming ancient admixture from O. eversi to O. sarasinorum was more supported than the models assuming no admixture; this supported the idea that the O. eversi-like mitochondrial haplotype in O. sarasinorum was introgressed from O. eversi. This study is the first to demonstrate ancient admixture of lacustrine or pond organisms in Sulawesi beyond 100 km. The complex geological history of this island enabled such island-wide admixture of lacustrine organisms, which usually experience limited migration.

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.

Immature symbiotic system between horizontally transmitted green algae and brown hydra.

Some strains of brown hydra (Hydra vulgaris) are able to harbor the green algae Chlorococcum in their endodermal epithelial cells as symbionts. However, the relationship between brown hydra and chlorococcum is considered to be incipient symbiosis because most artificially introduced symbionts are not stable and because symbiotic H. vulgaris strains are rare in the wild. In this study, we compared the gene expression levels of the newly established symbiotic hydra (strain 105G), the native symbiotic strain (J7), and their non-symbiotic polyps to determine what changes would occur at the early stage of the evolution of symbiosis. We found that both the 105G and J7 strains showed comparable expression patterns, exhibiting upregulation of lysosomal enzymes and downregulation of genes related to nematocyte development and function. Meanwhile, genes involved in translation and the respiratory chain were upregulated only in strain 105G. Furthermore, treatment with rapamycin, which inhibits translation activity, induced the degeneration of the symbiotic strains (105G and J7). This effect was severe in strain 105G. Our results suggested that evolving the ability to balance the cellular metabolism between the host and the symbiont is a key requirement for adapting to endosymbiosis with chlorococcum.

A sleep-like state in Hydra unravels conserved sleep mechanisms during the evolutionary development of the central nervous system.

Sleep behaviors are observed even in nematodes and arthropods, yet little is known about how sleep-regulatory mechanisms have emerged during evolution. Here, we report a sleep-like state in the cnidarian Hydra vulgaris with a primitive nervous organization. Hydra sleep was shaped by homeostasis and necessary for cell proliferation, but it lacked free-running circadian rhythms. Instead, we detected 4-hour rhythms that might be generated by ultradian oscillators underlying Hydra sleep. Microarray analysis in sleep-deprived Hydra revealed sleep-dependent expression of 212 genes, including cGMP-dependent protein kinase 1 (PRKG1) and ornithine aminotransferase. Sleep-promoting effects of melatonin, GABA, and PRKG1 were conserved in Hydra However, arousing dopamine unexpectedly induced Hydra sleep. Opposing effects of ornithine metabolism on sleep were also evident between Hydra and Drosophila, suggesting the evolutionary switch of their sleep-regulatory functions. Thus, sleep-relevant physiology and sleep-regulatory components may have already been acquired at molecular levels in a brain-less metazoan phylum and reprogrammed accordingly.

Evidence for sympatric speciation in a Wallacean ancient lake.

Sympatric speciation has been demonstrated in few empirical case studies, despite intense searches, because of difficulties in testing the criteria for this mode of speciation. Here, we report a possible case of sympatric speciation in ricefishes of the genus Oryzias on Sulawesi, an island of Wallacea. Three species of Oryzias are known to be endemic to Lake Poso, an ancient tectonic lake in central Sulawesi. Phylogenetic analyses using RAD-seq-derived single nucleotide polymorphisms (SNPs) revealed that these species are monophyletic. We also found that the three species are morphologically distinguishable and clearly separated by population-structure analyses based on the SNPs, suggesting that they are reproductively isolated from each other. A mitochondrial DNA chronogram suggested that their speciation events occurred after formation of the tectonic lake, and existence of a historical allopatric phase was not supported by coalescent-based demographic inference. Demographic inference also suggested introgressive hybridization from an outgroup population. However, differential admixture among the sympatric species was not supported by any statistical tests. These results all concur with criteria necessary to demonstrate sympatric speciation. Ricefishes in this Wallacean lake provide a promising new model system for the study of sympatric speciation.

Effects of gene duplication, epistasis, recombination and gene conversion on the fixation time of compensatory mutations.

Gene duplication is one of the major mechanisms of molecular evolution. Gene duplication enables copies of a gene to accumulate mutations through functional redundancy. If a gene encodes a specific protein that interacts with other proteins, RNA, or DNA, the relaxation of selective constraints caused by gene duplication might contribute to the fixation of compensatory mutations that occur at the interacting sites. In this study, we investigate the effect of gene duplication, epistasis among the duplicated copies and gene conversion on the fixation time of compensatory mutations by extending the original model of compensatory evolution proposed by Kimura in 1985. Our simulation results reveal that the time to fixation of compensatory mutations can be decreased remarkably by gene duplication if one of the duplicated loci can completely mask the deleterious effects of a mutation that occurs at the other locus. Conversely, the fixation time can be increased by gene duplication if such functional compensation is weak. We also show that the combination of the degree of functional compensation and the rate of gene conversion between duplicate loci have contrasting effects on the time to fixation of compensatory mutations.

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.

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.

Evolutionary Origin of OwlRep, a Megasatellite DNA Associated with Adaptation of Owl Monkeys to Nocturnal Lifestyle.

Rod cells of many nocturnal mammals have a "non-standard" nuclear architecture, which is called the inverted nuclear architecture. Heterochromatin localizes to the central region of the nucleus. This leads to an efficient light transmission to the outer segments of photoreceptors. Rod cells of diurnal mammals have the conventional nuclear architecture. Owl monkeys (genus Aotus) are the only taxon of simian primates that has a nocturnal or cathemeral lifestyle, and this adaptation is widely thought to be secondary. Their rod cells were shown to exhibit an intermediate chromatin distribution: a spherical heterochromatin block was found in the central region of the nucleus although it was less complete than that of typical nocturnal mammals. We recently demonstrated that the primary DNA component of this heterochromatin block was OwlRep, a megasatellite DNA consisting of 187-bp-long repeat units. However, the origin of OwlRep was not known. Here we show that OwlRep was derived from HSAT6, a simple repeat sequence found in the centromere regions of human chromosomes. HSAT6 occurs widely in primates, suggesting that it was already present in the last common ancestor of extant primates. Notably, Strepsirrhini and Tarsiformes apparently carry a single HSAT6 copy, whereas many species of Simiiformes contain multiple copies. Comparison of nucleotide sequences of these copies revealed the entire process of the OwlRep formation. HSAT6, with or without flanking sequences, was segmentally duplicated in New World monkeys. Then, in the owl monkey linage after its divergence from other New World monkeys, a copy of HSAT6 was tandemly amplified, eventually forming a megasatellite DNA.

The Effect of Two Amino acid Residue Substitutions via RNA Editing on Dark-operative Protochlorophyllide Oxidoreductase in the Black Pine Chloroplasts.

Dark-operative protochlorophyllide oxidoreductase (DPOR) is a key enzyme to produce chlorophyll in the dark. Among photosynthetic eukaryotes, all three subunits chlL, chlN, and chlB are encoded by plastid genomes. In some gymnosperms, two codons of chlB mRNA are changed by RNA editing to codons encoding evolutionarily conserved amino acid residues. However, the effect of these substitutions on DPOR activity remains unknown. We first prepared cyanobacterial ChlB variants with amino acid substitution(s) to mimic ChlB translated from pre-edited mRNA. Their activities were evaluated by measuring chlorophyll content of dark-grown transformants of a chlB-lacking mutant of the cyanobacterium Leptolyngbya boryana that was complemented with pre-edited mimic chlB variants. The chlorophyll content of the transformant cells expressing the ChlB variant from the fully pre-edited mRNA was only one-fourth of the control cells. Co-purification experiments of ChlB with Strep-ChlN suggested that a stable complex with ChlN is greatly impaired in the substituted ChlB variant. We then confirmed that RNA editing efficiency was markedly greater in the dark than in the light in cotyledons of the black pine Pinus thunbergii. These results indicate that RNA editing on chlB mRNA is important to maintain appropriate DPOR activity in black pine chloroplasts.

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.