Isolation of the floral morph-related genes in heterostylous flax (Linum grandiflorum): the genetic polymorphism and the transcriptional and post-transcriptional regulations of the S locus

Genetic Causes and Genomic Consequences of Breakdown of Distyly in Linum trigynum

Distyly is an iconic floral polymorphism governed by a supergene, which promotes efficient pollen transfer and outcrossing through reciprocal differences in the position of sexual organs in flowers, often coupled with heteromorphic self-incompatibility. Distyly has evolved convergently in multiple flowering plant lineages, but has also broken down repeatedly, often resulting in homostylous, self-compatible populations with elevated rates of self-fertilization. Here, we aimed to study the genetic causes and genomic consequences of the shift to homostyly in Linum trigynum, which is closely related to distylous Linum tenue. Building on a high-quality genome assembly, we show that L. trigynum harbors a genomic region homologous to the dominant haplotype of the distyly supergene conferring long stamens and short styles in L. tenue, suggesting that loss of distyly first occurred in a short-styled individual. In contrast to homostylous Primula and Fagopyrum, L. trigynum harbors no fixed loss-of-function mutations in coding sequences of S-linked distyly candidate genes. Instead, floral gene expression analyses and controlled crosses suggest that mutations downregulating the S-linked LtWDR-44 candidate gene for male self-incompatibility and/or anther height could underlie homostyly and self-compatibility in L. trigynum. Population genomic analyses of 224 whole-genome sequences further demonstrate that L. trigynum is highly self-fertilizing, exhibits significantly lower genetic diversity genome-wide, and is experiencing relaxed purifying selection and less frequent positive selection on nonsynonymous mutations relative to L. tenue. Our analyses shed light on the loss of distyly in L. trigynum, and advance our understanding of a common evolutionary transition in flowering plants.

PfPIN5 promotes style elongation by regulating cell length in Primula forbesii Franch

BACKGROUND AND AIMS

Style dimorphism is one of the polymorphic characteristics of flowers in heterostylous plants, which have two types of flowers: the pin morph, with long styles and shorter anthers, and the thrum morph, with short styles and longer anthers. The formation of dimorphic styles has received attention in the plant world. Previous studies showed that CYP734A50 in Primula determined style length and limited style elongation and that the brassinosteroid metabolic pathway was involved in regulation of style length. However, it is unknown whether there are other factors affecting the style length of Primula.

METHODS

Differentially expressed genes highly expressed in pin morph styles were screened based on Primula forbesii transcriptome data. Virus-induced gene silencing was used to silence these genes, and the style length and anatomical changes were observed 20 days after injection.

KEY RESULTS

PfPIN5 was highly expressed in pin morph styles. When PfPIN5 was silenced, the style length was shortened in pin and long-homostyle plants by shortening the length of style cells. Moreover, silencing CYP734A50 in thrum morph plants increased the expression level of PfPIN5 significantly, and the style length increased. The results indicated that PfPIN5, an auxin efflux transporter gene, contributed to regulation of style elongation in P. forbesii.

CONCLUSIONS

The results implied that the auxin pathway might also be involved in the formation of styles of P. forbesii, providing a new pathway for elucidating the molecular mechanism of style elongation in P. forbesii.

The genetic control of herkogamy

Herkogamy is the spatial separation of anthers and stigmas within complete flowers, and is a key floral trait that promotes outcrossing in many angiosperms. The degree of separation between pollen-producing anthers and receptive stigmas has been shown to influence rates of self-pollination amongst plants, with a reduction in herkogamy increasing rates of successful selfing in self-compatible species. Self-pollination is becoming a critical issue in horticultural crops grown in environments where biotic pollinators are limited, absent, or difficult to utilise. In these cases, poor pollination results in reduced yield and misshapen fruit. Whilst there is a growing body of work elucidating the genetic basis of floral organ development, the genetic and environmental control points regulating herkogamy are poorly understood. A better understanding of the developmental and regulatory pathways involved in establishing varying degrees of herkogamy is needed to provide insights into the production of flowers more adept at selfing to produce consistent, high-quality fruit. This review presents our current understanding of herkogamy from a genetics and hormonal perspective.

