ORIGINAL RESEARCH article

Front. Plant Sci., 27 November 2018

Sec. Plant Development and EvoDevo

Volume 9 - 2018 | https://doi.org/10.3389/fpls.2018.01570

Class XI Myosins Contribute to Auxin Response and Senescence-Induced Cell Death in Arabidopsis

  • 1. Department of Chemistry and Biotechnology, Tallinn University of Technology, Tallinn, Estonia

  • 2. Department of Botany and Plant Pathology, Oregon State University, Corvallis, OR, United States

Abstract

The integrity and dynamics of actin cytoskeleton is necessary not only for plant cell architecture but also for membrane trafficking-mediated processes such as polar auxin transport, senescence, and cell death. In Arabidopsis, the inactivation of actin-based molecular motors, class XI myosins, affects the membrane trafficking and integrity of actin cytoskeleton, and thus causes defective plant growth and morphology, altered lifespan and reduced fertility. To evaluate the potential contribution of class XI myosins to the auxin response, senescence and cell death, we followed the flower and leaf development in the triple gene knockout mutant xi1 xi2 xik (3KO) and in rescued line stably expressing myosin XI-K:YFP (3KOR). Assessing the development of primary inflorescence shoots we found that the 3KO plants produced more axillary branches. Exploiting the auxin-dependent reporters DR5::GUS and IAA2::GUS, a significant reduction in auxin responsiveness was found throughout the development of the 3KO plants. Examination of the flower development of the plants stably expressing the auxin transporter PIN1::PIN1-GFP revealed partial loss of PIN1 polarization in developing 3KO pistils. Surprisingly, the stable expression of PIN1::PIN1-GFP significantly enhanced the semi-sterile phenotype of the 3KO plants. Further we investigated the localization of myosin XI-K:YFP in the 3KOR floral organs and revealed its expression pattern in floral primordia, developing pistils, and anther filaments. Interestingly, the XI-K:YFP and PIN1::PIN1-GFP shared partially overlapping but distinct expression patterns throughout floral development. Assessing the foliar development of the 3KO plants revealed increased rosette leaf production with signs of premature yellowing. Symptoms of the premature senescence correlated with massive loss of chlorophyll, increased cell death, early plasmolysis of epidermal cells, and strong up-regulation of the stress-inducible senescence-associated gene SAG13 in 3KO plants. Simultaneously, the reduced auxin responsiveness and premature leaf senescence were accompanied by significant anthocyanin accumulation in 3KO tissues. Collectively, our results provide genetic evidences that Arabidopsis class XI myosins arrange the flower morphogenesis and leaf longevity via contributing to auxin responses, leaf senescence, and cell death.

Introduction

The dynamics of actin cytoskeleton, including filament assembly, disassembly and reorganization are regulated by diverse actin binding proteins (ABPs). Various genetic studies have shown that down-regulation of ABPs affects the organization and dynamics of the actin arrays in plant cells (Kandasamy et al., 2005; Henty-Ridilla et al., 2014; Li et al., 2014; Wu et al., 2015; Tang et al., 2017). Actin-dependent molecular motors, myosins, make up a one group of the ABPs that transport endomembranes and other macromolecular cargoes along actin filaments (AFs) in plant cells. In Arabidopsis, the cell elongation is affected significantly when multiple class XI myosins are simultaneously eliminated (Prokhnevsky et al., 2008; Peremyslov et al., 2010; Ueda et al., 2010; Ojangu et al., 2012; Cai et al., 2014; Madison et al., 2015; Scheuring et al., 2016; Talts et al., 2016; Abu-Abied et al., 2018). In particular, the stature of myosin XI single, double, triple, quadruple and quintuple knockout plants is decreased progressively along with the reduction in cell size (Ojangu et al., 2007, 2012; Peremyslov et al., 2008, 2010; Prokhnevsky et al., 2008; Ueda et al., 2010; Park and Nebenführ, 2013; Cai et al., 2014; Madison et al., 2015; Abu-Abied et al., 2018). At a subcellular level, myosin XI inactivation in triple knockout mutant xi1 xi2 xik (3KO) results in reorientation of prominent actin bundles, reduced dynamic behavior of actin arrays, changes in membrane trafficking and deceleration of cytoplasmic streaming (Peremyslov et al., 2010; Ueda et al., 2010; Park and Nebenführ, 2013; Cai et al., 2014; Scheuring et al., 2016). Thus, in the 3KO, cytoplasmic streaming is virtually arrested, the elongation of epidermal and other cell types is decreased and therefore plant size and fertility is affected (Peremyslov et al., 2010; Ueda et al., 2010; Ojangu et al., 2012; Cai et al., 2014; Scheuring et al., 2016).

In addition to being essential for cell integrity, actin cytoskeleton contributes to processes such as polar auxin transport (PAT) (Nick et al., 2009; Zhu and Geisler, 2015; Wu et al., 2015; Zhu et al., 2016; Eggenberger et al., 2017; Huang et al., 2017) and regulation of programmed cell death (PCD) (Kandasamy et al., 2005; Thomas et al., 2006; Breeze et al., 2011; Keech, 2011; Smertenko and Franklin-Tong, 2011).

Various genetic and pharmacological studies have revealed tight interplay between auxin signaling and actin cytoskeleton. On the one hand, the patterning of actin arrays is modulated by auxin; on the other hand, auxin transport depends on the organization and dynamics of microfilaments (Zhu and Geisler, 2015). Auxin regulates both expansion and polarity of individual cells, as well as initiation and patterning of organs. Transient auxin concentration gradients underlie developmental processes such as meristem initiation, organ primordia formation, embryo morphogenesis, lateral root formation, as well as regulation of phyllotaxy and vascular tissue differentiation, photo- and gravitropic responses (Berleth and Sachs, 2001; Vanneste and Friml, 2009; Cardarelli and Cecchetti, 2014; Adamowski and Friml, 2015). Transient auxin concentration gradients result from local biosynthesis and polar cell-to-cell transport of the hormone. PAT is mediated by specific auxin uptake permeases and efflux carrier proteins that localize to plasma membrane in an asymmetric manner. PIN-FORMED (PIN) and ATP-binding cassette transporters/P-glycoprotein (ABCB/PGP) families are principal auxin efflux carriers whereas members of Auxin-Resistant 1/LIKE-AUX1 (AUX1/LAX) family are major auxin uptake carriers (Okada et al., 1991; Bennett et al., 1996; Gälweiler et al., 1998; Luschnig et al., 1998; Marchant et al., 1999; Swarup et al., 2001; Geisler and Murphy, 2006; Krecek et al., 2009; Zažímalová et al., 2010; Swarup and Péret, 2012; Adamowski and Friml, 2015). The regulation of PAT depends mainly on the action of auxin efflux carriers of the PIN and ABCB/PGP families (Geisler and Murphy, 2006; Vanneste and Friml, 2009; Adamowski and Friml, 2015). Both auxin influx permease AUX1, as well as efflux carriers (PINs and ABCBs) cycle between the plasma membrane and endosomal compartments (Geldner et al., 2001; Kleine-Vehn et al., 2006, 2011; Cho et al., 2007; Kleine-Vehn and Friml, 2008; Titapiwatanakun et al., 2009; Swarup and Péret, 2012; Cho and Cho, 2013; Wang et al., 2013). The crosstalk between auxin and actin is confirmed by findings showing that both localization and recycling of auxin importers and exporters depends partially on AFs (Geldner et al., 2001; Kleine-Vehn et al., 2006; Dhonukshe et al., 2008; Wu et al., 2015; Zhu and Geisler, 2015). However, the exact role of the actin cytoskeleton in PAT is still unresolved. It has been suggested that interactions between auxin and actin are not only dose- and time-dependent, but also species- and organ-dependent (Zhu and Geisler, 2015).

Several recent findings have shown that alterations in the cytoskeleton polymerization status are also critical for triggering PCD in plants (Smertenko et al., 2003; Thomas et al., 2006; Keech et al., 2010; Breeze et al., 2011; Keech, 2011; Smertenko and Franklin-Tong, 2011; Chang et al., 2015). Senescence is the terminal phase in organ development that involves a programmed degradation of cellular components. The resulting degradation products of senescing tissues are reused to support the growth of newly forming organs like leaves, roots, tubers, shoots, flowers, fruits, and seeds (Himelblau and Amasino, 2001; Maillard et al., 2015). Leaf senescence, the most well studied type of organ senescence in Arabidopsis, is characterized by the elevated expression of the senescence-associated genes (SAGs), early auxin-responsive small auxin-up RNA genes (SAURs), loss of chlorophyll, degradation of organelles, autolysis distinguished by fragmentation of the tonoplast and subsequent removal of the cytoplasm (Lohman et al., 1994; Weaver et al., 1998; Lin and Wu, 2004; Zentgraf et al., 2004; Balazadeh et al., 2008; Hou et al., 2013; Watanabe et al., 2013; Woo et al., 2013; Ren and Gray, 2015; Kim et al., 2016). Many SAGs encode proteins that drive breakdown of cellular components, such as short-chain alcohol dehydrogenase SAG13, and cysteine protease SAG12 (Weaver et al., 1997; Buchanan-Wollaston et al., 2003; Lin and Wu, 2004; Zentgraf et al., 2004; Watanabe et al., 2013; Woo et al., 2013; Kim et al., 2016). It has been demonstrated that the SAG12 is specifically induced by developmental senescence, and SAG13 by a range of senescence-inducing stress-treatments such as detachment, hormonal treatment, darkness, drought, wounding and pathogen attack. Therefore, it is proposed that SAG12 could be the marker for age-related developmental senescence, and SAG13 for stress-induced senescence or general cell-death (Schippers et al., 2007).

Auxin involvement in senescence has been observed to a much longer than the role of cytoskeleton (Hodge and Sacher, 1975; Noodén and Noodén, 1985; Lim et al., 2010; Kim et al., 2011; Ren and Gray, 2015; Cha et al., 2016). Nevertheless, the precise functions of auxin in leaf senescence remain unclear due to controversial results reporting either negative or positive role of auxin in leaf senescence regulation (Lim et al., 2010; Kim et al., 2011; Hou et al., 2013; Jibran et al., 2013; Woo et al., 2013; Khan et al., 2014; Ren and Gray, 2015; Cha et al., 2016). However, it has been suggested that auxin may promote leaf senescence through the expression of SAUR36 gene in Arabidopsis (Hou et al., 2013).

Although the auxin impact on secondary metabolism in plants is not well understood, numerous investigations point to the auxin’s role in modulating flavonoid biosynthesis (Buer and Muday, 2004; Besseau et al., 2007; Lewis et al., 2011; Kuhn et al., 2016). Flavonoids, in turn, are considered to be endogenous regulators of auxin efflux carriers, suggesting a crosstalk between auxin- and flavonoid-dependent processes (Murphy et al., 2000; Brown et al., 2001; Buer and Muday, 2004; Peer et al., 2004; Besseau et al., 2007; Peer and Murphy, 2007; Santelia et al., 2008; Zažímalová et al., 2010). Flavonoid biosynthesis produces a variety of distinct flavonoid subclasses, including anthocyanins, a group of pink, red, purple or blue pigments widely produced in plants (Harborne and Baxter, 1999). In Arabidopsis, anthocyanins accumulate in variable amounts in leaves and stems, depending on light intensity and nutrition (Holton and Cornish, 1995; Gou et al., 2011). Accumulation of anthocyanins and acceleration of senescence are also well documented under certain stress conditions including reduced nitrogen levels and high light intensity (Feild et al., 2001; Peng et al., 2008; Sekhon et al., 2012; Misyura et al., 2013).

Even though myosin and AFs act in concert, the potential role of myosins in auxin responses is starting to be revealed, but the one in senescence and cell death largely remains addressed. It has been reported previously that myosin 3KO roots show moderate unresponsiveness to exogenous auxin treatment, exhibiting partially insensitive vacuoles (Scheuring et al., 2016). Because both the AF architecture and overall actin dynamics are altered in myosin 3KO cells (Peremyslov et al., 2010; Ueda et al., 2010; Cai et al., 2014; Scheuring et al., 2016), it is likely that the mutant cells are less responsive to physiological and developmental stimuli such as auxin signaling. Very recent findings of Abu-Abied et al. (2018) showed that the altered root architecture of the 3KO plants was in correlation with the reduced auxin gradient, and partial loss of PIN1 polarization in the stele cells. These results provide first evidence that PAT, at least partially, could be myosin-mediated process in Arabidopsis.

