Plant sexual reproduction during climate change: gene function in natura studied by ecological and evolutionary systems biology

Uncovering genes involved in pollinator-driven mating system shifts and selfing syndrome evolution in Brassica rapa

Shifts in pollinator occurrence and their pollen transport effectiveness drive the evolution of mating systems in flowering plants. Understanding the genomic basis of these changes is essential for predicting the persistence of a species under environmental changes. We investigated the genomic changes in Brassica rapa over nine generations of pollination by hoverflies associated with rapid morphological evolution toward the selfing syndrome. We combined a genotyping-by-sequencing (GBS) approach with a genome-wide association study (GWAS) to identify candidate genes, and assessed their functional role in the observed morphological changes by studying mutations of orthologous genes in the model plant Arabidopsis thaliana. We found 31 candidate genes involved in a wide range of functions from DNA/RNA binding to transport. Our functional assessment of orthologous genes in A. thaliana revealed that two of the identified genes in B. rapa are involved in regulating the size of floral organs. We found a protein kinase superfamily protein involved in petal width, an important trait in plant attractiveness to pollinators. Moreover, we found a histone lysine methyltransferase (HKMT) associated with stamen length. Altogether, our study shows that hoverfly pollination leads to rapid evolution toward the selfing syndrome mediated by polygenic changes.

Flowering time: From physiology, through genetics to mechanism

Plant species have evolved different requirements for environmental/endogenous cues to induce flowering. Originally, these varying requirements were thought to reflect the action of different molecular mechanisms. Thinking changed when genetic and molecular analysis in Arabidopsis thaliana revealed that a network of environmental and endogenous signaling input pathways converge to regulate a common set of "floral pathway integrators." Variation in the predominance of the different input pathways within a network can generate the diversity of requirements observed in different species. Many genes identified by flowering time mutants were found to encode general developmental and gene regulators, with their targets having a specific flowering function. Studies of natural variation in flowering were more successful at identifying genes acting as nodes in the network central to adaptation and domestication. Attention has now turned to mechanistic dissection of flowering time gene function and how that has changed during adaptation. This will inform breeding strategies for climate-proof crops and help define which genes act as critical flowering nodes in many other species.

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.

UV RESISTANCE LOCUS 8-Mediated UV-B Response Is Required Alongside CRYPTOCHROME 1 for Plant Survival in Sunlight under Field Conditions

As sessile, photoautotrophic organisms, plants are subjected to fluctuating sunlight that includes potentially detrimental ultraviolet-B (UV-B) radiation. Experiments under controlled conditions have shown that the UV-B photoreceptor UV RESISTANCE LOCUS 8 (UVR8) controls acclimation and tolerance to UV-B in Arabidopsis thaliana; however, its long-term impact on plant fitness under naturally fluctuating environments remain poorly understood. Here, we quantified the survival and reproduction of different Arabidopsis mutant genotypes under diverse field and laboratory conditions. We found that uvr8 mutants produced more fruits than wild type when grown in growth chambers under artificial low-UV-B conditions but not under natural field conditions, indicating a fitness cost in the absence of UV-B stress. Importantly, independent double mutants of UVR8 and the blue light photoreceptor gene CRYPTOCHROME 1 (CRY1) in two genetic backgrounds showed a drastic reduction in fitness in the field. Experiments with UV-B attenuation in the field and with supplemental UV-B in growth chambers demonstrated that UV-B caused the cry1 uvr8 conditional lethal phenotype. Using RNA-seq data of field-grown single and double mutants, we explicitly identified genes showing significant statistical interaction of UVR8 and CRY1 mutations in the presence of UV-B in the field. They were enriched in Gene Ontology categories related to oxidative stress, photoprotection and DNA damage repair in addition to UV-B response. Our study demonstrates the functional importance of the UVR8-mediated response across life stages in natura, which is partially redundant with that of cry1. Moreover, these data provide an integral picture of gene expression associated with plant responses under field conditions.

Seasonal pigment fluctuation in diploid and polyploid Arabidopsis revealed by machine learning-based phenotyping method PlantServation

Long-term field monitoring of leaf pigment content is informative for understanding plant responses to environments distinct from regulated chambers but is impractical by conventional destructive measurements. We developed PlantServation, a method incorporating robust image-acquisition hardware and deep learning-based software that extracts leaf color by detecting plant individuals automatically. As a case study, we applied PlantServation to examine environmental and genotypic effects on the pigment anthocyanin content estimated from leaf color. We processed >4 million images of small individuals of four Arabidopsis species in the field, where the plant shape, color, and background vary over months. Past radiation, coldness, and precipitation significantly affected the anthocyanin content. The synthetic allopolyploid A. kamchatica recapitulated the fluctuations of natural polyploids by integrating diploid responses. The data support a long-standing hypothesis stating that allopolyploids can inherit and combine the traits of progenitors. PlantServation facilitates the study of plant responses to complex environments termed "in natura".

Genome-wide association study of aphid abundance highlights a locus affecting plant growth and flowering in Arabidopsis thaliana

Plant life-history traits, such as size and flowering, contribute to shaping variation in herbivore abundance. Although plant genes involved in physical and chemical traits have been well studied, less is known about the loci linking plant life-history traits and herbivore abundance. Here, we conducted a genome-wide association study (GWAS) of aphid abundance in a field population of Arabidopsis thaliana. This GWAS of aphid abundance detected a relatively rare but significant variant on the third chromosome of A. thaliana, which was also suggestively but non-significantly associated with the presence or absence of inflorescence. Out of candidate genes near this significant variant, a mutant of a ribosomal gene (AT3G13882) exhibited slower growth and later flowering than a wild type under laboratory conditions. A no-choice assay with the turnip aphid, Lipaphis erysimi, found that aphids were unable to successfully establish on the mutant. Our GWAS of aphid abundance unexpectedly found a locus affecting plant growth and flowering.

Robustness and the generalist niche of polyploid species: Genome shock or gradual evolution?

The prevalence of polyploidy in wild and crop species has stimulated debate over its evolutionary advantages and disadvantages. Previous studies have focused on changes occurring at the polyploidization events, including genome-wide changes termed "genome shock," as well as ancient polyploidy. Recent bioinformatics advances and empirical studies of Arabidopsis and wheat relatives are filling a research gap: the functional evolutionary study of polyploid species using RNA-seq, DNA polymorphism, and epigenomics. Polyploid species can become generalists in natura through environmental robustness by inheriting and merging parental stress responses. Their evolvability is enhanced by mutational robustness working on inherited standing variation. The identification of key genes responsible for gradual adaptive evolution will encourage synthetic biological approaches to transfer polyploid advantages to other species.

ATTED-II v11: A Plant Gene Coexpression Database Using a Sample Balancing Technique by Subagging of Principal Components

