Vicki L. Friesen

Mechanisms of population differentiation in seabirds

Despite recent advances in population genetic theory and empirical research, the extent of genetic differentiation among natural populations of animals remains difficult to predict. We reviewed studies of geographic variation in mitochondrial DNA in seabirds to test the importance of various factors in generating population genetic and phylogeographic structure. The extent of population genetic and phylogeographic structure varies extensively among species. Species fragmented by land or ice invariably exhibit population genetic structure and most also have phylogeographic structure. However, many populations (26 of 37) display genetic structure in the absence of land, suggesting that other barriers to gene flow exist. In these populations, the extent of genetic structure is best explained by nonbreeding distribution: almost all species with two or more population‐specific nonbreeding areas (or seasons) have phylogeographic structure, and all species that are resident at or near breeding colonies year‐round have population genetic structure. Geographic distance between colonies and foraging range appeared to have a weak influence on the extent of population genetic structure, but little evidence was found for an effect of colony dispersion or population bottlenecks. In two species (Galapagos petrel, Pterodroma phaeopygia, and Xantuss murrelet, Synthliboramphus hypoleucus), population genetic structure, and even phylogeographic structure, exist in the absence of any recognizable physical or nonphysical barrier, suggesting that other selective or behavioural processes such as philopatry may limit gene flow. Retained ancestral variation may be masking barriers to dispersal in some species, especially at high latitudes. Allopatric speciation undoubtedly occurs in this group, but reproductive isolation also appears to have evolved through founder‐induced speciation, and there is strong evidence that parapatric and sympatric speciation occur. While many questions remain unanswered, results of the present review should aid conservation efforts by enabling managers to predict the extent of population differentiation in species that have not yet been studied using molecular markers, and, thus, enable the identification of management units and evolutionary significant units for conservation.

Intron variation in marbled murrelets detected using analyses of single-stranded conformational polymorphisms.

Combination of the targeted amplification of nuclear introns and the analysis of single‐stranded conformational polymorphisms has the potential to provide an inexpensive, rapid, versatile and sensitive genetic assay for evolutionary studies and conservation. We are developing primers and protocols to analyse nuclear introns in vertebrates, and are testing them in a population genetic study of marbled murrelets Brachyramphus marmoratus. Here we present protocols and results for introns for aldolase B, α‐enolase, glyceraldehyde‐3‐phosphate dehydrogenase and lamin A. Results suggest that this approach presents a potentially powerful method for detecting genetic variation within and among local populations and species of animals: (i) a variety of genes can be surveyed, including genes of special interest such as those involved in disease resistance; (ii) assays are rapid and relatively inexpensive; (iii) large numbers of genes can be assayed, enabling accurate estimation of variation in the total genome; (iv) almost any mutation can be detected in the genes amplified; (v) the exact nature of variation can be investigated by sequence analysis if desired; (vi) statistical methods previously developed for proteins and/or sequence data can be used; (vii) protocols can be easily transferred to other species and other laboratories; and (viii) assays can be performed on old or degraded samples, blood or museum skins, so that animals need not be killed. Results of analyses for murrelets support earlier evidence that North American and Asiatic subspecies represent reproductively isolated species, and that genetic differences exist among murrelets from different sites within North America.

A TEST OF THE GLACIAL REFUGIUM HYPOTHESIS USING PATTERNS OF MITOCHONDRIAL AND NUCLEAR DNA SEQUENCE VARIATION IN ROCK PTARMIGAN (LAGOPUS MUTUS)

The glacial refugium hypothesis (GRH) proposes that glaciers promoted differentiation and generation of intraspecific diversity by isolating populations in ice‐free refugia. We tested three predictions of this hypothesis for the evolutionary divergence of rock ptarmigan (Lagopus mutus) during the Wisconsin glaciation of the late Pleistocene. To do this, we examined subspecies distributions, population genetic structure, and phylogenetic relationships in 26 populations across North America and the Bering Sea region. First, we analyzed sequence variation in the mitochondrial control region, in a nuclear intron (Gapdh), and in an internal transcribed spacer (ITS1). Control region sequences of 154 rock ptarmigan revealed strong population and phylogeographic structure. Variation in intron sequences of 114 rock ptarmigan also revealed significant population structure compatible with results for the control region. Rock ptarmigan were invariant for ITS1. Second, we show that five known Nearctic refugia and an Icelandic refugium are concordant with the current distribution of morphologically distinct subspecies; five of these six refugia are geographically concordant with the distribution of closely related control region haplotypes. Third, our estimates of the time since phylogenetic lineages diverged predated the last glacial maximum for all but two lineages. In addition, all lines of evidence suggest that two unknown refugia in the Bering Sea region supported rock ptarmigan during the Wisconsin glaciation. Overall, our results are most consistent with the hypothesis that isolated populations of rock ptarmigan diverged in multiple refugia during the Wisconsin and that geographic variation reflects patterns of recolonization of the Nearctic after the ice receded. The GRH may therefore offer the most plausible explanation for similar biogeographic patterns in a variety of Nearctic vertebrates.

