Skip to content


  • Research article
  • Open Access

Arabidopsis thaliana ambient temperature responsive lncRNAs

BMC Plant Biology201818:145

  • Received: 31 October 2017
  • Accepted: 4 July 2018
  • Published:



Long non-coding RNAs (lncRNAs) have emerged as new class of regulatory molecules in animals where they regulate gene expression at transcriptional and post-transcriptional level. Recent studies also identified lncRNAs in plant genomes, revealing a new level of transcriptional complexity in plants. Thousands of lncRNAs have been predicted in the Arabidopsis thaliana genome, but only a few have been studied in depth.


Here we report the identification of Arabidopsis lncRNAs that are expressed during the vegetative stage of development in either the shoot apical meristem or in leaves. We found that hundreds of lncRNAs are expressed in these tissues, of which 50 show differential expression upon an increase in ambient temperature. One of these lncRNAs, FLINC, is down-regulated at higher ambient temperature and affects ambient temperature-mediated flowering in Arabidopsis.


A number of ambient temperature responsive lncRNAs were identified with potential roles in the regulation of temperature-dependent developmental changes, such as the transition from the vegetative to the reproductive (flowering) phase. The challenge for the future is to characterize the biological function and molecular mode of action of the large number of ambient temperature-regulated lncRNAs that have been identified in this study.


  • Long non-coding RNA (lncRNA)
  • Ambient temperature response
  • Flowering time
  • Arabidopsis thaliana


LncRNAs represent a new class of recently discovered regulatory molecules. A key role for lncRNAs in the transcriptional and post-transcriptional regulation of a plethora of biological processes is emerging. LncRNAs are defined as long RNA molecules (longer than 200 bp) that do not code for a protein. They have been divided in different categories based on their location in the genome relative to protein-coding genes: natural antisense transcripts (NATs) are lncRNAs transcribed on the opposite strand of a protein-coding gene; long intergenic RNAs (lincRNAs) are transcribed in intergenic regions; intronic lncRNAs (iLncRNA) are transcribed from an intron in a protein-coding gene and; promoter lncRNAs (pLncRNAs), are transcribed from the promoter region of a gene [1, 2].

In the past years, genome-wide studies have identified many lncRNAs in the Arabidopsis thaliana (Arabidopsis) genome, as well as in other plant species, including rice and maize [1]; however the biological function of most of these lncRNAs is not known [3]. Currently, only a small number of lncRNAs has been functionally characterized and the emerging picture shows that lncRNAs have several modes of action, ranging from transcriptional interference to chromatin looping and mRNA splicing [1, 2, 4]. Moreover, these lncRNAs act in a wide range of biological processes. Some examples include the lincRNA APOLO, which is involved in auxin response [5], ELENA1, which functions in plant immunity [6], Enod40, which plays a role in regulation of sucrose utilization in nodules [7] and IPS1, which acts as a mimic of a miRNA target in the control of Pi homeostasis [8].

Vernalization, a long period of cold that is required by many temperate plant species to make the transition from the vegetative to the reproductive phase, is regulated by lncRNAs. The first characterized lncRNAs in plant were the lncRNAs transcribed from the FLOWERING LOCUS C (FLC) locus: the NAT COOLAIR [9], the iLncRNA COLDAIR [10], and more recently the pLncRNA COLDWRAP [11] and NAT ASL [12]. These molecules regulate the epigenetic silencing of FLC during vernalization in response to winter cold. This, together with the fact that the key regulators of flowering time in the ambient temperature pathway, FLOWERING LOCUS M (FLM) [13, 14], MADS AFFECTING FLOWERING 2 (MAF2) [15], MAF3, MAF4 and MAF5 [16], all belong to the FLC clade [17], prompted us to speculate a possible role for lncRNAs in ambient temperature-mediated development.

How plants are able to respond in a fast and dynamic manner to continuous changes in ambient temperature is a long-standing question in plant biology research. It is also important, in light of global warming, to understand how temperature affects plant development, as temperature influences different aspects of plant biology and development, including hypocotyl elongation, leaf shape determination, flowering time control and pathogen responses [18, 19]. Current research on the role of ambient temperature in plant development is mainly focused on identification and characterization of protein coding genes involved in this process. For example, MADS-box transcription factor genes, the bHLH transcription factor PHYTOCHROME INTERACTING FACTOR 4 (PIF4) [2023], the TEMPRANILLO genes [24] and the MYB transcription factor EARLY FLOWERING MYB PROTEIN (EFM) [25] have been shown to regulate ambient temperature-mediated flowering. PIF4, together with the photoreceptor phytochrome B (phyB), also plays an important role in thermo-morphogenesis [2628]. Furthermore, there is an important role in ambient temperature-mediated responses for various circadian clock genes [2931].

Despite the clear role for lncRNAs in cold temperature signaling, their potential role as mediators of ambient temperature responses remains unexplored.


Identification of lncRNAs in Arabidopsis thaliana.

We used RNA-seq transcriptome analysis to identify lncRNAs in shoot apical meristem (SAM)-enriched tissue that might play a role in ambient temperature-mediated developmental processes. We focused on transcripts carrying a polyA tail and followed a TAIR10-guided assembly approach using cufflinks. This approach allowed the identification of both known lncRNAs and annotation of novel transcripts [32]. Strand-specific libraries were generated from shoot apical meristem (SAM)-enriched tissue of 5 weeks-old Arabidopsis Col-0 plants grown in short-day conditions, and therefore still in the vegetative stage. The SAM of plants grown in short-day conditions is larger than the SAM of those grown in long-day conditions [33], allowing the presence of a good proportion of SAM tissue in the sampled material. Using this approach we aimed to identify low-expressed transcripts in the SAM. To annotate the transcripts, the “cuffcompare” function of cufflinks was implemented, which classifies transcripts base on their relationship with the TAIR10 reference genome annotation. A similar approach was recently used successfully for the identification of lncRNAs in the tomato genome [34].

