Mitochondrial DNA haplotype distribution patterns in Pinus ponderosa (Pinaceae): range-wide evolutionary history and implications for conservation

New uses for ancient middens: bridging ecological and evolutionary perspectives

Rodent middens provide a fine-scale spatiotemporal record of plant and animal communities over the late Quaternary. In the Americas, middens have offered insight into biotic responses to past environmental changes and historical factors influencing the distribution and diversity of species. However, few studies have used middens to investigate genetic or ecosystem level responses. Integrating midden studies with neoecology and experimental evolution can help address these gaps and test mechanisms underlying eco-evolutionary patterns across biological and spatiotemporal scales. Fully realizing the potential of middens to answer cross-cutting ecological and evolutionary questions and inform conservation goals in the Anthropocene will require a collaborative research community to exploit existing midden archives and mount new campaigns to leverage midden records globally.

Identifying genetic variation associated with environmental gradients and drought-tolerance phenotypes in ponderosa pine

As climate changes, understanding the genetic basis of local adaptation in plants becomes an ever more pressing issue. Combining genotype-environment association (GEA) with genotype-phenotype association (GPA) analysis has an exciting potential to uncover the genetic basis of environmental responses. We use these approaches to identify genetic variants linked to local adaptation to drought in Pinus ponderosa. Over 4 million Single Nucleotide Polymorphisms (SNPs) were identified using 223 individuals from across the Sierra Nevada of California. 927,740 (22.3%) SNPs were retained after filtering for proximity to genes and used in our association analyses. We found 1374 associated with five major climate variables, with the largest number (1151) associated with April 1st snowpack. We also conducted a greenhouse study with various drought-tolerance traits measured in first-year seedlings of a subset of the genotyped trees grown in the greenhouse. 796 SNPs were associated with control-condition trait values, while 1149 were associated with responsiveness of these traits to drought. While no individual SNPs were associated with both the environmental variables and the measured traits, several annotated genes were associated with both, particularly those involved in cell wall formation, biotic and abiotic stress responses, and ubiquitination. However, the functions of many of the associated genes have not yet been determined due to the lack of gene annotation information for conifers. Future studies are needed to assess the developmental roles and ecological significance of these unknown genes.

Plant mitochondrial introns as genetic markers - conservation and variation

Plant genomes are comprised of nuclear, plastid and mitochondrial components characterized by different patterns of inheritance and evolution. Genetic markers from the three genomes provide complementary tools for investigations of inheritance, genetic relationships and phenotypic contributions. Plant mitochondrial genomes are challenging for universal marker development because they are highly variable in terms of size, gene order and intergenic sequences and highly conserved with respect to protein-coding sequences. PCR amplification of introns with primers that anneal to conserved, flanking exons is effective for the development of polymorphic nuclear genome markers. The potential for plant mitochondrial intron polymorphisms to distinguish between congeneric species or intraspecific varieties has not been systematically investigated and is possibly constrained by requirements for intron secondary structure and interactions with co-evolved organelle intron splicing factors. To explore the potential for broadly applicable plant mitochondrial intron markers, PCR primer sets based upon conserved sequences flanking 11 introns common to seven angiosperm species were tested across a range of plant orders. PCR-amplified introns were screened for indel polymorphisms among a group of cross-compatible Citrus species and relatives; two Raphanus sativus mitotypes; representatives of the two Phaseolus vulgaris gene pools; and congeneric pairs of Cynodon, Cenchrus, Solanum, and Vaccinium species. All introns were successfully amplified from each plant entry. Length polymorphisms distinguishable by gel electrophoresis were common among genera but infrequent within genera. Sequencing of three introns amplified from 16 entries identified additional short indel polymorphisms and nucleotide substitutions that separated Citrus, Cynodon, Cenchrus and Vaccinium congeners, but failed to distinguish Solanum congeners or representatives of the Phaseolus vulgaris major gene pools. The ability of primer sets to amplify a wider range of plant species' introns and the presence of intron polymorphisms that distinguish congeners was confirmed by in silico analysis. While mitochondrial intron variation is limited in comparison to nuclear introns, these exon-based primer sets provide robust tools for the amplification of mitochondrial introns across a wide range of plant species wherein useful polymorphisms can be identified.

