Systematics of the lizard family pygopodidae with implications for the diversification of Australian temperate biotas.

We conducted a phylogenetic study of pygopodid lizards, a group of 38 species endemic to Australia and New Guinea, with two major goals: to reconstruct a taxonomically complete and robustly supported phylogeny for the group and to use this information to gain insights into the tempo, mode, and timing of the pygopodid radiation. Phylogenetic analyses of mitochondrial DNA (mtDNA), nuclear DNA (nDNA), and previously published morphological data using parsimony, maximum likelihood, and Bayesian methods on the independent and combined three data sets yielded trees with similar and largely stable ingroup topologies. However, relationships among the six most inclusive and unambiguously supported clades (Aprasia, Delma, Lialis, Ophidiocephalus, Pletholax, and Pygopus) varied depending on data set analyzed. We used parametric bootstrapping to help us understand which of the three-branch schemes linking these six taxa was most plausible given our data. We conclude based on our results that the arrangement ((((Delma, Lialis)Pygopus)Pletholax)(Aprasia, Ophidiocephalus)) represents the best hypothesis of intergeneric relationships. A second major problem to arise in our study concerned the inability of our two outgroup taxa (Diplodactylus) to root trees properly; three different rooting locations were suggested depending upon analysis. This long-branch attraction problem was so severe that the outgroup branch also interfered with estimation of ingroup relationships. We therefore used the molecular clock method to root the pygopodid tree. Results of two independent molecular clock analyses (mtDNA and nDNA) converged upon the same root location (branch leading to Delma). We are confident that we have found the correct root because the possibility of our clock estimates agreeing by chance alone is remote given that there are 65 possible root locations (branches) on the pygopodid tree (approximately 1 in 65 odds). Our analysis also indicated that Delma fraseri is not monophyletic, a result supported by a parametric bootstrapping test. We elevated the Western Australian race, Delma f. petersoni, to species status (i.e., Delma petersoni) because hybridization and incomplete lineage sorting could be ruled out as potential causes of this paraphyletic gene tree and because D. grayii is broadly sympatric with its sister species D. fraseri. Climate changes over the past 23 million years, which transformed Australia from a wet, green continent to one that is largely dry and brown, have been suspected as playing a major role in the diversification of Australia's temperate biotas. Our phylogenetic analyses of pygopodid speciation and biogeography revealed four important findings consistent with this climate change diversification model: (1) our fossil-calibrated phylogeny shows that although some extant pygopodid lineages predate the onset of aridification, 28 of 33 pygopodid species included in our study seem to have originated in the last 23 million years; (2) relative cladogenesis tests suggest that several major clades underwent higher than expected rates of speciation; (3) our findings support earlier studies showing that speciation of mesic-adapted biotas in the southeastern and southwestern corners of Australia largely occurred within each of these regions between 12 and 23 million years ago as opposed to repeated dispersal between these regions; and (4) we have identified for the first time the existence of several pairs of sympatric sister species of lizards living in arid and semiarid ecosystems. These sympatric sister species seem to be younger than allopatric or parapatric sister-species pairs, which is not consistent with previous beliefs.

Jennings WB ，Pianka ER ，Donnellan S 《SYSTEMATIC BIOLOGY》

Is homoplasy or lineage sorting the source of incongruent mtdna and nuclear gene trees in the stiff-tailed ducks (Nomonyx-Oxyura)?

