Genome-Wide Identification and Expression Analysis of the MADS Gene Family in Tulips (Tulipa gesneriana).

To investigate the cold response mechanism and low temperature regulation of flowering in tulips, this study identified 32 MADS-box transcription factor family members in tulips based on full-length transcriptome sequencing, named TgMADS1-TgMADS32. Phylogenetic analysis revealed that these genes can be divided into two classes: type I and type II. Structural analysis showed that TgMADS genes from different subfamilies have a similar distribution of conserved motifs. Quantitative real-time PCR results demonstrated that some TgMADS genes (e.g., TgMADS3, TgMADS15, TgMADS16, and TgMADS19) were significantly upregulated in buds and stems under cold conditions, implying their potential involvement in the cold response of tulips. In summary, this study systematically identified MADS family members in tulips and elucidated their evolutionary relationships, gene structures, and cold-responsive expression patterns, laying the foundation for further elucidating the roles of these transcription factors in flowering and the cold adaptability of tulips.

Biological effects of gamma-ray radiation on tulip (Tulipa gesneriana L.).

Tulip, being an important ornamental plant, generally requires lengthy and laborious procedures to develop new varieties using traditional breeding methods requires. But ionizing radiation potentially accelerates the breeding process of ornamental plant species. The biological effects of γ-ray irradiation on tulip, therefore, were investigated through establishing an irradiation-mediated mutation breeding protocol to accelerate its breeding process. ISSR-PCR molecular marker technique was further used to identify the mutants of phenotypic variation plants. This study showed that low irradiation doses (5 Gy) stimulated bulb germination to improve the survival rate of tulip, while high irradiation doses (20 to 100 Gy) significantly (P < 0.05) inhibited its seed germination and growth, and decreased the flowering rate, petal number, flower stem length and flower diameter. More than 40 Gy significantly (P < 0.05) decreased the total chlorophyll content and increased the malondialdehyde (MDA) content in tulips. Interestingly, three types of both stigma variations and flower pattern variations, and four types of flower colour variations were observed. With increasing the irradiation dose from 5 to 100 Gy, the anthocyanin and flavonoid contents continuously decreased. Scanning electron microscopy (SEM) analysis evidenced that high irradiation doses altered the micromorphology of leaf stomata. Microscopic observations of tulip root apical mitosis further showed the abnormal chromosomal division behaviour occurring at different mitotic phases under irradiation treatment (80 Gy). Increasing the irradiation dose from 20 to 100 Gy enhanced the micronucleus rate. Moreover, the suspected genetic variation in tulips was evaluated by inter-simple sequence repeat (ISSR) analysis, and the percentage of polymorphic bands was 68%. Finally, this study concludes that that 80 Gy may be an appropriate radiation does to better enhance the efficiency of mutagenic breeds in tulip plants. Using γ-ray irradiation, therefore, is expected to offer a theoretical basis for mutation breeding in tulips.

Phylogeny and Comparative Analysis for the Plastid Genomes of Five Tulipa (Liliaceae).

Species of Tulipa (Liliaceae) are of great horticultural importance and are distributed across Europe, North Africa, and Asia. The Tien Shan Mountain is one of the primary diversity centres of Tulipa, but the molecular studies of Tulipa species from this location are lacking. In our study, we assembled four Tulipa plastid genomes from the Tien Shan Mountains, T. altaica, T. iliensis, T. patens, and T. thianschanica, combined with the plastid genome of T. sylvestris to compare against other Liliaceae plastid genomes. We focussed on the species diversity and evolution of their plastid genomes. The five Tulipa plastid genomes proved highly similar in overall size (151,691-152,088 bp), structure, gene order, and content. With comparative analysis, we chose 7 mononucleotide SSRs from the Tulipa species that could be used in further population studies. Phylogenetic analyses based on 24 plastid genomes robustly supported the monophyly of Tulipa and the sister relationship between Tulipa and Amana, Erythronium. T. iliensis, T. thianschanica, and T. altaica were clustered together, and T. patens was clustered with T. sylvestris, with our results clearly demonstrating the relationships between these five Tulipa species. Our results provide a more comprehensive understanding of the phylogenomics and comparative genomics of Tulipa.

