Intraspecific plant–soil feedback and intraspecific overyielding in Arabidopsis thaliana

Understanding the mechanisms of community coexistence and ecosystem functioning may help to counteract the current biodiversity loss and its potentially harmful consequences. In recent years, plant–soil feedback that can, for example, be caused by below‐ground microorganisms has been suggested to play a role in maintaining plant coexistence and to be a potential driver of the positive relationship between plant diversity and ecosystem functioning. Most of the studies addressing these topics have focused on the species level. However, in addition to interspecific interactions, intraspecific interactions might be important for the structure of natural communities. Here, we examine intraspecific coexistence and intraspecific diversity effects using 10 natural accessions of the model species Arabidopsis thaliana (L.) Heynh. We assessed morphological intraspecific diversity by measuring several above‐ and below‐ground traits. We performed a plant–soil feedback experiment that was based on these trait differences between the accessions in order to determine whether A. thaliana experiences feedback at intraspecific level as a result of trait differences. We also experimentally tested the diversity–productivity relationship at intraspecific level. We found strong differences in above‐ and below‐ground traits between the A. thaliana accessions. Overall, plant–soil feedback occurred at intraspecific level. However, accessions differed in the direction and strength of this feedback: Some accessions grew better on their own soils, some on soils from other accessions. Furthermore, we found positive diversity effects within A. thaliana: Accession mixtures produced a higher total above‐ground biomass than accession monocultures. Differences between accessions in their feedback response could not be explained by morphological traits. Therefore, we suggest that they might have been caused by accession‐specific accumulated soil communities, by root exudates, or by accession‐specific resource use based on genetic differences that are not expressed in morphological traits. Synthesis. Our results provide some of the first evidence for intraspecific plant–soil feedback and intraspecific overyielding. These findings may have wider implications for the maintenance of variation within species and the importance of this variation for ecosystem functioning. Our results highlight the need for an increased focus on intraspecific processes in plant diversity research to fully understand the mechanisms of coexistence and ecosystem functioning.     

Fisher's geometric model predicts the effects of random mutations when tested in the wild

Fisher's geometric model of adaptation (FGM) has been the conceptual foundation for studies investigating the genetic basis of adaptation since the onset of the neo Darwinian synthesis. FGM describes adaptation as the movement of a genotype toward a fitness optimum due to beneficial mutations. To date, one prediction of FGM, the probability of improvement is related to the distance from the optimum, has only been tested in microorganisms under laboratory conditions. There is reason to believe that results might differ under natural conditions where more mutations likely affect fitness, and where environmental variance may obscure the expected pattern. We chemically induced mutations into a set of 19 Arabidopsis thaliana accessions from across the native range of A. thaliana and planted them alongside the premutated founder lines in two habitats in the mid‐Atlantic region of the United States under field conditions. We show that FGM is able to predict the outcome of a set of random induced mutations on fitness in a set of A. thaliana accessions grown in the wild: mutations are more likely to be beneficial in relatively less fit genotypes. This finding suggests that FGM is an accurate approximation of the process of adaptation under more realistic ecological conditions.     

Population genetics structure of glyphosate‐resistant Johnsongrass (Sorghum halepense L. Pers) does not support a single origin of the resistance

Single sequence repeats (SSR) developed for Sorghum bicolor were used to characterize the genetic distance of 46 different Sorghum halepense (Johnsongrass) accessions from Argentina some of which have evolved toward glyphosate resistance. Since Johnsongrass is an allotetraploid and only one subgenome is homologous to cultivated sorghum, some SSR loci amplified up to two alleles while others (presumably more conserved loci) amplified up to four alleles. Twelve SSR providing information of 24 loci representative of Johnsongrass genome were selected for genetic distance characterization. All of them were highly polymorphic, which was evidenced by the number of different alleles found in the samples studied, in some of them up to 20. UPGMA and Mantel analysis showed that Johnsongrass glyphosate‐resistant accessions that belong to different geographic regions do not share similar genetic backgrounds. In contrast, they show closer similarity to their neighboring susceptible counterparts. Discriminant Analysis of Principal Components using the clusters identified by K‐means support the lack of a clear pattern of association among samples and resistance status or province of origin. Consequently, these results do not support a single genetic origin of glyphosate resistance. Nucleotide sequencing of the 5‐enolpyruvylshikimate‐3‐phosphate synthase (EPSPS) encoding gene from glyphosate‐resistant and susceptible accessions collected from different geographic origins showed that none presented expected mutations in aminoacid positions 101 and 106 which are diagnostic of target‐site resistance mechanism.     