In this study, we assessed the potential roles of class XI myosins in mediating auxin response and cell death during floral development and leaf senescence, respectively. We used well-characterized class XI myosin triple gene knockout mutant 3KO (Ojangu et al., 2012) as it has exhibited a prominent phenotype including stunted growth, partially impaired shoot development, and premature leaf yellowing, suggesting a possible connection between myosin function, auxin distribution and senescence signaling. Investigation of the genetically rescued 3KOR (Peremyslov et al., 2012) line confirmed that the observed defects have a myosin-dependent nature. Auxin-related processes in 3KO background were monitored through evaluating the activity of the auxin-responsive reporters DR5::GUS and IAA2::GUS, and comparing the expression patterns of the auxin efflux carrier PIN1::PIN1-GFP and of myosin XI-K:YFP in floral development. Senescence-related processes were analyzed by measuring the contents of chlorophylls and anthocyanins, and following the cell integrity, and patterning of AFs in senescent leaf cells. In addition, relative expression levels of auxin-responsive and senescence-related genes were evaluated. Collectively, our data imply that class XI myosins contribute significantly to auxin responses, stress-induced senescence, and cell death in Arabidopsis. At that, we provide first genetic evidences that actomyosin cytoskeleton mediates senescence-processes in Arabidopsis. Moreover, our results indicate that there is a mutual crosstalk between actomyosin cytoskeleton, auxin-regulated and senescence-dependent processes, and secondary metabolism.

Materials and Methods

Plant Material and Growth Conditions

Arabidopsis thaliana (ecotype Columbia-0) seeds of xi1 (Salk_022140; At1g17580), xi2 (Sail_632_D12; At5g43900), xik (Salk_067972; At5g20490) T-DNA mutant lines and PIN1::PIN1-GFP (ecotype Landsberg erecta, Ler) line (N23889) were obtained from the Nottingham Arabidopsis Stock Centre. Myosin mRNA levels of single mutant T-DNA lines were determined earlier by RT-qPCR (Talts et al., 2016). The myosin triple mutant line xi1 xi2 xik (3KO in this work) was published previously (Ojangu et al., 2012). PIN1::PIN1-GFP (Ler) and IAA2::GUS (Ler) lines were backcrossed four times to Columbia-0 (Columbia or Col in this work) prior to phenotypic analyses. The genetic background of the 3KO line transformed with the gene encoding YFP-tagged myosin XI-K, xi1 xi2 xik XI-K:YFP (3KOR in this work), was described earlier (Peremyslov et al., 2012). Reporter lines DR5::GUS (Col-0) and IAA2::GUS (Ler) were obtained from Malcolm Bennett’s lab, and seeds of the 35S::GFP-fABD2-GFP (Wang et al., 2008) line from Elison B. Blancaflor’s lab. The genes of DR5::GUS and IAA2::GUS auxin reporters, and the 35S::GFP-fABD2-GFP actin marker were introduced into the 3KO line by crossing the plant lines.

Vernalized seeds were held in water at 4°C for 1 day before sowing in the soil containing 50% (v/v) vermiculite. Plants were grown in growth chambers under 16 h light/8 h dark period at 22 ± 2°C and 60% of relative humidity. For the seedling analysis seeds were surface sterilized and grown on 0.5 × MS medium (Murashige and Skoog, 1962) supplemented with 1% sucrose in climate chambers as described above.

For 1-N-naphthylphthalamic acid (NPA) treatments, primary inflorescences were dipped twice (with 3-days interval) with 100 μm NPA (Sigma-Aldrich) and 0.01% (v/v) Silwet L-77 as adapted from Nemhauser et al. (2000). NPA was dissolved in dimethyl sulfoxide (DMSO), and mock treatments were performed with Milli-Q water containing 0.1% (v/v) DMSO and 0.01% Silwet L-77. For latrunculin B (LatB) (Abcam) treatments, 50 mg of 7-day-old seedlings grown on agar plates were incubated for 6 h in 5 ml of liquid MS medium supplemented with 0.5 μM LatB (dissolved in DMSO), and mock treatments were performed with liquid MS medium containing 0.025% (v/v) DMSO.

Quantitative Analysis of β-Glucuronidase (GUS) Activity in Plant Extracts

The activity of β-glucuronidase (GUS) can be determined in extracts of plant tissue using 4-methylumbelliferyl β-D-glucuronide (4-MUG) as a substrate. The 4-MUG fluorometric assay of plant extracts was performed to measure the GUS activity under the control of the DR5 promoter. For vegetative growth phase analysis 100 mg of 12-day-old seedlings and 23-day-old rosettes were collected and frozen. For generative growth phase analysis, 100 mg of primary inflorescence stems were collected when in height of 10–14 cm (floral transition stage) and 20–25 cm (silique formation stage). Protocol for quantitative GUS activity assay was adapted from Lewis and Muday (2009). Plant extract preparation was as follows: 150 μl GUS extraction buffer (50 mM sodium phosphate buffer pH 7.0; 10 mM EDTA pH 8.0; 0.1% SDS; 0.1% Triton X-100; 10 mM β-mercaptoethanol; 25 μg/μl PMSF) was added to the 100 mg of frozen tissues, homogenized at 28 Hz/min for 2 min (Qiagen TissueLyser), cell debris was removed by centrifugation (15 min 13,200 rpm at 4°C). The reaction was carried out by adding 50 μl of plant extract to the 450 μl pre-warmed reaction mixture (GUS extraction buffer supplemented with 2 mM 4-MUG), incubated 20 h at 37°C in darkness. The reaction was stopped by adding 100 μl of reaction mixture to 900 μl of ice-cold 0.2 M Na2CO3. 255 μl of the stopped reaction was loaded to the black 90-well microassay plates (Greiner) and protected from light. Released 4-methylumbelliferone (4-MU) (AppliChem) was excited at 340 nm and emission was measured at 492 nm using Tecan GENios Pro. Total protein content of extracts was measured following BioRad QuickStart 1 × Bradford protocol. Two to three experiments were performed with each tissue type, including 6–12 biological and two technical replicates per experiment. The GUS activity was expressed as picomoles of 4-MU produced per minute per milligram protein (pmol/min/mg).

Microscopy

The expression pattern analysis of Columbia and 3KO plants stably expressing PIN1::PIN1-GFP, XI-K:YFP or 35S::GFP-ABD2-GFP was performed with Carl Zeiss LSM 510 META confocal laser scanning microscope. GFP or YFP was excited at 488 nm with argon laser and fluorescence was detected with 505–550 nm band-pass filter; band-pass filter 575–615 IR was used for chlorophyll autofluorescence detection. Dissected flower parts, buds, and floral primordia were immersed in 50% glycerol or 95% perfluorodecalin (Sigma-Aldrich). Silicon spacer was applied in between glass slide and cover slip to prevent crushing tissues. Presented images show Z-stacks of confocal images combined into single image by maximal intensity projections using Zeiss LSM Image Browser software.

For histochemical analysis, DR5::GUS and IAA2::GUS plants were treated according to a standard protocol (Weigel and Glazebrook, 2002), immersed in Mowiol mounting medium, and analyzed under a light microscope. Images of DR5::GUS plants and trypan blue stained leaves were captured with digital camera Nikon D800E using a slide copying adapter ES-1. Images of IAA2::GUS inflorescences and flowers were captured with Zeiss SteREO Discovery.V8.

Scanning electron microscope (SEM) analysis was performed with Carl Zeiss EVO LS15. For flower architecture analysis in Arabidopsis, flower parts were dissected on double-sided tape using fine needle (27G) and SEM images of unfixed and uncoated flower tissues were captured using a reduced vacuum mode (100–200 Pa).

Adobe Photoshop CS6 was used to assemble photographs, indicate details, and measure the intensity of trypan blue staining on leaf photographs (as mean value of the blue channel). ImageJ2 software was used to analyze the PIN1-GFP fluorescence in single optical longitudinal sections (0.71 μm thick) of outer epidermal cells of stage 8 gynoecia. GFP fluorescence intensities of apical and lateral membranes of each cell was measured as quotients: mean gray values divided by values of selected areas. The ratios of the GFP fluorescence of apical membranes versus lateral membranes of each cell were used to evaluate the distribution of PIN1-GFP: the average ratio value of Columbia cells (3.3) was set as threshold to discriminate the cells with less polarized (<3.3) and more polarized (>3.3) PIN1-GFP distribution. The lower the ratio number was the less the fluorescence between apical and lateral membranes differentiated.

Quantification of Total Chlorophyll Content

For chlorophyll extraction, the fifth and sixth leaves from 23-day-old rosettes (before bolt formation) were dissected and frozen immediately. For dark-induced senescence, dissected leaves were placed in 30 ml of Milli-Q water and incubated for 3 days in growth chamber at darkness. One biological replicate contained six leaves – the fifth and sixth leaves from three different 23-day-old rosettes. Three independent experiments with six biological replicas in each were performed. Chlorophyll was extracted from the rosette leaves before and after dark-treatment as described previously (Hu et al., 2013; Misyura et al., 2013). Briefly, 1 ml of ice cold 80% acetone in 0.2 M Tris-buffer (pH 8.0) was added to frozen samples, homogenized with chrome-steal beads in TissueLyser 20 Hz/min for 1 min, and incubated 12 h in the dark at 4°C. Next day, samples were centrifuged 15 min 3000 rpm at 4°C. The extraction with 80% acetone solution was repeated three times, supernatants of every extraction were collected. 1 ml of extract was transferred to a disposable polymethyl methacrylate (PMMA) cuvette, and absorbance at 645 and 663 nm were measured. The total chlorophyll concentration in fresh weight (mg/g) were calculated using the equation:

Trypan Blue Staining

For trypan blue staining the fifth leaves from 21-day-old rosettes were selected. Fresh or dark treated leaves were boiled in lactophenol (10 ml of lactic acid, 10 ml of glycerol, 10 ml of liquid phenol, and 10 ml of distilled H2O) containing 10 mg of trypan blue for 1 min. Tissues were cleared in alcoholic lactophenol (2:1 95% ethanol:lactophenol) for 2 min, washed in 50% ethanol at room temperature, and stored in Milli-Q water. For the analysis, leaves were immersed in Mowiol mounting medium.

Quantification of Anthocyanin Content

Anthocyanin extraction from Arabidopsis seedlings was performed as described previously (Nakata and Ohme-Takagi, 2014). In brief, 100 mg of 7-day-old seedlings were frozen in liquid nitrogen, homogenized at 28 Hz/min for 2 min (TissueLyser), and suspended into five volumes of extraction buffer (45% methanol and 5% acetic acid). Cell debris was removed by centrifugation (20 min 13,200 rpm at room temperature). The relative anthocyanin content was calculated from the absorbance at 530 and 657 nm (Shimadzu Biospec Mini) using the equation:

RNA Extraction and Reverse Transcription – Quantitative Real-Time PCR (RT-qPCR)

Total RNA was extracted from 50 mg of plant material according to the method described by Oñate-Sánchez and Vicente-Carbajosa (2008). Buffer volumes were scaled up three times. Expression levels of AUX1, IAA2, PIN1, PIN3, PIN4, PIN7, SAUR36, SAG12, and SAG13 genes were analyzed. RNA was extracted from 7-day-old seedlings, 21-day-old leaves before and after dark-induced senescence, and mature inflorescences. cDNA was synthesized from 5 μg of DNase-treated RNA using Maxima Reverse Transcriptase (Thermo Scientific) and random hexamer primer. cDNAs were diluted twofold for qPCR. All RT-qPCR reactions were performed in 384-well plates on the LightCycler 480 instrument (Roche Applied Science). qPCR reactions were performed in duplicate and Cq values were averaged. Each 7 μl reaction contained 1.4 μl 5x HOT FIREPol® EvaGreen® qPCR Mix Plus (no ROX) (Solis Biodyne), 0.7 μl diluted cDNA and 3.5 pmol of each primer. qPCR conditions were as follows: initial denaturation at 95°C for 12 min, followed by 45 cycles of 95°C for 15 s, 59°C for 30 s and 72°C for 30 s. All primers used for qPCR experiments were designed for an annealing temperature of 60–62°C. Primers used for qPCR experiments are listed in Supplementary Table S1. Primers for reference genes were chosen as described previously (Czechowski et al., 2005; Supplementary Table S1). In all experiments, three reference genes were used for normalization: SAND, UBC, and expressed sequence EX70. Five to six biological replicates were analyzed in each experiment. Reference gene stability was analyzed using GeNorm M and coefficient of variation (CV) in qbasePLUS software (Hellemans et al., 2007). Statistical analysis (One-Way ANOVA) was performed with qbasePLUS software.

