Complete mitochondrial genome of the Galápagos sea lion, Zalophuswollebaeki (Carnivora, Otariidae): paratype specimen confirms separate species status.

The endangered Galápagos sea lion (Zalophuswollebaeki) inhabits the Galápagos Islands off the coast of Ecuador. We present a complete mitochondrial genome (16 465 bp) of a female paratype from the collections of the Natural History Museum Oslo, Norway, assembled from next-generation sequencing reads. It contains all canonical protein-coding, rRNA, tRNA genes, and the D-loop region. Sequence similarity is 99.93% to a previously published conspecific mitogenome sequence and 99.37% to the mitogenome sequence of the sister species Z.californianus. Sequence similarity of the D-loop region of the Z.wollebaeki paratype mitogenome is >99%, while the sequence difference to the Z.californianus sequences exceeds 2.5%. The paratype mitogenome sequence supports the taxonomic status of Z.wollebaeki as a separate species.

Insights into bear evolution from a Pleistocene polar bear genome.

The polar bear (Ursus maritimus) has become a symbol of the threat to biodiversity from climate change. Understanding polar bear evolutionary history may provide insights into apex carnivore responses and prospects during periods of extreme environmental perturbations. In recent years, genomic studies have examined bear speciation and population history, including evidence for ancient admixture between polar bears and brown bears (Ursus arctos). Here, we extend our earlier studies of a 130,000- to 115,000-y-old polar bear from the Svalbard Archipelago using a 10× coverage genome sequence and 10 new genomes of polar and brown bears from contemporary zones of overlap in northern Alaska. We demonstrate a dramatic decline in effective population size for this ancient polar bear's lineage, followed by a modest increase just before its demise. A slightly higher genetic diversity in the ancient polar bear suggests a severe genetic erosion over a prolonged bottleneck in modern polar bears. Statistical fitting of data to alternative admixture graph scenarios favors at least one ancient introgression event from brown bears into the ancestor of polar bears, possibly dating back over 150,000 y. Gene flow was likely bidirectional, but allelic transfer from brown into polar bear is the strongest detected signal, which contrasts with other published work. These findings may have implications for our understanding of climate change impacts: Polar bears, a specialist Arctic lineage, may not only have undergone severe genetic bottlenecks but also been the recipient of generalist, boreal genetic variants from brown bears during critical phases of Northern Hemisphere glacial oscillations.

High genomic diversity in the endangered East Greenland Svalbard Barents Sea stock of bowhead whales (Balaena mysticetus).

The East Greenland-Svalbard-Barents Sea (EGSB) bowhead whale stock (Balaena mysticetus) was hunted to near extinction and remains Endangered on the International Union of Conservation of Nature Red List. The intense, temporally extensive hunting pressure may have left the population vulnerable to other perturbations, such as environmental change. However, the lack of genomic baseline data renders it difficult to evaluate the impacts of various potential stressors on this stock. Twelve EGSB bowhead whales sampled in 2017/2018 were re-sequenced and mapped to a previously published draft genome. All individuals were unrelated and void of significant signs of inbreeding, with similar observed and expected homo- and heterozygosity levels. Despite the small population size, mean autosome-wide heterozygosity was 0.00102, which is higher than that of most mammals for which comparable estimates are calculated using the same parameters, and three times higher than a conspecific individual from the Eastern-Canada-West-Greenland bowhead whale stock. Demographic history analyses indicated a continual decrease of Ne from ca. 1.5 million to ca. 250,000 years ago, followed by a slight increase until ca. 100,000 years ago, followed by a rapid decrease in Ne between 50,000 and 10,000 years ago. These estimates are lower than previously suggested based on mitochondrial DNA, but suggested demographic patterns over time are similar.

Strong and lasting impacts of past global warming on baleen whales and their prey.

