Differential expression and evolutionary diversification of RNA helicases in Boechera sexual and apomictic reproduction.

In higher plants, sexual reproduction is characterized by meiosis of the first cells of the germlines, and double fertilization of the egg and central cell after gametogenesis. In contrast, in apomicts of the genus Boechera, meiosis is omitted or altered and only the central cell requires fertilization, while the embryo forms parthenogenetically from the egg cell. To deepen the understanding of the transcriptional basis underlying these differences, we applied RNA-seq to compare expression in reproductive tissues of different Boechera accessions. This confirmed previous evidence of an enrichment of RNA helicases in plant germlines. Furthermore, few RNA helicases were differentially expressed in female reproductive ovule tissues harboring mature gametophytes from apomictic and sexual accessions. For some of these genes, we further found evidence for a complex recent evolutionary history. This included a homolog of Arabidopsis thaliana FASCIATED STEM4 (FAS4). In contrast to AtFAS4, which is a single-copy gene, FAS4 is represented by three homologs in Boechera, suggesting a potential for subfunctionalization to modulate reproductive development. To gain first insights into functional roles of FAS4, we studied Arabidopsis lines carrying mutant alleles. This identified the crucial importance of AtFAS4 for reproduction, as we observed developmental defects and arrest during male and female gametogenesis.

Global Brassicaceae phylogeny based on filtering of 1,000-gene dataset.

The mustard family (Brassicaceae) is a scientifically and economically important family, containing the model plant Arabidopsis thaliana and numerous crop species that feed billions worldwide. Despite its relevance, most phylogenetic trees of the family are incompletely sampled and often contain poorly supported branches. Here, we present the most complete Brassicaceae genus-level family phylogenies to date (Brassicaceae Tree of Life or BrassiToL) based on nuclear (1,081 genes, 319 of the 349 genera; 57 of the 58 tribes) and plastome (60 genes, 265 genera; all tribes) data. We found cytonuclear discordance between the two, which is likely a result of rampant hybridization among closely and more distantly related lineages. To evaluate the impact of such hybridization on the nuclear phylogeny reconstruction, we performed five different gene sampling routines, which increasingly removed putatively paralog genes. Our cleaned subset of 297 genes revealed high support for the tribes, whereas support for the main lineages (supertribes) was moderate. Calibration based on the 20 most clock-like nuclear genes suggests a late Eocene to late Oligocene origin of the family. Finally, our results strongly support a recently published new family classification, dividing the family into two subfamilies (one with five supertribes), together representing 58 tribes. This includes five recently described or re-established tribes, including Arabidopsideae, a monogeneric tribe accommodating Arabidopsis without any close relatives. With a worldwide community of thousands of researchers working on Brassicaceae and its diverse members, our new genus-level family phylogeny will be an indispensable tool for studies on biodiversity and plant biology.

Transposition and duplication of MADS-domain transcription factor genes in annual and perennial Arabis species modulates flowering.

The timing of reproduction is an adaptive trait in many organisms. In plants, the timing, duration, and intensity of flowering differ between annual and perennial species. To identify interspecies variation in these traits, we studied introgression lines derived from hybridization of annual and perennial species, Arabis montbretiana and Arabis alpina, respectively. Recombination mapping identified two tandem A. montbretiana genes encoding MADS-domain transcription factors that confer extreme late flowering on A. alpina These genes are related to the MADS AFFECTING FLOWERING (MAF) cluster of floral repressors of other Brassicaceae species and were named A. montbretiana (Am) MAF-RELATED (MAR) genes. AmMAR1 but not AmMAR2 prevented floral induction at the shoot apex of A. alpina, strongly enhancing the effect of the MAF cluster, and MAR1 is absent from the genomes of all A. alpina accessions analyzed. Exposure of plants to cold (vernalization) represses AmMAR1 transcription and overcomes its inhibition of flowering. Assembly of the tandem arrays of MAR and MAF genes of six A. alpina accessions and three related species using PacBio long-sequence reads demonstrated that the MARs arose within the Arabis genus by interchromosomal transposition of a MAF1-like gene followed by tandem duplication. Time-resolved comparative RNA-sequencing (RNA-seq) suggested that AmMAR1 may be retained in A. montbretiana to enhance the effect of the AmMAF cluster and extend the duration of vernalization required for flowering. Our results demonstrate that MAF genes transposed independently in different Brassicaceae lineages and suggest that they were retained to modulate adaptive flowering responses that differ even among closely related species.