A pair of receptor‐like kinases is responsible for natural variation in shoot growth response to mannitol treatment in Arabidopsis thaliana

Growth is a complex trait that adapts to the prevailing conditions by integrating many internal and external signals. Understanding the molecular origin of this variation remains a challenging issue. In this study, natural variation of shoot growth under mannitol‐induced stress was analyzed by standard quantitative trait locus mapping methods in a recombinant inbred line population derived from a cross between the Col‐0 and Cvi‐0 Arabidopsis thaliana accessions. Cloning of a major QTL specific to mannitol‐induced stress condition led to identification of EGM1 and EGM2, a pair of tandem‐duplicated genes encoding receptor‐like kinases that are potentially involved in signaling of mannitol‐associated stress responses. Using various genetic approaches, we identified two non‐synonymous mutations in the EGM2[Cvi] allele that are shared by at least ten accessions from various origins and are probably responsible for a specific tolerance to mannitol. We have shown that the enhanced shoot growth phenotype contributed by the Cvi allele is not linked to generic osmotic properties but instead to a specific chemical property of mannitol itself. This result raises the question of the function of such a gene in A. thaliana, a species that does not synthesize mannitol. Our findings suggest that the receptor‐like kinases encoded by EGM genes may be activated by mannitol produced by pathogens such as fungi, and may contribute to plant defense responses whenever mannitol is present.     

Rapid identification of an Arabidopsis NLR gene as a candidate conferring susceptibility to Sclerotinia sclerotiorum using time‐resolved automated phenotyping

The broad host range necrotrophic fungus Sclerotinia sclerotiorum is a devastating pathogen of many oil and vegetable crops. Plant genes conferring complete resistance against S. sclerotiorum have not been reported. Instead, plant populations challenged by S. sclerotiorum exhibit a continuum of partial resistance designated as quantitative disease resistance (QDR). Because of their complex interplay and their small phenotypic effect, the functional characterization of QDR genes remains limited. How broad host range necrotrophic fungi manipulate plant programmed cell death is for instance largely unknown. Here, we designed a time‐resolved automated disease phenotyping pipeline enabling high‐throughput disease lesion measurement with high resolution, low footprint at low cost. We could accurately recover contrasted disease responses in several pathosystems using this system. We used our phenotyping pipeline to assess the kinetics of disease symptoms caused by seven S. sclerotiorum isolates on six A. thaliana natural accessions with unprecedented resolution. Large effect polymorphisms common to the most resistant A. thaliana accessions identified highly divergent alleles of the nucleotide‐binding site leucine‐rich repeat gene LAZ5 in the resistant accessions Rubezhnoe and Lip‐0. We show that impaired LAZ5 expression in laz5.1 mutant lines and in A. thaliana Rub natural accession correlate with enhanced QDR to S. sclerotiorum. These findings illustrate the value of time‐resolved image‐based phenotyping for unravelling the genetic bases of complex traits such as QDR. Our results suggest that S. sclerotiorum manipulates plant sphingolipid pathways guarded by LAZ5 to trigger programmed cell death and cause disease.     

Natural variation in tocochromanols content in Arabidopsis thaliana accessions – the effect of temperature and light intensity

In this study, 25 accessions of Arabidopsis thaliana originating from a variety of climate conditions were grown under controlled circumstances of different light intensity and temperature. The accessions were analyzed for prenyllipids content and composition, as well as expression of the genes involved in tocochromanol biosynthesis (vte1‐5). It was found that the applied conditions did not strongly affect total tocochromanols content and there was no apparent correlation of the tocochromanol content with the origin of the accessions. However, the presented results indicate that the temperature, more than the light intensity, affects the expression of the vte1‐5 genes and the content of some prenyllipids. An interesting observation was that under low growth temperature, the hydroxy‐plastochromanol (PC‐OH) to plastochromanol (PC) ratio was considerably increased regardless of the light intensity in most of the accessions. PC‐OH is known to be formed as a result of singlet oxygen stress, therefore this observation indicates that the singlet oxygen production is enhanced under low temperature. Unexpectedly, the highest increase in the PC‐OH/PC ratio was found for accessions originating from cold climate (Shigu, Krazo‐1 and Lov‐5), even though such plants could be expected to be more resistant to low temperature stress.     

