An expanded plastid DNA phylogeny of Orchidaceae and analysis of jackknife branch support strategy

Phylogenomic insights into the historical biogeography, character-state evolution, and species diversification rates of Cypripedioideae (Orchidaceae)

Cypripedioideae (slipper orchids; Orchidaceae) currently consist of ∼200 herbaceous species with a strikingly disjunctive distribution in tropical and temperate regions of both hemispheres. In this study, an updated phylogeny with representatives from all five cypripedioid genera was presented based on maximum likelihood and Bayesian inference of plastome and low-copy nuclear genes. Phylogenomic analyses indicated that each genus is monophyletic, but some relationships (e.g., those among Cypripedium sects. Acaulia, Arietinum, Bifolia, Flabellinervia, Obtusipetala and Palangshanensia) conflict with those in previous studies based on Sanger data. Cypripedioideae appeared to have arisen in South America and/or the adjacent Qinghai-Tibet Plateau and Hengduan Mountains ∼35 Mya. We inferred multiple dispersal events between East Asia and North America in Cypripedium, and between mainland Southeast Asia and the Malay Archipelago in Paphiopedilum. In the Americas, divergences among four genera (except Cypripedium) occurred around 31-20 Mya, long before the closure of the Isthmus of Panama, indicating the importance of long-distance dispersal. Evolutionary patterns between morphological and plastome character evolution suggested several traits, genome size and NDH genes, which are likely to have contributed to the success of slipper orchids in alpine floras and low-elevation forests. Species diversification rates were notably higher in epiphytic clades of Paphiopedilum than in other, terrestrial cypripedioids, paralleling similar accelerations associated with epiphytism in other groups. This study also suggested that sea-level fluctuations and mountain-building processes promoted the diversification of the largest genera, Paphiopedilum and Cypripedium.

Genus-level revision of the Heterozerconoidea (Parasitiformes: Mesostigmata)

The genera of Heterozerconoidea are revised based on a species-level analysis of relationships in the group. The family Discozerconidae in its current state may be paraphyletic. Diagnoses for the genera are updated, and a catalog of all described species is provided. As part of this re-analysis two new genera, Amyzozercon and Ecuazercon, and four new species are described, and a key to the genera is provided. Possible evolutionary implications of the proposed set of relationships in terms of biogeography and the evolution of podospermy are discussed.

Lineage and role in integrative taxonomy of a heterotrophic orchid complex

Lineage-based species definitions applying coalescent approaches to species delimitation have become increasingly popular. Yet, the application of these methods and the recognition of lineage-only definitions have recently been questioned. Species delimitation criteria that explicitly consider both lineages and evidence for ecological role shifts provide an opportunity to incorporate ecologically meaningful data from multiple sources in studies of species boundaries. Here, such criteria were applied to a problematic group of mycoheterotrophic orchids, the Corallorhiza striata complex, analysing genomic, morphological, phenological, reproductive-mode, niche, and fungal host data. A recently developed method for generating genomic polymorphism data-ISSRseq-demonstrates evidence for four distinct lineages, including a previously unidentified lineage in the Coast Ranges and Cascades of California and Oregon, USA. There is divergence in morphology, phenology, reproductive mode, and fungal associates among the four lineages. Integrative analyses, conducted in population assignment and redundancy analysis frameworks, provide evidence of distinct genomic lineages and a similar pattern of divergence in the extended data, albeit with weaker signal. However, none of the extended data sets fully satisfy the condition of a significant role shift, which requires evidence of fixed differences. The four lineages identified in the current study are recognized at the level of variety, short of comprising different species. This study represents the most comprehensive application of lineage + role to date and illustrates the advantages of such an approach.

Evolutionary Relationships and Range Evolution of Greenhood Orchids (Subtribe Pterostylidinae): Insights From Plastid Phylogenomics

Australia harbours a rich and highly endemic orchid flora with over 90% of native species found nowhere else. However, little is known about the assembly and evolution of Australia's orchid flora. Here, we used a phylogenomic approach to infer evolutionary relationships, divergence times and range evolution in Pterostylidinae (Orchidoideae), the second largest subtribe in the Australian orchid flora, comprising the genera Pterostylis and Achlydosa. Phylogenetic analysis of 75 plastid genes provided well-resolved and supported phylogenies. Intrageneric relationships in Pterostylis were clarified and monophyly of eight of 10 sections supported. Achlydosa was found to not form part of Pterostylidinae and instead merits recognition at subtribal level, as Achlydosinae. Pterostylidinae were inferred to have originated in eastern Australia in the early Oligocene, coinciding with the complete separation of Australia from Antarctica and the onset of the Antarctic Circumpolar Current, which led to profound changes in the world's climate. Divergence of all major lineages occurred during the Miocene, accompanied by increased aridification and seasonality of the Australian continent, resulting in strong vegetational changes from rainforest to more open sclerophyllous vegetation. The majority of extant species were inferred to have originated in the Quaternary, from the Pleistocene onwards. The rapid climatic oscillations during the Pleistocene may have acted as important driver of speciation in Pterostylidinae. The subtribe underwent lineage diversification mainly within its ancestral range, in eastern Australia. Long-distance dispersals to southwest Australia commenced from the late Miocene onwards, after the establishment of the Nullarbor Plain, which constitutes a strong edaphic barrier to mesic plants. Range expansions from the mesic into the arid zone of eastern Australia (Eremaean region) commenced from the early Pleistocene onwards. Extant distributions of Pterostylidinae in other Australasian regions, such as New Zealand and New Caledonia, are of more recent origin, resulting from long-distance dispersals from the Pliocene onwards. Temperate eastern Australia was identified as key source area for dispersals to other Australasian regions.

Target sequence capture in orchids: Developing a kit to sequence hundreds of single-copy loci

PREMISE

