Exploring the Effect of Asymmetric Mitochondrial DNA Introgression on Estimating Niche Divergence in Morphologically Cryptic Species

If potential morphologically cryptic species, identified based on differentiated mitochondrial DNA, express ecological divergence, this increases support for their treatment as distinct species. However, mitochondrial DNA introgression hampers the correct estimation of ecological divergence. We test the hypothesis that estimated niche divergence differs when considering nuclear DNA composition or mitochondrial DNA type as representing the true species range. We use empirical data of two crested newt species (Amphibia: Triturus) which possess introgressed mitochondrial DNA from a third species in part of their ranges. We analyze the data in environmental space by determining Fisher distances in a principal component analysis and in geographical space by determining geographical overlap of species distribution models. We find that under mtDNA guidance in one of the two study cases niche divergence is overestimated, whereas in the other it is underestimated. In the light of our results we discuss the role of estimated niche divergence in species delineation.     

Investigating niche and lineage diversification in widely distributed taxa: phylogeography and ecological niche modeling of the Peromyscus maniculatus species group

The relationship between lineage formation and variation in the ecological niche is a fundamental evolutionary question. Two prevailing hypotheses reflect this relationship: niche conservatism and niche divergence. Niche conservatism predicts a pattern where sister taxa will occupy similar niche spaces; whereas niche divergence predicts that sister taxa will occupy different niche spaces. Widely distributed species often show distinct phylogeographic structure, but little research has been conducted on how the environment may be related to these phylogenetic patterns. We investigated the relationship between lineage divergence and environmental space for the closely related species Peromyscus maniculatus and P. polionotus utilizing phylogenetic techniques and ecological niche modeling (ENM). We estimated the phylogenetic relationship among individuals based on complete cytochrome b sequences that represent individuals from a majority of the species ranges. Niche spaces that lineages occupy were estimated by using 12 environmental layers. Differences in niche space were tested using multivariate statistics based on location data, and ENMs were employed using maximum entropy algorithms. Two similarity indices estimated significant divergence in environmental space based on the ENM. Six geographically structured lineages were identified within P. maniculatus. Nested within P. maniculatus we found that P. polionotus recently diverged from a clade occupying central and western United States. We estimated that the majority of the genetic lineages occupy distinct environmental niches, which supports a pattern of niche divergence. Two sister taxa showed niche divergence and represent different ecomorphs, suggesting morphological, genetic and ecological divergence between the two lineages. Two other sister taxa were observed in the same environmental space based on multivariate statistics, suggesting niche conservatism. Overall our results indicate that a widely distributed species may exhibit both niche conservatism and niche divergence, and that most lineages seem to occupy distinct environmental niches.     

Niche position, but not niche breadth, differs in two coexisting amphibians having contrasting trends in Europe

Aim We explored general habitat‐related explanations (niche breadth and niche position) to the contrasting status of two amphibian species that have largely overlapping ranges and habitats – the rare and declining crested newt (Triturus cristatus), and the more common smooth newt (Triturus vulgaris). These closely related and ecologically similar species provide an excellent opportunity to study those methodologically challenging hypotheses, and this is the first such study on amphibians. Location Denmark. Methods We derived multivariate habitat models from 27 characteristics of 210 ponds and their surroundings, and their occupation by newts. In addition to the model performance, niche breadths were compared using the mean beta diversity of amphibian communities in the presence of each newt species. Results For each newt species, the habitat models comprised three variables and correctly classified 74% of observations. Diverse invertebrate fauna (prey base) and shorter distances to other ponds inhabited by conspecifics were positive for both species, while the surrounding habitat (notably dry grasslands with forests) was important for the crested newt and the sediment type of the pond for the smooth newt only. Beta diversity of the amphibian communities of occupied ponds did not differ between the two newt species. Hence, in an area of frequent coexistence, habitat requirements of the species differed in key variables, not in the extent of specialization. Main conclusions Our study supported the niche position rather than the niche‐breadth hypothesis of rarity. We suggest that the rarity and/or continuing decline of the crested newt is related to the degradation of (semi)natural terrestrial habitats around suitable water bodies in Europe. Consequently, special restoration of such habitats has a high potential for the recovery of this rare species, while general pond management appeared more beneficial for the common smooth newt.     

Morphological integration and ontogenetic niche shift: a study of crested newt limbs

