Genome-scale phylogenetics: inferring the plant tree of life from 18,896 gene trees

Phylogenetic analyses using genome-scale data sets must confront incongruence among gene trees, which in plants is exacerbated by frequent gene duplications and losses. Gene tree parsimony (GTP) is a phylogenetic optimization criterion in which a species tree that minimizes the number of gene duplications induced among a set of gene trees is selected. The run time performance of previous implementations has limited its use on large-scale data sets. We used new software that incorporates recent algorithmic advances to examine the performance of GTP on a plant data set consisting of 18,896 gene trees containing 510,922 protein sequences from 136 plant taxa (giving a combined alignment length of >2.9 million characters). The relationships inferred from the GTP analysis were largely consistent with previous large-scale studies of backbone plant phylogeny and resolved some controversial nodes. The placement of taxa that were present in few gene trees generally varied the most among GTP bootstrap replicates. Excluding these taxa either before or after the GTP analysis revealed high levels of phylogenetic support across plants. The analyses supported magnoliids sister to a eudicot + monocot clade and did not support the eurosid I and II clades. This study presents a nuclear genomic perspective on the broad-scale phylogenic relationships among plants, and it demonstrates that nuclear genes with a history of duplication and loss can be phylogenetically informative for resolving the plant tree of life.     

Reconciling gene and genome duplication events: using multiple nuclear gene families to infer the phylogeny of the aquatic plant family Pontederiaceae

Most plant phylogenetic inference has used DNA sequence data from the plastid genome. This genome represents a single genealogical sample with no recombination among genes, potentially limiting the resolution of evolutionary relationships in some contexts. In contrast, nuclear DNA is inherently more difficult to employ for phylogeny reconstruction because major mutational events in the genome, including polyploidization, gene duplication, and gene extinction can result in homologous gene copies that are difficult to identify as orthologs or paralogs. Gene tree parsimony (GTP) can be used to infer the rooted species tree by fitting gene genealogies to species trees while simultaneously minimizing the estimated number of duplications needed to reconcile conflicts among them. Here, we use GTP for five nuclear gene families and a previously published plastid data set to reconstruct the phylogenetic backbone of the aquatic plant family Pontederiaceae. Plastid-based phylogenetic studies strongly supported extensive paraphyly of Eichhornia (one of the four major genera) but also depicted considerable ambiguity concerning the true root placement for the family. Our results indicate that species trees inferred from the nuclear genes (alone and in combination with the plastid data) are highly congruent with gene trees inferred from plastid data alone. Consideration of optimal and suboptimal gene tree reconciliations place the root of the family at (or near) a branch leading to the rare and locally restricted E. meyeri. We also explore methods to incorporate uncertainty in individual gene trees during reconciliation by considering their individual bootstrap profiles and relate inferred excesses of gene duplication events on individual branches to whole-genome duplication events inferred for the same branches. Our study improves understanding of the phylogenetic history of Pontederiaceae and also demonstrates the utility of GTP for phylogenetic analysis.     

Phylogeny of γ‐proteobacteria: resolution of one branch of the universal tree?

The reconstruction of bacterial evolutionary relationships has proven to be a daunting task because variable mutation rates and horizontal gene transfer (HGT) among species can cause grave incongruities between phylogenetic trees based on single genes. Recently, a highly robust phylogenetic tree was constructed for 13 gamma-proteobacteria using the combined alignments of 205 conserved orthologous proteins.1 Only two proteins had incongruent tree topologies, which were attributed to HGT between Pseudomonas species and Vibrio cholerae or enterics. While the evolutionary relationships among these species appears to be resolved, further analysis suggests that HGT events with other bacterial partners likely occurred; this alters the implicit assumption of gamma-proteobacteria monophyly. Thus, any thorough reconstruction of bacterial evolution must not only choose a suitable set of molecular markers but also strive to reduce potential bias in the selection of species.     

A multilocus phylogeny of the Sulidae (Aves: Pelecaniformes)

