Transcriptome profiling reveals key genes in regulation of the tepal trichome development in Lilium pumilum D.C.

KEY MESSAGE: A number of potential genes and pathways involved in tepal trichome development were identified in a natural lily mutant by transcriptome analysis and were confirmed with trichome and trichomeless species. Trichome is a specialized structure found on the surface of the plant with an important function in survival against abiotic and biotic stress. It is also an important economic trait in crop breeding. Extensive research has investigated the foliar trichome in model plants (Arabidopsis and tomato). However, the developmental mechanism of tepal trichome remains elusive. Lilium pumilum is an edible ornamental bulb and a good breeding parent possessing cold and salt-alkali resistance. Here, we found a natural mutant of Lilium pumilum grown on a highland whose tepals are covered by trichomes. Our data indicate that trichomes of the mutant are multicellular and branchless. Notably, stomata are also developed on the tepal of the mutant as well, suggesting there may be a correlation between trichome and stomata regulation. Furthermore, we isolated 27 differentially expressed genes (DEGs) by comparing the transcriptome profiling between the natural mutant and the wild type. These 27 genes belong to 4 groups: epidermal cell cycle and division, trichome morphogenesis, stress response, and transcription factors. Quantitative real-time PCR in Lilium pumilum (natural mutant and the wild type) and other lily species (Lilium leichtlinii var. maximowiczii/trichome; Lilium davidii var. willmottiae/, trichomeless) confirmed the validation of RNA-seq data and identified several trichome-related genes.

3D mechanical characterization of single cells and small organisms using acoustic manipulation and force microscopy.

Quantitative micromechanical characterization of single cells and multicellular tissues or organisms is of fundamental importance to the study of cellular growth, morphogenesis, and cell-cell interactions. However, due to limited manipulation capabilities at the microscale, systems used for mechanical characterizations struggle to provide complete three-dimensional coverage of individual specimens. Here, we combine an acoustically driven manipulation device with a micro-force sensor to freely rotate biological samples and quantify mechanical properties at multiple regions of interest within a specimen. The versatility of this tool is demonstrated through the analysis of single Lilium longiflorum pollen grains, in combination with numerical simulations, and individual Caenorhabditis elegans nematodes. It reveals local variations in apparent stiffness for single specimens, providing previously inaccessible information and datasets on mechanical properties that serve as the basis for biophysical modelling and allow deeper insights into the biomechanics of these living systems.

The male gamete membrane protein DMP9/DAU2 is required for double fertilization in flowering plants.

All flowering plants exhibit a unique type of sexual reproduction called 'double fertilization' in which each pollen tube-delivered sperm cell fuses with an egg and a central cell. Proteins that localize to the plasma membrane of gametes regulate one-to-one gamete pairing and fusion between male and female gametes for successful double fertilization. Here, we have identified a membrane protein from Lilium longiflorum generative cells using proteomic analysis and have found that the protein is an ortholog of Arabidopsis DUF679 DOMAIN MEMBRANE PROTEIN 9 (DMP9)/DUO1-ACTIVATED UNKNOWN 2 (DAU2). The flowering plant DMP9 proteins analyzed in this study were predicted to have four transmembrane domains and be specifically expressed in both generative and sperm cells. Knockdown of DMP9 resulted in aborted seeds due to single fertilization of the central cell. Detailed imaging of DMP9-knockdown sperm cells during in vivo and semi-in vitro double fertilization revealed that DMP9 is involved in gamete interaction that leads to correct double fertilization.

Fountain streaming contributes to fast tip-growth through regulating the gradients of turgor pressure and concentration in pollen tubes.

Fountain streaming is a typical microfluidic pattern in plant cells, especially for cells with a high aspect ratio such as pollen tubes. Although it has been found that fountain streaming plays crucial roles in the transport of nutrients and metabolites, the positioning of organelles and the mixing of cytoplasms, its implications for the fast tip growth of pollen tubes remain a mystery. To address this, based on the observations of asiatic lily Lilium Casablanca, we developed physical models for reverse fountain streaming in pollen tubes and solved the hydrodynamics and advection-diffusion dynamics of viscous Stokes flow in the shank and apical region of pollen tubes. Theoretical and numerical results demonstrated that the gradients of turgor pressure and concentration of wall materials along the length of pollen tubes provide undamped driving force and high-efficiency materials supply, which are supposed to contribute to the fast tip-growth of pollen tubes. The sample experimental results show that the tip-growth will be abnormal when the gradients of turgor pressure change under osmotic stress induced by different concentrations of PEG-6000 (a dehydrant).

