The role of hybridization during ecological divergence of southwestern white pine ( Pinus strobiformis ) and limber pine ( P. flexilis )

Genomic insights into hybrid zone formation: The role of climate, landscape, and demography in the emergence of a novel hybrid lineage

Population demographic changes, alongside landscape, geographic and climate heterogeneity, can influence the timing, stability and extent of introgression where species hybridise. Thus, quantifying interactions across diverged lineages, and the relative contributions of interspecific genetic exchange and selection to divergence at the genome-wide level is needed to better understand the drivers of hybrid zone formation and maintenance. We used seven latitudinally arrayed transects to quantify the contributions of climate, geography and landscape features to broad patterns of genetic structure across the hybrid zone of Populus trichocarpa and P. balsamifera and evaluated the demographic context of hybridisation over time. We found genetic structure differed among the seven transects. While ancestry was structured by climate, landscape features influenced gene flow dynamics. Demographic models indicated a secondary contact event may have influenced contemporary hybrid zone formation with the origin of a putative hybrid lineage that inhabits regions with higher aridity than either of the ancestral groups. Phylogenetic relationships based on chloroplast genomes support the origin of this hybrid lineage inferred from demographic models based on the nuclear data. Our results point towards the importance of climate and landscape patterns in structuring the contact zones between P. trichocarpa and P. balsamifera and emphasise the value whole genome sequencing can have to advancing our understanding of how neutral processes influence divergence across space and time.

The role of introgressive hybridization in shaping the geographically isolated gene pools of wax palm populations (genus Ceroxylon)

The speciation continuum is the process by which genetic groups diverge until they reach reproductive isolation. It has become common in the literature to show that this process is gradual and flickering, with possibly many instances of secondary contact and introgression after divergence has started. The level of divergence might vary among genomic regions due to, among others, the different forces and roles of selection played by the shared regions. Through hybrid capture, we sequenced ca. 4,000 nuclear regions in populations of six species of wax palms, five of which form a monophyletic group (genus Ceroxylon, Arecaceae: Ceroxyloideae). We show that in this group, the different populations show varying degrees of introgressive hybridization, and two of them are backcrosses of the other three 'pure' species. This is particularly interesting because these three species are dioecious, have a shared main pollinator, and have slightly overlapping reproductive seasons but highly divergent morphologies. Our work supports shows wax palms diverge under positive and background selection in allopatry, and hybridize due to secondary contact and inefficient reproductive barriers, which sustain genetic diversity. Introgressed regions are generally not under positive selection. Peripheral populations are backcrosses of other species; thus, introgressive hybridization is likely modulated by demographic effects rather than selective pressures. In general, these species might function as an 'evolutionary syngameon' where expanding, peripheral, small, and isolated populations maintain diversity by crossing with available individuals of other wax palms. In the Andean context, species can benefit from gained variation from a second taxon or the enhancement of population sizes by recreating a common genetic pool.

Improving species delimitation for effective conservation: a case study in the endemic maple-leaf oak (Quercus acerifolia)

Species delimitation is challenging in lineages that exhibit both high plasticity and introgression. This challenge can be compounded by collection biases, which may downweight specimens morphologically intermediate between traditional species. Additionally, mismatch between named species and observable phenotypes can compromise species conservation. We studied the species boundaries of Quercus acerifolia, a tree endemic to Arkansas, U.S. We performed morphometric analyses of leaves and acorns from 527 field and 138 herbarium samples of Q. acerifolia and its close relatives, Q. shumardii and Q. rubra. We employed two novel approaches: sampling ex situ collections to detect phenotypic plasticity caused by environmental variation and comparing random field samples with historical herbarium samples to identify collection biases that might undermine species delimitation. To provide genetic evidence, we also performed molecular analyses on genome-wide SNPs. Quercus acerifolia shows distinctive morphological, ecological, and genomic characteristics, rejecting the hypothesis that Q. acerifolia is a phenotypic variant of Q. shumardii. We found mismatches between traditional taxonomy and phenotypic clusters. We detected underrepresentation of morphological intermediates in herbarium collections, which may bias species discovery and recognition. Rare species conservation requires considering and addressing taxonomic problems related to phenotypic plasticity, mismatch between taxonomy and morphological clusters, and collection biases.

Genomic clines across the species boundary between a hybrid pine and its progenitor in the eastern Tibetan Plateau