Glucosinolate profiles and phylogeny in Barbarea compared to other tribe Cardamineae (Brassicaceae) and Reseda (Resedaceae), based on a library of ion trap HPLC-MS/MS data of reference desulfoglucosinolates.

A library of ion trap MS2 spectra and HPLC retention times reported here allowed distinction in plants of at least 70 known glucosinolates (GSLs) and some additional proposed GSLs. We determined GSL profiles of selected members of the tribe Cardamineae (Brassicaceae) as well as Reseda (Resedaceae) used as outgroup in evolutionary studies. We included several accessions of each species and a range of organs, and paid attention to minor peaks and GSLs not detected. In this way, we obtained GSL profiles of Barbarea australis, Barbarea grayi, Planodes virginica selected for its apparent intermediacy between Barbarea and the remaining tribe and family, and Rorippa sylvestris and Nasturtium officinale, for which the presence of acyl derivatives of GSLs was previously untested. We also screened Armoracia rusticana, with a remarkably diverse GSL profile, the emerging model species Cardamine hirsuta, for which we discovered a GSL polymorphism, and Reseda luteola and Reseda odorata. The potential for aliphatic GSL biosynthesis in Barbarea vulgaris was of interest, and we subjected P-type and G-type B. vulgaris to several induction regimes in an attempt to induce aliphatic GSL. However, aliphatic GSLs were not detected in any of the B. vulgaris types. We characterized the investigated chemotypes phylogenetically, based on nuclear rDNA internal transcribed spacer (ITS) sequences, in order to understand their relation to the species B. vulgaris in general, and found them to be representative of the species as it occurs in Europe, as far as documented in available ITS-sequence repositories. In short, we provide GSL profiles of a wide variety of tribe Cardamineae plants and conclude aliphatic GSLs to be absent or below our limit of detection in two major evolutionary lines of B. vulgaris. Concerning analytical chemistry, we conclude that availability of authentic reference compounds or reference materials is critical for reliable GSL analysis and characterize two publicly available reference materials: seeds of P. virginica and N. officinale.

Comparison of glucosinolate diversity in the crucifer tribe Cardamineae and the remaining order Brassicales highlights repetitive evolutionary loss and gain of biosynthetic steps.

We review glucosinolate (GSL) diversity and analyze phylogeny in the crucifer tribe Cardamineae as well as selected species from Brassicaceae (tribe Brassiceae) and Resedaceae. Some GSLs occur widely, while there is a scattered distribution of many less common GSLs, tentatively sorted into three classes: ancient, intermediate and more recently evolved. The number of conclusively identified GSLs in the tribe (53 GSLs) constitute 60% of all GSLs known with certainty from any plant (89 GSLs) and apparently unique GSLs in the tribe constitute 10 of those GSLs conclusively identified (19%). Intraspecific, qualitative GSL polymorphism is known from at least four species in the tribe. The most ancient GSL biosynthesis in Brassicales probably involved biosynthesis from Phe, Val, Leu, Ile and possibly Trp, and hydroxylation at the ß-position. From a broad comparison of families in Brassicales and tribes in Brassicaceae, we estimate that a common ancestor of the tribe Cardamineae and the family Brassicaceae exhibited GSL biosynthesis from Phe, Val, Ile, Leu, possibly Tyr, Trp and homoPhe (ancient GSLs), as well as homologs of Met and possibly homoIle (intermediate age GSLs). From the comparison of phylogeny and GSL diversity, we also suggest that hydroxylation and subsequent methylation of indole GSLs and usual modifications of Met-derived GSLs (formation of sulfinyls, sulfonyls and alkenyls) occur due to conserved biochemical mechanisms and was present in a common ancestor of the family. Apparent loss of homologs of Met as biosynthetic precursors was deduced in the entire genus Barbarea and was frequent in Cardamine (e.g. C. pratensis, C. diphylla, C. concatenata, possibly C. amara). The loss was often associated with appearance of significant levels of unique or rare GSLs as well as recapitulation of ancient types of GSLs. Biosynthetic traits interpreted as de novo evolution included hydroxylation at rare positions, acylation at the thioglucose and use of dihomoIle and possibly homoIle as biosynthetic precursors. Biochemical aspects of the deduced evolution are discussed and testable hypotheses proposed. Biosyntheses from Val, Leu, Ile, Phe, Trp, homoPhe and homologs of Met are increasingly well understood, while GSL biosynthesis from mono- and dihomoIle is poorly understood. Overall, interpretation of known diversity suggests that evolution of GSL biosynthesis often seems to recapitulate ancient biosynthesis. In contrast, unprecedented GSL biosynthetic innovation seems to be rare.