LlFH1-mediated interaction between actin fringe and exocytic vesicles is involved in pollen tube tip growth.

Pollen tube tip growth is an extreme form of polarized cell growth, which requires polarized exocytosis based on a dynamic actin cytoskeleton. However, the molecular basis for the connection between actin filaments and exocytic vesicles is unclear. Here, we identified a Lilium longiflorum pollen-specific formin (LlFH1) and found that it localized at the apical vesicles and plasma membrane (PM). Overexpression of LlFH1 induced excessive actin cables in the tube tip region, and downregulation of LlFH1 eliminated the actin fringe. Fluorescence recovery after photobleaching (FRAP) analysis revealed that LlFH1-labeled exocytic vesicles exhibited an initial accumulation at the shoulder of the apex and coincided with the leading edge of the actin fringe. Time-lapse analysis suggested that nascent actin filaments followed the emergence of the apical vesicles, implying that LlFH1 at apical vesicles could initiate actin polymerization. Biochemical assays showed that LlFH1 FH1FH2 could nucleate actin polymerization, but then capped the actin filament at the barbed end and inhibited its elongation. However, in the presence of lily profilins, LlFH1 FH1FH2 could accelerate barbed-end actin elongation. In addition, LlFH1 FH1FH2 was able to bundle actin filaments. Thus, we propose that LlFH1 and profilin coordinate the interaction between the actin fringe and exocytic vesicle trafficking during pollen tube growth of lily.

High precision, localized proton gradients and fluxes generated by a microelectrode device induce differential growth behaviors of pollen tubes.

Pollen tubes are tip-growing plant cells that deliver the sperm cells to the ovules for double fertilization of the egg cell and the endosperm. Various directional cues can trigger the reorientation of pollen tube growth direction on their passage through the female tissues. Among the external stimuli, protons serve an important, regulatory role in the control of pollen tube growth. The generation of local guidance cues has been challenging when investigating the mechanisms of perception and processing of such directional triggers in pollen tubes. Here, we developed and characterized a microelectrode device to generate a local proton gradient and proton flux through water electrolysis. We confirmed that the cytoplasmic pH of pollen tubes varied with environmental pH change. Depending on the position of the pollen tube tip relative to the proton gradient, we observed alterations in the growth behavior, such as bursting at the tip, change in growth direction, or complete growth arrest. Bursting and growth arrest support the hypothesis that changes in the extracellular H+ concentration may interfere with cell wall integrity and actin polymerization at the growing tip. A change in growth direction for some pollen tubes implies that they can perceive the local proton gradient and respond to it. We also showed that the growth rate is directly correlated with the extracellular pH in the tip region. Our microelectrode approach provides a simple method to generate protons and investigate their effect on plant cell growth.

Ribosomal DNA locus variation and REMAP analysis of the diploid and triploid complexes of Lilium lancifolium.

Lilium lancifolium Thunb. (2n = 2x = 24) is a cytologically conspicuous species with both diploids and triploids in nature. Cytological and molecular genetic analyses were carried out in both diploids and triploids that were collected from 55 geographical locations in Korea, Japan, and China. While the 5S rRNA gene loci were located at duplicated loci on the long arm of chromosome 2, the 45S rRNA gene loci were present in chromosomes 1, 2, 4, 6, 7, and 11. While the loci on chromosomes 1 and 7 were constant, the loci on chromosomes 2, 4, 6, 7, and 11 were variable in some plants so that the L. lancifolium accessions were grouped into 7 cytotypes in diploids and 12 cytotypes in triploids. REMAP marker analysis revealed that the diploids were classified into seven clusters, and the triploids were classified into a large cluster. Geographic, cytological, and genetic differentiations were not related in both the diploid and triploid accessions of L. lancifolium. Thus, current genetic variations occurred prior to the geographic differentiation in both diploids and triploids, and the 45S rDNA cytotype variations occurred after geographic differentiation in the current habitats of L. lancifolium.

