Coalescent histories on phylogenetic networks and detection of hybridization despite incomplete lineage sorting

Analyses of the increasingly available genomic data continue to reveal the extent of hybridization and its role in the evolutionary diversification of various groups of species. We show, through extensive coalescent-based simulations of multilocus data sets on phylogenetic networks, how divergence times before and after hybridization events can result in incomplete lineage sorting with gene tree incongruence signatures identical to those exhibited by hybridization. Evolutionary analysis of such data under the assumption of a species tree model can miss all hybridization events, whereas analysis under the assumption of a species network model would grossly overestimate hybridization events. These issues necessitate a paradigm shift in evolutionary analysis under these scenarios, from a model that assumes a priori a single source of gene tree incongruence to one that integrates multiple sources in a unifying framework. We propose a framework of coalescence within the branches of a phylogenetic network and show how this framework can be used to detect hybridization despite incomplete lineage sorting. We apply the model to simulated data and show that the signature of hybridization can be revealed as long as the interval between the divergence times of the species involved in hybridization is not too small. We reanalyze a data set of 106 loci from 7 in-group Saccharomyces species for which a species tree with no hybridization has been reported in the literature. Our analysis supports the hypothesis that hybridization occurred during the evolution of this group, explaining a large amount of the incongruence in the data. Our findings show that an integrative approach to gene tree incongruence and its reconciliation is needed. Our framework will help in systematically analyzing genomic data for the occurrence of hybridization and elucidating its evolutionary role.     

A species-level phylogenetic study of the Verbena complex (Verbenaceae) indicates two independent intergeneric chloroplast transfers

Two major impediments to infer plant phylogenies at inter- or intra- species level include the lack of appropriate molecular markers and the gene tree/species tree discordance. Both of these problems require more extensive investigations. One of the foci of this study is examining the phylogenetic utility of a combined chloroplast DNA dataset (>5.0kb) of seven non-coding regions, in comparison with that of a large fragment (ca. 3.0kb) of a low-copy nuclear gene (waxy), in a recent, rapidly diversifying group, the Verbena complex. The complex includes three very closely related genera, Verbena (base chromosome number x=7), Glandularia (x=5), and Junellia (x=10), comprising some 150 species distributed predominantly in South and North America. Our results confirm the inadequacy of non-coding cpDNA in resolving relationships among closely related species due to lack of variation, and the great potential of low-copy nuclear gene as source of variation. However, this study suggests that when both cpDNA and nuclear DNA are employed in low-level phylogenetic studies, cpDNA might be very useful to infer organelle evolutionary history (e.g., chloroplast transfer) and more comprehensively understand the evolutionary history of organisms. The phylogenetic framework of the Verbena complex resulted from this study suggests that Junellia is paraphyletic and most ancestral among the three genera; both Glandularia and Verbena are monophyletic and have been derived from within Junellia. Implications of this phylogenetic framework to understand chromosome number evolution and biogeography are discussed. Most interestingly, the comparison of the cpDNA and nuclear DNA phylogenies indicates two independent intergeneric chloroplast transfers, both from Verbena to Glandularia. One is from a diploid North American Verbena species to a polyploid North American Glandularia species. The other is more ancient, from the South American Verbena group to the common ancestor of a major Glandularia lineage, which has radiated subsequently in both South and North America. The commonly assumed introgressive hybridization may not explain the chloroplast transfers reported here. The underlying mechanism remains uncertain.     

Miocene and Pliocene dominated diversification of the lichen-forming fungal genus Melanohalea (Parmeliaceae, Ascomycota) and Pleistocene population expansions