Molecular insights into self-incompatibility systems: From evolution to breeding

Plants have evolved diverse self-incompatibility (SI) systems for outcrossing. Since Darwin's time, considerable progress has been made toward elucidating this unrivaled reproductive innovation. Recent advances in interdisciplinary studies and applications of biotechnology have given rise to major breakthroughs in understanding the molecular pathways that lead to SI, particularly the strikingly different SI mechanisms that operate in Solanaceae, Papaveraceae, Brassicaceae, and Primulaceae. These best-understood SI systems, together with discoveries in other "nonmodel" SI taxa such as Poaceae, suggest a complex evolutionary trajectory of SI, with multiple independent origins and frequent and irreversible losses. Extensive exploration of self-/nonself-discrimination signaling cascades has revealed a comprehensive catalog of male and female identity genes and modifier factors that control SI. These findings also enable the characterization, validation, and manipulation of SI-related factors for crop improvement, helping to address the challenges associated with development of inbred lines. Here, we review current knowledge about the evolution of SI systems, summarize key achievements in the molecular basis of pollen‒pistil interactions, discuss potential prospects for breeding of SI crops, and raise several unresolved questions that require further investigation.

Haplotype-resolved genome assembly provides insights into the evolution of S-locus supergene in distylous Nymphoides indica

Distyly has evolved independently in numerous animal-pollinated angiosperm lineages. Understanding of its molecular basis has been restricted to a few species, primarily Primula. Here, we investigate the genetic architecture of the single diallelic locus (S-locus) supergene, a linkage group of functionally associated genes, and explore how it may have evolved in distylous Nymphoides indica, a lineage of flowering plants not previously investigated. We assembled haplotype-resolved genomes, used read-coverage-based genome-wide association study (rb-GWAS) to locate the S-locus supergene, co-expression network analysis to explore gene networks underpinning the development of distyly, and comparative genomic analyses to investigate the origins of the S-locus supergene. We identified three linked candidate S-locus genes - NinBAS1, NinKHZ2, and NinS1 - that were only evident in the short-styled morph and were hemizygous. Co-expression network analysis suggested that brassinosteroids contribute to dimorphic sex organs in the short-styled morph. Comparative genomic analyses indicated that the S-locus supergene likely evolved via stepwise duplications and has been affected by transposable element activities. Our study provides novel insight into the structure, regulation, and evolution of the supergene governing distyly in N. indica. It also provides high-quality genomic resources for future research on the molecular mechanisms underlying the striking evolutionary convergence in form and function across heterostylous taxa.

Two floral forms in the same species-distyly

MAIN CONCLUSION

This paper reviews the progress of research on the morphology, physiology and molecular biology of distyly in plants. It will help to elucidate the mysteries of distyly in plants. Distyly is a unique representative type of heterostyly in plants, primarily characterized by the presence of long style and short style within the flowers of the same species. This interesting trait has always fascinated researchers. With the rapid development of molecular biology, the molecular mechanism for the production of dimorphic styles in plants is also gaining ground. Researchers have been studying plant dimorphic styles from various perspectives. The researchers are gradually unravelling the mechanisms by which plants produce distyly traits. This paper reviews advances in the study of plant dimorphic style characteristics, mainly in terms of the morphology, physiology and molecular biology of plants with dimorphic styles. The aim is to provide a theoretical basis for the study of the mechanism of distyly formation in plants.

Comprehensive transcriptomic profiling reveals complex molecular mechanisms in the regulation of style-length dimorphism in Guettarda speciosa (Rubiaceae), a species with "anomalous" distyly

BACKGROUND

The evolution of heterostyly, a genetically controlled floral polymorphism, has been a hotspot of research since the 19th century. In recent years, studies on the molecular mechanism of distyly (the most common form of heterostyly) revealed an evolutionary convergence in genes for brassinosteroids (BR) degradation in different angiosperm groups. This floral polymorphism often exhibits considerable variability that some taxa have significant stylar dimorphism, but anther height differs less. This phenomenon has been termed "anomalous" distyly, which is usually regarded as a transitional stage in evolution. Compared to "typical" distyly, the genetic regulation of "anomalous" distyly is almost unknown, leaving a big gap in our understanding of this special floral adaptation strategy.