Integrating physiological and metabolites analysis to identify ethylene involvement in petal senescence in Tulipa gesneriana.

Flower senescence is classified into ethylene-dependent and ethylene-independent manners and determines the flower longevity which is valuable for ornamental plants. However, the manner of petal senescence in tulip is still less defined. In this study, we characterized the physiological indexes in the process of petal senescence, as well as metabolic and ethylene responses in tulip cultivar 'American Dream', and further identified the role of ethylene biosynthesis genes TgACS by transgenic and transient assays. Primary metabolites profiling revealed that sugars, amino acids and organic acids preferentially accumulated in senescent petals. Additionally, senescence-associated genes were identified and significantly up-regulated, coupled with increased ROS contents, rapid water loss and accelerated cell membrane breakdown. Moreover, ethylene production was stimulated as evidenced by increasing in ACS activity and ethylene biosynthesis-related genes expression. Exogenous treatment of cutting flowers with 1-MCP or ethephon resulted in delayed or enhanced petal senescence, respectively. Transient down-regulation of TgACS by VIGS assay in tulip petals delayed senescence, while over-expressed TgACS1 in tobacco promoted leaf senescence. Taken together, this study provides evidences to certify ethylene roles and TgACS functions during flower senescence in tulip.

One-Step Enzymatic Synthesis of 1-Tuliposide A Using Tuliposide-Converting Enzyme.

6-Tuliposides A (6-PosA) and B (6-PosB) are major secondary metabolites in tulip (Tulipa gesneriana), having an acyl group at the C-6 position of D-glucose. They serve as precursors of the antimicrobial α-methylene-γ-butyrolactones tulipalins A (PaA) and B (PaB). The conversions of 6-PosA/6-PosB to PaA/PaB are catalyzed by tuliposide-converting enzymes A and B (TCEA and TCEB), respectively. A minor Pos, 1-PosA, which has the acyl group at the C-1 position of D-glucose, has been identified in some wild tulip species, but availability of this compound is limited. Here, by using the TCEs, we established a facile enzymatic process for 1-PosA synthesis from the naturally occurring 1,6-diacyl-glucose type of Pos (PosD and PosF). We first discovered that TCEA and TCEB react preferentially with PosD and PosF, respectively, to form 1-PosA and the corresponding Pa derived from the 6-acyl group, demonstrating that the TCEs specifically acted on the 6-acyl group, but not the 1-acyl group, of the substrates. Using TCEB, 300 mg of PosF was completely converted to 1-PosA and PaB in 10 min at room temperature. Then, 160 mg of 1-PosA (75% molar yield) was purified by column chromatography. This one-step enzymatic process dramatically improves accessibility to 1-PosA.

Role of Tulipa gesneriana TEOSINTE BRANCHED1 (TgTB1) in the control of axillary bud outgrowth in bulbs.

KEY MESSAGE: Tulip vegetative reproduction. Tulips reproduce asexually by the outgrowth of their axillary meristems located in the axil of each bulb scale. The number of axillary meristems in one bulb is low, and not all of them grow out during the yearly growth cycle of the bulb. Since the degree of axillary bud outgrowth in tulip determines the success of their vegetative propagation, this study aimed at understanding the mechanism controlling the differential axillary bud activity. We used a combined physiological and "bottom-up" molecular approach to shed light on this process and found that first two inner located buds do not seem to experience dormancy during the growth cycle, while mid-located buds enter dormancy by the end of the growing season. Dormancy was assessed by weight increase and TgTB1 expression levels, a conserved TCP transcription factor and well-known master integrator of environmental and endogenous signals influencing axillary meristem outgrowth in plants. We showed that TgTB1 expression in tulip bulbs can be modulated by sucrose, cytokinin and strigolactone, just as it has been reported for other species. However, the limited growth of mid-located buds, even when their TgTB1 expression is downregulated, points at other factors, probably physical, inhibiting their growth. We conclude that the time of axillary bud initiation determines the degree of dormancy and the sink strength of the bud. Thus, development, apical dominance, sink strength, hormonal cross-talk, expression of TgTB1 and other possibly physical but unidentified players, all converge to determine the growth capacity of tulip axillary buds.

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.