We selected transcripts that belong to the class code “x”, “i” and “u” (Table 1). The class code “x” defines “Exonic overlap with reference on the opposite strand”, and includes, among others, the Natural Antisense Transcripts (NAT). The class code “i” defines “transcribed genomic fragments falling entirely within a reference intron”, which includes the intronic lncRNAs (iLncRNA). Finally, the class code “u” defines the “unknown intergenic transcript”, e.g. the long intergenic lncRNA (LincRNAs) and the promoter lncRNAs (pLncRNAs). In total, we identified 2132 new putative lncRNA that belong to one of these categories (Table 1 and Additional file 1: Table S1).
Table 1

putative lncRNAs



cufflink class



Natural antisense transcript




long intergenic ncRNA




intronic lncRNA







Only a few lncRNAs are annotated in the TAIR10 annotation. Therefore, we retrieved an extensive list of 14,954 RNAs annotated as lncRNAs from the Plant long non-coding RNA database (PlncDB) [35] and compared these with our identified putative lncRNA transcripts. In total 644 of the 2132 lncRNA identified in our annotation were also annotated in the PlncDB, 22 iLncRNA, 92 NATs and 530 lincRNAs (Additional file 1: Table S1).

Next, we compared the lncRNAs identified in our analysis with the ones reported in the CANTATAdb 2.0 [36]. We found 321 transcripts to be also present ion the CANTATA database, 10 iLncRNA, 38 NATs and 273 lincRNAs (Additional file 1: Table S1).

A new annotation of the Arabidopsis genome, Araport11, was released recently [37]. This new annotation is based on the analysis of a collection of RNA-seq experiments performed on different tissues from plants grown under different conditions. This approach allows us to identify new lncRNA transcripts expressed in specific conditions. We investigated the overlap between our lncRNA annotation and Araport11 transcripts and found 457 lncRNAs in both data sets (Additional file 1: Table S1).

Remarkably, 1323 lncRNAs were uniquely identified in our analysis, highlighting the power of our approach in isolating potentially new lncRNAs. Only a small number of these are likely to be falsely assigned due to sequencing or annotation errors. The identification of this relatively large number of novel transcripts is likely due to the fact that lncRNA are generally lowly expressed and therefore, their identification is dependent on the experimental setup used and tissues sampled for the RNA isolation.

Because one of the key developmental processes regulated by ambient temperature is flowering time control, we initially compared the genomic location of the lncRNAs in relation to the localization of key flowering time genes. A thorough investigation revealed that 22 lncRNAs are located nearby a flowering time gene [38], such as the lincRNA AtLnc1134 located in ABRE BINDING FACTOR 4 (ABF4), the AtLnc1860, AtLnc198 and AtLnc236, located respectively in the promoters of MADS AFFECTING FLOWERING 2 [15, 39, 40], FLAVIN-BINDING, KELCH REPEAT, F Box 1(FKF1) [41, 42] and EARLY FLOWERING IN SHORT DAYS (EFS) [43] (Additional file 2: Table S2).

LncRNA expression is modulated by ambient temperature

To identify lncRNAs that are responsive to an ambient temperature change, we compared the transcriptome of plants that were grown in short day condition at 16 °C and then moved to 25 °C with the transcriptome of plants grown continuously at 16 °C. For this purpose, tissue was collected at one, three and 5 days after the temperature increase, as described previously [44]. When grown in short-day conditions meristem size is increased, providing a relative higher amount of meristematic tissue to study the effect of a change in ambient growth temperature on the expression of meristematic lncRNAs. Using this approach, we found 50 lncRNAs (11 NATs, two intronic lncRNA and 37 lincRNAs) that were significantly differentially- expressed (log2 fold change> |1|, p-value < 0.01) in at least one time point upon the temperature change (Additional file 3: Table S3). Among the temperature responsive lncRNAs are AtLnc2 (Fig. 1a), an antisense RNA of CYTOCHROME P450 FAMILY 78, SUBFAMILY A, POLYPEPTIDE 8 (CYP78A8A), which plays a role in reproduction [45], AtLnc120 a LincRNA located nearby the MADS-box transcription factor AGL97 (Fig. 1a), and AtLnc1128 located in the intron of HOMOLOGUE OF CYANOBACTERIAL RBCX 1 (RBCX1), coding for a RUBISCO chaperon involved in cold response [46, 47]. We observed that expression of 28 of the differentially expressed lncRNAs (DElnc) was repressed upon an increase in growth temperature, while 22 lncRNAs showed an increase in expression (Fig. 1b-c). The majority of DElnc (36%) are differentially expressed at all time-points, while 16% of the DElnc are only differentially expressed 1 day after the temperature switch, showing a rapid and transient temperature response. 20% of the DElnc appeared to be differentially expressed at day five only, showing a late and probably indirect response to the temperature change (Fig. 1d).
Fig. 1
Fig. 1

Identification of ambient temperature-responsive lncRNAs. a. Example of a lncRNA identified in our strand-specific RNA-seq experiment. Protein coding genes and lncRNA are shown in grey and black respectively. In blue are reads mapping to the forward strand and in purple are reads mapping to the reverse strand from RNA-seq on plants growing at 16 °C in short day. b. Heat map showing lncRNA that are significantly (adjp-value < 0.01) differentially expressed (log2 fold change > |1|) in at least in one time point upon a temperature change from 16 °C to 25 °C. c. Histogram showing the number of lncRNAs whose expression is up-regulated and down-regulated at each time point after the temperature change. d. Venn diagram showing the overlap in differentially expressed lncRNA at each time point

