The Aquilegia genome provides insight into adaptive radiation and reveals an extraordinarily polymorphic chromosome with a unique history

Transcription Factors in Alkaloid Engineering

Plants produce a large variety of low-molecular-weight and specialized secondary compounds. Among them, nitrogen-containing alkaloids are the most biologically active and are often used in the pharmaceutical industry. Although alkaloid chemistry has been intensively investigated, characterization of alkaloid biosynthesis, including biosynthetic enzyme genes and their regulation, especially the transcription factors involved, has been relatively delayed, since only a limited number of plant species produce these specific types of alkaloids in a tissue/cell-specific or developmental-specific manner. Recent advances in molecular biology technologies, such as RNA sequencing, co-expression analysis of transcripts and metabolites, and functional characterization of genes using recombinant technology and cutting-edge technology for metabolite identification, have enabled a more detailed characterization of alkaloid pathways. Thus, transcriptional regulation of alkaloid biosynthesis by transcription factors, such as basic helix-loop-helix (bHLH), APETALA2/ethylene-responsive factor (AP2/ERF), and WRKY, is well elucidated. In addition, jasmonate signaling, an important cue in alkaloid biosynthesis, and its cascade, interaction of transcription factors, and post-transcriptional regulation are also characterized and show cell/tissue-specific or developmental regulation. Furthermore, current sequencing technology provides more information on the genome structure of alkaloid-producing plants with large and complex genomes, for genome-wide characterization. Based on the latest information, we discuss the application of transcription factors in alkaloid engineering.

Brassinosteroids regulate petal spur length in Aquilegia by controlling cell elongation

BACKGROUND AND AIMS

Aquilegia produce elongated, three-dimensional petal spurs that fill with nectar to attract pollinators. Previous studies have shown that the diversity of spur length across the Aquilegia genus is a key innovation that is tightly linked with its recent and rapid diversification into new ranges, and that evolution of increased spur lengths is achieved via anisotropic cell elongation. Previous work identified a brassinosteroid response transcription factor as being enriched in the early developing spur cup. Brassinosteroids are known to be important for cell elongation, suggesting that brassinosteroid-mediated response may be an important regulator of spur elongation and potentially a driver of spur length diversity in Aquilegia. In this study, we investigated the role of brassinosteroids in the development of the Aquilegia coerulea petal spur.

METHODS

We exogenously applied the biologically active brassinosteroid brassinolide to developing petal spurs to investigate spur growth under high hormone conditions. We used virus-induced gene silencing and gene expression experiments to understand the function of brassinosteroid-related transcription factors in A. coerulea petal spurs.

KEY RESULTS

We identified a total of three Aquilegia homologues of the BES1/BZR1 protein family and found that these genes are ubiquitously expressed in all floral tissues during development, yet, consistent with the previous RNAseq study, we found that two of these paralogues are enriched in early developing petals. Exogenously applied brassinosteroid increased petal spur length due to increased anisotropic cell elongation as well as cell division. We found that targeting of the AqBEH genes with virus-induced gene silencing resulted in shortened petals, a phenotype caused in part by a loss of cell anisotropy.

CONCLUSIONS

Collectively, our results support a role for brassinosteroids in anisotropic cell expansion in Aquilegia petal spurs and highlight the brassinosteroid pathway as a potential player in the diversification of petal spur length in Aquilegia.

Three chromosome-scale Papaver genomes reveal punctuated patchwork evolution of the morphinan and noscapine biosynthesis pathway

For millions of years, plants evolve plenty of structurally diverse secondary metabolites (SM) to support their sessile lifestyles through continuous biochemical pathway innovation. While new genes commonly drive the evolution of plant SM pathway, how a full biosynthetic pathway evolves remains poorly understood. The evolution of pathway involves recruiting new genes along the reaction cascade forwardly, backwardly, or in a patchwork manner. With three chromosome-scale Papaver genome assemblies, we here reveal whole-genome duplications (WGDs) apparently accelerate chromosomal rearrangements with a nonrandom distribution towards SM optimization. A burst of structural variants involving fusions, translocations and duplications within 7.7 million years have assembled nine genes into the benzylisoquinoline alkaloids gene cluster, following a punctuated patchwork model. Biosynthetic gene copies and their total expression matter to morphinan production. Our results demonstrate how new genes have been recruited from a WGD-induced repertoire of unregulated enzymes with promiscuous reactivities to innovate efficient metabolic pathways with spatiotemporal constraint.

Conserved and non-conserved triggers of 24-nucleotide reproductive phasiRNAs in eudicots

Small RNAs play important roles in plant growth and development by modulating expression of genes and transposons. In many flowering plant species, male reproductive organs, the anthers, produce abundant phased small interfering RNAs (phasiRNAs). Two classes of reproductive phasiRNAs are generally known, mostly from monocots: (i) pre-meiotic 21-nucleotide (nt) phasiRNAs triggered by miR2118 and (ii) meiotic 24-nt phasiRNAs triggered by miR2275. Here, we describe conserved and non-conserved triggers of 24-nt phasiRNAs in several eudicots. We found that the abundant 24-nt phasiRNAs in the basal eudicot columbine (Aquilegia coerulea) are produced by the canonical trigger miR2275, as well as by other non-canonical triggers, miR482/2118 and miR14051. These triggering microRNAs (miRNAs) are localized in microspore mother cells and tapetal cells of meiotic and post-meiotic stage anthers. Furthermore, we identified a lineage-specific trigger (miR11308) of 24-nt phasiRNAs and an expanded number of 24-PHAS loci in wild strawberry (Fragaria vesca). We validated the presence of the miR2275-derived 24-nt phasiRNA pathway in rose (Rosa chinensis). Finally, we evaluated all eudicots that have been validated for the presence of 24-nt phasiRNAs as possible model systems in which to study the biogenesis and function of 24-nt phasiRNAs. We conclude that columbine (Aquilegia coerulea) would be a strong model because of its extensive number of 24-PHAS loci and its diversity of trigger miRNAs.

Divergence in the Aquilegia ecalcarata complex is correlated with geography and climate oscillations: Evidence from plastid genome data

Quaternary climate oscillations and geographical heterogeneity play important roles in determining species and genetic diversity distribution patterns, but how these factors affect the migration and differentiation of East Asian plants species at the population level remains poorly understood. The Aquilegia ecalcarata complex, a group that originated in the Late Tertiary and is widely distributed throughout East Asia, displays high genetic variation that is suitable for studying elaborate phylogeographic patterns and demographic history related to the impact of Quaternary climate and geography. We used plastid genome data from 322 individuals in 60 populations of the A. ecalcarata complex to thoroughly explore the impact of Quaternary climate oscillations and geography on the phylogeographic patterns and demographic history of the A. ecalcarata complex through a series of phylogenetic, divergence time estimation, and demographic history analyses. The dry, cold climate and frequent climate oscillations that occurred during the early Pleistocene and the Mid-Pleistocene transition led to the differentiation of the A. ecalcarata complex, which was isolated in various areas. Geographically, the A. ecalcarata complex can be divided into Eastern and Western Clades and five subclades, which conform to the divergence of the East Asian flora. Our results clearly show the impact of Quaternary climate and geography on evolutionary history at the population level. These findings promote the understanding of the relationship between plant genetic differentiation and climate and geographical factors of East Asia at the population level.

