Genetic divergence, range expansion and possible homoploid hybrid speciation among pine species in Northeast China

Genetic differentiation across a steep and narrow environmental gradient: Quantitative genetic and genomic insights into Lake Superior populations of Quercus rubra

Adaptive differentiation of traits and underlying loci can occur at a small geographical scale if natural selection is stronger than countervailing gene flow and drift. We investigated this hypothesis using coupled quantitative genetic and genomic approaches for a wind-pollinated tree species, Quercus rubra, along the steep, narrow gradient of the Lake Superior coast that encompasses four USDA Hardiness Zones within 100 km. For the quantitative genetic component of this study, we examined phenotypic differentiation among eight populations in a common garden, measuring seed mass, germination, height, stem diameter, leaf number, specific leaf area and survival. For the genomic component, we quantified genetic differentiation for 26 populations from the same region using RAD-seq. Because hybridisation with Quercus ellipsoidalis occurs in other parts of the species' range, we included two populations of this congener for comparison. In the common garden study, we found a strong signal of population differentiation that was significantly associated with at least one climate factor for nine of 10 measured traits. In contrast, we found no evidence of genomic differentiation among populations based on FST or any other measures. However, both distance-based and genotype-environment association analyses identified loci showing the signature of selection, with one locus in common across five analyses. This locus was associated with the minimum temperature of the coldest month, a factor that defines the climate zones and was also significant in the common garden analyses. In addition, we documented introgression from Q. ellipsoidalis into Q. rubra, with rates of introgression correlated with the climate gradient. In sum, this study reveals signatures of selection at the quantitative trait and genomic level consistent with climate adaptation, a pattern that is more often documented at a much broader geographical scale, especially in long-lived wind-pollinated species.

Evolutionary radiation of the Eurasian Pinus species under pervasive gene flow

Evolutionary radiation, a pivotal aspect of macroevolution, offers valuable insights into evolutionary processes. The genus Pinus is the largest genus in conifers withc . $$ c. $$ 90% of the extant species emerged in the Miocene, which signifies a case of rapid diversification. Despite this remarkable history, our understanding of the mechanisms driving radiation within this expansive genus has remained limited. Using exome capture sequencing and a fossil-calibrated phylogeny, we investigated the divergence history, niche diversification, and introgression among 13 closely related Eurasian species spanning climate zones from the tropics to the boreal Arctic. We detected complex introgression among lineages in subsection Pinus at all stages of the phylogeny. Despite this widespread gene exchange, each species maintained its genetic identity and showed clear niche differentiation. Demographic analysis unveiled distinct population histories among these species, which further influenced the nucleotide diversity and efficacy of purifying and positive selection in each species. Our findings suggest that radiation in the Eurasian pines was likely fueled by interspecific recombination and further reinforced by their adaptation to distinct environments. Our study highlights the constraints and opportunities for evolutionary change, and the expectations of future adaptation in response to environmental changes in different lineages.

The Role of Hybridisation in the Making of the Species-Rich Arctic-Alpine Genus Saxifraga (Saxifragaceae)

Evolutionary processes fuelling rapid species diversification are not yet fully understood, although their major contribution to overall patterns of plant biodiversity is well established. Hybridisation is among the least understood of these processes, despite its multifaceted role in speciation processes being widely accepted. Species of the large arctic-alpine genus Saxifraga are notorious for their ability to hybridise; however, the overall role of hybridisation and polyploidisation for the diversification of this genus remains unknown. Here, we provide a comprehensive genus-wide review of hybridisation accounts and ploidy levels. We find that the sections of Saxifraga vary greatly in their propensity to hybridise. The majority of natural hybridisation accounts are from recent localised events (n = 71). Hybridisation hotspots were located in the Pyrenees and the European Alps, thus contrasting with the overall distribution of species richness in the genus. Hybrids or hybrid populations are often short-lived in Saxifraga due to a multitude of reproductive barriers, most commonly low F1 hybrid fertility. However, these barriers are not always fully effective, allowing for backcrossing and the formation of hybrid swarms. In addition, we find that the incidence of polyploidy varies widely across different sections of Saxifraga, with species-rich sections Porphyrion and Saxifraga showing divergent polyploidy proportions. Overall, we show that hybridisation and polyploidisation played differential roles in the diversification of this large genus. Nevertheless, a significant proportion of species are yet to be scrutinised, particularly among the Asian Saxifraga species, illustrating the need for systematic further study to fully unravel the role of hybridisation during the evolution of Saxifraga.

