Root branching under high salinity requires auxin-independent modulation of LATERAL ORGAN BOUNDARY DOMAIN 16 function.

Salinity stress constrains lateral root (LR) growth and severely affects plant growth. Auxin signaling regulates LR formation, but the molecular mechanism by which salinity affects root auxin signaling and whether salt induces other pathways that regulate LR development remains unknown. In Arabidopsis thaliana, the auxin-regulated transcription factor LATERAL ORGAN BOUNDARY DOMAIN 16 (LBD16) is an essential player in LR development under control conditions. Here, we show that under high-salt conditions, an alternative pathway regulates LBD16 expression. Salt represses auxin signaling but, in parallel, activates ZINC FINGER OF ARABIDOPSIS THALIANA 6 (ZAT6), a transcriptional activator of LBD16. ZAT6 activates LBD16 expression, thus contributing to downstream cell wall remodeling and promoting LR development under high-salt conditions. Our study thus shows that the integration of auxin-dependent repressive and salt-activated auxin-independent pathways converging on LBD16 modulates root branching under high-salt conditions.

A catalogue of recombination coldspots in interspecific tomato hybrids.

Increasing natural resistance and resilience in plants is key for ensuring food security within a changing climate. Breeders improve these traits by crossing cultivars with their wild relatives and introgressing specific alleles through meiotic recombination. However, some genomic regions are devoid of recombination especially in crosses between divergent genomes, limiting the combinations of desirable alleles. Here, we used pooled-pollen sequencing to build a map of recombinant and non-recombinant regions between tomato and five wild relatives commonly used for introgressive tomato breeding. We detected hybrid-specific recombination coldspots that underscore the role of structural variations in modifying recombination patterns and maintaining genetic linkage in interspecific crosses. Crossover regions and coldspots show strong association with specific TE superfamilies exhibiting differentially accessible chromatin between somatic and meiotic cells. About two-thirds of the genome are conserved coldspots, located mostly in the pericentromeres and enriched with retrotransposons. The coldspots also harbor genes associated with agronomic traits and stress resistance, revealing undesired consequences of linkage drag and possible barriers to breeding. We presented examples of linkage drag that can potentially be resolved by pairing tomato with other wild species. Overall, this catalogue will help breeders better understand crossover localization and make informed decisions on generating new tomato varieties.

Protein interaction mapping reveals widespread targeting of development-related host transcription factors by phytoplasma effectors.

Phytoplasmas are pathogenic bacteria that reprogram plant host development for their own benefit. Previous studies have characterized a few different phytoplasma effector proteins that destabilize specific plant transcription factors. However, these are only a small fraction of the potential effectors used by phytoplasmas; therefore, the molecular mechanisms through which phytoplasmas modulate their hosts require further investigation. To obtain further insights into the phytoplasma infection mechanisms, we generated a protein-protein interaction network between a broad set of phytoplasma effectors and a large, unbiased collection of Arabidopsis thaliana transcription factors and transcriptional regulators. We found widespread, but specific, interactions between phytoplasma effectors and host transcription factors, especially those related to host developmental processes. In particular, many unrelated effectors target specific sets of TCP transcription factors, which regulate plant development and immunity. Comparison with other host-pathogen protein interaction networks shows that phytoplasma effectors have unusual targets, indicating that phytoplasmas have evolved a unique and unusual infection strategy. This study contributes a rich and solid data source that guides further investigations of the functions of individual effectors, as demonstrated for some herein. Moreover, the dataset provides insights into the underlying molecular mechanisms of phytoplasma infection.

Interploidy Introgression Shaped Adaptation during the Origin and Domestication History of Brassica napus.

Polyploidy is recurrent across the tree of life and known as an evolutionary driving force in plant diversification and crop domestication. How polyploid plants adapt to various habitats has been a fundamental question that remained largely unanswered. Brassica napus is a major crop cultivated worldwide, resulting from allopolyploidy between unknown accessions of diploid B. rapa and B. oleracea. Here, we used whole-genome resequencing data of accessions representing the majority of morphotypes and ecotypes from the species B. rapa, B. oleracea, and B. napus to investigate the role of polyploidy during domestication. To do so, we first reconstructed the phylogenetic history of B. napus, which supported the hypothesis that the emergence of B. napus derived from the hybridization of European turnip of B. rapa and wild B. oleracea. These analyses also showed that morphotypes of swede and Siberian kale (used as vegetable and fodder) were domesticated before rapeseed (oil crop). We next observed that frequent interploidy introgressions from sympatric diploids were prominent throughout the domestication history of B. napus. Introgressed genomic regions were shown to increase the overall genetic diversity and tend to be localized in regions of high recombination. We detected numerous candidate adaptive introgressed regions and found evidence that some of the genes in these regions contributed to phenotypic diversification and adaptation of different morphotypes. Overall, our results shed light on the origin and domestication of B. napus and demonstrate interploidy introgression as an important mechanism that fuels rapid diversification in polyploid species.

