Mitochondrial haplotype distribution, seed dispersal and patterns of postglacial expansion of ponderosa pine

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.

Evolutionary History of Rhus chinensis (Anacardiaceae) From the Temperate and Subtropical Zones of China Based on cpDNA and Nuclear DNA Sequences and Ecological Niche Model

To explore the origin and evolution of local flora and vegetation, we examined the evolutionary history of Rhus chinensis, which is widely distributed in China’s temperate and subtropical zones, by sequencing three maternally inherited chloroplast DNAs (cpDNA: trnL-trnF, psbA-trnH, and rbcL) and the biparentally inherited nuclear DNA (nuDNA: LEAFY) from 19 natural populations of R. chinensis as well as the ecological niche modeling. In all, 23 chloroplast haplotypes (M1–M23) and 15 nuclear alleles (N1–N15) were detected. The estimation of divergence time showed that the most recent common ancestor dated at 4.2 ± 2.5 million years ago (Mya) from cpDNA, and the initial divergence of genotypes occurred at 4.8 ± 3.6 Mya for the nuDNA. Meanwhile, the multimodality mismatch distribution curves and positive Tajima’s D values indicated that R. chinensis did not experience population expansion after the last glacial maximum. Besides, our study was also consistent with the hypothesis that most refugia in the temperate and subtropical zones of China were in situ during the glaciation.

Genetic diversity and relationship between cultivated, weedy and wild rye species as revealed by chloroplast and mitochondrial DNA non-coding regions analysis

The genus Secale is small but very diverse. Despite the high economic importance, phylogenetic relationships of rye species have not been fully determined, and they are extremely important for the process of breeding of new cultivars that can be enriched with functional traits derived from wild rye species. The study analyzed the degree of relationship of 35 accessions of the genus Secale, representing 13 most often distinguished species and subspecies, originating from various seed collections in the world, based on the analysis of non-coding regions of the chloroplast (cpDNA) and mitochondrial genome (mtDNA), widely used in phylogenetic and population plant studies, because of a higher rate of evolution than the coding regions. There was no clear genetic structure between different species and subspecies, which may indicated the introgression between these taxa. The obtained data confirmed that S. vavilovii was very similar to S. cereale, which confirmed the assumption that they might share a common ancestor. The results also confirmed the divergence of S. sylvestre from other species and subspecies of rye. Areas that may be useful molecular markers in studies on closely related species of the genus Secale were also indicated.

Phylogeographic structure in three North American tent caterpillar species (Lepidoptera: Lasiocampidae): Malacosoma americana, M. californica, and M. disstria

While phylogeographic structure has been examined in many North American vertebrate species, insects have received much less attention despite their central ecological roles. The moth genus Malacosoma (Hübner, 1820), is an important group of forestry pests responsible for large-scale defoliation across much of the Nearctic and Palearctic. The present study uses sequence variation in the mitochondrial cytochrome c oxidase 1 (COI) gene to examine the population genetic structure of the three widespread Malacosoma species (M. americana, M. californica, and M. disstria). Populations of all three species showed highest diversity in the south, suggesting that modern populations derived from southern refugia with loss of variation as these lineages dispersed northwards. However, despite similar life histories and dispersal abilities, the extent of regional variation varied among the taxa. M. americana, a species restricted to eastern North America, showed much less genetic structure than the western M. californica or the widespread M. disstria. The regional differentiation in the latter reflects the likely derivation of modern lineages from several refugia, as well as taxonomic uncertainty in M. californica. In these respects, the three species of Malacosoma share phylogeographic patterns similar to those detected in vertebrates which are characterised by greater phylogeographic breaks in the western half of the continent and limited structure in the east.

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.

Pleistocene glaciations, demographic expansion and subsequent isolation promoted morphological heterogeneity: A phylogeographic study of the alpine Rosa sericea complex (Rosaceae)

