Ancient allopolyploid speciation in Geinae (Rosaceae): evidence from nuclear granule-bound starch synthase (GBSSI) gene sequences

A nuclear low-copy gene phylogeny provides strong evidence for the hybrid origin of seven polyploid species in Geinae (Rosaceae). In a gene tree, alleles at homologous loci in an allopolyploid species are expected to be sisters to orthologues in the ancestral taxa rather than to each other. Alleles at a duplicated locus in an autopolyploid, however, are expected to be more closely related to each other than they are to any orthologous copies in closely related species. We cloned and sequenced about 1.9 kilobases from the 5' end of the GBSSI-1 gene from two diploid, one tetraploid, and six hexaploid species. Each of the three loci in the hexaploid species forms a separate group, two of which are more closely related to copies in other species than they are to each other. This finding indicates that the hexaploid lineage evolved through two consecutive allopolyploidization events. Based on the GBSSI-1 gene tree, we hypothesized that there was an initial hybridization between a diploid species from the ancestral lineage of Coluria and Waldsteinia and an unknown diploid species to form the tetraploid Geum heterocarpum lineage. Backcrossing of G. heterocarpum with a representative of the unknown diploid lineage then resulted in a hexaploid lineage that has radiated considerably since its origin, comprising at least 40 extant species with various morphologies. A penalized likelihood analysis indicated that Geinae may be about 17 million years old, implying that the hypothesized allopolyploid speciation events are relatively ancient. Six of the 22 cloned Geinae GBSSI-1 copies in this study, which all are duplicate copies in polyploid taxa, may have become pseudogenes. We compared the GBSSI-1 phylogeny with one from chloroplast data and explored implications for the evolution of some fruit characters.     

Relationships between gene trees and species trees

It is well known that a phylogenetic tree (gene tree) constructed from DNA sequences for a genetic locus does not necessarily agree with the tree that represents the actual evolutionary pathway of the species involved (species tree). One of the important factors that cause this difference is genetic polymorphism in the ancestral species. Under the assumption of neutral mutations, this problem can be studied by evaluating the probability (P) that a gene tree has the same topology as that of the species tree. When one gene (allele) is used from each of the species involved, the probability can be expressed as a simple function of Ti = ti/(2N), where ti is the evolutionary time measured in generations for the ith internodal branch of the species tree and N is the effective population size. When any of the Ti's is less than 1, the probability P becomes considerably less than 1.0. This probability cannot be substantially increased by increasing the number of alleles sampled from a locus. To increase the probability, one has to use DNA sequences from many different loci that have evolved independently of each other.      

Chloroplast DNA phylogeny, reticulate evolution, and biogeography of Paeonia (Paeoniaceae)

The coding region of the mat K gene and two intergenic spacers, psb A-trn H and trn L(UAA)-trn F(GAA), of cpDNA were sequenced to study phylogenetic relationships of 32 Paeonia species. In the psb A-trn H intergenic spacer, short sequences bordered by long inverted repeats have undergone inversions that are often homoplasious mutations. Insertions/deletions found in the two intergenic spacers, mostly resulting from slipped-strand mispairing, provided relatively reliable phylogenetic information. The mat K coding region, evolving more rapidly than the trnL-trn F spacer and more slowly than the psb A-trn H spacer, produced the best resolved phylogenetic tree. The mat K phylogeny was compared with the phylogeny obtained from sequences of internal transcribed spacers (ITS) of nuclear ribosomal DNA. A refined hypothesis of species phylogeny of section Paeonia was proposed by considering the discordance between the nuclear and cpDNA phylogenies to be results of hybrid speciation followed by inheritance of cpDNA of one parent and fixation of ITS sequences of another parent. The Eurasian and western North American disjunct distribution of the genus may have resulted from interrruption of the continuous distribution of ancestral populations of extant peony species across the Bering land bridge during the Miocene. Pleistocene glaciation may have played an important role in triggering extensive reticulate evolution within section Paeonia and shifting distributional ranges of both parental and hybrid species.     

Evolution of parental ITS regions of nuclear rDNA in allopolyploid Aegilops (Poaceae) species