We evaluated the potential effects of homoplasy, ancestral polymorphism, and hybridization as obstacles to resolving phylogenetic relationships within Nomonyx-Oxyura stiff-tailed ducks (Oxyurinae; subtribe Oxyurina). Mitochondrial DNA (mtDNA) control region sequences from 94 individuals supported monophyly of mtDNA haplotypes for each of the six species and provided no evidence of extant incomplete lineage sorting or inter-specific hybridization. The ruddy ducks (O. j. jamaicensis,O. j. andina, O. j. ferruginea) are each others' closest relatives, but the lack of shared haplotypes between O. j. jamaicensis and O. j. ferruginea suggests long-standing historical isolation. In contrast, O. j. andina shares haplotypes with O. j. jamaicensis and O. j. ferruginea, which supports Todd's (1979) and Fjeldså's (1986) hypothesis that O. j. andina is an intergrade or hybrid subspecies of O. j. jamaicensis and O. j. ferruginea. Control region data and a much larger data set composed of approximately 8800 base pairs of mitochondrial and nuclear sequence for each species indicate that the two New World species, O. vittata and O. jamaicensis, branch basally within Oxyura. A clade of three Old World species (O. australis, O. maccoa, O. leucocephala) is well supported, but different loci and also different characters within the mtDNA data support three different resolutions of the Old World clade, yielding an essentially unresolved trichotomy. Fundamentally different factors limited the resolution of the mtDNA and nuclear gene trees. Gene trees for most nuclear loci were unresolved due to slow rates of mutation and a lack of informative variation, whereas uncertain resolution of the mtDNA gene tree was due to homoplasy. Within the mtDNA, approximately equal numbers of characters supported each of three possible resolutions. Parametric and nonparametric bootstrap analyses suggest that resolution of the mtDNA tree based on ~4300 bp per taxon is uncertain but that complete mtDNA sequences would yield a fully resolved gene tree. A short internode separating O. leucocephala from (O. australis, O. maccoa) in the best mtDNA tree combined with long terminal branches and substantial rate variation among nucleotide sites allowed the small number of changes occurring on the internode to be obscured by homoplasy in a significant portion of simulated data sets. Although most nuclear loci were uninformative, two loci supported a resolution of the Old World clade (O. maccoa, O. leucocephala) that is incongruent with the best mtDNA tree. Thus, incongruence between nuclear and mtDNA trees may be due to random sorting of ancestral lineages during the short internode, homoplasy in the mtDNA data, or both. The Oxyura trichotomy represents a difficult though likely common problem in molecular systematics. Given a short internode, the mtDNA tree has a greater chance of being congruent with the history of speciation because its effective population size (N(e)) is one-quarter that of any nuclear locus, but its resolution is more likely to be obscured by homoplasy. In contrast, gene trees for more slowly evolving nuclear loci will be difficult to resolve due to a lack of substitutions during the internode, and when resolved are more likely to be incongruent with the species history due to the stochastic effects of lineage sorting. We suggest that researchers consider first whether independent gene trees are adequately resolved and then whether those trees are congruent with the species history. In the case of Oxyura, the answer to both questions may be no. Complete mtDNA sequences combined with data from a very large number of nuclear loci may be the only way to resolve such trichotomies.

McCracken K ，Sorenson M 《SYSTEMATIC BIOLOGY》

Evolutionary relationships among the snakelike pygopodid lizards: a review of phylogenetic studies of an enigmatic Australian adaptive radiation.

Here, I review phylogenetic studies of the lizard family Pygopodidae, a group of 47 extant species that diversified in Australia and New Guinea. The goal of this study was to examine published phylogenetic and phylogenomic hypotheses on pygopodids to identify the strengths and weaknesses in our understanding of their phylogeny. Many parts of the pygopodid family tree are well established by multiple independent tree inferences including: (1) all multispecies genera (i.e., , , , , and ) are monophyletic groups; (2) the root of the pygopodid tree is located along the branch leading to the clade, thus showing that is the sister group to all other pygopodid genera; (3) the group, group, and several other groups of closely related species are demonstrated to be monophyletic entities; and (4) the monotypic is the sister lineage to the clade. Based on accumulated phylogenetic evidence, two taxonomic recommendations are given: merits generic status rather than being subsumed into as some earlier studies had suggested, and the monotypic should be recognized as a member of (following current practice) until future studies clarify its placement inside or outside the clade. One chronic problem with phylogenetic studies of pygopodids, which has limited the explanatory power of many tree hypotheses, concerns the undersampling of known species. Although the continual addition of newly described species, especially over the past two decades, has been a major reason for these taxon sampling gaps, deficits in species sampling for ingroups and/or outgroups in several studies of pygopodid species complexes has confounded the testing of some ingroup monophyly hypotheses. Ancient hybridization between non-sister lineages may also be confounding attempts to recover the relationships among pygopodids using molecular data. Indeed, such a phenomenon can explain at least five cases of mito-nuclear discordance and conflicts among trees based on nuclear DNA datasets. Another problem has been the lack of consensus on the relationships among most pygopodid genera, an issue that may stem from rapid diversification of these lineages early in the group's history. Despite current weaknesses in our understanding of pygopodid phylogeny, enough evidence exists to clarify many major and minor structural parts of their family tree. Accordingly, a composite tree for the Pygopodidae was able to be synthesized. This novel tree hypothesis contains all recognized pygopodid species and reveals that about half of the clades are corroborated by multiple independent tree hypotheses, while the remaining clades have less empirical support.

Jennings WB 《PeerJ》

Molecular phylogeography and systematics of the arid-zone members of the Egernia whitii (Lacertilia: Scincidae) species group.