Gene trees will often differ from the true species history, the species tree, as a result of processes such as incomplete lineage sorting. New methods such as Bayesian Estimation of the Species Tree (BEST) use the multispecies coalescent to model lineage sorting, and directly infer the species tree from multilocus DNA sequence data. The Sulidae (Aves: Pelecaniformes) is a family of ten booby and gannet species with a global distribution. We sequenced five nuclear intron loci and one mitochondrial locus to estimate a species tree for the Sulidae using both BEST and by concatenating nuclear loci. We also used fossil calibrated strict and relaxed molecular clocks in BEAST to estimate divergence times for major nodes in the sulid phylogeny. Individual gene trees showed little phylogenetic conflict but varied in resolution. With the exception of the mitochondrial gene tree, no gene tree was completely resolved. On the other hand, both the BEST and concatenated species trees were highly resolved, strongly supported, and topologically consistent with each other. The three sulid genera (Morus, Sula, Papasula) were monophyletic and the relationships within genera were mostly consistent with both a previously estimated mtDNA gene tree and the mtDNA gene tree estimated here. However, our species trees conflicted with the mtDNA gene trees in the relationships among the three genera. Most notably, we find that the endemic and endangered Abbott's booby (Papasula abbotti) is likely basal to all other members of the Sulidae and diverged from them approximately 22 million years ago.     

Gene trees, species trees, and morphology converge on a similar phylogeny of living gars (Actinopterygii: Holostei: Lepisosteidae), an ancient clade of ray-finned fishes

Extant gars represent the remaining members of a formerly diverse assemblage of ancient ray-finned fishes and have been the subject of multiple phylogenetic analyses using morphological data. Here, we present the first hypothesis of phylogenetic relationships among living gar species based on molecular data, through the examination of gene tree heterogeneity and coalescent species tree analyses of a portion of one mitochondrial (COI) and seven nuclear (ENC1, myh6, plagl2, S7 ribosomal protein intron 1, sreb2, tbr1, and zic1) genes. Individual gene trees displayed varying degrees of resolution with regards to species-level relationships, and the gene trees inferred from COI and the S7 intron were the only two that were completely resolved. Coalescent species tree analyses of nuclear genes resulted in a well-resolved and strongly supported phylogenetic tree of living gar species, for which Bayesian posterior node support was further improved by the inclusion of the mitochondrial gene. Species-level relationships among gars inferred from our molecular data set were highly congruent with previously published morphological phylogenies, with the exception of the placement of two species, Lepisosteus osseus and L. platostomus. Re-examination of the character coding used by previous authors provided partial resolution of this topological discordance, resulting in broad concordance in the phylogenies inferred from individual genes, the coalescent species tree analysis, and morphology. The completely resolved phylogeny inferred from the molecular data set with strong Bayesian posterior support at all nodes provided insights into the potential for introgressive hybridization and patterns of allopatric speciation in the evolutionary history of living gars, as well as a solid foundation for future examinations of functional diversification and evolutionary stasis in a "living fossil" lineage.     

基于COⅡ基因序列的斑腿蝗科部分亚科的分子系统学研究

采用PCR产物直接测序法测定了斑腿蝗科10个亚科16属22种的COⅡ基因585 bp的片段, 对序列的碱基组成进行了分析,并评估了数据集的系统发育信号;最后,以癞蝗科的肃南 短鼻蝗作为外群,采用NJ法、MP法、ML法以及贝叶斯推论法构建了系统树,以解决这些物种所代表的亚科之间的系统发育关系。结果表明：22种斑腿蝗科昆虫的COⅡ基因序列碱基组成表现强烈的A+T含量偏向性。对COⅡ基因585 bp序列片段构成的全数据组和根据密码子不同位点划分的密码子第一、第二和第三位点数据组的系统发育信号分析显示,所有数据组都具有一定的系统发育信息。在4种方法得到的合一树中发现： （1）星翅蝗亚科、刺胸蝗亚科、黑背蝗亚科、斑腿蝗亚科的亲缘关系较近;（2）卵翅蝗亚科与稻蝗亚科亲缘关系较近,建议卵翅蝗亚科似乎应归入稻蝗亚科中,板胸蝗亚科与这两个亚科的关系较近;（3）黑蝗亚科和秃蝗亚科似乎应合并为一个亚科;（4）切翅蝗亚科的4个属未聚在一起,表明这些属的区别较大,不是一个单系群;（5）黑蝗亚科和秃蝗亚科关系较近,且与本研究中其他几个亚科的亲缘关系相对较远。研究结果表明COⅡ基因在解决斑腿蝗科的亚科以下属种间的系统发育关系时是一个有效的分子标记。     

Estimating phylogenetic relationships despite discordant gene trees across loci: the species tree of a diverse species group of feather mites (Acari: Proctophyllodidae)

