Genome size diversity in orchids: consequences and evolution

Genome-Wide Identification of PEBP Gene Family in Two Dendrobium Species and Expression Patterns in Dendrobium chrysotoxum

The PEBP gene family plays a significant role in regulating flower development and formation. To understand its function in Dendrobium chrysotoxum and D. nobile flowering, we identified 22 PEBP genes (11 DchPEBPs and 11 DnoPEBPs) from both species. We conducted analyses on their conserved domains and motifs, phylogenetic relationships, chromosome distribution, collinear correlation, and cis elements. The classification results showed that the 22 PEBPs were mainly divided into three clades, as follows: FT, MFT, and TFL1. A sequence analysis showed that most PEBP proteins contained five conserved domains, while a gene structure analysis revealed that 77% of the total PEBP genes contained four exons and three introns. The promoter regions of the 22 PEBPs contained several cis elements related to hormone induction and light response. This suggests these PEBPs could play a role in regulating flower development by controlling photoperiod and hormone levels. Additionally, a collinearity analysis revealed three pairs of duplicate genes in the genomes of both D. chrysotoxum and D. nobile. Furthermore, RT-qPCR has found to influence the regulatory effect of DchPEBPs on the development of flower organs (sepals, petals, lip, ovary, and gynostemium) during the flowering process (bud, transparent stage, and initial bloom). The results obtained imply that DchPEBP8 and DchPEBP9 play a role in the initial bloom and that DchPEBP7 may inhibit flowering processes. Moreover, DchPEBP9 may potentially be involved in the development of reproductive functionality. PEBPs have regulatory functions that modulate flowering. FT initiates plant flowering by mediating photoperiod and temperature signals, while TFL1 inhibits flowering processes. These findings provide clues for future studies on flower development in Dendrobium.

Does the evolution of micromorphology accompany chromosomal changes on dysploid and polyploid levels in the Barnardia japonica complex (Hyacinthaceae)?

BACKGROUND

Chromosome number and genome size changes via dysploidy and polyploidy accompany plant diversification and speciation. Such changes often impact also morphological characters. An excellent system to address the questions of how extensive and structured chromosomal changes within one species complex affect the phenotype is the monocot species complex of Barnardia japonica. This taxon contains two well established and distinct diploid cytotypes differing in base chromosome numbers (AA: x = 8, BB: x = 9) and their allopolyploid derivatives on several ploidy levels (from 3x to 6x). This extensive and structured genomic variation, however, is not mirrored by gross morphological differentiation.

RESULTS

The current study aims to analyze the correlations between the changes of chromosome numbers and genome sizes with palynological and leaf micromorphological characters in diploids and selected allopolyploids of the B. japonica complex. The chromosome numbers varied from 2n = 16 and 18 (2n = 25 with the presence of supernumerary B chromosomes), and from 2n = 26 to 51 in polyploids on four different ploidy levels (3x, 4x, 5x, and 6x). Despite additive chromosome numbers compared to diploid parental cytotypes, all polyploid cytotypes have experienced genome downsizing. Analyses of leaf micromorphological characters did not reveal any diagnostic traits that could be specifically assigned to individual cytotypes. The variation of pollen grain sizes correlated positively with ploidy levels.

CONCLUSIONS

This study clearly demonstrates that karyotype and genome size differentiation does not have to be correlated with morphological differentiation of cytotypes.

Stoichiometry of carbon, nitrogen and phosphorus is closely linked to trophic modes in orchids

BACKGROUND

Mycorrhiza is a ubiquitous form of symbiosis based on the mutual, beneficial exchange of resources between roots of autotrophic (AT) plants and heterotrophic soil fungi throughout a complex network of fungal mycelium. Mycoheterotrophic (MH) and mixotrophic (MX) plants can parasitise this system, gaining all or some (respectively) required nutrients without known reciprocity to the fungus. We applied, for the first time, an ecological stoichiometry framework to test whether trophic mode of plants influences their elemental carbon (C), nitrogen (N), and phosphorus (P) composition and may provide clues about their biology and evolution within the framework of mycorrhizal network functioning.

RESULTS

We analysed C:N:P stoichiometry of 24 temperate orchid species and P concentration of 135 species from 45 plant families sampled throughout temperate and intertropical zones representing the three trophic modes (AT, MX and MH). Welch's one-way ANOVA and PERMANOVA were used to compare mean nutrient values and their proportions among trophic modes, phylogeny, and climate zones. Nutrient concentration and stoichiometry significantly differentiate trophic modes in orchids. Mean foliar C:N:P stoichiometry showed a gradual increase of N and P concentration and a decrease of C: nutrients ratio along the trophic gradient AT < MX < MH, with surprisingly high P requirements of MH orchids. Although P concentration in orchids showed the trophy-dependent pattern regardless of climatic zone, P concentration was not a universal indicator of trophic modes, as shown by ericaceous MH and MX plants.

CONCLUSION

The results imply that there are different evolutionary pathways of adaptation to mycoheterotrophic nutrient acquisition, and that the high nutrient requirements of MH orchids compared to MH plants from other families may represent a higher cost to the fungal partner and consequently lead to the high fungal specificity observed in MH orchids.

Bioinformatics Analysis of MSH1 Genes of Green Plants: Multiple Parallel Length Expansions, Intron Gains and Losses, Partial Gene Duplications, and Alternative Splicing

MutS homolog 1 (MSH1) is involved in the recombining and repairing of organelle genomes and is essential for maintaining their stability. Previous studies indicated that the length of the gene varied greatly among species and detected species-specific partial gene duplications in Physcomitrella patens. However, there are critical gaps in the understanding of the gene size expansion, and the extent of the partial gene duplication of MSH1 remains unclear. Here, we screened MSH1 genes in 85 selected species with genome sequences representing the main clades of green plants (Viridiplantae). We identified the MSH1 gene in all lineages of green plants, except for nine incomplete species, for bioinformatics analysis. The gene is a singleton gene in most of the selected species with conserved amino acids and protein domains. Gene length varies greatly among the species, ranging from 3234 bp in Ostreococcus tauri to 805,861 bp in Cycas panzhihuaensis. The expansion of MSH1 repeatedly occurred in multiple clades, especially in Gymnosperms, Orchidaceae, and Chloranthus spicatus. MSH1 has exceptionally long introns in certain species due to the gene length expansion, and the longest intron even reaches 101,025 bp. And the gene length is positively correlated with the proportion of the transposable elements (TEs) in the introns. In addition, gene structure analysis indicated that the MSH1 of green plants had undergone parallel intron gains and losses in all major lineages. However, the intron number of seed plants (gymnosperm and angiosperm) is relatively stable. All the selected gymnosperms contain 22 introns except for Gnetum montanum and Welwitschia mirabilis, while all the selected angiosperm species preserve 21 introns except for the ANA grade. Notably, the coding region of MSH1 in algae presents an exceptionally high GC content (47.7% to 75.5%). Moreover, over one-third of the selected species contain species-specific partial gene duplications of MSH1, except for the conserved mosses-specific partial gene duplication. Additionally, we found conserved alternatively spliced MSH1 transcripts in five species. The study of MSH1 sheds light on the evolution of the long genes of green plants.

Identification of three types of O-glycosylated flavonoids in Dendrobium loddigesii, Dendrobium primulinum, Dendrobium crepidatum, Dendrobium porphyrochilum, and Dendrobium hancockii using mass spectrometry

RATIONALE

Flavonoids, representing the pharmacologically active ingredients, are found widely in Dendrobium species. The biodiversity of Dendrobium makes the identification of its varieties all the more complicated. Previous studies showed that C-glycosylated flavones and a few O-glycosylated flavonols could be used in the identification of various Dendrobium species. Accordingly, this study further explores the significance of the identification of various types of O-glycosylated flavonoids in Dendrobium species.

METHODS

High-performance liquid chromatography coupled with electrospray ionization multistage tandem mass spectrometry (HPLC-ESI-MSⁿ ) was used to identify the chemical constituents in five types of Dendrobium: Dendrobium loddigesii, Dendrobium primulinum, Dendrobium crepidatum, Dendrobium porphyrochilum, and Dendrobium hancockii.

RESULTS

A total of 41 O-glycosylated flavonoids and 3 C-glycosylated flavones were identified, among which O-glycosylated dihydroflavones were the main flavonoids in D. loddigesii and D. primulinum, O-Glycosylated flavonols were rich in both D. crepidatum and D. porphyrochilum characterized by the main aglycone, substituted sugars, and their structural characteristics, and O-glycosylated flavones were the main constituents in D. hancockii.

CONCLUSIONS

In this study, three types of O-glycosylated flavonoids in the five Dendrobium species were determined to have certain significance. This also provides a reference for the identification of other O-glycosylated flavonoids in Chinese herbs.

Characterization of repetitive sequences in Dendrobium officinale and comparative chromosomal structures in Dendrobium species using FISH