A gene regulatory network critical for axillary bud dormancy directly controlled by Arabidopsis BRANCHED1.

The Arabidopsis thaliana transcription factor BRANCHED1 (BRC1) plays a pivotal role in the control of shoot branching as it integrates environmental and endogenous signals that influence axillary bud growth. Despite its remarkable activity as a growth inhibitor, the mechanisms by which BRC1 promotes bud dormancy are largely unknown. We determined the genome-wide BRC1 binding sites in vivo and combined these with transcriptomic data and gene co-expression analyses to identify bona fide BRC1 direct targets. Next, we integrated multi-omics data to infer the BRC1 gene regulatory network (GRN) and used graph theory techniques to find network motifs that control the GRN dynamics. We generated an open online tool to interrogate this network. A group of BRC1 target genes encoding transcription factors (BTFs) orchestrate this intricate transcriptional network enriched in abscisic acid-related components. Promoter::ß-GLUCURONIDASE transgenic lines confirmed that BTFs are expressed in axillary buds. Transient co-expression assays and studies in planta using mutant lines validated the role of BTFs in modulating the GRN and promoting bud dormancy. This knowledge provides access to the developmental mechanisms that regulate shoot branching and helps identify candidate genes to use as tools to adapt plant architecture and crop production to ever-changing environmental conditions.

Individual lipid transfer proteins from Tanacetum parthenium show different specificity for extracellular accumulation of sesquiterpenes.

KEY MESSAGE: A highly specialized function for individual LTPs for different products from the same terpenoid biosynthesis pathway is described and the function of an LTP GPI anchor is studied. Sequiterpenes produced in glandular trichomes of the medicinal plant Tanacetum parthenium (feverfew) accumulate in the subcuticular extracellular space. Transport of these compounds over the plasma membrane is presumably by specialized membrane transporters, but it is still not clear how these hydrophobic compounds are subsequently transported over the hydrophilic cell wall. Here we identified eight so-called non-specific Lipid transfer proteins (nsLTPs) genes that are expressed in feverfew trichomes. A putative function of these eight nsLTPs in transport of the lipophilic sesquiterpene lactones produced in feverfew trichomes, was tested in an in-planta transport assay using transient expression in Nicotiana benthamiana. Of eight feverfew nsLTP candidate genes analyzed, two (TpLTP1 and TpLTP2) can specifically improve extracellular accumulation of the sesquiterpene costunolide, while one nsLTP (TpLTP3) shows high specificity towards export of parthenolide. The specificity of the nsLTPs was also tested in an assay that test for the exclusion capacity of the nsLTP for influx of extracellular substrates. In such assay, TpLTP3 was identified as most effective in blocking influx of both costunolide and parthenolide, when these substrates are infiltrated into the apoplast. The TpLTP3 is special in having a GPI-anchor domain, which is essential for the export activity of TpLTP3. However, addition of the TpLTP3 GPI-anchor domain to TpLTP1 resulted in loss of TpLTP1 export activity. These novel export and exclusion assays thus provide new means to test functionality of plant nsLTPs.

Domestication Shapes Recombination Patterns in Tomato.

Meiotic recombination is a biological process of key importance in breeding, to generate genetic diversity and develop novel or agronomically relevant haplotypes. In crop tomato, recombination is curtailed as manifested by linkage disequilibrium decay over a longer distance and reduced diversity compared with wild relatives. Here, we compared domesticated and wild populations of tomato and found an overall conserved recombination landscape, with local changes in effective recombination rate in specific genomic regions. We also studied the dynamics of recombination hotspots resulting from domestication and found that loss of such hotspots is associated with selective sweeps, most notably in the pericentromeric heterochromatin. We detected footprints of genetic changes and structural variants, among them associated with transposable elements, linked with hotspot divergence during domestication, likely causing fine-scale alterations to recombination patterns and resulting in linkage drag.

