Steven J. Brunsfeld

INVESTIGATING THE EVOLUTIONARY HISTORY OF THE PACIFIC NORTHWEST MESIC FOREST ECOSYSTEM: HYPOTHESIS TESTING WITHIN A COMPARATIVE PHYLOGEOGRAPHIC FRAMEWORK

Abstract We examine the evolution of mesic forest ecosystems in the Pacific Northwest of North America using a statistical phylogeography approach in four animal and two plant lineages. Three a priori hypotheses, which explain the disjunction in the mesic forest ecosystem with either recent dispersal or ancient vicariance, are tested with phylogenetic and coalescent methods. We find strong support in three amphibian lineages (Ascaphus spp., and Dicampton spp., and Plethodon vandykei and P. idahoensis) for deep divergence between coastal and inland populations, as predicted by the ancient vicariance hypothesis. Unlike the amphibians, the disjunction in other Pacific Northwest lineages is likely due to recent dispersal along a northern route. Topological and population divergence tests support the northern dispersal hypothesis in the water vole (Microtus richardsoni) and northern dispersal has some support in both the dusky willow (Salix melanopsis) and whitebark pine (Pinus albicaulis). These analyses demonstrate that genetic data sampled from across an ecosystem can provide insight into the evolution of ecological communities and suggest that the advantages of a statistical phylogeographic approach are most pronounced in comparisons across multiple taxa in a particular ecosystem. Genetic patterns in organisms as diverse as willows and salamanders can be used to test general regional hypotheses, providing a consistent metric for comparison among members of an ecosystem with disparate life‐history traits.

Interspecific hybridization of plants and resistance to herbivores : hypotheses, genetics, and variable responses in a diverse herbivore community

We studied the morphology, molecular genetics, and hebivory of two species of willows (Salix sericea and S. eriocephala) and their interspecific hybrids to test four alternative hypotheses concerning the effects of hybridization on plant resistance. Individually marked plants were identified using morphological traits in the field and measurements of stipule and leaf pubescence were made and compared using Canonical Discriminant Function Analysis. DNA was extracted from the leaves of a sample of the marked plants and RAPD-PCR analysis was performed to establish the genetic status of parental and hybrid plants. RAPD band analysis generally verified the genetic status of parental plants. Hybrid plants were usually correctly identified in the field with a few exceptions. However, the hybrid plants were a heterogeneous group of plants made up of most plants that appear to be F1s and a few plants that appear to be backcrosses to S. sericea. Morphological variables were useful for distinguishing S. sericea from S. eriocephala and hybrids, but were not as dependable in distinguishing between S. eriocephala and hybrids. We compared the densities of 11 herbivore species and the infection by a leaf rust pathogen (Melampsora sp.) on the leaves and stems of two parents and the hybrids in the field. We found support for the Additive hypothesis (3 species), the Dominance hypothesis (2 species) and the Hybrid Susceptibility hypothesis (7 species, 6 herbivores and the Melampsora rust). We found no evidence for the Hybrid Resistance hypothesis. Guild membership was not a good predictor of similar responses of species to hybrid versus parental plants. A Canonical Discriminant Function Analysis showed discrete separation of the taxa based on herbivore densities, illustrating different community structures on hybrid and parental plants. This study demonstrates the diversity of responses of phytophages in response to interspecific hybridization.

Morphological and molecular evidence for hybridization and introgression in a willow ( Salix ) hybrid zone

Hybrid zones provide biologists with the opportunity to examine genetic and ecological interactions between differentiated populations. Accurate identification of hybrid genealogies is considered a necessary prerequisite to understanding observed patterns of hybridization–related phenomena. We analysed molecular and morphological data from individuals in a hybrid zone between two species of willows (Salix sericea Marshall and S. eriocephala Michaux) and report the use of randomly amplified polymorphic DNA (RAPD), chloroplast DNA (cpDNA), and ribosomal DNA (rDNA) markers, as well as vegetative morphology and foliar chemistry data to identify individuals in terms of hybrid genealogy and to infer the direction and extent of backcrossing and introgression within the hybrid zone. A novel version of a maximum likelihood estimate approach (developed for this study) was used to calculate hybrid index scores from RAPD marker data; this method produced results similar to those obtained using traditional arithmetic methods. Distribution of rDNA, cpDNA, and chemistry data were examined within the graphical context of RAPD–based hybrid index score histograms and principal component analyses (PCA) on RAPD and morphology data. Seven of the 21 plants classified as S. eriocephala in the field were possible introgressants. Another plant presented an unequivocal example of backcrossed S. sericea chemistry and RAPD markers. Inter– and intraspecific chloroplast diversity found within the hybrid zone suggests both historic introgression (perhaps in a glacial refugium), and contemporary hybridization. Patterns of inheritance and expression within the hybrid zone suggest that morphological characters are often not expressed in a simple additive fashion, and problems associated with both morphological and molecular data are considered.

DNA from bird-dispersed seed and wind-disseminated pollen provides insights into postglacial colonization and population genetic structure of whitebark pine (Pinus albicaulis)

