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.

First Report of Agroathelia rolfsii Causing Southern Blight on Lippia (Phyla canescens) in Guangzhou, China.

Lippia (Phyla canescens) is a fast-growing, mat-forming, and prostrate perennial plant well adapted to infertile, high-saline, and drought environments (Leigh, et al. 2004). It arrived in China from Japan as a flowering ground cover in 2001 (Cai, et al. 2004). In June 2022, southern blight appeared in our nursery of the Floriculture Research Institute of Guangdong Academy of Agricultural Sciences. High temperature and damp environment are major factors for this disease. The symptoms of top-layer plants were not easily detected, but they were slightly yellowed. A yellowish-brown water-soak lesion appeared on the stems and lowest leaves exposed to soil. White mycelium appeared in the middle stage. Finally, the surface plants showed water-soak decay, and a mass of beige to black-brown rapeseed-shaped sclerotia appeared on the residue and surrounding soil; these plants died. Sclerotia and mycelia were collected from disease tissue, and after surface sterilization, sclerotia was cultured on potato dextrose agar (PDA) at 28±2°C in an incubator without light. Eight fungal isolates with similar colony morphologies were consistently isolated by purifying from different sampling areas. The isolates exhibited obvious septa and a clamp connection structure within the white mycelium. The average growth rate was 26.86±0.06 mm/day. Numerous white granular sclerotia were produced on the mycelium 6 days later. The sclerotia with a diameter of 1.24±0.07mm (n=189) gradually changed from diage to yellow to brown. A typical strain B1 was selected for further identification, targeting its 18S rRNA and LSU rRNA sequences (Yang, et al. 2011; Xue, et al. 2019). Its 18S rRNA sequence (GenBank Accession No. OR517233, 1626 bp) is 99.63% and 99.57% identical to Athelia rolfsii (AY665774, 1179bp; KC670714, 1775bp; JF819726, 1781bp). Its LSU rRNA sequence (OR539570, 757 bp) is 99.87% identical to Agroathelia rolfsii (OR526537, 904 bp). For Athelia rolfsii, a synonym of Agroathelia rolfsii, by combining the morphological characteristics and molecular identification, the isolate pathogen B1 was confirmed to be Agroathelia rolfsii (the teleomorph of Sclerotium rolfsii). To fullfill Koch's postulates, we inoculated the mycelial plugs to healthy lippia stems and leaves which has grown for one year, with PDA plugs free of mycelium as the control. All the plants were kept in a greenhouse at 28±2°C with a 14-h photoperiod and 80% relative humidity. Each treatment was repeated thrice and vaccinated with 6 points. At 7 d following inoculation, all plants inoculated with B1 showed typical symptoms, but the control group was asymptomatic, and sclerotia appeared 17d after inoculation. Using the same protocol mentioned above, pathogenic fungal was reisolated only from treated groups, but not from the control group. Chose three of the pathogens for 18S rRNA and LSU rRNA sequencing, the results showed 100% identity to B1, the same as its microstructure. There are few reports about the disease on P. canescens. Sosa (2007) investigated the pathogens on P. canescens in Argentina, 16 fungi were found but no A. rolfsii. Sclerotium rolfsii were identified on P. nodiflora or P. lanceolata (Michaux) Greene in America (Farr, et al. 1989). To our knowledge, this is the first report in China. Because this pathogen has wide-ranging hosts and causes serious damage, the results from this study will offer guidance for the prevention and treatment of this disease.

Integrated proteomic, transcriptomic, and metabolomic profiling reveals that the gibberellin-abscisic acid hub runs flower development in the Chinese orchid Cymbidium sinense.

The seasonal flowering Chinese Cymbidium produce an axillary floral meristem and require a dormancy period during cold conditions for flower development. However, the bud activation mechanism remains elusive. This study evaluates the multi-omics across six stages of flower development, along with functional analysis of core genes to decipher the innate mechanism of floral bud initiation and outgrowth in the Chinese orchid Cymbidium sinense. Transcriptome and proteome analyses identified 10 modules with essential roles in floral bud dormancy and activation. Gene clusters in the early stages of flower development were mainly related to flowering time regulation and meristem determination, while the late stages were correlated with hormone signaling pathways. The metabolome identified 69 potential hormones in which gibberellin (GA) and abscisic acid (ABA) were the main regulatory hubs, and GA4 and GA53 exhibited a reciprocal loop. Extraneous GA application caused rapid elongation of flower buds and promoted the expression of flower development genes. Contrarily, exogenous ABA application extended the dormancy process and ABA inhibitors induced dormancy release. Moreover, CsAPETALA1 (CsAP1) was identified as the potential target of ABA for floral bud activation. Transformation of CsAP1 in Arabidopsis and its transient overexpression in C. sinense protoplasts not only affected flowering time and floral organ morphogenesis in Arabidopsis but also orchestrated the expression of flowering and hormone regulatory genes. The presence of ABA response elements in the CsAP1 promoter, rapid downregulation of CsAP1 after exogenous ABA application, and the activation of the floral bud after ABA inhibitor treatment suggest that ABA can control bud outgrowth through CsAP1.