The genus Aegilops comprises approximately 25 diploid, tetraploid and hexaploid species, in which the genome types of all allopolyploids involve either U or D genome, or both of them. The internal transcribed spacer (ITS) region of 18S-26S nuclear ribosomal DNA (rDNA) from 11 allopolyploid species and 7 related diploid species in the genus were directly sequenced by pooled PCR products. Phylogenetic analyses for tracing evolutionary patterns of parental rDNA in allopolyploid species were performed using the neighbor-joining method. The D genome involved tree included three clades (CC-DDCC, DDMM-DDMMSS-DDMMUU, and MM-MhMh-DDNN), but did not include Ae. squarrosa (DD). It indicated that the rDNA of ancestral D genome had been somewhat differentiated in allopolyploids. The U genome involved tree showed that the allopolyploids and their common ancestor, Ae. umbellulata, formed a clade, suggesting that rDNA in UUMM and UUSS genomes has been homogenizing toward that of ancestral U genome. The phylogenetic pattern of U genome based on ITS sequences also supported the "pivotal-differential" hypothesis.     

Boreotropical migration explains hybridization between geographically distant lineages in the pantropical clade Sideroxyleae (Sapotaceae)

To determine whether the fragmented pantropical distribution of present day Sideroxyleae primarily is the result of long-distance dispersals or represents the remnants of a once continuous distribution in the northern hemisphere, the boreotropical flora, we used phylogenetic analyses of chloroplast and nuclear ribosomal DNA data, Bayesian molecular dating, and Bayesian estimation of ancestral areas. Incongruence between the two data sets was examined with a nuclear low copy gene phylogeny to discover any occurrences of reticulate evolution. The Pacific clade Nesoluma was shown to have two distinct copies of the nuclear low copy gene AAT, one from an African and one from an American ancestral lineage, indicating that it is of allopolyploid origin. We conclude that Sideroxyleae, including the ancestral lineages of Nesoluma, were part of the boreotropical flora and entered the New World via the north Atlantic land bridge. We also suggest that the distribution of extant species resulted from the cooling climate at the end of the Eocene. Sideroxylon oxyacanthum is shown not to belong in the group, but in Chrysophylloideae. A classification reflecting phylogenetic relationships, as well as new combinations for the species in Nesoluma under Sideroxylon, is presented.     

Extracting species trees from complex gene trees: reconciled trees and vertebrate phylogeny

Paralogy is a pervasive problem in trying to use nuclear gene sequences to infer species phylogenies. One strategy for dealing with this problem is to infer species phylogenies from gene trees using reconciled trees, rather than directly from the sequences themselves. In this approach, the optimal species tree is the tree that requires the fewest gene duplications to be invoked. Because reconciled trees can identify orthologous from paralogous sequences, there is no need to do this prior to the analysis. Multiple gene trees can be analyzed simultaneously; however, the problem of nonuniform gene sampling raises practical problems which are discussed. In this paper the technique is applied to phylogenies for nine vertebrate genes (aldolase, alpha-fetoprotein, lactate dehydrogenase, prolactin, rhodopsin, trypsinogen, tyrosinase, vassopressin, and Wnt-7). The resulting species tree shows much similarity with currently accepted vertebrate relationships.     

Evolutionary dynamics of Waxy and the origin of hexaploid Spartina species (Poaceae)

We investigated the evolutionary dynamics of duplicated copies of the granule-bound starch synthase I gene (GBSSI or Waxy) within polyploid Spartina species. Molecular cloning, sequencing, and phylogenetic analyses revealed incongruences between the expected species phylogeny and the inferred gene trees. Some genes within species were more divergent than expected from ploidy level alone, suggesting the existence of paralogous sets of Waxy loci in Spartina. Phylogenetic analyses indicate that this paralogy originated from a duplication that occurred prior to the divergence of Spartina from other Chloridoideae. Gene tree topologies revealed three divergent homoeologous sequences in the hexaploid S. alterniflora that are consistent with the proposal of an allopolyploid origin of the hexaploid clade. Waxy sequences differ in insertion–deletion events in introns, which may be used to diagnose gene copies. Both paralogous and homoeologous coding regions appear to evolving under selective constraints.     