Hybridization and introgression during density-dependent range expansion: European wildcats as a case study

Interbreeding between historically allopatric species with incomplete reproductive barriers may result when species expand their range. The genetic consequences of such hybridization depend critically on the dynamics of the range expansion. Hybridization models during range expansion have been developed but assume dispersal to be independent from neighboring population densities. However, organisms may disperse because they are attracted by conspecifics or because they prefer depopulated areas. Here, through spatially explicit simulations, we assess the effect of various density-dependent dispersal modes on the introgression between two species. We find huge introgression from the local species into the invasive one with all dispersal modes investigated, even when the hybridization rate is relatively low. This represents a general expectation for neutral genes even if the dispersal modes differ in colonization times and amount of introgression. Invasive individuals attracted by conspecifics need more time to colonize the whole area and are more introgressed by local genes, whereas the opposite is found for solitary individuals. We applied our approach to a recent expansion of European wildcats in the Jura Mountains and the hybridization with domestic cats. We show that the simulations explained better the observed level of introgression at nuclear, mtDNA, and Y chromosome markers, when using solitary dispersal for wildcats instead of random or gregarious dispersal, in accordance with ecological knowledge. Thus, use of density-dependent dispersal models increases the predictive power of the approach.

Genetic analysis of admixture and hybrid patterns of Populus hopeiensis and P. tomentosa

Hybridization and introgression have resulted in reticulate evolution within the genus Populus. Consequently, the origin and evolutionary history of some hybrids has become blurred. P. hopeiensis and P. tomentosa are endemic to China, and there is still controversy about their origin. We employ phylogeny, Bayesian estimation of admixture, and approximate Bayesian computation to investigate their origin with 10 nuclear DNA and 6 cpDNA regions. The combined evidences firmly support the hypothesis that they are hybrids and dominated by F1s. P. hopeiensis was generated via hybridization between the paternal species P. alba and maternal species P. davidiana. Surprisingly, P. tomentosa was divided into two genetic types with different maternal parents. P. adenopoda hybridized with P. alba directly to generate the first genetic type (mb1) and hybridized with P. davidiana followed by P. alba to generate the second (mb2). In both genetic types, P. alba acted as the male parent. The maternal parent was P. adenopoda and P. davidiana for mb1 and mb2, respectively. Hybridization not only generated these hybrids but also resulted in a unidirectional gene flow from P. davidiana to P. adenopoda. The Populus species have maintained a delicate balance between their genetic integrity and gene exchange.

Effects of Geological and Environmental Events on the Diversity and Genetic Divergence of Four Closely Related Pines: Pinus koraiensis, P. armandii, P. griffithii, and P. pumila