While most temperate plants probably underwent glacial constriction to refugia and interglacial expansion, another type of interglacial refugia might have existed to maintain alpine plants during warm periods. To test this hypothesis, we applied phylogeographic methods to 763 individuals (62 populations) which belong to 7 taxonomically difficult species of the Rosa sericea complex distributed in alpine regions of the temperate and subtropical zones in eastern Asia. We used three chloroplast (cp) DNA fragments (trnL-trnF, ndhF-rpl32 and ndhJ-trnF) approximately 3,100 bp and nuclear microsatellite (nSSR) on eight sites to determine whether cold tolerant plants experienced expansion during the Pleistocene. The neutral test and mismatch distribution analysis (MDA) indicated that whole populations and major lineages of the Qinghai-Tibet Plateau (QTP) underwent expansion during the middle to late Pleistocene. Environmental niche modeling (ENM) indicates more suitable habitats during the Last Glacial Maximum (LGM) than at present. We concluded that the demographic history of R. sericea, which diverged in the middle Pleistocene, was mostly affected by climatic oscillations instead of by geographical barriers. The low genetic divergence, as well as the weak phylogenetic structure in the R. sericea complex both support treating this complex as a single taxon.

Molecular markers reveal limited population genetic structure in a North American corvid, Clark's nutcracker (Nucifraga columbiana)

The genetic impact of barriers and Pleistocene glaciations on high latitude resident species has not been widely investigated. The Clark’s nutcracker is an endemic North American corvid closely associated with Pinus-dominated forests. The nutcracker’s encompasses known barriers to dispersal for other species, and glaciated and unglaciated areas. Clark’s nutcrackers also irruptively disperse long distances in search of pine seed crops, creating the potential for gene flow among populations. Using the highly variable mitochondrial DNA control region, seven microsatellite loci, and species distribution modeling, we examined the effects of glaciations and dispersal barriers on population genetic patterns and population structure of nutcrackers. We sequenced 900 bp of mitochondrial control region for 169 individuals from 15 populations and analysed seven polymorphic microsatellite loci for 13 populations across the Clark’s nutcracker range. We used species distribution modeling and a range of phylogeographic analyses to examine evolutionary history. Clark’s nutcracker populations are not highly differentiated throughout their range, suggesting high levels of gene flow among populations, though we did find some evidence of isolation by distance and peripheral isolation. Our analyses suggested expansion from a single refugium after the last glacial maximum, but patterns of genetic diversity and paleodistribution modeling of suitable habitat were inconclusive as to the location of this refugium. Potential barriers to dispersal (e.g. mountain ranges) do not appear to restrict gene flow in Clark’s nutcracker, and postglacial expansion likely occurred quickly from a single refugium located south of the ice sheets.

Is ecological speciation a major trend in aphids? Insights from a molecular phylogeny of the conifer-feeding genus Cinara

Introduction

In the past decade ecological speciation has been recognized as having an important role in the diversification of plant-feeding insects. Aphids are host-specialised phytophagous insects that mate on their host plants and, as such, they are prone to experience reproductive isolation linked with host plant association that could ultimately lead to species formation. The generality of such a scenario remains to be tested through macroevolutionary studies. To explore the prevalence of host-driven speciation in the diversification of the aphid genus Cinara and to investigate alternative modes of speciation, we reconstructed a phylogeny of this genus based on mitochondrial, nuclear and Buchnera aphidicola DNA sequence fragments and applied a DNA-based method of species delimitation. Using a recent software (PhyloType), we explored evolutionary transitions in host-plant genera, feeding sites and geographic distributions in the diversification of Cinara and investigated how transitions in these characters have accompanied speciation events.

Results

The diversification of Cinara has been constrained by host fidelity to conifer genera sometimes followed by sequential colonization onto different host species and by feeding-site specialisation. Nevertheless, our analyses suggest that, at the most, only half of the speciation events were accompanied by ecological niche shifts. The contribution of geographical isolation in the speciation process is clearly apparent in the occurrence of species from two continents in the same clades in relatively terminal positions in our phylogeny. Furthermore, in agreement with predictions from scenarios in which geographic isolation accounts for speciation events, geographic overlap between species increased significantly with time elapsed since their separation.

Conclusions

The history of Cinara offers a different perspective on the mode of speciation of aphids than that provided by classic models such as the pea aphid. In this genus of aphids, the role of climate and landscape history has probably been as important as host-plant specialisation in having shaped present-day diversity.

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

PREMISE OF THE STUDY

Ponderosa pine (Pinus ponderosa Douglas ex P. Lawson & C. Lawson) exhibits complicated patterns of morphological and genetic variation across its range in western North America. This study aims to clarify P. ponderosa evolutionary history and phylogeography using a highly polymorphic mitochondrial DNA marker, with results offering insights into how geographical and climatological processes drove the modern evolutionary structure of tree species in the region.

METHODS

