Arabidopsis thaliana and its wild relatives: a model system for ecology and evolution

SNP-RFLP Markers for the Study of Arabidopsis lyrata

Arabidopsis lyrata has become a useful system for the study of comparative genomics, hybridization, polyploidization, and evolutionary transitions from outcrossing to selfing. Previous studies of its mating system have used microsatellite loci, but low allelic diversity, particularly in self-compatible populations characterized by low levels of outcrossing, reduces the utility of these markers for more detailed studies. Here, we aimed to develop population-level SNP markers for A. lyrata ssp. lyrata sampled from a self-compatible population at Rondeau Provincial Park, Ontario, Canada. We performed de novo SNP discovery and identified 6808 putative SNPs from genome-wide sequences of 22 individuals originating from a highly selfing population. Further filtering and marker validation enabled the development of 17 SNP marker loci that can be visualized using standard PCR-RFLP protocols. These markers had average minor-allele frequencies of 0.40 in the target population, and four of seven markers were variable in a small sample from nine other A. lyrata populations. These PCR-RFLP markers have the potential to be useful for the analysis of mating patterns within and beyond the inbred self-compatible populations of A. lyrata studied here and enable the continued development of A. lyrata as a model for studying evolutionary transitions from outcrossing to selfing.

Museum genomics reveals temporal genetic stasis and global genetic diversity in Arabidopsis thaliana

Global patterns of population genetic variation through time offer a window into evolutionary processes that maintain diversity. Over time, lineages may expand or contract their distribution, causing turnover in population genetic composition. At individual loci, migration, drift, environmental change (among other processes) may affect allele frequencies. Museum specimens of widely distributed species offer a unique window into the genetics of understudied populations and changes over time. Here, we sequenced genomes of 130 herbarium specimens and 91 new field collections of Arabidopsis thaliana and combined these with published genomes. We sought a broader view of genomic diversity across the species, and to test if population genomic composition is changing through time. We documented extensive and previously uncharacterized diversity in a range of populations in Africa, populations that are under threat from anthropogenic climate change. Through time, we did not find dramatic changes in genomic composition of populations. Instead, we found a pattern of genetic change every 100 years of the same magnitude seen when comparing Eurasian populations that are 185 km apart, potentially due to a combination of drift and changing selection. We found only mixed signals of polygenic adaptation at phenology and physiology QTL. We did find that genes conserved across eudicots show altered levels of directional allele frequency change, potentially due to variable purifying and background selection. Our study highlights how museum specimens can reveal new dimensions of population diversity and show how wild populations are evolving in recent history.

Butterfly eggs prime anti-herbivore defense in an annual but not perennial Arabidopsis species

MAIN CONCLUSION

Unlike Arabidopsis thaliana, defenses of Arabidopsis lyrata against Pieris brassicae larval feeding are not primable by P. brassicae eggs. Thus, egg primability of plant anti-herbivore defenses is not phylogenetically conserved in the genus Arabidopsis. While plant anti-herbivore defenses of the annual species Arabidopsis thaliana were shown to be primable by Pieris brassicae eggs, the primability of the phylogenetically closely related perennial Arabidopsis lyrata has not yet been investigated. Previous studies revealed that closely related wild Brassicaceae plant species, the annual Brassica nigra and the perennial B. oleracea, exhibit an egg-primable defense trait, even though they have different life spans. Here, we tested whether P. brassicae eggs prime anti-herbivore defenses of the perennial A. lyrata. We exposed A. lyrata to P. brassicae eggs and larval feeding and assessed their primability by (i) determining the biomass of P. brassicae larvae after feeding on plants with and without prior P. brassicae egg deposition and (ii) investigating the plant transcriptomic response after egg deposition and/or larval feeding. For comparison, these studies were also conducted with A. thaliana. Consistent with previous findings, A. thaliana's response to prior P. brassicae egg deposition negatively affected conspecific larvae feeding upon A. thaliana. However, this was not observed in A. lyrata. Arabidopsis thaliana responded to P. brassicae eggs with strong transcriptional reprogramming, whereas A. lyrata responses to eggs were negligible. In response to larval feeding, A. lyrata exhibited a greater transcriptome change compared to A. thaliana. Among the strongly feeding-induced A. lyrata genes were those that are egg-primed in feeding-induced A. thaliana, i.e., CAX3, PR1, PR5, and PDF1.4. These results suggest that A. lyrata has evolved a robust feeding response that is independent from prior egg exposure.

Insights into the evolution of herbivory from a leaf-mining fly

Herbivorous insects and their host plants comprise most known species on Earth. Illuminating how herbivory repeatedly evolved in insects from non-herbivorous lineages is critical to understanding how this biodiversity is created and maintained. We characterized the trophic niche of Scaptomyza flava, a representative of a lineage nested within the Drosophila that transitioned to herbivory ~10-15 million years ago. We used natural history studies to determine if S. flava is a true herbivore or a cryptic microbe-feeder, given that the ancestral character state for the family Drosophilidae is likely microbe-feeding. Specifically, we quantified oviposition substrate choice and larval viability across food-types, trophic-related morphological traits, and nitrogen isotope and sterol profiles across putatively herbivorous and non-herbivorous drosophilids. The results of these studies show that S. flava is an obligate herbivore of living plants. Paired with its genetic model host, Arabidopsis thaliana, S. flava is a novel and powerful system for exploring mechanisms underlying the evolution of herbivory, a complex trait that enabled the exceptional diversification of insects.

A Hyb-Seq phylogeny of Boechera and related genera using a combination of Angiosperms353 and Brassicaceae-specific bait sets

PREMISE

Although Boechera (Boechereae, Brassicaceae) has become a plant model system for both ecological genomics and evolutionary biology, all previous phylogenetic studies have had limited success in resolving species relationships within the genus. The recent effective application of sequence data from target enrichment approaches to resolve the evolutionary relationships of several other challenging plant groups prompted us to investigate their usefulness in Boechera and Boechereae.

METHODS

To resolve the phylogeny of Boechera and closely related genera, we utilized the Hybpiper pipeline to analyze two combined bait sets: Angiosperms353, with broad applicability across flowering plants; and a Brassicaceae-specific bait set designed for use in the mustard family. Relationships for 101 samples representing 81 currently recognized species were inferred from a total of 1114 low-copy nuclear genes using both supermatrix and species coalescence methods.

RESULTS

Our analyses resulted in a well-resolved and highly supported phylogeny of the tribe Boechereae. Boechereae is divided into two major clades, one comprising all western North American species of Boechera, the other encompassing the eight other genera of the tribe. Our understanding of relationships within Boechera is enhanced by the recognition of three core clades that are further subdivided into robust regional species complexes.