Evidence of Gene Conversion Events Between Paralogous Sequences Produced by Tetraploidization in Salmoninae Fish

We investigated the occurrence of gene conversions between paralogous sequences of Salmoninae derived from ancestral tetraploidization and their effect on the evolutionary history of DNA sequences. A microsatellite with long flanking regions (750 bp) including both coding and noncoding sequences was analyzed. Microsatellite size polymorphism was used to detect the alleles of both paralogous counterparts and infer linkage arrangement between loci. DNA sequencing of seven Salmoninae species revealed that paralogous sequences were highly differentiated within species, especially for noncoding regions. Ten gene conversion events between paralogous sequences were inferred. While these events appears to have homogenized regions of otherwise highly differential paralogous sequences, they amplified the differentiation among orthologous sequences. Their effects were larger on coding than on noncoding regions. As a consequence, noncoding sequences grouped by orthologous lineages in phylogenetic trees, whereas coding regions grouped by taxa. Based upon these results, we present a model showing how gene conversion events may also result in the PCR amplification of nonorthologous sequences in different taxa, with obvious complications for phylogenetic inferences, comparative mapping, and population genetic studies. Received: 11 October 2000 / Accepted: 18 September 2001     

Perils of paralogy: using HSP70 genes for inferring organismal phylogenies

Conserved genes have found their way into the mainstream of molecular systematics. Many of these genes are members of multigene families. A difficulty with using single genes of multigene families for phylogenetic inference is that genes from one species may be paralogous to those from another taxon. We focus attention on this problem using heat shock 70 (HSP70) genes. Using polymerase chain reaction techniques with genomic DNA, we isolated and sequenced 123 distinct sequences from 12 species of sharks. Phylogenetic analysis indicated that the sequences cluster with constituitively expressed cytoplasmic heat shock-like genes. Three highly divergent gene clades were sampled. A number of similar sequences were sampled from each species within each distinct gene clade. Comparison of published species trees with an HSP70 gene tree inferred using Bayesian phylogenetic analysis revealed several cases of gene duplication and differential sorting of gene lineages within this group of sharks. Gene tree parsimony based on the objective criteria of duplication and losses showed that previously published hypotheses of species relationships and two novel hypothesis based on Bayesian phylogenetics were concordant with the history of HSP70 gene duplication and loss. By contrast, two published hypotheses based on morphological data were not significantly different from the null hypothesis of a random association between species relatedness and the HSP70 gene tree. These results suggest that gene tree parsimony using data from multigene families can be used for inferring species relationships or testing published alternative hypotheses. More importantly, the results suggest that systematic studies relying on phylogenetic inferences from HSP70 genes may by plagued by unrecognized paralogy of sampled genes. Our results underscore the distinction between gene and species trees and highlight an underappreciated source of discordance between gene trees and organismal phylogeny, i.e., unrecognized paralogy of sampled genes.     

指环虫属的早期辐射及其与宿主鲤科鱼类的协同进化关系

本研究利用28SrDNAC1-D2区序列分析采自鲤科鱼类中6亚科宿主和寄生在花鲈、梅花鲈上的共17种指环虫的系统发育关系。同时,通过比较宿主鲤科鱼类与指环虫的系统发育树,检验指环虫与其宿主是否存在协同进化关系。结果表明:17种指环虫形成5个进化支(Clade),其中寄生在团头鲂(亚科)和鲢、鳙(鲢亚科)上的6种指环虫聚为一支(Clade1),而它们的宿主鱼类在系统发育分析中也表现为近缘关系;寄生在鲮鱼(野鲮亚科)上的D.quanfami(Clade5)位于系统树最基部,鲫鱼和鲤鱼(鲤亚科)的寄生指环虫处在系统树的次基部位置,而鲤亚科与野鲮亚科组成的姐妹群在宿主系统树上同样处在基部位置,寄生虫和宿主在进化上较为原始的地位得到了很好地相互印证。因而,本研究首次利用分子系统学手段分析指环虫属远缘物种间的系统关系,揭示了指环虫属与宿主鱼类之间存在协同进化关系。另外,本研究首次发现,野鲮亚科鱼类也可能是指环虫类的早期宿主,这与先前认为鲤亚科鱼类为指环虫类的祖先宿主的推测有所不同。     