Structural diversity, bioactivities, and biosynthesis of natural diterpenoid alkaloids

Covering: 2009 to 2018. Diterpenoid alkaloids, originating from the amination of natural tetracyclic diterpenes, are a diverse class of compounds having complex structural features with many stereocenters. The important pharmacological activities and structural complexity of the diterpenoid alkaloids have long interested scientists due to their medicinal uses, infamous toxicity, and unique biosynthesis. Since 2009, 373 diterpenoid alkaloids, assigned to 46 skeletons, have been isolated and identified from plants mostly in the Ranunculaceae family. The names, classes, molecular weight, molecular formula, NMR data, and plant sources of these diterpene alkaloids are collated here. This review will be a detailed update of the naturally occurring diterpene alkaloids reported from the plant kingdom from 2009-2018, providing an in-depth discussion of their diversity, biological activities, pharmacokinetics, toxicity, application, evolution, and biosynthesis.

Pre-meiotic 21-nucleotide reproductive phasiRNAs emerged in seed plants and diversified in flowering plants

Plant small RNAs are important regulatory elements that fine-tune gene expression and maintain genome integrity by silencing transposons. Reproductive organs of monocots produce abundant phased, small interfering RNAs (phasiRNAs). The 21-nt reproductive phasiRNAs triggered by miR2118 are highly enriched in pre-meiotic anthers, and have been found in multiple eudicot species, in contrast with prior reports of monocot specificity. The 24-nt reproductive phasiRNAs are triggered by miR2275, and are highly enriched during meiosis in many angiosperms. Here, we report the widespread presence of the 21-nt reproductive phasiRNA pathway in eudicots including canonical and non-canonical microRNA (miRNA) triggers of this pathway. In eudicots, these 21-nt phasiRNAs are enriched in pre-meiotic stages, a spatiotemporal distribution consistent with that of monocots and suggesting a role in anther development. Although this pathway is apparently absent in well-studied eudicot families including the Brassicaceae, Solanaceae and Fabaceae, our work in eudicots supports an earlier singular finding in spruce, a gymnosperm, indicating that the pathway of 21-nt reproductive phasiRNAs emerged in seed plants and was lost in some lineages.

Genetic architecture of floral traits in bee- and hummingbird-pollinated sister species of Aquilegia (columbine)

Interactions with animal pollinators have helped shape the stunning diversity of flower morphologies across the angiosperms. A common evolutionary consequence of these interactions is that some flowers have converged on suites of traits, or pollination syndromes, that attract and reward specific pollinator groups. Determining the genetic basis of these floral pollination syndromes can help us understand the processes that contributed to the diversification of the angiosperms. Here, we characterize the genetic architecture of a bee-to-hummingbird pollination shift in Aquilegia (columbine) using QTL mapping of 17 floral traits encompassing color, nectar composition, and organ morphology. In this system, we find that the genetic architectures underlying differences in floral color are quite complex, and we identify several likely candidate genes involved in anthocyanin and carotenoid floral pigmentation. Most morphological and nectar traits also have complex genetic underpinnings; however, one of the key floral morphological phenotypes, nectar spur curvature, is shaped by a single locus of large effect.

Analysis of the Coptis chinensis genome reveals the diversification of protoberberine-type alkaloids

Chinese goldthread (Coptis chinensis Franch.), a member of the Ranunculales, represents an important early-diverging eudicot lineage with diverse medicinal applications. Here, we present a high-quality chromosome-scale genome assembly and annotation of C. chinensis. Phylogenetic and comparative genomic analyses reveal the phylogenetic placement of this species and identify a single round of ancient whole-genome duplication (WGD) shared by the Ranunculaceae. We characterize genes involved in the biosynthesis of protoberberine-type alkaloids in C. chinensis. In particular, local genomic tandem duplications contribute to member amplification of a Ranunculales clade-specific gene family of the cytochrome P450 (CYP) 719. The functional versatility of a key CYP719 gene that encodes the (S)-canadine synthase enzyme involved in the berberine biosynthesis pathway may play critical roles in the diversification of other berberine-related alkaloids in C. chinensis. Our study provides insights into the genomic landscape of early-diverging eudicots and provides a valuable model genome for genetic and applied studies of Ranunculales.

The chromosome-level reference genome of Coptis chinensis provides insights into genomic evolution and berberine biosynthesis

Coptis chinensis Franch, a perennial herb, is mainly distributed in southeastern China. The rhizome of C. chinensis has been used as a traditional medicine for more than 2000 years in China and many other Asian countries. The pharmacological activities of C. chinensis have been validated by research. Here, we present a de novo high-quality genome of C. chinensis with a chromosome-level genome of ~958.20 Mb, a contig N50 of 1.58 Mb, and a scaffold N50 of 4.53 Mb. We found that the relatively large genome size of C. chinensis was caused by the amplification of long terminal repeat (LTR) retrotransposons. In addition, a whole-genome duplication event in ancestral Ranunculales was discovered. Comparative genomic analysis revealed that the tyrosine decarboxylase (TYDC) and (S)-norcoclaurine synthase (NCS) genes were expanded and that the aspartate aminotransferase gene (ASP5) was positively selected in the berberine metabolic pathway. Expression level and HPLC analyses showed that the berberine content was highest in the roots of C. chinensis in the third and fourth years. The chromosome-level reference genome of C. chinensis provides important genomic data for molecular-assisted breeding and active ingredient biosynthesis.

Datura genome reveals duplications of psychoactive alkaloid biosynthetic genes and high mutation rate following tissue culture

BACKGROUND

Datura stramonium (Jimsonweed) is a medicinally and pharmaceutically important plant in the nightshade family (Solanaceae) known for its production of various toxic, hallucinogenic, and therapeutic tropane alkaloids. Recently, we published a tissue-culture based transformation protocol for D. stramonium that enables more thorough functional genomics studies of this plant. However, the tissue culture process can lead to undesirable phenotypic and genomic consequences independent of the transgene used. Here, we have assembled and annotated a draft genome of D. stramonium with a focus on tropane alkaloid biosynthetic genes. We then use mRNA sequencing and genome resequencing of transformants to characterize changes following tissue culture.

RESULTS

Our draft assembly conforms to the expected 2 gigabasepair haploid genome size of this plant and achieved a BUSCO score of 94.7% complete, single-copy genes. The repetitive content of the genome is 61%, with Gypsy-type retrotransposons accounting for half of this. Our gene annotation estimates the number of protein-coding genes at 52,149 and shows evidence of duplications in two key alkaloid biosynthetic genes, tropinone reductase I and hyoscyamine 6 β-hydroxylase. Following tissue culture, we detected only 186 differentially expressed genes, but were unable to correlate these changes in expression with either polymorphisms from resequencing or positional effects of transposons.

CONCLUSIONS

We have assembled, annotated, and characterized the first draft genome for this important model plant species. Using this resource, we show duplications of genes leading to the synthesis of the medicinally important alkaloid, scopolamine. Our results also demonstrate that following tissue culture, mutation rates of transformed plants are quite high (1.16 × 10- 3 mutations per site), but do not have a drastic impact on gene expression.

Comparative chloroplast genomes and phylogenetic analysis of Aquilegia