Cytokinin-promoted secondary growth and nutrient storage in the perennial stem zone of Arabis alpina.

Perennial plants maintain their lifespan through several growth seasons. Arabis alpina serves as a model Brassicaceae species to study perennial traits. Lateral stems of A. alpina have a proximal vegetative zone with a dormant bud zone and a distal senescing seed-producing inflorescence zone. We addressed how this zonation is distinguished at the anatomical level, whether it is related to nutrient storage and which signals affect the zonation. We found that the vegetative zone exhibits secondary growth, which we termed the perennial growth zone (PZ). High-molecular-weight carbon compounds accumulate there in cambium and cambium derivatives. Neither vernalization nor flowering were requirements for secondary growth and the sequestration of storage compounds. The inflorescence zone with only primary growth, termed the annual growth zone (AZ), or roots exhibited different storage characteristics. Following cytokinin application cambium activity was enhanced and secondary phloem parenchyma was formed in the PZ and also in the AZ. In transcriptome analysis, cytokinin-related genes represented enriched gene ontology terms and were expressed at a higher level in the PZ than in the AZ. Thus, A. alpina primarily uses the vegetative PZ for nutrient storage, coupled to cytokinin-promoted secondary growth. This finding lays a foundation for future studies addressing signals for perennial growth.

The Arabidopsis V-ATPase is localized to the TGN/EE via a seed plant-specific motif.

The V-ATPase is a versatile proton-pump found in a range of endomembrane compartments yet the mechanisms governing its differential targeting remain to be determined. In Arabidopsis, VHA-a1 targets the V-ATPase to the TGN/EE whereas VHA-a2 and VHA-a3 are localized to the tonoplast. We report here that the VHA-a1 targeting domain serves as both an ER-exit and as a TGN/EE-retention motif and is conserved among seed plants. In contrast, Marchantia encodes a single VHA-isoform that localizes to the TGN/EE and the tonoplast in Arabidopsis. Analysis of CRISPR/Cas9 generated null alleles revealed that VHA-a1 has an essential function for male gametophyte development but acts redundantly with the tonoplast isoforms during vegetative growth. We propose that in the absence of VHA-a1, VHA-a3 is partially re-routed to the TGN/EE. Our findings contribute to understanding the evolutionary origin of V-ATPase targeting and provide a striking example that differential localization does not preclude functional redundancy.

Nested whole-genome duplications coincide with diversification and high morphological disparity in Brassicaceae.

Angiosperms have become the dominant terrestrial plant group by diversifying for ~145 million years into a broad range of environments. During the course of evolution, numerous morphological innovations arose, often preceded by whole genome duplications (WGD). The mustard family (Brassicaceae), a successful angiosperm clade with ~4000 species, has been diversifying into many evolutionary lineages for more than 30 million years. Here we develop a species inventory, analyze morphological variation, and present a maternal, plastome-based genus-level phylogeny. We show that increased morphological disparity, despite an apparent absence of clade-specific morphological innovations, is found in tribes with WGDs or diversification rate shifts. Both are important processes in Brassicaceae, resulting in an overall high net diversification rate. Character states show frequent and independent gain and loss, and form varying combinations. Therefore, Brassicaceae pave the way to concepts of phylogenetic genome-wide association studies to analyze the evolution of morphological form and function.

The Quaternary evolutionary history of Bristol rock cress (Arabis scabra, Brassicaceae), a Mediterranean element with an outpost in the north-western Atlantic region.

BACKGROUND AND AIMS: Bristol rock cress is among the few plant species in the British Isles considered to have a Mediterranean-montane element. Spatiotemporal patterns of colonization of the British Isles since the last interglacial and after the Last Glacial Maximum (LGM) from mainland Europe are underexplored and have not yet included such floristic elements. Here we shed light on the evolutionary history of a relic and outpost metapopulation of Bristol rock cress in the south-western UK. METHODS: Amplified fragment length polymorphisms (AFLPs) were used to identify distinct gene pools. Plastome assembly and respective phylogenetic analysis revealed the temporal context. Herbarium material was largely used to exemplify the value of collections to obtain a representative sampling covering the entire distribution range. KEY RESULTS: The AFLPs recognized two distinct gene pools, with the Iberian Peninsula as the primary centre of genetic diversity and the origin of lineages expanding before and after the LGM towards mountain areas in France and Switzerland. No present-day lineages are older than 51 ky, which is in sharp contrast to the species stem group age of nearly 2 My, indicating severe extinction and bottlenecks throughout the Pleistocene. The British Isles were colonized after the LGM and feature high genetic diversity. CONCLUSIONS: The short-lived perennial herb Arabis scabra, which is restricted to limestone, has expanded its distribution range after the LGM, following corridors within an open landscape, and may have reached the British Isles via the desiccated Celtic Sea at about 16 kya. This study may shed light on the origin of other rare and peculiar species co-occurring in limestone regions in the south-western British Isles.

