Low levels of genetic differentiation and structure in red fox populations in Eastern Canada and implications for Arctic fox rabies propagation potential.

Rabies is a lethal zoonosis present in most parts of the world which can be transmitted to humans through the bite from an infected mammalian reservoir host. The Arctic rabies virus variant (ARVV) persists mainly in populations of Arctic foxes (Vulpes lagopus), and to a lesser extent in red fox populations (Vulpes vulpes). Red foxes are thought to be responsible for sporadic southward movement waves of the ARVV outside the enzootic area of northern Canada. In this study, we wanted to investigate whether red foxes displayed notable levels of genetic structure across the Quebec-Labrador Peninsula, which includes portions of the provinces of Quebec and Newfoundland-Labrador in Canada, and is a region with a history of southward ARVV movement waves. We combined two datasets that were collected and genotyped using different protocols, totalling 675 red fox individuals across the whole region and genotyped across 13 microsatellite markers. We found two genetic clusters across the region, reflecting a latitudinal gradient, and characterized by low genetic differentiation. We also observed weak but significant isolation by distance, which seems to be marginally more important for females than for males. These findings suggest a general lack of resistance to movement in red fox populations across the Quebec-Labrador Peninsula, regardless of sex. Implications of these findings include additional support for the hypothesis of long-distance southward ARVV propagation through its red fox reservoir host.

Relationships between fox populations and rabies virus spread in northern Canada.

Rabies spreads in both Arctic (Vulpes lagopus) and red foxes (Vulpes vulpes) throughout the Canadian Arctic but limited wildlife disease surveillance, due to the extensive landmass of the Canadian north and its small widely scattered human population, undermines our knowledge of disease transmission patterns. This study has explored genetic population structure in both the rabies virus and its fox hosts to better understand factors that impact rabies spread. Phylogenetic analysis of 278 samples of the Arctic lineage of rabies virus recovered over 40 years identified four sub-lineages, A1 to A4. The A1 lineage has been restricted to southern regions of the Canadian province of Ontario. The A2 lineage, which predominates in Siberia, has also spread to northern Alaska while the A4 lineage was recovered from southern Alaska only. The A3 sub-lineage, which was also found in northern Alaska, has been responsible for virtually all cases across northern Canada and Greenland, where it further differentiated into 18 groups which have systematically evolved from a common predecessor since 1975. In areas of Arctic and red fox sympatry, viral groups appear to circulate in both hosts, but both mitochondrial DNA control region sequences and 9-locus microsatellite genotypes revealed contrasting phylogeographic patterns for the two fox species. Among 157 Arctic foxes, 33 mitochondrial control region haplotypes were identified but little genetic structure differentiating localities was detected. Among 162 red foxes, 18 control region haplotypes delineated three groups which discriminated among the Churchill region of Manitoba, northern Quebec and Labrador populations, and the coastal Labrador locality of Cartwright. Microsatellite analyses demonstrated some genetic heterogeneity among sampling localities of Arctic foxes but no obvious pattern, while two or three clusters of red foxes suggested some admixture between the Churchill and Quebec-Labrador regions but uniqueness of the Cartwright group. The limited population structure of Arctic foxes is consistent with the rapid spread of rabies virus subtypes throughout the north, while red fox population substructure suggests that disease spread in this host moves most readily down certain independent corridors such as the northeastern coast of Canada and the central interior. Interestingly the evidence suggests that these red fox populations have limited capacity to maintain the virus over the long term, but they may contribute to viral persistence in areas of red and Arctic fox sympatry.

Multi-host dispersal of known and novel carnivore amdoparvoviruses.