Molecular diversity of tuliposide B-converting enzyme in tulip (Tulipa gesneriana): identification of the root-specific isozyme.

6-Tuliposide B (PosB) is a glucose ester accumulated in tulip (Tulipa gesneriana) as a major secondary metabolite. PosB serves as the precursor of the antimicrobial lactone tulipalin B (PaB), which is formed by PosB-converting enzyme (TCEB). The gene TgTCEB1, encoding a TCEB, is transcribed in tulip pollen but scarcely transcribed in other tissues (e.g. roots) even though those tissues show high TCEB activity. This led to the prediction of the presence of a TCEB isozyme with distinct tissue specificity. Herein, we describe the identification of the TgTCEB-R gene from roots via native enzyme purification; this gene is a paralog of TgTCEB1. Recombinant enzyme characterization verified that TgTCEB-R encodes a TCEB. Moreover, TgTCEB-R was localized in tulip plastids, as found for pollen TgTCEB1. TgTCEB-R is transcribed almost exclusively in roots, indicating a tissue preference for the transcription of TCEB isozyme genes.

Molecular Regulation of Temperature-Dependent Floral Induction in Tulipa gesneriana.

The vegetative-to-reproductive phase change in tulip (Tulipa gesneriana) is promoted by increasing temperatures during spring. The warm winters of recent years interfere with this process and are calling for new adapted cultivars. A better understanding of the underlying molecular mechanisms would be of help, but unlike the model plant Arabidopsis (Arabidopsis thaliana), very little is known about floral induction in tulip. To shed light on the gene regulatory network controlling flowering in tulip, RNA sequencing was performed on meristem-enriched tissue collected under two contrasting temperature conditions, low and high. The start of reproductive development correlated with rounding of the shoot apical meristem and induction of TGSQA expression, a tulip gene with a high similarity to Arabidopsis APETALA1 Gene Ontology enrichment analysis of differentially expressed genes showed the overrepresentation of genes potentially involved in floral induction, bulb maturation, and dormancy establishment. Expression analysis revealed that TERMINAL FLOWER1 (TgTFL1) and SUPPRESSOR OF OVEREXPRESSION OF CONSTANS1-like1 (TgSOC1-like1) might be repressors, whereas TgSOC1-like2 likely is an activator, of flowering. Subsequently, the flowering time-associated expression of eight potential flowering time genes was confirmed in three tulip cultivars grown in the field. Additionally, heterologous functional analyses in Arabidopsis resulted in flowering time phenotypes in line with TgTFL1 being a floral repressor and TgSOC1-like2 being a floral activator in tulip. Taken together, we have shown that long before morphological changes occur in the shoot apical meristem, the expression of floral repressors in tulip is suppressed by increased ambient temperatures, leading either directly or indirectly to the activation of potential flowering activators shortly before the commencement of the phase change.

Elucidating and mining the Tulipa and Lilium transcriptomes.