This study deals with the ontogenetic and evolutionary aspects of integration patterns in the limbs of crested newt species, which, like most amphibians, have a biphasic life history with two morphologically distinct stages (larval vs. juvenile and adult) that occupy different environments (aquatic vs. terrestrial). We analyzed the structure and pattern of correlation between limb skeletal elements at three ontogenetic stages (larval, juvenile, and adult) of four closely related species that differ in their preferences of aquatic habitats (more terrestrial and more aquatic). We found dynamic changes in the pattern of morphological integration between successive ontogenetic stages, as well as changes over the course of crested newt phylogeny. Generally, equivalent ontogenetic stages of different species of crested newts show higher concordance in the correlation pattern than successive ontogenetic stages within species. Among species, two opposing correlation patterns were observed: in more terrestrial species, homologous limb elements are less correlated and within-limb elements are more correlated; in aquatic species, the reverse pattern occurs. These results indicate that the function seems to be the covariance-generating factor, which has shaped the patterns of morphological integration of crested newt limbs.     

Micro-habitat use within a guild of newt larvae (Trituridae) in an Alpine lake

Population and ecological parameters such as numbers of larvae, microhabitat use, niche breadth and niche overlap of three species of syntopic larval newts (Alpine newt Triturus alpestris, Italian crested newt T. carnifex, and common newt T. vulgaris) were studied for two years in a small pond at 1160 m a.s.l. in NE Slovenia. Differences in microhabitat partitioning among larval newts were small. The largest niche breadth was estimated for larval T. alpestris, and the narrowest estimate was for larval T. carnifex in both years. Ecological differences seem to be very small and quite variable among sites and years. It appears that the developmental stage and size of newt larvae are more important in explaining resource partitioning than the characteristics of each species. Because of the absence of potential invertebrate predators and adult newts in the second half of the breeding season, the injuries could only be caused by intra-and interspecific predation attempts.     

Adding more ecology into species delimitation: ecological niche models and phylogeography help define cryptic species in the black salamander (Aneides flavipunctatus)

Being able to efficiently and accurately delimit species is one of the most basic and important aspects of systematics because species are the fundamental unit of analysis in biogeography, ecology, and conservation. We present a rationale and approach for combining ecological niche modeling, spatially explicit analyses of environmental data, and phylogenetics in species delimitation, and we use our methodology in an empirical example focusing on Aneides flavipunctatus, the black salamander (Caudata: Plethodontidae), in California. We assess the relationships between genetic, environmental, and geographic distance among populations. We use 11 climatic variables and point locality data from public databases to create ecological niche models. The suitability of potential contact zones between parapatric lineages is also assessed using the data from ecological niche modeling. Phylogenetic analyses of portions of the mitochondrial genome reveal morphologically cryptic mitochondrial lineages in this species. In addition, we find that patterns of genetic divergence are strongly associated with divergence in the ecological niche. Our work demonstrates the ease and utility of using spatial analyses of environmental data and phylogenetics in species delimitation, especially for groups displaying fine-scaled endemism and cryptic species.     

Developmental life history is associated with variation in rates of climatic niche evolution in a salamander adaptive radiation*

Rates of climatic niche evolution vary widely across the tree of life and are strongly associated with rates of diversification among clades. However, why the climatic niche evolves more rapidly in some clades than others remains unclear. Variation in life history traits often plays a key role in determining the environmental conditions under which species can survive, and therefore, could impact the rate at which lineages can expand in available climatic niche space. Here, we explore the relationships among life-history variation, climatic niche breadth, and rates of climatic niche evolution. We reconstruct a phylogeny for the genus Desmognathus, an adaptive radiation of salamanders distributed across eastern North America, based on nuclear and mitochondrial genes. Using this phylogeny, we estimate rates of climatic niche evolution for species with long, short, and no aquatic larval stage. Rates of climatic niche evolution are unrelated to the mean climatic niche breadth of species with different life histories. Instead, we find that the evolution of a short larval period promotes greater exploration of climatic space, leading to increased rates of climatic niche evolution across species having this trait. We propose that morphological and physiological differences associated with variation in larval stage length underlie the heterogeneous ability of lineages to explore climatic niche space. Rapid rates of climatic niche evolution among species with short larval periods were an important dimension of the clade's adaptive radiation and likely contributed to the rapid rate of lineage accumulation following the evolution of an aquatic life history in this clade. Our results show how variation in a key life-history trait can constrain or promote divergence of the climatic niche, leading to variation in rates of climatic niche evolution among species.     