Results

The Activity of the Auxin-Responsive Reporters Is Reduced in 3KO Plants

Auxin signaling regulates all aspects of plant development, including determination of apical dominance during shoot growth, through finely tuned concentration gradients (Berleth and Sachs, 2001; Vernoux et al., 2010). Investigating the overall plant morphology, we noticed that the shoot development of 3KO plants varied significantly as mutant plants frequently displayed partially reduced apical dominance, and increased formation of axillary branches (Supplementary Figure S1). The average number of secondary inflorescences per cm of primary stems in 3KO plants was 1.8-fold, and the one of tertiary inflorescences 2.4-fold higher in comparison with Columbia control (Supplementary Figure S1). To assess if the developmental defects of 3KO plants have auxin-dependent nature, we used the auxin-responsive promoter-reporters DR5::GUS (β-glucuronidase), and IAA2::GUS which have been widely used as relevant markers to study endogenous auxin responses (Ulmasov et al., 1997; Shibasaki et al., 2009).

Fluorescent GUS assay using 4-methylumbelliferyl-β-D-galactopyranoside (4-MUG) fluorophore was used to measure DR5 promoter activity in plant extracts. Extracts prepared from 3KO seedlings stably expressing DR5::GUS (3KO DR5::GUS) showed 2.3-fold lower DR5 promoter activity when compared to the Columbia wild type control (Col DR5::GUS; Figure 1A). More specifically, histochemical GUS staining revealed reduced coloration in the roots of 3KO DR5::GUS plants (Figure 1B). GUS activity in rosette leaves of 3KO DR5::GUS plants was comparable with Columbia control (Figures 1C,D). Young 10–14 cm inflorescence stems and mature 20–25 cm inflorescence stems of the 3KO DR5::GUS plants showed 1.8-fold (Figure 1E) and 2.7-fold decrease (Figure 1G) in GUS activity, respectively, in comparison with Columbia control. Histochemical GUS assays revealed reduced GUS activity in the stems (Figure 1F) and siliques (Figure 1H) of 3KO DR5::GUS plants.

FIGURE 1

The IAA2::GUS reporter was used both for quantitative and histochemical analysis. Quantitative analysis of the IAA2::GUS activity in Columbia and 3KO extracts showed the similar reduction as the DR5::GUS reporter (data not shown). Histochemical staining of young inflorescences of Columbia and 3KO plants stably expressing IAA2::GUS showed similar staining patterns but different staining intensities (Figure 2). Throughout the flower development, the GUS staining in 3KO IAA2::GUS pistils (stigma, style and transmitting tract) (Figure 2B; st 9–13) and anthers (Figure 2B; st 10–11) was weaker than that in Columbia IAA2::GUS line. In addition, the GUS staining of 3KO stigmas at anthesis (Figure 2B; st 13) highlighted the retarded development of 3KO stigmas which was demonstrated previously (Ojangu et al., 2012).

FIGURE 2

We also analyzed the effect of NPA, the polar auxin transport inhibitor, on the IAA2::GUS expression. For this, inflorescences of Columbia IAA2::GUS and 3KO IAA2::GUS plants were dipped twice with 100 μM NPA. The histochemical examination of primary inflorescences revealed that NPA-treatment led to increased IAA2::GUS activity both in Columbia and 3KO inflorescences, when compared to DMSO-treated controls (Supplementary Figure S2). However, the responses to NPA were somewhat different in Columbia and 3KO. In NPA-treated Columbia, the strong GUS staining was spread all over the gynoecia, but it did not accumulate in valves of 3KO IAA2::GUS gynoecia (Supplementary Figure S2). Slightly weaker staining of pedicles, petals and sepals was noticeable in 3KO IAA2::GUS also (Supplementary Figure S2).

These results show that the simultaneous inactivation of three class XI myosins reduces the auxin responses of different Arabidopsis tissues.

Stable Expression of PIN1-GFP in 3KO Enhances Abnormalities in Flower Formation and Growth

In Arabidopsis, developmentally important auxin gradients are generated during PAT by modulating the organization and dynamics of actin cytoskeleton (Geldner et al., 2001; Ivakov and Persson, 2013; Zhu and Geisler, 2015). As the aberrant shoot development in 3KO plants correlated with the reduced auxin response, we further investigated if this was due to possible deviations in the distribution of auxin efflux carrier PIN1, since it plays an essential role in flower and inflorescence formation (Okada et al., 1991). For this, we examined the effect of stable expression of PIN1::PIN1-GFP in 3KO (3KO PIN1-GFP) and Columbia control plants (Col PIN1-GFP).

It is previously described that in 3KO plants, the fertility is slightly decreased as the pistil maturation (elongation of stigmas) partially delays (Ojangu et al., 2012). Interestingly, when examining the flowers at anthesis (stage 13/14) by using SEM, we found that the elongation of stamen filaments in parental 3KO line frequently delayed also (Figure 3L; white arrow). This indicated that during flower anthesis of the 3KO plants the rapid growth spurt of stamen filaments is partially retarded, and thus stamens often do not reach stigmas in time for proper pollination.

Even more surprising was the finding that previously described semi-sterile phenotype of 3KO plants was exacerbated in the 3KO PIN1-GFP line. SEM imaging revealed that the architecture of inflorescences, formation of floral primordia (Figure 4), and development of floral organs (Figure 3) were significantly affected in 3KO PIN1-GFP plants. The 3KO PIN1-GFP inflorescences displayed very irregular architecture being occasionally denser (Figure 4D and Supplementary Figure S3) or sparser (Figure 4E) in comparison with Columbia, 3KO or Col PIN1-GFP plants, indicating that the inflorescence meristem was partially disturbed. No matter if meristem produced more or less floral primordia (Figures 4D,E) their capability to develop into normal flowers was partially impaired. Some flowers showed completely arrested development (Figures 4D,E) and others exhibited a range of morphologies, from near-normal to severely deformed (Figure 3 and Supplementary Figure S3). Undeveloped flowers did not contain pistils, petals or anthers, and consisted only of one to three sepal-like structures (Figures 4D,E; white arrows). In extreme cases, the inflorescence meristem of the 3KO PIN1-GFP line produced only ∼20 near-normal flowers or floral buds. Thereafter, the emergence of new flowers stopped, since only needle-like floral primordia were formed (Figure 4E; asterisks). However, the most prevalent flower deformation was a significantly bent shape of pistils (Figures 3NP). Delayed apical closure of gynoecia (Figures 3M,N; white arrowheads) and pistils with decreased valves, and swollen style and stigma region (Figure 3M) were frequently observable during late stages of 3KO PIN1-GFP flower development. Besides serious aberrations, the flowers of 3KO PIN1-GFP plants exhibited also same deviations as parental 3KO line: retarded elongation of stamen filaments (anthesis stage 13/14) (Figure 3P; white arrow) and stigmatic papillae (pre-anthesis stage 11) (Figure 3N; asterisk). Modified architecture of 3KO PIN1-GFP inflorescences and flowers indicated that even the modest overexpression of PIN1 under native promoter may disturb auxin responses, and thus affect developmental decisions, when the actomyosin cytoskeleton is simultaneously affected.

FIGURE 3

FIGURE 4

Partial Loss of PIN1-GFP Polarization During 3KO Flower Development

Next, the distribution of auxin efflux carrier PIN1 in flower tissues of 3KO and Columbia plants was visualized with expression of the PIN1::PIN1-GFP. Examination of the PIN1-GFP distribution in flower tissues of Columbia plants revealed its presence in septum or valve margins of pistils, as well as in floral primordia (Figure 5A). The PIN1 patterning in the Columbia gynoecia was consistent with previously published data (Zúñiga-Mayo et al., 2014). The overall patterning of the PIN1-GFP in 3KO floral tissues was broadly the same with some exceptions. For instance, in heavily deformed 3KO pistils, such as valveless gynoecia, the PIN1-GFP signal was spread all over the gynoecium (Figure 5B; St 10; white bold arrow). Similarly, PIN-GFP pattern in developing 3KO floral primordia was occasionally aberrant (Supplementary Figure S4).

FIGURE 5

We examined the PIN1-GFP polarization in developing gynoecia at the cellular level because the significant pistil deformation was the most prominent phenotype of the 3KO PIN1::PIN1-GFP line. PIN1-GFP localization in longitudinal optical sections of the outer epidermal cell layer of stage 8 gynoecia was evaluated as ratio of apical membrane fluorescence to lateral membrane fluorescence of each cell. In outer epidermal cell layer of Columbia gynoecia, the PIN1-GFP was prominently apically localized (Figure 6A; arrowheads), the average ratio of GFP fluorescence of apical membranes versus lateral membranes was 3.3 ± 1.75 (n = 31). Epidermal cells of 3KO gynoecia displayed comparable PIN1-GFP fluorescence on both apical and lateral membranes (Figure 6B; asterisks), the average ratio of GFP fluorescence of apical membranes versus lateral membranes was 2.3 ± 1.33 [n = 42; p < 0.05 (Student’s t-test)]. To further validate the frequency of less and more polarized PIN1-GFP distribution of epidermal cells in Columbia and 3KO gynoecia, the ratio value of 3.3 of the GFP fluorescence was set as threshold (Figure 6C). The quantification showed that the ratio values less than 3.3 were observed in 86% of epidermal cells of developing 3KO PIN1-GFP gynoecia (Figure 6C) while in Col PIN1-GFP gynoecia, the frequency was 61% (Figure 6C). The ratio values higher than 3.3 were observed only in 14% of 3KO PIN1-GFP cells but in 39% of Col PIN1-GFP cells (Figure 6C). The gynoecium is the last organ to initiate from the floral meristem (Larsson et al., 2014), and apical domains such as the style and stigma are last structures which emerge during gynoecium development. The partial loss of PIN1 polarization in 3KO PIN1-GFP line is in accordance with severe pistil defects and indicated that the apical domains (style and stigma) of the developing gynoecium may not be sufficiently supplied with auxin in 3KO PIN1-GFP background. These results imply that class XI myosins, at least partially, contribute to the localization of PIN1 during floral development.

FIGURE 6

Myosin XI-K Is Expressed Throughout Floral Development

Despite the fact that class XI myosins share functional redundancy, it is known that the myosin XI-K plays important roles in such processes as membrane trafficking, cell expansion and division, plant growth, and fertility (Ojangu et al., 2007, 2012; Peremyslov et al., 2008, 2010, 2012; Avisar et al., 2012; Park and Nebenführ, 2013; Abu-Abied et al., 2018). Accordingly, to validate the XI-K role in flower development, we performed SEM analysis of the genetically rescued 3KOR line (Peremyslov et al., 2012), and found that the normal inflorescence development, including the elongation of stigmas and anther filaments, was restored in this plant line (Figures 3QT). Further, XI-K:YFP fluorescent signal in floral primordia, pistils and stamen filaments was examined (Figure 5C). The expression patterns of XI-K:YFP and PIN1::PIN1-GFP partially overlapped in the floral primordia, as well as in septum and valve margins of developing pistils (Figure 5C). At the same time, the XI-K:YFP showed distinct expression patterns in stigmas and anther filaments (Figure 5C), in comparison with PIN-GFP (Figures 5A,B). Both the expression pattern analysis, and flower architecture evaluation of the 3KOR line indicate that myosin XI-K contributes to the growth of floral organs, and thus to fertility.

The Senescence and Cell Death of 3KO Leaves Is Accelerated

Despite the delayed bolt formation and extended lifespan, described earlier (Peremyslov et al., 2010; Ojangu et al., 2012), we found that the 3KO plants displayed premature senescence of rosette leaves. Physical signs of early aging were particularly striking in the 12-day-old 3KO seedlings, which showed significant yellowing of cotyledons in comparison with Columbia or 3KOR line (Figure 7A). At the same time, 23-day-old mature rosettes (at bolt formation) of the 3KO plants produced 27% more leaves than those of Columbia control (Figure 7B). Given the very limited data on contributions of cytoskeleton to senescence and cell death in plant cells (Smertenko et al., 2003; Keech et al., 2010; Keech, 2011; Smertenko and Franklin-Tong, 2011), these observations prompted us to further explore myosin’s role in these processes.

FIGURE 7

First, to assess the senescence of 3KO plants, we measured chlorophyll content of detached rosette leaves before and after dark-induced senescence. Chlorophyll was extracted from the fifth and sixth leaves of 21-day-old rosettes. Quantification showed that the chlorophyll content in 3KO leaves before dark-treatment (Figure 8A; 3KO 0D) was only 19% of that in Columbia. After 3 days in darkness, the chlorophyll content of 3KO leaves was 9% of that of Columbia (Figure 8A; 3KO 3D). Photographs of detached leaves also show significant yellowing of the 3KO leaves after dark-treatment (Figure 8B; 3KO 3D).

