Geographic variation in the structure of oak hybrid zones provides insights into the dynamics of speciation

Does Asymmetric Reproductive Isolation Predict the Direction of Introgression in Plants?

The evolution of reproductive isolation (RI) results in the reduction of interspecific hybridization and the maintenance of species boundaries. Asymmetries in RI, where one species more frequently serves as the maternal or paternal parent in initial F1 hybrid formation, are commonly observed in plants. Asymmetric introgression, the predominantly unidirectional transfer of genetic material through hybridization and backcrossing, has also been frequently documented in hybridizing plant taxa as well. This study investigates whether asymmetries in total RI measured between species can predict the direction of introgression in naturally hybridizing plant taxa. A meta-analysis was conducted on 19 plant species pairs with published data on both asymmetric total RI, and asymmetric introgression. Species pairs that met these criteria were identified through a comprehensive literature review. A two-tailed binomial test was performed to evaluate whether asymmetric RI was associated with asymmetries in introgression. No significant relationship was found between asymmetries in total RI and the direction of introgression (p = 0.3593). Asymmetric RI largely does not predict the direction of introgression. Rather, introgression patterns may be better understood by examining F1 and later-generation hybrids in natural settings, focusing on their fitness, mating behaviors, and the ecological and demographic factors that shape hybrid zones.

Lipid and Amino Acid Pathway Metabolites Contribute to Cold Tolerance in Quercus wutaishanica

Cold is an important environmental stress affecting the growth, productivity, and geographic distribution of tree species. Oaks are important for environmental conservation and wood supplies. Oak metabolites respond to low temperatures (LTs). In this study, the physiological and metabolic responses of two oak species to cold stress were investigated and compared. The field observations and physiological responses showed that Quercus wutaishanica was more cold-tolerant than Q. acutissima. After frost, the one-year-old twigs of Q. wutaishanica had higher survival rates, accumulated more soluble sugar and protein, and exhibited higher superoxide dismutase (SOD) activity than those of Q. acutissima. Untargeted metabolomics identified 102 and 78 differentially accumulated metabolites in Q. acutissima and Q. wutaishanica, respectively, when the leaves were subjected to LTs (4 °C for 24 h). The carbohydrate and flavonoid metabolites contributed to the cold tolerance of both oak species. Succinate, an intermediate in the citric acid cycle, was significantly inhibited by LTs, a potential energy conservation strategy. Unlike Q. acutissima, Q. wutaishanica underwent metabolic reprogramming that significantly increased the contents of phosphatidylcholine, gallic acid, oxidized glutathione, shikimate, and phenylpyruvate under LTs. Our data provide a reference for characterizing the mechanisms involved in the response of oak species to cold temperatures and enhancing the cold tolerance of forest trees.

Genetic structure and trait variation within a maple hybrid zone underscore North China as an overlooked diversity hotspot

Tertiary relict flora in East Asia can be divided into northern and southern regions. North China is a diversity hotspot because it can be the secondary contact zone of ancient lineages from the two regions. To test the extent of ancient lineages hybridization and distinguish between the putative species pair Acer pictum subsp. mono and Acer truncatum, we conducted genetic and ecological studies within a maple hybrid zone in North China. Our results suggest that the two lineages of Acer coexist in the hybrid zone and that adult and offspring populations show typical bimodal genetic patterns. Hybrid individuals are established at intermediate altitudes between the two parental lineages. Flowering phenology is divergent between lineages, whereas the complex sexual system of Acer may ensure pollination among lineages. Leaf and fruit morphologies are different between the northern and southern origin lineages, corresponding to A. pictum subsp. mono and A. truncatum, respectively. Reduced gene flow between lineages suggests that they should be considered as two species. However, large morphological variations within each species and the existence of hybrids offer low reliability of species identification based solely on morphological traits. Our study underscores North China as an overlooked diversity hotspot that requires further study in the future.

Genomic Data Reveals Profound Genetic Structure and Multiple Glacial Refugia in Lonicera oblata (Caprifoliaceae), a Threatened Montane Shrub Endemic to North China

Characterizing genetic diversity and structure and identifying conservation units are both crucial for the conservation and management of threatened species. The development of high-throughput sequencing technology provides exciting opportunities for conservation genetics. Here, we employed the powerful SuperGBS method to identify 33, 758 high-quality single-nucleotide polymorphisms (SNP) from 134 individuals of a critically endangered montane shrub endemic to North China, Lonicera oblata. A low level of genetic diversity and a high degree of genetic differentiation among populations were observed based on the SNP data. Both principal component and phylogenetic analyses detected seven clusters, which correspond exactly to the seven geographic populations. Under the optimal K = 7, Admixture suggested the combination of the two small and geographically neighboring populations in the Taihang Mountains, Dongling Mountains, and Lijiazhuang, while the division of the big population of Jiankou Great Wall in the Yan Mountains into two clusters. High population genetic diversity and a large number of private alleles were detected in the four large populations, while low diversity and non-private alleles were observed for the remaining three small populations, implying the importance of these large populations as conservation units in priority. Demographic history inference suggested two drastic contractions of population size events that occurred after the Middle Pleistocene Transition and the Last Glacial Maximum, respectively. Combining our previous ecological niche modeling results with the present genomic data, there was a possible presence of glacial refugia in the Taihang and Yan Mountains, North China. This study provides valuable data for the conservation and management of L. oblata and broadens the understanding of the high biodiversity in the Taihang and Yan Mountains.

Influence of Pliocene and Pleistocene climates on hybridization patterns between two closely related oak species in China

BACKGROUND AND AIMS

Contemporary patterns of genetic admixture reflect imprints of both ancient and recent gene flow, which can provide us with valuable information on hybridization history in response to palaeoclimate change. Here, we examine the relationships between present admixture patterns and past climatic niche suitability of two East Asian Cerris oaks (Quercus acutissima and Q. chenii) to test the hypothesis that the mid-Pliocene warm climate promoted while the Pleistocene cool climate limited hybridization among local closely related taxa.

METHODS

We analyse genetic variation at seven nuclear microsatellites (1111 individuals) and three chloroplast intergenic spacers (576 individuals) to determine the present admixture pattern and ancient hybridization history. We apply an information-theoretic model selection approach to explore the associations of genetic admixture degree with past climatic niche suitability at multiple spatial scales.

KEY RESULTS

More than 70 % of the hybrids determined by Bayesian clustering analysis and more than 90 % of the individuals with locally shared chloroplast haplotypes are concentrated within a mid-Pliocene contact zone between ~30°N and 35°N. Climatic niche suitabilities for Q. chenii during the mid-Pliocene Warm Period [mPWP, ~3.264-3.025 million years ago (mya)] and during the Last Glacial Maximum (LGM, ~0.022 mya) best explain the admixture patterns across all Q. acutissima populations and across those within the ancient contact zone, respectively.

CONCLUSIONS

Our results highlight that palaeoclimate change shapes present admixture patterns by influencing the extent of historical range overlap. Specifically, the mid-Pliocene warm climate promoted ancient contact, allowing widespread hybridization throughout central China. In contrast, the Pleistocene cool climate caused the local extinction of Q. chenii, reducing the probability of interspecific gene flow in most areas except those sites having a high level of ecological stability.

Different population size change and migration histories created genetic diversity of three oaks in Tokai region, central Japan