Comparative chloroplast genomics, phylogenetic relationships and molecular markers development of Aglaonema commutatum and seven green cultivars of Aglaonema.

Aglaonema commutatum is a famous species in the Aglaonema genus, which has important ornamental and economic value. However, its chloroplast genome information and phylogenetic relationships among popular green cultivars of Aglaonema in southern China have not been reported. Herein, chloroplast genomes of one variety of A. commutatum and seven green cultivars of Aglaonema, namely, A. commutatum 'San Remo', 'Kai Sa', 'Pattaya Beauty', 'Sapphire', 'Silver Queen', 'Snow White', 'White Gem', and 'White Horse Prince', were sequenced and assembled for comparative analysis and phylogeny. These eight genomes possessed a typical quadripartite structure that consisted of a LSC region (90,799-91,486 bp), an SSC region (20,508-21,137 bp) and a pair of IR regions (26,661-26,750 bp). Each genome contained 112 different genes, comprising 79 protein-coding genes, 29 tRNA genes and 4 rRNA genes. The gene orders, GC contents, codon usage frequency, and IR/SC boundaries were highly conserved among these eight genomes. Long repeats, SSRs, SNPs and indels were analyzed among these eight genomes. Comparative analysis of 15 Aglaonema chloroplast genomes identified 7 highly variable regions, including trnH-GUG-exon1-psbA, trnS-GCU-trnG-UCC-exon1, trnY-GUA-trnE-UUC, psbC-trnS-UGA, trnF-GAA-ndhJ, ccsA-ndhD, and rps15-ycf1-D2. Reconstruction of the phylogenetic trees based on chloroplast genomes, strongly supported that Aglaonema was a sister to Anchomanes, and that the Aglaonema genus was classified into two sister clades including clade I and clade II, which corresponded to two sections, Aglaonema and Chamaecaulon, respectively. One variety and five cultivars, including A. commutatum 'San Remo', 'Kai Sa', 'Pattaya Beauty', 'Silver Queen', 'Snow White', and 'White Horse Prince', were classified into clade I; and the rest of the two cultivars, including 'Sapphire' and 'White Gem', were classified into clade II. Positive selection was observed in 34 protein-coding genes at the level of the amino acid sites among 77 chloroplast genomes of the Araceae family. Based on the highly variable regions and SSRs, 4 DNA markers were developed to differentiate the clade I and clade II in Aglaonema. In conclusion, this study provided chloroplast genomic resources for Aglaonema, which were useful for its classification and phylogeny.

Pigment Diversity in Leaves of Caladium × hortulanum Birdsey and Transcriptomic and Metabolic Comparisons between Red and White Leaves.