Interspecies association mapping links reduced CG to TG substitution rates to the loss of gene-body methylation.

Comparative genomics can unravel the genetic basis of species differences; however, successful reports on quantitative traits are still scarce. Here we present genome assemblies of 31 so-far unassembled Brassicaceae plant species and combine them with 16 previously published assemblies to establish the Brassicaceae Diversity Panel. Using a new interspecies association strategy for quantitative traits, we found a so-far unknown association between the unexpectedly high variation in CG to TG substitution rates in genes and the absence of CHROMOMETHYLASE3 (CMT3) orthologues. Low substitution rates were associated with the loss of CMT3, while species with conserved CMT3 orthologues showed high substitution rates. Species without CMT3 also lacked gene-body methylation (gbM), suggesting an evolutionary trade-off between the unknown function of gbM and low substitution rates in Brassicaceae, possibly due to low mutability of non-methylated cytosines.

A regulatory circuit conferring varied flowering response to cold in annual and perennial plants.

The reproductive strategies of plants are highly variable. Short-lived annuals flower abundantly soon after germination, whereas longer-lived perennials postpone and spatially restrict flowering. We used CRISPR/Cas9 and interspecies gene transfer to understand divergence in reproductive patterns between annual and perennial crucifers. We show that in perennial Arabis alpina, flowering in response to winter cold depends on the floral integrator SQUAMOSA PROMOTER BINDING PROTEIN-LIKE 15 (SPL15), whose activity is limited to older shoots and branches during cold exposure. In annuals, this regulatory system is conserved, but cold-induced flowering occurs in young shoots, without requirement for SPL15, through the photoperiodic pathway when plants return to warm. By reconstructing the annual response in perennials, we conclude that characteristic patterns of reproduction in annuals and perennials are conferred through variation in dependency on distinct flowering pathways acting in parallel.

Improving and correcting the contiguity of long-read genome assemblies of three plant species using optical mapping and chromosome conformation capture data.

Long-read sequencing can overcome the weaknesses of short reads in the assembly of eukaryotic genomes; however, at present additional scaffolding is needed to achieve chromosome-level assemblies. We generated Pacific Biosciences (PacBio) long-read data of the genomes of three relatives of the model plant Arabidopsis thaliana and assembled all three genomes into only a few hundred contigs. To improve the contiguities of these assemblies, we generated BioNano Genomics optical mapping and Dovetail Genomics chromosome conformation capture data for genome scaffolding. Despite their technical differences, optical mapping and chromosome conformation capture performed similarly and doubled N50 values. After improving both integration methods, assembly contiguity reached chromosome-arm-levels. We rigorously assessed the quality of contigs and scaffolds using Illumina mate-pair libraries and genetic map information. This showed that PacBio assemblies have high sequence accuracy but can contain several misassemblies, which join unlinked regions of the genome. Most, but not all, of these misjoints were removed during the integration of the optical mapping and chromosome conformation capture data. Even though none of the centromeres were fully assembled, the scaffolds revealed large parts of some centromeric regions, even including some of the heterochromatic regions, which are not present in gold standard reference sequences.

Pinpointing genes underlying annual/perennial transitions with comparative genomics.