Hydrogen peroxide affects ion channels in lily pollen grain protoplasts.

Ion homeostasis plays a central role in polarisation and polar growth. In several cell types ion channels are controlled by reactive oxygen species (ROS). One of the most important cells in the plant life cycle is the male gametophyte, which grows under the tight control of both ion fluxes and ROS balance. The precise relationship between these two factors in pollen tubes has not been completely elucidated, and in pollen grains it has never been studied to date. In the present study we used a simple model - protoplasts obtained from lily pollen grains at the early germination stage - to reveal the effect of H2 O2 on cation fluxes crucial for pollen germination. Here we present direct evidence for two ROS-sensitive currents on the pollen grain plasma membrane: the hyperpolarisation-activated calcium current, which is strongly enhanced by H2 O2 , and the outward potassium current, which is modestly enhanced by H2 O2 . We used low concentrations of H2 O2 that do not cause an intracellular oxidative burst and do not damage cells, as demonstrated with fluorescent staining.

Proteomic analysis reveals the differential histone programs between male germline cells and vegetative cells in Lilium davidii.

In flowering plants, male germline fate is determined after asymmetric division of the haploid microspore. Daughter cells have distinct fates: the generative cell (GC) undergoes further mitosis to generate sperm cells (SCs), and the vegetative cell (VC) terminally differentiates. However, our understanding of the mechanisms underlying germline development remains limited. Histone variants and modifications define chromatin states, and contribute to establishing and maintaining cell identities by affecting gene expression. Here, we constructed a lily protein database, then extracted and detailed histone entries into a comprehensive lily histone database. We isolated large amounts of nuclei from VCs, GCs and SCs from lily, and profiled histone variants of all five histone families in all three cell types using proteomics approaches. We revealed 92 identities representing 32 histone variants: six for H1, 11 for H2A, eight for H2B, five for H3 and two for H4. Nine variants, including five H1, two H2B, one H3 and one H4 variant, specifically accumulated in GCs and SCs. We also detected H3 modification patterns in the three cell types. GCs and SCs had almost identical histone profiles and similar H3 modification patterns, which were significantly different from those of VCs. Our study also revealed the presence of multiple isoforms, and differential expression patterns between isoforms of a variant. The results suggest that differential histone programs between the germline and companion VCs may be established following the asymmetric division, and are important for identity establishment and differentiation of the male germline as well as the VC.

Intra- and intertissular cytomictic interactions in mic-rosporogenesis of mono- and dicotyledonous plants.

A comparative cytological analysis of intra- and intertissular cytomictic interactions in early micro-sporogenesis of mono- and dicotyledonous plants was performed by the example of the two cellular systems - microsporocytes and tapetum. It is found that cytomixis is the component of intratissular interactions mainly. In the tapetum cells cytomixis is notable for structural and temporary taxon specific features. The nuclear migration in microsporocytes is confined mainly to zygotene-pachytene meiotic stages and characterized by a certain synchronism with cytomixis at the tapetum. Intertissular cytomictic interactions (tapetum - microsporocytes) were found in the monocot anthers only. Intertissular interactions are likely to reflect the intensification of competitive relations between the tapetum and microsporocytes for area in the process of anther tissue differentiation. Polyploid tapetum nucleus and syncytia being powerful acceptors are able to compete with microsporocytes and direct the chromatin translocation to their favor. The absence of intertissular interactions in dicots probably reflects a better balance between the processes of differentiation at somatic and generative tissues into microsporangium compared to monocots.

Automated cell-by-cell tissue imaging and single-cell analysis for targeted morphologies by laser ablation electrospray ionization mass spectrometry.

Mass spectrometry imaging (MSI) is an emerging technology for the mapping of molecular distributions in tissues. In most of the existing studies, imaging is performed by sampling on a predefined rectangular grid that does not reflect the natural cellular pattern of the tissue. Delivering laser pulses by a sharpened optical fiber in laser ablation electrospray ionization (LAESI) mass spectrometry (MS) has enabled the direct analysis of single cells and subcellular compartments. Cell-by-cell imaging had been demonstrated using LAESI-MS, where individual cells were manually selected to serve as natural pixels for tissue imaging. Here we describe a protocol for a novel cell-by-cell LAESI imaging approach that automates cell recognition and addressing for systematic ablation of individual cells. Cell types with particular morphologies can also be selected for analysis. First, the cells are recognized as objects in a microscope image. The coordinates of their centroids are used by a stage-control program to sequentially position the cells under the optical fiber tip for laser ablation. This approach increases the image acquisition efficiency and stability, and enables the investigation of extended or selected tissue areas. In the LAESI process, the ablation events result in mass spectra that represent the metabolite levels in the ablated cells. Peak intensities of selected ions are used to represent the metabolite distributions in the tissue with single-cell resolution.