Tandem repeats are one of the most conserved features in the eukaryote genomes. Dendrobium is the third largest genus in family Orchidaceae compromising over 1,200 species. However, the organization of repetitive sequences in Dendrobium species remains unclear. In this study, we performed the identification and characterization of the tandem repeats in D. officinale genome using graph-based clustering and Fluorescence in situ hybridization (FISH). Six major clusters including five satellite DNAs (DofSat1-5) and one 5S rDNA repeat (Dof5S) were identified as tandem repeats. The tandem organization of DofSat5 was verified by PCR amplification and southern blotting. The chromosomal locations of the repetitive DNAs in D. officinale were investigated by FISH using the tandem repeats and oligos probes. The results showed that each of the DofSat5, 5S and 45S rDNA had one pair of strong signals on D. officinale chromosomes. The distribution of repetitive DNAs along chromosomes was also investigated based on genomic in situ hybridization (GISH) among four Dendrobium species. The results suggested complex chromosomal fusion/segmentation and rearrangements during the evolution of Dendrobium species. In conclusion, the present study provides new landmarks for unequival differentiation of the Dendrobium chromosomes and facilitate the understanding the chromosome evolution in Dendrobium speceis.

Genome-wide identification and adaptive evolution of CesA/Csl superfamily among species with different life forms in Orchidaceae

Orchidaceae, with more than 25,000 species, is one of the largest flowering plant families that can successfully colonize wide ecological niches, such as land, trees, or rocks, and its members are divided into epiphytic, terrestrial, and saprophytic types according to their life forms. Cellulose synthase (CesA) and cellulose synthase-like (Csl) genes are key regulators in the synthesis of plant cell wall polysaccharides, which play an important role in the adaptation of orchids to resist abiotic stresses, such as drought and cold. In this study, nine whole-genome sequenced orchid species with three types of life forms were selected; the CesA/Csl gene family was identified; the evolutionary roles and expression patterns of CesA/Csl genes adapted to different life forms and abiotic stresses were investigated. The CesA/Csl genes of nine orchid species were divided into eight subfamilies: CesA and CslA/B/C/D/E/G/H, among which the CslD subfamily had the highest number of genes, followed by CesA, whereas CslB subfamily had the least number of genes. Expansion of the CesA/Csl gene family in orchids mainly occurred in the CslD and CslF subfamilies. Conserved domain analysis revealed that eight subfamilies were conserved with variations in orchids. In total, 17 pairs of CesA/Csl homologous genes underwent positive selection, of which 86%, 14%, and none belonged to the epiphytic, terrestrial, and saprophytic orchids, respectively. The inter-species collinearity analysis showed that the CslD genes expanded in epiphytic orchids. Compared with terrestrial and saprophytic orchids, epiphytic orchids experienced greater strength of positive selection, with expansion events mostly related to the CslD subfamily, which might have resulted in strong adaptability to stress in epiphytes. Experiments on stem expression changes under abiotic stress showed that the CslA might be a key subfamily in response to drought stress for orchids with different life forms, whereas the CslD might be a key subfamily in epiphytic and saprophytic orchids to adapt to freezing stress. This study provides the basic knowledge for the further systematic study of the adaptive evolution of the CesA/Csl superfamily in angiosperms with different life forms, and research on orchid-specific functional genes related to life-history trait evolution.

Are chromosome number and genome size associated with habit and environmental niche variables? Insights from the Neotropical orchids

BACKGROUND AND AIMS

The entangled relationship of chromosome number and genome size with species distribution has been the subject of study for almost a century, but remains an open question due to previous ecological and phylogenetic knowledge constraints. To better address this subject, we used the clade Maxillariinae, a widely distributed and karyotypically known orchid group, as a model system to infer such relationships in a robust methodological framework.

METHODS

Based on the literature and new data, we gathered the chromosome number and genome size for 93 and 64 species, respectively. We built a phylogenetic hypothesis and assessed the best macroevolutionary model for both genomic traits. Additionally, we collected together ecological data (preferences for bioclimatic variables, elevation and habit) used as explanatory variables in multivariate phylogenetic models explaining genomic traits. Finally, the impact of polyploidy was estimated by running the analyses with and without polyploids in the sample.

KEY RESULTS

The association between genomic and ecological data varied depending on whether polyploids were considered or not. Without polyploids, chromosome number failed to present consistent associations with ecological variables. With polyploids, there was a tendency to waive epiphytism and colonize new habitats outside humid forests. The genome size showed association with ecological variables: without polyploids, genome increase was associated with flexible habits, with higher elevation and with drier summers; with polyploids, genome size increase was associated with colonizing drier environments.

CONCLUSIONS

The chromosome number and genome size variations, essential but neglected traits in the ecological niche, are shaped in the Maxillariinae by both neutral and adaptive evolution. Both genomic traits are partially correlated to bioclimatic variables and elevation, even when controlling for phylogenetic constraints. While polyploidy was associated with shifts in the environmental niche, the genome size emerges as a central trait in orchid evolution by the association between small genome size and epiphytism, a key innovation to Neotropical orchid diversification.

Chromosome Number and Genome Size Evolution in Brasolia and Sobralia (Sobralieae, Orchidaceae)

Despite the clear circumscription of tribe Sobralieae (Orchidaceae), its internal relationships are still dubious. The recently delimited genus Brasolia, based on previous Sobralia species, is now assumed to be paraphyletic, with a third genus, Elleanthus, nested in it. The morphology of these three genera is significantly different, indicating the necessity of new data for a better genera delimitation. Though morphology and molecular data are available, cytogenetics data for Sobralieae is restricted to two Sobralia and one Elleanthus species. Aiming to evaluate the potential of cytogenetic data for Brasolia-Elleanthus-Sobralia genera delimitation, we present chromosome number and genome size data for 21 and 20 species, respectively, and used such data to infer the pattern of karyotype evolution in these genera. The analysis allowed us to infer x = 24 as the base chromosome number and genome size of average 1C-value of 5.0 pg for the common ancestor of Brasolia-Elleanthus-Sobralia. The recurrent descending dysploidy in Sobralieae and the punctual genome upsize suggest a recent diversification in Sobralieae but did not allow differing between Brasolia and Sobralia. However, the basal position of tribe Sobralieae in the subfamily Epidendroideae makes this tribe of interest to further studies clarifying the internal delimitation and pattern of karyotype evolution.

Genome-wide identification of Cymbidium sinense WRKY gene family and the importance of its Group III members in response to abiotic stress

Transcription factors (TFs) of the WRKY family play pivotal roles in defense responses and secondary metabolism of plants. Although WRKY TFs are well documented in numerous plant species, no study has performed a genome-wide investigation of the WRKY gene family in Cymbidium sinense. In the present work, we found 64 C. sinense WRKY (CsWRKY) TFs, and they were further divided into eight subgroups. Chromosomal distribution of CsWRKYs revealed that the majority of these genes were localized on 16 chromosomes, especially on Chromosome 2. Syntenic analysis implied that 13 (20.31%) genes were derived from segmental duplication events, and 17 orthologous gene pairs were identified between Arabidopsis thaliana WRKY (AtWRKY) and CsWRKY genes. Moreover, 55 of the 64 CsWRKYs were detectable in different plant tissues in response to exposure to plant hormones. Among them, Group III members were strongly induced in response to various hormone treatments, indicating their potential essential roles in hormone signaling. We subsequently analyzed the function of CsWRKY18 in Group III. The CsWRKY18 was localized in the nucleus. The constitutive expression of CsWRKY18 in Arabidopsis led to enhanced sensitivity to ABA-mediated seed germination and root growth and elevated plant tolerance to abiotic stress within the ABA-dependent pathway. Overall, our study represented the first genome-wide characterization and functional analysis of WRKY TFs in C. sinense, which could provide useful clues about the evolution and functional description of CsWRKY genes.

Cytological, Biochemical, and Transcriptomic Analyses of a Novel Yellow Leaf Variation in a Paphiopedilum (Orchidaceae) SCBG COP15

The genus Paphiopedilum, belonging to the Orchidaceae, has high ornamental value. Leaf variations can considerably improve the economic and horticultural value of the orchids. In the study, a yellow leaf mutant of a Paphiopedilum hybrid named P. SCBG COP15 was identified during the in vitro plant culture process; however, little is known about their molecular mechanisms. For this, RNA-seq libraries were created and used for the transcriptomic profiling of P. SCBG COP15 and the yellow mutant. The Chl a, Chl b, and carotenoid contents in the yellow leaves decreased by approximately 75.99%, 76.92%, and 56.83%, respectively, relative to the green leaves. Decreased chloroplasts per cell and abnormal chloroplast ultrastructure were observed by electron microscopic investigation in yellowing leaves; photosynthetic characteristics and Chl fluorescence parameters were also decreased in the mutant. Altogether, 34,492 unigenes were annotated by BLASTX; 1,835 DEGs were identified, consisting of 697 upregulated and 1138 downregulated DEGs. HEMA, CRD, CAO, and CHLE, involved in Chl biosynthesis, were predicted to be key genes responsible for leaf yellow coloration. Our findings provide an essential genetic resource for understanding the molecular mechanism of leaf color variation and breeding new varieties of Paphiopedilum with increased horticultural value.

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.

