Characterizing the physical and genetic structure of the lodgepole pine × jack pine hybrid zone: mosaic structure and differential introgression

Recurrent hybridization and gene flow shaped Norway and Siberian spruce evolutionary history over multiple glacial cycles

Most tree species underwent cycles of contraction and expansion during the Quaternary. These cycles led to an ancient and complex genetic structure that has since been affected by extensive gene flow and by strong local adaptation. The extent to which hybridization played a role in this multi-layered genetic structure is important to be investigated. To study the effect of hybridization on the joint population genetic structure of two dominant species of the Eurasian boreal forest, Picea abies and P. obovata, we used targeted resequencing and obtained around 480 K nuclear SNPs and 87 chloroplast SNPs in 542 individuals sampled across most of their distribution ranges. Despite extensive gene flow and a clear pattern of Isolation-by-Distance, distinct genetic clusters emerged, indicating the presence of barriers and corridors to migration. Two cryptic refugia located in the large hybrid zone between the two species played a critical role in shaping their current distributions. The two species repeatedly hybridized during the Pleistocene and the direction of introgression depended on latitude. Our study suggests that hybridization helped both species to overcome main shifts in their distribution ranges during glacial cycles and highlights the importance of considering whole species complex instead of separate entities to retrieve complex demographic histories.

Pest defences under weak selection exert a limited influence on the evolution of height growth and drought avoidance in marginal pine populations

While droughts, intensified by climate change, have been affecting forests worldwide, pest epidemics are a major source of uncertainty for assessing drought impacts on forest trees. Thus far, little information has documented the adaptability and evolvability of traits related to drought and pests simultaneously. We conducted common-garden experiments to investigate how several phenotypic traits (i.e. height growth, drought avoidance based on water-use efficiency inferred from δ¹³C and pest resistance based on defence traits) interact in five mature lodgepole pine populations established in four progeny trials in western Canada. The relevance of interpopulation variation in climate sensitivity highlighted that seed-source warm populations had greater adaptive capability than cold populations. In test sites, warming generated taller trees with higher δ¹³C and increased the evolutionary potential of height growth and δ¹³C across populations. We found, however, no pronounced gradient in defences and their evolutionary potential along populations or test sites. Response to selection was weak in defences across test sites, but high for height growth particularly at warm test sites. Response to the selection of δ¹³C varied depending on its selective strength relative to height growth. We conclude that warming could promote the adaptability and evolvability of growth response and drought avoidance with a limited evolutionary influence from pest (biotic) pressures.

Evidence of Coevolution Between Cronartium harknessii Lineages and Their Corresponding Hosts, Lodgepole Pine and Jack Pine