ABSTRACT: BACKGROUND: Factors promoting diversification in lichen symbioses remain largely unexplored. While Pleistocene events have been important for driving diversification and affecting distributions in many groups, recent estimates suggest that major radiations within some genera in the largest clade of macrolichens (Parmeliaceae, Ascomycota) vastly predate the Pleistocene. To better understand the temporal placement and sequence of diversification events in lichens, we estimated divergence times in a common lichen-forming fungal genus, Melanohalea, in the Northern Hemisphere. Divergence times were estimated using both concatenated gene tree and coalescence-based multilocus species tree approaches to assess the temporal context of major radiation events within Melanohalea. In order to complement our understanding of processes impacting genetic differentiation, we also evaluated the effects of Pleistocene glacial cycles on population demographics of distinct Melanohalea lineages, differing in reproductive strategies. RESULTS: We found that divergence estimates, from both concatenated gene tree and coalescence-based multilocus species tree approaches, suggest that diversification within Melanohalea occurred predominantly during the Miocene and Pliocene, although estimated of divergence times differed by up to 8.3 million years between the two methods. These results indicate that, in some cases, taxonomically diagnostic characters may be maintained among divergent lineages for millions of years. In other cases, similar phenotypic characters among non-sister taxa, including reproductive strategies, suggest the potential for convergent evolution due to similar selective pressures among distinct lineages. Our analyses provide evidence of population expansions predating the last glacial maximum in the sampled lineages. These results suggest that Pleistocene glaciations were not inherently unfavorable or restrictive for some Melanohalea species, albeit with apparently different demographic histories between sexually and vegetatively reproducing lineages. CONCLUSIONS: Our results contribute to the understanding of how major changes during the Miocene and Pliocene have been important in promoting diversification within common lichen-forming fungi in the northern Hemisphere. Additionally, we provide evidence that glacial oscillations have influenced current population structure of broadly distributed lichenized fungal species throughout the Holarctic.     

Complex evolution in Arundinarieae (Poaceae: Bambusoideae): incongruence between plastid and nuclear GBSSI gene phylogenies

The monophyly of tribe Arundinarieae (the temperate woody bamboos) has been unequivocally recovered in previous molecular phylogenetic studies. In a recent phylogenetic study, 10 major lineages in Arundinarieae were resolved based on eight non-coding plastid regions, which conflicted significantly with morphological classifications both at the subtribal and generic levels. Nevertheless, relationships among and within the 10 lineages remain unclear. In order to further unravel the evolutionary history of Arundinarieae, we used the nuclear GBSSI gene sequences along with those of eight plastid regions for phylogenetic reconstruction, with an emphasis on Chinese species. The results of the plastid analyses agreed with previous studies, whereas 13 primary clades revealed in the GBSSI phylogeny were better resolved at the generic level than the plastid phylogeny. Our analyses also revealed many inconsistencies between the plastid DNA and the nuclear GBSSI trees. These results implied that the nuclear genome and the plastid genome had different evolutionary trajectories. The patterns of incongruence suggested that lack of informative characters, incomplete lineage sorting, and/or hybridization (introgression) could be the causes. Seven putative hybrid species were hypothesized, four of which are discussed in detail on the basis of topological incongruence, chromosome numbers, morphology, and distribution patterns, and those taxa probably resulted from homoploid hybrid speciation. Overall, our study indicates that the tribe Arundinarieae has undergone a complex evolution.     

Dealing with incongruence in the quest for the species tree: a case study from the orchid genus Satyrium

A species tree was reconstructed for the mainly African terrestrial orchid genus Satyrium. Separate phylogenetic analysis of both plastid and ribosomal nuclear DNA sequences for 63 species, revealed extensive topological conflict. Here we describe a detailed protocol to deal with incongruence involving three steps: identifying incongruence and testing its significance, assessing the cause of incongruence, and reconstructing the species tree. The Incongruence Length Difference test revealed that many cases of incongruence were non-significant. For the remaining significant cases, results from taxon jack-knifing experiments and parametric bootstrap suggested that non-biological artefacts such as sparse taxon sampling and long-branch attraction could be excluded as causes for the observed incongruence. In order to evaluate biological causes, such as orthology/paralogy conflation, lineage sorting, and hybridization, the number of events was counted that needs to be invoked a-posteriori to explain the observed pattern. In most cases where incongruence was significant, this resulted in a similar number of events for each of these different causes. Only for the three species from south east Asia, that form a monophyletic clade, hybridization was favoured over the alternative causes. This conclusion is based on the large number of events that needs to be invoked, in order for either orthology/paralogy conflation or lineage sorting to have been the cause of the incongruence+morphological evidence. The final species tree presented here is the product of the combined analysis of plastid and ITS sequences for all non-incongruent species and a-posteriori grafting of the incongruent clades or accessions onto the tree.     

