Diploid and tetraploid genomes of Acorus and the evolution of monocots.

Monocots are a major taxon within flowering plants, have unique morphological traits, and show an extraordinary diversity in lifestyle. To improve our understanding of monocot origin and evolution, we generate chromosome-level reference genomes of the diploid Acorus gramineus and the tetraploid Ac. calamus, the only two accepted species from the family Acoraceae, which form a sister lineage to all other monocots. Comparing the genomes of Ac. gramineus and Ac. calamus, we suggest that Ac. gramineus is not a potential diploid progenitor of Ac. calamus, and Ac. calamus is an allotetraploid with two subgenomes A, and B, presenting asymmetric evolution and B subgenome dominance. Both the diploid genome of Ac. gramineus and the subgenomes A and B of Ac. calamus show clear evidence of whole-genome duplication (WGD), but Acoraceae does not seem to share an older WGD that is shared by most other monocots. We reconstruct an ancestral monocot karyotype and gene toolkit, and discuss scenarios that explain the complex history of the Acorus genome. Our analyses show that the ancestors of monocots exhibit mosaic genomic features, likely important for that appeared in early monocot evolution, providing fundamental insights into the origin, evolution, and diversification of monocots.

Genome-wide identification of YABBY genes in three Cymbidium species and expression patterns in C. ensifolium (Orchidaceae).

Members of the YABBY gene family play significant roles in lamina development in cotyledons, floral organs, and other lateral organs. The Orchidaceae family is one of the largest angiosperm groups. Some YABBYs have been reported in Orchidaceae. However, the function of YABBY genes in Cymbidium is currently unknown. In this study, 24 YABBY genes were identified in Cymbidium ensifolium, C. goeringii, and C. sinense. We analyzed the conserved domains and motifs, the phylogenetic relationships, chromosome distribution, collinear correlation, and cis-elements of these three species. We also analyzed expression patterns of C. ensifolium and C. goeringii. Phylogenetic relationships analysis indicated that 24 YABBY genes were clustered in four groups, INO, CRC/DL, YAB2, and YAB3/FIL. For most YABBY genes, the zinc finger domain was located near the N-terminus and the helix-loop-helix domain (YABBY domain) near the C-terminus. Chromosomal location analysis results suggested that only C. goeringii YABBY has tandem repeat genes. Almost all the YABBY genes displayed corresponding one-to-one relationships in the syntenic relationships analysis. Cis-elements analysis indicated that most elements were clustered in light-responsive elements, followed by MeJA-responsive elements. Expression patterns showed that YAB2 genes have high expression in floral organs. RT-qPCR analysis showed high expression of CeYAB3 in lip, petal, and in the gynostemium. CeCRC and CeYAB2.2 were highly expressed in gynostemium. These findings provide valuable information of YABBY genes in Cymbidium species and the function in Orchidaceae.

Genomes of leafy and leafless Platanthera orchids illuminate the evolution of mycoheterotrophy.

To improve our understanding of the origin and evolution of mycoheterotrophic plants, we here present the chromosome-scale genome assemblies of two sibling orchid species: partially mycoheterotrophic Platanthera zijinensis and holomycoheterotrophic Platanthera guangdongensis. Comparative analysis shows that mycoheterotrophy is associated with increased substitution rates and gene loss, and the deletion of most photoreceptor genes and auxin transporter genes might be linked to the unique phenotypes of fully mycoheterotrophic orchids. Conversely, trehalase genes that catalyse the conversion of trehalose into glucose have expanded in most sequenced orchids, in line with the fact that the germination of orchid non-endosperm seeds needs carbohydrates from fungi during the protocorm stage. We further show that the mature plant of P. guangdongensis, different from photosynthetic orchids, keeps expressing trehalase genes to hijack trehalose from fungi. Therefore, we propose that mycoheterotrophy in mature orchids is a continuation of the protocorm stage by sustaining the expression of trehalase genes. Our results shed light on the molecular mechanism underlying initial, partial and full mycoheterotrophy.

The Cymbidium genome reveals the evolution of unique morphological traits.

