Lu Wang, Fangqing Liu, Lan Zhao, Yaling Xu, Tianyu Zhang, Longying Wen. 2023: A test of genetic divergence of a bird existing in the Sichuan Basin and its surrounding mountain ranges. Avian Research, 14(1): 100144. DOI: 10.1016/j.avrs.2023.100144
Citation: Lu Wang, Fangqing Liu, Lan Zhao, Yaling Xu, Tianyu Zhang, Longying Wen. 2023: A test of genetic divergence of a bird existing in the Sichuan Basin and its surrounding mountain ranges. Avian Research, 14(1): 100144. DOI: 10.1016/j.avrs.2023.100144

A test of genetic divergence of a bird existing in the Sichuan Basin and its surrounding mountain ranges

Funds: 

the National Science Foundation of China 31372171

Leshan Normal University research grants 205210094

Leshan Normal University research grants ZZ201805

Leshan Normal University research grants CGZZ202002

Leshan Normal University research grants 205220114

Leshan Normal University research grants DGZZ202006

More Information
  • Corresponding author:

    E-mail address: lywen02@126.com (L. Wen)

  • Received Date: 26 Mar 2023
  • Rev Recd Date: 24 Sep 2023
  • Accepted Date: 06 Oct 2023
  • Available Online: 10 Jan 2024
  • Publish Date: 05 Nov 2023
  • Island ecosystems, serving as natural laboratories, facilitate geographical isolation, ecological specialization, and species divergence. The Sichuan Basin, surrounded by mountain ranges, represents a typical continental island due to its marked environmental spatial heterogeneity. This heterogeneity may contribute to geographical isolation and habitat heterogeneity, resulting in genetic divergence within populations. Therefore, we used the White-browed Laughingthrush (Garrulax sannio) as a model specimen to investigate the genetic divergence in the Sichuan Basin and its surrounding mountain ranges, given its presence in various habitats within and beyond this basin. Employing a RAD-seq dataset of 140 G. sannio individuals from 17 distinct ecological zones in the Sichuan Basin and its surrounding mountain ranges, we conducted PCA, population structure analysis, phylogenetic tree construction, and gene flow analysis to comprehensively analyze G. sannio groups. Additionally, in conjunction with geographical and ecological data, we performed isolation by distance, isolation by environment, PCA, and latent factor mixed model analysis to identify factors influencing the genetic divergence among these G. sannio groups. In summary, the 17 G. sannio groups were categorized into high-elevation, medium-elevation, and low-elevation groups. Genetic divergence in G. sannio may be attributed to both geographical distance and key ecological factors, particularly elevation and key climatic variables. Notably, the high-elevation group exhibited a greater number of SNPs and selected genes associated with the key ecological factors compared to the low-elevation group. The ADCY9 gene and several associated key pathways were identified as crucial elements driving ecological adaptation (elevation and key climatic variables) in the high-elevation group. Furthermore, climate changes during the glacial cycles may have facilitated gene flow among these groups residing in the Sichuan Basin and its surrounding mountain ranges. Our findings provide evidence of genetic divergence in G. sannio influenced by the geographical distance and key ecological factors between the Sichuan Basin and its surrounding mountain ranges. These results lay the groundwork for future research on the molecular systematics of continental islands.

  • In Pakistan, the Grey Francolin (Francolinus pondicerianus) is widely associated with the drier regions of the Indus plains and has penetrated the Thar desert in Sindh, as well as the Thal and Chohlistan deserts in Punjab. The species also occurs in the lower hills of the Makran and Lasbela districts in Balochistan, the Cherat and Kohat districts of Khyber Pakhtunkhwa Province, the salt range and agro-forestry tracks of the Pothwar Plateau in the Punjab and in the Margalla hills of Islamabad (Roberts, 1991; Mahmood et al., 1997; Mann and Chaudhry, 2000). The bird is normally found foraging in open cultivated tracks and grasslands intermixed with scrub forests and is rarely observed above an elevation of 1200 m in Pakistan. The calls of the Grey Francolin, uttered mostly at dawn and at dusk, are considered one of the most familiar bird songs of countryside (Roberts, 1991). In the Indian subcontinent, the Grey Francolin has remained a favorite game bird and continues to be hunted for food since it is considered a cheap source of meat to local inhabitants (Long, 1981). This bird is also used for cockfighting at village fairs (Javed et al., 1999). The species is a friend of farmers who believe this bird consumes large numbers of harmful insects, as well as their eggs and larvae and serves therefore as a biological control agent of insect pests in Pakistan (Beg and Qureshi, 1972; Qureshi, 1972; Mian, 1995).

    Although the bird is intensively hunted and captured in its native range in Pakistan, owing to which local populations could be declining, but the overall status of the species is regarded as stable (del Hoyo et al., 1994; Islam, 1999; Richard et al., 2002). Therefore, it is not believed that its populations are approaching the threshold criteria set for population decline in the IUCN Red List (i.e. declining more than 30% in ten years or in three generations) and thus the conservation status of this species has been rated as of Least Concern (IUCN, 2011). However, the Grey Francolin has a protected status in Punjab Province and falls in Schedule-I of the Punjab Wildlife Acts and Rules 1974. This law permits restricted hunting of the species during the non-breeding winter months, i.e. from 15 November to 15 January, with a limit to hunting not more than six birds per day per person (Shafiq, 2005). Such legal instruments and general public perception regarding current populations of the Grey Francolin have led us to believe that the species is under pressure, presumably due to illegal hunting, urbanization, expansion and intensification of agricultural practices and other pressures on scrub forests by local populations for their use in fodder, timber and fire wood needs.

    Rather than providing a complete picture of the field biology of this species in any particular habitat/area, most of the existing country level studies have reported different and exclusive life history traits (Faruqui et al., 1960; Ullah, 1991; Chaudhry and Bhatti, 1992; Mian and Wajid, 1994; Mahmood et al., 1997; Kamal, 2000; Mann and Chaudhry, 2000). However, a recent investigation by Khan (2010) has covered breeding and dietary aspects of a population inhabiting a protected area of desert ecology, about 500 km towards the southeast of the Pothwar Plateau. The Pothwar Plateau is one of the major distribution ranges of the Grey Francolin (Roberts, 1991). The plateau is characterized by broken terrain, with a total area of 1.82 million ha, out of which only 0.61 million ha is under cultivation and the remainder is occupied by scrub forests (Ahmad, 1990). The available literature provides limited information on the biology of the species in the Pothwar range of its distribution. Hence, in 2009, we explored some factors related to population, habitat use, diet composition and breeding potential of this elegant galliform in a representative site of the agro-ecosystem of the Pothwar Plateau.

    Population and breeding studies were conducted in the area of Jhungal village (33°26.754 N, 73°07.110 E), in the Gujar Khan district of Rawalpindi. The land utilization of Gujar Khan consists for 46% of cultivation and 26% is covered by natural vegetation and forest tracks. Its main crops are wheat, maize, millet, mustard and groundnut (Helders, 2007). The annual rain fall of the area ranges between 500 and 1000 mm. Major tree species, composed of wild flora, are Phulai (Acacia modesta), Kikar (Acacia nilotica), Toon (Cedrala toona), Dharek (Melia azederach), Ber (Zizyphus mausitiana) and Shisham (Dalbergia sissoo). A conspicuous species of grass called Dub (Desmostachia bipinnata) has been reported in abundance (Nizami et al., 2004). The landscape of the study area is represented by the typical characteristics of Pothwar's agro-ecology, consisting of undulated open land used for rain-fed agriculture. These tracks are inter-dispersed with gullies and low hills/ranges covered by dry sub-tropical forests. Generally, winter cultivation (November to April) consists of wheat with intercropping of grams (chickpeas), lentils and mustards whereas the summer (May to October) crops include sorghum, millet and groundnuts or peanuts (Beg et al., 1985; Ahmad, 1990).

    Estimation of the size of the Grey Francolin population was carried out using a line-transect method (Burnham et al., 1980). The study area (9 km2) included two representative sampling sites, i.e., cropland and scrub forest habitats. In both habitats three isolated permanent line transects, with a fixed width and length (500 m × 50 m), were established to monitor the direct and indirect evidence of the presence of the francolin species (Fig. 1). For direct sighting in each of the sampling sites, the observer walked at a moderate speed (2–3 km·hour–1) taking normal paces to record birds within the transect width. The observations were recorded once a month in the morning, the peak activity hours of the bird (Gould, 1966; Mahmood et al., 2010). Only direct sighting data were used for population assessment. Density (number of birds per ha) was calculated by dividing the number of sighted birds by the transect area (2.5 ha = 500 m × 50 m, length × width of the transect). Indirect evidence such as nests, droppings, foot prints, feathers and vocal calls were used to search the bird nests for an in-situ breeding study and to identify its roosting and feeding areas. Similarly, the vegetation associated with nesting sites were also observed and identified.

    Figure  1.  Maps of study site and transect lines established for monitoring Grey Francolin (Francolinus pondicerianus)populations in agro-ecosystems in the Rawalpindi district, Pakistan

    Assuming equal potential of both habitats for hosting and breeding activities of the Grey Francolin, investigations of the study location for nests were made during March, April and May. This period coincides with the main breeding season of this bird as reported by Roberts (1991) for March and April and by Khan (2010) for March through May. Their habitat were regularly visited and observed for any signs of nest construction before the start of the expected breeding season i.e., the first week of March 2009. Local hunters assisted in locating bird nests. Search of these nests was carried out by walking on existing pathways across a forest plantation, while randomly encountered female birds were followed to reach their nests. After confirming an active nest, it was marked and allotted a specific number. Nest location, general appearance and its structure were noted. Marked nests were visited twice a week with maximum precaution to minimize disturbance to the nesting birds. The frequency of observations/visits was increased (on alternate days) after egg laying in order to record incubation periods and time of hatching. The information recorded include dates of laying first and last eggs, the number of eggs, the date of hatching and the number of hatched eggs.

    In order to determine diet composition of the Grey Francolin, 13 (10 males and 3 females) freshly hunted birds were procured from local hunters from the Jhelum and Rawalpindi districts of the Pothwar Plateau. Simultaneously, surveys of the habitat of the hunted birds were carried out to collect reference vegetation, e.g. grasses, shrubs, trees and farm crops, for use to identify the gut contents. The crops and gizzards were dissected to observe their contents in situ and subsequently the contents were segregated into identifiable components. Seeds and plant parts were identified from the Weed Management Programme of the National Agricultural Research Center (NARC), Islamabad, while the insect samples were identified with technical assistance from the Biosystematics Laboratory, Department of Entomology, PMAS Arid Agriculture University Rawalpindi and the National Insect Museum, NARC, Islamabad. All the constituents of the crops and gizzards contents were quantified by measuring their dry weights.

    The transect data recorded during this nine-month study period are grouped into climatically distinctive seasons of the area, i.e., winter (mid-November through February), spring (March through April), summer including the monsoon season (May through mid-September) and autumn (mid-September through mid-November). Mean values mentioned in the text are given with their standard errors (SE). We used two-sample t-tests or one-way analyses of variance by employing Minitab-16 Statistical Software for analysis of differences between any two data sets.

    Records of bird sighting data from the line-transects revealed a mean perpendicular distance of birds from the observer in the cropland habitat of 5.7 ± 0.5 m (n = 34, range = 0.6–9.4 m) and in the forest habitat of 6.3 ± 0.5 m (n = 23, range = 0.6–9.1 m). Non-significant differences between these two data sets (t = 0.85, df = 55, p = 0.398) suggest that we can assume that the birds were equally detectable in both habitats. The observed flock sizes (n) with frequencies of their occurrence (given in parenthesis) are as follows: 1(20), 2(20), 3(7), 4(1), 5(5), 6(1), 7(2) and 9(1).

    The monthly/seasonal data of the Grey Francolin recorded by direct sightings from all the six transects (three each in the forest and agriculture habitats) are detailed in Tables 1 and 2. The density was estimated as 1.59 ± 0.39 (mean ± standard error) birds per ha in the cropland habitat and 0.87 ± 0.14 birds per ha in the forest track. Although a higher density was evident in the agricultural transects compared to the forest area, a comparison of the overall number of birds seen in these two habitats showed non-significant difference at the 0.05 level of significance (t = 1.82, df = 16, p = 0.098). Season-wise segregation of data indicated the highest population estimates in the cropland areas during the spring season, whereas in the forest habitat the population was (more or less) equally high during the winter and summer months. However, there were minimum sightings during the autumn season in both habitats (Fig. 2).

    Table  1.  Record of population densities of the Grey Francolin (Francolinus pondicerianus)in cropland habitats of the Gujar Khan area in the Rawalpindi district, Pakistan. Each transect has a length of 500 m and a width of 50 m with an area of 2.5 ha.
    Season Month Transect A Transect Transect Mean density (n per ha)(Mean ± SE)
    n Crop and vegetation n Crop and vegetation n Crop and vegetation Monthly Seasonal
    Autumn Oct 2 Lentil (Lens culinaris) 3 Mustard (Brassica campestris), wheat (Triticum aestivum) 1 Fallow fields 0.80 ± 0.23 0.80 ± 0.23
    Winter Nov 1 Fallow fields 1 Mustard, Wheat Fallow fields 0.27 ± 0.13 1.70 ± 0.61
    Dec 2 As above 15 As above 13 Mustard, wheat, lentil, tara mira (Eruca sativa) 4.00 ± 1.62
    Jan As above 9 Mustard, wheat, lentil 5 As above 1.87 ± 1.04
    Feb 2 As above As above 3 As above 0.67 ± 0.35
    Spring Mar 2 As above 5 As above 5 As above 1.60 ± 0.40 2.13 ± 0.40
    Apr 4 As above 7 Fallow fields 9 Fallow fields 2.67 ±0.58
    Summer May 2 As above 3 As above 7 As above 1.60 ± 0.61 1.20 ± 0.34
    Jun 1 Millet (Pennisetum americanum), sorghum (Sorghum bicolor) 2 As above 3 Millet, sorghum 0.80 ± 0.23
    Density (n per ha)Mean ± SE 0.71 ± 0.15 2.00 ± 0.63 2.04 ± 0.54 1.59 ± 0.39
     | Show Table
    DownLoad: CSV
    Table  2.  Record of population densities of the Grey Francolin (Francolinus pondicerianus)in scrub forest habitats of the Gujar Khan area in the Rawalpindi district, Pakistan. Each transect has a length of 500 m and a width of 50 m with an area of 2.5 ha.
    Season Month Transect A Transect B Transect C Mean density(n per ha)(Mean ± SE)
    n Crop and vegetation n Crop and vegetation n Crop and vegetation Monthly Seasonal
    Autumn Oct Acacia modesta (dominant), A. nilotica (rare) Acacia modesta (dominant), Zizyphus mausitiana, Typha angustata (common), Dalbergia sissoo (rare), Desmostachia bipinnata (scattered in patches) along a small wetland 5 Acacia modesta (dominant), Dalbergia sissoo, A. nilotica, Zizyphus mausitiana and Desmostachia bipinnata (common) 0.67 0.22 ± 0.22
    Winter Nov As above As above 5 As above 0.67 0.93 ± 0.45
    Dec As above As above 11 As above 1.47
    Jan 1 As above As above 9 As above 1.33 ± 1.14
    Feb As above As above 2 As above 0.27
    Spring Mar As above 3 As above 2 As above 0.67 ± 0.35 0.73 ± 0.32
    Apr 5 As above As above 1 As above 0.80 ± 0.61
    Summer May 2 As above As above 4 As above 0.80 ± 0.46 1.00 ± 0.35
    Jun As above 4 As above 5 As above 1.20 ± 0.61
    Density (n per ha)Mean ± SE 0.36 ± 0.23 0.31 ± 0.21 1.96 ± 0.44 0.87 ± 0.14
     | Show Table
    DownLoad: CSV
    Figure  2.  Seasonal variation in population densities of the Grey Francolin (Francolinus pondicerianus)in cropland and forest tracks of agro-ecosystem of the Pothwar Plateau, Pakistan

    We assume that the variable distribution of birds within the agriculture habitat (Table 1) is steered by the prevalence of cultivation. The agricultural fields occurring in transect A were recorded as uncultivated during most of the study period, except in October and June, whereas the crop fields of transects B and C remained under seasonal crops except during the periods of land preparation for sowing the next crop. Transect-wise densities of the Grey Francolin in these agricultural areas were recorded at 0.71 ± 0.15, 2.00 ± 0.63 and 2.04 ± 0.54 birds per ha in transects A, B and C, respectively. The breakup of data by season showed estimates in the following descending order: 2.13 ± 0.40 birds per ha during spring, 1.70 ± 0.61 in winter, followed by 1.20 ± 0.34 in the summer and the lowest estimates, i.e. 0.80 ± 0.23 birds per ha during the autumn season (Table 1).

    The record of population densities of the Grey Francolin in scrub forest tracks are presented in Table 2. The estimates of bird populations in the first two transects (A and B) were 0.36 ± 0.23 and 0.31 ± 0.21 birds per ha, respectively, whereas the third transect, C, had the maximum population i.e. 1.96 ± 0.44 birds per ha. The vegetation of transect A was that of a monoculture, i.e., Acacia modesta with a rare occurrence of A. nilotica. The flora of transects B and C had an almost similar inventory of plant species, except for the presence of a small wetland surrounded by Typha angustata along transect B. However, there were higher densities of Desmstachia bipinnata and Dalbergia sissoo on transect C where the highest population of Grey Francolins was recorded (Table 2). Although these plant species do not contribute to the food supply of the Grey Francolin (Table 5), the relationship can presumably be attributed to meet the roosting and sheltering needs of the species. The seasonal estimates of the Grey Francolin in the scrub forest habitat revealed occurrences of the birds at 0.22 ± 0.22, 0.93 ± 0.45, 0.73 ± 0.32 and 1.0 ± 0.35 birds per ha during the autumn, winter, spring and summer seasons, respectively.

    In total, six active nests were found, four in forest tracks A, B, C and D and two in wild vegetation on the boundaries of cultivated fields (E and F). All the nests were situated on the ground in vegetation mainly comprised of Desmostachia bipinnata, Acacia modesta, Imperata cylindrical, Zizipus jujuba and Euphorbia spp. (Table 3, Fig. 3). The mean egg laying period was estimated at 6 ± 0.36 days (range 5–7 days), which corresponds to the recorded average clutch size of 7 ± 0.36 (range of 6–8 eggs), indicating that this population was probably laying eggs daily.

    Table  3.  Location and nest material of the Grey Francolin (Francolinus pondicerianus)found in an agro-ecosystem of the Pothwar Plateau, Pakistan
    Nest No. (habitat) Date of egg laying Nest material Vegetation at nest location
    A (Forest) 23-03-09 Desmostachia bipinnata Acacia modesta
    B (Forest) 01-04-09 Twigs of Imperata cylindrical Desmostachia bipinnata, Acacia modesta, Euphorbia
    C (Forest) 01-05-09 Imperata cylindrical, Desmostachia bipinnata Desmostachia bipinnata, Typha angustata
    D (Forest) 07-05-09 Desmostachia bipinnata Imperata cylindrical, Desmostachia bipinnata, Acacia modesta
    E (Cropland) 18-04-09 Near field edge in soil with Desmostachia bipinnata Desmostachia bipinnata, Zizipus jujuba on field edge
    F (Cropland) 23-04-09 On ground in soil with plant stems/debris On field edge adjacent to fallow-land
     | Show Table
    DownLoad: CSV
    Figure  3.  A nest containing six eggs of the Grey Francolin (Francolinus pondicerianus)observed in an agro-ecosystem of the Pothwar Plateau, Pakistan

    The average incubation period recorded in the present study was of 20.6 ± 0.50 days (range 19–22 days). Out of a total of 42 eggs recorded, 32 hatched (75.39% success) with a mean number of 5.33 ± 1.22 eggs hatched per clutch. The mean number of fledglings was estimated as 3.83 ± 0.83 per clutch, indicating an overall success of 63.08% (Table 4).

    Table  4.  Breeding performance of the Grey Francolin (Francolinus pondicerianus)in an agro-ecosystem of the Pothwar Plateau, Pakistan
    Nest No. (habitat) Egg laying period (days) Clutch sizE Incubation period (days) Hatching success[n (%)] Survival or fledgling success [n (%)]
    A (Forest) 5 7 21 6 (85.7%) 4 (66%)
    B (Forest) 5 6 20 6 (100%) 6 (100%)
    C (Forest) a 7 7 0 0 (0%) 0 (0%)
    D (Forest) 6 6 19 4 (66.7%) 4 (100%)
    E (Cropland) 7 8 22 8 (100%) 5 (62.5%)
    F (Cropland) 6 8 21 8 (100%) 4 (50%)
    Mean ± SE (n= 6) 6 ± 0.36 7 ± 0.36 20.6 ± 0.50 5.33 ± 1.22 (75.39%) 3.83 ± 0.83 (63.08 %)
    a The female left the nest during incubation.
     | Show Table
    DownLoad: CSV

    The results of the gut content analysis suggest it an omnivorous and euryphagus species, with a limited range in food selection (Table 5). The overall diet of this population consisted of ten identifiable food items, consisting of seven plant species, two insects and one grit. However, not more than five identifiable food items were found in any of the individual gut samples. The main components of the diet were seeds of cultivated crops, wild vegetation and arthropods. In total, 42 incidences of occurrence of different identifiable food contents were recorded. Among these were 14 records (33.3%) of cultivated field crops, 10 records (23.8%) of wild vegetation, 7 records (16.7%) of insects and 11 records (26.2%) of grits. Proportions of different food items (by weight) were estimated as 45.6% plant seeds, 12.4% plant leaves, 20.5% insects, 15.7% grits and 5.7% of unidentified material. The proportion of food contents derived from cultivated crops and wild vegetation was found to be in the ratio of 61:39.

    Table  5.  Composition of different food items recovered from crops and gizzards of the Grey Francolin (Francolinus pondicerianus) inhabiting an agro-ecosystem of the Pothwar Plateau, Pakistan
    Sr. No. Timeof hunt Month (Season)a Habitat Body wt. (g) Dry wt. of food contents (g) Wt. of food contents per 100 g body wt. Millet (seeds) Sorghum(seeds) Vilayati Kikar (leaves) Mung Beans(seeds) Acaciasp. (seeds) Baru (seeds) Carthemus oxycantha(seeds) Termites Black Ant Grit Unidentified
    1 Evening Nov (W) Cropland 250 7.0 2.80 2.5 2.5 2.0
    2 Evening Dec (W) Forest 263 12.0 4.56 1.5 2.5 2.0 2.5 3.0 0.5
    3 Evening Dec (W) Cropland 210 9.0 4.29 1.7 3.2 1.3 1.0 1.5 0.3
    4 Evening Dec (W) Cropland 242 6.0 2.48 2.5 1.5 1.0 0.8 0.2
    5 Morning Jan (W) Forest 200 3.0 1.50 1.7 0.9 0.4
    6 Evening Jan (W) Cropland 267 13.0 4.87 5.0 5.5 1.0 1.0 0.5
    7 Evening Feb (W) Cropland 247 13.0 5.26 5.5 2.0 1.2 4.0 0.3
    8 Evening May (S) Cropland 118 1.0 0.85 0.6 0.25 0.15
    9 Morning May (S) Cropland 168 3.0 1.79 0.55 0.5 1.0 0.95
    10 Morning May (S) Cropland 155 2.0 1.29 0.8 0.5 0.7
    11 Evening July (S) Forest 257 8.5 3.31 1.5 5.5 1.0 0.5
    12 Evening July (S) Forest 240 7.5 3.13 2.0 3.8 1.5 0.2
    13 Morning July (S) Forest 235 2.0 0.85 1.2 0.5 0.3
    Frequency of occurrence of identifiable food items n (%) 5 (11.9) 4 (9.5) 5 (11.9) 5 (11.9) 3 (7.1) 1 (2.4) 1 (2.4) 4 (9.5) 3 (7.1) 11 (26.2) 12
    Mean ± SE 219.4 ± 12.9 6.7 ± 1.2 2.8 ± 0.4 1.4 ± 0.3 3.9 ± 0.8 2.2 ± 0.9 1.6 ± 0.5 1.6 ± 0.2 2.5 ± 0.0 1.7 ± 0.0 3.0 ± 1.2 2.0 ± 1.0 1.2 ± 0.2 0.4 ± 0.07
    Total weight (% proportion) b 7.1(8.2%) 15.5(17.8%) 10.8(12.4%) 8.2(9.4%) 4.7(5.4%) 2.5(2.9%) 1.7(2.0%) 11.8(13.6%) 6.0(6.9%) 13.7(15.7%) 5.0(5.7%)
    a W = winter; S = summer.
    b Totally, seeds = 45.6%; leaves = 12.4%; insects = 20.5%.
    Latin names of some of food items: Millet, Pennisetum typhoideum; Sorghum, Sorghum bicolor; Vilayati Kikar, Prosopis juliflora; Mung Beans, Phaseolus radiatus; Baru, Sorghum halepense; Termites, Coptotermes formosanus; Black Ant, Formica rufa.
     | Show Table
    DownLoad: CSV

    Adjusted mean weight (g per 100 g body weight) of food contents gleaned from crops and gizzards of hunted birds in the morning (1.36 ± 0.2, n = 4) was significantly less (F1, 12 = 8.83, p = 0.013) than that of birds hunted in the evening (3.51 ± 0.46, n = 9). The increased weight of food content in birds shot in the evening suggests that the birds fill their guts before sunset in order to meet the energy requirements during the non-feeding overnight period.

    Weights of the gut contents (g per 100 g body weight) of the birds collected from cropland and forest habitats, respectively, were recorded as 2.95 ± 0.59 g (n = 8) and 2.67 g ± 0.67 (n = 5). The difference between the two habitats was not significant (F1, 12 = 0.10, p = 0.76). These figures do not indicate any quantitative difference in the availability of food sources between both habitats.

    The pattern of habitat preference recorded for the Grey Francolin in the present study area is supported by a similar study conducted in a scrub habitat of the Pothwar area (Mahmood et al., 2010), suggesting that this species selects habitats dominated by trees and shrubs. The asymmetrical prevalence of the Grey Francolin in the study transects in both habitats (cropland and scrub forest) could be attributed to variable vegetation/crops providing food and shelter to the birds. The data from transect C in the forest habitat show that large numbers of birds in vegetation dominated by Acacia modesta, Dalbergia sissoo and D. bipinnata. Wijeyamohan et al. (2003) in Sri Lanka reported that low bushes and thorn scrubs, especially Acacia eburnean, provided the most favorite habitats of the Grey Francolin. The activity of the bird was lower in the area (transect B in the forest) along a small area of wetland, although it had a good diversity of natural vegetation containing Acacia modesta (dominant), Zizyphus mausitiana, Typha angustata (common), Dalbergia sissoo (rare) and Desmostachia bipinnata (scattered in patches). Some earlier reports have also shown that these birds do not prefer habitats/vegetation close to water sources and that they meet their water requirement from succulent vegetation, dew and/or insects (Yeatter, 1934; Mc Crow, 1982; Ullah, 1991; Mian and Wajid, 1994; Kamal, 2000).

    The relatively higher density estimated in the cropland habitat (2.13 ± 0.40 birds per ha) during the spring season, corresponding with the maturity of wheat and related crops, is supported by the report of Mian and Ghani (2007) that population density of the Grey Francolin increases with an increase in vegetation cover. The higher estimates in the cropland is also supported by Roberts (1991) who stated that the bird restricts its movement to areas rich in crop grains, weed seeds, succulent vegetation and grits. Relatively higher populations in the forest habitat during summer and winter months could be attributed to the availability of cover for protection from extreme weather conditions. However, populations of this species in the larger scrub-forests on the sub-mountainous tracks of Kheri Murat, Fateh Jang and Attock districts (Mann and Choudhry, 2000) and the Lehri Wildlife Park, in the Jhelum district of the Pothwar Plateau (Mahmood et al., 2010), were estimated at 0.05 and 0.47 bird per ha, respectively. The current population is higher than those reported (Kamal, 2000) in the different agro-ecological sub-habitats of the Mardan district in Khyber Pakhtunkhwa Province, Pakistan (i.e., 0.98 bird per ha in the scrub forest, 0.28 in an irrigated forest plantation, 1.3 in a cropland adjoining a scrub forest, 0.25 in a cropland, 0.1 in a cropland adjoining a wetland and 0.48 in an alkaline waste/cropland). It is also higher than that in a tropical thorn forest and irrigated forest of a protected area (Lal Suhanra National Park) located in the desert ecology of Chohlistan, Bahawalpur (Khan, 2010), i.e., 0.06 bird per ha.

    The higher population density of the Grey Francolin in our cropland area during the winter and spring could also be related to the abundance of food in this habitat, because at that time the crops of Brassica campestris (sarsoon), Triticum aestivum (wheat), Lens culinaris (masoor) and Eruca sativa (taramira)were at their different growth stages and had become sources of succulent leaves, small grains as well as hosting high populations of insects. The diet of the Grey Francolin during the winter months has reportedly been dominated by mustard seeds and leaves in Layyah, Punjab (Mian and Wajid, 1994). Our findings are also supported by Kaiser (1998), who stated that Grey Partridges (Perdix perdix) occur in large numbers on farmlands and fields with a high diversity of crops. Similarly, maize or oil seeds, winter cereals and beans have been reported as preferred habitat of the Grey Partridge in Europe (Bro et al., 2004).

    The current data on breeding biology provide only a limited scope from which to draw any conclusion but, as a first report on this aspect of the species in the Pothwar Plateau, it justifies its documentation. Our limited data is being reviewed in the light of published information, with the aim to consider it in designing any large-scale study.

    In the present study, all nests were found on the ground in vegetation mainly comprising Desmostachia bipinnata, Acacia modesta, Imperata cylindrical, Zizipus jujuba and Euphorbia spp. These records agree with the observations by Islam (1999) on the location of the Grey Francolin in North America and that of Johnsgard (1973) who reported the occurrence of francolin nests under some protective cover for temperature modulation. However, Bump and Bump (1964) did find a Grey Francolin nest in a stacked pile of sorghum < 1.5 m above the ground in India. According to Sharma (1983), grasslands and plowed fields, as well as Euphorbia spp. were selected as nesting sites by Grey Francolins during the breeding season in Rajasthan, India.

    The recorded egg-laying span (from 23 March to 7 May) of the current population coincided with earlier observations, i.e., for March to September (Ali and Ripley, 1969), March to April (Roberts, 1991) and mid-March to early May (Khan, 2010). The present study does not support the report by Waite (1948) from the salt-range part of the Pothwar Plateau suggesting that breeding may occur in December. The mean egg laying period of 6 ± 0.36 days and the average clutch size of 7 ± 0.36 eggs suggest that this population was probably laying eggs daily. The breeding behavior of this species studied in the Pakistan desert (Khan, 2010) has revealed that the number of days required to complete the clutch was 50% larger than clutch size, suggesting that the birds were not laying eggs daily. The average clutch size recorded in the present study is more or less within the range of clutch sizes reported by Baker (1921) of 4–10 eggs, Clark (1901), 8–10 eggs, Sharma (1983), 6–7 eggs. However, Ullah (1991) recorded clutch sizes in the range of 4–16 eggs and a study by Khan (2010) indicated an average clutch size of 7.22 ± 1.23 eggs (range 2–13, n = 76). A variation in clutch size from 4 to 10 eggs also has also been reported by Baker (1935), Long (1981) and Bump and Bump (1964). The maximum number of 8 eggs recorded in nests situated in wild vegetation on field boundaries of cropland habitats suggests that birds, living near food sources, have higher clutch sizes. The relationship between clutch size and food availability has also been documented for agricultural farms in Rajasthan, India (Sharma, 1983). The average incubation period (20.6 ± 0.50 days, range 19–22 days) recorded in the present study falls within the range of this duration as reported by Khan (2010) in a desert population (17.68 ± 2.10 days, range 16–21 days, n = 54). The current estimate also agrees with the values of this parameter in a range of 18–21 days, reported by Bump and Bump (1964), Ali and Ripley (1969) and Roberts (1991). The estimated hatching success (75.39%) is evidently higher than that (44.0 ± 3.36%, n = 85 clutches, 598 eggs) reported by Khan (2010) but any conclusion could not be drawn due to the low number of clutches available for this study. The same limitation may apply to the estimated overall breeding/fledging success (63.08%) of this study against the estimated (37.0 ± 3.25%) breeding success reported by Khan (2010). Khan (2010) has also reported a direct impact of rainfall on overall breeding parameters, indicating that the average clutch size and number of nestlings and fledgings increases with higher rainfall. The apparently higher breeding success in the Pothwar Plateau population could be attributed to the higher mean rainfall than that found in the desert of Bahawalpur, Pakistan.

    The present study could identify only seven plant species and two insect species in the gut contents of the Grey Francolin hunted in agro-ecosystems of the Pothwar Plateau. The small sample size of the current study does not allow a meaningful comparison with any other study. However, an overview suggests that the number of identifiable food species is less than that reported (35 plant species and 16 animal species) from the Lal Suhanra National Park, located in a desert ecosystem of Bahawalpur, Pakistan (Khan, 2010). These are also smaller than the wide food base of this species reported by other investigators (Faruqui et al., 1960; Ullah, 1991; Chaudhry and Bhatti, 1992; Mian and Wajid, 1994). The present analysis suggests that insects, constituting 20.5% of the diet, is the only animal based food consumed by this group of francolins. The occurrence of only insects in the present samples is supported by Mian and Wajid (1994), but contrasts with the reports by Ali and Riplay (1983), Roberts (1991) and Choudhry and Bhatti (1992) who have reported the presence of reptiles and earthworm along with insects.

    Our study shows that the Grey Francolin of the Pothwar agro-ecosystem is mainly associated with wild vegetation, inter-dispersed with cultivated tracks. The diet composition of the species suggests that it has little preference for any specific food species; rather, it consumes/prefers different food items depending upon their temporal and spatial availability in its habitat. Existing populations can be sustained provided the traditional system of cultivation and natural vegetation is maintained. The current rate of expansion in land use for urbanization and agriculture to meet fast growing human needs warrants extensive surveys of the Grey Francolin in the Pothwar ecology to assess its status, studying its feeding and breeding habit requirements, thus allowing the formulation of an effective conservation plan for the species.

  • Alström, P., Rasmussen, P.C., Zhao, C., Xu, J., Dalvi, S., Cai, T., et al., 2016. Integrative taxonomy of the Plain-backed Thrush (Zoothera mollissima) complex (Aves, Turdidae) reveals cryptic species, including a new species. Avian Res. 7, 1.
    Avise, J.C., 2000. Phylogeography: the History and Formation of Species. Harvard University Press, Harvard.
    Baird, N., Etter, P.D., Atwood, T.S., Currey, M.C., Shiver, A.L., Lewis, Z.A., et al., 2008. Rapid SNP discovery and genetic mapping using sequenced RAD markers. PLoS One 3, e3376.
    Bennett, K., Provan, J., 2008. What do we mean by ‘refugia’? Quat. Sci. Rev. 27, 2449–2455.
    Bryc, K., Auton, A., Nelson, M.R., Oksenberg, J.R., Hauser, S.L., Williams, S., et al., 2010. Genome-wide patterns of population structure and admixture in West Africans and African Americans. Proc. Natl. Acad. Sci. USA 107, 786–791.
    Camacho, C., Coulouris, G., Avagyan, V., Ma, N., Papadopoulos, J., Bealer, K., et al., 2009. BLAST+: architecture and applications. BMC Bioinf. 10, 421.
    Campagna, L., Toews, D.P.L., 2022. The genomics of adaptation in birds. Curr. Biol. 32, R1173–R1186.
    Chen, X.Y., He, F., 2009. Speciation and endemism under the model of island biogeography. Ecology 90, 39–45.
    Chen, Y., Li, R., Sun, J., Li, C., Xiao, H., Chen, S., 2022. Genome-wide population structure and selection signatures of Yunling Goat based on RAD-seq. Animals 12, 2401.
    Cheviron, Z.A., Whitehead, A., Brumfield, R.T., 2008. Transcriptomic variation and plasticity in rufous-collared sparrows (Zonotrichia capensis) along an altitudinal gradient. Mol. Ecol. 17, 4556–4569.
    Choi, J.Y., Dai, X., Alam, O., Peng, J.Z., Rughani, P., Hickey, S., et al., 2021. Ancestral polymorphisms shape the adaptive radiation of Metrosideros across the Hawaiian Islands. Proc. Natl. Acad. Sci. USA 118, e2023801118.
    Collar, N., Robson, C., de Juana, E., 2019. White-browed laughingthrush (Garrulax sannio). In: del Hoyo, J., Elliott, A., Sargatal, J., Christie, D.A., de Juana, E. (Eds.), Handbook of the Birds of the World Alive. Lynx, Barcelona.
    Cooper, E.A., Uy, J.A.C., 2017. Genomic evidence for convergent evolution of a key trait underlying divergence in island birds. Mol. Ecol. 26, 3760–3774.
    Cromie, G.A., Hyma, K.E., Ludlow, C.L., Garmendia-Torres, C., Gilbert, T.L., May, P., et al., 2013. Genomic sequence diversity and population structure of Saccharomyces cerevisiae assessed by RAD-seq. G3 Genes Genom. Genet. 3, 2163–2171.
    Cun, Y-Z., Wang, X-Q., 2010. Plant recolonization in the Himalaya from the southeastern Qinghai-Tibetan Plateau: geographical isolation contributed to high population differentiation. Mol. Phylogenet. Evol. 56, 972–982.
    Darwin, C., 1859. On the Origin of Species. John Murray, London.
    Dixon, P., Palmer, M.W., 2003. VEGAN, a package of R functions for community ecology. J. Veg. Sci. 14, 927–930.
    Eaton, D.A.R., Ree, R.H., 2013. Inferring phylogeny and introgression using RADseq data: an example from flowering plants (Pedicularis: Orobanchaceae). Syst. Biol. 62, 689–706.
    Emerson, B.C., 2002. Evolution on oceanic islands: molecular phylogenetic approaches to understanding pattern and process: phylogenetic investigation of island evolution. Mol. Ecol. 11, 951–966.
    Fick, S.E., Hijmans, R.J., 2017. WorldClim 2: new 1-km spatial resolution climate surfaces for global land areas. Int. J. Climatol. 37, 4302–4315.
    Gouin, N., Grandjean, F., Bouchon, D., Reynolds, J.D., Souty-Grosset, C., 2001. Population genetic structure of the endangered freshwater crayfish Austropotamobius pallipes, assessed using RAPD markers. Heredity 87, 80–87.
    Govindaraju, D.R., 1988. Relationship between dispersal ability and levels of gene flow in plants. Oikos 52, 31–35.
    Grada, A., Weinbrecht, K., 2013. Next-generation sequencing: methodology and application. J. Invest. Dermatol. 133, 1–4.
    Guillot, G., Santos, F., 2010. Using AFLP markers and the Geneland program for the inference of population genetic structure. Mol. Ecol. Resour. 10, 1082–1084.
    Hanes, M.M., Shell, S., Shimu, T., Crist, C., Machkour-M’Rabet, S., 2022. The phylogeographic history of Megistostegium (Malvaceae) in the dry, spiny thickets of southwestern Madagascar using RAD-seq data and ecological niche modeling. Ecol. Evol. 12.
    He, K., Jiang, X., 2014. Sky islands of southwest China. Ⅰ: an overview of phylogeographic patterns. Chin. Sci. Bull. 59, 585–597.
    He, K., Gutix errez, E.E., Heming, N.M., Koepfli, K., Wan, T., He, S., et al., 2019. Cryptic phylogeographic history sheds light on the generation of species diversity in skyisland mountains. J. Biogeogr. 46, 2232–2247.
    He, Y., Hou, X-Y., Li, C-X., Wang, Y., Ma, X-R., 2022. Soil microbial communities altered by titanium ions in different agroecosystems of pitaya and grape. Microbiol. Spectr. 10, e00907–e00921.
    Hijmans, R.J., Williams, E., Vennes, C., 2017. Geosphere: spherical trigonometry. R Package Version 1, 5–7.
    Jablonski, N.G., 1993. Quaternary environments and the evolution of primates in East Asia, with notes on two new specimens of fossil Cercopithecidae from China. IJFP 60, 118–132.
    Jiang, S., Luo, M.X., Gao, R.H., Zhang, W., Yang, Y.Z., Li, Y.J., et al., 2019. Isolation-byenvironment as a driver of genetic differentiation among populations of the only broad-leaved evergreen shrub Ammopiptanthus mongolicus in Asian temperate deserts. Sci. Rep. 9, 12008.
    Jin, X., He, M., Ferguson, B., Meng, Y., Ouyang, L., Ren, J., et al., 2012. Correction: an effort to use human-based exome capture methods to analyze chimpanzee and macaque exomes. PLoS One 7.
    Jordano, P., Godoy, J.A., 2000. RAPD variation and population genetic structure in Prunus mahaleb (Rosaceae), an animal-dispersed tree. Mol. Ecol. 9, 1293–1305.
    Kang, J., Ma, X., He, S., 2017. Population genetics analysis of the Nujiang catfish Creteuchiloglanis macropterus through a genome-wide single nucleotide polymorphisms resource generated by RAD-seq. Sci. Rep. 7, 2813.
    Kar, P.K., Vijayan, K., Mohandas, T.P., Nair, C.V., Saratchandra, B., Thangavelu, K., 2005. Genetic variability and genetic structure of wild and semi-domestic populations of Tasar Silkworm (Antheraea mylitta) Ecorace Daba as revealed through ISSR markers. Genetica 125, 173–183.
    Koenig, W.D., Knops, J.M.H., Pesendorfer, M.B., Zaya, D.N., Ashley, M.V., 2017. Drivers of synchrony of acorn production in the valley oak (Quercus lobata) at two spatial scales. Ecology 98, 3056–3062.
    Koppes, M., Gillespie, A.R., Burke, R.M., Thompson, S.C., Stone, J., 2008. Late quaternary glaciation in the Kyrgyz tien Shan. Quat. Sci. Rev. 27, 846–866.
    Kramer, I.M., Pfenninger, M., Feldmeyer, B., Dhimal, M., Gautam, I., Shreshta, P., et al., 2023. Genomic profiling of climate adaptation in Aedes aegypti along an altitudinal gradient in Nepal indicates nongradual expansion of the disease vector. Mol. Ecol. 32, 350–368.
    Lê, S., Josse, J., Husson, F., 2008. FactoMineR: an R package for multivariate analysis. J. Statist. Software 25, 1–18.
    Li, H., Durbin, R., 2009. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25, 1754–1760.
    (1000 Genome Project Data Processing Subgroup), Li, H., Handsaker, B., Wysoker, A., Fennell, T., Ruan, J., Homer, N., et al., 2009. The sequence alignment/map format and SAMtools. Bioinformatics 25, 2078–2079.
    Li, X., Ruhsam, M., Wang, Y., Zhang, H-Y., Fan, X-Y., Zhang, L., et al., 2022. Winddispersed seeds blur phylogeographic breaks: the complex evolutionary history of Populus lasiocarpa around the Sichuan Basin. Plant Divers. 45, 156–168.
    Liu, F., Guo, Q., Shi, H., Lv, F., Yu, Y., Lv, L., et al., 2014. Genetic diversity analysis of Perinereis aibuhitensis based on ISSR and SRAP markers of Chinese coast populations. Biochem. Syst. Ecol. 57, 262–269.
    Liu, C., Chen, H., Ren, Z., Zhang, C., Yang, X., 2019. Population genetic analysis of the domestic Bactrian camel in China by RAD-seq. Ecol. Evol. 9, 11232–11242.
    Liu, H., Wang, Z., Zhang, Y., Li, M., Wang, T., Su, Y., 2023. Geographic isolation and environmental heterogeneity contribute to genetic differentiation in Cephalotaxus oliveri. Ecol. Evol. 13, e9869.
    Losos, J.B., Ricklefs, R.E., 2009. Adaptation and diversification on islands. Nature 457, 80–83.
    Ma, J., Zhang, T., Wang, W., Chen, Y., Cai, W., Zhu, B., et al., 2022. Comparative transcriptome analysis of Gayal (Bos frontalis), Yak (Bos grunniens), and Cattle (Bos taurus) reveal the high-altitude adaptation. Front. Genet. 12, 778788.
    Manthey, J.D., Moyle, R.G., 2015. Isolation by environment in White-breasted Nuthatches (Sitta carolinensis) of the Madrean Archipelago sky islands: a landscape genomics approach. Mol. Ecol. 24, 3628–3638.
    Martin Cerezo, M.L., Kucka, M., Zub, K., Chan, Y.F., Bryk, J., 2020. Population structure of Apodemus flavicollis and comparison to Apodemus sylvaticus in northern Poland based on RAD-seq. BMC Genom. 21, 241.
    Martin, C.A., Armstrong, C., Illera, J.C., Emerson, B.C., Richardson, D.S., Spurgin, L.G., 2021. Genomic variation, population history and within-archipelago adaptation between island bird populations. R. Soc. Open Sci. 8, 201146.
    McCormack, J.E., Hird, S.M., Zellmer, A.J., Carstens, B.C., Brumfield, R.T., 2013. Applications of next-generation sequencing to phylogeography and phylogenetics. Mol. Phylogenet. Evol. 66, 526–538.
    Moysi, M., Ogolowa, B.O., Nikiforou, C., Dretakis, M., Purcell, J., Brelsford, A., et al., 2023. Genomic data reveal contrasting patterns of divergence among island and mainland birds of the Eastern Mediterranean. Ibis 165, 829–843.
    Murray, M.G., Thompson, W.F., 1980. Rapid isolation of high molecular weight plant DNA. Nucleic Acids Res. 8, 4321–4326.
    Naimi, B., Hamm, N.A.S., Groen, T.A., Skidmore, A.K., Toxopeus, A.G., 2014. Where is positional uncertainty a problem for species distribution modelling? Ecography 37, 191–203.
    Nanninga, G.B., Saenz-Agudelo, P., Manica, A., Berumen, M.L., 2014. Environmental gradients predict the genetic population structure of a coral reef fish in the Red Sea. Mol. Ecol. 23, 591–602.
    Neigel, J.E., Avise, J.C., 1986. Phylogenetic relationships of mitochondrial DNA under various demographic models of speciation. In: Nevo, E., Karlin, S. (Eds.), Evolutionary Processes and Theory. Academic Press, New York, pp. 513–534.
    O’Brien, K.A., Simonson, T.S., Murray, A.J., 2020. Metabolic adaptation to high altitude. Curr. Opin. Endocr. Metab. Res. 11, 33–41.
    Ortiz, D., Pekx ar, S., Dianat, M., 2022. Phylogenomics and loci dropout patterns of deeply diverged Zodarion ant-eating spiders suggest a high potential of RAD-seq for genuslevel spider phylogenetics. Cladistics 38, 320–334.
    Pedersen, T.L., 2022. Ggforce: Accelerating ‘ggplot2’. .
    Pickrell, J.K., Pritchard, J.K., 2012. Inference of population splits and mixtures from genome-wide allele frequency data. PLoS Genet. 8, e1002967.
    Pullin, A.S., 2002. Conservation Biology. Cambridge University Press, Cambridge.
    Qu, Y., Luo, X., Zhang, R., Song, G., Zou, F., Lei, F., 2011. Lineage diversification and historical demography of a montane bird Garrulax elliotii—implications for the Pleistocene evolutionary history of the eastern Himalayas. BMC Evol. Biol. 11, 174.
    Qu, Y., Chen, C., Chen, X., Hao, Y., She, H., Wang, M., et al., 2021. The evolution of ancestral and species-specific adaptations in snowfinches at the Qinghai-Tibet Plateau. Proc. Natl. Acad. Sci. USA 118, e2012398118.
    Reif, J.C., Hamrit, S., Heckenberger, M., Schipprack, W., Peter Maurer, H., Bohn, M., et al., 2005. Genetic structure and diversity of European flint maize populations determined with SSR analyses of individuals and bulks. Theor. Appl. Genet. 111, 906–913.
    Sazzini, M., Abondio, P., Sarno, S., Gnecchi-Ruscone, G.A., Ragno, M., Giuliani, C., et al., 2020. Genomic history of the Italian population recapitulates key evolutionary dynamics of both Continental and Southern Europeans. BMC Biol. 18, 51.
    Schuster, S.C., 2008. Next-generation sequencing transforms today’s biology. Nat. Methods 5, 16–18.
    Servedio, M.R., 2000. Reinforcement and the genetics of nonrandom mating. Evolution 54, 21–29.
    Stolpovsky, Y.A., Kol, N.V., Evsyukov, A.N., Ruzina, M.N., Shimiit, L.V., Sulimova, G.E., 2010. Analysis of the genetic structure of Tuvinian short-fat-tailed sheep populations with the use of the ISSR-PCR method. Russ. J. Genet. 46, 1462–1470.
    Stuart, S.N., Chanson, J.S., Cox, N.A., Young, B.E., Rodrigues, A.S., Fischman, D.L., et al., 2004. Status and trerds of amphibian declines and extinctions worldwide. Science 306, 1783–1786.
    Sun, Y., Buhler, J., 2007. Designing patterns for profile HMM search. Bioinformatics 23, e36–e43.
    Tang, Y., Wang, T., Zhang, A., Zhu, J., Zhou, T., Zhou, Y.L., et al., 2023. ADCY9 functions as a novel cancer suppressor gene in lung adenocarcinoma. J. Thorac. Dis. 15, 1018–1035.
    van Dijk, E.L., Auger, H., Jaszczyszyn, Y., Thermes, C., 2014. Ten years of nextgeneration sequencing technology. Trends Genet. 30, 418–426.
    Wang, Y.J., Cheng, H., Edwards, R.L., 2001. A high-resolution absolute-dated Late Pleistocene monsoon record from hulu cave, China. Science 294, 2345–2348.
    Wang, B., Xie, X., Liu, S., Wang, X., Pang, H., Liu, Y., 2017. Development and characterization of novel microsatellite markers for the Common Pheasant (Phasianus colchicus) using RAD-seq. Avian Res. 8, 4.
    Wang, Y., Feijo, x A., Cheng, J., Xia, L., Wen, Z., Ge, D., et al., 2021. Ring distribution patterns—diversification or speciation? Comparative phylogeography of two small mammals in the mountains surrounding the Sichuan Basin. Mol. Ecol. 30, 2641–2658.
    Wei, T.Y., Simko, V., 2021. R Package ‘corrplot’: Visualization of a Correlation Matrix Using Ggplot2. Version 0.92. .
    Wen, L., Wang, Y., Fu, Y., Dai, B., 2017a. The complete mitochondrial genome of the white-browed laughingthrush Garrulax sannio (Passeriformes: Leiothrichidae). Mitochond. DNA Part A 28, 23–24.
    Wen, Z., Wu, Y., Ge, D., Cheng, J., Chang, Y., Yang, Z., et al., 2017b. Heterogeneous distributional responses to climate warming: evidence from rodents along a subtropical elevational gradient. BMC Ecol. 17, 17.
    Wen, L., Antalffy, J.M., Messenger, K., Zhao, W., Liu, F., Zou, F., et al., 2020. Spatial variation in morphology of the White-browed Laughingthrush (Garrulax sannio) and its relationship with climate conditions in southern China. J. Ornithol. 161, 701–711.
    Wickham, H., 2016. ggplot2: Elegant Graphics for Data Analysis. Springer-Verlag, Berlin.
    Willing, E-M., Hoffmann, M., Klein, J.D., Weigel, D., Dreyer, C., 2011. Paired-end RADseq for de novo assembly and marker design without available reference. Bioinformatics 27, 2187–2193.
    Wright, S., 1978. Evolution and the genetic of population. In: Variability within and among Natural Populations, vol. 4. University of Chicago Press, Chicago.
    Wu, Y., DuBay, S.G., Colwell, R.K., Ran, J., Lei, F., 2017. Mobile hotspots and refugia of avian diversity in the mountains of south-west China under past and contemporary global climate change. J. Biogeogr. 44, 615–626.
    Xiang, L., Li, X-L., Wang, X-S., Yang, J., Lv, K., Xiong, Z-Q., et al., 2020. Genetic diversity and population structure of Distylium chinense revealed by ISSR and SRAP analysis in the three gorges reservoir region of the Yangtze river, China. Global Ecol. Conserv. 21, e00805.
    Yang, J., Lee, S.H., Goddard, M.E., Visscher, P.M., 2011. GCTA: a tool for genome-wide complex trait analysis. Am. J. Hum. Genet. 88, 76–82.
    Yang, W., Qi, Y., Fu, J., 2016. Genetic signals of high-altitude adaptation in amphibians: a comparative transcriptome analysis. BMC Genet. 17, 134.
    Yang, A., Chen, D., Wang, P., Fu, Y., Zhang, Z., 2017. Characterization of novel microsatellite markers of the Emei Shan Liocichla using restriction site-associated DNA sequencing. Avian Res. 8, 13.
    Yoichi, W., Takahashi, M., Nagano, A.J., Uehara, K., Abe, H., 2021. Evolutionary effects of geographic and climatic isolation between Rhododendron tsusiophyllum populations on the Izu Islands and mainland Honshu of Japan. Heredity 126, 859–868.
    Yu, G., Wang, L.G., Han, Y., He, Q.Y., 2012. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS 16, 284–287.
    Yu, L., Wang, G.D., Ruan, J., Chen, Y.B., Yang, C.P., Cao, X., et al., 2016. Genomic analysis of snub-nosed monkeys (Rhinopithecus) identifies genes and processes related to high-altitude adaptation. Nat. Genet. 48, 947–952.
    Zerbino, D.R., Birney, E., 2008. Velvet: algorithms for de novo short read assembly using de Bruijn graphs. Genome Res. 18, 821–829.
    Zhang, G., Li, C., Li, Q., Li, B., Larkin, D.M., Lee, C., et al., 2014. Comparative genomics reveals insights into avian genome evolution and adaptation. Science 346, 1311–1320.
    Zhang, S., Ding, J., Han, Z., Chen, S., Liu, Y., He, W., et al., 2022. Development of SSR markers and genetic diversity analysis based on RAD-seq technology among Chinese populations of Daphnia magna. Mol. Biol. Rep. 49, 4389–4397.
    Zhu, X., Zhou, S., Xu, X., Wang, J., Yu, Y., Yang, K., et al., 2017. Morphological differentiation in Asian honey bee (Apis cerana) populations in the basin and highlands of southwestern China. J. Apicult. Res. 56, 203–209.
  • Related Articles

Catalog

    Figures(5)

    Article Metrics

    Article views (5) PDF downloads (3) Cited by()

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return