To understand genetic diversity in focal species, it is important to consider the possibility of speciation with gene flow, especially in species with porous genomes such as oaks. We studied genetic diversity and structure in three oaks, Quercus mongolica var. mongolicoides (QM), Q. mongolica var. crispula (QC) and Q. serrata (QS), growing in the Tokai region, central Japan. QM is semi-endemic to the region while the others are common taxa. We also conducted demographic modeling to infer their population size change and migration histories using an approximate Bayesian computation (ABC) approach. The three taxa showed distinct genetic structures but there was genetic admixture among the taxa, especially between QM and QC. ABC analysis of population size change revealed that the population size of QM was stable during and after the last glacial period, while QC and QS showed population expansion after the last glacial maximum. ABC analysis of population divergence and migration revealed that continuous gene flow between QM and QC after their divergence was supported, while between QM and QS, and between QC and QS, secondary contact after sufficient isolation was supported. These historical migration patterns among the three taxa indicate that QM and QC are currently in the early stage or gray zone of speciation, whereas speciation of the other two taxon pairs is considered to have almost been established. Observed gene flow patterns and strength between QM and QC, and between QM and QS, were explained by both flowering patterns and historical distributions, but those between QC and QS were not.

Microsatellites and petal morphology reveal new patterns of admixture in Orchis hybrid zones

PREMISE

The genetic structure of hybrid zones provides insight into the potential for gene flow to occur between plant taxa. Four closely related European orchid species (Orchis anthropophora, O. militaris, O. purpurea, and O. simia) hybridize when they co-occur. We aimed to characterize patterns of hybridization in O. militaris-O. purpurea, O. purpurea-O. simia, and O. anthropophora-O. simia hybrid zones using molecular and morphological data.

METHODS

We used 11 newly isolated nuclear microsatellites to genotype 695 individuals collected from seven hybrid zones and six allopatric parental populations in France. Geometric morphometric analysis was conducted using 15 labellum landmarks to capture the main aspects of petal shape.

RESULTS

Backcrossing was asymmetric toward O. militaris in multiple O. militaris-O. purpurea hybrid zones. Hybrids in O. purpurea-O. simia and O. anthropophora-O. simia hybrid zones were largely limited to F1 and F2 generations, but further admixture had occurred. These patterns were reflected in labellum geometric morphometric data, which correlated strongly with nuclear microsatellite data in all three species combinations.

CONCLUSIONS

The coexistence of parental and admixed individuals in these Orchis hybrid zones implies they are likely to be tension zones being maintained by a balance between gene flow into the hybrid zone and selection acting against admixed individuals. The pattern of admixture in the three species combinations suggests intrinsic selection acting on the hybrids is weaker in more closely related taxa.

Recent demographic histories of temperate deciduous trees inferred from microsatellite markers

Background

Accurate inference of demographic histories for temperate tree species can aid our understanding of current climate change as a driver of evolution. Microsatellites are more suitable for inferring recent historical events due to their high mutation rates. However, most programs analyzing microsatellite data assume a strict stepwise mutation model (SMM), which could cause false detection of population shrinkage when microsatellite mutation does not follow SMM.

Results

This study aims to reconstruct the recent demographic histories of five cool-temperate tree species in Eastern Asia, Quercus mongolica, Q. liaotungensis, Juglans cathayensis, J. mandshurica and J. ailantifolia, by using 19 microsatellite markers with two methods considering generalized stepwise mutation model (GSM) (MIGRAINE and VarEff). Both programs revealed that all the five species experienced expansions after the Last Glacial Maximum (LGM). Within butternuts, J. cathayensis experienced a more serious bottleneck than the other species, and within oaks, Q. mongolica showed a moderate increase in population size and remained stable after the expansion. In addition, the point estimates of the multistep mutation proportion in the GSM model (p_GSM) for all five species were between 0.50 and 0.65, indicating that when inferring population demographic history of the cool-temperate forest species using microsatellite markers, it is better to assume a GSM rather than a SMM.

Conclusions

This study provides the first direct evidence that five cool-temperate tree species in East Asia have experienced expansions after the LGM with microsatellite data. Considering the mutation model of microsatellite has a vital influence on demographic inference, combining multiple programs such as MIGRAINE and VarEff can effectively reduce errors caused by inappropriate model selection and prior setting.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01805-w.

Unraveling the genetic diversity and structure of Quercus liaotungensis population through analysis of microsatellite markers

Background

Quercus liaotungensis Koidz. is an ecologically and economically important tree species widely distributed in Northern China. However, the effective assessment, utilization, and protection of Q. liaotungensis resources remain unexplored.

Methods

In total, 120 samples obtained from 12 Q. liaotungensis populations of Northern China were investigated for genetic diversity and structure using 19 simple sequence repeat (SSR) primer pairs.

Results

The total number of alleles detected was 293, the average number of effective allele (Ne) was 6.084, the genetic differentiation coefficient (Fst) was 0.033, and the mean observed heterozygosity (Ho) and expected heterozygosity (He) were 0.690 and 0.801, respectively. Moreover, analysis of molecular variance (AMOVA) showed a 5.5% genetic variation among 12 Q. liaotungensis populations, indicating that a high level of genetic diversity and a low degree of genetic differentiation among Q. liaotungensis populations. STRUCTURE and cluster analysis divided the 12 Q. liaotungensis populations into the following three subpopulations: Bashang Plateau subpopulation (SH), Liaodong Peninsula subpopulation (NC), and Loess Plateau subpopulation (other 10 populations). The cluster analysis based on 19 climatic factors was consistent with the genetic structure. A positive correlation was found between genetic distance and geographical distance (r = 0.638, p = 0.028) by the Mantel test, and two boundaries were found among the 12 Q. liaotungensis populations by the Barrier analysis, indicating that Q. liaotungensis populations existed isolated by geographical distance and physical barrier.

Conclusion

This study suggests that geographical isolation, physical barrier, climatic types, and natural hybridization promote the formation of genetic structures, which can contribute to future protection and genetic improvement of Q. liaotungensis.

Phylogenetic relationships in Chinese oaks (Fagaceae, Quercus): Evidence from plastid genome using low-coverage whole genome sequencing

China is a second center of oak diversity but with less intensively systematic studies. Here, with 49 species representing all four sections in China, we firstly gave insight into the comprehensive phylogenetic relationships of Chinese oaks based on 54 complete plastid genomes. Our results recovered a robust phylogenetic framework and provided strong support for most nodes. The phylogenetic tree supported Quercus section Ilex as not monophyletic, in which Quercus section Cyclobalanopsis and Quercus section Cerris were nested. Most likely, incomplete lineage sorting and/or introgression among ancestral lineages in these three sections resulted in this complex pattern. The current distribution, diversification and molecular differentiation of Q. sect. Ilex in China are likely consequences of local adaptation to the geographic and paleoclimatic changes, which were driven by the uplift of Tibetan Plateau, the Hengduan Mountains and the Himalayas.

Insights Into the Significance of the Chinense Loess Plateau for Preserving Biodiversity From the Phylogeography of Speranskia tuberculata (Euphorbiaceae)