Polymerase chain reaction (PCR) primers for the amplification of five nuclear introns in vertebrates

[Extract] Advancements in evolutionary genetics, as well as the conservation of biodiversity, increasingly require direct analyses of sequence variation in nuclear DNA. Recent studies indicate that nuclear introns have variabilities useful for both phylogenetics and population genetics (reviewed in Friesen 2000); however, use of introns is currently limited by the paucity of polymerase chain reaction (PCR) primers that have been demonstrated to have broad taxonomic utility (although several primers with less general or uncertain utilities have been published; reviewed in Friesen 2000). We have designed 30 general PCR primers for nuclear introns for vertebrates.

MECHANISMS OF POPULATION DIFFERENTIATION IN MARBLED MURRELETS: HISTORICAL VERSUS CONTEMPORARY PROCESSES

Abstract.— Mechanisms of population differentiation in highly vagile species such as seabirds are poorly understood. Previous studies of marbled murrelets (Brachyramphus marmoratus; Charadriiformes: Alcidae) found significant population genetic structure, but could not determine whether this structure is due to historical vicariance (e.g., due to Pleistocene glaciers), isolation by distance, drift or selection in peripheral populations, or nesting habitat selection. To discriminate among these possibilities, we analyzed sequence variation in nine nuclear introns from 120 marbled murrelets sampled from British Columbia to the western Aleutian Islands. Mismatch distributions indicated that murrelets underwent at least one population expansion during the Pleistocene and probably are not in genetic equilibrium. Maximum‐likelihood analysis of allele frequencies suggested that murrelets from “mainland” sites (from the Alaskan Peninsula east) are genetically different from those in the Aleutians and that these two lineages diverged prior to the last glaciation. Analyses of molecular variance, as well as estimates of gene flow derived using coalescent theory, indicate that population genetic structure is best explained by peripheral isolation of murrelets in the Aleutian Islands, rather than by selection associated with different nesting habitats. No isolation‐by‐distance effects could be detected. Our results are consistent with a rapid expansion of murrelets from a single refugium during the early–mid Pleistocene, subsequent isolation and divergence in two or more refugia during the final Pleistocene glacial advance, and secondary contact following retreat of the ice sheets. Population genetic structure now appears to be maintained by distance effects combined with small populations and a highly fragmented habitat in the Aleutian Islands.

POLYTOMIES AND THE POWER OF PHYLOGENETIC INFERENCE

Although phylogenetic hypotheses can provide insights into mechanisms of evolution, their utility is limited by our inability to differentiate simultaneous speciation events (hard polytomies) from rapid cladogenesis (soft polytomies). In the present paper, we tested the potential for statistical power analysis to differentiate between hard and soft polytomies in molecular phytogenies. Classical power analysis typically is used a priori to determine the sample size required to detect a particular effect size at a particular level of significance (a) with a certain power (1 – β). A posteriori, power analysis is used to infer whether failure to reject a null hypothesis results from lack of an effect or from insufficient data (i.e., low power). We adapted this approach to molecular data to infer whether polytomies result from simultaneous branching events or from insufficient sequence information. We then used this approach to determine the amount of sequence data (sample size) required to detect a positive branch length (effect size). A worked example is provided based on the auklets (Charadriiformes: Alcidae), a group of seabirds among which relationships are represented by a polytomy, despite analyses of over 3000 bp of sequence data. We demonstrate the calculation of effect sizes and sample sizes from sequence data using a normal curve test for difference of a proportion from an expected value and a t‐test for a difference of a mean from an expected value. Power analyses indicated that the data for the auklets should be sufficient to differentiate speciation events that occurred at least 100,000 yr apart (the duration of the shortest glacial and interglacial events of the Pleistocene), 2.6 million years ago.

Complex hybridization dynamics between golden‐winged and blue‐winged warblers (Vermivora chrysoptera and Vermivora pinus) revealed by AFLP, microsatellite, intron and mtDNA markers

Blue‐winged (Vermivora pinus) and golden‐winged warblers (Vermivora chrysoptera) have an extensive mosaic hybrid zone in eastern North America. Over the past century, the general trajectory has been a rapid replacement of chrysoptera by pinus in a broad, northwardly moving area of contact. Previous mtDNA‐based studies on these species’ hybridization dynamics have yielded variable results: asymmetric and rapid introgression from pinus into chrysoptera in some areas and bidirectional maternal gene flow in others. To further explore the hybridization genetics of this otherwise well‐studied complex, we surveyed variation in three nuclear DNA marker types — microsatellites, introns, and a panel of amplified fragment length polymorphisms (AFLPs) — with the goal of generating a multilocus assay of hybrid introgression. All markers were first tested on birds from phenotypically and mitochondrially pure parental‐type populations from outside the hybrid zone. Searches for private alleles and assignment test approaches found no combination of microsatellite or intron markers that could separate the parental populations, but seven AFLP characters exhibited significant frequency differences among them. We then used the AFLP markers to examine the extent and pattern of introgression in a population where pinus‐phenotype individuals have recently invaded a region that previously supported only a chrysoptera‐phenotype population. Despite the low frequency of phenotypic hybrids at this location, the AFLP data suggest that almost a third of the phenotypically pure chrysoptera have introgressed genotypes, indicating the presence of substantial cryptic hybridization in the history of this species. The evidence for extensive cryptic introgression, combined with the lack of differentiation at other nuclear loci, cautions against hybrid assessments based on single markers or on phenotypic traits that are likely to be determined by a small number of loci. Considered in concert, these results from four classes of molecular markers indicate that pinus and chrysoptera are surprisingly weakly differentiated and that far fewer genetically ‘pure’ populations of chrysoptera may exist than previously assumed, two findings with broad implications for the conservation of this rapidly declining taxon.