ATTED-II (https://atted.jp) is a gene coexpression database for nine plant species based on publicly available RNAseq and microarray data. One of the challenges in constructing condition-independent coexpression data based on publicly available gene expression data is managing the inherent sampling bias. Here, we report ATTED-II version 11, wherein we adopted a coexpression calculation methodology to balance the samples using principal component analysis and ensemble calculation. This approach has two advantages. First, omitting principal components with low contribution rates reduces the main contributors of noise. Second, balancing large differences in contribution rates enables considering various sample conditions entirely. In addition, based on RNAseq- and microarray-based coexpression data, we provide species-representative, integrated coexpression information to enhance the efficiency of interspecies comparison of the coexpression data. These coexpression data are provided as a standardized z-score to facilitate integrated analysis with different data sources. We believe that with these improvements, ATTED-II is more valuable and powerful for supporting interspecies comparative studies and integrated analyses using heterogeneous data.

Field transcriptome analysis reveals a molecular mechanism for cassava-flowering in a mountainous environment in Southeast Asia

The field survey in this article showed in 'KU50', a popular variety and late-branching type of cassava in Southeast Asia, that flowering rarely occurs in normal-field conditions in Southeast Asia but is strongly induced in the dry season in the mountainous region. Flowering time is correlated with the expression patterns of MeFT1 and homologs of Arabidopsis GI, PHYA, and NF-Ys. Cassava (Manihot esculenta Crantz) is a tropical crop that is propagated vegetatively rather than sexually by seed. Flowering rarely occurs in the erect-type variety grown in Southeast Asia, but it is known that cassava produces flowers every year in mountainous regions. Data pertaining to the effect of environmental factors on flowering time and gene expression in cassava, however, is limited. The aim of the present study was to determine the kinds of environmental conditions that regulate flowering time in cassava and the underlying molecular mechanisms. The flowering status of KU50, a popular variety in Southeast Asia and late-branching type of cassava, was monitored in six fields in Vietnam and Cambodia. At non-flowering and flowering field locations in North Vietnam, the two FLOWERING LOCUS T (FT)-like genes, MeFT1 and MeFT2, were characterized by qPCR, and the pattern of expression of flowering-related genes and genes responsive to environmental signals were analyzed by using RNA sequencing data from time-series samples. Results indicate that cassava flowering was induced in the dry season in the mountain region, and that flowering time was correlated with the expression of MeFT1, and homologs of Arabidopsis GI, PHYA, and NF-Ys. Based upon these data, we hypothesize that floral induction in cassava is triggered by some conditions present in the mountain regions during the dry season.

Fine-scale empirical data on niche divergence and homeolog expression patterns in an allopolyploid and its diploid progenitor species

Polyploidization is pervasive in plants, but little is known about the niche divergence of wild allopolyploids (species that harbor polyploid genomes originating from different diploid species) relative to their diploid progenitor species and the gene expression patterns that may underlie such ecological divergence. We conducted a fine-scale empirical study on habitat and gene expression of an allopolyploid and its diploid progenitors. We quantified soil properties and light availability of habitats of an allotetraploid Cardamine flexuosa and its diploid progenitors Cardamine amara and Cardamine hirsuta in two seasons. We analyzed expression patterns of genes and homeologs (homeologous gene copies in allopolyploids) using RNA sequencing. We detected niche divergence between the allopolyploid and its diploid progenitors along water availability gradient at a fine scale: the diploids in opposite extremes and the allopolyploid in a broader range between diploids, with limited overlap with diploids at both ends. Most of the genes whose homeolog expression ratio changed among habitats in C. flexuosa varied spatially and temporally. These findings provide empirical evidence for niche divergence between an allopolyploid and its diploid progenitor species at a fine scale and suggest that divergent expression patterns of homeologs in an allopolyploid may underlie its persistence in diverse habitats.

Experimental and Field Data Support Range Expansion in an Allopolyploid Arabidopsis Owing to Parental Legacy of Heavy Metal Hyperaccumulation

Empirical evidence is limited on whether allopolyploid species combine or merge parental adaptations to broaden habitats. The allopolyploid Arabidopsis kamchatica is a hybrid of the two diploid parents Arabidopsis halleri and Arabidopsis lyrata. A. halleri is a facultative heavy metal hyperaccumulator, and may be found in cadmium (Cd) and zinc (Zn) contaminated environments, as well as non-contaminated environments. A. lyrata is considered non-tolerant to these metals, but can be found in serpentine habitats. Therefore, the parents have adaptation to different environments. Here, we measured heavy metals in soils from native populations of A. kamchatica. We found that soil Zn concentration of nearly half of the sampled 40 sites was higher than the critical toxicity level. Many of the sites were near human construction, suggesting adaptation of A. kamchatica to artificially contaminated soils. Over half of the A. kamchatica populations had >1,000 μg g^–1 Zn in leaf tissues. Using hydroponic treatments, most genotypes accumulated >3,000 μg g^–1 Zn, with high variability among them, indicating substantial genetic variation in heavy metal accumulation. Genes involved in heavy metal hyperaccumulation showed an expression bias in the A. halleri-derived homeolog in widely distributed plant genotypes. We also found that two populations were found growing on serpentine soils. These data suggest that A. kamchatica can inhabit a range of both natural and artificial soil environments with high levels of ions that either of the parents specializes and that it can accumulate varying amount of heavy metals. Our field and experimental data provide a compelling example of combining genetic toolkits for soil adaptations to expand the habitat of an allopolyploid species.

Fermentation Ability of Gut Microbiota of Wild Japanese Macaques in the Highland and Lowland Yakushima: In Vitro Fermentation Assay and Genetic Analyses

Wild Japanese macaques (Macaca fuscata Blyth) living in the highland and lowland areas of Yakushima are known to have different diets, with highland individuals consuming more leaves. We aim to clarify whether and how these differences in diet are also reflected by gut microbial composition and fermentation ability. Therefore, we conduct an in vitro fermentation assay using fresh feces from macaques as inoculum and dry leaf powder of Eurya japonica Thunb. as a substrate. Fermentation activity was higher for feces collected in the highland, as evidenced by higher gas and butyric acid production and lower pH. Genetic analysis indicated separation of highland and lowland in terms of both community structure and function of the gut microbiota. Comparison of feces and suspension after fermentation indicated that the community structure changed during fermentation, and the change was larger for lowland samples. Analysis of the 16S rRNA V3-V4 barcoding region of the gut microbiota showed that community structure was clearly clustered between the two areas. Furthermore, metagenomic analysis indicated separation by gene and pathway abundance patterns. Two pathways (glycogen biosynthesis I and D-galacturonate degradation I) were enriched in lowland samples, possibly related to the fruit-eating lifestyle in the lowland. Overall, we demonstrated that the more leaf-eating highland Japanese macaques harbor gut microbiota with higher leaf fermentation ability compared with the more fruit-eating lowland ones. Broad, non-specific taxonomic and functional gut microbiome differences suggest that this pattern may be driven by a complex interplay between many taxa and pathways rather than single functional traits.

Adaptive reduction of male gamete number in the selfing plant Arabidopsis thaliana

The number of male gametes is critical for reproductive success and varies between and within species. The evolutionary reduction of the number of pollen grains encompassing the male gametes is widespread in selfing plants. Here, we employ genome-wide association study (GWAS) to identify underlying loci and to assess the molecular signatures of selection on pollen number-associated loci in the predominantly selfing plant Arabidopsis thaliana. Regions of strong association with pollen number are enriched for signatures of selection, indicating polygenic selection. We isolate the gene REDUCED POLLEN NUMBER1 (RDP1) at the locus with the strongest association. We validate its effect using a quantitative complementation test with CRISPR/Cas9-generated null mutants in nonstandard wild accessions. In contrast to pleiotropic null mutants, only pollen numbers are significantly affected by natural allelic variants. These data support theoretical predictions that reduced investment in male gametes is advantageous in predominantly selfing species.

Parallel evolution of dominant pistil-side self-incompatibility suppressors in Arabidopsis

Selfing is a frequent evolutionary trend in angiosperms, and is a suitable model for studying the recurrent patterns underlying adaptive evolution. Many plants avoid self-fertilization by physiological processes referred to as self-incompatibility (SI). In the Brassicaceae, direct and specific interactions between the male ligand SP11/SCR and the female receptor kinase SRK are required for the SI response. Although Arabidopsis thaliana acquired autogamy through loss of these genes, molecular evolution contributed to the spread of self-compatibility alleles requires further investigation. We show here that in this species, dominant SRK silencing genes have evolved at least twice. Different inverted repeat sequences were found in the relic SRK region of the Col-0 and C24 strains. Both types of inverted repeats suppress the functional SRK sequence in a dominant fashion with different target specificities. It is possible that these dominant suppressors of SI contributed to the rapid fixation of self-compatibility in A. thaliana.

In Brassicaceae, interaction between the pollen-derived peptide ligand SP11 and the pistil-expressed receptor kinase SRK leads to self-incompatibility. Here the authors provide evidence that in Arabidopsis dominant self-compatibility inducers evolved at least twice via insertion of inverted repeats in the SRK locus.

Transcriptional Variation in Glucosinolate Biosynthetic Genes and Inducible Responses to Aphid Herbivory on Field-Grown Arabidopsis thaliana

Recently, increasing attempts have been made to understand how plant genes function in natura. In this context, transcriptional profiles represent plant physiological status in response to environmental stimuli. Herein, we combined high-throughput RNA-Seq with insect survey data on 19 accessions of Arabidopsis thaliana grown at a field site in Switzerland. We found that genes with the gene ontology (GO) annotations of "glucosinolate biosynthetic process" and "response to insects" were most significantly enriched, and the expression of these genes was highly variable among plant accessions. Nearly half of the total expression variation in the glucosinolate biosynthetic genes (AOPs, ESM1, ESP, and TGG1) was explained by among-accession variation. Of these genes, the expression level of AOP3 differed among Col-0 accession individuals depending on the abundance of the mustard aphid (Lipaphis erysimi). We also found that the expression of the major cis-jasmone activated gene CYP81D11 was positively correlated with the number of flea beetles (Phyllotreta striolata and Phyllotreta atra). Combined with the field RNA-Seq data, bioassays confirmed that AOP3 was up-regulated in response to attack by mustard aphids. The combined results from RNA-Seq and our ecological survey illustrate the feasibility of using field transcriptomics to detect an inducible defense, providing a first step towards an in natura understanding of biotic interactions involving phenotypic plasticity.

Seasonal Stability and Dynamics of DNA Methylation in Plants in a Natural Environment

DNA methylation has been considered a stable epigenetic mark but may respond to fluctuating environments. However, it is unclear how they behave in natural environments. Here, we analyzed seasonal patterns of genome-wide DNA methylation in a single clone from a natural population of the perennial Arabidopsis halleri. The genome-wide pattern of DNA methylation was primarily stable, and most of the repetitive regions were methylated across the year. Although the proportion was small, we detected seasonally methylated cytosines (SeMCs) in the genome. SeMCs in the CHH context were detected predominantly at repetitive sequences in intergenic regions. In contrast, gene-body CG methylation (gbM) itself was generally stable across seasons, but the levels of gbM were positively associated with seasonal stability of RNA expression of the genes. These results suggest the existence of two distinct aspects of DNA methylation in natural environments: sources of epigenetic variation and epigenetic marks for stable gene expression.

A stigmatic gene confers interspecies incompatibility in the Brassicaceae

Pre-zygotic interspecies incompatibility in angiosperms is a male-female relationship that inhibits the formation of hybrids between two species. Here, we report on the identification of STIGMATIC PRIVACY 1 (SPRI1), an interspecies barrier gene in Arabidopsis thaliana. We show that the rejection activity of this stigma-specific plasma membrane protein is effective against distantly related Brassicaceae pollen tubes and is independent of self-incompatibility. Point-mutation experiments and functional tests of synthesized hypothetical ancestral forms of SPRI1 suggest evolutionary decay of SPRI1-controlled interspecies incompatibility in self-compatible A. thaliana. Hetero-pollination experiments indicate that SPRI1 ensures intraspecific fertilization in the pistil when pollen from other species are present. Our study supports the idea that SPRI1 functions as a barrier mechanism that permits entrance of pollen with an intrinsic signal from self species.

Plant trichomes and a single gene GLABRA1 contribute to insect community composition on field-grown Arabidopsis thaliana

BACKGROUND

Genetic variation in plants alters insect abundance and community structure in the field; however, little is known about the importance of a single gene among diverse plant genotypes. In this context, Arabidopsis trichomes provide an excellent system to discern the roles of natural variation and a key gene, GLABRA1, in shaping insect communities. In this study, we transplanted two independent glabrous mutants (gl1-1 and gl1-2) and 17 natural accessions of Arabidopsis thaliana to two localities in Switzerland and Japan.

RESULTS

Fifteen insect species inhabited the plant accessions, with the insect community composition significantly attributed to variations among plant accessions. The total abundance of leaf-chewing herbivores was negatively correlated with trichome density at both field sites, while glucosinolates had variable effects on leaf chewers between the sites. Interestingly, there was a parallel tendency for the abundance of leaf chewers to be higher on gl1-1 and gl1-2 than on their different parental accessions, Ler-1 and Col-0, respectively. Furthermore, the loss of function in the GLABRA1 gene significantly decreased the resistance of plants to the two predominant chewers; flea beetles and turnip sawflies.

CONCLUSIONS

Overall, our results indicate that insect community composition significantly varies among A. thaliana accessions across two distant field sites, with GLABRA1 playing a key role in altering the abundance of leaf-chewing herbivores. Given that such a trichome variation is widely observed in Brassicaceae plants, the present study exemplifies the community-wide effect of a single plant gene on crucifer-feeding insects in the field.

Evolutionary and ecological functional genomics, from lab to the wild

Plant phenotypes are the result of both genetic and environmental forces that act to modulate trait expression. Over the last few years, numerous approaches in functional genomics and systems biology have led to a greater understanding of plant phenotypic variation and plant responses to the environment. These approaches, and the questions that they can address, have been loosely termed evolutionary and ecological functional genomics (EEFG), and have been providing key insights on how plants adapt and evolve. In particular, by bringing these studies from the laboratory to the field, EEFG studies allow us to gain greater knowledge of how plants function in their natural contexts.

Morphological plasticity and adaptation level of distylous Primula nivalis in a heterogeneous alpine environment

Plant populations at high elevation face extreme climatic conditions and resource limitations. The existence of distylous species at different elevations can help us investigate their adaptation to high altitudes, the evolution of their morphological characteristics, as well as their responses to limited resources. Here, 17 populations of Primula nivalis at different elevations were evaluated regarding variations in plant morphological characteristics, biomass allocation, and morphological plasticity in a heterogeneous environment. Our results demonstrate that heterogeneous environments can affect plant morphological characteristics and resource allocation in each sexual morph of these plants. Moreover, environmental variations reduced morphological plasticity in the two plant morphs, and the plasticity of long style (LS) plants was greater than that of short style (SS) plants. There were significant negative correlations between morphological characteristics and elevation, rainfall, temperature, and sunshine, and these are the main variables that affect morphological characteristics and resource allocation of both morphs of P. nivalis plants in heterogeneous environments. The morphological characteristics of P. nivalis plants transplanted from high to lower elevations were not significantly different in either population. LS plants had greater morphological plasticity and adaptability in heterogeneous environments than SS plants. Elevational gradients and heterogeneous environments differentiated both morphs of P. nivalis plants with regards to morphology as well as adaptations. LS plants showed a higher level of adaptability than SS plants.

Molecular basis of flowering under natural long-day conditions in Arabidopsis

Plants sense light and temperature changes to regulate flowering time. Here, we show that expression of the Arabidopsis florigen gene, FLOWERING LOCUS T (FT), peaks in the morning during spring, a different pattern than we observe in the laboratory. Providing our laboratory growth conditions with a red/far-red light ratio similar to open-field conditions and daily temperature oscillation is sufficient to mimic the FT expression and flowering time in natural long days. Under the adjusted growth conditions, key light signalling components, such as phytochrome A and EARLY FLOWERING 3, play important roles in morning FT expression. These conditions stabilize CONSTANS protein, a major FT activator, in the morning, which is probably a critical mechanism for photoperiodic flowering in nature. Refining the parameters of our standard growth conditions to more precisely mimic plant responses in nature can provide a powerful method for improving our understanding of seasonal response.

Transcriptome response of the foundation plant Spartina alterniflora to the Deepwater Horizon oil spill

Despite the severe impacts of the Deepwater Horizon oil spill, the foundation plant species Spartina alterniflora proved resilient to heavy oiling, providing an opportunity to identify mechanisms of response to the anthropogenic stress of crude oil exposure. We assessed plants from oil-affected and unaffected populations using a custom DNA microarray to identify genomewide transcription patterns and gene expression networks that respond to crude oil exposure. In addition, we used T-DNA insertion lines of the model grass Brachypodium distachyon to assess the contribution of four novel candidate genes to crude oil response. Responses in S. alterniflora to hydrocarbon exposure across the transcriptome as well as xenobiotic specific response pathways had little overlap with those previously identified in the model plant Arabidopsis thaliana. Among T-DNA insertion lines of B. distachyon, we found additional support for two candidate genes, one (ATTPS21) involved in volatile production, and the other (SUVH5) involved in epigenetic regulation of gene expression, that may be important in the response to crude oil. The architecture of crude oil response in S. alterniflora is unique from that of the model species A. thaliana, suggesting that xenobiotic response may be highly variable across plant species. In addition, further investigations of regulatory networks may benefit from more information about epigenetic response pathways.

The Long-Term "In Natura" Study Sites of Arabidopsis halleri for Plant Transcription and Epigenetic Modification Analyses in Natural Environments

The majority of organismal phenomena show functional significance in the context of natural environments. However, we know little about how dynamic gene expression is controlled under natural complex conditions. One of the most attractive challenges in current biology is to understand organismal functions in natural environments. We established and have developed long-term "in natura" study sites of Arabidopsis halleri to evaluate precise control of gene expression in natural environments. At the sites, we monitored meteorological factors, recorded plant growth and phenology, and collected RNA and chromatin samples to investigate dynamics of transcription and epigenetic modifications. Here, we introduce the in natura study sites, especially with the emphasis on methodologies for setting up study sites in natural plant populations and collecting samples used in transcriptomics and epigenetics in natural environments. Although the methods introduced here need to be modified depending on situations of one's study systems, our case can be a model for planning new in natura studies.

Reference-guided de novo assembly approach improves genome reconstruction for related species

Background

The development of next-generation sequencing has made it possible to sequence whole genomes at a relatively low cost. However, de novo genome assemblies remain challenging due to short read length, missing data, repetitive regions, polymorphisms and sequencing errors. As more and more genomes are sequenced, reference-guided assembly approaches can be used to assist the assembly process. However, previous methods mostly focused on the assembly of other genotypes within the same species. We adapted and extended a reference-guided de novo assembly approach, which enables the usage of a related reference sequence to guide the genome assembly. In order to compare and evaluate de novo and our reference-guided de novo assembly approaches, we used a simulated data set of a repetitive and heterozygotic plant genome.

Results

The extended reference-guided de novo assembly approach almost always outperforms the corresponding de novo assembly program even when a reference of a different species is used. Similar improvements can be observed in high and low coverage situations. In addition, we show that a single evaluation metric, like the widely used N50 length, is not enough to properly rate assemblies as it not always points to the best assembly evaluated with other criteria. Therefore, we used the summed z-scores of 36 different statistics to evaluate the assemblies.

Conclusions

The combination of reference mapping and de novo assembly provides a powerful tool to improve genome reconstruction by integrating information of a related genome. Our extension of the reference-guided de novo assembly approach enables the application of this strategy not only within but also between related species. Finally, the evaluation of genome assemblies is often not straight forward, as the truth is not known. Thus one should always use a combination of evaluation metrics, which not only try to assess the continuity but also the accuracy of an assembly.

Electronic supplementary material

The online version of this article (10.1186/s12859-017-1911-6) contains supplementary material, which is available to authorized users.

Molecular phenology in plants: in natura systems biology for the comprehensive understanding of seasonal responses under natural environments

Phenology refers to the study of seasonal schedules of organisms. Molecular phenology is defined here as the study of the seasonal patterns of organisms captured by molecular biology techniques. The history of molecular phenology is reviewed briefly in relation to advances in the quantification technology of gene expression. High-resolution molecular phenology (HMP) data have enabled us to study phenology with an approach of in natura systems biology. I review recent analyses of FLOWERING LOCUS C (FLC), a temperature-responsive repressor of flowering, along the six steps in the typical flow of in natura systems biology. The extensive studies of the regulation of FLC have made this example a successful case in which a comprehensive understanding of gene functions has been progressing. The FLC-mediated long-term memory of past temperatures creates time lags with other seasonal signals, such as photoperiod and short-term temperature. Major signals that control flowering time have a phase lag between them under natural conditions, and hypothetical phase lag calendars are proposed as mechanisms of season detection in plants. Transcriptomic HMP brings a novel strategy to the study of molecular phenology, because it provides a comprehensive representation of plant functions. I discuss future perspectives of molecular phenology from the standpoints of molecular biology, evolutionary biology and ecology.

Evolution of Selfing: Recurrent Patterns in Molecular Adaptation

Selfing has evolved in animals, fungi, and plants, and since Darwin's pioneering study, it is considered one of the most frequent evolutionary trends in flowering plants. Generally, the evolution of selfing is characterized by a loss of self-incompatibility, the selfing syndrome, and changes in genome-wide polymorphism patterns. Recent interdisciplinary studies involving molecular functional experiments, genome-wide data, experimental evolution, and evolutionary ecology using Arabidopsis thaliana, Caenorhabditis elegans, and other species show that the evolution of selfing is not merely a degradation of outcrossing traits but a model for studying the recurrent patterns underlying adaptive molecular evolution. For example, in wild Arabidopsis relatives, self-compatibility evolved from mutations in the male specificity gene, S-LOCUS CYSTEINE-RICH PROTEIN/S-LOCUS PROTEIN 11 (SCR/SP11), rather than the female specificity gene, S-LOCUS RECEPTOR KINASE (SRK), supporting the theoretical prediction of sexual asymmetry. Prevalence of dominant self-compatible mutations is consistent with Haldane's sieve, which acts against recessive adaptive mutations. Time estimates based on genome-wide polymorphisms and self-incompatibility genes generally support the recent origin of selfing.

Multiple independent origins of auto-pollination in tropical orchids (Bulbophyllum) in light of the hypothesis of selfing as an evolutionary dead end

BACKGROUND

The transition from outcrossing to selfing has long been portrayed as an 'evolutionary dead end' because, first, reversals are unlikely and, second, selfing lineages suffer from higher rates of extinction owing to a reduced potential for adaptation and the accumulation of deleterious mutations. We tested these two predictions in a clade of Madagascan Bulbophyllum orchids (30 spp.), including eight species where auto-pollinating morphs (i.e., selfers, without a 'rostellum') co-exist with their pollinator-dependent conspecifics (i.e., outcrossers, possessing a rostellum). Specifically, we addressed this issue on the basis of a time-calibrated phylogeny by means of ancestral character reconstructions and within the state-dependent evolution framework of BiSSE (Binary State Speciation and Extinction), which allowed jointly estimating rates of transition, speciation, and extinction between outcrossing and selfing.

RESULTS

The eight species capable of selfing occurred in scattered positions across the phylogeny, with two likely originating in the Pliocene (ca. 4.4-3.1 Ma), one in the Early Pleistocene (ca. 2.4 Ma), and five since the mid-Pleistocene (ca. ≤ 1.3 Ma). We infer that this scattered phylogenetic distribution of selfing is best described by models including up to eight independent outcrossing-to-selfing transitions and very low rates of speciation (and either moderate or zero rates of extinction) associated with selfing.

CONCLUSIONS

The frequent and irreversible outcrossing-to-selfing transitions in Madagascan Bulbophyllum are clearly congruent with the first prediction of the dead end hypothesis. The inability of our study to conclusively reject or support the likewise predicted higher extinction rate in selfing lineages might be explained by a combination of methodological limitations (low statistical power of our BiSSE approach to reliably estimate extinction in small-sized trees) and evolutionary processes (insufficient time elapsed for selfers to go extinct). We suggest that, in these tropical orchids, a simple genetic basis of selfing (via loss of the 'rostellum') is needed to explain the strikingly recurrent transitions to selfing, perhaps reflecting rapid response to parallel and novel selective environments over Late Quaternary (≤ 1.3 Ma) time scales.

Genetic variation in the flowering and yield formation of timothy (Phleum pratense L.) accessions after different photoperiod and vernalization treatments

Timothy is a perennial forage grass grown commonly in Boreal regions. This study explored the effect of vernalization and photoperiod (PP) on flowering and growth characteristics and how this related to changes in expression of three flowering related genes in accessions from different geographic origin. Large variation was found in accessions in their vernalization and PP responses. In southern accessions vernalization response or requirement was not observed, the heading date remained unchanged, and plants flowered without vernalization. On the contrary, northern types had obligatory requirement for vernalization and long PP, but the tiller elongation did not require vernalization at 16-h PP. Longer vernalization or PP treatments reduced the genotypical differences in flowering. Moreover, the vernalization saturation progressed stepwise from main tiller to lateral tillers, and this process was more synchronized in southern accessions. The expression of PpVRN1 was associated with vernalization while PpVRN3 accumulated at long PP. A crucial role for PpVRN3 in the transition to flowering was supported as in southern accession the transcript accumulated in non-vernalized plants after transfer to 16-h PP, and the apices transformed to generative stage. Differences in vernalization requirements were associated with variation in expression levels of PpVRN1 and PpVRN3, with higher expression levels in southern type. Most divergent transcript accumulation of PpMADS10 was found under different vernalization conditions. These differences between accessions can be translated into agronomic traits, such as the tiller composition of canopy, which affects the forage yield. The southern types, with minimal vernalization response, have fast re-growth ability and rapidly decreasing nutritive value, whereas northern types grow slowly and have better quality. This information can be utilized in breeding for new cultivars for longer growing seasons at high latitudes.

Evolution of the Selfing Syndrome in Arabis alpina (Brassicaceae)

Introduction

The transition from cross-fertilisation (outcrossing) to self-fertilisation (selfing) frequently coincides with changes towards a floral morphology that optimises self-pollination, the selfing syndrome. Population genetic studies have reported the existence of both outcrossing and selfing populations in Arabis alpina (Brassicaceae), which is an emerging model species for studying the molecular basis of perenniality and local adaptation. It is unknown whether its selfing populations have evolved a selfing syndrome.

Methods

Using macro-photography, microscopy and automated cell counting, we compared floral syndromes (size, herkogamy, pollen and ovule numbers) between three outcrossing populations from the Apuan Alps and three selfing populations from the Western and Central Alps (Maritime Alps and Dolomites). In addition, we genotyped the plants for 12 microsatellite loci to confirm previous measures of diversity and inbreeding coefficients based on allozymes, and performed Bayesian clustering.

Results and Discussion

Plants from the three selfing populations had markedly smaller flowers, less herkogamy and lower pollen production than plants from the three outcrossing populations, whereas pistil length and ovule number have remained constant. Compared to allozymes, microsatellite variation was higher, but revealed similar patterns of low diversity and high Fis in selfing populations. Bayesian clustering revealed two clusters. The first cluster contained the three outcrossing populations from the Apuan Alps, the second contained the three selfing populations from the Maritime Alps and Dolomites.

Conclusion

We conclude that in comparison to three outcrossing populations, three populations with high selfing rates are characterised by a flower morphology that is closer to the selfing syndrome. The presence of outcrossing and selfing floral syndromes within a single species will facilitate unravelling the genetic basis of the selfing syndrome, and addressing which selective forces drive its evolution.

Female sterility associated with increased clonal propagation suggests a unique combination of androdioecy and asexual reproduction in populations of Cardamine amara (Brassicaceae)

BACKGROUND AND AIMS

The coexistence of hermaphrodites and female-sterile individuals, or androdioecy, has been documented in only a handful of plants and animals. This study reports its existence in the plant species Cardamine amara (Brassicaceae), in which female-sterile individuals have shorter pistils than seed-producing hermaphrodites.

METHODS

Morphological analysis, in situ manual pollination, microsatellite genotyping and differential gene expression analysis using Arabidopsis microarrays were used to delimit variation between female-sterile individuals and hermaphrodites.

KEY RESULTS

Female sterility in C. amara appears to be caused by disrupted ovule development. It was associated with a 2.4- to 2.9-fold increase in clonal propagation. This made the pollen number of female-sterile genets more than double that of hermaphrodite genets, which fulfils a condition of co-existence predicted by simple androdioecy theories. When female-sterile individuals were observed in wild androdioecious populations, their ramet frequencies ranged from 5 to 54 %; however, their genet frequencies ranged from 11 to 29 %, which is consistent with the theoretically predicted upper limit of 50 %.

CONCLUSIONS

The results suggest that a combination of sexual reproduction and increased asexual proliferation by female-sterile individuals probably explains the invasion and maintenance of female sterility in otherwise hermaphroditic populations. To our knowledge, this is the first report of the coexistence of female sterility and hermaphrodites in the Brassicaceae.

Ten years of transcriptomics in wild populations: what have we learned about their ecology and evolution?

Molecular ecology has moved beyond the use of a relatively small number of markers, often noncoding, and it is now possible to use whole-genome measures of gene expression with microarrays and RNAseq (i.e. transcriptomics) to capture molecular response to environmental challenges. While transcriptome studies are shedding light on the mechanistic basis of traits as complex as personality or physiological response to catastrophic events, these approaches are still challenging because of the required technical expertise, difficulties with analysis and cost. Still, we found that in the last 10 years, 575 studies used microarrays or RNAseq in ecology. These studies broadly address three questions that reflect the progression of the field: (i) How much variation in gene expression is there and how is it structured? (ii) How do environmental stimuli affect gene expression? (iii) How does gene expression affect phenotype? We discuss technical aspects of RNAseq and microarray technology, and a framework that leverages the advantages of both. Further, we highlight future directions of research, particularly related to moving beyond correlation and the development of additional annotation resources. Measuring gene expression across an array of taxa in ecological settings promises to enrich our understanding of ecology and genome function.

Gene duplication and genetic exchange drive the evolution of S-RNase-based self-incompatibility in Petunia

Self-incompatibility (SI) systems in flowering plants distinguish self- and non-self pollen to prevent inbreeding. While other SI systems rely on the self-recognition between specific male- and female-determinants, the Solanaceae family has a non-self recognition system resulting in the detoxification of female-determinants of S-ribonucleases (S-RNases), expressed in pistils, by multiple male-determinants of S-locus F-box proteins (SLFs), expressed in pollen. It is not known how many SLF components of this non-self recognition system there are in Solanaceae species, or how they evolved. We identified 16-20 SLFs in each S-haplotype in SI Petunia, from a total of 168 SLF sequences using large-scale next-generation sequencing and genomic polymerase chain reaction (PCR) techniques. We predicted the target S-RNases of SLFs by assuming that a particular S-allele must not have a conserved SLF that recognizes its own S-RNase, and validated these predictions by transformation experiments. A simple mathematical model confirmed that 16-20 SLF sequences would be adequate to recognize the vast majority of target S-RNases. We found evidence of gene conversion events, which we suggest are essential to the constitution of a non-self recognition system and also contribute to self-compatible mutations.

Mating system shifts and transposable element evolution in the plant genus Capsella

Background

Despite having predominately deleterious fitness effects, transposable elements (TEs) are major constituents of eukaryote genomes in general and of plant genomes in particular. Although the proportion of the genome made up of TEs varies at least four-fold across plants, the relative importance of the evolutionary forces shaping variation in TE abundance and distributions across taxa remains unclear. Under several theoretical models, mating system plays an important role in governing the evolutionary dynamics of TEs. Here, we use the recently sequenced Capsella rubella reference genome and short-read whole genome sequencing of multiple individuals to quantify abundance, genome distributions, and population frequencies of TEs in three recently diverged species of differing mating system, two self-compatible species (C. rubella and C. orientalis) and their self-incompatible outcrossing relative, C. grandiflora.

Results

We detect different dynamics of TE evolution in our two self-compatible species; C. rubella shows a small increase in transposon copy number, while C. orientalis shows a substantial decrease relative to C. grandiflora. The direction of this change in copy number is genome wide and consistent across transposon classes. For insertions near genes, however, we detect the highest abundances in C. grandiflora. Finally, we also find differences in the population frequency distributions across the three species.

Conclusion

Overall, our results suggest that the evolution of selfing may have different effects on TE evolution on a short and on a long timescale. Moreover, cross-species comparisons of transposon abundance are sensitive to reference genome bias, and efforts to control for this bias are key when making comparisons across species.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-602) contains supplementary material, which is available to authorized users.