The tomato cytochrome P450 CYP712G1 catalyses the double oxidation of orobanchol en route to the rhizosphere signalling strigolactone, solanacol.

Strigolactones (SLs) are rhizosphere signalling molecules and phytohormones. The biosynthetic pathway of SLs in tomato has been partially elucidated, but the structural diversity in tomato SLs predicts that additional biosynthetic steps are required. Here, root RNA-seq data and co-expression analysis were used for SL biosynthetic gene discovery. This strategy resulted in a candidate gene list containing several cytochrome P450s. Heterologous expression in Nicotiana benthamiana and yeast showed that one of these, CYP712G1, can catalyse the double oxidation of orobanchol, resulting in the formation of three didehydro-orobanchol (DDH) isomers. Virus-induced gene silencing and heterologous expression in yeast showed that one of these DDH isomers is converted to solanacol, one of the most abundant SLs in tomato root exudate. Protein modelling and substrate docking analysis suggest that hydroxy-orbanchol is the likely intermediate in the conversion from orobanchol to the DDH isomers. Phylogenetic analysis demonstrated the occurrence of CYP712G1 homologues in the Eudicots only, which fits with the reports on DDH isomers in that clade. Protein modelling and orobanchol docking of the putative tobacco CYP712G1 homologue suggest that it can convert orobanchol to similar DDH isomers as tomato.

Integrating structure-based machine learning and co-evolution to investigate specificity in plant sesquiterpene synthases.

Sesquiterpene synthases (STSs) catalyze the formation of a large class of plant volatiles called sesquiterpenes. While thousands of putative STS sequences from diverse plant species are available, only a small number of them have been functionally characterized. Sequence identity-based screening for desired enzymes, often used in biotechnological applications, is difficult to apply here as STS sequence similarity is strongly affected by species. This calls for more sophisticated computational methods for functionality prediction. We investigate the specificity of precursor cation formation in these elusive enzymes. By inspecting multi-product STSs, we demonstrate that STSs have a strong selectivity towards one precursor cation. We use a machine learning approach combining sequence and structure information to accurately predict precursor cation specificity for STSs across all plant species. We combine this with a co-evolutionary analysis on the wealth of uncharacterized putative STS sequences, to pinpoint residues and distant functional contacts influencing cation formation and reaction pathway selection. These structural factors can be used to predict and engineer enzymes with specific functions, as we demonstrate by predicting and characterizing two novel STSs from Citrus bergamia.

Meiotic recombination profiling of interspecific hybrid F1 tomato pollen by linked read sequencing.

Genome wide screening of pooled pollen samples from a single interspecific F1 hybrid obtained from a cross between tomato, Solanum lycopersicum and its wild relative, Solanum pimpinellifolium using linked read sequencing of the haploid nuclei, allowed profiling of the crossover (CO) and gene conversion (GC) landscape. We observed a striking overlap between cold regions of CO in the male gametes and our previously established F6 recombinant inbred lines (RILs) population. COs were overrepresented in non-coding regions in the gene promoter and 5'UTR regions of genes. Poly-A/T and AT rich motifs were found enriched in 1 kb promoter regions flanking the CO sites. Non-crossover associated allelic and ectopic GCs were detected in most chromosomes, confirming that besides CO, GC represents also a source for genetic diversity and genome plasticity in tomato. Furthermore, we identified processed break junctions pointing at the involvement of both homology directed and non-homology directed repair pathways, suggesting a recombination machinery in tomato that is more complex than currently anticipated.

Geometricus represents protein structures as shape-mers derived from moment invariants.

MOTIVATION: As the number of experimentally solved protein structures rises, it becomes increasingly appealing to use structural information for predictive tasks involving proteins. Due to the large variation in protein sizes, folds and topologies, an attractive approach is to embed protein structures into fixed-length vectors, which can be used in machine learning algorithms aimed at predicting and understanding functional and physical properties. Many existing embedding approaches are alignment based, which is both time-consuming and ineffective for distantly related proteins. On the other hand, library- or model-based approaches depend on a small library of fragments or require the use of a trained model, both of which may not generalize well. RESULTS: We present Geometricus, a novel and universally applicable approach to embedding proteins in a fixed-dimensional space. The approach is fast, accurate, and interpretable. Geometricus uses a set of 3D moment invariants to discretize fragments of protein structures into shape-mers, which are then counted to describe the full structure as a vector of counts. We demonstrate the applicability of this approach in various tasks, ranging from fast structure similarity search, unsupervised clustering and structure classification across proteins from different superfamilies as well as within the same family. AVAILABILITY AND IMPLEMENTATION: Python code available at https://git.wur.nl/durai001/geometricus.