Most species have clearly defined distribution ranges and ecological niches. The genetic and ecological causes of species differentiation and the mechanisms that maintain species boundaries between newly evolved taxa and their progenitors are, however, less clearly defined. This study investigated the genetic structure and clines in Pinus densata, a pine of hybrid origin on the southeastern Tibetan Plateau, to gain an understanding of the contemporary dynamics of species barriers. We analyzed genetic diversity in a range-wide collection of P. densata and representative populations of its progenitors, Pinus tabuliformis and Pinus yunnanensis, using exome capture sequencing. We detected four distinct genetic groups within P. densata that reflect its migration history and major gene-flow barriers across the landscape. The demographies of these genetic groups in the Pleistocene were associated with regional glaciation histories. Interestingly, population sizes rebounded rapidly during interglacial periods, suggesting persistence and resilience of the species during the Quaternary ice age. In the contact zone between P. densata and P. yunnanensis, 3.36% of the analyzed loci (57 849) showed exceptional patterns of introgression, suggesting their potential roles in either adaptive introgression or reproductive isolation. These outliers showed strong clines along critical climate gradients and enrichment in a number of biological processes relevant to high-altitude adaptation. This indicates that ecological selection played an important role in generating genomic heterogeneity and a genetic barrier across a zone of species transition. Our study highlights the forces that operate to maintain species boundaries and promote speciation in the Qinghai-Tibetan Plateau and other mountain systems.

Sequential hybridization may have facilitated ecological transitions in the Southwestern pinyon pine syngameon

Multispecies interbreeding networks, or syngameons, have been increasingly reported in natural systems. However, the formation, structure, and maintenance of syngameons have received little attention. Through gene flow, syngameons can increase genetic diversity, facilitate the colonization of new environments, and contribute to hybrid speciation. In this study, we evaluated the history, patterns, and consequences of hybridization in a pinyon pine syngameon using morphological and genomic data to assess genetic structure, demographic history, and geographic and climatic data to determine niche differentiation. We demonstrated that Pinus edulis, a dominant species in the Southwestern US and a barometer of climate change, is a core participant in the syngameon, involved in the formation of two drought-adapted hybrid lineages including the parapatric and taxonomically controversial fallax-type. We found that species remain morphologically and genetically distinct at range cores, maintaining species boundaries while undergoing extensive gene flow in areas of sympatry at range peripheries. Our study shows that sequential hybridization may have caused relatively rapid speciation and facilitated the colonization of different niches, resulting in the rapid formation of two new lineages. Participation in the syngameon may allow adaptive traits to be introgressed across species barriers and provide the changes needed to survive future climate scenarios.

A wide hybrid zone mediated by precipitation contributed to confused geographical structure of Scutiger boulengeri

Confused geographical structure of a population and mitonuclear discordance are shaped by a combination of rapid changes in population demographics and shifts in ecology. In this study, we generated a time-calibrated phylogeny of Scutiger boulengeri, an endemic Xizang alpine toad occurring in mountain streams on the Qinghai-Xizang (Tibet) Plateau (QTP). Based on three mitochondrial DNA (mtDNA) genes, eight clades were assigned to three deeply divergent lineages. Analysis of nuclear DNA (nuDNA) genes revealed three distinct clusters without geographic structure, indicating significantly high rates of gene flow. Coalescent theory framework analysis (approximate Bayesian computation model DIYABC and Migrate-N) suggested that divergence of the main intraspecific clusters was the result of hybridization after secondary contact in the Holocene around 0.59 million years ago (Ma). The ratio of mtDNA F _ST (fixation index) to nuDNA F _ST was 2.3, thus failing to show male-biased dispersal. Geographic cline analysis showed that a wide hybrid zone was initially established in southwestern China, without significant reproductive isolation but with strong introgression in S. boulengeri, suggesting high hybrid fitness. Furthermore, mtDNA genes exhibited isolation by distance (IBD) while nuDNA genes exhibited significant isolation by environment (IBE). Results suggested that mitonuclear discordance may have initially been caused by geographic isolation, followed by precipitation-mediated hybridization, producing a wide hybrid zone and geographic structure confusion of nuDNA genes in S. boulengeri. This study indicated that complicated historical processes may have led to specific genetic patterns, with a specific climate factor facilitating gene flow in the system.