BACKGROUND: Transitions between perennial and an annual life history occur often in plant lineages, but the genes that control whether a plant is an annual or perennial are largely unknown. To identify genes that confer differences between annuals and perennials we compared the gene content of four pairs of sister lineages (Arabidopsis thaliana/Arabidopsis lyrata, Arabis montbretiana/Arabis alpina, Arabis verna/Aubrieta parviflora and Draba nemorosa/Draba hispanica) in the Brassicaceae in which each pair contains one annual and one perennial, plus one extra annual species (Capsella rubella). RESULTS: After sorting all genes in all nine species into gene families, we identified five families in which well-annotated genes are present in the perennials A. lyrata and A. alpina, but are not present in any of the annual species. For the eleven genes in perennials in these families, an orthologous pseudogene or otherwise highly diverged gene was found in the syntenic region of the annual species in six cases. The five candidate families identified encode: a kinase, an oxidoreductase, a lactoylglutathione lyase, a F-box protein and a zinc finger protein. By comparing the active gene in the perennial to the pseudogene or heavily altered gene in the annual, dN and dS were calculated. The low dN/dS values in one kinase suggest that it became pseudogenized more recently, while the other kinase, F-box, oxidoreductase and zinc-finger became pseudogenized closer to the divergence between the annual-perennial pair. CONCLUSIONS: We identified five gene families that may be involved in the life history switch from perennial to annual. Considering the dN and dS data and whether syntenic pseudogenes were found and the potential functions of the genes, the F-box family is considered the most promising candidate for future functional studies to determine if it affects life history.

On the origin and evolutionary consequences of gene body DNA methylation.

In plants, CG DNA methylation is prevalent in the transcribed regions of many constitutively expressed genes (gene body methylation; gbM), but the origin and function of gbM remain unknown. Here we report the discovery that Eutrema salsugineum has lost gbM from its genome, to our knowledge the first instance for an angiosperm. Of all known DNA methyltransferases, only CHROMOMETHYLASE 3 (CMT3) is missing from E. salsugineum Identification of an additional angiosperm, Conringia planisiliqua, which independently lost CMT3 and gbM, supports that CMT3 is required for the establishment of gbM. Detailed analyses of gene expression, the histone variant H2A.Z, and various histone modifications in E. salsugineum and in Arabidopsis thaliana epigenetic recombinant inbred lines found no evidence in support of any role for gbM in regulating transcription or affecting the composition and modification of chromatin over evolutionary timescales.

Hybrid apomicts trapped in the ecological niches of their sexual ancestors.

Asexual reproduction is expected to reduce the adaptive potential to novel or changing environmental conditions, restricting or altering the ecological niche of asexual lineages. Asexual lineages of plants and animals are typically polyploid, an attribute that may influence their genetic variation, plasticity, adaptive potential, and niche breadth. The genus Boechera (Brassicaceae) represents an ideal model to test the relative ecological and biogeographic impacts of reproductive mode and ploidy because it is composed of diploid sexual and both diploid and polyploid asexual (i.e., apomictic) lineages. Here, we demonstrate a strong association between a transcriptionally conserved allele and apomictic seed formation. We then use this allele as a proxy apomixis marker in 1,649 accessions to demonstrate that apomixis is likely to be a common feature across the Boechera phylogeny. Phylogeographic analyses of these data demonstrate (i) species-specific niche differentiation in sexuals, (ii) extensive niche conservation between differing reproductive modes of the same species, (iii) ploidy-specific niche differentiation within and among species, and (iv) occasional niche drift between apomicts and their sexual ancestors. We conclude that ploidy is a substantially stronger and more common driver of niche divergence within and across Boechera species although variation in both traits may not necessarily lead to niche evolution on the species scale.

Genome expansion of Arabis alpina linked with retrotransposition and reduced symmetric DNA methylation.

Despite evolutionary conserved mechanisms to silence transposable element activity, there are drastic differences in the abundance of transposable elements even among closely related plant species. We conducted a de novo assembly for the 375 Mb genome of the perennial model plant, Arabis alpina. Analysing this genome revealed long-lasting and recent transposable element activity predominately driven by Gypsy long terminal repeat retrotransposons, which extended the low-recombining pericentromeres and transformed large formerly euchromatic regions into repeat-rich pericentromeric regions. This reduced capacity for long terminal repeat retrotransposon silencing and removal in A. alpina co-occurs with unexpectedly low levels of DNA methylation. Most remarkably, the striking reduction of symmetrical CG and CHG methylation suggests weakened DNA methylation maintenance in A. alpina compared with Arabidopsis thaliana. Phylogenetic analyses indicate a highly dynamic evolution of some components of methylation maintenance machinery that might be related to the unique methylation in A. alpina.

On the origin and evolution of apomixis in Boechera.

