Identification of putative orthologous genes for the phylogenetic reconstruction of temperate woody bamboos (Poaceae: Bambusoideae)

The temperate woody bamboos (Arundinarieae) are highly diverse in morphology but lack a substantial amount of genetic variation. The taxonomy of this lineage is intractable, and the relationships within the tribe have not been well resolved. Recent studies indicated that this tribe could have a complex evolutionary history. Although phylogenetic studies of the tribe have been carried out, most of these phylogenetic reconstructions were based on plastid data, which provide lower phylogenetic resolution compared with nuclear data. In this study, we intended to identify a set of desirable nuclear genes for resolving the phylogeny of the temperate woody bamboos. Using two different methodologies, we identified 209 and 916 genes, respectively, as putative single copy orthologous genes. A total of 112 genes was successfully amplified and sequenced by next‐generation sequencing technologies in five species sampled from the tribe. As most of the genes exhibited intra‐individual allele heterozygotes, we investigated phylogenetic utility by reconstructing the phylogeny based on individual genes. Discordance among gene trees was observed and, to resolve the conflict, we performed a range of analyses using BUCKy and HybTree. While caution should be taken when inferring a phylogeny from multiple conflicting genes, our analysis indicated that 74 of the 112 investigated genes are potential markers for resolving the phylogeny of the temperate woody bamboos.     

When is enough,enough in phylogenetics? A case in point from Hordeum (Poaceae)

Direct optimization was used to reconstruct the phylogeny of the 26 diploid taxa included in the genus Hordeum. The total data set was composed of 16 nucleotide sequence regions from the nuclear as well as the plastid genome. The nine nuclear regions were from single‐copy, protein coding genes located on six of the seven chromosome pairs in the diploid H. vulgare genome. The seven plastid regions comprise protein coding genes as well as intergenic regions. Studies of character congruence between data partitions showed no correlation between chromosomal location and congruence among the nuclear sequences and a level of congruence among the plastid sequences comparable with the level among the nuclear sequences. Combined analysis of all data resolved the phylogeny completely with most clades being robust and well supported. However, due to incongruence among data partitions some relationships are still and likely to remain ambiguously inferred. Rather than adding still more genes to the phylogenetic analyses, patterns of incongruence may be better explored by adding data from multiple specimens per taxon. For some species relationships the plastid data appear positively misleading, emphasizing the need for caution if plastid data are the only or dominant type of data used for phylogenetic reconstruction and subsequent re‐classification.
© The Willi Hennig Society 2011.     

Scanning electron microscopy of the leaf epidermis in Arundinarieae (Poaceae: Bambusoideae): evolutionary implications of selected micromorphological features

Arundinarieae (temperate woody bamboos) is a taxonomically challenging group and the taxa in this tribe show a low level of DNA sequence variation, especially for plastid markers. As a result of the scarcity of flowering collections, species identification mainly depends on vegetative features, and leaf epidermal micromorphology has proven to be useful in bamboo taxonomy. In this study, we used scanning electron microscopy to investigate the abaxial leaf epidermal micromorphology of 94 species in 24 genera of Arundinarieae and three species of Bambuseae. Leaf epidermal characteristics differed mainly in the number and distribution pattern of papillae. Seven papilla forms were assigned, which had little taxonomic value at the generic and subtribal levels. However, the papilla patterns combined with other leaf epidermal features were useful in the discrimination of several species. Most papilla types were shared between tribes Arundinarieae and Bambuseae. We assessed the phylogenetic implications of the leaf epidermal micromorphology based on recently published molecular phylogenetic analyses. With one exception, none of the papilla types corresponded to a particular clade in the plastid or nuclear phylogenetic trees. Based on these results, the papilla patterns have limited phylogenetic value. © 2014 The Linnean Society of London, Botanical Journal of the Linnean Society, 2014, 176 , 46–65.     

Plastid and nuclear DNA polymorphism reveals historical processes of isolation and reticulation in the olive tree complex (Olea europaea)

