Silene as a model system in ecology and evolution

Pas de deux: An Intricate Dance of Anther Smut and Its Host

The successful interaction between pathogen/parasite and host requires a delicate balance between fitness of the former and survival of the latter. To optimize fitness a parasite/pathogen must effectively create an environment conducive to reproductive success, while simultaneously avoiding or minimizing detrimental host defense response. The association between Microbotryum lychnidis-dioicae and its host Silene latifolia serves as an excellent model to examine such interactions. This fungus is part of a species complex that infects species of the Caryophyllaceae, replacing pollen with the fungal spores. In the current study, transcriptome analyses of the fungus and its host were conducted during discrete stages of bud development so as to identify changes in fungal gene expression that lead to spore development and to identify changes associated with infection in the host plant. In contrast to early biotrophic phase stages of infection for the fungus, the latter stages involve tissue necrosis and in the case of infected female flowers, further changes in the developmental program in which the ovary aborts and a pseudoanther is produced. Transcriptome analysis via Illumina RNA sequencing revealed enrichment of fungal genes encoding small secreted proteins, with hallmarks of effectors and genes found to be relatively unique to the Microbotryum species complex. Host gene expression analyses also identified interesting sets of genes up-regulated, including those involving stress response, host defense response, and several agamous-like MADS-box genes (AGL61 and AGL80), predicted to interact and be involved in male gametophyte development.

Paternal arbuscular mycorrhizal fungal status affects DNA methylation in seeds

Most land plants grow in association with arbuscular mycorrhizal fungi (AMF) in their roots and these fungi can cause transgenerational effects on plants' offspring. These may be caused by changes in DNA methylation of the offspring. In this study, we compared the amount of global DNA methylation in seeds of the gynodioecious plant Geranium sylvaticum in relation to the gender and the AMF status of the parents producing the seeds. The amount of DNA methylated was positively related to seed mass. Seeds produced by females had a similar proportion of methylated DNA regardless of the AMF status of the father siring the seed. By contrast, seeds from hermaphrodites had higher DNA methylation when sired by AMF fathers. We show to the best of our knowledge for the first time, that the AMF status of fathers can affect DNA methylation in seeds and that these changes in DNA methylation are further dependent on the gender of the mother producing the seeds.

Transmission dynamics: critical questions and challenges

This article overviews the dynamics of disease transmission in one-host-one-parasite systems. Transmission is the result of interacting host and pathogen processes, encapsulated with the environment in a 'transmission triangle'. Multiple transmission modes and their epidemiological consequences are often not understood because the direct measurement of transmission is difficult. However, its different components can be analysed using nonlinear transmission functions, contact matrices and networks. A particular challenge is to develop such functions for spatially extended systems. This is illustrated for vector transmission where a 'perception kernel' approach is developed that incorporates vector behaviour in response to host spacing. A major challenge is understanding the relative merits of the large number of approaches to quantifying transmission. The evolution of transmission mode itself has been a rather neglected topic, but is important in the context of understanding disease emergence and genetic variation in pathogens. Disease impacts many biological processes such as community stability, the evolution of sex and speciation, yet the importance of different transmission modes in these processes is not understood. Broader approaches and ideas to disease transmission are important in the public health realm for combating newly emerging infections.This article is part of the themed issue 'Opening the black box: re-examining the ecology and evolution of parasite transmission'.

The Evolution of Sex Chromosomes and Dosage Compensation in Plants

Plant sex chromosomes can be vastly different from those of the few historical animal model organisms from which most of our understanding of sex chromosome evolution is derived. Recently, we have seen several advancements from studies on green algae, brown algae, and land plants that are providing a broader understanding of the variable ways in which sex chromosomes can evolve in distant eukaryotic groups. Plant sex-determining genes are being identified and, as expected, are completely different from those in animals. Species with varying levels of differentiation between the X and Y have been found in plants, and these are hypothesized to be representing different stages of sex chromosome evolution. However, we are also finding that sex chromosomes can remain morphologically unchanged over extended periods of time. Where degeneration of the Y occurs, it appears to proceed similarly in plants and animals. Dosage compensation (a phenomenon that compensates for the consequent loss of expression from the Y) has now been documented in a plant system, its mechanism, however, remains unknown. Research has also begun on the role of sex chromosomes in sexual conflict resolution, and it appears that sex-biased genes evolve similarly in plants and animals, although the functions of these genes remain poorly studied. Because the difficulty in obtaining sex chromosome sequences is increasingly being overcome by methodological developments, there is great potential for further discovery within the field of plant sex chromosome evolution.

Evaluation of reference genes for reverse transcription quantitative real-time PCR (RT-qPCR) studies in Silene vulgaris considering the method of cDNA preparation

Accurate gene expression measurements are essential in studies of both crop and wild plants. Reverse transcription quantitative real-time PCR (RT-qPCR) has become a preferred tool for gene expression estimation. A selection of suitable reference genes for the normalization of transcript levels is an essential prerequisite of accurate RT-qPCR results. We evaluated the expression stability of eight candidate reference genes across roots, leaves, flower buds and pollen of Silene vulgaris (bladder campion), a model plant for the study of gynodioecy. As random priming of cDNA is recommended for the study of organellar transcripts and poly(A) selection is indicated for nuclear transcripts, we estimated gene expression with both random-primed and oligo(dT)-primed cDNA. Accordingly, we determined reference genes that perform well with oligo(dT)- and random-primed cDNA, making it possible to estimate levels of nucleus-derived transcripts in the same cDNA samples as used for organellar transcripts, a key benefit in studies of cyto-nuclear interactions. Gene expression variance was estimated by RefFinder, which integrates four different analytical tools. The SvACT and SvGAPDH genes were the most stable candidates across various organs of S. vulgaris, regardless of whether pollen was included or not.

Fungal Sex: The Basidiomycota