Not Only Severe Events: Moderate Dry Periods Impact the Hydraulic Functioning and Survival of Planted Ponderosa Pine Seedlings

Juvenile tree survival will increasingly shape the persistence of ponderosa pine forests in the western United States. In contrast to severe pulse disturbances that induce widespread adult and juvenile tree mortality, moderate periods of low rainfall and warm temperatures may reduce forest persistence by killing juvenile trees at the seedling stage. Intensification of these periods in a changing climate could therefore increasingly restrict both natural regeneration and artificial regeneration of planted seedlings. We conducted a controlled field experiment at a single site in the Front Range of Colorado, USA, to determine the responses and survival of 3 Colorado subpopulations of <1 year old potted ponderosa pines to moderately dry conditions, variation in small rainfall events based on observed patterns, and shaded and unshaded microsite environments. Near surface soil moisture increased slightly following small rainfall events, but declined over the 45-day experimental period. Seedling transpiration and associated canopy cooling declined after ∼13 days, and further declines in transpiration and canopy cooling suggest that the majority of trees in lower rainfall treatments experienced hydraulic dysfunction between days ∼20–30. After 45 days, mortality across all subpopulations and treatments, inferred by relative water loss, exceeded 90–95%. Despite some uncertainty pertaining to the stress tolerance of nursery grown versus naturally germinated conifers, our results show that planted ponderosa pine seedlings <1 year old are unlikely to survive moderate dry periods of 20+ days relying on small rainfall events. Although microsite conditions and soil moisture availability shaped tree hydraulic functioning early in the experiment (days 1–13), later functioning was shaped predominately by the legacy of rainfall treatments. Our results illustrate the importance of moderate dry events that occur consistently as part of seasonal variation in climate, and show how their intensification may constitute a sustained press that limits opportunities for natural and artificial regeneration.

Relative water content consistently predicts drought mortality risk in seedling populations with different morphology, physiology and times to death

Predicted increases in forest drought mortality highlight the need for predictors of incipient drought-induced mortality (DIM) risk that enable proactive large-scale management. Such predictors should be consistent across plants with varying morphology and physiology. Because of their integrative nature, indicators of water status are promising candidates for real-time monitoring of DIM, particularly if they standardize morphological differences among plants. We assessed the extent to which differences in morphology and physiology between Pinus ponderosa populations influence time to mortality and the predictive power of key indicators of DIM risk. Time to incipient mortality differed between populations but occurred at the same relative water content (RWC) and water potential (WP). RWC and WP were accurate predictors of drought mortality risk. These results highlight that variables related to water status capture critical thresholds during DIM and the associated dehydration processes. Both WP and RWC are promising candidates for large-scale assessments of DIM risk. RWC is of special interest because it allows comparisons across different morphologies and can be remotely sensed. Our results offer promise for real-time landscape-level monitoring of DIM and its global impacts in the near term.

Microhomologies Are Associated with Tandem Duplications and Structural Variation in Plant Mitochondrial Genomes

Short tandem repeats (STRs) contribute to structural variation in plant mitochondrial genomes, but the mechanisms underlying their formation and expansion are unclear. In this study, we detected high polymorphism in the nad7-1 region of the Pinus tabuliformis mitogenome caused by the rapid accumulation of STRs and rearrangements over a few million years ago. The STRs in nad7-1 have a 7-bp microhomology (TAG7) flanking the repeat array. We then scanned the mitogenomes of 136 seed plants to understand the role of microhomology in the formation of STR and mitogenome evolution. A total of 13,170 STRs were identified, and almost half of them were associated with microhomologies. A substantial amount (1,197) of microhomologies was long enough to mediate structural variation, and the length of microhomology is positively correlated with the length of tandem repeat unit. These results suggest that microhomology may be involved in the formation of tandem repeat via microhomology-mediated pathway, and the formation of longer duplicates required greater length of microhomology. We examined the abundance of these 1,197 microhomologies, and found 75% of them were enriched in the plant mitogenomes. Further analyses of the 400 prevalent microhomologies revealed that 175 of them showed differential enrichment between angiosperms and gymnosperms and 186 differed between angiosperms and conifers, indicating lineage-specific usage and expansion of microhomologies. Our study sheds light on the sources of structural variation in plant mitochondrial genomes and highlights the importance of microhomology in mitochondrial genome evolution.