Molecular phylogeny of Cotoneaster (Rosaceae) inferred from nuclear ITS and multiple chloroplast sequences

Cotoneaster Medik. (Rosaceae, Maloideae) is distributed in Europe, North Africa, and temperate areas of Asia except Japan. Members of the genus exhibit considerable morphological variation. The infrageneric classification is also obscured by polyploidy, hybridization, and apomixis. In this study, phylogenetic analyses were conducted to test infrageneric classifications of this genus using DNA sequence data from the nuclear ITS (nrITS) region and three chloroplast intergenic spacer regions. Maximum parsimony and Bayesian inference analyses of both datasets agreed with the two sections/subgenera of Koehne’s classification system, and suggested that four subsections (Microphylli, Chaenopetalum, Adpressi, and Cotoneaster) and the series of Koehne’s classification system were non-monophyletic. The incongruence length difference test indicated that the nrITS and cpDNA datasets were significantly incongruent (P = 0.001), and the placement of 14 species was discordant in phylogenetic trees derived from the two datasets. Within Cotoneaster, hybridization was indicated to be an important factor contributing to the incongruence between the nrITS and cpDNA data. By mapping nine morphological characters onto the combined nrITS–cpDNA phylogenetic tree, we inferred that a deciduous habit, glabrous fruit, white anthers, erect and light pink petals, and white filaments are plesiomorphic character states in Cotoneaster.     

Chloroplast DNA inheritance and variation in <Emphasis Type="Italic">Leucadendron</Emphasis> species (Proteaceae) as revealed by PCR-RFLP

The inheritance of chloroplast DNA (cpDNA) in Leucadendron species was studied by polymerase chain reaction–restriction fragment length polymorphism (PCR–RFLP) analysis. A total of 100 progeny from five interspecific crosses involving seven parental species were tested, and all progeny exhibited the cpDNA restriction fragment pattern of the female parent, indicating that cpDNA in Leucadendron is maternally inherited. PCR-RFLP was also employed to study cpDNA variation among 21 Leucadendron species. Parsimony analysis using a heuristic search resulted in a phylogenetic tree that showed limited agreement to the taxonomic classification of Leucadendron species, based on morphological characteristics. The incongruence between cpDNA phylogenetic and taxonomic groupings in Leucadendron may be due to reticulate evolution involving a combination of hybridization and introgression, convergent evolution and/or lineage sorting at the interspecific, intersubsectional and intersectional levels.     

Allopolyploid origin of Cardamine asarifolia (Brassicaceae): incongruence between plastid and nuclear ribosomal DNA sequences solved by a single-copy nuclear gene

Interspecific hybridization and polyploidization have played central roles in plant diversification. However, technical difficulties in the analyses of low-copy genes have limited the study of the origins of hybrid and polyploid plants. Here, we present a phylogenetic analysis of the hexaploid Cardamine asarifolia, distributed in the southern European Alps and northern Apennines. Our study included all relevant taxa of the genus found in Europe. A marked discrepancy was revealed between the trnL-trnF region of cpDNA and internal transcribed spacer (nrDNA ITS) sequences. To solve the incongruence, we sequenced a single-copy nuclear CHS gene (chalcone synthase) using a novel method to design homoeologue-specific PCR primers to bypass artefacts caused by artificial recombination of homoeologues during PCR and/or cloning. Three homoeologues were isolated from C. asarifolia, providing evidence for its allopolyploid origin. One homoeologue, showing the same phylogenetic position as the ITS sequences, most likely originated from an extinct parent. Furthermore, we documented recurrent polytopic hybridizations between C. asarifolia and diploid C. amara. The allohexaploidization and the following hybridization with a diploid species exemplify the ongoing dynamic processes of speciation in the genus Cardamine.     

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.     