Western Spruce Budworm Effects on Forest Resilience

Western spruce budworm (Choristoneura freemani Razowski) is the most destructive defoliator of forests in the western US. Forests in northern New Mexico experienced high levels of WSBW-caused defoliation and subsequent mortality between the 1980s and 2010s. The effects of severe western spruce budworm outbreaks on stand dynamics in the US Southwest are still relatively unknown, but understanding the impacts is important to the management and resilience of these forests. To begin addressing this knowledge gap, we conducted a study along two gradients: an elevational gradient from mixed-conifer to spruce-fir forests and a gradient of WSBW-caused defoliation intensity. We recorded overstory and understory stand conditions (size structure, species composition, damaging agents). Western spruce budworm was the primary damaging agent of host trees in all stands andcaused host tree mortality across all size classes, particularly in spruce-fir stands. Results indicate an unsustainable level of mortality in spruce-fir stands and a transition towards non-host species in mixed-conifer stands. Low levels of regeneration coupled with high overstory mortality rates indicate a potential lack of resilience in spruce-fir stands, whereas resilience to future western spruce budworm defoliation events may have increased in mixed-conifer stands affected by these outbreaks.

Reduced representation sequencing to understand the evolutionary history of Torrey pine (Pinus torreyana Parry) with implications for rare species conservation

Understanding the contribution of neutral and adaptive evolutionary processes to population differentiation is often necessary for better informed management and conservation of rare species. In this study, we focused on Pinus torreyana Parry (Torrey pine), one of the world's rarest pines, endemic to one island and one mainland population in California. Small population size, low genetic diversity, and susceptibility to abiotic and biotic stresses suggest Torrey pine may benefit from inter-population genetic rescue to preserve the species' evolutionary potential. We leveraged reduced representation sequencing to tease apart the respective contributions of stochastic and deterministic evolutionary processes to population differentiation. We applied these data to model spatial and temporal demographic changes in effective population sizes and genetic connectivity, to identify loci possibly under selection, and evaluate genetic rescue as a potential conservation strategy. Overall, we observed exceedingly low standing variation within both Torrey pine populations, reflecting consistently low effective population sizes across time, and limited genetic differentiation, suggesting maintenance of gene flow between populations following divergence. However, genome scans identified more than 2000 candidate SNPs potentially under divergent selection. Combined with previous observations indicating population phenotypic differentiation, this indicates natural selection has likely contributed to the evolution of population genetic differences. Thus, while reduced genetic diversity, small effective population size, and genetic connectivity between populations suggest genetic rescue could mitigate the adverse effects of rarity, evidence for adaptive differentiation suggests genetic mixing could disrupt adaptation. Further work evaluating the fitness consequences of inter-population admixture is necessary to empirically evaluate the trade-offs associated with genetic rescue in Torrey pine.

Comparative Association Mapping Reveals Conservation of Major Gene Resistance to White Pine Blister Rust in Southwestern White Pine (Pinus strobiformis) and Limber Pine (P. flexilis)

All native North American white pines are highly susceptible to white pine blister rust (WPBR) caused by Cronartium ribicola. Understanding genomic diversity and molecular mechanisms underlying genetic resistance to WPBR remains one of the great challenges in improvement of white pines. To compare major gene resistance (MGR) present in two species, southwestern white pine (Pinus strobiformis) Cr3 and limber pine (P. flexilis) Cr4, we performed association analyses of Cr3-controlled resistant traits using single nucleotide polymorphism (SNP) assays designed with Cr4-linked polymorphic genes. We found that ∼70% of P. flexilis SNPs were transferable to P. strobiformis. Furthermore, several Cr4-linked SNPs were significantly associated with the Cr3-controlled traits in P. strobiformis families. The most significantly associated SNP (M326511_1126R) almost colocalized with Cr4 on the Pinus consensus linkage group 8, suggesting that Cr3 and Cr4 might be the same R locus, or have localizations very close to each other in the syntenic region of the P. strobiformis and P. flexilis genomes. M326511_1126R was identified as a nonsynonymous SNP, causing amino acid change (Val₃₇₆Ile) in a putative pectin acetylesterase, with coding sequences identical between the two species. Moreover, top Cr3-associated SNPs were further developed as TaqMan genotyping assays, suggesting their usefulness as marker-assisted selection (MAS) tools to distinguish genotypes between quantitative resistance and MGR. This work demonstrates the successful transferability of SNP markers between two closely related white pine species in the hybrid zone, and the possibility for deployment of MAS tools to facilitate long-term WPBR management in P. strobiformis breeding and conservation.

