Population structure and connectivity among coastal and freshwater Kelp Gull (Larus dominicanus) populations from Patagonia

The genetic identification of evolutionary significant units and information on their connectivity can be used to design effective management and conservation plans for species of concern. Despite having high dispersal capacity, several seabird species show population structure due to both abiotic and biotic barriers to gene flow. The Kelp Gull is the most abundant species of gull in the southern hemisphere. In Argentina it reproduces in both marine and freshwater environments, with more than 100,000 breeding pairs following a metapopulation dynamic across 140 colonies in the Atlantic coast of Patagonia. However, little is known about the demography and connectivity of inland populations. We aim to provide information on the connectivity of the largest freshwater colonies (those from Nahuel Huapi Lake) with the closest Pacific and Atlantic populations to evaluate if these freshwater colonies are receiving immigrants from the larger coastal populations. We sampled three geographic regions (Nahuel Huapi Lake and the Atlantic and Pacific coasts) and employed a reduced-representation genomic approach to genotype individuals for single-nucleotide polymorphisms (SNPs). Using clustering and phylogenetic analyses we found three genetic groups, each corresponding to one of our sampled regions. Individuals from marine environments are more closely related to each other than to those from Nahuel Huapi Lake, indicating that the latter population constitutes the first freshwater Kelp Gull colony to be identified as an evolutionary significant unit in Patagonia.

Phylogenomics of the world's otters

Comparative whole-genome analyses hold great power to illuminate commonalities and differences in the evolution of related species that share similar ecologies. The mustelid subfamily Lutrinae includes 13 currently recognized extant species of otters,^1-5 a semiaquatic group whose evolutionary history is incompletely understood. We assembled a dataset comprising 24 genomes from all living otter species, 14 of which were newly sequenced. We used this dataset to infer phylogenetic relationships and divergence times, to characterize patterns of genome-wide genealogical discordance, and to investigate demographic history and current genomic diversity. We found that genera Lutra, Aonyx, Amblonyx, and Lutrogale form a coherent clade that should be synonymized under Lutra, simplifying the taxonomic structure of the subfamily. The poorly known tropical African Aonyx congicus and the more widespread Aonyx capensis were found to be reciprocally monophyletic (having diverged 440,000 years ago), supporting the validity of the former as a distinct species. We observed variable changes in effective population sizes over time among otters within and among continents, although several species showed similar trends of expansions and declines during the last 100,000 years. This has led to different levels of genomic diversity assessed by overall heterozygosity, genome-wide SNV density, and run of homozygosity burden. Interestingly, there were cases in which diversity metrics were consistent with the current threat status (mostly based on census size), highlighting the potential of genomic data for conservation assessment. Overall, our results shed light on otter evolutionary history and provide a framework for further in-depth comparative genomic studies targeting this group.

Genetics, Morphometrics and Health Characterization of Green Turtle Foraging Grounds in Mainland and Insular Chile

Simple Summary

Chilean waters constitute a foraging habitat for the endangered green turtle. Information about this species in the country has increased in recent years; nevertheless, little is known of its ecology and health status. Additionally, some populations have drastically decreased, probably due to human factors. Here, we studied the proportion of sex, age, morphological variation, genetic characteristics, origin, and health status of green turtles in mainland and insular Chile. We found that turtles from both regions are morphologically and genetically different. Individuals from the mainland territory are juveniles and probably originated from Galapagos. In contrast, the insular territory hosts juveniles and adults that probably originated from Galapagos and French Polynesia. We also found that turtles from both regions are facing numerous anthropic threats that must be controlled. We suggest the creation of protected areas for mainland foraging grounds, and strengthen the administrative plan of the insular region to ensure sea turtle population health.

Abstract

Two divergent genetic lineages have been described for the endangered green turtle in the Pacific Ocean, occurring sympatrically in some foraging grounds. Chile has seven known green turtle foraging grounds, hosting mainly juveniles of different lineages. Unfortunately, anthropic factors have led to the decline or disappearance of most foraging aggregations. We investigated age-class/sex structure, morphological variation, genetic diversity and structure, and health status of turtles from two mainland (Bahia Salado and Playa Chinchorro) and one insular (Easter Island) Chilean foraging grounds. Bahia Salado is composed of juveniles, and with Playa Chinchorro, exclusively harbors individuals of the north-central/eastern Pacific lineage, with Galapagos as the major genetic contributor. Conversely, Easter Island hosts juveniles and adults from both the eastern Pacific and French Polynesia. Morphological variation was found between lineages and foraging grounds, suggesting an underlying genetic component but also an environmental influence. Turtles from Easter Island, unlike Bahia Salado, exhibited injuries/alterations probably related to anthropic threats. Our findings point to establishing legal protection for mainland Chile’s foraging grounds, and to ensure that the administrative plan for Easter Island’s marine protected area maintains ecosystem health, turtle population viability, and related cultural and touristic activities.

Positive selection over the mitochondrial genome and its role in the diversification of gentoo penguins in response to adaptation in isolation