Variation in rate of infection and susceptibility of Pinus spp. to the fungus Cronartium harknessii (syn. Endocronartium harknessii), the causative agent of western gall rust, has been well documented. To test the hypothesis that there is a coevolutionary relationship between C. harknessii and its hosts, we examined genetic structure and virulence of C. harknessii associated with lodgepole pine (P. contorta var. latifolia), jack pine (P. banksiana), and their hybrids. A secondary objective was to improve assessment and diagnosis of infection in hosts. Using 18 microsatellites, we assessed genetic structure of C. harknessii from 90 sites within the ranges of lodgepole pine and jack pine. We identified two lineages (East and West, F_ST = 0.677) associated with host genetic structure (r = 0.81, P = 0.001), with East comprising three sublineages. In parallel, we conducted a factorial experiment in which lodgepole pine, jack pine, and hybrid seedlings were inoculated with spores from the two primary genetic lineages. With this experiment, we refined the phenotypic categories associated with infection and demonstrated that stem width can be used as a quantitative measure of host response to infection. Overall, each host responded differentially to the fungal lineages, with jack pine exhibiting more resiliency to infection than lodgepole pine and hybrids exhibiting intermediate resiliency. Taken together, the shared genetic structure between fungus and host species, and the differential interaction of the fungal species with the hosts, supports a coevolutionary relationship between host and pathogen.[Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

A holocenic and dynamic hybrid zone between two cactophilic Drosophila species in a coastal lowland plain of the Brazilian Atlantic Forest

Hybridization and introgression are processes that contribute to shaping biological diversity. The factors promoting the formation of these processes are multiples but poorly explored in a biogeographical and ecological context. In the southeast coastal plain of the Brazilian Atlantic Forest, a hybrid zone was described between two closely related cactophilic species, Drosophila antonietae and D. serido. Here, we revisited and analysed specimens from this hybrid zone to evaluate its temporal and spatial dynamic. We examined allopatric and sympatric populations of the flies using independent sources of data such as mitochondrial and nuclear sequences, microsatellite loci, morphometrics of wings and male genitalia, and climatic niche models. We also verified the emergence of the flies from necrotic tissues of collected cacti to verify the role of host association for the population dynamics. Our results support the existence of a hybrid zone due to secondary contact and limited to the localities where the two species are currently in contact. Furthermore, we detected asymmetric bidirectional introgression and the maintenance of the species integrity, ecological association and morphological characters, suggesting selection and limited introgression. Considering our paleomodels, probably this hybrid zone is recent and the contact occurred during the Holocene to the present day, favoured by range expansion of their populations due to expansion of open and dry areas in eastern South America during palaeoclimatic and geomorphological events.

Primary and Secondary Metabolite Profiles of Lodgepole Pine Trees Change with Elevation, but Not with Latitude

Climate change has a large influence on plant functional and phenotypic traits including plant primary and secondary metabolites. One well-established approach to investigating the variation in plant metabolites involves studying plant populations along elevation and latitude gradients. We considered how two space-for-time climate change gradients (elevation and latitude) influence carbohydrate reserves (soluble sugars, starches) and secondary metabolites (monoterpenes, diterpene resin acids) of lodgepole pine trees in western Canada. We were particularly interested in the relationship of terpenes and carbohydrates with a wide range of tree, site, and climatic factors. We found that only elevation had a strong influence on the expression of both terpenes and carbohydrates of trees. Specifically, as elevation increased, concentrations of monoterpenes and diterpenes generally increased and soluble sugars (glucose, sucrose, total sugars) decreased. In contrast, latitude had no impact on either of terpenes or carbohydrates. Furthermore, we found a positive relationship between concentrations of starch and total terpenes and diterpenes in the elevation study; whereas neither starches nor sugars were correlated to terpenes in the latitude study. Similarly, both terpenes and carbohydrates had a much greater number of significant correlations to site characteristics such as slope, basal area index, and sand basal area, in the elevational than in the latitude study. Overall, these results support the conclusion that both biotic and abiotic factors likely drive the patterns of primary and secondary metabolite profiles of lodgepole pine along geographical gradients. Also, presence of a positive relationship between terpenes and starches suggests an interaction between primary ad secondary metabolites of lodgepole pine trees.

The role of the environment in the spatial dynamics of an extensive hybrid zone between two neotropical cats

Identifying factors that create and maintain a hybrid zone is of great interest to ecology, evolution and, more recently, conservation biology. Here, we investigated the role of environmental features in shaping the spatial dynamics of a hybrid zone between the southern tigrina, Leopardus guttulus, and Geoffroy's cat, L. geoffroyi, testing for exogenous selection as the main force acting on its maintenance. These Neotropical felid species are mainly allopatric, with a restricted area of sympatry in the ecotone between the Atlantic Forest and Pampa biomes. As both biomes have experienced high rates of anthropogenic habitat alteration, we also analysed the influence of habitat conversion on the hybrid zone structure. To do this, we used 13 microsatellite loci to identify potential hybrids and generated ecological niche models for them and their parental species. We compared the influence of variables on parental species and hybrid occurrence and calculated the amount of niche overlap among them. Parental species showed different habitat requirements and predicted co-occurrence was restricted to the forest-grassland mosaic of the ecotone. However, hybrids were found beyond this area, mainly in the range of L. geoffroyi. Hybrids demonstrated higher tolerance to habitat alteration than parental types, with a probability of occurrence that was positively related with mosaics of cropland areas and remnants of natural vegetation. These results indicate that exogenous selection alone does not drive the dynamics of the hybrid zone, and that habitat conversion influences its structure, potentially favouring hybrids over parental species.

Linking genotype to phenotype to identify genetic variation relating to host susceptibility in the mountain pine beetle system

Identifying genetic variants responsible for phenotypic variation under selective pressure has the potential to enable productive gains in natural resource conservation and management. Despite this potential, identifying adaptive candidate loci is not trivial, and linking genotype to phenotype is a major challenge in contemporary genetics. Many of the population genetic approaches commonly used to identify adaptive candidates will simultaneously detect false positives, particularly in nonmodel species, where experimental evidence is seldom provided for putative roles of the adaptive candidates identified by outlier approaches. In this study, we use outcomes from population genetics, phenotype association, and gene expression analyses as multiple lines of evidence to validate candidate genes. Using lodgepole and jack pine as our nonmodel study species, we analyzed 17 adaptive candidate loci together with 78 putatively neutral loci at 58 locations across Canada (N > 800) to determine whether relationships could be established between these candidate loci and phenotype related to mountain pine beetle susceptibility. We identified two candidate loci that were significant across all population genetic tests, and demonstrated significant changes in transcript abundance in trees subjected to wounding or inoculation with the mountain pine beetle fungal associate Grosmannia clavigera. Both candidates are involved in central physiological processes that are likely to be invoked in a trees response to stress. One of these two candidate loci showed a significant association with mountain pine beetle attack status in lodgepole pine. The spatial distribution of the attack-associated allele further coincides with other indicators of susceptibility in lodgepole pine. These analyses, in which population genetics was combined with laboratory and field experimental validation approaches, represent first steps toward linking genetic variation to the phenotype of mountain pine beetle susceptibility in lodgepole and jack pine, and provide a roadmap for more comprehensive analyses.

Purifying selection does not drive signatures of convergent local adaptation of lodgepole pine and interior spruce

BACKGROUND

Lodgepole pine (Pinus contorta) and interior spruce (Picea glauca, Picea engelmannii, and their hybrids) are distantly related conifer species. Previous studies identified 47 genes containing variants associated with environmental variables in both species, providing evidence of convergent local adaptation. However, if the intensity of purifying selection varies with the environment, clines in nucleotide diversity could evolve through linked (background) selection that would yield allele frequency-environment signatures resembling local adaptation. If similar geographic patterns in the strength of purifying selection occur in these species, this could result in the convergent signatures of local adaptation, especially if the landscape of recombination is conserved. In the present study, we investigated whether spatially/environmentally varying purifying selection could give rise to the convergent signatures of local adaptation that had previously reported.

RESULTS

We analyzed 86 lodgepole pine and 50 interior spruce natural populations spanning heterogeneous environments in western Canada where previous analyses had found signatures of convergent local adaptation. We estimated nucleotide diversity and Tajima's D for each gene within each population and calculated the strength of correlations between nucleotide diversity and environmental variables. Overall, these estimates in the genes with previously identified convergent local adaptation signatures had no similar pattern between pine and spruce. Clines in nucleotide diversity along environmental variables were found for interior spruce, but not for lodgepole pine. In spruce, genes with convergent adaption signatures showed a higher strength of correlations than genes without convergent adaption signatures, but there was no such disparity in pine, which suggests the pattern in spruce may have arisen due to a combination of selection and hybridization.

CONCLUSIONS

The results rule out purifying/background selection as a driver of convergent local adaption signatures in lodgepole pine and interior spruce.

Phytochemicals as mediators for host range expansion of a native invasive forest insect herbivore

Mountain pine beetle (MPB) has recently invaded jack pine forests in western Canada. This invasion signifies a climate change-induced range expansion by a native insect. The mechanism underlying this invasion is unknown, but likely involves phytochemicals that play critical roles in MPB biology. Thus far, studies have investigated the compatibility of jack pine chemistry with beetles and their microbial symbionts. I have identified three phytochemical mechanisms that have likely facilitated the host range expansion of MPB. First, jack pine chemistry is overall similar to that of the historical hosts of MPB. In particular, jack pine chemistry is compatible with beetle pheromone production, aggregation on host trees and larval development. Furthermore, the compatibility of jack pine chemistry maintains beneficial interactions between MPB and its microbial symbionts. Second, compared with historical hosts, the novel host not only has lower concentrations of toxic and repellent defense chemicals, but also contains large concentrations of chemicals promoting host colonization by MPB. These patterns are especially pronounced when comparing novel hosts with well-defended historical hosts. Finally, before MPBs arrived in jack pine forests, they invaded a zone of hybrids of novel and historical hosts that likely improved beetle success on jack pine, as hybrids show chemical characteristics of both hosts. In conclusion, the phytochemistry of jack pine has likely facilitated the biological invasion of this novel host by MPB.

Growth gains from selective breeding in a spruce hybrid zone do not compromise local adaptation to climate

Hybrid zones contain extensive standing genetic variation that facilitates rapid responses to selection. The Picea glauca × Picea engelmannii hybrid zone in western Canada is the focus of tree breeding programs that annually produce ~90 million reforestation seedlings. Understanding the direct and indirect effects of selective breeding on adaptive variation is necessary to implement assisted gene flow (AGF) polices in Alberta and British Columbia that match these seedlings with future climates. We decomposed relationships among hybrid ancestry, adaptive traits, and climate to understand the implications of selective breeding for climate adaptations and AGF strategies. The effects of selection on associations among hybrid index estimated from ~6,500 SNPs, adaptive traits, and provenance climates were assessed for ~2,400 common garden seedlings. Hybrid index differences between natural and selected seedlings within breeding zones were small in Alberta (average +2%), but larger and more variable in BC (average −7%, range −24% to +1%), slightly favoring P. glauca ancestry. The average height growth gain of selected seedlings over natural seedlings within breeding zones was 36% (range 12%–86%). Clines in growth with temperature‐related variables were strong, but differed little between selected and natural populations. Seedling hybrid index and growth trait associations with evapotranspiration‐related climate variables were stronger in selected than in natural seedlings, indicating possible preadaptation to drier future climates. Associations among cold hardiness, hybrid ancestry, and cold‐related climate variables dominated signals of local adaptation and were preserved in breeding populations. Strong hybrid ancestry–phenotype–climate associations suggest that AGF will be necessary to match interior spruce breeding populations with shifting future climates. The absence of antagonistic selection responses among traits and maintenance of cold adaptation in selected seedlings suggests breeding populations can be safely redeployed using AGF prescriptions similar to those of natural populations.

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.

Change in soil fungal community structure driven by a decline in ectomycorrhizal fungi following a mountain pine beetle (Dendroctonus ponderosae) outbreak

Western North American landscapes are rapidly being transformed by forest die-off caused by mountain pine beetle (Dendroctonus ponderosae), with implications for plant and soil communities. The mechanisms that drive changes in soil community structure, particularly for the highly prevalent ectomycorrhizal fungi in pine forests, are complex and intertwined. Critical to enhancing understanding will be disentangling the relative importance of host tree mortality from changes in soil chemistry following tree death. Here, we used a recent bark beetle outbreak in lodgepole pine (Pinus contorta) forests of western Canada to test whether the effects of tree mortality altered the richness and composition of belowground fungal communities, including ectomycorrhizal and saprotrophic fungi. We also determined the effects of environmental factors (i.e. soil nutrients, moisture, and phenolics) and geographical distance, both of which can influence the richness and composition of soil fungi. The richness of both groups of soil fungi declined and the overall composition was altered by beetle-induced tree mortality. Soil nutrients, soil phenolics and geographical distance influenced the community structure of soil fungi; however, the relative importance of these factors differed between ectomycorrhizal and saprotrophic fungi. The independent effects of tree mortality, soil phenolics and geographical distance influenced the community composition of ectomycorrhizal fungi, while the community composition of saprotrophic fungi was weakly but significantly correlated with the geographical distance of plots. Taken together, our results indicate that both deterministic and stochastic processes structure soil fungal communities following landscape-scale insect outbreaks and reflect the independent roles tree mortality, soil chemistry and geographical distance play in regulating the community composition of soil fungi.

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.

Fitness of natural willow hybrids in a pioneer mosaic hybrid zone

Hybrid fitness is an important parameter to predict the evolutionary consequences of a hybridization event and to characterize hybrid zones. We studied fitness parameters of F1 and later-generation hybrids between the lowland species Salix purpurea and the alpine S. helvetica that have recently emerged during colonization of an alpine glacier forefield. Fruit production (number of capsules per catkin and fruit set) did not differ between hybrids and parents, but the number of seeds per capsule of F1 hybrids was slightly lower than that of later-generation hybrids and of the parents. Germination rates and seedling growth were tested on three substrates (pH 4.5, 7.0, and 8.0). Germination rates of seeds collected from F1 hybrids were lower on acid and neutral substrates, but equal at pH 8.0 compared to all other groups, while the seeds from later-generation hybrids performed as well as the parents on all three substrates. In seedling growth, the colonizer S. purpurea performed better than all other taxa on all three substrates, while hybrids resembled the subalpine species S. helvetica. Results suggest that endogenous selection acts against F1 hybrids, but favors fitter genotypes in later-generation hybrids. Exogenous selection via soil pH appears to be weak during seedling establishment. The pioneer vegetation on the glacier forefield may offer sufficient niche space for hybrid seedlings. Owing to the relatively high fitness of the hybrids and the scattered distribution of hybrids and parental individuals on the glacier forefield, this hybrid zone can be assigned to a mosaic model, probably facilitating gene flow and introgression between the parental species. As establishment of the hybrid zone appears to be linked to a colonization process, we propose to call it a pioneer mosaic hybrid zone.

Intervarietal and intravarietal genetic structure in Douglas-fir: nuclear SSRs bring novel insights into past population demographic processes, phylogeography, and intervarietal hybridization

Douglas-fir (Pseudotsuga menziesii) is one of numerous wide-range forest tree species represented by subspecies/varieties, which hybridize in contact zones. This study examined the genetic structure of this North American conifer and its two hybridizing varieties, coastal and Rocky Mountain, at intervarietal and intravarietal level. The genetic structure was subsequently associated with the Pleistocene refugial history, postglacial migration and intervarietal hybridization/introgression. Thirty-eight populations from the USA and Canada were genotyped for 13 nuclear SSRs and analyzed with simulations and traditional population genetic structuring methods. Eight genetic clusters were identified. The coastal clusters embodied five refugial populations originating from five distinct refugia. Four coastal refugial populations, three from California and one from western Canada, diverged during the Pleistocene (56.9–40.1 ka). The three Rocky Mountain clusters reflected distinct refugial populations of three glacial refugia. For Canada, ice covered during the Last Glacial Maximum, we present the following three findings. (1) One refugial population of each variety was revealed in the north of the distribution range. Additional research including paleodata is required to support and determine whether both northern populations originated from cryptic refugia situated south or north of the ice-covered area. (2) An interplay between intravarietal gene flow of different refugial populations and intervarietal gene flow by hybridization and introgression was identified. (3) The Canadian hybrid zone displayed predominantly introgressants of the Rocky Mountain into the coastal variety. This study provides new insights into the complex Quaternary dynamics of this conifer essential for understanding its evolution (outside and inside the native range), adaptation to future climates and for forest management.

Influence of water deficit on the molecular responses of Pinus contorta × Pinus banksiana mature trees to infection by the mountain pine beetle fungal associate, Grosmannia clavigera

Conifers exhibit a number of constitutive and induced mechanisms to defend against attack by pests and pathogens such as mountain pine beetle (Dendroctonus ponderosae Hopkins) and their fungal associates. Ecological studies have demonstrated that stressed trees are more susceptible to attack by mountain pine beetle than their healthy counterparts. In this study, we tested the hypothesis that water deficit affects constitutive and induced responses of mature lodgepole pine × jack pine hybrids (Pinus contorta Dougl. ex Loud. var. latifolia Engelm. ex S. Wats. × Pinus banksiana Lamb.) to inoculation with the mountain pine beetle fungal associate Grosmannia clavigera (Robinson-Jeffrey and Davidson) Zipfel, de Beer and Wingfield. The degree of stress induced by the imposed water-deficit treatment was sufficient to reduce photosynthesis. Grosmannia clavigera-induced lesions exhibited significantly reduced dimensions in water-deficit trees relative to well-watered trees at 5 weeks after inoculation. Treatment-associated cellular-level changes in secondary phloem were also observed. Quantitative RT-PCR was used to analyze transcript abundance profiles of 18 genes belonging to four families classically associated with biotic and abiotic stress responses: aquaporins (AQPs), dehydration-responsive element binding (DREB), terpene synthases (TPSs) and chitinases (CHIs). Transcript abundance profiles of a TIP2 AQP and a TINY-like DREB decreased significantly in fungus-inoculated trees, but not in response to water deficit. One TPS, Pcb(+)-3-carene synthase, and the Class II CHIs PcbCHI2.1 and PcbCHI2.2 showed increased expression under water-deficit conditions in the absence of fungal inoculation, while another TPS, Pcb(E)-β-farnesene synthase-like, and two CHIs, PcbCHI1.1 and PcbCHI4.1, showed attenuated expression under water-deficit conditions in the presence of fungal inoculation. The effects were observed both locally and systemically. These results demonstrate that both constitutive and induced carbon- and nitrogen-based defenses are affected by water deficit, suggesting potential consequences for mountain pine beetle dynamics, particularly in novel environments.

Cross-species outlier detection reveals different evolutionary pressures between sister species

Lodgepole pine (Pinus contorta var. latifolia) and jack pine (Pinus banksiana) hybridize in western Canada, an area of recent mountain pine beetle range expansion. Given the heterogeneity of the environment, and indications of local adaptation, there are many unknowns regarding the response of these forests to future outbreaks. To better understand this we aim to identify genetic regions that have adaptive potential. We used data collected on 472 single nucleotide polymorphism (SNP) loci from 576 tree samples collected across 13 lodgepole pine-dominated sites and four jack pine-dominated sites. We looked at the relationship of genetic diversity with the environment, and we identified candidate loci using both frequency-based (arlequin and bayescan) and correlation-based (matsam and bayenv) methods. We found contrasting relationships between environmental variation and genetic diversity for the species. While we identified a number of candidate outliers (34 in lodgepole pine, 25 in jack pine, and 43 interspecific loci), we did not find any loci in common between lodgepole and jack pine. Many of the outlier loci identified were correlated with environmental variation. Using rigorous criteria we have been able to identify potential outlier SNPs. We have also found evidence of contrasting environmental adaptations between lodgepole and jack pine which could have implications for beetle spread risk.

Mitochondrial DNA capture and divergence in Pinus provide new insights into the evolution of the genus

The evolution of the mitochondrial (mt) genome is far from being fully understood. Systematic investigations into the modes of inheritance, rates and patterns of recombination, nucleotide substitution, and structural changes in the mt genome are still lacking in many groups of plants. In this study, we sequenced >11kbp mtDNA segments from multiple accessions of 36 pine species to characterize the evolutionary patterns of mtDNA in the genus Pinus. We found extremely low substitution rates and complex repetitive sequences scattered across different genome regions, as well as chimeric structures that were probably generated by multiple intergenomic recombinations. The mtDNA-based phylogeny of the genus differed from that based on chloroplast and nuclear DNA in the placement of several groups of species. Such discordances suggest a series of mtDNA capture events during past range shifts of the pine species and that both vertical and horizontal inheritance are implicated in the evolution of mtDNA in Pinus. MtDNA dating revealed that most extant lineages of the genus originated during Oligocene-Miocene radiation and subgenus Strobus diversified earlier than subgenus Pinus. Our findings illustrate a reticular evolutionary pathway for the mt genome through capture and recombination in the genus Pinus, and provide new insights into the evolution of the genus.

Exploring the impact of climate variability during the Last Glacial Maximum on the pattern of human occupation of Iberia

The Last Glacial Maximum (LGM) was a global climate event, which had significant repercussions for the spatial distribution and demographic history of prehistoric populations. In Eurasia, the LGM coincides with a potential bottleneck for modern humans and may mark the divergence date for Asian and European populations (Keinan et al., 2007). In this research, the impact of climate variability on human populations in the Iberian Peninsula during the Last Glacial Maximum (LGM) is examined with the aid of downscaled high-resolution (16 × 16 km) numerical climate experiments. Human sensitivity to short time-scale (inter-annual) climate variability during this key time period, which follows the initial modern human colonisation of Eurasia and the extinction of the Neanderthals, is tested using the spatial distribution of archaeological sites. Results indicate that anatomically modern human populations responded to small-scale spatial patterning in climate variability, specifically inter-annual variability in precipitation levels as measured by the standard precipitation index. Climate variability at less than millennial scale, therefore, is shown to be an important component of ecological risk, one that played a role in regulating the spatial behaviour of prehistoric human populations and consequently affected their social networks.

Genetic introgression and species boundary of two geographically overlapping pine species revealed by molecular markers

Gene introgression and hybrid barriers have long been a major focus of studies of geographically overlapping species. Two pine species, Pinus massoniana and P. hwangshanensis, are frequently observed growing adjacent to each other, where they overlap in a narrow hybrid zone. As a consequence, these species constitute an ideal system for studying genetic introgression and reproductive barriers between naturally hybridizing, adjacently distributed species. In this study, we sampled 270 pine trees along an elevation gradient in Anhui Province, China and analyzed these samples using EST-SSR markers. The molecular data revealed that direct gene flow between the two species was fairly low, and that the majority of gene introgression was intermediated by backcrossing. On the basis of empirical observation, the on-site distribution of pines was divided into a P. massoniana zone, a hybrid zone, and a P. hwangshanensis zone. STRUCTURE analysis revealed the existence of a distinct species boundary between the two pine species. The genetic boundary of the hybrid zone, on the other hand, was indistinct owing to intensive backcrossing with parental species. Compared with P. massoniana, P. hwangshanensis was found to backcross with the hybrids more intensively, consistent with the observation that morphological and anatomical characteristics of trees in the contact zone were biased towards P. hwangshanensis. The introgression ability of amplified alleles varied across species, with some being completely blocked from interspecific introgression. Our study has provided a living example to help explain the persistence of adjacently distributed species coexisting with their interfertile hybrids.

Fitness dynamics within a poplar hybrid zone: I. Prezygotic and postzygotic barriers impacting a native poplar hybrid stand

Hybridization and introgression are pervasive evolutionary phenomena that provide insight into the selective forces that maintain species boundaries, permit gene flow, and control the direction of evolutionary change. Poplar trees (Populus L.) are well known for their ability to form viable hybrids and maintain their distinct species boundaries despite this interspecific gene flow. We sought to quantify the hybridization dynamics and postzygotic fitness within a hybrid stand of balsam poplar (Populus balsamifera L.), eastern cottonwood (P. deltoides Marsh.), and their natural hybrids to gain insight into the barriers maintaining this stable hybrid zone. We observed asymmetrical hybrid formation with P. deltoides acting as the seed parent, but with subsequent introgression biased toward P. balsamifera. Native hybrids expressed fitness traits intermediate to the parental species and were not universally unfit. That said, native hybrid seedlings were absent from the seedling population, which may indicate additional selective pressures controlling their recruitment. It is imperative that we understand the selective forces maintaining this native hybrid zone in order to quantify the impact of exotic poplar hybrids on this native system.

How the mountain pine beetle (Dendroctonus ponderosae) breached the Canadian Rocky Mountains

The mountain pine beetle (MPB; Dendroctonus ponderosae Hopkins), a major pine forest pest native to western North America, has extended its range north and eastward during an ongoing outbreak. Determining how the MPB has expanded its range to breach putative barriers, whether physical (nonforested prairie and high elevation of the Rocky Mountains) or climatic (extreme continental climate where temperatures can be below −40 °C), may contribute to our general understanding of range changes as well as management of the current epidemic. Here, we use a panel of 1,536 single nucleotide polymorphisms (SNPs) to assess population genetic structure, connectivity, and signals of selection within this MPB range expansion. Biallelic SNPs in MPB from southwestern Canada revealed higher genetic differentiation and lower genetic connectivity than in the northern part of its range. A total of 208 unique SNPs were identified using different outlier detection tests, of which 32 returned annotations for products with putative functions in cholesterol synthesis, actin filament contraction, and membrane transport. We suggest that MPB has been able to spread beyond its previous range by adjusting its cellular and metabolic functions, with genome scale differentiation enabling populations to better withstand cooler climates and facilitate longer dispersal distances. Our study is the first to assess landscape-wide selective adaptation in an insect. We have shown that interrogation of genomic resources can identify shifts in genetic diversity and putative adaptive signals in this forest pest species.

Comparison of lodgepole and jack pine resin chemistry: implications for range expansion by the mountain pine beetle, Dendroctonus ponderosae (Coleoptera: Curculionidae)

The mountain pine beetle, Dendroctonus ponderosae, is a significant pest of lodgepole pine in British Columbia (BC), where it has recently reached an unprecedented outbreak level. Although it is native to western North America, the beetle can now be viewed as a native invasive because for the first time in recorded history it has begun to reproduce in native jack pine stands within the North American boreal forest. The ability of jack pine trees to defend themselves against mass attack and their suitability for brood success will play a major role in the success of this insect in a putatively new geographic range and host. Lodgepole and jack pine were sampled along a transect extending from the beetle’s historic range (central BC) to the newly invaded area east of the Rocky Mountains in north-central Alberta (AB) in Canada for constitutive phloem resin terpene levels. In addition, two populations of lodgepole pine (BC) and one population of jack pine (AB) were sampled for levels of induced phloem terpenes. Phloem resin terpenes were identified and quantified using gas chromatography. Significant differences were found in constitutive levels of terpenes between the two species of pine. Constitutive α-pinene levels – a precursor in the biosynthesis of components of the aggregation and antiaggregation pheromones of mountain pine beetle – were significantly higher in jack pine. However, lower constitutive levels of compounds known to be toxic to bark beetles, e.g., 3-carene, in jack pine suggests that this species could be poorly defended. Differences in wounding-induced responses for phloem accumulation of five major terpenes were found between the two populations of lodgepole pine and between lodgepole and jack pine. The mountain pine beetle will face a different constitutive and induced phloem resin terpene environment when locating and colonizing jack pine in its new geographic range, and this may play a significant role in the ability of the insect to persist in this new host.

Effects of introgression on the genetic population structure of two ecologically and economically important conifer species: lodgepole pine (Pinus contorta var. latifolia) and jack pine (Pinus banksiana)

Forest trees exhibit a remarkable range of adaptations to their environment, but as a result of frequent and long-distance gene flow, populations are often only weakly differentiated. Lodgepole and jack pine hybridize in western Canada, which adds the opportunity for introgression through hybridization to contribute to population structure and (or) adaptive variation. Access to large sample size, high density SNP datasets for these species would improve our ability to resolve population structure, parameterize introgression, and separate the influence of demography from adaptation. To accomplish this, 454 transcriptome reads for lodgepole and jack pine were assembled using Newbler and MIRA, the assemblies mined for SNPs, and 1536 SNPs were selected for typing on lodgepole pine, jack pine, and their hybrids (N = 536). We identified population structure using both Bayesian clustering and discriminate analysis of principle components. Introgressed SNP loci were identified and their influence on observed population structure was assessed. We found that introgressed loci resulted in increased differentiation both within lodgepole and jack pine populations. These findings are timely given the recent mountain pine beetle population expansion in the hybrid zone, and will facilitate future studies of adaptive traits in these ecologically important species.

The lodgepole × jack pine hybrid zone in Alberta, Canada: a stepping stone for the mountain pine beetle on its journey East across the boreal forest?

Historical data show that outbreaks of the tree killing mountain pine beetle are often preceded by periods of drought. Global climate change impacts drought frequency and severity and is implicated in the range expansion of the mountain pine beetle into formerly unsuitable habitats. Its expanded range has recently reached the lodgepole × jack pine hybrid zone in central Alberta, Canada, which could act as a transition from its historical lodgepole pine host to a jack pine host present in the boreal forest. This field study tested the effects of water limitation on chemical defenses of mature trees against mountain pine beetle-associated microorganisms and on beetle brood success in lodgepole × jack pine hybrid trees. Tree chemical defenses as measured by monoterpene emission from tree boles and monoterpene concentration in needles were greater in trees that experienced water deficit compared to well-watered trees. Myrcene was identified as specific defensive compound, since it significantly increased upon inoculation with dead mountain pine beetles. Beetles reared in bolts from trees that experienced water deficit emerged with a higher fat content, demonstrating for the first time experimentally that drought conditions benefit mountain pine beetles. Further, our study demonstrated that volatile chemical emission from tree boles and phloem chemistry place the hybrid tree chemotype in-between lodgepole pine and jack pine, which might facilitate the host shift from lodgepole pine to jack pine.