Reticulate evolution and sea-level fluctuations together drove species diversification of slipper orchids (Paphiopedilum) in South-East Asia

To and fro in the archipelago: Repeated inter-island dispersal and New Guinea's orogeny affect diversification of Delias, the world's largest butterfly genus

The world's largest butterfly genus Delias, commonly known as Jezebels, comprises ca. 251 species found throughout Asia, Australia, and Melanesia. Most species are endemic to islands in the Indo-Australian Archipelago or to New Guinea and nearby islands in Melanesia, and many species are restricted to montane habitats over 1200 m. We inferred an extensively sampled and well-supported molecular phylogeny of the group to better understand the spatial and temporal dimensions of its diversification. The remarkable diversity of Delias evolved in just ca. 15-16 Myr (crown age). The most recent common ancestor of a clade with most of the species dispersed out of New Guinea ca. 14 Mya, but at least six subsequently diverging lineages dispersed back to the island. Diversification was associated with frequent dispersal of lineages among the islands of the Indo-Australian Archipelago, and the divergence of sister taxa on a single landmass was rare and occurred only on the largest islands, most notably on New Guinea. We conclude that frequent inter-island dispersal during the Neogene-likely facilitated by frequent sea level change-sparked much diversification during that period. Many extant New Guinea lineages started diversifying 5 Mya, suggesting that orogeny facilitated their diversification. Our results largely agree with the most recently proposed species group classification system, and we use our large taxon sample to extend this system to all described species. Finally, we summarize recent insights to speculate how wing pattern evolution, mimicry, and sexual selection might also contribute to these butterflies' rapid speciation and diversification.

Comparative and phylogenetic analysis of Chiloschista (Orchidaceae) species and DNA barcoding investigation based on plastid genomes

BACKGROUND

Chiloschista (Orchidaceae, Aeridinae) is an epiphytic leafless orchid that is mainly distributed in tropical or subtropical forest canopies. This rare and threatened orchid lacks molecular resources for phylogenetic and barcoding analysis. Therefore, we sequenced and assembled seven complete plastomes of Chiloschista to analyse the plastome characteristics and phylogenetic relationships and conduct a barcoding investigation.

RESULTS

We are the first to publish seven Chiloschista plastomes, which possessed the typical quadripartite structure and ranged from 143,233 bp to 145,463 bp in size. The plastomes all contained 120 genes, consisting of 74 protein-coding genes, 38 tRNA genes and eight rRNA genes. The ndh genes were pseudogenes or lost in the genus, and the genes petG and psbF were under positive selection. The seven Chiloschista plastomes displayed stable plastome structures with no large inversions or rearrangements. A total of 14 small inversions (SIs) were identified in the seven Chiloschista plastomes but were all similar within the genus. Six noncoding mutational hotspots (trnNGUU-rpl32 > rpoB-trnCGCA > psbK-psbI > psaC-rps15 > trnEUUC-trnTGGU > accD-psaI) and five coding sequences (ycf1 > rps15 > matK > psbK > ccsA) were selected as potential barcodes based on nucleotide diversity and species discrimination analysis, which suggested that the potential barcode ycf1 was most suitable for species discrimination. A total of 47-56 SSRs and 11-14 long repeats (> 20 bp) were identified in Chiloschista plastomes, and they were mostly located in the large single copy intergenic region. Phylogenetic analysis indicated that Chiloschista was monophyletic. It was clustered with Phalaenopsis and formed the basic clade of the subtribe Aeridinae with a moderate support value. The results also showed that seven Chiloschista species were divided into three major clades with full support.

CONCLUSION

This study was the first to analyse the plastome characteristics of the genus Chiloschista in Orchidaceae, and the results showed that Chiloschista plastomes have conserved plastome structures. Based on the plastome hotspots of nucleotide diversity, several genes and noncoding regions are suitable for phylogenetic and population studies. Chiloschista may provide an ideal system to investigate the dynamics of plastome evolution and DNA barcoding investigation for orchid studies.

The relationship between geographic range size and rates of species diversification

Range size is a universal characteristic of every biological species, and is often assumed to affect diversification rate. There are strong theoretical arguments that large-ranged species should have higher rates of diversification. On the other hand, the observation that small-ranged species are often phylogenetically clustered might indicate high diversification of small-ranged species. This discrepancy between theory and the data may be caused by the fact that typical methods of data analysis do not account for range size changes during speciation. Here we use a cladogenetic state-dependent diversification model applied to mammals to show that range size changes during speciation are ubiquitous and small-ranged species indeed diversify generally slower, as theoretically expected. However, both range size and diversification are strongly influenced by idiosyncratic and spatially localized events, such as colonization of an archipelago or a mountain system, which often override the general pattern of range size evolution.

Speciation across the Earth driven by global cooling in terrestrial orchids