The evolution of selfing from outcrossing ancestors in Brassicaceae: what have we learned from variation at the S-locus?

Evolutionary transitions between mating systems have occurred repetitively and independently in flowering plants. One of the most spectacular advances of the recent empirical literature in the field was the discovery of the underlying genetic machinery, which provides the opportunity to retrospectively document the scenario of the outcrossing to selfing transitions in a phylogenetic perspective. In this review, we explore the literature describing patterns of polymorphism and molecular evolution of the locus controlling self-incompatibility (S-locus) in selfing species of the Brassicaceae family in order to document the transition from outcrossing to selfing, a retrospective approach that we describe as the 'mating system genes approach'. The data point to strikingly contrasted scenarios of transition from outcrossing to selfing. We also perform original analyses of the fully sequenced genomes of four species showing self-compatibility, to compare the orthologous S-locus region with that of functional S-locus haplotypes. Phylogenetic analyses suggest that all species we investigated evolved independently towards loss of self-incompatibility, and in most cases almost intact sequences of either of the two S-locus genes suggest that these transitions occurred relatively recently. The S-locus region in Aethionema arabicum, representing the most basal lineage of Brassicaceae, showed unusual patterns so that our analysis could not determine whether self-incompatibility was lost secondarily, or evolved in the core Brassicaceae after the split with this basal lineage. Although the approach we detail can only be used when mating system genes have been identified in a clade, we suggest that its integration with phylogenetic and population genetic approaches should help determine the main routes of this predominant mating system shift in plants.