Frequent ploidy changes in Salicaceae indicates widespread sharing of the salicoid whole genome duplication by the relatives of Populus L. and Salix L

BACKGROUNDS

Populus and Salix belong to Salicaceae and are used as models to investigate woody plant physiology. The variation of karyotype and nuclear DNA content can partly reflect the evolutionary history of the whole genome, and can provide critical information for understanding, predicting, and potentially ameliorating the woody plant traits. Therefore, it is essential to study the chromosome number (CN) and genome size in detail to provide information for revealing the evolutionary process of Salicaceae.

RESULTS

In this study, we report the somatic CNs of seventeen species from eight genera in Salicaceae. Of these, CNs for twelve species and for five genera are reported for the first time. Among the three subfamilies of Salicaceae, the available data indicate CN in Samydoideae is n = 21, 22, 42. The only two genera, Dianyuea and Scyphostegia, in Scyphostegioideae respectively have n = 9 and 18. In Salicoideae, Populus, Salix and five genera closely related to them (Bennettiodendron, Idesia, Carrierea, Poliothyrsis, Itoa) are based on relatively high CNs from n = 19, 20, 21, 22 to n = 95 in Salix. However, the other genera of Salicoideae are mainly based on relatively low CNs of n = 9, 10, 11. The genome sizes of 35 taxa belonging to 14 genera of Salicaceae were estimated. Of these, the genome sizes of 12 genera and all taxa except Populus euphratica are first reported. Except for Dianyuea, Idesia and Bennettiodendron, all examined species have relatively small genome sizes of less than 1 pg, although polyploidization exists.

CONCLUSIONS

The variation of CN and genome size across Salicaceae indicates frequent ploidy changes and a widespread sharing of the salicoid whole genome duplication (WGD) by the relatives of Populus and Salix. The shrinkage of genome size after WGD indicates massive loss of genomic components. The phylogenetic asymmetry in clade of Populus, Salix, and their close relatives suggests that there is a lag-time for the subsequent radiations after the salicoid WGD event. Our results provide useful data for studying the evolutionary events of Salicaceae.

A chromosome-scale genome assembly and annotation of the spring orchid (Cymbidium goeringii)

Cymbidium goeringii, commonly known as the spring orchid, has long been favoured for horticultural purposes in Asian countries. It is a popular orchid with much demand for improvement and development for its valuable varieties. Until now, its reference genome has not been published despite its popularity and conservation efforts. Here, we report the de novo assembly of the C. goeringii genome, which is the largest among the orchids published to date, using a strategy that combines short- and long-read sequencing and chromosome conformation capture (Hi-C) information. The total length of all scaffolds is 3.99 Gb, with an N50 scaffold size of 178.2 Mb. A total of 29,556 protein-coding genes were annotated and 3.55 Gb (88.87% of genome) repetitive sequences were identified. We constructed pseudomolecular chromosomes using Hi-C, incorporating 89.4% of the scaffolds in 20 chromosomes. We identified 220 expanded and 106 contracted genes families in C. goeringii after divergence from its close relative. We also identified new gene families, resistance gene analogues and changes within the MADS-box genes, which control a diverse set of developmental processes during orchid evolution. Our high quality chromosomal-level assembly of C. goeringii can provide a platform for elucidating the genomic evolution of orchids, mining functional genes for agronomic traits and for developing molecular markers for accelerated breeding as well as accelerating conservation efforts.

Chromosome-scale assembly of the Dendrobium chrysotoxum genome enhances the understanding of orchid evolution

As one of the largest families of angiosperms, the Orchidaceae family is diverse. Dendrobium represents the second largest genus of the Orchidaceae. However, an assembled high-quality genome of species in this genus is lacking. Here, we report a chromosome-scale reference genome of Dendrobium chrysotoxum, an important ornamental and medicinal orchid species. The assembled genome size of D. chrysotoxum was 1.37 Gb, with a contig N50 value of 1.54 Mb. Of the sequences, 95.75% were anchored to 19 pseudochromosomes. There were 30,044 genes predicted in the D. chrysotoxum genome. Two whole-genome polyploidization events occurred in D. chrysotoxum. In terms of the second event, whole-genome duplication (WGD) was also found to have occurred in other Orchidaceae members, which diverged mainly via gene loss immediately after the WGD event occurred; the first duplication was found to have occurred in most monocots (tau event). We identified sugar transporter (SWEET) gene family expansion, which might be related to the abundant medicinal compounds and fleshy stems of D. chrysotoxum. MADS-box genes were identified in D. chrysotoxum, as well as members of TPS and Hsp90 gene families, which are associated with resistance, which may contribute to the adaptive evolution of orchids. We also investigated the interplay among carotenoid, ABA, and ethylene biosynthesis in D. chrysotoxum to elucidate the regulatory mechanisms of the short flowering period of orchids with yellow flowers. The reference D. chrysotoxum genome will provide important insights for further research on medicinal active ingredients and breeding and enhances the understanding of orchid evolution.

Genome Size Diversity in Rare, Endangered, and Protected Orchids in Poland

Orchidaceae is one of the largest and the most widespread plant families with many species threatened with extinction. However, only about 1.5% of orchids’ genome sizes have been known so far. The aim of this study was to estimate the genome size of 15 species and one infraspecific taxon of endangered and protected orchids growing wild in Poland to assess their variability and develop additional criterion useful in orchid species identification and characterization. Flow cytometric genome size estimation revealed that investigated orchid species possessed intermediate, large, and very large genomes. The smallest 2C DNA content possessed Liparis loeselii (14.15 pg), while the largest Cypripedium calceolus (82.10 pg). It was confirmed that the genome size is characteristic to the subfamily. Additionally, for four species Epipactis albensis, Ophrys insectifera, Orchis mascula, Orchis militaris and one infraspecific taxon, Epipactis purpurata f. chlorophylla the 2C DNA content has been estimated for the first time. Genome size estimation by flow cytometry proved to be a useful auxiliary method for quick orchid species identification and characterization.

Rapid radiation of angraecoids (Orchidaceae, Angraecinae) in tropical Africa characterised by multiple karyotypic shifts under major environmental instability

Angraecoid orchids present a remarkable diversity of chromosome numbers, which makes them a highly suitable system for exploring the impact of karyotypic changes on cladogenesis, diversification and morphological differentiation. We compiled an annotated cytotaxonomic checklist for 126 species of Angraecinae, which was utilised to reconstruct chromosomal evolution using a newly-produced, near-comprehensive phylogenetic tree that includes 245 angraecoid taxa. In tandem with this improved phylogenetic framework, using combined Bayesian, maximum likelihood and parsimony approaches on ITS-1 and five plastid markers, we propose a new cladistic nomenclature for the angraecoids, and we estimate a new timeframe for angraecoid radiation based on a secondary calibration, and calculate diversification rates using a Bayesian approach. Coincident divergence dates between clades with identical geographical distributions in the angraecoids and the pantropical orchid genus Bulbophyllum suggest that the same events may have intervened in the dispersal of these two epiphytic groups between Asia, continental Africa, Madagascar and the Neotropics. The major angraecoid lineages probably began to differentiate in the Middle Miocene, and most genera and species emerged respectively around the Late Miocene-Pliocene boundary and the Pleistocene. Ancestral state reconstruction using maximum likelihood estimation revealed an eventful karyotypic history dominated by descending dysploidy. Karyotypic shifts seem to have paralleled cladogenesis in continental tropical Africa, where approximately 90% of the species have descended from at least one inferred transition from n = 17-18 to n = 25 during the Middle Miocene Climatic Transition, followed by some clade-specific descending and ascending dysploidy from the Late Miocene to the Pleistocene. Conversely, detected polyploidy is restricted to a few species lineages mostly originating during the Pleistocene. No increases in net diversification could be related to chromosome number changes, and the apparent net diversification was found to be highest in Madagascar, where karyotypic stasis predominates. Finally, shifts in chromosome number appear to have paralleled the evolution of rostellum structure, leaflessness, and conspicuous changes in floral colour.

Effective double-digest RAD sequencing and genotyping despite large genome size

Obtaining informative data is the ambition of any genomic project, but in nonmodel species with very large genomes, pursuing such a goal requires surmounting a series of analytical challenges. Double-digest RAD sequencing is routinely used in nonmodel organisms and offers some control over the volume of data obtained. However, the volume of data recovered is not always an indication of the reliability of data sets, and quality checks are necessary to ensure that true and artefactual information is set apart. In the present study, we aim to fill the gap existing between the known applicability of RAD sequencing methods in plants with large genomes and the use of the retrieved loci for population genetic inference. By analysing two populations of Cypripedium calceolus, a nonmodel orchid species with a large genome size (1C ~ 31.6 Gbp), we provide a complete workflow from library preparation to bioinformatic filtering and inference of genetic diversity and differentiation. We show how filtering strategies to dismiss potentially misleading data need to be explored and adapted to data set-specific features. Moreover, we suggest that the occurrence of organellar sequences in libraries should not be neglected when planning the experiment and analysing the results. Finally, we explain how, in the absence of prior information about the genome of the species, seeking high standards of quality during library preparation and sequencing can provide an insurance against unpredicted technical or biological constraints.

