Identification and function analysis of BCL2 in immune response of Pteria penguin.

B-cell lymphoma/leukemia-2 (BCL2), an anti-apoptotic factor in the mitochondrial regulatory pathway of apoptosis, is critically important in immune defenses. In this study, a novel BCL2 gene was characterized from Pteria penguin (P. penguin). The PpBCL2 was 1482 bp long, containing an open reading frame (ORF) of 588 bp encoding 195 amino acids. Four highly conserved BCL-2 homology (BH) domains were found in PpBCL2. Amino acid alignment and phylogenetic tree showed that PpBCL2 had the highest similarity with BCL2 of Crassostrea gigas at 65.24 %. Tissue expression analysis showed that PpBCL2 had high constitutive expression in gill, digestive diverticulum and mantle, and was significantly increased 72 h of Vibrio parahaemolyticus (V. parahaemolyticus) challenge in these immune tissues. Furthermore, PpBCL2 silencing significantly inhibited antimicrobial activity of hemolymph supernatant by 1.4-fold, and significantly reduced the survival rate by 51.7 % at 72 h post infection in P. penguin. These data indicated that PpBCL2 played an important role in immune response of P. penguin against V. parahaemolyticus infection.

Spread of antibiotic resistance genes to Antarctica by migratory birds.

Recent studies have highlighted the presence of antibiotic resistance genes (ARGs) in Antarctica, which are typically indicative of human activity. However, these studies have concentrated in the Antarctic Peninsula region, and relatively less is known about ARG prevalence in East Antarctica, where human activity levels are lower compared to the Antarctic Peninsula. In addition, the mechanisms of ARG transmission to Antarctica through natural or anthropogenic pathways remain unclear. In this study, we analyzed the fecal samples of Adélie penguins and South polar skuas by using high-throughput sequencing and microfluidic quantitative PCR to detect potential pathogens and ARGs at their breeding colonies near Syowa Station in East Antarctica. These results revealed the presence of several potential pathogens in the fecal matter of both bird species. However, the HF183 marker, which indicates human fecal contamination, was absent in all samples, as well as seawater sampled near the breeding colonies. This suggests that the human fecal contamination was negligible in our study area. In addition to pathogens, we found a significant number of ARGs and metal resistance genes in the feces of both Adélie penguins and South polar skuas, with higher detection rates in skuas than in penguins. To better understand how these birds acquire and transmit these genes, we analyzed the migratory patterns of Adélie penguins and South polar skuas by geolocator tracking. We found that the skuas migrate to the tropical and subtropical regions of the Indian Ocean during the austral winter. On the other hand, Adélie penguins exhibited a more localized migration pattern, mainly staying within Antarctic waters. Because the Indian Ocean is considered one of the major reservoirs of ARGs, South polar skuas might be exposed to ARGs during their winter migration and transfer these genes to Antarctica.

Whole genome sequencing reveals stepping-stone dispersal buffered against founder effects in a range expanding seabird.

Many species are shifting their ranges in response to climate-driven environmental changes, particularly in high-latitude regions. However, the patterns of dispersal and colonization during range shifting events are not always clear. Understanding how populations are connected through space and time can reveal how species navigate a changing environment. Here, we present a fine-scale population genomics study of gentoo penguins (Pygoscelis papua), a presumed site-faithful colonial nesting species that has increased in population size and expanded its range south along the Western Antarctic Peninsula. Using whole genome sequencing, we analysed 129 gentoo penguin individuals across 12 colonies located at or near the southern range edge. Through a detailed examination of fine-scale population structure, admixture, and population divergence, we inferred that gentoo penguins historically dispersed rapidly in a stepping-stone pattern from the South Shetland Islands leading to the colonization of Anvers Island, and then the adjacent mainland Western Antarctica Peninsula. Recent southward expansion along the Western Antarctic Peninsula also followed a stepping-stone dispersal pattern coupled with limited post-divergence gene flow from colonies on Anvers Island. Genetic diversity appeared to be maintained across colonies during the historical dispersal process, and range-edge populations are still growing. This suggests large numbers of migrants may provide a buffer against founder effects at the beginning of colonization events to maintain genetic diversity similar to that of the source populations before migration ceases post-divergence. These results coupled with a continued increase in effective population size since approximately 500-800 years ago distinguish gentoo penguins as a robust species that is highly adaptable and resilient to changing climate.

The PENGUIN approach to reconstruct protein interactions at enhancer-promoter regions and its application to prostate cancer.