Genome-wide quantification of homeolog expression ratio revealed nonstochastic gene regulation in synthetic allopolyploid Arabidopsis

Genome duplication with hybridization, or allopolyploidization, occurs commonly in plants, and is considered to be a strong force for generating new species. However, genome-wide quantification of homeolog expression ratios was technically hindered because of the high homology between homeologous gene pairs. To quantify the homeolog expression ratio using RNA-seq obtained from polyploids, a new method named HomeoRoq was developed, in which the genomic origin of sequencing reads was estimated using mismatches between the read and each parental genome. To verify this method, we first assembled the two diploid parental genomes of Arabidopsis halleri subsp. gemmifera and Arabidopsis lyrata subsp. petraea (Arabidopsis petraea subsp. umbrosa), then generated a synthetic allotetraploid, mimicking the natural allopolyploid Arabidopsis kamchatica. The quantified ratios corresponded well to those obtained by Pyrosequencing. We found that the ratios of homeologs before and after cold stress treatment were highly correlated (r = 0.870). This highlights the presence of nonstochastic polyploid gene regulation despite previous research identifying stochastic variation in expression. Moreover, our new statistical test incorporating overdispersion identified 226 homeologs (1.11% of 20 369 expressed homeologs) with significant ratio changes, many of which were related to stress responses. HomeoRoq would contribute to the study of the genes responsible for polyploid-specific environmental responses.