The significance of the Chinese Loess Plateau (CLP) in maintaining biodiversity for northern China has rarely been shown, as previous phylogeographic studies are mostly woody species and they have revealed that Quaternary refugia are mainly located in mountain regions. We selected a drought-enduring endemic herb, Speranskia tuberculata (Euphorbiaceae), to determine its glacial refugia and postglacial demographic history. To this end, we sampled 423 individuals from 38 populations covering its entire geographic distribution. Three chloroplast DNA (cpDNA) fragments, two low-copy nuclear genes, and six nuclear microsatellites (nSSRs) were used and supplemented with ecological niche modeling (ENM) to infer the phylogeographic history of this species. Populations with private haplotypes and high haplotype diversity of cpDNA are mainly located in the CLP or scattered around northeastern China and the coastal region. Spatial expansion, detected using a neutrality test and mismatch distribution, may have resulted in a widely distributed ancestral cpDNA haplotype, especially outside of the CLP. For nuclear DNA, private haplotypes are also distributed mainly in the CLP. In nSSRs, STRUCTURE clustering identified two genetic clusters, which are distributed in the west (western cluster) and east (eastern cluster), respectively. Many populations belonged, with little to no admixture, to the western cluster while (hardly) pure populations of the eastern cluster were barely found. Genetic differentiation is significantly correlated with geographic distance, although genetic diversity is uniformly distributed. ENM suggests that the distribution of S. tuberculata has recently expanded northwards from the southern CLP, whereas it has experienced habitat loss in the south. Thus, S. tuberculata populations probably survived the last glacial maximum (LGM) in the southern CLP and experienced post-glacial expansion. Wind-dispersed pollen could bring the majority of genotypes to the front during spatial expansion, resulting in uniformly distributed genetic diversity. Based on evidence from molecular data and vegetation and climate changes since the LGM, we conclude that drought-enduring species, especially herbaceous species, are likely to have persisted in the CLP during the LGM and to have experienced expansion to other regions in northern China.

Looks can be deceiving: speciation dynamics of co-distributed Angophora (Myrtaceae) species in a varying landscape

Understanding the mechanisms underlying species divergence remains a central goal in evolutionary biology. Landscape genetics can be a powerful tool for examining evolutionary processes. We used genome-wide scans to genotype samples from populations of eight Angophora species. Angophora is a small genus within the eucalypts comprising common and rare species in a heterogeneous landscape, making it an appropriate group to study speciation. We found A. hispida was highly differentiated from the other species. Two subspecies of A. costata (subsp. costata and subsp. euryphylla) formed a group, while the third (subsp. leiocarpa, which is only distinguished by its smooth fruits and provenance) was supported as a distinct pseudocryptic species. Other species that are morphologically distinct could not be genetically differentiated (e.g., A. floribunda and A. subvelutina). Distribution and genetic differentiation within Angophora were strongly influenced by temperature and humidity, as well as biogeographic barriers, particularly rivers and higher elevation regions. While extensive introgression was found between many populations of some species (e.g., A. bakeri and A. floribunda), others only hybridized at certain locations. Overall, our findings suggest multiple mechanisms drove evolutionary diversification in Angophora and highlight how genome-wide analyses of related species in a diverse landscape can provide insights into speciation.

Salinity, not genetic incompatibilities, limits the establishment of the invasive hybrid cattail Typha × glauca in coastal wetlands

Hybrids of a single pair of parent species can be much more common in some geographical regions than in others. The reasons for this are not well understood, but could help explain processes such as species diversification or the range expansion of invasive hybrids. The widespread cattails Typha latifolia and T. angustifolia seldom hybridize in some parts of their range, but in other areas produce the dominant hybrid T. × glauca. We used a combination of field and greenhouse experiments to investigate why T. × glauca has invaded wetlands in the Laurentian Great Lakes region of southern Ontario, Canada, but is much less common in the coastal wetlands of Nova Scotia (NS) in eastern Canada. One potentially important environmental difference between these two regions is salinity. We therefore tested three hypotheses: (1) T. latifolia and T. angustifolia in NS are genetically incompatible; (2) the germination or growth of T. × glauca is reduced by salinity; and (3) T. latifolia, a main competitor of T. × glauca, is locally adapted to saline conditions in NS. Our experiments showed that NS T. latifolia and T. angustifolia are genetically compatible, and that saline conditions do not impede growth of hybrid plants. However, we also found that under conditions of high salinity, germination rates of hybrid seeds were substantially lower than those of NS T. latifolia. In addition, germination rates of NS T. latifolia were higher than those of Ontario T. latifolia, suggesting local adaptation to salinity in coastal wetlands. This study adds to the growing body of literature which identifies the important roles that local habitat and adaptation can play in the distributions and characteristics of hybrid zones.

Dispersal as a result of asymmetrical hybridization between two closely related oak species in China

Hybridization has played an important role in plant evolution. Less attention has been paid, however, to its role in dispersal. In this study, historical divergence and hybridization were investigated in two closely related Chinese oaks, Quercus mongolica and Q. liaotungensis, to estimate the role that hybridization played in their dispersal. We genotyped 27 Q. mongolica and Q. liaotungensis populations throughout the distributional range of the two oak species, using 14 single-copy nuclear genes and four noncoding chloroplast DNA regions. Bayesian cluster and population tree analyses indicated that there were three groups over all oak populations, namely, Q. mongolica, northwest-northern China (NW-NC) Q. liaotungensis, and northeastern China (NEC) Q. liaotungensis. Approximate Bayesian computation simulation supported an asymmetrical hybridization origin of NEC Q. liaotungensis, after a previous divergence between NW-NC Q. liaotungensis and Q. mongolica. IMa3 analyses suggested that Q. liaotungensis and Q. mongolica diverged in the NW-NC and NEC regions, respectively, and that NEC Q. liaotungensis arose from Q. mongolica, not from NW-NC Q. liaotungensis, and was greatly introgressed by NW-NC Q. liaotungensis. Oak populations in NW-NC and NEC regions held different chloroplast DNA haplotypes, and Q. liaotungensis in NEC shared most haplotypes with Q. mongolica populations, but none with NW-NC Q. liaotungensis populations, suggesting the maternal origin of NEC Q. liaotungensis from Q. mongolica. This study found clear signals of isolation divergence of Q. liaotungensis in NW-NC and Q. mongolica in NEC, and the results suggest that asymmetrical hybridization and introgression from Q. liaotungensis to Q. mongolica, mostly likely via pollen flow, facilitated Q. liaotungensis dispersal to NEC.

Framework Phylogeny, Evolution and Complex Diversification of Chinese Oaks

Oaks (Quercus L.) are ideal models to assess patterns of plant diversity. We integrated the sequence data of five chloroplast and two nuclear loci from 50 Chinese oaks to explore the phylogenetic framework, evolution and diversification patterns of the Chinese oak's lineage. The framework phylogeny strongly supports two subgenera Quercus and Cerris comprising four infrageneric sections Quercus, Cerris, Ilex and Cyclobalanopsis for the Chinese oaks. An evolutionary analysis suggests that the two subgenera probably split during the mid-Eocene, followed by intergroup divergence within the subgenus Cerris around the late Eocene. The initial diversification of sections in the subgenus Cerris was dated between the mid-Oligocene and the Oligocene-Miocene boundary, while a rapid species radiation in section Quercus started in the late Miocene. Diversification simulations indicate a potential evolutionary shift on section Quercus, while several phenotypic shifts likely occur among all sections. We found significant negative correlations between rates of the lineage diversification and phenotypic turnover, suggesting a complex interaction between the species evolution and morphological divergence in Chinese oaks. Our infrageneric phylogeny of Chinese oaks accords with the recently proposed classification of the genus Quercus. The results point to tectonic activity and climatic change during the Tertiary as possible drivers of evolution and diversification in the Chinese oak's lineage.

