Return of the Lemnaceae: duckweed as a model plant system in the genomics and postgenomics era.

The aquatic Lemnaceae family, commonly called duckweed, comprises some of the smallest and fastest growing angiosperms known on Earth. Their tiny size, rapid growth by clonal propagation, and facile uptake of labeled compounds from the media were attractive features that made them a well-known model for plant biology from 1950 to 1990. Interest in duckweed has steadily regained momentum over the past decade, driven in part by the growing need to identify alternative plants from traditional agricultural crops that can help tackle urgent societal challenges, such as climate change and rapid population expansion. Propelled by rapid advances in genomic technologies, recent studies with duckweed again highlight the potential of these small plants to enable discoveries in diverse fields from ecology to chronobiology. Building on established community resources, duckweed is reemerging as a platform to study plant processes at the systems level and to translate knowledge gained for field deployment to address some of society's pressing needs. This review details the anatomy, development, physiology, and molecular characteristics of the Lemnaceae to introduce them to the broader plant research community. We highlight recent research enabled by Lemnaceae to demonstrate how these plants can be used for quantitative studies of complex processes and for revealing potentially novel strategies in plant defense and genome maintenance.

Natural rubber reduces herbivory and alters the microbiome below ground.

Laticifers are hypothesized to mediate both plant-herbivore and plant-microbe interactions. However, there is little evidence for this dual function. We investigated whether the major constituent of natural rubber, cis-1,4-polyisoprene, a phylogenetically widespread and economically important latex polymer, alters plant resistance and the root microbiome of the Russian dandelion (Taraxacum koksaghyz) under attack of a root herbivore, the larva of the May cockchafer (Melolontha melolontha). Rubber-depleted transgenic plants lost more shoot and root biomass upon herbivory than normal rubber content near-isogenic lines. Melolontha melolontha preferred to feed on artificial diet supplemented with rubber-depleted rather than normal rubber content latex. Likewise, adding purified cis-1,4-polyisoprene in ecologically relevant concentrations to diet deterred larval feeding and reduced larval weight gain. Metagenomics and metabarcoding revealed that abolishing biosynthesis of natural rubber alters the structure but not the diversity of the rhizosphere and root microbiota (ecto- and endophytes) and that these changes depended on M. melolontha damage. However, the assumption that rubber reduces microbial colonization or pathogen load is contradicted by four lines of evidence. Taken together, our data demonstrate that natural rubber biosynthesis reduces herbivory and alters the plant microbiota, which highlights the role of plant-specialized metabolites and secretory structures in shaping multitrophic interactions.

Evolution of a Novel and Adaptive Floral Scent in Wild Tobacco.

Many plants emit diverse floral scents that mediate plant-environment interactions and attain reproductive success. However, how plants evolve novel and adaptive biosynthetic pathways for floral volatiles remains unclear. Here, we show that in the wild tobacco, Nicotiana attenuata, a dominant species-specific floral volatile (benzyl acetone, BA) that attracts pollinators and deters florivore is synthesized by phenylalanine ammonia-lyase 4 (NaPAL4), isoflavone reductase 3 (NaIFR3), and chalcone synthase 3 (NaCHAL3). Transient expression of NaFIR3 alone in N. attenuata leaves is sufficient and necessary for ectopic foliar BA emissions, and coexpressing NaIFR3 with NaPAL4 and NaCHAL3 increased the BA emission levels. Independent changes in transcription of NaPAL4 and NaCHAL3 contributed to intraspecific variations of floral BA emission. However, among species, the gain of expression of NaIFR3 resulted in the biosynthesis of BA, which was only found in N. attenuata. This study suggests that novel metabolic pathways associated with adaptation can arise via reconfigurations of gene expression.

Allelic differences of clustered terpene synthases contribute to correlated intraspecific variation of floral and herbivory-induced volatiles in a wild tobacco.

