Parallel genomic responses to historical climate change and high elevation in East Asian songbirds.

Parallel evolution can be expected among closely related taxa exposed to similar selective pressures. However, parallelism is typically stronger at the phenotypic level, while genetic solutions to achieve these phenotypic similarities may differ. For polygenic traits, the availability of standing genetic variation (i.e., heterozygosity) may influence such genetic nonparallelism. Here, we examine the extent to which high-elevation adaptation is parallel-and whether the level of parallelism is affected by heterozygosity-by analyzing genomes of 19 Paridae species distributed across East Asia with a dramatic east-west elevation gradient. We find that western highlands endemic parids have consistently lower levels of heterozygosity-likely the result of late-Pleistocene demographic contraction-than do parids found exclusively in eastern lowlands, which remained unglaciated during the late Pleistocene. Three widespread species (east to west) have high levels of heterozygosity similar to that observed in eastern species, although their western populations are less variable than eastern ones. Comparing genomic responses to extreme environments of the Qinghai-Tibet Plateau, we find that the most differentiated genomic regions between each high-elevation taxon and its low-elevation relative are significantly enriched for genes potentially related to the oxygen transport cascade and/or thermogenesis. Despite no parallelism at particular genes, high similarity in gene function is found among comparisons. Furthermore, parallelism is not higher in more heterozygous widespread parids than in highland endemics. Thus, in East Asian parids, parallel functional response to extreme elevation appears to rely on different genes, with differences in heterozygosity having no effect on the degree of genetic parallelism.

Genomic Consequences of Long-Term Population Decline in Brown Eared Pheasant.

Population genetic theory and empirical evidence indicate that deleterious alleles can be purged in small populations. However, this viewpoint remains controversial. It is unclear whether natural selection is powerful enough to purge deleterious mutations when wild populations continue to decline. Pheasants are terrestrial birds facing a long-term risk of extinction as a result of anthropogenic perturbations and exploitation. Nevertheless, there are scant genomics resources available for conservation management and planning. Here, we analyzed comparative population genomic data for the three extant isolated populations of Brown eared pheasant (Crossoptilon mantchuricum) in China. We showed that C. mantchuricum has low genome-wide diversity and a contracting effective population size because of persistent declines over the past 100,000 years. We compared genome-wide variation in C. mantchuricum with that of its closely related sister species, the Blue eared pheasant (C. auritum) for which the conservation concern is low. There were detrimental genetic consequences across all C. mantchuricum genomes including extended runs of homozygous sequences, slow rates of linkage disequilibrium decay, excessive loss-of-function mutations, and loss of adaptive genetic diversity at the major histocompatibility complex region. To the best of our knowledge, this study is the first to perform a comprehensive conservation genomic analysis on this threatened pheasant species. Moreover, we demonstrated that natural selection may not suffice to purge deleterious mutations in wild populations undergoing long-term decline. The findings of this study could facilitate conservation planning for threatened species and help recover their population size.

Not out of the woods yet: Signatures of the prolonged negative genetic consequences of a population bottleneck in a rapidly re-expanding wader, the black-faced spoonbill Platalea minor.

The long-term persistence of a population which has suffered a bottleneck partly depends on how historical demographic dynamics impacted its genetic diversity and the accumulation of deleterious mutations. Here we provide genomic evidence for the genetic effect of a recent population bottleneck in the endangered black-faced spoonbill (Platalea minor) after its rapid population recovery. Our data suggest that the bird's effective population size, Ne , had been relatively stable (7500-9000) since 22,000 years ago; however, a recent brief yet severe bottleneck (Ne  = 20) which we here estimated to occur around the 1940s wiped out >99% of its historical Ne in roughly three generations. Despite a >15-fold population recovery since 1988, we found that black-faced spoonbill population has higher levels of inbreeding (7.4 times more runs of homozygosity) than its sister species, the royal spoonbill (P. regia), which is not thought to have undergone a marked population contraction. Although the two spoonbills have similar levels of genome-wide genetic diversity, our results suggest that selection on more genes was relaxed in the black-faced spoonbill; moreover individual black-faced spoonbills carry more putatively deleterious mutations (Grantham's score > 50), and may therefore express more deleterious phenotypic effects than royal spoonbills. Here we demonstrate the value of using genomic indices to monitor levels of genetic erosion, inbreeding and mutation load in species with conservation concerns. To mitigate the prolonged negative genetic effect of a population bottleneck, we recommend that all possible measures should be employed to maintain population growth of a threatened species.