Although mitochondrial DNA has been widely used in phylogeography, evidence has emerged that factors such as climate, food availability, and environmental pressures that produce high levels of stress can exert a strong influence on mitochondrial genomes, to the point of promoting the persistence of certain genotypes in order to compensate for the metabolic requirements of the local environment. As recently discovered, the gentoo penguins (Pygoscelis papua) comprise four highly divergent lineages across their distribution spanning the Antarctic and sub-Antarctic regions. Gentoo penguins therefore represent a suitable animal model to study adaptive processes across divergent environments. Based on 62 mitogenomes that we obtained from nine locations spanning all four gentoo penguin lineages, we demonstrated lineage-specific nucleotide substitutions for various genes, but only lineage-specific amino acid replacements for the ND1 and ND5 protein-coding genes. Purifying selection (dN/dS < 1) is the main driving force in the protein-coding genes that shape the diversity of mitogenomes in gentoo penguins. Positive selection (dN/dS > 1) was mostly present in codons of the Complex I (NADH genes), supported by two different codon-based methods at the ND1 and ND4 in the most divergent lineages, the eastern gentoo penguin from Crozet and Marion Islands and the southern gentoo penguin from Antarctica respectively. Additionally, ND5 and ATP6 were under selection in the branches of the phylogeny involving all gentoo penguins except the eastern lineage. Our study suggests that local adaptation of gentoo penguins has emerged as a response to environmental variability promoting the fixation of mitochondrial haplotypes in a non-random manner. Mitogenome adaptation is thus likely to have been associated with gentoo penguin diversification across the Southern Ocean and to have promoted their survival in extreme environments such as Antarctica. Such selective processes on the mitochondrial genome may also be responsible for the discordance detected between nuclear- and mitochondrial-based phylogenies of gentoo penguin lineages.

Adaptation and cryptic pseudogenization in penguin Toll-like Receptors

Penguins (Sphenisciformes) are an iconic order of flightless, diving seabirds distributed across a large latitudinal range in the Southern Hemisphere. The extensive area over which penguins are endemic is likely to have fostered variation in pathogen pressure, which in turn will have imposed differential selective pressures on the penguin immune system. At the front line of pathogen detection and response, the Toll-like receptors (TLRs) provide insight into host evolution in the face of microbial challenge. TLRs respond to conserved pathogen-associated molecular patterns and are frequently found to be under positive selection, despite retaining specificity for defined agonist classes. We undertook a comparative immunogenetics analysis of TLRs for all penguin species and found evidence of adaptive evolution that was largely restricted to the cell surface-expressed TLRs, with evidence of positive selection at, or near, key agonist-binding sites in TLR1B, TLR4, and TLR5. Intriguingly, TLR15, which is activated by fungal products, appeared to have been pseudogenized multiple times in the Eudyptes spp., but a full-length form was present as a rare haplotype at the population level. However, in vitro analysis revealed that even the full-length form of Eudyptes TLR15 was nonfunctional, indicating an ancestral cryptic pseudogenization prior to its eventual disruption multiple times in the Eudyptes lineage. This unusual pseudogenization event could provide an insight into immune adaptation to fungal pathogens such as Aspergillus, which is responsible for significant mortality in wild and captive bird populations.

Origins and evolution of extreme life span in Pacific Ocean rockfishes

Pacific Ocean rockfishes (genus Sebastes) exhibit extreme variation in life span, with some species being among the most long-lived extant vertebrates. We de novo assembled the genomes of 88 rockfish species and from these identified repeated signatures of positive selection in DNA repair pathways in long-lived taxa and 137 longevity-associated genes with direct effects on life span through insulin signaling and with pleiotropic effects through size and environmental adaptations. A genome-wide screen of structural variation reveals copy number expansions in the immune modulatory butyrophilin gene family in long-lived species. The evolution of different rockfish life histories is coupled to genetic diversity and reshapes the mutational spectrum driving segregating CpG→TpG variants in long-lived species. These analyses highlight the genetic innovations that underlie life history trait adaptations and, in turn, how they shape genomic diversity.

Green, yellow or black? Genetic differentiation and adaptation signatures in a highly migratory marine turtle

Marine species may exhibit genetic structure accompanied by phenotypic differentiation related to adaptation despite their high mobility. Two shape-based morphotypes have been identified for the green turtle (Chelonia mydas) in the Pacific Ocean: the south-central/western or yellow turtle and north-central/eastern or black turtle. The genetic differentiation between these morphotypes and the adaptation of the black turtle to environmentally contrasting conditions of the eastern Pacific region has remained a mystery for decades. Here we addressed both questions using a reduced-representation genome approach (Dartseq; 9473 neutral SNPs) and identifying candidate outlier loci (67 outlier SNPs) of biological relevance between shape-based morphotypes from eight Pacific foraging grounds (n = 158). Our results support genetic divergence between morphotypes, probably arising from strong natal homing behaviour. Genes and enriched biological functions linked to thermoregulation, hypoxia, melanism, morphogenesis, osmoregulation, diet and reproduction were found to be outliers for differentiation, providing evidence for adaptation of C. mydas to the eastern Pacific region and suggesting independent evolutionary trajectories of the shape-based morphotypes. Our findings support the evolutionary distinctness of the enigmatic black turtle and contribute to the adaptive research and conservation genomics of a long-lived and highly mobile vertebrate.

Specialist versus generalist parasites: the interactions between host diversity, environment and geographic barriers in avian malaria

The specialist versus generalist strategies of hemoparasites in relation to their avian host, as well as environmental factors, can influence their prevalence, diversity and distribution. In this paper we investigated the influence of avian host species, as well as the environmental and geographical factors, on the strategies of Haemoproteus and Plasmodium hemoparasites. We determined prevalence and diversity by targeting their cytochrome b (Cytb) in a total of 2,590 passerine samples from 138 localities of Central and South America, and analysed biogeographic patterns and host-parasite relationships. We found a total prevalence of 23.2%. Haemoproteus presented a higher prevalence (15.3%) than Plasmodium (4.3%), as well as a higher diversity and host specificity. We determined that Plasmodium and Haemoproteus prevalences correlated positively with host diversity (Shannon index) and were significantly influenced by bird diversity, demonstrating a possible "amplification effect". We found an effect of locality and the avian family for prevalences of Haemoproteus and Plasmodium. These results suggest that Haemoproteus is more specialist than Plasmodium and could be mostly influenced by its avian host and the Andes Mountains.

