Fine-scale matrilineal population structure in the Galapagos fur seal and its implications for conservation management

Genomic evidence for homoploid hybrid speciation in a marine mammal apex predator

Hybridization is widespread and constitutes an important source of genetic variability and evolution. In animals, its role in generating novel and independent lineages (hybrid speciation) has been strongly debated, with only a few cases supported by genomic data. The South American fur seal (SAfs) Arctocephalus australis is a marine apex predator of Pacific and Atlantic waters, with a disjunct set of populations in Peru and Northern Chile [Peruvian fur seal (Pfs)] with controversial taxonomic status. We demonstrate, using complete genome and reduced representation sequencing, that the Pfs is a genetically distinct species with an admixed genome that originated from hybridization between the SAfs and the Galapagos fur seal (Arctocephalus galapagoensis) ~400,000 years ago. Our results strongly support the origin of Pfs by homoploid hybrid speciation over alternative introgression scenarios. This study highlights the role of hybridization in promoting species-level biodiversity in large vertebrates.

Population Genetics and Phylogeography of Galapagos Fur Seals

Pinnipeds found across islands provide an ideal opportunity to examine the evolutionary process of population subdivision affected by several mechanisms. Here, we report the genetic consequences of the geographic distribution of rookeries in Galapagos fur seals (GFS: Arctocephalus galapagoensis) in creating population structure. We show that rookeries across four islands (nine rookeries) are genetically structured into the following major groups: 1) a western cluster of individuals from Fernandina; 2) a central group from north and east Isabela, Santiago, and Pinta; and possibly, 3) a third cluster in the northeast from Pinta. Furthermore, asymmetric levels of gene flow obtained from eight microsatellites found migration from west Isabela to Fernandina islands (number of migrants Nm = 1), with imperceptible Nm in any other direction. Our findings suggest that the marked structuring of populations recovered in GFS is likely related to an interplay between long-term site fidelity and long-distance migration in both male and female individuals, probably influenced by varying degrees of marine productivity.

Declines and recovery in endangered Galapagos pinnipeds during the El Niño event

Currently, the Galapagos sea lion (GSL, Zalophus wollebaeki) and Galapagos fur seal (GFS, Arctocephalus galapagoensis) are among the most important endemic species for conservation in the Galapagos Archipelago. Both are classified as "Endangered" since their populations have undergone drastic declines over the last several decades. In this study we estimated the abundance of both otariids, and their population trends based using counts conducted between 2014 and 2018 in all their rookeries, and we analyzed the influence of environmental variability on pup production. The GSL population size in 2018 in the archipelago was estimated to be between 17,000 to 24,000 individuals and has increased at an average annual rate of 1% over the last five years after applying correction factors. The highest number of GSL counted in the archipelago was in 2014 followed by a population decline of 23.8% in 2015 that was associated with the El Niño event that occurred during that year. Following this event, the population increased mainly in the northern, central and southeastern rookeries. The GSL pup abundance showed a decreasing trend with the increase in intensity of the El Niño. The GFS population in 2018 was counted in 3,093 individuals and has increased at an annual rate of 3% from 2014 to 2018. A high number of GFS counted in 2014 was followed by a population decrease of 38% in 2015, mainly in the western rookeries. There was interannual population fluctuations and different growth trends among regions of the archipelago. GSL and GFS pup abundance has a strong decreasing tendency with the increase in the subthermocline temperature (ST) and the El Niño 1 + 2 index. Our results provide evidence that both species are highly vulnerable to periodic oceanographic-atmospheric events in the Galapagos Archipelago which impact prey abundance and the flow of energy in the unique Galapagos ecosystem.

High-quality whole-genome sequence of an abundant Holarctic odontocete, the harbour porpoise (Phocoena phocoena)

The harbour porpoise (Phocoena phocoena) is a highly mobile cetacean found across the Northern hemisphere. It occurs in coastal waters and inhabits basins that vary broadly in salinity, temperature and food availability. These diverse habitats could drive subtle differentiation among populations, but examination of this would be best conducted with a robust reference genome. Here, we report the first harbour porpoise genome, assembled de novo from an individual originating in the Kattegat Sea (Sweden). The genome is one of the most complete cetacean genomes currently available, with a total size of 2.39 Gb and 50% of the total length found in just 34 scaffolds. Using 122 of the longest scaffolds, we were able to show high levels of synteny with the genome of the domestic cattle (Bos taurus). Our draft annotation comprises 22,154 predicted genes, which we further annotated through matches to the NCBI nucleotide database, GO categorization and motif prediction. Within the predicted genes, we have confirmed the presence of >20 genes or gene families that have been associated with adaptive evolution in other cetaceans. Overall, this genome assembly and draft annotation represent a crucial addition to the genomic resources currently available for the study of porpoises and Phocoenidae evolution, phylogeny and conservation.

Population structuring of the ubiquitous stingless bee Tetragonisca angustula in southern Brazil as revealed by microsatellite and mitochondrial markers