Patterns of univariate and multivariate plasticity to elevated carbon dioxide in six European populations of Arabidopsis thaliana

The impact of elevated carbon dioxide on plants is a growing concern in evolutionary ecology and global change biology. Characterizing patterns of phenotypic integration and multivariate plasticity to elevated carbon dioxide can provide insights into ecological and evolutionary dynamics in future human‐altered environments. Here, we examined univariate and multivariate responses to carbon enrichment in six functional traits among six European accessions of Arabidopsis thaliana. We detected phenotypic plasticity in both univariate and multivariate phenotypes, but did not find significant variation in plasticity (genotype by environment interactions) within or among accessions. Eigenvector, eigenvalue variance, and common principal components analyses showed that elevated carbon dioxide altered patterns of trait covariance, reduced the strength of phenotypic integration, and decreased population‐level differentiation in the multivariate phenotype. Our data suggest that future carbon dioxide conditions may influence evolutionary dynamics in natural populations of A. thaliana.     

Length polymorphism and allele structure of trinucleotide microsatellites in natural accessions of Arabidopsis thaliana

 The objective of this work was to assess the degree of trinucleotide microsatellite length polymorphism in the selfing species Arabidopsis thaliana. PCR amplifications of 12 microsatellite loci among 49 natural populations revealed between one to eight length variants (alleles) for each locus. The average number of alleles per locus was four and the average genetic diversity index was 0.43. Divergence between length variants was investigated at the nucleotide level. Several observations emerge from the sequence data: (1) for most loci, length polymorphism results only from variations in the number of trinucleotide repeats; (2) for a few others, some variability was noted in the flanking sequences; (3) for compound and interrupted loci containing two arrays of trinucleotide repeats, length variations preferentially affect the longest one. Five of the Arabidopsis thaliana accessions were clearly composed of two sublines. In 2 other accessions, some heterozygous individual plants, probably resulting from recent outcrosses, were found. A phylogenetic tree constructed on the basis of trinucleotide microsatellite allelic diversity shows that genetic relationships among the accessions are not correlated with their geographic origin. Received: 4 November 1997 / Accepted: 3 March 1998     

The Pattern of Polymorphism in Arabidopsis thaliana

We resequenced 876 short fragments in a sample of 96 individuals of Arabidopsis thaliana that included stock center accessions as well as a hierarchical sample from natural populations. Although A. thaliana is a selfing weed, the pattern of polymorphism in general agrees with what is expected for a widely distributed, sexually reproducing species. Linkage disequilibrium decays rapidly, within 50 kb. Variation is shared worldwide, although population structure and isolation by distance are evident. The data fail to fit standard neutral models in several ways. There is a genome-wide excess of rare alleles, at least partially due to selection. There is too much variation between genomic regions in the level of polymorphism. The local level of polymorphism is negatively correlated with gene density and positively correlated with segmental duplications. Because the data do not fit theoretical null distributions, attempts to infer natural selection from polymorphism data will require genome-wide surveys of polymorphism in order to identify anomalous regions. Despite this, our data support the utility of A. thaliana as a model for evolutionary functional genomics.     

The pattern of polymorphism in Arabidopsis thaliana

We resequenced 876 short fragments in a sample of 96 individuals of Arabidopsis thaliana that included stock center accessions as well as a hierarchical sample from natural populations. Although A. thaliana is a selfing weed, the pattern of polymorphism in general agrees with what is expected for a widely distributed, sexually reproducing species. Linkage disequilibrium decays rapidly, within 50 kb. Variation is shared worldwide, although population structure and isolation by distance are evident. The data fail to fit standard neutral models in several ways. There is a genome-wide excess of rare alleles, at least partially due to selection. There is too much variation between genomic regions in the level of polymorphism. The local level of polymorphism is negatively correlated with gene density and positively correlated with segmental duplications. Because the data do not fit theoretical null distributions, attempts to infer natural selection from polymorphism data will require genome-wide surveys of polymorphism in order to identify anomalous regions. Despite this, our data support the utility of A. thaliana as a model for evolutionary functional genomics.     

