Pan-European study of genotypes and phenotypes in the Arabidopsis relative Cardamine hirsuta reveals how adaptation, demography, and development shape diversity patterns.

We study natural DNA polymorphisms and associated phenotypes in the Arabidopsis relative Cardamine hirsuta. We observed strong genetic differentiation among several ancestry groups and broader distribution of Iberian relict strains in European C. hirsuta compared to Arabidopsis. We found synchronization between vegetative and reproductive development and a pervasive role for heterochronic pathways in shaping C. hirsuta natural variation. A single, fast-cycling ChFRIGIDA allele evolved adaptively allowing range expansion from glacial refugia, unlike Arabidopsis where multiple FRIGIDA haplotypes were involved. The Azores islands, where Arabidopsis is scarce, are a hotspot for C. hirsuta diversity. We identified a quantitative trait locus (QTL) in the heterochronic SPL9 transcription factor as a determinant of an Azorean morphotype. This QTL shows evidence for positive selection, and its distribution mirrors a climate gradient that broadly shaped the Azorean flora. Overall, we establish a framework to explore how the interplay of adaptation, demography, and development shaped diversity patterns of 2 related plant species.

The trichome pattern diversity of Cardamine shares genetic mechanisms with Arabidopsis but differs in environmental drivers.

Natural variation in trichome pattern (amount and distribution) is prominent among populations of many angiosperms. However, the degree of parallelism in the genetic mechanisms underlying this diversity and its environmental drivers in different species remain unclear. To address these questions, we analyzed the genomic and environmental bases of leaf trichome pattern diversity in Cardamine hirsuta, a relative of Arabidopsis (Arabidopsis thaliana). We characterized 123 wild accessions for their genomic diversity, leaf trichome patterns at different temperatures, and environmental adjustments. Nucleotide diversities and biogeographical distribution models identified two major genetic lineages with distinct demographic and adaptive histories. Additionally, C. hirsuta showed substantial variation in trichome pattern and plasticity to temperature. Trichome amount in C. hirsuta correlated positively with spring precipitation but negatively with temperature, which is opposite to climatic patterns in A. thaliana. Contrastingly, genetic analysis of C. hirsuta glabrous accessions indicated that, like for A. thaliana, glabrousness is caused by null mutations in ChGLABRA1 (ChGL1). Phenotypic genome-wide association studies (GWAS) further identified a ChGL1 haplogroup associated with low trichome density and ChGL1 expression. Therefore, a ChGL1 series of null and partial loss-of-function alleles accounts for the parallel evolution of leaf trichome pattern in C. hirsuta and A. thaliana. Finally, GWAS also detected other candidate genes (e.g. ChETC3, ChCLE17) that might affect trichome pattern. Accordingly, the evolution of this trait in C. hirsuta and A. thaliana shows partially conserved genetic mechanisms but is likely involved in adaptation to different environments.

Seasonal Regulation of Petal Number.

Four petals characterize the flowers of most species in the Brassicaceae family, and this phenotype is generally robust to genetic and environmental variation. A variable petal number distinguishes the flowers of Cardamine hirsuta from those of its close relative Arabidopsis (Arabidopsis thaliana), and allelic variation at many loci contribute to this trait. However, it is less clear whether C. hirsuta petal number varies in response to seasonal changes in environment. To address this question, we assessed whether petal number responds to a suite of environmental and endogenous cues that regulate flowering time in C. hirsuta We found that petal number showed seasonal variation in C. hirsuta, such that spring flowering plants developed more petals than those flowering in summer. Conditions associated with spring flowering, including cool ambient temperature, short photoperiod, and vernalization, all increased petal number in C. hirsuta Cool temperature caused the strongest increase in petal number and lengthened the time interval over which floral meristems matured. We performed live imaging of early flower development and showed that floral buds developed more slowly at 15°C versus 20°C. This extended phase of floral meristem formation, coupled with slower growth of sepals at 15°C, produced larger intersepal regions with more space available for petal initiation. In summary, the growth and maturation of floral buds is associated with variable petal number in C. hirsuta and responds to seasonal changes in ambient temperature.

Heterochrony underpins natural variation in Cardamine hirsuta leaf form.