METHODS

Here we performed the first molecular-level study focusing on this floral polymorphism in Guettarda speciosa (Rubiaceae), a tropical tree with "anomalous" distyly. Comprehensive transcriptomic profiling was conducted to examine which genes and metabolic pathways were involved in the genetic control of style dimorphism and if they exhibit similar convergence with "typical" distylous species.

RESULTS

"Brassinosteroid homeostasis" and "plant hormone signal transduction" was the most significantly enriched GO term and KEGG pathway in the comparisons between L- and S-morph styles, respectively. Interestingly, homologs of all the reported S-locus genes either showed very similar expressions between L- and S-morph styles or no hits were found in G. speciosa. BKI1, a negative regulator of brassinosteroid signaling directly repressing BRI1 signal transduction, was identified as a potential gene regulating style length, which significantly up-regulated in the styles of S-morph.

DISCUSSION

These findings supported the hypothesis that style length in G. speciosa was regulated through a BR-related signaling network in which BKI1 may be one key gene. Our data suggested, in species with "anomalous" distyly, style length was regulated by gene differential expressions, instead of the "hemizygous" S-locus genes in "typical" distylous flowers such as Primula and Gelsemium, representing an "intermediate" stage in the evolution of distyly. Genome-level analysis and functional studies in more species with "typical" and "anomalous" distyly would further decipher this "most complex marriage arrangement" in angiosperms and improve our knowledge of floral evolution.

Annotation of the Turnera subulata (Passifloraceae) Draft Genome Reveals the S-Locus Evolved after the Divergence of Turneroideae from Passifloroideae in a Stepwise Manner

A majority of Turnera species (Passifloraceae) exhibit distyly, a reproductive system involving both self-incompatibility and reciprocal herkogamy. This system differs from self-incompatibility in Passiflora species. The genetic basis of distyly in Turnera is a supergene, restricted to the S-morph, and containing three S-genes. How supergenes and distyly evolved in Turnera, and the other Angiosperm families exhibiting distyly remain largely unknown. Unraveling the evolutionary origins in Turnera requires the generation of genomic resources and extensive phylogenetic analyses. Here, we present the annotated draft genome of the S-morph of distylous Turnera subulata. Our annotation allowed for phylogenetic analyses of the three S-genes' families across 56 plant species ranging from non-seed plants to eudicots. In addition to the phylogenetic analysis, we identified the three S-genes' closest paralogs in two species of Passiflora. Our analyses suggest that the S-locus evolved after the divergence of Passiflora and Turnera. Finally, to provide insights into the neofunctionalization of the S-genes, we compared expression patterns of the S-genes with close paralogs in Arabidopsis and Populus trichocarpa. The annotation of the T. subulata genome will provide a useful resource for future comparative work. Additionally, this work has provided insights into the convergent nature of distyly and the origin of supergenes.

Genomic evidence supports the genetic convergence of a supergene controlling the distylous floral syndrome

Heterostyly, a plant sexual polymorphism controlled by the S-locus supergene, has evolved numerous times among angiosperm lineages and represents a classic example of convergent evolution in form and function. Determining whether underlying molecular convergence occurs could provide insights on constraints to floral evolution. Here, we investigated S-locus genes in distylous Gelsemium (Gelsemiaceae) to determine whether there is evidence of molecular convergence with unrelated distylous species. We used several approaches, including anatomical measurements of sex-organ development and transcriptome and whole-genome sequencing, to identify components of the S-locus supergene. We also performed evolutionary analysis with candidate S-locus genes and compared them with those reported in Primula and Turnera. The candidate S-locus supergene of Gelsemium contained four genes, of which three appear to have originated from gene duplication events within Gelsemiaceae. The style-length genes GeCYP in Gelsemium and CYP734A50 in Primula likely arose from duplication of the same gene, CYP734A1. Three out of four S-locus genes in Gelsemium elegans were hemizygous, as previously reported in Primula and Turnera. We provide genomic evidence on the genetic convergence of the supergene underlying distyly among distantly related angiosperm lineages and help to illuminate the genetic architecture involved in the evolution of heterostyly.