Nucleotide diversity of vernalization and flowering-time-related genes in a germplasm collection of meadow fescue (Festuca pratensis Huds. syn. Lolium pratense (Huds.) Darbysh.)

In plant species, control of flowering time is an important factor for adaptation to local natural environments. The Vrn1,CO,FT1 and CK2α genes are key components in the flowering-specific signaling pathway of grass species. Meadow fescue is an agronomically important forage grass species, which is naturally distributed across Europe and Western Asia. In this study, meadow fescue flowering-time-related genes were resequenced to assess nucleotide diversity in European and Western Asian subpopulations. Identified sequence polymorphisms were then converted into PCR-based molecular genetic markers, and a meadow fescue germplasm collection was genotyped to investigate global allelic variation. Lower nucleotide diversities were observed for the Vrn1 and CO orthologs, while relatively higher values were observed for the FT1 and casein kinase II α-subunit (CK2α) orthologs. The nucleotide diversity for FT1 orthologs in the Western Asian subpopulation was significantly higher than those of the European subpopulation. Similarly, significant differences in nucleotide diversity for the remaining genes were observed between several combinations of subpopulation. The global allele distribution pattern was consistent with observed level of nucleotide diversity. These results suggested that the degree of purifying selection acting on the genes differs according to geographical location. As previously shown for model plant species, functional specificities of flowering-time-related genes may also vary according to environmental conditions.

