Landscape genomics: natural selection drives the evolution of mitogenome in penguins
© The Author(s). 2018
Received: 5 April 2017
Accepted: 28 December 2017
Published: 16 January 2018
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.
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.
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.
The mitochondrial genome (mtDNA) is under continuous selection since the 13 protein-coding genes produce polypeptide products that work in association with nuclear-encoded subunits of protein complexes involved in oxidative phosphorylation (OXPHOS) . Three types of hypotheses have been proposed to explain mitogenome evolution: 1) mitochondrial OXPHOS generates both energy and heat, and environmental temperature imposes a selective trade off between ATP production and heat [2–6]; 2) risk of starvation is associated with restricted food availability, since starvation can generate more efficient coupling of energy production by OXPHOS pathway [7–10]; 3) a more recent and debated hypothesis of immune response to pathogens [11–15]. Therefore, chlorophyll (Chl), as a measurement of primary productivity , and sea surface temperature (SST), are important environmental variables for marine organisms that are directly related to the efficiency of the OXPHOS and relevant to understanding mtDNA gene selection patterns. As OXPHOS is the terminal stage of cellular respiration, disruption of this process will be lethal for individual cells and if it occurs widely in an individual, it will adversely impact its fitness . Consequently, the mitochondrial genome on one hand evolves under continuous purifying selection (or negative selection) that removes deleterious variation from populations , while on the other hand, beneficial variation in extreme temperature environments may be under positive selection, increasing in frequency in a population or species. It is possible that selection may operate similarly among species in similar environments or it may drive adaptations in different genes.
Patterns of adaptation in the mtDNA of organisms have been evaluated mostly at the population level by associating each population or lineage with local climatic conditions [8, 18–21]; however, only relatively few have incorporated environmental data in the genetic analyses and generally only using nuclear genes or whole-genome data [22, 23]. Environmental temperature is the main variable that has been linked to selection in mtDNA genes [8, 18–20], along with altitude [24–26] and oxygen availability . Although most studies on the role of selection on mtDNA genome evolution in vertebrate species have focused on mammals at the intraspecific level [28–31], recently there has been an increased focus on birds [32–34]. Here, we provide one of the first species-level analyses in birds linking patterns of mtDNA selection with environmental variables (temperature, chlorophyll) which are associated with the evolution of this genome. Chlorophyll (Chl) is an excellent indirect measurement of food availability because primary productivity is associated with species distributions at trophic levels, including krill and penguins .
Penguins inhabit a broad latitudinal distribution in the Southern Hemisphere from the tropics to the extreme cold of the Antarctic. While all penguin species can live and forage at sea in cold, nutrient-rich water (including those breeding on the Galapagos Islands), they breed on land under varying climatic conditions. The northernmost species is the Galapagos penguin (Spheniscus mendiculus), while the southernmost distributed species, the Adélie penguin (P. adeliae) and Emperor penguin (Aptenodytes forsteri), breed in the Antarctic, south of the Antarctic Convergence (Polar Front). Emperor penguins breed on stable fast ice during the Antarctic winter when temperatures drop regularly to −50 °C. As the evolution, speciation and extinction of penguins are closely related to historical climatic changes [35–38], the broad distribution of penguins, which includes different and extreme climatic conditions, provides a unique opportunity to identify molecular genetic patterns associated with these climatic conditions.
Historical climate changes have been linked with penguin species diversification. Severe declines in Antarctic temperature that began approximately 12 mya coincide with the divergence of the four major penguin groups, around 11–16 mya . This penguin phylogeny is well characterized, making it possible to assess evidence of selection and adaptive evolution of the mtDNA and to explore possible links with environmental changes, including genes involved in penguin thermal adaptation. Sequences of entire mitochondrial genomes are available for five of the 18 extant species of penguins [47–50], including the Adélie, Chinstrap, African, Little, and Southern Rockhopper (E. chrysocome) penguins and more recently the entire genome became available for the Adélie and Emperor penguin . The mtDNA genome of P. antarcticus consists of 15,972 bp and has 94% identity with its sister species P. adeliae . Subramanian et al.  compared the mitochondrial genomes of modern and ancient (up to 44,000 years old) P. adeliae, documenting higher evolutionary rates than described earlier by phylogenetic studies. However, those studies evaluated the genome without a broader comparative approach to understand the evolution of the group and assess the evidence of signatures of selection on the mitochondrial genome.
We hypothesize that selection operates differentially on the mtDNA coding genes of penguin species in response to environmental variables. We evaluated variation in mitochondrial genomes at different taxonomic levels in seven species of penguins from South America and Africa (all four Spheniscus species) and Antarctica (all three Pygoscelis species). We then compared our results with two penguin species from temperate and sub-Antarctic regions (E. minor and E. chrysocome) and with A. forsteri, which survives and breeds during the most-extreme cold conditions of the Antarctic winter. We evaluated evolutionary rates and signatures of natural selection for all 13 mtDNA protein-coding genes and across the mtDNA genome phylogeny for the penguin species. Lastly, we assessed the genetic data for each gene against the spatial environmental data (sea surface temperature, chlorophyll and both variables combined) from remote sensing and GIS databases for each species along its distributional range.
Body Mass (kg)
Southern Rockhopper Penguin
Details of the individuals sequenced
Total genomic reads
Genome lenght (bp)
Antarctica, South Shetland Islands, Narebski base
Antarctic Peninsula, Kopaitic
Elephant Island, South Shetlands Islands
Antarctic Peninsula, Gabriel Gonzalez Videla base
Antarctica, South Shetland Islands, Arctowski base
Antarctica, Ross Sea, Cape Crozier
Galapagos Islands, Isabela Island
South America, Peru, Punta San Juan
South America, Chile, Cachagua
South America, Chile, Chiloé Island, Puñihuil
South America, Chile, Magdalena Island
South America, Argentina, Puerto Deseado
Genomic DNA was isolated from blood samples obtained from all penguins and preserved in ethanol using the salt method . DNA from six individuals of the three species was quantified and quality checked with Invitrogen’s PicoGreen®assay kit (fluorometry). Genomic sequencing in ABISOLiD™ 5500 XL was performed at Omics Solution, a Next Generation Sequencing facility. DNA was desalted then concentrated using standard EtOH/sodium acetate precipitation at 20 °C for 2 h, followed by two 70% EtOH washes. Then the pooled DNA was re-dissolved in low TE as per standard protocol for ABI SOLiD sequencing of genomic DNA fragment libraries. DNA samples were sheared in Covaris™S220 System, which fragments the input DNA, by sonication, into small fragments with a mean size of 160 bp.
The fragmented DNA was then purified with SOLiD Library Column Purification Kit, and libraries were prepared according to standard SOLiD protocols. Fragment libraries for the twelve penguins were prepared separately, P1 and P2 adapters were ligated, and each sample was tagged with a different barcode (a known sequence of ten bp in adapters). Prepared libraries were quantified by real-time PCR on a Light Cycler®Nano (Roche) using Invitrogen’s Quantification Kit for SOLiD. The double stranded library was added at a concentration of 0.2 pg/μL to the emulsion with 2400 million beads, according to manufacturers’ instructions. Thirty percent of the beads were P2 positive (contained amplified library fragments) before enrichment and 90% of the beads were P2 positive after enrichment, yielding 790 million beads deposited on the Flow Chip. Library beads were sequenced on a SOLiD™ 5500 XL using standard chemistry for paired-end fragment libraries and 75–35-bp read lengths.
SOLiD sequence alignment
The color-space reads (di-base encoded) were aligned with LifeScope™ (Applied Biosystems) using the genome assembly of P. adeliae  as a reference. The reference was translated into color-space with the aim of mapping the reads. The color-space reads helped to improve the quality of each base call, since each base was read twice during the sequencing step.
The consensus sequence was built from the binary alignment map (BAM) files obtained in the previous step. We used SAMtools  to obtain all bases mapped to each position, BCF tools to get the most probable genotype per position, and VCF utils to build the consensus sequence in FASTQ format. The FASTQ file was then converted to FASTA using SEQTK.
Mitochondrial sequence assembly
The mitochondrial genomes were assembled with the Denovo2 pipeline developed by Life Technologies using Velvet assembler version 1.2.09  in color-space mode. ASiD (Assembly Assistant for SOLiD) was used to fill gaps between contigs that formed scaffolds and to perform a color-space to base-space translation for the final assembly. Individual genes were selected from each of the penguins’ contigs using Sequencher software 5.1 (Gene Codes Corp., MI) and annotated by homology with genes of mitogenomes of the penguins previously published [48, 51, 56, 57]. Due to the low number of reads from one individual of P. papua (Table 2), this individual was not included in the data analyses. All 11 mitogenome were submitted to GenBank (accession numbers KM891593; KU361803-KU361806; KU356673-KU356677).
Molecular evolution and selection
Sequences were aligned and polymorphic sites were confirmed by visual review of the chromatograms using Sequencher v. 5.1. (Gene Codes, Ann Arbor, MI, USA) and multiple alignments were done using ClustalX v. 2.1 . For interspecific data analysis we compared the eleven mtDNA genomes with four other penguins genomes: S. demersus obtained from KwaZulu-Natal Province of South Africa (African penguin) , E. minor obtained from Nelson in New Zealand (Little blue penguin) , E. chrysocome (Rockhopper penguin)  and A. forsteri (Emperor penguin) .
The number of indels, polymorphic sites (S), nucleotide diversity (π), synonymous mutation rates (Ks), nonsynonymous mutation rates (Ka) and the Ka/Ks (or ω) ratio using the DNAsp v. 5  were estimated for all 13 protein-coding genes from the complete mtDNA. For analyses of the reading frame of the genes, such as Ka/Ks, only the ND3 gene was considered a different reading frame as described by Mindell et al. . All Ka and Ks estimates and corresponding ratios were obtained between pairwise species within Spheniscus and Pygoscelis genera, as well as between all five genera (Additional file 1 Table S1). Since we sequenced several species of Spheniscus and Pygoscelis genera, the values between genera were averages of the pairwise species Ka and Ks values for the genera comparison. We also calculated the pairwise distance between all 15 penguins studied here (10 species) for all 13 concatenated genes of the mtDNA genomes using the Maximun Likelihood model (Additional file 2 Table S2) with Mega v. 6.06 software [61, 62] and Arlequin software in R package [63, 64]. All positions containing gaps were eliminated, remaining a total of 10,874 positions in the final dataset.
Phylogenetic reconstruction was used to evaluate selection on each codon, gene and phylogenetic branch. We combined each of the 13 protein-coding genes from the mtDNA genome from the studied penguin species into a concatenated data set of 11,409 bp. The phylogeny was also estimated excluding the ND6 gene (a total of 10,886 bp), which is encoded by the light-strand  and has a significantly different base composition from the heavy-chain . Phylogenetic relationships were inferred using Bayesian inference methods (BA) with Markov chain Monte Carlo (MCMC) sampling as implemented in MrBayes v. 3.1.2 [66, 67]. The evolutionary model selected GTR + I + G, using jModelTest v. 0.1.1  and Akaike Information Criterion (AIC). The Mr. Bayes analyses were run for one million generations while sampling every 100 generations. At this point the standard deviation of split frequencies was <0.01, indicating that both runs had converged. Phylogenetic reconstructions were visualized using Figtree v. 1.2.2 . Although positive selection is frequently inferred by genetic data analysis as Ka/Ks > 1 (and purifying selection by Ka/Ks < 1) , a gene with Ka/Ks lower than 1 may show signature of positive selection at codon level. Evidence for signatures of selection at codon sites was evaluated for each gene with the concatenated data set using the BA phylogeny and the selected model of evolution without outgroups with the three algorithms available in HyPhy software [71–73]. For comparison, we used the Single Likelihood Ancestral Counting (SLAC), Fixed-Effects Likelihood (FEL) and the Random Effects Likelihood (REL) methods using the best fit nucleotide model . The single likelihood ancestor counting (SLAC) method is based on weighting across all possible ancestral reconstructions, or employs sampling from ancestral reconstructions. The fixed effects likelihood (FEL) method estimates nonsynonymous and synonymous substitution rates at each site. Random effects likelihood (REL) models variation in nonsynonymous and synonymous rates across sites based on a predefined distribution, where a Bayesian approach are used to infer selection pressure at an individual site . We also performed the branch Site REL  analysis in HyPhy software to estimate branch-specific episodes of positive selection along the phylogeny.
Environmental data and species distribution
Seasonal Chlorophyll concentrations (Chl) and Sea Surface Temperatures (SST) were acquired from Sea-viewing Wide Field-of-view Sensor (SeaWiFS) and Moderate Resolution Imaging Spectroradiometer (MODIS-Aqua) global seasonal databases that are based on a spatial resolution of 9 × 9 km . The standard SeaWiFS chlorophyll product, ChlOC4, is estimated based on an empirical maximal band-ratio algorithm that uses remote sensing reflectance measurements at 443, 490, 510 and 555 nm . This algorithm relates the maximum of three remote sensing reflectance band-ratios (443/555, 490/555 and 510/555) to field measures of sea chlorophyll concentration  SST is based on MODIS-calibrated mid- and far-infrared radiances (IR), estimated using an algorithm that exploits differences in atmospheric transmissivity in the different IR bands . The reliability of SST data is supported by its high correlation with sea surface temperatures measured during both oceanographic cruises and at local stations, with a slight positive bias at higher temperatures . Here we use Chl and SST as proxies for the environmental conditions along the distributional range of each penguin species.
Mantel test and GLM models
Mantel tests and GLM models (Generalized Linear Model) were implemented to explore the relationships between distance matrices of genetic parameters (π and ka/ks) and environmental parameters (Chl, SST and COM) as explanatory variables [81, 82]. To study the relationship between genetic patterns and the environmental parameters (the distances across Chl, SST and COM dimension), the Mantel matrix distances were divided into three sub-matrices (Example in Additional file 3 Table S3), each describing the patterns between pairs of species within a bounded interval of Chl, SST and COM distances . These bounded intervals consisted of a) short pairwise distances corresponding with the minimum to 33% of the data distance, b) middle distances corresponding with pairwise distances between 33% and 66% of the data distance, c) large distances consisting of the largest 33% of the pairwise distances and d) all pairwise data. For Mantel tests, p-values were obtained based on 10,000 simulations. Normality of the data (Ka/Ks and π) was tested applying the Shapiro-Wilk test. The GLMs p-values and AICs were calculated considering a Gaussian error structure for all data under normal distribution and gamma error structure for data without normal distribution. All statistical analyses were implemented using R packages .
MtDNA genomes of penguins
The number of reads for each of the 12 genomes ranged from 52,990,428 to 158,565,602. One individual of P. papua was excluded because of insufficient reads (7,341,640). The number of reads attributed to mtDNA genomes ranged from 5049 to 25,947 (1647 for P. papua) and the percent coverage varied from 16 to 82 (Table 2). We recovered the full sequence from all mtDNA regions except the repetitive region of the D-loop.
These are the first reported near-complete mtDNA genomes of S. humboldti, S. magellanicus, S. mendiculus and P. papua. The complete mitochondrial genome, without the D-loop, varied between 15,514 bp and 15,520 bp among the three Spheniscus ssp. and between 15,490 bp and 15,500 bp among the three Pygoscelis species. Including the partial D-loop, lengths varied between 16,367 bp and 16,660 bp among the three Spheniscus species, and between 16,653 bp and 16,765 bp among the three Pygoscelis species. These mt-genome lengths are comparable to the mtDNA genomes of P. antarcticus (15,972 bp), S. demersus (17,346 bp), E. minor (17,611 bp), E. chrysocome (16,930 bp), and P. adeliae (17,486 bp).
The ND6 and 8 tRNAs are transcribed from the light strand, while the other 12 protein coding genes and 14 tRNAs and two rRNAs are positioned on the heavy strand. This pattern of transcription has been described for P. adeliae, S. demersus and other bird species [49, 56, 83]. Indels, when comparing all ten species, were observed in five of the thirteen mtDNA genes (ATP8, CYTB, ND4, ND5 and ND6). All indels comprised three consecutive bases, such as for Pygoscelis and A. forsteri at ATP8; for E. chrysocome at CYTB; for A. forsteri at ND4 and ND6; and for E. minor at ND5.
Interspecific diversity across whole data set (all 15 penguins of 10 species) for all 13 protein-coding genes
Among all species
Gene Complex Length (bp)
Gene evolution and natural selection
Markedly low values of Ka/Ks ratios were found for within and between pairwise genera comparisons and were observed for COX1, followed by COX2, COX3 and ND1, suggesting negative selection (Fig. 4). In the intragenus comparisons, the Ka/Ks ratios were zero for COX1, ND3 and ND4L among all Spheniscus species and were also low for COX2, ND1 and COX3. Similarly, within Pygoscelis comparisons, Ka/Ks ratios were lowest for COX1, followed by COX3, COX2, ND1 and ATP6. Therefore, COX1 had a signature of purifying selection for amino acid sequence conservation for all pairwise comparisons among penguin genera and species.
Chlorophyll concentration (Fig. 2c) was high in the area occupied by S. humboldti, and varied highly among seasons (2.2 mg/m3 and 4.3 mg/m3 in winter and summer respectively). Chlorophyll concentrations in the distributional range of S. demersus were also high, but more stable (from 2.3 mg/m3 to 2.5 mg/m3). Highest SSTs were recorded throughout the range of S. mendiculus (32 °C and 37 °C during winter and summer respectively), followed by lower SSTs throughout areas occupied by S. demersus and S. humboldti (27.4–32.8 °C; and 23.6–29.6 °C during winter and summer respectively), with the lowest SST in areas where P. adelie, A. fosteri and P. antarcticus are distributed (0.5 °C and 5.2 °C during winter and summer respectively) (Fig. 2d).
Significant values of GLM and Mantel Test for Ka/Ks ratio for each gene and environmental data
Nucleotide diversity among penguin species was not correlated with the environmental variables (no significant value for Mantel test), nevertheless Ka/Ks ratios were significantly correlated. This supports the premise that broadly, genetic diversity of mtDNA genes are conserved, but that specific signals of selection related to environmental variables could be identified by comparative estimation of Ka/Ks ratios among genes.
Our analyses revealed strong and contrasting signals and selection patterns in the coding mtDNA genes of closely related penguin species, some of which correlated with patterns of sea surface water temperatures and chlorophyll concentrations. We identified a widespread signature of purifying selection (Ka/Ks < 1) across the mitochondrial genome, consistent with expected patterns if purifying selection were constraining the mitogenome evolution to maintain the OXPHOS proteins and functionality [20, 70]. These patterns included a signature of purifying selection in COX genes, especially COX1, which had Ka/Ks ratios close to zero for all comparison among species and genera. In contrast, one codon under positive selection was evidenced for ND4 and one for ATP8 (SLAC and FEL analysis respectively). Sea surface temperature (SST) was strongly correlated with the Ka/Ks ratio of ND4 and chlorophyll concentration was correlated with Ka/Ks ratios of COX1. The Ka/Ks ratio patterns for ATP8 varied widely among genera and there was only weak or no evidence of a correlation between genetic diversity of ATP8 and measured environmental factors. These results identify mtDNA candidate genes under selection which could be involved in broad-scale adaptations of penguins to their environment.
Evidence linking selection and gene function
In metazoans, mitogenome complexes I through V are physically positioned relative to each other in accordance with their function, with the complex II genes residing in the nucleus . A positively selected codon was detected by SLAC model on the mitogenome of penguins in the Complex I or NADH dehydrogenase gene (ND4). Complex I through IV generate the proton gradient in the intermembrane space of mitochondria. Complex I is the main entrance point of electrons and the largest complex of the respiratory chain. It catalyzes the transfer of two electrons from NADH to ubiquinone and translocates four protons across the inner mitochondrial membrane . Complex I consists of different protein subunits, with fourteen homologous subunits present in all organisms, and which constitute the enzyme.
ND1 and ND2 are situated at the connection between the peripheral and the membrane arm, whereas ND4 and ND5 are at the distal end [86, 87]. If ND2, ND4 and ND5 are the actual proton pumping devices , mutations in these genes may alter the efficiency of the proton-pumping process. The concentration of codons under positive selection in specific domains of the mitogenome are very likely to be related to protein function [20, 88] and a meta-analysis of 237 vertebrates showed that most of the genes under positive selection were located in Complex I . NADH dehydrogenase genes have also shown signatures of positive selection at a population level in birds (Parus spp.) in the ND2 gene , in lineages of an Australian bird (Eopsaltria australis)  in ND4 and ND4L. Similarly the Chinese snub-nosed monkey (Rinopithecus ssp.) showed evidence of adaptive variation related with physiological hypoxia and cold stress in altitude in ND2 and ND6 . Few studies have assessed if selection on mtDNA genes is correlated with environmental variables . In the Atlantic herring (Clupea harengus), significant positive selection was observed in ND2, ND4 and ND5. Correlations between genetic diversity and variation of environmental factors (temperature, salinity and latitude) were generally weak, although the authors addressed some correlation between temperature and genetic distances with ND3 and ND4L .
In our study, the Ka/Ks ratio for ATP8 was lower than 1, suggesting negative selection for the gene, although it was higher than any other value for mtDNA genes in penguins. However, one codon of ATP8 (Complex V) showed evidence of positive selection (under model FEL). Complex V (or ATP synthase) is the final enzyme in the OXPHOS pathway and is ubiquitous throughout all forms of life and biological functions. The enzyme uses the energy stored in a proton gradient across a membrane to drive the synthesis of ATP from ADP and phosphate. Complex V uses the proton gradient resulting from the activity of the respiratory chain enzymes complexes I, III and IV for ATP synthesis, generating 95% of cellular ATP. Complex V is a multi-subunit complex with two functional domains, F1 and F0, for which ATP8 encodes a core subunit of the F0 component. Across mammalian evolution, ATP6 is more conserved, in contrast to the ATP8 subunit which has some highly variable sites with extreme amino acid changes of properties, suggesting that there are differences in its regulatory function among species .
Cytochrome c oxidase
COXI, complex IV, showed purifying selection in penguins, as has been observed across all Aves . We estimated Ka/Ks ratios close to zero for all comparison among species and genera, and similar results were obtained using both the SLAC and FEL methods. These genes from complex IV represent the terminal complex of the respiratory chain. It catalyzes electron transfer from cytochrome c to molecular oxygen, reducing the latter to water. It has been postulated that Complex IV of the OXPHOS, in which COX1 is one of the two catalytic subunits, has the property of intrinsic uncoupling , through which it regulates and increases the coupling efficiency to produce ATP or heat. This mechanism is consistent with the negative selection observed for this gene, as adaptations can occur by changing the configuration of the protein, a mechanism that does not occur in complexes I and V.
Further studies are required to completely understand the factors constraining the evolutionary rate of mitogenomes in penguins. Evolution of protein sequences depends on functional constraints including differential gene expression, protein folding and connectivity, functional importance, and expression breadth among tissues . Here, we discussed the implications of functionality of different mtDNA genes on the evolutionary rate of a protein. Alternatively, Zang and Yang , suggest that this rate may be predominantly influenced by its expression level rather than functional importance.
Mitogenomics and ecological-evolutionary hypotheses to explain selection
There are two hypotheses generally offered to explain patterns of mitochondrial genome evolution. The first is that temperature is a primary driver of selection and the evolution of mtDNA genes. For penguins, speciation and adaptation are seemingly strongly correlated with temperature differences of ocean water between the Antarctic and the tropics [35–38]. The mitochondrial OXPHOS generates both energy and heat, with the balance of each depending on the coupling efficiency of the electron transport chain to the ATP synthase activity [2, 3]. Less efficient energy production and more heat production would be beneficial for survival in cold environmental conditions . A reduced proton leak would lead to incremental cold sensitivity that effectively constrains the biogeographic distribution of species with this pattern of mtDNA variation [4, 5]. Selection on mtDNA variation influencing mitochondrial bioenergetics would thus be plausible drivers of adaptation and speciation , and indeed we detected correlation between Ka/Ks from ND4 and sea suface water temperature which should be further investigated.
A second hypothesis is that the risk of starvation associated with restricted food availability (less abundant or less reliable food source) drives mtDNA gene evolution, since starvation can generate more efficient coupling of energy production by the OXPHOS pathway . Three mechanisms have been suggested for the evolution of starvation resistance and each may be influenced by mitochondrial diversity: (a) greater energy reserves (lipid stores); (b) reduced rate at which reserves are utilized; (c) and lowering of the minimal resources required for survival [7, 9]. Reducing energy requirements may occur by reducing ROS and heat production and thus increasing OXPHOS efficiency [8, 10]. Chl can be considered to be a measure or correlate of primary productivity and therefore resource availability in the Antarctic marine food web. Penguins as top predators, are dependent on the productivity of lower trophic levels. Therefore, the high variation in productivity of the habitats throughout the Spheniscus species distributional ranges may have led to the evolution of more efficient coupling in those species, compared to the stable environments of Emperor and Adélie penguins. For example, high variability in levels of Chl is observed within the range of S. humboldti, and becomes further exacerbated during El Niño events, leading to a higher probability of mortality due to starvation. In this circumstance, selection may favor OXPHOS efficiency rather than heat production. During El Niño events, the upwelling of nutrients is interrupted, reducing food availability and increasing water temperature in Chilean-Peruvian waters . This affects fish abundance, causing increased Galapagos and Humboldt penguin mortality, nest desertion and large-scale movements [40, 42, 96]. Although Chl could be an important variable for mtDNA evolution, only a weak correlation between Ka/Ks ratio of COXI and Chl across small geographic scales was detected here. Therefore, our results should be further investigated in penguins not only at a species level, but also in populations of each penguin species across a gradient of sea surface temperatures (e.g. across latitudes) and chlorophyll levels.
Although ATP8 was not correlated to SST and Chl here, the Ka/Ks ratio values were higher than in any other gene compared. Moreover, it showed a remarkable pattern for all comparisons among species distributed in extremely different environments (tropical and hot temperate climate versus polar climate). We found contrasting results for ATP8. A. forsteri and Spheniscus had the lowest values for ATP8 genes in pairwise comparisons among all genera. This may be indicative of different mechanisms of adaptation among genes and species depending on their individual habitat requirements. However, we cannot discount the possibility that this gene is associated with other selective variables, which in turn are correlated with climatic and productive variations observed in the distributional ranges of penguins. Indeed, it was recently suggested that evolution of mitochondrial genes might be linked with pathogen variability .
An increasing number of studies have implicated mitochondria with innate immune responses, including the control of apoptotic signaling, sterile inflammation, and antibacterial and antiviral response [12–15]. Although it is uncertain what direct role mitochondrial encoded genes have in immune response, it bears further evaluation. Antarctica and the tropics differ considerably in virulence, ensemble, abundance, development and diversity of pathogens. Just recently, the influenza virus was detected in penguins in Antarctica [97–99]. Antarctic wildlife is particularly susceptible to infectious agents because of its geographic isolation and extreme environment, which increases the likelihood that the fauna has evolved under low “pathogen pressure” and are thus immunologically naïve .
Implications for future adaptation to changing environments
Ongoing human-induced climatic changes are predicted to have large impacts on the world biota, including penguin species [44–46, 100]. Pygoscelis species appear to be changing their distributions. P. papua populations have been increasing along the Antarctic Peninsula in association with an abrupt increase of water temperatures, while populations of P. adeliae and P. antarcticus have been decreasing [44–46]. The decline and disappearance of an A. forsteri colony was associated with the rise in local mean annual air temperature and decline in seasonal sea ice duration . S. magellanicus populations have experienced increasing reproductive failure, measured in chicks’ mortality associated with starvation, rain, heat, and predation, which has been associated with climate change and an increase in frequency and intensity of storms . S. humboldti and S. mendiculus have been historically affected by strong El Niño events causing high mortality, nesting delay, nest flooding and breeding failure [40–42, 96, 103]. El Niño events are becoming more intense and frequent due to global climate change [104, 105]. The decreasing food availability for penguins in South America during El Niño events , and for penguins in Antarctica during the past 50 years  has caused not only direct mortality. Increased mortality has also been attributed to compromised immune systems, since adequate nutrition is fundamental to maintaining a functioning immune system . The expanding human activities in Antarctica also increase the risk of introduction of infectious diseases into this naïve wildlife community . Therefore, selective pressures on mtDNA genes by environmental factors such as temperature, food availability and novel pathogens could intensify significantly with climate change.
The effect of climate change on gene evolution does not have to be direct. As penguins are top marine predators, changes at all trophic levels will affect penguin population demography . However, understanding the genes involved in adaptation to climate by penguin species may permit us to predict how the species are physiologically limited and how they may be able to respond to environmental changes such as climate warming.
Here, we identify mtDNA candidate genes under selection which could be involved in broad-scale adaptations of penguins to their environment. Our results suggest a strong negative selection for COX1, and a positive selection at codon sites on ATP8 and ND4. We also found a high correlation between ND4 and sea surface temperature. Further studies on these genes at the intraspecific level are required to completely understand penguin adaptation across different environmental conditions. This is a novel study that illustrates the applicability of large-scale remote sensing as a useful approach to comprehend adaptation to the environment occurring at a molecular level. Integration of these different data provide insights into how penguins sepcies have adapted to their enviornments and therefore how they may respond to future, human-induced changes to their environment.
The research was funded by the Sea World and Busch Garden Conservation Fund, Fondecyt N° 11110060, 1150517 INACH N° G06–11, T12–13 and 12–14, UNAB DI-410-13/I, and CNPq 482501/2013–8 and 490403/2008–5, FAPESP 2009/08624. Sampling of Adélie penguins at Cape Crozier was supported through funding from the National Science Foundation to David G. Ainley (ANT 0944411) and logistics supplied by the U.S. Antarctic Program. MM was funded by a New Zealand’s Ministry of Science and Innovation grant (C01X1001). PP was funded by the Saint Louis Zoo and the Des Lee Collaborative Vision.
Availability of data and materials
All genome and sequences were submitted to GenBank (accession numbers KM891593; KU361803-KU361806; KU356673-KU356677).
BR and JAV conceived the entire study. DEL implemented the genomics bioinformatics. JAV, PGP, MM, DGA, GPMD contributed to the sample collection. BR, WEJ, GPMD, JAV analyzed and interpreted the genetic data. MDM conceived the acquisition of data, analysis and interpretation of spatial environmental data from remote sensing and GIS databases. All authors were involved in drafting or revising the content of the manuscript. All authors read and approved the final manuscript.
Ethics approval and consent to participate
All samples were collected with permits in accordance to Annex II, Article 3 of the Protocol on Environmental Protection to the Antarctic Treaty, and the regulation from the Scientific Committee on Antarctic Research (SCAR) provided by the Chilean Antarctic Authorities (permits INACH 44/2012, 45/2016, 46/2016) and the Brazilian Antarctic Authorities through the PROANTAR and Environmental Ministry (N°21/2011 and N 20/2012), and other permits from Subsecretaria de Pesca (Chile, 04/2008, 110/2012, 2086/2014) from other locations. The permits incorporated authorization for sample collection for all locations including special authorization to the studied protected areas. A bioethics permit was provided by Universidad de Concepción and Pontificia Universidad Católica de Chile, which was required to obtain INACH permits and Fondecyt project.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Boore JL. Animal mitochondrial genomes. Nucleic Acids Res. 1999;27:1767–80.View ArticlePubMedPubMed CentralGoogle Scholar
- Pörtner HO. Climate variability and the energetic pathways of evolution: the origin of Endothermy in mammals and birds. Physiol Biochem Zool. 2004;77:959–81.View ArticlePubMedGoogle Scholar
- Wallace DCA. Mitochondrial paradigm of metabolic and degenerative diseases, aging, and cancer: a Dawn for evolutionary medicine. Annu Rev Genet. 2005;39:359–407.View ArticlePubMedPubMed CentralGoogle Scholar
- Ballard JWO, James AC. Differential fitness of mitochondrial DNA in perturbation cage studies correlates with global abundance and population history in Drosophila Simulans. Proc R Soc B Biol Sci. 2004;271:1197–201.View ArticleGoogle Scholar
- Katewa SD, Ballard JWO. Sympatric Drosophila Simulans flies with distinct mtDNA show difference in mitochondrial respiration and electron transport. Insect Biochem Mol Biol. 2007;37:213–22.View ArticlePubMedGoogle Scholar
- Gershoni M, Templeton AR, Mishmar D. Mitochondrial bioenergetics as a major motive force of speciation. BioEssays. 2009;31:642–50.View ArticlePubMedGoogle Scholar
- Ballard JWO, Melvin RG. Linking the mitochondrial genotype to the organismal phenotype. Mol Ecol. 2010;19:1523–39.View ArticlePubMedGoogle Scholar
- Mishmar D, Ruiz-Pesini E, Golik P, Macaulay V, Clark AG, Hosseini S, et al. Natural selection shaped regional mtDNA variation in humans. Proc Natl Acad Sci. 2003;100:171–6.View ArticlePubMedGoogle Scholar
- Rion S, Kawecki TJ. Evolutionary biology of starvation resistance: what we have learned from Drosophila: starvation resistance in Drosophila. J Evol Biol. 2007;20:1655–64.View ArticlePubMedGoogle Scholar
- Ruiz-Pesini E, Mishmar D, Brandon M, Procaccio VV, Wallace DC. Effects of purifying and adaptive selection on regional variation in human mtDNA. Science. 2004;303:223–6.View ArticlePubMedGoogle Scholar
- Dobson A, Foufopoulos J. Emerging infectious pathogens of wildlife. Philos Trans R Soc B Biol Sci. 2001;356:1001–12.View ArticleGoogle Scholar
- Arnoult D, Carneiro L, Tattoli I, Girardin SE. The role of mitochondria in cellular defense against microbial infection. Semin Immunol. 2009;21:223–32.View ArticlePubMedGoogle Scholar
- Ohta A, Nishiyama Y. Mitochondria and viruses. Mitochondrion. 2011;11:1–12.View ArticlePubMedGoogle Scholar
- West AP, Shadel GS, Ghosh S. Mitochondria in innate immune responses. Nat Rev Immunol. 2011;11:389–402.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhou R, Yazdi AS, Menu P, Tschopp JA. Role for mitochondria in NLRP3 inflammasome activation. Nature. 2011;469:221–5.View ArticlePubMedGoogle Scholar
- Friedlaender AS, Johnston DW, Fraser WR, Burns J, Patrick NH, Costa DP. Ecological niche modeling of sympatric krill predators around Marguerite Bay, western Antarctic peninsula. Deep Sea Res Part II Top Stud Oceanogr. 2011;58:1729–40.View ArticleGoogle Scholar
- Bazin E. Population size does not influence mitochondrial genetic diversity in animals. Science. 2006;312:570–2.View ArticlePubMedGoogle Scholar
- Balloux F, Handley L-JL, Jombart T, Liu H, Manica A. Climate shaped the worldwide distribution of human mitochondrial DNA sequence variation. Proc R Soc B Biol Sci. 2009;276:3447–55.View ArticleGoogle Scholar
- Harrisson K, Pavlova A, Gan HM, Lee YP, Austin CM, Sunnucks P. Pleistocene divergence across a mountain range and the influence of selection on mitogenome evolution in threatened Australian freshwater cod species. Heredity. 2016;116:506–15.View ArticlePubMedPubMed CentralGoogle Scholar
- Morales HE, Pavlova A, Joseph L, Sunnucks P. Positive and purifying selection in mitochondrial genomes of a bird with mitonuclear discordance. Mol Ecol. 2015;24:2820–37.View ArticlePubMedGoogle Scholar
- Ribeiro ÂM, Lloyd P, Bowie RCK. A tight balance between natural selection and gene Flow in a southern African arid-zone endemic Bird: spatial heterogeneous environments and gene Flow. Evolution. 2011;65:3499–514.View ArticlePubMedGoogle Scholar
- Hancock AM, Clark VJ, Qian Y, Di Rienzo A. Population genetic analysis of the uncoupling proteins supports a role for UCP3 in human cold resistance. Mol Biol Evol. 2011;28:601–14.View ArticlePubMedGoogle Scholar
- Nunes VL, Beaumont MA, Butlin RK, Paulo OS. Multiple approaches to detect outliers in a genome scan for selection in ocellated lizards (Lacerta Lepida) along an environmental gradient: SELECTION IN OCELLATED LIZARDS. Mol Ecol. 2011;20:193–205.View ArticlePubMedGoogle Scholar
- Li Y, Ren Z, Shedlock AM, Wu J, Sang L, Tersing T, et al. High altitude adaptation of the schizothoracine fishes (Cyprinidae) revealed by the mitochondrial genome analyses. Gene. 2013;517:169–78.View ArticlePubMedGoogle Scholar
- Cheviron ZA, Brumfield RT. Migration-Selection Balance And Local Adaptation Of Mitochondrial Haplotypes In Rufous-Collared Sparrows ( Zonotrichia Capensis ) Along An Elevational Gradient. Evolution. 2009;63:1593–605.View ArticlePubMedGoogle Scholar
- Xu S, Luosang J, Hua S, He J, Ciren A, Wang W, et al. High altitude adaptation and phylogenetic analysis of Tibetan horse based on the mitochondrial genome. J Genet Genomics. 2007;34:720–9.View ArticlePubMedGoogle Scholar
- Zhuang X, CHC C. ND6 Gene “Lost” and Found: Evolution of Mitochondrial Gene Rearrangement in Antarctic Notothenioids. Mol Biol Evol. 2010;27:1391–403.View ArticlePubMedGoogle Scholar
- Menezes AN, Viana MC, Furtado C, Schrago CG, Seuánez HN. Positive selection along the evolution of primate mitogenomes. Mitochondrion. 2013;13:846–51.View ArticlePubMedGoogle Scholar
- Tomasco IH, Lessa EP. The evolution of mitochondrial genomes in subterranean caviomorph rodents: adaptation against a background of purifying selection. Mol Phylogenet Evol. 2011;61:64–70.View ArticlePubMedGoogle Scholar
- Foote AD, Morin PA, Durban JW, Pitman RL, Wade P, Willerslev E, et al. Positive selection on the killer whale mitogenome. Biol Lett. 2011;7:116–8.View ArticlePubMedGoogle Scholar
- Welch AJ, Bedoya-Reina OC, Carretero-Paulet L, Miller W, Rode KD, Lindqvist C. Polar bears exhibit genome-wide signatures of bioenergetic adaptation to life in the Arctic environment. Genome Biol Evol. 2014;6:433–50.View ArticlePubMedPubMed CentralGoogle Scholar
- Kerr KCR. Searching for evidence of selection in avian DNA barcodes: SEARCHING FOR SELECTION IN AVIAN DNA BARCODES. Mol Ecol Resour. 2011;11:1045–55.View ArticlePubMedGoogle Scholar
- Marshall HD, Baker AJ, Grant AR. Complete mitochondrial genomes from four subspecies of common chaffinch (Fringilla Coelebs): new inferences about mitochondrial rate heterogeneity, neutral theory, and phylogenetic relationships within the order Passeriformes. Gene. 2013;517:37–45.View ArticlePubMedGoogle Scholar
- Shen Y-Y, Shi P, Sun Y-B, Zhang Y-P. Relaxation of selective constraints on avian mitochondrial DNA following the degeneration of flight ability. Genome Res. 2009;19:1760–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Baker AJ, Pereira SL, Haddrath OP, Edge K-A. Multiple gene evidence for expansion of extant penguins out of Antarctica due to global cooling. Proc R Soc B Biol Sci. 2006;273:11–7.View ArticleGoogle Scholar
- Clarke JA, Ksepka DT, Stucchi M, Urbina M, Giannini N, Bertelli S, et al. Paleogene equatorial penguins challenge the proposed relationship between biogeography, diversity, and Cenozoic climate change. Proc Natl Acad Sci. 2007;104:11545–50.View ArticlePubMedPubMed CentralGoogle Scholar
- Ksepka DT, Bertelli S, Giannini NP. The phylogeny of the living and fossil Sphenisciformes (penguins). Cladistics. 2006;22:412–41.View ArticleGoogle Scholar
- Subramanian S, Beans-Picon G, Swaminathan SK, Millar CD, Lambert DM. Evidence for a recent origin of penguins. Biol Lett. 2013;9:20130748–8.Google Scholar
- Hays C. Effects of the el Niño 1982–83 on Humboldt penguin colonies in Perú. Biol Conserv. 1986;36:169–80.View ArticleGoogle Scholar
- Simeone A, Araya B, Bernal M, Diebold E, Grzybowski K, Michaels M, et al. Oceanographic and climatic factors influencing breeding and colony attendance patterns of Humboldt penguins Spheniscus Humboldti in central Chile. Mar Ecol Prog Ser. 2002;227:43–50.View ArticleGoogle Scholar
- Paredes R, Zavalaga CB, Battistini G, Majluf P, McGill P. Status of the Humboldt penguin in Peru, 1999-2000. Waterbirds. 2003;26:129.View ArticleGoogle Scholar
- Vianna JA, Cortes M, Ramos B, Sallaberry-Pincheira D, Gonzalez-Acuna D, Dantas GPM, et al. Changes in abundance and distribution of Humboldt penguin Spheniscus Humboldti. Mar Ornithol. 2014;42:153–9.Google Scholar
- Ainley DG. The Adélie penguin: bellwether of climate change. New York: Columbia University Press; 2002.View ArticleGoogle Scholar
- Forcada J, Trathan PN. Penguin responses to climate change in the Southern Ocean. Glob. Change Biol. 2009;15:1618–30.View ArticleGoogle Scholar
- Forcada J, Trathan PN, Reid K, Murphy EJ, Croxall JP. Contrasting population changes in sympatric penguin species in association with climate warming. Glob. Change Biol. 2006;12:411–23.View ArticleGoogle Scholar
- Lynch H, Fagan W, Naveen R, Trivelpiece S, Trivelpiece W. Differential advancement of breeding phenology in response to climate may alter staggered breeding among sympatric pygoscelid penguins. Mar Ecol Prog Ser. 2012;454:135–45.View ArticleGoogle Scholar
- Slack KE. Early penguin fossils, plus mitochondrial genomes, calibrate avian evolution. Mol Biol Evol. 2006;23:1144–55.View ArticlePubMedGoogle Scholar
- Slack KE, Janke A, Penny D, Arnason U. Two new avian mitochondrial genomes (penguin and goose) and a summary of bird and reptile mitogenomic features. Gene. 2003;302:43–52.View ArticlePubMedGoogle Scholar
- Subramanian S, Denver DR, Millar CD, Heupink T, Aschrafi A, Emslie SD, et al. High mitogenomic evolutionary rates and time dependency. Trends Genet. 2009;25:482–6.View ArticlePubMedGoogle Scholar
- Subramanian S, Lingala SG, Swaminathan S, Huynen L, Lambert D. Second generation DNA sequencing of the mitogenome of the chinstrap penguin and comparative genomics of Antarctic penguins. Mitochondrial DNA. 2014;25:271–2.View ArticlePubMedGoogle Scholar
- Li C, Zhang Y, Li J, Kong L, Hu H, Pan H, et al. Two Antarctic penguin genomes reveal insights into their evolutionary history and molecular changes related to the Antarctic environment. GigaScience [Internet]. 2014 [cited 2017 Mar 11];3. Available from: https://academic.oup.com/gigascience/article-lookup/.doi/10.1186/2047-217X-3-27Google Scholar
- Ritchie PA, Ancient DNA. Enables timing of the Pleistocene origin and Holocene expansion of two Adelie penguin lineages in Antarctica. Mol Biol Evol. 2003;21:240–8.View ArticlePubMedGoogle Scholar
- Aljanabi SM, Martinez I. Universal and rapid salt-extraction of high quality genomic DNA for PCR-based techniques. Nucleic Acids Res. 1997;25:4692–3.View ArticlePubMedPubMed CentralGoogle Scholar
- Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The sequence alignment/map format and SAMtools. Bioinformatics. 2009;25:2078–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Zerbino DR, Birney E. Velvet: algorithms for de novo short read assembly using de Bruijn graphs. Genome Res. 2008;18:821–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Labuschagne C, Kotzé A, Grobler JP, Dalton DL. The complete sequence of the mitochondrial genome of the African penguin (Spheniscus Demersus). Gene. 2014;534:113–8.View ArticlePubMedGoogle Scholar
- Watanabe M, Nikaido M, Tsuda TT, Kobayashi T, Mindell D, Cao Y, et al. New candidate species most closely related to penguins. Gene. 2006;378:65–73.View ArticlePubMedGoogle Scholar
- Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG. The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997;25:4876–82.View ArticlePubMedPubMed CentralGoogle Scholar
- Librado P, Rozas J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009;25:1451–2.View ArticlePubMedGoogle Scholar
- Mindell DP, Sorenson MD, Dimcheff DE. An extra nucleotide is not translated in mitochondrial ND3 of some birds and turtles. Mol Biol Evol. 1998;15:1568–71.View ArticlePubMedGoogle Scholar
- Tamura K, Nei M, Kumar S. Prospects for inferring very large phylogenies by using the neighbor-joining method. Proc Natl Acad Sci. 2004;101:11030–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: Molecular Evolutionary Genetics Analysis Version 6.0. Mol Biol Evol. 2013;30:2725–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Excoffier L, Lischer HEL. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and windows. Mol Ecol Resour. 2010;10:564–7.View ArticlePubMedGoogle Scholar
- R Core Team. R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2013. Available from: http://www.R-project.org.
- Gibson AA. Comprehensive analysis of mammalian mitochondrial Genome Base composition and improved phylogenetic methods. Mol Biol Evol. 2004;22:251–64.View ArticlePubMedGoogle Scholar
- Huelsenbeck JP, Ronquist F. MRBAYES: Bayesian inference of phylogenetic trees. Bioinforma Oxf Engl. 2001;17:754–5.View ArticleGoogle Scholar
- Ronquist F, Huelsenbeck JP. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics. 2003;19:1572–4.View ArticlePubMedGoogle Scholar
- Guindon S, Gascuel OA. Simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003;52:696–704.View ArticlePubMedGoogle Scholar
- Rambaut A. FIgTree 1.3.1. 2009. http://tree.bio.ed.ac.uk/software/figtree.
- Meiklejohn CD, Montooth KL, Rand DM. Positive and negative selection on the mitochondrial genome. Trends Genet. 2007;23:259–63.View ArticlePubMedGoogle Scholar
- Kosakovsky Pond SL, Frost SDW. Not so different after all: a comparison of methods for detecting amino acid sites under selection. Mol Biol Evol. 2005;22:1208–22.View ArticlePubMedGoogle Scholar
- Kosakovsky Pond SL, Murrell B, Fourment M, Frost SDW, Delport W, Scheffler KA. Random effects branch-site model for detecting episodic diversifying selection. Mol Biol Evol. 2011;28:3033–43.View ArticlePubMedPubMed CentralGoogle Scholar
- Kosakovsky Pond SLK, Frost SDW, Muse SV. HyPhy: hypothesis testing using phylogenies. Bioinformatics. 2005;21:676–9.View ArticleGoogle Scholar
- Feldman GC, McClain CR. Ocean Color Web, SeaWiFS Reprocessing # 4 NASA Goddard Space Flight Center. Kuring N, Bailey SW, Eds. 2014. http://oceancolor.gsfc.nasa.gov. Accessed 27 Nov 2014.
- Siegel DA, Behrenfeld MJ, Maritorena S, McClain CR, Antoine D, Bailey SW, et al. Regional to global assessments of phytoplankton dynamics from the SeaWiFS mission. Remote Sens Environ. 2013;135:77–91.View ArticleGoogle Scholar
- OceanColor Biology Processing Group. OceanColor [Internet]. Available from: http://oceancolor.gsfc.nasa.gov/REPROCESSING/R2009/ocv6/
- Parkinson CL, Greenstone R. EOS Data Products Handbook. V. 2. NASA Goddard Space Flight Center; Greenbelt, MD United States. 266p. 2000.Google Scholar
- Williams GN, Dogliotti AI, Zaidman P, Solis M, Narvarte MA, González RC, et al. Assessment of remotely-sensed sea-surface temperature and chlorophyll-a concentration in san Matías gulf (Patagonia, Argentina). Cont Shelf Res. 2013;52:159–71.View ArticleGoogle Scholar
- Borboroglu PG, Boersma PD, editors. Penguins: natural history and conservation. Seattle: University of Washington Press; 2013.Google Scholar
- QGIS Development Team. QGIS Geographic Information System. Open Source Geospatial Foundation [Internet]. 2009. Available from: http://qgis.osgeo.org
- Diniz-Filho JAF, Soares TN, Lima JS, Dobrovolski R, Landeiro VL, Telles MP d C, et al. Mantel test in population genetics. Genet Mol Biol. 2013;36:475–85.View ArticlePubMedPubMed CentralGoogle Scholar
- Storfer A, Murphy MA, Evans JS, Goldberg CS, Robinson S, Spear SF, et al. Putting the “landscape” in landscape genetics. Heredity. 2007;98:128–42.View ArticlePubMedGoogle Scholar
- Pereira SL. Mitochondrial genome organization and vertebrate phylogenetics. Genet Mol Biol. 2000;23:745–52.View ArticleGoogle Scholar
- Suzuki Y, Nei M. False-positive selection identified by ML-based methods: examples from the Sig1 gene of the diatom Thalassiosira Weissflogii and the tax gene of a human T-cell Lymphotropic virus. Mol Biol Evol. 2004;21:914–21.View ArticlePubMedGoogle Scholar
- Burger G, Lang BF, Reith M, Gray MW. Genes encoding the same three subunits of respiratory complex II are present in the mitochondrial DNA of two phylogenetically distant eukaryotes. Proc Natl Acad Sci. 1996;93:2328–32.View ArticlePubMedPubMed CentralGoogle Scholar
- Brandt U, Energy Converting NADH. Quinone Oxidoreductase (Complex I). Annu Rev Biochem. 2006;75:69–92.View ArticlePubMedGoogle Scholar
- Sazanov LA, Walker JE. Cryo-electron crystallography of two sub-complexes of bovine complex I reveals the relationship between the membrane and peripheral arms. J Mol Biol. 2000;302:455–64.View ArticlePubMedGoogle Scholar
- da Fonseca RR, Johnson WE, O’Brien SJ, Ramos M, Antunes A. The adaptive evolution of the mammalian mitochondrial genome. BMC Genomics. 2008;9:119.View ArticlePubMedPubMed CentralGoogle Scholar
- Garvin MR, Bielawski JP, Sazanov LA, Gharrett AJ. Review and meta-analysis of natural selection in mitochondrial complex I in metazoans. J Zool Syst Evol Res. 2015;53:1–17.View ArticleGoogle Scholar
- Zink RM. Natural selection on mitochondrial DNA in Parusand its relevance for phylogeographic studies. Proc R Soc B Biol Sci. 2005;272:71–8.View ArticleGoogle Scholar
- Yu L, Wang X, Ting N, Zhang Y. Mitogenomic analysis of Chinese snub-nosed monkeys: evidence of positive selection in NADH dehydrogenase genes in high-altitude adaptation. Mitochondrion. 2011;11:497–503.View ArticlePubMedGoogle Scholar
- Teacher AG, André C, Merilä J, Wheat CW. Whole mitochondrial genome scan for population structure and selection in the Atlantic herring. BMC Evol Biol. 2012;12:248.View ArticlePubMedPubMed CentralGoogle Scholar
- Kadenbach B. Intrinsic and extrinsic uncoupling of oxidative phosphorylation. Biochim Biophys Acta. 2003;1604:77–94.View ArticlePubMedGoogle Scholar
- Zhang J, Yang J-R. Determinants of the rate of protein sequence evolution. Nat Rev Genet. 2015;16:409–20.View ArticlePubMedPubMed CentralGoogle Scholar
- Thiel M, Macaya E, Acuna E, Arntz W, Bastias H, Brokordt K, et al. The Humboldt Current System of Northern and Central Chile: Oceanographic Processes, Ecological Interactions And Socioeconomic Feedback. In: Gibson R, Atkinson R, Gordon J, editors. Oceanogr. Mar. Biol. [Internet]. CRC Press; 2007 [cited 2017 Mar 11]. p. 195–344. Available from: http://www.crcnetbase.com/doi/abs/10.1201/9781420050943.ch6Google Scholar
- Vargas FH, Harrison S, Rea S, Macdonald DW. Biological effects of el Niño on the Galápagos penguin. Biol Conserv. 2006;127:107–14.View ArticleGoogle Scholar
- Barriga GP, Boric-Bargetto D, San Martin MC, Neira V, van Bakel H, Thompsom M, et al. Avian influenza virus H5 strain with north American and Eurasian lineage genes in an Antarctic penguin. Emerg Infect Dis. 2016;22:2221–3.View ArticlePubMedPubMed CentralGoogle Scholar
- Hurt AC, Vijaykrishna D, Butler J, Baas C, Maurer-Stroh S, Silva-de-la-Fuente MC, et al. Detection of Evolutionarily Distinct Avian Influenza A Viruses in Antarctica. mBio. 2014;5 e01098–14-e01098–14Google Scholar
- Hurt AC, Su YCF, Aban M, Peck H, Lau H, Baas C, et al. Evidence for the introduction, Reassortment, and persistence of diverse influenza a viruses in Antarctica. Schultz-cherry S, editor. J Virol 2016;90:9674–9682.Google Scholar
- Ainley D, Russell J, Jenouvrier S, Woehler E, Lyver PO, Fraser WR, et al. Antarctic penguin response to habitat change as Earth’s troposphere reaches 2°C above preindustrial levels. Ecol Monogr. 2010;80:49–66.View ArticleGoogle Scholar
- Trathan PN, Fretwell PT, Stonehouse B. First recorded loss of an emperor penguin Colony in the recent period of Antarctic regional warming: implications for other colonies. Briffa M, editor. PLoS One 2011;6:e14738.Google Scholar
- Boersma PD, Rebstock GA. Climate change increases reproductive failure in Magellanic penguins. Chiaradia a. PLoS One. 2014;9:e85602.View ArticlePubMedPubMed CentralGoogle Scholar
- Boersma PD. Population trends of the Galápagos penguin: impacts of el Niño and la Niña. condor. 1998;100:245–53.View ArticleGoogle Scholar
- Cobb KM, Charles CD, Cheng H, Edwards RL. El Niño/southern oscillation and tropical Pacific climate during the last millennium. Nature. 2003;424:271–6.View ArticlePubMedGoogle Scholar
- Latif M, Kleeman R, Eckert C. Greenhouse warming, decadal variability, or el Niño? An attempt to understand the anomalous 1990s. J Clim. 1997;10:2221–39.View ArticleGoogle Scholar
- Cimino MA, Fraser WR, Irwin AJ, Oliver MJ. Satellite data identify decadal trends in the quality of Pygoscelis penguin chick-rearing habitat. Glob Change Biol. 2013;19:136–48.View ArticleGoogle Scholar
- Kogut MH, Klasing K. An immunologist’s perspective on nutrition, immunity, and infectious diseases: introduction and overview. J Appl Poult Res. 2009;18:103–10.View ArticleGoogle Scholar
- Grimaldi WW, Seddon PJ, Lyver PO, Nakagawa S, Tompkins DM. Infectious diseases of Antarctic penguins: current status and future threats. Polar Biol. 2015;38:591–606.View ArticleGoogle Scholar
- McClintock J, Ducklow H, Fraser W. Ecological responses to climate change on the Antarctic peninsula. Am Sci. 2008;96:302.View ArticleGoogle Scholar
- Williams TD. The penguins: Spheniscidae. Oxford. New York: Oxford University Press; 1995.Google Scholar