Characterization of the Dual Subcellular Localization of Lilium LsGRP1, a Plant Class II Glycine-Rich Protein.

The defense-related gene LsGRP1 exhibits an increased level of expression in Lilium spp. after being infected by Botrytis elliptica, the fungal pathogen of lily leaf blight. In this study, the expression profile of the LsGRP1 protein (a plant class II glycine-rich protein) was characterized biochemically and its subcellular localization in lily leaves was evaluated using immunohistochemistry, enhanced green fluorescent protein (EGFP) imaging, and protein extraction analysis. Using an LsGRP1-specific antibody, LsGRP1 was found to be most abundant in epidermal cells and phloem tissues. Leaves from lily plants at different growth stages demonstrated similar levels of 14- and 16-kDa LsGRP1 and a decreased amount of 23-kDa LsGRP1 at the senescence stage. LsGRP1-EGFP imaging and protein extraction assays revealed that 14-kDa LsGRP1 was located in the plasma membrane whereas 16- and 23-kDa LsGRP1 was weakly bound to the cell wall. The time course analyses of LsGRP1 expression in response to salicylic acid treatment or B. elliptica infection showed an increased accumulation of 14- and 23-kDa LsGRP1 over time. Because 23-kDa LsGRP1 could be detected by an ubiquitin antibody, conversion of 14-kDa to 23-kDa LsGRP1 via mono-ubiquitination was presumed, which is a phenomenon that has not been reported for a plant class II glycine-rich protein.

Molecular phylogeny and genetic variation in the genus Lilium native to China based on the internal transcribed spacer sequences of nuclear ribosomal DNA.

We present a comprehensive phylogeny derived from nuclear ribosomal DNA (nrDNA) for 214 samples representing 98 species and five varieties, including 44 species and five varieties native to China. Our collection of 25 species and five varieties (44 samples) covering all five sections of the genus (Comber) distributed in China also were included in the internal transcribed spacer (ITS) database. This study incorporates previous research with an emphasis on Chinese species, including the controversial subsection, Sinomartagon 5c Comber. In the phylogenetic tree obtained by maximum parsimony (PAUP) and maximum likelihood (RAxML) analyses, the samples were divided into four major groups. Our results suggest that the subsection (subsect.) 5c Comber should be classified into the true subsect. 5c and the section (sect.) Lophophorum. And the latter was divided into three subsections (subsect. Lophophorum I, subsect. Lophophorum II, and subsect. Lophophorum III). Based on molecular phylogenetic analysis and fluorescence in situ hybridization, we report that L. henryi and L. rosthornii are closely related, and we propose their classification into subsect. Leucolirion 6a. Our results support Comber's subdivision of sect. Leucolirion, which was primarily based on bulb color. Chinese species were divided into five sections: sect. Martagon, sect. Archelirion, sect. Leucolirion, sect. Sinomartagon, and sect. Lophophorum. These findings contribute to our understanding of the phylogeny, origin, and classification of Lilium.

Comparative proteomic analysis of generative and sperm cells reveals molecular characteristics associated with sperm development and function specialization.