FIGURE 8

Second, we employed trypan blue staining to distinguish the extent of cell death in rosette leaves before and after dark incubation (Figure 8C). Before dark incubation, at day 21 of rosette growth, the fifth leaves of both WT and 3KO plants were weakly stained with trypan blue (Figure 8C; Col 0D, 3KO 0D). Measurement of blue channel intensity of leaf pictures showed that after dark incubation, the trypan blue staining of Columbia leaves was 1.6 times, and the one of 3KO leaves 2.1 times more intense in comparison with untreated Columbia (Figure 8C; Col 3D, 3KO 3D), indicating the prevalence of dying cells in mutant leaves. In addition, the measurement of total protein content revealed that after dark-treatment, the total protein concentration of 3KO leaves was 2.8-fold lower, and the one of Columbia was 1.6-fold lower in comparison with untreated controls (3KO 0D, Col 0D), respectively (Supplementary Figure S5).

Third, to monitor the cell integrity and architecture of actin arrays in abaxial epidermal cells of the fifth rosette leaf’s petiole, we used stable expression of an AF tracer GFP-ABD2-GFP under control of the 35S promoter (Figure 9). Epidermal cells of leaf petioles were selected for examination as the cell growth, and AF organization defects of 3KO plants are most pronounced in longest cells, such as root hairs and petiole cells, as demonstrated by Peremyslov et al. (2010). Since 3KO leaf petioles are 50% shorter than those of Columbia we selected three areas of leaf petioles for examination: near the leaf blade, in the middle, and near the proximal end (Figure 9). The cell integrity and the organization of the AFs was compared between leaves of 21- and 28-day-old rosettes of 3KO and Columbia plants. At day 21 of rosette growth, the petiole cells were intact, and the GFP-ABD2-GFP decorated thick longitudinal cables in epidermal cells of Columbia leaf petioles (Figure 9A; Col 21D), and remarkably thin and prominently transverse filaments in 3KO cells (Figure 9B; 3KO 21D). At 28 days, GFP-ABD2-GFP labeled prominently longitudinal cables in Columbia cells (Figure 9A; Col 28D) and only traces of AFs in the 3KO petiole cells (Figure 9B; 3KO 28D). This was due to massive plasmolysis and cell shape deformation of 3KO petioles. In Columbia, less pronounced changes in AF organization and cell shape were found: 66% of petioles (n = 12) showed some plasmolysis at the area near the leaf blade, and only 25% of petioles showed massive plasmolysis (Figure 9A; 28D). This indicated that in 3KO plants, at day 28 of rosette growth, the cell death of older rosette leaves was significantly progressed in comparison with Columbia. Taken together, the results of decreased chlorophyll content, increased trypan blue staining, and premature plasmolysis in 3KO mutant show that the loss of integrity of actomyosin cytoskeleton induces premature senescence and thus cell death in Arabidopsis leaves.

FIGURE 9

The Anthocyanin Pigments Accumulate in 3KO Tissues

Eventually, when investigating the overall 3KO phenotype, we often observed an accumulation of purple pigments in various tissues including root-shoot junctions of seedlings (Figure 7A; red framed magnification), floral buds and basal parts of inflorescence stems (Figure 10A; white empty arrows). Such increased purple pigmentation of plant tissues is associated primarily with anthocyanin accumulation (Gou et al., 2011; Misyura et al., 2013; Mushtaq et al., 2016). To validate this assumption, the anthocyanin content of 7-day-old seedlings was determined spectrophotometrically (Figure 10B). It was found that anthocyanin content in the 3KO seedlings was 1.9-fold higher than that in Columbia control (Figure 10B). The anthocyanin accumulation rate in the 3KO rosettes and inflorescence stems was similar to that in seedlings (data not shown). This excessive anthocyanin accumulation is likely an accompanying effect related to reduced auxin response and accelerated senescence in the 3KO plants.

FIGURE 10

Expression of the Auxin-Responsive and Senescence-Associated Genes Is Altered in 3KO

The main mechanism by which auxin- and senescence-responses are converted into cellular responses is via changes in transcription (Zentgraf et al., 2004; Paponov et al., 2008; Kim et al., 2016). To examine, if the elimination of multiple class XI myosins or the stable expression of PIN1::PIN1-GFP affect the regulation of the auxin-responsive and senescence-associated genes, we measured the relative expression levels of the selected mRNAs in seedlings, inflorescences, and rosette leaves of the Columbia, 3KO, 3KO PIN1-GFP, and Col PIN1-GFP plants using RT-qPCR. To assess auxin-related processes, we followed the expressions of auxin importer (AUX1) and exporter (PIN1, PIN3, PIN4, PIN7) genes whose expressions are related to flower development (Krecek et al., 2009; Lampugnani et al., 2013). We also measured the level of IAA2 as its expression is closely related to endogenous auxin (Shibasaki et al., 2009). To distinguish stress- and age-induced senescence-related processes, we followed the relative expressions of SAG13 and SAG12 genes, respectively (Swartzberg et al., 2006; Schippers et al., 2007; Hou et al., 2013). To assess mutual influences between auxin- and senescence-related processes the level of auxin-responsive SAUR36 was measured.

First, the expression levels of selected auxin- and senescence-associated genes in the 3KO seedlings, leaves and inflorescences were compared with Columbia control. In the 7-day-old 3KO seedlings, early senescence-associated gene SAG13 was 94-fold up-regulated when compared to the Columbia control (Figure 11A). The down-regulation of PIN7 in mutant seedlings was moderate, although statistically significant (Figure 11A).

FIGURE 11

The relative levels of selected genes expressed in rosette leaves of 3KO and Columbia plants were compared before and after dark treatment. Before dark-treatment, at 21 days of growth, 3KO rosette leaves (3KO 0D) showed a 34-fold up-regulation of SAG13 expression, when compared to untreated Columbia control (Col 0D) (Figure 11B). The expression of SAUR36 in untreated 3KO (3KO 0D) leaves was twofold increased and those of AUX1 and PIN1 were twofold decreased, in comparison with untreated control leaves (Col 0D) (Figure 11B). The dark-treatment of detached leaves significantly affected the expression levels of all selected auxin-responsive and senescence-associated genes both in 3KO (3KO 3D) and Columbia (Col 3D) plants when compared to untreated control (Col 0D) (Figure 11B). However, the changes in 3KO background were more pronounced than those in Columbia. Down-regulation of AUX1 (7-fold), IAA2 (10-fold), and PIN7 (110-fold) in dark-treated 3KO leaves was twice as high as in Columbia (Col 3D) (Figure 11B). It is noteworthy that the up-regulation of SAG13 in dark-treated Columbia (Col 3D) was comparable with untreated 3KO leaves (3KO 0D), whereas in dark-treated mutant (3KO 3D), the SAG13 was already 89-fold up-regulated (Figure 11B).

In 3KO inflorescences, levels of AUX1 and SAG12 were down-regulated twofold and those of IAA2 and PIN7 moderately, in comparison with Columbia control (Figure 11C). These results show that the onset of leaf senescence of 3KO plants is initiated at transcript level already in the 7-day-old seedlings, as demonstrated by the conspicuously high expression level of the SAG13 gene. The unchanged levels of SAG12 mRNA in seedlings and leaves indicates that the premature senescence of 3KO leaves most probably does not result from the activation of a developmental senescence program since the strong expression of SAG13 suggests the stress-induced senescence and cell death. The relative expression levels of auxin-responsive AUX1 and PIN1 in 3KO were affected the most in leaves and inflorescences, implying that defective development of 3KO shoots could be influenced by reduced auxin responses.

The PIN1 is the principal member of the PIN family involved in development of aerial organs and the lead player in floral development (Gälweiler et al., 1998; Benková et al., 2003; Scarpella et al., 2006; Adamowski and Friml, 2015). Therefore, even the modest overexpression of PIN1 under native promoter is thought to affect plant phenotype. Therefore, the levels of auxin- and senescence-responsive genes in PIN1::PIN1-GFP expressing lines were compared with untransformed Columbia control. As we expected, both in 3KO PIN1-GFP, as well as in Col PIN1-GFP seedlings the expression of PIN1 was about 1.5-fold up-regulated in comparison with Columbia control (Figure 11A). The level of IAA2 was three–fourfold down-regulated in seedlings of PIN1::PIN1-GFP expressing Columbia and 3KO (Figure 11A). In 3KO PIN1-GFP seedlings, the level of SAG13 was 92-fold up-regulated, similarly like in parental 3KO line (Figure 11A). Likewise, the level of SAUR36 was also twofold up-regulated in 3KO PIN1-GFP seedlings compared to untransformed Columbia and 3KO (Figure 11A). A moderate down-regulation in levels of AUX1, PIN3, and PIN7 in 3KO PIN1-GFP, and PIN3 in Col PIN1-GFP seedlings was also detectable (Figure 11A).

In 3KO PIN1-GFP inflorescences, the expression of PIN1 was up-regulated twofold, whereas in Col PIN1-GFP it was comparable with Columbia control (Figure 11C). The level of IAA2 both in 3KO PIN1-GFP and Col PIN1-GFP inflorescences was down-regulated six and threefold, respectively (Figure 11C). Relative expression levels of AUX1 and PIN7 were down-regulated moderately in both PIN1-GFP expressing lines (Figure 11C). Senescence-associated genes, SAG12 and SAG13, were four and twofold down-regulated in 3KO PIN1-GFP inflorescences, whereas down-regulation of SAG12 in Col PIN1-GFP shoots was marginal (Figure 11C). These data show that even the mild overexpression of the major non-redundant member of the auxin efflux carrier family, PIN1, inevitably affects the expression of auxin-responsive and senescence-associated genes not only in 3KO, but in Columbia background too. These results indicate also that there is a mutual crosstalk between myosin-mediated transport, auxin-signaling and senescence-related processes.

Third, we investigated how the treatment with AF destabilizing drug latrunculin B (LatB) affects the expression levels of AUX1, IAA2, PIN1, SAUR36, and SAG13 in 7-day-old Columbia seedlings in comparison with untreated Columbia and 3KOR seedlings (Figure 11D). Six hours after the application of 0.5 μM LatB the 1.9-fold up-regulation of SAG13, and moderate down-regulation of AUX1 in Columbia seedlings was detectable (Figure 11D). In 3KOR seedlings, the levels of AUX1, IAA2, PIN1, SAUR36, and SAG13 were comparable with untreated Columbia control (Figure 11D). These data show that the latrunculin B-mediated disruption of AFs activates the expression of stress-inducible SAG13 gene, indicating that the actin cytoskeleton may be necessary for delivering stress-responses in plant cells.

Discussion

Our results show that the three Arabidopsis class XI myosins, XI-1, XI-2, and XI-K, contribute to the auxin responses and cell death, and thereby affect developmental decisions during flower growth and leaf senescence. Using triple mutant line 3KO (Ojangu et al., 2012) and complemented line 3KOR (Peremyslov et al., 2012), we demonstrate that simultaneous depletion of these three myosins affects responsiveness of the auxin-dependent promoters and polarization of the PIN1 auxin efflux carrier, causes premature onset of senescence and cell death in leaves, elevates accumulation of anthocyanins, and changes the expression levels of genes related to these processes. We also show that the stable expression of myosin XI-K rescues the decreased fertility and prematurely senescent phenotype in 3KO background. This result implements that myosin XI-K plays important role not only in driving vegetative plant growth (Ojangu et al., 2007, 2012; Peremyslov et al., 2012) and gravitropic response (Talts et al., 2016), but contributes also to floral development and cell death.

Previous works on multiple gene knockout mutants have identified several myosin XI functions in plant development. In particular, the best studied phenotypes of 3KO plants included stunted growth, delayed bolting, incomplete development of stigmas and reduced fertility (Peremyslov et al., 2010; Ojangu et al., 2012). These phenotypic defects are explained with severe changes at the cellular level: disorganized and more static actin cytoskeleton, reduced membrane trafficking, nearly arrested cytoplasmic streaming (Peremyslov et al., 2010; Ueda et al., 2010; Cai et al., 2014), reduced sensitivity of vacuoles to exogenous auxin, and partial loss of PIN1 polarization in root cells (Scheuring et al., 2016; Abu-Abied et al., 2018).

