Incorporating fossils into the Pinaceae tree of life.

PREMISE OF THE STUDY: Pinaceae have a rich but enigmatic early fossil record, much of which is represented by permineralized seed cones. Our incomplete knowledge of morphology and anatomy in living and extinct species poses an important barrier to understanding their phylogenetic relationships and timing of diversification. METHODS: We expanded a morphology matrix to 46 fossil and 31 extant Pinaceae species, mainly adding characters from stem and leaf anatomy and seed cones. Using parsimony and Bayesian inference, we compared phylogenetic relationships for extant taxa with and without fossils from the morphology matrix combined with an alignment of plastid gene sequences. KEY RESULTS: Combined analysis of morphological and molecular characters resulted in a phylogeny of extant Pinaceae that was robust at all nodes except those relating to the interrelationships of Pinus, Picea, and Cathaya and the position of Cedrus. Simultaneous analysis of all fossil and extant species did not result in changes in the relationships among the extant species but did greatly reduce branch support. We found that the placement of most fossils was sensitive to the method of phylogenetic reconstruction when analyzing them singly with the extant species. CONCLUSIONS: A robust phylogenetic hypothesis for the main lineages of Pinaceae is emerging. Most Early Cretaceous fossils are stem or crown lineages of Pinus, but close relationships also were found between fossils and several other extant genera. The phylogenetic position of fossils broadly supports the existence of extant genera in the Lower Cretaceous.

Multi-locus phylogenetics, lineage sorting, and reticulation in Pinus subsection Australes.

PREMISE OF THE STUDY: Both incomplete lineage sorting and reticulation have been proposed as causes of phylogenetic incongruence. Disentangling these factors may be most difficult in long-lived, wind-pollinated plants with large population sizes and weak reproductive barriers. METHODS: We used solution hybridization for targeted enrichment and massive parallel sequencing to characterize low-copy-number nuclear genes and high-copy-number plastomes (Hyb-Seq) in 74 individuals of Pinus subsection Australes, a group of ~30 New World pine species of exceptional ecological and economic importance. We inferred relationships using methods that account for both incomplete lineage sorting and reticulation. KEY RESULTS: Concatenation- and coalescent-based trees inferred from nuclear genes mainly agreed with one another, but they contradicted the plastid DNA tree in recovering the Attenuatae (the California closed-cone pines) and Oocarpae (the egg-cone pines of Mexico and Central America) as monophyletic and the Australes sensu stricto (the southern yellow pines) as paraphyletic to the Oocarpae. The plastid tree featured some relationships that were discordant with morphological and geographic evidence and species limits. Incorporating gene flow into the coalescent analyses better fit the data, but evidence supporting the hypothesis that hybridization explains the non-monophyly of the Attenuatae in the plastid tree was equivocal. CONCLUSIONS: Our analyses document cytonuclear discordance in Pinus subsection Australes. We attribute this discordance to ancient and recent introgression and present a phylogenetic hypothesis in which mostly hierarchical relationships are overlain by gene flow.

Pinus ponderosa: A checkered past obscured four species.

PREMISE OF THE STUDY: Molecular genetic evidence can help delineate taxa in species complexes that lack diagnostic morphological characters. Pinus ponderosa (Pinaceae; subsection Ponderosae) is recognized as a problematic taxon: plastid phylogenies of exemplars were paraphyletic, and mitochondrial phylogeography suggested at least four subdivisions of P. ponderosa. These patterns have not been examined in the context of other Ponderosae species. We hypothesized that putative intraspecific subdivisions might each represent a separate taxon. METHODS: We genotyped six highly variable plastid simple sequence repeats in 1903 individuals from 88 populations of P. ponderosa and related Ponderosae (P. arizonica, P. engelmannii, and P. jeffreyi). We used multilocus haplotype networks and discriminant analysis of principal components to test clustering of individuals into genetically and geographically meaningful taxonomic units. KEY RESULTS: There are at least four distinct plastid clusters within P. ponderosa that roughly correspond to the geographic distribution of mitochondrial haplotypes. Some geographic regions have intermixed plastid lineages, and some mitochondrial and plastid boundaries do not coincide. Based on relative distances to other species of Ponderosae, these clusters diagnose four distinct taxa. CONCLUSIONS: Newly revealed geographic boundaries of four distinct taxa (P. benthamiana, P. brachyptera, P. scopulorum, and a narrowed concept of P. ponderosa) do not correspond completely with taxonomies. Further research is needed to understand their morphological and nuclear genetic makeup, but we suggest that resurrecting originally published species names would more appropriately reflect the taxonomy of this checkered classification than their current treatment as varieties of P. ponderosa.

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.

Estimating the species tree for Hawaiian Schiedea (Caryophyllaceae) from multiple loci in the presence of reticulate evolution.

Schiedea (Caryophyllaceae) is a monophyletic genus of 34 species, all endemic to the Hawaiian Islands, that arose from a single colonization, providing one of the best examples of adaptive radiation in Hawai'i. Species utilize a range of habitats and exhibit a variety of growth forms and transitions in breeding systems from hermaphroditism toward dimorphism or autogamy. Our study included the most thorough sampling to date: 2-5 individuals per species and 4 independent genetic partitions: eight plastid and three low-copy nuclear loci (9217bps), allowing a three-locus BEST species tree. Despite incomplete resolution at the tips, our results support monophyly for each extant species. Gene trees revealed several clear cases of cytonuclear incongruence, likely created by interspecific introgression. Conflict occurs at the divergence of section Alphaschiedea as well as at the tips. Ages inferred from a BEAST analysis allow an original colonization onto either Nihoa or Kauaì and inform some aspects of inter-island migrations. We suggest that several hard polytomies on the species tree are biologically realistic, signifying either nearly simultaneous speciation or historical introgressive hybridization. Based on inferred node ages that exceed expected coalescent times, we propose that undetected nuclear introgression may play a larger role than incomplete lineage sorting in sections Schiedea and Mononeura.