The effects of mountain uplift and environmental oscillations on nucleotide variability and species divergence remain largely unknown in East Asia. In this study, based on multiple nuclear DNA markers, we investigated the levels and patterns of nucleotide diversity and interspecific divergence in four closely related pines in China, i.e., Pinus koraiensis, P. armandii, P. griffithii, and P. pumila. The four pine taxa shared low levels of nucleotide polymorphisms at the species level. P. pumila had the highest silent nucleotide diversity (πsil = 0.00661) whereas P. griffithii had the lowest (πsil = 0.00175), while the levels of genetic polymorphism in P. armandii (πsil = 0.00508) and P. koraiensis (πsil = 0.00652) were intermediate between the other two species. Population genetic structure analysis showed that variations primarily existed within populations of the four pine species, presumably due to habitat fragmentation or the island-like distributions of Pinus species. Population divergence (FST) analysis showed that the genetic divergence between P. griffithii and P. koraiensis was much greater than that between P. koraiensis and the other two pines species. Isolation-with-migration analysis suggested that asymmetric gene flow had occurred between any two pairs of pine species. Phylogenetic analyses indicated that the four allied species split into two groups about 1.37 million years ago, where P. armandii and P. pumila were closer and clustered as sister species, whereas P. koraiensis and P. griffithii were clustered on another branch. Our results and those obtained in previous studies suggest that mountain uplift and geological climate oscillations may have led to the patterns of genetic divergence and nucleotide variations in these four pine species.

Hybridization as a facilitator of species range expansion

Explaining the evolution of species geographical ranges is fundamental to understanding how biodiversity is distributed and maintained. The solution to this classic problem in ecology and evolution remains elusive: we still do not fully know how species geographical ranges evolve and what factors fuel range expansions. Resolving this problem is now more crucial than ever with increasing biodiversity loss, global change and movement of species by humans. Here, we describe and evaluate the hypothesis that hybridization between species can contribute to species range expansion. We discuss how such a process can occur and the empirical data that are needed to test this hypothesis. We also examine how species can expand into new environments via hybridization with a resident species, and yet remain distinct species. Generally, hybridization may play an underappreciated role in influencing the evolution of species ranges. Whether-and to what extent-hybridization has such an effect requires further study across more diverse taxa.

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.

Importance of incomplete lineage sorting and introgression in the origin of shared genetic variation between two closely related pines with overlapping distributions

Genetic variation shared between closely related species may be due to retention of ancestral polymorphisms because of incomplete lineage sorting (ILS) and/or introgression following secondary contact. It is challenging to distinguish ILS and introgression because they generate similar patterns of shared genetic diversity, but this is nonetheless essential for inferring accurately the history of species with overlapping distributions. To address this issue, we sequenced 33 independent intron loci across the genome of two closely related pine species (Pinus massoniana Lamb. and Pinus hwangshanensis Hisa) from Southeast China. Population structure analyses revealed that the species showed slightly more admixture in parapatric populations than in allopatric populations. Levels of interspecific differentiation were lower in parapatry than in allopatry. Approximate Bayesian computation suggested that the most likely speciation scenario explaining this pattern was a long period of isolation followed by a secondary contact. Ecological niche modeling suggested that a gradual range expansion of P. hwangshanensis during the Pleistocene climatic oscillations could have been the cause of the overlap. Our study therefore suggests that secondary introgression, rather than ILS, explains most of the shared nuclear genomic variation between these two species and demonstrates the complementarity of population genetics and ecological niche modeling in understanding gene flow history. Finally, we discuss the importance of contrasting results from markers with different dynamics of migration, namely nuclear, chloroplast and mitochondrial DNA.

Transient hybridization, not homoploid hybrid speciation, between ancient and deeply divergent conifers

PREMISE OF THE STUDY

Homoploid hybrid speciation is receiving growing attention due the increasing recognition of its role in speciation. We investigate if individuals intermediate in morphology between the two species of the conifer genus Athrotaxis represent a homoploid hybrid species, A. laxifolia, or are spontaneous F1 hybrids.

METHODS

A total of 1055 individuals of Athrotaxis cupressoides and A. selaginoides, morphologically intermediate individuals, and two putative hybrid swarms were sampled across the range of the genus and genotyped with 13 microsatellites. We used simulations to test the power of our data to identify the pure species, F1s, F2s, and backcross generations.

KEY RESULTS

