Unlocking the Genetic Identity of Endangered Paphiopedilum Orchids: A DNA Barcoding Approach.

Orchids of the genus Paphiopedilum, also called slippers, are among the most valued representatives of the Orchidaceae family due to their aesthetic qualities. Due to overexploitation, deforestation, and illegal trade in these plants, especially in the vegetative phase, Paphiopedilum requires special protection. This genus is listed in Appendix I of the Convention on International Trade in Endangered Species of Wild Fauna and Flora. Their precise identification is of great importance for the preservation of genetic resources and biodiversity of the orchid family (Orchidaceae). Therefore, the main objective of the study was to investigate the usefulness of the DNA barcoding technique for the identification of endangered orchids of the genus Paphiopedilum and to determine the effectiveness of five loci: matK, rbcL, ITS2, atpF-atpH and trnH-psbA as potential molecular markers for species of this genus. Among single locus barcodes, matK was the most effective at identifying species (64%). Furthermore, matK, ITS2, matK + rbcL, and matK + trnH-psbA barcodes can be successfully used as a complementary tool to identify Paphiopedilum orchids while supporting morphological data provided by taxonomists.

Complete chloroplast genome structural characterization of two Aerides (Orchidaceae) species with a focus on phylogenetic position of Aerides flabellata.

BACKGROUND: The disputed phylogenetic position of Aerides flabellata Rolfe ex Downie, due to morphological overlaps with related species, was investigated based on evidence of complete chloroplast (cp) genomes. The structural characterization of complete cp genomes of A. flabellata and A. rosea Lodd. ex Lindl. & Paxton were analyzed and compared with those of six related species in "Vanda-Aerides alliance" to provide genomic information on taxonomy and phylogeny. RESULTS: The cp genomes of A. flabellata and A. rosea exhibited conserved quadripartite structures, 148,145 bp and 147,925 bp in length, with similar GC content (36.7 ~ 36.8%). Gene annotations revealed 110 single-copy genes, 18 duplicated in inverted regions, and ten with introns. Comparative analysis across related species confirmed stable sequence identity and higher variation in single-copy regions. However, there are notable differences in the IR regions between two Aerides Lour. species and the other six related species. The phylogenetic analysis based on CDS from complete cp genomes indicated that Aerides species except A. flabellata formed a monophyletic clade nested in the subtribe Aeridinae, being a sister group to Renanthera Lour., consistent with previous studies. Meanwhile, a separate clade consisted of A. flabellata and six Vanda R. Br. species was formed, as a sister taxon to Holcoglossum Schltr. CONCLUSIONS: This research was the first report on the complete cp genomes of A. flabellata. The results provided insights into understanding of plastome evolution and phylogenetic relationships of Aerides. The phylogenetic analysis based on complete cp genomes showed that A. flabellata should be placed in Vanda rather than in Aerides.

Physiological and transcriptomic analyses reveal the regulatory mechanisms of Anoectochilus roxburghii in response to high-temperature stress.

BACKGROUND: High temperatures significantly affect the growth, development, and yield of plants. Anoectochilus roxburghii prefers a cool and humid environment, intolerant of high temperatures. It is necessary to enhance the heat tolerance of A. roxburghii and breed heat-tolerant varieties. Therefore, we studied the physiological indexes and transcriptome of A. roxburghii under different times of high-temperature stress treatments. RESULTS: Under high-temperature stress, proline (Pro), H2O2 content increased, then decreased, then increased again, catalase (CAT) activity increased continuously, peroxidase (POD) activity decreased rapidly, then increased, then decreased again, superoxide dismutase (SOD) activity, malondialdehyde (MDA), and soluble sugars (SS) content all decreased, then increased, and chlorophyll and soluble proteins (SP) content increased, then decreased. Transcriptomic investigation indicated that a total of 2740 DEGs were identified and numerous DEGs were notably enriched for "Plant-pathogen interaction" and "Plant hormone signal transduction". We identified a total of 32 genes in these two pathways that may be the key genes for resistance to high-temperature stress in A. roxburghii. CONCLUSIONS: To sum up, the results of this study provide a reference for the molecular regulation of A. roxburghii's tolerance to high temperatures, which is useful for further cultivation of high-temperature-tolerant A. roxburghii varieties.