We amplified the mtDNA nad1 second intron minisatellite region for 3,100 trees representing 104 populations, and sequenced all length variants. We estimated population-level haplotypic diversity and determined diversity partitioning among varieties, races and populations. After aligning sequences of minisatellite repeat motifs, we evaluated evolutionary relationships among haplotypes.

KEY RESULTS

The geographical structuring of the 10 haplotypes corresponded with division between Pacific and Rocky Mountain varieties. Pacific haplotypes clustered with high bootstrap support, and appear to have descended from Rocky Mountain haplotypes. A greater proportion of diversity was partitioned between Rocky Mountain races than between Pacific races. Areas of highest haplotypic diversity were the southern Sierra Nevada mountain range in California, northwestern California, and southern Nevada.

CONCLUSIONS

Pinus ponderosa haplotype distribution patterns suggest a complex phylogeographic history not revealed by other genetic and morphological data, or by the sparse paleoecological record. The results appear consistent with long-term divergence between the Pacific and Rocky Mountain varieties, along with more recent divergences not well-associated with race. Pleistocene refugia may have existed in areas of high haplotypic diversity, as well as the Great Basin, Southwestern United States/northern Mexico, and the High Plains.

Development of genetic diversity, differentiation and structure over 500 years in four ponderosa pine populations