A key problem in biology is whether the same processes underlie morphological variation between and within species. Here, by using plant leaves as an example, we show that the causes of diversity at these two evolutionary scales can be divergent. Some species like the model plant Arabidopsis thaliana have simple leaves, whereas others like the A. thaliana relative Cardamine hirsuta bear complex leaves comprising leaflets. Previous work has shown that these interspecific differences result mostly from variation in local tissue growth and patterning. Now, by cloning and characterizing a quantitative trait locus (QTL) for C. hirsuta leaf shape, we find that a different process, age-dependent progression of leaf form, underlies variation in this trait within species. This QTL effect is caused by cis-regulatory variation in the floral repressor ChFLC, such that genotypes with low-expressing ChFLC alleles show both early flowering and accelerated age-dependent changes in leaf form, including faster leaflet production. We provide evidence that this mechanism coordinates leaf development with reproductive timing and may help to optimize resource allocation to the next generation.

The genetic architecture of petal number in Cardamine hirsuta.

Invariant petal number is a characteristic of most flowers and is generally robust to genetic and environmental variation. We took advantage of the natural variation found in Cardamine hirsuta petal number to investigate the genetic basis of this trait in a case where robustness was lost during evolution. We used quantitative trait locus (QTL) analysis to characterize the genetic architecture of petal number. Αverage petal number showed transgressive variation from zero to four petals in five C. hirsuta mapping populations, and this variation was highly heritable. We detected 15 QTL at which allelic variation affected petal number. The effects of these QTL were relatively small in comparison with alleles induced by mutagenesis, suggesting that natural selection may act to maintain petal number within its variable range below four. Petal number showed a temporal trend during plant ageing, as did sepal trichome number, and multi-trait QTL analysis revealed that these age-dependent traits share a common genetic basis. Our results demonstrate that petal number is determined by many genes of small effect, some of which are age-dependent, and suggests a mechanism of trait evolution via the release of cryptic variation.

Stochastic variation in Cardamine hirsuta petal number.

BACKGROUND AND AIMS: Floral development is remarkably robust in terms of the identity and number of floral organs in each whorl, whereas vegetative development can be quite plastic. This canalization of flower development prevents the phenotypic expression of cryptic genetic variation, even in fluctuating environments. A cruciform perianth with four petals is a hallmark of the Brassicaceae family, typified in the model species Arabidopsis thaliana However, variable petal loss is found in Cardamine hirsuta, a genetically tractable relative of A. thaliana Cardamine hirsuta petal number varies in response to stochastic, genetic and environmental perturbations, which makes it an interesting model to study mechanisms of decanalization and the expression of cryptic variation. METHODS: Multitrait quantitative trait locus (QTL) analysis in recombinant inbred lines (RILs) was used to identify whether the stochastic variation found in C. hirsuta petal number had a genetic basis. KEY RESULTS: Stochastic variation (standard error of the average petal number) was found to be a heritable phenotype, and four QTL that influenced this trait were identified. The sensitivity to detect these QTL effects was increased by accounting for the effect of ageing on petal number variation. All QTL had significant effects on both average petal number and its standard error, indicating that these two traits share a common genetic basis. However, for some QTL, a degree of independence was found between the age of the flowers where allelic effects were significant for each trait. CONCLUSIONS: Stochastic variation in C. hirsuta petal number has a genetic basis, and common QTL influence both average petal number and its standard error. Allelic variation at these QTL can, therefore, modify petal number in an age-specific manner via effects on the phenotypic mean and stochastic variation. These results are discussed in the context of trait evolution via a loss of robustness.

Cardamine hirsuta: a versatile genetic system for comparative studies.

A major goal in biology is to identify the genetic basis for phenotypic diversity. This goal underpins research in areas as diverse as evolutionary biology, plant breeding and human genetics. A limitation for this research is no longer the availability of sequence information but the development of functional genetic tools to understand the link between changes in sequence and phenotype. Here we describe Cardamine hirsuta, a close relative of the reference plant Arabidopsis thaliana, as an experimental system in which genetic and transgenic approaches can be deployed effectively for comparative studies. We present high-resolution genetic and cytogenetic maps for C. hirsuta and show that the genome structure of C. hirsuta closely resembles the eight chromosomes of the ancestral crucifer karyotype and provides a good reference point for comparative genome studies across the Brassicaceae. We compared morphological and physiological traits between C. hirsuta and A. thaliana and analysed natural variation in stamen number in which lateral stamen loss is a species characteristic of C. hirsuta. We constructed a set of recombinant inbred lines and detected eight quantitative trait loci that can explain stamen number variation in this population. We found clear phylogeographic structure to the genetic variation in C. hirsuta, thus providing a context within which to address questions about evolutionary changes that link genotype with phenotype and the environment.