The Role of Population and Half-Sib Family on Driving Suitable Functional Traits for Quercus suber L. Forest Restoration

Research Highlights: Seedlings of different Quercus suber L. populations and half-sib families differ in their response to multiple stressors, which may have consequences on the future distribution of this Mediterranean species. Background and Objectives: Global change will likely increase the frequency and severity of drought in drylands. Plant species’ distributions will largely depend on their ability to respond to the combined effect of drought and other environmental stressors. Genetic diversity in morpho-functional traits are key components of this response. Yet, information on the response to multiple stresses is scarce for many iconic species. The present study assessed the variability in the response of populations and half-sib families of a Mediterranean acidophilous tree, cork oak, to drought and changes in soil conditions. Materials and Methods: We sampled acorns of half-sib families from two cork oak populations genetically introgressed with the alkaline-tolerant species Quercus ilex L., and from a non-introgressed cork oak population located in its core habitat. We germinated the acorns and subjected seedlings to contrasted levels of water availability and additions of calcium and magnesium carbonate, and assessed their morpho-physiological response. Results: Response to drought and soil chemistry composition differed between populations and families. For some traits, introgressed populations responded similarly to drought than the non-introgressed population. Conversely, the response to soil chemistry was not clearly related to introgression. When considering half-sib families within populations, the population effect diminished, which revealed the importance of intra-population variation. However, relevant traits for water scarcity adaptations, such as specific leaf area and root:shoot ratio, remained significantly different at the population level, which highlights the relevance of these traits for management. Conclusions: Our study shows that the adaptive management and restoration of cork oak forests should consider not only geographic provenances, but also half-sib lines within populations.

The effect of phylogeographic history on species boundaries: a comparative framework in Hyla tree frogs

Because it is indicative of reproductive isolation, the amount of genetic introgression across secondary contact zones is increasingly considered in species delimitation. However, patterns of admixture at range margins can be skewed by the regional dynamics of hybrid zones. In this context, we posit an important role for phylogeographic history: hybrid zones located within glacial refugia (putatively formed during the Late-Pleistocene) should be better defined than those located in post-glacial or introduced ranges (putatively formed during the Holocene and the Anthropocene). We test this hypothesis in a speciation continuum of tree frogs from the Western Palearctic (Hyla), featuring ten identified contacts between species spanning Plio-Pleistocene to Miocene divergences. We review the rich phylogeographic literature of this group and examine the overlooked transition between H. arborea and H. molleri in Western France using a multilocus dataset. Our comparative analysis supports a trend that contacts zones resulting from post-glacial expansions and human translocations feature more extensive introgression than those established within refugial areas. Integrating the biogeographic history of incipient species, i.e. their age since first contact together with their genetic divergence, thus appears timely to draw sound evolutionary and taxonomic inferences from patterns of introgression across hybrid zones.

ddRAD Sequencing-Based Identification of Genomic Boundaries and Permeability in Quercus ilex and Q. suber Hybrids

Hybridization and its relevance is a hot topic in ecology and evolutionary biology. Interspecific gene flow may play a key role in species adaptation to environmental change, as well as in the survival of endangered populations. Despite the fact that hybridization is quite common in plants, many hybridizing species, such as Quercus spp., maintain their integrity, while precise determination of genomic boundaries between species remains elusive. Novel high throughput sequencing techniques have opened up new perspectives in the comparative analysis of genomes and in the study of historical and current interspecific gene flow. In this work, we applied ddRADseq technique and developed an ad hoc bioinformatics pipeline for the study of ongoing hybridization between two relevant Mediterranean oaks, Q. ilex and Q. suber. We adopted a local scale approach, analyzing adult hybrids (sensu lato) identified in a mixed stand and their open-pollinated progenies. We have identified up to 9,435 markers across the genome and have estimated individual introgression levels in adults and seedlings. Estimated contribution of Q. suber to the genome is higher, on average, in hybrid progenies than in hybrid adults, suggesting preferential backcrossing with this parental species, maybe followed by selection during juvenile stages against individuals with higher Q. suber genomic contribution. Most discriminating markers seem to be scattered throughout the genome, suggesting that a large number of small genomic regions underlie boundaries between these species. A noticeable proportion of the markers (26%) showed allelic frequencies in adult hybrids very similar to one of the parental species, and very different from the other; a finding that seems relevant for understanding the hybridization process and the occurrence of adaptive introgression. Candidate marker databases developed in this study constitute a valuable resource to design large scale re-sequencing experiments in Mediterranean sclerophyllous oak species and could provide insight in species boundaries and on adaptive introgression between Q. suber and Q. ilex.

A comparison of reproductive isolation between two closely related oak species in zones of recent and ancient secondary contact

BACKGROUND

Much of the debate over the evolutionary consequences of hybridization on genetic divergence and speciation results from the breakdown or reinforcement of reproductive barriers in secondary hybrid zones. Among hybrid populations established for different lengths of time following secondary contact, stronger reproductive barriers are generally expected to occur in zones with longer contact. However, in plants no detailed investigation of recent and ancient zones of secondary contact has been conducted despite the importance of such a comparative study. Here, we compare pre- and postzygotic reproductive barriers between two closely related oak species, Quercus mongolica and Q. liaotungensis, in such a situation.

RESULTS

The recorded flowering times of both species overlapped in both contact zones. The fruit set at 10 and 30 days after interspecific hand pollination was not significantly lower than that after intraspecific pollination whenever Q. mongolica or Q. liaotungensis comprised the maternal parents in both populations. These results indicated that neither prezygotic phenological barriers nor interspecific incompatibility could have resulted in the reproductive isolation between the two species in both hybrid zones. However, the proportion of hybrid seeds produced by both species in the ancient zone was significantly lower than that recorded in the recent zone of secondary contact. In addition, the proportion of hybrid seeds simulated to form, assuming both random mating and an absence of postpollination barriers, was significantly higher than that detected in the ancient contact zone but not in the recent contact zone. These results suggest stronger early-acting postzygotic isolation between the two oak species in the ancient relative to the recent contact zone.

CONCLUSIONS

Our comparative study demonstrated that postzygotic barriers during seed maturity were the main contributing factor to total reproductive isolation, particularly in the ancient contact zone, which aided species delimitation. In the recently formed secondary contact zone, pre- and postzygotic barriers were not well developed, and a high frequency of natural hybridization was evident. To our knowledge this study provides the first comparison of reproductive isolation between the ancient and recent secondary contact zones in plants and helps to clarify the evolutionary consequences of hybridization in a temporal context.

The Phytogeographic History of Common Walnut in China