We introduce Promoter-Enhancer-Guided Interaction Networks (PENGUIN), a method for studying protein-protein interaction (PPI) networks within enhancer-promoter interactions. PENGUIN integrates H3K27ac-HiChIP data with tissue-specific PPIs to define enhancer-promoter PPI networks (EPINs). We validated PENGUIN using cancer (LNCaP) and benign (LHSAR) prostate cell lines. Our analysis detected EPIN clusters enriched with the architectural protein CTCF, a regulator of enhancer-promoter interactions. CTCF presence was coupled with the prevalence of prostate cancer (PrCa) single nucleotide polymorphisms (SNPs) within the same EPIN clusters, suggesting functional implications in PrCa. Within the EPINs displaying enrichments in both CTCF and PrCa SNPs, we also show enrichment in oncogenes. We substantiated our identified SNPs through CRISPR/Cas9 knockout and RNAi screens experiments. Here we show that PENGUIN provides insights into the intricate interplay between enhancer-promoter interactions and PPI networks, which are crucial for identifying key genes and potential intervention targets. A dedicated server is available at https://penguin.life.bsc.es/ .

Corynebacterium antarcticum sp. nov., Corynebacterium marambiense sp. nov., Corynebacterium meridianum sp. nov., and Corynebacterium pygosceleis sp. nov., isolated from Adélie penguins (Pygoscelis adeliae).

A taxonomic study was conducted on 16 bacterial strains isolated from wild Adélie penguins (Pygoscelis adeliae) from Seymour (Marambio) Island and James Ross Island. An initial screening by repetitive sequence-based PCR fingerprinting divided the strains studied into four coherent groups. Phylogenetic analysis based on 16S rRNA gene sequences assigned all groups to the genus Corynebacterium and showed that Corynebacterium glyciniphilum and Corynebacterium terpenotabidum were the closest species with 16S rRNA gene sequence similarities between 95.4 % and 96.5 %. Further examination of the strains studied with ribotyping, MALDI-TOF mass spectrometry, comprehensive biotyping and calculation of average nucleotide identity and digital DNA-DNA hybridisation values confirmed the separation of the four groups from each other and from the other Corynebacterium species. Chemotaxonomically, the four strains P5828T, P5850T, P6136T, P7210T representing the studied groups were characterised by C16:0 and C18:1ω9c as the major fatty acids, by the presence of meso-diaminopimelic acid in the peptidoglycan, the presence of corynemycolic acids and a quinone system with the predominant menaquinone MK-9(H2). The results of this study show that the strains studied represent four new species of the genus Corynebacterium, for which the names Corynebacterium antarcticum sp. nov. (type strain P5850T = CCM 8835T = LMG 30620T), Corynebacterium marambiense sp. nov. (type strain P5828T = CCM 8864T = LMG 31626T), Corynebacterium meridianum sp. nov. (type strain P6136T = CCM 8863T = LMG 31628T) and Corynebacterium pygosceleis sp. nov. (type strain P7210T = CCM 8836T = LMG 30621T) are proposed.

Selection-driven adaptation to the extreme Antarctic environment in the Emperor penguin.

The eco-evolutionary history of penguins is characterised by shifting from temperate to cold environments. Breeding in Antarctica, the Emperor penguin appears as an extreme outcome of this process, with unique features related to insulation, heat production and energy management. However, whether this species actually diverged from a less cold-adapted ancestor, more ecologically similar to its sister species, the King penguin, is still an open question. As the Antarctic colonisation likely resulted in vast changes in selective pressure experienced by the Emperor penguin, the relative quantification of the genomic signatures of selection, unique to each sister species, could answer this question. Applying phylogeny-based selection tests on 7651 orthologous genes, we identified a more pervasive selection shift in the Emperor penguin than in the King penguin, supporting the hypothesis that its extreme cold adaptation is a derived state. Furthermore, among candidate genes under selection, four (TRPM8, LEPR, CRB1, and SFI1) were identified before in other cold-adapted homeotherms, like the woolly Mammoth, while other 161 genes can be assigned to biological functions relevant to cold adaptation identified in previous studies. Location and structural effects of TRPM8 substitutions in Emperor and King penguin lineages support their functional role with putative divergent effects on thermal adaptation. We conclude that extreme cold adaptation in the Emperor penguin largely involved unique genetic options which, however, affect metabolic and physiological traits common to other cold-adapted homeotherms.

Genomic insights into the secondary aquatic transition of penguins.