Chloroplast Genomes of Two Species of Cypripedium: Expanded Genome Size and Proliferation of AT-Biased Repeat Sequences

The size of the chloroplast genome (plastome) of autotrophic angiosperms is generally conserved. However, the chloroplast genomes of some lineages are greatly expanded, which may render assembling these genomes from short read sequencing data more challenging. Here, we present the sequencing, assembly, and annotation of the chloroplast genomes of Cypripedium tibeticum and Cypripedium subtropicum. We de novo assembled the chloroplast genomes of the two species with a combination of short-read Illumina data and long-read PacBio data. The plastomes of the two species are characterized by expanded genome size, proliferated AT-rich repeat sequences, low GC content and gene density, as well as low substitution rates of the coding genes. The plastomes of C. tibeticum (197,815 bp) and C. subtropicum (212,668 bp) are substantially larger than those of the three species sequenced in previous studies. The plastome of C. subtropicum is the longest one of Orchidaceae to date. Despite the increase in genome size, the gene order and gene number of the plastomes are conserved, with the exception of an ∼75 kb large inversion in the large single copy (LSC) region shared by the two species. The most striking is the record-setting low GC content in C. subtropicum (28.2%). Moreover, the plastome expansion of the two species is strongly correlated with the proliferation of AT-biased non-coding regions: the non-coding content of C. subtropicum is in excess of 57%. The genus provides a typical example of plastome expansion induced by the expansion of non-coding regions. Considering the pros and cons of different sequencing technologies, we recommend hybrid assembly based on long and short reads applied to the sequencing of plastomes with AT-biased base composition.

The Application of Flow Cytometry for Estimating Genome Size, Ploidy Level Endopolyploidy, and Reproductive Modes in Plants

Over the years, the amount of DNA in a nucleus (genome size) has been estimated using a variety of methods, but increasingly, flow cytometry (FCM) has become the method of choice. The popularity of this technique lies in the ease of sample preparation and in the large number of particles (i.e., nuclei) that can be analyzed in a very short period of time. This chapter presents a step-by-step guide to estimating the nuclear DNA content of plant nuclei using FCM. Attempting to serve as a tool for daily laboratory practice, we list, in detail, the equipment required, specific reagents and buffers needed, as well as the most frequently used protocols to carry out nuclei isolation. In addition, solutions to the most common problems that users may encounter when working with plant material and troubleshooting advice are provided. Finally, information about the correct terminology to use and the importance of obtaining chromosome counts to avoid cytological misinterpretations of the FCM data are discussed.

Morphological and Genetic Diversities of Habenaria radiata (Orchidaceae) in the Kinki Area, Japan

Floral organs have evolved from leaves for reproduction, and the morphological analyses help to understand the plant diversity and evolution. Habenaria radiata (syn. Pecteilis radiata) is a terrestrial orchid living in wetlands in Japan, Russia, South Korea, and China. The habitats of this plant in Japan have been reduced because of environmental destruction and overexploitation, and thus it is on the Red List of Japan as a Near Threatened species. One of the three petals of the H. radiata flower is called a lip or labellum, which resembles a flying white bird, egret, or white heron, with its proposed function being to attract pollinators. To understand the diversity of H. radiata plants in different areas, we examined the lip morphology and phylogeny of populations from eight habitats in the Kinki area, Japan. The complex shapes of the lips were quantified and presented as a radar chart, enabling characterization of the morphological difference among populations. Phylogenetic analysis with microsatellite markers that we generated showed the variation of genetic diversity among populations, suggesting the different degrees of inbreeding, outbreeding, and vegetative propagation. Our approach offers a basic method to characterize the morphological and genetic diversity in natural populations.

Why are there so many bee-orchid species? Adaptive radiation by intra-specific competition for mnesic pollinators

Adaptive radiations occur mostly in response to environmental variation through the evolution of key innovations that allow emerging species to occupy new ecological niches. Such biological innovations may play a major role in niche divergence when emerging species are engaged in reciprocal ecological interactions. To demonstrate coevolution is a difficult task; only a few studies have confirmed coevolution as driver of speciation and diversification. Herein we review current knowledge about bee orchid (Ophrys spp.) reproductive biology. We propose that the adaptive radiation of the Mediterranean orchid genus Ophrys, comprising several hundred species, is due to coevolutionary dynamics between these plants and their pollinators. We suggest that pollination by sexual swindling used by Ophrys orchids is the main driver of this coevolution. Flowers of each Ophrys species mimic a sexually receptive female of one particular insect species, mainly bees. Male bees are first attracted by pseudo-pheromones emitted by Ophrys flowers that are similar to the sexual pheromones of their females. Males then are lured by the flower shape, colour and hairiness, and attempt to copulate with the flower, which glues pollen onto their bodies. Pollen is later transferred to the stigma of another flower of the same Ophrys species during similar copulation attempts. In contrast to rewarding pollination strategies, Ophrys pollinators appear to be parasitized. Here we propose that this apparent parasitism is in fact a coevolutionary relationship between Ophrys and their pollinators. For plants, pollination by sexual swindling could ensure pollination efficiency and specificity, and gene flow among populations. For pollinators, pollination by sexual swindling could allow habitat matching and inbreeding avoidance. Pollinators might use the pseudo-pheromones emitted by Ophrys to locate suitable habitats from a distance within complex landscapes. In small populations, male pollinators would disperse once they have memorized the local diversity of sexual pseudo-pheromone bouquets or if all Ophrys flowers are fertilized and thus repel pollinators via production of repulsive pheromones that mimic those produced by fertilized female bees. We propose the following evolutionary scenario: Ophrys radiation is driven by strong intra-specific competition among Ophrys individuals for the attraction of species-specific pollinators, which is a consequence of the high cognitive abilities of pollinators. Male bees record the pheromone signatures of kin or of previously courted partners to avoid further copulation attempts, thereby inducing strong selection on Ophrys for variation in odour bouquets emitted by individual flowers. The resulting odour bouquets could by chance correspond to pseudo-pheromones of the females of another bee species, and thus attract a new pollinator. If such pollinator shifts occur simultaneously in several indivuals, pollen exchanges might occur and initiate speciation. To reinforce the attraction of the new pollinator and secure prezygotic isolation, the following step is directional selection on flower phenotypes (shape, colour and hairiness) towards a better match with the body of the pollinator's female. Pollinator shift and the resulting prezygotic isolation is adaptive for new Ophrys species because they may benefit from competitor-free space for limited pollinators. We end our review by proritizing several critical research avenues.

Molecular Identification and Evaluation of the Genetic Diversity of Dendrobium Species Collected in Southern Vietnam

Dendrobium has been widely used not only as ornamental plants but also as food and medicines. The identification and evaluation of the genetic diversity of Dendrobium species support the conservation of genetic resources of endemic Dendrobium species. Uniquely identifying Dendrobium species used as medicines helps avoid misuse of medicinal herbs. However, it is challenging to identify Dendrobium species morphologically during their immature stage. Based on the DNA barcoding method, it is now possible to efficiently identify species in a shorter time. In this study, the genetic diversity of 76 Dendrobium samples from Southern Vietnam was investigated based on the ITS (Internal transcribed spacer), ITS2, matK (Maturase_K), rbcL (ribulose-bisphosphate carboxylase large subunit) and trnH-psbA (the internal space of the gene coding histidine transfer RNA (trnH) and gene coding protein D1, a polypeptide of the photosystem I reaction center (psaB)) regions. The ITS region was found to have the best identification potential. Nineteen out of 24 Dendrobium species were identified based on phylogenetic tree and Indel information of this region. Among these, seven identified species were used as medicinal herbs. The results of this research contributed to the conservation, propagation, and hybridization of indigenous Dendrobium species in Southern Vietnam.

Volatile Organic Compounds from Orchids: From Synthesis and Function to Gene Regulation

Orchids are one of the most significant plants that have ecologically adapted to every habitat on earth. Orchids show a high level of variation in their floral morphologies, which makes them popular as ornamental plants in the global market. Floral scent and color are key traits for many floricultural crops. Volatile organic compounds (VOCs) play vital roles in pollinator attraction, defense, and interaction with the environment. Recent progress in omics technology has led to the isolation of genes encoding candidate enzymes responsible for the biosynthesis and regulatory circuits of plant VOCs. Uncovering the biosynthetic pathways and regulatory mechanisms underlying the production of floral scents is necessary not only for a better understanding of the function of relevant genes but also for the generation of new cultivars with desirable traits through molecular breeding approaches. However, little is known about the pathways responsible for floral scents in orchids because of their long life cycle as well as the complex and large genome; only partial terpenoid pathways have been reported in orchids. Here, we review the biosynthesis and regulation of floral volatile compounds in orchids. In particular, we focused on the genes responsible for volatile compounds in various tissues and developmental stages in Cymbidium orchids. We also described the emission of orchid floral volatiles and their function in pollination ecology. Taken together, this review will provide a broad scope for the study of orchid floral scents.

Transcriptome analysis to identify candidate genes associated with the yellow-leaf phenotype of a Cymbidium mutant generated by γ-irradiation