Paths to selection on life history loci in different natural environments across the native range of Arabidopsis thaliana.

Selection on quantitative trait loci (QTL) may vary among natural environments due to differences in the genetic architecture of traits, environment-specific allelic effects or changes in the direction and magnitude of selection on specific traits. To dissect the environmental differences in selection on life history QTL across climatic regions, we grew a panel of interconnected recombinant inbred lines (RILs) of Arabidopsis thaliana in four field sites across its native European range. For each environment, we mapped QTL for growth, reproductive timing and development. Several QTL were pleiotropic across environments, three colocalizing with known functional polymorphisms in flowering time genes (CRY2, FRI and MAF2-5), but major QTL differed across field sites, showing conditional neutrality. We used structural equation models to trace selection paths from QTL to lifetime fitness in each environment. Only three QTL directly affected fruit number, measuring fitness. Most QTL had an indirect effect on fitness through their effect on bolting time or leaf length. Influence of life history traits on fitness differed dramatically across sites, resulting in different patterns of selection on reproductive timing and underlying QTL. In two oceanic field sites with high prereproductive mortality, QTL alleles contributing to early reproduction resulted in greater fruit production, conferring selective advantage, whereas alleles contributing to later reproduction resulted in larger size and higher fitness in a continental site. This demonstrates how environmental variation leads to change in both QTL effect sizes and direction of selection on traits, justifying the persistence of allelic polymorphism at life history QTL across the species range.

Interspersed expression of CUP-SHAPED COTYLEDON2 and REDUCED COMPLEXITY shapes Cardamine hirsuta complex leaf form.

How genetically regulated growth shapes organ form is a key problem in developmental biology. Here, we investigate this problem using the leaflet-bearing complex leaves of Cardamine hirsuta as a model. Leaflet development requires the action of two growth-repressing transcription factors: REDUCED COMPLEXITY (RCO), a homeodomain protein, and CUP-SHAPED COTYLEDON2 (CUC2), a NAC-domain protein. However, how their respective growth-repressive actions are integrated in space and time to generate complex leaf forms remains unknown. By using live imaging, we show that CUC2 and RCO are expressed in an interspersed fashion along the leaf margin, creating a distinctive striped pattern. We find that this pattern is functionally important because forcing RCO expression in the CUC2 domain disrupts auxin-based marginal patterning and can abolish leaflet formation. By combining genetic perturbations with time-lapse imaging and cellular growth quantifications, we provide evidence that RCO-mediated growth repression occurs after auxin-based leaflet patterning and in association with the repression of cell proliferation. Additionally, through the use of genetic mosaics, we show that RCO is sufficient to repress both cellular growth and proliferation in a cell-autonomous manner. This mechanism of growth repression is different to that of CUC2, which occurs in proliferating cells. Our findings clarify how the two growth repressors RCO and CUC2 coordinate to subdivide developing leaf primordia into distinct leaflets and generate the complex leaf form. They also indicate different relationships between growth repression and cell proliferation in the patterning and post-patterning stages of organogenesis.

The cellular basis for synergy between RCO and KNOX1 homeobox genes in leaf shape diversity.

Leaves of seed plants provide an attractive system to study the development and evolution of form. Leaves show varying degrees of margin complexity ranging from simple, as in Arabidopsis thaliana, to fully dissected into leaflets in the closely related species Cardamine hirsuta. Leaflet formation requires actions of Class I KNOTTED1-LIKE HOMEOBOX (KNOX1) and REDUCED COMPLEXITY (RCO) homeobox genes, which are expressed in the leaves of C. hirsuta but not A. thaliana. Evolutionary studies indicate that diversification of KNOX1 and RCO genes was repeatedly associated with increased leaf complexity. However, whether this gene combination represents a developmentally favored avenue for leaflet formation remains unknown, and the cell-level events through which the combined action of these genes drives leaflet formation are also poorly understood. Here we show, through a genetic screen, that when a C. hirsuta RCO transgene is expressed in A. thaliana, then ectopic KNOX1 expression in leaves represents a preferred developmental path for leaflet formation. Using time-lapse growth analysis, we demonstrate that KNOX1 expression in the basal domain of leaves leads to prolonged and anisotropic cell growth. This KNOX1 action, in synergy with local growth repression by RCO, is instrumental in generating rachises and petiolules, the linear geometrical elements, that bear leaflets in complex leaves. Our results show how the combination of cell-level growth analyses and genetics can help us understand how evolutionary modifications in expression of developmentally important genes are translated into diverse leaf shapes.