Plant volatile emissions can recruit predators of herbivores for indirect defense and attract pollinators to aid in pollination. Although volatiles involved in defense and pollinator attraction are primarily emitted from leaves and flowers, respectively, they will co-evolve if their underlying genetic basis is intrinsically linked, due either to pleiotropy or to genetic linkage. However, direct evidence of co-evolving defense and floral traits is scarce. We characterized intraspecific variation of herbivory-induced plant volatiles (HIPVs), the key components of indirect defense against herbivores, and floral volatiles in wild tobacco Nicotiana attenuata. We found that variation of (E)-ß-ocimene and (E)-α-bergamotene contributed to the correlated changes in HIPVs and floral volatiles among N. attenuata natural accessions. Intraspecific variations of (E)-ß-ocimene and (E)-α-bergamotene emissions resulted from allelic variation of two genetically co-localized terpene synthase genes, NaTPS25 and NaTPS38, respectively. Analyzing haplotypes of NaTPS25 and NaTPS38 revealed that allelic variations of NaTPS25 and NaTPS38 resulted in correlated changes of (E)-ß-ocimene and (E)-α-bergamotene emission in HIPVs and floral volatiles in N. attenuata. Together, these results provide evidence that pleiotropy and genetic linkage result in correlated changes in defenses and floral signals in natural populations, and the evolution of plant volatiles is probably under diffuse selection.

Wild tobacco genomes reveal the evolution of nicotine biosynthesis.

Nicotine, the signature alkaloid of Nicotiana species responsible for the addictive properties of human tobacco smoking, functions as a defensive neurotoxin against attacking herbivores. However, the evolution of the genetic features that contributed to the assembly of the nicotine biosynthetic pathway remains unknown. We sequenced and assembled genomes of two wild tobaccos, Nicotiana attenuata (2.5 Gb) and Nicotiana obtusifolia (1.5 Gb), two ecological models for investigating adaptive traits in nature. We show that after the Solanaceae whole-genome triplication event, a repertoire of rapidly expanding transposable elements (TEs) bloated these Nicotiana genomes, promoted expression divergences among duplicated genes, and contributed to the evolution of herbivory-induced signaling and defenses, including nicotine biosynthesis. The biosynthetic machinery that allows for nicotine synthesis in the roots evolved from the stepwise duplications of two ancient primary metabolic pathways: the polyamine and nicotinamide adenine dinucleotide (NAD) pathways. In contrast to the duplication of the polyamine pathway that is shared among several solanaceous genera producing polyamine-derived tropane alkaloids, we found that lineage-specific duplications within the NAD pathway and the evolution of root-specific expression of the duplicated Solanaceae-specific ethylene response factor that activates the expression of all nicotine biosynthetic genes resulted in the innovative and efficient production of nicotine in the genus Nicotiana Transcription factor binding motifs derived from TEs may have contributed to the coexpression of nicotine biosynthetic pathway genes and coordinated the metabolic flux. Together, these results provide evidence that TEs and gene duplications facilitated the emergence of a key metabolic innovation relevant to plant fitness.

Efficient genetic transformation and CRISPR/Cas9-mediated genome editing in Lemna aequinoctialis.

The fast growth, ease of metabolic labelling and potential for feedstock and biofuels production make duckweeds not only an attractive model system for understanding plant biology, but also a potential future crop. However, current duckweed research is constrained by the lack of efficient genetic manipulation tools. Here, we report a case study on genome editing in a duckweed species, Lemna aequinoctialis, using a fast and efficient transformation and CRISPR/Cas9 tool. By optimizing currently available transformation protocols, we reduced the duration time of Agrobacterium-mediated transformation to 5-6 weeks with a success rate of over 94%. Based on the optimized transformation protocol, we generated 15 (14.3% success rate) biallelic LaPDS mutants that showed albino phenotype using a CRISPR/Cas9 system. Investigations on CRISPR/Cas9-mediated mutation spectrum among mutated L. aequinoctialis showed that most of mutations were short insertions and deletions. This study presents the first example of CRISPR/Cas9-mediated genome editing in duckweeds, which will open new research avenues in using duckweeds for both basic and applied research.