The evolutionary assembly of forest communities along environmental gradients: recent diversification or sorting of pre-adapted clades?

Recent studies have demonstrated that ecological processes that shape community structure and dynamics change along environmental gradients. However, much less is known about how the emergence of the gradients themselves shape the evolution of species that underlie community assembly. In this study, we address how the creation of novel environments leads to community assembly via two nonmutually exclusive processes: immigration and ecological sorting of pre-adapted clades (ISPC), and recent adaptive diversification (RAD). We study these processes in the context of the elevational gradient created by the uplift of the Central Andes. We develop a novel approach and method based on the decomposition of species turnover into within- and among-clade components, where clades correspond to lineages that originated before mountain uplift. Effects of ISPC and RAD can be inferred from how components of turnover change with elevation. We test our approach using data from over 500 Andean forest plots. We found that species turnover between communities at different elevations is dominated by the replacement of clades that originated before the uplift of the Central Andes. Our results suggest that immigration and sorting of clades pre-adapted to montane habitats is the primary mechanism shaping tree communities across elevations.

Multigenerational backcrossing and introgression between two woodrat species at an abrupt ecological transition

When organisms experience secondary contact after allopatric divergence, genomic regions can introgress differentially depending on their relationships with adaptation, reproductive isolation, recombination, and drift. Analyses of genome-wide patterns of divergence and introgression could provide insight into the outcomes of hybridization and the potential relationship between allopatric divergence and reproductive isolation. Here, we generate population genetic data (26,262 SNPs; 353 individuals) using a reduced-representation sequencing approach to quantify patterns of ancestry, differentiation, and introgression between a pair of ecologically distinct mammals-the desert woodrat (N. lepida) and Bryant's woodrat (N. bryanti)-that hybridize at a sharp ecotone in southern California. Individual ancestry estimates confirmed that hybrids were rare in this bimodal hybrid zone, and entirely consisted of a few F₁ individuals and a broad range of multigenerational backcrosses. Genomic cline analyses indicated more than half of loci had elevated introgression from one genomic background into the other. However, introgression was not associated with relative or absolute measures of divergence, and loci with extreme values for both were not typically found near detoxification enzymes previously implicated in dietary specialization for woodrats. The decoupling of differentiation and introgression suggests that processes other than adaptation, such as drift, may underlie the extreme clines at this contact zone.

Adaptive evolution in a conifer hybrid zone is driven by a mosaic of recently introgressed and background genetic variants

Extant conifer species may be susceptible to rapid environmental change owing to their long generation times, but could also be resilient due to high levels of standing genetic diversity. Hybridisation between closely related species can increase genetic diversity and generate novel allelic combinations capable of fuelling adaptive evolution. Our study unravelled the genetic architecture of adaptive evolution in a conifer hybrid zone formed between Pinus strobiformis and P. flexilis. Using a multifaceted approach emphasising the spatial and environmental patterns of linkage disequilibrium and ancestry enrichment, we identified recently introgressed and background genetic variants to be driving adaptive evolution along different environmental gradients. Specifically, recently introgressed variants from P. flexilis were favoured along freeze-related environmental gradients, while background variants were favoured along water availability-related gradients. We posit that such mosaics of allelic variants within conifer hybrid zones will confer upon them greater resilience to ongoing and future environmental change and can be a key resource for conservation efforts.

Field Translocation of Mountain Pine Beetles Suggests Phoretic Mite Communities Are Locally Adapted, and Mite Populations Respond Variably to Climate Warming

Simple Summary

Climate warming has significant effects on forest insect populations, particularly bark beetles, which cause millions of hectares of forest tree damage. Bark beetles live alongside a diverse host of other organisms which affect the success of beetle attacks on trees and are also affected by climate changes. Here, we explore climate effects on symbiotic mite communities associated with the mountain pine beetle (Dendroctonus ponderosae). We show that warming causes significant shifts in the abundance of mites. These effects were dependent on source population, suggesting mite populations are adapted to their local climates. Understanding beetle–mite patterns is important because mites can directly affect beetle reproduction by feeding on eggs, or indirectly affect beetle health by introducing fungi. Our results provide foundational information for understanding how climate change will affect beetle–mite associations; and serve to help determine how these shifting associations will affect the success of bark beetles in forest ecosystems.