Fungi of the Basidiomycota, representing major pathogen lineages and mushroom-forming species, exhibit diverse means to achieve sexual reproduction, with particularly varied mechanisms to determine compatibilities of haploid mating partners. For species that require mating between distinct genotypes, discrimination is usually based on both the reciprocal exchange of diffusible mating pheromones, rather than sexes, and the interactions of homeodomain protein signals after cell fusion. Both compatibility factors must be heterozygous in the product of mating, and genetic linkage relationships of the mating pheromone/receptor and homeodomain genes largely determine the complex patterns of mating-type variation. Independent segregation of the two compatibility factors can create four haploid mating genotypes from meiosis, referred to as tetrapolarity. This condition is thought to be ancestral to the basidiomycetes. Alternatively, cosegregation by linkage of the two mating factors, or in some cases the absence of the pheromone-based discrimination, yields only two mating types from meiosis, referred to as bipolarity. Several species are now known to have large and highly rearranged chromosomal regions linked to mating-type genes. At the population level, polymorphism of the mating-type genes is an exceptional aspect of some basidiomycete fungi, where selection under outcrossing for rare, intercompatible allelic variants is thought to be responsible for numbers of mating types that may reach several thousand. Advances in genome sequencing and assembly are yielding new insights by comparative approaches among and within basidiomycete species, with the promise to resolve the evolutionary origins and dynamics of mating compatibility genetics in this major eukaryotic lineage.

Divergence in style length and pollen size leads to a postmating-prezygotic reproductive barrier among populations of Silene latifolia

A central tenet of speciation research is the need to identify reproductive isolating barriers. One approach to this line of research is to identify the phenotypes that lead to reproductive isolation. Several studies on flowering plants have shown that differences in style length contribute to reproductive isolation between species, leading us to consider whether style length could act as a reproductive barrier among populations of a single species. This could occur if style length varied sufficiently and pollen size covaried with style length. Populations of Silene latifolia exhibit variation in flower size, including style length, that is negatively correlated with annual precipitation. We show that this divergence in style length has a genetic basis and acts as a reproductive barrier: males from small-flowered populations produced relatively small pollen grains that were poor at fertilizing ovules when crossed to females from large-flowered populations, leading to a significant reduction in seed production. Manipulating the distance pollen tubes had to travel revealed that this failure was purely mechanical and not the result of other incompatibilities. These results show that style length acts as a postmating-prezygotic reproductive barrier and indicate a potential link between ecotypic differentiation and reproductive isolation within a species.

Co-occurrence and hybridization of anther-smut pathogens specialized on Dianthus hosts

Host specialization has important consequences for the diversification and ecological interactions of obligate pathogens. The anther-smut disease of natural plant populations, caused by Microbotryum fungi, has been characterized by specialized host-pathogen interactions, which contribute in part to the isolation among these numerous fungal species. This study investigated the molecular variation of Microbotryum pathogens within the geographic and host-specific distributions on wild Dianthus species in southern European Alps. In contrast to prior studies on this pathogen genus, a range of overlapping host specificities was observed for four delineated Microbotryum lineages on Dianthus hosts, and their frequent co-occurrence within single-host populations was quantified at local and regional scales. In addition to potential consequences for direct pathogen competition, the sympatry of Microbotryum lineages led to hybridization between them in many populations, and these admixed genotypes suffered significant meiotic sterility. Therefore, this investigation of the anther-smut fungi reveals how variation in the degrees of host specificity can have major implications for ecological interactions and genetic integrity of differentiated pathogen lineages.

Non-coding RNA may be associated with cytoplasmic male sterility in Silene vulgaris

Cytoplasmic male sterility (CMS) is a widespread phenomenon in flowering plants caused by mitochondrial (mt) genes. CMS genes typically encode novel proteins that interfere with mt functions and can be silenced by nuclear fertility-restorer genes. Although the molecular basis of CMS is well established in a number of crop systems, our understanding of it in natural populations is far more limited. To identify CMS genes in a gynodioecious plant, Silene vulgaris, we constructed mt transcriptomes and compared transcript levels and RNA editing patterns in floral bud tissue from female and hermaphrodite full siblings. The transcriptomes from female and hermaphrodite individuals were very similar overall with respect to variation in levels of transcript abundance across the genome, the extent of RNA editing, and the order in which RNA editing and intron splicing events occurred. We found only a single genomic region that was highly overexpressed and differentially edited in females relative to hermaphrodites. This region is not located near any other transcribed elements and lacks an open-reading frame (ORF) of even moderate size. To our knowledge, this transcript would represent the first non-coding mt RNA associated with CMS in plants and is, therefore, an important target for future functional validation studies.

De novo transcriptome assembly of heavy metal tolerant Silene dioica

Silene dioica is a dioecious plant of the family Caryophyllaceae. In the present study, we used Illumina sequencing technology (MiSeq) to sequence, de novo assembly and annotate the transcriptomes of male and female copper tolerant S. dioica individuals. We sequenced the normalized mRNA of roots, shoots, flower buds and flowers for each sex. Raw reads of the transcriptome assembly project for S. dioica male and female individual have been deposited in NCBI's Sequence Read Archive (SRA) database with the accession number SRP094611. The Trinity and Detonate program was used to de novo assembly 92,347 transcripts for male and 94,757 transcripts for female transcriptome. The assembled transcriptome sequences for S. dioica male and female individuals can be accessed at NCBI with the following accession numbers: GFCG00000000 (male); GFCH00000000 (female). The obtained transcriptomic data will be useful for further studies focusing on copper tolerance, comparative transcriptome analysis with other Silene species and sex chromosomes evolution.

Differential adaptation drives ecological speciation in campions (Silene): evidence from a multi-site transplant experiment

In order to investigate the role of differential adaptation for the evolution of reproductive barriers, we conducted a multi-site transplant experiment with the dioecious sister species Silene dioica and S. latifolia and their hybrids. Crosses within species as well as reciprocal first-generation (F₁ ) and second-generation (F₂ ) interspecific hybrids were transplanted into six sites, three within each species' habitat. Survival and flowering were recorded over 4 yr. At all transplant sites, the local species outperformed the foreign species, reciprocal F₁ hybrids performed intermediately and F₂ hybrids underperformed in comparison to F₁ hybrids (hybrid breakdown). Females generally had slightly higher cumulative fitness than males in both within- and between-species crosses and we thus found little evidence for Haldane's rule acting on field performance. The strength of selection against F₁ and F₂ hybrids as well as hybrid breakdown increased with increasing strength of habitat adaptation (i.e. the relative fitness difference between the local and the foreign species) across sites. Our results suggest that differential habitat adaptation led to ecologically dependent post-zygotic reproductive barriers and drives divergence and speciation in this Silene system.

Development of highly variable microsatellite markers for the tetraploid Silene stellata (Caryophyllaceae)

PREMISE OF THE STUDY