Leaf color is an important agronomic trait in flowering plants, including orchids. However, factors underlying leaf phenotypes in plants remain largely unclear. A mutant displaying yellow leaves was obtained by the γ-ray-based mutagenesis of a Cymbidium orchid and characterized using RNA sequencing. A total of 144,918 unigenes obtained from over 25 million reads were assigned to 22 metabolic pathways in the Kyoto Encyclopedia of Genes and Genomes database. In addition, gene ontology was used to classify the predicted functions of transcripts into 73 functional groups. The RNA sequencing analysis identified 2,267 differentially expressed genes between wild-type and mutant Cymbidium sp. Genes involved in the chlorophyll biosynthesis and degradation, as well as ion transport, were identified and assayed for their expression levels in wild-type and mutant plants using quantitative real-time profiling. No critical expression changes were detected in genes involved in chlorophyll biosynthesis. In contrast, seven genes involved in ion transport, including two metal ion transporters, were down-regulated, and chlorophyllase 2, associated with chlorophyll degradation, was up-regulated. Together, these results suggest that alterations in chlorophyll metabolism and/or ion transport might contribute to leaf color in Cymbidium orchids.

Effects of the total dose and duration of γ-irradiation on the growth responses and induced SNPs of a Cymbidium hybrid

Purpose: Ionizing radiation has been used for developing new cultivars of diverse plant species, including Cymbidium orchid species. The effects of the total dose on mutation induction have been investigated; however, there is relatively little research on the influence of the dose rate or irradiation duration.Materials and methods: Thus, we analyzed the effects of the total dose and irradiation duration on the growth of Cymbidium hybrid RB001 protocorm-like bodies (PLBs). We completed a genotyping-by-sequencing analysis to compare the induced SNPs among five γ-irradiated populations with similar growth responses (LD₅₀) to γ-rays.Results: The optimal time to assess the effects of the γ-irradiation was at 6 months after the treatment. On the basis of the survival rate of γ-irradiated PLBs, the optimal doses (LD₅₀) for each irradiation duration were estimated: 1 h, 16.1 Gy; 4 h, 23.6 Gy; 8 h, 37.9 Gy; 16 h, 37.9 Gy; and 24 h, 40.0 Gy. The estimated optimal doses were duration-dependent at irradiation durations shorter than 8 h, but not at irradiation durations exceeding 8 h. A SNP comparison revealed a lack of significant differences among the mutations induced by γ-irradiations.Conclusions: These results indicate the irradiation duration affects PLB growth in response to γ-rays. Moreover, the mutations induced by a short-term treatment may be similar to those induced by a treatment over a longer period.

Diversity in genome size and GC content shows adaptive potential in orchids and is closely linked to partial endoreplication, plant life-history traits and climatic conditions

In angiosperms, genome size and nucleobase composition (GC content) exhibit pronounced variation with possible adaptive consequences. The hyperdiverse orchid family possessing the unique phenomenon of partial endoreplication (PE) provides a great opportunity to search for interactions of both genomic traits with the evolutionary history of the family. Using flow cytometry, we report values of both genomic traits and the type of endoreplication for 149 orchid species and compare these with a suite of life-history traits and climatic niche data using phylogeny-based statistics. The evolution of genomic traits was further studied using the Brownian motion (BM) and Ornstein-Uhlenbeck (OU) models to access their adaptive potential. Pronounced variation in genome size (341-54 878 Mb), and especially in GC content (23.9-50.5%), was detected among orchids. Diversity in both genomic traits was closely related to the type of endoreplication, plant growth form and climatic conditions. GC content was also associated with the type of dormancy. In all tested scenarios, OU models always outperformed BM models. Unparalleled GC content variation was discovered in orchids, setting new limits for plants. Our study indicates that diversity in both genome size and GC content has adaptive consequences and is tightly linked with evolutionary transitions to PE.

Iterative allogamy-autogamy transitions drive actual and incipient speciation during the ongoing evolutionary radiation within the orchid genus Epipactis (Orchidaceae)

BACKGROUND AND AIMS

The terrestrial orchid genus Epipactis has become a model system for the study of speciation via transitions from allogamy to autogamy, but close phylogenetic relationships have proven difficult to resolve through Sanger sequencing.

METHODS

We analysed with restriction site-associated sequencing (RAD-seq) 108 plants representing 29 named taxa that together span the genus, focusing on section Epipactis. Our filtered matrix of 12 543 single nucleotide polymorphisms was used to generate an unrooted network and a rooted, well-supported likelihood tree. We further inferred genetic structure through a co-ancestry heat map and admixture analysis, and estimated inbreeding coefficients per sample.

KEY RESULTS

The 27 named taxa of the ingroup were resolved as 11 genuine, geographically widespread species: four dominantly allogamous and seven dominantly autogamous. A single comparatively allogamous species, E. helleborine, is the direct ancestor of most of the remaining species, though one of the derived autogams has generated one further autogamous species. An assessment of shared ancestry suggested only sporadic hybridization between the re-circumscribed species. Taxa with the greatest inclination towards autogamy show less, if any, admixture, whereas the gene pools of more allogamous species contain a mixture alleles found in the autogams.

CONCLUSIONS

This clade is presently undergoing an evolutionary radiation driven by a wide spectrum of genotypic, phenotypic and environmental factors. Epipactis helleborine has also frequently generated many local variants showing inclinations toward autogamy (and occasionally cleistogamy), best viewed as incipient speciation from within the genetic background provided by E. helleborine, which thus becomes an example of a convincingly paraphyletic species. Autogams are often as widespread and ecologically successful as allogams.

RNA sequencing analysis of Cymbidium goeringii identifies floral scent biosynthesis related genes

BACKGROUND

Cymbidium goeringii belongs to the Orchidaceae, which is one of the most abundant angiosperm families. Cymbidium goeringii consist with high economic value and characteristics include fragrance and multiple flower colors. Floral scent is one of the important strategies for ensuring fertilization. However, limited genetic data is available in this non-model plant, and little known about the molecular mechanism responsible for floral scent in this orchid. Transcriptome and expression profiling data are needed to identify genes and better understand the biological mechanisms of floral scents in this species. Present transcriptomic data provides basic information on the genes and enzymes related to and pathways involved in flower secondary metabolism in this plant.

RESULTS

In this study, RNA sequencing analyses were performed to identify changes in gene expression and biological pathways related scent metabolism. Three cDNA libraries were obtained from three developmental floral stages: closed bud, half flowering stage and full flowering stage. Using Illumina technique 159,616,374 clean reads were obtained and were assembled into 85,868 final unigenes (average length 1194 nt), 33.85% of which were annotated in the NCBI non redundant protein database. Among this unigenes 36,082 were assigned to gene ontology and 23,164 were combined with COG groups. Total 33,417 unigenes were assigned in 127 pathways according to the Kyoto Encyclopedia of Genes and Genomes pathway database. According these transcriptomic data we identified number of candidates genes which differentially expressed in different developmental stages of flower related to fragrance biosynthesis. In q-RT-PCR most of the fragrance related genes highly expressed in half flowering stage.

CONCLUSIONS

RNA-seq and DEG data provided comprehensive gene expression information at the transcriptional level that could be facilitate the molecular mechanisms of floral biosynthesis pathways in three developmental phase's flowers in Cymbidium goeringii, moreover providing useful information for further analysis on C. goeringii, and other plants of genus Cymbidium.

Comparative transcriptomic analysis reveal the regulation mechanism underlying MeJA-induced accumulation of alkaloids in Dendrobium officinale

Dendrobium officinale is a traditional medicinal herb with a variety of bioactive components. Alkaloid is one of the major active ingredients of Dendrobium plants, and its immune regulatory effects have been well-studied. Although a number of genes involved in the biosynthetic pathway of alkaloids have been elucidated, the regulation mechanism underlying the methyl-jasmonate (MeJA)-induced accumulation of alkaloids in D. officinale is largely unknown. In our study, a total of 4,857 DEGs, including 2,943 up- and 1,932 down-regulated genes, were identified between the control and MeJA-treated groups. Kyoto Encyclopedia of Genes and Genomes annotation showed that a number of DEGs were associated with the putative alkaloid biosynthetic pathway in D. officinale. The main group of Dendrobium alkaloids are sesquiterpene alkaloids, which are the downstream products of mevalonate (MVA) and methylerythritol 4-phosphate (MEP) pathway. Several MVA and MEP pathway genes were significantly up-regulated by the MeJA treatment, suggesting an active precursor supply for the alkaloid biosynthesis under MeJA treatment. A number of MeJA-induced P450 family genes, aminotransferase genes and methyltransferase genes were identified, providing several important candidates to further elucidate the sesquiterpene alkaloid biosynthetic pathway of D. officinale. Furthermore, a large number of MeJA-induced transcript factor encoding genes were identified, suggesting a complex genetic network affecting the sesquiterpene alkaloid metabolism in D. officinale. Our data aids to reveal the regulation mechanism underlying the MeJA-induced accumulation of sesquiterpene alkaloids in D. officinale.