Coevolution-based prediction of protein-protein interactions in polyketide biosynthetic assembly lines.

MOTIVATION: Polyketide synthases (PKSs) are enzymes that generate diverse molecules of great pharmaceutical importance, including a range of clinically used antimicrobials and antitumor agents. Many polyketides are synthesized by cis-AT modular PKSs, which are organized in assembly lines, in which multiple enzymes line up in a specific order. This order is defined by specific protein-protein interactions (PPIs). The unique modular structure and catalyzing mechanism of these assembly lines makes their products predictable and also spurred combinatorial biosynthesis studies to produce novel polyketides using synthetic biology. However, predicting the interactions of PKSs, and thereby inferring the order of their assembly line, is still challenging, especially for cases in which this order is not reflected by the ordering of the PKS-encoding genes in the genome. RESULTS: Here, we introduce PKSpop, which uses a coevolution-based PPI algorithm to infer protein order in PKS assembly lines. Our method accurately predicts protein orders (93% accuracy). Additionally, we identify new residue pairs that are key in determining interaction specificity, and show that coevolution of N- and C-terminal docking domains of PKSs is significantly more predictive for PPIs than coevolution between ketosynthase and acyl carrier protein domains. AVAILABILITY AND IMPLEMENTATION: The code is available on http://www.bif.wur.nl/ (under 'Software'). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Comprehensive phenotyping reveals interactions and functions of Arabidopsis thaliana TCP genes in yield determination.

Members of the Arabidopsis thaliana TCP transcription factor (TF) family affect plant growth and development. We systematically quantified the effect of mutagenizing single or multiple TCP TFs and how altered vegetative growth or branching influences final seed yield. We monitored rosette growth over time and branching patterns and seed yield characteristics at the end of the lifecycle. Subsequently, an approach was developed to disentangle vegetative growth and to determine possible effects on seed yield. Analysis of growth parameters showed all investigated tcp mutants to be affected in certain growth aspects compared with wild-type plants, highlighting the importance of TCP TFs in plant development. Furthermore, we found evidence that all class II TCPs are involved in axillary branch outgrowth, either as inhibitors (BRANCHED-like genes) or enhancers (JAW- and TCP5-like genes). Comprehensive phenotyping of plants mutant for single or multiple TCP TFs reveals that the proposed opposite functions of class I and class II TCPs in plant growth needs revision and shows complex interactions between closely related TCP genes instead of full genetic redundancy. In various instances, the alterations in vegetative growth or in branching patterns result into negative trade-off effects on seed yield that were missed in previous studies, showing the importance of comprehensive and quantitative phenotyping.

Improved inference of intermolecular contacts through protein-protein interaction prediction using coevolutionary analysis.

MOTIVATION: Predicting residue-residue contacts between interacting proteins is an important problem in bioinformatics. The growing wealth of sequence data can be used to infer these contacts through correlated mutation analysis on multiple sequence alignments of interacting homologs of the proteins of interest. This requires correct identification of pairs of interacting proteins for many species, in order to avoid introducing noise (i.e. non-interacting sequences) in the analysis that will decrease predictive performance. RESULTS: We have designed Ouroboros, a novel algorithm to reduce such noise in intermolecular contact prediction. Our method iterates between weighting proteins according to how likely they are to interact based on the correlated mutations signal, and predicting correlated mutations based on the weighted sequence alignment. We show that this approach accurately discriminates between protein interaction versus non-interaction and simultaneously improves the prediction of intermolecular contact residues compared to a naive application of correlated mutation analysis. This requires no training labels concerning interactions or contacts. Furthermore, the method relaxes the assumption of one-to-one interaction of previous approaches, allowing for the study of many-to-many interactions. AVAILABILITY AND IMPLEMENTATION: Source code and test data are available at www.bif.wur.nl/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Novel routes towards bioplastics from plants: elucidation of the methylperillate biosynthesis pathway from Salvia dorisiana trichomes.