We designed and tested microsatellite markers for the North American native species Silene stellata (Caryophyllaceae) to investigate its population genetic structure and identify selection on floral design through male reproductive success.

METHODS AND RESULTS

A total of 153 candidate microsatellite loci were isolated based on next-generation sequencing. We identified 18 polymorphic microsatellite loci in three populations of S. stellata, with di- or trinucleotide repeats. Genotyping results showed the number of alleles per locus ranged from six to 45 and expected heterozygosity ranged from 0.511 to 0.951. Five of these loci were successfully amplified in S. virginica and S. caroliniana and were also polymorphic.

CONCLUSIONS

The microsatellite markers reported here provide a valuable tool for paternity analysis in S. stellata. They will also be useful for investigating the population genetic structures of S. stellata and related species.

Evolution of sex-biased gene expression in a dioecious plant

Separate sexes and sex-biased gene expression have repeatedly evolved in animals and plants, but the underlying changes in gene expression remain unknown. Here, we studied a pair of plant species, one in which separate sexes and sex chromosomes evolved recently and one which maintained hermaphrodite flowers resembling the ancestral state, to reconstruct expression changes associated with the evolution of dioecy. We found that sex-biased gene expression has evolved in autosomal and sex-linked genes in the dioecious species. Most expression changes relative to hermaphrodite flowers occurred in females rather than males, with higher and lower expression in females leading to female-biased and male-biased expression, respectively. Expression changes were more common in genes located on the sex chromosomes than the autosomes and led to feminization of the X chromosome and masculinization of the Y chromosome. Our results support a scenario in which sex-biased gene expression evolved during the evolution of dioecy to resolve intralocus sexual conflicts over the allocation of resources.

Strong phylogeographic co-structure between the anther-smut fungus and its white campion host

Although congruence between host and pathogen phylogenies has been extensively investigated, the congruence between host and pathogen genetic structures at the within-species level has received little attention. Using an unprecedented and comprehensive collection of associated plant-pathogen samples, we investigated the degree of congruence between the genetic structures across Europe of two evolutionary and ecological model organisms, the anther-smut pathogen Microbotryum lychnidis-dioicae and its host plant Silene latifolia. We demonstrated a significant and particularly strong level of host-pathogen co-structure, with three main genetic clusters displaying highly similar spatial ranges in Western Europe, Eastern Europe and Italy, respectively. Correcting for the geographical component of genetic variation, significant correlations were still found between the genetic distances of anther-smut and host populations. Inoculation experiments suggested plant local adaptation, at the cluster level, for resistance to pathogens. These findings indicate that the pathogen remained isolated in the same fragmented southern refugia as its host plant during the last glaciation, and that little long-distance dispersal has occurred since the recolonization of Europe for either the plant or the pathogen, despite their known ability to travel across continents. This, together with the inoculation results, suggests that coevolutionary and competitive processes may be drivers of host-pathogen co-structure.

Genetic variation in natural Melandrium album populations exposed to chronic ionizing radiation

The effect of radiation pollution on genetic variation in natural populations of Melandrium album was investigated at the head part of the East-Ural Radioactive Trace (EURT) and background areas. The highest genetic differentiation estimated using F _ST was revealed between compared pairs of the background and impact samples in populations of M. album. The highest rate of polymorphism was observed at the closest to nuclear accident, Impact-1 site. The unique alleles (Mdh-3¹⁰⁴, Pgi-2¹⁰⁶, Lap ¹⁰⁵, Mdh-2⁹⁶, and Dia ⁹⁴) were discovered at the EURT. Individuals from chronically low-level irradiated sites were genetically closer than to plants from background sites using Nadhdh locus. The increase of the frequency of unique homozygous and heterozygous genotypes was identified in populations of M. album growing under chronic radiation exposure conditions. The largest contribution to the group of unique heterozygous genotypes at the EURT was made by three loci - Lap, Pgi-2, and Nadhdh; the main role in interpopulation differentiation of samples was made by the alleles Sod-2¹¹⁵, Skdh ¹⁰⁰, and Nadhdh ¹⁰⁰. Our results provide evidence for the correlation between the increase of genetic variation other than the «genetic erosion» and chronic radiation exposure factor in natural plant populations.

Circadian rhythm of a Silene species favours nocturnal pollination and constrains diurnal visitation

Background and Aims Traits related to flower advertisement and reward sometimes vary in a circadian way, reflecting phenotypic specialization. However, specialized flowers are not necessarily restricted to specialized pollinators. This is the case of most Silene species, typically associated with diurnal or nocturnal syndromes of pollination but usually showing complex suites of pollinators. Methods A Silene species with mixed floral features between diurnal and nocturnal syndromes was used to test how petal opening, nectar production, scent emission and pollination success correlate in a circadian rhythm, and whether this is influenced by environmental conditions. The effect of diurnal and nocturnal visitation rates on plant reproductive success is also explored in three populations, including the effect of the pollinating seed predator Hadena sancta. Key Results The result showed that repeated petal opening at dusk was correlated with nectar secretion and higher scent production during the night. However, depending on environmental conditions, petals remain opened for a while in the morning, when nectar and pollen still were available. Pollen deposition was similarly effective at night and in the morning, but less effective in the afternoon. These results were consistent with field studies. Conclusions The circadian rhythm regulating floral attractiveness and reward in S. colorata is predominantly adapted to nocturnal flower visitors. However, favourable environmental conditions lengthen the optimal daily period of flower attraction and pollination towards morning. This allows the complementarity of day and night pollination. Diurnal pollination may help to compensate the plant reproductive success when nocturnal pollinators are scarce and when the net outcome of H. sancta shifts from mutualism to parasitism. These results suggest a functional mechanism explaining why the supposed nocturnal syndrome of many Silene species does not successfully predict their pollinator guilds.

Lower prevalence but similar fitness in a parasitic fungus at higher radiation levels near Chernobyl