CONCLUSIONS

This study presents the first broadly sampled, well-resolved phylogeny for most known sexual diploid Boechera. This effort provides the foundation for a new phylogenetically informed taxonomy of Boechera that is crucial for its continued use as a model system.

Plants cope with fluctuating light by frequency-dependent nonphotochemical quenching and cyclic electron transport

In natural environments, plants are exposed to rapidly changing light. Maintaining photosynthetic efficiency while avoiding photodamage requires equally rapid regulation of photoprotective mechanisms. We asked what the operation frequency range of regulation is in which plants can efficiently respond to varying light. Chlorophyll fluorescence, P700, plastocyanin, and ferredoxin responses of wild-types Arabidopsis thaliana were measured in oscillating light of various frequencies. We also investigated the npq1 mutant lacking violaxanthin de-epoxidase, the npq4 mutant lacking PsbS protein, and the mutants crr2-2, and pgrl1ab impaired in different pathways of the cyclic electron transport. The fastest was the PsbS-regulation responding to oscillation periods longer than 10 s. Processes involving violaxanthin de-epoxidase dampened changes in chlorophyll fluorescence in oscillation periods of 2 min or longer. Knocking out the PGR5/PGRL1 pathway strongly reduced variations of all monitored parameters, probably due to congestion in the electron transport. Incapacitating the NDH-like pathway only slightly changed the photosynthetic dynamics. Our observations are consistent with the hypothesis that nonphotochemical quenching in slow light oscillations involves violaxanthin de-epoxidase to produce, presumably, a largely stationary level of zeaxanthin. We interpret the observed dynamics of photosystem I components as being formed in slow light oscillations partially by thylakoid remodeling that modulates the redox rates.

Early lateral root formation in response to calcium and nickel shows variation within disjunct populations of Arabidopsis lyrata spp. lyrata

Root architecture is important in nutrient uptake and avoidance of toxic compounds within the soil. Arabidopsis lyrata spp. lyrata has widespread distribution in disjunct environments that encounter unique stressors starting at germination. Five populations of A. lyrata spp. lyrata show local adaptation to Nickel (Ni) but cross-tolerance to variations in Calcium (Ca) concentration within the soil. Differentiation among the populations begins early in development and appears to impact timing of lateral root formation; therefore the purpose of the study was to understand changes in root architecture and root exploration in response to Ca and Ni within the first three weeks of growth. Lateral root formation was first characterized under one concentration of Ca and Ni. Lateral root formation and tap root length were reduced in all five populations in response to Ni compared to Ca, with the least reduction in the three serpentine populations. When the populations were exposed to a gradient (either Ca or Ni) there were differences in population response based on the nature of the gradient. Start side was the greatest determinant of root exploration and lateral root formation under a Ca gradient, while population was the greatest determinant of root exploration and lateral root formation under a Ni gradient. All populations exhibited about the same frequency of root exploration under a Ca gradient, while the serpentine populations exhibited much higher levels of root exploration under a Ni gradient compared to the two non-serpentine populations. Differences among populations in response to Ca and Ni demonstrate the importance of stress responses early in development, particularly in species that have widespread distribution among disparate habitats.

Transcriptome analysis and differential expression in Arabidopsis thaliana in response to rohitukine (a chromone alkaloid) treatment

Rohitukine is a chromone alkaloid and precursor of potent anticancer drugs flavopiridol, P-276-00, and 2,6-dichloro-styryl derivative (11d) (IIIM-290). The metabolite is reported to possess anticancer, anti-inflammatory, antiadipogenic, immunomodulatory, gastroprotective, anti-implantation, antidyslipidemic, anti-arthritic, and anti-fertility properties. However, the physiological role of rohitukine in plant system is yet to be explored. Here, we studied the effect of rohitukine isolated from Dysoxylum gotadhora on Arabidopsis thaliana. The A. thaliana plants grown on a medium fortified with different rohitukine concentrations showed a significant effect on the growth and development. The root growth of A. thaliana seedlings showed considerable inhibition when grown on medium containing 1.0 mM of rohitukine. Transcriptomic analysis indicated the expression of 895 and 932 genes in control and treated samples respectively at a cut-off of FPKM ≥ 1 and P-value < 0.05. Gene ontology (GO) analysis revealed the upregulation of genes related to photosynthesis, membrane transport, antioxidation, xenobiotic degradation, and some transcription factors (TFs) in response to rohitukine. Conversely, rohitukine downregulated several genes including RNA helicases and those involved in nitrogen compound metabolism. The RNA-seq result was also validated by real-time qRT-PCR analysis. In light of these results, we discuss (i) likely ecological importance of rohitukine in parent plant as well as (ii) comparison between responses to rohitukine treatment in plants and mammals.

The selfing syndrome and beyond: diverse evolutionary consequences of mating system transitions in plants

The shift from outcrossing to self-fertilization (selfing) is considered one of the most prevalent evolutionary transitions in flowering plants. Selfing species tend to share similar reproductive traits in morphology and function, and such a set of traits is called the 'selfing syndrome'. Although the genetic basis of the selfing syndrome has been of great interest to evolutionary biologists, knowledge of the causative genes or mutations was limited until recently. Thanks to advances in population genomic methodologies combined with high-throughput sequencing technologies, several studies have successfully unravelled the molecular and genetic basis for evolution of the selfing syndrome in Capsella, Arabidopsis, Solanum and other genera. Here we first introduce recent research examples that have explored the loci, genes and mutations responsible for the selfing syndrome traits, such as reductions in petal size or in pollen production, that are mainly relevant to pre-pollination processes. Second, we review the relationship between the evolution of selfing and interspecific pollen transfer, highlighting the findings of post-pollination reproductive barriers at the molecular level. We then discuss the emerging view of patterns in evolution of the selfing syndrome, such as the pervasive involvement of loss-of-function mutations and the relative importance of selection versus neutral degradation. This article is part of the theme issue 'Genetic basis of adaptation and speciation: from loci to causative mutations'.

Globally Relaxed Selection and Local Adaptation in Boechera stricta

The strength of selection varies among populations and across the genome, but the determinants of efficacy of selection remain unclear. In this study, we used whole-genome sequencing data from 467 Boechera stricta accessions to quantify the strength of selection and characterize the pattern of local adaptation. We found low genetic diversity on 0-fold degenerate sites and conserved non-coding sites, indicating functional constraints on these regions. The estimated distribution of fitness effects and the proportion of fixed substitutions suggest relaxed negative and positive selection in B. stricta. Among the four population groups, the NOR and WES groups have smaller effective population size (Ne), higher proportions of effectively neutral sites, and lower rates of adaptive evolution compared with UTA and COL groups, reflecting the effect of Ne on the efficacy of natural selection. We also found weaker selection on GC-biased sites compared with GC-conservative (unbiased) sites, suggested that GC-biased gene conversion has affected the strength of selection in B. stricta. We found mixed evidence for the role of the recombination rate on the efficacy of selection. The positive and negative selection was stronger in high-recombination regions compared with low-recombination regions in COL but not in other groups. By scanning the genome, we found different subsets of selected genes suggesting differential adaptation among B. stricta groups. These results show that differences in effective population size, nucleotide composition, and recombination rate are important determinants of the efficacy of selection. This study enriches our understanding of the roles of natural selection and local adaptation in shaping genomic variation.