PREMISE

Aquilegia is an ideal taxon for studying the evolution of adaptive radiation. Current phylogenies of Aquilegia based on different molecular markers are inconsistent, and therefore a clear and accurate phylogeny remains uncertain. Analyzing the chloroplast genome, with its simple structure and low recombination rate, may help solve this problem.

METHODS

Next-generation sequencing data were generated or downloaded for Aquilegia species, enabling their chloroplast genomes to be assembled. The assemblies were used to estimate the genome characteristics and infer the phylogeny of Aquilegia.

RESULTS

In this study, chloroplast genome sequences were assembled for Aquilegia species distributed across Asia, North America, and Europe. Three of the genes analyzed (petG, rpl36, and atpB) were shown to be under positive selection and may be related to adaptation. The phylogenetic tree of Aquilegia showed that its member species formed two clades with high support, North American and European species, with the Asian species being paraphyletic; A. parviflora and A. amurensis clustered with the North American species, while the remaining Asian species were found in the European clade. In addition, A. oxysepala var. kansuensis should be considered as a separate species rather than a variety.

DISCUSSION

The complete chloroplast genomes of these Aquilegia species provide new insights into the reconstruction of the phylogeny of related species and contribute to the further study of this genus.

Genes and genome editing tools for breeding desirable phenotypes in ornamentals

We review the main genes underlying commercial traits in cut flower species and critically discuss the possibility to apply genome editing approaches to produce novel variation and phenotypes. Promoting flowering and flower longevity as well as creating novelty in flower structure, colour range and fragrances are major objectives of ornamental plant breeding. The novel genome editing techniques add new possibilities to study gene function and breed new varieties. The implementation of such techniques, however, relies on detailed information about structure and function of genomes and genes. Moreover, improved protocols for efficient delivery of editing reagents are required. Recent results of the application of genome editing techniques to elite ornamental crops are discussed in this review. Enabling technologies and genomic resources are reviewed in relation to the implementation of such approaches. Availability of the main gene sequences, underlying commercial traits and in vitro transformation protocols are provided for the world's best-selling cut flowers, namely rose, lily, chrysanthemum, lisianthus, tulip, gerbera, freesia, alstroemeria, carnation and hydrangea. Results obtained so far are described and their implications for the improvement of flowering, flower architecture, colour, scent and shelf-life are discussed.

A multilayered cross-species analysis of GRAS transcription factors uncovered their functional networks in plant adaptation to the environment

Introduction

Environmental stress is both a major force of natural selection and a prime factor affecting crop qualities and yields. The impact of the GRAS [gibberellic acid-insensitive (GAI), repressor of GA1-3 mutant (RGA), and scarecrow (SCR)] family on plant development and the potential to resist environmental stress needs much emphasis.

Objectives

This study aims to investigate the evolution, expansion, and adaptive mechanisms of GRASs of important representative plants during polyploidization.

Methods

We explored the evolutionary characteristics of GRASs in 15 representative plant species by systematic biological analysis of the genome, transcriptome, metabolite, protein complex map and phenotype.

Results

The GRAS family was systematically identified from 15 representative plant species of scientific and agricultural importance. The detection of gene duplication types of GRASs in all species showed that the widespread expansion of GRASs in these species was mainly contributed by polyploidization events. Evolutionary analysis reveals that most species experience independent genome-wide duplication (WGD) events and that interspecies GRAS functions may be broadly conserved. Polyploidy-related Chenopodium quinoa GRASs (CqGRASs) and Arabidopsis thaliana GRASs (AtGRASs) formed robust networks with flavonoid pathways by crosstalk with auxin and photosynthetic pathways. Furthermore, Arabidopsis thaliana population transcriptomes and the 1000 Plants (OneKP) project confirmed that GRASs are components of flavonoid biosynthesis, which enables plants to adapt to the environment by promoting flavonoid accumulation. More importantly, the GRASs of important species that may potentially improve important agronomic traits were mapped through TAIR and RARGE-II publicly available phenotypic data. Determining protein interactions and target genes contributes to determining GRAS functions.

Conclusion

The results of this study suggest that polyploidy-related GRASs in multiple species may be a target for improving plant growth, development, and environmental adaptation.

Comparative analysis using the draft genome sequence of California poppy (Eschscholzia californica) for exploring the candidate genes involved in benzylisoquinoline alkaloid biosynthesis

Genome characterization of California poppy (Eschscholzia californica cv. “Hitoezaki”), which produces pharmaceutically important benzylisoquinoline alkaloids (BIAs), was carried out using the draft genome sequence. The numbers of tRNA and rRNA genes were close to those of the other plant species tested, whereas the frequency of repetitive sequences was distinct from those species. Comparison of the predicted genes with those of Amborella trichopoda, Nelumbo nucifera, Solanum lycopersicum, and Arabidopsis thaliana, and analyses of gene ontology and Kyoto Encyclopedia of Genes and Genomes pathway indicated that the enzyme genes involved in BIA biosynthesis were highly enriched in the California poppy genome. Further comparative analysis using the genome information of Papaver somniferum and Aquilegia coerulea, both BIA-producing plants, revealed that many genes encoding BIA biosynthetic enzymes, transcription factors, transporters, and candidate proteins, possibly related to BIA biosynthesis, were specifically distributed in these plant species.

Chromosome-level genome assembly of Ophiorrhiza pumila reveals the evolution of camptothecin biosynthesis

Plant genomes remain highly fragmented and are often characterized by hundreds to thousands of assembly gaps. Here, we report chromosome-level reference and phased genome assembly of Ophiorrhiza pumila, a camptothecin-producing medicinal plant, through an ordered multi-scaffolding and experimental validation approach. With 21 assembly gaps and a contig N50 of 18.49 Mb, Ophiorrhiza genome is one of the most complete plant genomes assembled to date. We also report 273 nitrogen-containing metabolites, including diverse monoterpene indole alkaloids (MIAs). A comparative genomics approach identifies strictosidine biogenesis as the origin of MIA evolution. The emergence of strictosidine biosynthesis-catalyzing enzymes precede downstream enzymes' evolution post γ whole-genome triplication, which occurred approximately 110 Mya in O. pumila, and before the whole-genome duplication in Camptotheca acuminata identified here. Combining comparative genome analysis, multi-omics analysis, and metabolic gene-cluster analysis, we propose a working model for MIA evolution, and a pangenome for MIA biosynthesis, which will help in establishing a sustainable supply of camptothecin.

The RGF/GLV/CLEL Family of Short Peptides Evolved Through Lineage-Specific Losses and Diversification and Yet Conserves Its Signaling Role Between Vascular Plants and Bryophytes