Amdoparvoviruses (family Parvoviridae) are ssDNA viruses that cause an immune complex-mediated wasting syndrome in carnivores. They are multi-host pathogens and cross-species infection is facilitated by the fact that viral entry is mediated by cellular Fc receptors recognizing antibody-coated viruses. We developed a pan-amdoparvovirus PCR and screened tissue samples from 666 wild carnivores (families Felidae, Canidae, and Mustelidae) from Newfoundland or Labrador (Canada) and molecularly characterized the identified strains. Fifty-four out of 666 (8.1%) animals were amdoparvovirus-positive. Infection rate was the highest in American mink (34/47, 72.3%), followed by foxes (Arctic and red foxes, 13/311, 4.2%), lynx (2/58, 3.5%), and American martens (5/156, 3.4%). No virus was detected in samples from 87 coyotes and 17 ermines. Viruses from Newfoundland were classified as Aleutian mink disease virus (AMDV). Mink harvested near AMDV-affected fur farms had higher prevalence (24/24, 100%) than other mink (10/23, 43.5%; P < 0.001) and their viruses were phylogenetically closely related to those from farms, while most viruses from other mink were in other clades. Strains from three foxes and two lynx were highly related to mink strains. This proves that farms disperse AMDV that subsequently spreads among wild mink (maintenance host) and transmits to other spillover carnivore hosts. In Labrador two novel viruses were identified, Labrador amdoparvovirus 1 (LaAV-1) found in foxes (9/261, 3.5%) and martens (5/156, 3.4%), and LaAV-2 found in one fox (0.4%). LaAV-1 fulfills all requirements to be classified as a novel species. LaAV-1 was most similar to viruses of mink and skunks (AMDV and skunk amdoparvovirus (SKAV)) while LaAV-2 was more closely related to other viruses infecting canids. LaAV-1 capsid proteins were almost indistinguishable from those of AMDV in some regions, suggesting that LaAV-1 could be a virus of mustelids that can infect foxes. While intensive farming practices provide occasions for inter-species transmission in farms, niche overlap or predation could explain cross-species transmission in the wild, but competition among sympatric species reduces the chances of direct contacts, making this an infrequent event. Pan-amdoparvovirus detection methods in wide epidemiological investigations can play a crucial role in defining amdoparvoviral ecology and evolution and discovering novel viruses.

Phylogeographic mitogenomics of Atlantic cod Gadus morhua: Variation in and among trans-Atlantic, trans-Laurentian, Northern cod, and landlocked fjord populations.

The historical phylogeography, biogeography, and ecology of Atlantic cod (Gadus morhua) have been impacted by cyclic Pleistocene glaciations, where drops in sea temperatures led to sequestering of water in ice sheets, emergence of continental shelves, and changes to ocean currents. High-resolution, whole-genome mitogenomic phylogeography can help to elucidate this history. We identified eight major haplogroups among 153 fish from 14 populations by Bayesian, parsimony, and distance methods, including one that extends the species coalescent back to ca. 330 kya. Fish from the Barents and Baltic Seas tend to occur in basal haplogroups versus more recent distribution of fish in the Northwest Atlantic. There was significant differentiation in the majority of trans-Atlantic comparisons (ΦST = .029-.180), but little or none in pairwise comparisons within the Northwest Atlantic of individual populations (ΦST = .000-.060) or defined management stocks (ΦST = .000-.023). Monte Carlo randomization tests of population phylogeography showed significantly nonrandom trans-Atlantic phylogeography versus absence of such structure within various partitions of trans-Laurentian, Northern cod (NAFO 2J3KL) and other management stocks, and Flemish Cap populations. A landlocked meromictic fjord on Baffin Island comprised multiple identical or near-identical mitogenomes in two major polyphyletic clades, and was significantly differentiated from all other populations (ΦST = .153-.340). The phylogeography supports a hypothesis of an eastern origin of genetic diversity ca. 200-250 kya, rapid expansion of a western superhaplogroup comprising four haplogroups ca. 150 kya, and recent postglacial founder populations.

Patterns of MHC-DRB1 polymorphism in a post-glacial island canid, the Newfoundland red fox (Vulpes vulpes deletrix), suggest balancing selection at species and population timescales.