Characterizing genomic variation of Arabidopsis thaliana: the roles of geography and climate

Arabidopsis thaliana inhabits diverse climates and exhibits varied phenology across its range. Although A. thaliana is an extremely well‐studied model species, the relationship between geography, growing season climate and its genetic variation is poorly characterized. We used redundancy analysis (RDA) to quantify the association of genomic variation [214 051 single nucleotide polymorphisms (SNPs)] with geography and climate among 1003 accessions collected from 447 locations in Eurasia. We identified climate variables most correlated with genomic variation, which may be important selective gradients related to local adaptation across the species range. Climate variation among sites of origin explained slightly more genomic variation than geographical distance. Large‐scale spatial gradients and early spring temperatures explained the most genomic variation, while growing season and summer conditions explained the most after controlling for spatial structure. SNP variation in Scandinavia showed the greatest climate structure among regions, possibly because of relatively consistent phenology and life history of populations in this region. Climate variation explained more variation among nonsynonymous SNPs than expected by chance, suggesting that much of the climatic structure of SNP correlations is due to changes in coding sequence that may underlie local adaptation.     

The Arabidopsis STRESS RESPONSE SUPPRESSOR DEAD‐box RNA helicases are nucleolar‐ and chromocenter‐localized proteins that undergo stress‐mediated relocalization and are involved in epigenetic gene silencing

DEAD‐box RNA helicases are involved in many aspects of RNA metabolism and in diverse biological processes in plants. Arabidopsis thaliana mutants of two DEAD‐box RNA helicases, STRESS RESPONSE SUPPRESSOR1 (STRS1) and STRS2 were previously shown to exhibit tolerance to abiotic stresses and up‐regulated stress‐responsive gene expression. Here, we show that Arabidopsis STRS‐overexpressing lines displayed a less tolerant phenotype and reduced expression of stress‐induced genes confirming the STRSs as attenuators of Arabidopsis stress responses. GFP–STRS fusion proteins exhibited localization to the nucleolus, nucleoplasm and chromocenters and exhibited relocalization in response to abscisic acid (ABA) treatment and various stresses. This relocalization was reversed when stress treatments were removed. The STRS proteins displayed mis‐localization in specific gene‐silencing mutants and exhibited RNA‐dependent ATPase and RNA‐unwinding activities. In particular, STRS2 showed mis‐localization in three out of four mutants of the RNA‐directed DNA methylation (RdDM) pathway while STRS1 was mis‐localized in the hd2c mutant that is defective in histone deacetylase activity. Furthermore, heterochromatic RdDM target loci displayed reduced DNA methylation and increased expression in the strs mutants. Taken together, our findings suggest that the STRS proteins are involved in epigenetic silencing of gene expression to bring about suppression of the Arabidopsis stress response.     

Functional variation in a disease resistance gene in populations of Arabidopsis thaliana

Analyses of functional genetic diversity in natural populations may provide important new insights into gene function and are necessary to understand the evolutionary processes maintaining diversity itself. The importance of including diversity within and between local populations in such studies is often ignored although many of the processes affecting genetic diversity act on this scale. Here we examine the molecular diversity in RPW8 (Recognition of Powdery Mildew), a gene conferring broad‐spectrum resistance to powdery mildews in Arabidopsis thaliana stock‐center accessions. Our eight UK study populations of the weedy A. thaliana were from locations judged to be subject to a minimum of anthropogenic disturbance and potentially long established. The majority of populations comprised considerable variation both in disease phenotype and RPW8 genotype. Although resistant individuals shared a major RPW8 genotype, no single allele was uniquely associated with resistance. It is concluded that RPW8 is an essential component of resistance to powdery mildews in A. thaliana, but not the only genetic factor involved in this process. No signature of selection was detected at RPW8 with a microsatellite multilocus test using an empirical null model. Unlike many previous studies of this model plant species, we found high levels of genetic diversity and relatively low differentiation (F_ST = 0.31) between populations at 14 microsatellite markers. This is judged to be due to our sampling being aimed at potentially long established populations and highlights the importance of population choice for studies of genetic diversity within this species.     

