Next Article in Journal
Bajacanthon, a New Subgenus for the Mexican Deltochilini (Coleoptera: Scarabaeidae: Scarabaeinae) Fauna
Previous Article in Journal
The Nurse Plant Acacia spirorbis Enriches Ectomycorrhizal Community Composition of a Target Species: Tristaniopsis calobuxus
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Misael Daniel Mancilla-Morales
1,
Enriqueta Velarde
2,
Andres Aguilar
3,
Araceli Contreras-Rodríguez
4,
Exequiel Ezcurra
5,
Jesús A. Rosas-Rodríguez
6,
José G. Soñanez-Organis
6,* and
Enrico A. Ruiz
1,*
1
Departamento de Zoología, Escuela Nacional de Ciencias Biológicas, Instituto Politécnico Nacional, Colonia Santo Tomás, Ciudad de México 11340, Mexico
2
Instituto de Ciencias Marinas y Pesquerías, Universidad Veracruzana, Hidalgo 617, Colonia Río Jamapa, Boca del Río 94290, Mexico
3
Department of Biological Sciences, California State University, Los Angeles, CA 90032, USA
4
Departamento de Microbiología, Escuela Nacional de Ciencias Biológicas, Instituto Politécnico Nacional, Colonia Santo Tomás, Ciudad de México 11340, Mexico
5
Department of Botany and Plant Sciences, University of California, Riverside, CA 92521, USA
6
Departamento de Ciencias Químico-Biológicas y Agropecuarias, Universidad de Sonora, Colonia Francisco Villa, Navojoa 85880, Mexico
*
Authors to whom correspondence should be addressed.
Submission received: 29 December 2021 / Revised: 22 January 2022 / Accepted: 28 January 2022 / Published: 2 February 2022
(This article belongs to the Section Biodiversity Conservation)

Abstract

:
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.

Graphical Abstract

1. Introduction

Philopatry, sensu stricto, is the return of individuals to their natal site for breeding while dispersal is the opposite, and, in its original sense, a high degree of philopatry implies minimal gene flow [1,2,3]. In turn, gene flow means individuals moving into and out of a population, influencing the partitioning of genetic variation among populations [4]. The pattern of movement of individuals can profoundly affect the genetic structure of populations; as generations go by, the fact that individuals are being born and breed in the same places should result in the formation of discrete populations, suggesting that the stronger the philopatry, the stronger the genetic structure [3,5]. If so, genetic signatures thereof can persist among separated populations within widely distributed and highly mobile species [6,7,8]. Nonetheless, not all the birds are philopatric, and, by extension, not all show reduced gene flow, as this phenomenon can vary significantly per species, and even within species [9,10]. To better describe the connection between dispersal and genetic structure, it is advised the combination of methods such as mark-recapture and molecular markers [5,11,12].
It is a generality that, in birds, philopatry is biased towards males, while dispersal is biased towards females [1,13]. Direct observational methods allow one to estimate the sex ratio of a set of individuals returning to or leaving their natal site for breeding, leading to accurate estimates of dispersal if migrants have successful reproduction [14]. Philopatry or dispersal behavior reinforced by sexual selection may have some demographic, social, and genetic consequences, since it can determine the magnitude of gene flow, influence the social organization, and affect the genetic structure of populations [15,16,17,18]. For example, in the Magnificent Frigatebird (Fregata magnificens), a seabird with high dispersal capacity, geographic distance or physical barriers do not play an important role in determining genetic structure. Rather, their behavior traits (male-biased philopatry, female choice, and sequential monogamy) are more consistent with genetic differentiation among Eastern Pacific breeding colonies [19]. On the other hand, indirect methods that rely on molecular markers allow us to assess gene flow and genetic structure, making it possible to infer barriers among populations [20]. In this sense, two scenarios can be expected: (1) Weak genetic differentiation provided by the seabird’s ability to fly long distances, leading to panmixia [21,22]. For example, in the Elegant Tern (Thalasseus elegans), high connectivity among populations was revealed between the Gulf of California and southern California, with no significant genetic differentiation, even though their colonies are separated by land barriers and breed in different ocean regimes [23]. (2) Genetic structure can be noticed even in small areas due to philopatry and local selection pressures [5,24]. For example, the Least Storm-Petrel (Hydrobates microsoma) and the Black Storm-Petrel (Hydrobates melania) from Partida Norte Island showed genetic structure signatures nesting in the same island (possibly due to philopatry) [25].
Seabirds can exhibit genetic structure due to several factors, such as philopatry, local differences in selection pressure, and isolation by distance, which may constrain gene flow [5,19,26,27]. Therefore, investigating their genetic diversity can help us to elucidate processes and patterns of diversification [28,29]. In Charadriiformes, as in other long-distance migrants, ecological factors could better explain natal philopatry than genetic ones, although both factors may be interacting and even caused by sociality [30,31,32]. Moreover, ocean regimes and nesting habitats can exhibit differences in food availability and predation, both are important factors in reproductive decisions that influence fitness outcomes [33,34,35]. Thus, islands are important refuges because they bring protection from predators, while inter-colony segregation can occur if food resources are available in separated foraging areas [36,37,38]. Hence, local habitat differences may constrain the genetic exchange between populations, favoring ecological segregation by benefits of philopatry in high quality areas and discouraging movement to new environments through reduction of fitness by selection against dispersers [39,40,41,42].
Heermann’s Gull (Larus heermanni) reproductive ecology is coupled to global and local oceanographic conditions [43,44,45,46]. Historically, this gull has been reported nesting on several islands in the Pacific Ocean, although presently do so in only a few of them [47,48]. The two largest active colonies are in the Rasa (~130,000 pairs) and Cardonosa (~4000 pairs) islands, both located in the Midriff Islands Region [49]. Another important nesting site is Isabel Island (~400 pairs), located in the Mexican Province region [48]. Both regions have noticeable differences based on physical, chemical, biological oceanography, and other ecological criteria [50,51]. For example, the Gulf of California has deeper bathymetry, larger nutrient flux by upwellings in the Midriff Islands, and larger phytoplankton biomass [43,51,52,53]. Our aim was to infer the extent of Heermann’s Gull philopatry in its main breeding site (Rasa Island), as well as to examine the genetic structure between colonies near the northernmost and southernmost nesting limit of the species.

2. Materials and Methods

2.1. Study Area

The study was performed in Heermann’s Gull colonies from two biogeographic regions characterized by distinguishable physical, chemical, and ecological traits: The Midriff Islands Region and the Mexican Province [48,51,52,53] (Figure 1). First, to assess the dispersal behavior in Heermann’s Gull, we used data from Rasa Island banding in nine natal valleys between 1984 to 1993 (ES, Valle del Estero; IT, Vallecito del Estero; CA, Valle de la Casita; GE, Gran Estación Central; VG, Valle de los Gallitos; TV, Tapete Verde; CR, Casitas Viejas Arriba; CB, Casitas Viejas Abajo; and VW, Valle de la W), and monitoring in seven recapture valleys between 1989 and 2013 (see results). Next, to examine their genetic structure, cytochrome b sequences were obtained at twenty locations from Cardonosa (20 samples collected in 2011), Rasa (266 samples collected in 2011 and 2012), and Isabel (10 samples collected in 2017) islands (Figure 1).
Cardonosa Island is in the Midriff Islands Region of the Gulf of California (28°53′16″ N and 113°01′51″ W) [46,49]. It is an irregularly flattened island, 483 m long, 292 m wide, 18 m high, and ~0.14 km2 in area [43]. The vegetation consists mainly of the coastal Saltbush (Atriplex barclayana), the Baja California cholla (Cylindropuntia alcahes), and the cardon (Pachycereus pringlei) [25,54,55]. Here, Heermann’s Gulls nest in two places: up in an eroded flat area covered by guano and down on a boulder and pebble beach (Figure S1a,b). We sampled the second location in ‘Playa Rocosa’, a rugged area with scarce vegetation (Figure S1c).
Rasa Island is located 4 nautical miles SE of Cardonosa, in the same biogeographic region (28°49′29″ N and 112°58′50″ W) [48,49,50]. It is a young flattened volcanic island (~10,000 years), 1170 m long, 703 m wide, 30 m high, and ~0.68 km2 in area [43,54]. The relief consists of low flat areas covered by guano (valleys) and rocky slopes (hills), with relatively little change in the last fifty years [43,47,54] (Figure S2a–c). The vegetation consists of salt scrub in three esteros and desert scrub in the flats and hill areas, the last with a notable increase of chain-fruit cholla (Cylindropuntia fulgida) [47,54]. On Rasa Island, we sampled mainly in valley areas and surrounding rocky hills. Predation by the Yellow-footed Gull (Larus livens) (Figure S2d) and Common Raven (Corvus corax) has been observed in Heermann’s Gull chicks and eggs, and Peregrine Falcon (Falco peregrinus), also preys on adults and chicks of Heermann’s Gulls [43,47,55,56,57,58]. In normal years, initially, thousands of Heermann’s Gulls arrive at the nesting area, and, a few days later the Royal Tern (Thalasseus maximus) and the Elegant Tern (Thalasseus elegans) arrive in large numbers, displacing many gulls from their valley nesting areas and establishing their nests in a frequently mixed colony [44] (Figure S2a).
Isabel Island is located off the coast of Nayarit, in the Mexican Province (21°50′40″ N and 105°52′49″ W) [48,50]. This island is 1570 m long, 100–850 m wide, 30–85 m high, and ~1.06 km2 in area [43]. Isabel Island has cliffs with steep slopes, rocky and sandy beaches, and vegetation patches with dense grass and trees areas of low deciduous tropical forest with strong seasonality [59]. Here, we found an area with approximately fifty scattered nests of Heermann’s Gull near Brown Booby (Sula leucogaster) nests, over a mixed terrain of grass and solid tuff-breccia rock [59] (Figure S3a,b). Gull chicks were reared more successfully on nests among vegetation than in hard rock, and it is relatively common to find adult gulls stalking nests and consuming Blue-footed Booby (Sula nebouxii) eggs (Figure S3c), whereas the insular Milk-snake (Lampropeltis triangulum sinaloae) may also feed on gull chicks (Figure S3d).

2.2. Analyses of Philopatry and Dispersal by Banding Data