Genomic analyses of the Linum distyly supergene reveal convergent evolution at the molecular level

Supergenes govern multi-trait-balanced polymorphisms in a wide range of systems; however, our understanding of their origins and evolution remains incomplete. The reciprocal placement of stigmas and anthers in pin and thrum floral morphs of distylous species constitutes an iconic example of a balanced polymorphism governed by a supergene, the distyly S-locus. Recent studies have shown that the Primula and Turnera distyly supergenes are both hemizygous in thrums, but it remains unknown whether hemizygosity is pervasive among distyly S-loci. As hemizygosity has major consequences for supergene evolution and loss, clarifying whether this genetic architecture is shared among distylous species is critical. Here, we have characterized the genetic architecture and evolution of the distyly supergene in Linum by generating a chromosome-level genome assembly of Linum tenue, followed by the identification of the S-locus using population genomic data. We show that hemizygosity and thrum-specific expression of S-linked genes, including a pistil-expressed candidate gene for style length, are major features of the Linum S-locus. Structural variation is likely instrumental for recombination suppression, and although the non-recombining dominant haplotype has accumulated transposable elements, S-linked genes are not under relaxed purifying selection. Our findings reveal remarkable convergence in the genetic architecture and evolution of independently derived distyly supergenes, provide a counterexample to classic inversion-based supergenes, and shed new light on the origin and maintenance of an iconic floral polymorphism.

Variation in pollen sterility and gender specialization: an investigation with distylous species of Erythroxylum (Erythroxylaceae)

Differences in pollen sterility between morphs in distylous populations may represent the first step in the evolution of gender specialization. Theoretically, for partially or completely male sterile individuals to persist they must have some fitness advantage that compensates for the loss of male function. Gender specialization is considered a widespread process in Erythroxylum, but male sterility and the resource reallocation to female function have been investigated in few species and populations. In 18 populations of four distylous species of Erythroxylum, we quantified the levels of male sterility, estimated through pollen sterility, in short- and long-styled flowers to test if sterility is morph-biased. In one population per species, we also described the frequency of floral visitors, the production of flowers and fruits, and the quality of fruits and seeds of short- and long-styled plants to evaluate the expression of trade-offs in allocation to male and female function. In some populations of E. campestre and E. deciduum, short-styled flowers possessed higher levels of pollen sterility than long-styled flowers. Although most flowers of E. suberosum and E. tortuosum also expressed pollen sterility, the frequency of sterility was similar between morphs in all populations. Differences in reproductive output between morphs occurred only in populations of species with morph-biased sterility, but none of the variations reflect resource allocation to female fitness of short-styled plants. Differences in the level of sterility between morphs indicates the potential for gender specialization in populations of E. campestre and E. deciduum, despite the apparent lack of a trade-off in allocation.

A Draft Genome of the Ginger Species Alpinia nigra and New Insights into the Genetic Basis of Flexistyly

Angiosperms possess various strategies to ensure reproductive success, such as stylar polymorphisms that encourage outcrossing. Here, we investigate the genetic basis of one such dimorphism that combines both temporal and spatial separation of sexual function, termed flexistyly. It is a floral strategy characterised by the presence of two morphs that differ in the timing of stylar movement. We performed a de novo assembly of the genome of Alpinia nigra using high-depth genomic sequencing. We then used Pool-seq to identify candidate regions for flexistyly based on allele frequency or coverage differences between pools of anaflexistylous and cataflexistylous morphs. The final genome assembly size was 2 Gb, and showed no evidence of recent polyploidy. The Pool-seq did not reveal large regions with high FST values, suggesting large structural chromosomal polymorphisms are unlikely to underlie differences between morphs. Similarly, no region had a 1:2 mapping depth ratio which would be indicative of hemizygosity. We propose that flexistyly is governed by a small genomic region that might be difficult to detect with Pool-seq, or a complex genomic region that proved difficult to assemble. Our genome will be a valuable resource for future studies of gingers, and provides the first steps towards characterising this complex floral phenotype.

The Genomic Architecture and Evolutionary Fates of Supergenes