The genetic mechanisms causing seed development by gametophytic apomixis in plants are predominantly unknown. As apomixis is consistently associated with hybridity and polyploidy, these confounding factors may either (a) be the underlying mechanism for the expression of apomixis, or (b) obscure the genetic factors which cause apomixis. To distinguish between these hypotheses, we analyzed the population genetic patterns of diploid and triploid apomictic lineages and their sexual progenitors in the genus Boechera (Brassicaceae). We find that while triploid apomixis is associated with hybridization, the majority of diploid apomictic lineages are likely the product of intra-specific crosses. We then show that these diploid apomicts are more likely to sire triploid apomictic lineages than conspecific sexuals. Combined with flow cytometric seed screen phenotyping for male and female components of apomixis, our analyses demonstrate that hybridization is an indirect correlate of apomixis in Boechera.

A continental-wide perspective: the genepool of nuclear encoded ribosomal DNA and single-copy gene sequences in North American Boechera (Brassicaceae).

74 of the currently accepted 111 taxa of the North American genus Boechera (Brassicaceae) were subject to pyhlogenetic reconstruction and network analysis. The dataset comprised 911 accessions for which ITS sequences were analyzed. Phylogenetic analyses yielded largely unresolved trees. Together with the network analysis confirming this result this can be interpreted as an indication for multiple, independent, and rapid diversification events. Network analyses were superimposed with datasets describing i) geographical distribution, ii) taxonomy, iii) reproductive mode, and iv) distribution history based on phylogeographic evidence. Our results provide first direct evidence for enormous reticulate evolution in the entire genus and give further insights into the evolutionary history of this complex genus on a continental scale. In addition two novel single-copy gene markers, orthologues of the Arabidopsis thaliana genes At2g25920 and At3g18900, were analyzed for subsets of taxa and confirmed the findings obtained through the ITS data.

Systematics and evolution of Arctic-Alpine Arabis alpina (Brassicaceae) and its closest relatives in the eastern Mediterranean.

PREMISE OF THE STUDY: The high mountains in southern Anatolia and the eastern Mediterranean are assumed to play a major role as a primary center of genetic diversity and species richness in Eurasia. We tested this hypothesis by focusing on the widespread perennial arctic-alpine Arabis alpina and its sympatrically distributed closest relatives in the eastern Mediterranean. METHODS: Plastid (trnL intron, trnL-F intergenic spacer) and nuclear (ITS) DNA sequence analysis was used for phylogenetic reconstruction. Broad-scale plastid haplotype analyses were conducted to infer ancestral biogeographic patterns. KEY RESULTS: Five Arabis species, identified from the eastern Mediterranean (Turkey mainland and Cyprus), evolved directly and independently from A. alpina, leaving Arabis alpina as a paraphyletic taxon. These species are not affected by hybridization or introgression, and species divergence took place at the diploid level during the Pleistocene. CONCLUSIONS: Pleistocene climate fluctuations produced local altitudinal range-shifts among mountain glacial survival areas, resulting not only in the accumulation of intraspecific genotype diversity but also in the formation of five local species. We also show that the closest sister group of Arabis alpina consists exclusively of annuals/winter annuals and diverged prior to Pleistocene climatic fluctuations during the colonization of the lowland Mediterranean landscape. These findings highlight that Anatolia is not only a center of species richness but also a center for life-history diversification.

Aa TFL1 confers an age-dependent response to vernalization in perennial Arabis alpina.

Flowering of many plants is induced by environmental signals, but these responses can depend on the age of the plant. Exposure of Arabidopsis thaliana to vernalization (winter temperatures) at germination induces flowering, whereas a close perennial relative Arabis alpina only responds if exposed when at least 5 weeks old. We show that vernalization of these older A. alpina plants reduces expression of the floral repressor PEP1 and activates the orthologs of the Arabidopsis flowering genes SOC1 (Aa SOC1) and LFY (Aa LFY). By contrast, when younger plants are vernalized, PEP1 and Aa SOC1 mRNA levels change as in older plants, but Aa LFY is not expressed. We demonstrate that A. alpina TFL1 (Aa TFL1) blocks flowering and prevents Aa LFY expression when young plants are exposed to vernalization. In addition, in older plants, Aa TFL1 increases the duration of vernalization required for Aa LFY expression and flowering. Aa TFL1 has similar functions in axillary shoots, thus ensuring that following a flowering episode vegetative branches are maintained to continue the perennial life cycle. We propose that Aa TFL1 blocks flowering of young plants exposed to vernalization by setting a threshold for a flowering pathway that is increased in activity as the shoot ages, thus contributing to several perennial traits.