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

Cryptic diversity and rampant hybridization in annual gentians on the Qinghai-Tibet Plateau revealed by population genomic analysis

Understanding the evolutionary and ecological processes involved in population differentiation and speciation provides critical insights into biodiversity formation. In this study, we employed 29,865 single nucleotide polymorphisms (SNPs) and complete plastomes to examine genomic divergence and hybridization in Gentiana aristata, which is endemic to the Qinghai-Tibet Plateau (QTP) region. Genetic clustering revealed that G. aristata is characterized by geographic genetic structures with five clusters (West, East, Central, South and North). The West cluster has a specific morphological character (i.e., blue corolla) and higher values of FST compared to the remaining clusters, likely the result of the geological barrier formed by the Yangtze River. The West cluster diverged from the other clusters in the Early Pliocene; these remaining clusters diverged from one another in the Early Quaternary. Phylogenetic reconstructions based on SNPs and plastid data revealed substantial cyto-nuclear conflicts. Genetic clustering and D-statistics demonstrated rampant hybridization between the Central and North clusters, along the Bayankala Mountains, which form the geological barrier between the Central and North clusters. Species distribution modeling demonstrated the range of G. aristata expanded since the Last Interglacial period. Our findings provide genetic and morphological evidence of cryptic diversity in G. aristata, and identified rampant hybridization between genetic clusters along a geological barrier. These findings suggest that geological barriers and climatic fluctuations have an important role in triggering diversification as well as hybridization, indicating that cryptic diversity and hybridization are essential factors in biodiversity formation within the QTP region.

Molecular and acoustic evidence support the species status of Anthus rubescens rubescens and Anthus [rubescens] japonicus (Passeriformes: Motacillidae)

The Buff-bellied Pipit Anthus rubescens comprises two allopatric subspecies groups: A. r. rubescens and A. r. alticola in North America and A. [r.] japonicus in north-east Asia. Despite their great morphological resemblance in breeding plumage, most individuals can be assigned to one or the other subspecies group in non-breeding plumage. Allopatric distributions, morphological differentiation and previously reported molecular divergence suggested the need for additional taxonomic study to assess the rank of these two populations. To resolve the taxonomy of the Buff-bellied Pipit species complex we analysed i) two mitochondrial DNA (mtDNA) loci and ii) nine bioacoustic parameters across 69 sound recordings (338 flight calls) recovered from public databases using principal component analysis and Euclidean distance measures. By comparing our mtDNA and call divergence measures with similar values measured between long-recognised species pairs of the genus, we show that the level of mitochondrial and acoustic divergence between the two Buff-bellied Pipit subspecies groups is typical of species-level divergence in the genus Anthus. Therefore, we recommend splitting the Buff-bellied Pipit species complex into two species: Anthus rubescens (American Pipit) and Anthus japonicus (Siberian Pipit). Our results also suggest that the Water Pipit A. spinoletta deserves taxonomic reassessment as its lineages are highly divergent in acoustics and mtDNA, while mtDNA relationships suggest paraphyly relative to the Rock Pipit A. petrosus. Our work highlights the crucial importance of integrative approaches in taxonomy and the usefulness of bioacoustics in studying cryptic diversity.

Genomic Analyses Implicate the Amazon-Orinoco Plume as the Driver of Cryptic Speciation in a Swimming Crab

The Amazon-Orinoco plume (AOP) is the world's largest freshwater and sediment discharge into the ocean. Previous studies limited to mtDNA suggest that the swimming crab Callinectes ornatus Ordway, 1863 exists as two distinct genetic clusters separated by the AOP. However, questions concerning migration, diversification time, and species delimitation are unresolved. Densely sampling markers across the genome (SNPs) could elucidate the evolutionary processes within this species. Here, we combined mtDNA data and ddRAD-seq to explore the diversification patterns and processes within the swimming crab C. ornatus. We show great genetic differentiation between groups on the north and south sides of the plume but also signs of hybridization. Demographic modeling indicates the divergence between groups starting around 8 Mya following the AOP's formation. After a period of isolation, we detect two incidences of secondary contact with stronger migration in concordance with the North Brazil Current flow. Our results suggest speciation with gene flow explained by the interplay among the AOP, oceanographic currents, and long larval dispersal. This work represents the first investigation employing ddRAD-seq in a marine invertebrate species with distribution encompassing the north and south Atlantic and sheds light on the role of the AOP in the diversification of a marine species.