Although climate change has been implicated as a major catalyst of diversification, its effects are thought to be inconsistent and much less pervasive than localized climate or the accumulation of species with time. Focused analyses of highly speciose clades are needed in order to disentangle the consequences of climate change, geography, and time. Here, we show that global cooling shapes the biodiversity of terrestrial orchids. Using a phylogeny of 1,475 species of Orchidoideae, the largest terrestrial orchid subfamily, we find that speciation rate is dependent on historic global cooling, not time, tropical distributions, elevation, variation in chromosome number, or other types of historic climate change. Relative to the gradual accumulation of species with time, models specifying speciation driven by historic global cooling are over 700 times more likely. Evidence ratios estimated for 212 other plant and animal groups reveal that terrestrial orchids represent one of the best-supported cases of temperature-spurred speciation yet reported. Employing >2.5 million georeferenced records, we find that global cooling drove contemporaneous diversification in each of the seven major orchid bioregions of the Earth. With current emphasis on understanding and predicting the immediate impacts of global warming, our study provides a clear case study of the long-term impacts of global climate change on biodiversity.

Isolation of Tulasnella spp. from Cultivated Paphiopedilum Orchids and Screening of Germination-Enhancing Fungi

Ex situ conservation, an important way to increase the survival and sustainability of endangered species, is widely used in the conservation of endangered orchids. However, long-term ex situ conservation might affect the dominant group of orchid symbiotic fungi, which are crucial for orchid growth and reintroduction. This study investigated the culturable Tulasnella spp. associated with Paphiopedilum orchids after long-term greenhouse cultivation, and identified germination-enhancing isolates. A total of 44 Tulasnella isolates were obtained from the roots of 14 Paphiopedilum spp., and 29 of them were selected for phylogenetic analysis. They clustered mainly with Tulasnella deliquescens, Tulasnella calospora, Tulasnella bifrons, and Tulasnella irregularis, but included two potential new groups. Compared with published uncultured data, most of the isolates were grouped together with the reported types, and the dominant Tulasnella associated with P. armeniacum and P. micranthum could still be isolated after ten years of cultivation, most of which were the first isolation. In vitro symbiotic germination showed that certain root isolates could promote seed germination (e.g., parm152 isolated from P. armeniacum, Php12 from P. hirsutissimum, and prhi68 from P. rhizomatosum). These data indicated that the dominant Tulasnella types colonizing the roots of cultivated Paphiopedilum are stable over time, and germination-enhancing fungi colonizing the roots would benefit for seed reproduction after population reintroduction into the wild.

Origin and diversification of a Himalayan orchid genus Pleione

Pleione is an orchid endemically distributed in high mountain areas across the Hengduan Mountains (HDM), Himalayas, Southeast Asia and South of China. The unique flower shapes, rich colors and immense medicinal importance of Pleione are valuable ornamental and economic resources. However, the phylogenetic relationships and evolutionary history of the genus have not yet been comprehensively resolved. Here, the evolutionary history of Pleione was investigated using single-copy gene single nucleotide polymorphisms and chloroplast genome datasets. The data revealed that Pleione could be divided into five clades. Discordance in topology between the two phylogenetic trees and network and D-statistic analyses indicated the occurrence of reticulate evolution in the genus. The evolution could be attributed to introgression and incomplete lineage sorting. Ancestral area reconstruction suggested that Pleione was originated from the HDM. Uplifting of the HDM drove rapid diversification by creating conditions favoring rapid speciation. This coincided with two periods of consolidation of the Asian monsoon climate, which caused the first rapid diversification of Pleione from 8.87 to 7.83 Mya, and a second rapid diversification started at around 4.05 Mya to Pleistocene. The interaction between Pleione and climate changes, especially the monsoons, led to the current distribution pattern and shaped the dormancy characteristic of the different clades. In addition to revealing the evolutionary relationship of Pleione with orogeny and climate changes, the findings of this study provide insights into the speciation and diversification mechanisms of plants in the East Asian flora.

Unraveling the macroevolution of horseshoe bats (Chiroptera: Rhinolophidae: Rhinolophus)

Unraveling the diversification mechanisms of organisms is a fundamental and important macroevolutionary question regarding the diversity, ecological niche, and morphological divergence of life. However, many studies have only explored diversification mechanisms via isolated factors. Here, based on comparative phylogenetic analysis, we performed a macroevolutionary examination of horseshoe bats (Chiroptera: Rhinolophidae: Rhinolophus), to reveal the inter-relationships among diversification, intrinsic/extrinsic factors, and climatic ecological niche characteristics. Results showed a general slowing trajectory during diversification, with two dispersal events from Asia into Southeast Asia and Africa playing key roles in shaping regional heterogeneous diversity. Morphospace expansions of the investigated traits (e.g., body size, echolocation, and climate niche) revealed a decoupled pattern between diversification trajectory and trait divergence, suggesting that other factors (e.g., biotic interactions) potentially played a key role in recent diversification. Based on ancestral traits and pathway analyses, most Rhinolophus lineages belonging to the same region overlapped with each other geographically and were positively associated with the diversification rate, implying a competitive prelude to speciation. Overall, our study showed that multiple approaches need to be integrated to address diversification history. Rather than a single factor, the joint effects of multiple factors (biogeography, environmental drivers, and competition) are responsible for the current diversity patterns in horseshoe bats, and a corresponding multifaceted strategy is recommended to study these patterns in the future.