A biogeographic–ecological approach to disentangle reticulate evolution in the Triatoma phyllosoma species group (Heteroptera: Triatominae), vectors of Chagas disease

Introgression and incomplete lineage sorting (ILS) are two of the main sources of gene‐tree incongruence; both can confound the assessment of phylogenetic relationships among closely related species. The Triatoma phyllosoma species group is a clade of partially co‐distributed and cross‐fertile Chagas disease vectors. Despite previous efforts, the phylogeny of this group remains unresolved, largely because of substantial gene‐tree incongruence. Here, we sequentially address introgression and ILS to provide a robust phylogenetic hypothesis for the T. phyllosoma species group. To identify likely instances of introgression prior to molecular scrutiny, we assessed biogeographic data and information on fertility of inter‐specific crosses. We first derived a few explicit hybridization hypotheses by considering the degree of spatial overlap within each species pair. Then, we assessed the plausibility of these hypotheses in the light of each species pair's cross‐fertility. Using this contextual information, we evaluated mito‐nuclear (cyt b, ITS‐2) gene‐tree incongruence and found evidence suggesting introgression within two species pairs. Finally, we modeled ILS using a Bayesian multispecies coalescent approach and either (a) a “complete” dataset with all the specimens in our sample, or (b) a “filtered” dataset without putatively introgressed specimens. The “filtered tree” had higher posterior‐probability support, as well as more plausible topology and divergence times, than the “complete tree.” Detecting and filtering out introgression and modeling ILS allowed us to derive an improved phylogenetic hypothesis for the T. phyllosoma species group. Our results illustrate how biogeographic and ecological‐reproductive contextual information can help clarify the systematics and evolution of recently diverged taxa prone to introgression and ILS.     

Molecular phylogenetics of the sexually deceptive orchid genus Ophrys (Orchidaceae) based on nuclear and chloroplast DNA sequences.

We present a phylogenetic analysis of the major lineages of the sexually deceptive orchid genus Ophrys based on nuclear ribosomal (nr) DNA (internal transcribed spacer region) and noncoding chloroplast (cp) DNA (trnL-trnF region) sequences. Sequence divergence within and among major Ophrys lineages was low for both nrDNA and cpDNA sequences. Separate analyses resulted in similar but poorly resolved trees. An incongruence length difference test revealed that nrDNA and cpDNA data sets were not incongruent. A combined analysis resulted in a better-resolved phylogenetic hypothesis of relationships among the major Ophrys lineages. Our data strongly support a division of Ophrys into two groups. These groups do not correspond to the earlier proposed sections Euophrys and Pseudophrys and are thus in conflict with traditional classifications. Our results support a well-resolved monophyletic group that contains the geographically widespread O. bombyliflora, O. speculum, O. tenthredinifera, and the O. fusca-lutea lineage. Relationships in the other group are poorly resolved. Based on our observations that taxa with identical sequences at presumably rapidly evolving loci clearly differ in floral morphology, we hypothesize that the diversity in the genus Ophrys is the result of a recent radiation in this orchid lineage.     

A scalable and flexible approach for investigating the genomic landscapes of phylogenetic incongruence