With the increased availability of multilocus sequence data, the lack of concordance of gene trees estimated for independent loci has focused attention on both the biological processes producing the discord and the methodologies used to estimate phylogenetic relationships. What has emerged is a suite of new analytical tools for phylogenetic inference--species tree approaches. In contrast to traditional phylogenetic methods that are stymied by the idiosyncrasies of gene trees, approaches for estimating species trees explicitly take into account the cause of discord among loci and, in the process, provides a direct estimate of phylogenetic history (i.e. the history of species divergence, not divergence of specific loci). We illustrate the utility of species tree estimates with an analysis of a diverse group of feather mites, the pinnatus species group (genus Proctophyllodes). Discord among four sequenced nuclear loci is consistent with theoretical expectations, given the short time separating speciation events (as evident by short internodes relative to terminal branch lengths in the trees). Nevertheless, many of the relationships are well resolved in a Bayesian estimate of the species tree; the analysis also highlights ambiguous aspects of the phylogeny that require additional loci. The broad utility of species tree approaches is discussed, and specifically, their application to groups with high speciation rates--a history of diversification with particular prevalence in host/parasite systems where species interactions can drive rapid diversification.     

Protein coding gene trees in ciliates: comparison with rRNA-based phylogenies

We have reexamined the phylogeny of the ciliates using alpha-tubulin and phosphoglycerate kinase gene sequences. For alpha-tubulin, we have compared the amino acid and nucleotide sequences of 20 species representing seven of the nine classes of the phylum (Karyorelictea, Heterotrichea, Hypotrichea, Oligohymenophorea, Colpodea, Nassophorea, and Litostomatea). The phylogenetic tree resembles a bush from which three monophyletic lineages can be distinguished which correspond to the three classes Hypotrichea, Oligohymenophorea, and Litostomatea. For phosphoglycerate kinase, we have compared the amino acid sequences from 7 species representing three classes (Heterotrichea, Hypotrichea, and Oligohymenophorea). The branching pattern is resolved in three deeply separated branches with an early emergence of the heterotrich. Our comparative analysis shows that if alpha-tubulin phylogeny is not informative at the interclass level, the preliminary data from the phosphoglycerate kinase molecule appear more promising. Nevertheless, at low taxonomic level and at the class level, the resolved phylogenetic relationships inferred from both protein and rRNA sequence data are congruent.     

Nuclear gene sequences resolve species phylogeny and mitochondrial introgression in Leptocarabus beetles showing trans-species polymorphisms

We studied the phylogenetic relationships among Japanese Leptocarabus ground beetles, which show extensive trans-species polymorphisms in mitochondrial gene genealogies. Simultaneous analysis of combined nuclear data with partial sequences from the long-wavelength rhodopsin, wingless, phosphoenolpyruvate carboxykinase, and 28S rRNA genes resolved the relationships among the five species, although separate analyses of these genes provided topologies with low resolution. For both the nuclear gene tree resulting from the combined data from four genes and a mitochondrial cytochrome oxidase subunit I (COI) gene tree, we applied a Bayesian divergence time estimation using a common calibration method to identify mitochondrial introgression events that occurred after speciation. Three mitochondrial lineages shared by two or three species were likely subject to introgression due to interspecific hybridization because the coalescent times for these lineages were much shorter than the corresponding speciation times estimated from nuclear gene sequences. We demonstrated that when species phylogeny is fully resolved with nuclear gene sequence data, comparative analysis of nuclear and mitochondrial gene trees can be used to infer introgressive hybridization events that might cause trans-species polymorphisms in mitochondrial gene trees.     

Assessing phylogenetic resolution among mitochondrial, nuclear, and morphological datasets in Nothonotus darters (Teleostei: Percidae)

External morphological characters are the basis of our understanding of diversity and species relationships in many darter clades. The past decade has seen the publication of many studies utilizing mtDNA sequence data to investigate darter phylogenetics, but only recently have nuclear genes been used to investigate darter relationships. Despite a long tradition of use in darter systematics few studies have examined the phylogenetic utility of external morphological characters in estimating relationships among species in darter clades. We present DNA sequence data from the mitochondrial cytochrome b (cytb) gene, the nuclear encoded S7 intron 1, and discretely coded external morphological characters for all 20 species in the darter clade Nothonotus. Bayesian phylogenetic analyses result in phylogenies that are in broad agreement with previous studies. The cytb gene tree is well resolved, while the nuclear S7 gene tree lacks phylogenetic resolution, node support, and is characterized by a lack of reciprocal monophyly for many of the Nothonotus species. The phylogenies resulting from analysis of the morphological dataset lack resolution, but nodes present are found in the cytb and S7 gene trees. The highest resolution and node support is found in the Bayesian combined data phylogeny. Based on our results we propose continued exploration of the phylogenetic utility of external morphological characters in other darter clades. Given the extensive lack of reciprocal monophyly of species observed in the S7 gene tree we predict that nuclear gene sequences may have limited utility in intraspecific phylogeographic studies of Nothonotus darters.     

