Contrasting gene-level signatures of selection with reproductive fitness.

Selection leaves signatures in the DNA sequence of genes, with many test statistics devised to detect its action. While these statistics are frequently used to support hypotheses about the adaptive significance of particular genes, the effect these genes have on reproductive fitness is rarely quantified experimentally. Consequently, it is unclear how gene-level signatures of selection are associated with empirical estimates of gene effect on fitness. Eukaryotic data sets that permit this comparison are very limited. Using the model plant Arabidopsis thaliana, for which these resources are available, we calculated seven gene-level substitution and polymorphism-based statistics commonly used to infer selection (dN/dS, NI, DOS, Tajima's D, Fu and Li's D*, Fay and Wu's H, and Zeng's E) and, using knockout lines, compared these to gene-level estimates of effect on fitness. We found that consistent with expectations, essential genes were more likely to be classified as negatively selected. By contrast, using 379 Arabidopsis genes for which data was available, we found no evidence that genes predicted to be positively selected had a significantly different effect on fitness than genes evolving more neutrally. We discuss these results in the context of the analytic challenges posed by Arabidopsis, one of the only systems in which this study could be conducted, and advocate for examination in additional systems. These results are relevant to the evaluation of genome-wide studies across species where experimental fitness data is unavailable, as well as highlighting an increasing need for the latter.

Natural variation in Arabidopsis shoot branching plasticity in response to nitrate supply affects fitness.

The capacity of organisms to tune their development in response to environmental cues is pervasive in nature. This phenotypic plasticity is particularly striking in plants, enabled by their modular and continuous development. A good example is the activation of lateral shoot branches in Arabidopsis, which develop from axillary meristems at the base of leaves. The activity and elongation of lateral shoots depends on the integration of many signals both external (e.g. light, nutrient supply) and internal (e.g. the phytohormones auxin, strigolactone and cytokinin). Here, we characterise natural variation in plasticity of shoot branching in response to nitrate supply using two diverse panels of Arabidopsis lines. We find extensive variation in nitrate sensitivity across these lines, suggesting a genetic basis for variation in branching plasticity. High plasticity is associated with extreme branching phenotypes such that lines with the most branches on high nitrate have the fewest under nitrate deficient conditions. Conversely, low plasticity is associated with a constitutively moderate level of branching. Furthermore, variation in plasticity is associated with alternative life histories with the low plasticity lines flowering significantly earlier than high plasticity lines. In Arabidopsis, branching is highly correlated with fruit yield, and thus low plasticity lines produce more fruit than high plasticity lines under nitrate deficient conditions, whereas highly plastic lines produce more fruit under high nitrate conditions. Low and high plasticity, associated with early and late flowering respectively, can therefore be interpreted alternative escape vs mitigate strategies to low N environments. The genetic architecture of these traits appears to be highly complex, with only a small proportion of the estimated genetic variance detected in association mapping.

Transgressive segregation reveals mechanisms of Arabidopsis immunity to Brassica-infecting races of white rust (Albugo candida).

Arabidopsis thaliana accessions are universally resistant at the adult leaf stage to white rust (Albugo candida) races that infect the crop species Brassica juncea and Brassica oleracea We used transgressive segregation in recombinant inbred lines to test if this apparent species-wide (nonhost) resistance in A. thaliana is due to natural pyramiding of multiple Resistance (R) genes. We screened 593 inbred lines from an Arabidopsis multiparent advanced generation intercross (MAGIC) mapping population, derived from 19 resistant parental accessions, and identified two transgressive segregants that are susceptible to the pathogen. These were crossed to each MAGIC parent, and analysis of resulting F2 progeny followed by positional cloning showed that resistance to an isolate of A. candida race 2 (Ac2V) can be explained in each accession by at least one of four genes encoding nucleotide-binding, leucine-rich repeat (NLR) immune receptors. An additional gene was identified that confers resistance to an isolate of A. candida race 9 (AcBoT) that infects B. oleracea Thus, effector-triggered immunity conferred by distinct NLR-encoding genes in multiple A. thaliana accessions provides species-wide resistance to these crop pathogens.