Comparative transcriptomic analysis reveal the regulation mechanism underlying MeJA-induced accumulation of alkaloids in Dendrobium officinale

Dendrobium officinale is a traditional medicinal herb with a variety of bioactive components. Alkaloid is one of the major active ingredients of Dendrobium plants, and its immune regulatory effects have been well-studied. Although a number of genes involved in the biosynthetic pathway of alkaloids have been elucidated, the regulation mechanism underlying the methyl-jasmonate (MeJA)-induced accumulation of alkaloids in D. officinale is largely unknown. In our study, a total of 4,857 DEGs, including 2,943 up- and 1,932 down-regulated genes, were identified between the control and MeJA-treated groups. Kyoto Encyclopedia of Genes and Genomes annotation showed that a number of DEGs were associated with the putative alkaloid biosynthetic pathway in D. officinale. The main group of Dendrobium alkaloids are sesquiterpene alkaloids, which are the downstream products of mevalonate (MVA) and methylerythritol 4-phosphate (MEP) pathway. Several MVA and MEP pathway genes were significantly up-regulated by the MeJA treatment, suggesting an active precursor supply for the alkaloid biosynthesis under MeJA treatment. A number of MeJA-induced P450 family genes, aminotransferase genes and methyltransferase genes were identified, providing several important candidates to further elucidate the sesquiterpene alkaloid biosynthetic pathway of D. officinale. Furthermore, a large number of MeJA-induced transcript factor encoding genes were identified, suggesting a complex genetic network affecting the sesquiterpene alkaloid metabolism in D. officinale. Our data aids to reveal the regulation mechanism underlying the MeJA-induced accumulation of sesquiterpene alkaloids in D. officinale.

Detection and quantification of salep with real time PCR utilizing the nr-ITS2 region

BACKGROUND

Salep is a highly valued food product, susceptible to adulteration with less costly additives such as starch. The aim of this study is to develop a real-time polymerase chain reaction (PCR) method for the detection and quantification of authentic salep using primers designed from the nr-ITS2 region.

RESULTS

The designed primers were shown to amplify selectively with salep DNA and cloning/sequencing of end-point PCR products obtained with authentic salep (Orchideceae) confirmed amplimers to be copies of the intended target region. Detection of salep DNA using the designed primer set was possible at levels as low as 20 pg. After initial evaluation of the real-time PCR method with binary DNA mixtures of salep and cornflour, quantification of salep was conducted in binary mixtures of potato starch (a more likely adulterant) and salep, showing that quantification within the range 25%-100% could be achieved. The real-time assay was performed using the 18S rDNA region co-amplified as an internal standard and evaluated using the ΔΔCt method. The calibration curves obtained as a result of real-time PCR runs showed linearity, with R² values above 99%.

CONCLUSION

A real-time PCR method using SybreGreen intercalating dye was established for salep detection and quantification based on the specific amplification of salep DNA. Method validation was performed utilizing blind and commercial samples. © 2018 Society of Chemical Industry.

Morphology and genome size of Epipactis helleborine (L.) Crantz (Orchidaceae) growing in anthropogenic and natural habitats

BACKGROUND

The process of apophytism or spreading native species to human-made habitats is one of the main elements in the creation of plant cover in anthropogenic areas. Lately, an increase of anthropogenic localities with valuable flora has been observed. Apophytes are also members of the family Orchidaceae, especially from the genus Epipactis. The aim of the study was to (i) determine and compare the phenotypic variation of E. helleborine (L.) Crantz plants in anthropogenic and natural habitats, (ii) compare the genome size of plants growing in natural and anthropogenic habitats. The results reported in this study may indicate that a habitat influences morphological characteristics of plant species.

METHODS

Field studies were conducted on four native stands and four stands in anthropogenic areas of E. helleborine in Poland in years 2011-2013. Biometrical analyses were performed on shoots and flowers. The flowers were characterised by 25 biometric features and measured using a Nikon SMZ 800 binocular, microscopic Moticam-1SP cameras and the MIPlus07 programme (Conbest Co.). The nuclear DNA content was determined in fresh and young leaves of E. helleborine, collected from four natural and four anthropogenic populations.

RESULTS

We observed that in anthropogenic populations: (i) shoots were higher than shoots from natural populations, (ii) flowers differed significantly in terms of ten biometric features between habitats, (iii) the genome size of some population differed significantly between plants growing in natural and anthropogenic habitats.

DISCUSSION

According to some researchers, the presence of phenotypic variability and the occurrence of ecotypes are adaptation strategies of plants to environmental changes. In our opinion, in the case of the studied anthropogenic habitats (roadside) in which the E. helleborine populations grew, we can talk about ecofen due to the often repeated set of characteristic features, i.e., high shoots, long inflorescence and long, broad leaves. We agree, however, that it is difficult to isolate a taxonomic unit for ecofen due to the lack of experimental research.

Development of polymorphic microsatellite markers by using de novo transcriptome assembly of Calanthe masuca and C. sinica (Orchidaceae)

BACKGROUND

Calanthe masuca and C. sinica are two genetically closely related species in Orchidaceae. C. masuca is widely distributed in Asia, whereas C. sinica is restricted to Yunnan and Guangxi Provinces in southwest China. Both play important roles in horticulture and are under the pressure of population decline. Understanding their genetic background can greatly help us develop effective conservation strategies for these species. Simple sequence repeats (SSRs) are useful for genetic diversity analysis, presumably providing key information for the study and preservation of the wild populations of the two species we are interested in.

RESULTS

In this study, we performed RNA-seq analysis on the leaves of C. masuca and C. sinica, obtaining 40,916 and 71,618 unigenes for each species, respectively. In total, 2,019/3,865 primer pairs were successfully designed from 3,764/7,189 putative SSRs, among which 197 polymorphic SSRs were screened out according to orthologous gene pairs. After mononucleotide exclusion, a subset of 129 SSR primers were analysed, and 13 of them were found to have high polymorphism levels. Further analysis demonstrated that they were feasible and effective against C. masuca and C. sinica as well as transferable to another species in Calanthe. Molecular evolutionary analysis revealed functional pathways commonly enriched in unigenes with similar evolutionary rates in the two species, as well as pathways specific to each species, implicating species-specific adaptation. The divergence time between the two closely related species was tentatively determined to be 3.42 ± 1.86 Mya.

CONCLUSIONS

We completed and analysed the transcriptomes of C. masuca and C. sinica, assembling large numbers of unigenes and generating effective polymorphic SSR markers. This is the first report of the development of expressed sequence tag (EST)-SSR markers for Calanthe. In addition, our study could enable further genetic diversity analysis and functional and comparative genomic studies on Calanthe.

The complete plastid genomes of Ophrys iricolor and O. sphegodes (Orchidaceae) and comparative analyses with other orchids

Sexually deceptive orchids of the genus Ophrys may rapidly evolve by adaptation to pollinators. However, understanding of the genetic basis of potential changes and patterns of relationships is hampered by a lack of genomic information. We report the complete plastid genome sequences of Ophrys iricolor and O. sphegodes, representing the two most species-rich lineages of the genus Ophrys. Both plastomes are circular DNA molecules (146754 bp for O. sphegodes and 150177 bp for O. iricolor) with the typical quadripartite structure of plastid genomes and within the average size of photosynthetic orchids. 213 Simple Sequence Repeats (SSRs) (31.5% polymorphic between O. iricolor and O. sphegodes) were identified, with homopolymers and dipolymers as the most common repeat types. SSRs were mainly located in intergenic regions but SSRs located in coding regions were also found, mainly in ycf1 and rpoC2 genes. The Ophrys plastome is predicted to encode 107 distinct genes, 17 of which are completely duplicated in the Inverted Repeat regions. 83 and 87 putative RNA editing sites were detected in 25 plastid genes of the two Ophrys species, all occurring in the first or second codon position. Comparing the rate of nonsynonymous (dN) and synonymous (dS) substitutions, 24 genes (including rbcL and ycf1) display signature consistent with positive selection. When compared with other members of the orchid family, the Ophrys plastome has a complete set of 11 functional ndh plastid genes, with the exception of O. sphegodes that has a truncated ndhF gene. Comparative analysis showed a large co-linearity with other related Orchidinae. However, in contrast to O. iricolor and other Orchidinae, O. sphegodes has a shift of the junction between the Inverted Repeat and Small Single Copy regions associated with the loss of the partial duplicated gene ycf1 and the truncation of the ndhF gene. Data on relative genomic coverage and validation by PCR indicate the presence, with a different ratio, of the two plastome types (i.e. with and without ndhF deletion) in both Ophrys species, with a predominance of the deleted type in O. sphegodes. A search for this deleted plastid region in O. sphegodes nuclear genome shows that the deleted region is inserted in a retrotransposon nuclear sequence. The present study provides useful genomic tools for studying conservation and patterns of relationships of this rapidly radiating orchid genus.

Evolutionary and functional potential of ploidy increase within individual plants: somatic ploidy mapping of the complex labellum of sexually deceptive bee orchids

Background and Aims