Uniparentally inherited mitochondrial (mt)DNA and chloroplast (cp)DNA microsatellites (cpSSRs) were used to examine population genetic structure and biogeographic patterns of bird‐dispersed seed and wind‐disseminated pollen of whitebark pine (Pinus albicaulis Engelm.). Sampling was conducted from 41 populations throughout the range of the species. Analyses provide evidence for an ancestral haplotype and two derived mtDNA haplotypes with distinct regional distributions. An abrupt contact zone between mtDNA haplotypes in the Cascade Range suggests postglacial biogeographic movements. Among three cpSSR loci, 42 haplotypes were detected within 28 cpSSR sample populations that were aggregated into six regions. Analysis of molecular variance (amova) was used to determine the hierarchical genetic structure of cpSSRs. amova and population pairwise comparisons (FST) of cpSSR, and geographical distribution of mtDNA haplotypes provide insights into historical changes in biogeography. The genetic data suggest that whitebark pine has been intimately tied to climatic change and associated glaciation, which has led to range movements facilitated by seed dispersal by Clark’s nutcracker (Nucifraga columbiana Wilson). The two hypotheses proposed to explain the genetic structure are: (i) a northward expansion into Canada and the northern Cascades in the early Holocene; and (ii) historical gene flow between Idaho and the Oregon Cascades when more continuous habitat existed in Central Oregon during the late Pleistocene. Genetic structure and insights gained from historical seed movements provide a basis on which to develop recovery plans for a species that is at risk from multiple threats.

Multiple independent formations of Tragopogon tetraploids (Asteraceae): evidence from RAPD markers

Polyploidy is widely recognized as a significant force leading to the formation of new plant species. Estimates of the number of angiosperm species with polyploid origins are as high as ≈ 50%; however, in spite of this prevalence, many aspects of polyploid evolution remain poorly understood. Recent studies have suggested that recurrent origins of polyploid species are the rule rather than the exception. The present study is one of only a few designed to quantify the number of independent origins of a polyploid species. The two tetraploid species Tragopogon mirus and T. miscellus (Asteraceae) arose within the past 50 years in the Palouse region of eastern Washington and adjacent northern Idaho. Previous work using morphology, cpDNA and rDNA restriction site analyses, allozymes, cytology, and flavonoid chemistry established that T. mirus had arisen at least five times, and T. miscellus at least twice, on the Palouse. To assess the frequency of multiple origins of these species more rigorously, seven populations of T. mirus and three populations of T. miscellus that were indistinguishable based on previous markers were surveyed using random amplified polymorphic DNA (RAPD) markers; populations of the diploid progenitor species from the same sites were also analysed. Each tetraploid population had a unique RAPD marker profile, suggesting that each population surveyed originated independently of the other populations in the region. Only two of the tetraploid populations combined the RAPD marker profiles of the diploid progenitors occurring at the same site. Both polyploid species, whose ranges and numbers have greatly increased since their formation in the early part of the twentieth century, have formed repeatedly on a local geographical scale and during a short time frame. Furthermore, each tetraploid species is spreading not primarily by dispersal of propagules from a single population of origin, but through repeated, independent polyploidization events that recreate the polyploid taxa.

MOLECULAR EVIDENCE FOR POLYPLOID ORIGINS IN SAXIFRAGA (SAXIFRAGACEAE) : THE NARROW ARCTIC ENDEMIC S. SVALBARDENSIS AND ITS WIDESPREAD ALLIES

The recently described polyploid Saxifraga svalbardensis is endemic to the arctic archipelago of Svalbard. We investigated relationships among four closely related species of Saxifraga in Svalbard and tested three previously proposed hypotheses for the origin of S. svalbardensis: (1) differentiation from the morphologically and chromosomally variable polyploid S. cernua; (2) hybridization between the diploid S. hyperborea and S. cernua; and (3) hybridization between the tetraploid S. rivularis and S. cernua. Fifteen populations were analyzed using random amplified polymorphic DNAs (RAPDs) and nucleotide sequences of the chloroplast gene matK and the internal transcribed spacers (ITS) of nuclear ribosomal DNA (rDNA). RAPD and matK data suggest that S. svalbardensis has originated from a hybrid with S. rivularis as the maternal parent and S. cernua as the paternal parent, possibly a single time, whereas ITS data could not be used to discriminate among the hypotheses. The data also suggest that the diploid S. hyperborea is a progenitor of the tetraploid S. rivularis. The four populations examined of S. svalbardensis were virtually identical for RAPD and ITS markers, whereas S. cernua showed high levels of variation, suggesting that the latter polyploid either has formed recurrently or has undergone considerable differentiation since its origin.

Evolutionary patterns and processes in Salix sect. Longifoliae: evidence from chloroplast DNA