Plants produce a large variety of highly functionalized terpenoids. Functional groups such as partially unsaturated rings and carboxyl groups provide handles to use these compounds as feedstock for biobased commodity chemicals. For instance, methylperillate, a monoterpenoid found in Salvia dorisiana, may be used for this purpose, as it carries both an unsaturated ring and a methylated carboxyl group. The biosynthetic pathway of methylperillate in plants is still unclear. In this work, we identified glandular trichomes from S. dorisiana as the location of biosynthesis and storage of methylperillate. mRNA from purified trichomes was used to identify four genes that can encode the pathway from geranyl diphosphate towards methylperillate. This pathway includes a (-)-limonene synthase (SdLS), a limonene 7-hydroxylase (SdL7H, CYP71A76), and a perillyl alcohol dehydrogenase (SdPOHDH). We also identified a terpene acid methyltransferase, perillic acid O-methyltransferase (SdPAOMT), with homology to salicylic acid OMTs. Transient expression in Nicotiana benthamiana of these four genes, in combination with a geranyl diphosphate synthase to boost precursor formation, resulted in production of methylperillate. This demonstrates the potential of these enzymes for metabolic engineering of a feedstock for biobased commodity chemicals.

The santalene synthase from Cinnamomum camphora: Reconstruction of a sesquiterpene synthase from a monoterpene synthase.

Plant terpene synthases (TPSs) can mediate formation of a large variety of terpenes, and their diversification contributes to the specific chemical profiles of different plant species and chemotypes. Plant genomes often encode a number of related terpene synthases, which can produce very different terpenes. The relationship between TPS sequence and resulting terpene product is not completely understood. In this work we describe two TPSs from the Camphor tree Cinnamomum camphora (L.) Presl. One of these, CiCaMS, acts as a monoterpene synthase (monoTPS), and mediates the production of myrcene, while the other, CiCaSSy, acts as a sesquiterpene synthase (sesquiTPS), and catalyses the production of α-santalene, ß-santalene and trans-α-bergamotene. Interestingly, these enzymes share 97% DNA sequence identity and differ only in 22 amino acid residues out of 553. To understand which residues are essential for the catalysis of monoterpenes resp. sesquiterpenes, a number of hybrid synthases were prepared, and supplemented by a set of single-residue variants. These were tested for their ability to produce monoterpenes and sesquiterpenes by in vivo production of sesquiterpenes in E. coli, and by in vitro enzyme assays. This analysis pinpointed three residues in the sequence which could mediate the change in product specificity from a monoterpene synthase to a sesquiterpene synthase. Another set of three residues defined the sesquiterpene product profile, including the ratios between sesquiterpene products.

DNA sequence and shape are predictive for meiotic crossovers throughout the plant kingdom.

A better understanding of genomic features influencing the location of meiotic crossovers (COs) in plant species is both of fundamental importance and of practical relevance for plant breeding. Using CO positions with sufficiently high resolution from four plant species [Arabidopsis thaliana, Solanum lycopersicum (tomato), Zea mays (maize) and Oryza sativa (rice)] we have trained machine-learning models to predict the susceptibility to CO formation. Our results show that CO occurrence within various plant genomes can be predicted by DNA sequence and shape features. Several features related to genome content and to genomic accessibility were consistently either positively or negatively related to COs in all four species. Other features were found as predictive only in specific species. Gene annotation-related features were especially predictive for maize, whereas in tomato and Arabidopsis propeller twist and helical twist (DNA shape features) and AT/TA dinucleotides were found to be the most important. In rice, high roll (another DNA shape feature) and low CA dinucleotide frequency in particular were found to be associated with CO occurrence. The accuracy of our models was sufficient for Arabidopsis and rice (area under receiver operating characteristic curve, AUROC > 0.5), and was high for tomato and maize (AUROC â‰« 0.5), demonstrating that DNA sequence and shape are predictive for meiotic COs throughout the plant kingdom.

Transcription Factor-Mediated Control of Anthocyanin Biosynthesis in Vegetative Tissues.