Evolution of Olfactory Receptors Tuned to Mustard Oils in Herbivorous Drosophilidae

The diversity of herbivorous insects is attributed to their propensity to specialize on toxic plants. In an evolutionary twist, toxins betray the identity of their bearers when herbivores coopt them as cues for host-plant finding, but the evolutionary mechanisms underlying this phenomenon are poorly understood. We focused on Scaptomyza flava, an herbivorous drosophilid specialized on isothiocyanate (ITC)-producing (Brassicales) plants, and identified Or67b paralogs that were triplicated as mustard-specific herbivory evolved. Using in vivo heterologous systems for the expression of olfactory receptors, we found that S. flava Or67bs, but not the homologs from microbe-feeding relatives, responded selectively to ITCs, each paralog detecting different ITC subsets. Consistent with this, S. flava was attracted to ITCs, as was Drosophila melanogaster expressing S. flava Or67b3 in the homologous Or67b olfactory circuit. ITCs were likely coopted as olfactory attractants through gene duplication and functional specialization (neofunctionalization and subfunctionalization) in S. flava, a recently derived herbivore.

Cold stress and freezing tolerance negatively affect the fitness of Arabidopsis thaliana accessions under field and controlled conditions

Higher acclimated freezing tolerance improved winter survival, but reduced reproductive fitness of Arabidopsis thaliana accessions under field and controlled conditions. Low temperature is one of the most important abiotic factors influencing plant fitness and geographical distribution. In addition, cold stress is known to influence crop yield and is therefore of great economic importance. Increased freezing tolerance can be acquired by the process of cold acclimation, but this may be associated with a fitness cost. To assess the influence of cold stress on the fitness of plants, long-term field trials over 5 years were performed with six natural accessions of Arabidopsis thaliana ranging from very tolerant to very sensitive to freezing. Fitness parameters, as seed yield and 1000 seed mass, were measured and correlation analyses with temperature and freezing tolerance data performed. The results were compared with fitness parameters from controlled chamber experiments over 3 years with application of cold priming and triggering conditions. Winter survival and seed yield per plant were positively correlated with temperature in field experiments. In addition, winter survival and 1000 seed mass were correlated with the cold-acclimated freezing tolerance of the selected Arabidopsis accessions. The results provide strong evidence for a trade-off between higher freezing tolerance and reproductive fitness in A. thaliana, which might have ecological impacts in the context of global warming.

Transcriptome profiling of Arabidopsis thaliana roots in response to allelopathic effects of Conyza canadensis

The molecular mechanisms underlying allelopathy and their role in the interactions between invasive weeds and native species remain unclear. In this study, we aimed to explore the physiological and molecular response of plant roots of a native species to allelopathy from an invasive weed. We examined the growth and development of roots of native Arabidopsis thaliana for a 2-week period after being treated with aqueous extracts at different concentrations from invasive Conyza canadensis. Extracts with higher concentration in the Murashige and Skoog (MS) media (i.e., 4 mg of extract/mL of MS) significantly affected the root growth of A. thaliana. Roots of A. thaliana displayed weakened root tip activity and an accumulation of reactive oxygen species (ROS) in response to extracts from C. canadensis. The transcriptome analysis of A. thaliana roots exposed to phytotoxicity revealed differentially expressed genes (DEGs) involved in cell wall formation, abiotic stress, transporter genes and signal transduction. We found that genes associated with nutrient transport, such as major facilitator superfamily (MFS) and amino acid permease (AAP3) transporters as well as genes involved in stress response, including leucine-rich repeat receptor-like protein kinases (LRR-RLKs) were down-regulated. In addition, we found that many transcription factors associated with plant stress (such as APETALA2/ethylene response factors) were up-regulated while others (e.g., zinc-finger proteins) were down-regulated. Allelochemicals from C. canadensis also induced the up-regulation of detoxification (DTX) genes, ROS related genes, calcineurin B-like interacting protein kinases (CIPKs) and calmodulin. Overall, our findings provided insights into allelopathy in C. canadensis at the molecular level, and contributes to our understanding of invasion mechanisms of alien plant species. CLINICAL TRIALS REGISTRATION: This study does not contain any studies with clinical trials performed by any of the authors.

Evolution of seed mass associated with mating systems in multiple plant families

In flowering plants, the evolution of self‐fertilization (selfing) from obligate outcrossing is regarded as one of the most prevalent evolutionary transitions. The evolution of selfing is often accompanied by various changes in genomic, physiological and morphological properties. In particular, a set of reproductive traits observed typically in selfing species is called the “selfing syndrome”. A mathematical model based on the kinship theory of genetic imprinting predicted that seed mass should become smaller in selfing species compared with outcrossing congeners, as a consequence of the reduced conflict between maternally and paternally derived alleles in selfing plants. Here, we test this prediction by examining the association between mating system and seed mass across a wide range of taxa (642 species), considering potential confounding factors: phylogenetic relationships and growth form. We focused on three plant families—Solanaceae, Brassicaceae and Asteraceae—where information on mating systems is abundant, and the analysis was performed for each family separately. When phylogenetic relationships were controlled, we consistently observed that selfers (represented by self‐compatible species) tended to have a smaller seed mass compared with outcrossers (represented by self‐incompatible species) in these families. In summary, our analysis suggests that small seeds should also be considered a hallmark of the selfing syndrome, although we note that mating systems have relatively small effects on seed mass variation.

Genome-wide identification and expression analysis of U-box gene family in wild emmer wheat (Triticum turgidum L. ssp. dicoccoides)

In this study, 82 U-box genes were identified in wild emmer wheat (TdPUBs) through a genome-search method. Phylogenetic analysis classified them into seven groups and the genes belonging to the same group shared the similar exon-intron structure, motif organization and cis-element compositions. Synteny analysis of the U-box genes between different species revealed that segmental duplication and polyploidization mainly contributed to the expansion of TdPUBs. Furthermore, the genetic variations of U-box were investigated in wild emmer, domesticated emmer and durum wheat. Results showed that significant genetic bottleneck has occurred during domestication process of tetraploid emmer wheat. Meanwhile, 12 TdPUBs were co-located with known domestication related QTLs. Finally, the tissue-specific and stress-responsive TdPUB genes were identified through RNA-seq analysis. Combined with qPCR validation of 19 salt-responsive TdPUBs, the candidates involving in salt response were obtained. It lays the foundation to better understand the regulatory roles of U-box family in emmer wheat and beyond.