In flowering plants, two sperm cells (SCs) are generated from a generative cell (GC) in the developing pollen grain or growing pollen tube and are then delivered to the embryo sac to initiate double fertilization. SC development and function specialization involve the strict control of the protein (gene) expression program and coordination of diverse cellular processes. However, because methods for collecting a large amount of highly purified GCs and SCs for proteomic and transcriptomic studies from a plant are not available, molecular information about the program and the interconnections is lacking. Here, we describe a method for obtaining a large quantity of highly purified GCs and SCs from just-germinated lily pollen grains and growing pollen tubes for proteomic analysis. Our observation showed that SCs had less condensed chromatin and more vacuole-like structures than GCs and that mature SCs were arrested at the G2 phase. Comparison of SC and GC proteomes revealed 101 proteins differentially expressed in the two proteomes. These proteins are involved in diverse cellular and metabolic processes, with preferential involvement in metabolism, the cell cycle, signaling, the ubiquitin/proteasome pathway, and chromatin remodeling. Impressively, almost all proteins in SCF complex-mediated proteolysis and the cell cycle were up-regulated in SCs, whereas those in chromatin remodeling and stress response were down-regulated. Our data also reveal the coordination of SCF complex-mediated proteolysis, cell cycle progression, and DNA repair in SC development and function specialization. This study revealed for the first time a difference in protein profiles between GCs and SCs.

Basic rules for polarised cell growth.

Growth by cell elongation is a morphological process that transcends taxonomic kingdoms. Examples of this polarised growth form include hyphal tip growth in actinobacteria and filamentous fungi and pollen tube development. The biological processes required to produce polarisation in each of these examples are very different. However, commonality of the polarised growth habit suggests that certain "basic physical rules" of development are being followed. In this paper we are concerned with trying to further elucidate some of these basic rules. To this end, we focus on a simple and hence ubiquitous description of the polarised cell, its geometry, and using a mathematical model investigate how geometry and the deposition of new wall material could be related. We show that this simple model predicts both cell geometry and the location of maximal wall-deposition in a range of examples.

LlSR28 is involved in pollen germination by affecting filamentous actin dynamics.

Alternative splicing plays important roles in gene regulation and contributes to protein complexity. Previous studies suggest that alternative splicing exists in members of the villin/gelsolin/fragmin superfamily. In this study, a serine/argine-rich (SR) protein cDNA with 28 kDa protein (LlSR28) was isolated from a lily (Lilium longiflorum) expression library. Protein domain analysis showed that LlSR28 had similar structures to Arabidopsis SR45 (AtSR45), and LlSR28 could complement the phenotype of loss of AtSR45 function. Therefore, overexpression of LlSR28 and AtSR45 mutant (atsr45-1) were used in the following experiments. Overexpression of LlSR28 in Arabidopsis completely inhibited pollen germination. In contrast, the pollen germination of atsr45-1 was earlier than that of wild-type. In addition, pollen of atsr45-1 contained less F-actin at the corresponding hydration stage during pollen germination compared to that of wild-type. Alternative splicing analysis showed that Arabidopsis villin1 (AtVLN1) transcript encoding the full-length protein was increased, and that encoding the truncated protein was decreased in atst45-1. Moreover, the mRNA expression level of other actin-binding proteins (ABPs) abundant in Arabidopsis pollen was also changed in atsr45-1. In conclusion, we hypothesize that LlSR28 alters F-actin dynamics probably through its alternative splicing activities to affect directly or indirectly the alternative splicing of AtVLN1 and the expression of different ABPs, which then affects the pollen germination.

Generative cell-specific activation of the histone gH2A gene promoter of Lilium longiflorum in tobacco.

The Lilium longiflorum gH2A promoter is active exclusively in the generative cells of mature pollen in transgenic tobacco expressing the gH2A promoter::GUS (ß-glucuronidase) construct as a reporter gene. Temporal and spatial aspects of gH2A promoter activity examined during pollen development in transgenic tobacco reveal that GUS reporter activity was not detected until developing pollen entered the early bicellular developmental stage. Activity was first detected in generative cells at early-mid stages and gradually increased to maximum levels at mid-bicellular stages. The patterns of appearance and longevity of GUS activity in tobacco were very similar to those of gH2A mRNA during pollen development in Lilium. Exogenous treatment with colchicine, a well-known microtubule depolymerize, blocked microspore mitosis and inhibited generative cell differentiation. No GUS signal was detected in the resulting anomalous pollen, which lacked generative cell differentiation. These data strongly suggest that normal generative cell development is essential for activation of the gH2A promoter. Furthermore, these results indicate that common transcriptional activator(s) of the gH2A promoter may be present in both Lilium and Nicotiana, and that such putative factor(s) activates the gH2A promoter only when generative cells undergo normal development.

In vitro mutagenesis and identification of mutants via ISSR in lily (Lilium longiflorum).