Genome-Wide Identification and Expression Analysis of the GRAS Gene Family and Their Responses to Heat Stress in Cymbidium goeringii.

The GRAS gene family, responsible for encoding transcription factors, serves pivotal functions in plant development, growth, and responses to stress. The exploration of the GRAS gene family within the Orchidaceae has been comparatively limited, despite its identification and functional description in various plant species. This study aimed to conduct a thorough examination of the GRAS gene family in Cymbidum goeringii, focusing on its physicochemical attributes, phylogenetic associations, gene structure, cis-acting elements, and expression profiles under heat stress. The results show that a total of 54 CgGRASs were pinpointed from the genome repository and categorized into ten subfamilies via phylogenetic associations. Assessment of gene sequence and structure disclosed the prevalent existence of the VHIID domain in most CgGRASs, with around 57.41% (31/54) CgGRASs lacking introns. The Ka/Ks ratios of all CgGRASs were below one, indicating purifying selection across all CgGRASs. Examination of cis-acting elements unveiled the presence of numerous elements linked to light response, plant hormone signaling, and stress responsiveness. Furthermore, CgGRAS5 contained the highest quantity of cis-acting elements linked to stress response. Experimental results from RT-qPCR demonstrated notable variations in the expression levels of eight CgGRASs after heat stress conditions, particularly within the LAS, HAM, and SCL4/7 subfamilies. In conclusion, this study revealed the expression pattern of CgGRASs under heat stress, providing reference for further exploration into the roles of CgGRAS transcription factors in stress adaptation.

Phosphorylation of 399S at CsHsp70 of Cymbidium sinense is essential to maintain chlorophyll stability.

The Chinese orchids symbolise nobility and gentility in China, and the variation of leaf color makes Cymbidium sinense more diversified and valuable. However, its color variations especially at the protein level still remain largely unexplored. In this study, the proteomics and phosphoproteomics of Cymbidium sinense leaf color variation mutants were studied. A total of 1059 differentially abundant proteins (DAPs) and 1127 differentially abundant phosphorylation sites belonging to 644 phosphoproteins (DAPPs) were identified in the yellow section of leaf variegation mutant of Cymbidium sinense (MY) compared with the green section (MG). Moreover, 349 co-expressing proteins were found in both omics' datasets, while only 26 proteins showed the same expression patterns in the two omics. The interaction network analysis of kinases and phosphatases showed that DAPs and DAPPs in photosynthesis, response to hormones, pigment metabolic process, phosphorylation, glucose metabolic process, and dephosphorylation might contribute to leaf color variation. The abundance of 28 Hsps and 28 phosphorylation sites belonging to 10 Hsps showed significant differences between MG and MY. CsHsp70 was selected to explore the function in Cymbidium sinense leaf variegation. The results showed CsHsp70 is essential for maintaining photosynthetic pigment content and the 399S phosphorylation site is crucial to the function of CsHsp70. Collectively, our findings construct a comprehensive coverage of protein and protein phosphorylation in leaf variegation of C. sinense, providing valuable insights into its formation mechanisms.

Genome-Wide Identification and Expression Pattern Analysis of TIFY Family Genes Reveal Their Potential Roles in Phalaenopsis aphrodite Flower Opening.