Integrating morphological and genetic data at different spatial scales in a cosmopolitan marine turtle species: challenges for management and conservation

Patterns of genetic structure in highly mobile marine vertebrates may be accompanied by phenotypic variation. Most studies in marine turtles focused on population genetic structure have been performed at rookeries. We studied whether genetic and morphological variation of the endangered green turtle (Chelonia mydas) is consistent geographically, focusing on foraging grounds. An association between population genetic structure and body shape variation at broad (inter-lineage) and fine (foraging grounds) scales was predicted and analysed using mitochondrial DNA and geometric morphometrics. Although genetic and phenotypic differentiation patterns were congruent between lineages, no fine-scale association was found, suggesting adaptive divergence. Connectivity among Pacific foraging grounds found here suggests that temperatures of ocean surface currents may influence the genetic structure of C. mydas on a broad scale. Our results suggest that vicariance, dispersal, life-history traits and ecological conditions operating in foraging grounds have shaped the intraspecific morphology and genetic diversity of this species. Considering a range of geographic and temporal scales is useful when management strategies are required for cosmopolitan species. Integrating morphological and genetic tools at different spatial scales, conservation management is proposed based on protection of neutral and adaptive diversity. This approach opens new questions and challenges, especially regarding conservation genetics in cosmopolitan species.

Evidence of Pathogen-Induced Immunogenetic Selection across the Large Geographic Range of a Wild Seabird

Over evolutionary time, pathogen challenge shapes the immune phenotype of the host to better respond to an incipient threat. The extent and direction of this selection pressure depend on the local pathogen composition, which is in turn determined by biotic and abiotic features of the environment. However, little is known about adaptation to local pathogen threats in wild animals. The Gentoo penguin (Pygoscelis papua) is a species complex that lends itself to the study of immune adaptation because of its circumpolar distribution over a large latitudinal range, with little or no admixture between different clades. In this study, we examine the diversity in a key family of innate immune genes-the Toll-like receptors (TLRs)-across the range of the Gentoo penguin. The three TLRs that we investigated present varying levels of diversity, with TLR4 and TLR5 greatly exceeding the diversity of TLR7. We present evidence of positive selection in TLR4 and TLR5, which points to pathogen-driven adaptation to the local pathogen milieu. Finally, we demonstrate that two positively selected cosegregating sites in TLR5 are sufficient to alter the responsiveness of the receptor to its bacterial ligand, flagellin. Taken together, these results suggest that Gentoo penguins have experienced distinct pathogen-driven selection pressures in different environments, which may be important given the role of the Gentoo penguin as a sentinel species in some of the world's most rapidly changing environments.

Metapopulation dynamics and foraging plasticity in a highly vagile seabird, the southern rockhopper penguin

Population connectivity is driven by individual dispersal potential and modulated by natal philopatry. In seabirds, high vagility facilitates dispersal yet philopatry is also common, with foraging area overlap often correlated with population connectivity. We assess the interplay between these processes by studying past and current connectivity and foraging niche overlap among southern rockhopper penguin colonies of the coast of southern South America using genomic and stable isotope analyses. We found two distinct genetic clusters and detected low admixture between northern and southern colonies. Stable isotope analysis indicated niche variability between colonies, with Malvinas/Falklands colonies encompassing the species entire isotopic foraging niche, while the remaining colonies had smaller, nonoverlapping niches. A recently founded colony in continental Patagonia differed in isotopic niche width and position with Malvinas/Falklands colonies, its genetically identified founder population, suggesting the exploitation of novel foraging areas and/or prey items. Additionally, dispersing individuals found dead across the Patagonian shore in an unusual mortality event were also assigned to the northern cluster, suggesting northern individuals reach southern localities, but do not breed in these colonies. Facilitated by variability in foraging strategies, and especially during unfavorable conditions, the number of dispersing individuals may increase and enhance the probability of founding new colonies. Metapopulation demographic dynamics in seabirds should account for interannual variability in dispersal behavior and pay special attention to extreme climatic events, classically related to negative effects on population trends.

Latitudinal gradients of haemosporidian parasites: Prevalence, diversity and drivers of infection in the Thorn-tailed Rayadito (Aphrastura spinicauda)

Latitudinal gradients are well-suited systems that may be helpful explaining distribution of haemosporidian parasites and host susceptibility. We studied the prevalence, diversity and drivers of haemosporidian parasites (Leucocytozoon, Plasmodium and Haemoproteus) along a latitudinal gradient (30°–56° S), that encompass the total distribution (~3,000 km) of the Thorn-tailed Rayadito (Aphrastura spinicauda) in the South American temperate forests from Chile. We analyzed 516 individuals from 18 localities between 2010 and 2017 and observed an overall prevalence of 28.3% for haemosporidian parasites. Leucocytozoon was the most prevalent parasite (25.8%). We recorded 19 distinct lineages (13 for Leucocytozoon, five for Plasmodium, and one for Haemoproteus). Differences in haemosporidian prevalence and diversity by genus and type of habitat were observed in the latitudinal gradient. Further, we support the existence of a latitudinal associate distribution of Leucocytozoids in South America, where prevalence and diversity increase toward higher latitudes. Distribution of Leucocytozoon was associated with sub-antarctic habitat (higher latitude) and explained by cold temperature and high precipitation. On the other hand, we lacked to find a latitudinal associate pattern for Plasmodium and Haemoproteus, however low prevalence and high diversity were recorded in areas considered as a hotspot of biodiversity in Central Chile. Our findings confirmed the importance of habitat and climatic variables explaining prevalence, diversity and distribution of haemosporidian parasites in a huge latitudinal gradient, belonging the distribution of the Thorn-tailed Rayadito in the world's southernmost forests ecosystems.