Cryptic niche conservatism among evolutionary lineages of an invasive lizard

Aim There is increasing evidence that the quality and breadth of ecological niches vary among individuals, populations, evolutionary lineages and therefore also across the range of a species. Sufficient knowledge about niche divergence among clades might thus be crucial for predicting the invasion potential of species. We tested for the first time whether evolutionary lineages of an invasive species vary in their climate niches and invasive potential. Furthermore, we tested whether lineage‐specific models show a better performance than combined models. Location Europe. Methods We used species distribution models (SDMs) based on climatic information at native and invasive ranges to test for intra‐specific niche divergence among mitochondrial DNA (mtDNA) clades of the invasive wall lizard Podarcis muralis. Using DNA barcoding, we assigned 77 invasive populations in Central Europe to eight geographically distinct evolutionary lineages. Niche similarity among lineages was assessed and the predictive power of a combination of clade‐specific SDMs was compared with a combined SDM using the pooled records of all lineages. Results We recorded eight different invasive mtDNA clades in Central Europe. The analysed clades had rather similar realized niches in their native and invasive ranges, whereas inter‐clade niche differentiation was comparatively strong. However, we found only a weak correlation between geographic origin (i.e. mtDNA clade) and invasive occurrences. Clades with narrow realized niches still became successful invaders far outside their native range, most probably due to broader fundamental niches. The combined model using data for all invasive lineages achieved a much better prediction of the invasive potential. Conclusions Our results indicate that the observed niche differentiation among evolutionary lineages is mainly driven by niche realization and not by differences in the fundamental niches. Such cryptic niche conservatism might hamper the success of clade‐specific niche modelling. Cryptic niche conservatism may in general explain the invasion success of species in areas with apparently unsuitable climate.     

Niche differentiation in a postglacial colonizer,the bank vole Clethrionomys glareolus

Species‐level environmental niche modeling has been crucial in efforts to understand how species respond to climate variation and change. However, species often exhibit local adaptation and intraspecific niche differences that may be important to consider in predicting responses to climate. Here, we explore whether phylogeographic lineages of the bank vole originating from different glacial refugia (Carpathian, Western, Eastern, and Southern) show niche differentiation, which would suggest a role for local adaptation in biogeography of this widespread Eurasian small mammal. We first model the environmental requirements for the bank vole using species‐wide occurrences (210 filtered records) and then model each lineage separately to examine niche overlap and test for niche differentiation in geographic and environmental space. We then use the models to estimate past [Last Glacial Maximum (LGM) and mid‐Holocene] habitat suitability to compare with previously hypothesized glacial refugia for this species. Environmental niches are statistically significantly different from each other for all pairs of lineages in geographic and environmental space, and these differences cannot be explained by habitat availability within their respective ranges. Together with the inability of most of the lineages to correctly predict the distributions of other lineages, these results support intraspecific ecological differentiation in the bank vole. Model projections of habitat suitability during the LGM support glacial survival of the bank vole in the Mediterranean region and in central and western Europe. Niche differences between lineages and the resulting spatial segregation of habitat suitability suggest ecological differentiation has played a role in determining the present phylogeographic patterns in the bank vole. Our study illustrates that models pooling lineages within a species may obscure the potential for different responses to climate change among populations.     

One Species, Three Pleistocene Evolutionary Histories: Phylogeography of the Italian Crested Newt, Triturus carnifex

Phylogeographic patterns of temperate species from the Mediterranean peninsulas have been investigated intensively. Nevertheless, as more phylogeographies become available, either unique patterns or new lines of concordance continue to emerge, providing new insights on the evolution of regional biotas. Here, we investigated the phylogeography and evolutionary history of the Italian crested newt, Triturus carnifex, through phylogenetic, molecular dating and population structure analyses of two mitochondrial gene fragments (ND2 and ND4; overall 1273 bp). We found three main mtDNA lineages having parapatric distribution and estimated divergence times between Late Pliocene and Early Pleistocene. One lineage (S) was widespread south of the northern Apennine chain and was further geographically structured into five sublineages, likely of Middle Pleistocene origin. The second lineage (C) was widespread throughout the Padano-Venetian plain and did not show a clear phylogeographic structure. The third lineage (N) was observed in only two populations located on western Croatia/Slovenia. Results of analysis of molecular variance suggested that partitioning populations according to the geographic distribution of these lineages and sublineages explains 76% of the observed genetic variation. The phylogeographic structure observed within T. carnifex and divergence time estimates among its lineages, suggest that responses to Pleistocene environmental changes in this single species have been as diverse as those found previously among several codistributed temperate species combined. Consistent with the landscape heterogeneity, physiographic features, and palaeogeographical evolution of its distribution range, these responses encompass multiple refugia along the Apennine chain, lowland refugia in large peri-coastal plains, and a 'cryptic' northern refugium.     