Aim The olive tree is considered one of the best indicators of the Mediterranean climate. The species’ distribution is associated with geographical and bioclimatic factors, as well as being influenced by a long period of cultivation. Despite concerted efforts of different research groups, the origin of the Mediterranean olive tree still remains elusive. In the present study, relationships between taxa and populations covering the entire range of Olea europaea were investigated using both maternal (plastid genome) and biparental (nuclear genome) markers to disclose evolutionary patterns in the olive complex. Phylogenetic and phylogeographical results of the two‐genome analyses were interpreted in a biogeographical context. Location Mediterranean, temperate and subtropical floristic regions of the Old World. Methods Phylogeographical reconstructions of plastid DNA polymorphism were performed using microsatellites, restriction sites and indels on a wide sample of 185 representative trees across the Old World, including 28 herbarium specimens from remote areas. Additionally, the potential utility of one ITS‐1 pseudogene for phylogenetic analyses was explored using Bayesian and maximum parsimony approaches on a subsample of 38 olive trees. Results Forty plastid haplotypes were recognized and split into two lineages and seven sublineages. The analysis of ITS‐1 sequences also allowed the identification of seven well differentiated groups. Distribution of plastid and ribosomal DNA lineages was congruent, but particular cases of phylogenetic incongruence were disclosed (particularly in the Sahara and Madeira). Lastly, two divergent ITS‐1 copies were isolated from the same sample of four individuals of different subspecies. Main conclusions Phylogenetic congruence of both ITS‐1 and plastid lineages suggested an evolutionary scenario of predominant isolation during the Plio‐Pleistocene in Macaronesia, the Mediterranean, southern Africa, eastern Africa and Asia. The Saharan desert appeared to have played an important role of vicariant barrier between southern and northern African populations in early times. Incongruence of some plastid and nuclear results, as well as intermingled ITS‐1 copies of different lineages in single individuals, was interpreted as a result of recurrent reticulation events in the olive complex. We identified an ancient hybrid zone from the Sahara to north‐eastern African mountains, where divergent plastid and nuclear lineages still co‐exist. Results of this paper, and previous studies, suggest that the cultivated olive originated from a pre‐Quaternary Mediterranean ancestor, with no evidence for a recent hybrid origin. In contrast, a continuous process of olive domestication through local hybridization events of cultivated trees with natural populations may have brought about a remarkably high genomic diversity among cultivated trees across the Mediterranean.     

Phylogeny of Nolana (Solanaceae) of the Atacama and Peruvian deserts inferred from sequences of four plastid markers and the nuclear LEAFY second intron

The phylogeny of Nolana (Solanaceae), a genus primarily distributed in the coastal Atacama and Peruvian deserts with a few species in the Andes and one species endemic to the Galápagos Islands, was reconstructed using sequences of four plastid regions (ndhF, psbA-trnH, rps16-trnK and trnC-psbM) and the nuclear LEAFY second intron. The monophyly of Nolana was strongly supported by all molecular data. The LEAFY data suggested that the Chilean species, including Nolana sessiliflora, the N. acuminata group and at least some members of the Alona group, are basally diverged, supporting the Chilean origin of the genus. Three well-supported clades in the LEAFY tree were corroborated by the SINE (short interspersed elements) or SINE-like insertions. Taxa from Peru are grouped roughly into two clades. Nolana galapagensis from the Galápagos Island is most likely to have derived from a Peruvian ancestor. The monophyly of the morphologically well-diagnosed Nolana acuminata group (N. acuminata, N. baccata, N. paradoxa, N. parviflora, N. pterocarpa, N. rupicola and N. elegans) was supported by both plastid and LEAFY data. Incongruence between the plastid and the LEAFY data was detected concerning primarily the positions of N. sessiliflora, N. galapagensis, taxa of the Alona group and the two Peruvian clades. Such incongruence may be due to reticulate evolution or in some cases lineage sorting of plastid DNA. Incongruence between our previous GBSSI trees and the plastid-LEAFY trees was also detected concerning two well-supported major clades in the GBSSI tree. Duplication of the GBSSI gene may have contributed to this incongruence.     

Viburnum phylogeny: evidence from the duplicated nuclear gene GBSSI