•

Prevalence and diversity differed between haemosporidian genera and habitats in the latitudinal gradient.

•

Leucocytozoon was the most prevalent genus in the latitudinal distribution.

•

Leucocytozoon prevalence and diversity were associated to higher latitudes.

•

Distribution of Leucocytozoon was explained by cold temperature and high precipitation.

More than the eye can see: Genomic insights into the drivers of genetic differentiation in Royal/Macaroni penguins across the Southern Ocean

The study of systematics in wide-ranging seabirds can be challenging due to the vast geographic scales involved, as well as the possible discordance between molecular, morphological and behavioral data. In the Southern Ocean, macaroni penguins (Eudyptes chrysolophus) are distributed over a circumpolar range including populations in Antarctic and sub-Antarctic areas. Macquarie Island, in its relative isolation, is home to a closely related endemic taxon - the royal penguin (Eudyptes schlegeli), which is distinguishable from E. chrysolophus mainly by facial coloration. Although these sister taxa are widely accepted as representing distinct species based on morphological grounds, the extent of their genome-wide differentiation remains uncertain. In this study, we use genome-wide Single Nucleotide Polymorphisms to test genetic differentiation between these geographically isolated taxa and evaluate the main drivers of population structure among breeding colonies of macaroni/royal penguins. Genetic similarity observed between macaroni and royal penguins suggests they constitute a single evolutionary unit. Nevertheless, royal penguins exhibited a tendency to cluster only with macaroni individuals from Kerguelen Island, suggesting that dispersal occurs mainly between these neighboring colonies. A stepping stone model of differentiation of macaroni/royal populations was further supported by a strong pattern of isolation by distance detected across its whole distribution range, possibly driven by large geographic distances between colonies as well as natal philopatry. However, we also detected intraspecific genomic differentiation between Antarctic and sub-Antarctic populations of macaroni penguins, highlighting the role of environmental factors together with geographic distance in the processes of genetic differentiation between Antarctic and sub-Antarctic waters.

Uncovering population structure in the Humboldt penguin (Spheniscus humboldti) along the Pacific coast at South America

The upwelling hypothesis has been proposed to explain reduced or lack of population structure in seabird species specialized in food resources available at cold-water upwellings. However, population genetic structure may be challenging to detect in species with large population sizes, since variation in allele frequencies are more robust under genetic drift. High gene flow among populations, that can be constant or pulses of migration in a short period, may also decrease power of algorithms to detect genetic structure. Penguin species usually have large population sizes, high migratory ability but philopatric behavior, and recent investigations debate the existence of subtle population structure for some species not detected before. Previous study on Humboldt penguins found lack of population genetic structure for colonies of Punta San Juan and from South Chile. Here, we used mtDNA and nuclear markers (10 microsatellites and RAG1 intron) to evaluate population structure for 11 main breeding colonies of Humboldt penguins, covering the whole spatial distribution of this species. Although mtDNA failed to detect population structure, microsatellite loci and nuclear intron detected population structure along its latitudinal distribution. Microsatellite showed significant R_st values between most of pairwise locations (44 of 56 locations, R_st = 0.003 to 0.081) and 86% of individuals were assigned to their sampled colony, suggesting philopatry. STRUCTURE detected three main genetic clusters according to geographical locations: i) Peru; ii) North of Chile; and iii) Central-South of Chile. The Humboldt penguin shows signal population expansion after the Last Glacial Maximum (LGM), suggesting that the genetic structure of the species is a result of population dynamics and foraging colder water upwelling that favor gene flow and phylopatric rate. Our findings thus highlight that variable markers and wide sampling along the species distribution are crucial to better understand genetic population structure in animals with high dispersal ability.

Spatial distribution, prevalence and diversity of haemosporidians in the rufous-collared sparrow, Zonotrichia capensis

Background

Parasite prevalence and diversity are determined by the distribution of hosts and vectors and by the interplay among a suite of environmental factors. Distributions of parasite lineages vary based on host susceptibility and geographical barriers. Hemoparasites of the genera Haemoproteus and Plasmodium have wide distributions, and high prevalence and genetic diversity within perching birds (Order Passeriformes). The rufous-collared sparrow (Zonotrichia capensis) is widely distributed in Central and South America across an immense diversity of environments from sea level to more than 4000 meters above sea level. It therefore provides an excellent model to investigate whether altitudinal and latitudinal gradients influence the distribution, prevalence and diversity of haemosporidian parasites, their population structure and the biogeographical boundaries of distinct parasite lineages.

Results

We assembled samples from 1317 rufous-collared sparrows spanning 75 locales from across Central and South America (between 9.5°N and 54°S; 10–4655 meters above sea level). We used DNA sequence data from a fragment of the mitochondrial cytochrome b gene (cytb) of Haemoproteus and Plasmodium from 325 positive samples and found prevalences of 22 and 3%, respectively. Haemoproteus exhibited a higher prevalence than Plasmodium but with comparatively lower genetic diversity. We detected a relationship of Plasmodium and Haemoproteus prevalence with altitude and latitude; however, altitude and latitude did not influence parasite diversity.

Conclusions