Tetragonisca angustula is one of the most widespread stingless bees in the Neotropics. This species swarms frequently and is extremely successful in urban environments. In addition, it is one of the most popular stingless bee species for beekeeping in Latin America, so nest transportation and trading is common. Nest transportation can change the genetic structure of the host population, reducing inbreeding and increasing homogenization. Here, we evaluate the genetic structure of 17 geographic populations of T. angustula in southern Brazil to quantify the level of genetic differentiation between populations. Analyses were conducted on partially sequenced mitochondrial genes and 11 microsatellite loci of 1002 workers from 457 sites distributed on the mainland and on 3 islands. Our results show that T. angustula populations are highly differentiated as demonstrated by mitochondrial DNA (mtDNA) and microsatellite markers. Of 73 haplotypes, 67 were population-specific. MtDNA diversity was low in 9 populations but microsatellite diversity was moderate to high in all populations. Microsatellite data suggest 10 genetic clusters and low level of gene flow throughout the studied area. However, physical barriers, such as rivers and mountain ranges, or the presence or absence of forest appear to be unrelated to population clusters. Factors such as low dispersal, different ecological conditions, and isolation by distance are most likely shaping the population structure of this species. Thus far, nest transportation has not influenced the general population structure in the studied area. However, due to the genetic structure we found, we recommend that nest transportation should only occur within and between populations that are genetically similar.

Ancient female philopatry, asymmetric male gene flow, and synchronous population expansion support the influence of climatic oscillations on the evolution of South American sea lion (Otaria flavescens)

The South American sea lion (Otaria flavescens) is widely distributed along the southern Atlantic and Pacific coasts of South America with a history of significant commercial exploitation. We aimed to evaluate the population genetic structure and the evolutionary history of South American sea lion along its distribution by analyses of mitochondrial DNA (mtDNA) and 10 nuclear microsatellites loci. We analyzed 147 sequences of mtDNA control region and genotyped 111 individuals of South American sea lion for 10 microsatellite loci, representing six populations (Peru, Northern Chile, Southern Chile, Uruguay (Brazil), Argentina and Falkland (Malvinas) Islands) and covering the entire distribution of the species. The mtDNA phylogeny shows that haplotypes from the two oceans comprise two very divergent clades as observed in previous studies, suggesting a long period (>1 million years) of low inter-oceanic female gene flow. Bayesian analysis of bi-parental genetic diversity supports significant (but less pronounced than mitochondrial) genetic structure between Pacific and Atlantic populations, although also suggested some inter-oceanic gene flow mediated by males. Higher male migration rates were found in the intra-oceanic population comparisons, supporting very high female philopatry in the species. Demographic analyses showed that populations from both oceans went through a large population expansion ~10,000 years ago, suggesting a very similar influence of historical environmental factors, such as the last glacial cycle, on both regions. Our results support the proposition that the Pacific and Atlantic populations of the South American sea lion should be considered distinct evolutionarily significant units, with at least two managements units in each ocean.

Natal and breeding philopatry of female Steller sea lions in southeastern Alaska

Information on drivers of dispersal is critical for wildlife conservation but is rare for long-lived marine mammal species with large geographic ranges. We fit multi-state mark-recapture models to resighting data of 369 known-aged Steller sea lion (Eumetopias jubatus) females marked as pups on their natal rookeries in southeastern Alaska from 1994-2005 and monitored from 2001-15. We estimated probabilities of females being first observed parous at their natal site (natal philopatry), and of not moving breeding sites among years (breeding philopatry) at large (> 400 km, all five rookeries in southeastern Alaska) and small (< 4 km, all islands within the largest rookery, Forrester Island Complex, F) spatial scales. At the rookery scale, natal philopatry was moderately high (0.776-0.859) for most rookeries and breeding philopatry was nearly 1, with < 3% of females switching breeding rookeries between years. At more populous islands at F, natal philopatry was 0.500-0.684 versus 0.295-0.437 at less populous islands, and breeding philopatry was 0.919-0.926 versus 0.604-0.858. At both spatial scales, the probability of pupping at a non-natal site increased with population size of, and declined with distance from, the destination site. Natal philopatry of < 1 would increase gene flow, improve population resilience, and promote population recovery after decline in a heterogeneous environment. Very high breeding philopatry suggests that familiarity with neighboring females and knowledge of the breeding site (the topography of pupping sites and nearby foraging locations) may be a critical component to reproductive strategies of sea lions.

Population structure and historical demography of South American sea lions provide insights into the catastrophic decline of a marine mammal population

Understanding the causes of population decline is crucial for conservation management. We therefore used genetic analysis both to provide baseline data on population structure and to evaluate hypotheses for the catastrophic decline of the South American sea lion (Otaria flavescens) at the Falkland Islands (Malvinas) in the South Atlantic. We genotyped 259 animals from 23 colonies across the Falklands at 281 bp of the mitochondrial hypervariable region and 22 microsatellites. A weak signature of population structure was detected, genetic diversity was moderately high in comparison with other pinniped species, and no evidence was found for the decline being associated with a strong demographic bottleneck. By combining our mitochondrial data with published sequences from Argentina, Brazil, Chile and Peru, we also uncovered strong maternally directed population structure across the geographical range of the species. In particular, very few shared haplotypes were found between the Falklands and South America, and this was reflected in correspondingly low migration rate estimates. These findings do not support the prominent hypothesis that the decline was caused by migration to Argentina, where large-scale commercial harvesting operations claimed over half a million animals. Thus, our study not only provides baseline data for conservation management but also reveals the potential for genetic studies to shed light upon long-standing questions pertaining to the history and fate of natural populations.