Cathepsin L induced PC-12 cell apoptosis via activation of B-Myb and regulation of cell cycle proteins.

Cathepsin L (CTSL), a cysteine protease, is responsible for the degradation of a variety of proteins. It is known to participate in neuronal apoptosis associated with abnormal cell cycle. However, the mechanisms underlying CTSL-induced cell apoptosis remain largely unclear. We reported here that rotenone caused an activation of CTSL expression in PC-12 cells, while knockdown of CTSL by small interfering RNAs or its inhibitor reduced the rotenone-induced cell cycle arrest and apoptosis. Moreover, elevation of CTSL and increased-apoptosis were accompanied by induction of B-Myb, a crucial cell cycle regulator. We found that B-Myb was increased in rotenone-treated PC-12 cells and knockdown of B-Myb ameliorated rotenone-stimulated cell apoptosis. Further analysis demonstrated that CTSL influenced the expression of B-Myb as suppression of CTSL activity led to a decreased B-Myb expression, whereas overexpression of CTSL resulted in B-Myb induction. Reduction of B-Myb in CTSL-overexpressing cells revealed that regulation of cell cycle-related proteins, including cyclin A and cyclin B1, through CTSL was mediated by the transcription factor B-Myb. In addition, we demonstrated that the B-Myb target, Bim, and its regulator, Egr-1, which was also associated with CTSL closely, were both involved in rotenone-induced apoptosis in PC-12 cells. Our data not only revealed the role of CTSL in rotenone-induced neuronal apoptosis, but also indicated the involvement of B-Myb in CTSL-related cell cycle regulation.

SINS/CNS/GNSS Integrated Navigation Based on an Improved Federated Sage-Husa Adaptive Filter.

Among the methods of the multi-source navigation filter, as a distributed method, the federated filter has a small calculation amount with Gaussian state noise, and it is easy to achieve global optimization. However, when the state noise is time-varying or its initial estimation is not accurate, there will be a big difference with the true value in the result of the federated filter. For the systems with time-varying noise, adaptive filter is widely used for its remarkable advantages. Therefore, this paper proposes a federated Sage-Husa adaptive filter for multi-source navigation systems with time-varying or mis-estimated state noise. Because both the federated and the adaptive principles are different in updating the covariance of the state noise, it is required to weight the two updating methods to obtain a combined method with stability and adaptability. In addition, according to the characteristics of the system, the weighting coefficient is formed by the exponential function. This federated adaptive filter is applied to the SINS/CNS/GNSS integrated navigation, and the simulation results show that this method is effective.

SINS/Landmark Integrated Navigation Based on Landmark Attitude Determination.

Based on the situation that the traditional SINS (strapdown inertial navigation system)/CNS (celestial navigation system) integrated navigation system fails to realize all-day and all-weather navigation, this paper proposes a SINS/Landmark integrated navigation method based on landmark attitude determination to solve this problem. This integrated navigation system takes SINS as the basic scheme and uses landmark navigation to correct the error of SINS. The way of the attitude determination is to use the landmark information photographed by the landmark camera to complete feature matching. The principle of the landmark navigation and the process of attitude determination are discussed, and the feasibility of landmark attitude determination is analyzed, including the orthogonality of the attitude transform matrix, as well as the influences of the factors such as quantity and geometric position of landmarks. On this basis, the paper constructs the equations of the SINS/Landmark integrated navigation system, testifies the effectiveness of landmark attitude determination on the integrated navigation by Kalman filter, and improves the navigation precision of the system.

O-Acyl Sugars Protect a Wild Tobacco from Both Native Fungal Pathogens and a Specialist Herbivore.