Parasite lineages showed a phylogeographical boundary coincident with the Andes Mountains, although we also observed a north-south disjunction in Peru for Haemoproteus. Haemosporidian distribution was not homogeneous but differed based on latitude and altitude. This is most probably due to environmental factors that have influenced both vector distribution and abundance, as well as parasite development. Our study provides key insights on the distribution of haemoparasite lineages and parasite dynamics within hosts.

Electronic supplementary material

The online version of this article (10.1186/s13071-018-3243-4) contains supplementary material, which is available to authorized users.

Correction to: Chinstrap penguin population genetic structure: one or more populations along the Southern Ocean?

Correction to: BMC Evolutionary Biology (2018) 18:90 https://doi.org/10.1186/s12862-018-1207-0 .

Mercury Exposure in Humboldt (Spheniscus humboldti) and Chinstrap (Pygoscelis antarcticus) Penguins Throughout the Chilean Coast and Antarctica

Penguins are reliable sentinels for environmental assessments of mercury (Hg) due to their longevity, abundance, high trophic level, and relatively small foraging areas. We analyzed Hg concentrations from blood and feathers of adult Humboldt penguins (Spheniscus humboldti) and feathers of chinstrap penguins (Pygoscelis antarcticus) from different reproductive colonies with variable degrees of urbanization and industrialization along the Chilean and Antarctic coasts. We evaluated Hg concentration differences between species, sexes (Humboldt penguins), and localities. Our results showed significantly greater levels in Humboldt penguins than in chinstrap penguins and nonsignificant differences between sexes among Humboldts. Penguin Hg concentrations showed a latitudinal pattern, with greater values of the metal at lower latitudes, independent of the species. Both studied penguin species showed elevated Hg concentrations compared to their congeners, highlighting the necessity to investigate potential negative effects on their populations. Although differences between species are possibly due to variation in diet and trophic level, our results suggest an important effect of the degree of Hg pollution adjacent to foraging areas. Further research on Hg content in prey species and environmental samples, together with a larger overall sample size, and investigation on penguin's diet and trophic level are needed to elucidate Hg bioavailability in each location and the role of local Hg pollution levels. Likewise, it is important to monitor Hg and other heavy metals of ecotoxicological importance in penguin populations in vulnerable regions of Chile.

Chinstrap penguin population genetic structure: one or more populations along the Southern Ocean?

Background

Historical factors, demography, reproduction and dispersal are crucial in determining the genetic structure of seabirds. In the Antarctic marine environment, penguins are a major component of the avian biomass, dominant predators and important bioindicators of ecological change. Populations of chinstrap penguins have decreased in nearly all their breeding sites, and their range is expanding throughout the Antarctic Peninsula. Population genetic structure of this species has been studied in some colonies, but not between breeding colonies in the Antarctic Peninsula or at the species’ easternmost breeding colony (Bouvetøya).

Results

Connectivity, sex-biased dispersal, diversity, genetic structure and demographic history were studied using 12 microsatellite loci and a mitochondrial DNA region (HVRI) in 12 breeding colonies in the South Shetland Islands (SSI) and the Western Antarctic Peninsula (WAP), and one previously unstudied sub-Antarctic island, 3600 km away from the WAP (Bouvetøya). High genetic diversity, evidence of female bias-dispersal and a sign of population expansion after the last glacial maximum around 10,000 mya were detected. Limited population genetic structure and lack of isolation by distance throughout the region were found, along with no differentiation between the WAP and Bouvetøya (overall microsatellite F_ST = 0.002, p = 0.273; mtDNA F_ST = − 0.004, p = 0.766), indicating long distance dispersal. Therefore, genetic assignment tests could not assign individuals to their population(s) of origin. The most differentiated location was Georges Point, one of the southernmost breeding colonies of this species in the WAP.

Conclusions

The subtle differentiation found may be explained by some combination of low natal philopatric behavior, high rates of dispersal and/or generally high mobility among colonies of chinstrap penguins compared to other Pygoscelis species.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1207-0) contains supplementary material, which is available to authorized users.

Landscape genomics: natural selection drives the evolution of mitogenome in penguins

Background

Mitochondria play a key role in the balance of energy and heat production, and therefore the mitochondrial genome is under natural selection by environmental temperature and food availability, since starvation can generate more efficient coupling of energy production. However, selection over mitochondrial DNA (mtDNA) genes has usually been evaluated at the population level. We sequenced by NGS 12 mitogenomes and with four published genomes, assessed genetic variation in ten penguin species distributed from the equator to Antarctica. Signatures of selection of 13 mitochondrial protein-coding genes were evaluated by comparing among species within and among genera (Spheniscus, Pygoscelis, Eudyptula, Eudyptes and Aptenodytes). The genetic data were correlated with environmental data obtained through remote sensing (sea surface temperature [SST], chlorophyll levels [Chl] and a combination of SST and Chl [COM]) through the distribution of these species.

Results

We identified the complete mtDNA genomes of several penguin species, including ND6 and 8 tRNAs on the light strand and 12 protein coding genes, 14 tRNAs and two rRNAs positioned on the heavy strand. The highest diversity was found in NADH dehydrogenase genes and the lowest in COX genes. The lowest evolutionary divergence among species was between Humboldt (Spheniscus humboldti) and Galapagos (S. mendiculus) penguins (0.004), while the highest was observed between little penguin (Eudyptula minor) and Adélie penguin (Pygoscelis adeliae) (0.097). We identified a signature of purifying selection (Ka/Ks < 1) across the mitochondrial genome, which is consistent with the hypothesis that purifying selection is constraining mitogenome evolution to maintain Oxidative phosphorylation (OXPHOS) proteins and functionality. Pairwise species maximum-likelihood analyses of selection at codon sites suggest positive selection has occurred on ATP8 (Fixed-Effects Likelihood, FEL) and ND4 (Single Likelihood Ancestral Counting, SLAC) in all penguins. In contrast, COX1 had a signature of strong negative selection. ND4 Ka/Ks ratios were highly correlated with SST (Mantel, p-value: 0.0001; GLM, p-value: 0.00001) and thus may be related to climate adaptation throughout penguin speciation.