DNA sequencing studies of the granule-bound starch synthase gene (GBSSI) indicate the presence of two loci in Viburnum. Gene trees from separate and combined phylogenetic analyses of the GBSSI paralogues are generally congruent with each other and with trees from previous analyses, especially those of Donoghue et al. [Syst. Bot. 29 (2004) 188] based on nuclear ribosomal ITS and chloroplast trnK intron DNA sequences. Specifically, our GBSSI trees confirm (i) the monophyly of some and non-monophyly of other traditionally recognized taxonomic sections, (ii) the presence of three major supra-sectional lineages within Viburnum, and (iii) the resolution of many species relationships within the section-level clades. Analyses of GBSSI also provide greater resolution of relationships within the largest supra-sectional lineage. Relationships at the base of the Viburnum phylogeny remain uncertain; in particular, the position of the root, relationships among the supra-sectional clades, and the exact placement of several smaller groups (e.g., Viburnum clemensiae, Viburnum urceolatum, and section Pseudotinus). In two lineages each GBSSI paralogue is represented by two distinct sequences. The presence of additional copies appears to be correlated with polyploidy in these clades. Placement of the homoeologues in our gene trees suggests the possibility of a hybrid origin for these polyploids.     

Tracking ancient polyploids: a retroposon insertion reveals an extinct diploid ancestor in the polyploid origin of belladonna

Polyploidy is a prominent process in plant evolution and adaptation, but molecular phylogenetic studies of polyploids based on DNA sequences have often been confounded by their complex gene and genome histories. We report here a retroposon insertion in the nuclear gene granule-bound starch synthase I (GBSSI or "waxy") that clearly reveals the ancient hybrid history of the medically important polyploid species belladonna (Atropa belladonna) and resolves the controversy over the taxonomic group to which it belongs, the tribe Hyoscyameae (Solanaceae). Our inferences based on the pattern of presence or absence of the retroposon insertion are corroborated by phylogenetic analyses of the GBSSI gene sequences. This case may suggest that retroposons are promising molecular markers to study polyploid evolution.     

Subtribal and generic relationships of Maxillarieae (Orchidaceae) with emphasis on Stanhopeinae: combined molecular evidence

The monophyly of and phylogenetic relationships within the orchid tribe Maxillarieae Pfitzer were evaluated using parsimony analyses of combined nuclear ribosomal and plastid DNA sequence data of ITS 1 and 2, matK, and the trnL intron and the trnL-F intergene spacer. Each of the separate analyses produced highly congruent but weakly supported patterns (by the bootstrap), so these were combined in a single analysis. Analysis of 90 ingroup taxa (representing ～35% of currently recognized genera) and four outgroup taxa produced resolved and highly supported cladograms. Based on the cladograms, we recognize six subtribes: Eriopsidinae, Oncidiinae (including Pachyphyllinae, Ornithocephalinae, and Telipogoninae), Stanhopeinae, Coeliopsidinae, Maxillariinae (including Lycastinae and Bifrenariinae), and Zygopetalinae (including Cryptarrheninae, Dichaeinae, Huntleyinae, and Warreinae). Stanhopeinae were sampled most intensively; their generic relationships were highly resolved in the analysis and largely agree with currently accepted generic concepts based on morphology. Coeliopsidinae (Coeliopsis, Lycomormium, Peristeria) are sister to Stanhopeinae. Correlations are drawn among phylogeny, pollination mechanisms, and life history traits.     

Allopolyploidy in the Wintergreen Group of tribe Gaultherieae (Ericaceae) inferred from low‐copy nuclear genes

DNA sequence data from the low‐copy nuclear genes waxy (GBSSI) and leafy were compared with plastid and ITS sequence data from prior studies to reconstruct phylogenetic relationships in the Wintergreen Group of tribe Gaultherieae (Ericaceae). We conducted phylogenetic analysis with 102 species that includes representatives of all 15 major clades previously identified within the Wintergreen Group and that together span its circum‐Pacific distribution. Results yielded two distinct homeologous copies of waxy for two of the clades, each in widely separated parts of the tree. It also yielded two copies of leafy for one of the clades; only one copy of leafy was found for the other clade, but it was placed in the same major clade as its waxy counterpart and well away from its placement in a prior plastid analysis. A combined four‐locus (waxy, leafy, ITS and plastid data) phylogenetic analysis of all available relevant data placed the copies of each of the clades in two distinct positions in the phylogeny with strong overall statistical support. In combination with evidence from morphology, reproductive biology and cytology, the results suggest that these clades arose through allopolyploid hybridization between lineages deep in the phylogeny but relatively close geographically. This finding confirms previous assumptions that hybridization has played an important role in the evolution of the Gaultherieae.     