Analyses of DNA sequence datasets have repeatedly revealed inconsistencies in phylogenetic trees derived with different data. This is termed phylogenetic incongruence, and may arise from a methodological failure of the inference process or from biological processes, such as horizontal gene transfer, incomplete lineage sorting, and introgression. To better understand patterns of incongruence, we developed a method (PartFinder) that uses likelihood ratios applied to sliding windows for visualizing tree-support changes across genome-sequence alignments, allowing the comparative examination of complex phylogenetic scenarios among many species. As a pilot, we used PartFinder to investigate incongruence in the Homo-Pan-Gorilla group as well as Platyrrhini using high-quality bacterial artificial chromosome (BAC)-derived sequences as well as assembled whole-genome shotgun sequences. Our simulations verified the sensitivity of PartFinder, and our results were comparable to other studies of the Homo-Pan-Gorilla group. Analyses of the whole-genome alignments reveal significant associations between support for the accepted species relationship and specific characteristics of the genomic regions, such as GC-content, alignment score, exon content, and conservation. Finally, we analyzed sequence data generated for five platyrrhine species, and found incongruence that suggests a polytomy within Cebidae, in particular. Together, these studies demonstrate the utility of PartFinder for investigating the patterns of phylogenetic incongruence.     

Chloroplast DNA evidence for reticulate evolution in Eucalyptus (Myrtaceae)

Four highly differentiated chloroplast DNA (cpDNA) lineages were identified in the forest tree species Eucalyptus globulus Labill. (Myrtaceae) in Australia using restriction site polymorphisms from Southern analysis. The cpDNA variation did not conform with ssp. boundaries, yet there was a strong geographical pattern to the distribution of the lineages. One lineage (C) was geographically central and widespread, whereas the other three lineages were found in peripheral populations: Western (W), Northern (N) and Southern (S). Thirteen haplotypes were detected in E. globulus , seven of which belonged to clade C. At least three of the cpDNA lineages (C, N and S) were shared extensively with other species. On the east coast of the island of Tasmania, there was a major north–south difference in cpDNA in the virtually continuous distribution of E. globulus . Northern populations harboured haplotypes from clade C while southeastern populations harboured a single haplotype from clade S. This difference was also reflected in several co-occurring endemic species. It is argued that the extensive cpDNA differentiation within E. globulus is likely to originate from interspecific hybridization and 'chloroplast capture' from different species in different parts of its range. Superficially, this hybridization is not evident in taxonomic traits; however, large-scale common garden experiments have revealed a steep cline in quantitative genetic variation that coincides with the haplotype transition in Tasmania. Our cpDNA results provide the strongest evidence to date that hybridization has had a widespread impact on a eucalypt species and indicate that reticulate evolution may be occurring on an unappreciated scale in Eucalyptus .     

DIVERSIFICATION AND GENE FLOW IN NASCENT LINEAGES OF ISLAND AND MAINLAND NORTH AMERICAN TREE SQUIRRELS (TAMIASCIURUS)

Pleistocene climate cycles and glaciations had profound impacts on taxon diversification in the Boreal Forest Biome. Using population genetic analyses with multilocus data, we examined diversification, isolation, and hybridization in two sibling species of tree squirrels (Tamiasciurus douglasii and Tamiasciurus hudsonicus) with special attention to the geographically and genetically enigmatic population of T. hudsonicus on Vancouver Island, Canada. The two species differentiated only about 500,000 years ago, in the Late Pleistocene. The island population is phylogenetically nested within T. hudsonicus according to our nuclear analysis but within T. douglasii according to mitochondrial DNA. This conflict is more likely due to historical hybridization than to incomplete lineage sorting, and it appears that bidirectional gene flow occurred between the island population and both species on the mainland. This interpretation of our genetic analyses is consistent with our bioclimatic modeling, which demonstrates that both species were able to occupy this region throughout the Late Pleistocene. The divergence of the island population 40,000 years ago suggests that tree squirrels persisted in a refugium on Vancouver Island at the last glacial maximum, 20,000 years ago. Our observations demonstrate how Pleistocene climate change and habitat shifts have created incipient divergence in the presence of gene flow.     

Molecular phylogeny of Osmanthus(Oleaceae)based on non-coding chloroplast and nuclear ribosomal internal transcribed spacer regions