Supergenes are genomic regions containing sets of tightly linked loci that control multi-trait phenotypic polymorphisms under balancing selection. Recent advances in genomics have uncovered significant variation in both the genomic architecture as well as the mode of origin of supergenes across diverse organismal systems. Although the role of genomic architecture for the origin of supergenes has been much discussed, differences in the genomic architecture also subsequently affect the evolutionary trajectory of supergenes and the rate of degeneration of supergene haplotypes. In this review, we synthesize recent genomic work and historical models of supergene evolution, highlighting how the genomic architecture of supergenes affects their evolutionary fate. We discuss how recent findings on classic supergenes involved in governing ant colony social form, mimicry in butterflies, and heterostyly in flowering plants relate to theoretical expectations. Furthermore, we use forward simulations to demonstrate that differences in genomic architecture affect the degeneration of supergenes. Finally, we discuss implications of the evolution of supergene haplotypes for the long-term fate of balanced polymorphisms governed by supergenes.

Parallelism in Flower Evolution and Development

Flower evolution is characterized by widespread repetition, with adaptations to pollinator environment evolving in parallel. Recent studies have expanded our understanding of the developmental basis of adaptive floral novelties—petal fusion, bilateral symmetry, heterostyly, and floral dimensions. In this article, we describe patterns of trait evolution and review developmental genetic mechanisms underlying floral novelties. We discuss the diversity of mechanisms for parallel adaptation, the evidence for constraints on these mechanisms, and how constraints help explain observed macroevolutionary patterns. We describe parallel evolution resulting from similarities at multiple hierarchical levels—genetic, developmental, morphological, functional—which indicate general principles in floral evolution, including the central role of hormone signaling. An emerging pattern is mutational bias that may contribute to rapid patterns of parallel evolution, especially if the derived trait can result from simple degenerative mutations. We argue that such mutational bias may be less likely to govern the evolution of novelties patterned by complex developmental pathways.

Transcriptome and Network Analyses of Heterostyly in Turnera subulata Provide Mechanistic Insights: Are S-Loci a Red-Light for Pistil Elongation?

Heterostyly employs distinct hermaphroditic floral morphs to enforce outbreeding. Morphs differ structurally in stigma/anther positioning, promoting cross-pollination, and physiologically blocking self-fertilization. Heterostyly is controlled by a self-incompatibility (S)-locus of a small number of linked S-genes specific to short-styled morph genomes. Turnera possesses three S-genes, namely TsBAHD (controlling pistil characters), TsYUC6, and TsSPH1 (controlling stamen characters). Here, we compare pistil and stamen transcriptomes of floral morphs of T. subulata to investigate hypothesized S-gene function(s) and whether hormonal differences might contribute to physiological incompatibility. We then use network analyses to identify genetic networks underpinning heterostyly. We found a depletion of brassinosteroid-regulated genes in short styled (S)-morph pistils, consistent with hypothesized brassinosteroid-inactivating activity of TsBAHD. In S-morph anthers, auxin-regulated genes were enriched, consistent with hypothesized auxin biosynthesis activity of TsYUC6. Evidence was found for auxin elevation and brassinosteroid reduction in both pistils and stamens of S- relative to long styled (L)-morph flowers, consistent with reciprocal hormonal differences contributing to physiological incompatibility. Additional hormone pathways were also affected, however, suggesting S-gene activities intersect with a signaling hub. Interestingly, distinct S-genes controlling pistil length, from three species with independently evolved heterostyly, potentially intersect with phytochrome interacting factor (PIF) network hubs which mediate red/far-red light signaling. We propose that modification of the activities of PIF hubs by the S-locus could be a common theme in the evolution of heterostyly.

Buckwheat heteromorphic self-incompatibility: genetics, genomics and application to breeding