An efficient in vitro mutagenesis protocol for Lilium longiflorum Thunb. cv. White fox has been established. The effect of 6-BA and NAA on adventitious bud formation from the bulblet-scale thin cell layers was tested. Results showed that the optimal medium for adventitious bud induction is MS basal medium supplemented with 2.0 mg/l 6-BA and 0.1 mg/l NAA. The differentiation frequency and the average number of adventitious buds reached 95.55% and 3.00, respectively. Various doses (0.0, 0.5, 1.0, 1.5, 2.0, and 2.5 Gy) of gamma rays were applied to investigate the effect of radiation on adventitious bud formation from bulblet-scale thin cell layers. The forming capacity of the adventitious buds significantly decreased with the increase of radiation dose. The results suggested that the optimal irradiation dose is 1.0 Gy. Dose of 1.0 Gy treatment resulted in 55.33% survival of irradiated bulblet-scale thin cell layers and 39.27% mutagenesis rate. The genetic variations among the morphological mutants were evaluated by DNA fingerprinting using ISSR molecular marker. The genetic variation frequency reached 36.06% using seven ISSR primers. Out of the 50 mutant lines transferred to the greenhouse, 9 were observed to have significantly different morphological characters than those of the controls.

Calcium-regulated anion channels in the plasma membrane of Lilium longiflorum pollen protoplasts.

• Currents through anion channels in the plasma membrane of Lilium longiflorum pollen grain protoplasts were studied under conditions of symmetrical anionic concentrations by means of patch-clamp whole-cell configuration. • With Cl(-) -based intra- and extracellular solutions, three outward-rectifying anion conductances, I(Cl1) , I(Cl2) and I(Cl3) , were identified. These three activities were discriminated by differential rundown behaviour and sensitivity to 5-nitro-2-(phenylpropylamino)-benzoate (NPPB), which could not be attributed to one or more channel types. All shared strong outward rectification, activated instantaneously and displayed a slow time-dependent activation for positive potentials. All showed modulation by intracellular calcium ([Ca(2+) ](in) ), increasing intensity from 6.04 nM up to 0.5 mM (I(Cl1) ), or reaching a maximum value with 8.50 µM (I(Cl2) and I(Cl3) ). • After rundown, the anionic currents measured using NO(3) (-) -based solutions were indistinguishable, indicating that the permeabilities of the channels for Cl(-) and NO(3) (-) are similar. Additionally, unitary anionic currents were measured from outside-out excised patches, confirming the presence of individual anionic channels. • This study shows for the first time the presence of a large anionic conductance across the membrane of pollen protoplasts, resulting from the presence of Ca(2+) -regulated channels. A similar conductance was also found in germinated pollen. We hypothesize that these putative channels may be responsible for the large anionic fluxes previously detected by means of self-referencing vibrating probes.

Isobaric tags for relative and absolute quantification- based comparative proteomics reveals the features of plasma membrane-associated proteomes of pollen grains and pollen tubes from Lilium davidii.

Mature pollen grains (PGs) from most plant species are metabolically quiescent. However, once pollinated onto stigma, they quickly hydrate and germinate. A PG can give rise to a vegetative cell-derived polarized pollen tube (PT), which represents a specialized polar cell. The polarized PT grows by the tip and requires interaction of different signaling molecules localized in the apical plasma membrane and active membrane trafficking. The mechanisms underlying the interaction and membrane trafficking are not well understood. In this work, we purified PG and PT plasma-membrane vesicles from Lilium davidii Duch. using the aqueous two-phase partition technique, then enriched plasma membrane proteins by using Brij58 and KCl to remove loosely bound contaminants. We identified 223 integral and membrane-associated proteins in the plasma membrane of PGs and PTs by using isobaric tags for relative and absolute quantification (iTRAQ) and 2-D high-performance liquid chromatography-tandem mass spectrometry. More than 68% of the proteins have putative transmembrane domains and/or lipid-modified motifs. Proteins involved in signal transduction, membrane trafficking and transport are predominant in the plasma-membrane proteome. We revealed most components of the clathrin-dependent endocytosis pathway. Statistical analysis revealed 14 proteins differentially expressed in the two development stages: in PTs, six upregulated and eight downregulated are mainly involved in signaling, transport and membrane trafficking. These results provide novel insights into polarized PT growth.