The phylogenetic relationships of Osmanthus Lour.were investigated using the nuclear ribosomal internal transcribed spacer(ITS)regions and non-coding chloroplast regions(psbA-trnH,trnL-F).The two datasets support the conclusion that Osmanthus is polyphyletic.with some species of the subtribe Oleinae nested within Osmanthus.Osmanthus didymopetalus P S.Green is nested within the clade formed by species of section Osmanthus in two trees.Osmanthus attenuatus P.S.Green.O.Funnanensis P.S.Green,and O gracilinervis R.L.Lu of traditional section Osmanthus are clearly divergent from other accessions,and do not form a monophyletic group with other Osmanthus accessions.Osmanthus marginatus Hemsl.is embedded in the clade formed by species of section Osmanthus in the ITS tree.In cpDNA trees all species of section Osmanthus are placed in the large clade and all species of section Leiolea formed a group.The taxonomic incongruence among trees for ITS and cpDNA indicate hybridization.as introgression may have occurred among some species of sections Osmanthus and Leiolea.Phylogeny of Osmanthus is discussed in light of molecular and morphological data,and a revised infrageneric classification with three sections(Leiolea,Siphosmanthu,and Osmanthus)is presented.The section Linocieroides is abandoned and united with section Osmanthus.     

Detecting hybrid speciation in the presence of incomplete lineage sorting using gene tree incongruence: A model

The application of phylogenetic inference methods, to data for a set of independent genes sampled randomly throughout the genome, often results in substantial incongruence in the single-gene phylogenetic estimates. Among the processes known to produce discord between single-gene phylogenies, two of the best studied in a phylogenetic context are hybridization and incomplete lineage sorting. Much recent attention has focused on the development of methods for estimating species phylogenies in the presence of incomplete lineage sorting, but phylogenetic models that allow for hybridization have been more limited. Here we propose a model that allows incongruence in single-gene phylogenies to be due to both hybridization and incomplete lineage sorting, with the goal of determining the contribution of hybridization to observed gene tree incongruence in the presence of incomplete lineage sorting. Using our model, we propose methods for estimating the extent of the role of hybridization in both a likelihood and a Bayesian framework. The performance of our methods is examined using both simulated and empirical data.     

Testing hybridization hypotheses based on incongruent gene trees

Hybridization is an important evolutionary mechanism in plants and has been increasingly documented in animals. Difficulty in reconstruction of reticulate evolution, however, has been a long-standing problem in phylogenetics. Consequently, hybrid speciation may play a major role in causing topological incongruence between gene trees. The incongruence, in turn, offers an opportunity to detect hybrid speciation. Here we characterized certain distinctions between hybridization and other biological processes, including lineage sorting, paralogy, and lateral gene transfer, that are responsible for topological incongruence between gene trees. Consider two incongruent gene trees with three taxa, A, B, and C, where B is a sister group of A on gene tree 1 but a sister group of C on gene tree 2. With a theoretical model based on the molecular clock, we demonstrate that time of divergence of each gene between taxa A and C is nearly equal in the case of hybridization (B is a hybrid) or lateral gene transfer, but differs significantly in the case of lineage sorting or paralogy. After developing a bootstrap test to test these alternative hypotheses, we extended the model and test to account for incongruent gene trees with numerous taxa. Computer simulation studies supported the validity of the theoretical model and bootstrap test when each gene evolved at a constant rate. The computer simulation also suggested that the model remained valid as long as the rate heterogeneity was occurring proportionally in the same taxa for both genes. Although the model could not test hypotheses of hybridization versus lateral gene transfer as the cause of incongruence, these two processes may be distinguished by comparing phylogenies of multiple unlinked genes.     

Analysis of plastid and nuclear DNA data in plant phylogenetics——evaluation and improvement