The ecological adaptation of the unparalleled plastome character evolution in slipper orchids

Plastomes may have undergone adaptive evolution in the process of plant adaptation to diverse environments, whereby species may differ in plastome characters. Cypripedioideae successfully colonized distinct environments and could be an ideal group for studying the interspecific variation and adaptive evolution of plastomes. Comparative study of plastomes, ancestral state reconstruction, phylogenetic-based analysis, ecological niche modelling, and selective pressure analysis were conducted to reveal the evolutionary patterns of plastomes in Cypripedioideae and their relationship with environmental factors. The plastomes of the three evolved genera had reduced plastome size, increased GC content, and compacted gene content compared to the basal group. Variations in plastome size and GC content are proved to have clear relationships with climate regions. Furthermore, ecological niche modelling revealed that temperature and water factors are important climatic factors contributing to the distributional difference which is directly correlated with the climate regions. The temperature-sensitive genes ndh genes, infA, and rpl20 were found to be either lost/pseudogenized or under positive selection in the evolved groups. Unparalleled plastome character variations were discovered in slipper orchids. Our study indicates that variations in plastome characters have adaptive consequences and that temperature and water factors are important climatic factors that affect plastome evolution. This research highlights the expectation that plants can facilitate adaptation to different environmental conditions with the changes in plastome and has added critical insight for understanding the process of plastome evolution in plants.

Comparative Analyses of Chloroplast Genomes Provide Comprehensive Insights into the Adaptive Evolution of Paphiopedilum (Orchidaceae)

An elucidation of how the selection pressures caused by habitat environments affect plant plastid genomes and lead to the adaptive evolution of plants, is a very intense area of research in evolutionary biology. The genus Paphiopedilum is a predominant group of orchids that includes over 66 species with high horticultural and ornamental value. However, owing to the destructive exploitation and habitat deterioration of wild germplasm resources of Paphiopedilum, it needs more molecular genetic resources and studies on this genus. The chloroplast is cytoplasmically inherited and often used in evolutionary studies. Thus, for this study, we newly sequenced, assembled and annotated five chloroplast genomes of the Paphiopedilum species. The size of these genomes ranged from 155,886 bp (P. henryanum) to 160,503 bp (P. ‘GZSLKY’ Youyou) and they contained 121–122 genes, which consisted of 76 protein coding genes, eight ribosomal RNAs, and 37–38 transfer RNAs. Combined with the other 14 Paphiopedilum species, the characteristics of the repeat sequences, divergent hotspot regions, and the condo usage bias were evaluated and identified, respectively. The gene transfer analysis showed that some fragments of the ndh and ycf gene families were shared by both the chloroplast and nucleus. Although the genomic structure and gene content was conserved, there was a significant boundary shift caused by the inverted repeat (IR) expansion and small single copy (SSC) contraction. The lower GC content and loss of ndh genes could be the result of adaptive evolutionary responses to its unique habitats. The genes under positive selection, including accD, matK, psbM, rpl20, rps12, ycf1, and ycf2 might be regarded as potential candidate genes for further study, which significantly contribute to the adaptive evolution of Paphiopedilum.

Comparative Analysis of the Complete Chloroplast Genomes of Nine Paphiopedilum Species

Paphiopedilum is known as “lady’s or Venus” slipper orchids due to its prominent shoe-shaped labellum, with high ornamental value. Phylogenetic relationships among some species in Paphiopedilum genus cannot be effectively determined by morphological features alone or through the analysis of nuclear or chloroplast DNA fragments. In order to provide aid in understanding the evolutionary and phylogenetic relationship in Paphiopedilum at chloroplast (cp) genome-scale level, the complete cp genomes of six Paphiopedilum species were newly sequenced in this study, and three other published cp genome sequences of Paphiopedilum were included in the comparative analyses. The cp genomes of the six Paphiopedilum species ranged from 154,908 bp (P. hirsutissimum) to 161,300 bp (P. victoria-mariae) in size, all constituting four-part annular structures. Analyses of the nucleotide substitutions, insertions/deletions, and simple sequence repeats in the cp genomes were conducted. Ten highly variable regions that could serve as potential DNA barcodes or phylogenetic markers for this diverse genus were identified. Sequence variations in the non-coding regions were greater than that in the conserved protein-coding regions, as well as in the large single copy (LSC) and small single copy (SSC) regions than in the inverted repeat (IR) regions. Phylogenetic analysis revealed that all Paphiopedilum species clustered in one monophyletic clade in the Cypripedioideae subfamily and then subdivided into seven smaller branches corresponding to different subgenus or sections of the genus, with high bootstrap supports, indicate that cp genome sequencing can be an effective means in resolving the complex relationship in Paphiopedilum.