Abstract

Temperature is a key determining factor in the population dynamics of forest insects and their associated biota. Bark beetles, often considered key agents of change in forest ecosystems, are particularly affected by warming in their environment. Beetles associate with various phoretic mite species that have direct/indirect effects on beetle fitness and population dynamics, although there is limited knowledge of how temperature affects these communities. Here, we use a field reciprocal translocation experiment with the addition of a novel “warming” environment to represent future changes in local environment in two populations of a keystone bark beetle species (Dendroctonus ponderosae). We hypothesize that mite community abundances as carried by bark beetles are significantly altered when not in their native environments and when subjected to climate warming. We use multivariate generalized linear models based on species abundance data to show that mite community compositions significantly differ across different field climates; and that these patterns diverge between source populations, indicating local adaptation. Our study offers foundational information on the general effects of simulated climate-warming on the compositional shifts of common and abundant biotic associates of mountain pine beetles and may be used as a model system for other important insect–mite systems.

Effective double-digest RAD sequencing and genotyping despite large genome size

Obtaining informative data is the ambition of any genomic project, but in nonmodel species with very large genomes, pursuing such a goal requires surmounting a series of analytical challenges. Double-digest RAD sequencing is routinely used in nonmodel organisms and offers some control over the volume of data obtained. However, the volume of data recovered is not always an indication of the reliability of data sets, and quality checks are necessary to ensure that true and artefactual information is set apart. In the present study, we aim to fill the gap existing between the known applicability of RAD sequencing methods in plants with large genomes and the use of the retrieved loci for population genetic inference. By analysing two populations of Cypripedium calceolus, a nonmodel orchid species with a large genome size (1C ~ 31.6 Gbp), we provide a complete workflow from library preparation to bioinformatic filtering and inference of genetic diversity and differentiation. We show how filtering strategies to dismiss potentially misleading data need to be explored and adapted to data set-specific features. Moreover, we suggest that the occurrence of organellar sequences in libraries should not be neglected when planning the experiment and analysing the results. Finally, we explain how, in the absence of prior information about the genome of the species, seeking high standards of quality during library preparation and sequencing can provide an insurance against unpredicted technical or biological constraints.

RAD-sequencing reveals patterns of diversification and hybridization, and the accumulation of reproductive isolation in a clade of partially sympatric, tropical island trees

A common pattern observed in temperate tree clades is that species are often morphologically distinct and partially interfertile but maintain species cohesion despite ongoing hybridization where ranges overlap. Although closely related species commonly occur in sympatry in tropical ecosystems, little is known about patterns of hybridization within a clade over time, and the implications of this hybridization for the maintenance of species boundaries. In this study, we focused on a clade of sympatric trees in the genus Diospyros in the Mascarene islands and investigated whether species are genetically distinct, whether they hybridize, and how patterns of hybridization are related to the time since divergence among species. We sampled multiple populations from each of 12 Mascarene Diospyros species, generated genome-wide single nucleotide polymorphism data using 2bRADseq, and conducted population genomic and phylogenomic analyses. We found that Mascarene Diospyros species diverged millions of years ago and are today largely genetically distinct from one another. Although hybridization was observed between closely related species belonging to the same subclade, more distantly related species showed little evidence of interspecific hybridization. Phylogenomic analyses also suggested that introgression has occurred during the evolutionary history of the clade. This suggests that, as diversification progressed, interspecific hybridization occurred among species, but became infrequent as lineages diverged from one another and evolved reproductive barriers. Species now coexist in partial sympatry, and experience limited hybridization between close relatives. Additional research is needed to better understand the role that introgression may have played in adaptation and diversification of Mascarene Diospyros, and its relevance for conservation.

Morphological Differences in Pinus strobiformis Across Latitudinal and Elevational Gradients