Correct combination of plastid(cp)and nuclear(nr)DNA data for plant phylogenetic reconstructions is not a new issue,but with an increasing number of nrDNA loci being used,it is of ever greater practical concern.For accurately reconstructing the phylogeny and evolutionary history of plant groups,correct treatment of phylogenetic incongruence is a vital step in the proper analysis of cpDNA and nrDNA data.We first evaluated the current status of analyzing cpDNA and nrDNA data by searching all articles published in the journal Systematic Botany between 2005 and 2011.Many studies combining cpDNA and nrDNA data did not rigorously assess the combinability of the data sets,or did not address in detail possible reasons for incongruence between the two data sets.By reviewing various methods,we outline a procedure to more accurately analyze and/or combine cpDNA and nrDNA data,which includes four steps:identifying significant incongruence,determining conflicting taxa,providing possible interpretations for incongruence,and reconstructing the phylogeny after treating incongruence.Particular attention is given to explanation of the cause of incongruence.We hope that our procedure will help raise awareness of the importance of rigorous analysis and help identify the cause of incongruence before combining cpDNA and nrDNA data.     

Concordance between Phylogeographical and Biogeographical Patterns in the Brazilian Cerrado: Diversification of the Endemic Tree Dalbergia miscolobium (Fabaceae)

Few studies have addressed the phylogeography of species of the Cerrado, the largest savanna biome of South America. Here we aimed to investigate the phylogeographical structure of Dalbergia miscolobium, a widespread tree from the Cerrado, and to verify its concordance with plant phylogeographical and biogeographical patterns so far described. A total of 287 individuals from 32 populations were analyzed by sequencing the trnL intron of the chloroplast DNA and the internal transcribed spacer of the nuclear ribosomal DNA. Analysis of population structure and tests of population expansion were performed and the time of divergence of haplotypes was estimated. Twelve and 27 haplotypes were identified in the cpDNA and nrDNA data, respectively. The star-like network configuration and the mismatch distributions indicated a recent spatial and demographic expansion of the species. Consistent with previous tree phylogeographical studies of Cerrado trees, the cpDNA also suggested a recent expansion towards the southern Cerrado. The diversity of D. miscolobium was widespread but high levels of genetic diversity were found in the Central Eastern and in the southern portion of Central Western Cerrado. The combined analysis of cpDNA and nrDNA supported a phylogeographic structure into seven groups. The phylogeographical pattern showed many concordances with biogeographical and phylogeographical studies in the Cerrado, mainly with the Cerrado phytogeographic provinces superimposed to our sampling area. The data reinforced the uniqueness of Northeastern and Southeastern Cerrados and the differentiation between Eastern and Western Central Cerrados. The recent diversification of the species (estimated between the Pliocene and the Pleistocene) and the ‘genealogical concordances’ suggest that a shared and persistent pattern of species diversification might have been present in the Cerrado over time. This is the first time that an extensive ‘genealogical concordance’ between phylogeographic and phytogeographic patterns is shown for the Cerrado biome.     

Incongruence among mitochondrial,chloroplast and nuclear gene trees in <Emphasis Type="Italic">Pinus</Emphasis> subgenus <Emphasis Type="Italic">Strobus</Emphasis> (Pinaceae)

Introgression has been considered to be one of main factors leading to phylogenetic incongruence among different datasets at lower taxonomic levels. In the plants of Pinaceae, the mtDNA, cpDNA, and nuclear DNA (nrDNA) may have different evolutionary histories through introgression because they are inherited maternally, paternally and biparentally, respectively. We compared mtDNA, cpDNA, and two low-copy nrDNA phylogenetic trees in the genus Pinus subgenus Strobus, in order to detect unknown past introgression events in this group. nrDNA trees were mostly congruent with the cpDNA tree, and supported the recent sectional and subsectional classification system. In contrast, mtDNA trees split the members of sect. Quinquefoliae into two groups that were not observed in the other gene trees. The factors constituting incongruence may be divided into the following two categories: the different splits within subsect. Strobus, and the non-monophyly of subsect. Gerardianae. The former was hypothesized to have been caused by the past introgression of cpDNA, mtDNA or both between Eurasian and North American species through Beringia. The latter was likely caused by the chimeric structure of the mtDNA sequence of P. bungeana, which might have originated through past hybridization, or through a horizontal transfer event and subsequent recombination. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.