Target enrichment improves phylogenetic resolution in the genus Zanthoxylum (Rutaceae) and indicates both incomplete lineage sorting and hybridization events

Background and Aims Zanthoxylum is the only pantropical genus within Rutaceae, with a few species native to temperate eastern Asia and North America. Efforts using Sanger sequencing failed to resolve the backbone phylogeny of Zanthoxylum. In this study, we employed target-enrichment high-throughput sequencing to improve resolution. Gene trees were examined for concordance and sectional classifications of Zanthoxylum were evaluated. Off-target reads were investigated to identify putative single-copy markers for bait refinement, and low-copy markers for evidence of putative hybridization events.MethodsA custom bait set targeting 354 genes, with a median of 321 bp, was designed for Zanthoxylum and applied to 44 Zanthoxylum species and one Tetradium species as the outgroup. Illumina reads were processed via the HybPhyloMaker pipeline. Phylogenetic inferences were conducted using coalescent and maximum likelihood methods based on concatenated datasets. Concordance was assessed using quartet sampling. Additional phylogenetic analyses were performed on putative single and low-copy genes extracted from off-target reads.Key ResultsFour major clades are supported within Zanthoxylum: the African clade, the Z. asiaticum clade, the Asian–Pacific–Australian clade and the American–eastern Asian clade. While overall support has improved, regions of conflict are similar to those previously observed. Gene tree discordances indicate a hybridization event in the ancestor of the Hawaiian lineage, and incomplete lineage sorting in the American backbone. Off-target putative single-copy genes largely confirm on-target results, and putative low-copy genes provide additional evidence for hybridization in the Hawaiian lineage. Only two of the five sections of Zanthoxylum are resolved as monophyletic.ConclusionsTarget enrichment is suitable for assessing phylogenetic relationships in Zanthoxylum. Our phylogenetic analyses reveal that current sectional classifications need revision. Quartet tree concordance indicates several instances of reticulate evolution. Off-target reads are proven useful to identify additional phylogenetically informative regions for bait refinement or gene tree based approaches.     

Congruent mammalian trees from mitochondrial and nuclear genes using Bayesian methods

Analyses of mitochondrial and nuclear gene sequences have often produced different mammalian tree topologies, undermining confidence in the merit of molecular approaches with respect to "traditional" morphological classification. The recent sequencing of the complete mitochondrial genomes of two additional rodents (Spalax judaei and Jaculus jaculus) and one lagomorph (Ochotona princeps) has prompted us to reinvestigate the issue. Using Bayesian phylogenetics, we found phylogenetic relationships between mammalian species highly congruent with previous results based on nuclear genes. Our results show the existence of four primary lineages of placental mammals: Xenarthra, Afrotheria, Laurasiatheria, and Euarchontoglires. Relationships between and within these lineages strongly suggest that the gene trees may also be congruent with the underlying species phylogeny.     

Interspecific divergence, intrachromosomal recombination, and phylogenetic utility of Y-chromosomal genes in Drosophila

Reconstruction of phylogenetic relationships among recently diverged species is complicated by three general problems: segregation of polymorphisms that pre-date species divergence, gene flow during and after speciation, and intra-locus recombination. In light of these difficulties, the Y chromosome offers several important advantages over other genomic regions as a source of phylogenetic information. These advantages include the absence of recombination, rapid coalescence, and reduced opportunity for interspecific introgression due to hybrid male sterility. In this report, we test the phylogenetic utility of Y-chromosomal sequences in two groups of closely related and partially inter-fertile Drosophila species. In the D. bipectinata species complex, Y-chromosomal loci unambiguously recover the phylogeny most consistent with previous multi-locus analysis and with reproductive relationships, and show no evidence of either post-speciation gene flow or persisting ancestral polymorphisms. In the D. simulans species complex, the situation is complicated by the duplication of at least one Y-linked gene region, followed by intrachromosomal recombination between the duplicate genes that scrambles their genealogy. We suggest that Y-chromosomal sequences are a useful tool for resolving phylogenetic relationships among recently diverged species, especially in male-heterogametic organisms that conform to Haldane's rule. However, duplication of Y-linked genes may not be uncommon, and special care should be taken to distinguish between orthologous and paralogous sequences.     