Recent tissue-level observations made indirectly via flow cytometry suggest that endoreplication (duplication of the nuclear genome within the nuclear envelope in the absence of subsequent cell division) is widespread within the plant kingdom. Here, we also directly observe ploidy variation among cells within individual petals, relating size of nucleus to cell micromorphology and (more speculatively) to function.

Methods

We compared the labella (specialized pollinator-attracting petals) of two European orchid genera: Dactylorhiza has a known predisposition to organismal polyploidy, whereas Ophrys exhibits exceptionally complex epidermal patterning that aids pseudocopulatory pollination. Confocal microscopy using multiple staining techniques allowed us to observe directly both the sizes and the internal structures of individual nuclei across each labellum, while flow cytometry was used to test for progressively partial endoreplication.

Key Results

In Dactylorhiza, endoreplication was comparatively infrequent, reached only low levels, and appeared randomly located across the labellum, whereas in Ophrys endoreplication was commonplace, being most frequent in large peripheral trichomes. Endoreplicated nuclei reflected both endomitosis and endocycling, the latter reaching the third round of genome doubling (16C) to generate polytene nuclei. All Ophrys individuals studied exhibited progressively partial endoreplication.

Conclusions

Comparison of the two genera failed to demonstrate the hypothesized pattern of frequent polyploid speciation in genera showing extensive endoreplication. Endoreplication in Ophrys appears more strongly positively correlated with cell size/complexity than with cell location or secretory role. Epigenetic control of gene overexpression by localized induction of endoreplication within individual plant organs may represent a significant component of a plant's developmental programme, contributing substantially to organ plasticity.

The Role of Transposable Elements in Speciation

Understanding the phenotypic and molecular mechanisms that contribute to genetic diversity between and within species is fundamental in studying the evolution of species. In particular, identifying the interspecific differences that lead to the reduction or even cessation of gene flow between nascent species is one of the main goals of speciation genetic research. Transposable elements (TEs) are DNA sequences with the ability to move within genomes. TEs are ubiquitous throughout eukaryotic genomes and have been shown to alter regulatory networks, gene expression, and to rearrange genomes as a result of their transposition. However, no systematic effort has evaluated the role of TEs in speciation. We compiled the evidence for TEs as potential causes of reproductive isolation across a diversity of taxa. We find that TEs are often associated with hybrid defects that might preclude the fusion between species, but that the involvement of TEs in other barriers to gene flow different from postzygotic isolation is still relatively unknown. Finally, we list a series of guides and research avenues to disentangle the effects of TEs on the origin of new species.

Phylotranscriptomic analysis and genome evolution of the Cypripedioideae (Orchidaceae)

PREMISE OF THE STUDY

The slipper orchids (Cypripedioideae) are a morphologically distinct subfamily of Orchidaceae. They also have some of the largest genomes in the orchids, which may be due to polyploidy or some other mechanism of genome evolution. We generated 10 transcriptomes and incorporated existing RNA-seq data to infer a multilocus nuclear phylogeny of the Cypripedioideae and to determine whether a whole-genome duplication event (WGD) correlated with the large genome size of this subfamily. Knowing more about timing of ancient polyploidy events can help us understand the evolution of one of the most species-rich plant families.

METHODS

Transcriptome data were used to identify low-copy orthologous genes to infer a phylogeny of Orchidaceae and to identify paralogs to place any WGD events on the species tree.

KEY RESULTS

Our transcriptome phylogeny confirmed relationships published in previous studies that used fewer markers but incorporated more taxa. We did not find a WGD event at the base of the slipper orchids; however, we did identify one on the Orchidaceae stem lineage. We also confirmed the presence of a previously identified WGD event deeper in the monocot phylogeny.

CONCLUSIONS

Although WGD has played a role in the evolution of Orchidaceae, polyploidy does not appear to be responsible for the large genome size of slipper orchids. The conserved set of 775 largely single-copy nuclear genes identified in this study should prove useful in future studies of orchid evolution.

Integrating restriction site-associated DNA sequencing (RAD-seq) with morphological cladistic analysis clarifies evolutionary relationships among major species groups of bee orchids

BACKGROUND AND AIMS

Bee orchids (Ophrys) have become the most popular model system for studying reproduction via insect-mediated pseudo-copulation and for exploring the consequent, putatively adaptive, evolutionary radiations. However, despite intensive past research, both the phylogenetic structure and species diversity within the genus remain highly contentious. Here, we integrate next-generation sequencing and morphological cladistic techniques to clarify the phylogeny of the genus.

METHODS

At least two accessions of each of the ten species groups previously circumscribed from large-scale cloned nuclear ribosomal internal transcibed spacer (nrITS) sequencing were subjected to restriction site-associated sequencing (RAD-seq). The resulting matrix of 4159 single nucleotide polymorphisms (SNPs) for 34 accessions was used to construct an unrooted network and a rooted maximum likelihood phylogeny. A parallel morphological cladistic matrix of 43 characters generated both polymorphic and non-polymorphic sets of parsimony trees before being mapped across the RAD-seq topology.

KEY RESULTS

RAD-seq data strongly support the monophyly of nine out of ten groups previously circumscribed using nrITS and resolve three major clades; in contrast, supposed microspecies are barely distinguishable. Strong incongruence separated the RAD-seq trees from both the morphological trees and traditional classifications; mapping of the morphological characters across the RAD-seq topology rendered them far more homoplastic.

CONCLUSIONS

The comparatively high level of morphological homoplasy reflects extensive convergence, whereas the derived placement of the fusca group is attributed to paedomorphic simplification. The phenotype of the most recent common ancestor of the extant lineages is inferred, but it post-dates the majority of the character-state changes that typify the genus. RAD-seq may represent the high-water mark of the contribution of molecular phylogenetics to understanding evolution within Ophrys; further progress will require large-scale population-level studies that integrate phenotypic and genotypic data in a cogent conceptual framework.

Recent progress in whole genome sequencing, high-density linkage maps, and genomic databases of ornamental plants

Genome information is useful for functional analysis of genes, comparative genomic analysis, breeding of new varieties by marker-assisted selection, and map-based gene isolation. Genome-related research in ornamentals plants has been relatively slow to develop because of their heterozygosity or polyploidy. Advances in analytical instruments, such as next-generation sequencers and information processing technologies have revolutionized biology, and have been applied in a large number and variety of species, including ornamental plants. Recently, high-quality whole genome sequences have been reported in plant genetics and physiology studies of model ornamentals, such as those in genus Petunia and Japanese morning glory (Ipomoea nil). In this review, whole genome sequencing and construction of high-density genetic linkage maps based on SNP markers of ornamentals will be discussed. The databases that store this information for ornamentals are also described.

Post genomics era for orchid research

Among 300,000 species in angiosperms, Orchidaceae containing 30,000 species is one of the largest families. Almost every habitats on earth have orchid plants successfully colonized, and it indicates that orchids are among the plants with significant ecological and evolutionary importance. So far, four orchid genomes have been sequenced, including Phalaenopsis equestris, Dendrobium catenatum, Dendrobium officinale, and Apostaceae shengen. Here, we review the current progress and the direction of orchid research in the post genomics era. These include the orchid genome evolution, genome mapping (genome-wide association analysis, genetic map, physical map), comparative genomics (especially receptor-like kinase and terpene synthase), secondary metabolomics, and genome editing.

The Apostasia genome and the evolution of orchids

WebComparing the whole genome sequence of Apostasia shenzhenica with transcriptome and genome data from five orchid subfamilies permits the reconstruction of an ancestral gene toolkit, providing insight into orchid origins, evolution and diversification.

Around 10 per cent of flowering plant species are orchids, with a broad diversity in both morphology and lifestyle. Apostasia is one of the earliest-diverging genera of Orchidaceae. To study the evolution and diversity of Orchidaceae, Zhong-Jian Liu, Yves Van de Peer and colleagues sequenced the genome of Apostasia shenzhenica, a self-pollinating species found in southeast China. The authors also report improved genomes for two species of Epidendroideae, Phalaenopsis equestris and Dendrobium catenatum, as well as transcriptome analysis of representatives of subfamilies of Orchidaceae. Their analyses provide insights into orchid origins, genome evolution, adaptation and diversification.

Constituting approximately 10% of flowering plant species, orchids (Orchidaceae) display unique flower morphologies, possess an extraordinary diversity in lifestyle, and have successfully colonized almost every habitat on Earth^1,2,3. Here we report the draft genome sequence of Apostasia shenzhenica⁴, a representative of one of two genera that form a sister lineage to the rest of the Orchidaceae, providing a reference for inferring the genome content and structure of the most recent common ancestor of all extant orchids and improving our understanding of their origins and evolution. In addition, we present transcriptome data for representatives of Vanilloideae, Cypripedioideae and Orchidoideae, and novel third-generation genome data for two species of Epidendroideae, covering all five orchid subfamilies. A. shenzhenica shows clear evidence of a whole-genome duplication, which is shared by all orchids and occurred shortly before their divergence. Comparisons between A. shenzhenica and other orchids and angiosperms also permitted the reconstruction of an ancestral orchid gene toolkit. We identify new gene families, gene family expansions and contractions, and changes within MADS-box gene classes, which control a diverse suite of developmental processes, during orchid evolution. This study sheds new light on the genetic mechanisms underpinning key orchid innovations, including the development of the labellum and gynostemium, pollinia, and seeds without endosperm, as well as the evolution of epiphytism; reveals relationships between the Orchidaceae subfamilies; and helps clarify the evolutionary history of orchids within the angiosperms.

Karyological patterns in the European endemic genus Soldanella L.: Absolute genome size variation uncorrelated with cytotype chromosome numbers

PREMISE OF THE STUDY

Detailed knowledge about the karyological diversity of organisms undoubtedly represents one of the crucial steps toward a better understanding of their evolutionary trends and history. We investigated the cytotype and absolute genome size (AGS) patterns in the European mountain-dwelling genus Soldanella (Primulaceae) in light of its geographic distribution and ecological diversification.

METHODS

Our chromosome number survey was based on 34 newly determined and 125 previously published chromosome counts. AGS was estimated on the basis of propidium iodide (PI) flow cytometry (299 individuals, 110 populations).

KEY RESULTS

We confirmed the existence of two cytotypes with the same ploidy level, i.e., euploid 2n = 40 and dysploid 2n = 38. The overall infrageneric AGS variation ranged between 2.97 and 3.99 pg (25.6% variation). The 2n = 40 cytotype harbors a modest amount of continuous AGS variation. With regard to its distribution area and ecology, the cytotype is ubiquitous. By contrast, the 2n = 38 cytotype was detected only in six forest-dwelling taxa with AGS variation segregated into three discrete, geographically separated groups. The AGS variation of the 2n = 38 cytotype was strongly correlated with elevation and longitude.

CONCLUSIONS

Despite the apparent morphological and ecological variation, members of the genus Soldanella have not undergone any pronounced cytotype and AGS diversification during their evolutionary history. The lack of correlation between chromosome numbers and AGS indicates that the evolutionary mechanism behind the origin of the dysploid cytotype 2n = 38 was a chromosomal fusion.

Genome Size Diversity in Lilium (Liliaceae) Is Correlated with Karyotype and Environmental Traits

Genome size (GS) diversity is of fundamental biological importance. The occurrence of giant genomes in angiosperms is restricted to just a few lineages in the analyzed genome size of plant species so far. It is still an open question whether GS diversity is shaped by neutral or natural selection. The genus Lilium, with giant genomes, is phylogenetically and horticulturally important and is distributed throughout the northern hemisphere. GS diversity in Lilium and the underlying evolutionary mechanisms are poorly understood. We performed a comprehensive study involving phylogenetically independent analysis on 71 species to explore the diversity and evolution of GS and its correlation with karyological and environmental traits within Lilium (including Nomocharis). The strong phylogenetic signal detected for GS in the genus provides evidence consistent with that the repetitive DNA may be the primary contributors to the GS diversity, while the significant positive relationships detected between GS and the haploid chromosome length (HCL) provide insights into patterns of genome evolution. The relationships between GS and karyotypes indicate that ancestral karyotypes of Lilium are likely to have exhibited small genomes, low diversity in centromeric index (CVCI) values and relatively high relative variation in chromosome length (CVCL) values. Significant relationships identified between GS and annual temperature and between GS and annual precipitation suggest that adaptation to habitat strongly influences GS diversity. We conclude that GS in Lilium is shaped by both neutral (genetic drift) and adaptive evolution. These findings will have important consequences for understanding the evolution of giant plant genomes, and exploring the role of repetitive DNA fraction and chromosome changes in a plant group with large genomes and conservation of chromosome number.

DNA Remodeling by Strict Partial Endoreplication in Orchids, an Original Process in the Plant Kingdom

DNA remodeling during endoreplication appears to be a strong developmental characteristic in orchids. In this study, we analyzed DNA content and nuclei in 41 species of orchids to further map the genome evolution in this plant family. We demonstrate that the DNA remodeling observed in 36 out of 41 orchids studied corresponds to strict partial endoreplication. Such process is developmentally regulated in each wild species studied. Cytometry data analyses allowed us to propose a model where nuclear states 2C, 4E, 8E, etc. form a series comprising a fixed proportion, the euploid genome 2C, plus 2-32 additional copies of a complementary part of the genome. The fixed proportion ranged from 89% of the genome in Vanilla mexicana down to 19% in V. pompona, the lowest value for all 148 orchids reported. Insterspecific hybridization did not suppress this phenomenon. Interestingly, this process was not observed in mass-produced epiphytes. Nucleolar volumes grow with the number of endocopies present, coherent with high transcription activity in endoreplicated nuclei. Our analyses suggest species-specific chromatin rearrangement. Towards understanding endoreplication, V. planifolia constitutes a tractable system for isolating the genomic sequences that confer an advantage via endoreplication from those that apparently suffice at diploid level.

Karyotype diversity and genome size variation in Neotropical Maxillariinae orchids

Orchidaceae is a widely distributed plant family with very diverse vegetative and floral morphology, and such variability is also reflected in their karyotypes. However, since only a low proportion of Orchidaceae has been analysed for chromosome data, greater diversity may await to be unveiled. Here we analyse both genome size (GS) and karyotype in two subtribes recently included in the broadened Maxillariinea to detect how much chromosome and GS variation there is in these groups and to evaluate which genome rearrangements are involved in the species evolution. To do so, the GS (14 species), the karyotype - based on chromosome number, heterochromatic banding and 5S and 45S rDNA localisation (18 species) - was characterised and analysed along with published data using phylogenetic approaches. The GS presented a high phylogenetic correlation and it was related to morphological groups in Bifrenaria (larger plants - higher GS). The two largest GS found among genera were caused by different mechanisms: polyploidy in Bifrenaria tyrianthina and accumulation of repetitive DNA in Scuticaria hadwenii. The chromosome number variability was caused mainly through descending dysploidy, and x=20 was estimated as the base chromosome number. Combining GS and karyotype data with molecular phylogeny, our data provide a more complete scenario of the karyotype evolution in Maxillariinae orchids, allowing us to suggest, besides dysploidy, that inversions and transposable elements as two mechanisms involved in the karyotype evolution. Such karyotype modifications could be associated with niche changes that occurred during species evolution.

In situ morphometric survey elucidates the evolutionary systematics of the Eurasian Himantoglossum clade (Orchidaceae: Orchidinae)

Background and Aims

The charismatic Himantoglossum s.l. clade of Eurasian orchids contains an unusually large proportion of taxa that are of controversial circumscriptions and considerable conservation concern. Whereas our previously published study addressed the molecular phylogenetics and phylogeography of every named taxon within the clade, here we use detailed morphometric data obtained from the same populations to compare genotypes with associated phenotypes, in order to better explore taxonomic circumscription and character evolution within the clade.

Methods

Between one and 12 plants found in 25 populations that encompassed the entire distribution of the Himantoglossum s.l. clade were measured in situ for 51 morphological characters. Results for 45 of those characters were subjected to detailed multivariate and univariate analyses.

Key Results

Multivariate analyses readily separate subgenus Barlia and subgenus Comperia from subgenus Himantoglossum, and also the early-divergent H. formosum from the less divergent remainder of subgenus Himantoglossum. The sequence of divergence of these four lineages is confidently resolved. Our experimental approach to morphometric character analysis demonstrates clearly that phenotypic evolution within Himantoglossum is unusually multi-dimensional.

Conclusions

Degrees of divergence between taxa shown by morphological analyses approximate those previously shown using molecular analyses. Himantoglossum s.l. is readily divisible into three subgenera. The three sections of subgenus Himantoglossum—hircinum, caprinum and formosum—are arrayed from west to east with only limited geographical overlap. At this taxonomic level, their juxtaposition combines with conflict between contrasting datasets to complicate attempts to distinguish between clinal variation and the discontinuities that by definition separate bona fide species. All taxa achieve allogamy via food deceit and have only weak pollinator specificity. Artificial crossing demonstrates that intrinsic sterility barriers are weak. Although we have found evidence of gene flow among and within the three sections of subgenus Himantoglossum, reports of natural hybrids are surprisingly rare, probably because putative parents are sufficiently similar to questionably warrant the status of species. Phenological separation and increased xeromorphy characterise the origin of subgenus Barlia. Several individual morphological characters show evidence of parallel acquisition, and loss of features is especially frequent in floral markings among members of section caprinum. Detailed patterns of gain and loss demonstrate that several different categories of flower markings are inherited independently. Along with the dimensions of labellar lobes, these pigmentation characters have been over-emphasised in previous taxonomic treatments. Increased plant vigour was a crucial element of the origin of the genus, but vegetative characters underwent remarkably little subsequent evolution. Attempts to reconstruct hypothetical ancestors at internal nodes of the phylogeny are weakened by (a) uncertain placement of Steveniella as sister to Himantoglossum s.l. and (b) uncertain relationships among subtly different putative species within section caprinum. Nonetheless, heterochronic/allometric trends, ultimately limited by functional constraints, clearly dictate transitions between contrasting flower sizes and complex labellum shapes.