Validation of SNP allele frequencies determined by pooled next-generation sequencing in natural populations of a non-model plant species

Sequencing of pooled samples (Pool-Seq) using next-generation sequencing technologies has become increasingly popular, because it represents a rapid and cost-effective method to determine allele frequencies for single nucleotide polymorphisms (SNPs) in population pools. Validation of allele frequencies determined by Pool-Seq has been attempted using an individual genotyping approach, but these studies tend to use samples from existing model organism databases or DNA stores, and do not validate a realistic setup for sampling natural populations. Here we used pyrosequencing to validate allele frequencies determined by Pool-Seq in three natural populations of Arabidopsis halleri (Brassicaceae). The allele frequency estimates of the pooled population samples (consisting of 20 individual plant DNA samples) were determined after mapping Illumina reads to (i) the publicly available, high-quality reference genome of a closely related species (Arabidopsis thaliana) and (ii) our own de novo draft genome assembly of A. halleri. We then pyrosequenced nine selected SNPs using the same individuals from each population, resulting in a total of 540 samples. Our results show a highly significant and accurate relationship between pooled and individually determined allele frequencies, irrespective of the reference genome used. Allele frequencies differed on average by less than 4%. There was no tendency that either the Pool-Seq or the individual-based approach resulted in higher or lower estimates of allele frequencies. Moreover, the rather high coverage in the mapping to the two reference genomes, ranging from 55 to 284x, had no significant effect on the accuracy of the Pool-Seq. A resampling analysis showed that only very low coverage values (below 10-20x) would substantially reduce the precision of the method. We therefore conclude that a pooled re-sequencing approach is well suited for analyses of genetic variation in natural populations.