Climatic-niche evolution with key morphological innovations across clades within Scutiger boulengeri (Anura: Megophryidae)

The studies of climatic-niche shifts over evolutionary time accompanied by key morphological innovations have attracted the interest of many researchers recently. We applied ecological niche models (ENMs), ordination method (environment principal component analyses; PCA-env), combined phylogenetic comparative methods (PCMs), and phylogenetic generalized least squares (PGLS) regression methods to analyze the realized niche dynamics and correspondingly key morphological innovations across clades within Scutiger boulengeri throughout their distributions in Qinghai-Tibet Plateau (QTP) margins of China. Our results show there are six clades in S. boulengeri and obvious niche divergences caused by niche expansion in three clades. Moreover, in our system, niche expansion is more popular than niche unfilling into novel environmental conditions. Annual mean temperature, annual precipitation, and precipitation of driest month may contribute to such a shift. In addition, we identified several key climatic factors and morphological traits that tend to be associated with niche expansion in S. boulengeri clades correspondingly. We found phenotypic plasticity [i.e., length of lower arm and hand (LAHL), hind-limb length (HLL), and foot length (FL)] and evolutionary changes [i.e., snout-vent length (SVL)] may together contribute to niche expansion toward adapting novel niche, which provides us a potential pattern of how a colonizing toad might seed a novel habitat to begin the process of speciation and finally adaptive radiation. For these reasons, persistent phylogeographic divisions and accompanying divergences in niche occupancy and morphological adaption suggest that for future studies, distinct genetic structure and morphological changes corresponding to each genetic clade should be included in modeling niche evolution dynamics, but not just constructed at the species level.

Geographical disjunction and environmental conditions drive intraspecific differentiation in the chalk-hill blue butterfly

Drivers of evolution are often related to geographical isolation and/or diverging environmental conditions. Spatial variation in neutral genetic markers mostly reflects past geographical isolation, i.e. long-lasting allopatry, whereas morphology is often driven by local environmental conditions, resulting in more rapid evolution. In Europe, most thermophilic species persisted during the past glacial periods in geographically disjunct refugia, representing long-lasting isolates, frequently with diverging environmental conditions. This situation has driven the evolution of intraspecific signatures in species. Here, we analysed wing shape and wing pigmentation of the chalk-hill blue butterfly, Polyommatus coridon, across its entire distribution range restricted to the western Palaearctic. In addition, we compiled abiotic environmental parameters for each sampling site. Wing colour patterns differentiated a western and an eastern lineage. These lineages might represent two main Pleistocene refugia and differentiation centres, one located on the Italian Peninsula and the other in the Balkan region. The two lineages showed evidence of hybridization across Central Europe, from the Alps and across Germany. The intraspecific differentiation was strongest in the width of the brown band on the outer margin of the wings. The morphological structures obtained are in line with genetic signatures found in previous studies, but the latter are more fine-grained. Current environmental conditions, such as mean temperatures, were only marginally correlated with colour patterns. Our study underlines that Pleistocene range shifts, often resulting in allopatric isolation, shape intraspecific phenotypic structures within species; that pigmentation responds in a more sensitive manner to spatial disjunction than wing shape; and that morphometric and genetic structures in P. coridon provide concordant patterns and thus support identical biogeographical conclusions.

Genetic structure of coast redwood (Sequoia sempervirens [D. Don] Endl.) populations in and outside of the natural distribution range based on nuclear and chloroplast microsatellite markers