Common buckwheat (Fagopyrum esculentum Moench 2n = 2x = 16) is an outcrossing crop with heteromorphic self-incompatibility due to its distylous flowers, called pin and thrum. In pin plants, a long style is combined with short stamens and small pollen grains; in thrum plants, a short style is combined with long stamens and large pollen grains. Both the intra-morph self-incompatibility and flower morphology are controlled by a single genetic locus named the S locus; thrum plants are heterozygous (Ss) and pin plants are homozygous recessive (ss) at this locus. Self-incompatibility is an obstacle for establishing pure lines and fixation of agronomically useful genes. Elucidation of the molecular mechanism of heterostylous self-incompatibility of common buckwheat has continued for a quarter of a century. Recent advances in genomic and transcriptomic analyses using next-generation sequencing have made it possible to determine the genomic region harboring the buckwheat S locus and to identify novel genes at this locus. In this review, we summarize the current knowledge on buckwheat heterostyly gained from conventional and molecular genetics and genomics. We also discuss the application of these studies to breeding of common buckwheat.

The long and short of the S-locus in Turnera (Passifloraceae)

Distyly is an intriguing floral adaptation that increases pollen transfer precision and restricts inbreeding. It has been a model system in evolutionary biology since Darwin. Although the S-locus determines the long- and short-styled morphs, the genes were unknown in Turnera. We have now identified these genes. We used deletion mapping to identify, and then sequence, BAC clones and genome scaffolds to construct S/s haplotypes. We investigated candidate gene expression, hemizygosity, and used mutants, to explore gene function. The s-haplotype possessed 21 genes collinear with a region of chromosome 7 of grape. The S-haplotype possessed three additional genes and two inversions. TsSPH1 was expressed in filaments and anthers, TsYUC6 in anthers and TsBAHD in pistils. Long-homostyle mutants did not possess TsBAHD and a short-homostyle mutant did not express TsSPH1. Three hemizygous genes appear to determine S-morph characteristics in T. subulata. Hemizygosity is common to all distylous species investigated, yet the genes differ. The pistil candidate gene, TsBAHD, differs from that of Primula, but both may inactivate brassinosteroids causing short styles. TsYUC6 is involved in auxin synthesis and likely determines pollen characteristics. TsSPH1 is likely involved in filament elongation. We propose an incompatibility mechanism involving TsYUC6 and TsBAHD.

'A most complex marriage arrangement': recent advances on heterostyly and unresolved questions

Heterostylous genetic polymorphisms provide paradigmatic systems for investigating adaptation and natural selection. Populations are usually comprised of two (distyly) or three (tristyly) mating types, maintained by negative frequency-dependent selection resulting from disassortative mating. Theory predicts this mating system should result in equal style-morph ratios (isoplethy) at equilibrium. Here, I review recent advances on heterostyly, focusing on examples challenging stereotypical depictions of the polymorphism and unresolved questions. Comparative analyses indicate multiple origins of heterostyly, often within lineages. Ecological studies demonstrate that structural components of heterostyly are adaptations improving the proficiency of animal-mediated cross-pollination and reducing pollen wastage. Both neutral and selective processes cause deviations from isoplethy in heterostylous populations, and, under some ecological and demographic conditions, cause breakdown of the polymorphism, resulting in either the evolution of autogamy and mixed mating, or transitions to alternative outcrossing systems, including dioecy. Earlier ideas on the genetic architecture of the S-locus supergene governing distyly have recently been overturned by discovery that the dominant S-haplotype is a hemizygous region absent from the s-haplotype. Ecological, phylogenetic and molecular genetic data have validated some features of theoretical models on the selection of the polymorphism. Although heterostyly is the best-understood floral polymorphism in angiosperms, many unanswered questions remain.

Identification of a gene encoding polygalacturonase expressed specifically in short styles in distylous common buckwheat (Fagopyrum esculentum)

Common buckwheat (Fagopyrum esculentum) is a heteromorphic self-incompatible (SI) species with two types of floral architecture: thrum (short style) and pin (long style). The floral morphology and intra-morph incompatibility are controlled by a single genetic locus, S. However, the molecular mechanisms underlying the heteromorphic self-incompatibility of common buckwheat remain unclear. To identify these mechanisms, we performed proteomic, quantitative reverse-transcription PCR, and linkage analyses. Comparison of protein profiles between the long and short styles revealed a protein unique to the short style. Amino-acid sequencing revealed that it was a truncated form of polygalacturonase (PG); we designated the gene encoding this protein FePG1. Phylogenetic analysis classified FePG1 into the same clade as PGs that function in pollen development and floral morphology. FePG1 expression was significantly higher in short styles than in long styles. It was expressed in flowers of a short-homostyle line but not in flowers of a long-homostyle line. Linkage analysis indicated that FePG1 was not linked to the S locus; it could be a factor downstream of this locus. Our finding of a gene putatively working under the regulation of the S locus provides useful information for elucidation of the mechanism of heteromorphic self-incompatibility.

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.

GmMYB181, a Soybean R2R3-MYB Protein, Increases Branch Number in Transgenic Arabidopsis

Soybean (Glycine max) is an important economic crop that provides abundant oil and high quality protein for human beings. As the process of reproductive growth directly determines the crop seed yield and quality, we initiated studies to identify genes that regulate soybean floral organ development. One R2R3-MYB transcription factor gene, designated as GmMYB181, was found to be enriched in flowers based on microarray analysis and was further functionally investigated in transgenic Arabidopsis. GmMYB181 protein contains two MYB domains, which localized to the nucleus and displayed transcriptional activation in yeast hybrid system. Real-time quantitative PCR (qRT-PCR) results suggested GmMYB181 exclusively expressed in flower tissue. In Arabidopsis, overexpression of GmMYB181 altered the morphology of floral organs, fruit size and plant architecture, including outward curly sepals, smaller siliques, increased lateral branches and reduced plant height, indicating that GmMYB181 is involved in the development of reproductive organs and plays an important role in controlling plant architecture. Further, microarray analysis revealed that overexpressing GmMYB181 in Arabidopsis affected the expression of 3450 genes in mature flowers, including those involved in floral organ, seed/fruit development, and responded to different hormone signals.

Style polymorphism in Linum (Linaceae): a case of Mediterranean parallel evolution?

Heterostyly is a sex polymorphism that has challenged evolutionary biologists ever since Darwin. One of the lineages where heterostyly and related stylar conditions appear more frequently is Linum (Linaceae). This group is particularly suitable for testing competing hypotheses about ancestral and transitional stages on the evolutionary building up of heterostyly. We generated a phylogeny of Linum based on extensive sampling and plastid and nuclear DNA sequences, and used it to trace the evolution of character states of style polymorphism. We also revised available data on pollination, breeding systems, and polyploidy to analyse their associations. Our results supported former phylogenetic hypotheses: the paraphyly of Linum and the non-monophyly of current taxonomic sections. Heterostyly was common in the genus, but appeared concentrated in the Mediterranean Basin and the South African Cape. Ancestral character state reconstruction failed to determine a unique state as the most probable condition for style polymorphism in the genus. In contrast, approach herkogamy was resolved as ancestral state in some clades, giving support to recent hypotheses. Some traits putatively related to heterostyly, such as life history and polyploidy, did show marginal or non-significant phylogenetic correlation, respectively. Although pollinator data are limited, we suggest that beeflies are associated with specific cases of heterostyly. The consistent association between style polymorphism and heteromorphic incompatibility points to ecological factors as drivers of the multiple evolution of style polymorphism in Linum. Albeit based on limited evidence, we hypothesised that specialised pollinators and lack of mating opportunities drive evolution of style polymorphism and loss of the polymorphism, respectively.

A short story gets longer: recent insights into the molecular basis of heterostyly

Heterostyly is a fascinating adaptation to promote outbreeding and a classical paradigm of botany. In the most common type of heterostyly, plants either form flowers with long styles and short stamens, or short styles and long stamens. This reciprocal organ positioning reduces pollen wastage and promotes cross-pollination, thus increasing male fitness. In addition, in many heterostylous species selfing and the generation of unfit progeny due to inbreeding depression is limited by a self-incompatibility system, thus promoting female fitness. The two floral forms are genetically determined by the S locus as a complex supergene, namely a chromosomal region containing several individual genes that control the different traits, such as style or stamen length, and are held together by very tight linkage due to suppressed recombination. Recent molecular-genetic studies in several systems, including Turnera, Fagopyrum, Linum, and Primula have begun to identify and characterize the causal heterostyly genes residing at the S locus. An emerging theme from several families is that the dominant S haplotype represents a hemizygous region not present on the recessive s haplotype. This provides an explanation for the suppressed recombination and suggests a scenario for the chromosomal evolution of the S locus. In this review, we discuss the results from recent molecular-genetic analyses in light of the classical models on the genetics and evolution of heterostyly.

Comparative proteomic analysis of eggplant (Solanum melongena L.) heterostylous pistil development

Heterostyly is a common floral polymorphism, but the proteomic basis of this trait is still largely unexplored. In this study, self- and cross-pollination of L-morph and S-morph flowers and comparison of embryo sac development in eggplant (Solanum melongena L.) suggested that lower fruit set from S-morph flowers results from stigma-pollen incompatibility. To explore the molecular mechanism underlying heterostyly development, we conducted isobaric tags for relative and absolute quantification (iTRAQ) proteomic analysis of eggplant pistils for L- and S-morph flowers. A total of 5,259 distinct proteins were identified during heterostyly development. Compared S-morph flowers with L-morph, we discovered 57 and 184 differentially expressed proteins (DEPs) during flower development and maturity, respectively. Quantitative real time polymerase chain reactions were used for nine genes to verify DEPs from the iTRAQ approach. During flower development, DEPs were mainly involved in morphogenesis, biosynthetic processes, and metabolic pathways. At flower maturity, DEPs primarily participated in biosynthetic processes, metabolic pathways, and the formation of ribosomes and proteasomes. Additionally, some proteins associated with senescence and programmed cell death were found to be upregulated in S-morph pistils, which may lead to the lower fruit set in S-morph flowers. Although the exact roles of these related proteins are not yet known, this was the first attempt to use an iTRAQ approach to analyze proteomes of heterostylous eggplant flowers, and these results will provide insights into biochemical events taking place during the development of heterostyly.

Arabidopsis aspartic proteases A36 and A39 play roles in plant reproduction

Aspartic proteases (Aps, EC3.4.23) are one of the 4 major mechanistic classes of proteolytic enzymes with the conserved motifs Asp-Thr/Ser-Gly (DT/SG) at the active site and are activated at acidic pH. In Arabidopsis, 69 genes were identified as coding putative aspartic proteinases. However, little is known about most of these enzymes. Recently, we characterized 2 novel Arabidopsis Aps genes, A36 and A39, which encode 2 putative GPI-anchored pollen-high-expressed Aps. a36 a39 mutants display significant abortion. The pollen grains underwent apoptosis-like programmed cell death and the degeneration of female gametes was also appeared in the a36 a39 mutant. Besides, the pollen tube of a36 a39 has compromised micropylar guidance. A36 and A39 were membrane-anchored protein and co-localized with a reported GPI-anchored protein COBRA-LIKE 10 (COBL10). In apical region of a36 a39 pollen tubes cell wall, the abundance of highly methlyestered homogalacturonans and xyloglucans were significantly increased. These results indicated that A36 and A39 are vital factors involved in gametogenesis and pollen guidance in Arabidopsis.

De novo Sequencing and Comparative Transcriptomics of Floral Development of the Distylous Species Lithospermum multiflorum

Genes controlling the morphological, micromorphological, and physiological components of the breeding system distyly have been hypothesized, but many of the genes have not been investigated throughout development of the two floral morphs. To this end, the present study is an examination of comparative transcriptomes from three stages of development for the floral organs of the morphs of Lithospermum multiflorum. Transcriptomes of flowers of the two morphs, from various stages of development, were sequenced using an Illumina HiSeq 2000. The floral transcriptome of L. multiflorum was assembled, and differential gene expression (DE) was identified between morphs, throughout development. Additionally, Gene Ontology (GO) terms for DE genes were determined. Fewer genes were DE early in development compared to later in development, with more genes highly expressed in the gynoecium of the SS morph and the corolla and androecium of the LS morph. A reciprocal pattern was observed later in development, and many more genes were DE during this latter stage. During early development, DE genes appear to be involved in growth and floral development, and during later development, DE genes seem to affect physiological functions. Interestingly, many genes involved in response to stress were identified as DE between morphs.

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.

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.

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.