Postglacial species displacement in Triturus newts deduced from asymmetrically introgressed mitochondrial DNA and ecological niche models

ABSTRACT: BACKGROUND: If the geographical displacement of one species by another is accompanied by hybridization, mitochondrial DNA can introgress asymmetrically, from the outcompeted species into the invading species, over a large area. We explore this phenomenon using the two parapatric crested newt species, Triturus macedonicus and T. karelinii, distributed on the Balkan Peninsula in south-eastern Europe, as a model. RESULTS: We first delimit a ca. 54,000 km2 area in which T. macedonicus contains T. karelinii mitochondrial DNA. This introgression zone bisects the range of T. karelinii, cutting off a T. karelinii enclave. The high similarity of introgressed mitochondrial DNA haplotypes with those found in T. karelinii suggests a recent transfer across the species boundary. We then use ecological niche modeling to explore habitat suitability of the location of the present day introgression zone under current, mid-Holocene and Last Glacial Maximum conditions. This area was inhospitable during the Last Glacial Maximum for both species, but would have been habitable at the mid-Holocene. Since the mid-Holocene, habitat suitability generally increased for T. macedonicus, whereas it decreased for T. karelinii. CONCLUSION: The presence of a T. karelinii enclave suggests that T. karelinii was the first to colonize the area where the present day introgression zone is positioned after the Last Glacial Maximum. Subsequently, we propose T. karelinii was outcompeted by T. macedonicus, which captured T. karelinii mitochondrial DNA via introgressive hybridization in the process. Ecological niche modeling suggests that this replacement was likely facilitated by a shift in climate since the mid-Holocene. We suggest that the northwestern part of the current introgression zone was probably never inhabited by T. karelinii itself, and that T. karelinii mitochondrial DNA spread there through T. macedonicus exclusively. Considering the spatial distribution of the introgressed mitochondrial DNA and the signal derived from ecological niche modeling, we do not favor the hypothesis that foreign mitochondrial DNA was pulled into the T. macedonicus range by natural selection.     

Trophic specialization influences the rate of environmental niche evolution in damselfishes (Pomacentridae)

The rate of environmental niche evolution describes the capability of species to explore the available environmental space and is known to vary among species owing to lineage-specific factors. Trophic specialization is a main force driving species evolution and is responsible for classical examples of adaptive radiations in fishes. We investigate the effect of trophic specialization on the rate of environmental niche evolution in the damselfish, Pomacentridae, which is an important family of tropical reef fishes. First, phylogenetic niche conservatism is not detected in the family using a standard test of phylogenetic signal, and we demonstrate that the environmental niches of damselfishes that differ in trophic specialization are not equivalent while they still overlap at their mean values. Second, we estimate the relative rates of niche evolution on the phylogenetic tree and show the heterogeneity among rates of environmental niche evolution of the three trophic groups. We suggest that behavioural characteristics related to trophic specialization can constrain the evolution of the environmental niche and lead to conserved niches in specialist lineages. Our results show the extent of influence of several traits on the evolution of the environmental niche and shed new light on the evolution of damselfishes, which is a key lineage in current efforts to conserve biodiversity in coral reefs.     

Stream diatoms exhibit weak niche conservation along global environmental and climatic gradients