Global warming is affecting the population dynamics and trophic interactions across a wide range of ecosystems and habitats. Translating these real-time effects into their long-term consequences remains a challenge. The rapid and extreme warming period that occurred after the Last Glacial Maximum (LGM) during the Pleistocene-Holocene transition (7-12 thousand years ago) provides an opportunity to gain insights into the long-term responses of natural populations to periods with global warming. The effects of this post-LGM warming period have been assessed in many terrestrial taxa, whereas insights into the impacts of rapid global warming on marine taxa remain limited, especially for megafauna. In order to understand how large-scale climate fluctuations during the post-LGM affected baleen whales and their prey, we conducted an extensive, large-scale analysis of the long-term effects of the post-LGM warming on abundance and inter-ocean connectivity in eight baleen whale and seven prey (fish and invertebrates) species across the Southern and the North Atlantic Ocean; two ocean basins that differ in key oceanographic features. The analysis was based upon 7032 mitochondrial DNA sequences as well as genome-wide DNA sequence variation in 100 individuals. The estimated temporal changes in genetic diversity during the last 30,000 years indicated that most baleen whale populations underwent post-LGM expansions in both ocean basins. The increase in baleen whale abundance during the Holocene was associated with simultaneous changes in their prey and climate. Highly correlated, synchronized and exponential increases in abundance in both baleen whales and their prey in the Southern Ocean were indicative of a dramatic increase in ocean productivity. In contrast, the demographic fluctuations observed in baleen whales and their prey in the North Atlantic Ocean were subtle, varying across taxa and time. Perhaps most important was the observation that the ocean-wide expansions and decreases in abundance that were initiated by the post-LGM global warming, continued for millennia after global temperatures stabilized, reflecting persistent, long-lasting impacts of global warming on marine fauna.

Sea ice reduction drives genetic differentiation among Barents Sea polar bears.

Loss of Arctic sea ice owing to climate change is predicted to reduce both genetic diversity and gene flow in ice-dependent species, with potentially negative consequences for their long-term viability. Here, we tested for the population-genetic impacts of reduced sea ice cover on the polar bear (Ursus maritimus) sampled across two decades (1995-2016) from the Svalbard Archipelago, Norway, an area that is affected by rapid sea ice loss in the Arctic Barents Sea. We analysed genetic variation at 22 microsatellite loci for 626 polar bears from four sampling areas within the archipelago. Our results revealed a 3-10% loss of genetic diversity across the study period, accompanied by a near 200% increase in genetic differentiation across regions. These effects may best be explained by a decrease in gene flow caused by habitat fragmentation owing to the loss of sea ice coverage, resulting in increased inbreeding of local polar bears within the focal sampling areas in the Svalbard Archipelago. This study illustrates the importance of genetic monitoring for developing adaptive management strategies for polar bears and other ice-dependent species.

Genomes of Pleistocene Siberian Wolves Uncover Multiple Extinct Wolf Lineages.

Extant Canis lupus genetic diversity can be grouped into three phylogenetically distinct clades: Eurasian and American wolves and domestic dogs.1 Genetic studies have suggested these groups trace their origins to a wolf population that expanded during the last glacial maximum (LGM)1-3 and replaced local wolf populations.4 Moreover, ancient genomes from the Yana basin and the Taimyr peninsula provided evidence of at least one extinct wolf lineage that dwelled in Siberia during the Pleistocene.35 Previous studies have suggested that Pleistocene Siberian canids can be classified into two groups based on cranial morphology. Wolves in the first group are most similar to present-day populations, although those in the second group possess intermediate features between dogs and wolves.67 However, whether this morphological classification represents distinct genetic groups remains unknown. To investigate this question and the relationships between Pleistocene canids, present-day wolves, and dogs, we resequenced the genomes of four Pleistocene canids from Northeast Siberia dated between >50 and 14 ka old, including samples from the two morphological categories. We found these specimens cluster with the two previously sequenced Pleistocene wolves, which are genetically more similar to Eurasian wolves. Our results show that, though the four specimens represent extinct wolf lineages, they do not form a monophyletic group. Instead, each Pleistocene Siberian canid branched off the lineage that gave rise to present-day wolves and dogs. Finally, our results suggest the two previously described morphological groups could represent independent lineages similarly related to present-day wolves and dogs.

Arctic-adapted dogs emerged at the Pleistocene-Holocene transition.

