Cryptic speciation in gentoo penguins is driven by geographic isolation and regional marine conditions: Unforeseen vulnerabilities to global change

[Aim]: The conservation of biodiversity is hampered by data deficiencies, with many new species and subspecies awaiting description or reclassification. Population genomics and ecological niche modelling offer complementary new tools for uncovering functional units of phylogenetic diversity. We hypothesize that phylogenetically delineated lineages of gentoo penguins (Pygoscelis papua) distributed across Antarctica and sub-Antarctic Islands are subject to spatially explicit ecological conditions that have limited gene flow, facilitating genetic differentiation, and thereby speciation processes. [Location]: Antarctica and sub-Antarctic area. [Methods]: We identify divergent lineages for gentoo penguins using ddRAD-seq and mtDNA, and generated species distribution models (SDMs) based on terrestrial and marine parameters. Results: Analyses of our genomic data supports the existence of four major lineages of gentoo penguin: (i) spanning the sub-Antarctic archipelagos north of the Antarctic Polar Front (APF); (ii) Kerguelen Island; (iii) South America; and (iv) across maritime Antarctic and the Scotia Arc archipelagos. The APF, a major current system around Antarctica, acts as the most important barrier separating regional sister lineages. Our ecological analyses spanning both the terrestrial (breeding sites) and marine (feeding sites) realms recover limited niche overlap among the major lineages of gentoo penguin. We observe this pattern to correspond more closely with regional differentiation of marine conditions than to terrestrial macroenvironmental features. [Main conclusions]: Recognition of regional genetic lineages as discrete evolutionary entities that occupy distinct ecological niches and also differ morphologically should be considered a priority for conservation. Gentoo penguins provide a good example of how conservation policy can be directly impacted by new insights obtained through the integration of larger genomic datasets with novel approaches to ecological modelling. This is particularly pertinent to polar environments that are among the most rapidly changing environments on earth.Financial support was provided by INACH RT_12–14, Fondecyt Project 1150517, GAB (PIA CONICYT ACT172065), the CNPq (482501/2013-8 and 431463/2016-6) and the PROANTAR, IPEV prog 354 ETHOTAAF and by the Spanish Research Agency (CGL2004-01348, CGL2007-60369, POL2006-06635 and CTM2015-64720-R). LRP was contracted by the NICHEAPPS project (CGL2014-56416-P) and recipient of an Antarctic Science Bursary (2017-18 scheme). PA was supported by “Ramón y Cajal” contract (RYC-2011-07670, MINECO-MNCN)