The marvelously diverse Orchidaceae constitutes the largest family of angiosperms. The genus Cymbidium in Orchidaceae is well known for its unique vegetation, floral morphology, and flower scent traits. Here, a chromosome-scale assembly of the genome of Cymbidium ensifolium (Jianlan) is presented. Comparative genomic analysis showed that C. ensifolium has experienced two whole-genome duplication (WGD) events, the most recent of which was shared by all orchids, while the older event was the τ event shared by most monocots. The results of MADS-box genes analysis provided support for establishing a unique gene model of orchid flower development regulation, and flower shape mutations in C. ensifolium were shown to be associated with the abnormal expression of MADS-box genes. The most abundant floral scent components identified included methyl jasmonate, acacia alcohol and linalool, and the genes involved in the floral scent component network of C. ensifolium were determined. Furthermore, the decreased expression of photosynthesis-antennae and photosynthesis metabolic pathway genes in leaves was shown to result in colorful striped leaves, while the increased expression of MADS-box genes in leaves led to perianth-like leaves. Our results provide fundamental insights into orchid evolution and diversification.

The Euscaphis japonica genome and the evolution of malvids.

Malvids is one of the largest clades of rosids, includes 58 families and exhibits remarkable morphological and ecological diversity. Here, we report a high-quality chromosome-level genome assembly for Euscaphis japonica, an early-diverging species within malvids. Genome-based phylogenetic analysis suggests that the unstable phylogenetic position of E. japonica may result from incomplete lineage sorting and hybridization event during the diversification of the ancestral population of malvids. Euscaphis japonica experienced two polyploidization events: the ancient whole genome triplication event shared with most eudicots (commonly known as the γ event) and a more recent whole genome duplication event, unique to E. japonica. By resequencing 101 samples from 11 populations, we speculate that the temperature has led to the differentiation of the evergreen and deciduous of E. japonica and the completely different population histories of these two groups. In total, 1012 candidate positively selected genes in the evergreen were detected, some of which are involved in flower and fruit development. We found that reddening and dehiscence of the E. japonica pericarp and long fruit-hanging time promoted the reproduction of E. japonica populations, and revealed the expression patterns of genes related to fruit reddening, dehiscence and abscission. The key genes involved in pentacyclic triterpene synthesis in E. japonica were identified, and different expression patterns of these genes may contribute to pentacyclic triterpene diversification. Our work sheds light on the evolution of E. japonica and malvids, particularly on the diversification of E. japonica and the genetic basis for their fruit dehiscence and abscission.

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.

Wolfberry genomes and the evolution of Lycium (Solanaceae).

Wolfberry Lycium, an economically important genus of the Solanaceae family, contains approximately 80 species and shows a fragmented distribution pattern among the Northern and Southern Hemispheres. Although several herbaceous species of Solanaceae have been subjected to genome sequencing, thus far, no genome sequences of woody representatives have been available. Here, we sequenced the genomes of 13 perennial woody species of Lycium, with a focus on Lycium barbarum. Integration with other genomes provides clear evidence supporting a whole-genome triplication (WGT) event shared by all hitherto sequenced solanaceous plants, which occurred shortly after the divergence of Solanaceae and Convolvulaceae. We identified new gene families and gene family expansions and contractions that first appeared in Solanaceae. Based on the identification of self-incompatibility related-gene families, we inferred that hybridization hotspots are enriched for genes that might be functioning in gametophytic self-incompatibility pathways in wolfberry. Extremely low expression of LOCULE NUBER (LC) and COLORLESS NON-RIPENING (CNR) orthologous genes during Lycium fruit development and ripening processes suggests functional diversification of these two genes between Lycium and tomato. The existence of additional flowering locus C-like MADS-box genes might correlate with the perennial flowering cycle of Lycium. Differential gene expression involved in the lignin biosynthetic pathway between Lycium and tomato likely illustrates woody and herbaceous differentiation. We also provide evidence that Lycium migrated from Africa into Asia, and subsequently from Asia into North America. Our results provide functional insights into Solanaceae origins, evolution and diversification.

Genome-Wide Identification of the YABBY Gene Family in Seven Species of Magnoliids and Expression Analysis in Litsea.

The YABBY gene family, specific to seed plants, encodes a class of transcription factors in the lamina maintenance and development of lateral organs. Magnoliids are sisters to the clade-containing eudicots and monocots, which have rapidly diversified among the common ancestors of these three lineages. However, prior to this study, information on the function of the YABBY genes in magnoliids was extremely limited to the third major clades and the early diverging lineage of Mesangiospermae. In this study, the sum of 55 YABBY genes including five genes in INO, six in CRC, eight in YAB2, 22 in YAB5, and 14 in FIL clade were identified from seven magnoliid plants. Sequence analysis showed that all encoded YABBY protein sequences possess the highly conserved YABBY domain and C2C2 zinc-finger domain. Gene and protein structure analysis indicates that a certain number of exons were highly conserved and similar in the same class, and YABBY genes encode proteins of 71-392 amino acids and an open reading frame of 216-1179 bp in magnoliids. Additionally, the predicted molecular weight and isoelectric point of YABBY proteins in three species ranged from 7689.93 to 43578.13 and from 5.33 to 9.87, respectively. Meanwhile, the YABBY gene homolog expression of Litsea was detected at a temporal and spatial level during various developmental stages of leaf and reproductive tissues. This research could provide a brief overview of YABBY gene family evolution and its differential expression in magnoliids. Therefore, this comprehensive diversification analysis would provide a new insight into further understanding of the function of genes in seven magnoliids.

The Phoebe genome sheds light on the evolution of magnoliids.

Lauraceae includes the genus Phoebe, and the family is linked to the evolution of magnoliids. We sequenced the genome of Phoebe bournei Nanmu. The assembled genome size was 989.19 Mb, with a contig N50 value of 2.05 Mb. A total of 28,198 protein-coding genes were annotated in P. bournei. Whole-genome duplication (WGD) analysis showed that Lauraceae has experienced two WGD events; the older WGD event occurred just before the divergence of Lauraceae and Magnoliales, and the more recent WGD was shared by all lineages of Lauraceae. The phylogenetic tree showed that magnoliids form a sister clade to monocots and eudicots. We also identified 63 MADS-box genes, including AGL12-like genes that may be related to the regulation of P. bournei roots and FIN219-like genes encoding GH3 proteins, which are involved in photomorphogenesis. SAUR50-like genes involved in light signal-mediated pedicel or stem development were also identified. Four ATMYB46- and three PtrEPSP-homologous genes related to lignin biosynthesis were identified. These genes may be associated with the formation of straight trunks in P. bournei. Overall, the P. bournei reference genome provides insight into the origin, evolution, and diversification of Phoebe and other magnoliids.

The genome sequence of star fruit (Averrhoa carambola).

Oxalidaceae is one of the most important plant families in horticulture, and its key commercially relevant genus, Averrhoa, has diverse growth habits and fruit types. Here, we describe the assembly of a high-quality chromosome-scale genome sequence for Averrhoa carambola (star fruit). Ks distribution analysis showed that A. carambola underwent a whole-genome triplication event, i.e., the gamma event shared by most eudicots. Comparisons between A. carambola and other angiosperms also permitted the generation of Oxalidaceae gene annotations. We identified unique gene families and analyzed gene family expansion and contraction. This analysis revealed significant changes in MADS-box gene family content, which might be related to the cauliflory of A. carambola. In addition, we identified and analyzed a total of 204 nucleotide-binding site, leucine-rich repeat receptor (NLR) genes and 58 WRKY genes in the genome, which may be related to the defense response. Our results provide insights into the origin, evolution and diversification of star fruit.

Chromosome-scale assembly of the Kandelia obovata genome.

The mangrove Kandelia obovata (Rhizophoraceae) is an important coastal shelterbelt and landscape tree distributed in tropical and subtropical areas across East Asia and Southeast Asia. Herein, a chromosome-level reference genome of K. obovata based on PacBio, Illumina, and Hi-C data is reported. The high-quality assembled genome size is 177.99 Mb, with a contig N50 value of 5.74 Mb. A large number of contracted gene families and a small number of expanded gene families, as well as a small number of repeated sequences, may account for the small K. obovata genome. We found that K. obovata experienced two whole-genome polyploidization events: one whole-genome duplication shared with other Rhizophoreae and one shared with most eudicots (γ event). We confidently annotated 19,138 protein-coding genes in K. obovata and identified the MADS-box gene class and the RPW8 gene class, which might be related to flowering and resistance to powdery mildew in K. obovata and Rhizophora apiculata, respectively. The reference K. obovata genome described here will be very useful for further molecular elucidation of various traits, the breeding of this coastal shelterbelt species, and evolutionary studies with related taxa.

Transcriptome Analysis and Identification of Genes Associated with Starch Metabolism in Castanea henryi Seed (Fagaceae).

Starch is the most important form of carbohydrate storage and is the major energy reserve in some seeds, especially Castanea henryi. Seed germination is the beginning of the plant's life cycle, and starch metabolism is important for seed germination. As a complex metabolic pathway, the regulation of starch metabolism in C. henryi is still poorly understood. To explore the mechanism of starch metabolism during the germination of C. henryi, we conducted a comparative gene expression analysis at the transcriptional level using RNA-seq across four different germination stages, and analyzed the changes in the starch and soluble sugar contents. The results showed that the starch content increased in 0-10 days and decreased in 10-35 days, while the soluble sugar content continuously decreased in 0-30 days and increased in 30-35 days. We identified 49 candidate genes that may be associated with starch and sucrose metabolism. Three ADP-glucose pyrophosphorylase (AGPase) genes, two nucleotide pyrophosphatase/phosphodiesterases (NPPS) genes and three starch synthases (SS) genes may be related to starch accumulation. Quantitative real-time polymerase chain reaction (qRT-PCR) was used to validate the expression levels of these genes. Our study combined transcriptome data with physiological and biochemical data, revealing potential candidate genes that affect starch metabolism during seed germination, and provides important data about starch metabolism and seed germination in seed plants.

The complete chloroplast genome sequence of Monotropa uniflora (Ericaceae).

Monotropa uniflora is one of the representative plants of Ericaceae family, which was famous for entire translucent and 'ghostly' white. Also, unique lifestyle also attracts lots of researchers, which it obtains through fixed carbon from photosynthetic plants via a shared mycorrhizal network. In this study, the complete chloroplast (cp) genome of M. uniflora was assembled and annotated, its full-length is 26,913 bp. Plastid genome contains 31 genes, 14 protein-coding genes, 14 tRNA genes, and 3 rRNA genes. The phylogenetic analyses based on the complete chloroplast genome sequence provided solid evidence that M. uniflora has a close relationship M. odorata. The chloroplast genome presented here will help for further conservation of M. uniflora and other saprophytes.

Complete chloroplast genome of Chloranthus henryi (chloranthaceae).

The complete chloroplast genome of Chloranthus henryi, an important traditional Chinese herbal medicine, was sequenced and characterized in this study. The genome size is 157,990 bp in length with 37.3% GC content. Two inverted repeats of 26,151 bp are separated by a large single-copy region (87,148 bp), and a small single-copy region (18,569 bp). A total of 131 genes were identified, including 86 protein-coding genes, 37 tRNA genes and 8 rRNA genes. Eight plastome accessions from Chloranthales, Austrobaileyales, Nymphaeales, and Amborellales were selected to assess the phylogenetic placement of genus.

The complete chloroplast genome sequence of Euscaphis japonica (Staphyleaceae).

Euscaphis japonica is not only an ideal ornamental plant, but also a traditional medicinal plant, which is an extremely valuable species to study. We determined the complete chloroplast genome sequence for E. japonica using Illumina sequencing technology. The complete chloroplast sequence is 160,467 bp in length, including a large single-copy (LSC) region of 88,716 bp, a small single-copy (SSC) region of 18,614 bp, and a pair of invert repeats (IR) regions of 26,568 bp. Plastid genome contains 142 genes, 46 protein-coding genes, 39 tRNA genes, and eight rRNA genes. Phylogenetic analysis base on 14 chloroplast genomes indicates that E. japonica forms an isolated clade and sisters to Glycosmis-Gossypium clade with strong support.

The complete chloroplast genome sequence of Ludisia discolor from Hainan of China.

Ludisia discolor is one of the most important ornamental and medicinal orchids. To improve our understanding of the evolution of chloroplast, we resequenced complete chloroplast (cp) genome of L. discolor from Hainan, China. The cp genome sequence of L. discolor of Hainan was 153,324 bp in length, with a large single-copy (LSC) region of 82,922 bp, a small single-copy (SSC) region of 17,258 bp, and a pair of inverted repeats (IR) regions of 26,572 bp. Complete chloroplast genome contain 132 genes, there were 86 protein-coding genes, 38 tRNA genes and eight rRNA genes. The phylogenetic tree showed that L. discolor of Hainan is sister to L. discolor (unknown distributed region). Their cp genomes have same gene number but different in length of genome, indicating high conserved among them.