We found that Athrotaxis cupressoides and A. selaginoides are likely the most divergent congeneric conifers known, but the intermediates are F1 hybrids, sharing one allele each from A. cupressoides and A. selaginoides at six loci with completely species specific alleles. The hybrid swarms contain wide genetic variation with stronger affinities to the locally dominant species, A. selaginoides and A. selaginoides backcrosses outnumbering A. cupressoides backcrosses. In addition, we observed evidence for isolated advanced generation backcrosses within the range of the pure species.

CONCLUSIONS

We conclude that, even though they can be large and long-lived, Athrotaxis hybrid swarms are on a trajectory of decline and will eventually be reabsorbed by the parental species. However, this process may take millennia and fossil evidence suggests that such events have occurred repeatedly since the early Quaternary. Given this timeline, our study highlights the many obstacles to homoploid hybrid speciation.

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.

Phylogeny and biogeography of Primula sect. Armerina: implications for plant evolution under climate change and the uplift of the Qinghai-Tibet Plateau

BACKGROUND

The historical orogenesis and associated climatic changes of mountain areas have been suggested to partly account for the occurrence of high levels of biodiversity and endemism. However, their effects on dispersal, differentiation and evolution of many groups of plants are still unknown. In this study, we examined the detailed diversification history of Primula sect. Armerina, and used biogeographic analysis and macro-evolutionary modeling to investigate a series of different questions concerning the evolution of the geographical and ecological distribution of the species in this section.

RESULTS

We sequenced five chloroplast and one nuclear genes for species of Primula sect. Armerina. Neither chloroplast nor nuclear trees support the monophyly of the section. The major incongruences between the two trees occur among closely related species and may be explained by hybridization. Our dating analyses based on the chloroplast dataset suggest that this section began to diverge from its relatives around 3.55 million years ago, largely coinciding with the last major uplift of the Qinghai-Tibet Plateau (QTP). Biogeographic analysis supports the origin of the section in the Himalayan Mountains and dispersal from the Himalayas to Northeastern QTP, Western QTP and Hengduan Mountains. Furthermore, evolutionary models of ecological niches show that the two P. fasciculata clades have significantly different climatic niche optima and rates of niche evolution, indicating niche evolution under climatic changes and further providing evidence for explaining their biogeographic patterns.

CONCLUSION

Our results support the hypothesis that geologic and climatic events play important roles in driving biological diversification of organisms in the QTP area. The Pliocene uplift of the QTP and following climatic changes most likely promoted both the inter- and intraspecific divergence of Primula sect. Armerina. This study also illustrates how niche evolution under climatic changes influences biogeographic patterns.

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.

Phylogeographical analysis of a cold-temperate freshwater fish, the Amur sleeper (Perccottus glenii) in the Amur and Liaohe River basins of Northeast Asia

The Amur sleeper Perccottus glenii (Perciformes, Gobioidei, Odontobutidae) is well known as an invasive fish in the river basins of Eastern and Central Europe, but its genetic background is unavailable across its native habitats in northeast Asia. In this study, we used the mitochondrial cytochrome b gene by sampling 19 populations of P. glenii across its native distributional areas of Liaohe and Amur River basins to explore its evolutionary history. Phylogenetic analyses identified three major clades within P. glenii, among which Clade A and Clade B were co-distributed in the Liaohe and Amur River basins, and Clade C was restricted to the latter. Molecular dating showed that the splits of Clades A, B and C have happened in the late Early-early Middle Pleistocene and the most recent common ancestors of these clades have been presented in the late Middle-early Late Pleistocene. The P. glenii showed very high levels of genetic structure among populations (ΦST = 0.801), probably due to the characters of its life histories with very limited dispersal ability. The admixture of different clades in some populations of P. glenii probably reflects historical secondary contact. These findings indicate that Pleistocene climatic oscillation and river capture were major determinants for genetic variations and evolutionary history of the P. glenii.

Climatic adaptation and ecological divergence between two closely related pine species in Southeast China