Coast redwood (Sequoia sempervirens) naturally growing in southern Oregon and northern California is one of the few conifer tree species that are polyploid. Despite its unique ecological and economic importance, its population genetic structure is still insufficiently studied. To obtain additional data on its population genetic structure we genotyped 317 samples collected from populations in California (data set C) and 144 trees growing in a provenance trial in France (data set F) using 12 nuclear (five random nuclear genomic nSSRs and seven expressed sequence tag EST-SSRs) and six chloroplast (cpSSRs) microsatellite or simple sequence repeat (SSR) markers, respectively. These data sets were also used as reference to infer the origin of 147 coast redwood trees growing in Germany (data set G). Coast redwood was introduced to Europe, including Germany as an ornamental species, decades ago. Due to its fast growth and high timber quality, it could be considered as a potential commercial timber species, especially in perspective to climate warming that makes more regions in Germany suitable for its growing. The well performing trees in colder Germany could be potential frost resistant genotypes, but their genetic properties and origin are mostly unknown. Within the natural range in southern Oregon and northern California, only two relatively weak clusters were identified, one northern and one southern, separated by the San Francisco Bay. High genetic diversity, but low differentiation was found based on the 12 nuclear SSR markers for all three data sets F, C and G. We found that investigated 147 German trees represented only 37 different genotypes. They showed genetic diversity at the level less than diversity observed within the natural range in the northern or southern cluster, but more similar to the diversity observed in the southern cluster. It was difficult to assign German trees to the original single native populations using the six cpSSR markers, but rather to either the northern or southern cluster. The high number of haplotypes found in the data sets based on six cpSSR markers and low genetic differentiation based on 12 nuclear SSRs found in this study helps us study and better understand population genetic structure of this complex polyploid tree and supports the selection of potential genotypes for German forestry.

Ozone effects on plants in natural ecosystems

Tropospheric ozone (O₃ ) is an important stressor in natural ecosystems, with well-documented impacts on soils, biota and ecological processes. The effects of O₃ on individual plants and processes scale up through the ecosystem through effects on carbon, nutrient and hydrologic dynamics. Ozone effects on individual species and their associated microflora and fauna cascade through the ecosystem to the landscape level. Systematic injury surveys demonstrate that foliar injury occurs on sensitive species throughout the globe. However, deleterious impacts on plant carbon, water and nutrient balance can also occur without visible injury. Because sensitivity to O₃ may follow coarse physiognomic plant classes (in general, herbaceous crops are more sensitive than deciduous woody plants, grasses and conifers), the task still remains to use stomatal O₃ uptake to assess class and species' sensitivity. Investigations of the radial growth of mature trees, in combination with data from many controlled studies with seedlings, suggest that ambient O₃ reduces growth of mature trees in some locations. Models based on tree physiology and forest stand dynamics suggest that modest effects of O₃ on growth may accumulate over time, other stresses (prolonged drought, excess nitrogen deposition) may exacerbate the direct effects of O₃ on tree growth, and competitive interactions among species may be altered. Ozone exposure over decades may be altering the species composition of forests currently, and as fossil fuel combustion products generate more O₃ than deteriorates in the atmosphere, into the future as well.

Legacies of more frequent drought in ponderosa pine across the western United States

Despite widespread interest in drought legacies-multiyear impacts of drought on tree growth-the key implication of reported drought legacies remains unaddressed: as impaired growth and slow recovery associated with drought legacies are pervasive across forest ecosystems, what is the impact of more frequent drought conditions? We investigated the assumption that either multiple drought years occurring during a short period (multiyear droughts), or droughts occurring during the recovery period from previous drought (compounded droughts), are detrimental to subsequent growth. There is evidence that drought responses may vary among populations of widespread species, leading us to examine regional differences in responses of the conifer Pinus ponderosa to historic drought frequency in the western United States. More frequent drought conditions incurred additional growth declines and shifts in growth-climate sensitivities in the years following drought relative to single-drought events, with 'triple-droughts' being worse than 'double-droughts'. Notably, prediction skill was not strongly reduced when ignoring compounded droughts, a consequence of the temporally comprehensive formulation of our stochastic antecedent model that accounts for the climatic memory of tree growth. We argue that incorporating drought-induced temporal variability in tree growth sensitivities can aid inference gained from statistical models, where more simplistic models could overestimate the severity of drought legacies. We also found regional differences in response to repeated drought, and suggest plastic post-drought sensitivities and climatic memory may represent beneficial physiological adjustments in interior regions. Within-species variability may thus mediate forest responses to increasing drought frequency under future climate change, but experimental approaches using more species are necessary to improve our understanding of the mechanisms that underlie drought legacy effects on tree growth.