Nuclear disasters at Chernobyl and Fukushima provide examples of effects of acute ionizing radiation on mutations that can affect the fitness and distribution of species. Here, we investigated the prevalence of Microbotryum lychnidis-dioicae, a pollinator-transmitted fungal pathogen of plants causing anther-smut disease in Chernobyl, its viability, fertility and karyotype variation, and the accumulation of nonsynonymous mutations in its genome. We collected diseased flowers of Silene latifolia from locations ranging by more than two orders of magnitude in background radiation, from 0.05 to 21.03 μGy/h. Disease prevalence decreased significantly with increasing radiation level, possibly due to lower pollinator abundance and altered pollinator behaviour. Viability and fertility, measured as the budding rate of haploid sporidia following meiosis from the diploid teliospores, did not vary with increasing radiation levels and neither did karyotype overall structure and level of chromosomal size heterozygosity. We sequenced the genomes of twelve samples from Chernobyl and of four samples collected from uncontaminated areas and analysed alignments of 6068 predicted genes, corresponding to 1.04 × 10(7)  base pairs. We found no dose-dependent differences in substitution rates (neither dN, dS, nor dN/dS). Thus, we found no significant evidence of increased deleterious mutation rates at higher levels of background radiation in this plant pathogen. We even found lower levels of nonsynonymous substitution rates in contaminated areas compared to control regions, suggesting that purifying selection was stronger in contaminated than uncontaminated areas. We briefly discuss the possibilities for a mechanistic basis of radio resistance in this nonmelanized fungus.

Transcriptome and Biochemical Analysis of a Flower Color Polymorphism in Silene littorea (Caryophyllaceae)

Flower color polymorphisms are widely used as model traits from genetics to ecology, yet determining the biochemical and molecular basis can be challenging. Anthocyanin-based flower color variations can be caused by at least 12 structural and three regulatory genes in the anthocyanin biosynthetic pathway (ABP). We use mRNA-Seq to simultaneously sequence and estimate expression of these candidate genes in nine samples of Silene littorea representing three color morphs (dark pink, light pink and white) across three developmental stages in hopes of identifying the cause of flower color variation. We identified 29 putative paralogs for the 15 candidate genes in the ABP. We assembled complete coding sequences for 16 structural loci and nine of ten regulatory loci. Among these 29 putative paralogs, we identified 622 SNPs, yet only nine synonymous SNPs in Ans had allele frequencies that differentiated pigmented petals (dark pink and light pink) from white petals. These Ans allele frequency differences were further investigated with an expanded sequencing survey of 38 individuals, yet no SNPs consistently differentiated the color morphs. We also found one locus, F3h1, with strong differential expression between pigmented and white samples (>42x). This may be caused by decreased expression of Myb1a in white petal buds. Myb1a in S. littorea is a regulatory locus closely related to Subgroup 7 Mybs known to regulate F3h and other loci in the first half of the ABP in model species. We then compare the mRNA-Seq results with petal biochemistry which revealed cyanidin as the primary anthocyanin and five flavonoid intermediates. Concentrations of three of the flavonoid intermediates were significantly lower in white petals than in pigmented petals (rutin, quercetin and isovitexin). The biochemistry results for rutin, quercetin, luteolin and apigenin are consistent with the transcriptome results suggesting a blockage at F3h, possibly caused by downregulation of Myb1a.

The large-X effect in plants: increased species divergence and reduced gene flow on the Silene X-chromosome

The disproportionately large involvement of the X-chromosome in the isolation of closely related species (the large-X effect) has been reported for many animals, where X-linked genes are mostly hemizygous in the heterogametic sex. The expression of deleterious recessive mutations is thought to drive the frequent involvement of the X-chromosome in hybrid sterility, as well as to reduce interspecific gene flow for X-linked genes. Here, we evaluate the role of the X-chromosome in the speciation of two closely related plant species - the white and red campions (Silene latifolia and S. dioica) - that hybridize widely across Europe. The two species evolved separate sexes and sex chromosomes relatively recently (~10(7)  years), and unlike most animal species, most X-linked genes have intact Y-linked homologs. We demonstrate that the X-linked genes show a very small and insignificant amount of interspecific gene flow, while gene flow involving autosomal loci is significant and sufficient to homogenize the gene pools of the two species. These findings are consistent with the hypothesis of the large-X effect in Silene and comprise the first report of this effect in plants. Nonhemizygosity of many X-linked genes in Silene males indicates that exposure of recessive mutations to selection may not be essential for the occurrence of the large-X effect. Several possible causes of the large-X effect in Silene are discussed.

Evidence for asymmetrical hybridization despite pre- and post-pollination reproductive barriers between two Silene species

Interspecific hybridization is widespread among plants; nevertheless, pre- and post-zygotic isolating mechanisms may maintain species integrity for interfertile species in sympatry despite some gene flow. Interspecific hybridization and potential isolating barriers were evaluated between co-flowering Silene asclepiadea and Silene yunnanensis in an alpine community in southwest China. We investigated morphological and molecular (nuclear microsatellites and chloroplast gene sequence) variation in sympatric populations of S. asclepiadea and S. yunnanensis. Additionally, we analyzed pollinator behaviour and compared reproductive success between the putative hybrids and their parental species. Both the molecular and morphological data indicate that there were putative natural hybrids in the field, with S. asclepiadae the ovule parent and S. yunnanensis the pollen parent. Bumblebees were the primary visitors to S. asclepiadae and putative hybrids, while butterflies were the primary visitors to S. yunnanensis Pollen production and viability were significantly lower in putative hybrids than the parental species. The direction of hybridization is quite asymmetric from S. yunnanensis to S. asclepiadea Protandry combined with later peak flowering of S. yunnanensis, and pollinator preference may have contributed to the asymmetric pattern of hybridization, but putative hybrids were rare. Our results thus suggest that despite gene flow, S. asclepiadea and S. yunnanensis can maintain species boundaries, perhaps as a result of floral isolation and low fecundity of the hybrids.

Two complete chloroplast genomes of white campion (Silene latifolia) from male and female individuals

In this study, we assembled two individuals' complete chloroplast genomes with one male and one female from the dioecious plant species white campion (Silene latifolia). The two chloroplast genomes have an identical composition with each other as a circular molecule of 150 931 bp in length, with an overall GC content of 36.4%. The genomes consist of a pair of inverted repeats (IRs) of 25 503 bp, separated by a large single-copy (LSC) region and a small single-copy (SSC) region of 82 708 and 17 217 bp, respectively. The genomes contain 111 single copy genes, including 77 protein-coding genes, 4 ribosomal RNA genes and 30 transfer RNA genes. In these protein-coding genes, eight genes (rpl16, rpoC1, rps16, petD, petB, ndhB, ndhA and atpF) contain a single intron and three genes (rps12, clpP and ycf3) contain two introns. The two newly sequenced chloroplast genomes provide valuable information for detecting polymorphisms within species and between sexes.