Inflorescence photosynthetic contribution to fitness releases Arabidopsis thaliana plants from trade-off constraints on early flowering.

Leaves are thought to be the primary carbon source for reproduction in plants, so a positive relationship between vegetative size and reproductive output is expected, establishing a trade-off between time to reproduction and reproductive output. A common response to higher temperatures due to climate changes is the induction of earlier transition into reproduction. Thus, in annual plants, earlier transition into flowering can potentially constrain plant size and reduce seed production. However, trade-offs between early reproduction and fitness are not always observed, suggesting mechanisms to escape the constraints of early flowering do exist. Here, we test whether inflorescence photosynthesis contribution to the reproductive output of Arabidopsis thaliana can offset the cost of early reproduction. We followed the development, growth rate and fitness of 15 accessions, and removed all rosette leaves at flowering (prior to the completion of inflorescence development or any fruit production) in half of the plants to determine the ability of inflorescences to maintain fitness in the absence of leaves. Although leaf removal significantly reduced fruit number, seed weight and plant height, even the most severely impacted accessions maintained 35% of their fitness with the inflorescence as the sole photosynthetic organ; and some accessions experienced no reduction in fitness. Differences between accessions in their ability to maintain fitness after leaf removal is best explained by earlier flowering time and the ability to maintain as many or more branches after leaf removal as in the control treatment. Although earlier flowering does constrain plant vegetative size, we found that inflorescence photosynthesis can significantly contribute to seed production, explaining why early flowering plants can maintain high fitness despite a reduction in vegetative size. Thus, plants can be released from the usually assumed trade-offs associated with earlier reproduction, and selection on inflorescence traits can mediate the impact of climate change on phenology.

Engaging with primary schools: Supporting the delivery of the new curriculum in evolution and inheritance.

The official school regulator in England (OFSTED) recently reported that the delivery of science lessons has been significantly diminished in many primary schools. There is concern that the lack of good quality science in school can reduce the recruitment of young scientists, and the level of science literacy among the general public. We believe university scientists and undergraduate students can have a significant impact in the delivery of science in primary schools. However, a relatively small proportion of scientists engage with young children to improve curricular primary school science education. Here, we argue that long term engagement with primary schools can produce significant impact for the scientist's research, schools, and society. As an example, we describe our experience developing teaching materials for the topic of "Evolution and inheritance"; highlighting possible pitfalls and perceived benefits, in hope of encouraging and facilitating other scientists to engage with primary schools.

Genomic Rearrangements in Arabidopsis Considered as Quantitative Traits.

To understand the population genetics of structural variants and their effects on phenotypes, we developed an approach to mapping structural variants that segregate in a population sequenced at low coverage. We avoid calling structural variants directly. Instead, the evidence for a potential structural variant at a locus is indicated by variation in the counts of short-reads that map anomalously to that locus. These structural variant traits are treated as quantitative traits and mapped genetically, analogously to a gene expression study. Association between a structural variant trait at one locus, and genotypes at a distant locus indicate the origin and target of a transposition. Using ultra-low-coverage (0.3×) population sequence data from 488 recombinant inbred Arabidopsis thaliana genomes, we identified 6502 segregating structural variants. Remarkably, 25% of these were transpositions. While many structural variants cannot be delineated precisely, we validated 83% of 44 predicted transposition breakpoints by polymerase chain reaction. We show that specific structural variants may be causative for quantitative trait loci for germination and resistance to infection by the fungus Albugo laibachii, isolate Nc14. Further we show that the phenotypic heritability attributable to read-mapping anomalies differs from, and, in the case of time to germination and bolting, exceeds that due to standard genetic variation. Genes within structural variants are also more likely to be silenced or dysregulated. This approach complements the prevalent strategy of structural variant discovery in fewer individuals sequenced at high coverage. It is generally applicable to large populations sequenced at low-coverage, and is particularly suited to mapping transpositions.

Lineage-specific sequence evolution and exon edge conservation partially explain the relationship between evolutionary rate and expression level in A. thaliana.