Natural variation involving deletion alleles of FRIGIDA modulate temperature‐sensitive flowering responses in Arabidopsis thaliana

Ambient temperature is one of the major environmental factors that modulate plant growth and development. There is extensive natural genetic variation in thermal responses of plants exemplified by the variation exhibited by the accessions of Arabidopsis thaliana. In this work we have studied the enhanced temperature response in hypocotyl elongation and flowering shown by the Tsu‐0 accession in long days. Genetic mapping in the Col‐0 × Tsu‐0 recombinant inbred line (RIL) population identified several QTLs for thermal response including three major effect loci encompassing candidate genes FRIGIDA (FRI), FLOWERING LOCUS C (FLC) and FLOWERING LOCUS T (FT). We confirm and validate these QTLs. We show that the Tsu‐0 FRI allele, which is the same as FRI‐Ler is associated with late flowering but only at lower temperatures in long days. Using transgenic lines and accessions, we show that the FRI‐Ler allele confers temperature‐sensitive late flowering confirming a role for FRI in photoperiod‐dependent thermal response. Through quantitative complementation with heterogeneous inbred families, we further show that cis‐regulatory variation at FT contributes to the observed hypersensitivity of Tsu‐0 to ambient temperature. Overall our results suggest that multiple loci that interact epistatically govern photoperiod‐dependent thermal responses of A. thaliana.     

Incidence and post-pollination mechanisms of nonrandom mating in Arabidopsis thaliana

Compatible pollinations from many different taxa display nonrandom mating. Here we describe a system for examining questions of nonrandom mating in Arabidopsis thaliana. Using this system, we demonstrate that Arabidopsis thaliana displays nonrandom mating between distinct accessions. Statistical analysis of these data demonstrates aspects of both pollen competition and male–female complementarity in these matings. Cytological experiments implicate pollen germination and pollen tube growth rates as possible causal factors in these nonrandom mating efficiencies.     

Additive inheritance of histone modifications in Arabidopsis thaliana intra-specific hybrids

Plant genomes are earmarked with defined patterns of chromatin marks. Little is known about the stability of these epigenomes when related, but distinct genomes are brought together by intra‐species hybridization. Arabidopsis thaliana accessions and their reciprocal hybrids were used as a model system to investigate the dynamics of histone modification patterns. The genome‐wide distribution of histone modifications H3K4me2 and H3K27me3 in the inbred parental accessions Col‐0, C24 and Cvi and their hybrid offspring was compared by chromatin immunoprecipitation in combination with genome tiling array hybridization. The analysis revealed that, in addition to DNA sequence polymorphisms, chromatin modification variations exist among accessions of A. thaliana. The range of these variations was higher for H3K27me3 (typically a repressive mark) than for H3K4me2 (typically an active mark). H3K4me2 and H3K27me3 were rather stable in response to intra‐species hybridization, with mainly additive inheritance in hybrid offspring. In conclusion, intra‐species hybridization does not result in gross changes to chromatin modifications.     

Imputation of 3 million SNPs in the Arabidopsis regional mapping population

Natural variation has become a prime resource to identify genetic variants that contribute to phenotypic variation. The regional mapping (RegMap) population is one of the most important populations for studying natural variation in Arabidopsis thaliana, and has been used in a large number of association studies and in studies on climatic adaptation. However, only 413 RegMap accessions have been completely sequenced, as part of the 1001 Genomes (1001G) Project, while the remaining 894 accessions have only been genotyped with the Affymetrix 250k chip. As a consequence, most association studies involving the RegMap are either restricted to the sequenced accessions, reducing power, or rely on a limited set of SNPs. Here we impute millions of SNPs to the 894 accessions that are exclusive to the RegMap, using the 1135 accessions of the 1001G Project as the reference panel. We assess imputation accuracy using a novel cross‐validation scheme, which we show provides a more reliable measure of accuracy than existing methods. After filtering out low accuracy SNPs, we obtain high‐quality genotypic information for 2029 accessions and 3 million markers. To illustrate the benefits of these imputed data, we reconducted genome‐wide association studies on five stress‐related traits and could identify novel candidate genes.