O-Acyl sugars (O-AS) are abundant trichome-specific metabolites that function as indirect defenses against herbivores of the wild tobacco Nicotiana attenuata; whether they also function as generalized direct defenses against herbivores and pathogens remains unknown. We characterized natural variation in O-AS among 26 accessions and examined their influence on two native fungal pathogens, Fusarium brachygibbosum U4 and Alternaria sp. U10, and the specialist herbivore Manduca sexta At least 15 different O-AS structures belonging to three classes were found in N. attenuata leaves. A 3-fold quantitative variation in total leaf O-AS was found among the natural accessions. Experiments with natural accessions and crosses between high- and low-O-AS accessions revealed that total O-AS levels were associated with resistance against herbivores and pathogens. Removing O-AS from the leaf surface increased M. sexta growth rate and plant fungal susceptibility. O-AS supplementation in artificial diets and germination medium reduced M. sexta growth and fungal spore germination, respectively. Finally, silencing the expression of a putative branched-chain α-ketoacid dehydrogenase E1 ß-subunit-encoding gene (NaBCKDE1B) in the trichomes reduced total leaf O-AS by 20% to 30% and increased susceptibility to Fusarium pathogens. We conclude that O-AS function as direct defenses to protect plants from attack by both native pathogenic fungi and a specialist herbivore and infer that their diversification is likely shaped by the functional interactions among these biotic stresses.

Nicotiana attenuata Data Hub (NaDH): an integrative platform for exploring genomic, transcriptomic and metabolomic data in wild tobacco.

BACKGROUND: Nicotiana attenuata (coyote tobacco) is an ecological model for studying plant-environment interactions and plant gene function under real-world conditions. During the last decade, large amounts of genomic, transcriptomic and metabolomic data have been generated with this plant which has provided new insights into how native plants interact with herbivores, pollinators and microbes. However, an integrative and open access platform that allows for the efficient mining of these -omics data remained unavailable until now. DESCRIPTION: We present the Nicotiana attenuata Data Hub (NaDH) as a centralized platform for integrating and visualizing genomic, phylogenomic, transcriptomic and metabolomic data in N. attenuata. The NaDH currently hosts collections of predicted protein coding sequences of 11 plant species, including two recently sequenced Nicotiana species, and their functional annotations, 222 microarray datasets from 10 different experiments, a transcriptomic atlas based on 20 RNA-seq expression profiles and a metabolomic atlas based on 895 metabolite spectra analyzed by mass spectrometry. We implemented several visualization tools, including a modified version of the Electronic Fluorescent Pictograph (eFP) browser, co-expression networks and the Interactive Tree Of Life (iTOL) for studying gene expression divergence among duplicated homologous. In addition, the NaDH allows researchers to query phylogenetic trees of 16,305 gene families and provides tools for analyzing their evolutionary history. Furthermore, we also implemented tools to identify co-expressed genes and metabolites, which can be used for predicting the functions of genes. Using the transcription factor NaMYB8 as an example, we illustrate that the tools and data in NaDH can facilitate identification of candidate genes involved in the biosynthesis of specialized metabolites. CONCLUSION: The NaDH provides interactive visualization and data analysis tools that integrate the expression and evolutionary history of genes in Nicotiana, which can facilitate rapid gene discovery and comparative genomic analysis. Because N. attenuata shares many genome-wide features with other Nicotiana species including cultivated tobacco, and hence NaDH can be a resource for exploring the function and evolution of genes in Nicotiana species in general. The NaDH can be accessed at: http://nadh.ice.mpg.de/ .

Evidence of an evolutionary hourglass pattern in herbivory-induced transcriptomic responses.

