PHYLOGEOGRAPHY OF SPOTTED OWL (STRIX OCCIDENTALIS) POPULATIONS BASED ON MITOCHONDRIAL DNA SEQUENCES: GENE FLOW, GENETIC STRUCTURE, AND A NOVEL BIOGEOGRAPHIC PATTERN

Genomic Variation and Recent Population Histories of Spotted (Strix occidentalis) and Barred (Strix varia) Owls

Spotted owls (SOs, Strix occidentalis) are a flagship species inhabiting old-growth forests in western North America. In recent decades, their populations have declined due to ongoing reductions in suitable habitat caused by logging, wildfires, and competition with the congeneric barred owl (BO, Strix varia). The northern spotted owl (S. o. caurina) has been listed as "threatened" under the Endangered Species Act since 1990. Here, we use an updated SO genome assembly along with 51 high-coverage whole-genome sequences to examine population structure, hybridization, and recent changes in population size in SO and BO. We found that potential hybrids identified from intermediate plumage morphology were a mixture of pure BO, F1 hybrids, and F1 × BO backcrosses. Also, although SO underwent a population bottleneck around the time of the Pleistocene-Holocene transition, their population sizes rebounded and show no evidence of any historical (i.e., 100-10,000 years ago) population decline. This suggests that the current decrease in SO abundance is due to events in the past century. Finally, we estimate that western and eastern BOs have been genetically separated for thousands of years, instead of the previously assumed recent (i.e., <150 years) divergence. Although this result is surprising, it is unclear where the ancestors of western BO lived after the separation. In particular, although BO may have colonized western North America much earlier than the first recorded observations, it is also possible that the estimated divergence time reflects unsampled BO population structure within central or eastern North America.

Habitat Fragmentation Reduces Genetic Diversity and Connectivity of the Mexican Spotted Owl: A Simulation Study Using Empirical Resistance Models

We evaluated how differences between two empirical resistance models for the same geographic area affected predictions of gene flow processes and genetic diversity for the Mexican spotted owl (Strix occidentalis lucida). The two resistance models represented the landscape under low- and high-fragmentation parameters. Under low fragmentation, the landscape had larger but highly concentrated habitat patches, whereas under high fragmentation, the landscape had smaller habitat patches that scattered across a broader area. Overall habitat amount differed little between resistance models. We tested eight scenarios reflecting a factorial design of three factors: resistance model (low vs. high fragmentation), isolation hypothesis (isolation-by-distance, IBD, vs. isolation-by-resistance, IBR), and dispersal limit of species (200 km vs. 300 km). Higher dispersal limit generally had a positive but small influence on genetic diversity. Genetic distance increased with both geographic distance and landscape resistance, but landscape resistance displayed a stronger influence. Connectivity was positively related to genetic diversity under IBR but was less important under IBD. Fragmentation had a strong negative influence on the spatial patterns of genetic diversity and effective population size (Ns). Despite habitats being more concentrated and less widely distributed, the low-fragmentation landscape had greater genetic diversity than the high-fragmentation landscape, suggesting that highly concentrated but larger habitat patches may provide a genetic refuge for the Mexican spotted owl.

Complete mitochondrial genome sequences of the northern spotted owl (Strix occidentalis caurina) and the barred owl (Strix varia; Aves: Strigiformes: Strigidae) confirm the presence of a duplicated control region

We report here the successful assembly of the complete mitochondrial genomes of the northern spotted owl (Strix occidentalis caurina) and the barred owl (S. varia). We utilized sequence data from two sequencing methodologies, Illumina paired-end sequence data with insert lengths ranging from approximately 250 nucleotides (nt) to 9,600 nt and read lengths from 100–375 nt and Sanger-derived sequences. We employed multiple assemblers and alignment methods to generate the final assemblies. The circular genomes of S. o. caurina and S. varia are comprised of 19,948 nt and 18,975 nt, respectively. Both code for two rRNAs, twenty-two tRNAs, and thirteen polypeptides. They both have duplicated control region sequences with complex repeat structures. We were not able to assemble the control regions solely using Illumina paired-end sequence data. By fully spanning the control regions, Sanger-derived sequences enabled accurate and complete assembly of these mitochondrial genomes. These are the first complete mitochondrial genome sequences of owls (Aves: Strigiformes) possessing duplicated control regions. We searched the nuclear genome of S. o. caurina for copies of mitochondrial genes and found at least nine separate stretches of nuclear copies of gene sequences originating in the mitochondrial genome (Numts). The Numts ranged from 226–19,522 nt in length and included copies of all mitochondrial genes except tRNA^Pro, ND6, and tRNA^Glu. Strix occidentalis caurina and S. varia exhibited an average of 10.74% (8.68% uncorrected p-distance) divergence across the non-tRNA mitochondrial genes.

Northern Spotted Owl (Strix occidentalis caurina) Genome: Divergence with the Barred Owl (Strix varia) and Characterization of Light-Associated Genes