Rapidly evolving proteins can aid the identification of genes underlying phenotypic adaptation across taxa, but functional and structural elements of genes can also affect evolutionary rates. In plants, the 'edges' of exons, flanking intron junctions, are known to contain splice enhancers and to have a higher degree of conservation compared to the remainder of the coding region. However, the extent to which these regions may be masking indicators of positive selection or account for the relationship between dN/dS and other genomic parameters is unclear. We investigate the effects of exon edge conservation on the relationship of dN/dS to various sequence characteristics and gene expression parameters in the model plant Arabidopsis thaliana. We also obtain lineage-specific dN/dS estimates, making use of the recently sequenced genome of Thellungiella parvula, the second closest sequenced relative after the sister species Arabidopsis lyrata. Overall, we find that the effect of exon edge conservation, as well as the use of lineage-specific substitution estimates, upon dN/dS ratios partly explains the relationship between the rates of protein evolution and expression level. Furthermore, the removal of exon edges shifts dN/dS estimates upwards, increasing the proportion of genes potentially under adaptive selection. We conclude that lineage-specific substitutions and exon edge conservation have an important effect on dN/dS ratios and should be considered when assessing their relationship with other genomic parameters.

The genetic basis of natural variation in seed size and seed number and their trade-off using Arabidopsis thaliana MAGIC lines.

Offspring number and size are key traits determining an individual's fitness and a crop's yield. Yet, extensive natural variation within species is observed for these traits. Such variation is typically explained by trade-offs between fecundity and quality, for which an optimal solution is environmentally dependent. Understanding the genetic basis of seed size and number, as well as any possible genetic constraints preventing the maximization of both, is crucial from both an evolutionary and applied perspective. We investigated the genetic basis of natural variation in seed size and number using a set of Arabidopsis thaliana multiparent advanced generation intercross (MAGIC) lines. We also tested whether life history affects seed size, number, and their trade-off. We found that both seed size and seed number are affected by a large number of mostly nonoverlapping QTL, suggesting that seed size and seed number can evolve independently. The allele that increases seed size at most identified QTL is from the same natural accession, indicating past occurrence of directional selection for seed size. Although a significant trade-off between seed size and number is observed, its expression depends on life-history characteristics, and generally explains little variance. We conclude that the trade-off between seed size and number might have a minor role in explaining the maintenance of variation in seed size and number, and that seed size could be a valid target for selection.

Presence-absence variation in A. thaliana is primarily associated with genomic signatures consistent with relaxed selective constraints.

The sequencing of multiple genomes of the same plant species has revealed polymorphic gene and exon loss. Genes associated with disease resistance are overrepresented among those showing structural variations, suggesting an adaptive role for gene and exon presence-absence variation (PAV). To shed light on the possible functional relevance of polymorphic coding region loss and the mechanisms driving this process, we characterized genes that have lost entire exons or their whole coding regions in 17 fully sequenced Arabidopsis thaliana accessions. We found that although a significant enrichment in genes associated with certain functional categories is observed, PAV events are largely restricted to genes with signatures of reduced essentiality: PAV genes tend to be newer additions to the genome, tissue specific, and lowly expressed. In addition, PAV genes are located in regions of lower gene density and higher transposable element density. Partial coding region PAV events were associated with only a marginal reduction in gene expression level in the affected accession and occurred in genes with higher levels of alternative splicing in the Col-0 accession. Together, these results suggest that although adaptive scenarios cannot be ruled out, PAV events can be explained without invoking them.

Plant responses to elevated temperatures: a field study on phenological sensitivity and fitness responses to simulated climate warming.

Significant changes in plant phenology have been observed in response to increases in mean global temperatures. There are concerns that accelerated phenologies can negatively impact plant populations. However, the fitness consequence of changes in phenology in response to elevated temperature is not well understood, particularly under field conditions. We address this issue by exposing a set of recombinant inbred lines of Arabidopsis thaliana to a simulated global warming treatment in the field. We find that plants exposed to elevated temperatures flower earlier, as predicted by photothermal models. However, contrary to life-history trade-off expectations, they also flower at a larger vegetative size, suggesting that warming probably causes acceleration in vegetative development. Although warming increases mean fitness (fruit production) by ca. 25%, there is a significant genotype-by-environment interaction. Changes in fitness rank indicate that imminent climate change can cause populations to be maladapted in their new environment, if adaptive evolution is limited. Thus, changes in the genetic composition of populations are likely, depending on the species' generation time and the speed of temperature change. Interestingly, genotypes that show stronger phenological responses have higher fitness under elevated temperatures, suggesting that phenological sensitivity might be a good indicator of success under elevated temperature at the genotypic level as well as at the species level.