Cross-Species Comparisons of Nanoparticle Interactions with Innate Immune Systems: A Methodological Review

Many components of the innate immune system are evolutionarily conserved and shared across many living organisms, from plants and invertebrates to humans. Therefore, these shared features can allow the comparative study of potentially dangerous substances, such as engineered nanoparticles (NPs). However, differences of methodology and procedure between diverse species and models make comparison of innate immune responses to NPs between organisms difficult in many cases. To this aim, this review provides an overview of suitable methods and assays that can be used to measure NP immune interactions across species in a multidisciplinary approach. The first part of this review describes the main innate immune defense characteristics of the selected models that can be associated to NPs exposure. In the second part, the different modes of exposure to NPs across models (considering isolated cells or whole organisms) and the main endpoints measured are discussed. In this synergistic perspective, we provide an overview of the current state of important cross-disciplinary immunological models to study NP-immune interactions and identify future research needs. As such, this paper could be used as a methodological reference point for future nano-immunosafety studies.

DeepPod: a convolutional neural network based quantification of fruit number in Arabidopsis

Background

High-throughput phenotyping based on non-destructive imaging has great potential in plant biology and breeding programs. However, efficient feature extraction and quantification from image data remains a bottleneck that needs to be addressed. Advances in sensor technology have led to the increasing use of imaging to monitor and measure a range of plants including the model Arabidopsis thaliana. These extensive datasets contain diverse trait information, but feature extraction is often still implemented using approaches requiring substantial manual input.

Results

The computational detection and segmentation of individual fruits from images is a challenging task, for which we have developed DeepPod, a patch-based 2-phase deep learning framework. The associated manual annotation task is simple and cost-effective without the need for detailed segmentation or bounding boxes. Convolutional neural networks (CNNs) are used for classifying different parts of the plant inflorescence, including the tip, base, and body of the siliques and the stem inflorescence. In a post-processing step, different parts of the same silique are joined together for silique detection and localization, whilst taking into account possible overlapping among the siliques. The proposed framework is further validated on a separate test dataset of 2,408 images. Comparisons of the CNN-based prediction with manual counting (R² = 0.90) showed the desired capability of methods for estimating silique number.

Conclusions

The DeepPod framework provides a rapid and accurate estimate of fruit number in a model system widely used by biologists to investigate many fundemental processes underlying growth and reproduction

Let research on subterranean habitats resonate!

Whereas scientists interested in subterranean life typically insist that their research is exciting, adventurous, and important to answer general questions, this enthusiasm and potential often fade when the results are translated into scientific publications. This is because cave research is often written by cave scientists for cave scientists; thus, it rarely “leaves the cave”. However, the status quo is changing rapidly. We analysed 21,486 articles focused on subterranean ecosystems published over the last three decades and observed a recent, near-exponential increase in their annual citations and impact factor. Cave research is now more often published in non-specialized journals, thanks to a number of authors who are exploiting subterranean habitats as model systems for addressing important scientific questions. Encouraged by this positive trend, we here propose a few personal ideas for improving the generality of subterranean literature, including tips for framing broadly scoped research and making it accessible to a general audience, even when published in cave-specialized journals. Hopefully, this small contribution will succeed in condensing and broadcasting even further the collective effort taken by the subterranean biology community to bring their research “outside the cave”.

Probing natural variation of IRE1 expression and endoplasmic reticulum stress responses in Arabidopsis accessions

The environmental effects shape genetic changes in the individuals within plant populations, which in turn contribute to the enhanced genetic diversity of the population as a whole. Thus, individuals within the same species can acquire and accumulate genetic differences in their genomes depending on their local environment and evolutionary history. IRE1 is a universal endoplasmic reticulum (ER) stress sensor that activates an evolutionarily conserved signalling cascade in response to biotic and abiotic stresses. Here, we selected nine different Arabidopsis accessions along with the reference ecotype Columbia-0, based on their geographical origins and differential endogenous IRE1 expression under steady-state conditions to investigate the natural variation of ER stress responses. We cloned and analysed selected upstream regulatory regions of IRE1a and IRE1b, which revealed differential levels of their inducibility. We also subjected these accessions to an array of biotic and abiotic stresses including heat, ER stress-inducing chemical tunicamycin, phytohormone salicylic acid, and pathogen infection. We measured IRE1-mediated splicing of its evolutionarily conserved downstream client as well as transcript accumulation of ER-resident chaperones and co-chaperones. Collectively, our results illustrate the expression polymorphism of a major plant stress receptor and its relationship with molecular and physiological ER stress sensitivity.

Endopolyploidy is associated with leaf functional traits and climate variation in Arabidopsis thaliana

PREMISE

Endopolyploidy is widespread throughout the tree of life and is especially prevalent in herbaceous angiosperms. Its prevalence in this clade suggests that endopolyploidy may be adaptive, but its functional roles are poorly understood. To address this gap in knowledge, we explored whether endopolyploidy was associated with climatic factors and correlated with phenotypic traits related to growth.

METHODS

We sampled stem and leaf endopolyploidy in 56 geographically separated accessions of Arabidopsis thaliana grown in a common garden to explore species variation and to determine whether this variation was correlated with climatic variables and other plant traits.

RESULTS

Stem endopolyploidy was not associated with climate or other traits. However, leaf endopolyploidy was significantly higher in accessions from drier and colder environments. Moreover, leaf endopolyploidy was positively correlated with apparent chlorophyll content and leaf dry mass.

CONCLUSIONS

Endopolyploidy may have a functional role in the storage of chloroplasts and starch, and it may offer an adaptive avenue of tissue growth in cold and dry environments.

Arabidopsis species deploy distinct strategies to cope with drought stress

BACKGROUND AND AIMS

Water limitation is an important determinant of the distribution, abundance and diversity of plant species. Yet, little is known about how the response to limiting water supply changes among closely related plant species with distinct ecological preferences. Comparison of the model annual species Arabidopsis thaliana with its close perennial relatives A. lyrata and A. halleri, can help disentangle the molecular and physiological changes contributing to tolerance and avoidance mechanisms, because these species must maintain tolerance and avoidance mechanisms to increase long-term survival, but they are exposed to different levels of water stress and competition in their natural habitat.

METHODS