Reticulate evolution and incomplete lineage sorting among the ponderosa pines.

Interspecific gene flow via hybridization may play a major role in evolution by creating reticulate rather than hierarchical lineages in plant species. Occasional diploid pine hybrids indicate the potential for introgression, but reticulation is hard to detect because ancestral polymorphism is still shared across many groups of pine species. Nucleotide sequences for 53 accessions from 17 species in subsection Ponderosae (Pinus) provide evidence for reticulate evolution. Two discordant patterns among independent low-copy nuclear gene trees and a chloroplast haplotype are better explained by introgression than incomplete lineage sorting or other causes of incongruence. Conflicting resolution of three monophyletic Pinus coulteri accessions is best explained by ancient introgression followed by a genetic bottleneck. More recent hybridization transferred a chloroplast from P. jeffreyi to a sympatric P. washoensis individual. We conclude that incomplete lineage sorting could account for other examples of non-monophyly, and caution against any analysis based on single-accession or single-locus sampling in Pinus.

Chloroplast diversity in a putative hybrid swarm of Ponderosae (Pinaceae).

The Ponderosae subsection of the genus Pinus contains numerous taxa in disjunct mountain ranges of southern Arizona and New Mexico, differing for several leaf and cone traits, key among which is the number of leaf needles per fascicle. Trees with three needles are often found together with trees having five needles and mixed numbers. One taxonomic hypothesis is that there are swarms of hybrids between P. ponderosa and P. arizonica. A second hypothesis is that there are spatial mixtures of two separate taxa, five-needle P. arizonica and a "taxon X" containing three needle and mixed needle trees. We genotyped chloroplasts in one putative hybrid swarm on Mt. Lemmon using microsatellite markers and show that cpDNA is almost completely differentiated between two separate morphotypes corresponding to P. arizonica and "taxon X." Little if any introgression has occurred on Mt. Lemmon, and the simplest explanation is that little or no effective hybridization has occurred. Further results indicate that not only is taxon X not of hybrid origin, it is more closely related to nonregional Ponderosae other than P. ponderosa and P. arizonica. The results further suggest that other putative hybrid swarms in the region are also spatial mixtures of distinct taxa.

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.

Fossil calibration of molecular divergence infers a moderate mutation rate and recent radiations for pinus.

Silent mutation rate estimates for Pinus vary 50-fold, ranging from angiosperm-like to among the slowest reported for plants. These differences either reflect extraordinary genomic processes or inconsistent fossil calibration, and they have important consequences for population and biogeographical inferences. Here we estimate mutation rates from 4 Pinus species that represent the major lineages using 11 nuclear and 4 chloroplast loci. Calibration was tested at the divergence of Pinus subgenera with the oldest leaf fossil from subg. Strobus (Eocene; 45 MYA) or a recently published subg. Strobus wood fossil (Cretaceous; 85 MYA). These calibrations place the origin of Pinus 190-102 MYA and give absolute silent rate estimates of 0.70-1.31x10(-9) and 0.22-0.42x10(-9).site-1.year-1 for the nuclear and chloroplast genomes, respectively. These rates are approximately 4- to 20-fold slower than angiosperms, but unlike many previous estimates, they are more consistent with the high per-generation deleterious mutation rates observed in pines. Chronograms from nuclear and chloroplast genomes show that the divergence of subgenera accounts for about half of the time since Pinus diverged from Picea, with subsequent radiations occurring more recently. By extending the sampling to encompass the phylogenetic diversity of Pinus, we predict that most extant subsections diverged during the Miocene. Moreover, subsect. Australes, Ponderosae, and Contortae, containing over 50 extant species, radiated within a 5 Myr time span starting as recently as 18 MYA. An Eocene divergence of pine subgenera (using leaf fossils) does not conflict with fossil-based estimates of the Pinus-Picea split, but a Cretaceous divergence using wood fossils accommodates Oligocene fossils that may represent modern subsections. Because homoplasy and polarity of character states have not been tested for fossil pine assignments, the choice of fossil and calibration node represents a significant source of uncertainty. Based on several lines of evidence (including agreement with ages inferred using calibrations outside of Pinus), we conclude that the 85 MYA calibration at the divergence of pine subgenera provides a reasonable lower bound and that further refinements in age and mutation rate estimates will require a synthetic examination of pine fossil history.

Evolutionary relationships among Pinus (Pinaceae) subsections inferred from multiple low-copy nuclear loci.

Sequence data from nrITS and cpDNA have failed to fully resolve phylogenetic relationships among Pinus species. Four low-copy nuclear genes, developed from the screening of 73 mapped conifer anchor loci, were sequenced from 12 species representing all subsections. Individual loci do not uniformly support either the nrITS or cpDNA hypotheses and in some cases produce unique topologies. Combined analysis of low-copy nuclear loci produces a well-supported subsectional topology of two subgenera, each divided into two sections, congruent with prior hypotheses of deep divergence in Pinus. The placements of P. nelsonii, P. krempfii, and P. contorta have been of continued systematic interest. Results strongly support the placement of P. nelsonii as sister to the remaining members of sect. Parrya, suggest a moderately well-supported and consistent position of P. krempfii as sister to the remaining members of sect. Quinquefoliae, and are ambiguous about the placement of P. contorta. A successful phylogenetic strategy in Pinus will require many low-copy nuclear loci that include a high proportion of silent sites and derive from independent linkage groups. The locus screening and evaluation strategy presented here can be broadly applied to facilitate the development of phylogenetic markers from the increasing number of available genomic resources.