Heermann’s Gull nested synchronously in large numbers and was monitored on Rasa Island as part of other studies [44,45,46,60]. The banding was carried out in large flocks of fledglings moving together and restricted in chicken-wire enclosures in each natal valley (Figure 1). Inside the enclosure, chicks were captured by hand by four banders, banded, and immediately released to the outside of the enclosure. Disturbance was avoided by waiting till chicks had reached fledgling stage, since fledglings are already moving freely from parents and territory at that age, returning to their territory only to be fed by parents. We recaptured all the individuals included in the analysis, and band identification was based on numerical characters, not colors. The recaptures were performed using walk-in traps placed on their nests [60]. Sex assignment of recaptured Heermann’s Gull adults was carried out through mating behavior and somatic measurements (i.e., within members of each pair, males are heavier [60] and larger [61] than females).
To avoid duplicating records of individuals observed in more than one season, our calculations considered gulls that were observed breeding only for the first time they were recaptured. We organized a data matrix containing (a) the natal valley (where the banded bird was born), (b) the breeding valley (where the recovered adult was found nesting), (c) the distance between the natal and the breeding valleys (with reference to their location in Rasa Island; Figure 1), (d) the sex of the recaptured birds, (e) the re-sampling effort (the number of seasons where recaptures were performed in each breeding valley, (f) the number of birds from each specific natal valley recaptured in each different breeding valley, (g) the total number of birds recaptured in each breeding valley, and (h) the total number of birds banded in each specific natal valley (Supplementary Material S2).
If strong natal philopatry occurs in Heermann’s Gull, we would expect nesting banded birds to be recaptured in, or close to, their natal valley, i.e., the recovery of nesting banded birds would be a negative function of distance to the natal valley [3,10]. Because different number of birds were banded in each valley, we analyzed the proportion of birds recaptured in each breeding valley that had been originally banded in each natal valley. The analysis was done through a Generalized Linear Model for binomial proportions, with logit link and Chi-squared deviance. Using a stepwise additive model, we tried as predictors of the proportion of birds breeding in a given valley: (a) the distance between natal and breeding valleys, (b) the fixed effect of the individual’s sex, (c) the interaction between distance and sex, (d) the effect of re-sampling effort in the breeding valley, (e) the fixed of the natal and breeding valleys, and, lastly, (d) the factorial interaction between natal and breeding valleys.
Because at the time of banding the sex of chicks cannot be determined, we did not know the number of originally banded chicks in each sex in order to calculate sex-specific proportions of banded-to-recaptured. For this reason, we performed three different analyses under different assumptions to examine whether sexual bias is occurring (males prone to philopatry and females to dispersal). Model I: We first analyzed the proportion recaptured birds in each breeding valley in relation to the total number of recaptured birds from each natal valley. Model II: We assumed that the original sex ratio in the banded birds had been similar to the ratio found in the recaptured birds. Model III: Lastly, we also performed the analysis assuming that the original sex ratio in the banded birds had been 50%–50%. All analyses were done in the R package [62], using the function glm (Generalized Linear Model) with binomial error and logit link (see the logistic model script in Supplementary Material S3).

2.3. DNA Extraction, PCR Amplification, Sequencing, and Analyses of Genetic Diversity

DNA extraction, PCR amplification and sequencing was performed only with samples from Isabel Island (n = 10), which in turn were the most recently collected in 2017. All other sequences were those from a previous study by Ruiz et al. [49] and downloaded from GenBank (see below). Mitochondrial cytochrome b (Cyt-b) sequences were obtained from DNA isolated from ~60 µL of blood from the ulnar vein of Heermann’s Gulls using the DNeasy® Blood & Tissue kit (QIAGEN, Hilden, Germany). PCR reactions were performed in a final volume of 25 µL containing 50 ng of genomic DNA (1 µL), 10 mg mL−1 BSA (1 µL), DMSO 2% (1 µL), deionized water (10 µL), 2× Kapa Taq ReadyMix (Kapa Biosystems, R&D Cape Town, ZA) (10 µL), 10 mM primers L14841 (1 µL), and 10 mM primers H16065 (1µL) [49]. The PCR program and primers were previously described by Ruiz et al. [49], and the PCR products were outsourced for sequencing to Macrogen Inc. (Seoul, Korea). Then, the recovered sequences (1033 bp) were aligned and reviewed using Seaview v.4.5.4 [63]. Next, we identified the sequences as Cyt-b, by comparing them with previous GenBank sequences from Rasa and Cardonosa (accession numbers KT275317-KT275602; [49]). Finally, newly obtained sequences from Isabel Island were submitted to GenBank (accession numbers MT935687-MT935696). Jukes-Cantor 69 was selected as the best model fitting to molecular evolution by jModelTest v.2.1.7 [64]. Haplotype richness (h), number of segregating sites (S), number of nucleotide differences (k), haplotype diversity (Hd), and nucleotide diversity (π) were estimated using DnaSP v.6.12.01 [65]. Relationships among haplotypes from each island were inferred with a minimum spanning network by PopART v.1.7 [66]. To our knowledge, there is only another Cyt-b sequence for Larus heermanni from Grays Harbor, USA [67,68]. However, this was excluded from our analysis in the main text because it was shorter (290 bp) and the different diversity indices cannot be recovered from a single sequence. Instead, a phylogeny incorporating all the Cyt-b sequences reported to date is shown in the Supplementary Material (Table S2 and Figure S4). Phylogenetic relationships were assessed using the automatic selection of substitution model [69], subtree pruning and regrafting (SPR) tree searching [70], and aLRT SH-like branch support [71], all the tools implemented in PhyML 3.0 [72].

2.4. Genetic Differentiation

Our next aim was to explore the extent to which Heermann’s Gull colonies are genetically subdivided into subpopulations. Nonetheless, different statistical approaches vary in their power to detect genetic divergence due to sample size and genetic markers [73,74,75]. Thus, to better understand Heermann’s Gull population structure, different genetic differentiation tests were applied at different hierarchical geographic levels: regions, islands, clusters, and locations. Haplotype-based (χ2, [74]; HST [76]) and nucleotide sequence-based statistics (KST* and Z*, [76]; Snn [77]) were calculated using 1000 permutations through DnaSP v.6.12.01 [65]. Moreover, the genetic differentiation between populations at each geographic level was examined by statistics considering gene diversity (GST [76,78]), nucleotide diversity (NST [79]), and global partition from random mating departures (FST [80,81]). At the same time, the effective number of migrants per generation (Nm), was inferred for each statistic (GST, NST, and FST), based on the island model using DnaSP v.6.12.01 [65,80,82].
On the other hand, to assess their significance, the pairwise ΦST values among the twenty locations were examined with an exact test of population differentiation, with 100,000 steps of Markov chains and 10,000 permutations, using Arlequin v.3.5.2.2 [83].

2.5. Genetic Structure

To assess population structure in Heermann’s Gull, three approaches were followed. In all of them, the twenty locations sampled were regarded as separated populations (Figure 1). First, we assumed an a priori definition by analyzing molecular variance (AMOVA), which was used to test the statistical significance of the genetic partition within groupings inferred by geography and proximity (regions, islands, and clusters), with the Φ-statistics from haplotypic data using Arlequin v.3.5.2.2 [83]. Next, we used an a posteriori approach (SAMOVA) to search for discontinuity between the twenty locations; this analysis indirectly detects barriers and defines groups maximally differentiated [84]. With this second approach, the number of groups (K = n − 1) was identified using the highest FCT value. Significance tests were performed with 1000 permutations using SAMOVA v.2.0 [84]. Finally, a Bayesian analysis of population structure based on linked molecular information was carried out through clustering algorithms to assign individuals to clusters [73,85]. This third approach tests the population structure across populations employing a clustering model implemented in BAPS v.6 [86]. The clustering was implemented, running five independent replicates with linked loci, assuming a maximum number of groups (K = from 2 to 20) without specifying the sampling populations to geography.

2.6. Gene Flow

Migration rates between islands were estimated following a Bayesian approach implemented in Migrate v.4.4.3 [87]. The mutation scaled effective population size (Θ = xNeµ) and the estimate of migration per generation (M = m/µ) were the demographic parameters reported; both assume populations in migration-drift equilibrium [87]. This method (based on coalescence theory and a Monte Carlo Chains algorithm (MCMC)) considers the genealogical relationships of the samples to estimate demographic parameters, and provides reliable estimates of asymmetric gene flow [87]. The analysis was done using the Jukes-Cantor model, a constant mutation rate, and one long-chain of 10,000,000 steps (with three simultaneous replicates and 10% of the steps discarded as burn-in). To increase searching effectiveness, a static four-chain heating scheme was used with temperatures set to 1.0, 1.5, 3.0, and 1,000,000. Results were examined from independent runs, and the basic model of migration (migration free to vary between all the populations) was selected as the best model found having the highest probability, and each execution was accepted only when the estimates of effective sample sizes (ESS) were >1000. Pairwise comparisons among all locations could be problematic due to overparameterization and can inflate the effective population size estimates. Therefore, avoiding the violation of assumptions and to get more realistic gene flow estimates, we only made comparisons among islands considering each one as a population [87].

2.7. Correlation of Genetic Differentiation with Geographic Distances

To explore potential geographical isolation among populations, we used a Mantel test [88]. Under the assumption of isolation by distance (IBD), we expect genetic distance to increase with geographic distance. The Mantel test considers an island model, which assumes that the closer the individuals, the closer their relationship and, therefore, genetically more dissimilar from those further away [89]. Thus, the correlation between the matrix of linear geographic distances (natural logarithm of linear distances in kilometers) and the matrix of genetic distances (ΦST/1 − ΦST [83,90]) was used according to Rousset [89]. Each location in the three islands was considered as a different population, and 1000 permutations were used to assess its significance.

3. Results

3.1. Philopatry and Dispersal Behavior

A total of 998 Heermann’s Gulls (714 males and 284 females) were recaptured in 2721 occasions (males = 1917 recaptures, and females = 804 recaptures) breeding from one to fourteen times as adults during their lifetime at seven monitored valleys in Rasa Island (Table 1). In total, 68.34% of the individuals were faithful to their natal valley for breeding (n = 682 gulls), often returning to the same nesting site or only a few meters away (namely, natal philopatry [3]). The remaining 31.66% (n = 316 gulls) dispersed to breed in other areas on Rasa Island and repeatedly returned to that same valley for breeding (namely, breeding philopatry [3]).
In all three models (Table 2), the best predictor of the location of the breeding site was the distance to the natal site with a negative coefficient, i.e., nesting birds tend to nest in the same valley where they were born, or in the close neighboring valleys (r2 = 0.83 for Model I; 0.76 for Model II; and 0.72 for Model III; p < 0.0001 in all cases). Expectedly, the resampling effort had a small, but significant effect in all models (p ≤ 0.001 in all cases). Nonetheless, the interaction between sex and distance was also highly significant in all models, with a more pronounced negative slope for males than for females, indicating that males tend to be more philopatric. A significant interaction term between the natal and breeding valley suggests that some valleys tend to expel their offspring onto other valleys when they become adults, possibly due to differences in habitat quality. For example, Casitas Viejas Arriba is a significant (p < 0.001) expeller of offspring, many of which tend to choose Tapete Verde and Valle de la W as their preferred alternative for breeding as adults (see Table 1). Similarly, Valle de la W showed a very significant (p < 0.001) interaction term with itself, indicating that philopatry is higher here than in other natal sites, i.e., chicks raised in the valley tend to come back as adults to the same place avoiding other alternative sites. Finally, Model III, which assumed that banding at birth was done on a 50% ratio of males to females, indicated a very significant fixed effect of sex, suggesting that, if indeed sex ratio at banding was 50:50, then males are way more philopatric than females. Overall, the proportion of nesting adults that dispersed showed an inverse relationship with dispersal distance regardless of the model used, especially in males compared to females when assumed natal sex ratio is 1:1 (Model III) (Figure 2).

3.2. Genetic Diversity

Forty-two cytochrome b haplotypes were observed in 296 Heermann’s Gulls from the Midriff Islands Region and the Mexican Province (Table 3). Thirty-five haplotypes belonging to the Midriff Islands Region showed several segregating sites (S = 41), few nucleotide differences (k = 0.85), high haplotype (Hd = 0.61), and low nucleotide diversities (π = 0.00083). Conversely, only eight haplotypes were recovered from the Mexican Province, and showed fewer segregating sites (S = 12); nonetheless, there are higher values of nucleotide differences (k = 4.133), haplotype diversity (Hd = 0.93), and nucleotide diversity (π = 0.00401). In Cardonosa Island, we found only three haplotypes, two segregating sites (S = 2), and the lowest values of nucleotide differences (k = 0.54), haplotype diversity (Hd = 0.51), and nucleotide diversity (π = 0.00053). In Rasa Island, we found most of the haplotypes (h = 34), segregating sites (S = 40), and nucleotide differences (k = 0.88). In Isabel Island, we recovered eight haplotypes and higher haplotype and nucleotide diversity values than in Cardonosa and Rasa Islands. Despite sharing some haplotypes, we observed unique haplotypes on each island occurring in more than half of the locations.
The minimum spanning network indicates that only the haplotype H1 is shared among Cardonosa, Rasa, and Isabel Island (Figure 3), whereas the haplotype H2 was observed in both Cardonosa and Rasa Island, but not in Isabel Island. No other haplotype was observed in more than one island. Several low-frequency haplotypes with one to three mutational steps are derived from the two most common haplotypes (H1 and H2), and several low-frequency haplotypes with one to three mutational steps are derived. The occurrence of haplotype H1 in all three islands suggests that it is either very old and ancestral or a signal of genetic exchange even over long distances, or both. However, unique haplotypes suggest slight differentiation between islands, and, coupled with the resulting radial topology of relationships, there are signals of lineage sorting after a recent population expansion.

3.3. Genetic Differentiation

Comparisons between populations from different regions in all the tests revealed significant differentiation (Table 4). However, some tests revealed no significant differentiation regarding the lower geographic level. The Chi-square test (χ2) showing large values suggested that the genetic differences among populations are better explained by selection or drift. Heterogeneity of haplotype frequencies among populations was detected at region, island, and cluster levels. Nonetheless, the Chi-square test is less robust if there are too many unique haplotypes or if the number of each haplotype in each location is small [74,76,77]. In contrast, the HST statistic showed low values with statistical significance (p < 0.01), indicating considerable diversity within each sample analyzed, and slight movement of haplotypes [76,77,86]. The HST statistic only showed significant genetic differentiation between regions. This could be explained by rare haplotypes grouped together and reflects the effect of population structure on the haplotype diversity [76]. Similarly, the values of KST* near to zero indicated no differentiation; however, they were statistically significant (p < 0.01), supporting genetic differentiation [91,92]. The sequence-based statistic KST indicated differentiation at all the geographic levels. The Z* statistic was also statistically significant (p < 0.01), so the null hypothesis (no differentiation) was rejected [92]. This statistic revealed significant differences between regions and islands; however, it detected no considerable differentiation among clusters and locations. The sequence-based statistic Snn revealed robust differences among regions and islands but a high genetic exchange among clusters and locations. Overall, the higher the geographic level, the more significant and stronger the differentiation observed in Heermann’s Gull subpopulations analyzed.
The GST, NST, and FST statistics showed the strongest genetic differentiation at the regions and islands levels; however, it was lower at the finest geographic levels (clusters and locations) (Table 4). Conversely, the number of migrants (Nm) increased from higher to lower geographical level, a reflex of retention of homogeneity in the populations by lower exchange of individuals among regions and stronger exchange between locations. A widely accepted rule is that one migrant per generation will maintain a certain level of connectivity among subpopulations [93]. However, one to ten could be an appropriate number to minimize the loss of diversity within subpopulations [94,95]. In this way, despite the number of migrants per generation observed by nucleotide data (Nm2 = 1.71 to 7.73, and Nm3 = 1.71 to 7.75), it was not sufficient to avoid genetic differentiation among regions based on their nucleotide diversity (NST = 0.227) or random global departures of allele frequencies (FST = 0.226). In contrast, higher migration rates (Nm1 = 17 to 145) showed lower genetic differentiation (GST < 0.028).
Paired comparisons of genetic differentiation between locations revealed significant differences in few instances, mostly between both Rasa and Cardonosa Islands and Isabel Island (Table 5). Global FST or ΦST values in seabirds are usually lower than the maximal theoretical (FST or ΦST = 1) [5,75], and we interpreted their values as low (0 to 0.05), moderate (0.05 to 0.15), high (0.15 to 0.25), and very high (>0.25) [96]. Thus, the lower the value, the higher the gene flow. Low or negative ΦST values were interpreted as low genetic differentiation and high genetic exchange between populations [83,90], i.e., the ΦST comparisons were very high (ranging 0.21 to 0.32) for the most distant population, Costa Fragata, versus all other locations, although not all were statistically significant. Comparisons between locations from Cardonosa and Rasa are predominantly undifferentiated, mainly with negative or zero ΦST values. However, low to moderate genetic differentiation within Rasa Island was observed between Valle de la W versus Vallecito del Estero, Valle Norte, Valle del Fin, Casitas Viejas Abajo, and Casitas Viejas Rocas, with significant ΦST values ranging from 0.09 to 0.22. Only Valle de la W from Rasa showed significant differentiation versus Cardonosa. Overall, more prominent differentiation between some locations than others may be due to differences in migration between breeding places, as revealed by banding on Rasa Island (see Table 1).

3.4. Genetic Structure

The results of the AMOVA analysis are summarized in Table 6. Different patterns of structuring were found depending on groupings considered. For example, the inferred genetic sub-structuring among regions explained 48.57% of the variation, whereas comparisons among islands and clusters were 24.16% and 6.95%, respectively. However, the genetic differentiation was very high and significant among regions (ΦCT = 0.49, p = 0.05), although it was high, but not significant, between islands (ΦCT = 0.24, p = 0.10), and moderate, but not significant, among analysis clusters (ΦCT = 0.07, p = 0.12). Furthermore, from very high to moderate genetic differentiation was detected within populations (locations) and significant in all the groupings (regions, ΦST = 0.48; p = 0.0001; islands, ΦST = 0.23, p = 0.0001; and analysis clusters, ΦST = 0.05, p = 0.0001).
The SAMOVA analysis revealed two maximally differentiated groups from each other when considering the twenty Hermann’s Gull locations (Table S4). With this approach, we identified that location 1 (Costa Fragata, on Isabel Island), was grouped alone and against another group formed by all the other locations (K = 2, FCT = 0.486, p < 0.05). Similarly, the next maximally differentiated grouping was formed by the locations 7 (Valle de la W, Rasa Island), 20 (Costa Fragata, Isabel Island), and all other locations grouped (K = 3, FCT = 0.317, p < 0.05). All other K groupings showed lower, but significant, FCT values ranging from 0.105 to 0.177 (Table S4). Likewise, the Bayesian approach to estimate genetic structure implemented by BAPS also recovered two genetic groups (K = 2, log ML = −824.58, and KL divergence = 0.014). The larger one encompassed 283 individuals, mainly from gulls nesting on Cardonosa and Rasa Island (Midriff Islands Region), while the smaller one was formed by only 13 individuals, mostly from Isabel Island (Mexican Province region) (Figure 4a). Altogether, the three approaches are consistent in revealing strong genetic structure among Heermann’s Gull breeding colonies nesting separately between the Midriff Islands Region and the Mexican Province biogeographic regions. Such genetic structure may result from differences at local habitat, disrupting the gene flow among Heermann’s Gull subpopulations from both oceanic regimes.

3.5. Estimation of Gene Flow between Islands

The estimated effective population size was lower for Cardonosa (Θ = 0.004; 95% CI = 0.004–0.007) than for Rasa (Θ = 0.066; 95% CI = 0.063–0.083) and Isabel Island (Θ = 0.016; 95% CI = 0.005–0.024) (Figure 4b). The scaled gene flow was mostly non-symmetrical among islands. For example, higher migration was observed from Cardonosa to Rasa (MC→R = 962; 95% CI = 902–1000), while migrants from Rasa to Cardonosa were at least an order of magnitude smaller (MR→C = 60; 95% CI = 35–137). In addition, Cardonosa and Rasa Island showed high and similar migration toward Isabel Island (MR→I = 405; 95% CI = 33–807; MC→I = 641; 95% CI = 579–999), but reduced migration from Isabel toward Rasa and Cardonosa (MI→R = 29; 95% CI = 25–75; MI→C = 72; 95% CI = 39–152). Altogether, the scaled migration (M) was consistent with our other estimates of gene flow (Nm) and genetic differentiation (FST), as they also show a lower number of migrants and high differentiation among regions, indicating lower dispersal to remote regions (Table 4 and Table 5).

3.6. Isolation by Distance

The correlation between pairwise genetic distances (ΦST/1 − ΦST) and geographic distances was substantial and statistically significant (rM = 0.82, p = 0.03) (Figure 4c). Thus, an IBD model (isolation by distance), in addition to biogeographic barriers, could explain the mitochondrial genetic structure found in all the previous analyzes. Distance among locations ranged from a few meters to more than 1000 km (Table 5).

4. Discussion

This study provides the first evidence of population genetic structure among Heermann’s Gulls breeding on islands of different biogeographic regions [48]. As in other studies of seabirds that breed in two or more ocean regimes [5,19,21,26,27,28,29], our results support strong philopatry, local habitat, and distance as factors that reduce the gene flow among regions.

4.1. Philopatry and Dispersal Behavior

The results of our analysis confirm that Heermann’s Gulls are strongly philopatric tending to return as nesting adults not only to the island where they were born but specifically to the island valley where they were hatched and reared, or at a very close distance to it (Figure 2). Interestingly, the best predictor of the adult breeding ground was not the natal valley as a discrete factor but, rather, the quantitative distance to the place of birth (Table 2). A high proportion of Heermann’s Gull individuals were recaptured within 200 m of their hatching site, while a few moved greater distances elsewhere, as in the case of the Black-legged Kittiwake (Rissa tridactyla) [10]. This pattern may not be general in all species of birds, but there seems to be a strong tendency for philopatry among seabirds [5,75].
A second extremely consistent result, under any model hypothesis, is that males are more strongly philopatric than females (Table 2). Therefore, our study system is consistent with findings in other seabird species [1,10]. In seabirds, generally, the males will display higher philopatry than females, while the females will tend to nest farther away than males from their natal site [1,14]. Description of other patterns, however, have been found in land birds, such as the Brown Jay (Psilorhinus moio), in which philopatry is biased to females and dispersal to males [16]. Nonetheless, any interpretation should be taken with caution because any bias in the production of either sex or their survival (for which we have no control) may contribute to differences in dispersal among sexes. Such biases are not uncommon, as sex-biased production and mortality have been reported in other gulls [97,98,99]. For example, in the Black-headed Gull (Larus ridibundus), male-biased mortality has been reported: as males are larger, one possibility is that they are more prone to not fulfill their resource requirements [97]. Moreover, there is evidence that dimorphic birds overproduce the smaller sex under adverse conditions (usual females in Laridae) [61,98]. If so, females would be captured more than males because of this overproduction [98], which contrasts with our results.
We do not know the sex of the fledglings we banded initially. For this reason, we built three different models with varying fledgling sex ratio. Later, we recaptured more adult males than females, even though this is a slightly dimorphic species, with males slightly heavier and larger than females [51,61]. This can be unexpected because, as in the Black-headed Gull, males could have a higher mortality rate and, therefore, should be recaptured less [97]. However, selective factors in one species can be different from those of another, and effects may well vary [38]. Each species is in an environment with its own characteristics and natural selection pressures; in addition, each species may differ in the ways of reacting to these pressures, so we cannot expect to see the same results in two distinct species; on the contrary, what is interesting is to see the differences between them. A possible explanation for what we found is given by the Trivers-Willard hypothesis, which assumes differences in fitness returns with higher investment in the sex that gains more benefits [99]. In this way, an important part of the development of the natal population can be the more philopatric gender [100]. Then, as in the Banter See Common Tern (Sterna hirundo), male-biased natal recruitment in Heermann’s Gull suggests a greater influence of males in the natal population growth than their female equivalents [100]. The above suggests that more attention should be paid to gender bias in future studies of seabirds, as often only females are considered in population modeling and, in fact, males are left out.
Differences in the number of recaptured gulls with philopatry and dispersal behavior in each valley (see Table 1) are probably related to the quality of local habitat and interspecific interactions indirectly affecting the survival [38,39,101]. If some local habitats are more suited than others (i.e., with higher quality), then some valleys will show a higher proportion of gulls recaptured with natal philopatry returning every year (philopatry benefits) [38,39]. Conversely, dispersers (i.e., gulls coming from any other nesting area) should be captured less frequently (selection against immigrants) [42]. For instance, a recent study found evidence of selection against immigrants, which were found to have a lower fitness than residents in wild seabird populations of three long-lived Procellariiformes species (Wandering Albatross Diomedea exulans, Southern Fulmar Fulmarus glacialoides, and Snow Petrel Pagodroma nivea) [42]. Other interactions may include the deterring by successive occupation of other birds arriving in large numbers or predation (Figure S2), in addition to differences in egg sizes influenced by nest densities, habitat, and nest-site selection, as have been the cases for the Herring Gull (Larus argentatus) and the Common Gull (L. canus) [101,102]. Therefore, the nest-site selection must be a crucial decision because breeding synchronously and densely in certain areas might help to enhance offspring development and an effective defense against predation [30,31,32,33,34,35,36,37,38,44]. In Heermann’s Gull, as in many other species, producing offspring in the same environment experienced by the breeding adults at an early stage in life would be advantageous and favor the operation of group selection and local adaptation [26,30,31,32,103].
In conclusion, both natal and breeding place appear to play an important role with possible benefits for philopatric individuals, but more research is needed to corroborate this. At the end, we do not rule out any potential bias due to different return rates by mortality, our ability to recapture the birds, or even some bias in the behavior of birds with greater territoriality that allows them to be captured more easily. The latter is matter of future research.

4.2. Gene Diversity, Genetic Structure, and Gene Flow

In Heermann’s Gull, high haplotype diversity and low nucleotide diversity were observed in cytochrome b (Hd = 0.62 and π = 0.00098; Table 3). Both values of diversity indices are comparable, for instance, to mitochondrial ND2 and CR1 regions in the European Shag (Phalacrocorax aristotelis) (Hd = 0.95 and π = 0.006) [104], and to mitochondrial COI sequences in both the Black Storm-Petrel (Hd = 0.70 and π = 0.00291), and the Least-Storm Petrel (Hd = 0.95 and π = 0.00559) [25]. In seabirds with demographic expansions across a wide range, high haplotype diversity and low nucleotide diversities are probably the consequence of Pleistocene effects [49,105,106]. In our study system, unique haplotypes were found within each region, island, cluster, and even locations (Table 3, Figure 3); however, its distribution accounts for a small fraction of lineage sorting between the biogeographic regions Midriff Islands Region and Mexican Province. Overall, the diversity of cytochrome b haplotypes in Heermann’s Gull agrees with demographic expansion [49], and it could be explained by regional genetic differences [106].
Given the large distances between Rasa and Cardonosa Island versus Isabel Island (>1000 km), it was natural to expect more robust differentiation if distance plays a major role as a factor of genetic structuring. Conversely, if the migratory behavior of Heermann’s Gull individuals allows them to cover long distances during their lifespans (and there are no restrictions as to what places they can choose to nest and reproduce), then the genetic differentiation would be overturned by panmixia. Our results are consistent with the first scenario, as pairwise comparisons of ΦST showed very high values when comparing Heermann’s Gull sequences from locations on Rasa and Cardonosa against Isabel Island (ΦST = 0.26 to 0.32; Table 5). However, although the ΦST values were usually low when locations within Rasa and Cardonosa are compared (ΦST = 0), there were instances of very high values (ΦST = 0.22), indicating that not only geographic distance but other factors (e.g., behavior) are also playing a role in genetic structuring at fine scales. Our results are comparable to pairwise ΦST differences among Atlantic Kittiwake (Rissa t. tridactyla) colonies (ΦST = 0.10 to 0.40) [107] and between the Red-billed Tropicbirds (Phaethon aethereus) populations on islands from the Mexican Ocean Pacific (ΦST = 0.0 to 0.20) [22].
On the other hand, very high genetic structure was observed between Heermann’s Gull colonies from the Midriff Islands Region and the Mexican Province region. This was consistent through all three approaches implemented (AMOVA, ΦCT = 0.49, p < 0.05; SAMOVA, FCT = 0.486, p < 0.05; and BAPS, log ML = −824.58, and KL divergence = 0.014). This pattern of population genetic structure may be the result of differences between ocean regimes and intrinsic ecological barriers. For example, this was the case between Atlantic and Pacific populations of the Black-legged Kittiwakes (Rissa tridactyla) (AMOVA, global ΦST = 0.53, p < 0.05) [107], and between the Atlantic and Indian populations of the Thin-billed Prion (Pachyptila belcheri) (AMOVA, FST = 0.35) [21]. However, this is not always the case. For instance, low genetic structuring among groups was found by Rocha-Olivares et al. [19], between colonies of Mexican Frigatebirds (Fregata magnificens) from the Mexican Ocean Pacific and the Mexican Caribbean (AMOVA, ΦCT = 0.015 p = 0.05).
Although, in Heermann’s Gull, the number of migrants per generation was more than one (Nm = 2 to 8), this was insufficient to create a single panmictic and homogenized population according to the genetic structure observed (Table 6 and Table S3, Figure 4a). Moreover, the asymmetric gene flow observed between islands where Heermann’s Gull breeds showed a southward decrease in the exchange of individuals (Figure 4b). As in many other threatened seabirds, Heermann’s Gull are usually found in fragmented and partially isolated populations, in which cryptic speciation may be found [6,27,28,29,104]. Therefore, although the one migrant per generation rule simplifies assumptions about the connectivity and the loss of polymorphism among subpopulations, several factors may influence the genetic exchange between populations [5,93,94,95]. In this last sense, our genetic analyses showed that physical distance is an essential factor shaping genetic structure among Heermann’s Gull colonies breeding on islands of different ocean regimes (Figure 4c). Nonetheless, if the breeding success and reproductive decisions are occurring influenced by predation and multigenerational social interactions incentivized by inclusive fitness, then behavior may also play an important role in the formation of genetic structure in the Heermann’s Gull as well as in other species [15,32,103]. Altogether, our results agree with the rationale of Friesen [5], as geographically separated populations (both by physical and non-physical barriers) tend to be less connected by gene flow, preventing dispersal and leading to very high genetic structure.

5. Conclusions

In Heermann’s Gull, our banding data showed strong and male-biased natal philopatry in Rasa Island, while the genetic data revealed consistent genetic differentiation between colonies breeding in the Midriff Islands Region and the Mexican Province. Furthermore, despite their high dispersal potential, Heermann’s Gull displayed uneven gene flow and positive isolation by distance pattern between populations. The strong genetic structure observed among regions is supported by unique haplotypes in each subpopulation and may reflect philopatry as a crucial behavioral barrier reinforced by physical isolation (distance) and different local habitat conditions.

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/d14020108/s1; Figure S1. (a) Heermann’s Gull colony in Cardonosa on 19 June 2017 (up and down yellow arrows indicate the sampled nests). (b) Up, Heermann’s Gull nest in an elevated and flat area. (c) Down, Heermann’s Gull was trapped nesting in the boulder and pebble beach ‘Playa Rocosa’ on 23 April 2011. Photo credits: Misael Mancilla (a, b) and Julia Whitty (c); Figure S2. (a) Mixed nesting colonies of Heermann’s Gull (Larus heermanni), the Elegant Tern (Thalasseus elegans) and the Royal Tern (Thalasseus maximus) (drawn dotted line) in the nesting site Gran Estación Central of Isla Rasa on 26 May 2017. (b) Heermann’s Gull chicks at rocky hill area. (c) Adult gulls nesting on relatively soft soil formed by guano at valley area. (d) Predation of Heermann’s Gull chicks by Yellow-footed Gull (Larus livens). Photo credits: Misael Mancilla (a, b, c, d); Figure S3. (a) Heermann’s Gull small colony was observed on Isabel Island on 9 April 2017 nesting close by Brown Booby (Sula leucogaster) nests, (b) gull chicks in nest of grassland habitat, and (c) adult gull predating eggs of the Blue-footed Booby (Sula nebouxii). (d) Insular Milk-snake (Lampropeltis triangulum) on Isabel Island. Photo credits: Misael Mancilla (a, b, c, d); Figure S4. Phylogenetic relationships based on mitochondrial cytochrome b sequences of Heermann’s Gull from Table S2. Letters at the tip of each label indicates the geographical origin: Grays Harbor (_Grays), Cardonosa (_Car), Rasa (_Ras), and Isabel (_Isa). The complete cytochrome b from Rissa tridactyla was used as the outgroup (GenBank: DQ385229). The HKY+G model was selected as the best substitution model. For each branch, the aLTRT support is shown. The most likely tree is shown, and their log likelihood tree was −2287.969; Table S1. Number of Heermann’s Gull fledglings ringed in nine Rasa Island valleys between 1984 and 1993; Table S2. Comparison of sequences reported in the literature for Heermann’s Gull cytochrome b; Table S3. Accession numbers for each haplotype of Heermann’s Gull in the Gen Bank; Table S4. Fixation indexes for groups of populations maximally differentiated (K) and inferred by the Spatial Analysis of Molecular Variance (SAMOVA). The inferred groups based on Heermann’s Gull cytochrome b sequences from each location are described by the square brackets (Numbers between brackets correspond to the locations in Figure 1). Supplementary Material S2: Philopatry data. Supplementary Material S3: Logistic model script.

Author Contributions

Conceptualization, E.A.R., A.A., J.G.S.-O., E.V. and M.D.M.-M.; methodology, E.E., E.A.R., J.G.S.-O., A.A., E.V., A.C.-R. and M.D.M.-M.; software, E.E., E.A.R. and M.D.M.-M.; validation, E.E., E.A.R., E.V., A.C.-R. and A.A.; formal analysis, E.E., E.A.R., A.A., E.V. and M.D.M.-M.; investigation, E.A.R., E.V., A.A., J.G.S.-O. and M.D.M.-M.; resources, E.A.R., J.G.S.-O., A.A., E.V., A.C.-R. and J.A.R.-R.; data curation, E.A.R., E.V., A.A. and M.D.M.-M.; writing—original draft preparation, E.A.R., J.G.S.-O. and M.D.M.-M.; writing—review and editing, E.E., E.A.R., E.V., A.A., J.G.S.-O., J.A.R.-R., A.C.-R. and M.D.M.-M.; visualization, E.A.R., J.G.S.-O. and E.V.; supervision, E.A.R., J.G.S.-O. and E.V.; project administration, E.A.R., J.G.S.-O. and E.V.; funding acquisition, E.A.R., J.G.S.-O., A.A. and E.V. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Instituto Politécnico Nacional (IPN), projects SIP-20171835, SIP-20181044, and SIP-20195761 granted to E.A.R., and Universidad de Sonora grant 19-194 to J.G.S.-O. Consejo Nacional de Ciencia y Tecnología (CONACyT) supported the M.D.M.-M. doctoral degree (grant: 630058/577450), who also received support from BEIFI (projects 20200591, and 20210142) and COFAA-IPN (grant DEAE/LA0280/2018). Field work was funded by grants to E.V. between 1985 to 2018 by several donors, must importantly by the join fund from Fondo Mexicano para la Conservación de la Naturaleza/Lindblab Expeditions/National Geographic/Packard Foundation.

Institutional Review Board Statement

This study complies with Mexican regulations regarding the ethical treatment of research subjects. The Dirección General de Vida Silvestre, Subsecretaría de Gestión para la Protección Ambiental, Secretaría de Medio Ambiente y Recursos Naturales, granted the corresponding research and sample collection permits (SGPA/DGVS/02313/10, SGPA/DGVS/02017/11, SGPA/DGVS/01410/12, SGPA/DGVS/01522/12, and SGPA/DGVS/03382/12).

Data Availability Statement

The data presented in this study are available in Supplementary Materials. Also, can be found at the GenBank (https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/, accessed on 25 January 2022) or requested from the corresponding authors.

Acknowledgments

We thank Dirección General de Vida Silvestre from the Secretaría de Medio Ambiente y Recursos Naturales for research permits, Comisión Nacional de Áreas Naturales Protegidas, the staff of Bahía de los Ángeles’ office and José Arce for help and transportation to the islands. José Arce, Rosendo Olachea, Thor Morales, and Jorge Aguilar assisted in seabird data collection in the field. We also thank Gabriela De La Cruz Pino for assistance in band data curation, Benjamín López and Fabián Pérez for advice and assistance during genetic data analyses, and Jazmín García-Román for the illustration of graphical abstract. We are also grateful to the people of San Blas for assistance and logistic support, and to three anonymous reviewers for valuable comments on early versions of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Greenwood, P.J. Mating systems, philopatry and dispersal in birds and mammals. Anim. Behav. 1980, 28, 1140–1162. [Google Scholar] [CrossRef]
  2. Cury, P. Obstinate nature: An ecology of individuals. Thoughts on reproductive behavior and biodiversity. Can. J. Fish Aquat. Sci. 1994, 51, 1664–1673. [Google Scholar] [CrossRef]
  3. Pearce, J.M. Philopatry: A return to origins. Auk 2007, 124, 1085. [Google Scholar] [CrossRef]
  4. Hedgecock, D. Population genetics of marine organisms. Ocean Sci. 2019, 778–783. [Google Scholar] [CrossRef]
  5. Friesen, V.L. Speciation in seabirds: Why are there so many species… and why aren’t there more? J. Ornithol. 2015, 156, 27–39. [Google Scholar] [CrossRef]
  6. Cristofari, R.; Plaza, P.; Fernández, C.E.; Trucchi, E.; Gouin, N.; Le Bohec, C.; Luna-Jorquera, G. Unexpected population fragmentation in an endangered seabird: The case of the Peruvian diving-petrel. Sci. Rep. 2019, 9, 2021. [Google Scholar] [CrossRef]
  7. Dantas, G.P.M.; Oliveira, L.R.; Santos, A.M.; Flores, M.D.; De Melo, D.R.; Simeone, A.; Vianna, J.A. Uncovering population structure in the Humboldt penguin (Spheniscus humboldti) along the Pacific coast at South America. PLoS ONE 2019, 14, e0215293. [Google Scholar] [CrossRef]
  8. Milot, E.; Weimerskirch, H.; Bernatchez, L. The seabird paradox: Dispersal, genetic structure and population dynamics in a highly mobile, but philopatric albatross species. Mol. Ecol. 2008, 17, 1658–1673. [Google Scholar] [CrossRef]
  9. Cristofari, R.; Trucchi, E.; Whittington, J.D.; Vigetta, S.; Gachot-Neveu, H.; Stenseth, N.C.; Le Bohec, C. Spatial heterogeneity as a genetic mixing mechanism in highly philopatric colonial seabirds. PLoS ONE 2015, 10, e0117981. [Google Scholar] [CrossRef] [Green Version]
  10. Coulson, J.C.A. Review of Philopatry in Seabirds and Comparisons with Other Waterbird Species. Waterbirds 2016, 39, 229–240. [Google Scholar] [CrossRef] [Green Version]
  11. Peters, J.L.; Bolender, K.A.; Pearce, J.M. Behavioural vs. molecular sources of conflict between nuclear and mitochondrial DNA: The role of male-biased dispersal in a Holarctic sea duck. Mol. Ecol. 2012, 21, 3562–3575. [Google Scholar] [CrossRef] [PubMed]
  12. Schjørring, S. Ecologically determined natal philopatry within a colony of great cormorants. Behav. Ecol. 2001, 12, 287–294. [Google Scholar] [CrossRef] [Green Version]
  13. Griesser, M.; Halvarsson, P.; Sahlman, T.; Ekman, J. What are the strengths and limitations of direct and indirect assessment of dispersal? Insights from a long-term field study in a group-living bird species. Behav. Ecol. Sociobiol. 2014, 68, 485–497. [Google Scholar] [CrossRef] [Green Version]
  14. Mabry, K.E.; Shelley, E.L.; Davis, K.E.; Blumstein, D.T.; Van Vuren, D.H. Social Mating System and Sex-Biased Dispersal in Mammals and Birds: A Phylogenetic Analysis. PLoS ONE 2013, 8, e57980. [Google Scholar] [CrossRef] [Green Version]
  15. Ribeiro, A.M.; Lloyd, P.; Feldheim, K.A.; Bowie, R.C.K. Microgeographic socio-genetic structure of an African cooperative breeding passerine revealed: Integrating behavioural and genetic data. Mol. Ecol. 2012, 21, 662–672. [Google Scholar] [CrossRef]
  16. Beck, N.R.; Peakall, R.; Heinsohn, R. Social constraint and an absence of sex-biased dispersal drive fine-scale genetic structure in white-winged choughs. Mol. Ecol. 2008, 17, 4346–4358. [Google Scholar] [CrossRef]
  17. Berg, E.C.; Eadie, J.M.; Langen, T.A.; Russell, A.F. Reverse sex-biased philopatry in a cooperative bird: Genetic consequences and a social cause. Mol. Ecol. 2009, 18, 3486–3499. [Google Scholar] [CrossRef]
  18. Lecomte, N.; Gauthier, G.; Giroux, J.F.; Milot, E.; Bernatchez, L. Tug of war between continental gene flow and rearing site philopatry in a migratory bird: The sex-biased dispersal paradigm reconsidered. Mol. Ecol. 2009, 18, 593–602. [Google Scholar] [CrossRef]
  19. Rocha-Olivares, A.; González-Jaramillo, M. Population genetic structure of Mexican Magnificent Frigatebirds: An integrative analysis of the influence of reproductive behavior and sex-biased dispersal. Rev. Mex. Biodivers. 2014, 85, 532–545. [Google Scholar] [CrossRef] [Green Version]
  20. Clucas, G.V.; Younger, J.L.; Kao, D.; Emmerson, L.; Southwell, C.; Wienecke, B.; Hart, T. Comparative population genomics reveals key barriers to dispersal in Southern Ocean penguins. Mol. Ecol. 2018, 27, 4680–4697. [Google Scholar] [CrossRef] [Green Version]
  21. Quillfeldt, P.; Moodley, Y.; Weimerskirch, H.; Cherel, Y.; Delord, K.; Phillips, R.A.; Masello, J.F. Does genetic structure reflect differences in non-breeding movements? A case study in small, highly mobile seabirds. BMC Evol. Biol. 2017, 17, 1–11. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Castillo-Guerrero, J.A.; Piña-Ortiz, A.; Enríquez-Paredes, L.; van der Heiden, A.M.; Hernández-Vázquez, S.; Saavedra-Sotelo, N.C.; Fernández, G. Low genetic structure and diversity of Red-billed Tropicbirds in the Mexican Pacific. J. Field Ornithol. 2020, 91, 142–155. [Google Scholar] [CrossRef]
  23. Perez, G.S.; Goodenough, K.S.; Horn, M.H.; Patton, R.T.; Ruiz, E.A.; Velarde, E.; Aguilar, A. High Connectivity among Breeding Populations of the Elegant Tern (Thalasseus elegans) in Mexico and Southern California Revealed Through Population Genomic Analysis. Waterbirds 2020, 43, 17–27. [Google Scholar] [CrossRef]
  24. Ibarguchi, G.; Gaston, A.J.; Friesen, V.L. Philopatry, morphological divergence, and kin groups: Structuring in thick-billed murres Uria lomvia within a colony in Arctic Canada. J. Avian Biol. 2011, 42, 134–150. [Google Scholar] [CrossRef]
  25. Mancilla-Morales, M.D.; Romero-Fernández, S.; Contreras-Rodríguez, A.; Flores-Martínez, J.J.; Sánchez-Cordero, V.; Herrera, M.G.; López, M.F.; Ruiz, E.A. Diverging Genetic Structure of Coexisting Populations of the Black Storm-Petrel and the Least Storm-Petrel in the Gulf of California. Trop. Conserv. Sci. 2020, 13, 1–12. [Google Scholar] [CrossRef]
  26. Sexton, J.P.; Hangartner, S.B.; Hoffmann, A.A. Genetic isolation by environment or distance: Which pattern of gene flow is most common? Evolution 2014, 68, 1–15. [Google Scholar] [CrossRef] [Green Version]
  27. Pertierra, L.R.; Segovia, N.I.; Noll, D.; Martinez, P.A.; Pliscoff, P.; Barbosa, A.; Vianna, J.A. Cryptic speciation in gentoo penguins is driven by geographic isolation and regional marine conditions: Unforeseen vulnerabilities to global change. Divers. Distrib. 2020, 26, 958–975. [Google Scholar] [CrossRef]
  28. Nunes, G.T.; Mancini, P.L.; Bugoni, L. When Bergmann’s rule fails: Evidences of environmental selection pressures shaping phenotypic diversification in a widespread seabird. Ecography 2017, 40, 365–375. [Google Scholar] [CrossRef]
  29. Wallace, S.J.; Wolf, S.G.; Bradley, R.W.; Harvey, L.L.; Friesen, V.L. The influence of biogeographical barriers on the population genetic structure and gene flow in a coastal Pacific seabird. J. Biogeogr. 2015, 42, 390–400. [Google Scholar] [CrossRef]
  30. Martínez-Abraín, A.; Viedma, C.; Ramón, N.; Oro, D. A note on the potential role of philopatry and conspecific attraction as conservation tools in Audouin’s Gull Larus auduoii. Bird Conserv. Int. 2001, 11, 143–147. [Google Scholar] [CrossRef] [Green Version]
  31. Oro, D.; Tavecchia, G.; Genovart, M. Comparing demographic parameters for philopatric and immigrant individuals in a long-lived bird adapted to unstable habitats. Oecologia 2011, 165, 935–945. [Google Scholar] [CrossRef] [PubMed]
  32. Piper, W.H. Making habitat selection more “familiar”: A review. Behav. Ecol. Sociobiol. 2011, 65, 1329–1351. [Google Scholar] [CrossRef]
  33. Chalfoun, A.D.; Schmidt, K.A. Adaptive breeding-habitat selection: Is it for the birds? Auk 2012, 129, 589–599. [Google Scholar] [CrossRef] [Green Version]
  34. Oro, D.; Pradel, R.; Lebreton, J.D. Food availability and nest predation influence life history traits in Audouin’s gull, Larus audouinii. Oecologia 1999, 118, 438–445. [Google Scholar] [CrossRef] [PubMed]
  35. Wrege, P.H.; Shuford, W.D.; Winkler, D.W.; Jellison, R. Annual variation in numbers of breeding California Gulls at Mono Lake, California: The importance of natal philopatry and local and regional conditions. Condor 2006, 108, 82–96. [Google Scholar] [CrossRef]
  36. Ellis, J.C.; Shulman, M.J.; Jessop, H.; Suomala, R.; Morris, S.R.; Seng, V.; Mach, K. Impact of Raccoons on breeding success in large colonies of Great Black-backed Gulls and Herring Gulls. Waterbirds 2007, 30, 375–383. [Google Scholar] [CrossRef]
  37. Bolton, M.; Conolly, G.; Carroll, M.; Wakefield, E.D.; Caldow, R. A review of the occurrence of inter-colony segregation of seabird foraging areas and the implications for marine environmental impact assessment. Ibis 2019, 161, 241–259. [Google Scholar] [CrossRef] [Green Version]
  38. Antaky, C.C.; Young, L.; Ringma, J.; Price, M.R. Dispersal under the seabird paradox: Probability, life history, or spatial attributes? Mar. Ornithol. 2021, 49, 1–8. Available online: https://www.marineornithology.org/PDF/49_1/49_1_1-8.pdf (accessed on 16 June 2021).
  39. Stacey, P.B.; Ligon, J.D. The benefits-of-philopatry hypothesis for the evolution of cooperative breeding: Variation in territory quality and group size effects. Am. Nat. 1991, 137, 831–846. [Google Scholar] [CrossRef]
  40. Navarro, J.; Votier, S.C.; Aguzzi, J.; Chiesa, J.J.; Forero, M.G.; Phillips, R.A. Ecological Segregation in Space, Time and Trophic Niche of Sympatric Planktivorous Petrels. PLoS ONE 2013, 8, e62897. [Google Scholar] [CrossRef] [Green Version]
  41. O’Connor, E.A.; Hasselquist, D.; Nilsson, J.Å.; Westerdahl, H.; Cornwallis, C.K. Wetter climates select for higher immune gene diversity in resident, but not migratory, songbirds. Proc. R. Soc. B Biol. Sci. 2020, 287, 20192675. [Google Scholar] [CrossRef] [PubMed]
  42. Barbraud, C.; Delord, K. Selection against immigrants in wild seabird populations. Ecol. Lett. 2021, 24, 84–93. [Google Scholar] [CrossRef]
  43. Case, T.J.; Cody, M.L.; Ezcurra, E. A New Island Biogeography of the Sea of Cortés; Oxford University Press: Oxford, UK, 2002; 690p. [Google Scholar]
  44. Velarde, E. Breeding biology of Heermann’s Gulls on Isla Rasa, Gulf of California, Mexico. Auk 1999, 116, 513–519. [Google Scholar] [CrossRef]
  45. Velarde, E.; Ezcurra, E.; Cisneros-Mata, M.A.; Lavín, M.F. Seabird ecology, El Niño anomalies, and prediction of sardine fisheries in the Gulf of California. Ecol. Appl. 2004, 14, 607–615. [Google Scholar] [CrossRef] [Green Version]
  46. Velarde, E.; Ezcurra, E.; Anderson, D.W. Seabird diets provide early warning of sardine fishery declines in the Gulf of California. Sci. Rep. 2013, 3, 1332. [Google Scholar] [CrossRef]
  47. Boswall, J.; Barrett, M. Notes on the breeding birds of Isla Rasa, Baja California. West Birds 1978, 9, 93–108. [Google Scholar]
  48. Mellink, E. History and status of colonies of Heermann’s Gull in Mexico. Waterbirds 2001, 24, 188–194. [Google Scholar] [CrossRef]
  49. Ruiz, E.A.; Velarde, E.; Aguilar, A. Demographic history of Heermann’s Gull (Larus heermanni) from late Quaternary to present: Effects of past climate change in the Gulf of California. Auk 2017, 134, 308–316. [Google Scholar] [CrossRef] [Green Version]
  50. Muller-Karger, F.E.; Rueda-Roa, D.; Chavez, F.P.; Kavanaugh, M.T.; Roffer, M.A. Megaregions among the large marine ecosystems of the Americas. Environ. Dev. 2017, 22, 52–62. [Google Scholar] [CrossRef]
  51. Mercado-Santana, J.A.; Santamaría-del-Ángel, E.; González-Silvera, A.; Sánchez-Velasco, L.; Gracia-Escobar, M.F.; Millán-Núñez, R.; Torres-Navarrete, C. Productivity in the Gulf of California large marine ecosystem. Environ. Dev. 2017, 22, 18–29. [Google Scholar] [CrossRef]
  52. Lluch-Cota, S.E.; Parés-Sierra, A.; Magaña-Rueda, V.O.; Arreguín-Sánchez, F.; Bazzino, G.; Herrera-Cervantes, H.; Lluch-Belda, D. Changing climate in the Gulf of California. Prog. Oceanogr. 2010, 87, 114–126. [Google Scholar] [CrossRef]
  53. Espinosa-Carreón, T.L.; Escobedo-Urías, D. South region of the gulf of California large marine ecosystem upwelling, fluxes of CO2 and nutrients. Environ. Dev. 2017, 22, 42–51. [Google Scholar] [CrossRef]
  54. Velarde, E.; Wilder, B.T.; Felger, R.S.; Ezcurra, E. Floristic diversity and dynamics of Isla Rasa, Gulf of California—A globally important seabird island. Bot. Sci. 2014, 92, 89–101. [Google Scholar] [CrossRef]
  55. Flores-Martínez, J.J.; Herrera, L.G.; Arroyo-Cabrales, J.; Alarcón, I.; Ruiz, E.A. Seasonal dietary differences of the Yellow-footed Gull (Charadriiformes: Laridae) in Isla Partida Norte, Gulf of California, Mexico. Rev. Mex. Biodivers. 2015, 86, 412–418. [Google Scholar] [CrossRef] [Green Version]
  56. Velarde, E. Predation of Nesting Larids by Peregrine Falcons at Rasa Island, Gulf of California, Mexico. Condor 1993, 95, 706–708. [Google Scholar] [CrossRef]
  57. Velarde, E. Predation of Heermann’s Gull (Larus heermanni) Chicks by Yellow-footed Gulls (Larus livens) in Dense and scattered nesting sites. Colon. Waterbird 1992, 15, 8–13. [Google Scholar] [CrossRef]
  58. Lozano, E.; Carmona, R.; Brabata, G. Reproductive success of Heermann’s (Larus heermanni) and Yellow-footed (L. livens) gulls in the south of the Gulf of California, Mexico. Ornitol. Neotrop. 2004, 15, 237–246. Available online: https://sora.unm.edu/sites/default/files/journals/on/v015n02/p0237-p0246.pdf (accessed on 27 October 2020).
  59. Housh, T.B.; Aranda-Gómez, J.J.; Luhr, J.F. Isla Isabel (Nayarit, México): Quaternary alkalic basalts with mantle xenoliths erupted in the mouth of the Gulf of California. J. Volcanol. Geotherm. Res. 2010, 197, 85–107. [Google Scholar] [CrossRef]
  60. Vieyra, L.; Velarde, E.; Ezcurra, E. Effects of Parental Age and Food Availability on the Reproductive Success of Heermann’s Gulls in the Gulf of California. Ecology 2009, 90, 1084–1094. [Google Scholar] [CrossRef]
  61. Dwight, J. The gulls (Laridae) of the world; their plumages, moults, variations, relationships, and distribution. Bull. AMNH 1925, 52, 63–408. Available online: http://digitallibrary.amnh.org/handle/2246/1245 (accessed on 4 May 2021).
  62. R Core Team. R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna, Austria. Available online: https://www.R-project.org (accessed on 13 October 2021).
  63. Gouy, M.; Guindon, S.; Gascuel, O. Sea view version 4: A multiplatform graphical user interface for sequence alignment and phylogenetic tree building. Mol. Biol. Evol. 2010, 27, 221–224. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Darriba, D.; Taboada, G.L.; Doallo, R.; Posada, D. JModelTest 2: More models, new heuristics and parallel computing. Nat. Methods 2012, 9, 772. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Rozas, J.; Ferrer-Mata, A.; Sanchez-DelBarrio, J.C.; Guirao-Rico, S.; Librado, P.; Ramos-Onsins, S.E.; Sanchez-Gracia, A. DnaSP 6: DNA sequence polymorphism analysis of large data sets. Mol. Biol. Evol. 2017, 34, 3299–3302. [Google Scholar] [CrossRef] [PubMed]
  66. Leigh, J.W.; Bryant, D. POPART: Full-feature software for haplotype network construction. Methods Ecol. Evol. 2015, 6, 1110–1116. [Google Scholar] [CrossRef]
  67. Crochet, P.-A.; Desmaris, E. Slow rate of Evolution in the Mitochondrial Region of Gulls (Aves: Laridae). Mol. Biol. Evol. 2000, 17, 1797–1806. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  68. Crochet, P.-A.; Bonhomme, F.; Lebreton, J.-D. Molecular phylogeny and plumage evolution in gulls (Larini). J. Evol. Biol. 2000, 13, 47–57. Available online: http://gull-research.org/papers/16crochet.pdf (accessed on 25 January 2022).
  69. Lefort, V.; Longueville, J.E.; Gascuel, O. SMS: Smart Model Selection in PhyML. Mol. Biol. Evol. 2017, 34, 2422–2424. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  70. Hordijk, W.; Gascuel, O. Improving the efficiency of SPR moves in phylogenetic tree search methods based on maximum likelihood. Bioinformatics 2005, 21, 4338–4347. [Google Scholar] [CrossRef] [Green Version]
  71. Anisimova, M.; Gascuel, O. Approximate likelihood-ratio test for branches: A fast, accurate, and powerful alternative. Syst. Biol. 2006, 55, 539–552. [Google Scholar] [CrossRef]
  72. Guindon, S.; Dufayard, J.F.; Lefort, V.; Anisimova, M.; Hordijk, W.; Gascuel, O. New Algorithms and Methods to Estimate Maximim-Likelihood Phylogenies Assessing the Performance of PhyML 3.0. Syst. Biol. 2010, 59, 307–321. [Google Scholar] [CrossRef] [Green Version]
  73. Pearse, D.E.; Crandall, K.A. Beyond FST: Analysis of population genetic data for conservation. Conserv. Genet. 2004, 5, 585–602. Available online: https://0-link-springer-com.brum.beds.ac.uk/content/pdf/10.1007/s10592-004-1863-z.pdf (accessed on 25 March 2021). [CrossRef]
  74. Ryman, N.; Palm, S.; André, C.; Carvalho, G.R.; Dahlgren, T.G.; Jorde, P.E.; Ruzzante, D.E. Power for detecting genetic divergence: Differences between statistical methods and marker loci. Mol. Ecol. 2006, 15, 2031–2045. [Google Scholar] [CrossRef] [PubMed]
  75. Friesen, V.L.; Burg, T.M.; McCoy, K.D. Mechanisms of population differentiation in seabirds: Invited review. Mol. Ecol. 2007, 16, 1765–1785. [Google Scholar] [CrossRef] [PubMed]
  76. Hudson, R.R.; Boos, D.D.; Kaplan, N.L. A statistical test for detecting geographic subdivision. Mol. Biol. Evol. 1992, 9, 138–151. [Google Scholar] [CrossRef] [Green Version]
  77. Hudson, R.R. A new statistic for detecting genetic differentiation. Genetics 2000, 155, 2011–2014. [Google Scholar] [CrossRef] [PubMed]
  78. Nei, M. Analysis of gene diversity in subdivided populations. Proc. Natl. Acad. Sci. USA 1973, 70, 3321–3323. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  79. Lynch, M.; Crease, T.J. The analysis of population survey data on DNA sequence variation. Mol. Biol. Evol. 1990, 7, 377–394. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  80. Wright, S. The genetical structure of populations. Ann. Eugen. 1951, 15, 323–354. [Google Scholar] [CrossRef] [PubMed]
  81. Hudson, R.R.; Slatkin, M.; Madisson, W.P. Estimation of levels of gene flow from DNA sequence data. Genetics 1992, 132, 583–589. [Google Scholar] [CrossRef]
  82. Slatkin, M. A comparison of three indirect methods for estimating average levels of gene flow. Evolution 1989, 43, 1349–1368. [Google Scholar] [CrossRef]
  83. Excoffier, L.; Lischer, H.E.L. 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–567. [Google Scholar] [CrossRef]
  84. Dupanloup, I.; Schneider, S.; Excoffier, L. A simulated annealing approach to define the genetic structure of populations. Mol. Ecol. 2002, 11, 2571–2581. [Google Scholar] [CrossRef] [PubMed]
  85. Manel, S.; Gaggiotti, O.E.; Waples, R.S. Assignment methods: Matching biological questions with appropriate techniques. Trends Ecol. Evol. 2005, 20, 136–142. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  86. Corander, J.; Cheng, L.; Marttinen, P.; Tang, J. BAPS: Bayesian Analysis of Population Structure; Manual v 6.0; Department of Mathematics and Statistics, University of Helsinki: Helsinki, Finland, 2013; (Latest Update of This Manual on 14 February 2013). [Google Scholar]
  87. Beerli, P.; Palczewski, M. Unified framework to evaluate panmixia and migration direction among multiple sampling locations. Genetics 2010, 185, 313–326. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  88. Mantel, N. The detection of disease clustering and a generalized regression approach. Cancer Res. 1967, 27, 209–220. [Google Scholar]
  89. Rousset, F. Genetic differentiation and estimation of gene flow from F-statistics under isolation by distance. Genetics 1997, 145, 1219–1228. [Google Scholar] [CrossRef]
  90. Weir, B.S.; Cockerman, C.C. Estimating F-statistics for the analysis of population structure. Evolution 1984, 38, 1358–1370. [Google Scholar] [CrossRef]
  91. Lento, G.M.; Haddon, M.; Chambers, G.K.; Baker, C.S. Genetic Variation of Southern Hemisphere Fur Seals (Arctocephalus spp.): Investigation of population structure and species identity. J. Hered. 1997, 88, 202–208. [Google Scholar] [CrossRef] [Green Version]
  92. Hajizadeh, M.; Sokhandan-Bashir, N. Population genetic analysis of potato virus X based on the CP gene sequence. VirusDisease 2017, 28, 93–101. [Google Scholar] [CrossRef] [Green Version]
  93. Wang, J. Application of the one-migrant-per-generation rule to conservation and management. Conserv. Biol. 2004, 18, 332–343. [Google Scholar] [CrossRef]
  94. Mills, L.S.; Allendorf, F.W. The One-Migrant-per-Generation Rule in Conservation and Management. Conserv. Biol. 1996, 10, 1509–1518. [Google Scholar] [CrossRef]
  95. Nathan, L.R.; Kanno, Y.; Vokoun, J.C. Population demographics influence genetic responses to fragmentation: A demogenetic assessment of the ‘one migrant per generation’ rule of thumb. Biol. Conserv. 2017, 210, 261–272. [Google Scholar] [CrossRef]
  96. Ortiz-Ramírez, M.F.; Sánchez-González, L.A.; Castellanos-Morales, G.; Ornelas, J.F.; Navarro-Sigüenza, A.G. Concerted Pleistocene dispersal and genetic differentiation in passerine birds from the Tres Marías Archipelago, Mexico. Auk 2018, 135, 716–732. [Google Scholar] [CrossRef]
  97. Müller, W.; Groothuis, T.G.G.; Eising, C.M.; Dijkstra, C. An experimental study on the causes of sex-biased mortality in the black-headed gull—The possible role of testosterone. J. Anim. Ecol. 2005, 74, 735–741. [Google Scholar] [CrossRef] [Green Version]
  98. Kalmbach, E.; Nager, R.G.; Griffiths, R.; Furness, R.W. Increased reproductive effort results in male-biased offspring sex ratio: An experimental study in a species with reversed sexual size dimorphism. Proc. R. Soc. B 2001, 268, 2175–2179. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  99. Merkling, T.; Welcker, J.; Hewison, A.J.M.; Hatch, S.A.; Kitaysky, A.S.; Speakman, J.R.; Blanchard, P. Identifying the selective pressures underlying offspring sex-ratio adjustments: A case study in a wild seabird. Behav. Ecol. 2015, 26, 916–925. [Google Scholar] [CrossRef]
  100. Becker, P.H.; Ezard, T.H.G.; Ludwigs, J.D.; Sauer-Gürth, H.; Wink, M. Population sex ratio shift from fledging to recruitment: Consequences for demography in a philopatric seabird. Oikos 2008, 117, 60–68. [Google Scholar] [CrossRef]
  101. Becker, P.H.; Erdelen, M. Egg Size in Herring Gulls (Larus argentatus) on Mellum Island, North Sea, West Germany: The Influence of Nest Vegetation, Nest Density, and Colony Development. Colon. Waterbirds 1986, 9, 68–80. [Google Scholar] [CrossRef]
  102. Skórka, P.; Martyka, R.; Wójcik, J.D.; Babiarz, T.; Skórka, J. Habitat and nest site selection in the Common Gull Larus canus in southern Poland: Significance of man-made habitats for conservation of an endangered species. Acta Ornithol. 2006, 41, 137–144. [Google Scholar] [CrossRef]
  103. Kramer, J.; Meunier, J. Kin and multilevel selection in social evolution: A never-ending controversy? [version 1; peer review: 4 approved]. F1000Research 2016, 5, 1–13. [Google Scholar] [CrossRef]
  104. Humeau, L.; Le Corre, M.; Reynolds, S.J.; Wearn, C.; Hennicke, J.C.; Russell, J.C.; Gomard, Y.; Magalon, H.; Pinet, P.; Gélin, P.; et al. Genetic structuring among colonies of a pantropical seabird: Implication for subspecies validation and conservation. Ecol. Evol. 2020, 10, 11886–11905. [Google Scholar] [CrossRef]
  105. Thanou, E.; Sponza, S.; Nelson, E.J.; Perry, A.; Wanless, S.; Daunt, F.; Cavers, S. Genetic structure in the European endemic seabird, Phalacrocorax aristotelis, shaped by a complex interaction of historical and contemporary, physical and nonphysical drivers. Mol. Ecol. 2017, 26, 2796–2811. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  106. Lombal, A.J.; O’dwyer, J.E.; Friesen, V.; Woehler, E.J.; Burridge, C.P. Identifying mechanisms of genetic differentiation among populations in vagile species: Historical factors dominate genetic differentiation in seabirds. Biol. Rev. 2020, 95, 625–651. [Google Scholar] [CrossRef] [PubMed]
  107. Sauve, D.; Patirana, A.; Chardine, J.W.; Friesen, V.L. Mitochondrial DNA reveals population genetic structure within Atlantic but not Pacific populations of a holarctic seabird, the Black-legged Kittiwake Rissa tridactyla. Mar. Ornithol. 2019, 2, 199–208. Available online: https://www.marineornithology.org/content/get.cgi?rn=1318 (accessed on 21 January 2021).
Figure 1. Heermann’s Gull banding and monitoring were carried on Rasa Island between 1984 and 2013. Fledglings (blue dotted lines, n = 28,000; Table S1) were banded within their natal valleys (Orange continuous lines: ES, Valle del Estero; IT, Vallecito del Estero; CA, Valle de la Casita; GE, Gran Estación; VG, Valle de los Gallitos; TV, Tapete Verde; CR, Casitas Viejas Arriba; CB, Casitas Viejas Abajo; and VW, Valle de la W). Later, they were recaptured breeding as adults (n = 998). Mitochondrial cytochrome b sequences (1033 bp) were obtained from 296 gulls in nests from twenty locations distributed on Cardonosa, Rasa and Isabel Islands (dark circles, numbered from 1 to 20). Map modified from Mellink [51], showing the classic regions Midriff Islands Region and Mexican Province [54], as well as the Large Marine Ecosystems [53]. Photo credit: Misael Mancilla.
Figure 1. Heermann’s Gull banding and monitoring were carried on Rasa Island between 1984 and 2013. Fledglings (blue dotted lines, n = 28,000; Table S1) were banded within their natal valleys (Orange continuous lines: ES, Valle del Estero; IT, Vallecito del Estero; CA, Valle de la Casita; GE, Gran Estación; VG, Valle de los Gallitos; TV, Tapete Verde; CR, Casitas Viejas Arriba; CB, Casitas Viejas Abajo; and VW, Valle de la W). Later, they were recaptured breeding as adults (n = 998). Mitochondrial cytochrome b sequences (1033 bp) were obtained from 296 gulls in nests from twenty locations distributed on Cardonosa, Rasa and Isabel Islands (dark circles, numbered from 1 to 20). Map modified from Mellink [51], showing the classic regions Midriff Islands Region and Mexican Province [54], as well as the Large Marine Ecosystems [53]. Photo credit: Misael Mancilla.
Diversity 14 00108 g001
Figure 2. Functional response of banded individuals to distance (in meters) to their natal valley for the three GLM binomial models (A) Model I: Binomial ANOVA for recaptured/total recaptured. (B) Model II: Binomial ANOVA for recaptured/total banded (assuming sex ratios as in recaptured). (C) Model III: Binomial ANOVA for recaptured/total banded (assuming sex ratios to be 50:50%). See section Methods for more details.
Figure 2. Functional response of banded individuals to distance (in meters) to their natal valley for the three GLM binomial models (A) Model I: Binomial ANOVA for recaptured/total recaptured. (B) Model II: Binomial ANOVA for recaptured/total banded (assuming sex ratios as in recaptured). (C) Model III: Binomial ANOVA for recaptured/total banded (assuming sex ratios to be 50:50%). See section Methods for more details.
Diversity 14 00108 g002
Figure 3. Minimum spanning network based on 296 mitochondrial cytochrome b sequences of Heermann’s Gull. Circle sizes are proportional to haplotype frequency. Each color indicates the island of origin: Cardonosa (red), Rasa (yellow), and Isabel (green). Haplotypes differ by one single substitution, unless additional substitutions are present (defined by small perpendicular lines).
Figure 3. Minimum spanning network based on 296 mitochondrial cytochrome b sequences of Heermann’s Gull. Circle sizes are proportional to haplotype frequency. Each color indicates the island of origin: Cardonosa (red), Rasa (yellow), and Isabel (green). Haplotypes differ by one single substitution, unless additional substitutions are present (defined by small perpendicular lines).
Diversity 14 00108 g003
Figure 4. Genetic structure, scaled gene flow, and isolation by distance, calculated through 296 cytochrome b sequences of Heermann’s Gull. (a) Two clusters were inferred by the Bayesian approach BAPS (K = 2). An individual’s probability of belonging to one of the clusters is represented by different colors. (b) Estimates of gene flow obtained by MIGRATE. The effective population size (Θ) is shown within each island, while the scaled migration rate (M) is reported next to arrows showing the direction of migration. (c) Correlation plot among geographic (ln km) and genetic (ΦST /1 − ΦST) distances according to Rousset [89]. Association (rM) and probability (p) values of the Mantel test [88], along with the equation and Pearson coefficient (R2), are shown. A positive correlation among pairwise comparisons from Cardonosa (yellow), Rasa (red), and Isabel (green) is observed.
Figure 4. Genetic structure, scaled gene flow, and isolation by distance, calculated through 296 cytochrome b sequences of Heermann’s Gull. (a) Two clusters were inferred by the Bayesian approach BAPS (K = 2). An individual’s probability of belonging to one of the clusters is represented by different colors. (b) Estimates of gene flow obtained by MIGRATE. The effective population size (Θ) is shown within each island, while the scaled migration rate (M) is reported next to arrows showing the direction of migration. (c) Correlation plot among geographic (ln km) and genetic (ΦST /1 − ΦST) distances according to Rousset [89]. Association (rM) and probability (p) values of the Mantel test [88], along with the equation and Pearson coefficient (R2), are shown. A positive correlation among pairwise comparisons from Cardonosa (yellow), Rasa (red), and Isabel (green) is observed.
Diversity 14 00108 g004
Table 1. Heermann’s Gull individuals banded and recaptured on Rasa Island (n = 998). Individuals were banded at “Natal Valleys” and recaptured at either the same valley or at other “Recapture Valleys”: IT, Vallecito del Estero; CA, Valle de la Casita; GE, Gran Estación Central; TV, Tapete Verde; CR, Casitas Viejas Arriba; CB, Casitas Viejas Abajo; and VW, Valle de la W. Valle del Estero and Valle de los Gallitos valleys were not monitored after banding.
Table 1. Heermann’s Gull individuals banded and recaptured on Rasa Island (n = 998). Individuals were banded at “Natal Valleys” and recaptured at either the same valley or at other “Recapture Valleys”: IT, Vallecito del Estero; CA, Valle de la Casita; GE, Gran Estación Central; TV, Tapete Verde; CR, Casitas Viejas Arriba; CB, Casitas Viejas Abajo; and VW, Valle de la W. Valle del Estero and Valle de los Gallitos valleys were not monitored after banding.
Natal Valley Recapture Valley
SexITCAGETVCRCBVWTotal
ESmales71110009
females22000004
ITmales981700000116
females1752102027
CAmales6621102300282
females248210300110
GEmales24991010107
females0426141137
VGmales00100001
females02000103
TVmales10030110042
females1101591027
CRmales51023207157
females2001294229
CBmales00033230166
females1002137225
VWmales1010152634
females0132221222
All natal valleysmales18023410260674328714
females47973233401817284
Total227331134931076145998
Table 2. Analysis of deviance (binomial-error ANOVA) for recovered bands under three different model hypotheses (see Methods section for details).
Table 2. Analysis of deviance (binomial-error ANOVA) for recovered bands under three different model hypotheses (see Methods section for details).
Model 1: Binomial ANOVA for Recaptured/Total Recaptured
dfDeviancer2pSignif.
distance12564.10.8331<10−70***
effort19.80.00320.001717***
natal place893.00.03021.12 × 10−16***
breeding place595.20.03095.32 × 10−19***
natal × breed. pl.47208.40.06773.4 × 10−22***
sex10.10.00000.740144ns
sex × distance143.10.01405.23 × 10−11***
residuals6164.1
null model1253077.9
Model 2: Binomial ANOVA for Recaptured/Total Banded (Assuming Sex Ratios as in Recaptured)
dfDeviancer2pSignif.
distance12436.60.7626<10−70***
effort1140.70.04401.88 × 10−32***
natal place8330.00.10331.67 × 10−66***
breeding place552.80.01653.68 × 10−10***
natal × breed. pl.47141.00.04412.51 × 10−11***
sex11.20.00040.277356ns
sex × distance134.60.01084.13 × 10−09***
residuals6158.2
null model1253195.0
Model 3: Binomial ANOVA for Recaptured/Total Banded (Assuming Sex Ratios to be 50:50%)
dfDeviancer2pSignif.
distance12436.60.7180< 0−70***
effort1140.70.04151.88 × 10−32***
natal place8331.40.09768.54 × 10−67***
breeding place553.00.01563.32 × 10−10***
natal × breed. pl.47141.00.04152.51 × 10−11***
sex1196.40.05791.31 × 10−44***
sex × distance136.40.01071.63 × 10−09***
residuals6158.3
null model1253393.7
*** = p < 0.001. ns = non-significant.
Table 3. Nucleotide polymorphism for cytochrome b sequences (n = 296) of Heermann’s Gulls breeding in Cardonosa, Rasa, and Isabel Island, within the biogeographic regions Midriff Islands Region (MIR) and Mexican Province (MP). Coordinates of each location with the number of males (M = 132) and females (F = 164) sampled, respectively, are shown. Forty-two haplotypes were recovered (H1 to H42). The number of gulls for each haplotype is next, in brackets. Unique haplotypes are shown in bold.
Table 3. Nucleotide polymorphism for cytochrome b sequences (n = 296) of Heermann’s Gulls breeding in Cardonosa, Rasa, and Isabel Island, within the biogeographic regions Midriff Islands Region (MIR) and Mexican Province (MP). Coordinates of each location with the number of males (M = 132) and females (F = 164) sampled, respectively, are shown. Forty-two haplotypes were recovered (H1 to H42). The number of gulls for each haplotype is next, in brackets. Unique haplotypes are shown in bold.
RegionIslandClusterLocation Lat (N)/Long (W)n (M, F)hSkHdπHaplotypes
MIRCardonosaI128°53′16″/113°01′51″20 (8, 12)320.540.510.00053H1(13), H2(6), H3
RasaII228°49′29″/112°58′50″19 (9, 10)791.300.700.00126H1(10), H2(4), H4, H5, H6, H7, H8
328°49′27″/112°58′50″16 (8, 8)881.310.700.00127H1(9), H2, H8, H9, H10, H11, H12, H13
428°49′28″/112°58′55″12 (6, 6)210.490.490.00047H1(8), H2(4)
528°49′30″/112°58′50″18 (8, 10)430.650.580.00063H1(11), H2(5), H8, H14
628°49′26″/112°58′53″11 (5, 6)551.380.780.00134H1(5), H2 (2), H8, H15(2), H16
III728°49′16′″/112°58′37″13 (4, 9)650.770.770.00109H1(3), H2(6), H8, H14, H17, H18
IV828°49′24″/112°58′41″23(12, 11)781.010.650.00098H1(13), H2(5), H19, H20, H21, H22, H23
928°49′29″/112°58′46″8 (4, 4)320.820.680.00080H1(4), H2(3), H7
1028°49′32″/112°58′42″11 (3, 8)451.160.600.00113H1(7), H2(2), H15, H24
V1128°49′34″/112°58′50″19 (7, 12)320.520.490.00050H1(13), H2(5), H25
VI1228°49′27″/112°58′37″12 (4, 8)210.490.490.00047H1(8), H2(4)
1328°49′19″/112°58′35″20 (13, 7)650.880.680.00085H1(10), H2 (6), H8, H26, H27, H28
1428°49′22″/112°58′33″13 (7, 6)430.690.600.00067H1(8), H2 (3), H29, H30
VII1528°49′27″/112°58′31″14 (6, 8)540.990.590.00096H1(9), H2, H7, H8, H15(2)
VIII1628°49′17″/112°58′43″15 (11, 4)320.380.360.00037H1(12), H2(2), H31
1728°49′20″/112°58′46″16 (5, 11)540.860.680.00085H1(8), H2(5), H30, H32, H33
1828°49′20″/112°58′47″11 (1, 10)441.130.600.00019H1(7), H15(2), H26, H34
IX1928°49′17″/112°58′47″15 (7, 8)430.780.660.00076H1(7), H2(6), H8, H35
MPIsabelX2021°50′40″/105°52′49″10 (4, 6)8124.130.930.00401H1(3), H36, H37, H38, H39, H40, H41, H42
All 29642531.010.620.00098-
Sample size (n), number of haplotypes (h), number of segregating sites (S), an average of nucleotide differences (k), haplotype diversity (Hd), and nucleotide diversity (π). Clusters and locations as defined in Figure 1. Gene Bank accession numbers for each sequence of the forty-two haplotypes are listed in Table S3.
Table 4. Genetic differentiation and gene flow estimates of Heermann’s Gull at different geographic levels (regions, islands, clusters, and locations). Genetic differentiation was estimated by haplotype (χ2 and HST), and sequence-based statistics (KST*, Z*, and Snn), whereas the gene flow was calculated based on haplotype data (GST and Nm1), nucleotide data based on diversity (NST and Nm2), and nucleotide data based on global allele frequencies (FST and Nm3).
Table 4. Genetic differentiation and gene flow estimates of Heermann’s Gull at different geographic levels (regions, islands, clusters, and locations). Genetic differentiation was estimated by haplotype (χ2 and HST), and sequence-based statistics (KST*, Z*, and Snn), whereas the gene flow was calculated based on haplotype data (GST and Nm1), nucleotide data based on diversity (NST and Nm2), and nucleotide data based on global allele frequencies (FST and Nm3).
Levelχ2HSTKST*Z*SnnGSTNm1NSTNm2FSTNm3
Regions206 ***0.009 **0.047 ***9.795 **0.96 9 ***0.02817.560.2271.700.2261.71
Islands224 ***0.062 ns0.046 ***9.797 **0.842 ***0.01827.310.2061.930.2051.94
Clusters470 ***0.068 ns0.032 ***9.922 *0.156 *0.01241.770.1193.700.1183.70
Locations815 ns0.005 ns0.040 ***9.80 ns0.064 ***0.003145.240.0617.730.0607.75
* p < 0.05, ** p < 0.01, *** p < 0.001, and ns: non-significant.
Table 5. Exact test for pairwise ΦST comparisons performed (under diagonal), as well as geographic distances (km) among locations (above diagonal).
Table 5. Exact test for pairwise ΦST comparisons performed (under diagonal), as well as geographic distances (km) among locations (above diagonal).
1234567891011121314151617181920
1-8.578.628.488.548.579.088.838.618.608.438.819.069.018.918.978.848.838.891063.42
2−0.02-0.060.180.070.120.440.250.100.250.180.360.470.470.480.370.260.240.301054.96
30.01−0.01-0.140.100.070.450.280.170.300.220.410.500.510.530.360.230.220.251054.92
4−0.06−0.05−0.01-0.150.090.580.420.270.390.210.510.640.650.660.500.370.350.351055.07
5−0.04−0.03−0.02−0.06-0.140.510.310.120.240.120.370.540.530.520.440.320.310.361054.99
60.01−0.030.00−0.02−0.01-0.500.350.230.360.240.450.560.580.600.410.320.260.271054.98
70.090.040.100.040.080.02-0.230.440.510.600.330.130.230.370.130.220.250.301054.45
8−0.03−0.010.00−0.04−0.030.010.08-0.220.270.400.180.230.230.280.220.170.210.311054.69
90.01−0.040.04−0.030.01−0.03−0.050.00-0.140.180.250.450.420.400.400.290.300.371054.90
10−0.02−0.03−0.01−0.04−0.03−0.040.06−0.02−0.01-0.230.210.470.410.340.490.400.420.501054.90
11−0.04−0.020.00−0.06−0.040.010.12−0.030.03−0.02-0.440.630.600.560.560.430.430.481055.08
12−0.06−0.05−0.01−0.09−0.06−0.020.04−0.04−0.03−0.04−0.06-0.270.200.150.370.350.370.451054.68
13−0.03−0.030.00−0.06−0.03−0.010.03−0.01−0.03−0.02−0.02−0.06-0.110.270.240.310.350.411054.46
14−0.03−0.03−0.02−0.04−0.040.000.11−0.030.03−0.03−0.04−0.04−0.02-0.160.320.360.400.481054.49
15−0.01−0.04−0.01−0.04−0.03−0.060.07−0.01−0.02−0.05−0.01−0.04−0.02−0.02-0.450.450.490.571054.57
160.010.01−0.010.030.000.060.22−0.010.160.00−0.010.030.03−0.020.02-0.130.150.181054.57
17−0.03−0.020.02−0.06−0.02−0.010.02−0.01−0.05−0.02−0.02−0.06−0.03−0.03−0.010.05-0.040.131054.70
180.00−0.030.00−0.02−0.01−0.050.08−0.010.01−0.050.00−0.02−0.01−0.01−0.060.030.00-0.091054.72
19−0.03−0.030.01−0.06−0.03−0.030.00−0.01−0.06−0.02−0.01−0.06−0.04−0.01−0.020.07−0.040.00-1054.67
200.320.260.230.260.300.220.300.290.230.210.310.260.300.250.250.280.280.220.29-
Significant ΦST values are shown in bold (p ≤ 0.05). Locations of sampled nesting sites correspond to that in Figure 1 and Table 3: 1. Playa Rocosa (Cardonosa Island), 2. Valle de la Casita, 3. Vallecito del Estero, 4. Valle del Estero Norte, 5. Rocas Casita, 6. Valle del Estero Centro, 7. Valle de la W, 8. Gran Estación Central Sur, 9. Gran Estación Central Noroeste, 10. Gran Estación Central Rocas, 11. Valle Norte, 12. Valle de los Gallitos Norte, 13. Valle de los Gallitos Sur, 14. Valle de los Gallitos Este, 15. Valle del Fin, 16. Casitas Viejas Abajo, 17. Casitas Viejas Arriba, 18. Casitas Viejas Rocas, 19. Tapete Verde, and 20. Costa Fragata (Isabel Island).
Table 6. Analysis of molecular variance (AMOVA) of Heermann’s Gull cytochrome b, including groupings inferred a prioiri (based on regions, islands, and clusters). For each case, it is showed the source of variation, degrees of freedom (d.f.), the sum of squares (S.S.), percentage of variation (%V), and Φ-statistics. Significant values are shown in bold (* p < 0.05, and *** p < 0.0001).
Table 6. Analysis of molecular variance (AMOVA) of Heermann’s Gull cytochrome b, including groupings inferred a prioiri (based on regions, islands, and clusters). For each case, it is showed the source of variation, degrees of freedom (d.f.), the sum of squares (S.S.), percentage of variation (%V), and Φ-statistics. Significant values are shown in bold (* p < 0.05, and *** p < 0.0001).
GroupingSource of Variationd.f.S.S.%VΦ-Statistics
Regionsamong regions (ΦCT)19.1148.570.49 *
among populations within regions (ΦSC)186.72−0.80−0.02
within populations (ΦST)276133.6752.230.48 ***
Islandsamong islands (ΦCT)29.1824.160.24
among populations within islands (ΦSC)176.61−1.02−0.01
within populations (ΦST)276133.576.860.23 ***
Clustersamong clusters (ΦCT)912.166.950.070
among populations within clusters (ΦSC)103.64−1.65−0.018
within populations (ΦST)276133.594.710.053 ***
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mancilla-Morales, M.D.; Velarde, E.; Aguilar, A.; Contreras-Rodríguez, A.; Ezcurra, E.; Rosas-Rodríguez, J.A.; Soñanez-Organis, J.G.; Ruiz, E.A. Strong Philopatry, Isolation by Distance, and Local Habitat Have Promoted Genetic Structure in Heermann’s Gull. Diversity 2022, 14, 108. https://0-doi-org.brum.beds.ac.uk/10.3390/d14020108

AMA Style

Mancilla-Morales MD, Velarde E, Aguilar A, Contreras-Rodríguez A, Ezcurra E, Rosas-Rodríguez JA, Soñanez-Organis JG, Ruiz EA. Strong Philopatry, Isolation by Distance, and Local Habitat Have Promoted Genetic Structure in Heermann’s Gull. Diversity. 2022; 14(2):108. https://0-doi-org.brum.beds.ac.uk/10.3390/d14020108

Chicago/Turabian Style

Mancilla-Morales, Misael Daniel, Enriqueta Velarde, Andres Aguilar, Araceli Contreras-Rodríguez, Exequiel Ezcurra, Jesús A. Rosas-Rodríguez, José G. Soñanez-Organis, and Enrico A. Ruiz. 2022. "Strong Philopatry, Isolation by Distance, and Local Habitat Have Promoted Genetic Structure in Heermann’s Gull" Diversity 14, no. 2: 108. https://0-doi-org.brum.beds.ac.uk/10.3390/d14020108

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop