Molecular mechanisms of mutualistic and antagonistic interactions in a plant-pollinator association.

Many insects metamorphose from antagonistic larvae into mutualistic adult pollinators, with reciprocal adaptation leading to specialized insect-plant associations. It remains unknown how such interactions are established at molecular level. Here we assemble high-quality genomes of a fig species, Ficus pumila var. pumila, and its specific pollinating wasp, Wiebesia pumilae. We combine multi-omics with validation experiments to reveal molecular mechanisms underlying this specialized interaction. In the plant, we identify the specific compound attracting pollinators and validate the function of several key genes regulating its biosynthesis. In the pollinator, we find a highly reduced number of odorant-binding protein genes and an odorant-binding protein mainly binding the attractant. During antagonistic interaction, we find similar chemical profiles and turnovers throughout the development of galled ovules and seeds, and a significant contraction of detoxification-related gene families in the pollinator. Our study identifies some key genes bridging coevolved mutualists, establishing expectations for more diffuse insect-pollinator systems.

Primula vulgaris (primrose) genome assembly, annotation and gene expression, with comparative genomics on the heterostyly supergene.

Primula vulgaris (primrose) exhibits heterostyly: plants produce self-incompatible pin- or thrum-form flowers, with anthers and stigma at reciprocal heights. Darwin concluded that this arrangement promotes insect-mediated cross-pollination; later studies revealed control by a cluster of genes, or supergene, known as the S (Style length) locus. The P. vulgaris S locus is absent from pin plants and hemizygous in thrum plants (thrum-specific); mutation of S locus genes produces self-fertile homostyle flowers with anthers and stigma at equal heights. Here, we present a 411 Mb P. vulgaris genome assembly of a homozygous inbred long homostyle, representing ~87% of the genome. We annotate over 24,000 P. vulgaris genes, and reveal more genes up-regulated in thrum than pin flowers. We show reduced genomic read coverage across the S locus in other Primula species, including P. veris, where we define the conserved structure and expression of the S locus genes in thrum. Further analysis reveals the S locus has elevated repeat content (64%) compared to the wider genome (37%). Our studies suggest conservation of S locus genetic architecture in Primula, and provide a platform for identification and evolutionary analysis of the S locus and downstream targets that regulate heterostyly in diverse heterostylous species.

Agrobacterium-mediated transformation systems of Primula vulgaris.

BACKGROUND: Genetic transformation is a valuable tool and an important procedure in plant functional genomics contributing to gene discovery, allowing powerful insights into gene function and genetically controlled characteristics. Primulaceae species provide one of the best-known examples of heteromorphic flower development, a breeding system which has attracted considerable attention, including that of Charles Darwin. Molecular approaches, including plant transformation give the best opportunity to define and understand the role of genes involved in floral heteromorphy in the common primrose, Primula vulgaris, along with other Primula species. RESULTS: Two transformation systems have been developed in P. vulgaris. The first system, Agrobacterium-mediated vacuum infiltration of seedlings, enables the rapid testing of transgenes, transiently in planta. GUS expression was observed in the cotyledons, true leaves, and roots of Primula seedlings. The second system is based on Agrobacterium tumefaciens infection of pedicel explants with an average transformation efficiency of 4.6%. This transformation system, based on regeneration and selection of transformants within in vitro culture, demonstrates stable transgene integration and transmission to the next generation. CONCLUSION: The two transformation systems reported here will aid fundamental research into important traits in Primula. Although, stable integration of transgenes is the ultimate goal for such analyses, transient gene expression via Agrobacterium-mediated DNA transfer, offers a simple and fast method to analyse transgene functions. The second system describes, for the first time, stable Agrobacterium-mediated transformation of Primula vulgaris, which will be key to characterising the genes responsible for the control of floral heteromorphy.

Genetic architecture and evolution of the S locus supergene in Primula vulgaris.

Darwin's studies on heterostyly in Primula described two floral morphs, pin and thrum, with reciprocal anther and stigma heights that promote insect-mediated cross-pollination. This key innovation evolved independently in several angiosperm families. Subsequent studies on heterostyly in Primula contributed to the foundation of modern genetic theory and the neo-Darwinian synthesis. The established genetic model for Primula heterostyly involves a diallelic S locus comprising several genes, with rare recombination events that result in self-fertile homostyle flowers with anthers and stigma at the same height. Here we reveal the S locus supergene as a tightly linked cluster of thrum-specific genes that are absent in pins. We show that thrums are hemizygous not heterozygous for the S locus, which suggests that homostyles do not arise by recombination between S locus haplotypes as previously proposed. Duplication of a floral homeotic gene 51.7 million years (Myr) ago, followed by its neofunctionalization, created the current S locus assemblage which led to floral heteromorphy in Primula. Our findings provide new insights into the structure, function and evolution of this archetypal supergene.

On the origins of observations of heterostyly in Primula.

In 1862, Charles Darwin published his landmark study on the different forms of flower in Primula; he coined the term distyly and subsequently expanded his studies to other species, including those with tristyly. Darwin is widely recognized as the first to study pin and thrum flowers in Primula, and to provide an explanation for the functional significance of the two floral morphs. Our laboratory is pursuing the genes that underpin floral heteromorphy in Primula, work influenced by Darwin's observations. One day, while appreciating a print of Primula vulgaris from William Curtis' Flora Londinensis, I was struck by the fact that I was looking at images of dimorphic Primula flowers captured in a late-1700s copper-plate engraving that predated Darwin's observations by over 70 yr. This realization triggered a journey into archives of botanical texts, herbals and florilegia from the 16(th) to 19(th) Centuries, and correspondence archives, in search of earlier documents that could have influenced Darwin and the origins of an idea. Darwin was not the first to observe floral heteromorphy in Primula, but he was the first to realize the significance of the two floral morphs. Darwin's insight and exposition of purpose have underpinned all consequent work on the subject.

Molecular characterisation of four double-flowered mutants of Silene dioica representing four centuries of variation.

Records of double-flowered Silene dioica date from the late sixteenth century and four named varieties are grown today, as previously, for their horticultural interest. Although double-flowered mutants have been characterized in several plants, their study in dioecious species is of particular interest due to influences of the homeotic mutation on the different floral whorl configurations in males and females. We have analysed four double-flowered varieties of Silene dioica: Flore Pleno and Rosea Plena date back to the seventeenth and nineteenth centuries, Thelma Kay and Firefly were recognized in the latter part of the twentieth and early twenty-first centuries. We have analysed the floral structure of the four varieties, which have distinct floral architectures. Based on Y chromosome-specific PCR analysis we show that Firefly is male and that the other three varieties are female: Random Amplification of Polymorphic DNA (RAPD) analyses suggested a common origin for the three female varieties. The double-flowered phenotype in all four varieties is caused by mutation of the C-function MADS-box transcription factor gene SDM1. We show that Firefly carries a unique 44bp insertion into SDM1, revealing an independent origin for this variety. Comparative analysis of SDM1 cDNA and genomic sequences in Flore Pleno, Rosea Plena and Thelma Kay shows that all three are caused by the same 7bp insertion within SDM1 and therefore share a common origin. The three alleles also differ by several single nucleotide polymorphisms, which represent somatic mutations accumulated over four centuries of asexual propagation.

Integration of genetic and physical maps of the Primula vulgaris S locus and localization by chromosome in situ hybridization.

Heteromorphic flower development in Primula is controlled by the S locus. The S locus genes, which control anther position, pistil length and pollen size in pin and thrum flowers, have not yet been characterized. We have integrated S-linked genes, marker sequences and mutant phenotypes to create a map of the P. vulgaris S locus region that will facilitate the identification of key S locus genes. We have generated, sequenced and annotated BAC sequences spanning the S locus, and identified its chromosomal location. We have employed a combination of classical genetics and three-point crosses with molecular genetic analysis of recombinants to generate the map. We have characterized this region by Illumina sequencing and bioinformatic analysis, together with chromosome in situ hybridization. We present an integrated genetic and physical map across the P. vulgaris S locus flanked by phenotypic and DNA sequence markers. BAC contigs encompass a 1.5-Mb genomic region with 1 Mb of sequence containing 82 S-linked genes anchored to overlapping BACs. The S locus is located close to the centromere of the largest metacentric chromosome pair. These data will facilitate the identification of the genes that orchestrate heterostyly in Primula and enable evolutionary analyses of the S locus.