Fungal Infection Induces Sex-Specific Transcriptional Changes and Alters Sexual Dimorphism in the Dioecious Plant Silene latifolia

Sexual dimorphism, including differences in morphology, behavior and physiology between females and males, is widespread in animals and plants and is shaped by gene expression differences between the sexes. Such expression differences may also underlie sex-specific responses of hosts to pathogen infections, most notably when pathogens induce partial sex reversal in infected hosts. The genetic changes associated with sex-specific responses to pathogen infections on the one hand, and sexual dimorphism on the other hand, remain poorly understood. The dioecious White Campion (Silene latifolia) displays sexual dimorphism in floral traits and infection with the smut fungus Micobrotryum lychnidis-dioicae induces a partial sex reversal in females. We find strong sex-specific responses to pathogen infection and reduced sexual dimorphism in infected S. latifolia. This provides a direct link between pathogen-mediated changes in sex-biased gene expression and altered sexual dimorphism in the host. Expression changes following infection affected mainly genes with male-biased expression in healthy plants. In females, these genes were up-regulated, leading to a masculinization of the transcriptome. In contrast, infection in males was associated with down-regulation of these genes, leading to a demasculinization of the transcriptome. To a lesser extent, genes with female-biased expression in healthy plants were also affected in opposite directions in the two sexes. These genes were overall down-regulated in females and up-regulated in males, causing, respectively, a defeminization in infected females and a feminization of the transcriptome in infected males. Our results reveal strong sex-specific responses to pathogen infection in a dioecious plant and provide a link between pathogen-induced changes in sex-biased gene expression and sexual dimorphism.

Females and males differ from each other in many traits, including morphology, behavior and physiology. Differences in gene expression between the sexes, known as sex-biased gene expression, contribute to such sexual dimorphism. Here we characterize the responses of females and males of the dioecious plant Silene latifolia to infection with the anther smut fungus Micobrotryum lychnidis-dioicae. This fungus sterilizes the plant and induces a partial sex reversal in female hosts that form rudimentary stamens, thus allowing the fungus to transmit its spores via pollinators. Our comparisons of gene expression in healthy and infected plants reveal strong sex-specific responses to anther smut infection. Expression changes in females and males are in opposite directions and are associated with reduced sexual dimorphism between infected females and males. Our study reveals that infection with the anther smut fungus alters the extent of sex-biased gene expression in S. latifolia in a sex-specific manner and highlights how transcriptomic changes in females and males shape sexual dimorphism.

Rapid Y degeneration and dosage compensation in plant sex chromosomes

The nonrecombining regions of animal Y chromosomes are known to undergo genetic degeneration, but previous work has failed to reveal large-scale gene degeneration on plant Y chromosomes. Here, we uncover rapid and extensive degeneration of Y-linked genes in a plant species, Silene latifolia, that evolved sex chromosomes de novo in the last 10 million years. Previous transcriptome-based studies of this species missed unexpressed, degenerate Y-linked genes. To identify sex-linked genes, regardless of their expression, we sequenced male and female genomes of S. latifolia and integrated the genomic contigs with a high-density genetic map. This revealed that 45% of Y-linked genes are not expressed, and 23% are interrupted by premature stop codons. This contrasts with X-linked genes, in which only 1.3% of genes contained stop codons and 4.3% of genes were not expressed in males. Loss of functional Y-linked genes is partly compensated for by gene-specific up-regulation of X-linked genes. Our results demonstrate that the rate of genetic degeneration of Y-linked genes in S. latifolia is as fast as in animals, and that the evolutionary trajectories of sex chromosomes are similar in the two kingdoms.

Characterization of 14 microsatellite markers for Silene acaulis (Caryophyllaceae)(1)

PREMISE OF THE STUDY

Fifty candidate microsatellite markers, generated using 454 shotgun sequencing, were tested for the widespread arctic/alpine herb Silene acaulis (Caryophyllaceae).

METHODS AND RESULTS

Fourteen out of 50 markers resulted in polymorphic products with profiles that enabled interpretation. The numbers of alleles per locus ranged from two to six, and the expected heterozygosity per locus ranged from 0.06 to 0.68. Analysis of F0 and F1 samples proved that one allele was always inherited maternally. Four multiplex mixes have been developed.

CONCLUSIONS

Microsatellite markers for this species will be a valuable tool to study detailed small-scale genetic patterns in an arctic/alpine herb and to relate them to demographic parameters.

A glacial survivor of the alpine Mediterranean region: phylogenetic and phylogeographic insights into Silene ciliata Pourr. (Caryophyllaceae)

Silene ciliata Pourr. (Caryophyllaceae) is a species with a highly disjunct distribution which inhabits the alpine mountains of the Mediterranean Basin. We investigated the phylogeny and phylogeography of the species to (a) clarify the long-suggested division of S. ciliata into two subspecies, (b) evaluate its phylogenetic origin and (c) assess whether the species' diversification patterns were affected by the Mediterranean relief. For this purpose, we collected DNA from 25 populations of the species that inhabit the mountains of Portugal, Spain, France, Italy, former Yugoslav Republic of Macedonia, Bulgaria and Greece and studied the plastid regions rbcL, rps16 and trnL. Major intraspecific variation was supported by all analyses, while the possibility of the existence of more varieties or subspecies was not favoured. Plastid DNA (cpDNA) evidence was in accordance with the division of S. ciliata into the two subspecies, one spreading west (Iberian Peninsula and Central Massif) and the other east of the Alps region (Italian and Balkan Peninsula). This study proposes that the species' geographically disconnected distribution has probably derived from vicariance processes and from the Alps acting as a barrier to the species' dispersal. The monophyletic origin of the species is highly supported. cpDNA patterns were shown independent of the chromosome evolution in the populations and could have resulted from a combination of geographic factors providing links and barriers, climatic adversities and evolutionary processes that took place during Quaternary glaciations.

Identifying new sex-linked genes through BAC sequencing in the dioecious plant Silene latifolia

Background

Silene latifolia represents one of the best-studied plant sex chromosome systems. A new approach using RNA-seq data has recently identified hundreds of new sex-linked genes in this species. However, this approach is expected to miss genes that are either not expressed or are expressed at low levels in the tissue(s) used for RNA-seq. Therefore other independent approaches are needed to discover such sex-linked genes.

Results