Potential millennial-scale avian declines by humans in southern China.

Mounting observational records demonstrate human-caused faunal decline in recent decades, while accumulating archaeological evidence suggests an early biodiversity impact of human activities during the Holocene. A fundamental question arises concerning whether modern wildlife population declines began during early human disturbance. Here, we performed a population genomic analysis of six common forest birds in East Asia to address this question. For five of them, demographic history inference based on 25-33 genomes of each species revealed dramatic population declines by 4- to 48-fold over millennia (e.g. 2000-5000 thousand years ago). Nevertheless, summary statistics detected nonsignificant correlations between these population size trajectories and Holocene temperature variations, and ecological niche models explicitly predicted extensive range persistence during the Holocene, implying limited demographic consequence of Holocene climate change. Further analyses suggest high negative correlations between the reconstructed population declines and human disturbance intensities and indicate a potential driver of human activities. These findings provide a deep-time and large-scale insight into the recently recognized avifaunal decline and support an early origin hypothesis of human effects on biodiversity. Overall, our study sheds light on the current biodiversity crisis in the context of long-term human-environment interactions and offers a multi-evidential framework for quantitatively assessing the ecological consequences of human disturbance.

Standing genetic variation as the predominant source for adaptation of a songbird.

What kind of genetic variation contributes the most to adaptation is a fundamental question in evolutionary biology. By resequencing genomes of 80 individuals, we inferred the origin of genomic variants associated with a complex adaptive syndrome involving multiple quantitative traits, namely, adaptation between high and low altitudes, in the vinous-throated parrotbill (Sinosuthora webbiana) in Taiwan. By comparing these variants with those in the Asian mainland population, we revealed standing variation in 24 noncoding genomic regions to be the predominant genetic source of adaptation. Parrotbills at both high and low altitudes exhibited signatures of recent selection, suggesting that not only the front but also the trailing edges of postglacial expanding populations could be subjected to environmental stresses. This study verifies and quantifies the importance of standing variation in adaptation in a cohort of genes, illustrating that the evolutionary potential of a population depends significantly on its preexisting genetic diversity. These findings provide important context for understanding adaptation and conservation of species in the Anthropocene.

Strict allopatric speciation of sky island Pyrrhula erythaca species complex.

Increasing evidence of post-divergence gene flow between taxa is shifting our understanding on the mode of speciation. A fundamental question arises concerning the circumstances under which strict allopatric speciation still holds true. Sky island populations might undergo reduced gene flow by niche conservatism to highland habitats and follow divergence in an allopatric manner. In this study, we tested this hypothesis in the sky island Grey-headed Bullfinch (Pyrrhula erythaca) species complex via statistical analyses of both genetic and ecological data. Results of coalescent-based analysis of multiple nuclear loci suggested that P. e. owstoni likely colonized Taiwan island during the severe mid-Pleistocene glacial climate followed by strictly allopatric divergence from P. e. erythaca distributed in Himalayas-Hengduan mountains and central North China. Results of ecological niche modeling suggested that their speciation may be attributed to the niche conservatism of these birds and the lack of a suitable ecological corridor during subsequent milder glacial episodes. In addition, we delimited the traditionally defined P. erythaca into two full species, P. erythaca in the Asian mainland and P. owstoni on the island of Taiwan, based on both genetic and behavioural evidences. These results suggest that ecology can have a dynamic role in allowing highland populations to expand their ranges and isolated by habitat barriers to diversify in a strictly allopatric manner.

Disentangling the evolutionary history and biogeography of hill partridges (Phasianidae, Arborophila) from low coverage shotgun sequences.