Conclusions

These results identify mtDNA candidate genes under selection which could be involved in broad-scale adaptations of penguins to their environment. Such knowledge may be particularly useful for developing predictive models of how these species may respond to severe climatic changes in the future.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4424-9) contains supplementary material, which is available to authorized users.

Genetic structure of introduced American mink (Neovison vison) in Patagonia: colonisation insights and implications for control and management strategies

Context Biological invasions have caused dramatic changes in native biodiversity and ecosystem function. Studies of genetic variation and evolutionary changes are useful for understanding population dynamics during biological invasions, and shed light on management, prevention and restoration strategies. Aims This study aimed to investigate the structure and genetic variability of American mink (Neovison vison), an invasive species in southern South America, introduced for fur farming in the 1930s. Methods Samples from 153 mink were obtained from 12 locations in southern Chile to sequence the mitochondrial DNA (mtDNA) control region and to genotype 11 polymorphic microsatellite loci. Key results The highest mtDNA diversity was detected in Puerto Cisnes, suggesting multiple introductions and/or the most probable area where mink was first introduced. The latter is also supported by microsatellite data, because a high percentage of individuals from different locations were assigned to this location. All other locations showed low or no mtDNA diversity, possibly due to founder effect. The results also indicate marked population structure, with three genetic clusters coincident with the main historical introduction points, with low dispersal among them. Conclusions The results suggest that control strategies for American mink in southern Chile should be concentrated on these three genetically differentiated management units, and particularly on source populations and locations with low effective population size and restricted connectivity. Implications Genetic approaches have been used for the management of numerous alien species worldwide. Recommendations delivered here for American mink control could also be implemented in other regions and for other invasive species with similar genetic diversity distribution and connectivity.

Comparative genome-wide polymorphic microsatellite markers in Antarctic penguins through next generation sequencing

Microsatellites are valuable molecular markers for evolutionary and ecological studies. Next generation sequencing is responsible for the increasing number of microsatellites for non-model species. Penguins of the Pygoscelis genus are comprised of three species: Adélie (P. adeliae), Chinstrap (P. antarcticus) and Gentoo penguin (P. papua), all distributed around Antarctica and the sub-Antarctic. The species have been affected differently by climate change, and the use of microsatellite markers will be crucial to monitor population dynamics. We characterized a large set of genome-wide microsatellites and evaluated polymorphisms in all three species. SOLiD reads were generated from the libraries of each species, identifying a large amount of microsatellite loci: 33,677, 35,265 and 42,057 for P. adeliae, P. antarcticus and P. papua, respectively. A large number of dinucleotide (66,139), trinucleotide (29,490) and tetranucleotide (11,849) microsatellites are described. Microsatellite abundance, diversity and orthology were characterized in penguin genomes. We evaluated polymorphisms in 170 tetranucleotide loci, obtaining 34 polymorphic loci in at least one species and 15 polymorphic loci in all three species, which allow to perform comparative studies. Polymorphic markers presented here enable a number of ecological, population, individual identification, parentage and evolutionary studies of Pygoscelis, with potential use in other penguin species.

Marked phylogeographic structure of Gentoo penguin reveals an ongoing diversification process along the Southern Ocean

Two main hypotheses have been debated about the biogeography of the Southern Ocean: (1) the Antarctic Polar Front (APF), acting as a barrier between Antarctic and sub-Antarctic provinces, and (2) the Antarctic Circumpolar Current (ACC), promoting gene flow among sub-Antarctic areas. The Gentoo penguin is distributed throughout these two provinces, separated by the APF. We analyzed mtDNA (HVR1) and 12 microsatellite loci of 264 Gentoo penguins, Pygoscelis papua, from 12 colonies spanning from the Western Antarctic Peninsula and the South Shetland Islands (WAP) to the sub-Antarctic Islands (SAI). While low genetic structure was detected among WAP colonies (mtDNA Ф_ST=0.037-0.133; microsatellite F_ST=0.009-0.063), high differentiation was found between all SAI and WAP populations (mtDNA Ф_ST=0.678-0.930; microsatellite F_ST=0.110-0.290). These results suggest that contemporary dispersal around the Southern Ocean is very limited or absent. As predicted, the APF appears to be a significant biogeographical boundary for Gentoo penguin populations; however, the ACC does not promote connectivity in this species. Our data suggest demographic expansion in the WAP during the last glacial maximum (LGM, about 20kya), but stability in SAI. Phylogenetic analyses showed a deep divergence between populations from the WAP and those from the SAI. Therefore, taxonomy should be further revised. The Crozet Islands resulted as a basal clade (3.57Mya), followed by the Kerguelen Islands (2.32Mya) as well as a more recent divergence between the Falkland/Malvinas Islands and the WAP (1.27Mya). Historical isolation, local adaptation, and past climate scenarios of those Evolutionarily Significant Units may have led to different potentials to respond to climate changes.

Contrasting patterns of selection between MHC I and II across populations of Humboldt and Magellanic penguins