The right tool for the right question: contrasting biogeographic patterns in the notothenioid fish Harpagifer spp. along the Magellan Province

Molecular-based analysis has become a fundamental tool to understand the role of Quaternary glacial episodes. In the Magellan Province in southern South America, ice covering during the last glacial maximum (20 ka) radically altered the landscape/seascape, speciation rates and distribution of species. For the notothenioid fishes of the genus Harpagifer, in the area are described two nominal species. Nevertheless, this genus recently colonized South America from Antarctica, providing a short time for speciation processes. Combining DNA sequences and genotyping-by-sequencing SNPs, we evaluated the role of Quaternary glaciations over the patterns of genetic structure in Harpagifer across its distribution in the Magellan Province. DNA sequences showed low phylogeographic structure, with shared and dominant haplotypes between nominal species, suggesting a single evolutionary unit. SNPs identified contrastingly two groups in Patagonia and a third well-differentiated group in the Falkland/Malvinas Islands with limited and asymmetric gene flow. Linking the information of different markers allowed us to infer the relevance of postglacial colonization mediated by the general oceanographic circulation patterns. Contrasting rough- and fine-scale genetic patterns highlights the relevance of combined methodologies for species delimitation, which, depending on the question to be addressed, allows discrimination among phylogeographic structure, discarding incipient speciation, and contemporary spatial differentiation processes.

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.

Strong Philopatry, Isolation by Distance, and Local Habitat Have Promoted Genetic Structure in Heermann’s Gull

Philopatry can promote genetic differentiation among populations but remains undescribed in many seabirds. Hence, we explored such associations in Heermann’s Gull. Philopatry was observed monitoring 998 gulls in Rasa Island, while genetic differences were examined in the Cardonosa, Rasa, and Isabel islands using the cytochrome b of 296 gulls. Adults returned repeatedly to its natal valley or to a very close distance from it under different modelled hypotheses. Likewise, the interaction between sex and distance indicated significant male-biased philopatry. Besides, low to high genetic differentiation was observed between the Rasa and Cardonosa islands (ΦST = 0–0.22) (both in the Midriff Islands Region), but higher genetic differentiation against Isabel Island (ΦST > 0.25) (in the Mexican Province region). Consistently, genetic structure among regions was observed using different approaches (AMOVA: ΦCT = 0.49; SAMOVA: FCT = 0.49; and BAPS: K = 2). Similarly, a pattern of isolation by distance (rM = 0.82, p = 0.03), agrees with lower estimates of scaled migration rates between regions than among islands of the same region. Overall, it is suggested that the genetic structure found in Heermann’s Gull has been promoted by physical and behavioral barriers.

Morphometric and genetic evidence for four species of gentoo penguin

Gentoo penguins (Pygoscelis papua) are found across the Southern Ocean with a circumpolar distribution and notable genetic and morphological variation across their geographic range. Whether this geographic variation represents species-level diversity has yet to be investigated in an integrative taxonomic framework. Here, we show that four distinct populations of gentoo penguins (Iles Kerguelen, Falkland Islands, South Georgia, and South Shetlands/Western Antarctic Peninsula) are genetically and morphologically distinct from one another. We present here a revised taxonomic treatment including formal nomenclatural changes. We suggest the designation of four species of gentoo penguin: P. papua in the Falkland Islands, P. ellsworthi in the South Shetland Islands/Western Antarctic Peninsula, P. taeniata in Iles Kerguelen, and a new gentoo species P. poncetii, described herein, in South Georgia. These findings of cryptic diversity add to many other such findings across the avian tree of life in recent years. Our results further highlight the importance of reassessing species boundaries as methodological advances are made, particularly for taxa of conservation concern. We recommend reassessment by the IUCN of each species, particularly P. taeniata and P. poncetii, which both show evidence of decline.

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 (N _e ) 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 N _e 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.