The advent of the phylogenomic era has significantly improved our understanding of the evolutionary history and biogeography of Southeast Asia's diverse avian fauna. However, the taxonomy and phylogenetic relationships of many Southeast Asian birds remain poorly resolved, especially for those with large geographic ranges, which might have experienced both ancient and recent geological and environmental changes. In this study, we examined the evolutionary history and biogeography of the hill partridges (Galliformes: Phasianidae: Arborophila spp.), currently the second most speciose galliform genus, and thought to have colonized Southeast Asia from Africa. We present a well-resolved phylogeny of 14 Arborophila species inferred from ultra-conserved elements, exons, and mitochondrial genomes from both fresh and museum samples, which representing almost complete coverage of the genus. Our fossil-calibrated divergence time estimates and biogeographic modeling showed the ancestor of Arborophila arrived in Indochina during the early Miocene, but the initial divergence within Arborophila did not occur until ~10 Ma when global cooling intensified. Subsequent dispersal and diversification within Arborophila were driven by several tectonic and climatic events. In particular, we found evidence of rapid radiation in Indochinese Arborophila during the Pliocene global cooling and extensive dispersal and speciation of Sundaic Arborophila during the Pleistocene sea-level fluctuations. Taken together, these results suggest that the evolutionary history and biogeography of Arborophila were influenced by complex interactions among historical, geological and climatic events in Southeast Asia.

From the Himalayas to a continental Island: Integrative species delimitation in the Brownish-flanked Bush Warbler Horornis fortipes complex.

As species serve as basic units of study in many fields of biology, assessments of species limits are fundamental for such studies. Here, we used a multilocus dataset and different coalescent-based methods to analyze species delimitation and phylogenetic relationships in the Brownish-flanked Bush Warbler Horornis fortipes complex, which is widespread in the Sino-Himalayan region. We also examined the vocal and morphometric divergence within this complex. Our genetic results suggested that Horornis fortipes is composed of at least three independently evolving lineages, which diverged 1.1-1.8 million years ago. However, these lineages have hardly diverged in song or morphometrics and only very slightly in plumage. Our result indicate that there are three incipient species in Horonis fortipes complex diverged in central Himalayas and Hengduan Mountains, but not between the continent and Taiwan island.

Ice age unfrozen: severe effect of the last interglacial, not glacial, climate change on East Asian avifauna.

BACKGROUND: The glacial-interglacial cycles in the Pleistocene caused repeated range expansion and contraction of species in several regions in the world. However, it remains uncertain whether such climate oscillations had similar impact on East Asian biota, despite its widely recognized importance in global biodiversity. Here we use both molecular and ecological niche profiles on 11 East Asian avian species with various elevational ranges to reveal their response to the late Pleistocene climate changes. RESULTS: The ecological niche models (ENM) consistently showed that these avian species might substantially contract their ranges to the south during the Last Interglacial period (LIG) and expanded their northern range margins through the Last Glacial Maximum (LGM), leading to the LGM ranges observed for all 11 species. Consistently, coalescent simulations based on 25-30 nuclear genes retrieved signatures of significant population growth through the last glacial period across all species studied. Climate statistics suggested that high climatic variability during the LIG and a relatively mild climate at the LGM potentially explained the historical population dynamics of these birds. CONCLUSIONS: This is the first study based on multiple species and both lines of ecological niche profiles and genetic data to characterize the unique response of East Asian biota to late Pleistocene climate. The present study highlights regional differences in the evolutionary consequence of climate change during the last glacial cycle and implies that global warming might pose a great risk to species in this region given potentially higher climatic variation in the future analogous to that during the LIG.

Drastic population fluctuations explain the rapid extinction of the passenger pigeon.

To assess the role of human disturbances in species' extinction requires an understanding of the species population history before human impact. The passenger pigeon was once the most abundant bird in the world, with a population size estimated at 3-5 billion in the 1800s; its abrupt extinction in 1914 raises the question of how such an abundant bird could have been driven to extinction in mere decades. Although human exploitation is often blamed, the role of natural population dynamics in the passenger pigeon's extinction remains unexplored. Applying high-throughput sequencing technologies to obtain sequences from most of the genome, we calculated that the passenger pigeon's effective population size throughout the last million years was persistently about 1/10,000 of the 1800's estimated number of individuals, a ratio 1,000-times lower than typically found. This result suggests that the passenger pigeon was not always super abundant but experienced dramatic population fluctuations, resembling those of an "outbreak" species. Ecological niche models supported inference of drastic changes in the extent of its breeding range over the last glacial-interglacial cycle. An estimate of acorn-based carrying capacity during the past 21,000 y showed great year-to-year variations. Based on our results, we hypothesize that ecological conditions that dramatically reduced population size under natural conditions could have interacted with human exploitation in causing the passenger pigeon's rapid demise. Our study illustrates that even species as abundant as the passenger pigeon can be vulnerable to human threats if they are subject to dramatic population fluctuations, and provides a new perspective on the greatest human-caused extinction in recorded history.