Understanding relationships among orchid species and populations is of critical importance for orchid conservation. Target sequence capture has become a standard method for extracting hundreds of orthologous loci for phylogenomics. Up-front cost and time associated with design of bait sets makes this method prohibitively expensive for many researchers. Therefore, we designed a target capture kit to reliably sequence hundreds of orthologous loci across orchid lineages.

METHODS

We designed an Orchidaceae target capture bait set for 963 single-copy genes identified in published orchid genome sequences. The bait set was tested on 28 orchid species, with representatives of the subfamilies Cypripedioideae, Orchidoideae, and Epidendroideae.

RESULTS

Between 1,518,041 and 87,946,590 paired-end 150-base reads were generated for target-enriched genomic libraries. We assembled an average of 812 genes per library for Epidendroideae species and a mean of 501 genes for species in the subfamilies Orchidoideae and Cypripedioideae. Furthermore, libraries had on average 107 of the 254 genes that are included in the Angiosperms353 bait set, allowing for direct comparison of studies using either bait set.

DISCUSSION

The Orchidaceae963 kit will enable greater accessibility and utility of next-generation sequencing for orchid systematics, population genetics, and identification in the illegal orchid trade.

Molecular systematics and the evolution of mycoheterotrophy of tribe Neottieae (Orchidaceae, Epidendroideae)

Neottieae comprise about 150-200 species and are distributed mainly in temperate and subtropical zones of the northern hemisphere. Mycoheterotrophy is common in Neottieae. Based on three DNA markers and a broad sampling of Neottieae, these results indicate that Neottieae is strongly supported as monophyletic and Palmorchis is sister to the remaining genera of Neottieae. Holopogon and Neottia s.s. are deeply nested within Listera. The habit of leafless mycotrophy has independently evolved at least three times in Neottieae, one in Cephalanthera, another in Neottia s.l. and the third in the clade formed by Limodorum and Aphyllorchis.

Development and evolution of extreme synorganization in angiosperm flowers and diversity: a comparison of Apocynaceae and Orchidaceae

BACKGROUND AND AIMS

Apocynaceae and Orchidaceae are two angiosperm families with extreme flower synorganization. They are unrelated, the former in eudicots, the latter in monocots, but they converge in the formation of pollinia and pollinaria, which do not occur in any other angiosperm family, and for which extreme synorganization of floral organs is a precondition. In each family extensive studies on flower development and evolution have been performed; however, newer comparative studies focusing on flower synorganization and involving both families together are lacking.

SCOPE

For this study an extensive search through the morphological literature has been conducted. Based on this and my own studies on flowers in various Apocynaceae and Orchidaceae and complex flowers in other angiosperms with scanning electron microscopy and with microtome section series, a review on convergent floral traits in flower development and architecture in the two families is presented.

KEY FINDINGS

There is a tendency of protracted development of synorganized parts in Apocynaceae and Orchidaceae (development of synorganization of two or more organs begins earlier the more accentuated it is at anthesis). Synorganization (or complexity) also paves the way for novel structures. One of the most conspicuous such novel structures in Apocynaceae is the corona, which is not the product of synorganization of existing organs; however, it is probably enhanced by synorganization of other, existing, floral parts. In contrast to synorganized parts, the corona appears developmentally late.

CONCLUSIONS

Synorganization of floral organs may lead to a large number of convergences in clades that are only very distantly related. The convergences that have been highlighted in this comparative study should be developmentally investigated directly in parallel in future studies.

The effects of subsampling gene trees on coalescent methods applied to ancient divergences

Gene-tree-estimation error is a major concern for coalescent methods of phylogenetic inference. We sampled eight empirical studies of ancient lineages with diverse numbers of taxa and genes for which the original authors applied one or more coalescent methods. We found that the average pairwise congruence among gene trees varied greatly both between studies and also often within a study. We recommend that presenting plots of pairwise congruence among gene trees in a dataset be treated as a standard practice for empirical coalescent studies so that readers can readily assess the extent and distribution of incongruence among gene trees. ASTRAL-based coalescent analyses generally outperformed MP-EST and STAR with respect to both internal consistency (congruence between analyses of subsamples of genes with the complete dataset of all genes) and congruence with the concatenation-based topology. We evaluated the approach of subsampling gene trees that are, on average, more congruent with other gene trees as a method to reduce artifacts caused by gene-tree-estimation errors on coalescent analyses. We suggest that this method is well suited to testing whether gene-tree-estimation error is a primary cause of incongruence between concatenation- and coalescent-based results, to reconciling conflicting phylogenetic results based on different coalescent methods, and to identifying genes affected by artifacts that may then be targeted for reciprocal illumination. We provide scripts that automate the process of calculating pairwise gene-tree incongruence and subsampling trees while accounting for differential taxon sampling among genes. Finally, we assert that multiple tree-search replicates should be implemented as a standard practice for empirical coalescent studies that apply MP-EST.