Genome sequencing remains a challenge for species with large and complex genomes containing extensive repetitive sequences, of which the bulbous and monocotyledonous plants tulip and lily are examples. In such a case, sequencing of only the active part of the genome, represented by the transcriptome, is a good alternative to obtain information about gene content. In this study we aimed to generate a high quality transcriptome of tulip and lily and to make this data available as an open-access resource via a user-friendly web-based interface. The Illumina HiSeq 2000 platform was applied and the transcribed RNA was sequenced from a collection of different lily and tulip tissues, respectively. In order to obtain good transcriptome coverage and to facilitate effective data mining, assembly was done using different filtering parameters for clearing out contamination and noise of the RNAseq datasets. This analysis revealed limitations of commonly applied methods and parameter settings used in de novo transcriptome assembly. The final created transcriptomes are publicly available via a user friendly Transcriptome browser ( http://www.bioinformatics.nl/bulbs/db/species/index ). The usefulness of this resource has been exemplified by a search for all potential transcription factors in lily and tulip, with special focus on the TCP transcription factor family. This analysis and other quality parameters point out the quality of the transcriptomes, which can serve as a basis for further genomics studies in lily, tulip, and bulbous plants in general.

[Genetic polymorphism of Tulipa gesneriana L. evaluated on the basis of the ISSR marking data].

Using the method of ISSR analysis, the genetic diversity of 18 natural populations of Tulipa gesneriana L. from the north of the Lower Volga region was examined. The ten ISSR primers used in the study provided identification of 102 PCR fragments, of which 50 were polymorphic (49.0%). According to the proportion of polymorphic markers, two population groups were distinguished: (1) the populations in which the proportion of polymorphic markers ranged from 0.35 to 0.41; (2) the populations in which the proportion of polymorphic markers ranged from 0.64 to 0.85. UPGMA clustering analysis provided subdivision of the sample into two large clusters. The unrooted tree constructed using the Neighbor Joining algorithm had similar topology. The first cluster included slightly variable populations and the second cluster included highly variable populations. The AMOVA analysis showed statistically significant differences (F CT = 0.430; p = 0.000) between the two groups. Local populations are considerably genetically differentiated from each other (F ST = 0.632) and have almost no links via modern gene flow, as evidenced by the results of the Mantel test (r =­0.118; p = 0.819). It is suggested that the degree of genetic similarities and differences between the populations depends on the time and the species dispersal patterns on these territories.

[Morphologic and AFLP analysis of relationships between tulip species Tulipa biebersteiniana (Liliaceae)].

In populations of four species of tulips, (Tulipa biebersteiniana, T. patens, T. scytica and T. riparia) from the Volgograd, Kurgansk, Orenburg, and Chelyabinsk regions and the Republic of Bashkortostan, genetic diversity was studied by means of morphological and AFLP analysis. A morphological analysis of seven quantitative and two qualitative criteria was carried out. Three selective EcoRI/MseI primer pairs allowed one to genotype 81 individuals from 13 tulip populations with 87 loci. The low level of variability by AFLP loci were revealed in all species, including T. biebersteiniana (P = 20.41%, UH(e) = 0.075), T. patens (26.97%, 0.082), T. scytica (27.53%, 0.086), and T. riparia (27.72%, 0.096). According to the AMOVA results, the variability proportion that characterizes the differences between the four Tulip species was lower (F(CT) = 0.235) than between populations within species (F(ST) = 0.439). Tulipa patens is well differentiated by means of Nei's distances, coordination, and analysis in the STRUCTURE program. An analysis in the STRUCTURE revealed four genetic groups of tulips that are not completely in accordance with the analyzed species. This acknowledges the presence of complicated genetic process in the tulip population.

Molecular diversity of tuliposide A-converting enzyme in the tulip.

Tuliposide A-converting enzyme (TCEA) catalyzes the conversion of 6-tuliposide A to its lactonized aglycon, tulipalin A, in the tulip (Tulipa gesneriana). The TgTCEA gene, isolated previously from petals, was transcribed in all tulip tissues but not in the bulbs despite the presence of TCEA activity, which allowed prediction of the presence of a TgTCEA isozyme gene preferentially expressed in the bulbs. Here, the TgTCEA-b gene, the TgTCEA homolog, was identified in bulbs. TgTCEA-b polypeptides showed approximately 77% identity to the petal TgTCEA. Functional characterization of the recombinant enzyme verified that TgTCEA-b encoded the TCEA. Moreover, the TgTCEA-b was found to be localized to plastids, as found for the petal TgTCEA. Transcript analysis revealed that TgTCEA-b was functionally transcribed in the bulb scales, unlike the TgTCEA gene, whose transcripts were absent there. In contrast, TgTCEA-b transcripts were in the minority in other tissues where TgTCEA transcripts were dominant, indicating a tissue preference for the transcription of those isozyme genes.

Expression profiles of aquaporin homologues and petal movement during petal development in Tulipa gesneriana.

Previously, we have characterized two tonoplast intrinsic proteins (TIPs) and four plasma membrane intrinsic proteins (PIPs) from the 2-day-old petals of tulip (Tulipa gesneriana). In this study, we analyzed the development of tulip petals and stems, temperature-dependent petal movement, the amount of ³H2O transported into petals and stems during petal movement, and the transcript levels of two TIP (TgTIP1;1 and TgTIP1;2) and four TgPIP genes in petals and stems, from the first day of petal opening to day 12. The development of the petals and stems was completed by days 6 and 9, respectively, after the first day of petal opening. Temperature-dependent petal movement and the amount of ³H2O that was transported into petals could be detected at significant levels up to day 6 with petal movement reaching a peak at day 3. Real-time reverse transcription-polymerase chain reaction analysis revealed that TgTIP1;1 and TgTIP1;2 were expressed ubiquitously in petals, stems, leaves, bulbs and roots. However, the expression level of TgTIP1;2 was very low in bulbs. The expression of both TgTIP1 genes was upregulated in close association with the development of petals but not with that of the stem. The four TgPIP genes were expressed at almost the same level during the development of the petals and the stem. However, the levels of the TgTIP1 and TgPIP transcripts in petals decreased during the course of petal wilting from day 9 onwards. These results suggest that TgTIP1;1 and TgTIP1;2 may contribute to petal development.

Generation and analysis of expressed sequence tags in the extreme large genomes Lilium and Tulipa.

BACKGROUND: Bulbous flowers such as lily and tulip (Liliaceae family) are monocot perennial herbs that are economically very important ornamental plants worldwide. However, there are hardly any genetic studies performed and genomic resources are lacking. To build genomic resources and develop tools to speed up the breeding in both crops, next generation sequencing was implemented. We sequenced and assembled transcriptomes of four lily and five tulip genotypes using 454 pyro-sequencing technology. RESULTS: Successfully, we developed the first set of 81,791 contigs with an average length of 514 bp for tulip, and enriched the very limited number of 3,329 available ESTs (Expressed Sequence Tags) for lily with 52,172 contigs with an average length of 555 bp. The contigs together with singletons covered on average 37% of lily and 39% of tulip estimated transcriptome. Mining lily and tulip sequence data for SSRs (Simple Sequence Repeats) showed that di-nucleotide repeats were twice more abundant in UTRs (UnTranslated Regions) compared to coding regions, while tri-nucleotide repeats were equally spread over coding and UTR regions. Two sets of single nucleotide polymorphism (SNP) markers suitable for high throughput genotyping were developed. In the first set, no SNPs flanking the target SNP (50 bp on either side) were allowed. In the second set, one SNP in the flanking regions was allowed, which resulted in a 2 to 3 fold increase in SNP marker numbers compared with the first set. Orthologous groups between the two flower bulbs: lily and tulip (12,017 groups) and among the three monocot species: lily, tulip, and rice (6,900 groups) were determined using OrthoMCL. Orthologous groups were screened for common SNP markers and EST-SSRs to study synteny between lily and tulip, which resulted in 113 common SNP markers and 292 common EST-SSR. Lily and tulip contigs generated were annotated and described according to Gene Ontology terminology. CONCLUSIONS: Two transcriptome sets were built that are valuable resources for marker development, comparative genomic studies and candidate gene approaches. Next generation sequencing of leaf transcriptome is very effective; however, deeper sequencing and using more tissues and stages is advisable for extended comparative studies.

A novel lactone-forming carboxylesterase: molecular identification of a tuliposide A-converting enzyme in tulip.

Tuliposides, the glucose esters of 4-hydroxy-2-methylenebutanoate and 3,4-dihydroxy-2-methylenebutanoate, are major secondary metabolites in tulip (Tulipa gesneriana). Their lactonized aglycons, tulipalins, function as defensive chemicals due to their biological activities. We recently found that tuliposide-converting enzyme (TCE) purified from tulip bulbs catalyzed the conversion of tuliposides to tulipalins, but the possibility of the presence of several TCE isozymes was raised: TCE in tissues other than bulbs is different from bulb TCE. Here, to prove this hypothesis, TCE was purified from petals, which have the second highest TCE activity after bulbs. The purified enzyme, like the bulb enzyme, preferentially accepted tuliposides as substrates, with 6-tuliposide A the best substrate, which allowed naming the enzyme tuliposide A-converting enzyme (TCEA), but specific activity and molecular mass differed between the petal and bulb enzymes. After peptide sequencing, a novel cDNA (TgTCEA) encoding petal TCEA was isolated, and the functional characterization of the recombinant enzyme verified that TgTCEA catalyzes the conversion of 6-tuliposide A to tulipalin A. TgTCEA was transcribed in all tulip tissues but not in bulbs, indicating the presence of a bulb-specific TgTCEA, as suggested by the distinct enzymatic characters between the petal and bulb enzymes. Plastidial localization of TgTCEA enzyme was revealed, which allowed proposing a cytological mechanism of TgTCE-mediated tulipalin formation in the tulip defensive strategy. Site-directed mutagenesis of TgTCEA suggested that the oxyanion hole and catalytic triad characteristic of typical carboxylesterases are essential for the catalytic process of TgTCEA enzyme. To our knowledge, TgTCEA is the first identified member of the lactone-forming carboxylesterases, specifically catalyzing intramolecular transesterification.

Substitution of a single amino acid residue in the aromatic/arginine selectivity filter alters the transport profiles of tonoplast aquaporin homologs.

Aquaporins are integral membrane proteins that facilitate the transport of water and some small solutes across cellular membranes. X-ray crystallography of aquaporins indicates that four amino acids constitute an aromatic/arginine (ar/R) pore constriction known as the selectivity filter. On the basis of these four amino acids, tonoplast aquaporins called tonoplast intrinsic proteins (TIPs) are divided into three groups in Arabidopsis. Herein, we describe the characterization of two group I TIP1s (TgTIP1;1 and TgTIP1;2) from tulip (Tulipa gesneriana). TgTIP1;1 and TgTIP1;2 have a novel isoleucine in loop E (LE2 position) of the ar/R filter; the residue at LE2 is a valine in all group I TIPs from model plants. The homologs showed mercury-sensitive water channel activity in a fast kinetics swelling assay upon heterologous expression in Pichia pastoris. Heterologous expression of both homologs promoted the growth of P. pastoris on ammonium or urea as sole sources of nitrogen and decreased growth and survival in the presence of H(2)O(2). TgTIP1;1- and TgTIP1;2-mediated H(2)O(2) conductance was demonstrated further by a fluorescence assay. Substitutions in the ar/R selectivity filter of TgTIP1;1 showed that mutants that mimicked the ar/R constriction of group I TIPs could conduct the same substrates that were transported by wild-type TgTIP1;1. In contrast, mutants that mimicked group II TIPs showed no evidence of urea or H(2)O(2) conductance. These results suggest that the amino acid residue at LE2 position is critical for the transport selectivity of the TIP homologs and group I TIPs might have a broader spectrum of substrate selectivity than group II TIPs.

Micropropagation of tulip: production of virus-free stock plants.

We describe here a new tulip micropropagation method based on the cyclic shoot multiplication in presence of the thidiazuron (TDZ), which enables the production of virus-free stock plants, speeds up breeding, and provides new genotypes for the market. In our novel protocol, cyclic shoot multiplication can be performed for 2-3 years by using TDZ instead of other cytokinins, as 6-benzylaminopurine (BAP) and N(6)-(-isopentyl)adenine (2iP). It makes possible to produce 500-2,000 microbulbs from one healthy plant. There are six main stages of tulip micropropagation. Stage 0 is the selection of true-to-type and virus-free plants, confirmed by ELISA. Fragments of flower stems isolated from bulbs are used as initial explants. Shoot multiplication is based on the regeneration of adventitious shoots, which are sub-cultured every 8 weeks. In the Stage 3, the specially prepared shoots are induced by low temperature treatment to form bulbs which finally develop on a sucrose-rich medium at 20 degrees C. Bulbs are then dried for 6 weeks and rooted in vivo. The number of multiplication subcultures should be limited to 5-10 cycles in order to lower the risk of mutation. Virus indexing should be repeated 3-4 times, at the initial stage and then during shoot multiplication. Genetic stability of micropropagated shoots can be confirmed using molecular markers.

Alternative expression of vacuolar iron transporter and ferritin genes leads to blue/purple coloration of flowers in tulip cv. 'Murasakizuisho'.

Flowers of tulip cv. 'Murasakizuisho' have a purple perianth except for the bottom region, which is blue in color even though it has the same anthocyanin, delphinidin 3-O-rutinoside, as the entire perianth. The development of the blue coloration in the perianth bottom is due to complexation by anthocyanin, flavonol and iron (Fe), as well as a vacuolar iron transporter, TgVit1. Although transient expression of TgVit1 in the purple cells led to a color change to light blue, the coloration of the transformed cells did not coincide with the dark blue color of the cells of the perianth bottom. We thought that another factor is required for the blue coloration of the cells of perianth bottom. To examine the effect of ferritin (FER), an Fe storage protein, on blue color development, we cloned an FER gene (TgFER1) and performed expression analyses. TgFER1 transcripts were found in the cells located in the upper region of the petals along with purple color development by anthocyanin and were not found in the blue cells of the perianth bottom. This gene expression is in contrast to that of TgVit1, expressed only in the cells of the perianth bottom. Co-expression of TgVIT1 and TgFER-RNAi, constructed for suppressing endogenous TgFER1 by RNA interference (RNAi), changed the purple petal cells to a dark blue color similar to that of the natural perianth bottom. These results strongly suggest that TgVit1 expression and TgFER1 suppression are critical for the development of blue color in the perianth bottom.

Diploid endosperm formation in Tulipa spp. and identification of a 1:1 maternal-to-paternal genome ratio in endosperms of T. gesneriana L.

Most Liliaceae plants have the tetrasporic Fritillaria-type embryo sac and normally form diploid embryos and pentaploid endosperms derived from a 4:1 maternal-to-paternal genome ratio (4m:1p) after double fertilization. Here we characterize embryo sac and endosperm formation in Tulipa spp. of Liliaceae. Chromosome analysis using seeds derived from 2x x 2x crosses of Tulipa gesneriana (2n = 2x = 24) identified diploid chromosome number in the endosperm. Similarly, flow cytometric analysis confirmed diploid endosperm formation in T. gesneriana, T. fosteriana (2n = 2x = 24) and T. greigii (2n = 2x = 24). To further study the possible mechanism of diploid endosperm formation, we made interploidy crosses of triploid (2n = 3x = 36) x diploid in which aneuploid seeds with various chromosome numbers (2n = 25-36) were produced. Again, flow cytometric analysis confirmed the same ploidy level in both embryos and endosperms at all aneuploidy levels, suggesting that only a single haploid polar nucleus contributes to endosperm formation at fertilization. Histological observation further confirmed the physical separation of two polar nuclei by a large vacuole in the Fritillaria-type embryo sac of T. gesneriana that appeared to prevent the fusion of the two polar nuclei that originated at the micropylar and chalazal ends before fertilization. Taken together, these results indicate that diploid endosperms (1m:1p) are normally formed in Tulipa spp. by fusion of the micropylar polar nucleus (n) and a spermatid (n) but not by normal triple fusion. We also show that tulip endosperm partially overcomes the triploid block mechanism that occurs in interploidy crosses. Based on these observations, the possible role of triple nuclear fusion in double fertilization is discussed.