The TIFY gene family (formerly known as the zinc finger proteins expressed in inflorescence meristem (ZIM) family) not only functions in plant defense responses but also are widely involved in regulating plant growth and development. However, the identification and functional analysis of TIFY proteins remain unexplored in Orchidaceae. Here, we identified 19 putative TIFY genes in the Phalaenopsis aphrodite genome. The phylogenetic tree classified them into four subfamilies: 14 members from JAZ, 3 members from ZML, and 1 each from PPD and TIFY. Sequence analysis revealed that all Phalaenopsis TIFY proteins contained a TIFY domain. Exon-intron analysis showed that the intron number and length of Phalaenopsis TIFY genes varied, whereas the same subfamily and subgroup genes had similar exon or intron numbers and distributions. The most abundant cis-elements in the promoter regions of the 19 TIFY genes were associated with light responsiveness, followed by MeJA and ABA, indicating their potential regulation by light and phytohormones. The 13 candidate TIFY genes screened from the transcriptome data exhibited two types of expression trends, suggesting their different roles in cell proliferation and cell expansion of floral organ growth during Phalaenopsis flower opening. Overall, this study serves as a background for investigating the underlying roles of TIFY genes in floral organ growth in Phalaenopsis.

Assembly and comparative analysis of the complete multichromosomal mitochondrial genome of Cymbidium ensifolium, an orchid of high economic and ornamental value.

BACKGROUND: Orchidaceae is one of the largest groups of angiosperms, and most species have high economic value and scientific research value due to their ornamental and medicinal properties. In China, Chinese Cymbidium is a popular ornamental orchid with high economic value and a long history. However, to date, no detailed information on the mitochondrial genome of any species of Chinese Cymbidium has been published. RESULTS: Here, we present the complete assembly and annotation of the mitochondrial genome of Cymbidium ensifolium (L.) Sw. The mitogenome of C. ensifolium was 560,647 bp in length and consisted of 19 circular subgenomes ranging in size from 21,995 bp to 48,212 bp. The genome encoded 35 protein-coding genes, 36 tRNAs, 3 rRNAs, and 3405 ORFs. Repeat sequence analysis and prediction of RNA editing sites revealed a total of 915 dispersed repeats, 162 simple repeats, 45 tandem repeats, and 530 RNA editing sites. Analysis of codon usage showed a preference for codons ending in A/T. Interorganellar DNA transfer was identified in 13 of the 19 chromosomes, with plastid-derived DNA fragments representing 6.81% of the C. ensifolium mitochondrial genome. The homologous fragments of the mitochondrial genome and nuclear genome were also analysed. Comparative analysis showed that the GC content was conserved, but the size, structure, and gene content of the mitogenomes varied greatly among plants with multichromosomal mitogenome structure. Phylogenetic analysis based on the mitogenomes reflected the evolutionary and taxonomic statuses of C. ensifolium. Interestingly, compared with the mitogenomes of Cymbidium lancifolium Hook. and Cymbidium macrorhizon Lindl., the mitogenome of C. ensifolium lost 8 ribosomal protein-coding genes. CONCLUSION: In this study, we assembled and annotated the mitogenome of C. ensifolium and compared it with the mitogenomes of other Liliidae and plants with multichromosomal mitogenome structures. Our findings enrich the mitochondrial genome database of orchid plants and reveal the rapid structural evolution of Cymbidium mitochondrial genomes, highlighting the potential for mitochondrial genes to help decipher plant evolutionary history.

Phylogenomics and biogeographical diversification of Collabieae (Orchidaceae) and its implication in the reconstruction of the dynamic history of Asian evergreen broadleaved forests.

The tribe Collabieae (Epidendroideae, Orchidaceae) comprises approximately 500 species. Generic delimitation within Collabieae are confusing and phylogenetic interrelationships within the Collabieae have not been well resolved. Plastid genomes and nuclear internal transcribed spacer (ITS) sequences were used to estimate the phylogenetic relationships, ancestral ranges, and diversification rates of Collabieae. The results showed that Collabieae was subdivided into nine clades with high support. We proposed to combine Ancistrochilus and Pachystoma into Spathoglottis, merge Collabium and Chrysoglossum into Diglyphosa, and separate Pilophyllum and Hancockia as distinctive genera. The diversification of the nine clades of Collabieae might be associated with the uplift of the Himalayas during the Late Oligocene/Early Miocene. The enhanced East Asian summer monsoon in the Late Miocene may have promoted the rapid diversification of Collabieae at a sustained high diversification rate. The increased size of terrestrial pseudobulbs may be one of the drivers of Collabieae diversification. Our results suggest that the establishment and development of evergreen broadleaved forests facilitated the diversification of Collabieae.