Phylogenomics of Fargesia and Yushania reveals a history of reticulate evolution

Reticulate evolution is a common and important driving force in angiosperm evolution. In this study, we analyzed the phylogenetic signals of genomic regions with different inheritance patterns to understand the evolutionary process of organisms using species-rich Himalaya–Hengduan taxa of bamboos (Fargesia Franchet and Yushania Keng). We constructed phylogenetic trees using different sampling strategies and reconstruction methods based on genome skimming and double digest restriction-site-associated DNA sequencing data. We assessed the congruence of topologies generated from different datasets and employed several approaches to reveal the causes of phylogenetic incongruence, including the detection of hybridization and introgression using PhyloNetworks and the D-statistic test (ABBA-BABA test). We found that, in the plastome-based phylogeny, Fargesia bamboos can be clustered into three groups and Yushania was nested within one of them, which contradicts the nuclear–double digest restriction-site-associated DNA sequencing-based phylogeny. Moreover, the genetic variation of chloroplast DNA is significantly correlated with geographical distribution. The strong signal of incomplete lineage sorting, hybridization, introgression, and cytoplasmic gene flow found among genera and species suggests that reticulate evolution is the main cause for the phylogenetic incongruence between nuclear and chloroplast datasets. Our results add evidence that genomes with different inheritance patterns can reveal distinct evolutionary histories of species and suggest that reticulate evolution is prevalent in rapidly diversifying groups.     

Evolutionary processes in the genus Coreocarpus: insights from molecular phylogenetics

A molecular phylogenetic study of the plant genus Coreocarpus was conducted using nuclear (ITS) and plastid (rpl16 intron) DNA sequences, with phylogenies of the nuclear and plastid sequences highly congruent in defining a monophyletic group of six species (core Coreocarpus), although three other species often placed within the genus were excluded. Relationships within the genus are largely but not totally concordant with prior biosystematic studies. Despite strong molecular support, no morphological characters uniting the six species of core Coreocarpus have been identified; retention of plesiomorphic characters and the genetic lability of characters are two probable factors contributing to lack of consistent defining characters. The age of the core Coreocarpus is estimated at 1 million years because the basal species is endemic to a volcanic island that emerged in the past million years. Mapping the results of earlier breeding studies on the molecular phylogeny showed that use of cross-compatibility as a criterion for species delimitation would result in the recognition of paraphyletic species. Prior field, morphological, and biosystematic studies provided no indication of past hybridization in the evolution of Coreocarpus, and species in the genus appeared to be well defined morphologically. However, three instances of incongruence were observed. Two of these were between the nuclear and plastid partitions, and the third was between the morphological species assignment of one accession and the molecular data. If hybridization accounts for incongruence between the nuclear and plastid data, it occurred between species that now appear to be cross-incompatible and allopatric. The incongruence between morphological species assignment and the molecular data could be the result of parallel fixation of characters that have a simple genetic basis. This study suggests that the evolutionary history of Coreocarpus is much more complex than indicated from prior biosystematic investigations and that biosystematic and molecular phylogenetic studies may complement each other for elucidating the evolution and phylogeny of a group.     

Phylogeny of the tribe Aveneae (Pooideae, Poaceae) inferred from plastid trnT-F and nuclear ITS sequences

New insights into evolutionary trends in the economically important oat tribe (Aveneae) are presented. Plastid trnT-F and nuclear ribosomal ITS sequences were used to reconstruct the phylogeny of the Aveneae-Poeae-Seslerieae complex (Pooideae, Poaceae) through Bayesian- and maximum parsimony-based analyses, separately and in combination. The plastid data identified a strongly supported core Aveneae lineage that separated from other former Aveneae and Poeae groups. Koeleriinae, Aveninae, and Agrostidinae emerged as the main groups of this core Aveneae, which also included other minor subgroups with uncertain relationships and a few former Poeae members. Several former Aveneae representatives were also placed in independent sublineages in Poeae. Seslerieae resolved as close allies of Poeae or Aveneae in the plastid and nuclear topologies, respectively. Because of the intermingling of some Aveneae and Seslerieae lineages in Poeae and vice versa, we propose to expand Poeae to include all the aforementioned lineages. This best reflects our current understanding of the phylogeny of these important temperate grasses and sheds light on their evolutionary history.     