Subsequently, co-regulation was investigated for NATs and lincRNAs (Additional file 4: Table S4). In the case of NATs we compared their expression level to the expression of their corresponding sense transcript while for the lincRNAs to their two direct flanking protein coding genes. In general, we observed a positive Pearson correlation between sense and antisense transcript and we found the expression of the sense transcript was either not changed or changed in the same direction as the NAT (Additional file 4: Table S4). We also observed correlation (with prevalence of positive correlation) for most of the lincRNAs with at least one of the flanking genes, suggesting co-regulation between the lncRNAs and the protein coding genes located nearby (Additional file 4: Table S4). For example, AtLnc1444 showed a positive correlation with both the two direct flanking protein coding genes, C-REPEAT/DRE BINDING FACTOR 1 (CBF1) and CBF3, genes involved in response to low temperature [4851], while AtLnc488 showed a negative correlation with the transcription factor TEOSINTE BRANCHED1/CYCLOIDEA/PCF 15 (TCP15) [52, 53].

AtLnc428 plays a role in temperature-mediated flowering

We selected five ambient temperature responsive lncRNA for further functional characterization. We focused on lncRNAs for which a T-DNA insertion line is available in the predicted transcript or promoter region. We hypothesized that the promoter T-DNA insertions will affect the lncRNA expression level or pattern. We investigated the role of the temperature responsive AtLnc2, AtLnc120, AtLnc213, AtLnc428 and AtLnc1524 in the regulation of temperature-mediated flowering. Expression of AtLnc2 and AtLnc120 is up-regulated upon temperature increase while expression of AtLnc213, AtLnc428 and AtLnc1524 is down-regulated (Additional file 3: Table S3).

Wild-type and T-DNA insertion plants were grown for 3 weeks at 16 °C in long day conditions, and were phenotyped for their flowering time response after the switch in growth temperature to 25 °C and in comparison to plants that were maintained at 16 °C. Flowering time was quantified by counting rosette leaf number (RLN) and the number of days from sowing to the appearance of the main inflorescence (days after sowing; DAS). We measured the effect of temperature on flowering time as the ratio of RLN or DAS between mutant and Col-0 plants growing at 16 °C and 25 °C. Genotypes with complete temperature insensitivity will show a ratio of 1 [54]. Among the five lines tested, only the line carrying a T-DNA insertion in the AtLnc428 locus showed a significant change in temperature-mediated flowering compared to the wild-type (Additional file 5: Figure S1 and Fig. 2). We found AtLnc428 mutant plants to be significantly less sensitive to temperature-mediated flowering (Fig. 2a-c) and, for this reason, renamed this lncRNA ‘FLOWERING LONG INTERGENIC NON CODING RNA (FLINC)’.
Fig. 2
Fig. 2

FLINC plays a role in temperature-mediated flowering. a. Flowering time measured as rosette leaf number (RLN) for FLINC wild-type and mutant plants growing at 16 °C and 25 °C in long day conditions. The experiment was performed using four biological replicates with 13 plants per replicate for each genotype/condition. b. Flowering time measured as days after sowing (DAS) for FLINC wild-type (WT) and mutant plants growing at 16 °C and 25 °C in long day conditions. The experiment was performed using four biological replicates with 13 plants per replicate for each genotype/condition. c. Eight weeks-old FLINC WT and mutant plants grown at 16 °C and 25 °C in long day conditions. d. Flowering time in RLN for four independent T2 lines overexpressing FLINC (FLINC-OE) at 16 °C (blue) or 25 °C (red) in long day conditions. Box-plot showing the distribution of ca. 40 plants per replicate for each genotype/condition. FLINC-OE plants are more sensitive to temperature than the WT as shown by the higher RLN ratio. *** indicates significant differences at p-value < 0.0001, ** indicates significant differences at p-value < 0.001, and * indicates significant differences at p-value < 0.05 according to the Student t-test

FLINC is an intergenic lncRNA located in chromosome 1 between At1g56233, a gene coding for a defensin-like (DEFL) family protein and At1g56240, a gene coding for a phloem protein 2-B13. At1g56233 and At1g56240 are located respectively at 1.5 kb and 3.2 kb from FLINC (Additional file 6: Figure S2A). The effect of T-DNA insertion in FLINC on expression of FLINC was examined by qPCR by comparing wild-type and mutant plants in a segregating population (Additional file 6: Figure S2B). No FLINC expression was observed in plants with the T-DNA insertion, suggesting that the flinc mutant is a full knock-out. We also examined the effect of the FLINC T-DNA insertion on the expression of the genes flanking FLINC in all above ground tissues of 2 weeks-old plants in the vegetative stage of development. At1g56233 expression could not be detected in either wild type Col-0 or flinc plants, and we did not observe any significant change in the expression of At1g56240 (Additional file 6: Figure S2C). This result indicates that the T-DNA insertion does not affect the expression of the genes flanking FLINC in the investigated material.

Our results indicate a role for FLINC in the regulation of ambient temperature-mediated flowering. To confirm this observation, we generated plants overexpressing the lncRNA under control of the constitutive CaMV35S promoter and analyzed their flowering behavior at different ambient temperature conditions (Fig. 2d). In contrast to the flinc T-DNA insertion line, we found that plants overexpressing the lncRNA transcript are more sensitive to a temperature change than wild-type plants, and providing additional evidence for the role of FLINC in the control of ambient temperature-mediated flowering.

In our RNA-seq experiment, FLINC was down-regulated upon the temperature change under short day conditions (Additional file 3: Table S3). Because flowering time investigations were performed under long day conditions, we examined the effect of an increased ambient temperature on FLINC expression under these conditions. This analysis revealed a similar response upon an increase of temperature from 16 °C to 25 °C as was observed in the initial experiment under short day conditions (Additional file 6: Figure S2E), suggesting that regulation of FLINC expression by temperature is not day-length dependent.