Although sled dogs are one of the most specialized groups of dogs, their origin and evolution has received much less attention than many other dog groups. We applied a genomic approach to investigate their spatiotemporal emergence by sequencing the genomes of 10 modern Greenland sled dogs, an ~9500-year-old Siberian dog associated with archaeological evidence for sled technology, and an ~33,000-year-old Siberian wolf. We found noteworthy genetic similarity between the ancient dog and modern sled dogs. We detected gene flow from Pleistocene Siberian wolves, but not modern American wolves, to present-day sled dogs. The results indicate that the major ancestry of modern sled dogs traces back to Siberia, where sled dog-specific haplotypes of genes that potentially relate to Arctic adaptation were established by 9500 years ago.

Specialized sledge dogs accompanied Inuit dispersal across the North American Arctic.

Domestic dogs have been central to life in the North American Arctic for millennia. The ancestors of the Inuit were the first to introduce the widespread usage of dog sledge transportation technology to the Americas, but whether the Inuit adopted local Palaeo-Inuit dogs or introduced a new dog population to the region remains unknown. To test these hypotheses, we generated mitochondrial DNA and geometric morphometric data of skull and dental elements from a total of 922 North American Arctic dogs and wolves spanning over 4500 years. Our analyses revealed that dogs from Inuit sites dating from 2000 BP possess morphological and genetic signatures that distinguish them from earlier Palaeo-Inuit dogs, and identified a novel mitochondrial clade in eastern Siberia and Alaska. The genetic legacy of these Inuit dogs survives today in modern Arctic sledge dogs despite phenotypic differences between archaeological and modern Arctic dogs. Together, our data reveal that Inuit dogs derive from a secondary pre-contact migration of dogs distinct from Palaeo-Inuit dogs, and probably aided the Inuit expansion across the North American Arctic beginning around 1000 BP.

Evolutionary Implications of the microRNA- and piRNA Complement of Lepidodermella squamata (Gastrotricha).

Gastrotrichs-'hairy bellies'-are microscopic free-living animals inhabiting marine and freshwater habitats. Based on morphological and early molecular analyses, gastrotrichs were placed close to nematodes, but recent phylogenomic analyses have suggested their close relationship to flatworms (Platyhelminthes) within Spiralia. Small non-coding RNA data on e.g., microRNAs (miRNAs) and PIWI-interacting RNAs (piRNA) may help to resolve this long-standing question. MiRNAs are short post-transcriptional gene regulators that together with piRNAs play key roles in development. In a 'multi-omics' approach we here used small-RNA sequencing, available transcriptome and genomic data to unravel the miRNA- and piRNA complements along with the RNAi (RNA interference) protein machinery of Lepidodermella squamata (Gastrotricha, Chaetonotida). We identified 52 miRNA genes representing 35 highly conserved miRNA families specific to Eumetazoa, Bilateria, Protostomia, and Spiralia, respectively, with overall high similarities to platyhelminth miRNA complements. In addition, we found four large piRNA clusters that also resemble flatworm piRNAs but not those earlier described for nematodes. Congruently, transcriptomic annotation revealed that the Lepidodermella protein machinery is highly similar to flatworms, too. Taken together, miRNA, piRNA, and protein data support a close relationship of gastrotrichs and flatworms.

Mitochondrial genomes of ancient bowhead whales (Balaena mysticetus) from Svalbard.

The endangered Spitsbergen stock of bowhead whales (Balaena mysticetus) has once been large with up to estimated 100,000 individuals. Genetic diversity of the extant Spitsbergen stock is unknown. We present 10 complete mitochondrial genomes of heterochronous ancient bowhead whale samples from Svalbard (14C age estimate range: 215-8885 years) obtained via NGS of total genomic DNA extracts. The ten mitogenomes differed by nucleotide substitutions and/or indels, and there was a total of 160 variable positions. The average nucleotide diversity was π = 0.0029. There was no statistically significant correlation between genetic divergence and time.

Population genomics of grey wolves and wolf-like canids in North America.