Short secreted plant peptides act as key signaling molecules and control a plethora of developmental and physiological processes. The ROOT GROWTH FACTOR (RGF)/GOLVEN (GLV)/CLE-Like (CLEL) family of peptides was discovered to be involved in root development in Arabidopsis thaliana. In contrast to active research efforts, which have been revealing receptors and downstream signaling components, little attention has been paid to evolutionary processes that shaped the RGF signaling system as we know it in angiosperms today. As a first step toward understanding how RGF signaling emerged and evolved, this study aimed to elucidate the phylogenetic distribution and functional conservation of RGF-like sequences. Using publicly available, genome and transcriptome data, RGF-like sequences were searched in 27 liverworts, 22 mosses, 8 hornworts, 23 lycophytes, 23 ferns, 38 gymnosperms, and 8 angiosperms. This led to the identification of more than four hundreds of RGF-like sequences in all major extant land plant lineages except for hornworts. Sequence comparisons within and between taxonomic groups identified lineage-specific characters. Notably, one of the two major RGF subgroups, represented by A. thaliana RGF6/GLV1/CLEL6, was found only in vascular plants. This subgroup, therefore, likely emerged in a common ancestor of vascular plants after its divergence from bryophytes. In bryophytes, our results infer independent losses of RGF-like sequences in mosses and hornworts. On the other hand, a single, highly similar RGF-like sequence is conserved in liverworts, including Marchantia polymorpha, a genetically tractable model species. When constitutively expressed, the M. polymorpha RGF-like sequence (MpRGF) affected plant development and growth both in A. thaliana and M. polymorpha. This suggests that MpRGF can exert known RGF-like effects and that MpRGF is under transcriptional control so that its potent activities are precisely controlled. These data suggest that RGFs are conserved as signaling molecules in both vascular plants and bryophytes and that lineage-specific diversification has increased sequence variations of RGFs. All together, our findings form a basis for further studies into RGF peptides and their receptors, which will contribute to our understandings of how peptide signaling pathways evolve.

Genomic and transcriptomic resources for candidate gene discovery in the Ranunculids

PREMISE

Multiple transitions from insect to wind pollination are associated with polyploidy and unisexual flowers in Thalictrum (Ranunculaceae), yet the underlying genetics remains unknown. We generated a draft genome of Thalictrum thalictroides, a representative of a clade with ancestral floral traits (diploid, hermaphrodite, and insect pollinated) and a model for functional studies. Floral transcriptomes of T. thalictroides and of wind-pollinated, andromonoecious T. hernandezii are presented as a resource to facilitate candidate gene discovery in flowers with different sexual and pollination systems.

METHODS

A draft genome of T. thalictroides and two floral transcriptomes of T. thalictroides and T. hernandezii were obtained from HiSeq 2000 Illumina sequencing and de novo assembly.

RESULTS

The T. thalictroides de novo draft genome assembly consisted of 44,860 contigs (N50 = 12,761 bp, 243 Mbp total length) and contained 84.5% conserved embryophyte single-copy genes. Floral transcriptomes contained representatives of most eukaryotic core genes, and most of their genes formed orthogroups.

DISCUSSION

To validate the utility of these resources, potential candidate genes were identified for the different floral morphologies using stepwise data set comparisons. Single-copy gene analysis and simple sequence repeat markers were also generated as a resource for population-level and phylogenetic studies.

The Evolution of Chromosome Numbers: Mechanistic Models and Experimental Approaches

Chromosome numbers have been widely used to describe the most fundamental genomic attribute of an organism or a lineage. Although providing strong phylogenetic signal, chromosome numbers vary remarkably among eukaryotes at all levels of taxonomic resolution. Changes in chromosome numbers regularly serve as indication of major genomic events, most notably polyploidy and dysploidy. Here, we review recent advancements in our ability to make inferences regarding historical events that led to alterations in the number of chromosomes of a lineage. We first describe the mechanistic processes underlying changes in chromosome numbers, focusing on structural chromosomal rearrangements. Then, we focus on experimental procedures, encompassing comparative cytogenomics and genomics approaches, and on computational methodologies that are based on explicit models of chromosome-number evolution. Together, these tools offer valuable predictions regarding historical events that have changed chromosome numbers and genome structures, as well as their phylogenetic and temporal placements.

POPOVICH, encoding a C2H2 zinc-finger transcription factor, plays a central role in the development of a key innovation, floral nectar spurs, in Aquilegia

The evolution of novel features, such as eyes or wings, that allow organisms to exploit their environment in new ways can lead to increased diversification rates. Therefore, understanding the genetic and developmental mechanisms involved in the origin of these key innovations has long been of interest to evolutionary biologists. In flowering plants, floral nectar spurs are a prime example of a key innovation, with the independent evolution of spurs associated with increased diversification rates in multiple angiosperm lineages due to their ability to promote reproductive isolation via pollinator specialization. As none of the traditional plant model taxa have nectar spurs, little is known about the genetic and developmental basis of this trait. Nectar spurs are a defining feature of the columbine genus Aquilegia (Ranunculaceae), a lineage that has experienced a relatively recent and rapid radiation. We use a combination of genetic mapping, gene expression analyses, and functional assays to identify a gene crucial for nectar spur development, POPOVICH (POP), which encodes a C2H2 zinc-finger transcription factor. POP plays a central role in regulating cell proliferation in the Aquilegia petal during the early phase (phase I) of spur development and also appears to be necessary for the subsequent development of nectaries. The identification of POP opens up numerous avenues for continued scientific exploration, including further elucidating of the genetic pathway of which it is a part, determining its role in the initial evolution of the Aquilegia nectar spur, and examining its potential role in the subsequent evolution of diverse spur morphologies across the genus.

Identification of the target genes of AqAPETALA3-3 (AqAP3-3) in Aquilegia coerulea (Ranunculaceae) helps understand the molecular bases of the conserved and nonconserved features of petals

Identification and comparison of the conserved and variable downstream genes of floral organ identity regulators are critical to understanding the mechanisms underlying the commonalities and peculiarities of floral organs. Yet, because of the lack of studies in nonmodel species, a general picture of the regulatory evolution between floral organ identity genes and their targets is still lacking. Here, by conducting extensive chromatin immunoprecipitation followed by high-throughput sequencing (ChIP-seq), electrophoretic mobility shift assay and bioinformatic analyses, we identify and predict the target genes of a petal identity gene, AqAPETALA3-3 (AqAP3-3), in Aquilegia coerulea (Ranunculaceae) and compare them with those of its counterpart in Arabidopsis thaliana, AP3. In total, 7049 direct target genes are identified for AqAP3-3, of which 2394 are highly confident and 1085 are shared with AP3. Gene Ontology enrichment analyses further indicate that conserved targets are largely involved in the formation of identity-related features, whereas nonconserved targets are mostly required for the formation of species-specific features. These results not only help understand the molecular bases of the conserved and nonconserved features of petals, but also pave the way to studying the regulatory evolution between floral organ identity genes and their targets.

A chromosome-scale reference genome of Aquilegia oxysepala var. kansuensis

The genus Aquilegia (Ranunculaceae) has been cultivated as ornamental and medicinal plants for centuries. With petal spurs of strikingly diverse size and shape, Aquilegia has also been recognized as an excellent system for evolutionary studies. Pollinator-mediated selection for longer spurs is believed to have shaped the evolution of this genus, especially the North American taxa. Recently, however, an opposite evolutionary trend was reported in an Asian lineage, where multiple origins of mini- or even nonspurred morphs have occurred. Interesting as it is, the lack of genomic resources has limited our ability to decipher the molecular and evolutionary mechanisms underlying spur reduction in this special lineage. Using long-read sequencing (PacBio Sequel), in combination with optical maps (BioNano DLS) and Hi-C, we assembled a high-quality reference genome of A. oxysepala var. kansuensis, a sister species to the nonspurred taxon. The final assembly is approximately 293.2 Mb, 94.6% (277.4 Mb) of which has been anchored to 7 pseudochromosomes. A total of 25,571 protein-coding genes were predicted, with 97.2% being functionally annotated. When comparing this genome with that of A. coerulea, we detected a large rearrangement between Chr1 and Chr4, which might have caused the Chr4 of A. oxysepala var. kansuensis to partly deviate from the "decaying" path that was taken before the split of Aquilegia and Semiaquilegia. This high-quality reference genome is fundamental to further investigations on the development and evolution of petal spurs and provides a strong foundation for the breeding of new horticultural Aquilegia cultivars.

