Genetic structure in the European endemic seabird, Phalacrocorax aristotelis, shaped by a complex interaction of historical and contemporary, physical and nonphysical drivers

Hybridization of Atlantic puffins in the Arctic coincides with 20th-century climate change

The Arctic is experiencing the fastest rates of global warming, leading to shifts in the distribution of its biota and increasing the potential for hybridization. However, genomic evidence of recent hybridization events in the Arctic remains unexpectedly rare. Here, we use whole-genome sequencing of contemporary and 122-year-old historical specimens to investigate the origin of an Arctic hybrid population of Atlantic puffins (Fratercula arctica) on Bjørnøya, Norway. We show that the hybridization between the High Arctic, large-bodied subspecies F. a. naumanni and the temperate, smaller-sized subspecies F. a. arctica began as recently as six generations ago due to an unexpected southward range expansion of F. a. naumanni. Moreover, we find a significant temporal loss of genetic diversity across Arctic and temperate puffin populations. Our observations provide compelling genomic evidence of the impacts of recent distributional shifts and loss of diversity in Arctic communities during the 20th century.

Geography vs. past climate: the drivers of population genetic structure of the Himalayan langur

Background

Contemporary species distribution, genetic diversity and evolutionary history in many taxa are shaped by both historical and current climate as well as topography. The Himalayas show a huge variation in topography and climatic conditions across its entire range, and have experienced major climatic fluctuations in the past. However, very little is known regarding how this heterogenous landscape has moulded the distribution of Himalayan fauna. A recent study examined the effect of these historical events on the genetic diversity of the Himalayan langurs in Nepal Himalaya. However, this study did not include the samples from the Indian Himalayan region (IHR). Therefore, here we revisit the questions addressed in the previous study with a near complete sampling from the IHR, along with the samples from the Nepal Himalaya. We used the mitochondrial Cytochrome-b (Cyt-b, 746 bp) region combined with multiple phylogeographic analyses and palaeodistribution modelling.

Results

Our dataset contained 144 sequences from the IHR as well as the Nepal Himalaya. Phylogenetic analysis showed a low divergent western clade nested within high divergent group of eastern lineages and in the network analysis we identified 22 haplotypes over the entire distribution range of the Himalayan langurs. Samples from the Nepal Himalaya showed geographically structured haplotypes corresponding to different river barriers, whereas samples from IHR showed star-like topology with no structure. Our statistical phylogeography analysis using diyABC supported the model of east to west colonisation of these langurs with founder event during colonisation. Analysis of demographic history showed that the effective population size of the Himalayan langurs decreased at the onset of last glacial maximum (LGM) and started increasing post LGM. The palaeodistribution modelling showed that the extent of suitable habitat shifted from low elevation central Nepal, and adjoining parts of north India, during LGM to the western Himalaya at present.

Conclusion

The current genetic diversity and distribution of Himalayan langurs in the Nepal Himalaya has been shaped by river barriers, whereas the rivers in the IHR had relatively less time to act as a strong genetic barrier after the recent colonisation event. Further, the post LGM expansion could have had confounding effect on Himalayan langur population structure in both Nepal Himalaya and IHR.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02054-1.

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.

Geolocator tagging links distributions in the non-breeding season to population genetic structure in a sentinel North Pacific seabird

We tested the hypothesis that segregation in wintering areas is associated with population differentiation in a sentinel North Pacific seabird, the rhinoceros auklet (Cerorhinca monocerata). We collected tissue samples for genetic analyses on five breeding colonies in the western Pacific Ocean (Japan) and on 13 colonies in the eastern Pacific Ocean (California to Alaska), and deployed light-level geolocator tags on 12 eastern Pacific colonies to delineate wintering areas. Geolocator tags were deployed previously on one colony in Japan. There was strong genetic differentiation between populations in the eastern vs. western Pacific Ocean, likely due to two factors. First, glaciation over the North Pacific in the late Pleistocene might have forced a southward range shift that historically isolated the eastern and western populations. And second, deep-ocean habitat along the northern continental shelf appears to act as a barrier to movement; abundant on both sides of the North Pacific, the rhinoceros auklet is virtually absent as a breeder in the Aleutian Islands and Bering Sea, and no tagged birds crossed the North Pacific in the non-breeding season. While genetic differentiation was strongest between the eastern vs. western Pacific, there was also extensive differentiation within both regional groups. In pairwise comparisons among the eastern Pacific colonies, the standardized measure of genetic differentiation (Fꞌ_ST) was negatively correlated with the extent of spatial overlap in wintering areas. That result supports the hypothesis that segregation in the non-breeding season is linked to genetic structure. Philopatry and a neritic foraging habit probably also contribute to the structuring. Widely distributed, vulnerable to anthropogenic stressors, and exhibiting extensive genetic structure, the rhinoceros auklet is fully indicative of the scope of the conservation challenges posed by seabirds.