Coupled ecohydrology and plant hydraulics modeling predicts ponderosa pine seedling mortality and lower treeline in the US Northern Rocky Mountains

We modeled hydraulic stress in ponderosa pine seedlings at multiple scales to examine its influence on mortality and forest extent at the lower treeline in the northern Rockies. We combined a mechanistic ecohydrologic model with a vegetation dynamic stress index incorporating intensity, duration and frequency of hydraulic stress events, to examine mortality from loss of hydraulic conductivity. We calibrated our model using a glasshouse dry-down experiment and tested it using in situ monitoring data on seedling mortality from reforestation efforts. We then simulated hydraulic stress and mortality in seedlings within the Bitterroot River watershed of Montana. We show that cumulative hydraulic stress, its legacy and its consequences for mortality are predictable and can be modeled at local to landscape scales. We demonstrate that topographic controls on the distribution and availability of water and energy drive spatial patterns of hydraulic stress. Low-elevation, south-facing, nonconvergent locations with limited upslope water subsidies experienced the highest rates of modeled mortality. Simulated mortality in seedlings from 2001 to 2015 correlated with the current distribution of forest cover near the lower treeline, suggesting that hydraulic stress limits recruitment and ultimately constrains the low-elevation extent of conifer forests within the region.

Genome-wide markers reveal a complex evolutionary history involving divergence and introgression in the Abert's squirrel (Sciurus aberti) species group

Background

Genetic introgression between divergent lineages is now considered more common than previously appreciated, with potentially important consequences for adaptation and speciation. Introgression is often asymmetric between populations and patterns can vary for different types of loci (nuclear vs. organellar), complicating phylogeographic reconstruction. The taxonomy of the ecologically specialized Abert’s squirrel species group has been controversial, and previous studies based on mitochondrial data have not fully resolved the evolutionary relationships among populations. Moreover, while these studies identified potential areas of secondary contact between divergent lineages, the possibility for introgression has not been tested.

Results

We used RAD-seq to unravel the complex evolutionary history of the Abert’s squirrel species group. Although some of our findings reinforce inferences based on mitochondrial data, we also find significant areas of discordance. Discordant signals generally arise from previously undetected introgression between divergent populations that differentially affected variation at mitochondrial and nuclear loci. Most notably, our results support earlier claims (disputed by mitochondrial data) that S. aberti kaibabensis, found only on the north rim of the Grand Canyon, is highly divergent from other populations. However, we also detected introgression of S. aberti kaibabensis DNA into other S. aberti populations, which likely accounts for the previously inferred close genetic relationship between this population and those south of the Grand Canyon.

Conclusions

Overall, the evolutionary history of Abert’s squirrels appears to be shaped largely by divergence during periods of habitat isolation. However, we also found evidence for interbreeding during periods of secondary contact resulting in introgression, with variable effects on mitochondrial and nuclear markers. Our results support the emerging view that populations often diversify under scenarios involving both divergence in isolation and gene flow during secondary contact, and highlight the value of genome-wide datasets for resolving such complex evolutionary histories.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1248-4) contains supplementary material, which is available to authorized users.

Cascading speciation among mutualists and antagonists in a tree-beetle-fungi interaction