Genome Size and Labellum Epidermal Cell Size Are Evolutionarily Correlated With Floral Longevity in Paphiopedilum Species

Genome size is known to influence phenotypic traits in leaves and seeds. Although genome size is closely related to cellular and developmental traits across biological kingdoms, floral longevity is a floral trait with important fitness consequence, but less is known about the link between floral longevity and sizes of genomes and cells. In this study, we examined evolutionary coordination between genome size, floral longevity, and epidermal cell size in flowers and leaves in 13 Paphiopedilum species. We found that, across all the study species, the genome size was positively correlated with floral longevity but negatively associated with labellum epidermal cell size, and a negative relationship was found between floral longevity and labellum epidermal cell size. This suggested that genome size is potentially correlated with floral longevity, and genome size has an important impact on life-history trait. In addition, genome size was positively correlated with leaf epidermal cell size, which was different from the relationship in flower due to different selective pressures they experienced or different functions they performed. Therefore, genome size constraints floral longevity, and it is a strong predictor of cell size. The impact of genome size on reproduction might have more implications for the evolution of flowering plants and pollination ecology.

Phylogenetic incongruence in Cymbidium orchids

Cymbidium, which includes approximately 80 species, is one of the most ornamental and cultivated orchid genera. However, a lack of markers and sparse sampling have posed great challenges to resolving the phylogenetic relationships within the genus. In the present study, we reconstructed the phylogenetic relationships by utilizing one nuclear DNA (nrITS) and seven plastid genes (rbcL, trnS, trnG, matK, trnL, psbA, and atpI) from 70 species (varieties) in Cymbidium. We also examined the occurrence of phylogenetic conflict between nuclear (nrITS) and plastid loci and investigated how phylogenetic conflict bears on taxonomic classification within the genus. We found that phylogenetic conflict and low support values may be explained by hybridization and a lack of informative characteristics. Our results do not support previous classification of the subgenera and sections within Cymbidium. Discordance between gene trees and network analysis indicate that reticulate evolution occurred in the genus Cymbidium. Overall, our study indicates that Cymbidium has undergone a complex evolution.

Species Phylogeny versus Gene Trees: A Case Study of an Incongruent Data Matrix Based on Paphiopedilum Pfitz. (Orchidaceae)

The phylogeny of the genus Paphiopedilum based on the plastome is consistent with morphological analysis. However, to date, none of the analyzed nuclear markers has confirmed this. Topology incongruence among the trees of different nuclear markers concerns entire sections of the subgenus Paphiopedilum. The low-copy nuclear protein-coding gene PHYC was obtained for 22 species representing all sections and subgenera of Paphiopedilum. The nuclear-based phylogeny is supported by morphological characteristics and plastid data analysis. We assumed that an incongruence in nuclear gene trees is caused by ancestral homoploid hybridization. We present a model for inferring the phylogeny of the species despite the incongruence of the different tree topologies. Our analysis, based on six low-copy nuclear genes, is congruent with plastome phylogeny and has been confirmed by phylogenetic network analysis.

The chloroplast genome evolution of Venus slipper (Paphiopedilum): IR expansion, SSC contraction, and highly rearranged SSC regions

BACKGROUND

Paphiopedilum is the largest genus of slipper orchids. Previous studies showed that the phylogenetic relationships of this genus are not well resolved, and sparse taxon sampling documented inverted repeat (IR) expansion and small single copy (SSC) contraction of the chloroplast genomes of Paphiopedilum.

RESULTS

Here, we sequenced, assembled, and annotated 77 plastomes of Paphiopedilum species (size range of 152,130 - 164,092 bp). The phylogeny based on the plastome resolved the relationships of the genus except for the phylogenetic position of two unstable species. We used phylogenetic and comparative genomic approaches to elucidate the plastome evolution of Paphiopedilum. The plastomes of Paphiopedilum have a conserved genome structure and gene content except in the SSC region. The large single copy/inverted repeat (LSC/IR) boundaries are relatively stable, while the boundaries of the inverted repeat and small single copy region (IR/SSC) varied among species. Corresponding to the IR/SSC boundary shifts, the chloroplast genomes of the genus experienced IR expansion and SSC contraction. The IR region incorporated one to six genes of the SSC region. Unexpectedly, great variation in the size, gene order, and gene content of the SSC regions was found, especially in the subg. Parvisepalum. Furthermore, Paphiopedilum provides evidence for the ongoing degradation of the ndh genes in the photoautotrophic plants. The estimated substitution rates of the protein coding genes show accelerated rates of evolution in clpP, psbH, and psbZ. Genes transferred to the IR region due to the boundary shift also have higher substitution rates.

CONCLUSIONS

We found IR expansion and SSC contraction in the chloroplast genomes of Paphiopedilum with dense sampling, and the genus shows variation in the size, gene order, and gene content of the SSC region. This genus provides an ideal system to investigate the dynamics of plastome evolution.

The chloroplast genome evolution of Venus slipper (Paphiopedilum): IR expansion, SSC contraction, and highly rearranged SSC regions