Glacial vicariance and historical biogeography of rock ptarmigan (Lagopus mutus) in the Bering region

In this paper, we address alternative hypotheses for the evolution of subspecies of rock ptarmigan (Lagopus mutus) endemic to the Aleutian Archipelago. To do this we examined patterns of genetic differentiation among populations of rock ptarmigan in the Aleutian Islands and parts of both Alaska and Siberia. Variation in mitochondrial control region sequences of 105 rock ptarmigan from 10 subspecies within the Bering region revealed three major phylogenetic lineages, two of which are endemic to the Aleutian Islands. Accordingly, haplotype and nucleotide diversities of rock ptarmigan within the archipelago are much higher than within mainland Alaska or Siberia. For Aleutian rock ptarmigan, analyses of molecular variance indicated significant genetic structuring and low estimates of gene flow among populations, despite small interisland distances within the archipelago. However, isolation by distance did not describe the pattern of gene flow or differentiation at this scale. Our estimates of divergence times of lineages suggest that Aleutian rock ptarmigan became isolated prior to the most recent Pleistocene glaciation event (late Wisconsin Stade) and that current patterns of genetic variation reflect the postglacial redistribution of divergent lineages and subsequent limited gene flow. In addition, genetic divergence among lineages was concordant with the distribution of plumage types among subspecies. The patterns of genetic variation described here for rock ptarmigan provide evidence for the role of glacial vicariance in contributing to genetic diversity within this and other Bering region species.

The Isthmus of Panama: a major physical barrier to gene flow in a highly mobile pantropical seabird

To further test the hypothesis that the Isthmus of Panama is a major barrier to gene flow in pantropical seabirds, we applied phylogeographic methods to mitochondrial control sequence variation in masked booby (Sula dactylatra) populations on either side of the Isthmus of Panama and the southern tip of Africa. In accord with Steeves et al. (2003) , we found that all Caribbean masked boobies with the ‘secondary contact’ cytochrome b haplotype (m‐B) shared a control region haplotype (Sd_100), which grouped with Indian–Pacific haplotypes and not Caribbean–Atlantic haplotypes. In addition, Sd_100 was more closely related to control region haplotypes in the Indian Ocean than in the Pacific. We also found that the ‘secondary contact’ birds diverged more recently from extant populations in the Indian Ocean than in the Pacific. Thus, it appears that these masked boobies did not breach the Isthmus of Panama. Rather, birds likely dispersed around the southern tip of Africa during favourable oceanographic conditions in the Pleistocene.

Appraisal of the consequences of the DDT-induced bottleneck on the level and geographic distribution of neutral genetic variation in Canadian peregrine falcons, Falco peregrinus

Peregrine falcon populations underwent devastating declines in the mid‐20th century due to the bioaccumulation of organochlorine contaminants, becoming essentially extirpated east of the Great Plains and significantly reduced elsewhere in North America. Extensive re‐introduction programs and restrictions on pesticide use in Canada and the United States have returned many populations to predecline sizes. A proper population genetic appraisal of the consequences of this decline requires an appropriate context defined by (i) meaningful demographic entities; and (ii) suitable reference populations. Here we explore the validity of currently recognized subspecies designations using data from the mitochondrial control region and 11 polymorphic microsatellite loci taken from 184 contemporary individuals from across the breeding range, and compare patterns of population genetic structure with historical patterns inferred from 95 museum specimens. Of the three North American subspecies, the west coast marine subspecies Falco peregrinus pealei is well differentiated genetically in both time periods using nuclear loci. In contrast, the partitioning of continental Falco peregrinus anatum and arctic Falco peregrinus tundrius subspecies is not substantiated, as individuals from these subspecies are historically indistinguishable genetically. Bayesian clustering analyses demonstrate that contemporary genetic differentiation between these two subspecies is mainly due to changes within F. p. anatum (specifically the southern F. p. anatum populations). Despite expectations and a variety of tests, no genetic bottleneck signature is found in the identified populations; in fact, many contemporary indices of diversity are higher than historical values. These results are rationalized by the promptness of the recovery and the possible introduction of new genetic material.