Cascading speciation is predicted to occur when multiple interacting species diverge in parallel as a result of divergence in one species promoting adaptive differentiation in other species. However, there are few examples where ecological interactions among taxa have been shown to result in speciation that cascades across multiple trophic levels. Here, we test for cascading speciation occurring among the western pine beetle (Dendroctonus brevicomis), its primary host tree (Pinus ponderosa), and the beetle's fungal mutualists (Ceratocystiopsis brevicomi and Entomocorticium sp. B). We assembled genomes for the beetle and a fungal symbiont and then generated reduced representation genomic data (RADseq) from range-wide samples of these three interacting species. Combined with published data for the host tree, we present clear evidence that the tree, the beetle, and the fungal symbionts are all genetically structured into at least two distinct groups that have strongly codiverged with geographical isolation. We then combine our genomic results with diverse population and laboratory-based data to show evidence for reproductive isolation at each level of the cascade and for coevolution of both antagonistic and mutualistic species interactions within this complex network.

Distributions of ectomycorrhizal and foliar endophytic fungal communities associated with Pinus ponderosa along a spatially constrained elevation gradient

PREMISE OF THE STUDY

Understanding distributions of plant-symbiotic fungi is important for projecting responses to environmental change. Many coniferous trees host ectomycorrhizal fungi (EM) in association with roots and foliar endophytic fungi (FE) in leaves. We examined how EM and FE associated with Pinus ponderosa each vary in abundance, diversity, and community structure over a spatially constrained elevation gradient that traverses four plant communities, 4°C in mean annual temperature, and 15 cm in mean annual precipitation.

METHODS

We sampled 63 individuals of Pinus ponderosa in 10 sites along a 635 m elevation gradient that encompassed a geographic distance of 9.8 km. We used standard methods to characterize each fungal group (amplified and sequenced EM from root tips; isolated and sequenced FE from leaves).

KEY RESULTS

Abundance and diversity of EM were similar across sites, but community composition and distributions of the most common EM differed with elevation (i.e., with climate, soil chemistry, and plant communities). Abundance and composition of FE did not differ with elevation, but diversity peaked in mid-to-high elevations.

CONCLUSIONS

Our results suggest relatively tight linkages between EM and climate, soil chemistry, and plant communities. That FE appear less linked with these factors may speak to limitations of a culture-based approach, but more likely reflects the small spatial scale encompassed by our study. Future work should consider comparable methods for characterizing these functional groups, and additional transects to understand relationships of EM and FE to environmental factors that are likely to shift as a function of climate change.

Reproductive isolation and environmental adaptation shape the phylogeography of mountain pine beetle (Dendroctonus ponderosae)

Chromosomal rearrangement can be an important mechanism driving population differentiation and incipient speciation. In the mountain pine beetle (MPB, Dendroctonus ponderosae), deletions on the Y chromosome that are polymorphic among populations are associated with reproductive incompatibility. Here, we used RAD sequencing across the entire MPB range in western North America to reveal the extent of the phylogeographic differences between Y haplotypes compared to autosomal and X-linked loci. Clustering and geneflow analyses revealed three distinct Y haplogroups geographically positioned within and on either side of the Great Basin Desert. Despite close geographic proximity between populations on the boundaries of each Y haplogroup, there was extremely low Y haplogroup mixing among populations, and gene flow on the autosomes was reduced across Y haplogroup boundaries. These results are consistent with a previous study suggesting that independent degradation of a recently evolved neo-Y chromosome in previously isolated populations causes male sterility or inviability among Y haplotype lineages. Phylogeographic results supported historic contraction of MPB into three separate Pleistocene glacial refugia followed by postglacial range expansion and secondary contact. Distinct sets of SNPs were statistically associated with environmental data among the most genetically distinct sets of geographic populations. This finding suggests that the process of adaptation to local climatic conditions is influenced by population genetic structure, with evidence for largely independent evolution in the most genetically isolated Y haplogroup.

Exploring Climate Niches of Ponderosa Pine (Pinus ponderosa Douglas ex Lawson) Haplotypes in the Western United States: Implications for Evolutionary History and Conservation