Since a number of genes involved in the ambient temperature pathway are circadian-regulated or code for circadian clock components [55], we monitored the expression of FLINC in a 24 h-time course. We found that FLINC expression is not influenced by the circadian rhythm (Additional file 6: Figure S2F).

We also examined the FLINC expression profile during Arabidopsis development and in different Arabidopsis tissues (Additional file 6: Figure S2G and Figure S2H). We found that FLINC is expressed broadly in the plant, suggesting that FLINC plays a more extensive role during development than just the control of flowering.

To shed some light on the molecular basis of the flinc phenotype, we investigated the expression of selected flowering-related genes in wild type and 3 weeks-old flinc plants grown in long day conditions at 16 °C. These include the flowering repressors FLC and FLM, the florigen FLOWERING LOCUS T (FT) and the floral integrators APETALA 1 (AP1), SUPRESSOR OF COSTANS 1 (SOC1), LEAFY (LFY) and TERMINAL FLOWER 1 (TFL1) (Additional file 7: Figure S3A). We detected higher expression of FT in flinc than in the Col-0 wild type, in agreement with the early flowering phenotype of flinc (Fig. 2a).

LncRNAs have been show to exert their effect in numerous ways, including nucleotide pairing with mRNAs [1] of protein-coding genes or miRNAs. Therefore, we searched for regions homologous to FLINC within the Arabidopsis genome. We found ten regions with partial sequence homology to FLINC, of which five are located within protein coding genes. However, for none of these five genes differential expression was found in flinc in comparison to Col-0 wild type (Additional file 8: Table S5).

To obtain insight in conservation of FLINC, we searched for FLINC-related lncRNAs in genomes of other plant species using Blast and found various similar genes (e-value < 0.001), e.g. in the close Arabidopsis relatives Arabis alpina and Brassica rapa, but also in less-related species such as Solanum pennelli and Vitis vinifera (Additional file 7: Figure S3C). These results point to positive selection on FLINC. Further research is needed to confirm the proposed role of FLINC in ambient temperature-mediated flowering time control in Arabidopsis and to explore the function of the closely related genes in other species.


The discovery of lncRNAs as regulatory molecules revealed a new layer of complexity in gene expression regulation. So far, only few of the thousands of identified lncRNAs in plants have been functionally characterized. Nevertheless, these studies revealed that lncRNAs play a role in many different biological processes, from lateral root development to photomorphogenesis [1, 2]. We identified new Arabidopsis lncRNAs in SAM-enriched tissue and show ambient temperature-sensitive expression for a small subset of these lncRNAs. Furthermore, we showed that one of these lncRNAs, FLINC, regulates temperature-mediated flowering in Arabidopsis. Previous studies showed that expression of the lncRNAs transcribed from the FLC locus, COOLAIR [9], COLDAIR [10], ASL [12] and COLDWRAP [11], is regulated by a prolonged period of cold (4 °C). By contrast, FLINC expression is affected rapidly by a milder change of temperature. FLINC is an intergenic lncRNA and at this moment it is not clear which genes are regulated by FLINC. Thus, further research is needed to elucidate FLINC’s molecular mode of action. The observation that FT is higher expressed at low temperature in flinc mutant plants, suggests that FLINC acts prior to the integration of the environmental signals by FT. Next, the two transcription factors with a known role in temperature-mediated flowering time control, SOC1 [56] and AGAMOUS-LIKE 15 (AGL15) [57], bind to the FLINC locus suggesting a regulatory interaction. Moreover, the presence of sequence conservation in different plant species points to positive selection and a potential conserved role for FLINC in ambient temperature-mediated flowering time control.

Advances in RNA sequencing technology allow a more precise identification of lncRNA transcripts. Thanks to the use of strand specific RNA-seq data, it is now possible to identify transcripts derived from the opposite stand and therefore, get insight into expression profiles of NATs. The current sequence depth that can be obtained facilitates identification of low abundant transcripts, including as numerous lncRNAs. LncRNA expression appeared to be heavily responsive to environmental conditions, including temperature (this study) or drought [58], suggesting a role for these molecules in plant adaptation to the environment. The tissue-specific and environmental regulation of lncRNAs illustrates the need for defined experiments to identify all lncRNA-encoding genes in the genome.


The discovery of lncRNA as regulatory molecules revealed a new layer of complexity in the regulation of gene expression, and functional studies in plants revealed roles for lncRNAs in various biological processes. Here we show that lncRNA expression is influenced by ambient temperature changes and identified a significant role in ambient temperature-mediated flowering time the FLINC lncRNA. The challenge for the future is to decipher the biological function and molecular mode of action of ambient temperature responsive lncRNAs.


RNA-seq experiments

RNA-seq data were generated previously by Pajoro et al., 2017 [44]. In brief, plants were grown under short day conditions (8 h light, 16 h dark) on rock-wool in growth cabinets with LED lamps with light intensity of 200 μmol m-2 s-1 and 75% relative humidity at 16 °C for 5 weeks. Nutrients were supplied by sub-irrigation with Hyponex. After 5 weeks, plants were either left at the same temperature or moved to 25 °C (6 h after lights on). Three biological samples were generated for each condition. For each sample, tissue from ca. 10 plants was collected 4 hours after light on. Using jeweler’s forceps, leaves were removed to obtain SAM-enriched tissue. Total RNA was extracted using the Invitek InviTrap Spin Plant RNA Mini Kit (REF: 1064100300) according to the manufacturer’s protocol. DNase treatment was performed to remove genomic DNA. DNase I digestion was performed on total RNA using Turbo DNase from Ambion according to the manufacturer’s protocol. RNA integrity was checked on a by agarose gel electrophoresis after DNase I treatment. Samples were prepared for Illumina sequencing using the Illumina TruSeq Stranded mRNA Sample Prep kit (REF: 15032613) according to the manufacturer’s protocol, which generates polyadenylated transcripts. Libraries were analyzed on the Bioanalyzer and quantified using a Qubit fluorometer before pooling for sequencing on an Illumina HiSeq2500. Two lanes were used on a 125 bp paired-end (PE) run.