Niche conservatism (NC) describes the scenario in which species retain similar characteristics or traits over time and space, and thus has potentially important implications for understanding their biogeographic distributions. Evidence consistent with NC includes similar niche properties across geographically distant regions. We investigated whether NC was evident in stream diatom morphospecies by modeling species responses to environmental and climatic variables in a set of calibration sites (from the US) and then evaluated the models with test sets (from France, Finland, New Zealand, Antilles and La Réunion). We also examined whether diatom species showed congruency in environmental niche optima and niche breadths between the study regions, and whether species occupancy and functional traits influenced the observed patterns. We used boosted regression tree models with local environmental variables and climatic variables as predictors. We detected low NC in both environmental and climate models and a lack of consistent differences in NC between widely distributed and regionally rare species and among functional groups. For all species, diatom environmental and climatic optima varied clearly between the regions but showed some positive relationships especially for pH and total phosphorus. Diatom niche breadths were only weakly correlated between the US and the other regions. We demonstrated that diatoms showed overall relatively little NC globally, and NC was especially low for climatic variables. Collectively, these findings suggest that there may exist locally adapted lineages within the diatom morphospecies or diatoms possess some adaptation potential for differences in temperature. We argue that in diatoms, environmental and especially climate models may not be transferrable in space globally but need regional diatom data for calibration because species niches seem to differ among geographical regions.     

DOES NICHE DIVERGENCE ACCOMPANY ALLOPATRIC DIVERGENCE IN APHELOCOMA JAYS AS PREDICTED UNDER ECOLOGICAL SPECIATION?: INSIGHTS FROM TESTS WITH NICHE MODELS

The role of ecology in the origin of species has been the subject of long‐standing interest to evolutionary biologists. New sources of spatially explicit ecological data allow for large‐scale tests of whether speciation is associated with niche divergence or whether closely related species tend to be similar ecologically (niche conservatism). Because of the confounding effects of spatial autocorrelation of environmental variables, we generate null expectations for niche divergence for both an ecological‐niche modeling and a multivariate approach to address the question: do allopatrically distributed taxa occupy similar niches? In a classic system for the study of niche evolution—the Aphelocoma jays—we show that there is little evidence for niche divergence among Mexican Jay (A. ultramarina) lineages in the process of speciation, contrary to previous results. In contrast, Aphelocoma species that exist in partial sympatry in some regions show evidence for niche divergence. Our approach is widely applicable to the many cases of allopatric lineages in the beginning stages of speciation. These results do not support an ecological speciation model for Mexican Jay lineages because, in most cases, the allopatric environments they occupy are not significantly more divergent than expected under a null model.     

Measuring ecological niche overlap from occurrence and spatial environmental data

Aim Concerns over how global change will influence species distributions, in conjunction with increased emphasis on understanding niche dynamics in evolutionary and community contexts, highlight the growing need for robust methods to quantify niche differences between or within taxa. We propose a statistical framework to describe and compare environmental niches from occurrence and spatial environmental data. Location Europe, North America and South America. Methods The framework applies kernel smoothers to densities of species occurrence in gridded environmental space to calculate metrics of niche overlap and test hypotheses regarding niche conservatism. We use this framework and simulated species with pre‐defined distributions and amounts of niche overlap to evaluate several ordination and species distribution modelling techniques for quantifying niche overlap. We illustrate the approach with data on two well‐studied invasive species. Results We show that niche overlap can be accurately detected with the framework when variables driving the distributions are known. The method is robust to known and previously undocumented biases related to the dependence of species occurrences on the frequency of environmental conditions that occur across geographical space. The use of a kernel smoother makes the process of moving from geographical space to multivariate environmental space independent of both sampling effort and arbitrary choice of resolution in environmental space. However, the use of ordination and species distribution model techniques for selecting, combining and weighting variables on which niche overlap is calculated provide contrasting results. Main conclusions The framework meets the increasing need for robust methods to quantify niche differences. It is appropriate for studying niche differences between species, subspecies or intra‐specific lineages that differ in their geographical distributions. Alternatively, it can be used to measure the degree to which the environmental niche of a species or intra‐specific lineage has changed over time.     

Cryptic diversity in European bats

Different species of bat can be morphologically very similar. In order to estimate the amount of cryptic diversity among European bats we screened the intra- and interspecific genetic variation in 26 European vespertilionid bat species. We sequenced the DNA of subunit 1 of the mitochondrial protein NADH dehydrogenase (ND1) from several individuals of a species, which were sampled in a variety of geographical regions. A phylogeny based on the mitochondrial (mt) DNA data is in good agreement with the current classification in the family. Highly divergent mitochondrial lineages were found in two taxa, which differed in at least 11% of their ND1 sequence. The two mtDNA lineages in Plecotus austriacus correlated with the two subspecies Plecotus austriacus austriacus and Plecotus austriacus kolombatovici. The two mtDNA lineages in Myotis mystacinus were partitioned among two morphotypes. The evidence for two new bat species within Europe is discussed. Convergent adaptive evolution might have contributed to the morphological similarity among distantly related species if they occupy similar ecological niches. Closely related species may differ in their ecology but not necessarily in their morphology. On the other hand, two morphologically clearly different species (Eptesicus serotinus and Eptesicus nilssonii) were found to be genetically very similar. Neither morphological nor mitochondrial DNA sequence analysis alone can be guaranteed to identify species.     