Polyploidization events shaped the transcription factor repertoires in legumes (Fabaceae)

Transcription factors (TFs) are essential for plant growth and development. Several legumes (e.g. soybean) are rich sources of protein and oil and have great economic importance. Here we report a phylogenomic analysis of TF families in legumes and their potential association with important traits (e.g. nitrogen fixation). We used TF DNA-binding domains to systematically screen the genomes of 15 leguminous and five non-leguminous species. Transcription factor orthologous groups (OGs) were used to estimate OG sizes in ancestral nodes using a gene birth-death model, which allowed the identification of lineage-specific expansions. The OG analysis and rate of synonymous substitutions show that major TF expansions are strongly associated with whole-genome duplication (WGD) events in the legume (approximately 58 million years ago) and Glycine (approximately 13 million years ago) lineages, which account for a large fraction of the Phaseolus vulgaris and Glycine max TF repertoires. Of the 3407 G. max TFs, 1808 and 676 have homeologs within single syntenic regions in Phaseolus vulgaris and Vitis vinifera, respectively. We found a trend for TFs expanded in legumes to be preferentially transcribed in roots and nodules, supporting their recruitment early in the evolution of nodulation in the legume clade. Some families also showed count differences between G. max and the wild soybean Glycine soja, including genes located within important quantitative trait loci. Our findings strongly support the roles of two WGDs in shaping the TF repertoires in the legume and Glycine lineages, and these are probably related to important aspects of legume and soybean biology.

Evolutionary Understanding of Aquaporin Transport System in the Basal Eudicot Model Species Aquilegia coerulea

Aquaporins (AQPs) play a pivotal role in the cellular transport of water and many other small solutes, influencing many physiological and developmental processes in plants. In the present study, extensive bioinformatics analysis of AQPs was performed in Aquilegia coerulea L., a model species belonging to basal eudicots, with a particular focus on understanding the AQPs role in the developing petal nectar spur. A total of 29 AQPs were identified in Aquilegia, and their phylogenetic analysis performed with previously reported AQPs from rice, poplar and Arabidopsis depicted five distinct subfamilies of AQPs. Interestingly, comparative analysis revealed the loss of an uncharacterized intrinsic protein II (XIP-II) group in Aquilegia. The absence of the entire XIP subfamily has been reported in several previous studies, however, the loss of a single clade within the XIP family has not been characterized. Furthermore, protein structure analysis of AQPs was performed to understand pore diversity, which is helpful for the prediction of solute specificity. Similarly, an AQP AqcNIP2-1 was identified in Aquilegia, predicted as a silicon influx transporter based on the presence of features such as the G-S-G-R aromatic arginine selectivity filter, the spacing between asparagine-proline-alanine (NPA) motifs and pore morphology. RNA-seq analysis showed a high expression of tonoplast intrinsic proteins (TIPs) and plasma membrane intrinsic proteins (PIPs) in the developing petal spur. The results presented here will be helpful in understanding the AQP evolution in Aquilegia and their expression regulation, particularly during floral development.

Over 100 Million Years of Enzyme Evolution Underpinning the Production of Morphine in the Papaveraceae Family of Flowering Plants

Phylogenomic analysis of whole genome sequences of five benzylisoquinoline alkaloid (BIA)-producing species from the Ranunculales and Proteales orders of flowering plants revealed the sequence and timing of evolutionary events leading to the diversification of these compounds. (S)-Reticuline is a pivotal intermediate in the synthesis of many BIAs and our analyses revealed parallel evolution between the two orders, which diverged ∼122 million years ago (MYA). Berberine is present in species across the entire Ranunculales, and we found co-evolution of genes essential for production of the protoberberine class. The benzophenanthridine class, which includes the antimicrobial compound sanguinarine, is specific to the Papaveraceae family of Ranunculales, and biosynthetic genes emerged after the split with the Ranunculaceae family ∼110 MYA but before the split of the three Papaveraceae species used in this study at ∼77 MYA. The phthalideisoquinoline noscapine and morphinan class of BIAs are exclusive to the opium poppy lineage. K_s estimation of paralogous pairs indicates that morphine biosynthesis evolved more recently than 18 MYA in the Papaver genus. In the preceding 100 million years gene duplication, neofunctionalization and recruitment of additional enzyme classes, combined with gene clustering, gene fusion, and gene amplification, resulted in emergence of medicinally valuable BIAs including morphine and noscapine.

Nanopore sequencing-based genome assembly and evolutionary genomics of circum-basmati rice

Background

The circum-basmati group of cultivated Asian rice (Oryza sativa) contains many iconic varieties and is widespread in the Indian subcontinent. Despite its economic and cultural importance, a high-quality reference genome is currently lacking, and the group’s evolutionary history is not fully resolved. To address these gaps, we use long-read nanopore sequencing and assemble the genomes of two circum-basmati rice varieties.

Results

We generate two high-quality, chromosome-level reference genomes that represent the 12 chromosomes of Oryza. The assemblies show a contig N50 of 6.32 Mb and 10.53 Mb for Basmati 334 and Dom Sufid, respectively. Using our highly contiguous assemblies, we characterize structural variations segregating across circum-basmati genomes. We discover repeat expansions not observed in japonica—the rice group most closely related to circum-basmati—as well as the presence and absence variants of over 20 Mb, one of which is a circum-basmati-specific deletion of a gene regulating awn length. We further detect strong evidence of admixture between the circum-basmati and circum-aus groups. This gene flow has its greatest effect on chromosome 10, causing both structural variation and single-nucleotide polymorphism to deviate from genome-wide history. Lastly, population genomic analysis of 78 circum-basmati varieties shows three major geographically structured genetic groups: Bhutan/Nepal, India/Bangladesh/Myanmar, and Iran/Pakistan.

Conclusion

The availability of high-quality reference genomes allows functional and evolutionary genomic analyses providing genome-wide evidence for gene flow between circum-aus and circum-basmati, describes the nature of circum-basmati structural variation, and reveals the presence/absence variation in this important and iconic rice variety group.

The complete chloroplast genome of Aquilegia barnebyi, a basal eudicot species

Aquilegia barnebyi, belonging to the genus Aquilegia (Ranunculaceae), is a member of basal eudicot species. In this study, we obtained the complete chloroplast (cp) genome of A. barnebyi. The genome size is 161,954 bp with a GC content of 38.98%. A total of 113 unique genes including 79 protein-coding genes, 30 tRNA genes, four rRNA genes were annotated. The large single-copy region and small single-copy region contains 91,250 bp and 17,359 bp, respectively. The inverted repeat regions are 26,671 bp in length. The phyologenetic analysis indicated that A. barnebyi had a close relationship with A. coerulea. And four species in genus Aquilegia formed a monophyletic group with high support value. The availability of A. barnebyi cp genomic resources will greatly helpful for taxonomy, phylogeny and conservation genetic studies of basal eudicot specie.