Oakleaf: an S locus-linked mutation of Primula vulgaris that affects leaf and flower development.

In Primula vulgaris outcrossing is promoted through reciprocal herkogamy with insect-mediated cross-pollination between pin and thrum form flowers. Development of heteromorphic flowers is coordinated by genes at the S locus. To underpin construction of a genetic map facilitating isolation of these S locus genes, we have characterised Oakleaf, a novel S locus-linked mutant phenotype. We combine phenotypic observation of flower and leaf development, with classical genetic analysis and next-generation sequencing to address the molecular basis of Oakleaf. Oakleaf is a dominant mutation that affects both leaf and flower development; plants produce distinctive lobed leaves, with occasional ectopic meristems on the veins. This phenotype is reminiscent of overexpression of Class I KNOX-homeodomain transcription factors. We describe the structure and expression of all eight P. vulgaris PvKNOX genes in both wild-type and Oakleaf plants, and present comparative transcriptome analysis of leaves and flowers from Oakleaf and wild-type plants. Oakleaf provides a new phenotypic marker for genetic analysis of the Primula S locus. We show that none of the Class I PvKNOX genes are strongly upregulated in Oakleaf leaves and flowers, and identify cohorts of 507 upregulated and 314 downregulated genes in the Oakleaf mutant.

Floral heteromorphy in Primula vulgaris: progress towards isolation and characterization of the S locus.

BACKGROUND: The common primrose, Primula vulgaris, along with many other species of the Primulaceae, exhibits floral heteromorphy in which different individuals develop one of two possible forms of flower, known as pin and thrum. Both flower types are hermaphrodite and exhibit reciprocal positions of male and female reproductive structures, which together with a sporophytic incompatibility system, prevent self-pollination and promote out-crossing. The development of the two different forms of flower is controlled by a co-adapted linkage group of genes known as the S locus. SCOPE: Here progress towards identification and characterization of these genes is described to provide a molecular genetic explanation of the different floral characteristics that define heterostyly in Primula as observed and described by Charles Darwin. Previous work to identify and characterize developmental mutations linked to the P. vulgaris S locus, together with the isolation of S locus-linked genes and polymorphic DNA sequences markers, is summarized. The development of tools are described which will facilitate isolation and characterization of the S locus and its environs, including the creation of two expressed sequence tag libraries from pin and thrum flowers, as well as the construction and screening of two bacterial artificial chromosome (BAC) libraries containing thrum genomic DNA. Screening of these libraries with four S locus-linked sequences has enabled us to assemble four BAC contigs representing over 40 individual overlapping BAC clones which represent over 2·2 Mb of S locus-linked genomic sequence. PCR-based approaches for identification of the allelic origin of these BACs are described as well as identification of an additional 14 S locus-linked genes within BAC-end sequences. CONCLUSIONS: On-going work to assemble the four S locus-linked contigs into one contiguous sequence spanning the S locus is outlined in preparation for sequence analysis and characterization of the genes located within this region.

Hose in Hose, an S locus-linked mutant of Primula vulgaris, is caused by an unstable mutation at the Globosa locus.

Hose in Hose mutants of primrose and cowslip have been cultivated since the early 17th century and show dominant homeotic conversion of sepals to petals. The phenotype shows variable penetrance and expressivity and is linked to the S locus, which controls floral heteromorphy in Primula species. Here we demonstrate that the homeotic conversion of sepals to petals in Hose in Hose is associated with up-regulation of both Primula B-function MADS box genes PvDef and PvGlo in the first floral whorl. We have defined a restriction fragment length polymorphism associated with PvGlo that cosegregates with the Hose in Hose phenotype and have also identified and characterized a retrotransposon insertion in the PvGlo promoter which is associated with the up-regulated expression of PvGlo. Excision of this retrotransposon, associated with epigenetic changes at the locus, causes reversion toward normal calyces and restores wild-type flower development. These data define the molecular basis of the Hose in Hose mutation and provide an explanation for its long-documented phenotypic instability.

Wind-borne insects mediate directional pollen transfer between desert fig trees 160 kilometers apart.

The question of how far pollen can move between plants has implications for topics as diverse as habitat fragmentation, conservation management, and the containment of genetically modified crops. The monoecious African fig tree Ficus sycomorus L. relies on the small, short-lived, night-flying, host-specific fig wasp Ceratosolen arabicus Mayr for pollination. We used microsatellite markers to characterize a geographically isolated riparian population of F. sycomorus growing along the Ugab River in the Namib Desert, Namibia, together with paternity analysis of seedlings from known mothers, to map pollen movement within this population. In this way we tracked insect movements between individually recognizable trees by means of their pollen cargo and documented the movement of C. arabicus between known trees separated by more than 160 km, with a mean distance for confirmed successful pollination events of 88.6 km. The predominant observed movement of pollinators was in a westerly direction, toward the sea, reflecting seasonal nighttime wind direction and the wind-borne dispersal of fig wasps. Our results suggest the existence of an extensive panmictic population of trees that are well suited to overcome the effects of geographical isolation.

From gene expression to gene regulatory networks in Arabidopsis thaliana.

BACKGROUND: The elucidation of networks from a compendium of gene expression data is one of the goals of systems biology and can be a valuable source of new hypotheses for experimental researchers. For Arabidopsis, there exist several thousand microarrays which form a valuable resource from which to learn. RESULTS: A novel Bayesian network-based algorithm to infer gene regulatory networks from gene expression data is introduced and applied to learn parts of the transcriptomic network in Arabidopsis thaliana from a large number (thousands) of separate microarray experiments. Starting from an initial set of genes of interest, a network is grown by iterative addition to the model of the gene, from another defined set of genes, which gives the 'best' learned network structure. The gene set for iterative growth can be as large as the entire genome. A number of networks are inferred and analysed; these show (i) an agreement with the current literature on the circadian clock network, (ii) the ability to model other networks, and (iii) that the learned network hypotheses can suggest new roles for poorly characterized genes, through addition of relevant genes from an unconstrained list of over 15,000 possible genes. To demonstrate the latter point, the method is used to suggest that particular GATA transcription factors are regulators of photosynthetic genes. Additionally, the performance in recovering a known network from different amounts of synthetically generated data is evaluated. CONCLUSION: Our results show that plausible regulatory networks can be learned from such gene expression data alone. This work demonstrates that network hypotheses can be generated from existing gene expression data for use by experimental biologists.

The S locus-linked Primula homeotic mutant sepaloid shows characteristics of a B-function mutant but does not result from mutation in a B-function gene.

Floral homeotic and flower development mutants of Primula, including double, Hose in Hose, Jack in the Green and Split Perianth, have been cultivated since the late 1500s as ornamental plants but until recently have attracted limited scientific attention. Here we describe the characterization of a new mutant phenotype, sepaloid, that produces flowers comprising only sepals and carpels. The sepaloid mutation is recessive, and is linked to the S locus that controls floral heteromorphy. The phenotype shows developmental variability, with flowers containing three whorls of sepals surrounding fertile carpels, two whorls of sepals with a diminished third whorl of sepals surrounding a fourth whorl of carpels, or three whorls of sepals surrounding abnormal carpels. In some respects, these phenotypes resemble the Arabidopsis and Antirrhinum homeotic B-function mutants apetala3/deficiens (ap3/def) and pistillata/globosa (pi/glo). We have isolated the Primula vulgaris B-function genes PvDEFICIENS (PvDEF) and PvGLOBOSA (PvGLO), expression of both of which is affected in the sepaloid mutant. PvGLO, like sepaloid, is linked to the S locus, whereas PvDEF is not. However, our analyses reveal that sepaloid and PvGLO represent different genes. We conclude that SEPALOID is an S-linked independent regulator of floral organ identity genes including PvDEF and PvGLO.

Identification and characterization of pin and thrum alleles of two genes that co-segregate with the Primula S locus.

The study of heteromorphy in Primula over the past 140 years has established the reproductive significance of this breeding system. Plants produce either thrum or pin flowers that demonstrate reciprocal herkogamy. Thrums have short styles and produce large pollen from anthers at the mouth of the flower; pins have long styles and produce small pollen from anthers located within the corolla tube. The control of heteromorphy is orchestrated by the S locus with dominant (S) and recessive (s) alleles that comprise a co-adapted linkage group of genes. Thrum plants are heterozygous (Ss) and pin plants are homozygous (ss). Reciprocal crosses between the two forms are required for fertilization; within-morph crosses are impeded by a sporophytic self-incompatibility system. Rare recombination events within the S locus produce self-fertile homostyles. As a first step towards identifying genes located at the S locus, we used fluorescent differential display to screen for differential gene expression in pin and thrum flowers. Rather than only detecting differentially regulated genes, we identified two S locus linked genes by virtue of allelic variation between pin and thrum transcripts. Analysis of pin and thrum plants together with homostyle recombinant reveals that one gene flanks the locus, whereas the other shows complete linkage. One gene is related to Arabidopsis flower-timing genes Col9 and Col10; the other encodes a small predicted membrane protein of unknown function. Notwithstanding the diallelic behaviour of the Primula S locus, analysis of pin and thrum plants reveal three alleles for each gene: two pin and one thrum.

Conservation, convergence, and divergence of light-responsive, circadian-regulated, and tissue-specific expression patterns during evolution of the Arabidopsis GATA gene family.

In vitro analyses of plant GATA transcription factors have implicated some proteins in light-mediated and circadian-regulated gene expression, and, more recently, the analysis of mutants has uncovered further diverse roles for plant GATA factors. To facilitate function discovery for the 29 GATA genes in Arabidopsis (Arabidopsis thaliana), we have experimentally verified gene structures and determined expression patterns of all family members across adult tissues and suspension cell cultures, as well as in response to light and signals from the circadian clock. These analyses have identified two genes that are strongly developmentally light regulated, expressed predominantly in photosynthetic tissue, and with transcript abundance peaking before dawn. In contrast, several GATA factor genes are light down-regulated. The products of these light-regulated genes are candidates for those proteins previously implicated in light-regulated transcription. Coexpression of these genes with well-characterized light-responsive transcripts across a large microarray data set supports these predictions. Other genes show additional tissue-specific expression patterns suggesting novel and unpredicted roles. Genome-wide analysis using coexpression scatter plots for paralogous gene pairs reveals unexpected differences in cocorrelated gene expression profiles. Clustering the Arabidopsis GATA factor gene family by similarity of expression patterns reveals that genes of recent descent do not uniformly show conserved current expression profiles, yet some genes showing more distant evolutionary origins have acquired common expression patterns. In addition to defining developmental and environmental dynamics of GATA transcript abundance, these analyses offer new insights into the evolution of gene expression profiles following gene duplication events.

Analysis of late stage flower development in Primula vulgaris reveals novel differences in cell morphology and temporal aspects of floral heteromorphy.

Heterostyly in Primula is characterized by the development of long-styled pin and short-styled thrum flowers, with anthers midway down the corolla tube in pin flowers, and at its mouth in thrum flowers. Other differences include pollen size and stigmatic papillae length. Several linked genes at the S locus control these differences. In this study we have analyzed pin and thrum flowers through the temporal development of heteromorphy.These studies indicate that the S locus linked genes that orchestrate heteromorphic flower development act in coordination, but with different temporal and spatial dynamics. Style length is differentiated by longer style cells in pin than thrum. However, our studies on cell shape and size within the corolla tube show that a different mechanism mediates the dissimilar elevation of anthers between pin and thrum types. These studies have also revealed that upper corolla tube cells in thrum flowers are wider than those in pin flowers. This results in a larger corolla tube mouth in thrum flowers and represents a new and previously undocumented heteromorphic variation between pin and thrum flowers.

Arabidopsis Co-expression Tool (ACT): web server tools for microarray-based gene expression analysis.

The Arabidopsis Co-expression Tool, ACT, ranks the genes across a large microarray dataset according to how closely their expression follows the expression of a query gene. A database stores pre-calculated co-expression results for approximately 21,800 genes based on data from over 300 arrays. These results can be corroborated by calculation of co-expression results for user-defined sub-sets of arrays or experiments from the NASC/GARNet array dataset. Clique Finder (CF) identifies groups of genes which are consistently co-expressed with each other across a user-defined co-expression list. The parameters can be altered easily to adjust cluster size and the output examined for optimal inclusion of genes with known biological roles. Alternatively, a Scatter Plot tool displays the correlation coefficients for all genes against two user-selected queries on a scatter plot which can be useful for visual identification of clusters of genes with similar r-values. User-input groups of genes can be highlighted on the scatter plots. Inclusion of genes with known biology in sets of genes identified using CF and Scatter Plot tools allows inferences to be made about the roles of the other genes in the set and both tools can therefore be used to generate short lists of genes for further characterization. ACT is freely available at www.Arabidopsis.leeds.ac.uk/ACT.

The Arabidopsis co-expression tool (ACT): a WWW-based tool and database for microarray-based gene expression analysis.

We present a new WWW-based tool for plant gene analysis, the Arabidopsis Co-Expression Tool (ACT), based on a large Arabidopsis thaliana microarray data set obtained from the Nottingham Arabidopsis Stock Centre. The co-expression analysis tool allows users to identify genes whose expression patterns are correlated across selected experiments or the complete data set. Results are accompanied by estimates of the statistical significance of the correlation relationships, expressed as probability (P) and expectation (E) values. Additionally, highly ranked genes on a correlation list can be examined using the novel clique finder tool to determine the sets of genes most likely to be regulated in a similar manner. In combination, these tools offer three levels of analysis: creation of correlation lists of co-expressed genes, refinement of these lists using two-dimensional scatter plots, and dissection into cliques of co-regulated genes. We illustrate the applications of the software by analysing genes encoding functionally related proteins, as well as pathways involved in plant responses to environmental stimuli. These analyses demonstrate novel biological relationships underlying the observed gene co-expression patterns. To demonstrate the ability of the software to develop testable hypotheses on gene function within a defined biological process we have used the example of cell wall biosynthesis genes. The resource is freely available at http://www.arabidopsis.leeds.ac.uk/ACT/

Sex ratio strategies and the evolution of cue use.

Quantitative tests of sex allocation theory have often indicated that organism strategies deviate from model predictions. In pollinating fig wasps, Lipporrhopalum tentacularis, whole fig (brood) sex ratios are generally more female-biased than predicted by local mate competition (LMC) theory where females (foundresses) use density as a cue to assess potential LMC. We use microsatellite markers to investigate foundress sex ratios in L. tentacularis and show that they actually use their clutch size as a cue, with strategies closely approximating the predictions of a new model we develop of these conditions. We then provide evidence that the use of clutch size as a cue is common among species experiencing LMC, and given the other predictions of our model argue that this is because their ecologies mean it provides sufficiently accurate information about potential LMC that the use of other more costly cues has not evolved. We further argue that the use of these more costly cues by other species is due to the effect that ecological differences have on cue accuracy. This implies that deviations from earlier theoretical predictions often indicate that the cues used to assess environmental conditions differ from those assumed by models, rather than limits on the ability of natural selection to produce "perfect" organisms.

Analysis of members of the Silene latifolia Cys2/His2 zinc-finger transcription factor family during dioecious flower development and in a novel stamen-defective mutant ssf1.

Sex determination in dioecious Silene latifolia Poir. is governed by the inheritance of heteromorphic sex chromosomes. In male plants the Y chromosome influences two aspects of male organogenesis, the continued differentiation of stamen primordia and male fertility, and one aspect of female organogenesis, the arrest of development of the pistil. S. latifolia is susceptible to infection by the parasitic smut fungus Ustilago violacea, which induces stamen development in genetically female plants. Here we describe the identification and characterisation of a novel male mutant, short stamen filaments 1 (ssf1), defective in stamen differentiation. Although several independent studies have identified genes expressed during sex-determination in S. latifolia, analyses suggest that none of these encode regulatory proteins involved in the control of sex determination. We therefore isolated six S. latifolia cDNAs encoding members of a family of transcriptional regulators, the ZPT-type Cys2/His2 zinc-finger proteins that had previously shown to be co-ordinately regulated during stamen development in Petunia x hybrida hort. Vilm.-Andr. We have analysed the genomic organisation of these genes in male and female plants and their expression dynamics in male and female plants, in smut-infected female plants and in the ssf1 mutant. Our studies reveal expression patterns during development of the androecium that suggest a possible role for SlZPT2-1 in filament elongation and SlZPT4-1 in aspects of male fertility during stamen differentiation.