Long Huang, Lishi Zhang, Dan Li, Rongfei Yan, Weiping Shang, Yunlei Jiang, Shi Li. 2022: Molecular evidence of introgressive hybridization between related species Jankowski's Bunting (Emberiza jankowskii) and Meadow Bunting (Emberiza cioides) (Aves: Passeriformes). Avian Research, 13(1): 100035. DOI: 10.1016/j.avrs.2022.100035
Citation:
Long Huang, Lishi Zhang, Dan Li, Rongfei Yan, Weiping Shang, Yunlei Jiang, Shi Li. 2022: Molecular evidence of introgressive hybridization between related species Jankowski's Bunting (Emberiza jankowskii) and Meadow Bunting (Emberiza cioides) (Aves: Passeriformes). Avian Research, 13(1): 100035. DOI: 10.1016/j.avrs.2022.100035
Long Huang, Lishi Zhang, Dan Li, Rongfei Yan, Weiping Shang, Yunlei Jiang, Shi Li. 2022: Molecular evidence of introgressive hybridization between related species Jankowski's Bunting (Emberiza jankowskii) and Meadow Bunting (Emberiza cioides) (Aves: Passeriformes). Avian Research, 13(1): 100035. DOI: 10.1016/j.avrs.2022.100035
Citation:
Long Huang, Lishi Zhang, Dan Li, Rongfei Yan, Weiping Shang, Yunlei Jiang, Shi Li. 2022: Molecular evidence of introgressive hybridization between related species Jankowski's Bunting (Emberiza jankowskii) and Meadow Bunting (Emberiza cioides) (Aves: Passeriformes). Avian Research, 13(1): 100035. DOI: 10.1016/j.avrs.2022.100035
Molecular evidence of introgressive hybridization between related species Jankowski's Bunting (Emberiza jankowskii) and Meadow Bunting (Emberiza cioides) (Aves: Passeriformes)
Natural hybridization, which often occurs between closely related species exhibiting sympatric or parapatric distributions, is an important source of genetic variation within populations. The closely related Jankowski's Bunting (Emberiza jankowskii) and Meadow Bunting (E. cioides) are similar in morphology and genetics, occupy overlapping niches, and are sympatric in eastern Inner Mongolia. Previous studies have reported trans-species polymorphisms of alleles between the two species, as well as an unexpectedly high genetic diversity of the endangered E. jankowskii. We speculate that introgressive hybridization has occurred between the two species and contributed to the additional unexpected variation to E. jankowskii. We used mitochondrial NADH dehydrogenase subunit 2 (ND2) gene and 15 nuclear microsatellite markers to compare the genetic diversity of E. jankowskii and E. cioides, and inferred the origin of trans-species polymorphisms between the two species by phylogenetic reconstruction and Bayesian cluster analysis. The two species could be clearly distinguished by population cluster analysis. Despite the large number of mutational differences, we still detected sharing of major haplotypes and the presence of hybrids between the two species. Our study confirmed that weak introgressive hybridization has occurred between sympatric E. jankowskii and E. cioides, which may be mediated by female E. cioides individuals, and that interspecific introgression has contributed to the maintenance of high genetic diversity in E. jankowskii. While being wary of the potential negative effects of introgressive hybridization, we suggest that expanding the habitat of E. jankowskii remains the most effective conservation strategy at present.
Birds have developed a suite of behavioural, physiological and even morphological responses to cope with heat stress (e.g., Carroll et al., 2015; Cook et al., 2020; Pessato et al., 2020; Fernandez et al., 2021; Kim et al., 2021; Playà-Montmany et al., 2021), which is particularly important for ground-dwelling species inhabiting extreme environments (Talbot et al., 2017). The absence or presence of tegumentary structures such as hairs or feathers and their features play an important thermoregulatory role in animals (Mota-Rojas et al., 2021). To date, however, the discussion about the purpose of bare facial and neck skin patches as well as of the pigmentation of these and other fleshy parts in birds has been centred on breeding phenology, social signalling and sexual selection (e.g., Van Niekerk, 1985, 2009; Kimball et al., 2008; Iverson and Karubian, 2017; Stuart-Fox et al., 2017), which is deemed to influence the evolution of morphological traits, particularly when sexual dimorphism is overt (Moore, 1990). Nevertheless, functional traits may evolve under different selection pressures, resulting into multifunctionality (Friedman et al., 2019).
Although the involvement of bare head and neck parts in sexual selection has been largely documented in birds, evidence points towards an important role of these structures also in thermoregulation through evaporative cooling and convection (Crowe and Withers, 1979; Stettenheim, 2000; Eduardo et al., 2010; Gerson et al., 2014; Galván et al., 2017). A pioneering study showed that the surface of the naked head and neck is adjusted to dissipate body heat in the American Black Vulture (Coragyps atratus) (Larochelle et al., 1982). Likewise, it is known that the Southern Ground Hornbill (Bucorvus leadbeatri) dissipates heat via the beak and unfeathered skin (Van Vuuren et al., 2020). Buchholz (1996) tested the heat-dissipation theory by insulating the bare heads and necks of Wild Turkeys (Meleagris gallopavo) with a hood, as if they were feathered, observing increased oxygen consumption, higher body temperature and decreased cooling capacity. However, little research has previously been undertaken to determine the role of bare facial colours in thermoregulation (Crowe, 1979; Walsberg, 1982; Kimball et al., 2008). Peculiar morphological features of head and neck may also have a thermoregulatory function. In a recent study, the long debated function of the peculiar helmet (or casque) of the Cassowary (Casuarius casuarius) has been indicated as primarily associated with thermoregulation, with this spongy bone acting as a thermal window (Eastick et al., 2019). It is possible that morphological traits originally developed to serve an antipredatory function evolved under different environmental pressures leading to their involvement also into homeostasis maintenance to cope with fluctuating ambient temperatures. This could be the case of the Helmeted Guineafowl (Numida meleagris), whose sponge-filled helmet was deemed to have originally evolved to protect eyes and skull of birds running away from predators through the scrub savannah (Crowe, 1978; Van Niekerk, unpublished). In some subspecies, it could be argued that this structure is also implied in thermoregulation, but supporting evidence has yet to be produced.
Interestingly, the Helmeted Guineafowl also displays bare throat skin, often hanging like a loose sack, that varies greatly in size and colour across its range. This raises the question about the adaptive meaning of this morphological trait and why different phenotypes evolved across morphological subspecies ranging in different habitats, which possibly points towards the effect of environmental factors. The implication of the throat sack in thermoregulation seems plausible since the neck hosts blood vessels involved in brain cooling (Crowe and Withers, 1979). By implanting thermocouples in the forebrain of Helmeted Guineafowl it was demonstrated that when the ambient temperature increases, their brain temperature also does, ranging between 36 and 42 ℃ in a healthy bird. Furthermore, heat radiation and convection experiments show that this species can lower its brain temperature through convective cooling of head and neck (Withers and Crowe, 1980) by stretching out the bare neck to dissipate heat (Crowe and Withers, 1979) while retracting it in the collar feathers when it is cold. Crowe (1979) studied morphological variation in Helmeted Guineafowl and showed how it relates to environmental variables, but no mention was made of the sack role in regulating the brain temperature. Nevertheless, in light of its polytypicity and distribution over a wide latitudinal range, the Helmeted Guineafowl is an excellent model to explore the alleged adaptive role of phenotypic variation (i.e., morphology and pigmentation) in the sack with a special focus on thermoregulation. A number of additional elements make the Helmeted Guineafowl a species of choice among other available ground-dwelling birds with bare throat skin to address this research topic. While the ground hornbills consist of just a couple of taxa, spurfowls (Pternistis spp.), which include as many as 24 species (Madge and McGowan, 2002), display bright yellow, orange or red skin patches which are not deemed to be related to heat absorption unlike darker colours (Stuart-Fox et al., 2017), being rather involved in sexual signalling (Van Niekerk, 1979, 2021). Moreover, spurfowl bare throats are inflatable (Van Niekerk, 2010) unlike those of Helmeted Guineafowl (Van Niekerk, 1983). Additional features that make this galliform a suitable model species are its sedentariety and high site-fidelity (with home ranges stretching over 1–2 km2) along the entire lifetime, which was inferred by means of mark-recapture experiments (Van Niekerk et al., 2009).
The Helmeted Guineafowl ranges across tropical and subtropical areas of Sub-Saharan Africa with eight extant morphological subspecies: N. m. coronata (southern Africa), N. m. papillosus (South West Africa), N. m. marungensis (Zambia, southern Democratic Republic of Congo, central and eastern Angola), N. m. mitrata (eastern and southeastern Africa), N. m. reichenowi (Tanzania and Kenya), N. m meleagris (Uganda, Sudan, Central African Republic, Chad, Ethiopia and Eritrea), N. m. somaliensis (Somalia) and N. m. galeatus (West-Central and Western Africa) (Fig. 1). One disjunct subspecies (N. m. sabyi) occurred in Morocco (Crowe, 1979) before it went extinct during the 1950s (Madge and McGowan, 2002; Birdfinding.info/helmeted-guineafowl/). For the sake of clarity, some authors assume the range of N. m. mitrata to include Limpopo Province of South Africa (Clancey, 1980; Chittenden et al., 2012), but we opted to rather follow Crowe (1978), who made a thorough investigation based on 1833 museum skins, and other authors (Madge and McGowan, 2002; Little, 2016).
Figure
1.
Distribution of Helmeted Guineafowl subspecies across Africa. The map was compiled using distributions from Crowe (1978). For each taxon, a close-up of the head and an image of the main habitat are provided. The coloured lines delineate the range of each subspecies. The photographs of guineafowl in this figure were obtained either from the Macaulay Library or from private photographers, while those of the habitats were downloaded from various websites upon permission. Credits for all the photographs are listed in the Acknowledgements.
The evolution process leading to present-day guineafowl across African biomes has been addressed in a seminal work based on multidisciplinary evidence pointing to the bare throat as a derived morphological feature (Crowe, 1978). In the early Pliocene, proto-Agelastes guineafowl with bare head and neck radiated into forests, which are now inhabited by species with fully feathered throats. Subsequently, the proto-Acryllium guineafowl moved into desert steppes, while proto-Guterra guineafowl moved into forest edge biotopes and proto-Numida guineafowl to savannahs. Each group speciated in this suite of different habitats under the selective pressures associated with different climatic conditions and landscape features, with several taxa developing specific neck and head structures presumably as adaptations to control brain temperature in arid and hot environments. For the proto-Numida group, this diversification process occurred in allopatry due to forest expansion isolating savannah patches. According to Crowe (1978), the typical morphological features of Numida, including the horny helmet, filoplumes (a likely protection against UV radiation) on the differently-coloured hindneck and the throat sack, are likely physioecological adaptations of former forest birds to cope with thermoregulation in open dry habitats but could also play a role in social interactions. Noteworthy, homeostasis maintenance and social signalling are crucial for the Helmeted Guineafowl, since this species forage over many hours (often up to 9 h per day) under direct sunlight. Since their food items are dry and often with poor nutritional value, guineafowl ingest about 25 g of biomass to obtain enough energy for roosting in trees at night, when temperature may drop to around ℃ in winter (Van Niekerk, unpublished). At the same time, the long time spent on the open ground for feeding exposes these birds to a high predation risk. For this reason, guineafowl form flocks of 15–40 birds that forage in closely knitted formations as an antipredator strategy (Crowe and Elbin, 1987), which shows why social signalling mediated by colourful bare skin patches is so important in this species.
In this study, we combined unpublished 40-year-old field data with material outsourced from a large online digital photographic archive and information from online climatic databases to explore the Helmeted Guineafowl range-wide sack variation and pigmentation in relation to thermoregulation and sexual signalling. For this purpose, we explored the involvement of environmental factors as possible drivers influencing the evolution of structural and colour variation in the Helmeted Guineafowl throat sack across a spatial gradient in Africa. Specifically, our purposes were to: (1) determine if sack size covaries with deep body temperature in the field; (2) correlate this information with sack area and colour across subspecies in response to different environmental variables; and (3) propose a scientifically plausible explanation of its role in thermoregulation.
2.
Materials and methods
2.1
Fieldwork
From 1 March to 1 July 1982, 56 Helmeted Guineafowl (N. m. coronata) were trapped and marked with numbered metal rings in Krugersdorp Game Reserve (KGR: Gauteng, South Africa; 27°42ʹ41.14ʺ E, 26°5ʹ13.11ʺ S; Appendix Table S1; Van Niekerk, 1982). Birds were caught with a funnel walk-in trap (Van Niekerk, 2010), with the observers monitoring them 200 m away using binoculars to see when they entered it. Once in, the guineafowl kept feeding normally as if they were undisturbed. The trap was then approached as unobtrusively and quickly as possible and two people entered it, catching the birds to measure throat sack size and take body weight readings with a Persola spring balance scale (nearest gram). Guineafowl ring number, age (n = 56), mass (n = 28) and sex (n = 37) were recorded. We also measured the deep body temperature by means of a Hanna digital thermometer fitted with a thermocouple (needle probe) and placed about 5 cm deep into the vent. Measurements were recorded once the values stabilised on the digital screen of the thermometer, which took approximately 30 s (n = 15: young birds were excluded; see below) (Appendix Table S1). The surface temperature on the throat sack (Fig. 2A) was also recorded by enveloping it around the tip of the probe (about 10 mm long) with a self-made wooden forceps (to prevent heat convection). Multiple vent or skin temperatures per bird were not recorded to lessen stress (Nilsson et al., 2016). The surface of the throat sack was measured by unwrapping it gently as to assume its typical triangular shape against a white sizeable thin card box (shoebox). The sack outline was drawn with a pencil onto the paper, which allowed to determine its area (mm2) as that of a triangle calculated as half base multiplied by its height once linear measurements to the nearest mm had been taken with a vernier calliper. Only one side of the sack was measured and outlined (Fig. 2A). We also estimated its elasticity (i.e., sack variation across time) by sizing its surface up in adult birds both before (n = 18) and after (n = 13) handling (i.e., fitting individual rings, taking body temperature and sexing), which was an invariable source of stress. Importantly, the measurements pre- and post-handling were taken on different individuals and on different days. For the sake of clarity, the sack cannot be inflated or deflated: it is essentially a 1–2 mm thick flap which is relaxed when the weather is warm but tensed up around the neck when it is cold. At the same time, this structure is not deemed to be affected by gular fluttering, which is produced by pulsating the skin in the lower jaw with the mouth open (cf. Nassar et al., 2001). When Black-backed Jackals (Lupulella mesomelas) or predatory birds were observed near the trap, the guineafowl were not handled and measured as this would almost certainly further increased their stress levels, thus resulting in biased sack size measures. Field data collection took 6–10 min per bird (Appendix Table S1).
Figure
2.
(A) Dimensions of the head of Helmeted Guineafowl measured with a calliper (mm) to determine the surface area of the sack (throat skin) and of the side of the neck. a × b is equal to the area of one side of the neck, while c × d is equal to the area of one side of the sack. Only the sack surface was measured in birds trapped in KGR, while both surfaces were measured directly from the photographs. (B) Relationship between deep body temperature recorded by a thermocouple and sack size in N. m. coronata free-ranging individuals trapped in KGR. (C) Relationship between sack size and body mass in adult N. m. coronata male free-ranging individuals trapped in KGR.
The photographs of 432 Helmeted Guineafowl (some of which captured in the same shot) stored in the Macaulay Library (Cornell Lab of Ornithology: www.macaulaylibrary.org) (Appendix Table S2) were included to show the variation in sack size along a spatial gradient (Fig. 1). The pictorial data of this repository are being used in other bird studies (e.g., Pyle, 2021), even though caution is needed, especially with digital images (Stevens et al., 2007 and references therein). We strived to select photographs shot as evenly as possible throughout the year in order to account for possible seasonal variation.
The inspection of this material also enabled to determine the sex of adult guineafowl (Appendix Table S3). Males have longer bony and uniformly apricot casques, slightly reddish at the base, which are often backward arching (Siegfried, 1966; Van Niekerk, unpublished), while in females casques — of the same colour — are shorter. Conversely, in first year birds of both sexes this morphological feature is markedly different, being fleshy and with reddish and blackish blotches (Van Niekerk, unpublished), with males starting to develop larger casques at the age of about 27 weeks (Siegfried, 1966). Males in breeding pairs often stand upright to spot predators (Van Niekerk, 1979; 2014; Prinsloo et al., 2005) and assume a low or high intensity hunchback posture when they threaten or chase each other (Crowe and Elbin, 1987; Van Niekerk, 1980), while females spend most of the time in a horizontal position.
No photographic material was found for the extinct N. m. sabyi from Morocco. The photos were selected as to have the lateral side of the head portrayed in the two-dimensional space (i.e., height and length) and minimise obliqueness (i.e., depth) around the edges, which could potentially bias measurements. Moreover, pictures of birds in alert posture, stampede mode or running were excluded since sack areas might have been affected by their stress condition (Van Niekerk, unpublished). Likewise, pictures of birds engaging in gular fluttering (for criteria cf. Nassar et al., 2001) were avoided because of the likely negligible yet possible effects on sack size. Photographs of N. m. coronata individuals west of its natural distribution in South Africa (Fig. 1) were avoided since these birds, introduced to the Western Cape over the past 150 years, are deemed to be admixed with domestic guineafowl (Clancey, 1980).
The sack invariably revealed itself as either a triangular or rectangular (most frequently observed in N. m. papillosus) shape (Fig. 2A), and a different approach was followed to obtain the surface depending on each case once linear measurements to the nearest mm had been taken with a vernier calliper. For triangles, the surface was calculated as the product of half base and height, while for rectangles as the product of their lengths and widths. In both cases, the side of the bare neck was calculated as a rectangle. For the sake of clarity, sack sizes were expressed as the percentage of the total bare skin visible in the picture (sack plus neck; only one side of the head). Measurements were taken from pictures of only seven extant subspecies since we could not find enough pertinent material dealing with N. m. somaliensis, for which, however, information on facial colours was obtained based on a couple of pictures (see acknowledgments).
In order to correlate monthly variation in ambient temperatures with monthly variation in sack size across subspecies, we gathered data from online climatic databases of meteorological stations from across their range (climate-data.com) and determined the mean monthly temperatures for each region. For the sake of clarity, the localities at issue included Skukuza and Bela Bela (South Africa) for N. m. coronata, Sossusvlei and Otjiwarongo (Namibia) for N. m. papilossus, Sofala (Mozambique) and Harare (Zimbabwe) for N. m. mitrata, Lusaka (Zambia) for N. m. marungensis, Nairobi and Tsavo West (Kenya) for N. m. reichenowi, Libreville (Gabon), Ziguinchor (Senegal) and Yaounde (Cameroon) for N. m. galeatus, and Addis Ababa (Ethiopia) for N. m. meleagris. Furthermore, we used D' Agostino test to ascertain that the frequency in the number of pictures selected per month and subspecies followed an unskewed distribution, which ensured an even representation of pictorial data across the year.
A wind factor was determined by multiplying the average monthly wind speed (km/h) by the average number of days per month during which wind was recorded (weather-and-climate.com/average). Humidity data were obtained from weather-and-climate.com/average and given as percentage. Annual direct normal irradiation (kWh/m2) (hereafter referred to as just annual solar irradiation) was sourced from World Bank (2019). Vegetation biomass values (Mg/ha) based on cover, canopy height and tree allometry information from Hanan et al. (2020) were implemented along with the precipitation values outsourced at https://commons.wikimedia.org/wiki/File:Africa_Precipitation_Map-sr.svg#file filelinks.
2.3
Statistical methods
The sack size of males and female Helmeted Guineafowl trapped in KGR were compared using Kruskal-Wallis test (SciStatCalc.com). The same statistics was employed to assess possible variation in sack size as measured before and after bird handling. Spearman rank correlation (Social Science Statisticks.com) was used to explore the relation between sack size and deep body temperature of trapped birds.
We fitted the pictorial data on the sack size of five subspecies with a generalized linear model (GLM) with Gamma residual distribution and log as linking function. We included subspecies and sex of the birds as factor predictors in this model, as well as their interaction, to investigate whether significant differences in sack size between males and females could be generalized across seven N. meleagris subspecies (we had no data for N. m. somaliensis).
We also performed Mann-Whitney test comparisons of the sack size between subspecies pairs and calculated a correlation matrix with the ambient variables including annual solar irradiation (kWh/m2), biomass (Mg/ha), humidity (%) and wind factor (km/h). We confirmed the existence of a strong autocorrelation. Therefore, we summarised these variables in a single environmental factor (see results) by means of a principal component analysis (PCA). We employed a generalized mixed model with sack size as response, the environmental principal factor as fixed predictor and subspecies affiliation as random factor.
We also tested whether average sack size per subspecies covaried with monthly average minimum and maximum temperatures. Hence, we used a GLM where subspecies, maximum average temperature (ATmax), and minimum average temperature (ATmin) were set as predictors along with the interactions subspecies × ATmax and subspecies × ATmin to test if these relationships differed across taxa.
3.
Results
3.1
Sack elasticity and development
Sack size differed significantly between birds measured upon arrival at the trap and those handled later (i.e., 6–10 min). The average sack size before and after handling was 255.05 mm2 (n = 18, SD = 77.77, range = 93–400) and 314.23 mm2 (n = 13, SD = 77.12, range = 204–493), respectively; the difference was significant (n = 31, H = 3.92, P = 0.047). Six birds were captured twice and other two four times, each one on different days, revealing a variation in sack size (Appendix Fig. S1). We found a significantly negative correlation between sack size and deep body temperature (n = 15; rs = −0.067, P = 0.005; Fig. 2B, Appendix Table S1). In some adult birds, sex could not be determined, neither was it possible to determine the sack sizes or sex of birds aged 9–16 weeks (Appendix Table S1). Of the 56 Helmeted Guineafowl trapped in KGR, 15 were immature; of these, one — trapped in March — was a 9-week-old chick which, as such, possessed no throat sack, while head and facial parts were completely covered with natal down (cf. Siegfried, 1966), consistent with what is expected in the hatching season (January to March, namely the hot summer months in the southern hemisphere).
We also trapped seven chicks that were 14 to 16-week-old between March and June when, under cooler ambient conditions, the dark facial skins start appearing (cf. Siegfried, 1966), but no sack was present yet. However, this feature — first visible at the age of about six months — was present in some of the seven older juveniles trapped from June to July (the coldest months), and measured on average 186.9 mm2 (juveniles in Appendix Table S1: SD = 51.29, range 96–285.5). The sudden appearance of the throat sack in both sexes coincided with the loss of facial feathers in place of blue bare skin. No statistical difference in sack size emerged between males and females (n = 36, H = 0.01, P = 0.919; Appendix Fig. S2). Adult males weighting about 1150 g displayed the largest sacks, but their size receded as they got older and heavier (Fig. 2C).
3.2
Bare neck and facial features
The nine recognised subspecies fell in two main groups based on throat sack colour. The subspecies from the southern half of Africa, including N. m. coronata, N. m. papillosus, N. m. mitrata and N. m. marungensis, possess blue sacks. On the other hand, those from Central and East Africa, including N. m. reichenowi, N. m. galeatus and N. m. meleagris, possess black sacks. Finally, the disjunct N. m. sabyi in Morocco displayed mostly bluish-white facial parts (Madge and McGowan, 2002, Fig. 1). The representatives of the first group have long thin hanging wattles with red tips as opposed to their counterpart in which these are thick, rounded and either red or wholly blue (see Fig. 1). All the subspecies display blue bare facial colours around the eyes.
3.3
Comparisons of sack size among subspecies
We found significant differences in sack size of N. m. meleagris when compared to N. m. papillosus and N. m. coronata, as well as in N. m. reichenowi when compared to N. m. papilossus and N. m. coronata (Fig. 3; Table 1). Furthermore, when inspecting the photographic material we could not see any throat sack (meaning that it was completely retracted around the neck) in some individuals of N. m. meleagris (29.09%), N. m. reichenowi (12.66%) and N. m. mitrata (1.2%), while it was always visible in N. m. papillosus and N. m. coronata.
Figure
3.
Box plot showing sack size variation in Helmeted Guineafowl subspecies. Boxes and whiskers are used to indicate mean ± standard error and 0.95 confidence interval, respectively.
Table
1.
Summary statistics derived from a two-way GLM test conducted on sack size in male and female Helmeted Guineafowl as inferred from photographic material (see Appendix Table S2 and Table S4 for details).
We determined the sex and the likely related sack size differences in 266 out of 432 photographed birds encompassing 5 subspecies (Appendix Table S2). We found significant sexual dimorphism (SSD) in sack size of N. m. meleagris and N. m. mitrata (Fig. 4; Table 2).
Figure
4.
Comparisons of sack sizes (measured as ratios of the bare skin in the head and neck) of male and female individuals of Helmeted Guineafowl subspecies.
Table
2.
Mann-Whitney comparisons of sack size data between pairs of Helmeted Guineafowl subspecies (see Appendix Table S4 for sample size and standard deviation).
Pair
Mean Rank difference
Z
SE
Critical value
P-value
N. m. coronata–N. m. meleagris
179.488
5.973
30.049
91.292
2.328E-9
N. m. coronata–N. m. mitrata
45.805
1.631
28.085
85.325
0.103
N. m. coronata–N. m. papillosus
50.636
1.845
27.439
83.361
0.065
N. m. coronata–N. m. reichenowi
143.226
5.346
26.789
81.387
8.971E-8
N. m. coronata–N. m. marungensis
33.877
0.763
44.377
134.822
0.445
N. m. coronata–N. m. galeatus
142.851
4.317
33.093
100.539
1.584E-5
N. m. mitrata–N. m. meleagris
-133.683
4.204
31.803
96.62
2.628E-5
N. m. papillosus–N. m. meleagris
-128.852
4.125
31.234
94.89
3.7E-5
N. m. meleagris–N. m. reichenowi
-36.262
1.183
30.664
93.161
0.237
N. m. marungensis–N. m. meleagris
-145.611
3.110
46.819
142.239
0.002
N. m. meleagris–N. m. galeatus
-36.637
1.009
36.302
110.287
0.313
N. m. mitrata–N. m. papillosus
4.831
0.165
29.349
89.163
0.869
N. m. mitrata–N. m. reichenowi
97.421
3.39
28.742
87.321
0.001
N. m. mitrata–N. m. marungensis
-11.928
0.262
45.583
138.484
0.794
N. m. mitrata–N. m. galeatus
97.046
2.797
34.693
105.4
0.005
N. m. papillosus–N. m. reichenowi
92.590
3.294
28.111
85.403
0.001
N. m. papillosus–N. m. marungensis
-16.759
0.371
45.188
137.283
0.711
N. m. papillosus–N. m. galeatus
92.215
2.699
34.172
103.817
0.007
N. m. marungensis–N. m. reichenowi
-109.349
2.441
44.796
136.094
0.015
N. m. reichenowi–N. m. galeatus
-0.375
0.011
33.653
102.238
0.991
N. m. marungensis–N. m. galeatus
108.974
2.232
48.828
148.343
0.026
SE: standard error; Z: Z-scores; PBonferroni: P-value after Bonferroni correction (0.002).
The single principal factor derived from the environmental data (eigenvalue: 3.67) explained 73.36% of the initial variance in environmental variables. PC1 was positively correlated with humidity (0.86), rainfall (0.98) and biomass (0.80), whereas it was negatively correlated with annual solar irradiation (−0.92). Thus, when PC1 increased, humidity, rainfall and biomass did too as opposed to annual solar irradiation, which instead decreased. The general mixed model indicated a significant and negative relationship between sack size and the environment PC1 (β = −0.32 ± 0.04; Chi1,430 = 51.87; P < 0.001). Thus, birds with larger throat sacks are distributed in drier regions of Africa, with less vegetation and with higher annual solar irradiation (Fig. 5).
Figure
5.
Relationship between ambient variables and sack size in Helmeted Guineafowl across Africa.
Concerning the monthly correlation between ATmin-ATmax and average sack size across subspecies, there was a significant effect of the subspecies by temperatures interaction (subspecies*ATmin: F4,44 = 5.52, P = 0.001; subspecies*ATmax = F4,44 = 2.64, P = 0.046). Hence, the three southernmost subspecies (i.e., N. m. mitrata, N. m. papillosus and N. m. coronata; see Fig. 1) showed a significantly positive correlation between sack size and both high and low ambient temperatures. However, the effects were stronger in relation with ATmin (Appendix Fig. S3, Table S3). On the other hand, no correlation was found for ATmin and ATmax in the northernmost black-throated subspecies.
4.
Discussion
The multifunctionality of bare facial and neck skin patches in birds has been addressed in several studies. However, while its role in sexual selection and social signalling has been widely explored, its involvement in thermoregulation is still scarcely investigated. In this context, wide-ranging polytypic species may represent insightful models to address important questions. In this study, we explored sack size variation in Helmeted Guineafowl subspecies encompassing different habitats and monthly varying conditions across Africa combining empirical data collected in the field with photographic material and environmental information stored in online repositories to elucidate the alleged adaptive significance of this trait. To the best of our knowledge, our work is among the first ones to rely on this novel approach which bears huge potential for future comparative studies. Noteworthy, the photographic repository allowed us to access a spatially representative and larger sample than the one we could have obtained based on museum skins, which, as noted by Crowe (1978), suffer from alterations of the original colours of naked body parts.
Our results are in line with the hypothesis that the heat dissipated through the sack of birds in the southernmost third of the African continent likely turns into lowering the temperatures of the supra-laryngeal, oesophageal, pharyngeal and vagus arteries situated just behind the sack around the throat, thus ultimately assisting cooling of brain and body (Crowe and Crowe, 1979). This function was specifically supported by the negative correlation of sack size with deep body temperature. At the same time, the blue-throated subspecies from this African region are known to withdraw the bare skin of the sack into the neck feathers to prevent heat loss from the brain when experiencing low ambient temperatures (Crowe and Withers, 1979). Concordantly, birds engage in sunning soon after coming from their roost trees during cold winter mornings with their bare necks and sacks exposed to the sunrays (Van Niekerk, unpublished), which suggests a possible role of these structures in transferring warm blood to the brain. This behaviour may also occur in the northern black-throated subspecies, in which a large sack is observed during cold months, with the dark pigments of the sack maximising the absorption of solar radiation (cf. Galván et al., 2018).
Immature individuals displayed no bare facial parts at hatching (Siegfried, 1966), with the natal down or feathers on their head and neck protecting them from solar radiation during summer months. Birds of this age are frequently observed to shelter under dense vegetation cover (Van Niekerk, 2010). The appearance of bare facial parts and sack in young birds during the height of winter invariably causes them to retract their neck into the collar feathers in the mornings and late afternoons to minimise the exposition of fleshy areas and retain body heat. On the other hand, adult birds are often observed standing with their necks outstretched and the sack rolled out when basking in the morning, supporting the idea that this behaviour maximises solar radiation absorbance to warm up the brain (Van Niekerk, unpublished). Young birds still have a seemingly limited ability to regulate their body temperature when ambient temperature fluctuates as observed in hatchlings of other species (Pereyra and Morton, 2001; Remeš and Martin, 2002; Nord and Giroud, 2020). Interestingly, the tendency of the throat sack to regress in older and heavier males can be imputed to a slower metabolism and, hence, an energy-preservation strategy similar to what was observed in the secondary sexual characters across other bird species (e.g., Møller and de Lope, 1999; Balbontín et al., 2011; Cornwallis et al., 2014) and in line with the reduction of thermogenic capacity in aging mammals (Florez-Duquet and Mcdonald, 1998; Grosiak et al., 2020). This phenomenon, presumably occurring also in older females, will require further elucidation.
Bare body parts may also work as social signals in birds by facilitating interactions between mates (Iverson and Karubian, 2017). Accordingly, we found significant SSD in the throat sack of N. m. meleagris and N. m. mitrata. However, this morphological trait was not sexually dimorphic in N. m. papillosus and N. m. coronata. In this respect, it should be noted that neither did we find a difference between the average sack size of N. m. coronata males and females trapped and measured in KGR, essentially showing consistency between fieldwork and photographic measurements. Admittedly, we could not test this SSD hypothesis in N. m. galeatus, N. m. marungensis, N. m. somaliensis and N. m. sabyi. The degree of sexual dimorphism in a species (i.e., phenotypic difference between sexes), is widely associated with the strength of sexual selection (Stuart-Fox and Ord, 2004; Fairbairn et al., 2008; Hoops et al., 2022). We have no evidence in this study to conclude that the latter varies across Helmeted Guineafowl subspecies, but it is arguable that different selective pressures may turn into lowering the degree of sexual dimorphism in some taxa.
The significant association among subspecies affiliation, sex and sack size supported the involvement of Helmeted Guineafowl throat sack in sexual signalling. This structure may have more than one function, assisting these birds in thermoregulation and mate choice (e.g., Iverson and Karubian, 2017; Tattersall et al., 2017; Van Jaarsveld et al., 2021). However, the lack of SSD in the other subspecies suggests that in some of the habitats occupied by N. meleagris across Africa the environmental pressures (e.g., climate) may exert a greater influence on sack size evolution than mate choice (Buchholz, 1996; cf. Galván et al., 2017).
The taxon-dependent multifunctionality of the sack is concordant with a multilocus galliform phylogeny postulating that bare facial parts play a role in heat loss but is also a product of sexual selection (Kimball et al., 2008). It appears that the SSD, occurring in some taxa only, developed as a secondary trait as opposed to the thermoregulatory function occurring in all subspecies. This interpretation is fully consistent with the scenario hypothesised by previous authors postulating that the plesiomorphic fleshy structures in galliforms were primarily involved in heat loss, possibly ending up with playing a role in sexual selection only at a later stage (Kimball et al., 2008, 2011).
Several lines of evidence suggest that the sack has a likely thermoregulatory function: (1) field data on N. m. coronata evidenced that its size correlates negatively with deep body temperature, which most likely suggests heat dissipation through a larger surface; (2) a positive correlation was found between sack size and monthly ambient temperature variation in the southern blue-throated subspecies, exposed to a sharper thermal variation across seasons; (3) even if such correlation was not found in the northern black-throated subspecies, the sack was larger during winter, which is likely to maximise heat absorbance; (4) the biggest sacks in the southern blue-throated subspecies were found in open and drier habitats, probably since larger skin surfaces assist evaporative cooling when less shade is available in the surrounding habitat.
Likewise, the pigmentation of throat sack as it manifests in all subspecies seems related to thermoregulation, with heat absorption facilitated by the dark colour of northern taxa. For instance, the interplay of black sack and less blue facial parts observed in N. m. galeatus aligns with Gloger's rule (Delhey, 2017), stating that birds inhabiting densely vegetated and humid habitats display dark colours (see also Crowe, 1979), which is also functional to camouflage. Black patches are thought to increase heat dissipation through evaporative cooling in humid tropical environments (Galván et al., 2018), but given their reduced size in N. m. galeatus, their role in regulating brain temperature must be limited and possibly constrained by other factors. Our contention is therefore that the lighter blue sacks of southern subspecies make them better adapted to sunny and windy open conditions by deflecting solar radiation while promoting evaporative cooling as compared to black throats, absorbing heat but constraining this process in their northern counterpart. The invariant occurrence of light blue colours around the eyes across all subspecies might be an adaptation to deflect the sunlight and prevent interference with vision (Jeren et al., 2019). Indeed, the limited vascularisation of this area would impair thermoregulation and, consequently, an increase of local temperature might be harmful.
The Helmeted Guineafowl dispersed throughout African open savannah during late Pliocene (Crowe, 1978). However, when forest expansion took place during early Pleistocene, some populations adapted to more densely vegetated areas around remaining forest patches (after a second spell of contraction) in the tropics, originating N. m. galeatus, while others settled in less vegetated regions, diverging into N. m. reichenowi and N. m. meleagris. Other subspecies remained in drier open habitats either in southern (N. m. coronata and N. m. papillosus) or northern (N. m. sabyi) Africa. In these environments, the bare facial skin and the throat sack with its thermoregulatory function can assist these birds in accessing and exploiting open areas away from vegetation patches. In other words, the Helmeted Guineafowl can feed in sunlight and remain cool headed. This also happens during winter, when the weather is dry, and food becomes scarce. In the southern third of Africa, temperature during winter afternoons can reach up to 25 ℃, so the sack would also serve to dissipate heat. Indeed, field observations in KGR showed that the sack is rolled out also during times of heat stress in winter or retracted during cold spells in summer after hailstorms (JHvN personal observations).
The Helmeted Guineafowl subspecies had hitherto never been classified into a southern blue-throated and a northern black-throated group (cf. Crowe, 1978), but it is strikingly evident how several homologous morphological characteristics and environmental correlations are associated with each one of these categories. The northern black-throated group is characterised by black facial parts, rounded wattles and smaller throat sacks. Their body temperature balance is seemingly less affected by ambient temperature, which varies more moderately in and around tropical Africa, than that of their counterpart. Conversely, the southern blue-throated subspecies from more open habitats have developed larger sacks to facilitate heat dissipation and the maintenance of homeostasis, with N. m. coronata and N. m. papillosus showing — unlike other subspecies — no evidence of involvement in sexual signalling.
The distinctiveness of northern and southern group suggests that they might have originated from separate ancestral lineages, and that their divergence was essentially influenced by vegetation biomass and other environmental variables among which ambient temperature played an important role (Crowe, 1979; Cooney et al., 2016). The discrete spatial patterns of phenotypic features such as wattle shape and colours could conceivably suggest that they diverged to minimise interbreeding between neighbouring subspecies, as in the case of N. m. galeatus and N. m. meleagris, but supporting genetic evidence has yet to be produced.
After the work by Crowe (1978), addressing the evolution of guineafowl as a whole, the only study attempting to make a reappraisal of the Helmeted Guineafowl intraspecific phylogeny was actually focused on elucidating the origin of domestic guineafowl (Shen et al., 2021). When using the Vulturine Guineafowl (Acryllium vulturinum) and the Crested Guineafowl (Guterra edouardi) as outgroups, the clade of wild eastern African Helmeted Guineafowl (Sudan and Kenya) diverged first, followed by wild conspecifics from Nigeria and, subsequently, the domestic guineafowl stock originated in West Africa between 5500 and 1300 years ago. However, an in-depth genomic study encompassing present-day N. meleagris subspecies as a whole is needed to make conclusive statements on the species diversification process and taxonomy. In this respect, the genome assembly that was recently produced for the domestic N. meleagris line (Vignal et al., 2019) may represent a useful resource, even though a wild individual would work better as reference.
5.
Conclusions
Our results pointed towards a primary role of Helmeted Guineafowl sack morphology and pigmentation in thermoregulation, with homeostasis being putatively maintained by both evaporative cooling and deflection of solar radiation. On the other hand, observational data suggest that throat sack can assist the birds to gain heat during sunning in cold mornings. This process may be facilitated by the dense arteries of the sack and neck feeding blood to the brain and the rest of the body. We suggest that evaporative cooling is maximised by lighter pigments of the subspecies occurring in the southern third of Africa and inhabiting more open areas with higher wind factors. Conversely, the black-throated subspecies further north seems to be less impacted by ambient temperature, which, at those latitudes, is less prone to sharp fluctuations; this is the likely reason why a large sack is not required in these taxa. Our assertion is that the colour combinations, in conjunction with sack size, is functional to the modulation of heat absorption across the bare neck and facial parts, thus optimising the physiological functions of each subspecies in its own specific habitat.
The alleged thermoregulatory role, however, does not exclude an additional implication of sack size and colour with social signalling. Concordantly, sack size was sexually dimorphic in N. m. meleagris and N. m. mitrata, pointing to its implication in sexual selection (Mori et al., 2017). On the other hand, the absence of SSD of the sack sizes in the two most southerly subspecies, N. m. coronata and N. m. papillosus, could be likely ascribed to its overriding thermoregulatory role in regions where solar radiation and draught are more prominent - and homeostasis maintenance more energetically expensive - than elsewhere. Admittedly, the availability of real-time versus average monthly environmental data other than field information for the entire subspecies set would have provided a more convincing evidence supporting our results, which should be corroborated by further empirical research like the one conducted in KGR. Likewise, the comparison of local environmental climate data versus time-energy budget (i.e., individual metabolic rate) would further elucidate the function of the throat sack in guineafowl. Additional studies are needed to deepen our knowledge of the poorly explored field of thermoregulation in avian research.
Authors' contribution
JVN: conceptualisation, data curation, formal analysis, funding acquisition, investigation, writing – original draft, writing – review & editing. RMP: formal analysis, methodology, writing – review & editing. GF: formal analysis, methodology, writing – review & editing, funding acquisition, visualisation. All authors read and approved the final manuscript.
Ethics statement
All operations were conducted with permits issued by the Department of Nature Conservation of the former Transvaal Provincial Authority. The Parks Board of Krugersdorp Municipality (currently Mogale City Parks and Recreation) provided permission to conduct research in KGR. JHvN was approved as a B-grade (no B 717) ringer by SAFRING, Cape Town University at the time of fieldwork.
Declaration of competing interest
The authors declare that they have no competing interests.
Abbott, R., Albach, D., Ansell, S., Arntzen, J.W., Baird, S.J.E., Bierne, N., et al., 2013. Hybridization and speciation. J. Evol. Biol. 26, 229-246. https://onlinelibrary.wiley.com/doi/abs/10.1111/j.1420-9101.2012.02599.x
Ait Belkacem, A., Gast, O., Stuckas, H., Canal, D., LoValvo, M., Giacalone, G., et al., 2016. North African hybrid sparrows (Passer domesticus, P. hispaniolensis) back from oblivion - ecological segregation and asymmetric mitochondrial introgression between parental species. Ecol. Evol. 6, 5190-5206.
Anderson, E.C., 2008. Bayesian inference of species hybrids using multilocus dominant genetic markers. Phil. Trans. R. Soc. B. 363, 2841-2850. https://www.ncbi.nlm.nih.gov/pubmed/18508754
Anderson, E.C., Thompson, E.A., 2002. A model-based method for identifying species hybrids using multilocus genetic data. Genetics 160, 1217-1229
Arantes, L.S., Ferreira, L.C.L., Driller, M., Repinaldo Filho, F.P.M., Mazzoni, C.J., Santos, F.R., 2020. Genomic evidence of recent hybridization between sea turtles at Abrolhos Archipelago and its association to low reproductive output. Sci. Rep. 10, 12847
Araya-Anchetta, A., Scoles, G.A., Giles, J., Busch, J.D., Wagner, D.M., 2013. Hybridization in natural sympatric populations of Dermacentor ticks in northwestern North America. Ecol. Evol. 3, 714-724
Ardren, W.R., Borer, S., Thrower, F., Joyce, J.E., Kapuscinski, A.R., 1999. Inheritance of 12 microsatellite loci in Oncorhynchus mykiss. J. Hered. 90, 529-536
Arnold, M., 1997. Natural Hybridization and Evolution. Oxford University Press, Oxford
Bandelt, H.J., Forster, P., Rohl, A., 1999. Median-joining networks for inferring intraspecific phylogenies. Mol. Biol. Evol. 16, 37-48
Barbujani, G., Pilastro, A., De Domenico, S., Renfrew, C., 1994. Genetic variation in North Africa and Eurasia: neolithic demic diffusion vs. Paleolithic colonisation. Am. J. Phys. Anthropol. 95, 137-154
Barton, N.H., 2001. The role of hybridization in evolution. Mol. Ecol. 10, 551-568
Beatty, G.E., Philipp, M., Provan, J., 2010. Unidirectional hybridization at a species’ range boundary: implications for habitat tracking. Divers. Distrib. 16, 1-9.
Benson, J.F., Hostetler, J.A., Onorato, D.P., Johnson, W.E., Roelke, M.E., O'Brien, S.J., et al., 2011. Intentional genetic introgression influences survival of adults and subadults in a small, inbred felid population. J. Anim. Ecol. 80, 958-967
Brower, A.V., 2013. Introgression of wing pattern alleles and speciation via homoploid hybridization in Heliconius butterflies: a review of evidence from the genome. Proc. R. Soc. B. 280, 20122302
Cai, T., Wu, G., Sun, L., Zhang, Y., Peng, Z., Guo, Y., et al., 2021. Biogeography and diversification of Old World buntings (Aves: Emberizidae): radiation in open habitats. J. Avian Biol. 52, 6.
Chan, W.Y., Hoffmann, A.A., van Oppen, M.J.H., 2019. Hybridization as a conservation management tool. Conserv. Lett. 12, e12652.
Chapuis, M.P., Estoup, A., 2007. Microsatellite null alleles and estimation of population differentiation. Mol. Biol. Evol. 24, 621-631
Charles, K.M., Stehlik, I., 2021. Assisted species migration and hybridization to conserve cold-adapted plants under climate change. Conserv. Biol. 35, 559-566
Chesser, R.K., Baker, R.J., 1996. Effective sizes and dynamics of uniparentally and diparentally inherited genes. Genetics 144, 1225-1235
Dabrowski, A., Fraser, R., Confer, J.L., Lovette, I.J., 2005. Geographic variability in mitochondrial introgression among hybridizing populations of Golden-winged (Vermivora chrysoptera) and Blue-winged (V. opinus) Warblers. Conserv. Genet. 6, 843-853.
Delph, L.F., Demuth, J.P., 2016. Haldane's rule: genetic bases and their empirical support. J. Hered. 107, 383-391
Evanno, G., Regnaut, S., Goudet, J., 2005. Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol. Ecol. 14, 2611-2620
Excoffier, L., Lischer, H.E., 2010. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol. Ecol. Resour. 10, 564-567
Frankel, O.H., Soule, M.E., 1981. Conservation and Evolution. Cambridge University Press, Cambridge
Frankowski, J., Lubke, K., Coke, M., Weyl, O.L.F., 2020. Genetic variability and demographic history of Anguilla mossambica (Peters, 1852) from continental Africa and Madagascar. J. Fish Biol. 96, 1251-1259. https://www.ncbi.nlm.nih.gov/pubmed/31777080
Gainsford, A., Jones, G.P., Hobbs, J-P.A., Heindler, F.M., van Herwerden, L., 2020. Species integrity, introgression, and genetic variation across a coral reef fish hybrid zone. Ecol. Evol. 10, 11998-12014. https://onlinelibrary.wiley.com/doi/abs/10.1002/ece3.6769
Goudet, J., 1995. FSTAT (Version 1.2): a computer program to calculate F-statistics. J. Hered. 86, 485-486
Grant, P.R., Grant, B.R., 2019. Hybridization increases population variation during adaptive radiation. Proc. Natl. Acad. Sci. U.S.A. 116, 23216-23224
Grant, P.R., Grant, B.R., 2014. Synergism of natural selection and introgression in the origin of a new species. Am. Nat. 183, 671-681
Hall, T.A., 1999. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp. Ser. 41, 95-98
Han, Z., Zhang, L-S., Qin, B.O., Wang, L.I.N., Liu, Y.U., Fu, V.W.K., et al., 2017. Updated breeding distribution and population status of Jankowski’s Bunting Emberiza jankowskii in China. Bird Conserv. Int. 28, 643-652
Harrisson, K.A., Pavlova, A., Goncalves da Silva, A., Rose, R., Bull, J.K., Lancaster, M.L., et al., 2016. Scope for genetic rescue of an endangered subspecies though re-establishing natural gene flow with another subspecies. Mol. Ecol. 25, 1242-1258
Hoban, S.M., McCleary, T.S., Schlarbaum, S.E., Romero-Severson, J., 2009. Geographically extensive hybridization between the forest trees American butternut and Japanese walnut. Biol. Lett. 5, 324-327. https://pubmed.ncbi.nlm.nih.gov/19324631
Hoelzer, G.A., 1997. Inferring phylogenies from mtDNA variation: mitochondrial-gene trees versus nuclear-gene trees revisited. Evolution 51, 622-626
Irwin, D.E., Rubtsov, A.S., Panov, E.N., 2009. Mitochondrial introgression and replacement between yellowhammers (Emberiza citrinella) and pine buntings (Emberiza leucocephalos) (Aves: Passeriformes). Biol. J. Linn. Soc. 98, 422-438
Jacquemyn, H., Brys, R., Honnay, O., Roldan-Ruiz, I., 2012. Asymmetric gene introgression in two closely related Orchis species: evidence from morphometric and genetic analyses. BMC Evol. Biol. 12, 178
Jangjoo, M., Matter, S.F., Roland, J., Keyghobadi, N., 2016. Connectivity rescues genetic diversity after a demographic bottleneck in a butterfly population network. Proc. Natl. Acad. Sci. U.S.A. 113, 10914-10919
Kutschera, V.E., Bidon, T., Hailer, F., Rodi, J.L., Fain, S.R., Janke, A., 2014. Bears in a forest of gene trees: phylogenetic inference is complicated by incomplete lineage sorting and gene flow. Mol. Biol. Evol. 31, 2004-2017
Lade, J.A., Murray, N.D., Marks, C.A., Robinson, N.A., 1996. Microsatellite differentiation between Phillip Island and mainland Australian populations of the red fox Vulpes vulpes. Mol. Ecol. 5, 81-87
Langin, K.M., Sillett, T.S., Funk, W.C., Morrison, S.A., Ghalambor, C.K., 2017. Partial support for the central-marginal hypothesis within a population: reduced genetic diversity but not increased differentiation at the range edge of an island endemic bird. Heredity 119, 8-15.
Lepais, O., Petit, R.J., Guichoux, E., Lavabre, J.E., Alberto, F., Kremer, A., et al., 2009. Species relative abundance and direction of introgression in oaks. Mol. Ecol. 18, 2228-2242
Lexer, C., Fay, M.F., Joseph, J.A., Nica, M.S., Heinze, B., 2005. Barrier to gene flow between two ecologically divergent Populus species, P. alba (white poplar) and P. tremula (European aspen): the role of ecology and life history in gene introgression. Mol. Ecol. 14, 1045-1057
Li, D., Sun, K., Zhao, Y., Lin, A., Li, S., Jiang, Y., et al., 2017. Polymorphism in the major histocompatibility complex (MHC class II B) genes of the Rufous-backed Bunting (Emberiza jankowskii). PeerJ 5, e2917
Li, S., Li, D., Zhang, L.S., Shang, W.P., Qin, B., Jiang, Y.L., 2021. High levels of genetic diversity and an absence of genetic structure among breeding populations of the endangered Rufous-backed Bunting in Inner Mongolia, China: implications for conservation. Avian Res. 12, 7
Liang, G., Li, T., Yin, Z.H., Lei, F.M., 2008. Molecular phylogenetic analysis of some Fringillidae species based on mitochondrial coi gene sequences. Zool. Res. 29, 465-475
Lima-Rezende, C.A., Rocha, A.V., Junior, A.F.C., Martins, E.S., Vasconcelos, V., Caparroz, R., 2019. Late Pleistocene climatic changes promoted demographic expansion and population reconnection of a Neotropical savanna-adapted bird, Neothraupis fasciata (Aves: Thraupidae). PLoS ONE 14, e0212876
Mallet, J., 2005. Hybridization as an invasion of the genome. Trends Ecol. Evol. 20, 229-237. https://www.ncbi.nlm.nih.gov/pubmed/16701374
McCracken, K.G., Wilson, R.E., 2011. Gene flow and hybridization between numerically imbalanced populations of two duck species in the Falkland Islands. PLoS ONE 6, e23173
Murphy, C., Burnett, S., Conroy, G.C., Howland, B.W.A., Lamont, R.W., Sumner, J., et al., 2019. Genetic diversity and structure of the threatened striped legless lizard, Delma impar: management implications for the species and a translocated population. Conserv. Genet. 20, 245-257.
Newhouse, D.J., Balakrishnan, C.N., 2015. High major histocompatibility complex class I polymorphism despite bottlenecks in wild and domesticated populations of the zebra finch (Taeniopygia guttata). BMC Evol. Biol. 15, 265
Nolte, A.W., Tautz, D., 2010. Understanding the onset of hybrid speciation. Trends Genet. 26, 54-58
Oosterhout, C.V., Hutchinson, W.F., Wills, D.P.M., Shipley, P., 2010. Micro-checker: software for identifying and correcting genotyping errors in microsatellite data. Mol. Ecol. Resour. 4, 535-538
Packert, M., Ait Belkacem, A., Wolfgramm, H., Gast, O., Canal, D., Giacalone, G., et al., 2019. Genetic admixture despite ecological segregation in a North African sparrow hybrid zone (Aves, Passeriformes, Passer domesticus × Passer hispaniolensis). Ecol. Evol. 9, 12710-12726.
Packert, M., Favre, A., Schnitzler, J., Martens, J., Sun, Y-H., Tietze, D.T., et al., 2020. "Into and Out of" the Qinghai-Tibet Plateau and the Himalayas: centers of origin and diversification across five clades of Eurasian montane and alpine passerine birds. Ecol. Evol. 10, 9283-9300. https://pubmed.ncbi.nlm.nih.gov/32953061
Packert, M., Sun, Y-H., Strutzenberger, P., Valchuk, O., Tietze, T., Martens, J., 2015. Phylogenetic relationships of endemic bunting species (Aves, Passeriformes, Emberizidae, Emberiza koslowi ) from the eastern Qinghai-Tibet Plateau. Vertebr. Zool. 65, 135-150
Pamilo, P., Nei, M., 1988. Relationships between gene trees and species trees. Mol. Biol. Evol. 5, 568-583
Paxton, R.J., Thoren, P.A., Tengo, J., Estoup, A., Pamilo, P., 1996. Mating structure and nestmate relatedness in a communal bee, Andrena jacobi (Hymenoptera, Andrenidae), using microsatellites. Mol. Ecol. 5, 511-519
Peakall, R., Smouse, P.E., 2012. GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and research--an update. Bioinformatics 28, 2537-2539
Peng, X., Wu, X., Zhang, Y., Zhang, H., Tang, J., Kang, B., 2019. Genetic diversity and metapopulation structures of two intertidal species along the coast of Zhejiang, China: implications for conservation. Mitochondrial DNA A DNA Mapp. Seq. Anal. 30, 674-681
Peters, K.J., Myers, S.A., Dudaniec, R.Y., O'Connor, J.A., Kleindorfer, S., 2017. Females drive asymmetrical introgression from rare to common species in Darwin's tree finches. J. Evol. Biol. 30, 1940-1952
Pierce, A.A., Gutierrez, R., Rice, A.M., Pfennig, K.S., 2017. Genetic variation during range expansion: effects of habitat novelty and hybridization. Proc. R. Soc. B. 284, 20170007
Pohjoismaki, J.L.O., Michell, C., Levanen, R., Smith, S., 2021. Hybridization with mountain hares increases the functional allelic repertoire in brown hares. Sci. Rep. 11, 15771
Price, T.D., Bouvier, M.M., 2002. The evolution of F1 postzygotic incompatibilities in birds. Evolution 56, 2083-2089
Pritchard, J.K., Stephens, M., Donnelly, P., 2000. Inference of population structure using multilocus genotype data. Genetics 155, 945-959
Qu, Y., Zhang, R., Quan, Q., Song, G., Li, S.H., Lei, F., 2012. Incomplete lineage sorting or secondary admixture: disentangling historical divergence from recent gene flow in the Vinous-throated parrotbill (Paradoxornis webbianus). Mol. Ecol. 21, 6117-6133
Quilodran, C.S., Currat, M., Montoya-Burgos, J.I., 2018. Effect of hybridization with genome exclusion on extinction risk. Conserv. Biol. 32, 1139-1149
Ravinet, M., Elgvin, T.O., Trier, C., Aliabadian, M., Gavrilov, A., Saetre, G.P., 2018. Signatures of human-commensalism in the house sparrow genome. Proc. R. Soc. B. 285, 20181246
Raymond, M., Rousset, F., 1995. GENEPOP (Version 1.2): population genetics software for exact tests and ecumenicism. J. Hered. 86, 248-249.
Ronka, N., Pakanen, V-M., Blomqvist, D., Degtyaryev, V., Golovatin, M., Isakov, G., et al., 2019. Near panmixia at the distribution-wide scale but evidence of genetic differentiation in a geographically isolated population of the Terek Sandpiper Xenus cinereus. Ibis 161, 632-647.
Rousset, F., 1997. Genetic differentiation and estimation of gene flow from F-statistics under isolation by distance. Genetics 145, 1219-1228
Rozas, J., Ferrer-Mata, A., Sanchez-DelBarrio, J.C., Guirao-Rico, S., Librado, P., Ramos-Onsins, S.E., et al., 2017. DnaSP 6: DNA sequence polymorphism analysis of large data sets. Mol. Biol. Evol. 34, 3299-3302
Toews, D.P., Brelsford, A., 2012. The biogeography of mitochondrial and nuclear discordance in animals. Mol. Ecol. 21, 3907-3930. https://www.ncbi.nlm.nih.gov/pubmed/22738314
Tompkins, D.M., Mitchell, R.A., Bryant, D.M., 2006. Hybridization increases measures of innate and cell-mediated immunity in an endangered bird species. J. Anim. Ecol. 75, 559-564.
Vaha, J.P., Primmer, C.R., 2006. Efficiency of model-based Bayesian methods for detecting hybrid individuals under different hybridization scenarios and with different numbers of loci. Mol. Ecol. 15, 63-72
Wang, H.T., Jiang, Y.L., Gao, W., 2010. Jankowski’s Bunting (Emberiza jankowskii): current status and conservation. Chinese Birds 1, 251-258
Wang, W.J., Dai, C.Y., Alstrom, P., Zhang, C.L., Qu, Y.H., Li, S.H., et al., 2014. Past hybridization between two East Asian long-tailed tits (Aegithalos bonvaloti and A. fuliginosus). Front. Zool. 11, 40
Wesselmann, M., Gonzalez-Wanguemert, M., Serrao, E.A., Engelen, A.H., Renault, L., Garcia-March, J.R., et al., 2018. Genetic and oceanographic tools reveal high population connectivity and diversity in the endangered pen shell Pinna nobilis. Sci. Rep. 8, 4770. https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5859023/
Wolfgramm, H., Martens, J., Topfer, T., Vamberger, M., Pathak, A., Stuckas, H., et al., 2021. Asymmetric allelic introgression across a hybrid zone of the coal tit (Periparus ater) in the central Himalayas. Ecol. Evol. 11, 17332-17351.
Yoshida, K., Miyagi, R., Mori, S., Takahashi, A., Makino, T., Toyoda, A., et al., 2016. Whole-genome sequencing reveals small genomic regions of introgression in an introduced crater lake population of threespine stickleback. Ecol. Evol. 6, 2190-2204.
Zhai, D.D., Li, W.J., Liu, H.Z., Cao, W.X., Gao, X., 2019. Genetic diversity and temporal changes of an endemic cyprinid fish species, Ancherythroculter nigrocauda, from the upper reaches of Yangtze River. Zool. Res. 40, 427-438. https://www.ncbi.nlm.nih.gov/pubmed/31111694
Zheng, G., 2005. A Checklist on the Classification and Distribution of the Birds of China. Science Press, Beijing
Zink, R.M., 2005. Natural selection on mitochondrial DNA in Parus and its relevance for phylogeographic studies. Proc. R. Soc. B. 272, 71-78
Table
1.
Summary statistics derived from a two-way GLM test conducted on sack size in male and female Helmeted Guineafowl as inferred from photographic material (see Appendix Table S2 and Table S4 for details).
Table
2.
Mann-Whitney comparisons of sack size data between pairs of Helmeted Guineafowl subspecies (see Appendix Table S4 for sample size and standard deviation).
Pair
Mean Rank difference
Z
SE
Critical value
P-value
N. m. coronata–N. m. meleagris
179.488
5.973
30.049
91.292
2.328E-9
N. m. coronata–N. m. mitrata
45.805
1.631
28.085
85.325
0.103
N. m. coronata–N. m. papillosus
50.636
1.845
27.439
83.361
0.065
N. m. coronata–N. m. reichenowi
143.226
5.346
26.789
81.387
8.971E-8
N. m. coronata–N. m. marungensis
33.877
0.763
44.377
134.822
0.445
N. m. coronata–N. m. galeatus
142.851
4.317
33.093
100.539
1.584E-5
N. m. mitrata–N. m. meleagris
-133.683
4.204
31.803
96.62
2.628E-5
N. m. papillosus–N. m. meleagris
-128.852
4.125
31.234
94.89
3.7E-5
N. m. meleagris–N. m. reichenowi
-36.262
1.183
30.664
93.161
0.237
N. m. marungensis–N. m. meleagris
-145.611
3.110
46.819
142.239
0.002
N. m. meleagris–N. m. galeatus
-36.637
1.009
36.302
110.287
0.313
N. m. mitrata–N. m. papillosus
4.831
0.165
29.349
89.163
0.869
N. m. mitrata–N. m. reichenowi
97.421
3.39
28.742
87.321
0.001
N. m. mitrata–N. m. marungensis
-11.928
0.262
45.583
138.484
0.794
N. m. mitrata–N. m. galeatus
97.046
2.797
34.693
105.4
0.005
N. m. papillosus–N. m. reichenowi
92.590
3.294
28.111
85.403
0.001
N. m. papillosus–N. m. marungensis
-16.759
0.371
45.188
137.283
0.711
N. m. papillosus–N. m. galeatus
92.215
2.699
34.172
103.817
0.007
N. m. marungensis–N. m. reichenowi
-109.349
2.441
44.796
136.094
0.015
N. m. reichenowi–N. m. galeatus
-0.375
0.011
33.653
102.238
0.991
N. m. marungensis–N. m. galeatus
108.974
2.232
48.828
148.343
0.026
SE: standard error; Z: Z-scores; PBonferroni: P-value after Bonferroni correction (0.002).
DownLoad:
Email Alerts
RSS Feeds