The phenotype of trees is determined by the relationships and interactions among genetic and environmental influences. Understanding the patterns and processes that are responsible for phenotypic variation is facilitated by studying the relationships between phenotype and the environment among many individuals across broad ecological and climatic gradients. We used Pinus strobiformis, which has a wide latitudinal distribution, as a model species to: (a) estimate the relative importance of different environmental factors in predicting these morphological traits and (b) characterize the spatial patterns of standing phenotypic variation of cone and seed traits across the species' range. A large portion of the total variation in morphological characteristics was explained by ecological, climatic and geographical variables (54.7% collectively). The three climate, vegetation and geographical variable groups, each had similar total ability to explain morphological variation (43.4%, 43.8%, 51.5%, respectively), while the topographical variable group had somewhat lower total explanatory power (36.9%). The largest component of explained variance (33.6%) was the four-way interaction of all variable sets, suggesting that there is strong covariation in environmental, climate and geographical variables in their relationship to morphological traits of southwest white pine across its range. The regression results showed that populations in more humid and warmer climates expressed greater cone length and seed size. This may in part facilitate populations of P. strobiformis in warmer and wetter portions of its range growing in dense, shady forest stands, because larger seeds provide greater resources to germinants at the time of germination. Our models provide accurate predictions of morphological traits and important insights regarding the factors that contribute to their expression. Our results indicate that managers should be conservative during reforestation efforts to ensure match between ecotypic variation in seed source populations. However, we also note that given projected large range shift due to climate change, managers will have to balance the match between current ecotypic variation and expected range shift and changes in local adaptive optima under future climate conditions.

Patterns of hybridization and cryptic introgression among one- and four-needled pinyon pines

BACKGROUND AND AIMS

Pinyon pine hybridization is widely acknowledged, but the frequency of and contributors to such interspecific mating remain largely unstudied. Pinus quadrifolia has three to four needles per fascicle, suggesting that it is a result of hybridization between the five-needled P. juarezensis and the single-needled P. monophylla. In this study we address the taxonomic validity of P. juarezensis, the hybrid origin of P. quadrifolia and the presence of hybridization and intermediate morphology as a result of interspecific hybridization in this complex.

METHODS

We address these questions by combining a genomic and morphological approach. We generated 1868 single nucleotide polymorphisms (SNPs) to detect genetic clusters using principal co-ordinates analyis, discriminant analysis of principal components, fastSTRUCTURE and ADMIXTURE analyses, and performed a morphological analysis of the leaves.

KEY RESULTS

We found that the five-needled pinyons did not differ genetically from the four-needled P. quadrifolia, reducing the status of P. juarezensis to P. quadrifolia. We also found no evidence that P. quadrifolia is of hybrid origin from P. juarezensis × P. monophylla but is instead a genetically distinct species with natural needle number variation that has yet to be explained. Hybridization does occur in this complex, but mostly between P. quadrifolia and P. californiarum, and less commonly between P. quadrifolia and P. monophylla. Interestingly, some hybrid derivatives were detected between both single-needled taxa, P. monophylla and P. californiarum, a hybrid combination that has not yet been proposed. Hybrids have intermediate morphology when they have similar genetic contributions from both parental species; however, when one parent contributes more, hybrid derivatives resemble the parent with higher genetic contribution, resulting in cryptic introgression.

CONCLUSIONS

Our detailed sampling across the distribution of this complex allows us to describe the patterns of hybridization among these taxa, resolves an ancient taxonomic conflict and provides insights into the challenges of exclusively using morphological traits when identifying these taxa with cryptic hybridization and variable morphology.

Genotyping-by-sequencing and ecological niche modeling illuminate phylogeography, admixture, and Pleistocene range dynamics in quaking aspen (Populus tremuloides)

Populus tremuloides is the widest-ranging tree species in North America and an ecologically important component of mesic forest ecosystems displaced by the Pleistocene glaciations. Using phylogeographic analyses of genome-wide SNPs (34,796 SNPs, 183 individuals) and ecological niche modeling, we inferred population structure, ploidy levels, admixture, and Pleistocene range dynamics of P. tremuloides, and tested several historical biogeographical hypotheses. We found three genetic lineages located mainly in coastal-Cascades (cluster 1), east-slope Cascades-Sierra Nevadas-Northern Rockies (cluster 2), and U.S. Rocky Mountains through southern Canadian (cluster 3) regions of the P. tremuloides range, with tree graph relationships of the form ((cluster 1, cluster 2), cluster 3). Populations consisted mainly of diploids (86%) but also small numbers of triploids (12%) and tetraploids (1%), and ploidy did not adversely affect our genetic inferences. The main vector of admixture was from cluster 3 into cluster 2, with the admixture zone trending northwest through the Rocky Mountains along a recognized phenotypic cline (Utah to Idaho). Clusters 1 and 2 provided strong support for the "stable-edge hypothesis" that unglaciated southwestern populations persisted in situ since the last glaciation. By contrast, despite a lack of clinal genetic variation, cluster 3 exhibited "trailing-edge" dynamics from niche suitability predictions signifying complete northward postglacial expansion. Results were also consistent with the "inland dispersal hypothesis" predicting postglacial assembly of Pacific Northwestern forest ecosystems, but rejected the hypothesis that Pacific-coastal populations were colonized during outburst flooding from glacial Lake Missoula. Overall, congruent patterns between our phylogeographic and ecological niche modeling results and fossil pollen data demonstrate complex mixtures of stable-edge, refugial locations, and postglacial expansion within P. tremuloides. These findings confirm and refine previous genetic studies, while strongly supporting a distinct Pacific-coastal genetic lineage of quaking aspen.