We assembled a molecular phylogeny for the arid-zone members of the Egernia whitii species group to test Pianka's [Zoogeography and speciation of Australian desert lizards: an ecological perspective, Copeia (1972) 127-145] hypothesis that habitat specificity to the three major arid-zone vegetation communities is the primary cause of lizard speciation within the arid interior of Australia. This hypothesis predicts that species should exhibit phylogeographic structuring concordant with the major arid-zone vegetation types. Sequence data were obtained from four of the five arid-zone members of the E. whitii species group, and from across the ranges of the ecologically generalized E. inornata and E. multiscutata and the more specialized E. striata. We targeted a fragment (696 base pair (bp)) of the mitochondrial genome comprising the 3' half of the ND4 gene. We analysed the data using parsimony, maximum likelihood and Bayesian methods. Our phylogeny confirms the monophyly of the arid-zone members of the species group, although the phylogenetic relationships among species were not fully resolved. Although our topology does not support the recognition of the existing subspecies within E. multiscutata, there is a substantial phylogeographic break between South Australian/Victorian (Clade 1) and Western Australian (Clade 2) populations. We found considerable phylogeographic structure within E. inornata, with six major clades identified. However, these clades were not concordant with the distribution of habitat types in the arid-zone. Phylogeographic structure was also observed in the more specialized E. striata, although our analysis revealed close phylogenetic affinities between the sympatric species E. striata and E. kintorei. Shimodaira-Hasegawa topology tests were equivocal in regard to whether the phylogeographic structure within E. striata was in accordance with Pianka's predictions. Although our data failed to provide strong support for the suggestion that ecological and habitat factors are responsible for the diversification of arid-zone lizards, most E. inornata and E. striata populations had similar habitats, indicating that adaptation to particular habitats may have some role in the speciation of lizards in the Australian arid-zone.

Chapple DG ，Keogh JS ，Hutchinson MN 《MOLECULAR PHYLOGENETICS AND EVOLUTION》

Molecular phylogeny of the carnivora (mammalia): assessing the impact of increased sampling on resolving enigmatic relationships.

This study analyzed 76 species of Carnivora using a concatenated sequence of 6243 bp from six genes (nuclear TR-i-I, TBG, and IRBP; mitochondrial ND2, CYTB, and 12S rRNA), representing the most comprehensive sampling yet undertaken for reconstructing the phylogeny of this clade. Maximum parsimony and Bayesian methods were remarkably congruent in topologies observed and in nodal support measures. We recovered all of the higher level carnivoran clades that had been robustly supported in previous analyses (by analyses of morphological and molecular data), including the monophyly of Caniformia, Feliformia, Arctoidea, Pinnipedia, Musteloidea, Procyonidae + Mustelidae sensu stricto, and a clade of (Hyaenidae + (Herpestidae + Malagasy carnivorans)). All of the traditional "families," with the exception of Viverridae and Mustelidae, were robustly supported as monophyletic groups. We further have determined the relative positions of the major lineages within the Caniformia, which previous studies could not resolve, including the first robust support for the phylogenetic position of marine carnivorans (Pinnipedia) within the Arctoidea (as the sister-group to musteloids [sensu lato], with ursids as their sister group). Within the pinnipeds, Odobenidae (walrus) was more closely allied with otariids (sea lions/fur seals) than with phocids ("true" seals). In addition, we recovered a monophyletic clade of skunks and stink badgers (Mephitidae) and resolved the topology of musteloid interrelationships as: Ailurus (Mephitidae (Procyonidae, Mustelidae [sensu stricto])). This pattern of interrelationships of living caniforms suggests a novel inference that large body size may have been the primitive condition for Arctoidea, with secondary size reduction evolving later in some musteloids. Within Mustelidae, Bayesian analyses are unambiguous in supporting otter monophyly (Lutrinae), and in both MP and Bayesian analyses Martes is paraphyletic with respect to Gulo and Eira, as has been observed in some previous molecular studies. Within Feliformia, we have confirmed that Nandinia is the outgroup to all other extant feliforms, and that the Malagasy Carnivora are a monophyletic clade closely allied with the mongooses (Herpestidae [sensu stricto]). Although the monophyly of each of the three major feliform clades (Viverridae sensu stricto, Felidae, and the clade of Hyaenidae + (Herpestidae + Malagasy carnivorans)) is robust in all of our analyses, the relative phylogenetic positions of these three lineages is not resolvable at present. Our analyses document the monophyly of the "social mongooses," strengthening evidence for a single origin of eusociality within the Herpestidae. For a single caniform node, the position of pinnipeds relative to Ursidae and Musteloidea, parsimony analyses of data for the entire Carnivora did not replicate the robust support observed for both parsimony and Bayesian analyses of the caniform ingroup alone. More detailed analyses and these results demonstrate that outgroup choice can have a considerable effect on the strength of support for a particular topology. Therefore, the use of exemplar taxa as proxies for entire clades with diverse evolutionary histories should be approached with caution. The Bayesian analysis likelihood functions generally were better able to reconstruct phylogenetic relationships (increased resolution and more robust support for various nodes) than parsimony analyses when incompletely sampled taxa were included. Bayesian analyses were not immune, however, to the effects of missing data; lower resolution and support in those analyses likely arise from non-overlap of gene sequence data among less well-sampled taxa. These issues are a concern for similar studies, in which different gene sequences are concatenated in an effort to increase resolving power.

Flynn JJ ，Finarelli JA ，Zehr S ，Hsu J ，Nedbal MA ... -  《SYSTEMATIC BIOLOGY》