Phenotypic variation seems not to be associated with the genetic profile in Zygopetalum (Orchidaceae): a case study of a high-elevation rocky complex.

BACKGROUND: Hybridization associated with polyploidy studies is rare in the tropics. The genus Zygopetalum (Orchidaceae) was investigated here as a case study of Neotropical plants. In the rocky highlands of the Ibitipoca State Park (ISP), southeast Brazil, individuals with intermediate colors and forms between the species Z. maculatum and Z. triste were commonly identified. METHODS AND RESULTS: Chromosomal analysis and DNA quantity showed a uniform population. Regardless of the aspects related to the color and shape of floral structures, all individuals showed 2n = 96 chromosomes and an average of 14.05 pg of DNA. Irregularities in meiosis associated with chromosome number and C value suggest the occurrence of polyploidy. The genetic distance estimated using ISSR molecular markers revealed the existence of genetic variability not related to morphological clusters. Morphometric measurements of the flower pieces revealed that Z. maculatum shows higher variation than Z. triste although lacking a defined circumscription. CONCLUSION: The observed variation can be explained by the polyploid and phenotypic plasticity resulting from the interaction of the genotypes with the heterogeneous environments observed in this habitat.

The Complete Chloroplast Genomes of Bulbophyllum (Orchidaceae) Species: Insight into Genome Structure Divergence and Phylogenetic Analysis.

Bulbophyllum is one of the largest genera and presents some of the most intricate taxonomic problems in the family Orchidaceae, including species of ornamental and medical importance. The lack of knowledge regarding the characterization of Bulbophyllum chloroplast (cp) genomes has imposed current limitations on our study. Here, we report the complete cp genomes of seven Bulbophyllum species, including B. ambrosia, B. crassipes, B. farreri, B. hamatum, B. shanicum, B. triste, and B. violaceolabellum, and compared with related taxa to provide a better understanding of their genomic information on taxonomy and phylogeny. A total of 28 Bulbophyllum cp genomes exhibit typical quadripartite structures with lengths ranging from 145,092 bp to 165,812 bp and a GC content of 36.60% to 38.04%. Each genome contained 125-132 genes, encompassing 74-86 protein-coding genes, 38 tRNA genes, and eight rRNA genes. The genome arrangements, gene contents, and length were similar, with differences observed in ndh gene composition. It is worth noting that there were exogenous fragment insertions in the IR regions of B. crassipes. A total of 18-49 long repeats and 38-80 simple sequence repeats (SSRs) were detected and the single nucleotide (A/T) was dominant in Bulbophyllum cp genomes, with an obvious A/T preference. An analysis of relative synonymous codon usage (RSCU) revealed that leucine (Leu) was the most frequently used codon, while cysteine (Cys) was the least used. Six highly variable regions (rpl32-trnLUAG > trnTUGU-trnLUAA > trnFGAA-ndhJ > rps15-ycf1 > rbcL-accD > psbI-trnSGCU) and five coding sequences (ycf1 > rps12 > matK > psbK > rps15) were identified as potential DNA markers based on nucleotide diversity. Additionally, 31,641 molecular diagnostic characters (MDCs) were identified in complete cp genomes. A phylogenetic analysis based on the complete cp genome sequences and 68 protein-coding genes strongly supported that 28 Bulbophyllum species can be divided into four branches, sects. Brachyantha, Cirrhopetalum, and Leopardinae, defined by morphology, were non-monophyly. Our results enriched the genetic resources of Bulbophyllum, providing valuable information to illustrate the complicated taxonomy, phylogeny, and evolution process of the genus.