A Phylogenomic Analysis of the Floral Transcriptomes of Sexually Deceptive and Rewarding European Orchids, Ophrys and Gymnadenia

The orchids (Orchidaceae) constitute one of the largest and most diverse families of flowering plants. They have evolved a great variety of adaptations to achieve pollination by a diverse group of pollinators. Many orchids reward their pollinators, typically with nectar, but the family is also well-known for employing deceptive pollination strategies in which there is no reward for the pollinator, in the most extreme case by mimicking sexual signals of pollinators. In the European flora, two examples of these different pollination strategies are the sexually deceptive genus Ophrys and the rewarding genus Gymnadenia, which differ in their level of pollinator specialization; Ophrys is typically pollinated by pseudo-copulation of males of a single insect species, whilst Gymnadenia attracts a broad range of floral visitors. Here, we present and describe the annotated floral transcriptome of Ophrys iricolor, an Andrena-pollinated representative of the genus Ophrys that is widespread throughout the Aegean. Furthermore, we present additional floral transcriptomes of both sexually deceptive and rewarding orchids, specifically the deceptive Ophrys insectifera, Ophrys aymoninii, and an updated floral transcriptome of Ophrys sphegodes, as well as the floral transcriptomes of the rewarding orchids Gymnadenia conopsea, Gymnadenia densiflora, Gymnadenia odoratissima, and Gymnadenia rhellicani (syn. Nigritella rhellicani). Comparisons of these novel floral transcriptomes reveal few annotation differences between deceptive and rewarding orchids. Since together, these transcriptomes provide a representative sample of the genus-wide taxonomic diversity within Ophrys and Gymnadenia (Orchidoideae: Orchidinae), we employ a phylogenomic approach to address open questions of phylogenetic relationships within the genera. Specifically, this includes the controversial placement of O. insectifera within the Ophrys phylogeny and the placement of "Nigritella"-type morphologies within the phylogeny of Gymnadenia. Whereas in Gymnadenia, several conflicting topologies are supported by a similar number of gene trees, a majority of Ophrys gene topologies clearly supports a placement of O. insectifera as sister to a clade containing O. sphegodes.

The Aquilegia genome reveals a hybrid origin of core eudicots

Background

Whole-genome duplications (WGDs) have dominated the evolutionary history of plants. One consequence of WGD is a dramatic restructuring of the genome as it undergoes diploidization, a process under which deletions and rearrangements of various sizes scramble the genetic material, leading to a repacking of the genome and eventual return to diploidy. Here, we investigate the history of WGD in the columbine genus Aquilegia, a basal eudicot, and use it to illuminate the origins of the core eudicots.

Results

Within-genome synteny confirms that columbines are ancient tetraploids, and comparison with the grape genome reveals that this tetraploidy appears to be shared with the core eudicots. Thus, the ancient gamma hexaploidy found in all core eudicots must have involved a two-step process: first, tetraploidy in the ancestry of all eudicots, then hexaploidy in the ancestry of core eudicots. Furthermore, the precise pattern of synteny sharing suggests that the latter involved allopolyploidization and that core eudicots thus have a hybrid origin.

Conclusions

Novel analyses of synteny sharing together with the well-preserved structure of the columbine genome reveal that the gamma hexaploidy at the root of core eudicots is likely a result of hybridization between a tetraploid and a diploid species.

Chromosome Painting Facilitates Anchoring Reference Genome Sequence to Chromosomes In Situ and Integrated Karyotyping in Banana (Musa Spp.)

Oligo painting FISH was established to identify all chromosomes in banana (Musa spp.) and to anchor pseudomolecules of reference genome sequence of Musa acuminata spp. malaccensis "DH Pahang" to individual chromosomes in situ. A total of 19 chromosome/chromosome-arm specific oligo painting probes were developed and were shown to be suitable for molecular cytogenetic studies in genus Musa. For the first time, molecular karyotypes of diploid M. acuminata spp. malaccensis (A genome), M. balbisiana (B genome), and M. schizocarpa (S genome) from the Eumusa section of Musa, which contributed to the evolution of edible banana cultivars, were established. This was achieved after a combined use of oligo painting probes and a set of previously developed banana cytogenetic markers. The density of oligo painting probes was sufficient to study chromosomal rearrangements on mitotic as well as on meiotic pachytene chromosomes. This advance will enable comparative FISH mapping and identification of chromosomal translocations which accompanied genome evolution and speciation in the family Musaceae.

Homologs of LEAFY and UNUSUAL FLORAL ORGANS Promote the Transition From Inflorescence to Floral Meristem Identity in the Cymose Aquilegia coerulea

Homologs of the transcription factor LEAFY (LFY) and the F-box family member UNUSUAL FLORAL ORGANS (UFO) have been found to promote floral meristem identity across diverse dicot model systems. The lower eudicot model Aquilegia produces cymose inflorescences that are independently evolved from the well-studied cymose models Petunia and tomato. We have previously characterized the expression pattern of the Aquilegia homolog AqLFY but in the current study, we add expression data on the two UFO homologs, AqUFO1 and 2, and conduct virus-induced gene silencing of all the loci. Down-regulation of AqLFY or AqUFO1 and 2 does not eliminate floral meristem identity but, instead, causes the transition from inflorescence to floral identity to become gradual rather than discrete. Inflorescences in down-regulated plants generate several nodes of bract/sepal chimeras and, once floral development does commence, flowers initiate several whorls of sepals before finally producing the wildtype floral whorls. In addition, silencing of AqUFO1/2 appears to specifically impact petal identity and/or the initiation of petal and stamen whorls. In general, however, there is no evidence for an essential role of AqLFY or AqUFO1/2 in transcriptional activation of the B or C gene homologs. These findings highlight differences between deeply divergent dicot lineages in the functional conservation of the floral meristem identity program.

The story of promiscuous crucifers: origin and genome evolution of an invasive species, Cardamine occulta (Brassicaceae), and its relatives

BACKGROUND AND AIMS

Cardamine occulta (Brassicaceae) is an octoploid weedy species (2n = 8x = 64) originated in Eastern Asia. It has been introduced to other continents including Europe and considered to be an invasive species. Despite its wide distribution, the polyploid origin of C. occulta remained unexplored. The feasibility of comparative chromosome painting (CCP) in crucifers allowed us to elucidate the origin and genome evolution in Cardamine species. We aimed to investigate the genome structure of C. occulta in comparison with its tetraploid (2n = 4x = 32, C. kokaiensis and C. scutata) and octoploid (2n = 8x = 64, C. dentipetala) relatives.

METHODS

Genomic in situ hybridization (GISH) and large-scale CCP were applied to uncover the parental genomes and chromosome composition of the investigated Cardamine species.

KEY RESULTS