Ponderosa pine (Pinus ponderosa Douglas ex Lawson) occupies montane environments throughout western North America, where it is both an ecologically and economically important tree species. A recent study using mitochondrial DNA analysis demonstrated substantial genetic variation among ponderosa pine populations in the western U.S., identifying 10 haplotypes with unique evolutionary lineages that generally correspond spatially with distributions of the Pacific (P. p. var. ponderosa) and Rocky Mountain (P. p. var. scopulorum) varieties. To elucidate the role of climate in shaping the phylogeographic history of ponderosa pine, we used nonparametric multiplicative regression to develop predictive climate niche models for two varieties and 10 haplotypes and to hindcast potential distribution of the varieties during the last glacial maximum (LGM), ~22,000 yr BP. Our climate niche models performed well for the varieties, but haplotype models were constrained in some cases by small datasets and unmeasured microclimate influences. The models suggest strong relationships between genetic lineages and climate. Particularly evident was the role of seasonal precipitation balance in most models, with winter- and summer-dominated precipitation regimes strongly associated with P. p. vars. ponderosa and scopulorum, respectively. Indeed, where present-day climate niches overlap between the varieties, introgression of two haplotypes also occurs along a steep clinal divide in western Montana. Reconstructed climate niches for the LGM suggest potentially suitable climate existed for the Pacific variety in the California Floristic province, the Great Basin, and Arizona highlands, while suitable climate for the Rocky Mountain variety may have existed across the southwestern interior highlands. These findings underscore potentially unique phylogeographic origins of modern ponderosa pine evolutionary lineages, including potential adaptations to Pleistocene climates associated with discrete temporary glacial refugia. Our predictive climate niche models may inform strategies for further genetic research (e.g., sampling design) and conservation that promotes haplotype compatibility with projected changes in future climate.

A genomic investigation of the putative contact zone between divergent Brown Creeper (Certhia americana) lineages: chromosomal patterns of genetic differentiation

Sky islands, or montane forest separated by different lowland habitats, are highly fragmented regions that potentially limit gene flow between isolated populations. In the sky islands of the Madrean Archipelago (Arizona, USA), various taxa display different phylogeographic patterns, from unrestricted gene flow among sky islands to complex patterns with multiple distinct lineages. Using genomic-level approaches allows the investigation of differential patterns of gene flow, selection, and genetic differentiation among chromosomes and specific genomic regions between sky island populations. Here, we used thousands of SNPs to investigate the putative contact zone of divergent Brown Creeper (Certhia americana) lineages in the Madrean Archipelago sky islands. We found the two lineages to be completely allopatric (during the breeding season) with a lack of hybridization and gene flow between lineages and no genetic structure among sky islands within lineages. Additionally, the two lineages inhabit different climatic and ecosystem conditions and have many local primary song dialects in the southern Arizona mountain ranges. We identified a positive relationship between genetic differentiation and chromosome size, but the sex chromosome (Z) was not found to be an outlier. Differential patterns of genetic differentiation per chromosome may be explained by genetic drift--possibly in conjunction with non-random mating and non-random gene flow--due to variance in recombination rates among chromosomes.

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.

A set of plastid loci for use in multiplex fragment length genotyping for intraspecific variation in Pinus (Pinaceae)

PREMISE OF THE STUDY

Recently released Pinus plastome sequences support characterization of 15 plastid simple sequence repeat (cpSSR) loci originally published for P. contorta and P. thunbergii. This allows selection of loci for single-tube PCR multiplexed genotyping in any subsection of the genus. •

METHODS

Unique placement of primers and primer conservation across the genus were investigated, and a set of six loci were selected for single-tube multiplexing. We compared interspecific variation between cpSSRs and nucleotide sequences of ycf1 and tested intraspecific variation for cpSSRs using 911 samples in the P. ponderosa species complex. •

RESULTS

The cpSSR loci contain mononucleotide and complex repeats with additional length variation in flanking regions. They are not located in hypervariable regions, and most primers are conserved across the genus. A single PCR per sample multiplexed for six loci yielded 45 alleles in 911 samples. •

DISCUSSION

The protocol allows efficient genotyping of many samples. The cpSSR loci are too variable for Pinus phylogenies but are useful for the study of genetic structure within and among populations. The multiplex method could easily be extended to other plant groups by choosing primers for cpSSR loci in a plastome alignment for the target group.