Organelle Genomes of Epipogium roseum Provide Insight into the Evolution of Mycoheterotrophic Orchids.

Epipogium roseum, commonly known as one of the ghost orchids due to its rarity and almost transparent color, is a non-photosynthetic and fully mycoheterotrophic plant. Given its special nutritional strategies and evolutionary significance, the mitogenome was first characterized, and three plastomes sampled from Asia were assembled. The plastomes were found to be the smallest among Orchidaceae, with lengths ranging from 18,339 to 19,047 bp, and exhibited high sequence variety. For the mitogenome, a total of 414,552 bp in length, comprising 26 circular chromosomes, were identified. A total of 54 genes, including 38 protein-coding genes, 13 tRNA genes, and 3 rRNA genes, were annotated. Multiple repeat sequences spanning a length of 203,423 bp (45.47%) were discovered. Intriguingly, six plastid regions via intracellular gene transfer and four plastid regions via horizontal gene transfer to the mitogenome were observed. The phylogenomics, incorporating 90 plastomes and 56 mitogenomes, consistently revealed the sister relationship of Epipogium and Gastrodia, with a bootstrap percentage of 100%. These findings shed light on the organelle evolution of Orchidaceae and non-photosynthetic plants.

Identification and Analysis of PEPC Gene Family Reveals Functional Diversification in Orchidaceae and the Regulation of Bacterial-Type PEPC.

Phosphoenolpyruvate carboxylase (PEPC) gene family plays a crucial role in both plant growth and response to abiotic stress. Approximately half of the Orchidaceae species are estimated to perform CAM pathway, and the availability of sequenced orchid genomes makes them ideal subjects for investigating the PEPC gene family in CAM plants. In this study, a total of 33 PEPC genes were identified across 15 orchids. Specifically, one PEPC gene was found in Cymbidium goeringii and Platanthera guangdongensis; two in Apostasia shenzhenica, Dendrobium chrysotoxum, D. huoshanense, Gastrodia elata, G. menghaiensis, Phalaenopsis aphrodite, Ph. equestris, and Pl. zijinensis; three in C. ensifolium, C. sinense, D. catenatum, D. nobile, and Vanilla planifolia. These PEPC genes were categorized into four subgroups, namely PEPC-i, PEPC-ii, and PEPC-iii (PTPC), and PEPC-iv (BTPC), supported by the comprehensive analyses of their physicochemical properties, motif, and gene structures. Remarkably, PEPC-iv contained a heretofore unreported orchid PEPC gene, identified as VpPEPC4. Differences in the number of PEPC homolog genes among these species were attributed to segmental duplication, whole-genome duplication (WGD), or gene loss events. Cis-elements identified in promoter regions were predominantly associated with light responsiveness, and circadian-related elements were observed in each PEPC-i and PEPC-ii gene. The expression levels of recruited BTPC, VpPEPC4, exhibited a lower expression level than other VpPEPCs in the tested tissues. The expression analyses and RT-qPCR results revealed diverse expression patterns in orchid PEPC genes. Duplicated genes exhibited distinct expression patterns, suggesting functional divergence. This study offered a comprehensive analysis to unveil the evolution and function of PEPC genes in Orchidaceae.

First report of the complete mitochondrial genome of 3 beetles (Coleoptera: Scarabaeidae) harming Gastrodia elata (Asparagales: Orchidaceae).