As the only native insular Newfoundland canid between the extinction of the wolf in the 1930s and the recent arrival of coyotes, the red fox (Vulpes vulpes deletrix Bangs 1898) poses interesting questions about genetic distinctiveness and the post-glacial colonization history of the island's depauperate mammalian fauna. Here, we characterized genetic variability at the major histocompatibility complex (MHC) class II DR ß1 domain (DRB1) locus in 28 red foxes from six sampling localities island-wide and compared it with mitochondrial control region (CR) diversity and DRB1 diversity in other canids. Our goals were to describe novel DRB1 alleles in a new canid population and to make inferences about the role of selection in maintaining their diversity. As in numerous studies of vertebrates, we found an order-of-magnitude higher nucleotide diversity at the DRB1 locus compared with the CR and significantly positive nonsynonymous-to-synonymous substitution ratios, indicative of selection in the distant past. Although the evidence is weaker, the Ewens-Watterson test of neutrality and the geographical distribution of variation compared with the CR suggest a role for selection over the evolutionary timescale of populations. We report the first genetic data from the DRB1 locus in the red fox and establish baseline information regarding immunogenetic variation in this island canid population which should inform continued investigations of population demography, adaptive genetic diversity, and wildlife disease in red foxes and related species.

Quantitative Phylogenomics of Within-Species Mitogenome Variation: Monte Carlo and Non-Parametric Analysis of Phylogeographic Structure among Discrete Transatlantic Breeding Areas of Harp Seals (Pagophilus groenlandicus).

Phylogenomic analysis of highly-resolved intraspecific phylogenies obtained from complete mitochondrial DNA genomes has had great success in clarifying relationships within and among human populations, but has found limited application in other wild species. Analytical challenges include assessment of random versus non-random phylogeographic distributions, and quantification of differences in tree topologies among populations. Harp Seals (Pagophilus groenlandicus Erxleben, 1777) have a biogeographic distribution based on four discrete trans-Atlantic breeding and whelping populations located on "fast ice" attached to land in the White Sea, Greenland Sea, the Labrador ice Front, and Southern Gulf of St Lawrence. This East to West distribution provides a set of a priori phylogeographic hypotheses. Outstanding biogeographic questions include the degree of genetic distinctiveness among these populations, in particular between the Greenland Sea and White Sea grounds. We obtained complete coding-region DNA sequences (15,825 bp) for 53 seals. Each seal has a unique mtDNA genome sequence, which differ by 6 ~ 107 substitutions. Six major clades / groups are detectable by parsimony, neighbor-joining, and Bayesian methods, all of which are found in breeding populations on either side of the Atlantic. The species coalescent is at 180 KYA; the most recent clade, which accounts for 66% of the diversity, reflects an expansion during the mid-Wisconsinan glaciation 40~60 KYA. FST is significant only between the White Sea and Greenland Sea or Ice Front populations. Hierarchal AMOVA of 2-, 3-, or 4-island models identifies small but significant ΦSC among populations within groups, but not among groups. A novel Monte-Carlo simulation indicates that the observed distribution of individuals within breeding populations over the phylogenetic tree requires significantly fewer dispersal events than random expectation, consistent with island or a priori East to West 2- or 3-stepping-stone biogeographic models, but not a simple 1-step trans-Atlantic model. Plots of the cumulative pairwise sequence difference curves among seals in each of the four populations provide continuous proxies for phylogenetic diversification within each. Non-parametric Kolmogorov-Smirnov (K-S) tests of maximum pairwise differences between these curves indicates that the Greenland Sea population has a markedly younger phylogenetic structure than either the White Sea population or the two Northwest Atlantic populations, which are of intermediate age and homogeneous structure. The Monte Carlo and K-S assessments provide sensitive quantitative tests of within-species mitogenomic phylogeography. This is the first study to indicate that the White Sea and Greenland Sea populations have different population genetic histories. The analysis supports the hypothesis that Harp Seals comprises three genetically distinguishable breeding populations, in the White Sea, Greenland Sea, and Northwest Atlantic. Implications for an ice-dependent species during ongoing climate change are discussed.

Hereditary hyperplastic gingivitis in North American farmed silver fox (Vulpes vulpes).

Hereditary hyperplastic gingivitis is a progressive growth of gingival tissues in foxes resulting in dental encapsulation. It is an autosomal recessive condition displaying a gender-biased penetrance, with an association with superior fur quality. This disease has been primarily described in European farmed foxes. Here we document its emergence in Canada.