Climate envelope modelling reveals intraspecific relationships among flowering phenology, niche breadth and potential range size in Arabidopsis thaliana.

Species often harbour large amounts of phenotypic variation in ecologically important traits, and some of this variation is genetically based. Understanding how this genetic variation is spatially structured can help to understand species' ecological tolerances and range limits. We modelled the climate envelopes of Arabidopsis thaliana genotypes, ranging from early- to late-flowering, as a function of several climatic variables. We found that genotypes with contrasting alleles at individual flowering time loci differed significantly in potential range size and niche breadth. We also found that later flowering genotypes had more restricted range potentials and narrower niche breadths than earlier flowering genotypes, indicating that local selection on flowering can constrain or enhance the ability of populations to colonise other areas. Our study demonstrates how climate envelope models that incorporate ecologically important genetic variation can provide insights into the macroecology of a species, which is important to understand its responses to changing environments.

Natural variation in GA1 associates with floral morphology in Arabidopsis thaliana.

• The genetic architecture of floral traits is evolutionarily important due to the fitness consequences of quantitative variation in floral morphology. Yet, little is known about the genes underlying these traits in natural populations. Using Arabidopsis thaliana, we examine molecular variation at GIBBERELLIC ACID REQUIRING 1 (GA1) and test for associations with floral morphology. • We examined full-length sequence in 32 accessions and describe two haplotypes (comprising four nonsynonymous polymorphisms) in GA1 that segregate at intermediate frequencies. In 133 A. thaliana accessions, we test for genotype-phenotype associations and corroborate these findings in segregating progenies. • The two common GA1 haplotypes were associated with the length of petals, stamens, and to a lesser extent style-stigma length. Associations were confirmed in a segregating progeny developed from 19 accessions. We find analogous results in recombinant inbred lines of the Bayreuth × Shahdara cross, which differ only at one of 4 SNPs, suggesting that this SNP may contribute to the observed association. • Assuming GA1 causally affects floral organ size, it is interesting that adjacent petal and stamen whorls are most strongly affected. This pattern suggests that GA1 could contribute to the greater strength of petal-stamen correlations relative to other floral-length correlations observed in some Brassicaceous species.

Mapping the genetic basis of ecologically and evolutionarily relevant traits in Arabidopsis thaliana.

There has been a long standing interest in the relationship between genetic and phenotypic variation in natural populations, in order to understand the genetic basis of adaptation and to discover natural alleles to improve crops. Here we review recent developments in mapping approaches that have significantly improved our ability to identify causal polymorphism explaining natural variation in ecological and evolutionarily relevant traits. However, challenges in interpreting these discoveries remain. In particular, we need more detailed transcriptomic, epigenomic, and gene network data to help understand the mechanisms behind identified associations. Also, more studies need to be performed under field conditions or using experimental evolution to determine whether polymorphisms identified in the lab are relevant for adaptation and improvement under natural conditions.

Multiple reference genomes and transcriptomes for Arabidopsis thaliana.

Genetic differences between Arabidopsis thaliana accessions underlie the plant's extensive phenotypic variation, and until now these have been interpreted largely in the context of the annotated reference accession Col-0. Here we report the sequencing, assembly and annotation of the genomes of 18 natural A. thaliana accessions, and their transcriptomes. When assessed on the basis of the reference annotation, one-third of protein-coding genes are predicted to be disrupted in at least one accession. However, re-annotation of each genome revealed that alternative gene models often restore coding potential. Gene expression in seedlings differed for nearly half of expressed genes and was frequently associated with cis variants within 5 kilobases, as were intron retention alternative splicing events. Sequence and expression variation is most pronounced in genes that respond to the biotic environment. Our data further promote evolutionary and functional studies in A. thaliana, especially the MAGIC genetic reference population descended from these accessions.

Functional genetics of intraspecific ecological interactions in Arabidopsis thaliana.

Studying the genetic basis of traits involved in ecological interactions is a fundamental part of elucidating the connections between evolutionary and ecological processes. Such knowledge allows one to link genetic models of trait evolution with ecological models describing interactions within and between species. Previous work has shown that connections between genetic and ecological processes in Arabidopsis thaliana may be mediated by the fact that quantitative trait loci (QTL) with 'direct' effects on traits of individuals also have pleiotropic 'indirect' effects on traits expressed in neighbouring plants. Here, we further explore these connections by examining functional relationships between traits affected directly and indirectly by the same QTL. We develop a novel approach using structural equation models (SEMs) to determine whether observed pleiotropic effects result from traits directly affected by the QTL in focal individuals causing the changes in the neighbours' phenotypes. This hypothesis was assessed using SEMs to test whether focal plant phenotypes appear to mediate the connection between the focal plants' genotypes and the phenotypes of their neighbours, or alternatively, whether the connection between the focal plants' genotypes and the neighbours' phenotypes is mediated by unmeasured traits. We implement this analysis using a QTL of major effect that maps to the well-characterized flowering locus, FRIGIDA. The SEMs support the hypothesis that the pleiotropic indirect effects of this locus arise from size and developmental timing-related traits in focal plants affecting the expression of developmental traits in their neighbours. Our findings provide empirical insights into the genetics and nature of intraspecific ecological interactions. Our technique holds promise in directing future work into the genetic basis and functional relationship of traits mediating and responding to ecological interactions.

Paternal effects in Arabidopsis indicate that offspring can influence their own size.

The existence of genetic variation in offspring size in plants and animals is puzzling because offspring size is often strongly associated with fitness and expected to be under stabilizing selection. An explanation for variation in seed size is conflict between parents and between parents and offspring. However, for this hypothesis to be true, it must be shown that the offspring genotype can affect its own size. The existence of paternal effects would support this hypothesis, but these have rarely been shown. Using a diallel cross among four natural accessions of Arabidopsis thaliana we show that maternal, paternal and positional effects jointly influence seed size, number and the frequency of seed abortion. We found that seed abortion (%) depends on the combination of maternal and paternal genotypes, suggesting the existence of mate choice or epistatic incompatibility among accessions of A. thaliana. In addition, since paternal genotype explains approximately 10 per cent of the variation in seed size, we propose that A. thaliana's offspring must influence the amount of resources allocated to themselves. Identification of paternal effects in Arabidopsis should facilitate dissection of the genetic mechanisms involved in paternal effects.

A Multiparent Advanced Generation Inter-Cross to fine-map quantitative traits in Arabidopsis thaliana.

Identifying natural allelic variation that underlies quantitative trait variation remains a fundamental problem in genetics. Most studies have employed either simple synthetic populations with restricted allelic variation or performed association mapping on a sample of naturally occurring haplotypes. Both of these approaches have some limitations, therefore alternative resources for the genetic dissection of complex traits continue to be sought. Here we describe one such alternative, the Multiparent Advanced Generation Inter-Cross (MAGIC). This approach is expected to improve the precision with which QTL can be mapped, improving the outlook for QTL cloning. Here, we present the first panel of MAGIC lines developed: a set of 527 recombinant inbred lines (RILs) descended from a heterogeneous stock of 19 intermated accessions of the plant Arabidopsis thaliana. These lines and the 19 founders were genotyped with 1,260 single nucleotide polymorphisms and phenotyped for development-related traits. Analytical methods were developed to fine-map quantitative trait loci (QTL) in the MAGIC lines by reconstructing the genome of each line as a mosaic of the founders. We show by simulation that QTL explaining 10% of the phenotypic variance will be detected in most situations with an average mapping error of about 300 kb, and that if the number of lines were doubled the mapping error would be under 200 kb. We also show how the power to detect a QTL and the mapping accuracy vary, depending on QTL location. We demonstrate the utility of this new mapping population by mapping several known QTL with high precision and by finding novel QTL for germination data and bolting time. Our results provide strong support for similar ongoing efforts to produce MAGIC lines in other organisms.