BACKGROUND

Paphiopedilum is the largest genus of slipper orchids. Previous studies showed that the phylogenetic relationships of this genus are not well resolved, and sparse taxon sampling documented inverted repeat (IR) expansion and small single copy (SSC) contraction of the chloroplast genomes of Paphiopedilum.

RESULTS

Here, we sequenced, assembled, and annotated 77 plastomes of Paphiopedilum species (size range of 152,130 - 164,092 bp). The phylogeny based on the plastome resolved the relationships of the genus except for the phylogenetic position of two unstable species. We used phylogenetic and comparative genomic approaches to elucidate the plastome evolution of Paphiopedilum. The plastomes of Paphiopedilum have a conserved genome structure and gene content except in the SSC region. The large single copy/inverted repeat (LSC/IR) boundaries are relatively stable, while the boundaries of the inverted repeat and small single copy region (IR/SSC) varied among species. Corresponding to the IR/SSC boundary shifts, the chloroplast genomes of the genus experienced IR expansion and SSC contraction. The IR region incorporated one to six genes of the SSC region. Unexpectedly, great variation in the size, gene order, and gene content of the SSC regions was found, especially in the subg. Parvisepalum. Furthermore, Paphiopedilum provides evidence for the ongoing degradation of the ndh genes in the photoautotrophic plants. The estimated substitution rates of the protein coding genes show accelerated rates of evolution in clpP, psbH, and psbZ. Genes transferred to the IR region due to the boundary shift also have higher substitution rates.

CONCLUSIONS

We found IR expansion and SSC contraction in the chloroplast genomes of Paphiopedilum with dense sampling, and the genus shows variation in the size, gene order, and gene content of the SSC region. This genus provides an ideal system to investigate the dynamics of plastome evolution.

A comparative plastomics approach reveals available molecular markers for the phylogeographic study of Dendrobium huoshanense, an endangered orchid with extremely small populations

Comparative plastomics approaches have been used to identify available molecular markers for different taxonomic level studies of orchid species. However, the adoption of such methods has been largely limited in phylogeographic studies. Therefore, in this study, Dendrobium huoshanense, an endangered species with extremely small populations, was used as a model system to test whether the comparative plastomic approaches could screen available molecular markers for the phylogeographic study. We sequenced two more plastomes of D. huoshanense and compared them with our previously published one. A total of 27 mutational hotspot regions and six polymorphic cpSSRs have been screened for the phylogeographic studies of D. huoshanense. The cpDNA haplotype data revealed that the existence of haplotype distribution center was located in Dabieshan Mts. (Huoshan). The genetic diversity and phylogenetic analyses showed that the populations of D. huoshanense have been isolated and evolved independently for long period. On the contrary, based on cpSSR data, the genetic structure analysis revealed a mixed structure among the populations in Anhui and Jiangxi province, which suggested that the hybridization or introgression events have occurred among the populations of D. huoshanense. These results indicated that human activities have played key roles in shaping the genetic diversity and distributional patterns of D. huoshanense. According to our results, both two markers showed a high resolution for the phylogeographic studies of D. huoshanense. Therefore, we put forth that comparative plastomic approaches could revealed available molecular markers for phylogeographic study, especially for the species with extremely small populations.

Phylogeny and Historical Biogeography of Paphiopedilum Pfitzer (Orchidaceae) Based on Nuclear and Plastid DNA

The phylogeny and biogeography of the genus Paphiopedilum were evaluated by using phylogenetic trees derived from analysis of nuclear ribosomal internal transcribed spacer (ITS) sequences, the plastid trnL intron, the trnL-F spacer, and the atpB-rbcL spacer. This genus was divided into three subgenera: Parvisepalum, Brachypetalum, and Paphiopedilum. Each of them is monophyletic with high bootstrap supports according to the highly resolved phylogenetic tree reconstructed by combined sequences. There are five sections within the subgenus Paphiopedilum, including Coryopedilum, Pardalopetalum, Cochlopetalum, Paphiopedilum, and Barbata. The subgenus Parvisepalum is phylogenetic basal, which suggesting that Parvisepalum is comprising more ancestral characters than other subgenera. The evolutionary trend of genus Paphiopedilum was deduced based on the maximum likelihood (ML) tree and Bayesian Evolutionary Analysis Sampling Trees (BEAST). Reconstruct Ancestral State in Phylogenies (RASP) analyses based on the combined sequence data. The biogeographic analysis indicates that Paphiopedilum species were firstly derived in Southern China and Southeast Asia, subsequently dispersed into the Southeast Asian archipelagoes. The subgenera Paphiopedilum was likely derived after these historical dispersals and vicariance events. Our research reveals the relevance of the differentiation of Paphiopedilum in Southeast Asia and geological history. Moreover, the biogeographic analysis explains that the significant evolutionary hotspots of these orchids in the Sundaland and Wallacea might be attributed to repeated migration and isolation events between the south-eastern Asia mainland and the Sunda Super Islands.