Incipient speciation with gene flow on a continental island: Species delimitation of the Hainan Hwamei (Leucodioptron canorum owstoni, Passeriformes, Aves).

Because of their isolation, continental islands (e.g., Madagascar) are often thought of as ideal systems to study allopatric speciation. However, many such islands have been connected intermittently to their neighboring continent during recent periods of glaciation, which may cause frequent contact between the diverging populations on the island and continent. As a result, the speciation processes on continental islands may not meet the prerequisites for strictly allopatric speciation. We used multiple lines of evidence to re-evaluate the taxonomic status of the Hainan Hwamei (Leucodioptron canorum owstoni), which is endemic to Hainan, the largest continental island in the South China Sea. Our analysis of mitochondrial DNA and twelve nuclear loci suggests that the Hainan Hwamei can be regarded as an independent species (L. owstoni); the morphological traits of the Hainan Hwamei also showed significant divergence from those of their mainland sister taxon, the Chinese Hwamei (L. canorum). We also inferred the divergence history of the Hainan and Chinese Hwamei to see whether their divergence was consistent with a strictly allopatric model. Our results suggest that the two Hwameis split only 0.2 million years ago with limited asymmetrical post-divergence gene flow. This implies that the Hainan Hwamei is an incipient species and that speciation occurred through ecologically divergent selection and/or assortative mating rather than a strictly allopatric process.

Within-island speciation with an exceptional case of distinct separation between two sibling lizard species divided by a narrow stream.

Delimitation of genetic and geographic boundaries between species is a focus of evolutionary biology. In this study, we demonstrated fine-scale differentiation of Takydromus formosanus species complex comprising four insular endemics on Taiwan Island. Phylogeny and ancestral range reconstruction based on mitochondrial DNA sequences of 430 Takydromus lizards (405 lizards of this complex throughout their distribution range, and 25 lizards from 11 other species) indicated that the major branching process occurred within Taiwan, which represented a solid evidence of within-island speciation on this small island. We further demonstrated an exceptional case of a pair of sister species, T. viridipunctatus and T. luyeanus, that were separated by a narrow stream with a width of only 15m. This pattern might be one of the narrowest contact zones ever documented in terrestrial vertebrates. To evaluate the level of genetic introgression between these sister species, a fine-scale collection of another 382 lizards was conducted along a transect line across the stream. A total of 13 microsatellite markers and mtDNA genotyping was used to detect a low proportion of hybrids (5.7-9.9% from STRUCTURE, and 2.3% from DAPC). Our results indicated that the two clades are highly differentiated across this extremely short distance.

Was the exposed continental shelf a long-distance colonization route in the ice age? The Southeast Asia origin of Hainan and Taiwan partridges.

Research on island biotas has greatly contributed to the development of modern evolutionary and biogeographic theories. Until now, most studies have suggested that continental islands received their biotas directly from the adjacent mainland. However, only a few studies have indicated that species on continental islands might originate from other distantly non-adjacent regions. Here, we used the hill partridges (genus Arborophila) that are widely distributed in the southwest and southeast China mainland, Indochina, Hainan and Taiwan islands to test whether species on continental islands might originate from distant regions rather than the adjacent mainland. Based on molecular phylogenies inferred from three mitochondrial fragments and three nuclear introns, together with ancestral area reconstruction, we found that the ancestors of the endemic Hainan and Taiwan partridges (A. ardens and A. crudigularis) likely originated from Indochina, rather than the nearby southeast China mainland. The divergence time estimates demonstrate that their ancestors likely colonized Hainan and Taiwan islands using the long exposed continental shelf between Indochina, Hainan and Taiwan islands during glacial periods, which had not been demonstrated before. Thus, integrating distribution data with phylogenetic information can shed new lights on the historical biogeography of continental islands and surrounding mainland regions.

Climate niche differentiation between two passerines despite ongoing gene flow.