Penguins lost the ability to fly more than 60 million years ago, subsequently evolving a hyper-specialized marine body plan. Within the framework of a genome-scale, fossil-inclusive phylogeny, we identify key geological events that shaped penguin diversification and genomic signatures consistent with widespread refugia/recolonization during major climate oscillations. We further identify a suite of genes potentially underpinning adaptations related to thermoregulation, oxygenation, diving, vision, diet, immunity and body size, which might have facilitated their remarkable secondary transition to an aquatic ecology. Our analyses indicate that penguins and their sister group (Procellariiformes) have the lowest evolutionary rates yet detected in birds. Together, these findings help improve our understanding of how penguins have transitioned to the marine environment, successfully colonizing some of the most extreme environments on Earth.

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

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

Penguin: A tool for predicting pseudouridine sites in direct RNA nanopore sequencing data.

Pseudouridine is one of the most abundant RNA modifications, occurring when uridines are catalyzed by Pseudouridine synthase proteins. It plays an important role in many biological processes and has been reported to have application in drug development. Recently, the single-molecule sequencing techniques such as the direct RNA sequencing platform offered by Oxford Nanopore technologies have enabled direct detection of RNA modifications on the molecule being sequenced. In this study, we introduce a tool called Penguin that integrates several machine learning (ML) models to identify RNA Pseudouridine sites on Nanopore direct RNA sequencing reads. Pseudouridine sites were identified on single molecule sequencing data collected from direct RNA sequencing resulting in 723 K reads in Hek293 and 500 K reads in Hela cell lines. Penguin extracts a set of features from the raw signal measured by the Oxford Nanopore and the corresponding basecalled k-mer. Those features are used to train the predictors included in Penguin, which in turn, can predict whether the signal is modified by the presence of Pseudouridine sites in the testing phase. We have included various predictors in Penguin, including Support vector machines (SVM), Random Forest (RF), and Neural network (NN). The results on the two benchmark data sets for Hek293 and Hela cell lines show outstanding performance of Penguin either in random split testing or in independent validation testing. In random split testing, Penguin has been able to identify Pseudouridine sites with a high accuracy of 93.38% by applying SVM to Hek293 benchmark dataset. In independent validation testing, Penguin achieves an accuracy of 92.61% by training SVM with Hek293 benchmark dataset and testing it for identifying Pseudouridine sites on Hela benchmark dataset. Thus, Penguin outperforms the existing Pseudouridine predictors in the literature by 16 % higher accuracy than those predictors using independent validation testing. Employing penguin to predict Pseudouridine sites revealed a significant enrichment of "regulation of mRNA 3'-end processing" in Hek293 cell line and 'positive regulation of transcription from RNA polymerase II promoter involved in cellular response to chemical stimulus' in Hela cell line. Penguin software and models are available on GitHub at https://github.com/Janga-Lab/Penguin and can be readily employed for predicting Ψ sites from Nanopore direct RNA-sequencing datasets.

Adaptation and Cryptic Pseudogenization in Penguin Toll-Like Receptors.

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

Genetic diversity and structure of captive gentoo penguin populations in Japan.

Until the last decade, gentoo penguins were usually split into two subspecies, northern gentoo penguins (Pygoscelis papua papua) breeding in the Falkland Islands, South Georgia, and other subantarctic islands and southern gentoo penguins (P. papua ellsworthi) breeding in the South Sandwich, South Orkney and South Shetland islands, and Antarctic Peninsula. Recent genetics research, however, suggests that the population at South Georgia is much more closely related to those further south and should be included in P. papua ellsworthi. In Japanese zoos and aquariums, captive breeding of gentoo penguins is conducted separately in three populations: "Captive-South Georgia," originating from South Georgia, "Captive-South Shetlands," originating from South Shetlands, and "Captive-Unknown," originating from at least one founder of unknown subspecies. The aims of the present study were to investigate the genetic diversity and differentiation of these captive populations using microsatellite analysis. Genetic diversity in each captive population was similar to that found in the wild, although they had much lower contemporary effective population sizes. Pairwise genetic differentiation indexes (FST ) among the three captive populations were as follows: 0.0309 ("Captive-South Georgia" and "Captive-Unknown"), 0.1094 ("Captive-South Georgia" and "Captive-South Shetlands"), and 0.1214 ("Captive-South Shetlands" and "Captive-Unknown"). Using Bayesian clustering, there was relatively high genetic differentiation between the "Captive-South Shetlands" group, which formed a distinct cluster, and individuals of the "Captive-Unknown" group, which were assigned to clusters in common with "Captive-South Georgia." The results from the present study are useful for future management of captive gentoo penguin populations in Japan.