Gingivite hyperplasique héréditaire chez le renard argenté d'élevage d'Amérique du Nord(Vulpes vulpes). La gingivite hyperplasique héréditaire est une croissance progressive des tissus gingivaux chez les renards qui produit une encapsulation dentaire. Il s'agit d'une affection récessive autosomique qui manifeste une pénétration privilégiant un sexe et qui présente une association avec une qualité de fourrure supérieure. Cette maladie a été principalement décrite chez les renards d'élevage européen. Nous documentons ici son émergence au Canada.(Traduit par Isabelle Vallières).

Genome-wide expression analysis of hereditary hyperplastic gingivitis in silver foxes (Vulpes vulpes) using canine microarrays.

Hereditary hyperplastic gingivitis (HHG) is an autosomal recessive condition found predominantly in farmed silver foxes, first documented in Europe in the 1940s. Hereditary gingival fibromatosis (HGF) is an analogous condition occurring in humans. HGF has a heterogeneous aetiology with emphasis placed on the autosomal dominant forms of inheritance for which there are three known loci: HGF1, HGF2, and HGF3. Among these, only one causative mutation has been determined, in the Son of sevenless homolog 1 (SOS1) gene. The goal of this study was to explore potential molecular or cellular mechanisms underlying HHG by analysis of global gene expression patterns from Affymetrix Canine 2.0 microarrays cross-referenced against candidate genes within the human loci. We conclude that the SOS1 gene involved in HGF1 is not significantly up-regulated in HHG. However, the structurally and functionally similar SOS2 gene is up-regulated in affected foxes, and we propose this as a candidate gene for HHG. At HGF2 we identify RASA1 (rat sarcoma viral p21 protein activator 1) as a candidate gene for HHG, as it is up-regulated in affected foxes and is involved in MAPK signalling. From comparison to the genes within the HGF3 locus, we find evidence for a role of androgens in HHG phenotype severity by differential up-regulation of SRD5A2 in HHG-affected foxes. We hypothesize that the putative mutation occurs upstream of RAS in the extracellular signal-regulated kinase component of MAPK signalling.

Complete mitochondrial genomes from four subspecies of common chaffinch (Fringilla coelebs): new inferences about mitochondrial rate heterogeneity, neutral theory, and phylogenetic relationships within the order Passeriformes.