Leaf color is a key ornamental characteristic of cultivated caladium (Caladium × hortulanum Birdsey), a plant with diverse leaf colors. However, the genetic improvement of leaf color in cultivated caladium is hindered by the limited understanding of leaf color diversity and regulation. In this study, the chlorophyll and anthocyanin content of 137 germplasm resources were measured to explore the diversity and mechanism of leaf color formation in cultivated caladium. Association analysis of EST-SSR markers and pigment traits was performed, as well as metabolomics and transcriptomics analysis of a red leaf variety and its white leaf mutant. We found significant differences in chlorophyll and anthocyanin content among different color groups of cultivated caladium, and identified three, eight, three, and seven EST-SSR loci significantly associated with chlorophyll-a, chlorophyll-b, total chlorophyll and total anthocyanins content, respectively. The results further revealed that the white leaf mutation was caused by the down-regulation of various anthocyanins (such as cyanidin-3-O-rutinoside, quercetin-3-O-glucoside, and others). This change in concentration is likely due to the down-regulation of key genes (four PAL, four CHS, six CHI, eight F3H, one F3'H, one FLS, one LAR, four DFR, one ANS and two UFGT) involved in anthocyanin biosynthesis. Concurrently, the up-regulation of certain genes (one FLS and one LAR) that divert the anthocyanin precursors to other pathways was noted. Additionally, a significant change in the expression of numerous transcription factors (12 NAC, 12 bZIP, 23 ERF, 23 bHLH, 19 MYB_related, etc.) was observed. These results revealed the genetic and metabolic basis of leaf color diversity and change in cultivated caladium, and provided valuable information for molecular marker-assisted selection and breeding of leaf color in this ornamental plant.

Thirteen complete chloroplast genomes of the costaceae family: insights into genome structure, selective pressure and phylogenetic relationships.

BACKGROUND: Costaceae, commonly known as the spiral ginger family, consists of approximately 120 species distributed in the tropical regions of South America, Africa, and Southeast Asia, of which some species have important ornamental, medicinal and ecological values. Previous studies on the phylogenetic and taxonomic of Costaceae by using nuclear internal transcribed spacer (ITS) and chloroplast genome fragments data had low resolutions. Additionally, the structures, variations and molecular evolution of complete chloroplast genomes in Costaceae still remain unclear. Herein, a total of 13 complete chloroplast genomes of Costaceae including 8 newly sequenced and 5 from the NCBI GenBank database, representing all three distribution regions of this family, were comprehensively analyzed for comparative genomics and phylogenetic relationships. RESULT: The 13 complete chloroplast genomes of Costaceae possessed typical quadripartite structures with lengths from 166,360 to 168,966 bp, comprising a large single copy (LSC, 90,802 - 92,189 bp), a small single copy (SSC, 18,363 - 20,124 bp) and a pair of inverted repeats (IRs, 27,982 - 29,203 bp). These genomes coded 111 - 113 different genes, including 79 protein-coding genes, 4 rRNA genes and 28 - 30 tRNAs genes. The gene orders, gene contents, amino acid frequencies and codon usage within Costaceae were highly conservative, but several variations in intron loss, long repeats, simple sequence repeats (SSRs) and gene expansion on the IR/SC boundaries were also found among these 13 genomes. Comparative genomics within Costaceae identified five highly divergent regions including ndhF, ycf1-D2, ccsA-ndhD, rps15-ycf1-D2 and rpl16-exon2-rpl16-exon1. Five combined DNA regions (ycf1-D2 + ndhF, ccsA-ndhD + rps15-ycf1-D2, rps15-ycf1-D2 + rpl16-exon2-rpl16-exon1, ccsA-ndhD + rpl16-exon2-rpl16-exon1, and ccsA-ndhD + rps15-ycf1-D2 + rpl16-exon2-rpl16-exon1) could be used as potential markers for future phylogenetic analyses and species identification in Costaceae. Positive selection was found in eight protein-coding genes, including cemA, clpP, ndhA, ndhF, petB, psbD, rps12 and ycf1. Maximum likelihood and Bayesian phylogenetic trees using chloroplast genome sequences consistently revealed identical tree topologies with high supports between species of Costaceae. Three clades were divided within Costaceae, including the Asian clade, Costus clade and South American clade. Tapeinochilos was a sister of Hellenia, and Parahellenia was a sister to the cluster of Tapeinochilos + Hellenia with strong support in the Asian clade. The results of molecular dating showed that the crown age of Costaceae was about 30.5 Mya (95% HPD: 14.9 - 49.3 Mya), and then started to diverge into the Costus clade and Asian clade around 23.8 Mya (95% HPD: 10.1 - 41.5 Mya). The Asian clade diverged into Hellenia and Parahellenia at approximately 10.7 Mya (95% HPD: 3.5 - 25.1 Mya). CONCLUSION: The complete chloroplast genomes can resolve the phylogenetic relationships of Costaceae and provide new insights into genome structures, variations and evolution. The identified DNA divergent regions would be useful for species identification and phylogenetic inference in Costaceae.

Genome-wide association analysis identified molecular markers and candidate genes for flower traits in Chinese orchid (Cymbidium sinense).

The orchid, the champagne of flowers, brings luxury, elegance, and novelty to nature. Cymbidium sinense is a symbol of gigantic floral variability on account of wavering shapes and sizes of floral organs, although marker-trait association (MTA) has not been studied for its floral traits. We evaluated markers associated with 14 floral traits of C. sinense through a genome-wide association study (GWAS) of 195 accessions. A total of 65 318 522 single-nucleotide polymorphisms (SNPs) and 3 906 176 insertion/deletion (InDel) events were identified through genotyping-by-sequencing. Among these, 4694 potential SNPs and 477 InDels were identified as MTAs at -log10 P > 5. The genes related to these SNPs and InDels were largely associated with floral regulators, hormonal pathways, cell division, and metabolism, playing essential roles in tailoring floral morphology. Moreover, 20 candidate SNPs/InDels linked to 11 genes were verified, 8 of which were situated on exons, one was located in the 5'-UTR and two were positioned in introns. Here, the multitepal trait-related gene RABBIT EARS (RBE) was found to be the most crucial gene. We analyzed the role of CsRBE in the regulation of flower-related genes via efficient transient overexpression in C. sinense protoplasts, and found that the floral homeotic genes CsAP3 and CsPI, as well as organ boundary regulators, including CsCUC and CsTCP genes, were regulated by CsRBE. Thus, we obtained key gene loci for important ornamental traits of orchids using genome-wide association analysis of populations with natural variation. The findings of this study can do a great deal to expedite orchid breeding programs for shape variability.

Genome-Wide Identification, Expression, and Molecular Characterization of the CONSTANS-like Gene Family in Seven Orchid Species.

The orchid is one of the most distinctive and highly valued flowering plants. Nevertheless, the CONSTANS-like (COL) gene family plays significant roles in the control of flowering, and its functions in Orchidaceae have been minimally explored. This research identified 68 potential COL genes within seven orchids' complete genome, divided into three groups (groups I, II, and III) via a phylogenetic tree. The modeled three-dimensional structure and the conserved domains exhibited a high degree of similarity among the orchid COL proteins. The selection pressure analysis showed that all orchid COLs suffered a strong purifying selection. Furthermore, the orchid COL genes exhibited functional and structural heterogeneity in terms of collinearity, gene structure, cis-acting elements within their promoters, and expression patterns. Moreover, we identified 50 genes in orchids with a homology to those involved in the COL transcriptional regulatory network in Arabidopsis. Additionally, the first overexpression of CsiCOL05 and CsiCOL09 in Cymbidium sinense protoplasts suggests that they may antagonize the regulation of flowering time and gynostemium development. Our study will undoubtedly provide new resources, ideas, and values for the modern breeding of orchids and other plants.

Comparative Chloroplast Genomics of 21 Species in Zingiberales with Implications for Their Phylogenetic Relationships and Molecular Dating.

Zingiberales includes eight families and more than 2600 species, with many species having important economic and ecological value. However, the backbone phylogenetic relationships of Zingiberales still remain controversial, as demonstrated in previous studies, and molecular dating based on chloroplast genomes has not been comprehensively studied for the whole order. Herein, 22 complete chloroplast genomes from 21 species in Zingiberales were sequenced, assembled, and analyzed. These 22 genomes displayed typical quadripartite structures, which ranged from 161,303 bp to 163,979 bp in length and contained 111-112 different genes. The genome structures, gene contents, simple sequence repeats, long repeats, and codon usage were highly conserved, with slight differences among these genomes. Further comparative analysis of the 111 complete chloroplast genomes of Zingiberales, including 22 newly sequenced ones and the remaining ones from the national center for biotechnology information (NCBI) database, identified three highly divergent regions comprising ccsA, psaC, and psaC-ndhE. Maximum likelihood and Bayesian inference phylogenetic analyses based on chloroplast genome sequences found identical topological structures and identified a strongly supported backbone of phylogenetic relationships. Cannaceae was sister to Marantaceae, forming a clade that was collectively sister to the clade of (Costaceae, Zingiberaceae) with strong support (bootstrap (BS) = 100%, and posterior probability (PP) = 0.99-1.0); Heliconiaceae was sister to the clade of (Lowiaceae, Strelitziaceae), then collectively sister to Musaceae with strong support (BS = 94-100%, and PP = 0.93-1.0); the clade of ((Cannaceae, Marantaceae), (Costaceae, Zingiberaceae)) was sister to the clade of (Musaceae, (Heliconiaceae, (Lowiaceae, Strelitziaceae))) with robust support (BS = 100%, and PP = 1.0). The results of divergence time estimation of Zingiberales indicated that the crown node of Zingiberales occurred approximately 85.0 Mya (95% highest posterior density (HPD) = 81.6-89.3 million years ago (Mya)), with major family-level lineages becoming from 46.8 to 80.5 Mya. These findings proved that chloroplast genomes could contribute to the study of phylogenetic relationships and molecular dating in Zingiberales, as well as provide potential molecular markers for further taxonomic and phylogenetic studies of Zingiberales.

Cross-Compatibility in Interspecific Hybridization of Different Curcuma Accessions.

Curcuma is extensively cultivated as a medicinal and ornamental plant in tropical and subtropical regions. Due to the bright bract color, distinctive inflorescence and long blooming period, it has become a new favorite in terms of the urban landscape, potted flowers and cut flowers. However, little research on breeding new cultivars using traditional plant breeding methods is available on the genus Curcuma. In the present study, pollen viability and stigma receptivity evaluation were performed, and the genetic relationship of 38 Curcuma accessions was evaluated, then 5 C. alismatifolia Gagnep. (Ca), 2 C. hybrid (Ch), 2 C. sparganiifolia Gagnep. cultivars and 4 Curcuma native species were selected as parents for subsequent interspecific cross-breeding. A total of 132 reciprocal crosses were carried out for interspecific hybridization, including 70 obverse and 62 inverse crosses. Obvious discrepancies among fruit-setting rates were manifested in different combinations and in reciprocal crosses. Results showed that the highest fruit-setting rate (87.5%) was observed in the Ca combinations. There were 87 combinations with a fruit-setting rate of 0%, which meant nearly 65.9% was incompatible. We concluded that C. alismatifolia 'Siam Shadow' (Ch34) was suitable as a male parent and C. petiolata Roxb. (Cpet) was suitable as a female parent to improve the fruit-setting rates. The maximum number of seeds per fruit (45.4) was obtained when C. alismatifolia 'Chiang Mai Pink' (Ca01) was used as a female parent followed by C. attenuata Wall. ex Baker (Catt) (42.8) and C. alismatifolia 'Splash' (Ca63) (39.6) as male parents. The highest germination rate was observed for the Ca group followed by Catt and C. sparganiifolia 'Maetang Sunrise' (Csms). The germination rates of Ca accessions ranged from 58.2% (C. alismatifolia 'Siam Scarlet' (Ca06) as a male parent) to 89.3% (C. alismatifolia 'Sitone' (Ca10) as a male parent) with an average value of 74.0%. Based on the results of hybrid identification, all the individuals from the four combinations exhibited paternal-specific bands, indicating that the true hybrid rates of crossings were 100%. Our results would facilitate the interspecific hybridization and introduction of genetic variation from wild species into the cultivars in Curcuma in the future, which could be helpful in realizing the sustainable application in urban green areas.

Functional conservation and divergence of SEPALLATA-like genes in floral development in Cymbidium sinense.

Cymbidium sinense is one of the most important traditional Chinese Orchids due to its unique and highly ornamental floral organs. Although the ABCDE model for flower development is well-established in model plant species, the precise roles of these genes in C. sinense are not yet fully understood. In this study, four SEPALLATA-like genes were isolated and identified from C. sinense. CsSEP1 and CsSEP3 were grouped into the AGL9 clade, while CsSEP2 and CsSEP4 were included in the AGL2/3/4 clade. The expression pattern of CsSEP genes showed that they were significantly accumulated in reproductive tissues and expressed during flower bud development but only mildly detected or even undetected in vegetative organs. Subcellular localization revealed that CsSEP1 and CsSEP4 were localized to the nucleus, while CsSEP2 and CsSEP3 were located at the nuclear membrane. Promoter sequence analysis predicted that CsSEP genes contained a number of hormone response elements (HREs) and MADS-box binding sites. The early flowering phenotype observed in transgenic Arabidopsis plants expressing four CsSEP genes, along with the expression profiles of endogenous genes, such as SOC1, LFY, AG, FT, SEP3 and TCPs, in both transgenic Arabidopsis and C. sinense protoplasts, suggested that the CsSEP genes played a regulatory role in the flowering transition by influencing downstream genes related to flowering. However, only transgenic plants overexpressing CsSEP3 and CsSEP4 caused abnormal phenotypes of floral organs, while CsSEP1 and CsSEP2 had no effect on floral organs. Protein-protein interaction assays indicated that CsSEPs formed a protein complex with B-class CsAP3-2 and CsSOC1 proteins, affecting downstream genes to regulate floral organs and flowering time. Our findings highlighted both the functional conservation and divergence of SEPALLATA-like genes in C. sinense floral development. These results provided a valuable foundation for future studies of the molecular network underlying floral development in C. sinense.

The Integrated mRNA and miRNA Approach Reveals Potential Regulators of Flowering Time in Arundina graminifolia.

Orchids are among the most precious flowers in the world. Regulation of flowering time is one of the most important targets to enhance their ornamental value. The beauty of Arundina graminifolia is its year-round flowering, although the molecular mechanism of this flowering ability remains masked. Therefore, we performed a comprehensive assessment to integrate transcriptome and miRNA sequencing to disentangle the genetic regulation of flowering in this valuable species. Clustering analyses provided a set of molecular regulators of floral transition and floral morphogenesis. We mined candidate floral homeotic genes, including FCA, FPA, GI, FT, FLC, AP2, SOC1, SVP, GI, TCP, and CO, which were targeted by a variety of miRNAs. MiR11091 targeted the highest number of genes, including candidate regulators of phase transition and hormonal control. The conserved miR156-miR172 pathway of floral time regulation was evident in our data, and we found important targets of these miRNAs in the transcriptome. Moreover, endogenous hormone levels were determined to decipher the hormonal control of floral buds in A. graminifolia. The qRT-PCR analysis of floral and hormonal integrators validated the transcriptome expression. Therefore, miRNA-mediated mining of candidate genes with hormonal regulation forms the basis for comprehending the complex regulatory network of perpetual flowering in precious orchids. The findings of this study can do a great deal to broaden the breeding programs for flowering time manipulation of orchids.

Diversity, classification, and EST-SSR-based association analysis of caladium ornamental traits.

Caladium (Caladium × Hortulanum Birdsey) is a popular ornamental plant with a wide range of vibrant leaf color among Araceae. Even after years of breeding, creating new caladium leaf color variations is extremely difficult. Molecular marker-assisted selection is an effective approach for accelerating breeding, but few studies on the molecular markers associated with caladium traits have been performed. In the current study, 144 caladium accessions were used to examine 12 phenotypic characteristics. The coefficient of variation for four numerical characters ranged from 23.94% to 43.22%, and the Shannon-Wiener indexes for eight descriptive characters ranged from 0.13 to 1.52. Based on L*, a*, b*, C, h° values determined by a colorimeter and hierarchical cluster analysis, the leaf color can be divided into four groups: pale green, green, light pink, and red. Furthermore, 7708 new SSR loci were identified by transcriptome sequencing, and 26 SSR markers with high polymorphism and reproducibility were screened. Genetic structure, NJ clustering, and PCoA analysis revealed that 144 accessions could be divided into three clusters, with genetic structure being closely related to germplasm origin. An association analysis revealed that the SSR markers 2, 1, 1, 1, 1, and 1 were mainly associated with petiole color, main vein color, blade upperside glossiness, and C, b*, and L* of leaf color (p < 0.01). These findings will serve as a valuable reference for evaluating germplasm resources and caladium molecular marker-assisted breeding.

Complete Chloroplast Genome Sequences of Four Species in the Caladium Genus: Comparative and Phylogenetic Analyses.

Caladiums are promising colorful foliage plants due to their dazzling colors of the leaves, veins, stripes, and patches, which are often cultivated in pots or gardens as decorations. Four wild species, including C. bicolor, C. humboldtii, C. praetermissum, and C. lindenii, were employed in this study, where their chloroplast (cp) genomes were sequenced, assembled, and annotated via high-throughput sequencing. The whole cp genome size ranged from 162,776 bp to 168,888 bp, and the GC contents ranged from 35.09% to 35.91%. Compared with the single large copy (LSC) and single small copy (SSC) regions, more conserved sequences were identified in the inverted repeat regions (IR). We further analyzed the different region borders of nine species of Araceae and found the expansion or contraction of IR/SSC regions might account for the cp genome size variation. Totally, 131 genes were annotated in the cp genomes, including 86 protein-coding genes (PCGs), 37 tRNAs, and eight rRNAs. The effective number of codons (ENC) values and neutrality plot analyses provided the foundation that the natural selection pressure could greatly affect the codon preference. The GC3 content was significantly lower than that of GC1 and GC2, and codons ending with A/U had higher usage preferences. Finally, we conducted phylogenetic relationship analysis based on the chloroplast genomes of twelve species of Araceae, in which C. bicolor and C. humboldtii were grouped together, and C. lindenii was furthest from the other three Caladium species occupying a separate branch. These results will provide a basis for the identification, development, and utilization of Caladium germplasm.

The genomic and bulked segregant analysis of Curcuma alismatifolia revealed its diverse bract pigmentation.

Compared with most flowers where the showy part comprises specialized leaves (petals) directly subtending the reproductive structures, most Zingiberaceae species produce showy "flowers" through modifications of leaves (bracts) subtending the true flowers throughout an inflorescence. Curcuma alismatifolia, belonging to the Zingiberaceae family, a plant species originating from Southeast Asia, has become increasingly popular in the flower market worldwide because of its varied and esthetically pleasing bracts produced in different cultivars. Here, we present the chromosome-scale genome assembly of C. alismatifolia "Chiang Mai Pink" and explore the underlying mechanisms of bract pigmentation. Comparative genomic analysis revealed C. alismatifolia contains a residual signal of whole-genome duplication. Duplicated genes, including pigment-related genes, exhibit functional and structural differentiation resulting in diverse bract colors among C. alismatifolia cultivars. In addition, we identified the key genes that produce different colored bracts in C. alismatifolia, such as F3'5'H, DFR, ANS and several transcription factors for anthocyanin synthesis, as well as chlH and CAO in the chlorophyll synthesis pathway by conducting transcriptomic analysis, bulked segregant analysis using both DNA and RNA data, and population genomic analysis. This work provides data for understanding the mechanism of bract pigmentation and will accelerate breeding in developing novel cultivars with richly colored bracts in C. alismatifolia and related species. It is also important to understand the variation in the evolution of the Zingiberaceae family. Supplementary Information: The online version contains supplementary material available at 10.1007/s42994-022-00081-6.

Comparative chloroplast genomes and phylogenetic relationships of Aglaonema modestum and five variegated cultivars of Aglaonema.

Aglaonema, commonly called Chinese evergreens, are widely used for ornamental purposes. However, attempts to identify Aglaonema species and cultivars based on leaf morphology have been challenging. In the present study, chloroplast sequences were used to elucidate the phylogenetic relationships of cultivated Aglaonema in South China. The chloroplast genomes of one green species and five variegated cultivars of Aglaonema, Aglaonema modestum, 'Red Valentine', 'Lady Valentine', 'Hong Yan', 'Hong Jian', and 'Red Vein', were sequenced for comparative and phylogenetic analyses. The six chloroplast genomes of Aglaonema had typical quadripartite structures, comprising a large single copy (LSC) region (91,092-91,769 bp), a small single copy (SSC) region (20,816-26,501 bp), and a pair of inverted repeat (IR) regions (21,703-26,732 bp). The genomes contained 112 different genes, including 79-80 protein coding genes, 28-29 tRNAs and 4 rRNAs. The molecular structure, gene order, content, codon usage, long repeats, and simple sequence repeats (SSRs) were generally conserved among the six sequenced genomes, but the IR-SSC boundary regions were significantly different, and 'Red Vein' had a distinct long repeat number and type frequency. For comparative and phylogenetic analyses, Aglaonema costatum was included; it was obtained from the GenBank database. Single-nucleotide polymorphisms (SNPs) and insertions/deletions (indels) were determined among the seven Aglaonema genomes studied. Nine divergent hotspots were identified: trnH-GUG-CDS1_psbA, trnS-GCU_trnS-CGA-CDS1, rps4-trnT-UGU, trnF-GAA-ndhJ, petD-CDS2-rpoA, ycf1-ndhF, rps15-ycf1-D2, ccsA-ndhD, and trnY-GUA-trnE-UUC. Additionally, positive selection was found for rpl2, rps2, rps3, ycf1 and ycf2 based on the analyses of Ka/Ks ratios among 16 Araceae chloroplast genomes. The phylogenetic tree based on whole chloroplast genomes strongly supported monophyletic Aglaonema and clear relationships among Aroideae, Lasioideae, Lemnoideae, Monsteroideae, Orontioideae, Pothoideae and Zamioculcadoideae in the family Araceae. By contrast, protein coding gene phylogenies were poorly to strongly supported and incongruent with the whole chloroplast genome phylogenetic tree. This study provided valuable genome resources and helped identify Aglaonema species and cultivars.

Transcriptome Analysis Reveals Endogenous Hormone Changes during Spike Development in Phalaenopsis.

Phalaenopsis orchids are popular worldwide due to their high ornamental and economic value; the spike and inflorescence formation of their flowers could be efficiently controlled under proper conditions. In this study, transcriptomic profiles and endogenous hormone changes were investigated to better understand the spike formation of Phalaenopsis. Morphological observations revealed four spike initiation statuses (i.e., S0: the status refers to axillary buds remaining dormant in the leaf axils; S1: the status refers to the 0.5 cm-long initial spike; S2: the status refers to the 1 cm-long spike; S3: the status refers to the 3 cm-long spike) during the process of spike development, while anatomical observations revealed four related statuses of inflorescence primordium differentiation. A total of 4080 differentially expressed genes were identified based on pairwise comparisons of the transcriptomic data obtained from the S0 to S3 samples; high levels of differential gene expression were mostly observed in S1 vs. S2, followed by S0 vs. S1. Then, the contents of 12 endogenous hormones (e.g., irindole-3-acetic acid (IAA), salicylic acid (SA), abscisic acid (ABA), gibberellins, and cytokinins) were measured. The results showed that the ABA content was decreased from S0 to S1, while the gibberellic acid 1 (GA1) content exhibited an opposite trend, indicating the reduction in ABA levels combined with the increase in GA1 levels in S0 promoted the axillary bud dormancy breaking, preparing for the following spike initiation. The GA20 oxidase and ABA 8'-hydroxylase genes, which are involved in endogenous hormone metabolism and signaling pathways, displayed similar expression patterns, suggesting they were probably the key genes participating in the GA and ABA regulation. Taken together, the findings of this study indicate that GA and ABA may be the key endogenous hormones breaking the dormancy and promoting the germination of axillary buds in Phalaenopsis.

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.

Complete chloroplast genomes and comparative analyses of Hippeastrum 'milady', Hippeastrum albertii and Hippeastrum reticulatum (Amaryllidaceae).

Hippeastrum is a genus of ornamental plants with large, brightly colored flowers. Due to the very high seed-setting rate of the hybridization of Hippeastrum, the large population of hybrid progeny and the existence of superparent inheritance, it is difficult to trace the origin of the varieties collected from the market during breeding. In this study, we analyzed the chloroplast genomes of Hippeastrum 'Milady', H. alberti, and H. reticulatum using the Illumina NovaSeq sequencing platform and generated full-length sequences of 158,067, 158,067, and 158,522 bp, respectively. All three genomes had the typical tetrad structure. The large single copy, small single copy, and inverted repeat regions of H. reticulatum were observed to be respectively 277, 138, and 20 bp longer than the corresponding regions of H. 'Milady' and H. alberti. The results of comparative analysis of simple sequence repeats (SSRs), Ka/Ks ratios, codon preferences, and complete sequences of chloroplasts of these three taxa and 14 other plant species were as follows. First, the chloroplast genomes of H. 'Milady', H. alberti, and H. reticulatum contain 209, 209, and 211 SSR sites, respectively, most of which (123, 123, and 122, respectively) are single nucleotide repeats. Second, leucine, arginine, and serine are the most frequently used amino acids in the three chloroplast genomes. Third, H. 'Milady', H. alberti, and H. reticulatum are more closely related to Lycoris and Narcissus than to Allium and Agapanthus. Our results will provide information on the study of origins or relatedness of native species, and the identification of cultivars.

The Transcriptome Profiling of Flavonoids and Bibenzyls Reveals Medicinal Importance of Rare Orchid Arundina graminifolia.

Orchids are very important flowering plants that spend long juvenile phases before flowering. Along with aesthetic importance, they are rich sources of medicinal components. However, their long reproductive cycle is the major hurdle to study the medicinal efficacy. Arundina graminifolia is a rare orchid that grows fast, unlike other orchids, and this characteristic makes it an ideal plant to study the medicinal enrichment of orchids. Therefore, this study presents the identification of important medicinal components in various parts of A. graminifolia. Transcriptome analysis was performed for five stages (FD1-FD5) of flower development and four tissue types (mature flower, silique, root, and leaf) to ascertain genetic regulators of flavonoids and bibenzyls. Most of the genes showed the highest expression in roots as compared with other tissues. Weighted gene coexpression network analysis (WGCNA) was performed to identify the coexpression modules and the candidate genes involving biosynthesis pathways of these chemicals. MEyellow module contained the highly coexpressed genes. Moreover, the concentrations of phenylpropanoid, bibenzyls, and flavone were ascertained through high-performance liquid chromatography-tandem mass spectrometry (HPLC-MS/MS). Phenylpropanoid and bibenzyl were comparatively high in the leaf, while flavone showed a high concentration in the stem. The selected candidate genes [bibenzyl biosynthesis (BIBSY212), CYP84A1, CYP73A4, 4CLL7, UGT88B1, UGT73C3, anthocyanin synthase (ANS), phenylalanine ammonia-lyase (PAL), flavanone synthase FLS, and CHS8] were validated through quantitative real-time PCR (qRT-PCR). Most of these genes showed high expression in leaf and root as compared with other tissue. Therefore, the presence of bibenzyls and flavonoids in different parts of A. graminifolia and their molecular regulators can provide a quick source to decipher the medicinal efficacy of orchids.