Molecular epidemiology of aspergillosis in Magellanic penguins and susceptibility patterns of clinical isolates.

Aspergillus section Fumigati is reported in up to 99% of aspergillosis cases in penguins. So far, no data regarding molecular epidemiology and azole resistance are available for A. fumigatus isolates collected from Magellanic penguins. The aim of this work was to perform molecular identification of Aspergillus section Fumigati at species level, to genotype those isolates using microsatellite markers, to evaluate the in vitro susceptibility patterns of A. fumigatus sensu stricto, and to characterize the cyp51A gene in clinical A. fumigatus strains isolated from Magellanic penguins with proven aspergillosis. All 34 isolates included in the study were identified as A. fumigatus sensu stricto. Analyzing the genetic diversity of the isolates of A. fumigatus sensu stricto, we identified two possible outbreaks in the rehabilitation center and we also observed the maintenance of clonal strains through the years. One A. fumigatus sensu stricto isolate was resistant to posaconazole, but the mutations found in the cyp51A gene of this isolate have not been described as conferring phenotypic resistance, suggesting that other mechanisms of resistance could be involved in the resistance of this isolate. With this study, we were able to understand the molecular diversity of Aspergillus fumigatus isolates collected from Magellanic penguins, to characterize them and to associate them with the described global population of Aspergillus fumigatus.

A. fumigatus sensu stricto is of great importance in penguins' aspergillosis. We could identify two outbreaks in the rehabilitation center and the maintenance of clonal strains through the years. Regarding antifungal prophylaxis, it may proceed, but preferably with surveillance for azole resistance.

Genetic evidence of hybridization between Magellanic (Sphensicus magellanicus) and Humboldt (Spheniscus humboldti) penguins in the wild.

The process of hybridization between closely related species plays an important role in defining the genetic integrity and overall genetic diversity of species. The distribution range of Magellanic (Spheniscus magellanicus) and Humboldt (Spheniscus humboldti) penguins is predominantly allopatric; however, the species share a region of sympatry where they may hybridize. We analyzed four types of genetic markers (including nuclear and mitochondrial markers) to assess their utility in detecting hybridization events between Magellanic and Humboldt penguins. Genetic assessment of non-introgressed reference samples allowed us to identify three types of informative markers (microsatellites, major histocompatibility complex, and mitochondrial DNA) and detect positive evidence of introgressive hybridization in the wild. Four out of six putative hybrids showed positive evidence of hybridization, revealed by the detection of Humboldt mitochondrial DNA and Magellanic species-specific alleles from nuclear markers. Bayesian Structure analysis, including samples from the sympatric region of the species in the southern Pacific Ocean, confirmed the use of nuclear markers for detecting hybridization and genetic admixture of putative hybrids, but revealed relatively low levels of genetic introgression at the population level. These findings provide insights into the role of hybridization in regions of species sympatry and its potential consequences on the levels of genetic introgression, genetic diversity, and conservation of these penguin species.

Genome-wide analyses reveal drivers of penguin diversification.

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

Foster rather than biological parental telomere length predicts offspring survival and telomere length in king penguins.

Because telomere length and dynamics relate to individual growth, reproductive investment and survival, telomeres have emerged as possible markers of individual quality. Here, we tested the hypothesis that, in species with parental care, parental telomere length can be a marker of parental quality that predicts offspring phenotype and survival. In king penguins (Aptenodytes patagonicus), we experimentally swapped the single egg of 66 breeding pairs just after egg laying to disentangle the contribution of prelaying parental quality (e.g., genetics, investment in the egg) and/or postlaying parental quality (e.g., incubation, postnatal feeding rate) on offspring growth, telomere length and survival. Parental quality was estimated through the joint effects of biological and foster parent telomere length on offspring traits, both soon after hatching (day 10) and at the end of the prewinter growth period (day 105). We expected that offspring traits would be mostly related to the telomere lengths (i.e., quality) of biological parents at day 10 and to the telomere lengths of foster parents at day 105. Results show that chick survival up to 10 days was negatively related to biological fathers' telomere length, whereas survival up to 105 days was positively related to foster fathers' telomere lengths. Chick growth was not related to either biological or foster parents' telomere length. Chick telomere length was positively related to foster mothers' telomere length at both 10 and 105 days. Overall, our study shows that, in a species with biparental care, parents' telomere length is foremost a proxy of postlaying parental care quality, supporting the "telomere - parental quality hypothesis."