Genetic Diversity and Identification of Vietnamese Paphiopedilum Species Using DNA Sequences

Paphiopedilum is among the most popular ornamental orchid genera due to its unique slipper flowers and attractive leaf coloration. Most of the Paphiopedilum species are in critical danger due to over-exploitation. They were listed in Appendix I of the Convention on International Trade in Endangered Species of Wild Fauna and Flora, which prevents their being traded across borders. While most Paphiopedilum species are distinctive, owing to their respective flowers, their vegetative features are more similar and undistinguished. Hence, the conservation of these species is challenging, as most traded specimins are immature and non-flowered. An urgent need exists for effective identification methods to prevent further illegal trading of Paphiopedilum species. DNA barcoding is a rapid and sensitive method for species identification, at any developmental stage, using short DNA sequences. In this study, eight loci, i.e., ITS, LEAFY, ACO, matK, trnL, rpoB, rpoC1, and trnH-psbA, were screened for potential barcode sequences on the Vietnamese Paphiopedilum species. In total, 17 out of 22 Paphiopedilum species were well identified. The studied DNA sequences were deposited to GenBank, in which Paphiopedilum dalatense accessions were introduced for the first time. ACO, LEAFY, and trnH-psbA were limited in amplification rate for Paphiopedilum. ITS was the best single barcode. Single ITS could be used along with nucleotide polymorphism characteristics for species discrimination. The combination of ITS + matK was the most efficient identification barcode for Vietnamese Paphiopedilum species. This barcode also succeeded in recognizing misidentified or wrongly-named traded samples. Different bioinformatics programs and algorithms for establishing phylogenetic trees were also compared in the study to propose quick, simple, and effective tools for practical use. It was proved that both the Bayesian Inference method in the MRBAYES program and the neighbor-joining method in the MEGA software met the criteria. Our study provides a barcoding database of Vietnamese Paphiopedilum which may significantly contribute to the control and conservation of these valuable species.

Satellite DNA in Paphiopedilum subgenus Parvisepalum as revealed by high-throughput sequencing and fluorescent in situ hybridization

BACKGROUND

Satellite DNA is a rapidly diverging, largely repetitive DNA component of many eukaryotic genomes. Here we analyse the evolutionary dynamics of a satellite DNA repeat in the genomes of a group of Asian subtropical lady slipper orchids (Paphiopedilum subgenus Parvisepalum and representative species in the other subgenera/sections across the genus). A new satellite repeat in Paphiopedilum subgenus Parvisepalum, SatA, was identified and characterized using the RepeatExplorer pipeline in HiSeq Illumina reads from P. armeniacum (2n = 26). Reconstructed monomers were used to design a satellite-specific fluorescent in situ hybridization (FISH) probe. The data were also analysed within a phylogenetic framework built using the internal transcribed spacer (ITS) sequences of 45S nuclear ribosomal DNA.

RESULTS

SatA comprises c. 14.5% of the P. armeniacum genome and is specific to subgenus Parvisepalum. It is composed of four primary monomers that range from 230 to 359 bp and contains multiple inverted repeat regions with hairpin-loop motifs. A new karyotype of P. vietnamense (2n = 28) is presented and shows that the chromosome number in subgenus Parvisepalum is not conserved at 2n = 26, as previously reported. The physical locations of SatA sequences were visualised on the chromosomes of all seven Paphiopedilum species of subgenus Parvisepalum (2n = 26-28), together with the 5S and 45S rDNA loci using FISH. The SatA repeats were predominantly localisedin the centromeric, peri-centromeric and sub-telocentric chromosome regions, but the exact distribution pattern was species-specific.

CONCLUSIONS

We conclude that the newly discovered, highly abundant and rapidly evolving satellite sequence SatA is specific to Paphiopedilum subgenus Parvisepalum. SatA and rDNA chromosomal distributions are characteristic of species, and comparisons between species reveal that the distribution patterns generate a strong phylogenetic signal. We also conclude that the ancestral chromosome number of subgenus Parvisepalum and indeed of all Paphiopedilum could be either 2n = 26 or 28, if P. vietnamense is sister to all species in the subgenus as suggested by the ITS data.

In Silico Study on Molecular Sequences for Identification of Paphiopedilum Species

Our study searched all available sequences of Paphiopedilum from NCBI (National Center for Biotechnology Information) and tested for their species resolution capability in single as well as in combination forms. A total of 28 loci were applied for analyses in the study. From the nuclear genome, the highest resolution was of LFY, followed by ACO, DEF4, and RAD51. These 4 loci were found to be even better than the popular region ITS for Paphiopedilum identification. Among the chloroplast regions, the intergenic spacer atpB-rbcL gave the highest species resolution (76.7%), followed by matK, trnL, rpoC2, and ycf1. The divergence of CHS, XDH, 18S, Nad1, ccsA, rbcL, and ycf2 was very low and should not be used as identifying markers for Paphiopedilum. In addition, 2-locus combinations could improve significantly the resolving capability for the genus, in which 14/36 data sets could be resolved completely (100%) with interspecies relationships. The indel information was also effective supporting data for molecular discrimination of species.