Common walnut (Juglans regia L.) is an economically important hardwood tree species cultivated worldwide for its high quality wood and edible nuts. It is generally accepted that after the last glaciation J. regia survived and grew in almost completely isolated stands in Asia, and that ancient humans dispersed walnuts across Asia and into new habitats via trade and cultural expansion. The history of common walnut in China is a matter of debate, however. We estimated the genetic diversity and spatial genetic structure of 31 walnut populations sampled across its Chinese range using 22 microsatellite markers (13 neutral and 9 non-neutral). Using historical data and population genetic analysis, including approximate Bayesian analysis (ABC), we reconstructed the demographic history of J. regia in China. The genetic data indicated the likely presence of J. regia in glacial refugia in the Xinjiang province (Northwest China), Northeastern China (Beijing, Shandong, and Changbai Mountains), Central China (Qinling and Baishan Mountains and Xi'an), and Southwestern China (Tibet, Yunnan, Guizhou, and Sichuan provinces). Based on DIY-ABC analysis, we identified three ancient lineages of J. regia in China. Two lineages (subpopulation A and subpopulation B+C) diverged about 2.79 Mya, while Southwestern China, and Qinling and Baishan Mountains lineages diverged during the Quaternary glaciations (about 1.13 Mya). Remnants of these once-distinct genetic clusters of J. regia may warrant ecological management if they are to be retained as in situ resources. A population size expansion in Northeastern China was detected in the last five centuries. The present distribution of walnut in China resulted from the combined effects of expansion/contraction from multiple refugia after the Last Glacial Maximum and later human exploitation.

Species Boundaries Between Three Sympatric Oak Species: Quercus aliena, Q. dentata, and Q. variabilis at the Northern Edge of Their Distribution in China

Oaks are important timber trees with wide distributions in China, but few genetic studies have been conducted on a fine scale. In this study, we seek to investigate the genetic diversity and differentiation of three sympatric oak species (Quercus aliena Blume, Quercus dentata Thunb. ex Murray, and Quercus variabilis Blume) in their northern distribution in China using 17 bi-parentally inherited nSSRs markers and five maternally inherited chloroplast DNA (cpDNA) fragments. Both the cpDNA and the nSSRs show a high level of genetic differentiation between different oak sections. The chloroplast haplotypes are clustered into two lineages. Clear species boundaries are detected between Q. variabilis and either Q. aliena or Q. dentata. The sharing of chloroplast haplotype H1 between Q. aliena and Q. dentata suggests very recent speciation and incomplete lineage sorting or introgression of H1 from one species to another. The nSSRs data indicate a complete fixation of variation within sites for all three oak species, and that extensive gene flow occurs within species whereas only limited gene flow is detected between Q. aliena and Q. dentata and nearly no gene flow can be detected between Q. aliena and Q. variabilis and between Q. dentata and Q. variabilis. Prezygotic isolation may have contributed to the species boundaries of these three sympatric oak species.

Plastid Genome Comparative and Phylogenetic Analyses of the Key Genera in Fagaceae: Highlighting the Effect of Codon Composition Bias in Phylogenetic Inference

Fagaceae is one of the largest and economically important taxa within Fagales. Considering the incongruence among inferences from plastid and nuclear genes in the previous Fagaceae phylogeny studies, we assess the performance of plastid phylogenomics in this complex family. We sequenced and assembled four complete plastid genomes (Fagus engleriana, Quercus spinosa, Quercus aquifolioides, and Quercus glauca) using reference-guided assembly approach. All of the other 12 published plastid genomes in Fagaceae were retrieved for genomic analyses (including repeats, sequence divergence and codon usage) and phylogenetic inference. The genomic analyses reveal that plastid genomes in Fagaceae are conserved. Comparing the phylogenetic relationships of the key genera in Fagaceae inferred from different codon positions and gene function datasets, we found that the first two codon sites dataset recovered nearly all relationships and received high support. Thus, the result suggested that codon composition bias had great influence on Fagaceae phylogenetic inference. Our study not only provides basic understanding of Fagaceae plastid genomes, but also illuminates the effectiveness of plastid phylogenomics in resolving relationships of this intractable family.

Climatic and Soil Factors Shape the Demographical History and Genetic Diversity of a Deciduous Oak (Quercus liaotungensis) in Northern China

Past and current climatic changes have affected the demography, patterns of genetic diversity, and genetic structure of extant species. The study of these processes provides valuable information to forecast evolutionary changes and to identify conservation priorities. Here, we sequenced two functional nuclear genes and four chloroplast DNA regions for 105 samples from 21 populations of Quercus liaotungensis across its distribution range. Coalescent-based Bayesian analysis, approximate Bayesian computation (ABC), and ecological niche modeling (ENM) were integrated to investigate the genetic patterns and demographical history of this species. Association estimates including Mantel tests and multiple linear regressions were used to infer the effects of geographical and ecological factors on temporal genetic variation and diversity of this oak species. Based on multiple loci, Q. liaotungensis populations clustered into two phylogenetic groups; this grouping pattern could be the result of adaptation to habitats with different temperature and precipitation seasonality conditions. Demographical reconstructions and ENMs suggest an expansion decline trend of this species during the Quaternary climatic oscillations. Association analyses based on nuclear data indicated that intraspecific genetic differentiation of Q. liaotungensis was clearly correlated with ecological distance; specifically, the genetic diversity of this species was significantly correlated with temperature seasonality and soil pH, but negatively correlated with precipitation. Our study highlights the impact of Pleistocene climate oscillations on the demographic history of a tree species in Northern China, and suggests that climatic and soil conditions are the major factors shaping the genetic diversity and population structure of Q. liaotungensis.

Asymmetrical natural hybridization varies among hybrid swarms between two diploid Rhododendron species

Background and Aims

The extent to which hybridization leads to gene flow between plant species depends on the structure of hybrid populations. However, if this varies between locations, species barriers might prove permeable in some locations but not in others. To assess possible variation in hybrid population structure, the magnitude and direction of natural hybridization between two Chinese endemic species, Rhododendron spiciferum and Rhododendron spinuliferum , were evaluated.

Methods

Thirteen nuclear microsatellite markers were employed to characterize 566 individuals collected from 15 non-allopatric populations and nine allopatric parental populations. Chloroplast DNA (cpDNA) sequences were obtained from a subset of samples. Genetic structure and direction of gene flow was determined using a combination of STRUCTURE and NEWHYBRIDS analysis.

Key Results

Nuclear analysis revealed that parental taxa formed two genetically distinct clusters and hybrids shared the genetic background of both parents and did not form a separate genetic lineage. Overall, hybrid swarms were dominated by early- and later-generation hybrids, with a significantly higher proportion of hybrids (59·6 %) possessing >50 % R. spiciferum-like nuclear germplasm. The cpDNA analysis further indicated that a significantly greater proportion of hybrids (61·1 %) possessed the R. spiciferum cpDNA haplotype.

Conclusions

Gene flow between R. spiciferum and R. spinuliferum was found to be bidirectional in 14 of the 15 hybrid swarms and asymmetrical in six hybrid swarms. Asymmetrical gene flow was evident for only nuclear DNA (nDNA) in two populations, for only cpDNA in three populations, and for both nDNA and cpDNA in one population. Collectively, the variation in genetic structure found among the 15 hybrid swarms suggests that introgression rather than hybrid speciation is a more likely outcome of hybridization between these hybridizing taxa.

Development of Chloroplast and Nuclear DNA Markers for Chinese Oaks (Quercus Subgenus Quercus) and Assessment of Their Utility as DNA Barcodes

Chloroplast DNA (cpDNA) is frequently used for species demography, evolution, and species discrimination of plants. However, the lack of efficient and universal markers often brings particular challenges for genetic studies across different plant groups. In this study, chloroplast genomes from two closely related species (Quercus rubra and Castanea mollissima) in Fagaceae were compared to explore universal cpDNA markers for the Chinese oak species in Quercus subgenus Quercus, a diverse species group without sufficient molecular differentiation. With the comparison, nine and 14 plastid markers were selected as barcoding and phylogeographic candidates for the Chinese oaks. Five (psbA-trnH, matK-trnK, ycf3-trnS, matK, and ycf1) of the nine plastid candidate barcodes, with the addition of newly designed ITS and a single-copy nuclear gene (SAP), were then tested on 35 Chinese oak species employing four different barcoding approaches (genetic distance-, BLAST-, character-, and tree-based methods). The four methods showed different species identification powers with character-based method performing the best. Of the seven barcodes tested, a barcoding gap was absent in all of them across the Chinese oaks, while ITS and psbA-trnH provided the highest species resolution (30.30%) with the character- and BLAST-based methods, respectively. The six-marker combination (psbA-trnH + matK-trnK + matK + ycf1 + ITS + SAP) showed the best species resolution (84.85%) using the character-based method for barcoding the Chinese oaks. The barcoding results provided additional implications for taxonomy of the Chinese oaks in subg. Quercus, basically identifying three major infrageneric clades of the Chinese oaks (corresponding to Groups Quercus, Cerris, and Ilex) referenced to previous phylogenetic classification of Quercus. While the morphology-based allocations proposed for the Chinese oaks in subg. Quercus were challenged. A low variation rate of the chloroplast genome, and complex speciation patterns involving incomplete lineage sorting, interspecific hybridization and introgression, possibly have negative impacts on the species assignment and phylogeny of oak species.

Impacts of human-induced environmental disturbances on hybridization between two ecologically differentiated Californian oak species

Natural hybridization, which can be involved in local adaptation and in speciation processes, has been linked to different sources of anthropogenic disturbance. Here, we use genotypic data to study range-wide patterns of genetic admixture between the serpentine-soil specialist leather oak (Quercus durata) and the widespread Californian scrub oak (Quercus berberidifolia). First, we estimated hybridization rates and the direction of gene flow. Second, we tested the hypothesis that genetic admixture increases with different sources of environmental disturbance, namely anthropogenic destruction of natural habitats and wildfire frequency estimated from long-term records of fire occurrence. Our analyses indicate considerable rates of hybridization (> 25%), asymmetric gene flow from Q. durata into Q. berberidifolia, and a higher occurrence of hybrids in areas where both species live in close parapatry. In accordance with the environmental disturbance hypothesis, we found that genetic admixture increases with wildfire frequency, but we did not find a significant effect of other sources of human-induced habitat alteration (urbanization, land clearing for agriculture) or a suite of ecological factors (climate, elevation, soil type). Our findings highlight that wildfires constitute an important source of environmental disturbance, promoting hybridization between two ecologically well-differentiated native species.

Phylogeography Study of the Siberian Apricot (Prunus sibirica L.) in Northern China Assessed by Chloroplast Microsatellite and DNA Makers

There is evidence that a band of dry climate separated plants in East Asia into distinct northern and southern groups. However, few studies have focused on the arid belt in this region, especially with regard to plants. We analyzed genetic variation in 22 populations of Siberian apricot (Prunus sibirica L.), a temperate deciduous species distributed in this arid belt, using two chloroplast DNA (cpDNA) sequences, seven chloroplast microsatellite loci (cpSSRs), and 31 nuclear microsatellite loci (nSSRs), to study its phylogeography. Chloroplast data showed the complete fixation of two different genetic groups: the eastern and western groups. Genetic differentiation between the two groups was significant (F_ST = 0.90925, p < 0.01). This pronounced phylogeographic break was also indicated by nSSR data, but there were disparities regarding individual populations. An asymmetric gene flow via pollen and seeds likely resulted in discordance between the present-day geography of nuclear and chloroplast lineages. There was a distinct boundary between the two large groups, which were fixed for two of the most ancestral chlorotypes. Two populations with the highest chloroplast genetic diversity were located in the Yanshan Mountains and Jinzhou, considered to be the glacial refugia. The split of chloroplasts between the eastern and western groups was estimated to have occurred ~0.1795 Ma, whereas nuclear divergence occurred approximately 13,260 years ago. Linear regression analysis showed that climatic factors (annual precipitation and annual mean temperature) had a significant correlation with mean ancestry value (P < 0.05) indicated that they were potential factors for the formation of the two groups. In addition, this boundary was a contact zone between two groups from different refugia.

Phylogeography of the endangered orchid Dendrobium moniliforme in East Asia inferred from chloroplast DNA sequences

The aim of the current study was to elucidate the phylogeographic history of Dendrobium moniliforme, an endangered orchid species, based on two chloroplast DNA (cpDNA) markers (trnC-petN and trnE-trnT). One hundred and thirty-five samples were collected from 18 natural populations of D. moniliforme covering the entire range of the Sino-Japanese Floristic Region (SJFR) of East Asia. A total of 35 distinct cpDNA haplotypes were identified in these populations, of which 23 haplotypes were each present in only one sample and thus restricted to a single population. The significantly larger N_ST value (0.586) than G_ST (0.328) (p < 0.05) demonstrated the presence of strong phylogeographic structure. Phylogenetic analyses indicated that all haplotypes were clustered into two lineages. The genetic diversity of D. moniliforme was high at the species level, reflected in its haplotype diversity (H_d=0.8862), nucleotide diversity (P_i=0.00361), total genetic diversity (H_T=0.9011), and significant differentiation (Φ_ST=0.5482). Based on mismatch distribution analysis and neutrality tests, population expansion was evident in all sampled populations and also in all populations sampled in mainland China. Three refuge areas were identified, one each in southwestern China, central-southeastern China, and the CKJ (Taiwan, Japan and Korea) Islands. The results supported the hypothesis that glacial refugia were maintained on different spatial-temporal scales in the SJFR during the last glacial maximum or earlier cold periods, suggesting that Quaternary refugial isolation promoted allopatric speciation of D. moniliforme in East Asia.

Environmental (in)dependence of a hybrid zone: Insights from molecular markers and ecological niche modeling in a hybrid zone of Origanum (Lamiaceae) on the island of Crete

The role of environment and the relative significance of endogenous versus exogenous selection in shaping hybrid zones have been crucial issues in the studies of hybridization. Recent advances in ecological niche modeling (ENM) offer new methodological tools, especially in combination with the genotyping of individuals in the hybrid zone. Here, we study the hybrid zone between the widely known spices Origanum onites and Origanum vulgare ssp. hirtum in Crete. We analyze the genetic structure of both parental taxa and their hybrid Origanum × intercendens using AFLP markers on 15 sympatric and 12 allopatric populations and employ ecological niche modeling and niche similarity tests to study their niche patterns. We complement these analyses with seed viability measurements. Our study revealed that the hybridizing taxa O. onites and O. vulgare ssp. hirtum and the resulting genotypic classes showed geographical and environmental niche similarities based on the predictions of ENMs and the subsequent similarity tests. The occurrence of the hybrid zone is not directly dependent on environmental factors which favor the fitness of the hybrid compared to the parental taxa, but rather on aspects such as historical factors and management practices, which may contribute to the localization and maintenance of the contact zone between parental species. Our results suggest that if a minimum required niche differentiation between genotypic classes is not achieved, environmental dependence might not have a prominent role on the outcome of the hybridization.

Phylogeography of postglacial range expansion in Juglans mandshurica (Juglandaceae) reveals no evidence of bottleneck, loss of genetic diversity, or isolation by distance in the leading-edge populations

The past studies of postglacial recolonization patterns in high latitude regions have revealed a significant role of dispersal capacity in shaping the genetic diversity and population structure of temperate trees. However, most of these studies have focused on species with long-distance dispersal followed by exponential population growth and were therefore unable to reveal the patterns in the case of a gradual expansion. Here we studied the impacts of postglacial range expansions on the distribution of genetic diversity in the Manchurian walnut (Juglans mandshurica), a common tree of East Asian cool-temperate deciduous forests that apparently lacks long-distance seed dispersal ability. The genetic diversity and structure of 19 natural walnut populations in Northeast China and the Korean Peninsula were examined using 17 nuclear simple sequence repeat (SSR) loci. Potential habitats under current and past climatic conditions were predicted using the ecological niche modelling (ENM) method. Bayesian clustering analysis revealed three groups, which were inferred to have diverged through multiple glacial-interglacial cycles in multiple refugia during the Quaternary Period. ENM estimated a southward range shift at the LGM, but high suitability scores still occurred in the western parts of the Changbai Mountains (Northeast China), the Korean peninsula and the exposed seafloor of the Yellow Sea. In contrast to most other cool-temperate trees co-occurring in the same region, the Manchurian walnut did not show any evidence of a population bottleneck, loss of genetic diversity or isolation by distance during the postglacial expansion. Our study clearly indicates that current northern populations originated from one glacial lineage and recolonization via a gradually advancing front due to the lack of a long-distance seed dispersal mechanism led to no latitudinal decrease in genetic diversity.

Contrasting genetic patterns between two coexisting Eleutherococcus species in northern China

Climate oscillations are the key factors to understand the patterns in modern biodiversity. East Asia harbors the most diverse temperate flora, largely because an extensive terrestrial ice cap was absent during repeated Pleistocene glaciation-interglacial cycles. Comparing the demographic histories of species that are codistributed and are close relatives may provide insight into how the process of climate change influences species ranges. In this study, we compared the spatial genetic structure and demographic histories of two coexisting Eleutherococcus species, Eleutherococcus senticosus and E. sessiliflorus. Both species are distributed in northern China, regions that are generally considered to be sensitive to climatic fluctuations. These regions once hosted temperate forest, but this temperate forest was replaced by tundra and taiga forest during the Last Glacial Maximum (LGM), according to pollen records. Using three chloroplast DNA fragments, we assessed the genetic structure of 20 and 9 natural populations of E. senticosus and E. sessiliflorus, respectively. Extremely contrasting genetic patterns were found between the two species; E. sessiliflorus had little genetic variation, whereas E. senticosus had considerably higher levels of genetic variation (15 haplotypes). We speculated that a recent severe bottleneck may have resulted in the extremely low genetic diversity in E. sessiliflorus. In E. senticosus, populations in Northeast China (NEC) harbored all of the haplotypes found in this species and included private haplotypes. The populations in NEC had higher levels of genetic diversity than did those from North China (NC). Therefore, we suggest that both the NC and NEC regions can sustain LGM refugia and that lineage admixture from multiple refugia took place after the LGM elevated the local genetic diversity in NEC. In NEC, multiple genetic hot spots were found in the Changbai Mountains and the Xiaoxing'an Range, which implied that multiple locations in NEC may sustain LGM refugia, even in the Xiaoxing'an Range.

Genetic Structure and Evolutionary History of Three Alpine Sclerophyllous Oaks in East Himalaya-Hengduan Mountains and Adjacent Regions

The East Himalaya-Hengduan Mountains (EH-HM) region has a high biodiversity and harbors numerous endemic alpine plants. This is probably the result of combined orographic and climate oscillations occurring since late Tertiary. Here, we determined the genetic structure and evolutionary history of alpine oak species (including Quercus spinosa, Quercus aquifolioides, and Quercus rehderiana) using both cytoplasmic-nuclear markers and ecological niche models (ENMs), and elucidated the impacts of climate oscillations and environmental heterogeneity on their population demography. Our results indicate there were mixed genetic structure and asymmetric contemporary gene flow within them. The ENMs revealed a similar demographic history for the three species expanded their ranges from the last interglacial (LIG) to the last glacial maximum (LGM), which was consistent with effective population sizes changes. Effects of genetic drift and fragmentation of habitats were responsible for the high differentiation and the lack of phylogeographic structure. Our results support that geological and climatic factors since Miocene triggered the differentiation, evolutionary origin and range shifts of the three oak species in the studied area and also emphasize that a multidisciplinary approach combining molecular markers, ENMs and population genetics can yield deep insights into diversification and evolutionary dynamics of species.

Genetic Structure and Evolutionary History of Three Alpine Sclerophyllous Oaks in East Himalaya-Hengduan Mountains and Adjacent Regions

The East Himalaya-Hengduan Mountains (EH-HM) region has a high biodiversity and harbors numerous endemic alpine plants. This is probably the result of combined orographic and climate oscillations occurring since late Tertiary. Here, we determined the genetic structure and evolutionary history of alpine oak species (including Quercus spinosa, Quercus aquifolioides, and Quercus rehderiana) using both cytoplasmic-nuclear markers and ecological niche models (ENMs), and elucidated the impacts of climate oscillations and environmental heterogeneity on their population demography. Our results indicate there were mixed genetic structure and asymmetric contemporary gene flow within them. The ENMs revealed a similar demographic history for the three species expanded their ranges from the last interglacial (LIG) to the last glacial maximum (LGM), which was consistent with effective population sizes changes. Effects of genetic drift and fragmentation of habitats were responsible for the high differentiation and the lack of phylogeographic structure. Our results support that geological and climatic factors since Miocene triggered the differentiation, evolutionary origin and range shifts of the three oak species in the studied area and also emphasize that a multidisciplinary approach combining molecular markers, ENMs and population genetics can yield deep insights into diversification and evolutionary dynamics of species.

Contributions of multiple refugia during the last glacial period to current mainland populations of Korean pine (Pinus koraiensis)

The northern microrefugia that existed during the Last Glacial Maximum (LGM) are a key factor in the demographic history of species. Pinus koraiensis has a unique distribution in northeast Asia. The Changbai Mountains and the Korean peninsula (CM/KP) are usually considered to be the LGM refugia for P. koraiensis. However, the Xiaoxingan Range (XR), at the northern part of this species' distribution, is another possible refugium. We used chloroplast sequencing and ten nuclear single-copy gene loci to calculate the genetic diversity pattern of P. koraiensis. The probabilities of a single LGM refugium and of multiple LGM refugia were calculated based on approximate Bayesian computation. The effect of the latitudinal gradient on genetic diversity was not significant. However, unique alleles occurred at low frequencies in CM/KP and XR. A conservative estimate of the coalescence time between CM/KP and XR is 0.4 million years ago, a time prior to the LGM. Gene flow between CM/KP and XR was estimated to be more than one in per generation, an amount that may be sufficient to limit genetic divergence between the regions. Our study strongly supports the hypothesis that XR was another LGM refugium in addition to CM/KP.

Influence of niche similarity on hybridization between Myriophyllum sibiricum and M. spicatum

The impact of ecological factors on natural hybridization is of widespread interest. Here, we asked whether climate niche influences hybridization between the two closely related plant species Myriophyllum sibiricum and M. spicatum. Eight microsatellite loci and two chloroplast fragments were used to investigate the occurrence of hybridization between these two species in two co-occurring regions: north-east China (NEC) and the Qinghai-Tibetan Plateau (QTP). The climate niches of the species were quantified by principal component analysis with bioclimatic data, and niche comparisons were performed between the two species in each region. Reciprocal hybridization was observed, and M. sibiricum was favoured as the maternal species. Furthermore, hybrids were rare in NEC but common in the QTP. Accordingly, in NEC, the two species were climatically distinct, and hybrids only occurred in the narrow geographical or ecological transition zone, whereas in the QTP, obvious niche overlaps were found for the two species, and hybrids occurred in multiple contact zones. This association between hybridization pattern and climate niche similarity suggests that the level of hybridization was promoted by niche overlap. Compared with the parental species, similar climate niches were found for the hybrid populations in the QTP, indicating that other environmental factors rather than climate were important for hybrid persistence. Our findings highlight the significance of climate niche with respect to hybridization patterns in plants.

Living on the edge: the role of geography and environment in structuring genetic variation in the southernmost populations of a tropical oak

Understanding the factors determining genetic diversity and structure in peripheral populations is a long-standing goal of evolutionary biogeography, yet little empirical information is available for tropical species. In this study, we combine information from nuclear microsatellite markers and niche modelling to analyse the factors structuring genetic variation across the southernmost populations of the tropical oak Quercus segoviensis. First, we tested the hypothesis that genetic variability decreases with population isolation and increases with local habitat suitability and stability since the Last Glacial Maximum (LGM). Second, we employed a recently developed multiple matrix regression with randomisation (MMRR) approach to study the factors associated with genetic divergence among the studied populations and test the relative contribution of environmental and geographic isolation to contemporary patterns of genetic differentiation. We found that genetic diversity was negatively correlated with average genetic differentiation with other populations, indicating that isolation and limited gene flow have contributed to erode genetic variability in some populations. Considering the relatively small size of the study area (<120 km), analyses of genetic structure indicate a remarkable inter-population genetic differentiation. Environmental dissimilarity and differences in current and past climate niche suitability and their additive effects were not associated with genetic differentiation after controlling for geographic distance, indicating that local climate does not contribute to explain spatial patterns of genetic structure. Overall, our data indicate that geographic isolation, but not current or past climate, is the main factor determining contemporary patterns of genetic diversity and structure within the southernmost peripheral populations of this tropical oak.

Variation in the frequency and extent of hybridization between Leucosceptrum japonicum and L. stellipilum (Lamiaceae) in the Central Japanese Mainland

Variations in the frequency and extent of hybridization among mixed populations located in the same contact zone provide natural laboratories for the study of extrinsic reproductive isolation maintaining species integrity. In this study, we examined the pattern of hybridization between L. japonicum and L. stellipilum among mixed populations in different localities of a contact zone. The genetic structures from three sympatric populations and six mixed populations in the hybrid zone, and five reference populations far from the contact zone, were characterized using 10 neutral nuclear microsatellite markers. Evidence from genetic distance-based clustering analysis, the frequency distribution of admixture proportion values, and the hybrid category assignment approaches indicated that the frequency and extent of hybridization varied considerably among populations in the contact zone between L. japonicum and L. stellipilum. One likely explanation is that variation in exogenous (ecological) selection among populations might contribute to differences in frequency and extent of hybridization. The present study will facilitate future research exploring the evolution of reproductive isolation between L. japonicum and L. stellipilum.

Evolutionary history of a widespread tree species Acer mono in East Asia

East Asia has the most diverse temperate flora in the world primarily due to the lack of Pleistocene glaciation and the geographic heterogeneity. Although increasing phylogeography studies in this region provided more proofs in this issue, discrepancies and uncertainty still exist, especially in northern temperate deciduous broad-leaved and coniferous mixed forest region (II). And a widespread plant species could reduce the complexity to infer the relationship between diversity and physiographical pattern. Hence, we studied the evolution history of a widespread temperate tree, Acer mono, populations in region II and the influence of physiographic patterns on intraspecific genetic diversity. Analyses of chloroplast sequences and nuclear microsatellites indicated high levels of genetic diversity. The diversity distribution was spatially heterogeneous and a latitudinal cline existed in both markers. The spatial distribution pattern between genetic diversity within A. mono and the diversity at species level was generally consistent. Western subtropical evergreen broad-leaved forest subregion (IVb) had a unique ancient chloroplast clade (CP3) and a nuclear gene pool (GP5) with dominance indicating the critical role of this area in species diversification. Genetic data and ecological niche model results both suggested that populations in region II disappeared during the last glacial maximum (LGM) and recovered from south of Changbai Mt. and the Korean Peninsula. Two distribution centers were likely during the LGM, one in the north edge of warm temperate deciduous broad-leaved forest region (III) and another in the south edge of region III. This was reflected by the genetic pattern with two spatially independent genetic groups. This study highlights the key role of region III in sustaining genetic diversity in the northern range and connecting diversity between southern and northern range. We elucidated the diversity relationship between vegetation regions which could facilitate the understanding of biodiversity origin and maintenance in East Asia.

Bioclimatic, ecological, and phenotypic intermediacy and high genetic admixture in a natural hybrid of octoploid strawberries

PREMISE OF THE STUDY

Hybrid zones provide "natural laboratories" for understanding the processes of selection, reinforcement, and speciation. We sought to gain insight into the degree of introgression and the extent of ecological-phenotypic intermediacy in the natural hybrid strawberry, Fragaria × ananassa subsp. cuneifolia. •

METHODS

We used whole-plastome sequencing to identify parental species-specific (Fragaria chiloensis and F. virginiana) chloroplast single-nucleotide polymorphisms and combined the use of these with nuclear microsatellite markers to genetically characterize the hybrid zone. We assessed the potential role of selection in the observed geographic patterns by bioclimatically characterizing the niche of the hybrid populations and phenotypically characterizing hybrid individuals of known genomic constitution. •

KEY RESULTS

Significant admixture and little overall maternal bias in chloroplast or nuclear genomes suggest a high degree of interfertility among the parental and hybrid species and point to a long history of backcrossing and genetic mixing in the hybrid zone. Even though hybrids were phenotypically intermediate to the parental species, there was a discernible fingerprint of the parental genotype within hybrid individuals. Thus, although the pattern of introgression observed suggests geographic limitations to gene flow, it may be reinforced by selection for specific parental traits in the bioclimatically intermediate habitat occupied by the hybrid. •

CONCLUSIONS

This work uncovered the genetic complexity underlying the hybrid zone of the wild relatives of the cultivated strawberry. It lays the foundation for experimental dissection of the causes of genomic introgression and nuclear-cytoplasmic disassociation, and for understanding other parts of Fragaria evolutionary history.

Reinforcement in plants

A major goal of evolutionary biology is to understand how diverging populations become species. The evolution of reproductive isolation (RI) halts the genomic homogenization caused by gene flow and recombination, and enables differentiation and local adaptations to become fixed between newly forming species. Selection can favor the strengthening of RI through a process termed reinforcement. Reinforcement occurs when selection favors traits that decrease mating between two incipient species in response to costly mating or the production of maladapted hybrids. Although this process has been investigated more frequently in animals, there is also evidence of reinforcement in plants. There are three strategies for the investigation of the process of reinforcement: case studies of species or diverging taxa; experimental evolution studies; and comparative studies. Here, I discuss how all three strategies find evidence consistent with reinforcement occurring in plants. I focus largely on case studies, and use research on Phlox drummondii to illustrate the importance of testing alternative hypotheses. Although the existing evidence suggests that reinforcement can occur, further investigations, particularly using large-scale comparative studies, are needed to determine the importance of reinforcement in plant speciation.