Integration of conflict into integrative taxonomy: fitting hybridization in species delimitation of Mesocarabus (Coleoptera: Carabidae)

In species differentiation, characters may not diverge synchronously, and there are also processes that shuffle character states in lineages descendant from a common ancestor. Species are thus expected to show some degree of incongruence among characters; therefore, taxonomic delimitation can benefit from integrative approaches and objective strategies that account for character conflict. We illustrate the potential of exploiting conflict for species delimitation in a study case of ground beetles of the subgenus Carabus (Mesocarabus), where traditional taxonomy does not accurately delimit species. The molecular phylogenies of four mitochondrial and three nuclear genes, cladistic analysis of the aedeagus, ecological niche divergence and morphometry of pronotal shape in more than 500 specimens of Mesocarabus show that these characters are not fully congruent. For these data, a three‐step operational strategy is proposed for species delimitation by (i) delineating candidate species based on the integration of incongruence among conclusive lines of evidence, (ii) corroborating candidate species with inconclusive lines of evidence and (iii) refining a final species proposal based on an integrated characterization of candidate species based on the evolutionary analysis of incongruence. This procedure provided a general understanding of the reticulate process of hybridization and introgression acting on Mesocarabus and generated the hypothesis of seven Mesocarabus species, including two putative hybrid lineages. Our work emphasizes the importance of incorporating critical analyses of character and phylogenetic conflict to infer both the evolutionary history and species boundaries through an integrative taxonomic approach.     

Molecular phylogeny of the critically endangered Indochinese box turtle (Cuora galbinifrons)

The Indochinese box turtle Cuora galbinifrons is a polytypic, critically endangered species from Vietnam, Laos, and Hainan Island, China. We analyze up to 1790bp of mitochondrial DNA under maximum parsimony and maximum likelihood criteria to test if the five historically recognized subspecies represent evolutionary lineages, and to elucidate the relationship of C. galbinifrons to other Cuora. C. galbinifrons is composed of three major mitochondrial DNA clades corresponding to the three subspecies galbinifrons, bourreti, and picturata. These three lineages are also morphologically diagnosable, and consequently we recommend elevating each to full species. Cuora galbinifrons hainanensis nests within the galbinifrons clade, and we retain it as a synonym of galbinifrons, as supported by morphology. Cuora "serrata" is known to be a hybrid of male Cuora mouhotii and female C. galbinifrons, and our findings show that C. "serrata" originates from both female galbinifrons and bourreti. Little or no mitochondrial DNA variation was found among the morphologically distinct species Cuora aurocapitata, Cuora pani, and Cuora trifasciata, for which hypotheses are proposed. Recognizing galbinifrons, bourreti, and picturata as separate species has consequences for ongoing ex situ captive breeding programs and prioritization of in situ conservation activities, particularly in Vietnam.     

Phylogeny and diversification of the gallopheasants (Aves: Galliformes): Testing roles of sexual selection and environmental niche divergence

The gallopheasants comprise a clade of 22 species including some of the most elaborately plumaged and highly ornamented birds in the world. They also occupy a remarkable breath of environments and habitats, ranging from lowland rainforests to high grasslands and steppes of the Tibetan plateau. Here, we provide the first well-resolved species phylogeny of this charismatic group, inferred from ultraconserved elements, nuclear introns and mitochondrial DNA sequences. Unlike previous studies which found unresolvable relationships and suggested a rapid initial burst of diversification, we identified a well-resolved phylogeny supported in both concatenated and coalescent analytical frameworks, and a steady accrual of lineages through time. Morphological trait reconstructions demonstrated strong phylogenetic signal, not only for highly ornamented males, but also in more cryptically plumaged females. Environmental niche similarly exhibited strong phylogenetic signal. Moreover, evolution of male traits, female traits and environmental niche were all significantly correlated, making it difficult to disentangle their individual roles in gallopheasant diversification.