Snapdragon (Antirrhinum majus L.), a member of the Plantaginaceae family, is an important model for plant genetics and molecular studies on plant growth and development, transposon biology and self-incompatibility. Here we report a near-complete genome assembly of A. majus cultivar JI7 (A. majus cv.JI7) comprising 510 Megabases (Mb) of genomic sequence and containing 37,714 annotated protein-coding genes. Scaffolds covering 97.12% of the assembled genome were anchored on eight chromosomes. Comparative and evolutionary analyses revealed that a whole-genome duplication event occurred in the Plantaginaceae around 46–49 million years ago (Ma). We also uncovered the genetic architectures associated with complex traits such as flower asymmetry and self-incompatibility, identifying a unique duplication of TCP family genes dated to around 46–49 Ma and reconstructing a near-complete ψS-locus of roughly 2 Mb. The genome sequence obtained in this study not only provides a representative genome sequenced from the Plantaginaceae but also brings the popular plant model system of Antirrhinum into the genomic age.

The self-incompatibility (SI) system with the broadest taxonomic distribution in angiosperms is based on multiple S-locus F-box genes (SLFs) tightly linked to an S-RNase termed type-1. Multiple SLFs collaborate to detoxify nonself S-RNases while being unable to detoxify self S-RNases. However, it is unclear how such a system evolved, because in an ancestral system with a single SLF, many nonself S-RNases would not be detoxified, giving low cross-fertilization rates. In addition, how the system has been maintained in the face of whole-genome duplications (WGDs) or lost in other lineages remains unclear. Here we show that SLFs from a broad range of species can detoxify S-RNases from Petunia with a high detoxification probability, suggestive of an ancestral feature enabling cross-fertilization and subsequently modified as additional SLFs evolved. We further show, based on its genomic signatures, that type-1 was likely maintained in many lineages, despite WGD, through deletion of duplicate S-loci. In other lineages, SI was lost either through S-locus deletions or by retaining duplications. Two deletion lineages regained SI through type-2 (Brassicaceae) or type-4 (Primulaceae), and one duplication lineage through type-3 (Papaveraceae) mechanisms. Thus, our results reveal a highly dynamic process behind the origin, maintenance, loss, and regain of SI.

Angiosperms have evolved the metabolic capacity to synthesize p-hydroxyphenyl, guaiacyl (G), and syringyl (S) lignin subunits in their cell walls to better adapt to the harsh terrestrial environment. The structural characteristics of lignin subunits are essentially determined by three cytochrome P450-catalzyed reactions. NADPH-dependent cytochrome P450 oxidoreductase (CPR) is commonly regarded as the electron carrier for P450-catalyzed reactions during monolignol biosynthesis. Here, we show that cytochrome b 5 isoform D (CB5D) is an indispensable electron shuttle protein specific for S-lignin biosynthesis. Arabidopsis (Arabidopsis thaliana) CB5D localizes to the endoplasmic reticulum membrane and physically associates with monolignol P450 enzymes. Disrupting CB5D in Arabidopsis resulted in a >60% reduction in S-lignin subunit levels but no impairment in G-lignin formation compared with the wild type, which sharply contrasts with the impaired G-and S-lignin synthesis observed after disrupting ATR2, encoding Arabidopsis CPR. The defective S-lignin synthesis in cb5d mutants was rescued by the expression of the gene encoding CB5D but not with mutant CB5D devoid of its electron shuttle properties. Disrupting ATR2 suppressed the catalytic activity of both cinnamic acid 4-hydroxylase and ferulate 5-hydroxylase (F5H), but eliminating CB5D specifically depleted the latter's activity. Therefore, CB5D functions as an obligate electron shuttle intermediate that specifically augments F5H-catalyzed reactions, thereby controlling S-lignin biosynthesis.