Niche evolution underpins the generation and maintenance of biological diversity, but niche conservatism, in which niches remain little changed over time in closely related taxa, and the role of ecology in niche evolution are continually debated. To test whether climate niches are conserved in two closely related passerines in East Asia - the vinous-throated (Paradoxornis webbianus) and ashy-throated (P. alphonsianus) parrotbills - we established their potential allopatric and sympatric regions using ecological niche models and compared differences in their climate niches using niche overlap indices in background tests and multivariate statistical analyses. We also used polymorphism data on 44 nuclear genes to infer their divergence demography. We found that these two parrotbills occupy different climate niches, in both their allopatric and potential sympatric regions. Because the potential sympatric region is the area predicted to be suitable for both parrotbills based on the ecological niche models, it can serve as a natural common garden. Therefore, their observed niche differences in this potential sympatry were not simply rendered by phenotypic plasticity and probably had a genetic basis. Our genetic analyses revealed that the two parrotbills are not evolutionarily independent for the most recent part of their divergence history. The two parrotbills diverged c. 856,000 years ago and have had substantial gene flow since a presumed secondary contact c. 290,000 years ago. This study provides an empirical case demonstrating that climate niches may not be homogenized in nascent species in spite of substantial, ongoing gene flow, which in turn suggests a role for ecology in promoting and maintaining diversification among incipient species.

Inferring the geographic mode of speciation by contrasting autosomal and sex-linked genetic diversity.

When geographic isolation drives speciation, concurrent termination of gene flow among genomic regions will occur immediately after the formation of the barrier between diverging populations. Alternatively, if speciation is driven by ecologically divergent selection, gene flow of selectively neutral genomic regions may go on between diverging populations until the completion of reproductive isolation. It may also lead to an unsynchronized termination of gene flow between genomic regions with different roles in the speciation process. Here, we developed a novel Approximate Bayesian Computation pipeline to infer the geographic mode of speciation by testing for a lack of postdivergence gene flow and a concurrent termination of gene flow in autosomal and sex-linked markers jointly. We applied this approach to infer the geographic mode of speciation for two allopatric highland rosefinches, the vinaceous rosefinch Carpodacus vinaceus and the Taiwan rosefinch C. formosanus from DNA polymorphisms of both autosomal and Z-linked loci. Our results suggest that the two rosefinch species diverged allopatrically approximately 0.5 Ma. Our approach allowed us further to infer that female effective population sizes are about five times larger than those of males, an estimate potentially useful when comparing the intensity of sexual selection across species.

Testing hypotheses of mitochondrial gene-tree paraphyly: unravelling mitochondrial capture of the Streak-breasted Scimitar Babbler (Pomatorhinus ruficollis) by the Taiwan Scimitar Babbler (Pomatorhinus musicus).

Species-level paraphyly inferred from mitochondrial gene trees is a prevalent phenomenon in taxonomy and systematics, but there are several potential causes that are not easily explained by currently used methods. This study aimed to test the underlying causes behind the observed paraphyly of Streak-breasted Scimitar Babbler (Pomatorhinus ruficollis) via statistical analyses of four mitochondrial (mtDNA) and nine nuclear (nuDNA) genes. Mitochondrial gene trees show paraphyly of P. ruficollis with respect to the Taiwan Scimitar Babbler (Pomatorhinus musicus), but nuclear genealogies support a sister-group relationship. Predictive coalescent simulations imply several hypothetical explanations, the most likely being mitochondrial capture of P. ruficollis by P. musicus for the observed cyto-nuclear incongruence. Further approximate Bayesian computation suggests a unidirectional introgression model with substantial level of gene flow from P. ruficollis to P. musicus during their initial divergence during the Late Pleistocene. This specific observation frames several potential causes for incongruent outcomes of mitochondrial and nuclear introgression in general, and on the whole, our results underscore the strength of multiple independent loci for species delimitation and importance of testing hypotheses that explain disparate causes of mitochondrial gene-tree paraphyly.

Influence of gene flow on divergence dating - implications for the speciation history of Takydromus grass lizards.

Dating the time of divergence and understanding speciation processes are central to the study of the evolutionary history of organisms but are notoriously difficult. The difficulty is largely rooted in variations in the ancestral population size or in the genealogy variation across loci. To depict the speciation processes and divergence histories of three monophyletic Takydromus species endemic to Taiwan, we sequenced 20 nuclear loci and combined with one mitochondrial locus published in GenBank. They were analysed by a multispecies coalescent approach within a Bayesian framework. Divergence dating based on the gene tree approach showed high variation among loci, and the divergence was estimated at an earlier date than when derived by the species-tree approach. To test whether variations in the ancestral population size accounted for the majority of this variation, we conducted computer inferences using isolation-with-migration (IM) and approximate Bayesian computation (ABC) frameworks. The results revealed that gene flow during the early stage of speciation was strongly favoured over the isolation model, and the initiation of the speciation process was far earlier than the dates estimated by gene- and species-based divergence dating. Due to their limited dispersal ability, it is suggested that geographical isolation may have played a major role in the divergence of these Takydromus species. Nevertheless, this study reveals a more complex situation and demonstrates that gene flow during the speciation process cannot be overlooked and may have a great impact on divergence dating. By using multilocus data and incorporating Bayesian coalescence approaches, we provide a more biologically realistic framework for delineating the divergence history of Takydromus.