Sex and the Catasetinae (Darwin's favourite orchids)

Two sexual systems are predominant in Catasetinae (Orchidaceae), namely protandry (which has evolved in other orchid lineages as well) and environmental sex determination (ESD) being a unique trait among Orchidaceae. Yet, the lack of a robust phylogenetic framework for Catasetinae has hampered deeper insights in origin and evolution of sexual systems. To investigate the origins of protandry and ESD in Catasetinae, we sequenced nuclear and chloroplast loci from 77 species, providing the most extensive data matrix of Catasetinae available so far with all major lineages represented. We used Maximum Parsimony, Maximum Likelihood and Bayesian methods to infer phylogenetic relationships and evolution of sexual systems. Irrespectively of the methods used, Catasetinae were monophyletic in molecular phylogenies, with all established generic lineages and their relationships resolved and highly supported. According to comparative reconstruction approaches, the last common ancestor of Catasetinae was inferred as having bisexual flowers (i.e., lacking protandry and ESD as well), and protandry originated once in core Catasetinae (comprising Catasetum, Clowesia, Cycnoches, Dressleria and Mormodes). In addition, three independent gains of ESD are reliably inferred, linked to corresponding loss of protandry within core Catasetinae. Thus, prior gain of protandry appears as the necessary prerequisite for gain of ESD in orchids. Our results contribute to a comprehensive evolutionary scenario for sexual systems in Catasetinae and more generally in orchids as well.

Mining from transcriptomes: 315 single-copy orthologous genes concatenated for the phylogenetic analyses of Orchidaceae

Phylogenetic relationships are hotspots for orchid studies with controversial standpoints. Traditionally, the phylogenies of orchids are based on morphology and subjective factors. Although more reliable than classic phylogenic analyses, the current methods are based on a few gene markers and PCR amplification, which are labor intensive and cannot identify the placement of some species with degenerated plastid genomes. Therefore, a more efficient, labor-saving and reliable method is needed for phylogenic analysis. Here, we present a method of orchid phylogeny construction using transcriptomes. Ten representative species covering five subfamilies of Orchidaceae were selected, and 315 single-copy orthologous genes extracted from the transcriptomes of these organisms were applied to reconstruct a more robust phylogeny of orchids. This approach provided a rapid and reliable method of phylogeny construction for Orchidaceae, one of the most diversified family of angiosperms. We also showed the rigorous systematic position of holomycotrophic species, which has previously been difficult to determine because of the degenerated plastid genome. We concluded that the method presented in this study is more efficient and reliable than methods based on a few gene markers for phylogenic analyses, especially for the holomycotrophic species or those whose DNA sequences have been difficult to amplify. Meanwhile, a total of 315 single-copy orthologous genes of orchids are offered and more informative loci could be used in the future orchid phylogenetic studies.

Phylogenetic placement and taxonomy of the genus Hederorkis (Orchidaceae)

Three plastid regions, matK, rpl32-trnL and rpl16 intron and the ITS1-5.8S-ITS2 nuclear ribosomal DNA were used to demonstrate a phylogenetic placement of the genus Hederorkis (Orchidaceae) for the first time. The taxonomic position of this genus has been unclear thus far. The phylogenetic and morphological relations of Hederorkis to the most closely related genera Sirhookera, Adrorhizon, Bromheadia and Polystachya are also discussed. A hypothesis concerning an origin and evolution of Hederorkis is proposed. Hederorkis is an epiphytic two-leaved orchid genus with lateral inflorescence, non-resupinate flowers, elongate gynostemium and rudimentary column foot. It is native to the Indian Ocean Islands. Two species of Hederorkis are recognized worldwide, H. scandens endemic to Mauritius and Réunion and H. seychellensis endemic to Seychelles. For each of the species treated a full synonymy, detailed description and illustration are included. The distribution map and dichotomous keys to the species have also been provided.

A new myco-heterotrophic genus, Yunorchis, and the molecular phylogenetic relationships of the tribe Calypsoeae (Epidendroideae, Orchidaceae) inferred from plastid and nuclear DNA sequences

We identified a new holomycotrophic orchid that is related to the myco-heterotrophic Calypsoeae. Because chloroplast genes are primarily lacking or are highly divergent, key morphological characters are either reduced or lost from many myco-heterotrophs, and the phylogenetic relationships of weakly supported paraphyletic Calypsoeae within Epidendroideae have been poorly understood in previous molecular systematic studies. Using chloroplast rbcL, psaB, and matK and nuclear Xdh and ITS sequences, we determined the circumscription and systematic positions of the new orchid and the tribe. The results indicate that the epidendroid taxa include most of the clades that are successively sister to the grade of clades representing previously recognized tribes. Calypsoeae comprising four well-supported clades with 12 genera (except for the previous temporarily placed Wullschlaegelia) is supported as a monophyletic and sister clade to Epidendreae (excluding Coeliinae). The new orchid is nested in Calypsoeae and is a sister to Dactylostalix and/or Calypso. This new holomycotrophic orchid presents a subumbel inflorescence that grows underground, and flower with a long pedicel reputing the ground to open and two fragments at the base of the hook, which are obviously morphologically different from those of Calypsoeae. To accommodate this species in the current generic circumscription, a new genus Yunorchis was created.

Phylogenetic relationships in Epidendroideae (Orchidaceae), one of the great flowering plant radiations: progressive specialization and diversification

BACKGROUND AND AIMS

The largest subfamily of orchids, Epidendroideae, represents one of the most significant diversifications among flowering plants in terms of pollination strategy, vegetative adaptation and number of species. Although many groups in the subfamily have been resolved, significant relationships in the tree remain unclear, limiting conclusions about diversification and creating uncertainty in the classification. This study brings together DNA sequences from nuclear, plastid and mitochrondrial genomes in order to clarify relationships, to test associations of key characters with diversification and to improve the classification.

METHODS

Sequences from seven loci were concatenated in a supermatrix analysis for 312 genera representing most of epidendroid diversity. Maximum-likelihood and parsimony analyses were performed on this matrix and on subsets of the data to generate trees and to investigate the effect of missing values. Statistical character-associated diversification analyses were performed.

KEY RESULTS

Likelihood and parsimony analyses yielded highly resolved trees that are in strong agreement and show significant support for many key clades. Many previously proposed relationships among tribes and subtribes are supported, and some new relationships are revealed. Analyses of subsets of the data suggest that the relatively high number of missing data for the full analysis is not problematic. Diversification analyses show that epiphytism is most strongly associated with diversification among epidendroids, followed by expansion into the New World and anther characters that are involved with pollinator specificity, namely early anther inflexion, cellular pollinium stalks and the superposed pollinium arrangement.

CONCLUSIONS

All tested characters show significant association with speciation in Epidendroideae, suggesting that no single character accounts for the success of this group. Rather, it appears that a succession of key features appeared that have contributed to diversification, sometimes in parallel.

The velamen protects photosynthetic orchid roots against UV-B damage, and a large dated phylogeny implies multiple gains and losses of this function during the Cenozoic

UV-B radiation damage in leaves is prevented by epidermal UV-screening compounds that can be modulated throughout ontogeny. In epiphytic orchids, roots need to be protected against UV-B because they photosynthesize, sometimes even replacing the leaves. How orchid roots, which are covered by a dead tissue called velamen, avoid UV-B radiation is currently unknown. We tested for a UV-B protective function of the velamen using gene expression analyses, mass spectrometry, histochemistry, and chlorophyll fluorescence in Phalaenopsis × hybrida roots. We also investigated its evolution using comparative phylogenetic methods. Our data show that two paralogues of the chalcone synthase (CHS) gene family are UV-B-induced in orchid root tips, triggering the accumulation of two UV-B-absorbing flavonoids and resulting in effective protection of the photosynthetic root cortex. Phylogenetic and dating analyses imply that the two CHS lineages duplicated c. 100 million yr before the rise of epiphytic orchids. These findings indicate an additional role for the epiphytic orchid velamen previously thought to function solely in absorbing water and nutrients. This new function, which fundamentally differs from the mechanism of UV-B avoidance in leaves, arose following an ancient duplication of CHS, and has probably contributed to the family's expansion into the canopy during the Cenozoic.

Taxonomy and Biogeography of Apomixis in Angiosperms and Associated Biodiversity Characteristics

Apomixis in angiosperms is asexual reproduction from seed. Its importance to angiospermous evolution and biodiversity has been difficult to assess mainly because of insufficient taxonomic documentation. Thus, we assembled literature reporting apomixis occurrences among angiosperms and transferred the information to an internet database (http://www.apomixis.uni-goettingen.de). We then searched for correlations between apomixis occurrences and well-established measures of taxonomic diversity and biogeography. Apomixis was found to be taxonomically widespread with no clear tendency to specific groups and to occur with sexuality at all taxonomic levels. Adventitious embryony was the most frequent form (148 genera) followed by apospory (110) and diplospory (68). All three forms are phylogenetically scattered, but this scattering is strongly associated with measures of biodiversity. Across apomictic-containing orders and families, numbers of apomict-containing genera were positively correlated with total numbers of genera. In general, apomict-containing orders, families, and subfamilies of Asteraceae, Poaceae, and Orchidaceae were larger, i.e., they possessed more families or genera, than non-apomict-containing orders, families or subfamilies. Furthermore, many apomict-containing genera were found to be highly cosmopolitan. In this respect, 62% occupy multiple geographic zones. Numbers of genera containing sporophytic or gametophytic apomicts decreased from the tropics to the arctic, a trend that parallels general biodiversity. While angiosperms appear to be predisposed to shift from sex to apomixis, there is also evidence of reversions to sexuality. Such reversions may result from genetic or epigenetic destabilization events accompanying hybridization, polyploidy, or other cytogenetic alterations. Because of increased within-plant genetic and genomic heterogeneity, range expansions and diversifications at the species and genus levels may occur more rapidly upon reversion to sexuality. The significantly-enriched representations of apomicts among highly diverse and geographically-extensive taxa, from genera to orders, support this conclusion.

Comparative chloroplast genomes of photosynthetic orchids: insights into evolution of the Orchidaceae and development of molecular markers for phylogenetic applications

The orchid family Orchidaceae is one of the largest angiosperm families, including many species of important economic value. While chloroplast genomes are very informative for systematics and species identification, there is very limited information available on chloroplast genomes in the Orchidaceae. Here, we report the complete chloroplast genomes of the medicinal plant Dendrobium officinale and the ornamental orchid Cypripedium macranthos, demonstrating their gene content and order and potential RNA editing sites. The chloroplast genomes of the above two species and five known photosynthetic orchids showed similarities in structure as well as gene order and content, but differences in the organization of the inverted repeat/small single-copy junction and ndh genes. The organization of the inverted repeat/small single-copy junctions in the chloroplast genomes of these orchids was classified into four types; we propose that inverted repeats flanking the small single-copy region underwent expansion or contraction among Orchidaceae. The AT-rich regions of the ycf1 gene in orchids could be linked to the recombination of inverted repeat/small single-copy junctions. Relative species in orchids displayed similar patterns of variation in ndh gene contents. Furthermore, fifteen highly divergent protein-coding genes were identified, which are useful for phylogenetic analyses in orchids. To test the efficiency of these genes serving as markers in phylogenetic analyses, coding regions of four genes (accD, ccsA, matK, and ycf1) were used as a case study to construct phylogenetic trees in the subfamily Epidendroideae. High support was obtained for placement of previously unlocated subtribes Collabiinae and Dendrobiinae in the subfamily Epidendroideae. Our findings expand understanding of the diversity of orchid chloroplast genomes and provide a reference for study of the molecular systematics of this family.

Phylogenetics of tribe Collabieae (Orchidaceae, Epidendroideae) based on four chloroplast genes with morphological appraisal

Collabieae (Orchidaceae) is a long neglected tribe with confusing tribal and generic delimitation and little-understood phylogenetic relationships. Using plastid matK, psaB, rbcL, and trnH-psbA DNA sequences and morphological evidence, the phylogenetic relationships within the tribe Collabieae were assessed as a basis for revising their tribal and generic delimitation. Collabieae (including the previously misplaced mycoheterotrophic Risleya) is supported as monophyletic and nested within a superclade that also includes Epidendreae, Podochileae, Cymbidieae and Vandeae. Risleya is nested in Collabiinae and sister to Chrysoglossum, a relationship which, despite their great vegetative differences, is supported by floral characters. Ania is a distinct genus supported by both morphological and molecular evidence, while redefined Tainia includes Nephelaphyllum and Mischobulbum. Calanthe is paraphyletic and consists four clades; the genera Gastrorchis, Phaius and Cephalantheropsis should be subsumed within Calanthe. Calanthe sect. Ghiesbreghtia is nested within sect. Calanthe, to which the disputed Calanthe delavayi belongs as well. Our results indicate that, in Collabieae, habit evolved from being epiphytic to terrestrial.

Resolving ancient radiations: can complete plastid gene sets elucidate deep relationships among the tropical gingers (Zingiberales)?

BACKGROUND AND AIMS

Zingiberales comprise a clade of eight tropical monocot families including approx. 2500 species and are hypothesized to have undergone an ancient, rapid radiation during the Cretaceous. Zingiberales display substantial variation in floral morphology, and several members are ecologically and economically important. Deep phylogenetic relationships among primary lineages of Zingiberales have proved difficult to resolve in previous studies, representing a key region of uncertainty in the monocot tree of life.

METHODS

Next-generation sequencing was used to construct complete plastid gene sets for nine taxa of Zingiberales, which were added to five previously sequenced sets in an attempt to resolve deep relationships among families in the order. Variation in taxon sampling, process partition inclusion and partition model parameters were examined to assess their effects on topology and support.

KEY RESULTS

Codon-based likelihood analysis identified a strongly supported clade of ((Cannaceae, Marantaceae), (Costaceae, Zingiberaceae)), sister to (Musaceae, (Lowiaceae, Strelitziaceae)), collectively sister to Heliconiaceae. However, the deepest divergences in this phylogenetic analysis comprised short branches with weak support. Additionally, manipulation of matrices resulted in differing deep topologies in an unpredictable fashion. Alternative topology testing allowed statistical rejection of some of the topologies. Saturation fails to explain observed topological uncertainty and low support at the base of Zingiberales. Evidence for conflict among the plastid data was based on a support metric that accounts for conflicting resampled topologies.

CONCLUSIONS

Many relationships were resolved with robust support, but the paucity of character information supporting the deepest nodes and the existence of conflict suggest that plastid coding regions are insufficient to resolve and support the earliest divergences among families of Zingiberales. Whole plastomes will continue to be highly useful in plant phylogenetics, but the current study adds to a growing body of literature suggesting that they may not provide enough character information for resolving ancient, rapid radiations.

DNA barcoding of Orchidaceae in Korea

Species of Orchidaceae are under severe threat of extinction mainly due to overcollection and habitat destruction; accurate identification of orchid species is critical in conservation biology and sustainable utilization of orchids as plant resources. We examined 647 sequences of the cpDNA regions rbcL, matK, atpF-atpH IGS, psbK-psbI IGS and trnH-psbA IGS from 89 orchid species (95 taxa) and four outgroup taxa to develop an efficient DNA barcode for Orchidaceae in Korea. The five cpDNA barcode regions were successfully amplified and sequenced for all chlorophyllous taxa, but the amplification and sequencing of the same regions in achlorophyllous taxa produced variable results. psbK-psbI IGS showed the highest mean interspecific K2P distance (0.1192), followed by matK (0.0803), atpF-atpH IGS (0.0648), trnH-psbA IGS (0.0460) and rbcL (0.0248). The degree of species resolution for individual barcode regions ranged from 60.5% (rbcL) to 83.5% (trnH-psbA IGS). The degree of species resolution was significantly enhanced in multiregion combinations of the five barcode regions. Of the 26 possible combinations of the five regions, six provided the highest degree of species resolution (98.8%). Among these, a combination of atpF-atpH IGS, psbK-psbI IGS and trnH-psbA IGS, which comprises the least number of DNA regions, is the best option for barcoding of the Korean orchid species.

Complete chloroplast genome of the genus Cymbidium: lights into the species identification, phylogenetic implications and population genetic analyses

BACKGROUND

Cymbidium orchids, including some 50 species, are the famous flowers, and they possess high commercial value in the floricultural industry. Furthermore, the values of different orchids are great differences. However, species identification is very difficult. To a certain degree, chloroplast DNA sequence data are a versatile tool for species identification and phylogenetic implications in plants. Different chloroplast loci have been utilized for evaluating phylogenetic relationships at each classification level among plant species, including at the interspecies and intraspecies levels. However, there is no evidence that a short sequence can distinguish all plant species from each other in order to infer phylogenetic relationships. Molecular markers derived from the complete chloroplast genome can provide effective tools for species identification and phylogenetic resolution.

RESULTS

The complete nucleotide sequences of eight individuals from a total of five Cymbidium species' chloroplast (cp) genomes were determined using Illumina sequencing technology of the total DNA via a combination of de novo and reference-guided assembly. The length of the Cymbidium cp genome is about 155 kb. The cp genomes contain 123 unique genes, and the IR regions contain 24 duplicates. Although the genomes, including genome structure, gene order and orientation, are similar to those of other orchids, they are not evolutionarily conservative. The cp genome of Cymbidium evolved moderately with more than 3% sequence divergence, which could provide enough information for phylogeny. Rapidly evolving chloroplast genome regions were identified and 11 new divergence hotspot regions were disclosed for further phylogenetic study and species identification in Orchidaceae.

CONCLUSIONS

Phylogenomic analyses were conducted using 10 complete chloroplast genomes from seven orchid species. These data accurately identified the individuals and established the phylogenetic relationships between the species. The results reveal that phylogenomics based on organelle genome sequencing lights the species identification-organelle-scale "barcodes", and is also an effective approach for studying whole populations and phylogenetic characteristics of Cymbidium.

Evolution and biogeography of the slipper orchids: Eocene vicariance of the conduplicate genera in the Old and New World Tropics

Intercontinental disjunctions between tropical regions, which harbor two-thirds of the flowering plants, have drawn great interest from biologists and biogeographers. Most previous studies on these distribution patterns focused on woody plants, and paid little attention to herbs. The Orchidaceae is one of the largest families of angiosperms, with a herbaceous habit and a high species diversity in the Tropics. Here we investigate the evolutionary and biogeographical history of the slipper orchids, which represents a monophyletic subfamily (Cypripedioideae) of the orchid family and comprises five genera that are disjunctly distributed in tropical to temperate regions. A relatively well-resolved and highly supported phylogeny of slipper orchids was reconstructed based on sequence analyses of six maternally inherited chloroplast and two low-copy nuclear genes (LFY and ACO). We found that the genus Cypripedium with a wide distribution in the northern temperate and subtropical zones diverged first, followed by Selenipedium endemic to South America, and finally conduplicate-leaved genera in the Tropics. Mexipedium and Phragmipedium from the neotropics are most closely related, and form a clade sister to Paphiopedilum from tropical Asia. According to molecular clock estimates, the genus Selenipedium originated in Palaeocene, while the most recent common ancestor of conduplicate-leaved slipper orchids could be dated back to the Eocene. Ancestral area reconstruction indicates that vicariance is responsible for the disjunct distribution of conduplicate slipper orchids in palaeotropical and neotropical regions. Our study sheds some light on mechanisms underlying generic and species diversification in the orchid family and tropical disjunctions of herbaceous plant groups. In addition, we suggest that the biogeographical study should sample both regional endemics and their widespread relatives.

Assessment of phylogenetic inter-relationships in the genus Cymbidium (Orchidaceae) based on internal transcribed spacer region of rDNA

Sequence data obtained from nrITS region were used to assess phylogenetic inter-relationships and infrageneric classification of ten Cymbidium species collected from north-east India. The final aligned data matrix of combined ITS 1, 5.8S and ITS 2 yielded 684 characters. The ITS 1 and ITS 2 regions showed variable sequence lengths and G+C content (%). The 5.8S region was found to be more conserved (98.71%) followed by ITS 1 (86.12%) and ITS 2 (69.40%). ITS 2 recorded highest percentage of parsimony informative sites (7.46%), high sequence divergence with indels (24.63%), high number of transitions and transversions. ITS sequence data determined the phylogeny of Asiatic Cymbidiums with high bootstrap values. All three proposed subgenera could be distinguished clearly by all four (MP, ML, NJ, and BI) phylogenetic methods. This study validates the utility of ITS rDNA region as a reliable indicator of phylogenetic relationships, especially ITS 2 as probable DNA barcode at higher levels and can serve as an additional approach for identification of broader range of plant taxa especially orchids.

Radical instability and spurious branch support by likelihood when applied to matrices with non-random distributions of missing data

Non-random distributions of missing data are a general problem for likelihood-based statistical analyses, including those in a phylogenetic context. Extensive non-randomly distributed missing data are particularly problematic in supermatrix analyses that include many terminals and/or loci. It has been widely reported that missing data can lead to loss of resolution, but only very rarely create misleading or otherwise unsupported results in a parsimony context. Yet this does not hold for all parametric-based analyses because of their assumption of homogeneity across characters and lineages, which can lead to both long-branch attraction and long-branch repulsion. Contrived examples were used to demonstrate that non-random distributions of missing data, even without rate heterogeneity among characters and a well fitting model, can provide misleading likelihood-based topologies and branch-support values that are radically unstable based on slight modifications to character sampling. The same can occur despite complete absence of parsimony-informative characters. Otherwise unsupported resolution and high branch support for these clades were found to occur frequently in 22 empirical examples derived from a published supermatrix. Partitioning characters based on the distribution of missing data helped to decrease, but did not eliminate, these artifacts. These artifacts were exacerbated by low quality tree searches, particularly when holding only a single optimal tree that must be fully resolved.

Paraholcoglossum and Tsiorchis, two new orchid genera established by molecular and morphological analyses of the Holcoglossum alliance

BACKGROUND

Holcoglossum is a small orchid genus of 12 species ranging from SW China to Thailand and NE India. Although molecular and morphological analyses have been performed to establish the phylogenetic relationships within this genus, the interspecific relations and its relations with allied genera, such as Rhynchostylis, Aerides and Vanda, remain unclear.

METHODOLOGY/PRINCIPAL FINDINGS

In addition to morphological analysis, maximum parsimony, maximum likelihood, and Bayesian inference analyses were performed based on fragments of the nuclear ITS and chloroplast trnL-F and matK genes of 31 taxa (15 Holcoglossum, 14 Aeridinae, 2 outgroups) representing all major clades of the Holcoglossum alliance. The results suggest that Holcoglossum is triphyletic, comprising three clades: the Holcoglossum clade, its sister clade, and a distant clade more closely related to Rhynchostylis, Aerides, and Vanda than to the Holcoglossum clade. The Holcoglossum clade is further divided into three subclades; the genetic distances between these three subclades also support this delimitation. The molecular conclusion is consistent with their distinct morphological characters.

CONCLUSIONS

We propose that the latter two clades comprise two new genera, Paraholcoglossum and Tsiorchis, and Holcoglossum clade divides into three sections. In addition, a new section, Holcoglossum sect. Nujiangensia, and a new species, Holcoglossum linearifolium, are proposed. Some new combinations are made, and a new scheme is provided for the classification of all species of Holcoglossum, Paraholcoglossum, and Tsiorchis.

Rate accelerations in nuclear 18S rDNA of mycoheterotrophic and parasitic angiosperms

Rate variation in genes from all three genomes has been observed frequently in plant lineages with a parasitic and mycoheterotrophic mode of life. While the loss of photosynthetic ability leads to a relaxation of evolutionary constraints in genes involved in the photosynthetic apparatus, it remains to be determined how prevalent increased substitution rates are in nuclear DNA of non-photosynthetic angiosperms. In this study we infer rates of molecular evolution of 18S rDNA of all parasitic and mycoheterotorphic plant families (except Lauraceae and Polygalaceae) using relative rate tests. In several holoparasitic and mycoheterotrophic plant lineages extremely high substitution rates are observed compared to other photosynthetic angiosperms. The position and frequency of these substitutions have been identified to understand the mutation dynamics of 18S rRNA in achlorophyllous plants. Despite the presence of significantly elevated substitution rates, very few mutations occur in major functional and structural regions of the small ribosomal molecule, providing evidence that the efficiency of the translational apparatus in non-photosynthetic plants has not been affected.

DNA barcode sequence identification incorporating taxonomic hierarchy and within taxon variability

For DNA barcoding to succeed as a scientific endeavor an accurate and expeditious query sequence identification method is needed. Although a global multiple-sequence alignment can be generated for some barcoding markers (e.g. COI, rbcL), not all barcoding markers are as structurally conserved (e.g. matK). Thus, algorithms that depend on global multiple-sequence alignments are not universally applicable. Some sequence identification methods that use local pairwise alignments (e.g. BLAST) are unable to accurately differentiate between highly similar sequences and are not designed to cope with hierarchic phylogenetic relationships or within taxon variability. Here, I present a novel alignment-free sequence identification algorithm--BRONX--that accounts for observed within taxon variability and hierarchic relationships among taxa. BRONX identifies short variable segments and corresponding invariant flanking regions in reference sequences. These flanking regions are used to score variable regions in the query sequence without the production of a global multiple-sequence alignment. By incorporating observed within taxon variability into the scoring procedure, misidentifications arising from shared alleles/haplotypes are minimized. An explicit treatment of more inclusive terminals allows for separate identifications to be made for each taxonomic level and/or for user-defined terminals. BRONX performs better than all other methods when there is imperfect overlap between query and reference sequences (e.g. mini-barcode queries against a full-length barcode database). BRONX consistently produced better identifications at the genus-level for all query types.

Complete chloroplast genome of Oncidium Gower Ramsey and evaluation of molecular markers for identification and breeding in Oncidiinae

BACKGROUND

Oncidium spp. produce commercially important orchid cut flowers. However, they are amenable to intergeneric and inter-specific crossing making phylogenetic identification very difficult. Molecular markers derived from the chloroplast genome can provide useful tools for phylogenetic resolution.

RESULTS

The complete chloroplast genome of the economically important Oncidium variety Onc. Gower Ramsey (Accession no. GQ324949) was determined using a polymerase chain reaction (PCR) and Sanger based ABI sequencing. The length of the Oncidium chloroplast genome is 146,484 bp. Genome structure, gene order and orientation are similar to Phalaenopsis, but differ from typical Poaceae, other monocots for which there are several published chloroplast (cp) genome. The Onc. Gower Ramsey chloroplast-encoded NADH dehydrogenase (ndh) genes, except ndhE, lack apparent functions. Deletion and other types of mutations were also found in the ndh genes of 15 other economically important Oncidiinae varieties, except ndhE in some species. The positions of some species in the evolution and taxonomy of Oncidiinae are difficult to identify. To identify the relationships between the 15 Oncidiinae hybrids, eight regions of the Onc. Gower Ramsey chloroplast genome were amplified by PCR for phylogenetic analysis. A total of 7042 bp derived from the eight regions could identify the relationships at the species level, which were supported by high bootstrap values. One particular 1846 bp region, derived from two PCR products (trnHGUG -psbA and trnFGAA-ndhJ) was adequate for correct phylogenetic placement of 13 of the 15 varieties (with the exception of Degarmoara Flying High and Odontoglossum Violetta von Holm). Thus the chloroplast genome provides a useful molecular marker for species identifications.

CONCLUSION

In this report, we used Phalaenopsis. aphrodite as a prototype for primer design to complete the Onc. Gower Ramsey genome sequence. Gene annotation showed that most of the ndh genes inOncidiinae, with the exception of ndhE, are non-functional. This phenomenon was observed in all of the Oncidiinae species tested. The genes and chloroplast DNA regions that would be the most useful for phylogenetic analysis were determined to be the trnHGUG-psbA and the trnFGAA-ndhJ regions. We conclude that complete chloroplast genome information is useful for plant phylogenetic and evolutionary studies in Oncidium with applications for breeding and variety identification.

Evidence of transoceanic dispersion of the genus Vanilla based on plastid DNA phylogenetic analysis

The phylogeny and the biogeographical history of the genus Vanilla was investigated using four chloroplastic genes (psbB, psbC; psaB and rbcL), on 47 accessions of Vanilla chosen from the ex situ CIRAD collection maintained in Reunion Island and additional sequences from GenBank. Bayesian methods provided a fairly well supported reconstruction of the phylogeny of the Vanilloideae sub-family and more particularly of the genus Vanilla. Three major phylogenetic groups in the genus Vanilla were differentiated, which is in disagreement with the actual classification in two sections (Foliosae and Aphyllae) based on morphological traits. Recent Bayesian relaxed molecular clock methods allowed to test the two main hypotheses of the phylogeography of the genus Vanilla. Early radiation of the Vanilla genus and diversification by vicariance consecutive to the break-up of Gondwana, 95 million years ago (Mya), was incompatible with the admitted age of origin of Angiosperm. Based on the Vanilloideae age recently estimated to 71 million years ago (Mya), we conclude that the genus Vanilla would have appeared approximately 34 Mya in South America, when continents were already separated. Nevertheless, whatever the two extreme scenarios tested, at least three long distance migration events are needed to explain the present distribution of Vanilla species in tropical areas. These transoceanic dispersions could have occurred via transoceanic passageway such as the Rio Grande Ridge and the involvement of floating vegetation mats and migratory birds.

A phylogenetic study of Laeliinae (Orchidaceae) based on combined nuclear and plastid DNA sequences

BACKGROUND AND AIMS

Laeliinae are a neotropical orchid subtribe with approx. 1500 species in 50 genera. In this study, an attempt is made to assess generic alliances based on molecular phylogenetic analysis of DNA sequence data.

METHODS

Six DNA datasets were gathered: plastid trnL intron, trnL-F spacer, matK gene and trnK introns upstream and dowstream from matK and nuclear ITS rDNA. Data were analysed with maximum parsimony (MP) and Bayesian analysis with mixed models (BA).

KEY RESULTS

Although relationships between Laeliinae and outgroups are well supported, within the subtribe sequence variation is low considering the broad taxonomic range covered. Localized incongruence between the ITS and plastid trees was found. A combined tree followed the ITS trees more closely, but the levels of support obtained with MP were low. The Bayesian analysis recovered more well-supported nodes. The trees from combined MP and BA allowed eight generic alliances to be recognized within Laeliinae, all of which show trends in morphological characters but lack unambiguous synapomorphies.

CONCLUSIONS

By using combined plastid and nuclear DNA data in conjunction with mixed-models Bayesian inference, it is possible to delimit smaller groups within Laeliinae and discuss general patterns of pollination and hybridization compatibility. Furthermore, these small groups can now be used for further detailed studies to explain morphological evolution and diversification patterns within the subtribe.

Molecular phylogenetics of Prescottiinae s.l. and their close allies (Orchidaceae, Cranichideae) inferred from plastid and nuclear ribosomal DNA sequences

The Andes are a cradle of orchid evolution, but most phylogenetic studies of Orchidaceae in this biodiversity hotspot have dealt with epiphytic epidendroid lineages. Here we present a study on neotropical, terrestrial, orchidoid taxa of Prescottiinae s.l. (8 genera, ∼100 species), which are adapted to some of the highest elevation habitats on earth that support orchids. They are currently included within an expanded concept of Cranichidinae in the tribe Cranichideae, but DNA sequence data show that neither Prescottiinae s.l. nor Cranichidinae s.s. are monophyletic. Prescottiinae s.l. consist of two strongly supported lineages: the Altensteinia and Prescottia clades, which have closer affinities to Spiranthinae than to Cranichidinae. The Prescottia clade comprises two well-supported subclades, one including most sampled species of Prescottia and a second one with Pseudocranichis thysanochila sister to Prescottia tubulosa. As a group, they are sister to Spiranthinae. Sister to this pair is the Altensteinia clade comprised of six genera, whose intergeneric relationships are well resolved. Finally, Cranichidinae s.s. is sister to all three of these clades. Morphological and ecological features distinguishing the major groups are discussed, as are potential synapomorphies to define them. The reconstructed phylogeny indicates that the classification of Cranichideae needs to be reexamined.

Positive selection and ancient duplications in the evolution of class B floral homeotic genes of orchids and grasses

BACKGROUND

Positive selection is recognized as the prevalence of nonsynonymous over synonymous substitutions in a gene. Models of the functional evolution of duplicated genes consider neofunctionalization as key to the retention of paralogues. For instance, duplicate transcription factors are specifically retained in plant and animal genomes and both positive selection and transcriptional divergence appear to have played a role in their diversification. However, the relative impact of these two factors has not been systematically evaluated. Class B MADS-box genes, comprising DEF-like and GLO-like genes, encode developmental transcription factors essential for establishment of perianth and male organ identity in the flowers of angiosperms. Here, we contrast the role of positive selection and the known divergence in expression patterns of genes encoding class B-like MADS-box transcription factors from monocots, with emphasis on the family Orchidaceae and the order Poales. Although in the monocots these two groups are highly diverse and have a strongly canalized floral morphology, there is no information on the role of positive selection in the evolution of their distinctive flower morphologies. Published research shows that in Poales, class B-like genes are expressed in stamens and in lodicules, the perianth organs whose identity might also be specified by class B-like genes, like the identity of the inner tepals of their lily-like relatives. In orchids, however, the number and pattern of expression of class B-like genes have greatly diverged.

RESULTS

The DEF-like genes from Orchidaceae form four well-supported, ancient clades of orthologues. In contrast, orchid GLO-like genes form a single clade of ancient orthologues and recent paralogues. DEF-like genes from orchid clade 2 (OMADS3-like genes) are under less stringent purifying selection than the other orchid DEF-like and GLO-like genes. In comparison with orchids, purifying selection was less stringent in DEF-like and GLO-like genes from Poales. Most importantly, positive selection took place before the major organ reduction and losses in the floral axis that eventually yielded the zygomorphic grass floret.

CONCLUSION

In DEF-like genes of Poales, positive selection on the region mediating interactions with other proteins or DNA could have triggered the evolution of the regulatory mechanisms behind the development of grass-specific reproductive structures. Orchidaceae show a different trend, where gene duplication and transcriptional divergence appear to have played a major role in the canalization and modularization of perianth development.

MADS about the evolution of orchid flowers

Orchids have unique flowers involving three types of perianth organs: outer tepals, lateral inner tepals, and a lip. Expression studies indicate that the identity of these organs is specified by the combinatorial interaction of four different DEFICIENS-like MADS-box genes. We suggest that clarifying the evolution of these genes provides a rational framework for reconstructing the enigmatic origin and unique diversification of the orchid flower. For example, two rounds of gene duplications during early orchid evolution might have generated the genes that were probably recruited to distinguish the different types of perianth organs. This hypothesis suggests intriguing, experimentally testable mechanisms by which gene duplications followed by sub- and neo-functionalization events might have contributed to the evolutionary origin of morphological novelties in orchids - and well beyond.