Here we used 10 well-characterized S. latifolia sex-linked genes and their homologs in Silene vulgaris, a species without sex chromosomes, to screen BAC libraries of both species. We isolated and sequenced 4 Mb of BAC clones of S. latifolia X and Y and S. vulgaris genomic regions, which yielded 59 new sex-linked genes (with S. vulgaris homologs for some of them). We assembled sequences that we believe represent the tip of the Xq arm. These sequences are clearly not pseudoautosomal, so we infer that the S. latifolia X has a single pseudoautosomal region (PAR) on the Xp arm. The estimated mean gene density in X BACs is 2.2 times lower than that in S. vulgaris BACs, agreeing with the genome size difference between these species. Gene density was estimated to be extremely low in the Y BAC clones. We compared our BAC-located genes with the sex-linked genes identified in previous RNA-seq studies, and found that about half of them (those with low expression in flower buds) were not identified as sex-linked in previous RNA-seq studies. We compiled a set of ~70 validated X/Y genes and X-hemizygous genes (without Y copies) from the literature, and used these genes to show that X-hemizygous genes have a higher probability of being undetected by the RNA-seq approach, compared with X/Y genes; we used this to estimate that about 30 % of our BAC-located genes must be X-hemizygous. The estimate is similar when we use BAC-located genes that have S. vulgaris homologs, which excludes genes that were gained by the X chromosome.

Conclusions

Our BAC sequencing identified 59 new sex-linked genes, and our analysis of these BAC-located genes, in combination with RNA-seq data suggests that gene losses from the S. latifolia Y chromosome could be as high as 30 %, higher than previous estimates of 10-20 %.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1698-7) contains supplementary material, which is available to authorized users.

Diversity of sexual systems within different lineages of the genus Silene

Species and populations can be categorized by their sexual systems, depending on the spatial distribution of female and male reproductive structures within and among plants. Although a high diversity of sexual systems exists in Silene, their relative frequency at the genus and infrageneric level is unknown. Here, we carried out an extensive literature search for direct or indirect descriptions of sexual systems in Silene species. We found descriptions of sexual systems for 98 Silene species, where 63 and 35 correspond to the phylogenetically supported subgenera Silene and Behenantha, respectively. Hermaphroditism was the commonest sexual system (58.2 %), followed by dioecy (14.3 %), gynodioecy (13.3 %) and gynodioecy-gynomonoecy (i.e. hermaphroditic, female and gynomonoecious plants coexisting in the same population; 12.2 %). The presence of these sexual systems in both subgenera suggests their multiple origins. In 17 species, the description of sexual systems varied, and in most cases these differences corresponded to variations within or among populations. Interestingly, the poorly studied gynodioecy-gynomonoecy sexual system showed similar frequency to dioecy and gynodioecy in both subgenera. In addition, the incidence of gynodioecy-gynomonoecy was analysed in the species of section Psammophilae (Silene littorea, S. psammitis, S. adscendens and S. cambessedesii), in a survey of 26 populations across the distribution area of the species. The four species showed gynomonoecy-gynodioecy in most populations. Hermaphrodites were the most frequent morph, with a low number of females and gynomonoecious plants in all populations. The frequency of sexual morphs varied significantly among the studied populations but not among species. Female plants generally produced smaller numbers of flowers than hermaphroditic or gynomonoecious plants, and the percentages of female flowers per population were low. All these findings suggest that the gynodioecious-gynomonoecious sexual system in section Psammophilae is closer to hermaphroditism or gynomonoecy than gynodioecy.

Genetics of dioecy and causal sex chromosomes in plants

Dioecy (separate male and female individuals) ensures outcrossing and is more prevalent in animals than in plants. Although it is common in bryophytes and gymnosperms, only 5% of angiosperms are dioecious. In dioecious higher plants, flowers borne on male and female individuals are, respectively deficient in functional gynoecium and androecium. Dioecy is inherited via three sex chromosome systems: XX/XY, XX/X0 and WZ/ZZ, such that XX or WZ is female and XY, X0 or ZZ are males. The XX/XY system generates the rarer XX/X0 and WZ/ZZ systems. An autosome pair begets XY chromosomes. A recessive loss-of-androecium mutation (ana) creates X chromosome and a dominant gynoecium-suppressing (GYS) mutation creates Y chromosome. The ana/ANA and gys/GYS loci are in the sex-determining region (SDR) of the XY pair. Accumulation of inversions, deleterious mutations and repeat elements, especially transposons, in the SDR of Y suppresses recombination between X and Y in SDR, making Y labile and increasingly degenerate and heteromorphic from X. Continued recombination between X and Y in their pseudoautosomal region located at the ends of chromosomal arms allows survival of the degenerated Y and of the species. Dioecy is presumably a component of the evolutionary cycle for the origin of new species. Inbred hermaphrodite species assume dioecy. Later they suffer degenerate-Y-led population regression. Cross-hybridization between such extinguishing species and heterologous species, followed by genome duplication of segregants from hybrids, give rise to new species.

Recent and massive expansion of the mating-type-specific region in the smut fungus Microbotryum

The presence of large genomic regions with suppressed recombination (SR) is a key shared property of some sex- and mating-type determining (mat) chromosomes identified to date in animals, plants, and fungi. Why such regions form and how they evolve remain central questions in evolutionary genetics. The smut fungus Microbotryum lychnis-dioicae is a basidiomycete fungus in which dimorphic mat chromosomes have been reported, but the size, age, and evolutionary dynamics of the SR region remains unresolved. To identify the SR region in M. lychnis-dioicae and to study its evolution, we sequenced 12 genomes (6 per mating type) of this species and identified the genomic contigs that show fixed sequence differences between the mating types. We report that the SR region spans more than half of the mat chromosome (>2.3 Mbp) and that it is of very recent origin (∼2 × 10(6) years) as the average sequence divergence between mating types was only 2% in the SR region. This contrasts with a much higher divergence in and around the mating-type determining pheromone receptor locus in the SR, suggesting a recent and massive expansion of the SR region. Our results comprise the first reported case of recent massive SR expansion documented in a basidiomycete fungus.

Fragaria: a genus with deep historical roots and ripe for evolutionary and ecological insights

