- Research article
- Open Access
The tomato WV gene encoding a thioredoxin protein is essential for chloroplast development at low temperature and high light intensity
BMC Plant Biologyvolume 19, Article number: 265 (2019)
Chloroplast biogenesis, a complex process in higher plants, is the key to photoautotrophic growth in plants. White virescent (wv) mutants have been used to unfold the molecular mechanisms underlying the regulation of chloroplast development and chloroplast gene expression in plants. However, most of genes controlling white virescent phenotype still remain unknown.
In this study, we identified a temperature- and light intensity-sensitive mutant, named as wv. The content of chlorophyll was dramatically decreased in the immature leaves of wv mutant under the conditions of low temperature and high-light intensity. TEM observation showed that the chloroplasts in the young leaves of wv mutant lacked an organized thylakoid membrane, whereas crescent-shaped chloroplasts with well-developed stromal and stacked grana thylakoids in the mature leaves were developed. Immunoblot analyses suggested that proteins of photosynthetic complexes were decreased substantially in wv mutants. Based on map-based cloning and transgenic analysis, we determined that the wv phenotype was caused by single base mutation in the first intron of WV gene, which encoded a thioredoxin protein with 365 amino acids. qRT-PCR analysis revealed that the expression of WV gene was significantly down-regulated in wv mutant. In addition, knockdown of WV gene through RNAi also resulted in white virescent young leaves, suggesting that the mutation possibly blocks the differentiation of chloroplasts through inhibiting the expression of WV gene. Furthermore, the expression of WV peaked in apical buds and gradually decreased along with the developmental stage, which was consistent with the wv mutant phenotype. Expression analysis of chloroplast-encoded genes by qRT-PCR showed that the wv mutation affected the expression pattern of chloroplast-encoded PEP dependent genes.
Our results suggested that wv mutant was sensitive to low temperature and light intensity. WV gene was essential for chloroplast differentiation. A single base mutation in the first intron resulted in down-regulation of WV gene expression, which inhibited the expression of chloroplast-encoded genes, thereby blocking chloroplast formation and chlorophyll synthesis.
Chloroplasts are important cellular organelles of photoautotrophic eukaryotes, which differentiate from plastids in plant cells . Chloroplasts are coordinately controlled by the proteins encoded in both plastid and nuclear genomes . These proteins play important roles in protein translation, foldingand transport during chloroplast development . Chloroplast development is tightly coupled with the leaf cell division and elongation. The later developmental process requires illumination, including chlorophyll accumulation and construction of the photosynthetic machinery . Most of photosynthesis and some metabolic processes occurred in chloroplasts, such as synthesis of amino acids, fatty acids, carotenoids, vitamins, and a range of specialized metabolites .
Chlorophyll is critical for photosynthesis, which can capture light energy and drive photosynthetic reactions . Chlorophyll biosynthesis is cooperatively regulated by 16 enzymes in Arabidopsis . Glutamyl tRNA reductases (GluTR) encoded by HEMA1, HEMA2, and HEMA3, are important limiting factors of chlorophyll synthesis . 5-aminolevulinic acid (ALA) is converted into uroporphyrinogen III and then oxidized into Proto by protoporphyrinogen oxidase. Three enzymes, including MgProtoMe cyclase, NADPH-Pchlide oxidoreductase, and Chl synthase, catalyze Proto into chlorophyll . Block of each step can cause changes in chlorophyll content and leaf color. Different mutations caused diverse leaf colors, like albino, yellow, chlorina, light green, dark green, etc. In addition, leaf color patterns are tremendously varied, such as mono-colored, spotted, striped, blotched, zebra, and variegated . Recently, many other leaf color mutants have been reported in rice [11, 12], wheat , maize , and foxtail millet .
Nowadays, many genes involved in chlorophyll biosynthesis or chloroplast formation have been characterized in many plants [7, 16,17,18]. YGL1, encoding a chlorophyll synthase, functions in the esterification of chlorophyllide and phytol in rice . NYC encodes a chlorophyll b reductase, the mutation of which impaired chlorophyll degradation during senescence . Chloroplast formation is essential for chlorophyll biosynthesis. Therefore, mutations in many genes required for chloroplast formation also caused leaf color changes. AtSIG6, a nuclear-encoded sigma factor in Arabidopsis, is critical for light-dependent chloroplast development. AtSIG6 null mutant exhibited pale green phenotype in cotyledons . EMB1303, a chloroplast-localized protein, controls chloroplast development in Arabidopsis. The emb1303–1 mutant displayed albino rosette leaves . In addition, overexpression of a RanBP2-type zinc finger transcription factor, SlRBZ, possibly impaired the biosynthesis of chlorophyll and resulted in chlorosis phenotype through blocking chloroplast development in tomato .
White or yellow virescent leaf mutants widely exist in nature, which have accelerated the characterization of molecular mechanisms underlying chlorophyll synthesis and chloroplast development. Chlorina-1 mutant exhibited yellowish-green leaf phenotype under normal growth conditions at seedling stage, and then turned into normal green . virescent-2 mutant developed chlorotic leaves at restrictive temperature (20 °C), whereas developed green leaves at permissive temperature (30 °C) . Two rice temperature-conditional mutants, virescent3 (v3) and stripe1 (st1), produced bleached leaves at a constant temperature of 20 °C or 30 °C, and green leaves under diurnal 30 °C/20 °C conditions . In addition, one mutation in NUS1 gene caused chlorotic rice leaves under low-temperature conditions . Tu, a chloroplast protein synthesis elongation factor, was involved in chloroplast development. Loss-of-function mutation in Tu resulted in green-revertible albino phenotype in rice .
In this study, we identified a chlorophyll deficient mutant, named as white virescent (wv), which showed white virescent apical buds and immature leaves. We developed insertion-deletion (InDel) markers for genotyping of an F2 segregating population. wv was delimited in an approximate 94 kb region, which contained 18 putative candidates. Transgenic analysis indicated that the wv phenotype was caused by a single base mutation in the first intron of WV gene, which encodes a thioredoxin protein of 365 amino acids. qRT-PCR analysis and GUS activity revealed that the expression level of WV peaked in apical buds, and gradually decreased along with the developmental stage. Knockdown of WV gene by RNAi also caused white virescent young leaves. We thus speculate that the single base mutation in WV gene possibly blocked chloroplast differentiation by inhibiting the transcription of WV gene.
wv mutant exhibits a thermo-sensitive phenotype
wv mutant LA1526 displayed white/yellow leaves at 16 °C compared with AC plants (Fig. 1a). The yellow leaves gradually became virescent with temperature rising and exhibited green leaves at 30 °C. In addition, wv leaves gradually turned green from top to base (Fig. 1b). Interestingly, LA1526 showed the same leaf color as AC plants at 16 °C and low-light intensity of 50 μmol m− 2 s− 1 (Fig. 1c). These results suggested that wv mutant LA1526 may be sensitive to temperature and light. We further measured chlorophyll content and found that chlorophyll a, chlorophyll b and total chlorophyll drastically decreased in the immature leaves of LA1526 compared to AC plants at 16 °C and high-light intensity (250 μmol m− 2 s− 1) (Fig. 1d). However, both LA1526 and AC plants showed similar chlorophyll content when growing at 30 °C and 250 μmol m− 2 s− 1 intense luminosity or at 16 °C and intense luminosity (50 μmol m− 2 s− 1), respectively (Fig. 1e, f). The chlorophyll content in mature leaves of LA1526 and AC showed no significant difference regardless of growth situation (Fig. 1d, e, f). These results indicated that LA1526 exhibited the white virescent leaves that could be enhanced by low temperature and high light intensity.
Chloroplast development was impaired due to the defects in photosynthesis-related proteins accumulation in wv mutant
Chloroplast development has an important impact on chlorophyll content, the defect of which impaired chlorophyll synthesis . We observed the chloroplast ultrastructure by TEM and found that chloroplasts in AC immature, mature leaves, and LA1526 mature green leaves were crescent-shaped with well-developed stromal and stacked grana thylakoids (Fig. 2a, b, d). Additionally, starch grains were also observed in these plastids. However, the chloroplasts in LA1526 immature yellow leaves showed much smaller sizes and lacked organized membrane structures, including thylakoid membrane system and lamellar layer system of thylakoid. Moreover, chloroplasts contained no starch grains (Fig. 2c). These results demonstrated that wv possibly plays an important role in chloroplast biogenesis at an early stage. Analysis of total proteins by Coomassie blue staining assay after SDS-PAGE showed that many proteins, particularly the proteins of approximately 50 kDa in size, were decreased in the wv mutants (Fig. 2e). In addition, the accumulation of photosynthetic proteins in wv mutant was also detected by immunoblot analysis with the corresponding antibodies. The results showed that PsaD, PsbA, and rbcL, were sharply reduced in wv mutant (Fig. 2f; Additional file 1: Figure S1), suggesting that the accumulation of photosynthetic-related proteins was obviously inhibited in wv immature yellow leaves.
wv was delimited to a 94 kb physical region on chromosome 2
Genetic analysis indicated that the white virescent phenotype in wv mutant LA1526 was controlled by a single recessive gene (Additional file 2: Table S4). The high-density molecular linkage map showed that wv was mapped to an interval flanked by TG14 and TG454 on the long arm of chromosome 2 . These two markers were located at approximate 41.48 and 45.32 Mb on chromosome 2, respectively. To further determine the position of wv gene, we generated an F2 mapping population from the cross between LA1526 and IL2–3. Five markers, including 2–3-9, wv-c12, wv-c13, wv-c24, and 2–3-15 (Additional file 2: Table S5), were developed to screen 186 F2 recessive individuals. The linkage analysis showed that wv was delimited to an interval between wv-c13 and wv-c24, which were at a distance of approximate 0.27 and 1.09 cM, respectively (Fig. 3a). To further narrow the interval spanning the target locus, we conducted fine mapping based on 1602 F2 recessive plants. Forty-eight individuals displayed recombination events between two InDel markers, wv-c13 and wv-c24 (Fig. 3b). We thus developed six InDel molecular markers within this interval, including wv-c39, wv-c47, wc-c65, wv-c53, wv-c28, and wv-c75 (Additional file 2: Table S5). We investigated the genotypes of these recombinants by using these newly developed markers. We finally mapped wv gene in the region between wv-c53 and wv-c75, both of which were closely linked to wv with six recombination event identified. One marker wv-c28 was co-segregated with wv gene (Fig. 3b). We confirmed that wv gene was mapped to an approximately 94 kb fragment based on the tomato reference genome sequence (Fig. 3c).
Solyc02g079730 was the candidate gene of wv
Based on the tomato genome annotation [ITAG Release 2.5 predicted coding sequence (CDS)], 18 putative ORFs were predicted in the target region (Fig. 3c). Analysis of the ~ 94 kb sequence by GENESCAN and FGENESH showed identical result. The best hits of these ORFs include receptor like protein kinases (ORF1- ORF13), an EF-hand type centrin, thioredoxin protein, U-box domain-containing protein, Flavoprotein wrbA and bHLH transcription factor (ORF14-ORF18), respectively (Additional file 2: Table S6). To further determine the possible candidate, we amplified and sequenced the genomic and cDNA sequences of these candidates from wv mutant LA1526, AC and M82. Sequence alignment showed that ORF15 (Solyc02g079730) in wv mutant LA1526 contains a single base mutation (T to G) (Fig. 3d; Additional file 1: Figure S2). We thus considered Solyc02g079730 as a good candidate for wv gene. Solyc02g079730 encoded a thioredoxin family protein, which is comprised of 365 amino acids. Interestingly, this nucleotide substitution did not result in amino acid substitution, implying that it possibly caused transcriptional change (Fig. 3d). We examined the expression of this gene in young leaves of AC, LA1526, F1 and F2 individuals from the cross IL2–3 × LA1526 (three genotypes, WV/WV, WV/wv and wv/wv) by using qRT-PCR. The expression of Solyc02g079730 was significantly repressed in LA1526 and F2 plants with wv/wv genotype (Fig. 3e). Overall, we inferred that the single base mutation possibly decreased the expression of Solyc02g079730 and led to white virescent leaves in LA1526.
Mutation in Solyc02g079730 responsible for the white virescent phenotype in wv mutant
In order to investigate the biological function of Solyc02g079730, we prepared an expression construct by placing the full-length Solyc02g079730 from Ailsa Craig into pHellesgate 8 under control of its native promoter and introduced this construct into wv mutant LA1526 by stable transformation. We obtained three independent transgenic lines. qRT-PCR analysis showed that the expression of Solyc02g079730 was evidently up-regulated in transgenic plants (Fig. 4b; Additional file 1: Figure S3). As the sensitivity of white virescent phenotype to low temperature, we analyzed the phenotypic change of transgenic plants when growing at 16 °C. We found that the white virescent phenotype was rescued by the fragment from Ailsa Craig (Fig. 4a). Furthermore, we also generated RNAi transgenic lines (Ri) by using AC as recipient and found that they showed white virescent phenotype at 16 °C or 30 °C when Solyc02g079730 expression was decreased (Fig. 4a; Fig. 4b; Additional file 1: Figure S4). These results suggested that Solyc02g079730 is the candidate for wv gene. Through genomic and cDNA sequence comparison of wv, we determined that there were four exons and three introns. In addition, we conducted nucleotide and protein BLAST search against the National Center for Biotechnology Information and found that the homologs of WV gene widely exist in many different species, except for Synechocystis sp. PCC 6803 and Clamydomonas reinhardtii. Multi-alignment of the amino acid sequences encoded by homologues from other species indicated that they contained a highly conserved domain (Fig. 5a). WV showed a high amino acid sequence identity with a thioredoxin protein in cucumber (60%), Arabidopsis (63%), rice (67%), Asparagus officinalis (67%) and maize (69%), respectively., The homolog of WV in Arabidopsis had the highest homology and the same gene structure as SVR4/MRL7, which was involved in the regulation of early chloroplasts biogenesis [30, 31]. The phylogenetic tree demonstrated that these genes can be divided into two subgroups, in accordance with the biological taxonomy, monocots and dicots (Fig. 5b). These results suggested that an ancestral WV gene diverged prior to the split between monocots and dicots, and WV might be the functional ortholog of SVR4.
Expression pattern of WV
To examine the expression pattern of WV gene, we conducted qRT-PCR analysis with total RNA extracted from several tissues of AC and LA1526, including roots, stems, apical buds, immature leaves and mature leaves. WV was expressed in all these tested tissues, being highest in apical buds of AC plants (Fig. 6a). However, the expression of WV gene was dramatically decreased in several tissues of LA1526 compared to AC, such as stems, apical buds, immature leaves, and mature leaves (Fig. 6a). It has been suggested that some putative cis-elements in its promoter play important roles in the tissue-specific gene expression. Therefore, we detected cis-element of the promoter sequence of WV by using the Plant CARE and the PLACE program. We identified five types of putative cis-elements, including light response elements (3-AF1 binding site, ACE, Box 4, G-Box, Gap-box), MYB binding site (MBS), heat stress response element (HSE), defense and stress response element (TC-rich repeats), gibberellin-responsive element (GARE-motif). We further conducted WV promoter-driven GUS transformation to analyze the spatial expression pattern of WV gene. In accordance with our previous results, GUS was highly expressed in apical buds (Fig. 6b). Furthermore, we identified that the expression of WV was significantly down-regulated in both wv mutant and AC plants at 16 °C (Fig. 6c). Taken together, the expression pattern of WV is consistent with the target phenotype.
WV affected the expression of chloroplast-encoded genes
The expression of chloroplast-encoded genes was closely related to chloroplast developmental status, the expression of which was relied on collaboration between plastid-encoded polymerase (PEP) and nuclear-encoded polymerase (NEP) . Therefore, we investigated the transcript profiles of various chloroplast-encoded genes in 4-week-old AC and wv mutant (Fig. 7). The results showed that the PEP-dependent genes, like PsaB, PsaA, PsbA, and PsbB, were decreased significantly in wv mutant compared to AC plants. In contrast, the expression level of NEP-dependent genes, like AccD, Ycf2, RpoA, and RpoB, were increased in wv mutant. Among PEP&NEP dependent genes, the abundance of ClpP transcripts in wv mutant was about 5 times as much as that in AC plants, whereas the abundance of AtpE and AtpB transcripts in wv mutant was the same as that in AC plants. In addition, the expression of 16SrRNA was sharply reduced in wv mutant compared to AC plants. Furthermore, the expression level of nuclear-encoded genes destined for chloroplasts, such as CAO and Psbw, were much lower in wv mutant than that in AC plants. However, the transcripts of PetC and PsbO in wv mutant were almost identical to that in AC plants. These results indicated that the accumulation of many chloroplast-encoded genes was affected in wv mutants, which was probably the consequence of defective chloroplast biogenesis.
Plant leaves play essential roles in photosynthesis and then crop yield and quality. Chlorophyll biosynthesis and chloroplast formation are co-regulated by nuclear and plastid-encoded genes . Chlorophyll-deficient mutants, like albino, variegated and virescent mutants, have been used to characterize the molecular mechanisms underlying the chloroplast development and chlorophyll synthesis in plants. Recently, many genes responsible for chlorophyll-deficient phenotypes have been identified in several plants, like VAR2 in Arabidopsis and V2 in rice [34, 35]. In this study, we identified a white virescent mutant. It has been suggested that many leaf-color mutants are sensitive to low temperature. For example, mutation in Virescent 1 (V1) gene caused chlorotic leaves and sensitivity to low temperature . Maize virescent 16 mutant was unable to accumulate chlorophyll and exhibited chlorotic phenotype when grown at temperature lower than 25 °C . Another maize temperature-sensitive mutant, M-11, was not able to accumulate chlorophyll below 17 °C . Accordingly, we confirmed that wv mutant developed white leaves under low temperature condition, whereas exhibited green leaves under high temperature condition. Furthermore, we presented the evidence that the content of chlorophyll significantly decreased in the immature leaves of LA1526 at low temperature compared to that at high temperature. We thus speculated that the expression of WV gene responsible for wv phenotype may be sensitive to low temperature. In accordance with this reference, the expression of WV was significantly repressed in both wv mutant and AC plants at 16 °C compared to that at 30 °C. We further confirmed that the wv phenotype can be enhanced by high-light intensity, suggesting that WV may respond to light signal. It was further supported by the existence of temperature and light responsive cis-elements in the promoter of WV.
Since chlorophyll is synthesized in chloroplasts, impairment of chloroplast formation causes the reduction of chlorophyll content . MDA1 gene plays an important role in chloroplast formation, the mutation of which resulted in pigmentation reduction in Arabidopsis . SVR4 regulates early chloroplast biogenesis in Arabidopsis and svr4 mutant exhibit the virescent phenotype . In the present study, we observed the ultrastructure of chloroplasts through TEM and found that the chloroplasts in the young leaves of wv mutant lacked mature thylakoids and starch grains. Moreover, the expression of chloroplast-encoded PEP-dependent genes was obviously inhibited in immature leaves of wv mutant LA1526 compared to AC, including PsaB, PsaA, PsbA and PsbB. Consistently, a homolog of WV in Arabidopsis, SVR4, showed a similar regulatory role in the expression of chloroplast-encoded genes . Interestingly, WV and SVR4 had the same gene structure and conserved motif. Furthermore, the expression of chloroplast-encoded NEP dependent genes was significantly increased in immature leaves of wv mutant LA1526 compared to AC, such as AccD, Ycf2, RpoA and RpoB. And the expression of some nucleus-encoded photosynthetic genes was inhibited in immature leaves of LA1526, such as CAO and PsbW. Taken together, the data suggest that WV may be the functional ortholog of SVR4, which regulates chloroplast formation and chlorophyll accumulation through controlling the transcription of chloroplast-encoded PEP dependent genes. In addition, the up-regulated expression of NEP-dependent chloroplast genes and down-regulated expression of some nucleus-encoded photosynthetic genes were mainly influenced by the defects in the transcription of chloroplast-encoded PEP-dependent genes .
Based on the analysis of transgenic plants, we determined that the wv phenotype was caused by a single base mutation in WV gene, which encoded a thioredoxin family protein. Thioredoxins are small redox-active proteins, which are widely distributed among living organisms . Previous studies showed that some of thioredoxins were involved in the regulation of chlorophyll metabolism and other processes in chloroplasts . Therefore, knock-out mutations in some of these proteins usually impaired chloroplast formation and caused a severe reduction in chlorophyll accumulation . Interestingly, we found that the mutation in the first intron of WV gene did not lead to amino acid substitution. We thus speculated that the white virescent leaves of wv mutant was possibly caused by the reduced level of WV transcripts. Accordingly, the expression of WV gene was significantly decreased in wv mutant LA1526. Furthermore, knockdown of WV gene through RNAi also resulted in similar phenotype to wv mutant, implying that the single base mutation inhibited the expression of WV gene, and then led to white virescent leaves. It has been reported that several plant introns probably contain enhancer elements and can enhance gene expression . For example, the high-locule-number tomato fruit was possibly caused by an insertion in the first intron of fas gene, which repressed the expression of fas gene during early fruit development . Additionally, small mutations in the second intron of the Arabidopsis AG gene could inhibit its expression . Thus, it can be speculated that the first intron of WV gene possibly plays similar roles in the regulation of gene expression. Here, one pair of 11 bp reverse repeated sequences was found in the first intron of WV. This sequence-structure may be involved in intron-mediated enhancement. Coincidently, the nucleotide substitution occurred in the second repeat sequence, which possibly impaired the regulatory structure and inhibited the expression of WV gene. Moreover, the secondary structure of wv pre-mRNA may be altered by the single base mutation, which can reduce mRNA splicing efficiency and cause low abundance of WV gene transcripts. However, there is no more evidence to support this speculation. Therefore, it will be interesting to figure out the regulatory role of the first intron in the WV gene expression in the future.
In this study, we identified a temperature- and light intensity-sensitive mutant, named as wv mutant. Based on positional cloning and transgenic analysis, we speculated that the single base mutation in the first intron of WV gene caused significant decreased expression of WV gene and blocked the differentiation of chloroplasts and synthesis of chlorophyll through inhibiting the transcription of chloroplast-encoded PEP dependent genes, which resulted in white virescent phenotype. Our results will help us to understand the regulatory role of thioredoxin protein in chloroplast differentiation. Moreover, further researches are necessary to elucidate the intron-mediated enhancement of gene expression during chloroplast formation and chlorophyll synthesis.
Plant materials and growth conditions
One introgression line (IL) 2–3 and wv mutant LA1526 were provided by the Tomato Genetics Resource Center (TGRC). IL2–3 contains the middle fragment of chromosome 2 of S. pennellii LA0716 genome in M82 background . F2 mapping population segregated for the target phenotype was derived from a cross between IL2–3 and LA1526. All plants were grown under standard glasshouse conditions at 16 °C with 16 h day/8 h night cycle and 250 μmol m− 2 s− 1 intense luminosity. F2 individuals with white virescent leaves were used for mapping analysis. To detect whether wv mutant LA1526 was sensitive to temperature, plants were grown under the condition as follows: 16 h day/8 h night cycle, 24 °C, and 250 μmol m− 2 s− 1 intense luminosity. Two-week-old plants were divided into two groups, each of which was grown under 16 °C, and 30 °C conditions, respectively. Ailsa Craig (AC) plants as control were grown under the same conditions as described above. In addition, to detect the effect of light on wv mutant, we compared the second group with the third group, which was grown under 16 °C, 50 μmol m− 2 s− 1 intense luminosity.
Chlorophyll and carotenoid content determination
Immature yellow leaves (the second leaves) and mature green leaves (the fourth leaves) were collected from 4-week-old LA1526 and AC plants and treated with 80% acetone in the dark. After grinding and high speed centrifugation, supernatant liquor containing pigments were immediately spectrophotometrically assayed at specific absorption coefficients using a microplate reader as described by a previous study .
Transmission electron microscopy (TEM) observation
Samples (the second leaves and fourth leaves) collected from 4-week-old LA1526 and AC plants were cut into small pieces (~ 1 mm2) and inbubated in 3.5% (v/v) glutaraldehyde solution overnight. The fixed samples were then washed with 0.1 M phosphate buffer, post-fixed with 1% osmiophilic tetroxide, and dehydrated with a series of alcohol solutions, and then infiltrated and embedded with resin. The ultra-thin sections were prepared according to the methods described by Fan et al., 2016 . Finally, the samples were observed and photographed using TEM (Hitachi H-7650, Tokyo, Japan) and Gatan 832 digital imaging system.
Immunological detection of photosynthetic proteins
Total proteins were extracted from the top two leaves of 30-day-old AC plants and wv mutants LA1526 as described previous . Approximately 1 g of each sample leaves was grinded into powder in liquid nitrogen, then lysed with ice-cold plant lysis buffer (50 mM Tris-HCl (pH 7.5), 150 mM NaCl, 0.5% TritonX-100) with complete protease inhibitor cocktail tablets (Roche, Basel, Switzerland), respectively. After incubating on ice for 30 min, the samples were centrifuged at 13000 g for 20 min at 4 °C to obtain total proteins in the supernatants. The supernatants transferred into a new tube were mixed well with equal amounts of 5x loading buffer and then boiled for 10 min.
Approximately 30 μg of total protein per sample were separated on 12% (m/v) sodium dodecyl sulfate PAGE (SDS-PAGE). After electrophoresis, the contents of the total proteins were detected by coomassie blue staining according to the manufacturer’s instructions. For immunodetection, the gels were blotted onto polyvinylidene fluoride (PVDF) membranes (Trans-Blot® Turbo TM Mini PVDF Transfer, Bio-Rad). After blocking with 5% skim milk for 2 h at room temperature, the membranes were washed twice in tris buffered saline with Tween 20 (TBST) and then incubated overnight at 4 °C with anti- PsaD, PsbA and rbcL rabbit primary antibody diluted 1000 times, respectively. Subsequently, the membranes were washed five times in TBST and incubated with secondary antibodies antibody against rabbit IgG at 1:3000 dilutions for 1 h at room temperature. Finally, the proteins were detected using chemiluminescence (ECL) method (Millipore, Burlington, MA). For all samples, the protein level of actin was used as a loading control which was detected with mouse anti-actin (AT3G12110; actin-11) primary antibody (manufacturer details) at 1:300 dilution. All polyclonal antibodies used in this study were obtained from ABclonal Technology. The common protocols and the manufacturers’ manuals for electrophoresis, semi-dry blotting and western detection using the ECL western blotting substrate were followed in this study.
Positional cloning of wv gene
Plant genomic DNA was isolated from the young leaves of 3 to 4 weeks old seedlings by using a cetyltrimethyl ammonium bromide protocol . The concentration of each DNA sample was determined by a NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific) and adjusted to final 100 ng/μL. Primers for InDel markers were designed by using Primer 5.0 software based on the genomic sequences of M82 and LA0716 . 186 F2 recessive individuals from the cross IL2–3 × LA1526 were initially used for rough mapping of wv. Then, a larger population of 1602 F2 recessive individuals was used for fine mapping of wv gene. Genotype data of each individual were adopted for linkage analysis by using the MAPMAKER/EXP 3.0 program .
Molecular cloning and sequencing
We downloaded the genome sequence of Heinz 1706 in the target fragment based on the fine mapping results. The predicated genes in the target region were further analyzed by FGENESH (http://linux1.softberry.com/berry.phtml?topic=fgenesh&group=programs&subgroup=gfind) and GENESCAN (http://hollywood.mit.edu/GENSCAN.html). We amplified the candidate genes by using genomic DNA and cDNA of LA1526 and AC as templates. The primers are listed in Additional file 2: Table S1. PCR was performed in 20 μL mixture containing 100 ng of gDNA, 1 μL of 10 × Taq buffer, 0.2 μL of 10 mM dNTPs, 1 U of Taq DNA polymerase (Invitrogen, USA), and 2 μm of each primer. The amplification program was performed at 94 °C for 2 min, followed by 35 cycles at 94 °C for 30 s, 56 °C for 60 s, and 72 °C for 1 min, with a final extension at 72 °C for 5 min. PCR products were mixed with 2 μL of loading buffer and analyzed on a 2% agarose gel.
Construction of expression vectors and generation of transgenic lines
To repress the expression of the Solyc02g079730.2.1, RNAi fragment was amplified with the primers listed in Additional file 2: Table S2 and inserted into the vector pHellsgate2 by using clonase BP reaction (Invitrogen). Furthermore, the full length genomic DNA of WV gene together with its 2.5 kb promoter sequence was cloned into pMV2 vector as described by a previous study . The construct proWV::WV was transformed into wv mutant for functional complementation. In addition, we amplified an approximate 2.5 kb endogenous promoter fragment upstream of the start codon of WV gene. The fusion construct proWV::GUS was acquired by inserting the promoter fragment in front of the GUS coding region as described by Gao et al., 2017 . The proWV::WV fusion construct was transformed into LA1526, whereas other constructs were transformed into AC mediated by Agrobacterium tumefaciens strain GV3101. The empty vectors of pHellsgate2 and pMV2 were also transformed as negative controls.
Cis-element analysis of the WV gene promoter and histochemical GUS assay
An approximate 2.5 kb promoter fragment upstream the start codon of the WV gene was analyzed by using both NewPLACE (https://www.dna.affrc.go.jp/PLACE/?action=newplace) and plant CARE (http://bioinformatics.psb.ugent.be/webtools/plantcare/html/). Expression pattern of WV gene was reflected by GUS activity. In order to evaluate GUS activity, the 10 days old transgenic plants transformed with the proWV::GUS fusion construct were immersed into a GUS staining solution and incubated overnight at 37 °C after vacuum infiltrating for 5 min. Then, the staining buffer was removed and chlorophyll was removed by incubating all samples in 70% ethanol. The GUS activity were observed and photographed by a microscope.
Sequence and phylogenetic analysis
In order to identify homologs of WV gene in different species, pBLAST was conducted with WV amino acids at National Center for Biotechnology Information (NCBI) (https://www.ncbi.nlm.nih.gov/). Then, multiple sequence alignment was carried out by the Clustal Omega program (http://www.ebi.ac.uk/Tools/msa/clustalo/) and GENEDOC. Furthermore, the phylogenetic tree was constructed and tested by MEGA6.0 software based on the neighbor joining method.
RNA extraction and qRT-PCR analysis
Total RNA was isolated using Trizol reagent (Invitrogen, USA) and treated by DNase I to remove any genomic DNA contamination. cDNA was synthesized by using reverse transcriptase (Toyobo, Japan) according to the manufacturer’s protocol. And the quality of cDNA was detected by PCR using Ef1a primers as follows: Ef1a forward (5′-GGCCACAGGGATTTCATCAAG-3′) and reverse (5′-GTCCCTTGTACCAGTCGAGGTTG-3′). The concentration of each cDNA sample was adjusted to 100 ng/μL using a NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific). qRT-PCR was performed in 10 μL reactions, which consisted of 5 μL TransStart Green qRT-PCR SuperMix (TransGen Biotech), 0.5 μL of each primer, and 4 μL of first-strand cDNAs. The PCR reactions were carried out on the Roche LightCycler 480 II. The program was described as follows: preheated at 94 °C for 3 min, followed by 40 cycles of amplification (94 °C for 20 s, 58 °C for 20 s, and 72 °C for 30 s), and stopped by an extension (72 °C for 10 min). Actin was used as the internal control to normalize the relative expression level of each gene. The primers of actin were listed as follows: actin forward (5′-GTCCTCTTCCAGCCATCCAT-3′) and actin reverse (5′-ACCACTGAGCACAATGTTACCG-3′). For all qRT-PCR experiments were conducted in triplicate. All qRT-PCR primers of genes related to chlorophyll biosynthesis and photosynthesis are listed in Additional file 2: Table S3.
Glutamyl tRNA reductase
National Center for Biotechnology Information
Nuclear-encoded RNA polymerase
Open reading frame
Plastid-encoded RNA polymerase
Quantitative real time polymerase chain reaction
sodium dodecyl sulfate PAGE
tris buffered saline with Tween 20
Transmission electron microscopy
Tomato Genetics Resource Center
Gruissem W. Chloroplast gene expression: how plants turn their plastids on. Cell. 1989;56(2):161–70.
Chen X, Schnell DJ. Protein import into chloroplasts. Trends Cell Biol. 1999;9:222–7.
Fuks B, Schnell DJ. Mechanism of protein transport across the chloroplast envelope. Plant Physiol. 1997;114(2):405–10.
Kusumi K, Iba K. Establishment of the chloroplast genetic system in rice during early leaf development and at low temperatures. Front Plant Sci. 2014;5:386.
Armbruster U, Pesaresi P, Pribil M, Hertle A, Leister D. Update on chloroplast research: new tools, new topics, and new trends. Mol Plant. 2011;04(1):1–16.
Katz JJ, Norris JR, Shipman LL, Thurnauer MC, Wasielewski MR. Chlorophyll function in the photosynthetic reaction center. Annu Rev Biophys Bioeng. 1978;7(1):393–434.
Eckhardt U, Grimm B, Hörtensteiner S. Recent advances in chlorophyll biosynthesis and breakdown in higher plants. Plant Mol Biol. 2004;56(1):1–14.
McCormac AC, Fischer A, Kumar AM, Söll D, Terry MJ. Regulation of HEMA1 expression by phytochrome and a plastid signal during de-etiolation in Arabidopsis thaliana. Plant J. 2001;25(5):549–61.
Reinbothe C, Bakkouri ME, Buhr F, Muraki N, Nomata J, Kurisu G, Fujita Y, Reinbothe S. Chlorophyll biosynthesis: spotlight on protochlorophyllide reduction. Trends Plant Sci. 2010;15(11):614–24.
Wilson-Sánchez D, Rubio-Díaz S, Muñoz-Viana R, Pérez-Pérez JM, Jover-Gil S, Ponce MR, Micol JL. Leaf phenomics: a systematic reverse genetic screen for Arabidopsis leaf mutants. Plant J. 2014;79(5):878–91.
Dong H, Fei GL, Wu CY, Wu FQ, Sun YY, Chen MJ, Ren YL, Zhou KN, Cheng ZJ, Wang JL, Jiang L, Zhang X, Guo XP, Lei CL, Su N, Wang H, Wan JM. A rice virescent-yellow leaf mutant reveals new insights into the role and assembly of plastid caseinolytic protease in higher plants. Plant Physiol. 2013;162(4):1867–80.
Deng XJ, Zhang HQ, Wang Y, He F, Liu JL, Xiao X, Wang GL. Mapped clone and functional analysis of leaf-color gene ygl7 in a rice hybrid (Oryza sativa L. ssp. indica). PLoS One. 2014;9(6):e99564.
Falbel TG, Meehl JB, Staehelin JA. Severity of mutant phenotype in a series of chlorophyll-deficient wheat mutants depends on light intensity and the severity of the block in chlorophyll synthesis. Plant Physiol. 1996;112(2):821–32.
Sawers RJ, Viney J, Farmer PR, Bussey RR, Olsefski G, Anufrikova K, Hunter CN, Brutnell TP. The maize oil yellow1 (Oy1) gene encodes the I subunit of magnesium chelatase. Plant Mol Biol. 2006;60(1):95–106.
Li W, Tang S, Zhang S, Shan J, Tang C, Chen Q, Jia G, Han Y, Zhi H, Diao X. Gene mapping and functional analysis of the novel leaf color gene siygl1 in foxtail millet [setaria italica (l.) p. beauv]. Physiol Plant. 2016;157(1):24–37.
Huang X, Zhang X, Yang S. A novel chloroplast-localized protein EMB1303 is required for chloroplast development in Arabidopsis. Cell Res. 2009;19(10):1205–16.
Huang Y, Li H. Arabidopsis CHLI2 can substitute for CHLI1. Plant Physiol. 2009;150(2):636–45.
Grossman AR, Lohr M, Im CS. Chlamydomonas reinhardtii in the landscape of pigments. Annu Rev Genet. 2004;38(38):119–73.
Wu Z, Zhang X, He B, Diao L, Sheng S, Wang J, Guo X, Su N, Wang L, Jiang L, Wang C, Zhai H, Wan J. A chlorophyll-deficient rice mutant with impaired chlorophyllide esterification in chlorophyll biosynthesis. Plant Physiol. 2007;145(1):29–40.
Kusaba M, Ito H, Morita R, Iida S, Sato Y, Fujimoto M, Kawasaki S, Tanaka R, Hirochika H, Nishimura M, Tanaka A. Rice NON-YELLOW COLORING1 is involved in light-harvesting complex II and grana degradation during leaf senescence. Plant Cell. 2007;19(4):1362–75.
Ishizaki Y, Tsunoyama Y, Hatano K, Ando K, Kato K, Shinmyo A, Kobori M, Takeba G, Nakahira Y, Shiina T. A nuclear-encoded sigma factor, Arabidopsis SIG6, recognizes sigma- 70 type chloroplast promoters and regulates early chloroplast development in cotyledons. Plant J. 2005;42(2):133–44.
Fan M, Gao S, Ren J, Yang Q, Li H, Yang C, Ye Z. Overexpression of SIRBZ results in chlorosis and dwarfism through impairing chlorophyll, carotenoid, and gibberellin biosynthesis in tomato. Front Plant Sci. 2016;7:907.
Zhang H, Li J, Yoo JH, Yoo SC, Cho SH, Koh HJ, Seo HS, Paek NC. Rice Chlorina-1 and Chlorina-9 encode ChlD and ChlI subunits of mg-chelatase, a key enzyme for chlorophyll synthesis and chloroplast development. Plant Mol Biol. 2006;62(3):325–37.
Sugimoto H, Kusumi K, Tozawa Y, Yazaki J, Kishimoto N, Kikuchi S, Iba K. The virescent-2 mutation inhibits translation of plastid transcripts for the plastid genetic system at an early stage of chloroplast differentiation. Plant Cell Physiol. 2004;45(8):985–96.
Yoo SC, Cho SH, Sugimoto H, Li J, Kusumi K, Koh HJ, Iba K, Paek NC. Rice Virescent3 and Stripe1 encoding the large and small subunits of ribonucleotide reductase are required for chloroplast biogenesis during early leaf development. Plant Physiol. 2009;150(1):388–401.
Kusumi K, Sakata C, Nakamura T, Kawasaki S, Yoshimura A, Iba K. A plastid protein NUS1 is essential for build-up of the genetic system for early chloroplast development under cold stress conditions. Plant J. 2011;68(6):1039–50.
Chen T, Zhang Y, Zhao L, Zhu Z, Lin J, Zhang S, Wang C. Fine mapping and candidate gene analysis of a green-revertible albino gene gra(t) in rice. J Genet Genomics. 2009;36(2):117–23.
Sakamoto W, Uno Y, Zhang Q, Miura E, Kato Y. Sodmergen. Arrested differentiation of proplastids into chloroplasts in variegated leaves characterized by plastid ultrastructure and nucleoid morphology. Plant Cell Physiol. 2009;50(12):2069–83.
Tanksley SD, Ganal MW, Prince JP, de-Vicente MC, Bonierbale MW, Broun P, Fulton TM, Giovannoni JJ, Grandillo S, Martin GB, Messeguer R, Miller JC, Miller L, Paterson AH, Pineda O, Roder MS, Wing RA, Wu W, Young ND. High density molecular linkage maps of the tomato and potato genomes. Genetics. 1992;132(4):1141–60.
Yu Q, Ma Q, Kong M, Zhao T, Zhang X, Zhou Q, Huang C, Chong K, Yang Z. AtECB1/MRL7, a Thioredoxin-like fold protein with disulfide reductase activity, regulates chloroplast gene expression and chloroplast biogenesis in Arabidopsis thaliana. Mol Plant. 2014;7(1):206–17.
Yu F, Park SS, Liu X, Foudree A, Fu A, Powikrowska M, Khrouchtchova A, Jensen PE, Kriger JN, Gray GR, Rodermel SR. Suppressor of variegation4, a new var2 Suppressor locus, encodes a pioneer protein that is required for chloroplast biogenesis. Mol Plant. 2011;4(2):229–40.
Yu Q, Jiang Y, Chong K, Yang Z. AtECB2, a pentatricopeptide repeat protein, is required for chloroplast transcript accD RNA editing and early chloroplast biogenesis in Arabidopsis thaliana. Plant J. 2009;59(6):1011–23.
Taylor WC. Regulatory interactions between nuclear and plastid genomes. Annu Rev Plant Physiol Plant Mol Biol. 1989;40(40):211–33.
Chen M, Choi Y, Voytas DF, Rodermel S. Mutations in the Arabidopsis VAR2 locus cause leaf variegation due to the loss of a chloroplast FtsH protease. Plant J. 2000;22(4):303–13.
Sugimoto H, Kusumi K, Noguchi K, Yano M, Yoshimura A, Iba K. The rice nuclear gene, VIRESCENT 2, is essential for chloroplast development and encodes a novel type of guanylate kinase targeted to plastids and mitochondria. Plant J. 2007;52(3):512–27.
Hopkins WG, Walden DB. Temperature sensitivity of virescent mutants of maize. J Hered. 1977;68(5):283–6.
Millerd A, Mcwilliam JR. Studies on a maize mutant sensitive to low temperature I. influence of temperature and light on the production of chloroplast pigments. Plant Physiol. 1968;43(12):1967–72.
Xing S, Miao J, Li S, Qin G, Tang S, Li H, Gu H, Qu L. Disruption of the 1-deoxy-D-xylulose-5-phosphate reductoisomerase (DXR) gene results in albino, dwarf and defects in trichome initiation and stomata closure in Arabidopsis. Cell Res. 2010;20(6):688–700.
Robles P, Micol JL, Quesada V. Arabidopsis MDA1, a nuclear-encoded protein, functions in chloroplast development and abiotic stress responses. PLoS One. 2012;7(8):e42924.
Buchanan BB, Schürmann P, Wolosiuk RA, Jacquot JP. The ferredoxin/thioredoxin system: from discovery to molecular structures and beyond. Photosynth Res. 2002;73(1–3):215–22.
Bartsch S, Monnet J, Selbach K, Quigley F, Gray J, von Wettstein D, Reinbothe S, Reinbothe C. Three thioredoxin targets in the inner envelope membrane of chloroplasts function in protein import and chlorophyll metabolism. Proc Natl Acad Sci U S A. 2008;105(12):4933–8.
Cao Z, Yu Q, Sun Y, Lu Y, Cui Y, Yang Z. A point mutation in the pentatricopeptide repeat motif of the AtECB2 protein causes delayed chloroplast development. J Integr Plant Biol. 2011;53(4):258–69.
Rose AB. Intron-mediated regulation of gene expression. Curr Top Microbiol Immunol. 2008;326:277–90.
Cong B, Barrero LS, Tanksley SD. Regulatory change in YABBY-like transcription factor led to evolution of extreme fruit size during tomato domestication. Nat Genet. 2008;40(6):800–4.
Sridhar VV, Surendrarao A, Liu Z. APETALA1 and SEPALLATA3 interact with SEUSS to mediate transcription repression during flower development. Development. 2006;133(16):3159–66.
Eshed Y. ZamirD. A genomic library of Lycopersicon pennellii in L. esculentum: a tool for fine mapping of genes. Euphytica. 1994;79(3):175–9.
Inskeep WP, Bloom PR. Extinction coefficients of chlorophyll a and b in N, N-Dimethylformamide and 80% acetone. Plant Physiol. 1985;77(2):483–5.
Saghai-Maroof MA, Soliman KM, Jorgensen RA, Allard RW. Ribosomal DNA spacer-length polymorphisms in barley: mendelian inheritance, chromosomal location, and population dynamics. Proc Natl Acad Sci U S A. 1984;81(24):8014–8.
Lin T, Zhu G, Zhang J, Xu X, Yu Q, Zheng Z, Zhang Z, Lun Y, Li S, Wang X, Huang Z, Li J, Zhang C, Wang T, Zhang Y, Wang A, Zhang Y, Lin K, Li C, Xiong G, Xue Y, Mazzucato A, Causse M, Fei Z, Giovannoni JJ, Chetelat RT, Zamir D, Stadler T, Li J, Ye Z, Du Y, Huang S. Genomic analyses provide insights into the history of tomato breeding. Nat Genet. 2014;46(11):1220–6.
Lander ES, Green P, Abrahamson J, Barlow A, Daly MJ, Lincoln SE, Newburg L. Mapmaker: an interactive computer package for constructing primary genetic linkage maps of experimental and natural populations. Genomics. 1987;1(2):174–81.
Gao S, Gao Y, Xiong C, Yu G, Chang J, Yang Q, Yang C, Ye Z. The tomato B type cyclin gene SlCycB2, plays key roles in reproductive organ development, trichome initiation, terpenoids biosynthesis and Prodenia litura defense. Plant Sci. 2017;262:103–14.
We thank the Tomato Genetics Resource Center for providing the seed stocks of the mutant LA1526 and the introgression line IL2-3.
This work was supported by grants of the Natural Science Foundation of China (31872122 and 31572128), the National Key Research and Development Program (2017YFD0101902), the Foundation for the Author of National Excellent Doctoral Dissertation of China (No. 201468), and the China Agricultural Research System (CARS-23-A-03). The funding body played no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. Immunoblot analysis of chloroplast proteins in wv mutant and AC plants. These proteins were PsaD (the photosystem I subunits); PsbA (photosystem II reaction center subunit; 39 kDa); rbcL (the large subunit of Rubisco enzyme). M represented marker. The α-actin was used as a loading control. Figure S2. Alignment of CDS and genomic sequences of wv from AC and LA1526. Figure S3. qRT-PCR analysis of complemented WV in transgenic plants. T1–1,T1–2 and T1–3 represented three independent functional complemented transgenic lines. The expression level of WV in apical buds of LA1526 plants at 16 °C with 250 μmolm− 2 s− 1 light intensity was provided as controls.Values were mean ± SD of three technical replicates. Asterisks indicated statistical significance at P < 0.01. Figure S4. qRT-PCR analysis of Solyc02g079730.2.1 RNAi transgenic plants. The expression level of Solyc02g079730.2.1 in AC plants was provided as controls. Ri-1, Ri-2 and Ri-3 represented three independent RNAi transgenic lines. Each value represented the mean ± SE of three replicates. Asterisks indicated statistical significance at P < 0.01. (DOC 3295 kb)
Table S1. Primers for amplifying CDS and promoters of candidate genes. Table S2. Primers for Solyc02g079730.2.1 fusion vectors construction. Table S3. Primers for qRT-PCR of chloroplast and nucleus encoded photosynthetic-related genes. Table S4. Genetic segregation analysis of wv mutants in different generations. Table S5. Details of the markers for definitive mapping of wv. Table S6. Predicted genes between marker wv-c53 and wv-c75. (DOC 131 kb)