We report here the assembly of a northern spotted owl (Strix occidentalis caurina) genome. We generated Illumina paired-end sequence data at 90× coverage using nine libraries with insert lengths ranging from ∼250 to 9,600 nt and read lengths from 100 to 375 nt. The genome assembly is comprised of 8,108 scaffolds totaling 1.26 × 109 nt in length with an N50 length of 3.98 × 106 nt. We calculated the genome-wide fixation index (FST) of S. o. caurina with the closely related barred owl (Strix varia) as 0.819. We examined 19 genes that encode proteins with light-dependent functions in our genome assembly as well as in that of the barn owl (Tyto alba). We present genomic evidence for loss of three of these in S. o. caurina and four in T. alba. We suggest that most light-associated gene functions have been maintained in owls and their loss has not proceeded to the same extent as in other dim-light-adapted vertebrates.

Genetic differentiation and inferred dynamics of a hybrid zone between Northern Spotted Owls (Strix occidentalis caurina) and California Spotted Owls (S. o. occidentalis) in northern California

Genetic differentiation among Spotted Owl (Strix occidentalis) subspecies has been established in prior studies. These investigations also provided evidence for introgression and hybridization among taxa but were limited by a lack of samples from geographic regions where subspecies came into close contact. We analyzed new sets of samples from Northern Spotted Owls (NSO: S. o. caurina) and California Spotted Owls (CSO: S. o. occidentalis) in northern California using mitochondrial DNA sequences (mtDNA) and 10 nuclear microsatellite loci to obtain a clearer depiction of genetic differentiation and hybridization in the region. Our analyses revealed that a NSO population close to the northern edge of the CSO range in northern California (the NSO Contact Zone population) is highly differentiated relative to other NSO populations throughout the remainder of their range. Phylogenetic analyses identified a unique lineage of mtDNA in the NSO Contact Zone, and Bayesian clustering analyses of the microsatellite data identified the Contact Zone as a third distinct population that is differentiated from CSO and NSO found in the remainder of the subspecies' range. Hybridization between NSO and CSO was readily detected in the NSO Contact Zone, with over 50% of individuals showing evidence of hybrid ancestry. Hybridization was also identified among 14% of CSO samples, which were dispersed across the subspecies' range in the Sierra Nevada Mountains. The asymmetry of hybridization suggested that the hybrid zone may be dynamic and moving. Although evidence of hybridization existed, we identified no F1 generation hybrid individuals. We instead found evidence for F2 or backcrossed individuals among our samples. The absence of F1 hybrids may indicate that (1) our 10 microsatellites were unable to distinguish hybrid types, (2) primary interactions between subspecies are occurring elsewhere on the landscape, or (3) dispersal between the subspecies' ranges is reduced relative to historical levels, potentially as a consequence of recent regional fires.

Prioritizing populations for conservation using phylogenetic networks

In the face of inevitable future losses to biodiversity, ranking species by conservation priority seems more than prudent. Setting conservation priorities within species (i.e., at the population level) may be critical as species ranges become fragmented and connectivity declines. However, existing approaches to prioritization (e.g., scoring organisms by their expected genetic contribution) are based on phylogenetic trees, which may be poor representations of differentiation below the species level. In this paper we extend evolutionary isolation indices used in conservation planning from phylogenetic trees to phylogenetic networks. Such networks better represent population differentiation, and our extension allows populations to be ranked in order of their expected contribution to the set. We illustrate the approach using data from two imperiled species: the spotted owl Strix occidentalis in North America and the mountain pygmy-possum Burramys parvus in Australia. Using previously published mitochondrial and microsatellite data, we construct phylogenetic networks and score each population by its relative genetic distinctiveness. In both cases, our phylogenetic networks capture the geographic structure of each species: geographically peripheral populations harbor less-redundant genetic information, increasing their conservation rankings. We note that our approach can be used with all conservation-relevant distances (e.g., those based on whole-genome, ecological, or adaptive variation) and suggest it be added to the assortment of tools available to wildlife managers for allocating effort among threatened populations.

Introgression and dispersal among spotted owl (Strix occidentalis) subspecies

Population genetics plays an increasingly important role in the conservation and management of declining species, particularly for defining taxonomic units. Subspecies are recognized by several conservation organizations and countries and receive legal protection under the US Endangered Species Act (ESA). Two subspecies of spotted owls, northern (Strix occidentalis caurina) and Mexican (S. o. lucida) spotted owls, are ESA-listed as threatened, but the California (S. o. occidentalis) spotted owl is not listed. Thus, determining the boundaries of these subspecies is critical for effective enforcement of the ESA. We tested the validity of previously recognized spotted owl subspecies by analysing 394 spotted owls at 10 microsatellite loci. We also tested whether northern and California spotted owls hybridize as suggested by previous mitochondrial DNA studies. Our results supported current recognition of three subspecies. We also found bi-directional hybridization and dispersal between northern and California spotted owls centered in southern Oregon and northern California. Surprisingly, we also detected introgression of Mexican spotted owls into the range of northern spotted owls, primarily in the northern part of the subspecies' range in Washington, indicating long-distance dispersal of Mexican spotted owls. We conclude with a discussion of the conservation implications of our study.