Chloroplast DNA (cpDNA) restriction site analysis was used to explore evolutionary relationships within a taxonomically difficult group of diploid willows comprising seven taxa (Salix sect. Longifoliae). Forty-five population samples from throughout North America were analyzed with up to 23 restriction endonucleases. Clades revealed by cpDNA data are generally congruent with groups delineated by enzyme electrophoresis. Both sets of molecular data disagree with traditional taxonomic concepts. Morphological characters such as stigma length and pubescence, traditionally viewed as taxonomically important, are not correlated with major genetic divisions in the section. High levels of intraspecific cpDNA variability were detected, mostly representing geographic pat- terns of divergence. Salix interior, for example, contains six distinct chloroplast genomes, and has a mean sequence divergence of 0.09%. Conversely, S. sessilifolia, S. fluviatilis, and populations of S. exigua from the Pacific Northwest all possess identical chloroplast genomes. Hybridization appears to be an important process in the section. Five cases of hybridization and subsequent cytoplasmic capture were documented. Chloroplast DNA data suggest a southern origin for the section and a northward radiation into geographically isolated montane and arid lowland habitats. The high incidence of cytoplasmic capture and intraspecific cpDNA variability detected in this study under- scores the importance of analyzing multiple cpDNA samples of a taxon in combination with nuclear markers in molecular systematic studies.

Interspecific and temporal variation in herbivore responses to hybrid willows

We studied herbivory of two species of willows (Salix sericea and S. eriocephala) and their interspecific hybrids to test alternative hypotheses concerning the effects of hybridization on plant resistance. Individually marked plants were identified using morphological traits in the field and random amplified polymorphic DNA (RAPD) band analysis was used to verify the genetic status of many parental and hybrid plants. The desities of 12 herbivore species on plants in the field were compared between two parents and their F2-type hybrids. We found about equal support for the additive, dominance, and hybrid susceptibility hypotheses over 4 years. In one year, one species supported the hybrid resistance hypothesis. Guild membership was not a good predictor of similar responses of species to hybrid versus parental plants. There were marked differences in support for particular hypotheses among years for four herbivore species. This study demonstrates the diversity of responses of phytophages in response to interspecific hybridization, and indicates that year-to-year variation in relative resistance of hybrid plants can be important.

Pollen dispersal and genetic structure of the tropical tree Dipteryx panamensis in a fragmented Costa Rican landscape

In the face of widespread deforestation, the conservation of rainforest trees relies increasingly on their ability to maintain reproductive processes in fragmented landscapes. Here, we analysed nine microsatellite loci for 218 adults and 325 progeny of the tree Dipteryx panamensis in Costa Rica. Pollen dispersal distances, genetic diversity, genetic structure and spatial autocorrelation were determined for populations in four habitats: continuous forest, forest fragments, pastures adjacent to fragments and isolated pastures. We predicted longer but less frequent pollen movements among increasingly isolated trees. This pattern would lead to lower outcrossing rates for pasture trees, as well as lower genetic diversity and increased structure and spatial autocorrelation among their progeny. Results generally followed these expectations, with the shortest pollen dispersal among continuous forest trees (240 m), moderate distances for fragment (343 m) and adjacent pasture (317 m) populations, and distances of up to 2.3 km in isolated pastures (mean: 557 m). Variance around pollen dispersal estimates also increased with fragmentation, suggesting altered pollination conditions. Outcrossing rates were lower for pasture trees and we found greater spatial autocorrelation and genetic structure among their progeny, as well as a trend towards lower heterozygosity. Paternal reproductive dominance, the pollen contributions from individual fathers, did not vary among habitats, but we did document asymmetric pollen flow between pasture and adjacent fragment populations. We conclude that long‐distance pollen dispersal helps maintain gene flow for D. panamensis in this fragmented landscape, but pasture and isolated pasture populations are still at risk of long‐term genetic erosion.

Genetic Variation in Resistance of Hybrid Willows to Herbivores

We measured herbivore abundance on two species of willows (Salix sericea and S. eriocephala) and their interspecific hybrids in two common garden experiments to test alternative hypotheses concerning the genetic effects of hybridization on plant resistance. The first experiment used potted clones of naturally occurring parental and hybrid plants. Purity of parental plants was determined using 20 RAPD markers (pure S. eriocephala had a score of 1.0 and pure S. sericea had a score of 0). The cloned hybrid plants had hybrid scores between 0.40 and 0.60, indicating they were intermediate hybrids and probably had recombinant (F 2 -type) genotypes. The second experiment used potted one-year-old interspecific F 1 progeny and intraspecific progeny from crosses between genetically pure parents. We counted numbers of herbivores on plants after exposure to attack for about 2 months. Herbivore abundance data were analyzed with ANOVA and a priori contrasts were used to test fit of the data to 5 hypotheses (No Difference, Additive, Dominance, Hybrid Susceptibility, or Hybrid Resistance). We found genetic effects of hybridization on susceptibility to some herbivores in both experiments, with the additive, dominance, and hybrid susceptibility hypotheses supported by different herbivore species. There was also significant genetic variation among clones and genotypes within hybrid and parental taxa for herbivore resistance. For several herbivore species we found no difference in the susceptibility among hybrids and parents, which indicates genetic effects of hybridization on resistance were absent for these species. Differences between the hypotheses supported from this study and data from four years of prior field studies suggest that environmental variation is an important factor affecting relative hybrid and parental susceptibility in the field.