Gastrodia elata Blume, a valuable traditional Chinese medicine with significant clinical and nutritional importance, is a fungal heterotrophic orchid. We present the first report of the mitochondrial genome structure and characteristics of 3 Scarabaeidae pests affecting G. elata: Sophrops peronosporus Gu & Zhang, Anomala rufiventris Kollar & Redtenbacher, and Callistethus plagiicollis Fairmaire. Each mitogenome contained 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), 2 ribosomal RNAs (rRNAs), and a control region, with no gene rearrangements observed. All 21 tRNAs, except trnS1 that lacks a dihydrouridine, had a stable cloverleaf secondary structure. Maximum likelihood and Bayesian inference analyses based on the 13 PCGs produced 2 topologically similar phylogenetic trees, both of with high nodal support. Larvae of these Scarabaeidae pests cause substantial damage by gnawing on the tubers and roots of G. elata, leading to reduced yield and compromised quality. These findings contribute to phylogenetic studies of Scarabaeidae, expand knowledge of G. elata pests, and offer valuable reference materials for their identification and control.

Mechanism underlying the rapid growth of Phalaenopsis equestris induced by 60Co-γ-ray irradiation.

Gamma (γ)-ray irradiation is one of the important modern breeding methods. Gamma-ray irradiation can affect the growth rate and other characteristics of plants. Plant growth rate is crucial for plants. In horticultural crops, the growth rate of plants is closely related to the growth of leaves and flowering time, both of which have important ornamental value. In this study, 60Co-γ-ray was used to treat P. equestris plants. After irradiation, the plant's leaf growth rate increased, and sugar content and antioxidant enzyme activity increased. Therefore, we used RNA-seq technology to analyze the differential gene expression and pathways of control leaves and irradiated leaves. Through transcriptome analysis, we investigated the reasons for the rapid growth of P. equestris leaves after irradiation. In the analysis, genes related to cell wall relaxation and glucose metabolism showed differential expression. In addition, the expression level of genes encoding ROS scavenging enzyme synthesis regulatory genes increased after irradiation. We identified two genes related to P. equestris leaf growth using VIGS technology: PeNGA and PeEXPA10. The expression of PeEXPA10, a gene related to cell wall expansion, was down-regulated, cell wall expansion ability decreased, cell size decreased, and leaf growth rate slowed down. The TCP-NGATHA (NGA) molecular regulatory module plays a crucial role in cell proliferation. When the expression of the PeNGA gene decreases, the leaf growth rate increases, and the number of cells increases. After irradiation, PeNGA and PeEXPA10 affect the growth of P. equestris leaves by influencing cell proliferation and cell expansion, respectively. In addition, many genes in the plant hormone signaling pathway show differential expression after irradiation, indicating the crucial role of plant hormones in plant leaf growth. This provides a theoretical basis for future research on leaf development and biological breeding.

Breeding of ornamental orchids with focus on Phalaenopsis: current approaches, tools, and challenges for this century.

Ornamental orchid breeding programs have been conducted to develop commercially valuable cultivars with improved characteristics of commercial interest, such as size, flower color, pattern, shape, and resistance to pathogens. Conventional breeding, including sexual hybridization followed by selection of desirable characteristics in plants, has so far been the main method for ornamental breeding, but other techniques, including mutation induction by polyploidization and gamma irradiation, and biotechnological techniques, such as genetic transformation, have also been studied and used in ornamental breeding programs. Orchids are one of the most commercially important families in floriculture industry, having very particular reproductive biology characteristics and being a well-studied group of ornamentals in terms of genetic improvement. The present review focuses on the conventional and biotechnological techniques and approaches specially employed in breeding Phalaenopsis orchids, the genus with highest worldwide importance as an ornamental orchid, highlighting the main limitations and strengths of the approaches. Furthermore, new opportunities and future prospects for ornamental breeding in the CRISPR/Cas9 genome editing era are also discussed. We conclude that conventional hybridization remains the most used method to obtain new cultivars in orchids. However, the emergence of the first biotechnology-derived cultivars, as well as the new biotechnological tools available, such as CRISPR-Cas9, rekindled the full potential of biotechnology approaches and their importance for improve ornamental orchid breeding programs.

The origin and speciation of orchids.