Climate is one of the most important drivers for adaptive evolution in forest trees. Climatic selection contributes greatly to local adaptation and intraspecific differentiation, but this kind of selection could also have promoted interspecific divergence through ecological speciation. To test this hypothesis, we examined intra- and interspecific genetic variation at 25 climate-related candidate genes and 12 reference loci in two closely related pine species, Pinus massoniana Lamb. and Pinus hwangshanensis Hisa, using population genetic and landscape genetic approaches. These two species occur in Southeast China but have contrasting ecological preferences in terms of several environmental variables, notably altitude, although hybrids form where their distributions overlap. One or more robust tests detected signals of recent and/or ancient selection at two-thirds (17) of the 25 candidate genes, at varying evolutionary timescales, but only three of the 12 reference loci. The signals of recent selection were species specific, but signals of ancient selection were mostly shared by the two species likely because of the shared evolutionary history. FST outlier analysis identified six SNPs in five climate-related candidate genes under divergent selection between the two species. In addition, a total of 24 candidate SNPs representing nine candidate genes showed significant correlation with altitudinal divergence in the two species based on the covariance matrix of population history derived from reference SNPs. Genetic differentiation between these two species was higher at the candidate genes than at the reference loci. Moreover, analysis using the isolation-with-migration model indicated that gene flow between the species has been more restricted for climate-related candidate genes than the reference loci, in both directions. Taken together, our results suggest that species-specific and divergent climatic selection at the candidate genes might have counteracted interspecific gene flow and played a key role in the ecological divergence of these two closely related pine species.

Homoploid hybrid origin of Yucca gloriosa: intersectional hybrid speciation in Yucca (Agavoideae, Asparagaceae)

There is a growing appreciation for the importance of hybrid speciation in angiosperm evolution. Here, we show that Yucca gloriosa (Asparagaceae: Agavoideae) is the product of intersectional hybridization between Y. aloifolia and Y. filamentosa. These species, all named by Carl Linnaeus, exist in sympatry along the southeastern Atlantic coast of the United States. Yucca gloriosa was found to share a chloroplast haplotype with Y. aloifolia in all populations sampled. In contrast, nuclear gene-based microsatellite markers in Y. gloriosa are shared with both parents. The hybrid origin of Y. gloriosa is supported by multilocus analyses of the nuclear microsatellite markers including principal coordinates analysis (PCO), maximum-likelihood hybrid index scoring (HINDEX), and Bayesian cluster analysis (STRUCTURE). The putative parental species share only one allele at a single locus, suggesting there is little to no introgressive gene flow occurring between these species and Y. gloriosa. At the same time, diagnostic markers are segregating in Y. gloriosa populations. Lack of variation in the chloroplast of Y. aloifolia, the putative maternal parent, makes it difficult to rule out multiple hybrid origins of Y. gloriosa, but allelic variation at nuclear loci can be explained by a single hybrid origin of Y. gloriosa. Overall, these data provide strong support for the homoploid hybrid origin of Y. gloriosa.

Genetic architecture of maize kernel composition in the nested association mapping and inbred association panels

The maize (Zea mays) kernel plays a critical role in feeding humans and livestock around the world and in a wide array of industrial applications. An understanding of the regulation of kernel starch, protein, and oil is needed in order to manipulate composition to meet future needs. We conducted joint-linkage quantitative trait locus mapping and genome-wide association studies (GWAS) for kernel starch, protein, and oil in the maize nested association mapping population, composed of 25 recombinant inbred line families derived from diverse inbred lines. Joint-linkage mapping revealed that the genetic architecture of kernel composition traits is controlled by 21-26 quantitative trait loci. Numerous GWAS associations were detected, including several oil and starch associations in acyl-CoA:diacylglycerol acyltransferase1-2, a gene that regulates oil composition and quantity. Results from nested association mapping were verified in a 282 inbred association panel using both GWAS and candidate gene association approaches. We identified many beneficial alleles that will be useful for improving kernel starch, protein, and oil content.