Population history plays an important role in shaping contemporary levels of genetic variation and geographic structure. This is especially true in small, isolated range-margin populations, where effects of inbreeding, genetic drift and gene flow may be more pronounced than in large continuous populations. Effects of landscape fragmentation and isolation distance may have implications for persistence of range-margin populations if they are demographic sinks. We studied four small, disjunct populations of ponderosa pine over a 500-year period. We coupled demographic data obtained through dendroecological methods with microsatellite data to discern how and when contemporary levels of allelic diversity, among and within-population levels of differentiation, and geographic structure, arose. Alleles accumulated rapidly following initial colonization, demonstrating proportionally high levels of gene flow into the populations. At population sizes of approximately 100 individuals, allele accumulation saturated. Levels of genetic differentiation among populations (F(ST) and Jost's D(est)) and diversity within populations (F(IS)) remained stable through time. There was no evidence of geographic genetic structure at any time in the populations' history. Proportionally, high gene flow in the early stages of population growth resulted in rapid accumulation of alleles and quickly created relatively homogenous genetic patterns among populations. Our study demonstrates that contemporary levels of genetic diversity were formed quickly and early in population development. How contemporary genetic diversity accumulates over time is a key facet of understanding population growth and development. This is especially relevant given the extent and speed at which species ranges are predicted to shift in the coming century.

Interspecific phylogenetic analysis enhances intraspecific phylogeographical inference: a case study in Pinus lambertiana

Pinus lambertiana (sugar pine) is an economically and ecologically important conifer with a 1600-km latitudinal range extending from Oregon, USA, to northern Baja California, Mexico. Like all North American white pines (subsect. Strobus), sugar pine is highly susceptible to white pine blister rust, a disease caused by the fungus Cronartium ribicola. We conducted a chloroplast DNA (cpDNA) survey of Pinus subsect. Strobus with comprehensive geographical sampling of P. lambertiana. Sequence analysis of 12 sugar pine individuals revealed strong geographical differentiation for two chloroplast haplotypes. A diagnostic restriction site survey of an additional 72 individuals demarcated a narrow 150-km contact zone in northeastern California. In the contact zone, maternal (megagametophtye) and paternal (embryo) haplotypes were identified in 31 single seeds, demonstrating bidirectional pollen flow extending beyond the range of maternal haplotypes. The frequencies of the Cr1 allele for white pine blister rust major gene resistance, previously determined for 41 seed zones, differ significantly among seed zones that are fixed for the alternate haplotypes, or contain a mixture of both haplotypes. Interspecific phylogenetic analysis reveals that the northern sugar pine haplotype belongs to a clade that includes Pinus albicaulis (whitebark pine) and all of the East Asian white pines. Furthermore, there is little cpDNA divergence between northern sugar pine and whitebark pine (dS = 0.00058). These results are consistent with a Pleistocene migration of whitebark pine into North America and subsequent chloroplast introgression from whitebark pine to sugar pine. This study demonstrates the importance of placing phylogeographical results in a broader phylogenetic context.

Phylogeography of a northeast Asian spruce, Picea jezoensis, inferred from genetic variation observed in organelle DNA markers

Range-wide genetic variation of the widespread cold-temperate spruce Picea jezoensis was studied throughout northeast Asia using maternally inherited mitochondrial DNA and paternally inherited chloroplast DNA markers. This study assessed 33 natural populations including three varieties of the species in Japan, Russia, China, and South Korea. We depicted sharp suture zones in straits around Japan in the geographical distribution pattern of mitochondrial haplotypes (GST=0.901; NST=0.934). In contrast, we detected possible extensive pollen flow without seed flow across the straits around Japan during the past population history in the distribution pattern of chloroplast haplotypes (GST=0.233; NST=0.333). The analysis of isolation by distance of the species implied that by acting as a barrier for the movement of seeds and pollen, the sharp suture zones contributed considerably to the level of genetic differentiation between populations. Constructed networks of mitochondrial haplotypes allowed inference of the phylogeographical history of the species. We deduced that the disjunction with Kamchatka populations reflects range expansion and contraction to the north of the current distribution. Within Japan, we detected phylogeographically different types of P. jezoensis between Hokkaido and Honshu islands; P. jezoensis in Honshu Island may have colonized this region from the Asian continent via the Korean peninsula and the species in Hokkaido Island is likely to have spread from the Asian continent via Sakhalin through land bridges. Japanese endemism of mitochondrial haplotypes in Hokkaido and Honshu islands might have been promoted by separation of these islands from each other and from the Asian continent by the straits during the late Quaternary.

Colonization patterns of the invasive Brazilian peppertree, Schinus terebinthifolius, in Florida

Invasive species are believed to spread through a process of stratified dispersal consisting of short-distance diffusive spread around established foci and human mediated long-distance jumps. Brazilian peppertree (Schinus terebinthifolius), native to South America, was introduced twice as an ornamental plant into Florida, USA, just over 100 years ago. A previous study indicated that these two introductions were from genetically differentiated source populations in the native range. In this study, we took advantage of these contrasting genetic signatures to study the spatial spread of Brazilian peppertree across its entire range in Florida. A combination of spatial genetic and geostatistical analyses using chloroplast and nuclear microsatellite markers revealed evidence for both diffusive dispersal and long-distance jumps. Chloroplast DNA haplotype distributions and extensive bands of intra-specific hybridization revealed extensive dispersal by both introduced populations across the state. The strong genetic signature around the original introduction points, the presence of a general southeast to northwest genetic cline, and evidence for short-distance genetic spatial autocorrelation provided evidence of diffusive dispersal from an advancing front, probably by birds and small mammals. In the northernmost part of the range, there were patches having a high degree of ancestry from each introduction, suggesting long-distance jump dispersal, probably by the movement of humans. The evidence for extensive movement throughout the state suggests that Brazilian peppertree will be capable of rapidly recolonizing areas from which it has been eradicated. Concerted eradication efforts over large areas or the successful establishment of effective biocontrol agents over a wide area will be needed to suppress this species.

Molecular phylogeny and biogeography of Picea (Pinaceae): Implications for phylogeographical studies using cytoplasmic haplotypes

The center of diversity is not necessarily the place of origin, as has been established by many plant molecular phylogenies. Picea is a complicated but very important genus in coniferous forests of the Northern Hemisphere, with a high species diversity in Asia. Its phylogeny and biogeography were investigated here using sequence analysis of the paternally inherited chloroplast trnC-trnD and trnT-trnF regions and the maternally inherited mitochondrial nad5 intron 1. We found that the North American P. breweriana and P. sitchensis were basal to the other spruces that were further divided into three clades in the cpDNA phylogeny, and that the New World species harbored four of five mitotypes detected, including two ancestral ones and three endemics. These results, combined with biogeographic analyses using DIVA and MacClade and fossil evidence, suggest that Picea originated in North America, and that its present distribution could stem from two times of dispersal from North America to Asia by the Beringian land bridge, and then from Asia to Europe. Most of the northeastern Asian species and the European P. abies could arise from a recent radiation given the very low interspecific genetic differentiation and pure mitotype of them. Considering that the ancestral mtDNA polymorphism can be preserved in many descendant species, even distantly related ones, we suggest that more species, at least the closely related ones, should be sampled in the phylogeographical study using cytoplasmic haplotypes if possible. In addition, we also discussed the evolution and phylogenetic utility of morphological characters in Picea.

Effects of Pleistocene glaciations on population structure of North American chestnut-backed chickadees

The postglacial recolonization of northern North America was heavily influenced by the Pleistocene glaciation. In the Pacific Northwest, there are two disjunct regions of mesic temperate forest, one coastal and the other interior. The chestnut-backed chickadee is one of the species associated with this distinctive ecosystem. Using seven microsatellite markers we found evidence of population structure among nine populations of chestnut-backed chickadees. High levels of allelic variation were found in each of the populations. Northern British Columbia and central Alaska populations contained a large number of private alleles compared to other populations, including those from unglaciated regions. The disjunct population in the interior was genetically distinct from the coastal population. Genetic and historical records indicate that the interior population originated from postglacial inland dispersal. Population structuring was found within the continuous coastal population, among which the peripheral populations, specifically those on the Queen Charlotte Islands and the central Alaska mainland, were genetically distinct. The pattern of population structure among contemporary chickadee populations is consistent with a pioneer model of recolonization. The persistence of genetic structure in western North American chestnut-backed chickadees may be aided by their sedentary behaviour, linear distribution, and dependence on cedar-hemlock forests.

Distribution of chloroplast DNA haplotypes in the contact zone of Fagus crenata in the southwest of Kanto District, Japan

We have investigated cpDNA haplotype distribution in 24 populations of Fagus crenata in the southwest of Kanto District, Japan, and clarified the extent of intermixing of haplotypes in the contact zone by additional fine-scale analysis of two areas. Two cpDNA haplotypes belonging to different lineages were detected, and their distribution had geographical structure. Intermixed populations with the two haplotypes were limited to a narrow area. The geographical boundary between the haplotypes extended from Hakone to the west of the Kanto Mountains through the northern foot of Mt Fuji. No relationship was observed between the boundary location and the current topography of the southwest of Kanto District.

Tension versus ecological zones in a two-locus system

Previous theories show that tension and ecological zones are indistinguishable in terms of gene frequency clines. Here I analytically show that these two types of zones can be distinguished in terms of genetic statistics other than gene frequency. A two-locus cline model is examined with the assumptions of random mating, weak selection, no drift, no mutation, and multiplicative viabilities. The genetic statistics for distinguishing the two types of zones are the deviations of one- or two-locus genotypic frequencies from Hardy-Weinberg equilibrium (HWE) or from random association of gametes (RAG), and the deviations of additive and dominance variances from the values at HWE. These deviations have a discontinuous distribution in space and different extents of interruptions in the ecological zone with a sharp boundary, but exhibit a continuous distribution in the tension zone. Linkage disequilibrium enhances the difference between the deviations from HWE and from RAG for any two-locus genotypic frequency.

Recolonization and radiation in Larix (Pinaceae): evidence from nuclear ribosomal DNA paralogues

Gene paralogy frequently causes the conflict between gene tree and species tree, but sometimes the coexistence of a few paralogous copies could provide more markers for tracing the phylogeographical process of some organisms. In the present study, nrDNA ITS paralogues were cloned from all but one species of Larix, an Eocene genus having two sections, Larix and Multiserialis, with a huge circumboreal distribution and an Eastern Asia-Western North America disjunction, respectively. A total of 96 distinct clones, excluding five putative pseudogenes or recombinants, were obtained and used in the gene genealogy analysis. The clones from all Eurasian species of section Larix are mixed together, suggesting that recolonization and recent morphological differentiation could have played important roles in the evolution of this section. In contrast, the species diversification of the Eurasian section Multiserialis may result from radiation in the east Himalayas and its vicinity, considering extensive nrDNA founder effects in this group. Our study also suggests that the distribution pattern analysis of members of multiple gene family would be very useful in tracking the evolutionary history of some taxa with recent origin or rapid radiation that cannot be resolved by other molecular markers.

Impact of Ice Ages on the genetic structure of trees and shrubs

Data on the genetic structure of tree and shrub populations on the continental scale have accumulated dramatically over the past decade. However, our ability to make inferences on the impact of the last ice age still depends crucially on the availability of informative palaeoecological data. This is well illustrated by the results from a recent project, during which new pollen fossil maps were established and the variation in chloroplast DNA was studied in 22 European species of trees and shrubs. Species exhibit very different levels of genetic variation between and within populations, and obviously went through very different histories after Ice Ages. However, when palaeoecological data are non-informative, inferences on past history are difficult to draw from entirely genetic data. On the other hand, as illustrated by a study in ponderosa pine, when we can infer the species' history with some certainty, coalescent simulations can be used and new hypotheses can be tested.