All investigated species descended from a common ancestral Cardamine genome (n = 8), structurally resembling the Ancestral Crucifer Karyotype (n = 8), but differentiated by a translocation between chromosomes AK6 and AK8. Allotetraploid C. scutata originated by hybridization between two diploid species, C. parviflora and C. amara (2n = 2x = 16). By contrast, C. kokaiensis has an autotetraploid origin from a parental genome related to C. parviflora. Interestingly, octoploid C. occulta probably originated through hybridization between the tetraploids C. scutata and C. kokaiensis. The octoploid genome of C. dentipetala probably originated from C. scutata via autopolyploidization. Except for five species-specific centromere repositionings and one pericentric inversion post-dating the polyploidization events, the parental subgenomes remained stable in the tetra- and octoploids.

CONCLUSIONS

Comparative genome structure, origin and evolutionary history was reconstructed in C. occulta and related species. For the first time, whole-genome cytogenomic maps were established for octoploid plants. Post-polyploid evolution in Asian Cardamine polyploids has not been associated with descending dysploidy and intergenomic rearrangements. The combination of different parental (sub)genomes adapted to distinct habitats provides an evolutionary advantage to newly formed polyploids by occupying new ecological niches.

Developmental and Molecular Changes Underlying the Vernalization-Induced Transition to Flowering in Aquilegia coerulea (James)

Reproductive success in plants is dependent on many factors but the precise timing of flowering is certainly among the most crucial. Perennial plants often have a vernalization or over-wintering requirement in order to successfully flower in the spring. The shoot apical meristem undergoes drastic developmental and molecular changes as it transitions into inflorescence meristem (IM) identity, which then gives rise to floral meristems (FMs). In this study, we have examined the developmental and gene expression changes underlying the transition from the vegetative to reproductive phases in the basal eudicot Aquilegia coerulea, which has evolved a vernalization response independently relative to other established model systems. Results from both our histology and scanning electron studies demonstrate that developmental changes in the meristem occur gradually during the third and fourth weeks of vernalization. Based on RNAseq data and cluster analysis, several known flowering time loci, including AqFT and AqFL1, exhibit dramatic changes in expression during the fourth week. Further consideration of candidate gene homologs as well as unexpected loci of interest creates a framework in which we can begin to explore the genetic basis of the flowering time transition in Aquilegia.

Fluorescence in situ hybridization in plants: recent developments and future applications

Fluorescence in situ hybridization (FISH) was developed more than 30 years ago and has been the most paradigm-changing technique in cytogenetic research. FISH has been used to answer questions related to structure, mutation, and evolution of not only individual chromosomes but also entire genomes. FISH has served as an important tool for chromosome identification in many plant species. This review intends to summarize and discuss key technical development and applications of FISH in plants since 2006. The most significant recent advance of FISH is the development and application of probes based on synthetic oligonucleotides (oligos). Oligos specific to a repetitive DNA sequence, to a specific chromosomal region, or to an entire chromosome can be computationally identified, synthesized in parallel, and fluorescently labeled. Oligo probes designed from conserved DNA sequences from one species can be used among genetically related species, allowing comparative cytogenetic mapping of these species. The advances with synthetic oligo probes will significantly expand the applications of FISH especially in non-model plant species. Recent achievements and future applications of FISH and oligo-FISH are discussed.

Comparative transcriptomics of early petal development across four diverse species of Aquilegia reveal few genes consistently associated with nectar spur development

BACKGROUND

Petal nectar spurs, which facilitate pollination through animal attraction and pollen placement, represent a key innovation promoting diversification in the genus Aquilegia (Ranunculaceae). Identifying the genetic components that contribute to the development of these three-dimensional structures will inform our understanding of the number and types of genetic changes that are involved in the evolution of novel traits. In a prior study, gene expression between two regions of developing petals, the laminar blade and the spur cup, was compared at two developmental stages in the horticultural variety A. coerulea 'Origami'. Several hundred genes were differentially expressed (DE) between the blade and spur at both developmental stages. In order to narrow in on a set of genes crucial to early spur formation, the current study uses RNA sequencing (RNAseq) to conduct comparative expression analyses of petals from five developmental stages between four Aquilegia species, three with morphologically variable nectar spurs, A. sibirica, A. formosa, and A. chrysantha, and one that lacks nectar spurs, A. ecalcarata.

RESULTS

Petal morphology differed increasingly between taxa across the developmental stages assessed, with petals from all four taxa being indistinguishable pre-spur formation at developmental stage 1 (DS1) and highly differentiated by developmental stage 5 (DS5). In all four taxa, genes involved in mitosis were down-regulated over the course of the assessed developmental stages, however, many genes involved in mitotic processes remained expressed at higher levels later in development in the spurred taxa. A total of 690 genes were identified that were consistently DE between the spurred taxa and A. ecalcarata at all five developmental stages. By comparing these genes with those identified as DE between spur and blade tissue in A. coerulea 'Origami', a set of only 35 genes was identified that shows consistent DE between petal samples containing spur tissue versus those without spur tissue.

CONCLUSIONS

The results of this study suggest that expression differences in very few loci are associated with the presence and absence of spurs. In general, it appears that the spurless petals of A. ecalcarata cease cell divisions and enter the cell differentiation phase at an earlier developmental time point than those that produce spurs. This much more tractable list of 35 candidates genes will greatly facilitate targeted functional studies to assess the genetic control and evolution of petal spurs in Aquilegia.

Rapid Divergence Followed by Adaptation to Contrasting Ecological Niches of Two Closely Related Columbine Species Aquilegia japonica and A. oxysepala

Elucidating the mechanisms underlying the genetic divergence between closely related species is crucial to understanding the origin and evolution of biodiversity. The genus Aquilegia L. has undergone rapid adaptive radiation, generating about 70 well-recognized species that are specialized to distinct habitats and pollinators. In this study, to address the underlying evolutionary mechanisms that drive the genetic divergence, we analyzed the whole genomes of two ecologically isolated Aquilegia species, A. oxysepala and A. japonica as well as their putative hybrid. Our comparative genomic analyses reveal that while the two species diverged only recently and experienced recurrent gene flow, a high level of genetic divergence is observed in their nuclear genomes. In particular, candidate genomic regions that show signature of selection differ dramatically between the two species. Given that the splitting time of the two species is broadly matched with the decrease in effective population sizes, we propose that allopatric isolation together with natural selection have preceded the interspecific gene flow in the process of speciation. The observed high genetic divergence is likely an outcome of combined effects of natural selection, genetic drift and divergent sorting of ancestral polymorphisms. Our study provides a genome-wide view of how genetic divergence has evolved between closely related species.

Fluorescence in situ hybridization in plants: recent developments and future applications

Fluorescence in situ hybridization (FISH) was developed more than 30 years ago and has been the most paradigm-changing technique in cytogenetic research. FISH has been used to answer questions related to structure, mutation, and evolution of not only individual chromosomes but also entire genomes. FISH has served as an important tool for chromosome identification in many plant species. This review intends to summarize and discuss key technical development and applications of FISH in plants since 2006. The most significant recent advance of FISH is the development and application of probes based on synthetic oligonucleotides (oligos). Oligos specific to a repetitive DNA sequence, to a specific chromosomal region, or to an entire chromosome can be computationally identified, synthesized in parallel, and fluorescently labeled. Oligo probes designed from conserved DNA sequences from one species can be used among genetically related species, allowing comparative cytogenetic mapping of these species. The advances with synthetic oligo probes will significantly expand the applications of FISH especially in non-model plant species. Recent achievements and future applications of FISH and oligo-FISH are discussed.