First, we used flower development of 3KO plants as a model for examining myosin-dependent auxin-responsive processes in Arabidopsis. We show that the apical dominance and branching architecture of 3KO shoots is partially affected as the mutant plants produce more axillary branches on primary inflorescence stem. Investigating inflorescence development, four traits are usually evaluated: bolting time, length of the reproductive phase, number of rosette leaves at bolting, number of axillary branches and fruits (Ungerer et al., 2002; Pouteau and Albertini, 2009). In myosin 3KO mutant, all these developmental aspects are partly affected: bolting time delays, reproductive phase is expanded, more rosette leaves are formed, branching architecture is affected, and silique size is decreased. The branching architecture of inflorescence shoots is mainly regulated by auxin, cytokinin, and strigolactone, which control initiation and outgrowth of axillary meristem (Domagalska and Leyser, 2011). Basipetal auxin transport from shoot apex toward the base suppresses axillary meristem outgrowth, and leads to apical dominance (Davies et al., 1966). In Arabidopsis and tomato, the axillary meristem initiation is characterized by preparative auxin depletion and the subsequent meristem emergence by a local auxin accumulation (Wang et al., 2014). It has been showed that local auxin gradients necessary during phyllotactic patterning in leaf and inflorescence meristems are regulated by auxin importers together with PIN1 exporter (Bainbridge et al., 2008). Consistent with this, our results show that the auxin responsiveness is reduced throughout the development of the 3KO plants, with the most pronounced effects in seedlings and inflorescence stems. In addition, as the PIN1 polarization in epidermal cells of developing 3KO gynoecia is partially disturbed, and two major auxin transporter genes, AUX1 and PIN1, are down-regulated in 3KO leaves and inflorescences, a correlation between these processes can be found, indicating that PAT-related processes are disturbed in 3KO background. Interestingly, our results complement recently published data which demonstrated that the production of lateral and adventitious roots in 3KO plants was increased (Abu-Abied et al., 2018). The abnormal formation of the lateral roots was related to changed auxin gradient, which was in correlation with partial loss of PIN1 polarization in stele cells of 3KO plants. Moreover, authors revealed the myosin XI-K role in cell division in both the root and shoot meristem (Abu-Abied et al., 2018). The reduced auxin responsiveness and loss of myosin XI-K function in 3KO meristem was used to explain the altered branching architecture of mutant roots (Abu-Abied et al., 2018). Although these two studies (current work; Abu-Abied et al., 2018) used different models (flowers and leaves versus roots) and different experimental approaches, both highlight myosin-auxin connections.

In this work, we demonstrate for the first time that the myosin XI-K:YFP is expressed in floral primordia and in developing flowers, indicating its role in floral development. Moreover, correlation between the expression pattern of XI-K:YFP and the delayed elongation of stigmas and anther filaments in 3KO can be found, though we do not exclude the roles of myosin XI-1 and XI-2 in these processes.

It is broadly assumed that development of the stamens and gynoecia, the most modified floral organs, require enhanced production of hormones, including auxin (Okada et al., 1991; Sessions et al., 1997; Nemhauser et al., 2000; Aloni et al., 2006; Cheng et al., 2006; Cecchetti et al., 2008; Sundberg and Østergaard, 2009; Hawkins and Liu, 2014). It is proposed that the finely tuned auxin gradient is necessary for gynoecium patterning in general and style and stigma development in particular (Nemhauser et al., 2000). To complete the development of valves, style and stigma, the medial, lateral and apical domains of the growing gynoecium are provided with auxin from the base. This has been demonstrated by Larsson et al. (2014), who showed that the PIN1-GFP localization in outer epidermal cells of stage 7 gynoecia was prominently apical. Interestingly, we found that stable expression of the PIN1::PIN1-GFP augmented a semi-sterile phenotype of the 3KO plants (Peremyslov et al., 2010; Ojangu et al., 2012). Our results indicate that the ability of 3KO PIN1-GFP inflorescence meristem to produce flowers with determined number and shape of inner whorl organs is affected to a variable degree, whereas the gynoecium development is affected the most. Similar gynoecium phenotypes have been described upon direct genetic or pharmaceutical disruption of auxin signaling (Sessions et al., 1997; Nemhauser et al., 2000; Cheng et al., 2006; Zúñiga-Mayo et al., 2014). The gynoecium is the last organ which emerges from the floral meristem (Larsson et al., 2014), and the style and stigma are last structures which emerge from the growing gynoecium. The partial loss of PIN1 polarization in developing 3KO PIN1-GFP gynoecia may indicate that the auxin transport from the base may not be sufficient to complete the growth of style, stigma or valves of 3KO PIN1-GFP pistils correctly. Here we demonstrate also that the late phase of stamen development, pre-anthesis filament elongation, often delays both in the parental line as well as in 3KO PIN1-GFP line. Our results indicate that the defective floral development of the 3KO PIN1-GFP line could be due to inability of the disrupted actomyosin cytoskeleton to properly allocate an excess of PIN1 which in turn, via feedback signaling pathway, may affect the expression of auxin-responsive genes. All these results together indicate the possible cooperation between myosin-mediated trafficking and auxin responses in floral development.

Second, we used leaf development of 3KO plants as a model for examining myosin-dependent senescence-associated processes in Arabidopsis. It is well documented that late flowering plants have more rosette leaves, and that there is a strong correlation between leaf number and bolting time (Pouteau and Albertini, 2009; Schmalenbach et al., 2014). In accordance with this, rosettes of 3KO plants create more leaves, but at the same time show premature signs of leaf senescence. Premature senescence in 3KO leaves is confirmed by massive loss of chlorophyll, enhanced cell death, premature rupture of petiole epidermal cells, and abnormal accumulation of anthocyanins. Our finding that the SAG13 is significantly up-regulated prior to observable signs of leaf yellowing both in young 3KO seedlings as well as in rosette leaves, indicates that the premature senescence in this mutant line is highly accelerated at transcript level. In addition, the expression of late senescence-associated gene SAG12 is not different in the 3KO seedlings and leaves relative to the Columbia. Moreover, the accumulation of anthocyanins in 3KO plants is in accordance with the fact that increased flavonoid production (e.g., anthocyanins) has been also associated with stress responses (Solfanelli et al., 2006). As we detect the anthocyanin accumulation, and SAG13 up-regulation already in young 3KO seedlings, this supports our assumption that the early leaf senescence of 3KO mutants is not typical developmental senescence but could be caused by cellular stresses such as changes in membrane trafficking, AF rearrangement, and auxin responses. For example, when wild type Arabidopsis seedlings are suffering from salt and methyl jasmonate co-stress, the leaf senescence is prematurely initiated as the SAG13 of transcript level is significantly up-regulated (Chen et al., 2017). In Arabidopsis accelerated-cell-death11 (acd11) mutant, the SAG13 mRNA is strongly expressed, whereas SAG12 mRNA does not accumulate. The authors concluded that the cell death in acd11 does not result from activation of a senescence program, but from PCD as suggested by conspicuously strong SAG13 expression (Brodersen et al., 2002). It is tempting to speculate that in 3KO cells, the loss of integrity of actomyosin cytoskeleton activates the stress-signaling pathway, and as a consequence of this the senescence and cell death are prematurely initiated.

It has been shown that the integrity of actin cytoskeleton contributes to senescence and cell death through reorganization of the endoplasmic reticulum (ER) network (Boevink et al., 1998; Kasaras et al., 2012; Xu et al., 2012; Carvalho et al., 2014; Chang et al., 2015) and transcription factor activation (Hao and August, 2005; Breeze et al., 2011; Smertenko and Franklin-Tong, 2011). Smertenko et al. (2003) suggested that the filamentous actin is crucial for the developmentally programmed cell death in the Norway spruce embryos as AFs disappeared significantly later than the microtubules. It has been proposed also that during leaf senescence, a functional actin cytoskeleton could be essential for maintaining the primary metabolism until the cell death (Keech, 2011). Because the integrity of actin arrays, and thus membrane trafficking is affected in 3KO epidermal cells (Peremyslov et al., 2010; Cai et al., 2014), a causal connection between these processes can be suggested. Myosin’s potential role in senescence is supported by the facts that ER-streaming is dramatically suppressed in 3KO plants (Ueda et al., 2010), and a portion of the XI-K:YFP is aligned and co-fractionated with a motile ER subdomain (Peremyslov et al., 2012). Intriguingly, myosins XI-K and XI-2 that are inactivated in 3KO are not only broadly expressed throughout plant tissues, but show particularly high expression level in the senescent leaves too (Peremyslov et al., 2011).

In general, the role of auxin in regulating leaf senescence still remains elusive. It has been implied that auxin may promote leaf senescence via SAUR36 gene since the auxin-inducible transcript accumulation promotes premature senescence of young leaves, whereas in saur36 null mutant plants senescence delays substantially (Hou et al., 2013). Our findings show that the premature leaf senescence of 3KO leaves is indeed accompanied by up-regulation of SAUR36, although it is not clear how the early aging is influenced by reduced auxin responses in this mutant. Up-regulation of SAG13 and SAUR36 during leaf senescence, regardless of whether it is induced naturally or by darkness suggests that the senescence mechanism under these two different conditions share common features. It is well known that the dark-induced senescence of detached leaves is radical intervention as major changes in leaf physiology inevitably affect the gene expression. Although, the dark-treatment of detached 3KO leaves does not reflect the situation of natural senescence, the results obtained from these experiments, provide additional clues that myosins contribute to the auxin- and senescence-responsive processes, and indicate that these processes are mutually interconnected.

Accumulation of anthocyanins in plant tissues is related both with auxin signaling and senescence (Brown et al., 2001; Feild et al., 2001; Buer and Muday, 2004; Besseau et al., 2007; Schippers et al., 2007; Falcone Ferreyra et al., 2012). In particular, auxin transport is elevated in inflorescences, hypocotyls, and roots of plants of the flavonoids-deficient mutants (Murphy et al., 2000; Brown et al., 2001; Lewis et al., 2011). During senescence, anthocyanins reduce the risk of photo-oxidative damage in leaf cells and thereby help retrieve nutrients from senescing tissues (Feild et al., 2001; Falcone Ferreyra et al., 2012). We propose that the accumulation of anthocyanins in 3KO tissues is an accompanying effect influenced both by reduced auxin responses as well as by premature onset of senescence. Collectively, this work highlights the importance of actomyosin cytoskeleton in auxin responsiveness and senescence, and leads to further questions as to how reorganization of metabolism is achieved in cells having a drastically affected actomyosin cytoskeleton. The next critical step in understanding these networks is investigation of the mechanistic contributions of myosins to the function of these networks. It is anticipated that such contributions could involve specific role of myosin-dependent transport in PIN and AUX targeting, overall reduction in cytoplasmic streaming that could affect auxin transport or auxin diffusion, as well as myosin-dependent changes in AF organization. Whether these myosins contribute to auxin responses through shaping actin cytoskeleton or via reorganization of the ER remains to be elucidated. Important aspect of the future studies will be identifying the myosin cargoes that affect auxin signaling and senescence. These cargoes could include MyoB myosin receptors that appear to drive cytoplasmic streaming (Peremyslov et al., 2013, 2015) or newly identified myosin adaptors of MadA and MadB families that presumably mediate more specialized myosin-dependent processes (Kurth et al., 2017).

Taken together, our results provide the genetic evidence that the integrity of actomyosin cytoskeleton, signaling of auxin and senescence, as well as secondary metabolism are functionally intertwined in a finely tuned network as illustrated in a tentative model (Figure 12). According to this, remodeling of actomyosin cytoskeleton affects not only distribution of auxin exporter PIN1, but also the expression of auxin-responsive and senescence-associated genes. Most likely, the levels of auxin-responsive genes are influenced indirectly, through feedback signaling mechanism, mediated by changes in auxin responses. The actomyosin cytoskeleton-mediated signaling of premature leaf senescence involves the control of expression of stress-inducible gene SAG13. Senescence and auxin signaling could be mutually regulated by secondary metabolism (e.g., anthocyanin accumulation) and by expression of auxin-responsive SAUR genes (e.g., SAUR36). All these processes are mutually interconnected, and if one process is unbalanced the others are affected too. It is tempting to speculate that changes in auxin and senescence responses influence reproduction of 3KO plants in opposite directions (Figure 12). From the one hand, the reduced auxin response could disturb flower development; on the other hand, it seems that the premature leaf senescence could constitute a rescue mechanism for supporting the production of inflorescence biomass in conditions where flower development has been compromised.

FIGURE 12

Statements

Author contributions