Mass flowering of the tropical tree Shorea beccariana was preceded by expression changes in flowering and drought-responsive genes

Community-level mass flowering, known as general flowering, which occurs in South-East Asia at supra-annual irregular intervals, is considered a particularly spectacular phenomenon in tropical ecology. Recent studies have proposed several proximate factors inducing general flowering, such as drought and falls in minimum temperature. However, limited empirical data on the developmental and physiological processes have been available to test the significance of such factors. To overcome this limitation and test the hypotheses that general flowering is triggered by the proposed factors, we conducted an 'ecological transcriptome' study of a mass flowering species, Shorea beccariana, comparing meteorological data with genome-wide expression patterns obtained using next-generation sequencing. Among the 98 flowering-related genes identified, the homologs of a floral pathway integrator, SbFT, and a floral repressor, SbSVP, showed dramatic transcriptional changes before flowering, and their flowering functions were confirmed using transgenic Arabidopsis thaliana. Expression in drought-responsive and sucrose-induced genes also changed before flowering. All these expression changes occurred when the flowering-inducing level of drought was reached, as estimated using data from the preceding 10 years. These genome-wide expression data support the hypothesis that drought is a trigger for general flowering.

Effects of pollen availability and the mutation bias on the fixation of mutations disabling the male specificity of self-incompatibility

The evolution of self-compatibility (SC) by the loss of self-incompatibility (SI) is regarded as one of the most frequent transitions in flowering plants. SI systems are generally characterized by specific interactions between the male and female specificity genes encoded at the S-locus. Recent empirical studies have revealed that the evolution of SC is often driven by male SC-conferring mutations at the S-locus rather than by female mutations. In this study, using a forward simulation model, we compared the fixation probabilities of male vs. female SC-conferring mutations at the S-locus. We explicitly considered the effects of pollen availability in the population and bias in the occurrence of SC-conferring mutations on the male and female specificity genes. We found that male SC-conferring mutations were indeed more likely to be fixed than were female SC-conferring mutations in a wide range of parameters. This pattern was particularly strong when pollen availability was relatively high. Under such a condition, even if the occurrence of mutations was biased strongly towards the female specificity gene, male SC-conferring mutations were much more often fixed. Our study demonstrates that fixation probabilities of those two types of mutation vary strongly depending on ecological and genetic conditions, although both types result in the same evolutionary consequence-the loss of SI.