A dry-down experiment was conducted to mimic a period of missing precipitation. The covariation of a progressive decrease in soil water content (SWC) with various physiological and morphological plant traits across a set of representative genotypes in A. thaliana, A. lyrata and A. halleri was quantified. Transcriptome changes to soil dry-down were further monitored.

KEY RESULTS

The analysis of trait covariation demonstrates that the three species differ in the strategies they deploy to respond to drought stress. Arabidopsis thaliana showed a drought avoidance reaction but failed to survive wilting. Arabidopsis lyrata efficiently combined avoidance and tolerance mechanisms. In contrast, A. halleri showed some degree of tolerance to wilting but it did not seem to protect itself from the stress imposed by drought. Transcriptome data collected just before plant wilting and after recovery corroborated the phenotypic analysis, with A. lyrata and A. halleri showing a stronger activation of recovery- and stress-related genes, respectively.

CONCLUSIONS

The response of the three Arabidopsis species to soil dry-down reveals that they have evolved distinct strategies to face drought stress. These strategic differences are in agreement with the distinct ecological priorities of the stress-tolerant A. lyrata, the competitive A. halleri and the ruderal A. thaliana.

Ancient polymorphisms contribute to genome-wide variation by long-term balancing selection and divergent sorting in Boechera stricta

BACKGROUND

Genomic variation is widespread, and both neutral and selective processes can generate similar patterns in the genome. These processes are not mutually exclusive, so it is difficult to infer the evolutionary mechanisms that govern population and species divergence. Boechera stricta is a perennial relative of Arabidopsis thaliana native to largely undisturbed habitats with two geographic and ecologically divergent subspecies. Here, we delineate the evolutionary processes driving the genetic diversity and population differentiation in this species.

RESULTS

Using whole-genome re-sequencing data from 517 B. stricta accessions, we identify four genetic groups that diverged around 30-180 thousand years ago, with long-term small effective population sizes and recent population expansion after the Last Glacial Maximum. We find three genomic regions with elevated nucleotide diversity, totaling about 10% of the genome. These three regions of elevated nucleotide diversity show excess of intermediate-frequency alleles, higher absolute divergence (dXY), and lower relative divergence (FST) than genomic background, and significant enrichment in immune-related genes, reflecting long-term balancing selection. Scattered across the genome, we also find regions with both high FST and dXY among the groups, termed FST-islands. Population genetic signatures indicate that FST-islands with elevated divergence, which have experienced directional selection, are derived from divergent sorting of ancient polymorphisms.

CONCLUSIONS

Our results suggest that long-term balancing selection on disease resistance genes may have maintained ancestral haplotypes across different geographical lineages, and unequal sorting of balanced polymorphisms may have generated genomic regions with elevated divergence. This study highlights the importance of ancestral balanced polymorphisms as crucial components of genome-wide variation.

Cucumber mosaic virus infection as a potential selective pressure on Arabidopsis thaliana populations

It has been proposed that in wild ecosystems viruses are often plant mutualists, whereas agroecosystems favour pathogenicity. We seek evidence for virus pathogenicity in wild ecosystems through the analysis of plant-virus coevolution, which requires a negative effect of infection on the host fitness. We focus on the interaction between Arabidopsis thaliana and Cucumber mosaic virus (CMV), which is significant in nature. We studied the genetic diversity of A. thaliana for two defence traits, resistance and tolerance, to CMV. A set of 185 individuals collected in 76 A. thaliana Iberian wild populations were inoculated with different CMV strains. Resistance was estimated from the level of virus multiplication in infected plants, and tolerance from the effect of infection on host progeny production. Resistance and tolerance to CMV showed substantial genetic variation within and between host populations, and depended on the virus x host genotype interaction, two conditions for coevolution. Resistance and tolerance were co-occurring independent traits that have evolved independently from related life-history traits involved in adaptation to climate. The comparison of the genetic structure for resistance and tolerance with that for neutral traits (QST/FST analyses) indicated that both defence traits are likely under uniform selection. These results strongly suggest that CMV infection selects for defence on A. thaliana populations, and support plant-virus coevolution. Thus, we propose that CMV infection reduces host fitness under the field conditions of the wild A. thaliana populations studied.

Genome of Crucihimalaya himalaica, a close relative of Arabidopsis, shows ecological adaptation to high altitude

Crucihimalaya himalaica, a close relative of Arabidopsis and Capsella, grows on the Qinghai-Tibet Plateau (QTP) about 4,000 m above sea level and represents an attractive model system for studying speciation and ecological adaptation in extreme environments. We assembled a draft genome sequence of 234.72 Mb encoding 27,019 genes and investigated its origin and adaptive evolutionary mechanisms. Phylogenomic analyses based on 4,586 single-copy genes revealed that C. himalaica is most closely related to Capsella (estimated divergence 8.8 to 12.2 Mya), whereas both species form a sister clade to Arabidopsis thaliana and Arabidopsis lyrata, from which they diverged between 12.7 and 17.2 Mya. LTR retrotransposons in C. himalaica proliferated shortly after the dramatic uplift and climatic change of the Himalayas from the Late Pliocene to Pleistocene. Compared with closely related species, C. himalaica showed significant contraction and pseudogenization in gene families associated with disease resistance and also significant expansion in gene families associated with ubiquitin-mediated proteolysis and DNA repair. We identified hundreds of genes involved in DNA repair, ubiquitin-mediated proteolysis, and reproductive processes with signs of positive selection. Gene families showing dramatic changes in size and genes showing signs of positive selection are likely candidates for C. himalaica's adaptation to intense radiation, low temperature, and pathogen-depauperate environments in the QTP. Loss of function at the S-locus, the reason for the transition to self-fertilization of C. himalaica, might have enabled its QTP occupation. Overall, the genome sequence of C. himalaica provides insights into the mechanisms of plant adaptation to extreme environments.

Biology in Bloom: A Primer on the Arabidopsis thaliana Model System

Arabidopsis thaliana could have easily escaped human scrutiny. Instead, Arabidopsis has become the most widely studied plant in modern biology despite its absence from the dinner table. Pairing diminutive stature and genome with prodigious resources and tools, Arabidopsis offers a window into the molecular, cellular, and developmental mechanisms underlying life as a multicellular photoautotroph. Many basic discoveries made using this plant have spawned new research areas, even beyond the verdant fields of plant biology. With a suite of resources and tools unmatched among plants and rivaling other model systems, Arabidopsis research continues to offer novel insights and deepen our understanding of fundamental biological processes.

Wild tomato endosperm transcriptomes reveal common roles of genomic imprinting in both nuclear and cellular endosperm