Herbivory-induced defenses are specific and activated in plants when elicitors, frequently found in the herbivores' oral secretions, are introduced into wounds during attack. While complex signaling cascades are known to be involved, it remains largely unclear how natural selection has shaped the evolution of these induced defenses. We analyzed herbivory-induced transcriptomic responses in wild tobacco, Nicotiana attenuata, using a phylotranscriptomic approach that measures the origin and sequence divergence of herbivory-induced genes. Highly conserved and evolutionarily ancient genes of primary metabolism were activated at intermediate time points (2-6 h) after elicitation, while less constrained and young genes associated with defense signaling and biosynthesis of specialized metabolites were activated at early (before 2 h) and late (after 6 h) stages of the induced response, respectively - a pattern resembling the evolutionary hourglass pattern observed during embryogenesis in animals and the developmental process in plants and fungi. The hourglass patterns found in herbivory-induced defense responses and developmental process are both likely to be a result of signaling modularization and differential evolutionary constraints on the modules involved in the signaling cascade.

Catechol, a major component of smoke, influences primary root growth and root hair elongation through reactive oxygen species-mediated redox signaling.

Nicotiana attenuata germinates from long-lived seedbanks in native soils after fires. Although smoke signals have been known to break seed dormancy, whether they also affect seedling establishment and root development remains unclear. In order to test this, seedlings were treated with smoke solutions. Seedlings responded in a dose-dependent manner with significantly increased primary root lengths, due mainly to longitudinal cell elongation, increased numbers of lateral roots and impaired root hair development. Bioassay-driven fractionations and NMR were used to identify catechol as the main active compound for the smoke-induced root phenotype. The transcriptome analysis revealed that mainly genes related to auxin biosynthesis and redox homeostasis were altered after catechol treatment. However, histochemical analyses of reactive oxygen species (ROS) and the inability of auxin applications to rescue the phenotype clearly indicated that highly localized changes in the root's redox-status, rather than in levels of auxin, are the primary effector. Moreover, H2 O2 application rescued the phenotype in a dose-dependent manner. Chemical cues in smoke not only initiate seed germination, but also influence seedling root growth; understanding how these cues work provides new insights into the molecular mechanisms by which plants adapt to post-fire environments.

Molecular mechanisms of adaptation and speciation: why do we need an integrative approach?

Understanding divergent adaptation and ecological speciation requires the synthesis of multiple approaches, including phenotypic characterization, genetics and genomics, realistic assessment of fitness and population genetic modelling. Current research in this field often approaches this problem from one of two directions: either a mechanistic approach-seeking to link phenotype, genotype and fitness, or a genomic approach-searching for signatures of divergence or selection across the genome. In most cases, these two approaches are not synthesized, and as a result, our understanding is incomplete. We argue that research in adaptation and evolutionary genetics needs to integrate these approaches for multiple reasons, including progress towards understanding the architecture and evolutionary history of adaptation and speciation loci, the ability to untangle linkage and pleiotropy, increased knowledge of mechanisms of genomic evolution and insights into parallel evolutionary events. Identifying the genetic underpinnings of adaptation and ecological speciation is not necessarily the end goal of research, but it is an integral part of understanding the evolutionary process. As a result, it is critical to utilize both genetic and genomic approaches. Challenges remain, particularly in nonmodel organisms and in our ability to synthesize results from multiple experimental systems. Nonetheless, advances in genetic and genomic techniques are increasingly available in a diverse array of systems, and the time is ripe to exploit the synthesis of these two approaches to increase our understanding of evolution.

Auxin Is Rapidly Induced by Herbivore Attack and Regulates a Subset of Systemic, Jasmonate-Dependent Defenses.

Plant responses to herbivore attack are regulated by phytohormonal networks. To date, the role of the auxin indole-3-acetic acid (IAA) in this context is not well understood. We quantified and manipulated the spatiotemporal patterns of IAA accumulation in herbivore-attacked Nicotiana attenuata plants to unravel its role in the regulation of plant secondary metabolism. We found that IAA is strongly, rapidly, and specifically induced by herbivore attack. IAA is elicited by herbivore oral secretions and fatty acid conjugate elicitors and is accompanied by a rapid transcriptional increase of auxin biosynthetic YUCCA-like genes. IAA accumulation starts 30 to 60 s after local induction and peaks within 5 min after induction, thereby preceding the jasmonate (JA) burst. IAA accumulation does not require JA signaling and spreads rapidly from the wound site to systemic tissues. Complementation and transport inhibition experiments reveal that IAA is required for the herbivore-specific, JA-dependent accumulation of anthocyanins and phenolamides in the stems. In contrast, IAA does not affect the accumulation of nicotine or 7-hydroxygeranyllinalool diterpene glycosides in the same tissue. Taken together, our results uncover IAA as a rapid and specific signal that regulates a subset of systemic, JA-dependent secondary metabolites in herbivore-attacked plants.