Plants accumulate secondary metabolites to adapt to environmental conditions. These compounds, here exemplified by the purple-colored anthocyanins, are accumulated upon high temperatures, UV-light, drought, and nutrient deficiencies, and may contribute to tolerance to these stresses. Producing compounds is often part of a more broad response of the plant to changes in the environment. Here we investigate how a transcription-factor-mediated program for controlling anthocyanin biosynthesis also has effects on formation of specialized cell structures and changes in the plant root architecture. A systems biology approach was developed in tomato (Solanum lycopersicum) for coordinated induction of biosynthesis of anthocyanins, in a tissue- and development-independent manner. A transcription factor couple from Antirrhinum that is known to control anthocyanin biosynthesis was introduced in tomato under control of a dexamethasone-inducible promoter. By application of dexamethasone, anthocyanin formation was induced within 24 h in vegetative tissues and in undifferentiated cells. Profiles of metabolites and gene expression were analyzed in several tomato tissues. Changes in concentration of anthocyanins and other phenolic compounds were observed in all tested tissues, accompanied by induction of the biosynthetic pathways leading from Glc to anthocyanins. A number of pathways that are not known to be involved in anthocyanin biosynthesis were observed to be regulated. Anthocyanin-producing plants displayed profound physiological and architectural changes, depending on the tissue, including root branching, root epithelial cell morphology, seed germination, and leaf conductance. The inducible anthocyanin-production system reveals a range of phenomena that accompanies anthocyanin biosynthesis in tomato, including adaptions of the plants architecture and physiology.

Distribution, position and genomic characteristics of crossovers in tomato recombinant inbred lines derived from an interspecific cross between Solanum lycopersicum and Solanum pimpinellifolium.

We determined the crossover (CO) distribution, frequency and genomic sequences involved in interspecies meiotic recombination by using parent-assigned variants of 52 F6 recombinant inbred lines obtained from a cross between tomato, Solanum lycopersicum, and its wild relative, Solanum pimpinellifolium. The interspecific CO frequency was 80% lower than reported for intraspecific tomato crosses. We detected regions showing a relatively high and low CO frequency, so-called hot and cold regions. Cold regions coincide to a large extent with the heterochromatin, although we found a limited number of smaller cold regions in the euchromatin. The CO frequency was higher at the distal ends of chromosomes than in pericentromeric regions and higher in short arm euchromatin. Hot regions of CO were detected in euchromatin, and COs were more often located in non-coding regions near the 5' untranslated region of genes than expected by chance. Besides overrepresented CCN repeats, we detected poly-A/T and AT-rich motifs enriched in 1-kb promoter regions flanking the CO sites. The most abundant sequence motifs at CO sites share weak similarity to transcription factor-binding sites, such as for the C2H2 zinc finger factors class and MADS box factors, while InterPro scans detected enrichment for genes possibly involved in the repair of DNA breaks.

Tulipa gesneriana and Lilium longiflorum PEBP Genes and Their Putative Roles in Flowering Time Control.

Floral induction in Tulipa gesneriana and Lilium longiflorum is triggered by contrasting temperature conditions, high and low temperature, respectively. In Arabidopsis, the floral integrator FLOWERING LOCUS T (FT), a member of the PEBP (phosphatidyl ethanolamine-binding protein) gene family, is a key player in flowering time control. In this study, one PEBP gene was identified and characterized in lily (LlFT) and three PEBP genes were isolated from tulip (TgFT1, TgFT2 and TgFT3). Overexpression of these genes in Arabidopsis thaliana resulted in an early flowering phenotype for LlFT and TgFT2, but a late flowering phenotype for TgFT1 and TgFT3. Overexpression of LlFT in L. longiflorum also resulted in an early flowering phenotype, confirming its proposed role as a flowering time-controlling gene. The tulip PEBP genes TgFT2 and TgFT3 have a similar expression pattern in tulip, but show opposite effects on the timing of flowering in Arabidopsis. Therefore, the difference between these two proteins was further investigated by interchanging amino acids thought to be important for the FT function. This resulted in the conversion of phenotypes in Arabidopsis upon overexpressing the substituted TgFT2 and TgFT3 genes, revealing the importance of these interchanged amino acid residues. Based on all obtained results, we hypothesize that LlFT is involved in creating meristem competence to flowering-related cues in lily, and TgFT2 is considered to act as a florigen involved in the floral induction in tulip. The function of TgFT3 remains unclear, but, based on our observations and phylogenetic analysis, we propose a bulb-specific function for this gene.