History of Diversification and Adaptation from North to South Revealed by Genomic Data: Guanacos from the Desert to Sub-Antarctica.

The increased availability of quality genomic data has greatly improved the scope and resolution of our understanding of the recent evolutionary history of wild species adapted to extreme environments and their susceptibility to anthropogenic impacts. The guanaco (Lama guanicoe), the largest wild ungulate in South America, is a good example. The guanaco is well adapted to a wide range of habitats, including the Sechura Desert, the high Andes Mountains to the north, and the extreme temperatures and conditions of Navarino Island to the south. Guanacos also have a long history of overexploitation by humans. To assess the evolutionary impact of these challenging habitats on the genomic diversity, we analyzed 38 genomes (∼10 to 16×) throughout their extensive latitudinal distribution from the Sechura and Atacama Desert to southward into Tierra del Fuego Island. These included analyses of patterns of unique differentiation in the north and geographic region further south with admixture among L. g. cacsilensis and L. g. guanicoe. Our findings provide new insights on the divergence of the subspecies ∼800,000 yr BP and document two divergent demographic trajectories and to the initial expansion of guanaco into the more southern portions of the Atacama Desert. Patagonian guanacos have experienced contemporary reductions in effective population sizes, likely the consequence of anthropogenic impacts. The lowest levels of genetic diversity corresponded to their northern and western limits of distribution and some varying degrees of genetic differentiation. Adaptive genomic diversity was strongly linked with environmental variables and was linked with colonization toward the south followed by adaptation.

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.

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.

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.

The systematics of Tachymenini (Serpentes, Dipsadidae): an updated classification based on molecular and morphological evidence

Despite the recent advances in the systematics of snakes, the diversity of several Neotropical groups of species remains poorly understood. The lack of studies focused on the phylogenetic relationship within most of the 20 tribes of Dipsadidae precludes a better understanding of the evolution of this diverse family. Here, we present a comprehensive phylogenetic analysis of Tachymenini, a heterogeneous tribe that comprises 36 viviparous species of dipsadids, grouped in seven genera based on morphological similarities. The tribe is widely distributed throughout South America presenting very distinctive phenotypes, habitats, and behaviours. The phylogenetic relationship among tachymenins is a well-recognized challenge regarding the systematics of the South American dipsadids. The similar morphotype of some generalist species, combined with the very derived morphology of some strict specialists, creates a complexity of traits that has prevented the comprehension of the systematics of the group. To address such a challenge, we combine molecular (six loci) and morphological (70 characters) datasets in an integrative phylogenetic approach. The resultant phylogenetic trees indicate, with strong support, that three of the seven current recognized genera (Tachymenis, Tomodon, and Thamnodynastes) are non-monophyletic and, consequently, we propose a new systematic arrangement for Tachymenini. We revalidate two genera, Dryophylax and Mesotes, and we describe three additional monotypic genera, Apographon gen. n., Tachymenoides gen. n., and Zonateres gen. n. to accommodate Tomodon orestes, Tachymenis affinis, and Thamnodynastes lanei, respectively. We also include Tomodon ocellatus and Pseudotomodon trigonatus in Tachymenis and describe a new genus, Galvarinus gen. n., to accommodate the Tachymenis chilensis species group. Furthermore, we also provide an evolutionary scenario for the speciation events based on a time-calibrated tree, commenting on the diversification and origin of the tribe, and on the probable existence of undescribed species of Mesotes and Dryophylax.

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.

Genome-wide analyses reveal drivers of penguin diversification.

Penguins are the only extant family of flightless diving birds. They currently comprise at least 18 species, distributed from polar to tropical environments in the Southern Hemisphere. The history of their diversification and adaptation to these diverse environments remains controversial. We used 22 new genomes from 18 penguin species to reconstruct the order, timing, and location of their diversification, to track changes in their thermal niches through time, and to test for associated adaptation across the genome. Our results indicate that the penguin crown-group originated during the Miocene in New Zealand and Australia, not in Antarctica as previously thought, and that Aptenodytes is the sister group to all other extant penguin species. We show that lineage diversification in penguins was largely driven by changing climatic conditions and by the opening of the Drake Passage and associated intensification of the Antarctic Circumpolar Current (ACC). Penguin species have introgressed throughout much of their evolutionary history, following the direction of the ACC, which might have promoted dispersal and admixture. Changes in thermal niches were accompanied by adaptations in genes that govern thermoregulation and oxygen metabolism. Estimates of ancestral effective population sizes (Ne ) confirm that penguins are sensitive to climate shifts, as represented by three different demographic trajectories in deeper time, the most common (in 11 of 18 penguin species) being an increased Ne between 40 and 70 kya, followed by a precipitous decline during the Last Glacial Maximum. The latter effect is most likely a consequence of the overall decline in marine productivity following the last glaciation.

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.

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.

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.

Receding ice drove parallel expansions in Southern Ocean penguins.

Climate shifts are key drivers of ecosystem change. Despite the critical importance of Antarctica and the Southern Ocean for global climate, the extent of climate-driven ecological change in this region remains controversial. In particular, the biological effects of changing sea ice conditions are poorly understood. We hypothesize that rapid postglacial reductions in sea ice drove biological shifts across multiple widespread Southern Ocean species. We test for demographic shifts driven by climate events over recent millennia by analyzing population genomic datasets spanning 3 penguin genera (Eudyptes, Pygoscelis, and Aptenodytes). Demographic analyses for multiple species (macaroni/royal, eastern rockhopper, Adélie, gentoo, king, and emperor) currently inhabiting southern coastlines affected by heavy sea ice conditions during the Last Glacial Maximum (LGM) yielded genetic signatures of near-simultaneous population expansions associated with postglacial warming. Populations of the ice-adapted emperor penguin are inferred to have expanded slightly earlier than those of species requiring ice-free terrain. These concerted high-latitude expansion events contrast with relatively stable or declining demographic histories inferred for 4 penguin species (northern rockhopper, western rockhopper, Fiordland crested, and Snares crested) that apparently persisted throughout the LGM in ice-free habitats. Limited genetic structure detected in all ice-affected species across the vast Southern Ocean may reflect both rapid postglacial colonization of subantarctic and Antarctic shores, in addition to recent genetic exchange among populations. Together, these analyses highlight dramatic, ecosystem-wide responses to past Southern Ocean climate change and suggest potential for further shifts as warming continues.

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 Rst values between most of pairwise locations (44 of 56 locations, Rst = 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.

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 .

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.

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.

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.