NaMYB8 regulates distinct, optimally distributed herbivore defense traits.

When herbivores attack, plants specifically reconfigure their metabolism. Herbivory on the wild tobacco Nicotiana attenuata strongly induces the R2R3 MYB transcriptional activator MYB8, which was reported to specifically regulate the accumulation of phenolamides (PAs). We discovered that transcriptional regulation of trypsin protease inhibitors (TPIs) and a threonine deaminase (TD) also depend on MYB8 expression. Induced distributions of PAs, TD and TPIs all meet predictions of optimal defense theory: their leaf concentrations increase with the fitness value and the probability of attack of the tissue. Therefore, we suggest that these defensive compounds have evolved to be co-regulated by MYB8.

Insect herbivory elicits genome-wide alternative splicing responses in Nicotiana attenuata.

Changes in gene expression and alternative splicing (AS) are involved in many responses to abiotic and biotic stresses in eukaryotic organisms. In response to attack and oviposition by insect herbivores, plants elicit rapid changes in gene expression which are essential for the activation of plant defenses; however, the herbivory-induced changes in AS remain unstudied. Using mRNA sequencing, we performed a genome-wide analysis on tobacco hornworm (Manduca sexta) feeding-induced AS in both leaves and roots of Nicotiana attenuata. Feeding by M. sexta for 5 h reduced total AS events by 7.3% in leaves but increased them in roots by 8.0% and significantly changed AS patterns in leaves and roots of existing AS genes. Feeding by M. sexta also resulted in increased (in roots) and decreased (in leaves) transcript levels of the serine/arginine-rich (SR) proteins that are involved in the AS machinery of plants and induced changes in SR gene expression that were jasmonic acid (JA)-independent in leaves but JA-dependent in roots. Changes in AS and gene expression elicited by M. sexta feeding were regulated independently in both tissues. This study provides genome-wide evidence that insect herbivory induces changes not only in the levels of gene expression but also in their splicing, which might contribute to defense against and/or tolerance of herbivory.

Herbivore associated elicitor-induced defences are highly specific among closely related Nicotiana species.

BACKGROUND: Herbivore-induced defence responses are often specific - different herbivores induce different defence responses in plants - and their specificity is largely mediated by chemical cues (herbivore-associated elicitors: HAEs) in insect oral or oviposition secretions. However, the specificity and the mechanisms of HAE-induced defence have not been investigated in the context of the evolutionary relationships among plant species. Here we compare the responses of six closely related Nicotiana species to a synthetic elicitor, N-linolenoyl-glutamic acid (C18:3-Glu) and HAE of two insect herbivores (the Solanaceae specialist Manduca sexta and generalist Spodoptera littoralis). RESULTS: HAE-induced defences are highly specific among closely related Nicotiana species at three perspectives. 1) A single Nicotiana species can elicit distinct responses to different HAEs. N. pauciflora elicited increased levels of JA and trypsin proteinase inhibitors (TPI) in response to C18:3-Glu and the oral secretions of M. sexta (OS Ms ) but not to oral secretions of S. littoralis (OS Sl ). In contrast, N. miersii only responded to OS Sl but not to the other two HAEs. The specific responses to different HAEs in Nicotiana species are likely due to the perception by the plant of each specific component of the HAE. 2) One HAE can induce different defence responses among closely related Nicotiana species. OS Ms and C18:3-Glu induced JA and TPI accumulations in N. linearis, N. attenuata, N. acuminata and N. pauciflora, but not in N. miersii and N. obtusifolia. 3) The effect of HAE-induced defences differ for the Solanaceae specialist M. sexta and the generalist S. littoralis. Among the four tested Nicotiana species, while the growth rate of M. sexta was only reduced by the induced defences elicited by C18:3-Glu; the growth rate of S. littoralis can be reduced by the induced defences elicited by all three HAEs. This is likely due to differences in the susceptibility of the specialist M. sexta and generalist S. littoralis to induced defences of their host. CONCLUSIONS: Closely related Nicotiana species elicit highly specific defence responses to herbivore associated elicitors and provide an ideal framework for investigating the molecular mechanisms and evolutionary divergence of induced resistance in plants.