Molecular evidence of reticulate evolution in the subgenus Plantago (Plantaginaceae)

Polyploidization is a frequent evolutionary event in plants that has a large influence on speciation and evolution of the genome. Molecular phylogenetic analyses of the taxonomically complex subgenus Plantago were conducted to elucidate intrasubgeneric phylogenetic relationships. A nuclear-encoding single-copy gene, SUC1 (1.0-1.8 kb), was sequenced in 24 taxa representing all five sections of the subgenus Plantago and two taxa from subgenus Coronopus as the outgroup. Fifteen known polyploids and one putative polyploid were sampled to examine polyploid origins and occurrence of reticulate evolution by cloning and sequence analysis of SUC1. Phylogenetic relationships were estimated using maximum parsimony, neighbor-joining, and Bayesian analyses. For the first time, our analysis provides a highly resolved phylogenetic tree. Subgenus Plantago formed a well-supported monophyletic clade. In contrast, alleles from polyploid species were scattered across the whole SUC1 phylogenetic tree, and some independent allopolyploids originated from hybridization between distant lineages. One reason for this taxonomic complexity can be attributed to reticulate evolution within the subgenus Plantago. Our results also suggest the possibility of two independent long-distance dispersals between the northern and southern hemispheres.     

Reconciling gene trees of eastern North American minnows

Most eastern North American cyprinid fishes belong to a clade known as the "open posterior myodome" (OPM) minnows, but phylogenetic relationships within this clade have been difficult to ascertain. Previous attempts to resolve relationships among the generally benthic "chubs" and the more pelagic "shiners", that constitute the majority of OPM minnows, have led to highly discordant phylogenetic hypotheses. To further examine relationships among the OPM minnows, we utilized both a concatenated Bayesian approach and a coalescent-based species tree method to analyze data from six protein coding nuclear loci (Enc1, Ptr, Ryr3, Sh3px3, Tbr1, and Zic1), as well as the mitochondrial locus (Cytb). We focused our analyses on the chub-like genus Phenacobius, a group that has drifted topologically between other benthic chubs and the more pelagic shiners, and also included exemplar taxa from 11 other OPM lineages. Individual gene trees were highly discordant regarding relationships within Phenacobius and across the OPM clade. The concatenated Bayesian analysis and coalescent-based species tree reconstruction recovered slightly different phylogenetic topologies. Additionally, the posterior support values for clades using the coalescent-based approach were consistently lower than the concatenated analysis. However, Phenacobius was resolved as monophyletic and as the sister lineage to Erimystax regardless of the combined data approach taken. Furthermore, Phenacobius+Erimystax was recovered as more closely related to the shiners we examined than to other chubs. Relationships within Phenacobius varied depending on the combined phylogenetic method utilized. Our results highlight the importance of multi-locus, coalescent-based approaches for resolving the phylogeny of diverse clades like the eastern North American OPM minnows.     

Molecular systematics of the deep-sea hydrothermal vent endemic Brachyuran family Bythograeidae: a comparison of three Bayesian species tree methods

Brachyuran crabs of the family Bythograeidae are endemic to deep-sea hydrothermal vents and represent one of the most successful groups of macroinvertebrates that have colonized this extreme environment. Occurring worldwide, the family includes six genera (Allograea, Austinograea, Bythograea, Cyanagraea, Gandalfus, and Segonzacia) and fourteen formally described species. To investigate their evolutionary relationships, we conducted Maximum Likelihood and Bayesian molecular phylogenetic analyses, based on DNA sequences from fragments of three mitochondrial genes (16S rDNA, Cytochrome oxidase I, and Cytochrome b) and three nuclear genes (28S rDNA, the sodium-potassium ATPase a-subunit 'NaK', and Histone H3A). We employed traditional concatenated (i.e., supermatrix) phylogenetic methods, as well as three recently developed Bayesian multilocus methods aimed at inferring species trees from potentially discordant gene trees. We found strong support for two main clades within Bythograeidae: one comprising the members of the genus Bythograea; and the other comprising the remaining genera. Relationships within each of these two clades were partially resolved. We compare our results with an earlier hypothesis on the phylogenetic relationships among bythograeid genera based on morphology. We also discuss the biogeography of the family in the light of our results. Our species tree analyses reveal differences in how each of the three methods weighs conflicting phylogenetic signal from different gene partitions and how limits on the number of outgroup taxa may affect the results.     