Tall Pinus luzmariae trees with genes from P. herrerae

CONTEXT

Pinus herrerae and P. luzmariae are endemic to western Mexico, where they cover an area of more than 1 million hectares. Pinus herrerae is also cultivated in field trials in South Africa and South America, because of its considerable economic importance as a source of timber and resin. Seed quality, afforestation success and desirable traits may all be influenced by the presence of hybrid trees in seed stands.

AIMS

We aimed to determine the degree of hybridization between P. herrerae and P. luzmariae in seed stands of each species located in the Sierra Madre Occidental, Durango, Mexico.

METHODS

AFLP molecular markers from samples of 171 trees across five populations were analyzed with STRUCTURE and NewHybrids software to determine the degree of introgressive hybridization. The accuracy of STRUCTURE and NewHybrids in detecting hybrids was quantified using the software Hybridlab 1.0. Morphological analysis of 131 samples from two populations of P. herrerae and two populations of P. luzmariae was also conducted by Random Forest classification. The data were compared by Principal Coordinate Analysis (PCoA) in GenAlex 6.501.

RESULTS

Hybridization between Pinus herrerae and P. luzmariae was observed in all seed stands under study and resulted in enhancement of desirable silvicultural traits in the latter species. In P. luzmariae, only about 16% molecularly detected hybrids correspond to those identified on a morphological basis. However, the morphology of P. herrerae is not consistent with the molecularly identified hybrids from one population and is only consistent with 3.3 of those from the other population.

CONCLUSIONS

This is the first report of hybrid vigour (heterosis) in Mexican pines. Information about hybridization and introgression is essential for developing effective future breeding programs, successful establishment of plantations and management of natural forest stands. Understanding how natural hybridization may influence the evolution and adaptation of pines to climate change is a cornerstone to sustainable forest management including adaptive silviculture.

Tracing the footprints of a moving hybrid zone under a demographic history of speciation with gene flow

A lack of optimal gene combinations, as well as low levels of genetic diversity, is often associated with the formation of species range margins. Conservation efforts rely on predictive modelling using abiotic variables and assessments of genetic diversity to determine target species and populations for controlled breeding, germplasm conservation and assisted migration. Biotic factors such as interspecific competition and hybridization, however, are largely ignored, despite their prevalence across diverse taxa and their role as key evolutionary forces. Hybridization between species with well-developed barriers to reproductive isolation often results in the production of offspring with lower fitness. Generation of novel allelic combinations through hybridization, however, can also generate positive fitness consequences. Despite this possibility, hybridization-mediated introgression is often considered a threat to biodiversity as it can blur species boundaries. The contribution of hybridization towards increasing genetic diversity of populations at range margins has only recently gathered attention in conservation studies. We assessed the extent to which hybridization contributes towards range dynamics by tracking spatio-temporal changes in the central location of a hybrid zone between two recently diverged species of pines: Pinus strobiformis and P. flexilis. By comparing geographic cline centre estimates for global admixture coefficient with morphological traits associated with reproductive output, we demonstrate a northward shift in the hybrid zone. Using a combination of spatially explicit, individual-based simulations and linkage disequilibrium variance partitioning, we note a significant contribution of adaptive introgression towards this northward movement, despite the potential for differences in regional population size to aid hybrid zone movement. Overall, our study demonstrates that hybridization between recently diverged species can increase genetic diversity and generate novel allelic combinations. These novel combinations may allow range margin populations to track favourable climatic conditions or facilitate adaptive evolution to ongoing and future climate change.