RNA-seq data analysis

RNA-seq reads were mapped against the Arabidopsis genome version TAIR10 ( using TopHat2 [59]. To identify un-annotated transcripts, reference based full length transcript reconstruction was performed separately for each sample using Cufflinks with strand-specific awareness [32]. Cuffmerge, which is part of the cufflinks package, was finally used for merging the individual cufflinks results into an overall set of full-length transcripts. The function “Cuffcompared” was used to compare the newly annotate transcripts with the TAIR10 reference annotation. We selected transcripts that belong to the class code “x”, “i” and “u”. The class code “x” defines “Exonic overlap with reference on the opposite strand”, to this class belong e.g. the Natural Antisense Transcripts (NAT). The class code “i” defines “transcribed genomic fragments falling entirely within a reference intron”, to this class belong the intronic lncRNAs (iLncRNA). Finally, the class code “u” defines the “unknown intergenic transcript”, to this class belong e.g. the long intergenic lncRNA (LincRNAs). To retrieve differentially expressed transcripts upon the temperature changes, the number of fragments mapping to the newly annotated transcripts was determined using HTseq count [60]. Differentially expressed genes were detected using DESeq2 [61].

Flowering time assay

Seeds of T-DNA insertion lines (Table 2), were obtained from the Nottingham Arabidopsis Stock Center (NASC). Plants were genotyped using the primers listed in Table 2. The progeny of mutant and wild-type plants from the segregating population was used in the flowering time experiments for AtLnc213 and AtLnc428. Since AtLnc2, AtLnc120 and AtLnc1524 homozygous T-DNA insertion lines were obtained from the Stock Center, Col-8 wild-type plants were used as comparison. For flowering time analyses, plants were randomly arranged in trays and grown in growth cabinets, as described above, but in long day-conditions (16 h light, 8 h dark). After 3 weeks of growth at 16 °C, half of the plants were moved to the 25 °C growth cabinet with the same day length, light intensity and humidity conditions. Flowering time was quantified by determining the time until the macroscopic appearance of the first flower bud (days after sowing, DAS; screening was done every day), and by counting rosette leaf numbers (RLN) per plant after bolting of all plants in the tray. Four biological replicates were used, with 13 plants per replicate for each genotype/condition. Data are shown as the average of ratios between the replicates. Student’s t-tests were calculated with GraphPad QuickCalcs.
Table 2

lncRNAs T-DNA insertion lines


T-DNA line

primer LP

primer RP





















Expression profile by qRT-PCR

Plant were grown under long day conditions (16 h light, 8 h dark) at 21 °C as described above. Three biological samples were generated for each condition tested, unless stated otherwise. For each sample, tissue from ca. 10 plants was collected using jeweler’s forceps. Total RNA was extracted and DNAse-treated as described above. cDNA was synthetized using the iScript kit from BIORAD (Cat. No. 1708890) using one μg of total RNA. Gene expression was measured by quantitative RT-PCR (qPCR) using the iQ™ SYBR® Green Supermix from BIORAD (Cat. No. 1708885). The qPCR primers are listed in Additional file 9: Table S6. Relative expression was calculated using TIP41 as reference gene for normalization and the formula 2-(ΔCttarget – ΔCtTIP41). Normalized fold expression was calculated by setting the wild-type (WT) value to one.

Generation of FLINC overexpression lines

Gateway cloning was used to generate the FLINC overexpression construct. The FLINC transcript was amplified by PCR from cDNA and cloned into the pDONR201 entry vector using a BP reaction. The entry vector was then recombined with the destination vector pB7FWG2 carrying the CaMV35S promoter. The resulting expression vector with the p35S:FLINC construct was transformed into Arabidopsis Col-0 plants by floral dip. Transformed plants were selected on media with 10 μg/ml of phosphinothricin.

Identification of protein coding genes with sequences homologous to FLINC

We used the BLASTN software [62] with default parameter to identify Arabidopsis protein coding genes with sequence similarity to FLINC. The genome coordinates that matched the FLINC sequence were compared to the coordinates of the Arabidopsis gene models using the Bedtools intersect software [63] and the gene coordinates as reference and –F 1 parameter. To identify sequences similar to FLINC in other plant species, we used Blastn software (NCBI) [64] with the parameter optimized for ‘somewhat similar sequences’. Only species with an e-value < 0.001 were considered.











days after sowing














intronic lncRNAs




long intergenic RNAs


Long non-coding RNAs




natural antisense transcripts


Next Generation Sequence.




Plant long non-coding RNA database




rosette leaf number









The authors would like to thank Kim Boutilier for critical proofreading of the manuscript.

Availability of data and material

The datasets analysed during the current study are available in the Gene Expression Omnibus (GEO) database, Series GSE85281.


This work was supported by an NWO grant (849.13.005) for the project ERACAPS13.012, FLOWPLAST. The RNA-seq experiments were supported by a ZonMw Enabling Technology Hotel project (40–43500–98-054). Funding body did not play any roles in the design of the study and collection, analysis, or interpretation of data neither in writing the manuscript.

Authors’ contributions

PA designed the study, performed the experiments, analyzed the data and wrote the paper. ES and LF analized the data. SJ, VF, JB and MB preformed experiments. YZ and FB generated material. GCA and RGHI conceived the study and revised the manuscript. All the authors read and approved the final manuscript.

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Authors’ Affiliations

Max Planck Institute for Plant Breeding Research, 50829 Köln, Germany
Laboratory of Phytopathology, Wageningen University and Research, 6708PB Wageningen, The Netherlands
Laboratory of Molecular Biology, Wageningen University and Research, 6708PB, Wageningen, The Netherlands
Bioscience, Wageningen University and Research, 6708PB Wageningen, The Netherlands


  1. Ariel F, Romero-Barrios N, Jegu T, Benhamed M, Crespi M. Battles and hijacks: noncoding transcription in plants. Trends Plant Sci. 2015;20(6):362–71.View ArticlePubMedGoogle Scholar
  2. Chekanova JA. Long non-coding RNAs and their functions in plants. Curr Opin Plant Biol. 2015;27:207–16.View ArticlePubMedGoogle Scholar
  3. Liu X, Hao LL, Li DY, Zhu LH. Hu SN: long non-coding RNAs and their biological roles in plants. Genom Proteom Bioinf. 2015;13(3):137–47.View ArticleGoogle Scholar
  4. Kopp F, Mendell JT. Functional classification and experimental dissection of long noncoding RNAs. Cell. 2018;172(3):393–407.View ArticlePubMedGoogle Scholar
  5. Ariel F, Jegu T, Latrasse D, Romero-Barrios N, Christ A, Benhamed M, Crespi M. Noncoding transcription by alternative RNA polymerases dynamically regulates an auxin-driven chromatin loop. Mol Cell. 2014;55(3):383–96.View ArticlePubMedGoogle Scholar
  6. Seo JS, Sun HX, Park BS, Huang CH, Yeh SD, Jung C, Chua NH. ELF18-INDUCED LONG-NONCODING RNA associates with mediator to enhance expression of innate immune response genes in Arabidopsis. Plant Cell. 2017;29(5):1024–38.View ArticlePubMedPubMed CentralGoogle Scholar
  7. Rohrig H, Schmidt J, Miklashevichs E, Schell J, John M. Soybean ENOD40 encodes two peptides that bind to sucrose synthase. Proc Natl Acad Sci U S A. 2002;99(4):1915–20.View ArticlePubMedPubMed CentralGoogle Scholar
  8. Franco-Zorrilla JM, Valli A, Todesco M, Mateos I, Puga MI, Rubio-Somoza I, Leyva A, Weigel D, Garcia JA, Paz-Ares J. Target mimicry provides a new mechanism for regulation of microRNA activity. Nat Genet. 2007;39(8):1033–7.View ArticlePubMedGoogle Scholar
  9. Swiezewski S, Liu FQ, Magusin A, Dean C. Cold-induced silencing by long antisense transcripts of an Arabidopsis Polycomb target. Nature. 2009;462(7274):799–U122.View ArticlePubMedGoogle Scholar
  10. Heo JB, Sung S. Vernalization-mediated epigenetic silencing by a long Intronic noncoding RNA. Science. 2011;331(6013):76–9.View ArticlePubMedGoogle Scholar
  11. Kim DH, Sung S. Vernalization-triggered intragenic chromatin loop formation by long noncoding RNAs. Dev Cell. 2017;40(3):302–12. e304View ArticlePubMedPubMed CentralGoogle Scholar
  12. Shin JH, Chekanova JA. Arabidopsis RRP6L1 and RRP6L2 Function in FLOWERING LOCUS C Silencing via Regulation of Antisense RNA Synthesis. PLoS Genet. 2014;10(9):e1004612.View ArticlePubMedPubMed CentralGoogle Scholar
  13. Lee JH, Ryu HS, Chung KS, Pose D, Kim S, Schmid M, Ahn JH. Regulation of temperature-responsive flowering by MADS-box transcription factor repressors. Science. 2013;342(6158):628–32.View ArticlePubMedGoogle Scholar
  14. Pose D, Verhage L, Ott F, Yant L, Mathieu J, Angenent GC, Immink RGH, Schmid M. Temperature-dependent regulation of flowering by antagonistic FLM variants. Nature. 2013;503(7476):414.View ArticlePubMedGoogle Scholar
  15. Airoldi CA, Mckay M, Davies B. MAF2 is regulated by temperature-dependent splicing and represses flowering at low temperatures in parallel with FLM. PLoS One. 2015;10(5):e0126516.View ArticlePubMedPubMed CentralGoogle Scholar
  16. Gu XF, Le C, Wang YZ, Li ZC, Jiang DH, Wang YQ, He YH. Arabidopsis FLC clade members form flowering-repressor complexes coordinating responses to endogenous and environmental cues. Nat Commun. 2013;4:1947.View ArticlePubMedPubMed CentralGoogle Scholar
  17. Verhage L, Angenent GC, Immink RG. Research on floral timing by ambient temperature comes into blossom. Trends Plant Sci. 2014;19(9):583–91.View ArticlePubMedGoogle Scholar
  18. Quint M, Delker C, Franklin KA, Wigge PA, Halliday KJ, van Zanten M. Molecular and genetic control of plant thermomorphogenesis. Nature plants. 2016;2:15190.View ArticlePubMedGoogle Scholar
  19. Wigge PA. Ambient temperature signalling in plants. Curr Opin Plant Biol. 2013;16(5):661–6.View ArticlePubMedGoogle Scholar
  20. Fernandez V, Takahashi Y, Le Gourrierec J, Coupland G. Photoperiodic and thermosensory pathways interact through CONSTANS to promote flowering at high temperature under short days. Plant J. 2016;86(5):426–40.View ArticlePubMedGoogle Scholar
  21. Galvao VC, Collani S, Horrer D, Schmid M. Gibberellic acid signaling is required for ambient temperature-mediated induction of flowering in Arabidopsis thaliana. Plant J. 2015;84(5):949–62.View ArticlePubMedGoogle Scholar
  22. Kumar SV, Lucyshyn D, Jaeger KE, Alos E, Alvey E, Harberd NP, Wigge PA. Transcription factor PIF4 controls the thermosensory activation of flowering. Nature. 2012;484(7393):242–U127.View ArticlePubMedPubMed CentralGoogle Scholar
  23. Press MO, Lanctot A, Queitsch C. PIF4 and ELF3 act independently in Arabidopsis thaliana Thermoresponsive flowering. PLoS One. 2016;11(8):e0161791.View ArticlePubMedPubMed CentralGoogle Scholar
  24. Marin-Gonzalez E, Matias-Hernandez L, Aguilar-Jaramillo AE, Lee JH, Ahn JH, Suarez-Lopez P, Pelaz S. SHORT VEGETATIVE PHASE Up-Regulates TEMPRANILLO2 floral repressor at low ambient temperatures. Plant Physiol. 2015;169(2):1214–24.View ArticlePubMedPubMed CentralGoogle Scholar
  25. Yan YY, Shen LS, Chen Y, Bao SJ, Thong ZH, Yu H. A MYB-domain protein EFM mediates flowering responses to environmental cues in Arabidopsis. Dev Cell. 2014;30(4):437–48.View ArticlePubMedGoogle Scholar
  26. Delker C, van Zanten M, Quint M, Enlightened T. Trends Plant Sci. 2017;22(3):185–7.View ArticlePubMedGoogle Scholar
  27. Jung JH, Domijan M, Klose C, Biswas S, Ezer D, Gao MJ, Khattak AK, Box MS, Charoensawan V, Cortijo S, et al. Phytochromes function as thermosensors in Arabidopsis. Science. 2016;354(6314):886–9.View ArticlePubMedGoogle Scholar
  28. Legris M, Klose C, Burgie ES, Rojas CC, Neme M, Hiltbrunner A, Wigge PA, Schafer E, Vierstra RD, Casal JJ. Phytochrome B integrates light and temperature signals in Arabidopsis. Science. 2016;354(6314):897–900.View ArticlePubMedGoogle Scholar
  29. Box MS, Huang BE, Domijan M, Jaeger KE, Khattak AK, Yoo SJ, Sedivy EL, Jones DM, Hearn TJ, AAR W, et al. ELF3 controls Thermoresponsive growth in Arabidopsis. Curr Biol. 2015;25(2):194–9.View ArticlePubMedGoogle Scholar
  30. Ezer D, Jung JH, Lan H, Biswas S, Gregoire L, Box MS, Charoensawan V, Cortijo S, Lai XL, Stockle D, et al. The evening complex coordinates environmental and endogenous signals in Arabidopsis. Nature plants. 2017;3(7):17087.View ArticlePubMedPubMed CentralGoogle Scholar
  31. Seo PJ, Park MJ, Lim MH, Kim SG, Lee M, Baldwin IT, Park CM. A self-regulatory circuit of CIRCADIAN CLOCK-ASSOCIATED1 underlies the circadian clock regulation of temperature responses in Arabidopsis. Plant Cell. 2012;24(6):2427–42.View ArticlePubMedPubMed CentralGoogle Scholar
  32. Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, Salzberg SL, Wold BJ, Pachter L. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol. 2010;28(5):511–U174.View ArticlePubMedPubMed CentralGoogle Scholar
  33. Landrein B, Refahi Y, Besnard F, Hervieux N, Mirabet V, Boudaoud A, Vernoux T, Hamant O. Meristem size contributes to the robustness of phyllotaxis in Arabidopsis. J Exp Bot. 2015;66(5):1317–24.View ArticlePubMedGoogle Scholar
  34. Scarano D, Rao R, Corrado G. In Silico identification and annotation of noncoding RNAs by RNA-seq and De Novo assembly of the transcriptome of Tomato Fruits. PLoS One. 2017;12(2):e0171504.View ArticlePubMedPubMed CentralGoogle Scholar
  35. Jin JJ, Liu J, Wang H, Wong L, Chua NH. PLncDB: plant long non-coding RNA database. Bioinformatics. 2013;29(8):1068–71.View ArticlePubMedPubMed CentralGoogle Scholar
  36. Szczesniak MW, Rosikiewicz W, Makalowska I. CANTATAdb: a collection of plant long non-coding RNAs. Plant Cell Physiol. 2016;57(1):e8.View ArticlePubMedGoogle Scholar
  37. Krishnakumar V, Hanlon MR, Contrino S, Ferlanti ES, Karamycheva S, Kim M, Rosen BD, Cheng CY, Moreira W, Mock SA, et al. Araport: the Arabidopsis information portal. Nucleic Acids Res. 2015;43(D1):D1003–9.View ArticlePubMedGoogle Scholar
  38. Bouche F, Lobet G, Tocquin P, Perilleux C. FLOR-ID: an interactive database of flowering-time gene networks in Arabidopsis thaliana. Nucleic Acids Res. 2016;44(D1):D1167–71.View ArticlePubMedGoogle Scholar
  39. Ratcliffe OJ, Kumimoto RW, Wong BJ, Riechmann JL. Analysis of the Arabidopsis MADS AFFECTING FLOWERING gene family: MAF2 prevents vernalization by short periods of cold. Plant Cell. 2003;15(5):1159–69.View ArticlePubMedPubMed CentralGoogle Scholar
  40. Rosloski SM, Singh A, Jali SS, Balasubramanian S, Weigel D, Grbic V. Functional analysis of splice variant expression of MADS AFFECTING FLOWERING 2 of Arabidopsis thaliana. Plant Mol Biol. 2013;81(1–2):57–69.View ArticlePubMedGoogle Scholar
  41. Song YH, Estrada DA, Johnson RS, Kim SK, Lee SY, MacCoss MJ, Imaizumi T. Distinct roles of FKF1, GIGANTEA, and ZEITLUPE proteins in the regulation of CONSTANS stability in Arabidopsis photoperiodic flowering. Proc Natl Acad Sci U S A. 2014;111(49):17672–7.View ArticlePubMedPubMed CentralGoogle Scholar
  42. Song YH, Smith RW, BJ T, Millar AJ, Imaizumi T. FKF1 conveys timing information for CONSTANS stabilization in photoperiodic flowering. Science. 2012;336(6084):1045–9.View ArticlePubMedPubMed CentralGoogle Scholar
  43. Kim SY, He YH, Jacob Y, Noh YS, Michaels S, Amasino R. Establishment of the vernalization-responsive, winter-annual habit in Arabidopsis requires a putative histone H3 methyl transferase. Plant Cell. 2005;17(12):3301–10.View ArticlePubMedPubMed CentralGoogle Scholar
  44. Pajoro A, Severing E, Angenent GC, Immink RGH. Histone H3 lysine 36 methylation affects temperature-induced alternative splicing and flowering in plants. Genome Biol. 2017;18(1):102.View ArticlePubMedPubMed CentralGoogle Scholar
  45. Sotelo-Silveira M, Cucinotta M, Chauvin AL, Montes RAC, Colombo L, Marsch-Martinez N, de Folter S. Cytochrome P450 CYP78A9 is involved in Arabidopsis reproductive development. Plant Physiol. 2013;162(2):779–99.View ArticlePubMedPubMed CentralGoogle Scholar
  46. Kolesinski P, Golik P, Grudnik P, Piechota J, Markiewicz M, Tarnawski M, Dubin G, Szczepaniak A. Insights into eukaryotic Rubisco assembly - crystal structures of RbcX chaperones from Arabidopsis thaliana. Bba-Gen Subjects. 2013;1830(4):2899–906.View ArticleGoogle Scholar
  47. Kolesinski P, Piechota J, Szczepaniak A. Initial characteristics of RbcX proteins from Arabidopsis thaliana. Plant Mol Biol. 2011;77(4–5):447–59.View ArticlePubMedPubMed CentralGoogle Scholar
  48. Cook D, Fowler S, Fiehn O, Thomashow MF. A prominent role for the CBF cold response pathway in configuring the low-temperature metabolome of Arabidopsis. Proc Natl Acad Sci U S A. 2004;101(42):15243–8.View ArticlePubMedPubMed CentralGoogle Scholar
  49. Lee CM, Thomashow MF. Photoperiodic regulation of the C-repeat binding factor (CBF) cold acclimation pathway and freezing tolerance in Arabidopsis thaliana. Proc Natl Acad Sci U S A. 2012;109(37):15054–9.View ArticlePubMedPubMed CentralGoogle Scholar
  50. Medina J, Catala R, Salinas J. The CBFs: three arabidopsis transcription factors to cold acclimate. Plant Sci. 2011;180(1):3–11.View ArticlePubMedGoogle Scholar
  51. Zhou MQ, Chen H, Wei DH, Ma H, Lin J. Arabidopsis CBF3 and DELLAs positively regulate each other in response to low temperature. Sci Rep-Uk. 2017;7:39819.View ArticleGoogle Scholar
  52. Daviere JM, Wild M, Regnault T, Baumberger N, Eisler H, Genschik P, Achard P. Class I TCP-DELLA interactions in inflorescence shoot apex determine plant height. Curr Biol. 2014;24(16):1923–8.View ArticlePubMedGoogle Scholar
  53. Kieffer M, Master V, Waites R, Davies B. TCP14 and TCP15 affect internode length and leaf shape in Arabidopsis. Plant J. 2011;68(1):147–58.View ArticlePubMedPubMed CentralGoogle Scholar
  54. Pajoro A, Verhage L, Immink RGH. Plasticity versus adaptation of ambient-temperature flowering response. Trends Plant Sci. 2016;21(1):6–8.View ArticlePubMedGoogle Scholar
  55. Capovilla G, Pajoro A, Immink RG, Schmid M. Role of alternative pre-mRNA splicing in temperature signaling. Curr Opin Plant Biol. 2015;27:97–103.View ArticlePubMedGoogle Scholar
  56. Immink RGH, Pose D, Ferrario S, Ott F, Kaufmann K, Valentim FL, de Folter S, van der Wal F, van Dijk ADJ, Schmid M, et al. Characterization of SOC1's central role in flowering by the identification of its upstream and downstream regulators. Plant Physiol. 2012;160(1):433–49.View ArticlePubMedPubMed CentralGoogle Scholar
  57. Zheng YM, Ren N, Wang H, Stromberg AJ, Perry SE. Global identification of targets of the Arabidopsis MADS domain protein AGAMOUS-Like15. Plant Cell. 2009;21(9):2563–77.View ArticlePubMedPubMed CentralGoogle Scholar
  58. Liu J, Jung C, Xu J, Wang H, Deng SL, Bernad L, Arenas-Huertero C, Chua NH. Genome-wide analysis uncovers regulation of long intergenic noncoding RNAs in Arabidopsis. Plant Cell. 2012;24(11):4333–45.View ArticlePubMedPubMed CentralGoogle Scholar
  59. Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10(3):R25.View ArticlePubMedPubMed CentralGoogle Scholar
  60. Anders S, Pyl PT, Huber W. HTSeq-a python framework to work with high-throughput sequencing data. Bioinformatics. 2015;31(2):166–9.View ArticlePubMedGoogle Scholar
  61. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.View ArticlePubMedPubMed CentralGoogle Scholar
  62. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215(3):403–10.View ArticlePubMedGoogle Scholar
  63. Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010;26(6):841–2.View ArticlePubMedPubMed CentralGoogle Scholar
  64. Agarwala R, Barrett T, Beck J, Benson DA, Bollin C, Bolton E, Bourexis D, Brister JR, Bryant SH, Lanese K, et al. Database resources of the National Center for biotechnology information. Nucleic Acids Res. 2016;44(D1):D7–D19.View ArticleGoogle Scholar


© The Author(s). 2018