E-LO, BI, KTal, HP, ET, and VD contributed to design of the study and interpretation of data. E-LO, BI, KTal, KTan, and EI conducted experiments and analyzed data. KTan performed anthocyanin and chlorophyll measurements, and trypan blue stainings. BI performed RT-qPCR. EI performed histochemical GUS stainings. KTal helped with protocol optimization for quantitative MUG assay. HP helped with SEM and confocal imaging, image processing, and data analysis. VD generated 3KOR line. E-LO, BI, KTal, HP, VD, and ET drafted the manuscript and revised it critically. All authors read and approved the final manuscript, and they agreed to be accountable for all aspects of the work.

Funding

This research was supported by institutional research grant IUT 193 from the Estonian Ministry of Education and Research, and Estonian Research Council grant ETF8604. Each of the funding bodies granted the funds based on a research proposal. They had no influence over the experimental design, data analysis, interpretation, or writing the manuscript.

Acknowledgments

We thank prof. Malcolm Bennett from the University of Nottingham for providing DR5::GUS line and prof. Elison B. Blancaflor from The Samuel Roberts Noble Foundation for providing 35S::GFP-ABD2-GFP seeds. We thank Cecilia Sarmiento for thoroughly editing manuscript; Kairi Kärblane for helping with Adobe Photoshop CS6; Signe Nõu for taking care of plants. We sincerely thank the reviewers for their detailed comments, suggestions and constructive criticism for the manuscript.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2018.01570/full#supplementary-material