On the origin of orphan hybrids between Aquilegia formosa and Aquilegia flavescens

We report the investigation of an Aquilegia flavescens × A. formosa population in British Columbia that is disjunct from its parents-the latter species is present locally but ecologically separated, while the former is entirely absent. To confirm hybridity, we used multivariate analysis of floral characters of field-sampled populations to ordinate phenotypes of putative hybrids in relation to those of the parental species. Microsatellite genotypes at 11 loci from 72 parental-type and putative hybrid individuals were analysed to assess evidence for admixture. Maternally inherited plastid sequences were analysed to infer the direction of hybridization and test hypotheses on the origin of the orphan hybrid population. Plants from the orphan hybrid population are on average intermediate between typical A. formosa and A. flavescens for most phenotypes examined and show evidence of genetic admixture. This population lies beyond the range of A. flavescens, but within the range of A. formosa. No pure A. flavescens individuals were observed in the vicinity, nor is this species known to occur within 200 km of the site. The hybrids share a plastid haplotype with local A. formosa populations. Alternative explanations for this pattern are evaluated. While we cannot rule out long-distance pollen dispersal followed by proliferation of hybrid genotypes, we consider the spread of an A. formosa plastid during genetic swamping of a historical A. flavescens population to be more parsimonious.

Genome structure and evolution of Antirrhinum majus L

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.

Genome structure and evolution of Antirrhinum majus L

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.

Genome sequence of Malania oleifera, a tree with great value for nervonic acid production

BACKGROUND

Malania oleifera, a member of the Olacaceae family, is an IUCN red listed tree, endemic and restricted to the Karst region of southwest China. This tree's seed is valued for its high content of precious fatty acids (especially nervonic acid). However, studies on its genetic makeup and fatty acid biogenesis are severely hampered by a lack of molecular and genetic tools.

FINDINGS

We generated 51 Gb and 135 Gb of raw DNA sequences, using Pacific Biosciences (PacBio) single-molecule real-time and 10× Genomics sequencing, respectively. A final genome assembly, with a scaffold N50 size of 4.65 Mb and a total length of 1.51 Gb, was obtained by primary assembly based on PacBio long reads plus scaffolding with 10× Genomics reads. Identified repeats constituted ∼82% of the genome, and 24,064 protein-coding genes were predicted with high support. The genome has low heterozygosity and shows no evidence for recent whole genome duplication. Metabolic pathway genes relating to the accumulation of long-chain fatty acid were identified and studied in detail.

CONCLUSIONS

Here, we provide the first genome assembly and gene annotation for M. oleifera. The availability of these resources will be of great importance for conservation biology and for the functional genomics of nervonic acid biosynthesis.

Chromosome Painting Facilitates Anchoring Reference Genome Sequence to Chromosomes In Situ and Integrated Karyotyping in Banana (Musa Spp.)

Oligo painting FISH was established to identify all chromosomes in banana (Musa spp.) and to anchor pseudomolecules of reference genome sequence of Musa acuminata spp. malaccensis "DH Pahang" to individual chromosomes in situ. A total of 19 chromosome/chromosome-arm specific oligo painting probes were developed and were shown to be suitable for molecular cytogenetic studies in genus Musa. For the first time, molecular karyotypes of diploid M. acuminata spp. malaccensis (A genome), M. balbisiana (B genome), and M. schizocarpa (S genome) from the Eumusa section of Musa, which contributed to the evolution of edible banana cultivars, were established. This was achieved after a combined use of oligo painting probes and a set of previously developed banana cytogenetic markers. The density of oligo painting probes was sufficient to study chromosomal rearrangements on mitotic as well as on meiotic pachytene chromosomes. This advance will enable comparative FISH mapping and identification of chromosomal translocations which accompanied genome evolution and speciation in the family Musaceae.

The Aquilegia genome provides insight into adaptive radiation and reveals an extraordinarily polymorphic chromosome with a unique history

The columbine genus Aquilegia is a classic example of an adaptive radiation, involving a wide variety of pollinators and habitats. Here we present the genome assembly of A. coerulea 'Goldsmith', complemented by high-coverage sequencing data from 10 wild species covering the world-wide distribution. Our analyses reveal extensive allele sharing among species and demonstrate that introgression and selection played a role in the Aquilegia radiation. We also present the remarkable discovery that the evolutionary history of an entire chromosome differs from that of the rest of the genome - a phenomenon that we do not fully understand, but which highlights the need to consider chromosomes in an evolutionary context.

Transcriptome Profiling of Two Ornamental and Medicinal Papaver Herbs

The Papaver spp. (Papaver rhoeas (Corn poppy) and Papaver nudicaule (Iceland poppy)) genera are ornamental and medicinal plants that are used for the isolation of alkaloid drugs. In this study, we generated 700 Mb of transcriptome sequences with the PacBio platform. They were assembled into 120,926 contigs, and 1185 (82.2%) of the benchmarking universal single-copy orthologs (BUSCO) core genes were completely present in our assembled transcriptome. Furthermore, using 128 Gb of Illumina sequences, the transcript expression was assessed at three stages of Papaver plant development (30, 60, and 90 days), from which we identified 137 differentially expressed transcripts. Furthermore, three co-occurrence heat maps are generated from 51 different plant genomes along with the Papaver transcriptome, i.e., secondary metabolite biosynthesis, isoquinoline alkaloid biosynthesis (BIA) pathway, and cytochrome. Sixty-nine transcripts in the BIA pathway along with 22 different alkaloids (quantified with LC-QTOF-MS/MS) were mapped into the BIA KEGG map (map00950). Finally, we identified 39 full-length cytochrome transcripts and compared them with other genomes. Collectively, this transcriptome data, along with the expression and quantitative metabolite profiles, provides an initial recording of secondary metabolites and their expression related to Papaver plant development. Moreover, these profiles could help to further detail the functional characterization of the various secondary metabolite biosynthesis and Papaver plant development associated problems.

Chromosome painting in meiosis reveals pairing of specific chromosomes in polyploid Solanum species

Analysis of chromosome pairing has been an important tool to assess the genetic similarity of homologous and homoeologous chromosomes in polyploids. However, it is technically challenging to monitor the pairing of specific chromosomes in polyploid species, especially for plant species with a large number of small chromosomes. We developed oligonucleotide-based painting probes for four different potato chromosomes. We demonstrate that these probes are robust enough to monitor a single chromosome throughout the prophase I of meiosis in polyploid Solanum species. Cultivated potato (Solanum tuberosum, 2n = 4x = 48) is an autotetraploid. We demonstrate that the four copies of each potato chromosome pair as a quadrivalent in 66-78% of the meiotic cells at the pachytene stage. Solanum demissum (2n = 6x = 72) is a hexaploid and has been controversial regarding its nature as an autopolyploid or allopolyploid. Interestingly, no hexavalent pairing was observed in meiosis. Instead, we observed three independent bivalents in 83-98% of the meiotic cells at late diakinesis and early metaphase I for the four chromosomes. These results suggest that S. demissum has evolved into a cytologically stable state with predominantly bivalent pairing in meiosis.