Are phylogenies resolved at the genus level appropriate for studies on phylogenetic structure of species assemblages?

Phylogenies are essential to studies investigating the effect of evolutionary history on assembly of species in ecological communities and geographical and ecological patterns of phylogenetic structure of species assemblages. Because phylogenies well resolved at the species level are lacking for many major groups of organisms such as vascular plants, researchers often generate a species-level phylogenies using a phylogeny well resolved at the genus level as a backbone and attaching species to their respective genera in the phylogeny as polytomies or by using a megaphylogeny well resolved at the genus level as a backbone and adding additional species to the megaphylogeny as polytomies of their respective genera. However, whether the result of a study using species-level phylogenies generated in these ways is robust, compared to that based on phylogenies fully resolved at the species level, has not been assessed. Here, we use 1093 angiosperm tree assemblages (each in a 110 × 110 km quadrat) in North America as a model system to address this question, by examining six commonly used metrics of phylogenetic structure (phylogenetic diversity and phylogenetic relatedness) and six climate variables commonly used in ecology. Our results showed that (1) the scores of phylogenetic metrics derived from species-level phylogenies resolved at the genus level with species being attached to their respective genera as polytomies are very strongly or perfectly correlated to those derived from a phylogeny fully resolved at the species level (the mean of correlation coefficients is 0.973), and (2) the relationships between the scores of phylogenetic metrics and climate variables are consistent between the two sets of analyses based on the two types of phylogeny. Our study suggests that using species-level phylogenies resolved at the genus level with species being attached to their genera as polytomies is appropriate in studies exploring patterns of phylogenetic structure of species in ecological communities across geographical and ecological gradients.     

Phylogenetic analysis of RAD‐seq data: examining the influence of gene genealogy conflict on analysis of concatenated data

One of the major issues in phylogenetic analysis is that gene genealogies from different gene regions may not reflect the true species tree or history of speciation. This has led to considerable debate about whether concatenation of loci is the best approach for phylogenetic analysis. The application of Next‐generation sequencing techniques such as RAD‐seq generates thousands of relatively short sequence reads from across the genomes of the sampled taxa. These data sets are typically concatenated for phylogenetic analysis leading to data sets that contain millions of base pairs per taxon. The influence of gene region conflict among so many loci in determining the phylogenetic relationships among taxa is unclear. We simulated RAD‐seq data by sampling 100 and 500 base pairs from alignments of over 6000 coding regions that each produce one of three highly supported alternative phylogenies of seven species of Drosophila. We conducted phylogenetic analyses on different sets of these regions to vary the sampling of loci with alternative gene trees to examine the effect on detecting the species tree. Irrespective of sequence length sampled per region and which subset of regions was used, phylogenetic analyses of the concatenated data always recovered the species tree. The results suggest that concatenated alignments of Next‐generation data that consist of many short sequences are robust to gene tree/species tree conflict when the goal is to determine the phylogenetic relationships among taxa.     

Mitochondrial DNA,morphology, and the phylogenetic relationships of Antarctic icefishes (Notothenioidei: Channichthyidae)