Candidate gene association mapping of Arabidopsis flowering time.

The pathways responsible for flowering time in Arabidopsis thaliana comprise one of the best characterized genetic networks in plants. We harness this extensive molecular genetic knowledge to identify potential flowering time quantitative trait genes (QTGs) through candidate gene association mapping using 51 flowering time loci. We genotyped common single nucleotide polymorphisms (SNPs) at these genes in 275 A. thaliana accessions that were also phenotyped for flowering time and rosette leaf number in long and short days. Using structured association techniques, we find that haplotype-tagging SNPs in 27 flowering time genes show significant associations in various trait/environment combinations. After correction for multiple testing, between 2 and 10 genes remain significantly associated with flowering time, with CO arguably possessing the most promising associations. We also genotyped a subset of these flowering time gene SNPs in an independent recombinant inbred line population derived from the intercrossing of 19 accessions. Approximately one-third of significant polymorphisms that were associated with flowering time in the accessions and genotyped in the outbred population were replicated in both mapping populations, including SNPs at the CO, FLC, VIN3, PHYD, and GA1 loci, and coding region deletions at the FRI gene. We conservatively estimate that approximately 4-14% of known flowering time genes may harbor common alleles that contribute to natural variation in this life history trait.

Standing genetic variation in FRIGIDA mediates experimental evolution of flowering time in Arabidopsis.

The role of standing genetic variation in adaptive evolution remains unclear. Although there has been much progress in identifying candidate genes that underlie adaptive traits, we still lack direct evidence that natural allelic variation in these genes can actually mediate adaptive evolution. In this study, we investigate the role of natural allelic variation in two candidate flowering time genes, in response to selection for early flowering in Arabidopsis thaliana: FRIGIDA (FRI) and FLOWERING LOCUS C (FLC). We performed artificial selection for early flowering under 'spring-' and 'winter-annual' growth conditions using an outbred population of A. thaliana produced by intermating 19 natural accessions. FRI and FLC are involved in A. thaliana's response to winter conditions, and nonfunctional and weak alleles at these loci are know to reduce flowering time, particularly under spring-annual conditions. Our results provide direct evidence that natural allelic variation in FRI can provide rapid and predictable adaptive evolution in flowering time under spring-annual conditions. We observed a strong response to selection, in terms of reducing flowering time, in both growth conditions (approximately 2 standard deviation reduction). Concomitantly, the frequency of functional FRI alleles under spring-annual conditions was reduced by 68%, in agreement with predicted changes. No significant changes in allele frequencies were observed in FRI in the winter-annual growth condition or in FLC for either growth conditions. These results indicate that changes in flowering time are mediated by different genetic factors under spring- and winter-annual growth conditions, and that other loci must also be contributing to the response to selection.

Antagonistic pleiotropic effects reduce the potential adaptive value of the FRIGIDA locus.

Although the occurrence of epistasis and pleiotropy is widely accepted at the molecular level, its effect on the adaptive value of fitness-related genes is rarely investigated in plants. Knowledge of these features of a gene is critical to understand the molecular basis of adaptive evolution. Here we investigate the importance of pleiotropy and epistasis in determining the adaptive value of a candidate gene using the gene FRI (FRIGIDA), which is thought to be the major gene controlling flowering time variation in Arabidopsis thaliana. The effect of FRI on flowering time was analyzed in an outbred population created by randomly mating 19 natural accessions of A. thaliana. This unique population allows the estimation of FRI effects independent of any linkage association with other loci due to demographic processes or to coadapted genes. It also allows for the estimation of pleiotropic effects of FRI on fitness and inflorescence architecture. We found that FRI explains less variation in flowering time than previously observed among natural accessions, and interacts epistatically with the FLC locus. Although early flowering plants produce more fruits under spring conditions, and nonfunctional alleles of FRI were associated with early flowering, variation at FRI was not associated with fitness. We show that nonfunctional FRI alleles have negative pleiotropic effects on fitness by reducing the numbers of nodes and branches on the inflorescence. We propose that these antagonistic pleiotropic effects reduce the adaptive value of FRI, and helps explain the maintenance of alternative life history strategies across natural populations of A. thaliana.