We describe whole mitochondrial genome sequences from four subspecies of the common chaffinch (Fringilla coelebs), and compare them to 31 publicly available mitochondrial genome sequences from other Passeriformes. Rates and patterns of mitochondrial gene evolution are analyzed at different taxonomic levels within this avian order, and evidence is adduced for and against the nearly neutral theory of molecular evolution and the role of positive selection in shaping genetic variation of this small but critical genome. We find evidence of mitochondrial rate heterogeneity in birds as in other vertebrates, likely due to differences in mutational pressure across the genome. Unlike in gadine fish and some of the human mitochondrial work we do not observe strong support for the nearly neutral theory of molecular evolution; instead evidence from molecular clocks, distribution of dN/dS ratios at different levels of the taxonomic hierarchy and in different lineages, McDonald-Kreitman tests within Fringillidae, and site-specific tests of selection within Passeriformes, all point to a role for positive selection, especially for the complex I NADH dehydrogenase genes. The protein-coding mitogenome phylogeny of the order Passeriformes is broadly consistent with previously-reported molecular findings, but provides support for a sister relationship between the superfamilies Muscicapoidea and Passeroidea on a short basal internode of the Passerida where relationships have been difficult to resolve. An unexpected placement of the Paridae (represented by Hume's groundpecker) within the Muscicapoidea was observed. Consistent with other molecular studies the mtDNA phylogeny reveals paraphyly within the Muscicapoidea and a sister relationship of Fringilla with Carduelis rather than Emberiza.

Sequence analysis of three pigmentation genes in the Newfoundland population of Canis latrans links the Golden Retriever Mc1r variant to white coat color in coyotes.

Three genes, Mc1r, Agouti, and CBD103, interact in a type-switching process that controls much of the pigmentation variation observed in mammals. A deletion in the CBD103 gene is responsible for dominant black color in dogs, while the white-phased black bear ("spirit bear") of British Columbia, Canada, is the lightest documented color variant caused by a mutation in Mc1r. Rare all-white animals have recently been discovered in a new northeastern population of the coyote in insular Newfoundland and Labrador, Canada. To investigate the causative gene and mutation of white coat in coyotes, we sequenced the three type-switching genes in white and dark-phased animals from Newfoundland. The only sequence variants unambiguously associated with white color were in Mc1r, and one of these variants causes the amino acid variant R306Ter, a premature stop codon also linked to coat color in Golden Retrievers and other dogs with yellow/red coats. The allele carrying R306Ter in coyotes matches that in the Golden Retriever at other variable amino acid sites and hence may have originated in these dogs. Coyotes experienced introgression with wolves and dogs as they colonized northeastern North America, and coyote/Golden Retriever interactions have been observed in Newfoundland. We speculate that natural selection, with or without a founder effect, may contribute to the observed frequency of white coyotes in Newfoundland, as it has contributed to the high frequency of white bears, and of a domestic dog-derived CBD allele in gray wolves.

Near neutrality, rate heterogeneity, and linkage govern mitochondrial genome evolution in Atlantic cod (Gadus morhua) and other gadine fish.

The mitochondrial DNA (mtDNA) genome figures prominently in evolutionary investigations of vertebrate animals due to a suite of characteristics that include absence of Darwinian selection, high mutation rate, and inheritance as a single linkage group. Given complete linkage and selective neutrality, mtDNA gene trees are expected to correspond to intraspecific phylogenies, and mtDNA diversity will reflect population size. The validity of these assumptions is, however, rarely tested on a genome-wide scale. Here, we analyze rates and patterns of molecular evolution among 32 whole mitochondrial genomes of Atlantic Cod (Gadus morhua) as compared with its sister taxon, the walleye pollock (Gadus [Theragra] chalcogrammus), and genomes of seven other gadine codfish. We evaluate selection within G. morhua, between sister species, and among species and intraspecific measures of linkage disequilibrium and recombination within G. morhua. Strong rate heterogeneity occurs among sites and genes at all levels of hierarchical comparison, consistent with variation in mutation rates across the genome. Neutrality indices (dN/dS) are significantly greater than unity among G. morhua genomes and between sister species, which suggests that polymorphisms within species are slightly deleterious, as expected under the nearly neutral theory of molecular evolution. Among species of gadines, dN/dS ratios are heterogeneous among genes, consistent with purifying selection and variation in functional constraint among genes rather than positive selection. The dN/dS ratio for ND4L is anomalously high across all hierarchical levels. There is no evidence for recombination within G. morhua. These patterns contrast strongly with those reported for humans: genome-wide patterns in other vertebrates should be investigated to elucidate the complex patterns of mtDNA molecular evolution.

Intraspecific phylogeographic genomics from multiple complete mtDNA genomes in Atlantic cod (Gadus morhua): origins of the "codmother," transatlantic vicariance and midglacial population expansion.

On the basis of multiple complete mitochondrial DNA genome sequences, we describe the temporal phylogeography of Atlantic cod (Gadus morhua), a lineage that has undergone a complex pattern of vicariant evolution, postglacial demographic shifts, and historic sharp population declines due to fishing and/or environmental shifts. Each of 32 fish from four spawning aggregations from the northwest Atlantic and Norway has a unique mtDNA sequence, which differs by 6-60 substitutions. Phylogenetic analysis identifies six major haplogroups that range in age from 37 to 75 KYA. The widespread haplotype identified by previous single-locus analyses at the center of a "star phylogeny" is shown to be a paraphyletic assemblage of genome lineages. The coalescent that includes all cod occurs 162 KYA. The most basal clade comprises two fish from the western Atlantic. The most recent superclade that includes all fish examined from Norway, and which includes 84% of all fish examined, dates to 128 KYA at the Sangamon/Würm interglacial, when ocean depths on continental shelves would have favored transcontinental movement. The pairwise mismatch distribution dates population expansion of this superclade to the middle of the Wisconsinan/Weichsel glaciation 59 KYA, rather than to a postglacial emergence from a marine refugium 12 KYA, or to more recent historic events. We discuss alternative scenarios for the expansion and distribution of the descendants of the "codmother" in the North Atlantic. Mitochondrial phylogenomic analyses generate highly resolved trees that enable fine-scale tests of temporal hypotheses with an accuracy not possible with single-locus methods.

Phylogeographic analysis of complete mtDNA genomes from Walleye Pollock (Gadus chalcogrammus Pallas, 1811) shows an ancient origin of genetic biodiversity.

Ursvik et al. compared the complete mitochondrial DNA (mtDNA) genome sequences of Walleye Pollock (Gadus ( = Theragra) chalcogrammus) from the Pacific Ocean with a pair of fish from an isolated population of Norwegian pollock in the Barents Sea. They concluded that the Norwegian population was recently introduced from the Pacific. We test this hypothesis within a temporal framework provided by a phylogeographic analysis of complete genomes from the pollocks' sister species, Atlantic Cod (Gadus morhua), and their divergence 3.5 mya. Pollock have a coalescent ancestor 189 +/- 25 kya. The two Norwegian fish have a common ancestor 87 +/- 7 kya, which suggests an ancient origin rather than a recent human-mediated introduction. Mitochondrial genomic biodiversity in pollock antedates the most recent glacial cycle. The clade structure of the whole-genome tree indicates that previously described single-locus mtDNA haplotypes and haplogroups are typically paraphyletic.

Phylogeographic genomics of mitochondrial DNA: Highly-resolved patterns of intraspecific evolution and a multi-species, microarray-based DNA sequencing strategy for biodiversity studies.

Phylogeographic genomics, based on multiple complete mtDNA genome sequences from within individual vertebrate species, provides highly-resolved intraspecific trees for the detailed study of evolutionary biology. We describe new biogeographic and historical insights from our studies of the genomes of codfish, wolffish, and harp seal populations in the Northwest Atlantic, and from the descendants of the founding human population of Newfoundland. Population genomics by conventional sequencing methods remains laborious. A new biotechnology, iterative DNA "re-sequencing", uses a DNA microarray to recover 30-300 kb of contiguous DNA sequence in a single experiment. Experiments with a single-species mtDNA microarray show that the method is accurate and efficient, and sufficiently species-specific to discriminate mtDNA genomes of moderately-divergent taxa. Experiments with a multi-species DNA microarray (the "ArkChip") show that simultaneous sequencing of species in different orders and classes detects SNPs within each taxon with equal accuracy as single-species-specific experiments. Iterative DNA sequencing offers a practical method for high-throughput biodiversity genomics that will enable standardized, coordinated investigation of multiple species of interest to Species at Risk and conservation biologists.

Mitochondrial genomics of gadine fishes: implications for taxonomy and biogeographic origins from whole-genome data sets.

Phylogenetic analysis of 13 substantially complete mitochondrial DNA genome sequences (14,036 bp) from 10 taxa of gadine codfishes and pollock provides highly corroborated resolution of outstanding questions on their biogeographic evolution. Of 6 resolvable nodes among species, 4 were supported by >95% of bootstrap replications in parsimony, distance, likelihood, and similarly high posterior probabilities in bayesian analyses, one by 85%-95% according to the method of analysis, and one by 99% by one method and a majority of the other two. The endemic Pacific species, walleye pollock (Theragra chalcogramma), is more closely related to the endemic Atlantic species, Atlantic cod (Gadus macrocephalus), than either is to a second Pacific endemic, Pacific cod (Gadus macrocephalus). The walleye pollock should thus be referred to the genus Gadus as originally described (Gadus chalcogrammus Pallas 1811). Arcto-Atlantic Greenland cod, previously regarded as a distinct species (G. ogac), are a genomically distinguishable subspecies within pan-Pacific G. macrocephalus. Of the 2 endemic Arctic Ocean genera, Polar cod (Boreogadus) as the outgroup to Arctic cod (Arctogadus) and Gadus sensu lato is more strongly supported than a pairing of Boreogadus and Arctogadus as sister taxa. Taking into consideration historical patterns of hydrogeography, we outline a hypothesis of the origin of the 2 endemic Pacific species as independent but simultaneous invasions through the Bering Strait from an Arcto-Atlantic ancestral lineage. In contrast to the genome data, the complete proteome sequence (3830 amino acids) resolved only 3 nodes with >95% confidence, and placed Alaska pollock outside the Gadus clade owing to reversal mutations in the ND5 locus that restore ancestral, non-Gadus, amino acid residues in that species.