North America is currently home to a number of grey wolf (Canis lupus) and wolf-like canid populations, including the coyote (Canis latrans) and the taxonomically controversial red, Eastern timber and Great Lakes wolves. We explored their population structure and regional gene flow using a dataset of 40 full genome sequences that represent the extant diversity of North American wolves and wolf-like canid populations. This included 15 new genomes (13 North American grey wolves, 1 red wolf and 1 Eastern timber/Great Lakes wolf), ranging from 0.4 to 15x coverage. In addition to providing full genome support for the previously proposed coyote-wolf admixture origin for the taxonomically controversial red, Eastern timber and Great Lakes wolves, the discriminatory power offered by our dataset suggests all North American grey wolves, including the Mexican form, are monophyletic, and thus share a common ancestor to the exclusion of all other wolves. Furthermore, we identify three distinct populations in the high arctic, one being a previously unidentified "Polar wolf" population endemic to Ellesmere Island and Greenland. Genetic diversity analyses reveal particularly high inbreeding and low heterozygosity in these Polar wolves, consistent with long-term isolation from the other North American wolves.

Interspecific Gene Flow Shaped the Evolution of the Genus Canis.

The evolutionary history of the wolf-like canids of the genus Canis has been heavily debated, especially regarding the number of distinct species and their relationships at the population and species level [1-6]. We assembled a dataset of 48 resequenced genomes spanning all members of the genus Canis except the black-backed and side-striped jackals, encompassing the global diversity of seven extant canid lineages. This includes eight new genomes, including the first resequenced Ethiopian wolf (Canis simensis), one dhole (Cuon alpinus), two East African hunting dogs (Lycaon pictus), two Eurasian golden jackals (Canis aureus), and two Middle Eastern gray wolves (Canis lupus). The relationships between the Ethiopian wolf, African golden wolf, and golden jackal were resolved. We highlight the role of interspecific hybridization in the evolution of this charismatic group. Specifically, we find gene flow between the ancestors of the dhole and African hunting dog and admixture between the gray wolf, coyote (Canis latrans), golden jackal, and African golden wolf. Additionally, we report gene flow from gray and Ethiopian wolves to the African golden wolf, suggesting that the African golden wolf originated through hybridization between these species. Finally, we hypothesize that coyotes and gray wolves carry genetic material derived from a "ghost" basal canid lineage.

Finding Evolutionary Processes Hidden in Cryptic Species.

Cryptic species could represent a substantial fraction of biodiversity. However, inconsistent definitions and taxonomic treatment of cryptic species prevent informed estimates of their contribution to biodiversity and impede our understanding of their evolutionary and ecological significance. We propose a conceptual framework that recognizes cryptic species based on their low levels of phenotypic (morphological) disparity relative to their degree of genetic differentiation and divergence times as compared with non-cryptic species. We discuss how application of a more rigorous definition of cryptic species in taxonomic practice will lead to more accurate estimates of their prevalence in nature, better understanding of their distribution patterns on the tree of life, and increased abilities to resolve the processes underlying their evolution.

Atp8 is in the ground pattern of flatworm mitochondrial genomes.

BACKGROUND: To date, mitochondrial genomes of more than one hundred flatworms (Platyhelminthes) have been sequenced. They show a high degree of similarity and a strong taxonomic bias towards parasitic lineages. The mitochondrial gene atp8 has not been confidently annotated in any flatworm sequenced to date. However, sampling of free-living flatworm lineages is incomplete. We addressed this by sequencing the mitochondrial genomes of the two small-bodied (about 1 mm in length) free-living flatworms Stenostomum sthenum and Macrostomum lignano as the first representatives of the earliest branching flatworm taxa Catenulida and Macrostomorpha respectively. RESULTS: We have used high-throughput DNA and RNA sequence data and PCR to establish the mitochondrial genome sequences and gene orders of S. sthenum and M. lignano. The mitochondrial genome of S. sthenum is 16,944 bp long and includes a 1,884 bp long inverted repeat region containing the complete sequences of nad3, rrnS, and nine tRNA genes. The model flatworm M. lignano has the smallest known mitochondrial genome among free-living flatworms, with a length of 14,193 bp. The mitochondrial genome of M. lignano lacks duplicated genes, however, tandem repeats were detected in a non-coding region. Mitochondrial gene order is poorly conserved in flatworms, only a single pair of adjacent ribosomal or protein-coding genes - nad4l-nad4 - was found in S. sthenum and M. lignano that also occurs in other published flatworm mitochondrial genomes. Unexpectedly, we unambiguously identified the full metazoan mitochondrial protein-coding gene complement including atp8 in S. sthenum and M. lignano. A subsequent search detected atp8 in all mitochondrial genomes of polyclad flatworms published to date, although the gene wasn't previously annotated in these species. CONCLUSIONS: Manual, but not automated genome annotation revealed the presence of atp8 in basally branching free-living flatworms, signifying both the importance of manual data curation and of diverse taxon sampling. We conclude that the loss of atp8 within flatworms is restricted to the parasitic taxon Neodermata.