Development and characterization of polymorphic EST-SSR markers for Paphiopedilum henryanum (Orchidaceae)

PREMISE OF THE STUDY

Microsatellite primers were developed for Paphiopedilum henryanum (Orchidaceae), a species threatened with extinction, to assess genetic diversity and population genetic structure.

METHODS AND RESULTS

Based on the transcriptome data of P. henryanum, 34 novel polymorphic microsatellite expressed sequence tag-simple sequence repeat markers were developed and characterized in 33 individuals from two P. henryanum populations. The results showed the number of alleles per locus ranged from two to four, and levels of observed and expected heterozygosity per locus ranged from 0.000 to 1.000 and from 0.000 to 0.7333, respectively. Of these markers, some primers showed good amplification results in seven other Paphiopedilum species.

CONCLUSIONS

The developed microsatellite markers will be useful in exploring the genetic diversity and structure of P. henryanum. Furthermore, most loci were successfully cross-amplified in seven other species of Paphiopedilum, indicating that they will be of great value for genetic study across this genus.

Taxonomic placement of Paphiopedilum rungsuriyanum (Cypripedioideae; Orchidaceae) based on morphological, cytological and molecular analyses

BACKGROUND

Paphiopedilum rungsuriyanum from Northern Laos was discovered and described in 2014. It is characterized by having miniature tessellated leaves, a flower having a helmet shaped lip with a V-shaped neckline, and a semi-lunate, 3-dentate staminode with an umbo. These morphological features distinguish P. rungsuriyanum from the other known sections/subgenera of Paphiopedilum, making it difficult to group with existing infrageneric units.

RESULTS

Paphiopedilum rungsuriyanum has chromosome number of 2n = 26. Fluorescence in situ hybridization study demonstrates that there are two 45S rDNA signals in the telomeric region of chromosomes, and more than 20 5S rDNA signals dispersed signals in the pericentromeric and centromeric regions. Phylogenetic analyses based on four nuclear (i.e. ITS, ACO, DEF4 and RAD51) and four plastid (i.e. atpI-atpH, matK, trnS-trnfM and ycf1) gene regions indicate that P. rungsuriyanum is nested in subgenus Paphiopedilum and is a sister to section Paphiopedilum.

CONCLUSIONS

The results in combination with karyomorphological, rDNA FISH patterns, morphological and phylogenetic analyses suggest a new section Laosianum to accommodate this species in the current sectional circumscription of subgenus Paphiopedilum.

The molecular basis of differential morphology and bleaching thresholds in two morphs of the coral Pocillopora acuta

Processes of cnidarian evolution, including hybridization and phenotypic plasticity, have complicated the clear diagnosis of species boundaries within the phylum. Pocillopora acuta, a species of scleractinian coral that was recently split from the widespread Pocillopora damicornis species complex, occurs in at least two distinct morphs on the Great Barrier Reef. Contrasting morphology combined with evidence of differential bleaching thresholds among sympatrically distributed colonies suggest that the taxonomy of this recently described species is not fully resolved and may represent its own species complex. To examine the basis of sympatric differentiation between the two morphs, we combined analyses of micro- and macro-skeletal morphology with genome wide sequencing of the coral host, as well as ITS2 genotyping of the associated Symbiodinium communities. We found consistent differences between morphs on both the macro- and micro-skeletal scale. In addition, we identified 18 candidate functional genes that relate to skeletal formation and morphology that may explain how the two morphs regulate growth to achieve their distinct growth forms. With inconclusive results in endosymbiotic algal community diversity between the two morphs, we propose that colony morphology may be linked to bleaching susceptibility. We conclude that cryptic speciation may be in the early stages within the species P. acuta.

Iteration expansion and regional evolution: phylogeography of Dendrobium officinale and four related taxa in southern China

The genus Dendrobium was used as a case study to elucidate the evolutionary history of Orchidaceae in the Sino-Japanese Floristic Region (SJFR) and Southeast Asia region. These evolutionary histories remain largely unknown, including the temporal and spatial distribution of the evolutionary events. The present study used nuclear and plastid DNA to determine the phylogeography of Dendrobium officinale and four closely related taxa. Plastid DNA haplotype and nuclear data were shown to be discordant, suggesting reticulate evolution drove the species' diversification. Rapid radiation and genetic drift appeared to drive the evolution of D. tosaense and D. flexicaule, whereas introgression or hybridization might have been involved in the evolution of D. scoriarum and D. shixingense. The phylogeographical structure of D. officinale revealed that core natural distribution regions might have served as its glacial refuges. In recent years, human disturbances caused its artificial migration and population extinction. The five taxa may have originated from the Nanling Mountains and the Yungui Plateau and then migrated northward or eastward. After the initial iteration expansion, D. officinale populations appeared to experience the regional evolutionary patterns in different regions and follow the sequential or rapid decline in gene exchange.