The cultivated strawberry, Fragaria ×ananassa, is one of the youngest domesticated plants. Its 18th century origin via hybridization in Europe between the North American F. virginiana and the South American F. chiloensis was documented by the botanist Antoine Nicolas Duchesne. His 1766 "Natural History of Strawberries" is an extraordinary work that integrates fundamental discoveries on the biology, ecology, and phylogeny of Fragaria with applied information on cultivation and ethnobotanical uses, serving as an inspiration for current research in the genus. Fragaria species exhibit the full range of sexual systems in the gynodioecy pathway from hermaphroditism to dioecy (and back again), as well as variation in self-compatibility, and evidence of sex chromosomes with female heterogamety. The genus is also characterized by interspecific hybridization and polyploidy, with a natural range of ploidy levels from diploids to decaploids. This biological diversity, combined with the availability of genomic resources and the ease of growing and experimenting with the plants, makes Fragaria a very attractive system for ecological and evolutionary genomics. The goal of this review is to introduce Fragaria as a model genus and to provide a roadmap for future integrative research. These research directions will deepen our understanding of the ecological and evolutionary context that shaped the ancestors of the cultivated strawberry, not only providing information that can be applied to efforts to shape the future of this important fruit crop but also our understanding of key transitions in plant evolution.

Marginal likelihood estimate comparisons to obtain optimal species delimitations in Silene sect. Cryptoneurae (Caryophyllaceae)

Coalescent-based inference of phylogenetic relationships among species takes into account gene tree incongruence due to incomplete lineage sorting, but for such methods to make sense species have to be correctly delimited. Because alternative assignments of individuals to species result in different parametric models, model selection methods can be applied to optimise model of species classification. In a Bayesian framework, Bayes factors (BF), based on marginal likelihood estimates, can be used to test a range of possible classifications for the group under study. Here, we explore BF and the Akaike Information Criterion (AIC) to discriminate between different species classifications in the flowering plant lineage Silene sect. Cryptoneurae (Caryophyllaceae). We estimated marginal likelihoods for different species classification models via the Path Sampling (PS), Stepping Stone sampling (SS), and Harmonic Mean Estimator (HME) methods implemented in BEAST. To select among alternative species classification models a posterior simulation-based analog of the AIC through Markov chain Monte Carlo analysis (AICM) was also performed. The results are compared to outcomes from the software BP&P. Our results agree with another recent study that marginal likelihood estimates from PS and SS methods are useful for comparing different species classifications, and strongly support the recognition of the newly described species S. ertekinii.

Determinants of genetic structure in a nonequilibrium metapopulation of the plant Silene latifolia

Population genetic differentiation will be influenced by the demographic history of populations, opportunities for migration among neighboring demes and founder effects associated with repeated extinction and recolonization. In natural populations, these factors are expected to interact with each other and their magnitudes will vary depending on the spatial distribution and age structure of local demes. Although each of these effects has been individually identified as important in structuring genetic variance, their relative magnitude is seldom estimated in nature. We conducted a population genetic analysis in a metapopulation of the angiosperm, Silene latifolia, from which we had more than 20 years of data on the spatial distribution, demographic history, and extinction and colonization of demes. We used hierarchical Bayesian methods to disentangle which features of the populations contributed to among population variation in allele frequencies, including the magnitude and direction of their effects. We show that population age, long-term size and degree of connectivity all combine to affect the distribution of genetic variance; small, recently-founded, isolated populations contributed most to increase F_ST in the metapopulation. However, the effects of population size and population age are best understood as being modulated through the effects of connectivity to other extant populations, i.e. F_ST diminishes as populations age, but at a rate that depends how isolated the population is. These spatial and temporal correlates of population structure give insight into how migration, founder effect and within-deme genetic drift have combined to enhance and restrict genetic divergence in a natural metapopulation.

Isolating fungal pathogens from a dynamic disease outbreak in a native plant population to establish plant-pathogen bioassays for the ecological model plant Nicotiana attenuata

The wild tobacco species Nicotiana attenuata has been intensively used as a model plant to study its interaction with insect herbivores and pollinators in nature, however very little is known about its native pathogen community. We describe a fungal disease outbreak in a native N. attenuata population comprising 873 plants growing in an area of about 1500 m2. The population was divided into 14 subpopulations and disease symptom development in the subpopulations was monitored for 16 days, revealing a waxing and waning of visible disease symptoms with some diseased plants recovering fully. Native fungal N. attenuata pathogens were isolated from diseased plants, characterized genetically, chemotaxonomically and morphologically, revealing several isolates of the ascomycete genera Fusarium and Alternaria, that differed in the type and strength of the disease symptoms they caused in bioassays on either detached leaves or intact soil-grown plants. These isolates and the bioassays will empower the study of N. attenuata-pathogen interactions in a realistic ecological context.

Diversity of Secondary Metabolites in the Genus Silene L. (Caryophyllaceae)—Structures, Distribution, and Biological Properties

The genus Silene (family Caryophyllaceae) comprises more than 700 species, which are widely distributed in temperate zones of the Northern Hemisphere, but are also present in Africa and have been introduced in other continents. Silene produces a high diversity of secondary metabolites and many of them show interesting biological and pharmacological activities. More than 450 compounds have been isolated; important classes include phytoecdysteroids (which mimic insect molting hormones), triterpene saponins (with detergent properties), volatiles, other terpenoids and phenolics. This review focusses on the phytochemical diversity, distribution of Silene secondary metabolites and their biological activities.

Interactions between a pollinating seed predator and its host plant: the role of environmental context within a population

Plant-insect interactions often are important for plant reproduction, but the outcome of these interactions may vary with environmental context. Pollinating seed predators have positive and negative effects on host plant reproduction, and the interaction outcome is predicted to vary with density or abundance of the partners. We studied the interaction between Silene stellata, an herbaceous perennial, and Hadena ectypa, its specialized pollinating seed predator. Silene stellata is only facultatively dependent upon H. ectypa for pollination because other nocturnal moth co-pollinators are equally effective at pollen transfer. We hypothesized that for plants without conspecific neighbors, H. ectypa would have higher visitation rates compared to co-pollinators, and the plants would experience lower levels of H. ectypa pollen deposition. We predicted similar oviposition throughout the study site but greater H. ectypa predation in the area without conspecific neighbors compared to plants embedded in a naturally high density area. We found that H. ectypa had consistently higher visitation than moth co-pollinators in all host plant contexts. However, H. ectypa pollinator importance declined in areas with low conspecific density because of reduced pollen deposition, resulting in lower seed set. Conversely, oviposition was similar across the study site independent of host plant density. Greater likelihood of very high fruit predation combined with lower pollination by H. ectypa resulted in reduced S. stellata female reproductive success in areas with low conspecific density. Our results demonstrate local context dependency of the outcomes of pollinating seed predator interactions with conspecific host plant density within a population.