The mitochondrial genome of the egg-laying flatworm Aglaiogyrodactylus forficulatus (Platyhelminthes: Monogenoidea).

BACKGROUND: The rather species-poor oviparous gyrodactylids are restricted to South America. It was suggested that they have a basal position within the otherwise viviparous Gyrodactylidae. Accordingly, it was proposed that the species-rich viviparous gyrodactylids diversified and dispersed from there. METHODS: The mitochondrial genome of Aglaiogyrodactylus forficulatus was bioinformatically assembled from next-generation illumina MiSeq sequencing reads, annotated, and compared to previously published mitochondrial genomes of other monogenoidean flatworm species. RESULTS: The mitochondrial genome of A. forficulatus consists of 14,371 bp with an average A + T content of 75.12 %. All expected 12 protein coding, 22 tRNA, and 2 rRNA genes were identified. Furthermore, there were two repetitive non-coding regions essentially consisting of 88 bp and 233 bp repeats, respectively. Maximum Likelihood analyses placed the mitochondrial genome of A. forficulatus in a well-supported clade together with the viviparous Gyrodactylidae species. The gene order differs in comparison to that of other monogenoidean species, with rearrangements mainly affecting tRNA genes. In comparison to Paragyrodactylus variegatus, four gene order rearrangements, i.e. three transpositions and one complex tandem-duplication-random-loss event, were detected. CONCLUSION: Mitochondrial genome sequence analyses support a basal position of the oviparous A. forficulatus within Gyrodactylidae, and a sister group relationship of the oviparous and viviparous forms.

Genetic diversity of historical Atlantic walruses (Odobenus rosmarus rosmarus) from Bjørnøya and Håøya (Tusenøyane), Svalbard, Norway.

BACKGROUND: The population size of Atlantic walruses (Odobenus rosmarus rosmarus) is depleted relative to historical abundance levels. In Svalbard, centuries of over-exploitation brought the walrus herds to the verge of extinction, and such bottlenecks may have caused loss of genetic variation. To address this for Svalbard walruses, mitochondrial haplotypes of historical walruses from two major haul-out sites, Bjørnøya and Håøya, within the Archipelago were explored using bone samples from animals killed during the peak period of harvesting. RESULTS: Using ancient DNA methodologies, the mitochondrial NADH dehydrogenase 1 (ND1) gene, the cytochrome c oxidase 1 (COI) gene, and the control region (CR) were targeted for 15 specimens from Bjørnøya (of which five were entirely negative) and 9 specimens from Håøya (of which one was entirely negative). While ND1 and COI sequences were obtained for only a few samples, the CR delivered the most comprehensive data set, and the average genetic distance among historic Svalbard samples was 0.0028 (SD = 0.0023). CONCLUSIONS: The CR sequences from the historical samples appear to be nested among contemporary Atlantic walruses, and no distinct mitochondrial haplogroups were identified in the historical samples that may have been lost during the periods of extensive hunting. However, given the low sample size and poor phylogenetic resolution it cannot be excluded that such haplogroups existed.

Mitogenomes of contemporary Spitsbergen stock bowhead whales (Balaena mysticetus).

The Spitsbergen stock of bowhead whales (Balaena mysticetus) is considered Critically Endangered by IUCN. Over recent decades, there have been only a few sightings, and very few biological samples are available for molecular analyses. Hence, genetic diversity of the extant Spitsbergen stock is unknown. Here, we present mitochondrial genomes from eight skin biopsy samples that were collected in 2006 and 2010. There were three different haplotypes, two of which have not previously been detected. Six samples shared the same haplotype, indicating that they were obtained from closely related whales, or possibly the same individual. Average nucleotide diversity was π = 0.0037, with a total of 93 variable positions among the haplotypes.