Genomic imprinting is a conspicuous feature of the endosperm, a triploid tissue nurturing the embryo and synchronizing angiosperm seed development. An unknown subset of imprinted genes (IGs) is critical for successful seed development and should have highly conserved functions. Recent genome-wide studies have found limited conservation of IGs among distantly related species, but there is a paucity of data from closely related lineages. Moreover, most studies focused on model plants with nuclear endosperm development, and comparisons with properties of IGs in cellular-type endosperm development are lacking. Using laser-assisted microdissection, we characterized parent-specific expression in the cellular endosperm of three wild tomato lineages (Solanum section Lycopersicon). We identified 1025 candidate IGs and 167 with putative homologs previously identified as imprinted in distantly related taxa with nuclear-type endosperm. Forty-two maternally expressed genes (MEGs) and 17 paternally expressed genes (PEGs) exhibited conserved imprinting status across all three lineages, but differences in power to assess imprinted expression imply that the actual degree of conservation might be higher than that directly estimated (20.7% for PEGs and 10.4% for MEGs). Regardless, the level of shared imprinting status was higher for PEGs than for MEGs, indicating dissimilar evolutionary trajectories. Expression-level data suggest distinct epigenetic modulation of MEGs and PEGs, and gene ontology analyses revealed MEGs and PEGs to be enriched for different functions. Importantly, our data provide evidence that MEGs and PEGs interact in modulating both gene expression and the endosperm cell cycle, and uncovered conserved cellular functions of IGs uniting taxa with cellular- and nuclear-type endosperm.

High altitude population of Arabidopsis thaliana is more plastic and adaptive under common garden than controlled condition

BACKGROUND

Population differentiation and their adaptation to a particular environment depend on their ability to respond to a new environment. This, in turn is governed to an extent, by the degree of phenotypic plasticity exhibited by the populations. The populations of same species inhabiting different climatic conditions may differ in their phenotypic plasticity. Himalayan populations of Arabidopsis thaliana originating from a steep altitude are exposed to different climatic conditions ranging from sub-tropical to temperate. Thus they might have experienced different selection pressures during evolution and may respond differently under common environmental condition.

RESULTS

Phenotypic plasticity and differentiation of natural populations of A. thaliana grown under common garden and controlled conditions were determined. A total of seventeen morphological traits, their plasticity, association between traits and environment were performed using 45 accessions from three populations. Plants from different altitudes differed in phenotypes, their selection and fitness under two conditions. Under both the conditions lower altitude population was characterized by higher leaf count and larger silique than higher and middle altitude population. Flowering time of high altitude population increased while that of low and medium altitude decreased under controlled condition compared to open field. An increase in seed weight and germination was observed for all the population under open field than controlled. Rosette area was under divergent selection in both the condition. The heritability of lower altitude population was the highest under both the conditions, where as it was the least for higher altitude population further indicating that the high altitude populations are more responsive towards phenotypic changes under new environmental conditions. Ninety-nine percent of variability in traits and their plasticity co-varied with the altitude of their origin. The population of high altitude was more plastic and differentiated as compared to the lower altitude one.

CONCLUSIONS

Arabidopsis thaliana population native to different altitudes of the west Himalaya responds differently when grown under common environments. The success of high altitude population is more in common garden than the controlled conditions. The significant variability in phenotype and its association with altitude of origin predicts for non-random genetic differentiation among the populations.

Characterization of Leaf Transcriptome in Banksia hookeriana

Banksia is a significant element in vegetation of southwestern Australia, a biodiversity hotspot with global significance. In particular, Banksia hookeriana represents a species with significant economic and ecological importance in the region. For better conservation and management, we reported an overview of transcriptome of B. hookeriana using RNA-seq and de novo assembly. We have generated a total of 202.7 million reads (18.91 billion of nucleotides) from four leaf samples in four plants of B. hookeriana, and assembled 59,063 unigenes (average size = 1098 bp) through de novotranscriptome assembly. Among them, 39,686 unigenes were annotated against the Swiss-Prot, Clusters of Orthologous Groups (COG), and NCBI non-redundant (NR) protein databases. We showed that there was approximately one single nucleotide polymorphism (SNP) per 5.6–7.1 kb in the transcriptome, and the ratio of transitional to transversional polymorphisms was approximately 1.82. We compared unigenes of B. hookeriana to those of Arabidopsis thaliana and Nelumbo nucifera through sequence homology, Gene Ontology (GO) annotation, and KEGG pathway analyses. The comparative analysis revealed that unigenes of B. hookeriana were closely related to those of N. nucifera. B. hookeriana, N. nucifera, and A. thaliana shared similar GO annotations but different distributions in KEGG pathways, indicating that B. hookeriana has adapted to dry-Mediterranean type shrublands via regulating expression of specific genes. In total 1927 potential simple sequence repeat (SSR) markers were discovered, which could be used in the genotype and genetic diversity studies of the Banksia genus. Our results provide valuable sequence resource for further study in Banksia.

Multiple alleles at a single locus control seed dormancy in Swedish Arabidopsis

Seed dormancy is a complex life history trait that determines the timing of germination and is crucial for local adaptation. Genetic studies of dormancy are challenging, because the trait is highly plastic and strongly influenced by the maternal environment. Using a combination of statistical and experimental approaches, we show that multiple alleles at the previously identified dormancy locus DELAY OF GERMINATION1 jointly explain as much as 57% of the variation observed in Swedish Arabidopsis thaliana, but give rise to spurious associations that seriously mislead genome-wide association studies unless modeled correctly. Field experiments confirm that the major alleles affect germination as well as survival under natural conditions, and demonstrate that locally adaptive traits can sometimes be dissected genetically.

Tissue Culture as a Source of Replicates in Nonmodel Plants: Variation in Cold Response in Arabidopsis lyrata ssp. petraea

While genotype-environment interaction is increasingly receiving attention by ecologists and evolutionary biologists, such studies need genetically homogeneous replicates-a challenging hurdle in outcrossing plants. This could be potentially overcome by using tissue culture techniques. However, plants regenerated from tissue culture may show aberrant phenotypes and "somaclonal" variation. Here, we examined somaclonal variation due to tissue culturing using the response to cold treatment of photosynthetic efficiency (chlorophyll fluorescence measurements for Fv/Fm, Fv'/Fm', and ΦPSII, representing maximum efficiency of photosynthesis for dark- and light-adapted leaves, and the actual electron transport operating efficiency, respectively, which are reliable indicators of photoinhibition and damage to the photosynthetic electron transport system). We compared this to variation among half-sibling seedlings from three different families of Arabidopsis lyrata ssp. petraea Somaclonal variation was limited, and we could detect within-family variation in change in chlorophyll fluorescence due to cold shock successfully with the help of tissue-culture derived replicates. Icelandic and Norwegian families exhibited higher chlorophyll fluorescence, suggesting higher performance after cold shock, than a Swedish family. Although the main effect of tissue culture on Fv/Fm, Fv'/Fm', and ΦPSII was small, there were significant interactions between tissue culture and family, suggesting that the effect of tissue culture is genotype-specific. Tissue-cultured plantlets were less affected by cold treatment than seedlings, but to a different extent in each family. These interactive effects, however, were comparable to, or much smaller than the single effect of family. These results suggest that tissue culture is a useful method for obtaining genetically homogenous replicates for studying genotype-environment interaction related to adaptively-relevant phenotypes, such as cold response, in nonmodel outcrossing plants.