Hybrid genera in Liatrinae (Asteraceae: Eupatorieae)

Liatrinae is a small subtribe of Eupatorieae that occurs in North America with a center of generic-level diversity in the southeastern United States. Molecular phylogenetic data were sought to assess whether two monotypic genera, Garberia and Hartwrightia, are accurately placed in the subtribe, and to resolve questions of the generic-level classification of Carphephorus. Phylogenetic analyses of nuclear ITS/ETS and plastid DNA data indicated that Garberia is the basalmost diverging lineage, and that Hartwrightia is phylogenetically embedded in the subtribe. There was significant incongruence between the ITS/ETS and plastid DNA datasets in the placement of Hartwrightia and another monotypic genus, Litrisa, suggesting that both are of original hybrid origin. The results also showed that Carphephorus s.l. is not monophyletic, and even after removal of the two species of Trilisa, it is still paraphyletic to Liatris. The apparent hybrid origin of Hartwrightia, which is morphologically transgressive relative to its inferred parental lineages, suggests that reticulation between phylogenetically distinct lineages may be a recurrent problem for phylogenetic estimation in Asteraceae.     

Molecular phylogeny and historical biogeography of Caltha (Ranunculaceae) based on analyses of multiple nuclear and plastid sequences

Caltha is a widely distributed genus in the buttercup family (Ranunculaceae) showing interesting distribution patterns in both hemispheres. Evolutionary history ofCaltha was examined by means of phylogenetic, molecular dating, and historical biogeographic analyses with a more comprehensive sampling than previous studies. The internal transcribed spacer from the nuclear genome and trnL-F and atpB-rbcL regions from the plastid genome were used and analyzed using parsimony and Bayesian methods. Divergence time was estimated using Bayesian dating analyses with multiple fossil calibrations. Historical biogeography was inferred using the Bayes-DIVA method implemented in RASP. We obtained a well-resolved and well-supported phylogeny within the Caltha lineage. Caltha natansPall. diverged first from the genus and the other species grouped into two clades. Our expanded sampling scheme revealed a complicated evolutionary pattern in theC. palustris complex. Caltha sinogracilis W. T. Wang was resolved to be a member of the C. palustris complex, rather than closely related to C. scaposa Hook. f. & Thomson. Caltha rubriflora B. L. Burtt & Lauener was also revealed to be not just a red-flower form of C. sinogracilis. The diversification of the genus began at 50.5 mya (95% high posterior density: 37.1–63.9 mya), and its ancestral range was very probably in the Northern Hemisphere. The South American species may derive from western North American ancestors that dispersed along the western American Cordillera during the Cenozoic era. The vicariance model of the Southern Hemisphere species proposed by a previous study was rejected in this study.     

Phylogenetic utility of nuclear nitrate reductase: A multi-locus comparison of nuclear and chloroplast sequence data for inference of relationships among American Lycieae (Solanaceae)

There has been considerable interest and research into the benefits of multiple low-copy nuclear regions for phylogenetic studies at low taxonomic levels. In this study, the phylogenetic utility of DNA sequence data from two low-copy nuclear genes, nitrate reductase (NIA) and granule-bound starch synthase I (GBSSI), was compared with data from nuclear ribosomal ITS and the cpDNA spacers trnT–trnF and trnD–trnT across 33 closely related taxa in tribe Lycieae (Solanaceae). The NIA data (introns 1–2) had the most parsimony-informative characters, with over twice the number provided by GBSSI, but NIA also had greater homoplasy. Although gene trees were generally concordant across the four regions, there were some notable areas of incongruence, suggesting both incomplete lineage sorting as well as possible reticulate origins.     

Genome-scale data, angiosperm relationships, and "ending incongruence": a cautionary tale in phylogenetics

As systematists grapple with assembling the Tree of Life, recent studies have encouraged a genomic-scale approach, obtaining DNA sequence data for entire nuclear, plastid or mitochondrial genomes for a few exemplar taxa. Some have proclaimed that this comparative genomic strategy heralds the end of incongruence in phylogeny reconstruction. Although we applaud the use of many genes to resolve phylogenetic patterns, there is a significant caveat. In spite of, or even because of, the abundant data per taxon, whole-genome sequencing for a few exemplars can provide completely resolved and strongly supported, but incorrect, evolutionary reconstructions. We provide a conspicuous example that includes Amborella, the putative sister of all other extant angiosperms, highlighting the limits of phylogenetics when whole genomes are used but taxon sampling is poor.     