Taxonomy and systematics are key to biological information: Arabidopsis, Eutrema (Thellungiella), Noccaea and Schrenkiella (Brassicaceae) as examples

Taxonomy and systematics provide the names and evolutionary framework for any biological study. Without these names there is no access to a biological context of the evolutionary processes which gave rise to a given taxon: close relatives and sister species (hybridization), more distantly related taxa (ancestral states), for example. This is not only true for the single species a research project is focusing on, but also for its relatives, which might be selected for comparative approaches and future research. Nevertheless, taxonomical and systematic knowledge is rarely fully explored and considered across biological disciplines. One would expect the situation to be more developed with model organisms such as Noccaea, Arabidopsis, Schrenkiella and Eutrema (Thellungiella). However, we show the reverse. Using Arabidopsis halleri and Noccaea caerulescens, two model species among metal accumulating taxa, we summarize and reflect past taxonomy and systematics of Arabidopsis and Noccaea and provide a modern synthesis of taxonomic, systematic and evolutionary perspectives. The same is presented for several species of Eutrema s. l. and Schrenkiella recently appeared as models for studying stress tolerance in plants and widely known under the name Thellungiella.

Recent loss of self-incompatibility by degradation of the male component in allotetraploid Arabidopsis kamchatica

The evolutionary transition from outcrossing to self-fertilization (selfing) through the loss of self-incompatibility (SI) is one of the most prevalent events in flowering plants, and its genetic basis has been a major focus in evolutionary biology. In the Brassicaceae, the SI system consists of male and female specificity genes at the S-locus and of genes involved in the female downstream signaling pathway. During recent decades, much attention has been paid in particular to clarifying the genes responsible for the loss of SI. Here, we investigated the pattern of polymorphism and functionality of the female specificity gene, the S-locus receptor kinase (SRK), in allotetraploid Arabidopsis kamchatica. While its parental species, A. lyrata and A. halleri, are reported to be diploid and mainly self-incompatible, A. kamchatica is self-compatible. We identified five highly diverged SRK haplogroups, found their disomic inheritance and, for the first time in a wild allotetraploid species, surveyed the geographic distribution of SRK at the two homeologous S-loci across the species range. We found intact full-length SRK sequences in many accessions. Through interspecific crosses with the self-incompatible and diploid congener A. halleri, we found that the female components of the SI system, including SRK and the female downstream signaling pathway, are still functional in these accessions. Given the tight linkage and very rare recombination of the male and female components on the S-locus, this result suggests that the degradation of male components was responsible for the loss of SI in A. kamchatica. Recent extensive studies in multiple Brassicaceae species demonstrate that the loss of SI is often derived from mutations in the male component in wild populations, in contrast to cultivated populations. This is consistent with theoretical predictions that mutations disabling male specificity are expected to be more strongly selected than mutations disabling female specificity, or the female downstream signaling pathway.

The predominantly selfing plant Arabidopsis thaliana experienced a recent reduction in transposable element abundance compared to its outcrossing relative Arabidopsis lyrata

Background

Transposable elements (TEs) are major contributors to genome evolution. One factor that influences their evolutionary dynamics is whether their host reproduces through selfing or through outcrossing. According to the recombinational spreading hypothesis, for instance, TEs can spread more easily in outcrossing species through recombination, and should thus be less abundant in selfing species. We here studied the distribution and evolutionary dynamics of TE families in the predominantly selfing plant Arabidopsis thaliana and its close outcrossing relative Arabidopsis lyrata on a genome-wide scale. We characterized differences in TE abundance between them and asked which, if any, existing hypotheses about TE abundances may explain these differences.

Results

We identified 1,819 TE families representing all known classes of TEs in both species, and found three times more copies in the outcrossing A. lyrata than in the predominantly selfing A. thaliana, as well as ten times more TE families unique to A. lyrata. On average, elements in A. lyrata are younger than elements in A. thaliana. In particular, A. thaliana shows a marked decrease in element number that occurred during the most recent 10% of the time interval since A. thaliana split from A. lyrata. This most recent period in the evolution of A. thaliana started approximately 500,000 years ago, assuming a splitting time of 5 million years ago, and coincides with the time at which predominant selfing originated.

Conclusions

Our results indicate that the mating system may be important for determining TE copy number, and that selfing species are likely to have fewer TEs.

New perspectives on the evolution of plant mating systems

BACKGROUND

The remarkable diversity of mating patterns and sexual systems in flowering plants has fascinated evolutionary biologists for more than a century. Enduring questions about this topic include why sexual polymorphisms have evolved independently in over 100 plant families, and why proportions of self- and cross-fertilization often vary dramatically within and among populations. Important new insights concerning the evolutionary dynamics of plant mating systems have built upon a strong foundation of theoretical models and innovative field and laboratory experiments. However, as the pace of advancement in this field has accelerated, it has become increasingly difficult for researchers to follow developments outside their primary area of research expertise.

SCOPE

In this Viewpoint paper we highlight three important themes that span and integrate different subdisciplines: the changes in morphology, phenology, and physiology that accompany the transition to selfing; the evolutionary consequences of pollen pool diversity in flowering plants; and the evolutionary dynamics of sexual polymorphisms. We also highlight recent developments in molecular techniques that will facilitate more efficient and cost-effective study of mating patterns in large natural populations, research on the dynamics of pollen transport, and investigations on the genetic basis of sexual polymorphisms. This Viewpoint also serves as the introduction to a Special Issue on the Evolution of Plant Mating Systems. The 15 papers in this special issue provide inspiring examples of recent discoveries, and glimpses of exciting developments yet to come.

Evolutionary systems biology: historical and philosophical perspectives on an emerging synthesis

Systems biology (SB) is at least a decade old now and maturing rapidly. A more recent field, evolutionary systems biology (ESB), is in the process of further developing system-level approaches through the expansion of their explanatory and potentially predictive scope. This chapter will outline the varieties of ESB existing today by tracing the diverse roots and fusions that make up this integrative project. My approach is philosophical and historical. As well as examining the recent origins of ESB, I will reflect on its central features and the different clusters of research it comprises. In its broadest interpretation, ESB consists of five overlapping approaches: comparative and correlational ESB; network architecture ESB; network property ESB; population genetics ESB; and finally, standard evolutionary questions answered with SB methods. After outlining each approach with examples, I will examine some strong general claims about ESB, particularly that it can be viewed as the next step toward a fuller modern synthesis of evolutionary biology (EB), and that it is also the way forward for evolutionary and systems medicine. I will conclude with a discussion of whether the emerging field of ESB has the capacity to combine an even broader scope of research aims and efforts than it presently does.