Effect of Drought on Herbivore-Induced Plant Gene Expression: Population Comparison for Range Limit Inferences

Low elevation “trailing edge” range margin populations typically face increases in both abiotic and biotic stressors that may contribute to range limit development. We hypothesize that selection may act on ABA and JA signaling pathways for more stable expression needed for range expansion, but that antagonistic crosstalk prevents their simultaneous co-option. To test this hypothesis, we compared high and low elevation populations of Boechera stricta that have diverged with respect to constitutive levels of glucosinolate defenses and root:shoot ratios; neither population has high levels of both traits. If constraints imposed by antagonistic signaling underlie this divergence, one would predict that high constitutive levels of traits would coincide with lower plasticity. To test this prediction, we compared the genetically diverged populations in a double challenge drought-herbivory growth chamber experiment. Although a glucosinolate defense response to the generalist insect herbivore Spodoptera exigua was attenuated under drought conditions, the plastic defense response did not differ significantly between populations. Similarly, although several potential drought tolerance traits were measured, only stomatal aperture behavior, as measured by carbon isotope ratios, was less plastic as predicted in the high elevation population. However, RNAseq results on a small subset of plants indicated differential expression of relevant genes between populations as predicted. We suggest that the ambiguity in our results stems from a weaker link between the pathways and the functional traits compared to transcripts.

Transcriptome sequencing of Crucihimalaya himalaica (Brassicaceae) reveals how Arabidopsis close relative adapt to the Qinghai-Tibet Plateau

The extreme environment of the Qinghai-Tibet Plateau (QTP) provides an ideal natural laboratory for studies on adaptive evolution. Few genome/transcriptome based studies have been conducted on how plants adapt to the environments of QTP compared to numerous studies on vertebrates. Crucihimalaya himalaica is a close relative of Arabidopsis with typical QTP distribution, and is hoped to be a new model system to study speciation and ecological adaptation in extreme environment. In this study, we de novo generated a transcriptome sequence of C. himalaica, with a total of 49,438 unigenes. Compared to five relatives, 10,487 orthogroups were shared by all six species, and 4,286 orthogroups contain putative single copy gene. Further analysis identified 487 extremely significantly positively selected genes (PSGs) in C. himalaica transcriptome. Theses PSGs were enriched in functions related to specific adaptation traits, such as response to radiation, DNA repair, nitrogen metabolism, and stabilization of membrane. These functions are responsible for the adaptation of C. himalaica to the high radiation, soil depletion and low temperature environments on QTP. Our findings indicate that C. himalaica has evolved complex strategies for adapting to the extreme environments on QTP and provide novel insights into genetic mechanisms of highland adaptation in plants.

Epigenetics of drought-induced trans-generational plasticity: consequences for range limit development

Genetic variation gives plants the potential to adapt to stressful environments that often exist beyond their geographic range limits. However, various genetic, physiological or developmental constraints might prevent the process of adaptation. Alternatively, environmentally induced epigenetic changes might sustain populations for several generations in stressful areas across range boundaries, but previous work on Boechera stricta, an upland mustard closely related to Arabidopsis, documented a drought-induced trans-generational plastic trade-off that could contribute to range limit development. Offspring of parents who were drought treated had higher drought tolerance, but lower levels of glucosinolate toxins. Both drought tolerance and defence are thought to be needed to expand the range to lower elevations. Here, we used methylation-sensitive amplified fragment length polymorphisms to determine whether environmentally induced DNA methylation and thus epigenetics could be a mechanism involved in the observed trans-generational plastic trade-off. We compared 110 offspring from the same self-fertilizing lineages whose parents were exposed to experimental drought stress treatments in the laboratory. Using three primer combinations, 643 polymorphic epi-loci were detected. Discriminant function analysis (DFA) on the amount of methylation detected resulted in significant combinations of epi-loci that distinguished the parent drought treatments in the offspring. Principal component (PC) and univariate association analyses also detected the significant differences, even after controlling for lineage, planting flat, developmental differences and multiple testing. Univariate tests also indicated significant associations between the amount of methylation and drought tolerance or glucosinolate toxin concentration. One epi-locus that was implicated in DFA, PC and univariate association analysis may be directly involved in the trade-off because increased methylation at this site on the genome decreased drought tolerance, but increased glucosinolate concentration.

Genetic diversity and population structure of Arabidopsis thaliana along an altitudinal gradient

The natural genetic variation within a plant species is primarily a consequence of its phylogeography and evolutionary history. This variation largely determines its present-day population structure. Arabidopsis thaliana, as a model plant, has been studied in great detail including its probable origin, local as well as global genetic diversity pattern, population structure, adaptation, etc. However, no such studies have so far been reported from the Indian Himalayan region. Here, we describe a comprehensive study on the genetic diversity and population structure of A. thaliana from an altitudinal range of 700-3400 m above mean sea level the highest altitudinal range reported so far. We also compare these populations with previously reported worldwide populations. A total of 48 accessions representing six populations were analysed using 19 microsatellites and 11 chloroplast markers. Genetic diversity analysis indicated populations to be highly diverse and comparable with worldwide populations. STRUCTURE, principal coordinate and isolation by distance (IBD) analyses showed that genetic variation in different populations is structured at geographical and altitudinal level. Further analyses indicate that these populations are genetically distinct from the rest of the world populations. Different parameters of the demographic expansion model support a rapid expansion. Based on mismatch distribution, the initial time of expansion of west Himalayan populations was found to be about 130 000 years. Bayesian analysis of divergence time indicated that these populations have a long evolutionary history in this region. Based on the results of genetic diversity parameters, demographic expansion and divergence time estimation, it appears that west Himalayan populations may be the source of the west-east expansion model.

Increased heterosis in selfing populations of a perennial forb

Quantifying the importance of random genetic drift in natural populations is central to understanding the potential limits to natural selection. One approach is to estimate the magnitude of heterosis, the increased fitness of progeny derived from crosses between populations relative to crosses within populations caused by the heterozygous masking of deleterious recessive or nearly recessive alleles that have been fixed by drift within populations. Self-fertilization is expected to reduce the effective population size by half relative to outcrossing, and population bottlenecks may be common during the transition to selfing. Therefore, chance fixation of deleterious alleles due to drift in selfing populations should increase heterosis between populations. Increased homozygosity due to fixation or loss of alleles should also decrease inbreeding depression within populations. Most populations of the perennial herb Arabidopsis lyrata ssp. lyrata are self-incompatible (SI), but several have evolved self-compatibility and are highly selfing. We quantified heterosis and inbreeding depression in two predominantly self-compatible (SC) and seven SI populations in a field common garden experiment within the species' native range and examined the correlation between these metrics to gauge the similarity in their genetic basis. We measured proportion germination in the lab, and survival and fecundity (flower and seed production) for 2 years in the field, and calculated estimates of cumulative fitness. We found 7.2-fold greater heterosis in SC compared with SI populations, despite substantial heterosis in SI populations (56 %). Inbreeding depression was >61 %, and not significantly different between SC and SI populations. There was no correlation between population estimates of heterosis and inbreeding depression, suggesting that they have somewhat different genetic bases. Combined with other sources of information, our results suggest a history of bottlenecks in all of these populations. The bottlenecks in SC populations may have been severe, but their strong inbreeding depression remains enigmatic.

Timing of shoot development transitions affects degree of perenniality in Arabidopsis lyrata (Brassicaceae)

BACKGROUND

Perenniality is best understood in quantitative terms, involving the relationship between production vs. turnover of meristems, biomass, or energy reserves. Previous quantitative trait locus (QTL) studies using divergent populations of the perennial rock cress Arabidopsis lyrata have shown that trade-offs in vegetative growth vs. reproduction are due to cascading effects of differences in early vegetative development, which contribute to local adaptation. However, details of the developmental differences and how they affect perenniality remained unclear. In this study, we investigated in detail the developmental differences in perenniality between populations. A. lyrata from Norway and North Carolina populations, representing contrasting environments and degrees of perenniality, were grown under controlled conditions, and data were collected on plant phenology and shoot-level development. We tested hypotheses that differences in perenniality involve strict allocation of lateral meristems to vegetative vs. reproductive fates, or alternatively quantitative effects of pre-reproductive vegetative development.

RESULTS

The two populations showed large differences in the degree of vegetative development on individual shoots prior to reproductive transitions. The number of leaves produced on shoots prior to bolting, and not strict meristem allocation or variation in apical dominance, was able to explain variation in the number of inflorescences on individual plants. These results suggested that allocation of time to shoot vegetative vs. reproductive development could be a major factor in resource allocation differences between the populations.

CONCLUSIONS

Based on these results and those of previous QTL studies, we propose a model in which the degree of shoot vegetative development shapes the developmental context for reproduction and subsequent vegetative growth in different environments. Climate-specific effects of shoot development patterns on reproductive output and survival may result in divergent evolutionary trajectories along a perenniality continuum, which may have broader relevance for plant life history 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.

Genomic variation in Arabidopsis: tools and insights from next-generation sequencing

The release of a reference genome for Arabidopsis thaliana in 2000 has been an enormous boon for the study of plant genetics. Less than a decade later, however, a revolution in sequencing technology had enabled rapid and inexpensive re-sequencing of whole A. thaliana genomes. Large-scale efforts to characterize natural genomic variation in A. thaliana have revealed remarkable intra-specific variation in this species, ranging from single-nucleotide differences to large structural rearrangements. The partitioning of this variation by geography and local adaptation has been described using powerful new methods and tools. Simultaneously, an ambitious research agenda has emerged to sequence 1001 A. thaliana lines from around the world, while sequencing of related species is enabling powerful evolutionary genomic analyses. In this review, I summarize recent progress in genomic analysis of natural variation in A. thaliana and its close relatives. This progress has set the stage for the emergence of Arabidopsis as a model genus for evolutionary and functional genomics.

Natural variation in abiotic stress responsive gene expression and local adaptation to climate in Arabidopsis thaliana

Gene expression varies widely in natural populations, yet the proximate and ultimate causes of this variation are poorly known. Understanding how variation in gene expression affects abiotic stress tolerance, fitness, and adaptation is central to the field of evolutionary genetics. We tested the hypothesis that genes with natural genetic variation in their expression responses to abiotic stress are likely to be involved in local adaptation to climate in Arabidopsis thaliana. Specifically, we compared genes with consistent expression responses to environmental stress (expression stress responsive, "eSR") to genes with genetically variable responses to abiotic stress (expression genotype-by-environment interaction, "eGEI"). We found that on average genes that exhibited eGEI in response to drought or cold had greater polymorphism in promoter regions and stronger associations with climate than those of eSR genes or genomic controls. We also found that transcription factor binding sites known to respond to environmental stressors, especially abscisic acid responsive elements, showed significantly higher polymorphism in drought eGEI genes in comparison to eSR genes. By contrast, eSR genes tended to exhibit relatively greater pairwise haplotype sharing, lower promoter diversity, and fewer nonsynonymous polymorphisms, suggesting purifying selection or selective sweeps. Our results indicate that cis-regulatory evolution and genetic variation in stress responsive gene expression may be important mechanisms of local adaptation to climatic selective gradients.

Abiotic and biotic controls on local spatial distribution and performance of Boechera stricta

This study investigates the relative influence of biotic and abiotic factors on community dynamics using an integrated approach and highlights the influence of space on genotypic and phenotypic traits in plant community structure. We examined the relative influence of topography, environment, spatial distance, and intra- and interspecific interactions on spatial distribution and performance of Boechera stricta (rockcress), a close perennial relative of model plant Arabidopsis. First, using Bayesian kriging, we mapped the topography and environmental gradients and explored the spatial distribution of naturally occurring rockcress plants and two neighbors, Taraxacum officinale (dandelion) and Solidago missouriensis (goldenrod) found in close proximity within a typical diverse meadow community across topographic and environmental gradients. We then evaluated direct and indirect relationships among variables using Mantel path analysis and developed a network displaying abiotic and biotic interactions in this community. We found significant spatial autocorrelation among rockcress individuals, either because of common microhabitats as displayed by high density of individuals at lower elevation and high soil moisture area, or limited dispersal as shown by significant spatial autocorrelation of naturally occurring inbred lines, or a combination of both. Goldenrod and dandelion density around rockcress does not show any direct relationship with rockcress fecundity, possibly due to spatial segregation of resources. However, dandelion density around rockcress shows an indirect negative influence on rockcress fecundity via herbivory, indicating interspecific competition. Overall, we suggest that common microhabitat preference and limited dispersal are the main drivers for spatial distribution. However, intra-specific interactions and insect herbivory are the main drivers of rockcress performance in the meadow community.