CRISPR/Cas9-Mediated Mutagenesis of RCO in Cardamine hirsuta.

The small crucifer Cardamine hirsuta bears complex leaves divided into leaflets. This is in contrast to its relative, the reference plant Arabidopsis thaliana, which has simple leaves. Comparative studies between these species provide attractive opportunities to study the diversification of form. Here, we report on the implementation of the CRISPR/Cas9 genome editing methodology in C. hirsuta and with it the generation of novel alleles in the RCO gene, which was previously shown to play a major role in the diversification of form between the two species. Thus, genome editing can now be deployed in C. hirsuta, thereby increasing its versatility as a model system to study gene function and evolution.

Demographic analyses of a new sample of haploid genomes from a Swedish population of Drosophila melanogaster.

European and African natural populations of Drosophila melanogaster have been the focus of several studies aiming at inferring demographic and adaptive processes based on genetic variation data. However, in these analyses little attention has been given to gene flow between African and European samples. Here we present a dataset consisting of 14 fully sequenced haploid genomes sampled from a natural population from the northern species range (Umeå, Sweden). We co-analyzed this new data with an African population to compare the likelihood of several competing demographic scenarios for European and African populations and show that gene flow improves the fit of demographic models to data.

The role of APETALA1 in petal number robustness.

Invariant floral forms are important for reproductive success and robust to natural perturbations. Petal number, for example, is invariant in Arabidopsis thaliana flowers. However, petal number varies in the closely related species Cardamine hirsuta, and the genetic basis for this difference between species is unknown. Here we show that divergence in the pleiotropic floral regulator APETALA1 (AP1) can account for the species-specific difference in petal number robustness. This large effect of AP1 is explained by epistatic interactions: A. thaliana AP1 confers robustness by masking the phenotypic expression of quantitative trait loci controlling petal number in C. hirsuta. We show that C. hirsuta AP1 fails to complement this function of A. thaliana AP1, conferring variable petal number, and that upstream regulatory regions of AP1 contribute to this divergence. Moreover, variable petal number is maintained in C. hirsuta despite sufficient standing genetic variation in natural accessions to produce plants with four-petalled flowers.

The Cardamine hirsuta genome offers insight into the evolution of morphological diversity.

Finding causal relationships between genotypic and phenotypic variation is a key focus of evolutionary biology, human genetics and plant breeding. To identify genome-wide patterns underlying trait diversity, we assembled a high-quality reference genome of Cardamine hirsuta, a close relative of the model plant Arabidopsis thaliana. We combined comparative genome and transcriptome analyses with the experimental tools available in C. hirsuta to investigate gene function and phenotypic diversification. Our findings highlight the prevalent role of transcription factors and tandem gene duplications in morphological evolution. We identified a specific role for the transcriptional regulators PLETHORA5/7 in shaping leaf diversity and link tandem gene duplication with differential gene expression in the explosive seed pod of C. hirsuta. Our work highlights the value of comparative approaches in genetically tractable species to understand the genetic basis for evolutionary change.

Construction of a genetic linkage map of Thlaspi caerulescens and quantitative trait loci analysis of zinc accumulation.

Zinc (Zn) hyperaccumulation seems to be a constitutive species-level trait in Thlaspi caerulescens. When compared under conditions of equal Zn availability, considerable variation in the degree of hyperaccumulation is observed among accessions originating from different soil types. This variation offers an excellent opportunity for further dissection of the genetics of this trait. A T. caerulescens intraspecific cross was made between a plant from a nonmetallicolous accession [Lellingen (LE)], characterized by relatively high Zn accumulation, and a plant from a calamine accession [La Calamine (LC)], characterized by relatively low Zn accumulation. Zinc accumulation in roots and shoots segregated in the F3 population. This population was used to construct an LE/LC amplified fragment length polymorphism (AFLP)-based genetic linkage map and to map quantitative trait loci (QTL) for Zn accumulation. Two QTL were identified for root Zn accumulation, with the trait-enhancing alleles being derived from each of the parents, and explaining 21.7 and 16.6% of the phenotypic variation observed in the mapping population. Future development of more markers, based on Arabidopsis orthologous genes localized in the QTL regions, will allow fine-mapping and map-based cloning of the genes underlying the QTL.