Do flower color and floral scent of silene species affect host preference of Hadena bicruris, a seed-eating pollinator, under field conditions?

Specialization in plant–insect interactions is an important driver of evolutionary divergence; yet, plant traits mediating such interactions are poorly understood. In this study, we investigated how flower color and floral scent are related to seed predation by a seed-eating pollinator. We used field-transplanted recombinant F₂ hybrids between Silene latifolia and S. dioica that are the preferred and alternative hosts of the moth Hadena bicruris and crosses within these species for comparison. We scored seed predation and flower color and analyzed floral scent. Pinker S. dioica-like flowers and emission of α-pinene decreased the odds of seed predation while emission of benzyl acetate and 6-methyl-5-hepten-2-one increased the odds of seed predation. Emission of these compounds did not differ significantly between the two Silene species. Our results suggest that flower color plays an important role in the specific interaction of H. bicruris with its preferred host S. latifolia. The compounds α-pinene, benzyl acetate and 6-methyl-5-hepten-2-one could represent non-specific deterrents and attractants to ovipositing moths. Alternatively, emission of these compounds could be related to herbivory or pathogen attack and act as a signal for host quality. This would weaken the predictability of the plant's costs and benefits of the interaction and act to maintain an imperfect degree of specialization.

Identification of internal reference genes for gene expression normalization between the two sexes in dioecious white Campion

Quantitative real time (qRT)-PCR is a precise and efficient method for studying gene expression changes between two states of interest, and is frequently used for validating interesting gene expression patterns in candidate genes initially identified in genome-wide expression analyses, such as RNA-seq experiments. For an adequate normalisation of qRT-PCR data, it is essential to have reference genes available whose expression intensities are constant among the different states of interest. In this study we present and validate a catalogue of traditional and newly identified reference genes that were selected from RNA-seq data from multiple individuals from the dioecious plant Silene latifolia with the aim of studying gene expression differences between the two sexes in both reproductive and vegetative tissues. The catalogue contains more than 15 reference genes with both stable expression intensities and a range of expression intensities in flower buds and leaf tissues. These reference genes were used to normalize expression differences between reproductive and vegetative tissues in eight candidate genes with sex-biased expression. Our results suggest a trend towards a reduced sex-bias in sex-linked gene expression in vegetative tissues. In this study, we report on the systematic identification and validation of internal reference genes for adequate normalization of qRT-PCR-based analyses of gene expression differences between the two sexes in S. latifolia. We also show how RNA-seq data can be used efficiently to identify suitable reference genes in a wide diversity of species.

A quantitative genetic signature of senescence in a short-lived perennial plant

The evolution of senescence (the physiological decline of organisms with age) poses an apparent paradox because it represents a failure of natural selection to increase the survival and reproductive performance of organisms. The paradox can be resolved if natural selection becomes less effective with age, because the death of postreproductive individuals should have diminished effects on Darwinian fitness [1, 2]. A substantial body of empirical work is consistent with this prediction for animals, which transmit their genes to progeny via an immortal germline. However, such evidence is still lacking in plants, which lack a germline and whose reproduction is diffuse and modular across the soma. Here, we provide experimental evidence for a genetic basis of senescence in the short-lived perennial plant Silene latifolia. Our pedigree-based analysis revealed a marked increase with age in the additive genetic variance of traits closely associated with fitness. This result thus extends to plants the quantitative genetic support for the evolutionary theory of senescence.

The roles of female and hermaphroditic flowers in the gynodioecious-gynomonoecious Silene littorea: insights into the phenology of sex expression

Some gynodioecious species have intermediate individuals that bear both female and hermaphroditic flowers. This phenomenon is known as a gynodioecious-gynomonoecious sexual system. Gender expression in such species has received little attention in the past, and the phenologies of male and female functions have also yet to be explored. In this study, we examined variations in gender patterns, their effects on female reproductive success and sex expression in depth throughout the flowering period in two populations. The studied populations of Silene littorea contained mostly gynomonoecious plants and the number of pure females was very low. The gynomonoecious plants showed high variability in the total proportion of female flowers. In addition, the proportion of female flowers in each plant varied widely across the flowering season. Although there was a trend towards maleness, our measures of functional gender suggested that most plants transmit their genes via both pollen and ovules. Fruit set and seed set were not significantly different among populations; in contrast, flower production significantly varied between the two populations - and among plants - with consequent variation in total seed production. Conversely, gender and sex expression were similar in both populations. Plants with higher phenotypic femaleness did not have higher fruit set, seed set or total female fecundity. The mating environment fluctuated little across the flowering period, but fluctuations were higher in the population with low flower production. We therefore conclude that the high proportion of gynomonoecious individuals in our studied populations of S. littorea may be advantageous for the species, providing the benefits of both hermaphroditic and female flowers.

Patterns of cyto-nuclear linkage disequilibrium in Silene latifolia: genomic heterogeneity and temporal stability

Non-random association of alleles in the nucleus and cytoplasmic organelles, or cyto-nuclear linkage disequilibrium (LD), is both an important component of a number of evolutionary processes and a statistical indicator of others. The evolutionary significance of cyto-nuclear LD will depend on both its magnitude and how stable those associations are through time. Here, we use a longitudinal population genetic data set to explore the magnitude and temporal dynamics of cyto-nuclear disequilibria through time. We genotyped 135 and 170 individuals from 16 and 17 patches of the plant species Silene latifolia in Southwestern VA, sampled in 1993 and 2008, respectively. Individuals were genotyped at 14 highly polymorphic microsatellite markers and a single-nucleotide polymorphism (SNP) in the mitochondrial gene, atp1. Normalized LD (D') between nuclear and cytoplasmic loci varied considerably depending on which nuclear locus was considered (ranging from 0.005-0.632). Four of the 14 cyto-nuclear associations showed a statistically significant shift over approximately seven generations. However, the overall magnitude of this disequilibrium was largely stable over time. The observed origin and stability of cyto-nuclear LD is most likely caused by the slow admixture between anciently diverged lineages within the species' newly invaded range, and the local spatial structure and metapopulation dynamics that are known to structure genetic variation in this system.