The evolutionary and adaptive potential of populations or species facing an emerging infectious disease depends on their genetic diversity in genes, such as the major histocompatibility complex (MHC). In birds, MHC class I deals predominantly with intracellular infections (e.g., viruses) and MHC class II with extracellular infections (e.g., bacteria). Therefore, patterns of MHC I and II diversity may differ between species and across populations of species depending on the relative effect of local and global environmental selective pressures, genetic drift, and gene flow. We hypothesize that high gene flow among populations of Humboldt and Magellanic penguins limits local adaptation in MHC I and MHC II, and signatures of selection differ between markers, locations, and species. We evaluated the MHC I and II diversity using 454 next-generation sequencing of 100 Humboldt and 75 Magellanic penguins from seven different breeding colonies. Higher genetic diversity was observed in MHC I than MHC II for both species, explained by more than one MHC I loci identified. Large population sizes, high gene flow, and/or similar selection pressures maintain diversity but limit local adaptation in MHC I. A pattern of isolation by distance was observed for MHC II for Humboldt penguin suggesting local adaptation, mainly on the northernmost studied locality. Furthermore, trans-species alleles were found due to a recent speciation for the genus or convergent evolution. High MHC I and MHC II gene diversity described is extremely advantageous for the long-term survival of the species.

Molecular Epidemiology of Avian Malaria in Wild Breeding Colonies of Humboldt and Magellanic Penguins in South America

Avian malaria is a disease caused by species of the genera Haemoproteus, Leucocytozoon, and Plasmodium. It affects hundreds of bird species, causing varied clinical signs depending on the susceptibility of the host species. Although high mortality has been reported in captive penguins, limited epidemiological studies have been conducted in wild colonies, and isolated records of avian malaria have been reported mostly from individuals referred to rehabilitation centers. For this epidemiological study, we obtained blood samples from 501 adult Humboldt and 360 adult Magellanic penguins from 13 colonies throughout South America. To identify malaria parasitaemia, we amplified the mtDNA cytochrome b for all three parasite genera. Avian malaria was absent in most of the analyzed colonies, with exception of the Punta San Juan Humboldt penguin colony, in Peru, where we detected at least two new Haemoproteus lineages in three positive samples, resulting in a prevalence of 0.6% for the species. The low prevalence of avian malaria detected in wild penguins could be due to two possible causes: A low incidence, with high morbidity and mortality in wild penguins or alternatively, penguins sampled in the chronic stage of the disease (during which parasitaemia in peripheral blood samples is unlikely) would be detected as false negatives.

Comparative phylogeography of co-distributed Phrygilus species (Aves, Thraupidae) from the Central Andes

The Neotropical ecoregion has been an important place of avian diversification where dispersal and allopatric events coupled with periods of active orogeny and climate change (Late Pliocene-Pleistocene) have shaped the biogeography of the region. In the Neotropics, avian population structure has been sculpted not only by geographical barriers, but also by non-allopatric factors such as natural selection and local adaptation. We analyzed the genetic variation of six co-distributed Phrygilus species from the Central Andes, based on mitochondrial and nuclear markers in conjunction with morphological differentiation. We examined if Phrygilus species share patterns of population structure and historical demography, and reviewed the intraspecific taxonomy in part of their geographic range. Our results showed different phylogeographic patterns between species, even among those belonging to the same phylogenetic clade. P. alaudinus, P. atriceps, and P. unicolor showed genetic differentiation mediated by allopatric mechanisms in response to specific geographic barriers; P. gayi showed sympatric lineages in northern Chile, while P. plebejus and P. fruticeti showed a single genetic group. We found no relationship between geographic range size and genetic structure. Additionally, a signature of expansion was found in three species related to the expansion of paleolakes in the Altiplano region and the drying phase of the Atacama Desert. Morphological analysis showed congruence with molecular data and intraspecific taxonomy in most species. While we detected genetic and phenotypic patterns that could be related to natural selection and local adaptation, our results indicate that allopatric events acted as a major factor in the population differentiation of Phrygilus species.

Cross-amplification of nonspecific microsatellites markers: a useful tool to study endangered/vulnerable species of southern Andes deer

Thirty-nine microsatellite loci that are highly conserved in red deer, sika deer, reindeer, Soay sheep, and other artiodactyls were tested in two vulnerable and endangered Neotropical deer (pudu: Pudu puda and huemul: Hippocamelus bisulcus) with the aim of producing a standardized set of markers that can be used successfully in noninvasive samples from these species. We also compared these nonspecific loci against eight polymorphic loci that were recently developed for huemul to determine whether the nonspecific markers could reflect the huemul's genetic variation that was observed with the specific loci. We identified 10 suitable loci, six of which constitute a standardized set for the two species and can be used to identify them in the absence of phenotypic data. The expected heterozygosity per locus for the panel of six loci ranged from 0.461 to 0.889 (average 0.665), and the maximum probability of identity value was 6.9x10(-6) and 3.2x10(-4) in pudu and huemul, respectively. This set of loci has potential applications in evolutionary, ecological, forensic, and conservation studies in pudu and huemul.

Have historical climate changes affected Gentoo penguin (Pygoscelis papua) populations in Antarctica?

