Herman L Mays Jr, Bailey D McKay, Dieter Thomas Tietze, Cheng-Te Yao, Lindsey N Miller, Kathleen N Moreland, Fumin Lei. 2015: A multilocus molecular phylogeny for the avian genus Liocichla (Passeriformes: Leiothrichidae: Liocichla). Avian Research, 6(1): 17. DOI: 10.1186/s40657-015-0025-y
Citation: Herman L Mays Jr, Bailey D McKay, Dieter Thomas Tietze, Cheng-Te Yao, Lindsey N Miller, Kathleen N Moreland, Fumin Lei. 2015: A multilocus molecular phylogeny for the avian genus Liocichla (Passeriformes: Leiothrichidae: Liocichla). Avian Research, 6(1): 17. DOI: 10.1186/s40657-015-0025-y

A multilocus molecular phylogeny for the avian genus Liocichla (Passeriformes: Leiothrichidae: Liocichla)

More Information
  • Corresponding author:

    Herman L Mays Jr, maysh@marshall.edu

  • Received Date: 17 Feb 2015
  • Accepted Date: 20 Jul 2015
  • Available Online: 24 Apr 2022
  • Publish Date: 04 Aug 2015
  • Background 

    Historically the babblers have been assigned to the family Timaliidae but several recent studies have attempted to rest the taxonomy of this diverse passerine assemblage on a more firm evolutionary footing. The result has been a major rearrangement of the group. A well-supported and comprehensive phylogeny for this widespread avian group is an important part of testing evolutionary and biogeographic hypotheses, especially in Asia where the babblers are a key component of many forest ecosystems. However, the genus Liocichla is poorly represented in these prior studies of babbler systematics.

    Methods 

    We used a multilocus molecular genetic approach to generate a phylogenetic hypothesis for all five currently recognized species in the avian genus Liocichla. Multilocus DNA sequence data was used to construct individual gene trees using maximum likelihood and species trees were estimated from gene trees using Bayesian analyses. Divergence dates were obtained using a molecular clock approach.

    Results 

    Molecular data estimate a probable window of time for the origin for the Liocichla from the mid to late Miocene, between 5.55 and 12.87 Ma. Despite plumage similarities between the insular Taiwan endemic, L. steerii, and the continental L. bugunorum and L. omeiensis, molecular data suggest that L. steerii is the sister taxon to all continental Liocichla. The continental Liocichla are comprised of two lineages; a lineage containing L. omeiensis and L. bugunorum and a lineage comprised of L. phoenicea and L. ripponi. The comparatively early divergence of L. steerii within the Liocichla may be illusory due to extinct and therefore unsampled lineages. L. ripponi and L. phoenicea are parapatric with a Pleistocene split (0.07-1.88 Ma) occurring between an Eastern Himalayan L. phoenicea and a Northern Indochina distributed L. ripponi. L. bugunorum and L. omeiensis underwent a similar split between the Eastern Himalaya (L. bugunorum) and Central China (L. omeiensis) divided by the Hengduan Mountains.

    Conclusions 

    This study supports an origin of the Liocichla occurring sometime prior to the Miocene-Pliocene boundary, a period of significant climatic upheaval in Asia. The biogeographical patterns within the Liocichla mirror those of other birds in the region and allude to common geological and climatic drivers of avian diversification in Asia.

  • Communal roosting, commonly seen in birds, is the aggregation of more than two conspecifics spending the diurnal or nocturnal resting periods together (Bijleveld et al. 2010; Laughlin et al. 2014). However, the driving forces behind communal roosting and its adaptive value are still not fully understood. Several hypotheses that explain the benefits of communal roosting have been proposed, including reduced thermoregulation demand, lower predation risk, and greater foraging efficiency (Eiserer 1984; Beauchamp 1999).

    An increase in foraging efficiency is often considered as the main benefit of communal roosting in birds. Roost sites can act as information centres, whereby unsuccessful foragers can learn about the locations of food sources by following successful foragers (Ward and Zahavi 1973; Chou 2015). Locations where birds form roosts are often thought to be influenced by the proximity of food resources. Alternatively, the patch-sitting hypothesis proposes that roosts tend to be close to supplemental feeding areas, reducing the commuting cost from distant foraging sites to roost sites (Caccamise and Morrison 1986, 1988).

    The thermoregulation hypothesis proposes that the presence of nearby companion birds can reduce high energy costs during inactivity phases (McKechnie and Lovegrove 2002; McKechnie et al. 2006), and/or the physical structure of the roost can protect the roosting birds from detrimental weather conditions (Broom et al. 1976). A considerable body of studies found that the microhabitat temperature of roost sites was higher than that of non-roost sites (Clergeau and Quenot 2007; Erickson 2015), which emphasised that roost site selection in cold weather may be constrained by thermoregulatory considerations (Buttemer 1985).

    Many studies have demonstrated that the roosting behaviours of urban birds have strong associations with human activities, such as night light and land use (Khadraoui and Toews 2015; Daoud-Opit and Jones 2016; Saiyad et al. 2017). This could be attributed to predator avoidance, and specifically that roost sites found in urban areas with large open space and good illumination can aid the roosting birds in detecting their predators (Gorenzel and Salmon 1995; Peh and Sodhi 2002; Jaggard et al. 2015).

    There has been a myriad of research aimed to find out the environmental factors influencing communal roosting by comparing the characteristics between roosts and non-roosts (Clergeau and Simonnet 1996; Peh and Sodhi 2002; Yap et al. 2002; Erickson 2015; Jaggard et al. 2015). However, some of these studies offered contradictory findings. For instance, both positive and negative impacts of illumination on roosting have been recorded (Erickson 2015; Jaggard et al. 2015). These inconsistent findings indicate that variations among roosts may be significant, but have been insufficiently studied. Only a few studies have emphasised the seasonal variation of roosts (Barrows 1981; Buehler et al. 1991; Gorenzel and Salmon 1995). Despite the enhanced understanding of roost site selection, knowledge regarding seasonal switching remains limited.

    Yellow-crested Cockatoo (Cacatua sulphurea) is a critically endangered species endemic to Indonesia and Timor-Leste (BirdLife International 2018). As a result of the international pet trade, Yellow-crested Cockatoo was introduced to Hong Kong, China. A released and/or escaped population (c. 100‒150 individuals) has inhabited the remnant green patches in highly urbanised areas (Leven and Corlett 2004; Mok et al. 2012). We observed that as many as 102 cockatoos formed communally roosting flocks throughout the year (Wang 2020). The birds would first come to a pre-roost gathering which is defined as an aggregation of individuals that occurs before departure to the nocturnal roosts (Moore and Switzer 1998). But while the pre-roost gathering was always in the same place, the locations of communal roosts changed seasonally. To some extent, this pattern supported the hypothesis that the function of pre-roost gathering is to provide individual birds with the knowledge about the changeable roost locations (Zahavi 1971; Sonerud et al. 2002). However, the causes of this seasonal variation of roost sites are not understood.

    We hypothesised that the seasonal switching of roost sites was influenced by multiple environmental factors in an urban environment, and that roost sites would have different microhabitat characteristics. To test this, we assessed the microhabitat characteristics of roosts from three aspects: land use composition surrounding roost sites, microclimate temperature and human disturbance. We attempted to find out the crucial factors distinguishing different seasonal roosts and explore the possible causes of these associations. We hope that correlations between these environmental factors and roosts would contribute to the understanding of roost switching.

    The study was conducted in Hong Kong, which is located on the south coast of China (22° 09′‒22° 33′ N, 113° 50′‒114° 26′ E). It is a compact and highly populated city, with the majority of 7.48 million population residing on only 24.4% of 1110 km2 total land (Planning Department 2020). Hong Kong has a humid subtropical monsoon climate, tending towards temperate for nearly half the year. Due to monsoons and typhoons, 80% of the precipitation is concentrated between May and September (Hong Kong Observatory 2020). The mean daily minimum temperature varies monthly, ranging from 14.6 to 26.9 ℃. In contrast, the native habitats of Yellow-crested Cockatoo in Indonesia are almost entirely tropical. The mean daily minimum temperature remains fairly constant throughout the year, averaging 23.2 ℃. The minimum temperature variations between Hong Kong and Indonesia are relatively greater between December and March than other months, with Hong Kong being 5.6 ℃ to 8.6 ℃ colder during this season (Additional file 1: Table S1). This indicates a possible night-time cold stress of Yellow-crested Cockatoo in the winter of Hong Kong.

    The studied roosting flock, which can be regularly observed and easily tracked, assembles on the northern coast of Hong Kong Island. This district is a prime financial and commercial centre of the city, characterised by dense and tall buildings. It also contains two large urban parks, i.e. Hong Kong Park (8.2 ha) and Hong Kong Zoological and Botanical Gardens (5.6 ha), which are the major feeding grounds for Yellow-crested Cockatoo, providing abundant and various plant food sources (Fig. 1).

    Figure  1.  Map showing the seven roost sites (yellow circles with roost names) among those high-rises and the two adjacent urban parks (Map source: Google Earth image © 2019)

    The roosting population was counted once a week in clear evenings for two years (Year 2014/15 and 2015/16) from March 2014 to March 2016. As the evening progressed, birds in pairs or small flocks successively flew into the pre-roost gathering sites in Hong Kong Park from all directions. The observation lasted from late afternoon to early evening, beginning when the individuals flew to the pre-roost aggregation sites until no conspicuous flying and calling came from the flock at the final roost site (Davis 1955). We located the birds visually using binoculars (Nikon Monarch 5, 10 × 42), standing at a 30-m high observation tower in Hong Kong Park (Fig. 1). The observation tower afforded an unobstructed view to track the birds as they entered the pre-roost gathering sites and then settled into the roost sites, with a few exceptional locations where the birds could not be seen but their presence could later be identified by faecal droppings (Gorenzel and Salmon 1995). The number of cockatoos attending the pre-roost aggregation was counted every 5 min until the cockatoos flew into their roost site(s). The roosting flock size was counted at the roost site; otherwise, it was estimated from the final number obtained at the pre-roost site when not all individuals could be distinguished at the roost site. The pre-roosting aggregation would sometimes split up to occupy two to three roost sites, that the one greater in number was regarded as the major roost while the others were minor roosts.

    In total, seven roost sites were identified to be used by Yellow-crested Cockatoo during the 105 observational days. Five sites were located on trees (BT, CT, GT, PT and RT), one on a building roof (CR), and the other on a cluster of lamp posts (LP). Given the occupied months and night minimum air temperature of occupied days of each roost (Additional file 2: Fig. S1), we classified the roost sites into three distinct groups, i.e. spring roosts (GT and PT), summer roosts (BT, CR and LP), and winter roosts (CT and RT). As roost CR was not accessible, we excluded it from all microhabitat-related analyses.

    To investigate whether land use around roosts affects roost site selection, we measured the proportion of land use types within a 50-m radius plot around each roost delineated by Google Earth Pro (Google Inc.). We modified the classification of land uses defined by Yap et al. (2002), whereby "tree-dominated area" refers to lands covered by tree canopies, "built-up area" represents developed areas with buildings, and "open space" consists of areas covered by roads, low-growing vegetation or other vacant lands.

    To explore the level of human disturbance of each roost site, we quantified human disturbance using multiple anthropogenic components (Gorenzel and Salmon 1995; Peh and Sodhi 2002; Yap et al. 2002; Jaggard et al. 2015). Specifically, we collected data on (1) pedestrians: night-time occurrence of pedestrians passing by the roosts in 10-min durations; (2) traffic: night-time car volume on the nearest roads to the roosts in 10-min durations; (3) night-time noise level at the roost sites; (4) night-time light intensity at the roost sites; (5) mean roosting height above ground; and (6) distance to the nearest main road, building, streetlight and tree. We measured the mean number of pedestrians and number of vehicles of each roost by conducting a 10-min survey soon after the roosting flock settled down, once a month from July 2015 to February 2016. Similarly, we measured the noise levels and light intensity once a month, using an extendable sound meter and light meter, respectively. We held the instruments at the roosting height at four directions around each roost, soon after the roosting flock settled down. We recorded the light intensity three times per sampling evening and the average decibel for a total of 5-min per sampling evening. We measured the mean roosting height using a rangefinder (Nikon Forestry Pro). We estimated the distance to the nearest main road, building, streetlight, and tree on the Hong Kong GEOINFO MAP (www.map.gov.hk/gm/) since the measurements could not be accurately made on the ground. The mean values of the anthropogenic disturbance variables are given in Additional file 1: Table S1.

    As described above, the night-time temperature in the winter of Hong Kong is much lower than that in the native habitats of Yellow-crested Cockatoo. We measured the microclimate of each roost site when the winter roosts were occupied, and recorded the microhabitat temperature at each roost site by positioning a thermometer on an extendable pole at the roosting height. The temperature was recorded for 5 min soon after the roosting flock settled down on each roost on 12 randomly selected days from late November 2015 to January 2016, from which the mean values (Tmicro) were calculated.

    We compared the roosting flock size of each roost site using Kruskal-Wallis test followed by Wilcoxon rank sum tests. We conducted Pearson's chi-squared test to determine whether roost sites and land use types were independent of each other. We used Pearson residuals to measure the discrepancy between observed and expected values (Friendly 1994; Friendly and Meyer 2015). The formula is:

    Pearson residual=(actualexpected)/expected

    Cut-off points of Pearson residuals at ± 2 and ± 4 implied that the residuals were significant at α=0.05 and α=0.0001 levels respectively (Meyer et al. 2006). We performed a hierarchical clustering analysis using Ward's method to investigate how the roosts were grouped by human disturbance variables. Thereafter, we performed a partial least squares determinant analysis (PLS-DA) to maximise the variation between seasonal roost groups and identify the type of human disturbance important in roost site selection. PLS regression is designed explicitly for analytical situations where predictor variables are highly correlated and/or the sample size is smaller than the number of observations (Carrascal et al. 2009). The special case of PLS-DA is a supervised method where the response is a categorical variable (Pérez-Enciso and Tenenhaus 2003), like the seasonal roost pattern in the present study. This method extracts components (or latent variables) from predictor variables and indicates its importance in explaining the response variable. The results of the PLS-DA model include R2X and R2Y scores which indicate the explained variance and Q2Y score which indicates the predictive variance of model fit (Wu and Guo 2018). The performance of PLS-DA model was assessed via overall misclassification error rate and significance of R2Y and Q2Y by leave-one-out cross-validation and permutation testing (100 cycles) respectively (Thévenot et al. 2015; Rohart et al. 2017). The predictor variables were standardised before PLS-DA regression. We performed a one-way repeated measures ANOVA followed by paired t-tests to determine whether there was a significant difference of microclimate temperature between winter roosts and other roosts. We conducted the Mauchly's test of sphericity to test if the variances of each paired roost site were equal and assumption of sphericity was met. Prior to the analysis, Tmicro values were standardised to eliminate the variation among sampling times.

    We used R version 3.5.1 (R Core Team 2018) to conduct the statistical analyses. Prior to statistical analyses, we checked the data normality by Shapiro-Wilk normality tests and homogeneity of variance by F-tests. We used the car package (Fox and Weisberg 2019) for one-way repeated measures ANOVA, vcd package (Meyer et al. 2017) for Pearson chi-squared test and mosaic plot (Friendly 2013), ggtern package (Hamilton and Ferry 2018) for ternary diagram, cluster and factoextra packages (Kassambara et al. 2017; Maechler et al. 2019) for clustering analysis, and ropls and mixOmics packages (Thévenot et al. 2015; Rohart et al. 2017) for PLS-DA analysis.

    BT was the most frequently used roost (44.8% of total observation days), followed by CT (23.8%), PT (16.2%) and RT (15.2%); LP, CR, and GT were the least frequently used (Fig. 2a). Winter roosts CT and RT were usually occupied between November and March. Specifically, CT was mainly occupied in the first year, while RT was mostly used in the second year. Likewise, spring roost PT was primarily used in 2014/5 whereas GT was used more in 2015/6. PT was often used between March and June, but GT was only occupied between April and May. The occupancy of summer roost BT lasted for 8 months, from April to November. In contrast, summer roosts CR and LP were usually occupied between June and September, and on most occasions, they were used as the minor roosts for BT. At the time of seasonal roost switching, the roosts of the last season and next season were sometimes used interchangeably. This implied that the environmental factors affecting roost switching had been changing but were not stable, and cockatoos were sensitive to these changes by choosing the optimal roost site on a daily basis.

    Figure  2.  Bar plots of roosting days (a) and violin plots of roosting flock size (b) of each roost site. Different letters in (b) denote significant differences of roosting flock size between roost sites (Kruskal-Wallis test followed by pairwise Wilcoxon rank sum-tests; p adjusted by Benjamini-Hochberg method)

    In terms of roosting flock size, the differences among BT, CT, GT, PT and RT were not significant, with a mean flock size ranging from 63 to 74. However, the mean roosting flock size of CR and LP was 21 and 22, respectively, which were significantly smaller than those of other roosts (Kruskal-Wallis test, χ2=37.008, df=6, p < 0.001). This is because these sites were artificial perching structures with limited space. Violin plots showed the data distribution of flock size of each roost, where the wider sections of the violin plots represented a higher probability of that flock size and vice versa (Fig. 2b). With the exception of BT, the distribution of the data was similar within seasons, but different across seasons. For example, the flock size was around the median for spring roosts GT and PT, while the distribution of flock size of summer roosts CR and LP had two noticeable peaks. It indicated that variations in flock size were larger for roosts used in different seasons than those used in the same season.

    Land use types were significantly associated with roost sites (Pearson's chi-squared test, χ2=130.82, df=10, p < 0.001). As indicated by the Pearson residuals, tree-dominated areas were strongly positively associated with spring roosts GT and PT (Fig. 3a), with a proportion of 50% and 63% respectively (Fig. 3b). In contrast, tree-dominated areas were significantly negatively related to summer roost BT (18%) and LP (14%), and winter roost CT (15%). The proportion of tree-dominated areas separated the roosts used during spring from the other roosts (Fig. 3b). This indicated that spring roosts were closer to green patches than other roosts. Overall, the amount of built-up area varied considerably; for example, the percentage of built-up area was extremely low surrounding spring roost PT (3%) but was much larger for summer roost BT (44%).

    Figure  3.  A mosaic plot (a) with residual-based shading of land use types versus roost sites, and a ternary diagram (b) of the proportion of land use types within a 50-m circular plot of roost sites. In the mosaic plot, blue colour indicates the observed value is higher than the expected value, while red colour denotes the observed value is lower than the expected value

    The seasonal roost categories were well distinguished by the cluster analysis using anthropogenic disturbance parameters (Fig. 4a). This suggests that human disturbance was associated with seasonal switching among roosts. Moreover, two components (Comp 1 and Comp 2) extracted from the PLS-DA regression were sufficient to achieve the best classification performance of roosts by cross-validation (Overall error rate=0.000) and permutation testing (p=0.04). They captured 73.6% of the human disturbance variables and explained 98.4% of the variance of seasonal roost types (Table 1). Performance of the prediction was good with a cumulative predict variance (Q2Y) of 89.1%. Comp 1 discriminated spring roosts (PT and GT) from other roosts (Fig. 4b), where the distance to the nearest building, streetlight and road were positively correlated with spring roosts while traffic, noise level, pedestrians and light intensity had negative loadings on spring roosts (Fig. 4c, d). Consistent with the analyses on land use types, spring roosts were much closer to other trees. These indicated that spring roosts were distant from all human activities. Comp 2 further discriminated winter and summer roosts (Fig. 4b). Both roost groups were close to building, light and road. However, winter roosts were positively associated with roosting height, pedestrians, and noise level, where they were usually near bus stops, zebra crossings, and a tourist service centre. Summer roosts were positively related to light intensity (Fig. 4c, e).

    Figure  4.  Hierarchical dendrogram using Ward's clustering method (a), and PLS-DA regression score plot (b), correlation circle plot (c) and loading plots (d, e) of anthropogenic disturbance parameters on seasonal roost patterns. The score plot visualises the component scores, and prediction areas of each roost category. The correlation circle plot displays the Pearson correlation between the variables and latent components. The loading plots show the loading weights (contributions) of each variable on each component. The most important variables (according to coefficient values) are ordered from bottom to top. Colours indicate the seasonal roost pattern where the mean value of this variable is maximal. Roost Height: roosting height; Light: light intensity; Noise: noise level; DistRoad/DistTree/DistLight/DistBuild: distance to the nearest road/tree/light/building
    Table  1.  Percentage of explained variance in the predictor variables (R2X) and the response variable (R2Y) per component and their cumulative values (R2X cum and R2Y cum) derived from PLS-DA regression
    Component R2X R2X cum R2Y R2Y cum Q2Y cum Error rate
    Comp 1 0.491 0.491 0.491 0.491 0.165 0.833
    Comp 2 0.245 0.736 0.494 0.984 0.891 0.000
    Q2Y cum is the cumulative predictive variance of fit by cross-validation. Error rate using maximum distance in Comp 2 includes the performance of the PLS-DA model in Comp 1
     | Show Table
    DownLoad: CSV

    Roosts were categorised into two distinct groups by Tmicro, viz. the high and low Tmicro group, that exhibited a highly significant difference (one-way repeated measures ANOVA, sphericity assumed, F (5, 55)=37.6, p < 0.001; Fig. 5). The temperature difference between groups attained as much as 1.1 ℃, where the high Tmicro group included CT, RT and LP, and the low Tmicro group comprised of BT, GT and PT. Winter roosts were all categorised into the high Tmicro group, which were significantly warmer than the spring roosts and the major summer roost. It might be because winter roosts are close to the main road with heavy traffic and pedestrians (Additional file 1: Table S2), from which the external source of heat was produced.

    Figure  5.  Boxplots of standardised Tmicro of each roost site. Different letters denote significant differences of Tmicro between roost sites (one-way repeated measures ANOVA followed by paired pairwise t-tests; p adjusted by Benjamini-Hochberg method)

    Seasonal switching of roost sites was associated with multiple environmental factors, wherein some microhabitat characteristics showed consistency among roost sites used in the same season but significant differences among sites used in different seasons.

    Spring roosts were distinguished from other roosts by the high proportion of canopy cover at and adjacent to the roost sites. Besides, spring roosts were farther away from the main roads, streetlights, and buildings. This indicates that spring roosts were closer to green patches than other roosts. In addition, these related green patches are urban parks which serve as the daily feeding grounds and diurnal activity centre for Yellow-crested Cockatoo (Wang et al. 2018). A handful of research had shed light on the importance of food resources in roost site selection and shown that roost sites are usually located in the most visited foraging areas (Jenni 1991; Yap et al. 2002; Saiyad et al. 2017). Caccamise and Morrison (1988) found an increase in the number of European Starling using major roosts coincided with seasonal changes in foraging-substrate selection and diet composition. However, the information centre hypothesis may not apply to this study as the locations of food sources (e.g. plants) did not change dramatically in the short-term.

    Summer roosts were often located in less vegetated areas but surrounded by a larger proportion of built-up and open area. Besides, night light had a positive correlation on roosting in summer for Yellow-crested Cockatoo. Rainbow Lorikeet (Trichoglossus haematodus) and Carnaby's Black Cockatoo (Calyptorhynchus latirostris) had similar preference for roost sites in parking lots and sports grounds that received more artificial light at night (Berry 2008; Jaggard et al. 2015; Daoud-Opit and Jones 2016; Le Roux 2017). The roosts of House Crow (Corvus splendens) in mid-July to mid-October in metropolitan Singapore had better illumination, shorter distance to buildings, and less surrounding vegetation (Peh and Sodhi 2002), a result which was very similar to the urban setting of the summer roosts of Yellow-crested Cockatoo in this study. It is generally acknowledged that light is associated with predator avoidance because parrot species have relatively poor vision in low-light conditions (Birkhead 2012). New-borns of Yellow-crested Cockatoo were often observed in the roost sites. Indeed, roost switching in summer was probably related to the protection of juveniles because they were likely to be vulnerable to some commonly seen urban predators (e.g. wild cats, masked palm civets).

    Winter roosts were characterised by higher microhabitat temperatures, which suggests a preference for relatively warmer roost sites in cold weather. Similar studies found that the minimum temperature was usually 0.5 ℃ to 1.4 ℃ greater at roosts compared to non-roosts for American Crow (Corvus brachyrhynchos) and European Starling (Sturnus vulgaris) in an urban environment (Gorenzel and Salmon 1995; Clergeau and Quenot 2007; Erikson 2015). The mean temperature of Hong Kong is 6 ℃ to 9 ℃ lower than that in Indonesia between November to April, which is an indication of their behavioural flexibility that Yellow-crested Cockatoo coped with low-temperature stress in a non-native urban environment. Switches between roosts over a short period also highlight that Yellow-crested Cockatoo is sensitive to abrupt temperature changes. During late November to early December 2015, the cockatoos started to roost on the winter roost RT when the ambient temperature dropped to 15.9 ℃, but they returned quickly to the summer roost BT as soon as the ambient temperature increased to 22.0 ℃. After a few days, they switched back to the winter roost again when cold weather prevailed. The warm roosts were usually associated with intensive human activities because the concrete substrate beneath and urban heat island effect is an external source of heat for the roosts (Clergeau and Quenot 2007; Le Roux 2017). Due to urban heat island effect, the temperature in commercial areas is usually greater than residential areas and urban parks (Lee et al. 2009), which is consistent with land uses composition for roost sites in this study. For example, Asian Rose-necked Parakeet (Alexandrinus manillensis) roosted in the warmest places in cities where traffic was highest (Braun et al. 2018). American Crow preferred roosting in commercial areas and along streets, where the interior canopy temperatures of the roost trees were higher (Gorenzel and Salmon 1995). Likewise, the winter roosts of Yellow-crested Cockatoo were characterised by high levels of noise and pedestrian flow under the roost trees, as well as a short distance to the main road with heavy traffic.

    Our results confirmed that differences in microhabitat characteristics among seasonal roosts could be significant. Consequently, it was not surprising that some studies, focusing on variations between roosts and non-roosts, have obtained inconsistent results with ours. For example, Erickson (2015) found that the night-time light intensity was lower at roosts than that at non-roosts probably because they only investigated the winter roosts. Clergeau and Quenot (2007) found ambient temperature but not illumination significantly higher in the roost site, perhaps as they only studied the roosts from November to February. In contrast, our results indicated that illumination was possibly just not the major selective factor for birds during this period. Despite these inconsistencies with previous work, which may be the result of habitat differences at the species level, our study provides an alternative explanation for variation among roost sites between which birds can seasonally switch.

    This study was limited by the small sample size. Cockatoos (c. 20-30 individuals) outside the study area were not surveyed due to their small flock size and unstable roosting behaviour observed. Another limitation was that the results based on correlation analyses were not robust to support any hypothesis related to the benefits and driving forces behind communal roosting. Sampling size could be increased by identifying other relatively stable roosting flocks in Hong Kong or studying urban populations in other countries such as Singapore (Neo 2012). Besides, future studies can be designed to explain the specific requirements for roost sites in different seasons.

    Roost site switching is the result of the selection of different environmental factors in different seasons. Spring roosts were close to green patches, and distant from the main roads, streetlights and buildings. Summer roosts usually received more artificial light in less vegetated areas. Winter roosts had higher microhabitat temperatures, usually being isolated tall trees located in the most crowded and noisiest urban areas. The study provides evidence that the variations among roost sites were associated with differences in microhabitat characteristics and suggests that special attention should be paid to the seasonal variation among roost sites in future studies.

    The online version contains supplementary material available at https://doi.org/10.1186/s40657-021-00270-9.

    Additional file 1: Table S1. The difference of mean daily minimum temperature (Tmin) between Hong Kong (Tmin, HK) and five locations of Indonesia (− Tmin), respectively. Table S2. Mean values and standard deviations (if any) of anthropogenic disturbance variables of roost sites of Yellow-crested Cockatoo in urban areas of Hong Kong.

    Additional file 2: Figure S1. Boxplots of the Tambient (night minimum air temperature of occupied days) of each roost site during the study period. Different letters denote significant differences between roost sites (Kruskal–Wallis test followed by pairwise Wilcoxon rank sum tests; p adjusted by Benjamini–Hochberg method)

    We thank Mr. Ben Yeung and Mr. Banson Leung for their technical assistance.

    LMC and SW conceived the idea. SW collected and analysed the data, and leaded the writing. LMC revised and polished the manuscript. Both authors read and approved the final manuscript.

    This study was supported by the Hong Kong Environment and Conservation Fund (ECF Project 02/2014) and the United College Endowment Fund of The Chinese University of Hong Kong.

    The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

    Not applicable.

    Not applicable.

    The authors declare that they have no competing interests.

  • Collar NJ, Robson C (2007) Family Timaliidae (Babblers). In: del Hoyo J, Elliott A, Christie D (eds) Handbook of birds of the World. Lynx Edicions, Barcelona, pp 70-291
    Pfeifer R (2013) What is a babbler? Fascinating diversity of life strategies and challenge for systematists. Vogelwarte 51:117-126
    Alström P, Ericson PGP, Olsson U, Sundberg P (2006) Phylogeny and classification of the avian superfamily Sylvioidea. Mol Phylogenet Evol 38:381-397
    Alström P, Olsson U, Lei FM (2013) A review of the recent advances in the systematics of the avian superfamily Sylvioidea. Chin Birds 4:99-131
    Alström P, Hopper DM, Liu Y, Olsson U, Mohan D, Gelang M et al (2014) Discovery of a relict lineage and monotypic family of passerine birds. Biol Lett 10:20131067
    Cibois A (2003) Mitochondrial DNA phylogeny of babblers (Timaliidae). Auk 120:35-54
    Cibois A, Kalyakin MV, Han L-X, Pasquet E (2002) Molecular phylogenetics of babblers (Timaliidae): revaluation of the genera Yuhina and Stachyris. J Avian Biol 33:380-390
    Dong F, Li S-H, Yang X-J (2010) Molecular systematics and diversification of the Asian scimitar babblers (Timaliidae, Aves) based on mitochondrial and nuclear DNA sequences. Mol Phylogenet Evol 57:1268-1275
    Dong F, Wu F, Liu LM, Yang XJ (2010) Molecular phylogeny of the barwings (Aves: Timaliidae: Actinodura), a paraphyletic group, and its taxonomic implications. Zool Stud 49:703-709
    Gelang M, Cibois A, Pasquet E, Olsson U, Alström P, Ericson PGP (2009) Phylogeny of babblers (Aves, Passeriformes): major lineages, family limits and classification. Zool Scr 38:225-236
    Moyle RG, Andersen MJ, Oliveros CH, Steinheimer FD, Reddy S (2012) Phylogeny and biogeography of the core babblers (Aves: Timaliidae). Syst Biol 61:631-651
    Reddy S (2008) Systematics and biogeography of the shrike-babblers (Pteruthius): species limits, molecular phylogenetics, and diversification patterns across southern Asia. Mol Phylogenet Evol 47:54-72
    Reddy S, Cracraft J (2007) Old World Shrike-babblers (Pteruthius) belong with New World Vireos (Vireonidae). Mol Phylogenet Evol 44:1352-1357
    Oliveros CH, Reddy S, Moyle RG (2012) The phylogenetic position of some Philippine "babblers" spans the muscicapoid and sylvioid bird radiations. Mol Phylogenet Evol 65:799-804
    Beecher WJ (1953) A phylogeny of the Oscines. Auk 70:270-333
    Collar NJ (2006) A partial revision of the Asian babblers (Timaliidae). Forktail 22:85-112
    Delacour J (1946) Les timaliinés. L'Oiseau 16:7-36
    Mayr E, Amadon D (1951) A classification of recent birds. Am Mus Novit 1496:1-42
    Pasquet E, Bourdon E, Kalyakin MV, Cibois A (2006) The fulvettas (Alcippe, Timaliidae, Aves): a polyphyletic group. Zool Scr 35:559-566
    Barker FK, Cibois A, Schikler P, Feinstein J, Cracraft J (2004) Phylogeny and diversification of the largest avian radiation. Proc Natl Acad Sci USA 101:11040-11045
    Fregin S, Haase M, Olsson U, Alström P (2012) New insights into family relationships within the avian superfamily Sylvioidea (Passeriformes) based on seven molecular markers. BMC Evol Biol 12:157
    Gill F, Donsker D (eds) (2015) IOC World Bird List (v 5.2). doi:
    Collar NJ (2011) Taxonomic notes on some Asian babblers (Timaliidae). Forktail 27:100-102
    Dickinson EC, Christidis L (2014) The Howard and Moore complete checklist of the birds of the world, vol 2: Passerines, 4th edn. Christopher Helm, London
    Athreya R (2006) A new species of Liocichla (Aves: Timaliidae) from Eaglenest Wildlife Sanctuary, Arunachal Pradesh, India. Indian Birds 2:82-94
    Fu Y, Dowell SD, Zhang Z (2013) Emei Shan Liocichla: population, behavior and conservation. Chin Birds 4:260-264
    Peterson AT, Papeş M (2006) Potential geographic distribution of the Bugun Liocichla Liocichla bugunorum, a poorly-known species from north-eastern India. Indian Birds 2:146-149
    McKay BD, Mays HL Jr, Peng Y-W, Kozak KH, Yao C-T, Yuan H-W et al (2010) Recent range-wide demographic expansion in a Taiwan endemic montane bird, Steere's Liocichla (Liocichla steerii). BMC Evol Biol 10:71
    Dickinson EC (2003) The Howard and Moore complete checklist of the birds of the world, 3rd edn. Christopher Helm, London
    Luo X, Qu YH, Han LX, Li S-H, Lei FM (2009) A phylogenetic analysis of laughingthrushes (Timaliidae: Garrulax) and allies based on mitochondrial and nuclear DNA sequences. Zool Scr 38:9-22
    Päckert M, Martens J, Sun Y-H, Severinghaus LL, Nazarenko AA, Ting J et al (2011) Horizontal and elevational phylogeographic patterns of Himalayan and Southeast Asian forest passerines (Aves: Passeriformes). J Biogeogr 39:556-573
    de Bruyn M, Stelbrink B, Morley RJ, Hall R, Carvalho GR, Cannon CH et al (2014) Borneo and Indochina are major evolutionary hotspots for Southeast Asian biodiversity. Syst Biol 63:879-901
    Price TD, Hooper DM, Buchanan CD, Johansson US, Tietze DT, Alström P et al (2014) Niche filling slows the diversification of Himalayan songbirds. Nature 509:222-225
    Backström N, Fagerberg S, Ellegren H (2008) Genomics of natural bird populations: a gene-based set of reference markers evenly spread across the avian genome. Mol Ecol 17:964-980
    Primmer CR, Borge T, Lindell J, Sætre GP (2002) Single-nucleotide polymorphism characterization in species with limited available sequence information: high nucleotide diversity revealed in the avian genome. Mol Ecol 11:603-612
    Marini MÂ, Hackett SJ (2002) A multifaceted approach to the characterization of an intergeneric hybrid manakin (Pipridae) from Brazil. Auk 119:1114-1120
    Moonsamy PV, Williams T, Bonella P, Holcomb CL, Höglund BN, Hillman G et al (2013) High throughput HLA genotyping using 454 sequencing and the Fluidigm Access Array™ System for simplified amplicon library preparation. Tissue Antigens 81:141-149
    Sorenson MD, Ast JC, Dimcheff DE, Yuri T, Mindell DP (1999) Primers for a PCR-based approach to mitochondrial genome sequencing in birds and other vertebrates. Mol Phylogenet Evol 12(2):105-114
    Chesser RT (1999) Molecular systematics of the Rhinocryptid genus Pteroptochos. Condor 101:439-446
    Groth JG (1998) Molecular phylogenetics of finches and sparrows: consequences of character state removal in cytochrome b sequences. Mol Phylogenet Evol 10(3):377-390
    Reddy S, Moyle RG (2011) Systematics of the scimitar babblers (Pomatorhinus: Timaliidae): Phylogeny, biogeography and species-limits of four species complexes. Biol J Linn Soc 102:846-869
    Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S et al (2012) Geneious basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics 28:1647-1649
    Stephens M, Smith NJ, Donnelly P (2001) A new statistical method for haplotype reconstruction from population data. Am J Hum Genet 68:978-989
    Librado P, Rozas J (2009) DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics 25:1451-1452
    Bruen TC, Philippe H, Bryant D (2006) A simple and robust statistical test for detecting the presence of recombination. Genetics 172:2665-2681
    Huson DH, Bryant D (2006) Application of phylogenetic networks in evolutionary studies. Mol Biol Evol 23:254-267
    Hudson RR, Kreitman M, Aguadé M (1987) A test of neutral molecular evolution based on nucleotide data. Genetics 116:153-159
    Darriba D, Taboada GL, Ramón D, Posada D (2012) JModelTest 2: more models, new hueristics and parallel computing. Nat Methods 9:772
    Guindon S, Dufayard JF, Lefort V, Anisimova M, Hordijk W, Gascuel O (2010) New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. Syst Biol 59:307-321
    Guindon S, Gascuel O (2003) A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol 52:696-704
    Swofford DL (2003) PAUP* phylogenetic analysis using parsimony (*and Other Methods) v 4. Sinauer Associates, Sunderland
    Bouckaert R, Heled J, Kühnert D, Vaughan T, Wu C-H, Xie D et al (2014) BEAST 2: a software platform for Bayesian evolutionary analysis. PLoS Comput Biol 10:e1003537
    Heled J, Drummond AJ (2010) Bayesian inference of species trees from multilocus data. Mol Biol Evol 27:570-580
    Bouckaert R (2015) bModelTest: Bayesian site model selection for nucleotide data. bioRxiv. . Accessed 28 July 2015
    Drummond AJ, Ho SYW, Phillips MJ, Rambaut A (2006) Relaxed phylogenetics and dating with confidence. PLoS Biol 4:e88
    Weir JT, Schluter D (2008) Calibrating the avian molecular clock. Mol Ecol 17:2321-2328
    Pereira SL, Baker AJ (2006) A mitogenomic timescale for birds detects variable phylogenetic molecular evolution and refutes the standard molecular clock. Mol Biol Evo 23:1731-1740
    Fleischer RC, McIntosh CE, Tarr CL (1998) Evolution on a volcanic conveyor belt: using phylogeographic reconstructions and K-Ar-based ages of the Hawaiian Islands to estimate molecular evolutionary rates. Mol Ecol 7:533-545
    Rambaut A, Suchard M, Xie D, Drummond AJ (2014) Tracer v 1.6. . Accessed 28 July 2015
    Baele G, Lemey P, Bedford T, Rambaut A, Suchard MA, Alekseyenko AV (2012) Improving the accuracy of demographic and molecular clock model comparison while accommodating phylogenetic uncertainty. Mol Biol Evol 29:2157-2167
    Zink RM, Barrowclough GF (2008) Mitochondrial DNA under siege in avian phylogeography. Mol Ecol 17(9):2107-2121
    Mays HL Jr, Doucet SM, Yao C-T, Yuan H-W (2006) Sexual dimorphism and dichromatism in Steere's Liocichla (Liocichla steerii). J Field Ornithol 77:437-443
    Dong F, Li S-H, Zou FS, Lei FM, Liang W, Yang JX et al (2014) Molecular systematics and plumage coloration evolution of an enigmatic babbler (Pomatorhinus ruficollis) in East Asia. Mol Phylogenet Evol 70:76-83
    McKay BD, Mays HL Jr, Wu YC, Li H, Yao C-T, Nishiumi I et al (2013) An empirical comparison of character-based and coalescent-based approaches to species delimitation in a young avian complex. Mol Ecol 22:4943-4957
    Reddy S, Sharief S, Yohe LR, Witkowski J, Hosner PA, Nyári ÁS et al (2015) Untangling taxonomic confusion and diversification patterns of the Streak-breasted Scimitar Babblers (Timaliidae: Pomatorhinus ruficollis complex) in southern Asia. Mol Phylogenet Evol 82:183-192
    Stoddard MC, Prum RO (2011) How colorful are birds? Evolution of the avian plumage color gamut. Behav Ecol 22:1042-1052
    Friedman NR, McGraw KJ, Omland KE (2014) Evolution of carotenoid pigmentation in caciques and meadowlarks (Icteridae): repeated gains of red plumage coloration by carotenoid C4-oxygenation. Evolution 68:791-801
    Omland KE, Lanyon SM (2000) Reconstructing plumage evolution in orioles (Icterus): repeated convergence and reversal in patterns. Evolution 54:2119-2133
    Molnar P (2005) Mio-Pliocene growth of the Tibetan plateau and evolution of East Asian climate. Palaeontol Electron 8:2A
    Zhang R, Jiang DB, Zhang ZS, Yu E (2015) The impact of regional uplift of the Tibetan Plateau on the Asian monsoon climate. Palaeogeogr Palaeoclimatol Palaeoecol 417:137-150
    Bird MI, Taylor D, Hunt C (2005) Palaeoenvironments of insular Southeast Asia during the Last Glacial Period: a savanna corridor in Sundaland? Quat Sci Rev 24:2228-2242
    Cerling TE, Harris JM, MacFadden BJ, Leakey MG, Quade J, Eisenmann V et al (1997) Global vegetation change through the Miocene/Pliocene boundary. Nature 389:153-158
    McKay BD (2011) A new timeframe for the diversification of Japan's mammals. J Biogeogr 39:1134-1143
    Johansson US, Alström P, Olsson U, Ericson PGP, Sundberg P, Price TD (2007) Build-up of the Himalayan avifauna through immigration: a biogeographical analysis of the Phylloscopus and Seicercus warblers. Evolution 61:324-333
    Qu YH, Song G, Gao B, Quan Q, Ericson PGP, Lei FM (2014) The influence of geological events on the endemism of East Asian birds studied through comparative phylogeography. J Biogeogr 42:179-192
    Che J, Zhou WW, Hu JS, Yan F, Papenfuss TJ, Wake DB et al (2010) Spiny frogs (Paini) illuminate the history of the Himalayan region and Southeast Asia. Proc Natl Acad Sci USA 107:13765-13770
    Qiu YX, Guan BC, Fu CX, Comes HP (2009) Did glacials and/or interglacials promote allopatric incipient speciation in East Asian temperate plants? Phylogeographic and coalescent analysis on refugial isolation and divergence in Dysosma versipellis. Mol Phylogenet Evol 51:281-293
    Wang L, Schneider H, Zhang XC, Xiang QP (2012) The rise of the Himalaya enforced the diversification of SE Asian ferns by altering the monsoon regimes. BMC Plant Biol 12:210
    Wang WJ, McKay BD, Dai CY, Zhao N, Zhang RY, Qu YH et al (2013) Glacial expansion and diversification of an East Asian montane bird, the green-backed tit (Parus monticolus). J Biogeogr 40:1156-1169
    Song G, Qu YH, Yin ZH, Li S-H, Liu N, Lei FM (2009) Phylogeography of the Alcippe morrisonia (Aves: Timaliidae): long population history beyond late Pleistocene glaciations. BMC Evol Biol 9:143
    Zou FS, Lim H-C, Marks BD, Moyle RG, Sheldon FH (2007) Molecular phylogenetic analysis of the Grey-cheeked Fulvetta (Alcippe morrisonia) of China and Indochina: a case of remarkable genetic divergence in a "species". Mol Phylogenet Evol 44:165-174
    Lovette IJ (2004) Mitochondrial dating and mixed support for the "2% Rule" in birds. Auk 121:1-6
    Bromham L, Penny D (2003) The modern molecular clock. Nat Rev Genet 4(3):216-224
    Welch JJ, Bromham L (2005) Molecular dating when rates vary. Trends Ecol Evol 20:320-327
    Teng LS (1990) Geotectonic evolution of late Cenozoic arc-continent collision in Taiwan. Tectonophysics 183:57-76
    Huang C-Y, Yuan PB, Tsao S-J (2006) Temporal and spatial records of active arc-continent collision in Taiwan: a synthesis. Geol Soc Am Bull 118:274-288
    Lee Y-H, Chen C-C, Liu T-K, Ho H-C, Lu H-Y, Lo W (2006) Mountain building mechanisms in the Southern Central Range of the Taiwan Orogenic Belt—From accretionary wedge deformation to arc—continental collision. Earth Planet Sci Lett 252:413-422
    Päckert M, Martens J, Severinghaus LL (2008) The Taiwan Firecrest (Regulus goodfellowi) belongs to the Goldcrest assemblage (Regulus regulus s. l.): evidence from mitochondrial DNA and the territorial song of the Regulidae. J Ornithol 150:205-220
    Wu HC, Lin RC, Hung H-Y, Yeh C-F, Chu J-H, Yang XJ et al (2011) Molecular and morphological evidences reveal a cryptic species in the Vinaceous rosefinch Carpodacus vinaceus (Fringillidae; Aves). Zool Scr 40:468-478
    Tietze DT, Päckert M, Martens J, Lehmann H, Sun Y-H (2013) Complete phylogeny and historical biogeography of true rosefinches (Aves: Carpodacus). Zool J Linn Soc 169:215-234
    McKay BD, Mays HL Jr, Yao C-T, Wan D, Higuchi H, Nishiumi I (2014) Incorporating color into integrative taxonomy: analysis of the varied tit (Sittiparus varius) complex in East Asia. Syst Biol 63:505-517
    Chou Y-W, Thomas PI, Ge X-J, LePage BA, Wang C-N (2011) Refugia and phylogeography of Taiwania in East Asia. J Biogeogr 38:1992-2005
    Jang-Liaw N-H, Lee T-H, Chou W-H (2008) Phylogeography of Sylvirana latouchii (Anura, Ranidae) in Taiwan. Zool Sci 25:68-79
    Jang-Liaw N-H, Chou W-H (2011) Phylogeography of the fanged dicroglossine frog, Limnonectes fujianensis (Anura, Ranidae), in Taiwan. Zool Sci 28:254-263
    Huang J-P, Lin C-P (2011) Lineage-specific late pleistocene expansion of an endemic subtropical gossamer-wing damselfly, Euphaea formosa, in Taiwan. BMC Evol Biol 11:94
    Lei FM, Qu YH, Lu JL, Liu Y, Yin ZH (2003) Conservation on diversity and distribution patterns of endemic birds in China. Biodivers Conserv 12:239-254
    Favre A, Päckert M, Pauls SU, Jähnig SC, Uhl D, Michalak I et al (2014) The role of the uplift of the Qinghai-Tibetan Plateau for the evolution of Tibetn biotas. Biol Rev Camb Philos Soc 90:236-253
    Li Y, Zhai SN, Qiu YX, Guo YP, Ge XJ, Comes HP (2011) Glacial survival east and west of the 'Mekong-Salween Divide' in the Himalaya-Hengduan Mountains region as revealed by AFLPs and cpDNA sequence variation in Sinopodophyllum hexandrum (Berberidaceae). Mol Phylogenet Evol 59:412-424
    Kingdon-Ward F (1921) The Mekong-Salween divide as a geographical barrier. Geogr J 58:49-56
  • Related Articles

  • Cited by

    Periodical cited type(4)

    1. Sajeewa S. N. Maharachchikumbura, Yanpeng Chen, Hiran A. Ariyawansa, et al. Integrative approaches for species delimitation in Ascomycota. Fungal Diversity, 2021, 109(1): 155. DOI:10.1007/s13225-021-00486-6
    2. Alice Cibois, Magnus Gelang, Per Alström, et al. Comprehensive phylogeny of the laughingthrushes and allies (Aves, Leiothrichidae) and a proposal for a revised taxonomy. Zoologica Scripta, 2018, 47(4): 428. DOI:10.1111/zsc.12296
    3. Jiekun He, Zefang Gao, Yuanyuan Su, et al. Geographical and temporal origins of terrestrial vertebrates endemic to Taiwan. Journal of Biogeography, 2018, 45(11): 2458. DOI:10.1111/jbi.13438
    4. Pengcheng Wang, Hongyan Yao, Kadeem J. Gilbert, et al. Glaciation-based isolation contributed to speciation in a Palearctic alpine biodiversity hotspot: Evidence from endemic species. Molecular Phylogenetics and Evolution, 2018, 129: 315. DOI:10.1016/j.ympev.2018.09.006

    Other cited types(0)

Catalog

    Figures(5)  /  Tables(3)

    Article Metrics

    Article views (439) PDF downloads (11) Cited by(4)

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return