On the persistence of reproductive barriers in Eucalyptus: the bridging of mechanical barriers to zygote formation by F1 hybrids is counteracted by intrinsic post-zygotic incompatibilities

BACKGROUND AND AIMS

Many previous studies conclude that pre-zygotic barriers such as mechanical isolation account for most reproductive isolation between pairs of taxa. However, the inheritance and persistence of barriers such as these after the first generation of hybridization is rarely quantified, even though it is a vital consideration in understanding gene flow potential. There is an asymmetrical pre-zygotic mechanical barrier to hybridization between Eucalyptus nitens and Eucalyptus globulus, which completely prevents small-flowered E. nitens pollen from mating with large E. globulus flowers, while the reverse cross is possible. We aimed to determine the relative importance of pre- and post-zygotic barriers in preventing gene flow following secondary contact between E. nitens and E. globulus, including the inheritance of barriers in advanced-generation hybrids.

METHODS

Experimental crossing was used to produce outcrossed E. nitens, E. globulus and their F1, F2, BCg and BCn hybrids. The strength and inheritance of a suite of pre- and post-zygotic barriers were assessed, including 20-year survival, growth and reproductive capacity.

KEY RESULTS

The mechanical barrier to hybridization was lost or greatly reduced in the F1 hybrid. In contrast, intrinsic post-zygotic barriers were strong and persistent. Line-cross analysis indicated that the outbreeding depression in the hybrids was best explained by epistatic loss.

CONCLUSIONS

The removal of strong mechanical barriers between E. nitens and E. globulus allows F1 hybrids to act as a bridge for bi-directional gene flow between these species. However, strong and persistent post-zygotic barriers exist, meaning that wherever F1 hybridization does occur, intrinsic post-zygotic barriers will be responsible for most reproductive isolation in this system. This potential transient nature of mechanical barriers to zygote formation due to additive inheritance in hybrids appears under-appreciated, and highlights the often important role that intrinsic post-mating barriers play in maintaining species boundaries at zones of secondary contact.

Promise and Challenge of DNA Barcoding in Venus Slipper (Paphiopedilum)

Orchidaceae are one of the largest families of flowering plants, with over 27,000 species described and all orchids are listed in CITES. Moreover, the seedlings of orchid species from the same genus are similar. The objective of DNA barcoding is rapid, accurate, and automated species identification, which may be used to identify illegally traded endangered species from vegetative specimens of Paphiopedilum (Venus slipper), a flagship group for plant conservation with high ornamental and commercial values. Here, we selected eight chloroplast barcodes and nrITS to evaluate their suitability in Venus slippers. The results indicate that all tested barcodes had no barcoding gap and the core plant barcodes showed low resolution for the identification of Venus slippers (18.86%). Of the single-locus barcodes, nrITS is the most efficient for the species identification of the genus (52.27%), whereas matK + atpF-atpH is the most efficient multi-locus combination (28.97%). Therefore, we recommend the combination of matK + atpF-atpH + ITS as a barcode for Venus slippers. Furthermore, there is an upper limit of resolution of the candidate barcodes, and only half of the taxa with multiple samples were identified successfully. The low efficiency of these candidate barcodes in Venus slippers may be caused by relatively recent speciation, the upper limit of the barcodes, and/or the sampling density. Although the discriminatory power is relatively low, DNA barcoding may be a promising tool to identify species involved in illegal trade, which has broad applications and is valuable for orchid conservation.

Comparative Transcriptome and Chloroplast Genome Analyses of Two Related Dipteronia Species

Dipteronia (order Sapindales) is an endangered genus endemic to China and has two living species, D.sinensis and D. dyeriana. The plants are closely related to the genus Acer, which is also classified in the order Sapindales. Evolutionary studies on Dipteronia have been hindered by the paucity of information on their genomes and plastids. Here, we used next generation sequencing to characterize the transcriptomes and complete chloroplast genomes of both Dipteronia species. A comparison of the transcriptomes of both species identified a total of 7814 orthologs. Estimation of selection pressures using Ka/Ks ratios showed that only 30 of 5435 orthologous pairs had a ratio significantly >1, i.e., showing positive selection. However, 4041 orthologs had a Ka/Ks < 0.5 (p < 0.05), suggesting that most genes had likely undergone purifying selection. Based on orthologous unigenes, 314 single copy nuclear genes (SCNGs) were identified. Through a combination of de novo and reference guided assembly, plastid genomes were obtained; that of D. sinensis was 157,080 bp and that of D. dyeriana was 157,071 bp. Both plastid genomes encoded 87 protein coding genes, 40 tRNAs, and 8 rRNAs; no significant differences were detected in the size, gene content, and organization of the two plastomes. We used the whole chloroplast genomes to determine the phylogeny of D. sinensis and D. dyeriana and confirmed that the two species were highly divergent. Overall, our study provides comprehensive transcriptomic and chloroplast genomic resources, which will be valuable for future evolutionary studies of Dipteronia.