The West Antarctic Peninsula (WAP) has been suffering an increase in its atmospheric temperature during the last 50 years, mainly associated with global warming. This increment of temperature trend associated with changes in sea-ice dynamics has an impact on organisms, affecting their phenology, physiology and distribution range. For instance, rapid demographic changes in Pygoscelis penguins have been reported over the last 50 years in WAP, resulting in population expansion of sub-Antarctic Gentoo penguin (P. papua) and retreat of Antarctic Adelie penguin (P. adeliae). Current global warming has been mainly associated with human activities; however these climate trends are framed in a historical context of climate changes, particularly during the Pleistocene, characterized by an alternation between glacial and interglacial periods. During the last maximal glacial (LGM∼21,000 BP) the ice sheet cover reached its maximum extension on the West Antarctic Peninsula (WAP), causing local extinction of Antarctic taxa, migration to lower latitudes and/or survival in glacial refugia. We studied the HRVI of mtDNA and the nuclear intron βfibint7 of 150 individuals of the WAP to understand the demographic history and population structure of P. papua. We found high genetic diversity, reduced population genetic structure and a signature of population expansion estimated around 13,000 BP, much before the first paleocolony fossil records (∼1,100 BP). Our results suggest that the species may have survived in peri-Antarctic refugia such as South Georgia and North Sandwich islands and recolonized the Antarctic Peninsula and South Shetland Islands after the ice sheet retreat.

Phylogeography and demographic history of the neotropical otter (Lontra longicaudis)

The Neotropical otter (Lontra longicaudis) is a medium-sized semiaquatic carnivore with a broad distribution in the Neotropical region. Despite being apparently common in many areas, it is one of the least known otters, and genetic studies on this species are scarce. Here, we have investigated its genetic diversity, population structure, and demographic history across a large portion of its geographic range by analyzing 1471 base pairs (bp) of mitochondrial DNA from 52 individuals. Our results indicate that L. longicaudis presents high levels of genetic diversity and a consistent phylogeographic pattern, suggesting the existence of at least 4 distinct evolutionary lineages in South America. The observed phylogeographic partitions are partially congruent with the subspecies classification previously proposed for this species. Coalescence-based analyses indicate that Neotropical otter mitochondrial DNA lineages have shared a rather recent common ancestor, approximately 0.5 Ma, and have subsequently diversified into the observed phylogroups. A consistent scenario of recent population expansion was identified in Eastern South America based on several complementary analyses of historical demography. The results obtained here provide novel insights on the evolutionary history of this largely unknown Neotropical mustelid and should be useful to design conservation and management policies on behalf of this species and its habitats.

Priorities for the conservation of the pudu (Pudu puda) in southern South America

The southern pudu (Pudu puda) is a threatened deer that is endemic to the South American temperate forests. Despite its assumed threatened status, there is relatively little understanding on the ecology and conservation of this species. Considering this situation and the fact that there are some research groups currently working on this species, we organised a symposium to discuss research and management priorities – as well as to coordinate efforts – to move forward on the conservation of the pudu. We agreed that main research priorities should be to increase the understanding of the threats that jeopardize the viability of pudu populations, with a strong emphasis on research questions that will provide information for the management of these threats. The main management recommendations were to implement monitoring of pudu populations at least in protected areas, to implement specific actions to remove threats from protected areas and to start following internationally-accepted guidelines for the management of rescued and confiscated animals.

Phylogeography of the Marine Otter (Lontra felina): historical and contemporary factors determining its distribution

The evolutionary history of a species can be revealed by phylogeographical analysis; nevertheless, not only historical but also contemporary processes can imprint on the distribution of genetic diversity. We report on the phylogeny of Lontra ssp. in South America, and the role of spatial heterogeneity in shaping the distribution and population structure of the endangered marine otter, Lontra felina. Analyzing a total of 2261 bp of mitochondrial DNA (mtDNA) revealed the recent divergence of L. felina from L. provocax. A strong population structure (Φ(st) = 0.83, P < 0.0001) and a significant pattern of isolation by distance were described for L. felina (n = 168) across a wide geographical distribution (13°53'S to 43°36'S). Lontra felina mtDNA phylogeny is composed of 2 main clades: a clade from Peru and another composed of Chilean haplotypes. Northern populations show different divergent lineages and higher genetic diversity when compared with more recently colonized southern populations. Furthermore, long sandy beaches seem to act as barriers to dispersal, creating 2 evolutionary significant units in agreement with subspecies previous description, and at least 5 different management units (MUs). At a fine spatial scale, the size of rocky seashore patches, the distance between patches and anthropogenic factors also play important roles in species gene flow.

Phylogeography, phylogeny and hybridization in trichechid sirenians: implications for manatee conservation

Abstract The three living species of manatees, West Indian (Trichechus manatus), Amazonian (Trichechus inunguis) and West African (Trichechus senegalensis), are distributed across the shallow tropical and subtropical waters of America and the western coast of Africa. We have sequenced the mitochondrial DNA control region in 330 Trichechus to compare their phylogeographic patterns. In T. manatus we observed a marked population structure with the identification of three haplotype clusters showing a distinct spatial distribution. A geographic barrier represented by the continuity of the Lesser Antilles to Trinidad Island, near the mouth of the Orinoco River in Venezuela, appears to have restricted the gene flow historically in T. manatus. However, for T. inunguis we observed a single expanding population cluster, with a high diversity of very closely related haplotypes. A marked geographic population structure is likely present in T. senegalensis with at least two distinct clusters. Phylogenetic analyses with the mtDNA cytochrome b gene suggest a clade of the marine Trichechus species, with T. inunguis as the most basal trichechid. This is in agreement with previous morphological analyses. Mitochondrial DNA, autosomal microsatellites and cytogenetic analyses revealed the presence of hybrids between the T. manatus and T. inunguis species at the mouth of the Amazon River in Brazil, extending to the Guyanas and probably as far as the mouth of the Orinoco River. Future conservation strategies should consider the distinct population structure of manatee species, as well as the historical barriers to gene flow and the likely occurrence of interspecific hybridization.