Genetic structuring among colonies of a pantropical seabird: Implication for subspecies validation and conservation

Investigations of the genetic structure of populations over the entire range of a species yield valuable information about connectivity among populations. Seabirds are an intriguing taxon in this regard because they move extensively when not breeding, facilitating intermixing of populations, but breed consistently on the same isolated islands, restricting gene flow among populations. The degree of genetic structuring of populations varies extensively among seabird species but they have been understudied in their tropical ranges. Here, we address this across a broad spatial scale by using microsatellite and mitochondrial data to explore the population connectivity of 13 breeding populations representing the six subspecies of the white-tailed tropicbird (Phaethon lepturus) in the Atlantic, Indian, and Pacific Oceans. Our primary aim was to identify appropriate conservation units for this little known species. Three morphometric characters were also examined in the subspecies. We found a clear pattern of population structuring with four genetic groups. The most ancient and the most isolated group was in the northwestern Atlantic Ocean. The South Atlantic populations and South Mozambique Channel population on Europa were genetically isolated and may have had a common ancestor. Birds from the Indo-Pacific region showed unclear and weak genetic differentiation. This structuring was most well defined from nuclear and mtDNA markers but was less well resolved by morphological data. The validity of classifying white-tailed tropicbirds into six distinct subspecies is discussed in light of our new findings. From a conservation standpoint our results highlight that the three most threatened conservation units for this species are the two subspecies of the tropical North and South Atlantic Oceans and that of Europa Island in the Indian Ocean.

Riverine barrier effects on population genetic structure of the Hanuman langur (Semnopithecus entellus) in the Nepal Himalaya

BACKGROUND

Past climatological events and contemporary geophysical barriers shape the distribution, population genetic structure, and evolutionary history of many organisms. The Himalayan region, frequently referred to as the third pole of the Earth, has experienced large-scale climatic oscillations in the past and bears unique geographic, topographic, and climatic areas. The influences of the Pleistocene climatic fluctuations and present-day geographical barriers such as rivers in shaping the demographic history and population genetic structure of organisms in the Nepal Himalaya have not yet been documented. Hence, we examined the effects of late-Quaternary glacial-interglacial cycles and riverine barriers on the genetic composition of Hanuman langurs (Semnopithecus entellus), a colobine primate with a wide range of altitudinal distribution across the Nepalese Himalaya, using the mitochondrial DNA control region (CR, 1090 bp) and cytochrome B (CYTB, 1140 bp) sequences combined with paleodistribution modeling.

RESULTS

DNA sequences were successfully retrieved from 67 non-invasively collected fecal samples belonging to 18 wild Hanuman langur troops covering the entire distribution range of the species in Nepal. We identified 37 haplotypes from the concatenated CR + CYTB (2230 bp) sequences, with haplotype and nucleotide diversities of 0.958 ± 0.015 and 0.0237 ± 0.0008, respectively. The troops were clustered into six major clades corresponding to their river-isolated spatial distribution, with the significantly high genetic variation among these clades confirming the barrier effects of the snow-fed Himalayan rivers on genetic structuring. Analysis of demographic history projected a decrease in population size with the onset of the last glacial maximum (LGM); and, in accordance with the molecular analyses, paleodistribution modeling revealed a range shift in its suitable habitat downward/southward during the LGM. The complex genetic structure among the populations of central Nepal, and the stable optimal habitat through the last interglacial period to the present suggest that the central mid-hills of Nepal served as glacial refugia for the Hanuman langur.

CONCLUSIONS

Hanuman langurs of the Nepal Himalaya region exhibit high genetic diversity, with their population genetic structure is strongly shaped by riverine barrier effects beyond isolation by distance; hence, this species demands detailed future phylogenetic study.