Further insights into the phylogeny of Arabidopsis (Brassicaceae) from nuclear Atmyb2 flanking sequence

Arabidopsis thaliana is the preeminent plant model organism. However, significant advances in evolution and ecology are being made by expanding the scope of research beyond this single species into the broader genus Arabidopsis. Surprisingly, few studies have rigorously investigated phylogenetic relationships between the nine Arabidopsis species, and this study evaluates both these and hypotheses related to two instances of intra-generic hybridization. DNA sequences from the 5' flanking region of the nuclear Atmyb2 gene from 12 of the 14 Arabidopsis taxa were used to reconstruct the generic phylogeny. The strict consensus tree was highly concordant with previous studies, identifying lineages corresponding to widespread species but exhibiting a large basal polytomy. Our data indicates that the paternal parent of the allopolyploid A. suecica is A. neglecta rather than A. arenosa s.l., although the need for a detailed phylogeographical study of these three species is noted. Finally, our data provided additional phylogenetic evidence of hybridization between Arabidopsis lyrata s.l. and A. halleri s.l. Taken together, the well-defined lineages within the genus and the potential for hybridization between them highlight Arabidopsis as a promising group for comparative and experimental studies of hybridization.     

Intrageneric phylogeny of Capsella (Brassicaceae) and the origin of the tetraploid C. bursa-pastoris based on chloroplast and nuclear DNA sequences

Polyploidization, often accompanied by hybridization, has been of major importance in flowering plant evolution. Here we investigate the importance of these processes for the evolution of the tetraploid crucifer Capsella bursa-pastoris using DNA sequences from two chloroplast loci as well as from three nuclear low-copy genes. The near-absence of variation at the C. bursa-pastoris chloroplast markers suggests a single and recent origin of the tetraploid. However, despite supporting a single phylogeny, chloroplast data indicate that neither of the extant Capsella diploids is the maternal parent of the tetraploid. Combined with data from the three nuclear loci, our results do not lend support to previous hypotheses on the origin of C. bursa-pastoris as an allopolyploid between the diploids C. grandiflora and C. rubella or an autopolyploid of C. grandiflora. Nevertheless, for each locus, some of the C. bursa-pastoris accessions harbored C. rubella alleles, indicating that C. rubella contributed to the gene pool of C. bursa-pastoris, either through allopolyploid speciation or, more likely, through hybridization and introgression. To our knowledge, this study is the first of a wild, nonmodel plant genus that uses a combination of chloroplast and multiple low-copy nuclear loci for phylogenetic inference of polyploid evolution.     

Incongruence of mitochondrial and nuclear gene trees in the Carabid beetles Ohomopterus

We studied the molecular phylogeny of the carabid subgenus Ohomopterus (genus Carabus), using two mitochondrial (mt) DNA regions (16SrRNA and NADH dehydrogenase subunit 5) and three nuclear DNA regions (wingless, phosphoenolpyruvate carboxykinase, and an anonymous locus). We revisited the previously reported incongruence between the distribution of mtDNA markers and morphologically defined species (Su et al., 1996; J. Mol. Evol. 43:662-671), which those authors attributed to "type switching", a concerted change in many morphological characters that results in the repeated evolution of a particular morphological type. Our mtDNA gene tree obtained from 44 individuals representing all 15 currently recognized species of Ohomopterus revealed that haplotypes isolated from individuals of a single "species" were frequently separated into distant clades, confirming the previous report. The three nuclear markers generally conformed better-with the morphologically defined species than did the mitochondrial markers. The phylogenetic signal in mtDNA and nuclear DNA data differed strongly, and these two partitions were significantly incongruent with each other according to the incongruence length difference test of Farris et al. (1994; Cladistics 10:315-320), although the three nuclear partitions were not homogeneous either. Our results did not support the type-switching hypothesis that had been proposed to fit the morphological data to the mitochondrial gene tree: The incongruence of the mtDNA tree with other nuclear markers indicates that the mtDNA-based tree does not reflect species history any better than the morphological data do. Incongruence of gene trees in Ohomopterus may have been promoted by the complex processes of geographic isolation and hybridization in the Japanese Archipelago that have led to occasional gene flow and recombination between separated entities. The occurrence of reticulate patterns in this group is intriguing, because species of Ohomopterus exhibit extremely divergent genitalic structures that represent a highly efficient reproductive isolation mechanism.     