Chloroplast phylogenomics of liverworts: a reappraisal of the backbone phylogeny of liverworts with emphasis on Ptilidiales

As one of the four main lineages diverging from the early diversification of land plants, the phylogeny of liverworts holds the information about nearly 500 Myr of independent adaptation to changing environments. Thus, resolving the phylogenetic history of liverworts will provide unique insights into the successful diversification of early land plants in terrestrial ecosystems. However, the deep diverging events of this group remain incompletely resolved, such as the definite position of Ptilidiales. Here, we aimed to reconstruct the backbone relationships of liverworts using 84 protein-coding chloroplast genes, a dataset comprising 35 representatives from all major lineages of liverworts, and three phylogenetic analyses, namely maximum parsimony, maximum likelihood and Bayesian inference. To test the impact of composition biases, the phylogenetic analyses were carried out using three alignments representing the same dataset either as: (i) nucleotides, (ii) amino acids, or (iii) recoded nucleotides applying ambiguity base code. Chloroplast genome data consistently supported the monophyletic origin of three major lineages in liverworts, as well as the majority of backbone relationships. Ptilidiales were found to be sister to Jungermanniales. The rapid accumulation of G/C tracks as a consequence of increased GC content is an important cause for the long branches inferred in this group. Our study not only provides empirical evidence to support the significance of plastid genome sequencing to reconstruct the phylogeny of this important plant lineage, but also suggests that the GC content has played a critical role in the evolutionary dynamics of plastid genomes in land plants.     

Genome skimming by shotgun sequencing helps resolve the phylogeny of a pantropical tree family

Whole genome sequencing is helping generate robust phylogenetic hypotheses for a range of taxonomic groups that were previously recalcitrant to classical molecular phylogenetic approaches. As a case study, we performed a shallow shotgun sequencing of eight species in the tropical tree family Chrysobalanaceae to retrieve large fragments of high‐copy number DNA regions and test the potential of these regions for phylogeny reconstruction. We were able to assemble the nuclear ribosomal cluster (nrDNA), the complete plastid genome (ptDNA) and a large fraction of the mitochondrial genome (mtDNA) with approximately 1000×, 450× and 120× sequencing depth respectively. The phylogenetic tree obtained with ptDNA resolved five of the seven internal nodes. In contrast, the tree obtained with mtDNA and nrDNA data were largely unresolved. This study demonstrates that genome skimming is a cost‐effective approach and shows potential in plant molecular systematics within Chrysobalanaceae and other under‐studied groups.     

Using nuclear genes to reconstruct angiosperm phylogeny at the species level: A casestudy with Brassicaceae species

Angiosperm phylogeny has been investigated extensively using organellar sequences; recent efforts using nuclear genes have also been successful in reconstructing angiosperm phylogenies at family or deeper levels. However, it is not clear whether nuclear genes are also effective in understanding relationships between species in a genus. Here we present a case study of phylogeny at generic and specific levels with nuclear genes, using Brassicaceae taxa as examples. Brassicaceae includes various crops and the model plant Arabidopsis thaliana. A recent study showed that nuclear genes can provide well-resolved relationships between tribes and larger lineages in Brassicaceae, but few species were included in any given genus. We present a phylogeny with multiple species in each of five genera within Brassicaceae for a total of 65 taxa, using three protein-coding nuclear genes, MLH1, SMC2, and MCM5, with up to approximately 10 200 base pairs (in both exons and introns). Maximum likelihood and Bayesian analyses of the separate gene regions and combined data reveal high resolution at various phylogenetic depths. The relationships between genera here were largely congruent with previous results, with further resolution at the species level. Also, we report for the first time the affinity of Cardamine rockii with tribe Camelineae instead of other Cardamine members. In addition, we report sequence divergence at three levels: across angiosperms, among Brassicaceae species, and between Arabidopsis ecotypes. Our results provide a robust species-level phylogeny for a number of Brassicaceae members and support an optimistic perspective on the phylogenetic utility of conserved nuclear data for relatively recent clades.