Comparative Genomics Supports Ecologically Induced Selection as a Putative Driver of Banded Penguin Diversification.

The relative importance of genetic drift and local adaptation in facilitating speciation remains unclear. This is particularly true for seabirds, which can disperse over large geographic distances, providing opportunities for intermittent gene flow among distant colonies that span the temperature and salinity gradients of the oceans. Here, we delve into the genomic basis of adaptation and speciation of banded penguins, Galápagos (Spheniscus mendiculus), Humboldt (Spheniscus humboldti), Magellanic (Spheniscus magellanicus), and African penguins (Spheniscus demersus), by analyzing 114 genomes from the main 16 breeding colonies. We aim to identify the molecular mechanism and genomic adaptive traits that have facilitated their diversifications. Through positive selection and gene family expansion analyses, we identified candidate genes that may be related to reproductive isolation processes mediated by ecological thermal niche divergence. We recover signals of positive selection on key loci associated with spermatogenesis, especially during the recent peripatric divergence of the Galápagos penguin from the Humboldt penguin. High temperatures in tropical habitats may have favored selection on loci associated with spermatogenesis to maintain sperm viability, leading to reproductive isolation among young species. Our results suggest that genome-wide selection on loci associated with molecular pathways that underpin thermoregulation, osmoregulation, hypoxia, and social behavior appears to have been crucial in local adaptation of banded penguins. Overall, these results contribute to our understanding of how the complexity of biotic, but especially abiotic, factors, along with the high dispersal capabilities of these marine species, may promote both neutral and adaptive lineage divergence even in the presence of gene flow.

Climate change and population persistence in a hibernating marsupial.

Climate change has physiological consequences on organisms, ecosystems and human societies, surpassing the pace of organismal adaptation. Hibernating mammals are particularly vulnerable as winter survival is determined by short-term physiological changes triggered by temperature. In these animals, winter temperatures cannot surpass a certain threshold, above which hibernators arouse from torpor, increasing several fold their energy needs when food is unavailable. Here, we parameterized a numerical model predicting energy consumption in heterothermic species and modelled winter survival at different climate change scenarios. As a model species, we used the arboreal marsupial monito del monte (genus Dromiciops), which is recognized as one of the few South American hibernators. We modelled four climate change scenarios (from optimistic to pessimistic) based on IPCC projections, predicting that northern and coastal populations (Dromiciops bozinovici) will decline because the minimum number of cold days needed to survive the winter will not be attained. These populations are also the most affected by habitat fragmentation and changes in land use. Conversely, Andean and other highland populations, in cooler environments, are predicted to persist and thrive. Given the widespread presence of hibernating mammals around the world, models based on simple physiological parameters, such as this one, are becoming essential for predicting species responses to warming in the short term.

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.

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.

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.

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.

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.

Widespread Infection with Hemotropic Mycoplasmas in Free-Ranging Dogs and Wild Foxes Across Six Bioclimatic Regions of Chile.

Blood samples of 626 rural dogs, 140 Andean foxes (Lycalopex culpaeus), and 83 South American grey foxes (L. griseus) from six bioregions of Chile spanning 3000 km were screened for Mycoplasma DNA by conventional PCR and sequencing. Risk factors of infection were inferred using Generalized Linear Mixed Models and genetic structure by network analyses. Overall, Mycoplasma haemocanis/Mycoplasma haemofelis (Mhc/Mhf) and Candidatus Mycoplasma haematoparvum (CMhp) observed prevalence was 23.8% and 12.8% in dogs, 20.1% and 7.2% in Andean foxes, and 26.5% and 8.4% in grey foxes, respectively. Both hemoplasmas were confirmed in all the bioregions, with higher prevalence in those where ticks from the Rhipicephalus sanguineus species group were absent. Candidatus M. haematominutum and a Mycoplasma sp. previously found in South American carnivores were detected in one fox each. Although the most prevalent Mhc/Mhf and CMhp sequence types were shared between dogs and foxes, network analysis revealed genetic structure of Mhc/Mhf between hosts in some regions. Male sex was associated with a higher risk of Mhc/Mhf and CMhp infection in dogs, and adult age with CMhp infection, suggesting that direct transmission is relevant. No risk factor was identified in foxes. Our study provides novel information about canine hemoplasmas with relevance in distribution, transmission routes, and cross-species transmission.

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.

Utility of genetic variation in coat color genes to distinguish wild, domestic and hybrid South American camelids for forensic and judicial applications.

A molecular genetic protocol for distinguishing pure and hybrid South American camelids was developed to provide strong, quantifiable, and unbiased species identification. We detail the application of the approach in the context of a criminal case in the Andes Mountains of central Chile where the defendants were alleged to have illegally hunted three wild guanacos (Lama guanicoe), as opposed to hybrid domestic llama (Lama glama)/wild guanaco crosses, which are unregulated. We describe a workflow that differentiates among wild, domestic and hybrid South American camelids (Lama versus Vicugna) based on mitochondrial cytochrome b genetic variation (to distinguish between Lama and Vicugna), and MC1R and exon 4 variation of the ASIP gene (to differentiate wild from domestic species). Additionally, we infer the population origin and sex of each of the three individuals from a panel of 15 autosomal microsatellite loci and the presence or absence of the SRY gene. Our analyses strongly supported the inference that the confiscated carcasses corresponded with 2 male and 1 female guanacos that were hunted illegally. Statistical power analyses suggested that there was an extremely low probability of misidentifying domestic camelids as wild camelids (an estimated 0 % Type I error rate), or using more conservative approached a 1.17 % chance of misidentification of wild species as domestic camelids (Type II error). Our case report and methodological and analytical protocols demonstrate the power of genetic variation in coat color genes to identify hybrids between wild and domestic camelid species and highlight the utility of the approach to help combat illegal wildlife hunting and trafficking.

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.

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.

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

Abstract 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.