References

  • 1

    Abu-AbiedM.BelausovE.HagayS.PeremyslovV.DoljaV.SadotE. (2018). Myosin XI-K is involved in root organogenesis, polar auxin transport and cell division.J. Exp. Bot.6928692881. 10.1093/jxb/ery112

  • 2

    AdamowskiM.FrimlJ. (2015). PIN-dependent auxin transport: action, regulation, and evolution.Plant Cell272032. 10.1105/tpc.114.134874

  • 3

    AloniR.AloniE.LanghansM.UllrichC. I. (2006). Role of auxin in regulating Arabidopsis flower development.Planta223315328. 10.1007/s00425-005-0088-9

  • 4

    AvisarD.Abu-AbiedM.BelausovE.SadotE. (2012). Myosin XIK is a major player in cytoplasm dynamics and is regulated by two amino acids in its tail.J. Exp. Bot.63241249. 10.1093/jxb/err265

  • 5

    BainbridgeK.Guyomarc’hS.BayerE.SwarupR.BennettM.MandelT.et al (2008). Auxin influx carriers stabilize phyllotactic patterning.Genes Dev.22810823. 10.1101/gad.462608

  • 6

    BalazadehS.ParlitzS.Mueller-RoeberB.MeyerR. C. (2008). Natural developmental variations in leaf and plant senescence in Arabidopsis thaliana.Plant Biol.10(Suppl. 1), 136147. 10.1111/j.1438-8677.2008.00108.x

  • 7

    BenkováE.MichniewiczM.SauerM.TeichmannT.SeifertováD.JürgensG.et al (2003). Local, efflux-dependent auxin gradients as a common module for plant organ formation.Cell115591602. 10.1016/S0092-8674(03)00924-3

  • 8

    BennettM. J.MarchantA.GreenH. G.MayS. T.WardS. P.MillnerP. A.et al (1996). Arabidopsis AUX1 gene: a permease-like regulator of root gravitropism.Science273948950. 10.1126/science.273.5277.948

  • 9

    BerlethT.SachsT. (2001). Plant morphogenesis: long-distance coordination and local patterning.Curr. Opin. Plant Biol.45762. 10.1016/S1369-5266(00)00136-9

  • 10

    BesseauS.HoffmannL.GeoffroyP.LapierreC.PolletB.LegrandM. (2007). Flavonoid accumulation in Arabidopsis repressed in lignin synthesis affects auxin transport and plant growth.Plant Cell19148162. 10.1105/tpc.106.044495

  • 11

    BoevinkP.OparkaK.CruzS. S.MartinB.BetteridgeA.HawesC. (1998). Stacks on tracks: the plant Golgi apparatus traffics on an actin/ER network.Plant J.15441447. 10.1046/j.1365-313X.1998.00208.x

  • 12

    BreezeE.HarrisonE.McHattieS.HughesL.HickmanR.HillC.et al (2011). High-resolution temporal profiling of transcripts during Arabidopsis leaf senescence reveals a distinct chronology of processes and regulation.Plant Cell23873894. 10.1105/tpc.111.083345

  • 13

    BrodersenP.PetersenM.PikeH. M.OlszakB.SkovS.OdumN.et al (2002). Knockout of Arabidopsis accelerated-cell-death11 encoding a sphingosine transfer protein causes activation of programmed cell death and defense.Genes Dev.16490502. 10.1101/gad.218202

  • 14

    BrownD. E.RashotteA. M.MurphyA. S.NormanlyJ.TagueB. W.PeerW. A.et al (2001). Flavonoids act as negative regulators of auxin transport in vivo in Arabidopsis.Plant Physiol.126524535. 10.1104/pp.126.2.524

  • 15

    Buchanan-WollastonV.EarlS.HarrisonE.MathasE.NavabpourS.PageT.et al (2003). The molecular analysis of leaf senescence – A genomics approach.Plant Biotechnol. J.1322. 10.1046/j.1467-7652.2003.00004.x

  • 16

    BuerC. S.MudayG. K. (2004). The transparent testa4 mutation prevents flavonoid synthesis and alters auxin transport and the response of Arabidopsis roots to gravity and light.Plant Cell1611911205. 10.1105/tpc.020313

  • 17

    CaiC.Henty-RidillaJ. L.SzymanskiD. B.StaigerC. J. (2014). Arabidopsis myosin XI: a motor rules the tracks.Plant Physiol.16613591370. 10.1104/pp.114.244335

  • 18

    CardarelliM.CecchettiV. (2014). Auxin polar transport in stamen formation and development: how many actors?Front. Plant Sci.5:333. 10.3389/fpls.2014.00333

  • 19

    CarvalhoH. H.SilvaP. A.MendesG. C.BrustoliniO. J. B.PimentaM. R.GouveiaB. C.et al (2014). The endoplasmic reticulum binding protein BiP displays dual function in modulating cell death events.Plant Physiol.164654670. 10.1104/pp.113.231928

  • 20

    CecchettiV.AltamuraM. M.FalascaG.CostantinoP.CardarelliM. (2008). Auxin regulates Arabidopsis anther dehiscence, pollen maturation, and filament elongation.Plant Cell2017601774. 10.1105/tpc.107.057570

  • 21

    ChaJ.-Y.KimM. R.JungI. J.KangS. B.ParkH. J.KimM. G.et al (2016). The thiol reductase activity of YUCCA6 mediates delayed leaf senescence by regulating genes involved in auxin redistribution.Front. Plant Sci.7:626. 10.3389/fpls.2016.00626

  • 22

    ChangX.RiemannM.LiuQ.NickP. (2015). Actin as deathly switch? How auxin can suppress cell-death related defence.PLoS One10:e0125498. 10.1371/journal.pone.0125498

  • 23

    ChenY.WangY.HuangJ.ZhengC.CaiC.WangQ.et al (2017). Salt and methyl jasmonate aggravate growth inhibition and senescence in Arabidopsis seedlings via the JA signaling pathway.Plant Sci.26119. 10.1016/j.plantsci.2017.05.005

  • 24

    ChengY.DaiX.ZhaoY. (2006). Auxin biosynthesis by the YUCCA flavin monooxygenases controls the formation of floral organs and vascular tissues in Arabidopsis.Genes Dev.2017901799. 10.1101/gad.1415106

  • 25

    ChoM.ChoH.-T. (2013). The function of ABCB transporters in auxin transport.Plant Signal. Behav.8:e22990. 10.4161/psb.22990

  • 26

    ChoM.LeeS. H.ChoH.-T. (2007). P-glycoprotein4 displays auxin efflux transporter–like action in Arabidopsis root hair cells and tobacco cells.Plant Cell1939303943. 10.1105/tpc.107.054288

  • 27

    CzechowskiT.StittM.AltmannT.UdvardiM. K.ScheibleW.-R. (2005). Genome-wide identification and testing of superior reference genes for transcript normalization in Arabidopsis.Plant Physiol.139517. 10.1104/pp.105.063743

  • 28

    DaviesC. R.SethA. K.WareingP. F. (1966). Auxin and kinetin interaction in apical dominance.Science151468469. 10.1126/science.151.3709.468

  • 29

    DhonuksheP.GrigorievI.FischerR.TominagaM.RobinsonD. G.HasekJ.et al (2008). Auxin transport inhibitors impair vesicle motility and actin cytoskeleton dynamics in diverse eukaryotes.Proc. Natl. Acad. Sci. U.S.A.10544894494. 10.1073/pnas.0711414105

  • 30

    DomagalskaM. A.LeyserO. (2011). Signal integration in the control of shoot branching.Nat. Rev. Mol. Cell Biol.12211221. 10.1038/nrm3088

  • 31

    EggenbergerK.SanyalP.HundtS.WadhwaniP.UlrichA. S.NickP. (2017). Challenge integrity: the cell-penetrating peptide BP100 interferes with the auxin-actin oscillator.Plant Cell Physiol.587185. 10.1093/pcp/pcw161

  • 32

    Falcone FerreyraM. L.RiusS. P.CasatiP. (2012). Flavonoids: biosynthesis, biological functions, and biotechnological applications.Front. Plant Sci.3:222. 10.3389/fpls.2012.00222

  • 33

    FeildT. S.LeeD. W.HolbrookN. M. (2001). Why leaves turn red in autumn. The role of anthocyanins in senescing leaves of red-osier dogwood.Plant Physiol.127566574. 10.1104/pp.010063

  • 34

    GälweilerL.GuanC.MüllerA.WismanE.MendgenK.YephremovA.et al (1998). Regulation of polar auxin transport by AtPIN1 in Arabidopsis vascular tissue.Science28222262230. 10.1126/science.282.5397.2226

  • 35

    GeislerM.MurphyA. S. (2006). The ABC of auxin transport: the role of p-glycoproteins in plant development.FEBS Lett.58010941102. 10.1016/j.febslet.2005.11.054

  • 36

    GeldnerN.FrimlJ.StierhofY.-D.JürgensG.PalmeK. (2001). Auxin transport inhibitors block PIN1 cycling and vesicle trafficking.Nature413425428. 10.1038/35096571

  • 37

    GouJ.-Y.FelippesF. F.LiuC.-J.WeigelD.WangJ.-W. (2011). Negative regulation of anthocyanin biosynthesis in Arabidopsis by a miR156-targeted SPL transcription factor.Plant Cell2315121522. 10.1105/tpc.111.084525

  • 38

    HaoS.AugustA. (2005). Actin depolymerization transduces the strength of B-cell receptor stimulation.Mol. Biol. Cell1622752284. 10.1091/mbc.e04-10-0881

  • 39

    HarborneJ. B.BaxterH.(eds.) (1999). The Handbook of Natural Flavonoids,Vol. 1 and 2. Chichester: John Wiley and Sons, 889.

  • 40

    HawkinsC.LiuZ. (2014). A model for an early role of auxin in Arabidopsis gynoecium morphogenesis.Front. Plant Sci.5:327. 10.3389/fpls.2014.00327

  • 41

    HellemansJ.MortierG.De PaepeA.SpelemanF.VandesompeleJ. (2007). qBase relative quantification framework and software for management and automated analysis of real-time quantitative PCR data.Genome Biol.8:R19. 10.1186/gb-2007-8-2-r19

  • 42

    Henty-RidillaJ. L.LiJ.DayB.StaigerC. J. (2014). ACTIN DEPOLYMERIZING FACTOR4 regulates actin dynamics during innate immune signaling in Arabidopsis.Plant Cell26340352. 10.1105/tpc.113.122499

  • 43

    HimelblauE.AmasinoR. M. (2001). Nutrients mobilized from leaves of Arabidopsis thaliana during leaf senescence.J. Plant Physiol.15813171323. 10.1078/0176-1617-00608

  • 44

    HodgeE. T.SacherJ. A. (1975). Effects of kinetin, auxin and abscisic acid on ribonuclease and acid phosphatase during senescence of leaf tissue.Biochem. Physiol. Pflanz.168433441. 10.1016/S0015-3796(17)30148-8

  • 45

    HoltonT.CornishE. (1995). Genetics and biochemistry of anthocyanin biosynthesis.Plant Cell710711083. 10.1105/tpc.7.7.1071

  • 46

    HouK.WuW.GanS.-S. (2013). SAUR36, a small auxin up RNA gene, is involved in the promotion of leaf senescence in Arabidopsis.Plant Physiol.16110021009. 10.1104/pp.112.212787

  • 47

    HuX.TanakaA.TanakaR. (2013). Simple extraction methods that prevent the artifactual conversion of chlorophyll to chlorophyllide during pigment isolation from leaf samples.Plant Methods9:19. 10.1186/1746-4811-9-19

  • 48

    HuangX.MaischJ.NickP. (2017). Sensory role of actin in auxin-dependent responses of tobacco BY-2.J. Plant Physiol.218615. 10.1016/j.jplph.2017.07.011

  • 49

    IvakovA.PerssonS. (2013). Plant cell shape: modulators and measurements.Front. Plant Sci.4:439. 10.3389/fpls.2013.00439

  • 50

    JibranR.A HunterD.P DijkwelP. (2013). Hormonal regulation of leaf senescence through integration of developmental and stress signals.Plant Mol. Biol.82547561. 10.1007/s11103-013-0043-2

  • 51

    KandasamyM. K.McKinneyE. C.DealR. B.MeagherR. B. (2005). Arabidopsis ARP7 is an essential actin-related protein required for normal embryogenesis, plant architecture, and floral organ abscission.Plant Physiol.13820192032. 10.1104/pp.105.065326

  • 52

    KasarasA.MelzerM.KunzeR. (2012). Arabidopsis senescence-associated protein DMP1 is involved in membrane remodeling of the ER and tonoplast.BMC Plant Biol.12:54. 10.1186/1471-2229-12-54

  • 53

    KeechO. (2011). The conserved mobility of mitochondria during leaf senescence reflects differential regulation of the cytoskeletal components in Arabidopsis thaliana.Plant Signal. Behav.6147150. 10.4161/psb.6.1.14307

  • 54

    KeechO.PesquetE.GutierrezL.AhadA.BelliniC.SmithS. M.et al (2010). Leaf senescence is accompanied by an early disruption of the microtubule network in Arabidopsis.Plant Physiol.15417101720. 10.1104/pp.110.163402

  • 55

    KhanM.RozhonW.PoppenbergerB. (2014). The role of hormones in the aging of plants – A mini-review.Gerontology604955. 10.1159/000354334

  • 56

    KimJ.WooH. R.NamH. G. (2016). Toward systems understanding of leaf senescence: an integrated multi-omics perspective on leaf senescence research.Mol. Plant9813825. 10.1016/j.molp.2016.04.017

  • 57

    KimJ. I.MurphyA. S.BaekD.LeeS.-W.YunD.-J.BressanR. A.et al (2011). YUCCA6 over-expression demonstrates auxin function in delaying leaf senescence in Arabidopsis thaliana.J. Exp. Bot.6239813992. 10.1093/jxb/err094

  • 58

    Kleine-VehnJ.DhonuksheP.SwarupR.BennettM.FrimlJ. (2006). Subcellular trafficking of the Arabidopsis auxin influx carrier AUX1 uses a novel pathway distinct from PIN1.Plant Cell1831713181. 10.1105/tpc.106.042770

  • 59

    Kleine-VehnJ.FrimlJ. (2008). Polar targeting and endocytic recycling in auxin-dependent plant development.Annu. Rev. Cell Dev. Biol.24447473. 10.1146/annurev.cellbio.24.110707.175254

  • 60

    Kleine-VehnJ.WabnikK.MartinièreA.ŁangowskiŁ.WilligK.NaramotoS.et al (2011). Recycling, clustering, and endocytosis jointly maintain PIN auxin carrier polarity at the plasma membrane.Mol. Syst. Biol.7:540. 10.1038/msb.2011.72

  • 61

    KrecekP.SkupaP.LibusJ.NaramotoS.TejosR.FrimlJ.et al (2009). The PIN-FORMED (PIN) protein family of auxin transporters.Genome Biol.10:249. 10.1186/gb-2009-10-12-249

  • 62

    KuhnB. M.ErrafiS.BucherR.DobrevP.GeislerM.BiglerL.et al (2016). 7-Rhamnosylated flavonols modulate homeostasis of the plant hormone auxin and affect plant development.J. Biol. Chem.29153855395. 10.1074/jbc.M115.701565

  • 63

    KurthE. G.PeremyslovV. V.TurnerH. L.MakarovaK. S.IranzoJ.MekhedovS. L.et al (2017). Myosin-driven transport network in plants.Proc. Natl. Acad. Sci. U.S.A.114E1385E1394. 10.1073/pnas.1620577114

  • 64

    LampugnaniE. R.KilincA.SmythD. R. (2013). Auxin controls petal initiation in Arabidopsis.Development140185194. 10.1242/dev.084582

  • 65

    LarssonE.RobertsC. J.ClaesA. R.FranksR. G.SundbergE. (2014). Polar auxin transport is essential for medial versus lateral tissue specification and vascular-deiated valve outgrowth in Arabidopsis gynoecia.Plant Physiol.16619982012. 10.1104/pp.114.245951

  • 66

    LewisD. R.MudayG. K. (2009). Measurement of auxin transport in Arabidopsis thaliana.Nat. Protoc.4437451. 10.1038/nprot.2009.1

  • 67

    LewisD. R.RamirezM. V.MillerN. D.VallabhaneniP.RayW. K.HelmR. F.et al (2011). Auxin and ethylene induce flavonol accumulation through distinct transcriptional networks.Plant Physiol.156144164. 10.1104/pp.111.172502

  • 68

    LiG.LiangW.ZhangX.RenH.HuJ.BennettM. J.et al (2014). Rice actin-binding protein RMD is a key link in the auxin–actin regulatory loop that controls cell growth.Proc. Natl. Acad. Sci. U.S.A.1111037710382. 10.1073/pnas.1401680111

  • 69

    LimP. O.LeeI. C.KimJ.KimH. J.RyuJ. S.WooH. R.et al (2010). Auxin response factor 2 (ARF2) plays a major role in regulating auxin-mediated leaf longevity.J. Exp. Bot.6114191430. 10.1093/jxb/erq010

  • 70

    LinJ.-F.WuS.-H. (2004). Molecular events in senescing Arabidopsis leaves.Plant J.39612628. 10.1111/j.1365-313X.2004.02160.x

  • 71

    LohmanK. N.GanS.JohnM. C.AmasinoR. M. (1994). Molecular analysis of natural leaf senescence in Arabidopsis thaliana.Physiol. Plant.92322328. 10.1111/j.1399-3054.1994.tb05343.x

  • 72

    LuschnigC.GaxiolaR. A.GrisafiP.FinkG. R. (1998). EIR1, a root-specific protein involved in auxin transport, is required for gravitropism in Arabidopsis thaliana.Genes Dev.1221752187. 10.1101/gad.12.14.2175

  • 73

    MadisonS. L.BuchananM. L.GlassJ. D.McClainT. F.ParkE.NebenführA. (2015). Class XI myosins move specific organelles in pollen tubes and are required for normal fertility and pollen tube growth in Arabidopsis.Plant Physiol.16919461960. 10.1104/pp.15.01161

  • 74

    MaillardA.DiquélouS.BillardV.LaînéP.GarnicaM.PrudentM.et al (2015). Leaf mineral nutrient remobilization during leaf senescence and modulation by nutrient deficiency.Front. Plant Sci.6:317. 10.3389/fpls.2015.00317

  • 75

    MarchantA.KargulJ.MayS. T.MullerP.DelbarreA.Perrot-RechenmannC.et al (1999). AUX1 regulates root gravitropism in Arabidopsis by facilitating auxin uptake within root apical tissues.EMBO J.1820662073. 10.1093/emboj/18.8.2066

  • 76

    MisyuraM.ColasantiJ.RothsteinS. J. (2013). Physiological and genetic analysis of Arabidopsis thaliana anthocyanin biosynthesis mutants under chronic adverse environmental conditions.J. Exp. Bot.64229240. 10.1093/jxb/ers328

  • 77

    MurashigeT.SkoogF. (1962). A revised medium for rapid growth and bio assays with tobacco tissue cultures.Physiol. Plant.15473497. 10.1111/j.1399-3054.1962.tb08052.x

  • 78

    MurphyA.PeerW. A.TaizL. (2000). Regulation of auxin transport by aminopeptidases and endogenous flavonoids.Planta211315324. 10.1007/s004250000300

  • 79

    MushtaqM. A.PanQ.ChenD.ZhangQ.GeX.LiZ. (2016). Comparative leaves transcriptome analysis emphasizing on accumulation of anthocyanins in Brassica: molecular regulation and potential interaction with photosynthesis.Front. Plant Sci.7:311. 10.3389/fpls.2016.00311

  • 80

    NakataM.Ohme-TakagiM. (2014). Quantification of anthocyanin content.BIO-Protoc.4:e1098. 10.21769/BioProtoc.1098

  • 81

    NemhauserJ. L.FeldmanL. J.ZambryskiP. C. (2000). Auxin and ETTIN in Arabidopsis gynoecium morphogenesis.Development12738773888.

  • 82

    NickP.HanM.-J.AnG. (2009). Auxin stimulates its own transport by shaping actin filaments.Plant Physiol.151155167. 10.1104/pp.109.140111

  • 83

    NoodénL. D.NoodénS. M. (1985). Effects of morphactin and other auxin transport inhibitors on soybean senescence and pod development.Plant Physiol.78263266. 10.1104/pp.78.2.263

  • 84

    OjanguE.-L.JärveK.PavesH.TruveE. (2007). Arabidopsis thaliana myosin XIK is involved in root hair as well as trichome morphogenesis on stems and leaves.Protoplasma230193202. 10.1007/s00709-006-0233-8

  • 85

    OjanguE.-L.TannerK.PataP.JärveK.HolwegC. L.TruveE.et al (2012). Myosins XI-K, XI-1, and XI-2 are required for development of pavement cells, trichomes, and stigmatic papillae in Arabidopsis.BMC Plant Biol.12:81. 10.1186/1471-2229-12-81

  • 86

    OkadaK.UedaJ.KomakiM.BellC.ShimuraY. (1991). Requirement of the auxin polar transport system in early stages of Arabidopsis floral bud formation.Plant Cell3677684. 10.1105/tpc.3.7.677

  • 87

    Oñate-SánchezL.Vicente-CarbajosaJ. (2008). DNA-free RNA isolation protocols for Arabidopsis thaliana, including seeds and siliques.BMC Res. Notes1:93. 10.1186/1756-0500-1-93

  • 88

    PaponovI. A.PaponovM.TealeW.MengesM.ChakraborteeS.MurrayJ. A. H.et al (2008). Comprehensive transcriptome analysis of auxin responses in Arabidopsis.Mol. Plant1321337. 10.1093/mp/ssm021

  • 89

    ParkE.NebenführA. (2013). Myosin XIK of Arabidopsis thaliana accumulates at the root hair tip and is required for fast root hair growth.PLoS One8:e76745. 10.1371/journal.pone.0076745

  • 90

    PeerW. A.BandyopadhyayA.BlakesleeJ. J.MakamS. N.ChenR. J.MassonP. H.et al (2004). Variation in expression and protein localization of the PIN family of auxin efflux facilitator proteins in flavonoid mutants with altered auxin transport in Arabidopsis thaliana.Plant Cell1618981911. 10.1105/tpc.021501

  • 91

    PeerW. A.MurphyA. S. (2007). Flavonoids and auxin transport: modulators or regulators?Trends Plant Sci.12556563. 10.1016/j.tplants.2007.10.003

  • 92

    PengM.HudsonD.SchofieldA.TsaoR.YangR.GuH.et al (2008). Adaptation of Arabidopsis to nitrogen limitation involves induction of anthocyanin synthesis which is controlled by the NLA gene.J. Exp. Bot.5929332944. 10.1093/jxb/ern148

  • 93

    PeremyslovV. V.ColeR. A.FowlerJ. E.DoljaV. V. (2015). Myosin-powered membrane compartment drives cytoplasmic streaming, cell expansion and plant development.PLoS One10:e0139331. 10.1371/journal.pone.0139331

  • 94

    PeremyslovV. V.KlockoA. L.FowlerJ. E.DoljaV. V. (2012). Arabidopsis myosin XI-K localizes to the motile endomembrane vesicles associated with F-actin.Front. Plant Sci.3:184. 10.3389/fpls.2012.00184

  • 95

    PeremyslovV. V.MocklerT. C.FilichkinS. A.FoxS. E.JaiswalP.MakarovaK. S.et al (2011). Expression, splicing, and evolution of the myosin gene family in plants.Plant Physiol.15511911204. 10.1104/pp.110.170720

  • 96

    PeremyslovV. V.MorgunE. A.KurthE. G.MakarovaK. S.KooninE. V.DoljaV. V. (2013). Identification of myosin XI receptors in Arabidopsis defines a distinct class of transport vesicles.Plant Cell2530223038. 10.1105/tpc.113.113704

  • 97

    PeremyslovV. V.ProkhnevskyA. I.AvisarD.DoljaV. V. (2008). Two class XI myosins function in organelle trafficking and root hair development in Arabidopsis.Plant Physiol.14611091116. 10.1104/pp.107.113654

  • 98

    PeremyslovV. V.ProkhnevskyA. I.DoljaV. V. (2010). Class XI myosins are required for development, cell expansion, and F-actin organization in Arabidopsis.Plant Cell2218831897. 10.1105/tpc.110.076315

  • 99

    PouteauS.AlbertiniC. (2009). The significance of bolting and floral transitions as indicators of reproductive phase change in Arabidopsis.J. Exp. Bot.6033673377. 10.1093/jxb/erp173

  • 100

    ProkhnevskyA. I.PeremyslovV. V.DoljaV. V. (2008). Overlapping functions of the four class XI myosins in Arabidopsis growth, root hair elongation, and organelle motility.Proc. Natl. Acad. Sci. U.S.A.1051974419749. 10.1073/pnas.0810730105

  • 101

    RenH.GrayW. M. (2015). SAUR proteins as effectors of hormonal and environmental signals in plant growth.Mol. Plant811531164. 10.1016/j.molp.2015.05.003

  • 102

    SanteliaD.HenrichsS.VincenzettiV.SauerM.BiglerL.KleinM.et al (2008). Flavonoids redirect PIN-mediated polar auxin fluxes during root gravitropic responses.J. Biol. Chem.2833121831226. 10.1074/jbc.M710122200

  • 103

    ScarpellaE.MarcosD.FrimlJ.BerlethT. (2006). Control of leaf vascular patterning by polar auxin transport.Genes Dev.2010151027. 10.1101/gad.1402406

  • 104

    ScheuringD.LöfkeC.KrügerF.KittelmannM.EisaA.HughesL.et al (2016). Actin-dependent vacuolar occupancy of the cell determines auxin-induced growth repression.Proc. Natl. Acad. Sci. U.S.A.113452457. 10.1073/pnas.1517445113

  • 105

    SchippersJ.JingH.-C.HilleJ.DijkwelP. (2007). Developmental and hormonal control of leaf senescence.Annu. Plant Rev.26145170. 10.1002/9780470988855.ch7

  • 106

    SchmalenbachI.ZhangL.ReymondM.Jiménez-GómezJ. M. (2014). The relationship between flowering time and growth responses to drought in the Arabidopsis Landsberg erecta x Antwerp-1 population.Front. Plant Sci.5:609. 10.3389/fpls.2014.00609

  • 107

    SekhonR. S.ChildsK. L.SantoroN.FosterC. E.BuellC. R.de LeonN.et al (2012). Transcriptional and metabolic analysis of senescence induced by preventing pollination in maize.Plant Physiol.15917301744. 10.1104/pp.112.199224

  • 108

    SessionsA.NemhauserJ. L.McCollA.RoeJ. L.FeldmannK. A.ZambryskiP. C. (1997). ETTIN patterns the Arabidopsis floral meristem and reproductive organs.Development12444814491.

  • 109

    ShibasakiK.UemuraM.TsurumiS.RahmanA. (2009). Auxin response in Arabidopsis under cold stress: underlying molecular mechanisms.Plant Cell2138233838. 10.1105/tpc.109.069906

  • 110

    SmertenkoA.Franklin-TongV. E. (2011). Organisation and regulation of the cytoskeleton in plant programmed cell death.Cell Death Differ.1812631270. 10.1038/cdd.2011.39

  • 111

    SmertenkoA. P.BozhkovP. V.FilonovaL. H.von ArnoldS.HusseyP. J. (2003). Re-organisation of the cytoskeleton during developmental programmed cell death in Picea abies embryos.Plant J.33813824. 10.1046/j.1365-313X.2003.01670.x

  • 112

    SolfanelliC.PoggiA.LoretiE.AlpiA.PerataP. (2006). Sucrose-specific induction of theanthocyanin biosynthetic pathway in Arabidopsis.Plant Physiol.140637646. 10.1104/pp.105.072579

  • 113

    SundbergE.ØstergaardL. (2009). Distinct and dynamic auxin activities during reproductive development.Cold Spring Harb. Perspect. Biol.1:a001628. 10.1101/cshperspect.a001628

  • 114

    SwartzbergD.DaiN.GanS.AmasinoR.GranotD. (2006). Effects of cytokinin production under two SAG promoters on senescence and development of tomato plants.Plant Biol.8579586. 10.1055/s-2006-924240

  • 115

    SwarupR.FrimlJ.MarchantA.LjungK.SandbergG.PalmeK.et al (2001). Localization of the auxin permease AUX1 suggests two functionally distinct hormone transport pathways operate in the Arabidopsis root apex.Genes Dev.1526482653. 10.1101/gad.210501

  • 116

    SwarupR.PéretB. (2012). AUX/LAX family of auxin influx carriers — An overview.Front. Plant Sci.3:225. 10.3389/fpls.2012.00225

  • 117

    TaltsK.IlauB.OjanguE.-L.TannerK.PeremyslovV. V.DoljaV. V.et al (2016). Arabidopsis myosins XI1, XI2, and XIK are crucial for gravity-induced bending of inflorescence stems.Front. Plant Sci.7:1932. 10.3389/fpls.2016.01932

  • 118

    TangL. P.LiX. M.DongY. X.ZhangX. S.SuY. H. (2017). Microfilament depolymerization is a pre-requisite for stem cell formation during in vitro shoot regeneration in Arabidopsis.Front. Plant Sci.8:158. 10.3389/fpls.2017.00158

  • 119

    ThomasS. G.HuangS.LiS.StaigerC. J.Franklin-TongV. E. (2006). Actin depolymerization is sufficient to induce programmed cell death in self-incompatible pollen.J. Cell Biol.174221229. 10.1083/jcb.200604011

  • 120

    TitapiwatanakunB.BlakesleeJ. J.BandyopadhyayA.YangH.MravecJ.SauerM.et al (2009). ABCB19/PGP19 stabilises PIN1 in membrane microdomains in Arabidopsis.Plant J.572744. 10.1111/j.1365-313X.2008.03668.x

  • 121

    UedaH.YokotaE.KutsunaN.ShimadaT.TamuraK.ShimmenT.et al (2010). Myosin-dependent endoplasmic reticulum motility and F-actin organization in plant cells.Proc. Natl. Acad. Sci. U.S.A.10768946899. 10.1073/pnas.0911482107

  • 122

    UlmasovT.MurfettJ.HagenG.GuilfoyleT. J. (1997). Aux/IAA proteins repress expression of reporter genes containing natural and highly active synthetic auxin response elements.Plant Cell919631971. 10.1105/tpc.9.11.1963

  • 123

    UngererM. C.HalldorsdottirS. S.ModliszewskiJ. L.MackayT. F. C.PuruggananM. D. (2002). Quantitative trait loci for inflorescence development in Arabidopsis thaliana.Genetics16011331151.

  • 124

    VannesteS.FrimlJ. (2009). Auxin: a trigger for change in plant development.Cell13610051016. 10.1016/j.cell.2009.03.001

  • 125

    VernouxT.BesnardF.TraasJ. (2010). Auxin at the shoot apical meristem.Cold Spring Harb. Perspect. Biol.2:a001487. 10.1101/cshperspect.a001487

  • 126

    WangC.YanX.ChenQ.JiangN.FuW.MaB.et al (2013). Clathrin light chains regulate clathrin-mediated trafficking, auxin signaling, and development in Arabidopsis.Plant Cell25499516. 10.1105/tpc.112.108373

  • 127

    WangQ.KohlenW.RossmannS.VernouxT.TheresK. (2014). Auxin depletion from the leaf axil conditions competence for axillary meristem formation in Arabidopsis and tomato.Plant Cell2620682079. 10.1105/tpc.114.123059

  • 128

    WangY.-S.YooC.-M.BlancaflorE. B. (2008). Improved imaging of actin filaments in transgenic Arabidopsis plants expressing a green fluorescent protein fusion to the C- and N-termini of the fimbrin actin-binding domain 2.New Phytol.177525536. 10.1111/j.1469-8137.2007.02261.x

  • 129

    WatanabeM.BalazadehS.TohgeT.ErbanA.GiavaliscoP.KopkaJ.et al (2013). Comprehensive dissection of spatiotemporal metabolic shifts in primary, secondary, and lipid metabolism during developmental senescence in Arabidopsis.Plant Physiol.16212901310. 10.1104/pp.113.217380

  • 130

    WeaverL. M.GanS. S.QuirinoB.AmasinoR. M. (1998). A comparison of the expression patterns of several senescence-associated genes in response to stress and hormone treatment.Plant Mol. Biol.37455469. 10.1023/A:1005934428906

  • 131

    WeaverL. M.HimelblauE.AmasinoR. M. (1997). Leaf senescence: gene expression and regulation.Genetic Eng.19215234. 10.1007/978-1-4615-5925-2_12

  • 132

    WeigelD.GlazebrookJ. (2002). Arabidopsis: A Laboratory Manual. New York, NY: CSHL Press, 355.

  • 133

    WooH. R.KimH. J.NamH. G.LimP. O. (2013). Plant leaf senescence and death – Regulation by multiple layers of control and implications for aging in general.J. Cell Sci.12648234833. 10.1242/jcs.109116

  • 134

    WuS.XieY.ZhangJ.RenY.ZhangX.WangJ.et al (2015). VLN2 regulates plant architecture by affecting microfilament dynamics and polar auxin transport in rice.Plant Cell2728292845. 10.1105/tpc.15.00581

  • 135

    XuG.LiS.XieK.ZhangQ.WangY.TangY.et al (2012). Plant ERD2-like proteins function as endoplasmic reticulum luminal protein receptors and participate in programmed cell death during innate immunity.Plant J.725769. 10.1111/j.1365-313X.2012.05053.x

  • 136

    ZažímalováE.MurphyA. S.YangH.HoyerováK.HošekP. (2010). Auxin transporters — Why so many?Cold Spring Harb. Perspect. Biol.2:a001552. 10.1101/cshperspect.a001552

  • 137

    ZentgrafU.JobstJ.KolbD.RentschD. (2004). Senescence-related gene expression profiles of rosette leaves of Arabidopsis thaliana: leaf age versus plant age.Plant Biol.6178183. 10.1055/s-2004-815735

  • 138

    ZhuJ.BaillyA.ZwiewkaM.SoveroV.Di DonatoM.GeP.et al (2016). TWISTED DWARF1 mediates the action of auxin transport inhibitors on actin cytoskeleton dynamics.Plant Cell28930948. 10.1105/tpc.15.00726

  • 139

    ZhuJ.GeislerM. (2015). Keeping it all together: auxin-actin crosstalk in plant development.J. Exp. Bot.6649834998. 10.1093/jxb/erv308

  • 140

    Zúñiga-MayoV. M.Reyes-OlaldeJ. I.Marsch-MartinezN.de FolterS. (2014). Cytokinin treatments affect the apical-basal patterning of the Arabidopsis gynoecium and resemble the effects of polar auxin transport inhibition.Front. Plant Sci.5:191. 10.3389/fpls.2014.00191

Summary

Keywords

Arabidopsis, myosin XI, auxin response, PIN1, flower development, senescence, SAG13, anthocyanins

Citation

Ojangu E-L, Ilau B, Tanner K, Talts K, Ihoma E, Dolja VV, Paves H and Truve E (2018) Class XI Myosins Contribute to Auxin Response and Senescence-Induced Cell Death in Arabidopsis. Front. Plant Sci. 9:1570. doi: 10.3389/fpls.2018.01570

Received

29 May 2018

Accepted

08 October 2018

Published

27 November 2018

Volume

9 - 2018

Edited by

Jian Xu, National University of Singapore, Singapore

Reviewed by

Andreas Nebenführ, The University of Tennessee, Knoxville, United States; Jie Le, Institute of Botany (CAS), China

Updates

Copyright

*Correspondence: Eve-Ly Ojangu,

This article was submitted to Plant Evolution and Development, a section of the journal Frontiers in Plant Science

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics