The genus Liriodendron belongs to the family Magnoliaceae, which resides within the magnoliids, an early diverging lineage of the Mesangiospermae. However, the phylogenetic relationship of magnoliids with eudicots and monocots has not been conclusively resolved and thus remains to be determined^1–6. Liriodendron is a relict lineage from the Tertiary with two distinct species—one East Asian (L. chinense (Hemsley) Sargent) and one eastern North American (L. tulipifera Linn)—identified as a vicariad species pair. However, the genetic divergence and evolutionary trajectories of these species remain to be elucidated at the whole-genome level⁷. Here, we report the first de novo genome assembly of a plant in the Magnoliaceae, L. chinense. Phylogenetic analyses suggest that magnoliids are sister to the clade consisting of eudicots and monocots, with rapid diversification occurring in the common ancestor of these three lineages. Analyses of population genetic structure indicate that L. chinense has diverged into two lineages—the eastern and western groups—in China. While L. tulipifera in North America is genetically positioned between the two L. chinense groups, it is closer to the eastern group. This result is consistent with phenotypic observations that suggest that the eastern and western groups of China may have diverged long ago, possibly before the intercontinental differentiation between L. chinense and L. tulipifera. Genetic diversity analyses show that L. chinense has tenfold higher genetic diversity than L. tulipifera, suggesting that the complicated regions comprising east–west-orientated mountains and the Yangtze river basin (especially near 30° N latitude) in East Asia offered more successful refugia than the south–north-orientated mountain valleys in eastern North America during the Quaternary glacial period.

Metasequoia glyptostroboides Hu et Cheng is the only species in the genus Metasequoia Miki ex Hu et Cheng, which belongs to the Cupressaceae family. There were around 10 species in the Metasequoia genus, which were widely spread across the Northern Hemisphere during the Cretaceous of the Mesozoic and in the Cenozoic. M. glyptostroboides is the only remaining representative of this genus. Here, we report the complete chloroplast (cp) genome sequence and the cp genomic features of M. glyptostroboides. The M. glyptostroboides cp genome is 131,887 bp in length, with a total of 117 genes comprised of 82 protein-coding genes, 31 tRNA genes and four rRNA genes. In this genome, 11 forward repeats, nine palindromic repeats, and 15 tandem repeats were detected. A total of 188 perfect microsatellites were detected through simple sequence repeat (SSR) analysis and these were distributed unevenly within the cp genome. Comparison of the cp genome structure and gene order to those of several other land plants indicated that a copy of the inverted repeat (IR) region, which was found to be IR region A (IRA), was lost in the M. glyptostroboides cp genome. The five most divergent and five most conserved genes were determined and further phylogenetic analysis was performed among plant species, especially for related species in conifers. Finally, phylogenetic analysis demonstrated that M. glyptostroboides is a sister species to Cryptomeria japonica (L. F.) D. Don and to Taiwania cryptomerioides Hayata. The complete cp genome sequence information of M. glyptostroboides will be great helpful for further investigations of this endemic relict woody plant and for in-depth understanding of the evolutionary history of the coniferous cp genomes, especially for the position of M. glyptostroboides in plant systematics and evolution.

The plant elicitor peptides (Peps), a family of damage/danger-associated molecular patterns (DAMPs), are perceived by two receptors, PEPR1 and PEPR2, and contribute to plant defense against pathogen attack and abiotic stress. Here we show that the Peps-PEPR signaling pathway functions in stomatal immunity by activating guard cell anion channels in Arabidopsis thaliana. The mutant plants lacking both PEPR1 and PEPR2 (pepr1 pepr2) displayed enhanced bacterial growth after being sprayed with Pseudomonas syringae pv. tomato (Pst) DC3000, but not after pathogen infiltration into leaves, implicating PEPR function in stomatal immunity. Indeed, synthetic Arabidopsis Peps (AtPeps) effectively induced stomatal closure in wild-type but not pepr1 pepr2 mutant leaves, suggesting that the AtPeps-PEPR signaling pathway triggers stomatal closure. Consistent with this finding, patch-clamp recording revealed AtPep1-induced activation of anion channels in the guard cells of wild-type but not pepr1 pepr2 mutant plants. We further identified two guard cell-expressed anion channels SLOW ANION CHANNEL 1 (SLAC1) and its homologue SLAH3, as functionally overlapping components responsible for AtPep1-induced stomatal closure. The slac1 slah3 double mutant, but not slac1 or slah3 single mutants, failed to respond to AtPep1 in stomatal closure assays. Interestingly, disruption of OPEN STOMATA 1 (OST1), an essential gene for abscisic acid-triggered stomatal closure, did not affect the AtPep1-induced anion channel activity and stomatal response. Together, these results illustrate a DAMP-triggered signaling pathway that, unlike the flagellin22-FLAGELLIN-SENSITIVE 2 pathway, triggers stomata immunity through an OST1-independent mechanism.

Knowledge of the proteome of the early gymnosperm embryo could provide important information for optimizing plant cloning procedures and for establishing platforms for research into plant development/regulation and in vitro transgenic studies. Compared with angiosperms, it is more difficult to induce somatic embryogenesis in gymnosperms; success in this endeavour could be increased, however, if proteomic information was available on the complex, dynamic, and multistage processes of gymnosperm embryogenesis in vivo. A proteomic analysis of Chinese fir seeds in six developmental stages was carried out during early embryogenesis. Proteins were extracted from seeds dissected from immature cones and separated by two-dimensional difference gel electrophoresis. Analysis with DeCyder 6.5 software revealed 136 spots that differed in kinetics of appearance. Analysis by liquid chromatography coupled to tandem mass spectrometry and MALDI-TOF mass spectrometry identified proteins represented by 71 of the spots. Functional annotation of these seed proteins revealed their involvement in programmed cell death and chromatin modification, indicating that the proteins may play a central role in determining the number of zygotic embryos generated and controlling embryo patterning and shape remodelling. The analysis also revealed other proteins involved in carbon metabolism, methionine metabolism, energy production, protein storage, synthesis and stabilization, disease/defence, the cytoskeleton, and embryo development. The comprehensive protein expression profiles generated by our study provide new insights into the complex developmental processes in the seeds of the Chinese fir.

To date, several studies have indicated a major role for microRNAs (miRNAs) in regulating plant development, but miRNA-mediated regulation of the developing somatic embryo is poorly understood, especially during early stages of somatic embryogenesis in hardwood plants. In this study, Solexa sequencing and miRNA microfluidic chips were used to discover conserved and species-specific miRNAs during somatic embryogenesis of hybrid yellow poplar (Liriodendron tulipifera×L. chinense). A total of 17,214,153 reads representing 7,421,623 distinct sequences were obtained from a short RNA library generated from small RNAs extracted from all stages of somatic embryos. Through a combination of deep sequencing and bioinformatic analyses, we discovered 83 sequences with perfect matches to known miRNAs from 33 conserved miRNA families and 273 species-specific candidate miRNAs. MicroRNA microarray results demonstrated that many conserved and species-specific miRNAs were expressed in hybrid yellow poplar embryos. In addition, the microarray also detected another 149 potential miRNAs, belonging to 29 conserved families, which were not discovered by deep sequencing analysis. The biological processes and molecular functions of the targets of these miRNAs were predicted by carrying out BLAST search against Arabidopsis thaliana GenBank sequences and then analyzing the results with Gene Ontology. Solexa sequencing and microarray hybridization were used to discover 232 candidate conserved miRNAs from 61 miRNA families and 273 candidate species-specific miRNAs in hybrid yellow poplar. In these predicted miRNAs, 64 conserved miRNAs and 177 species-specific miRNAs were detected by both sequencing and microarray hybridization. Our results suggest that miRNAs have wide-ranging characteristics and important roles during all stages of somatic embryogenesis in this economically important species.

Liriodendron tulipifera, also known as tuliptree, is a popular ornamental horticultural plant with extraordinary tulip-shaped flowers characterized by an orange band near their base. The mechanisms underlying petal band-specific pigmentation during L. tulipifera flower development are unclear. Here, we combined nontargeted and targeted metabolomics and transcriptomics to identify a pathway cascade leading to carotenoid biosynthesis that is specifically activated in the petal band. The comparative analysis of carotenoid metabolites between L. tulipifera and Liriodendron hybrids indicates that γ-carotene, a rare carotene in plants, is the most likely orange pigment responsible for the coloration of the petal band. Phenotypic and transcriptomic analyses of developing petals reveal that the band area is first predefined by the loss of green color. Later, the band is maintained by locally activating and repressing carotenoid and chlorophyll biosynthesis genes, respectively. Two rate-limiting genes of carotene biosynthesis, carotenoid isomerase (CRTISO) and epsilon lycopene cyclase (ε-LCY), encode the core enzymes responsible for petal band-specific orange pigmentation in L. tulipifera. In particular, a putative additional ε-LCY copy specific to L. tulipifera may contribute to the distinct petal coloration pattern, compared with L. chinense. Taken together, our work provides a first glimpse of the metabolome and transcriptome dynamics in tuliptree flower coloration and provides a valuable resource for flower breeding or metabolic engineering as well as for understanding flower evolution in an early woody angiosperm.