Multiomics analysis reveals a distinct mechanism of oleaginousness in the emerging model alga Chromochloris zofingiensis

Chromochloris zofingiensis, featured due to its capability to simultaneously synthesize triacylglycerol (TAG) and astaxanthin, is emerging as a leading candidate alga for production uses. To better understand the oleaginous mechanism of this alga, we conducted a multiomics analysis by systematically integrating time-resolved transcriptomes, lipidomes and metabolomes in response to nitrogen deprivation. The data analysis unraveled the distinct mechanism of TAG accumulation, which involved coordinated stimulation of multiple biological processes including supply of energy and reductants, carbon reallocation from protein and starch, and 'pushing' and 'pulling' carbon to TAG synthesis. Unlike the model alga Chlamydomonas, de novo fatty acid synthesis in C. zofingiensis was promoted, together with enhanced turnover of both glycolipids and phospholipids, supporting the drastic need of acyls for TAG assembly. Moreover, genomewide analysis identified many key functional enzymes and transcription factors that had engineering potential for TAG modulation. Two genes encoding glycerol-3-phosphate acyltransferase (GPAT), the first committed enzyme for TAG assembly, were found in the C. zofingiensis genome; in vivo functional characterization revealed that extrachloroplastic GPAT instead of chloroplastic GPAT played a central role in TAG synthesis. These findings illuminate distinct oleaginousness mechanisms in C. zofingiensis and pave the way towards rational manipulation of this alga to becone an emerging model for trait improvements.

Holocene Population Decline and Conservation Implication for the Western Hercules Beetle, Dynastes grantii (Coleoptera, Scarabaeidae)

The Western Hercules beetle (Dynastes grantii) is endemic to the highland forest habitats of southwestern United States and northern Mexico. The habitats harbor many endemic species, but are being threatened by rapid climate change and urban development. In this study, the genetic structure of D. grantii populations from southwestern United States was investigated. Specifically, genomic data from double-digest restriction-site-associated DNA sequencing libraries were utilized to test whether geographically distant populations from the Mogollon Rim (Arizona [N = 12 individuals] and New Mexico [N = 10 individuals]) are genetically structured. The study also estimated the effective population size of the Mogollon Rim populations based on genetic diversity. The results indicated that the 2 geographic populations from the Mogollon Rim were not genetically structured. A population size reduction was detected since the end of the last glacial period, which coincided with a reduction of forest habitat in the study area. The results implied that the connectivity and the size of highland forest habitats in the Mogollon Rim could have been the major factors shaping the population genetic structure and demographic history of D. grantii. The Western Hercules beetle could be a useful flagship species for local natural history education and to promote the conservation of highland forest habitats.

Combining mitochondrial and nuclear genome analyses to dissect the effects of colonization, environment, and geography on population structure in Pinus tabuliformis

The phylogeographic histories of plants in East Asia are complex and shaped by both past large-scale climatic oscillations and dramatic tectonic events. The impact of these historic events, as well as ecological adaptation, on the distribution of biodiversity remains to be elucidated. Pinus tabuliformis is the dominant coniferous tree in northern China, with a large distribution across wide environmental gradients. We examined genetic variation in this species using genotyping-by-sequencing and mitochondrial (mt) DNA markers. We found population structure on both nuclear and mt genomes with a geographic pattern that corresponds well with the landscape of northern China. To understand the contributions of environment, geography, and colonization history to the observed population structure, we performed ecological niche modeling and partitioned the among-population genomic variance into isolation by environment (IBE), isolation by distance (IBD), and isolation by colonization (IBC). We used mtDNA, which is transmitted by seeds in pine, to reflect colonization. We found little impact of IBE, IBD, and IBC on variation in neutral SNPs, but significant impact of IBE on a group of outlier loci. The lack of IBC illustrates that the maternal history can be quickly eroded from the nuclear genome by high rates of gene flow. Our results suggest that genomic variation in P. tabuliformis is largely affected by neutral and stochastic processes, and the signature of local adaptation is visible only at robust outlier loci. This study enriches our understanding on the complex evolutionary forces that shape the distribution of genetic variation in plant taxa in northern China, and guides breeding, conservation, and reforestation programs for P. tabuliformis.