Molecular systematics and plumage coloration evolution of an enigmatic babbler (Pomatorhinus ruficollis) in East Asia.

The streak-breasted scimitar babbler, Pomatorhinus ruficollis, is a polytypic and taxonomically enigmatic babbler common in southern, eastern, and southeastern Asia. To infer the phylogeny of the P. ruficollis, we examined the sequences of two complete mitochondrial genes (2184 bp in total) from fourteen of the fifteen known subspecies, and an additional five nuclear genes (2657 bp in total) from ten subspecies. The mitochondrial phylogeny indicates four major clades with large geographical identity in P. ruficollis and paraphyly of the P. ruficollis species complex, with the inclusion of the olivaceus group of congeneric P. schisticeps. Together with their interbreeding in northern Indochina, we propose to lump this group into P. ruficollis. Analysis of both multilocus networks and species-tree inference recovered poor phylogenetic structure among mainland/ Hainan subspecies and exclusive groupings of the Taiwanese subspecies, consistent with the recent taxonomic revision of its species status. Our analyses also suggest strong incongruence between the morphological-based classification and molecular systematics, implying the strength of multilocus data for taxonomy.

Past hybridization between two East Asian long-tailed tits (Aegithalos bonvaloti and A. fuliginosus).

INTRODUCTION: Incomplete lineage sorting and hybridization are two major nonexclusive causes of haplotype sharing between species. Distinguishing between these two processes is notoriously difficult as they can generate similar genetic signatures. Previous studies revealed that the mitochondrial DNA (mtDNA) differentiation between two East Asian long-tailed tits (Aegithalos bonvaloti and A. fuliginosus) was extremely low, even lower than intraspecific differentiation in some other long-tailed tits. Using a combination of multilocus and coalescent analyses, we explored the causes of the anomalous lack of mtDNA differentiation between the two species. RESULTS: The mtDNA divergence between the two species was shallow, while the nuclear DNA (nuDNA) divergence was considerably deeper. The IMa analyses based on the mtDNA dataset suggested relatively high gene flow from A. fuliginosus to A. bonvaloti, while negligible gene flow in the opposite direction. In contrast to mtDNA, the migration rates at autosomal and Z-linked nuDNA loci were negligible or much lower. The NEWHYBRIDS analysis assigned all individuals except one to pure parental species with high posterior probability. The Bayesian skyline plot showed that both species underwent population expansions during the Last Glacial Maximum (LGM), and the ecological niche modelling suggested that their ranges overlapped more during the LGM than at present. CONCLUSIONS: We suggest that historical hybridization, in combination with selective sweep and/or genetic drift might be the main causes of the extremely low mtDNA differentiation between the two species. The hybridization probably occurred mainly between A. fuliginosus females and A. bonvaloti males. The LGM distribution expansion might have facilitated hybridization, while the post-LGM distribution contraction could have facilitated some mtDNA sorting. Ongoing hybridization between the two species might be very limited, but further studies with more samples from the contact zone are needed to test this conclusion.

Convergent evolution of kidney sizes and supraorbital salt glands for birds living in saline habitats.

Only a small number of avian species inhabit salty environments. To understand how they adapted, we examined the evolution of kidney sizes, supraorbital salt glands (SSGs), and the utilization of salty habitats across 230 species spanning 25 avian orders. Phylogenetic analysis indicates that SSGs, large kidneys, and thriving in salty habitats emerged convergently in birds. Transition rate analysis reveals that species possessing SSGs and large kidneys tended to move from low-to high-salinity environments, while others moved in the opposite direction. However, habitat salinity also influenced kidney evolution; lineages residing in high-salinity environments tended to develop larger kidneys than those in low-salinity environments. Our findings suggest that SSGs and large kidneys may have evolved through adaptation to high salinity. Overall, habitat conditions and physiological traits influenced avian adaptation to salty environments in a reciprocal manner. These results shed the new light on the evolutionary mechanisms underlying functional diversity in birds.