Orchids constitute one of the most spectacular radiations of flowering plants. However, their origin, spread across the globe, and hotspots of speciation remain uncertain due to the lack of an up-to-date phylogeographic analysis. We present a new Orchidaceae phylogeny based on combined high-throughput and Sanger sequencing data, covering all five subfamilies, 17/22 tribes, 40/49 subtribes, 285/736 genera, and c. 7% (1921) of the 29 524 accepted species, and use it to infer geographic range evolution, diversity, and speciation patterns by adding curated geographical distributions from the World Checklist of Vascular Plants. The orchids' most recent common ancestor is inferred to have lived in Late Cretaceous Laurasia. The modern range of Apostasioideae, which comprises two genera with 16 species from India to northern Australia, is interpreted as relictual, similar to that of numerous other groups that went extinct at higher latitudes following the global climate cooling during the Oligocene. Despite their ancient origin, modern orchid species diversity mainly originated over the last 5 Ma, with the highest speciation rates in Panama and Costa Rica. These results alter our understanding of the geographic origin of orchids, previously proposed as Australian, and pinpoint Central America as a region of recent, explosive speciation.

Comparative chloroplast genomics of 24 species shed light on the genome evolution and phylogeny of subtribe Coelogyninae (Orchidaceae).

BACKGROUND: The orchids of the subtribe Coelogyninae are among the most morphologically diverse and economically important groups within the subfamily Epidendroideae. Previous molecular studies have revealed that Coelogyninae is an unambiguously monophyletic group. However, intergeneric and infrageneric relationships within Coelogyninae are largely unresolved. There has been long controversy over the classification among the genera within the subtribe. RESULTS: The complete chloroplast (cp.) genomes of 15 species in the subtribe Coelogyninae were newly sequenced and assembled. Together with nine available cp. genomes in GenBank from representative clades of the subtribe, we compared and elucidated the characteristics of 24 Coelogyninae cp. genomes. The results showed that all cp. genomes shared highly conserved structure and contained 135 genes arranged in the same order, including 89 protein-coding genes, 38 tRNAs, and eight rRNAs. Nevertheless, structural variations in relation to particular genes at the IR/SC boundary regions were identified. The diversification pattern of the cp. genomes showed high consistency with the phylogenetic placement of Coelogyninae. The number of different types of SSRs and long repeats exhibited significant differences in the 24 Coelogyninae cp. genomes, wherein mononucleotide repeats (A/T), and palindromic repeats were the most abundant. Four mutation hotspot regions (ycf1a, ndhF-rp132, psaC-ndhE, and rp132-trnL) were determined, which could serve as effective molecular markers. Selection pressure analysis revealed that three genes (ycf1a, rpoC2 and ycf2 genes) might have experienced apparent positive selection during the evolution. Using the alignments of whole cp. genomes and protein-coding sequences, this study presents a well-resolved phylogenetic framework of Coelogyninae. CONCLUSION: The inclusion of 55 plastid genome data from a nearly complete generic-level sampling provide a comprehensive view of the phylogenetic relationships among genera and species in subtribe Coelogyninae and illustrate the diverse genetic variation patterns of plastid genomes in this species-rich plant group. The inferred relationships and informally recognized major clades within the subtribe are presented. The genetic markers identified here will facilitate future studies on the genetics and phylogeny of subtribe Coelogyninae.

The absence of bumblebees on an oceanic island blurs the species boundary of two closely related orchids.