Sustained RNA virome diversity in Antarctic penguins and their ticks.

Despite its isolation and extreme climate, Antarctica is home to diverse fauna and associated microorganisms. It has been proposed that the most iconic Antarctic animal, the penguin, experiences low pathogen pressure, accounting for their disease susceptibility in foreign environments. There is, however, a limited understanding of virome diversity in Antarctic species, the extent of in situ virus evolution, or how it relates to that in other geographic regions. To assess whether penguins have limited microbial diversity we determined the RNA viromes of three species of penguins and their ticks sampled on the Antarctic peninsula. Using total RNA sequencing we identified 107 viral species, comprising likely penguin associated viruses (n = 13), penguin diet and microbiome associated viruses (n = 82), and tick viruses (n = 8), two of which may have the potential to infect penguins. Notably, the level of virome diversity revealed in penguins is comparable to that seen in Australian waterbirds, including many of the same viral families. These data run counter to the idea that penguins are subject to lower pathogen pressure. The repeated detection of specific viruses in Antarctic penguins also suggests that rather than being simply spill-over hosts, these animals may act as key virus reservoirs.

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

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

The role of demographic history and selection in shaping genetic diversity of the Galápagos penguin (Spheniscus mendiculus).

Although many studies have documented the effects of demographic bottlenecks on the genetic diversity of natural populations, there is conflicting evidence of the roles that genetic drift and selection may play in driving changes in genetic variation at adaptive loci. We analyzed genetic variation at microsatellite and mitochondrial loci in conjunction with an adaptive MHC class II locus in the Galápagos penguin (Spheniscus mendiculus), a species that has undergone serial demographic bottlenecks associated with El Niño events through its evolutionary history. We compared levels of variation in the Galápagos penguin to those of its congener, the Magellanic penguin (Spheniscus magellanicus), which has consistently maintained a large population size and thus was used as a non-bottlenecked control. The comparison of neutral and adaptive markers in these two demographically distinct species allowed assessment of the potential role of balancing selection in maintaining levels of MHC variation during bottleneck events. Our analysis suggests that the lack of genetic diversity at both neutral and adaptive loci in the Galápagos penguin likely resulted from its restricted range, relatively low abundance, and history of demographic bottlenecks. The Galápagos penguin revealed two MHC alleles, one mitochondrial haplotype, and six alleles across five microsatellite loci, which represents only a small fraction of the diversity detected in Magellanic penguins. Despite the decreased genetic diversity in the Galápagos penguin, results revealed signals of balancing selection at the MHC, which suggest that selection can mitigate some of the effects of genetic drift during bottleneck events. Although Galápagos penguin populations have persisted for a long time, increased frequency of El Niño events due to global climate change, as well as the low diversity exhibited at immunological loci, may put this species at further risk of extinction.

Phylogeography, Population Structure, and Species Delimitation in Rockhopper Penguins (Eudyptes chrysocome and Eudyptes moseleyi).

Rockhopper penguins are delimited as 2 species, the northern rockhopper (Eudyptes moseleyi) and the southern rockhopper (Eudyptes chrysocome), with the latter comprising 2 subspecies, the western rockhopper (Eudyptes chrysocome chrysocome) and the eastern rockhopper (Eudyptes chrysocome filholi). We conducted a phylogeographic study using multilocus data from 114 individuals sampled across 12 colonies from the entire range of the northern/southern rockhopper complex to assess potential population structure, gene flow, and species limits. Bayesian and likelihood methods with nuclear and mitochondrial DNA, including model testing and heuristic approaches, support E. moseleyi and E. chrysocome as distinct species lineages with a divergence time of 0.97 Ma. However, these analyses also indicated the presence of gene flow between these species. Among southern rockhopper subspecies, we found evidence of significant gene flow and heuristic approaches to species delimitation based on the genealogical diversity index failed to delimit them as species. The best-supported population models for the southern rockhoppers were those where E. c. chrysocome and E. c. filholi were combined into a single lineage or 2 lineages with bidirectional gene flow. Additionally, we found that E. c. filholi has the highest effective population size while E. c. chrysocome showed similar effective population size to that of the endangered E. moseleyi. We suggest that the current taxonomic definitions within rockhopper penguins be upheld and that E. chrysocome populations, all found south of the subtropical front, should be treated as a single taxon with distinct management units for E. c. chrysocome and E. c. filholi.