The Channichthyidae is a lineage of 16 species in the Notothenioidei, a clade of fishes that dominate Antarctic near-shore marine ecosystems with respect to both diversity and biomass. Among four published studies investigating channichthyid phylogeny, no two have produced the same tree topology, and no published study has investigated the degree of phylogenetic incongruence between existing molecular and morphological datasets. In this investigation we present an analysis of channichthyid phylogeny using complete gene sequences from two mitochondrial genes (ND2 and 16S) sampled from all recognized species in the clade. In addition, we have scored all 58 unique morphological characters used in three previous analyses of channichthyid phylogenetic relationships. Data partitions were analyzed separately to assess the amount of phylogenetic resolution provided by each dataset, and phylogenetic incongruence among data partitions was investigated using incongruence length difference (ILD) tests. We utilized a parsimony-based version of the Shimodaira-Hasegawa test to determine if alternative tree topologies are significantly different from trees resulting from maximum parsimony analysis of the combined partition dataset. Our results demonstrate that the greatest phylogenetic resolution is achieved when all molecular and morphological data partitions are combined into a single maximum parsimony analysis. Also, marginal to insignificant incongruence was detected among data partitions using the ILD. Maximum parsimony analysis of all data partitions combined results in a single tree, and is a unique hypothesis of phylogenetic relationships in the Channichthyidae. In particular, this hypothesis resolves the phylogenetic relationships of at least two species (Channichthys rhinoceratus and Chaenocephalus aceratus), for which there was no consensus among the previous phylogenetic hypotheses. The combined data partition dataset provides substantial statistical power to discriminate among alternative hypotheses of channichthyid relationships. These findings suggest the optimal strategy for investigating the phylogenetic relationships of channichthyids is one that uses all available phylogenetic data in analyses of combined data partitions.     

PHYLOGENETIC RELATIONSHIPS WITHIN THE COLONIAL VOLVOCALES (CHLOROPHYTA) INFERRED FROM rbcL GENE SEQUENCE DATA

The chloroplast-encoded large subunit of the ribulose-1, 5-bisphosphate carboxylase / oxygenase (rbcL) gene was sequenced from 20 species of the colonial Volvocales (the Volvacaceae, Goniaceae, and Tetrabaenaceae) in order to elucidate phylogenetic relationships within the colonial Volvocales. Eleven hundred twenty-eight base pairs in the coding regions of the (rbcL) gene were analyzed by the neighbor-joining (NJ) method using three kinds of distance estimations, as well as by the maximum parsimony (MP) method. A large group comprising all the anisogamous and oogamous volvocacean species was resolved in the MP tree as well as in the NJ trees based on overall and synonymous substitutions. In all the trees constructed, Basichlamys and Tetrabaena (Tetrabaenaceae) constituted a very robust phylogenetic group. Although not supported by high bootstrap values, the MP tree and the NJ tree based on nonsynonymous substitutions indicated that the Tetrabaenaceae is the sister group to the large group comprising the Volvocaceae and the Goniaceae. In addition, the present analysis strongly suggested that Pandorina and Astrephomene are monophyletic genera whereas Eudorina is nonmonophyletic. These results are essentially consistent with the results of the recent cladistic analyses of morphological data. However, the monophyly of the Volvocaceae previously supported by four morphological synapomorphies is found only in the NJ tree based on nonsynonymous substitutions (with very low bootstrap values). The genus Volvox was clearly resolved as a polyphyletic group with V. rousseletii Pocock separated from other species of Volvox in the rbcL gene comparisons, although this genus represents a monophyletic group in the previous morphological analyses. Furthermore, none of the rbcL gene trees supported the monophyly of the Goniaceae; Astrephomene was placed in various phylogenetic positions .     

Are phylogenies derived from family‐level supertrees robust for studies on macroecological patterns along environmental gradients?

Ecologists frequently use a supertree method to generate phylogenies in ecological studies. However, the robustness of research results based on phylogenies generated with a supertree method has not been well evaluated. Here, we use the angiosperm tree flora of North America as a model system to test the robustness of phylogenies generated with a supertree method for studies on the relationship between phylogenetic properties and environment, by comparing the relationship between phylogenetic metrics and environmental variables derived from a phylogeny reconstructed with a supertree method to that derived from a phylogeny resolved at species level. North America was divided into equal area quadrats of 12 100 km². Nine indices of phylogenetic structure were calculated for angiosperm tree assemblages in each quadrat using two phylogenies resolved at different levels (one resolved at the family level and the other resolved at the species level). Scores of phylogenetic indices were related to two major climatic variables (temperature and precipitation) using correlation and regression analyses. Scores of phylogenetic indices resulting from the two phylogenies are perfectly or nearly perfectly correlated. On average, there is no difference in the variation explained by the two climatic variables between scores of phylogenetic indices derived from the two phylogenies. Our study suggests that a phylogeny derived from a well resolved family-level supertree as backbone with genera and species attached to the backbone as polytomies is robust for studies investigating the relationship between phylogenetic structure and environment in biological assemblages at a broad spatial scale.