A simple algorithm to infer gene duplication and speciation events on a gene tree

MOTIVATION: When analyzing protein sequences using sequence similarity searches, orthologous sequences (that diverged by speciation) are more reliable predictors of a new protein's function than paralogous sequences (that diverged by gene duplication), because duplication enables functional diversification. The utility of phylogenetic information in high-throughput genome annotation ('phylogenomics') is widely recognized, but existing approaches are either manual or indirect (e.g. not based on phylogenetic trees). Our goal is to automate phylogenomics using explicit phylogenetic inference. A necessary component is an algorithm to infer speciation and duplication events in a given gene tree. RESULTS: We give an algorithm to infer speciation and duplication events on a gene tree by comparison to a trusted species tree. This algorithm has a worst-case running time of O(n(2)) which is inferior to two previous algorithms that are approximately O(n) for a gene tree of sequences. However, our algorithm is extremely simple, and its asymptotic worst case behavior is only realized on pathological data sets. We show empirically, using 1750 gene trees constructed from the Pfam protein family database, that it appears to be a practical (and often superior) algorithm for analyzing real gene trees. AVAILABILITY: http://www.genetics.wustl.edu/eddy/forester.     

Species trees for the tree swallows (Genus Tachycineta): An alternative phylogenetic hypothesis to the mitochondrial gene tree

The New World swallow genus Tachycineta comprises nine species that collectively have a wide geographic distribution and remarkable variation both within- and among-species in ecologically important traits. Existing phylogenetic hypotheses for Tachycineta are based on mitochondrial DNA sequences, thus they provide estimates of a single gene tree. In this study we sequenced multiple individuals from each species at 16 nuclear intron loci. We used gene concatenated approaches (Bayesian and maximum likelihood) as well as coalescent-based species tree inference to reconstruct phylogenetic relationships of the genus. We examined the concordance and conflict between the nuclear and mitochondrial trees and between concatenated and coalescent-based inferences. Our results provide an alternative phylogenetic hypothesis to the existing mitochondrial DNA estimate of phylogeny. This new hypothesis provides a more accurate framework in which to explore trait evolution and examine the evolution of the mitochondrial genome in this group.     

Phylogeny of PgiC gene in Shorea and its closely related genera (Dipterocarpaceae), the dominant trees in Southeast Asian tropical rain forests

Dipterocarpaceae, trees that dominate tropical rain forests in Southeast Asia consist of many economically and ecologically important species. We determined partial sequences of the PgiC gene from species of Shorea, Hopea, Neobalanocarpus, and Parashorea to elucidate phylogenetic relationships among the species of these genera, which have been regarded as interrelated. The sequences generated a gene tree with better resolution than previous cpDNA trees. The PgiC tree is essentially consistent with cpDNA trees, except for the placement of Neobalanocarpus. The PgiC tree shows that Neobalanocarpus is nested within White Meranti of Shorea, whereas this genus forms a clade with Hopea in cpDNA trees. This conflict suggests that Neobalanocarpus is derived via hybridization between White Meranti of Shorea and Hopea. Species belonging to each of three timber groups (Yellow Meranti, Balau, and Red Meranti) within Shorea are monophyletic. Together they form a monophyletic clade distinct from White Meranti. Botanical sections within Red Meranti appear not to be monophyletic. An extensive number of shared polymorphisms among species and consequential lack of monophyly of intraspecific haplotypes are found in Red Meranti. Potential causes of this phenomenon, including persistence of ancestral polymorphisms and gene flow via interspecific hybridization, are discussed.