The genetic basis of pollinator adaptation in a sexually deceptive orchid.

In plants, pollinator adaptation is considered to be a major driving force for floral diversification and speciation. However, the genetic basis of pollinator adaptation is poorly understood. The orchid genus Ophrys mimics its pollinators' mating signals and is pollinated by male insects during mating attempts. In many species of this genus, chemical mimicry of the pollinators' pheromones, especially of alkenes with different double-bond positions, plays a key role for specific pollinator attraction. Thus, different alkenes produced in different species are probably a consequence of pollinator adaptation. In this study, we identify genes that are likely involved in alkene biosynthesis, encoding stearoyl-acyl carrier protein (ACP) desaturases (SAD), in three closely related Ophrys species, O. garganica, O. sphegodes, and O. exaltata. Combining floral odor and gene expression analyses, two SAD homologs (SAD1/2) showed significant association with the production of (Z)-9- and (Z)-12-alkenes that were abundant in O. garganica and O. sphegodes, supporting previous biochemical data. In contrast, two other newly identified homologs (SAD5/6) were significantly associated with (Z)-7-alkenes that were highly abundant only in O. exaltata. Both molecular evolutionary analyses and pollinator preference tests suggest that the alkenes associated with SAD1/2 and SAD5/6 are under pollinator-mediated divergent selection among species. The expression patterns of these genes in F1 hybrids indicate that species-specific expression differences in SAD1/2 are likely due to cis-regulation, while changes in SAD5/6 are likely due to trans-regulation. Taken together, we report a genetic mechanism for pollinator-mediated divergent selection that drives adaptive changes in floral alkene biosynthesis involved in reproductive isolation among Ophrys species.

Stearoyl-acyl carrier protein desaturases are associated with floral isolation in sexually deceptive orchids.

The orchids Ophrys sphegodes and O. exaltata are reproductively isolated from each other by the attraction of two different, highly specific pollinator species. For pollinator attraction, flowers chemically mimic the pollinators' sex pheromones, the key components of which are alkenes with different double-bond positions. This study identifies genes likely involved in alkene biosynthesis, encoding stearoyl-acyl carrier protein (ACP) desaturase (SAD) homologs. The expression of two isoforms, SAD1 and SAD2, is flower-specific and broadly parallels alkene production during flower development. SAD2 shows a significant association with alkene production, and in vitro assays show that O. sphegodes SAD2 has activity both as an 18:0-ACP Δ(9) and a 16:0-ACP Δ(4) desaturase. Downstream metabolism of the SAD2 reaction products would give rise to alkenes with double-bonds at position 9 or position 12, matching double-bond positions observed in alkenes in the odor bouquet of O. sphegodes. SAD1 and SAD2 show evidence of purifying selection before, and positive or relaxed purifying selection after gene duplication. By contributing to the production of species-specific alkene bouquets, SAD2 is suggested to contribute to differential pollinator attraction and reproductive isolation among these species. Taken together, these data are consistent with the hypothesis that SAD2 is a florally expressed barrier gene of large phenotypic effect and, possibly, a genic target of pollinator-mediated selection.