Oceanic islands offer valuable natural laboratories for studying evolution. The Izu Islands, with their recent geological origin, provide an exceptional opportunity to explore the initial evolution on oceanic islands. Another noteworthy aspect is the absence of bumblebee species on most Izu Islands. We used ecological, morphological, and molecular data to investigate the impact of bumblebee absence on the evolution of two closely related orchid species, Goodyera henryi and Goodyera similis, focusing on Kozu Island, the Izu Islands. Our investigation revealed that while G. henryi exclusively relies on a bumblebee species for pollination on the mainland, G. similis is pollinated by scoliid wasps on both the mainland and the island. Intriguingly, all specimens initially categorized as G. henryi on Kozu Island are hybrids of G. henryi and G. similis, leading to the absence of pure G. henryi distribution on the island. These hybrids are pollinated by the scoliid wasp species that also pollinates G. similis on the island. The absence of bumblebees might result in sporadic and inefficient pollination of G. henryi by scoliid wasps, consequently promoting hybrid proliferation on the island. Our findings suggest that the absence of bumblebees can blur plant species boundaries.

Metagenomic Analyses of Viruses in the Orchid Mycorrhizal Interaction Using Improved Assemble Tools.

In nature, mycorrhizal association with soil-borne fungi is indispensable for orchid species. Compatible mycorrhizal fungi form endo-mycorrhizal structures in orchid cells, and the fungal structures are digested in orchid cells to be supplied to orchids as nutrition. Because orchid seeds lack the reserves for germination, they keep receiving nutrition through mycorrhizal formation from seed germination until nonphotosynthetic young seedlings develop leaves and become photoautotrophic. Seeds of all orchids germinate with the help of their own fungal partners, and therefore, specific partnership has been acquired in a long evolutionary history between orchids and fungi. Assuming that horizontal transmission of viruses may occur in such a close relationship, we are focusing on viruses that infect orchids and their mycorrhizal fungi. We prepared aseptically germinated orchid plants (i.e., fungi-free plants) together with pure-cultured fungal isolates and conducted transcriptome analyses (RNA-seq) by next-generation sequencing (NGS) approach. To reconstruct virus-related sequences that would have been present in the RNA sample of interest, de novo assembly process is required using short read sequences obtained from RNA-seq. In the previous version of our protocol (see Viral Metagenomics, first edition 2018), virus searches were conducted using contig sets constructed by a single assembler, but this time we devised a method to construct more reliable contigs using multiple assemblers and again reinvestigated that viruses could be detected. Because the virus detection efficiency and number of detected virus species clearly differed depending on the assembly pipeline and the number of the input data, multiple methods should be used to identify viral infection, if possible.

Sex-Related Gene Network Revealed by Transcriptome Differentiation of Bisexual and Unisexual Flowers of Orchid Cymbidium tortisepalum.

Despite extensive research on orchid reproductive strategies, the genetic studies of sex differentiation in the orchid family are still lacking. In this study, we compared three sexual phenotypes of Cymbidium tortisepalum bisexual flowers as well as female and male unisexual mutants. Through comparative transcriptomes, we analyzed the sex-biased differentially expressed genes (DEGs) and gene co-expression networks of sex organs (gynostemium and ovary) among them, identified the candidate genes of sex differentiation, and validated their expression by qRT-PCR. The C. tortisepalum unisexual mutants with degenerated phenotypes were compared to the bisexual plants with respect to both the flower organs and plant morphologies. Totally, 12,145, 10,789, and 14,447 genes were uniquely expressed in the female, male, and hermaphrodite sex organs, respectively. A total of 4291 sex-biased DEGs were detected among them, with 871, 2867, and 1937 DEGs in the comparisons of bisexual vs. female, bisexual vs. male, and male vs. female flowers, respectively. Two co-expressed network modules, with 81 and 419 genes were tightly correlated with female sexual traits, while two others with 265 and 135 genes were highly correlated with male sexual traits. Two female-biased hub genes (CtSDR3b and CtSDR3b-like) nested in the female modules, the homologs of maize sex determinant tasselseed2, may control the feminization of C. tortisepalum. At the same time, two male-biased hub genes (CtYAB2 and CtYAB5) nested in the male modules, the homologs of grape sex determinant VviYABBY3, may control the androphany of C. tortisepalum. This study discovered the molecular regulation networks and proposed a model for orchid sex differentiation, therefore providing for the first time the genetic basis for the sex separation in the orchid family.