Plant Cell Advance Publication. Published on May 18, 2017, doi:10.1105/tpc.16.00484
RESEARCH ARTICLE
Carbonic Anhydrases Function in Anther Cell Differentiation Downstream of the Receptor-like Kinase EMS1 Jian Huang,a Zhiyong Li,a Gabriel Biener,b Erhui Xiong,a,c Shikha Malik,a Nathan Eaton,a Catherine Z. Zhao,d Valerica Raicu,a,b Hongzhi Kong,e and Dazhong Zhaoa,f,1 a
Department of Biological Sciences, University of Wisconsin-Milwaukee, Milwaukee, WI 53211, USA Department of Physics, University of Wisconsin-Milwaukee, Milwaukee, WI 53211, USA c College of Life Science, Henan Agricultural University, Zhenzhou 450002, China d Whitefish Bay High School, Whitefish Bay, WI 53217, USA e State Key Laboratory of Systematic and Evolutionary Botany, Institute of Botany, Chinese Academy of Sciences, Beijing 100093, China f College of Life Science, Shandong Normal University, Jinan 250014, China b
1
Corresponding Author:
[email protected].
Short title: Signaling Role of Carbonic Anhydrases One-sentence summary: β-carbonic anhydrases are post-translationally modified by EMS1 and act as direct downstream effectors of this receptor-like kinase, highlighting their crucial role in cell differentiation. The author responsible for distribution of materials integral to the findings presented in this article in accordance with the policy described in the Instructions for Authors (www.plantcell.org) is: Dazhong Zhao (
[email protected]).
ABSTRACT Plants extensively employ leucine-rich repeat receptor-like kinases (LRR-RLKs), the largest family of RLKs, to control a wide range of growth and developmental processes as well as defense responses. To date, only a few direct downstream effectors for LRR-RLKs have been identified. We previously showed that the LRR-RLK EMS1 (EXCESS MICROSPOROCYTES1) and its ligand TPD1 (TAPETUM DETERMINANT1) are required for the differentiation of somatic tapetal cells and reproductive microsporocytes during early anther development in Arabidopsis thaliana. Here, we report the identification of β-carbonic anhydrases (βCAs) as the direct downstream targets of EMS1. EMS1 biochemically interacts with βCA proteins. Loss-of-function of βCA genes caused defective tapetal cell differentiation, while overexpression of βCA1 led to the formation of extra tapetal cells. EMS1 phosphorylates βCA1 at four sites, resulting in increased βCA1 activity. Furthermore, phosphorylation-blocking mutations impaired the function of βCA1 in tapetal cell differentiation; however, a phosphorylation mimic mutation promoted the formation of tapetal cells. βCAs are also involved in pH regulation in tapetal cells. Our findings highlight the role of βCA in controlling cell differentiation and provide insights into the post-translational modification of carbonic anhydrases via receptor-like kinase-mediated phosphorylation.
1 ©2017 American Society of Plant Biologists. All Rights Reserved
1
INTRODUCTION
2 3
In flowering plants, successful sexual reproduction depends on the normal development of
4
anthers, which produce the male gametophyte, pollen. Within each of the four lobes
5
(microsporangia) of a mature anther, the central reproductive microsporocytes (or pollen mother
6
cells) are surrounded by four concentrically organized somatic cell layers: the epidermis,
7
endothecium, middle layer, and tapetum (from the outside to inside) (Goldberg et al., 1993; Scott
8
et al., 2004; Zhao, 2009; Walbot and Egger, 2016). Microsporocytes give rise to pollen via
9
meiosis, while the somatic cell layers, particularly the tapetum, are required for pollen
10
development and release. Due to the central importance of anthers for plant yield and breeding, it
11
is imperative to obtain an in-depth understanding of anther cell differentiation.
12
In Arabidopsis thaliana, the tapetum consists of a single layer of endopolyploid cells,
13
which enclose microsporocytes, tetrads, microspores, and developing pollen grains as anther
14
development progresses (Goldberg et al., 1993; Scott et al., 2004; Walbot and Egger, 2016).
15
Tapetum development comprises three stages: differentiation, maturation, and programmed cell
16
death (PCD). Early differentiated tapetal cells secrete enzymes required for the release of haploid
17
microspores from meiotic tetrads (Pacini et al., 1985; Clement and Pacini, 2001; Ishiguro et al.,
18
2001; Hsieh and Huang, 2007; Parish, 2010). During endoreduplication, mature tapetal cells,
19
characterized by multiple nuclei, non-photosynthetic plastids (elaioplasts), and tapetosomes, are
20
highly metabolically active, and thereby serve as a nutritional source to provide energy and
21
materials (polysaccharides, lipids, and proteins) for the development of microspores and pollen
22
grains. In the later stage, tapetal cells are degenerated via PCD and the resulting materials are
23
deposited onto pollen to form the pollen coat (tryphine). Mutants lacking tapetum or plants with
24
genetically ablated tapetum fail to produce pollen grains (Mariani et al., 1990; Zhao et al., 2002;
25
Yang et al., 2003; Huang et al., 2016b). In addition, precocious or delayed tapetum degeneration
26
causes abnormal pollen development (Zhang et al., 2006; Zhang and Yang, 2014). Furthermore,
27
stress-induced male sterility is mainly ascribed to alterations in tapetum development (Parish et
28
al., 2012; De Storme and Geelen, 2014; Smith and Zhao, 2016). Although the tapetum is
29
essential for pollen development, the molecular mechanisms underlying tapetal cell
30
differentiation and functional maintenance are not clear.
31
We previously demonstrated that the EMS1 (EXCESS MICROSPOROCYTES1) 2
32
leucine-rich repeat receptor-like kinase (LRR-RLK)-linked signaling pathway plays a
33
fundamental role in the differentiation of somatic tapetal cells and reproductive microsporocytes
34
during early anther development. In Arabidopsis, anthers in ems1 [also known as exs (extra
35
sporogenous cells)] and tpd1 (tapetum determinant1) mutants have no tapetal cells; instead, they
36
produce excess microsporocytes (Canales et al., 2002; Zhao et al., 2002; Yang et al., 2003). We
37
identified TPD1 as a putative small protein ligand of EMS1 (Jia et al., 2008). Our results further
38
showed that TPD1 is secreted from precursors of microsporocytes and then activates EMS1,
39
which is localized at the plasma membrane of tapetal precursor cells/tapetal cells. The
40
TPD1-EMS1 signaling pathway initially promotes periclinal division of parietal cells to form
41
tapetal precursor cells, and later determines and maintains the fate of functional tapetal cells.
42
Additionally, tapetal cells suppress microsporocyte proliferation (Huang et al., 2016c). We
43
recently also found that the SERK1/2 (SOMATIC EMBRYOGENESIS RECEPTOR-LIKE
44
KINASE) LRR-RLK acts as a potential co-receptor of EMS1 (Li et al., 2017). However, the
45
downstream signaling molecules of TPD1-EMS1/SERK1/2 are currently unclear.
46
In this study, we identified β-carbonic anhydrases (βCAs) as EMS1-interacting proteins.
47
CAs catalyze the rapid, reversible reaction CO2+H2O↔HCO3‾+H+ (Tripp et al., 2001; Alterio et
48
al., 2009; Becker et al., 2011). In animals, CAs are important for various processes in normal and
49
especially pathological states, such as gluconeogenesis, lipogenesis, ureagenesis, and tumor
50
formation (Chegwidden et al., 2000; Supuran, 2008; Davis et al., 2010). In plants, CAs
51
principally belong to three independently evolved families: α, β, and γ, with the β family being
52
predominant (Badger, 2003; Rowlett, 2014). βCAs function in photosynthesis mainly by
53
concentrating CO2 and facilitating the movement of CO2 and HCO3–. Arabidopsis has six βCAs
54
(βCA1 to βCA6), and each βCA has various splice forms (Wang et al., 2014). βCA1
55
loss-of-function mutants show reduced seedling establishment (Ferreira et al., 2008). βCA1 and
56
βCA4 coordinate with EPIDERMAL PATTERNING FACTORS2 (EPF2) and CO2 RESPONSE
57
SECRETED PROTEASE (CRSP) to control stomatal development and movement (Hu et al.,
58
2010; Engineer et al., 2014). Furthermore, βCA4 and its interactor protein PIP2;1 aquaporin are
59
required for CO2-induced stomatal movement (Wang et al., 2016a). A recent study demonstrated
60
that βCA2 and βCA4 are important for optimal plant growth via affecting amino acid
61
biosynthesis but not photosynthesis (DiMario et al., 2016). Thus, βCAs may possess many
62
unidentified non-photosynthetic roles in plant growth and development. 3
63
In this study, we showed that βCA1, βCA2 and βCA4 biochemically interact with EMS1.
64
Loss-of-function of βCA genes led to abnormal tapetal cell differentiation, whereas
65
overexpression of βCA1 caused the formation of extra tapetal cells. EMS1 phosphorylates βCA1,
66
βCA2 and βCA4. The phosphorylation of βCA1 significantly increases its activity. Moreover,
67
phosphorylation-blocking mutations caused the failure of βCA1 to recover tapetal cell
68
differentiation in the βca1 βca2 βca4 mutant; however, a phosphorylation mimic mutation
69
promoted the formation of tapetal cells. βCAs likely regulate tapetal cell pH. Our results suggest
70
that βCAs serve as the direct downstream targets of EMS1, which highlights a role for βCAs in
71
controlling cell differentiation and provides a paradigm for LRR-RLK-linked signal transduction
72
pathways. Furthermore, our research sheds light on the post-translational modification of CA via
73
receptor-like kinase-mediated phosphorylation.
74 75 76
RESULTS
77 78
EMS1 Interacts With βCAs
79 80
To identify the downstream targets of EMS1, we performed yeast two-hybrid (Y2H) screening
81
for potential EMS1-interacting proteins using the EMS1 kinase domain (852-1192 amino acids)
82
as the bait (Figure 1). We generated a cDNA library using young buds with stage-5 and -6
83
anthers, in which EMS1 expression is at its peak. Out of 129 selected positive clones, 79 cDNAs
84
encoded β-carbonic anhydrase 1 (βCA1, At3g01500). Three splice variants for βCA1 were found
85
in The Arabidopsis Information Resource (TAIR). Our Y2H screening only identified βCA1.3
86
(16/79) and a new splice variant (63/79) which lacks the first and last exon, named βCA1.4
87
(Supplemental Figure 1). Y2H tests showed that βCA1.3 and βCA1.4 interacted with EMS1 at a
88
higher level than βCA1.1 and βCA1.2 (Supplemental Figure 2). Phylogenetic analysis
89
demonstrated that βCA1 and βCA2 share the highest similarity and are closely related to βCA3
90
and βCA4, whereas βCA5 and βCA6 form a different clade (Supplemental Figure 3,
91
Supplemental File 1). All βCA genes except βCA3 have different numbers of splice variants
92
(TAIR). Further Y2H assays showed that βCA2.2, βCA3, βCA4.1, βCA5.2, and βCA6.1, which
93
have the closest structures to βCA1.4 compared with other splice forms, also interacted with 4
94 95
EMS1, possibly at different levels (Figure 1A). We then examined interactions between EMS1 and βCAs by performing a bimolecular
96
fluorescence complementation (BiFC) assay with Arabidopsis mesophyll protoplasts. EMS1 and
97
βCAs were fused to the amino- and carboxy-terminal halves of Enhanced Yellow Fluorescent
98
Protein (EYFP), respectively. EYFP signals were observed at the plasma membrane when
99
EMS1-cEYFP was co-transfected with βCA1.4-nEYFP (Figure 1B), βCA2.2-nEYFP (Figure
100
1C), and βCA4.1-nEYFP (Figure 1E); however, no EYFP signal was observed at the plasma
101
membrane when EMS1-cEYFP was co-transfected with βCA3-nEYFP (Figure 1D),
102
βCA5.2-nEYFP (Figure 1F), or βCA6.1-nEYFP (Figure 1G). Protein gel blot analysis showed
103
EMS1-cEYFP and βCAs-nEYFP were expressed at similar levels in the BiFC assay
104
(Supplemental Figure 4). As controls, no EYFP signal was observed from EMS1-cEYFP and
105
αCA1.1-nEYFP or BRI1-cEYFP and βCA1.4-nEYFP combinations (Figure 1H and 1I). Our
106
results support the notion that EMS1 specifically interacts with βCA1.4, βCA2.2 and βCA4.1 at
107
the plasma membrane in living cells.
108
We also performed a Förster Resonance Energy Transfer (FRET) experiment to quantify
109
the EMS1-βCA1.4 interaction using a two-photon micro-spectroscope (Raicu et al., 2009; Biener
110
et al., 2014). FRET efficiencies were calculated using images generated by unmixing measured
111
spectra of the donor (EMS1-CFP) at 860 nm (Figure 1J) and acceptor (βCA1.4-EYFP) following
112
excitation at 860 nm (Figure 1K), as well as 960 nm (Figure 1L). A FRET efficiency of 13±4%
113
(Supplemental Table 1) was obtained from 29 regions of interest at the plasma membrane,
114
indicating that the EMS1-βCA1.4 interaction occurs at the plasma membrane.
115
Finally, we carried out co-immunoprecipitation (co-IP) experiments to examine whether
116
EMS1 interacts with βCAs by representatively testing βCA1 in planta. Membrane proteins
117
extracted from young buds of ProEMS1:EMS1-4xcMyc ProβCA1:βCA1-Flag double transgenic
118
plants, where both EMS1 and βCA1 were expressed in early anthers, were immunoprecipitated
119
with an anti-cMyc antibody. Protein gel blot analysis using anti-cMyc and anti-Flag antibodies
120
against the immunoprecipitated proteins showed that βCA1-Flag was co-immunoprecipitated by
121
EMS1-4xcMyc (Figure 1M), indicating that EMS1 interacts with βCA1 in planta. To test
122
whether the EMS1 ligand TPD1 affects the EMS1 and βCA1 interaction, we performed co-IP
123
using a protoplast transient expression system. βCA1.4-Flag was expressed in protoplasts
124
prepared from the leaves of Pro35S:EMS1-cMyc and Pro35S:EMS1-cMyc 5
125
Pro35S:TPD1p-GFP-∆TPD1 transgenic plants (Huang et al., 2016c). Protein gel blot analysis
126
using an anti-Flag antibody against membrane proteins immunoprecipitated by anti-cMyc
127
detected a similar level of βCA1.4-Flag in the presence and absence of TPD1 (Supplemental
128
Figure 5), suggesting that TPD1 is not essential for the interaction between βCA1.4 and EMS1.
129
Collectively, our results show that EMS1 biochemically interacts with βCA1, βCA2 and
130
βCA4, suggesting that βCA1, βCA2 and βCA4 may act as downstream signaling effectors of
131
EMS1.
132 133 134
βCAs Are Expressed in Tapetal Cells
135 136
To investigate whether βCAs serve as EMS1 downstream signaling molecules, we examined the
137
expression of βCAs in anthers and tapetal cells (Figure 2). The βCA1 gene has 10 exons, which
138
leads to four alternative splice variants (Supplemental Figure 1). Our RT-PCR results showed
139
that the expression of βCA1.3 was dominant in leaf, seedling, silique and stem tissues, while
140
βCA1.4 was only weakly expressed in seedlings (Supplemental Figure 6). All βCA1 variants
141
except βCA1.1 were expressed in young buds (Figure 2A); however, only βCA1.3 and βCA1.4
142
were detected in wild-type stage-5/6 anthers (Figure 2A). In ems1 stage-5/6 anthers, the
143
expression of βCA1.3 was not altered compared to wild type. Conversely, the expression of
144
βCA1.4 was not detected (Figure 2A), suggesting that βCA1.4 might be a major downstream
145
molecule of EMS1. We also detected the expression of βCA2 and βCA4 in wild-type stage-5/6
146
anthers, and their expression was reduced in ems1 stage-5/6 anthers (Figure 2B). We did not
147
observe the expression of βCA3 in wild-type or ems1 stage-5/6 anthers (Figure 2B). Our results
148
suggest that βCA1 (mainly βCA1.4), together with βCA2 and βCA4, are potential downstream
149
molecules of EMS1.
150
Using the βCA1 1.4 kb promoter region and transcribed region without the last intron and
151
exon, we generated ProβCA1:βCA1-GFP transgenic plants expressing βCA1.3-GFP and
152
βCA1.4-GFP proteins. We observed a weak GFP signal in epidermal cells in stage-4 anthers
153
(Figure 2C). At stage 5 and stage 6, two key stages for tapetal cell differentiation, GFP signals
154
strongly accumulated in tapetal cells (Figure 2D and 2E); however, GFP signals became
155
gradually reduced in tapetal cells of stage-7 (Figure 2F) and stage-8 (Figure 2G) anthers. GFP 6
156
signals were not detectable in tapetal cells at stage 10, when tapetal cell degeneration initiates
157
(Figure 2H). Furthermore, we found that βCA1 was localized at the plasma membrane and in the
158
cytoplasm of tapetal cells in stage-5 anthers (Figure 2I to 2K). We also generated
159
ProβCA2:βCA2-GFP, ProβCA3:βCA3-GFP and ProβCA4:βCA4-GFP transgenic plants. Our
160
results showed that instead of βCA3 (Figure 2M), βCA2 (Figure 2L) and βCA4 (Figure 2N) were
161
detected in both tapetal cells and microsporocytes of stage-5 anthers.
162
Taken together, similar to EMS1 (Huang et al., 2016c), βCA1, βCA2 and βCA4 are
163
localized in tapetal cells, supporting the notion that βCA1, βCA2, and βCA4 are downstream
164
signaling partners of EMS1.
165 166 167
βCAs Are Required for Tapetal Cell Differentiation
168 169
To investigate the function of βCAs in anther cell differentiation, we analyzed the phenotypes of
170
βCAs loss-of-function mutants (Figure 3). We did not detect mutant phenotypes in βca1
171
(SALK_106570), βca2 (CS303346: identified in this study, Supplemental Figure 7), βca3
172
(Salk_144106) or βca4 (CS859392) single mutant anthers, nor in βca1 βca2, βca1 βca3, βca1
173
βca4, βca2 βca4 double or βca1 βca2 βca3 triple mutant anthers. Compared to wild-type plants
174
(Figure 3A and 3H), βca1 βca2 βca4 triple mutant plants were smaller and did not produce
175
pollen grains in anthers (Figure 3B and 3I). Employing the artificial miRNA approach (Schwab
176
et al., 2006), we used the 35S promoter to knock down βCA1 to βCA4 genes. Among the 80
177
Pro35S:amirβCA1-4 transgenic plants examined, 52.5% (42/80) of plants were small and only
178
formed a few dead pollen grains (Figure 3C and 3J). To rule out the possibility that the pollen
179
production defect was caused by abnormal vegetative growth, we specifically knocked down
180
βCA1 to βCA4 in tapetal cells using the tapetum-specific promoter A9 (Paul et al., 1992; Feng
181
and Dickinson, 2010). Among the 90 examined ProA9:amirβCA1-4 transgenic plants, all of
182
which showed normal vegetative growth (Figure 3D), 68.9% (62/90) of plants produced
183
completely empty anthers (Figure 3K), indicating that βCAs are required for pollen formation.
184
To confirm that βCAs are responsible for pollen development, we conducted
185
complementation experiments. Our results showed that 64.0% (16/25) of ProβCA1:βCA1/βca1
186
βca2 βca4 plants (Figure 3E), 60.6% (20/33) of ProβCA2:βCA2/βca1 βca2 βca4 plants 7
187
(Supplemental Figure 8A) and 55.0% (22/40) of ProβCA4:βCA4/βca1 βca2 βca4 (Supplemental
188
Figure 8B) plants had normal growth and development. Although some plants were still smaller
189
than the wild type, 80.0% (20/25) of ProβCA1:βCA1/βca1 βca2 βca4 plants (Figure 3L), 75.8%
190
(25/33) of ProβCA2:βCA2/βca1 βca2 βca4 plants (Supplemental Figure 8C) and 82.5% (33/40)
191
of ProβCA4:βCA4/βca1 βca2 βca4 plants (Supplemental Figure 8D) produced normal pollen
192
grains. Moreover, although 70.0% (14/20) of ProA9:βCA1.4/βca1 βca2 βca4 plants were similar
193
to βca1 βca2 βca4 plants in terms of vegetative growth, their seed production was restored
194
(Figure 3F). Further analysis revealed that these plants produced normal pollen grains (Figure
195
3M). By contrast, 100% (22/22) of ProA9:βCA1.3/βca1 βca2 βca4 plants exhibited short siliques
196
(Figure 3G) and had no pollen grains (Figure 3N), suggesting that βCA1.4 but not βCA1.3 is
197
mainly responsible for early anther development. In conclusion, our results support the notion
198
that the disruption of βCA1, βCA2 and βCA4 caused the failure of pollen formation.
199
To further investigate the function of βCAs in anther development, we analyzed anther
200
cell differentiation in semi-thin sections (Figure 4). Our results showed that at stage 6, wild-type
201
anther lobes contained four somatic cell layers (epidermis, endothecium, the middle layer and
202
tapetum) and microsporocytes in the center (Figure 4A); however, tapetal-like cells were
203
vacuolated in βca1 βca2 βca4 mutant anthers (Figure 4B). Similar defects were observed in
204
anthers of Pro35S:amirβCA1-4 (Figure 4C) and ProA9:amirβCA1-4 (Figure 4D) plants in which
205
βCA1 to βCA4 were knocked down, with defects observed throughout the plant and specifically
206
in the tapetum, respectively. At stage-7 wild-type anthers, tetrads had formed (Figure 4F). By
207
contrast, in βca1 βca2 βca4 (Figure 4G), Pro35S:amirβCA1-4 (Figure 4H) and
208
ProA9:amirβCA1-4 (Figure 4I) anthers, tetrads had not formed. Instead, tapetal-like cells
209
continuously expanded and microsporocytes were degenerating. In stage-9 wild-type anthers,
210
tapetal cells were still present and the microspore wall was becoming thickened, indicating
211
normal pollen development (Figure 4K). Conversely, in βca1 βca2 βca4 (Figure 4L),
212
Pro35S:amirβCA1-4 (Figure 4M) and ProA9:amirβCA1-4 (Figure 4N) anthers, both tapetal-like
213
cells and microsporocytes were degenerated, resulting in empty anther lobes. In
214
ProA9:βCA1.4/βca1 βca2 βca4 plants, anther cell differentiation was the same as that of
215
wild-type plants (Figure 4E, 4J and 4O).
216 217
We then examined the expression of A9, a tapetum-specific marker gene, via in situ hybridization and qRT-PCR. At stage 6, the A9 gene was strongly expressed in tapetal cells in 8
218
wild-type anthers (Figure 4P and 4Y), but the expression levels of A9 were significantly reduced
219
in tapetal-like cells in βca1 βca2 βca4 (Figure 4Q and 4Y), Pro35S:amirβCA1-4 (Figure 4R and
220
4Y) and ProA9:amirβCA1-4 (Figure 4S and 4Y) anthers. As expected, A9 expression was
221
normal in ProA9:βCA1.4/βca1 βca2 βca4 tapetal cells (Figure 4T).
222
Elaioplasts (non-photosynthetic plastids) and tapetosomes, which are tapetum-specific
223
lipid-accumulating organelles, serve as nutrient sources to provide energy and materials
224
(polysaccharides, lipids, and proteins) for the development of microspores and pollen grains
225
(Dickinson, 1973; Owen and Makaroff, 1995; Clement and Pacini, 2001; Hsieh and Huang,
226
2007). Thus, we introduced the elaioplast marker FIB1a-GFP and the tapetosome marker
227
GRP17-GFP (Suzuki et al., 2013) into ProA9:amirβCA1-4 plants. Confocal microscopy analysis
228
revealed strong FIB1a-GFP signals in tapetal cells of stage-8 FIB1a-GFP (wild-type) anthers
229
(Figure 4U), but in the anthers from the 12 ProA9:amirβCA1-4 FIB1a-GFP plants examined, the
230
number and size of FIB1a-GFP signals were strongly reduced (Figure 4V). Similarly, compared
231
to the wild type (Figure 4W), GRP17-GFP signals were strongly reduced in stage-10 anthers
232
from the 17 analyzed ProA9:amirβCA1-4 GRP17-GFP plants (Figure 4X). Our results indicate
233
that βCAs are involved in the formation of elaioplasts and tapetosomes during tapetal cell
234
differentiation.
235
In summary, tapetal cell differentiation is impaired in βCA loss-of-function mutants,
236
resulting in the early degeneration of tapetal cells and consequently failed pollen formation. Thus,
237
our results support the notion that βCAs are required for tapetal cell differentiation.
238 239 240
Overexpression of βCA1 Leads to Extra Tapetal Cells
241 242
To further investigate whether βCAs are essential for tapetal cell differentiation, we
243
overexpressed the βCA1 gene by generating Pro4x35S-βCA1:βCA1 transgenic plants using four
244
cauliflower mosaic virus (CaMV) 35S enhancers (Weigel et al., 2000) and the βCA1 genomic
245
fragment (Figure 5). Among the 178 transgenic plants analyzed, 25.8% (46/178) of plants
246
showed short siliques (Supplemental Figure 9A and 9B), and a majority of pollen grains were
247
dead in the anthers (Supplemental Figure 9D and 9E). Moreover, 10.1% (18/178) of plants were
248
completely sterile and did not produce any pollen grains (Supplemental Figure 9C and 9F). 9
249
Further analysis of semi-thin sections showed that the 4x35S-ProβCA1:βCA1 anthers formed
250
extra tapetal-like cells at stage 6 (Figure 5A, 5B and 5O). At stage 7, extra vacuolated
251
tapetal-like cells were observed in Pro4x35S-βCA1:βCA1 anthers (Figure 5C, 5D and 5P).
252
However, when we introduced Pro4x35S-βCA1:βCA1 into the ems1 mutant,
253
Pro4x35S-βCA1:βCA1/ems1 anthers had the same phenotype as ems1 mutant anthers at stage 6
254
(Figure 5E and 5F), suggesting that the βCA1 gain-of-function effect is dependent on the normal
255
functioning of EMS1.
256
To confirm the cellular identity of the observed extra tapetal-like cells, we introduced
257
ProA9:mGFP5er, a tapetal cell marker, into Pro4x35S-βCA1:βCA1 plants. Confocal microscopy
258
analysis showed that the GFP signal was observed in a single layer of tapetal cells in wild-type
259
anthers at stage 6 (Figure 5G) and 7 (Figure 5I); however, in Pro4x35S-βCA1:βCA1 anthers,
260
additional cells expressed ProA9:mGFP5er at stage 6 (Figure 5H) and 7 (Figure 5J). Moreover,
261
the in situ hybridization results, which agreed with the results of confocal microscopy, showed
262
stronger and expanded expression of A9 in Pro4x35S-βCA1:βCA1 anthers (Figure 5L and 5N)
263
compared with wild type (Figure 5K and 5M). Furthermore, qRT-PCR showed that the
264
expression of A9 (Figure 5Q) and another tapetum-specific gene, ATA7 (Figure 5R), was
265
significantly higher in Pro4x35S-βCA1:βCA1 anthers compared to wild type. In conclusion, our
266
results suggest that overexpression of βCA1 promotes tapetal cell formation.
267 268 269
βCA1 is Phosphorylated by EMS1
270 271
Based on the well-established paradigm for the action of receptor-like kinases, EMS1 should
272
phosphorylate its interacting proteins. Therefore, we tested whether βCA1 can be phosphorylated
273
by EMS1 (Figure 6). Our in vitro phosphorylation assay showed that βCA1.4 was strongly
274
phosphorylated by EMS1-KD-GST (kinase domain) (Figure 6A), as was βCA2.2 and βCA4.1
275
(Supplemental Figure 10). We then performed protein gel blot analysis to investigate the in vivo
276
phosphorylation of βCA1. Membrane proteins extracted from young ProβCA1:βCA1-GFP and
277
ProβCA1:βCA1-GFP/ems1 buds were immunoprecipitated with an anti-GFP antibody. We
278
detected a weak signal from young ProβCA1:βCA1-GFP/ems1 buds but a strong signal from
279
young ProβCA1:βCA1-GFP buds using an anti-phospho-(Ser/Thr) antibody (Figure 6B). To 10
280
further confirm whether EMS1 can phosphorylate βCA1, we transiently expressed βCA1.4-EYFP
281
in protoplasts prepared from wild-type and Pro35S:EMS1 Pro35S:TPD1 leaves (Huang et al.,
282
2016c). The phosphorylation of βCA1.4 was detected in the presence of EMS1 and TPD1
283
(Figure 6C). Collectively, our data suggest that βCA1 can be phosphorylated by EMS1 both in
284
vitro and in vivo.
285 286 287
Phosphorylation of βCA1 by EMS1 Enhances Its Activity
288 289
We examined how phosphorylation mediated by EMS1 affects βCA enzyme activity (Figure 7).
290
Four phosphorylation sites (Thr35, Thr54, Thr69 and Ser189) were identified by mass
291
spectrometry (Figure 7A; Supplemental Figure 11; Supplemental Table 2). Thr35, Thr54 and
292
Thr69 are upstream of the β_CA_cladeB domain (NCBI Accession: cd00884; the CA domain for
293
short), while only Ser189 is located in the CA domain (Figure 7A). Although βCA2.2 and
294
βCA4.1 can also be phosphorylated by EMS1, only Thr35 is conserved (Supplemental Figure
295
11).
296
We found that EMS1 significantly enhanced βCA1.4 activity (Figure 7B to 7E). To
297
further analyze the significance of identified phosphorylation sites for the βCA1 function, we
298
produced mutated βCA1.4 in which Thr35, Thr54, Thr69 and Ser189 were substituted by alanine
299
(A) to block Thr (T)/Ser (S) phosphorylation and by aspartate (D) to functionally mimic
300
phosphorylated Thr (T)/Ser (S). βCA1.4T35A/D lost CA enzyme activity, and EMS1 treatment did
301
not increase its activity (Figure 7B). By contrast, the activities of βCA1.4T54A/D and βCA1.4T69A/D
302
were similar to that of βCA1.4, which were significantly enhanced by EMS1 treatment, although
303
they were lower than that of EMS1-treated wild-type βCA1.4 (Figure 7C and 7D). The activity
304
of βCA1.4S189A was similar to that of βCA1.4, but EMS1 treatment did not increase its activity
305
(Figure 7E). Interestingly, βCA1.4S189D activity was significantly higher than that of wild type,
306
and EMS1 treatment increased its activity more strongly than that of EMS1-treated βCA1.4
307
(Figure 7E).
308
To test the in vivo effect of EMS1 signaling on CA activity, we measured CA activity in
309
Pro35S:EMS1 Pro35S:TPD1 plants. Our results showed that CA activity was significantly
310
higher in both leaves and young buds of Pro35S:EMS1 Pro35S:TPD1 plants compared to wild 11
311
type (Figure 7F).
312
In summary, our results suggest that the phosphorylation of βCA1 by the kinase EMS1
313
increases its enzyme activity. Phosphorylation of Thr54, Thr69 and Ser189 likely contribute to
314
the enhancement of βCA1 activity, although Ser189 is more critical.
315 316 317
Phosphorylation of βCA1 Is Not Required for Its Subcellular Localization, Dimerization or
318
Interaction with EMS1
319 320
We investigated whether the phosphorylation of βCA1 affects its subcellular localization,
321
dimerization or interaction with EMS1 (Figure 8). The subcellular localization of a protein is
322
strongly associated with its function. βCA1 is localized to chloroplasts and the vicinity of the
323
plasma membrane (Fabre et al., 2007; Hu et al., 2010; Hu et al., 2015). We transiently expressed
324
βCA1.3-EYFP in Arabidopsis leaf protoplasts and detected EYFP signals in chloroplasts and at
325
the plasma membrane (Figure 8A and 8B), whereas βCA1.4-EYFP was localized at the plasma
326
membrane and in the cytoplasm (Figure 8C). We did not detect changes in the localization of
327
βCA1.4T35A (Figure 8D) or βCA1.4S189A (Figure 8E).
328
Studies on the structure of βCAs have revealed that the βCA dimer serves as the
329
functional unit of βCAs (Kimber and Pai, 2000; Strop et al., 2001). Our BiFC assay results
330
support the formation of homodimers of βCA1.4 at the plasma membrane (Figure 8F). Similarly,
331
we found dimerization between βCA1.4 and βCA1.4T35A (Figure 8G), βCA1.4T35A and
332
βCA1.4T35A (Figure 8H), βCA1.4 and βCA1.4 S189A (Figure 8I), as well as βCA1.4S189A and
333
βCA1.4 S189A (Figure 8J). In addition, βCA1.4T35A and βCA1.4S189A were still able to interact with
334
EMS1 (Figure 8K and 8L). Together, our results indicate that mutation of Thr35A or Ser189A
335
does not affect βCA1.4 localization, dimerization, or interaction with EMS1. Therefore, the
336
phosphorylation of βCA1 by EMS1 appears to mainly affect the enzyme activity of βCA1, which
337
might be required for its function in tapetal cell differentiation.
338 339 340
Phosphorylation of βCA1 by EMS1 Is Important for Tapetal Cell Differentiation
341
12
342
To investigate the functional significance of βCA1 phosphorylation in anther development, we
343
used ProA9:βCA1.4T35A, ProA9:βCA1.4T54A, ProA9:βCA1.4T69A and ProA9:βCA1.4S189A to
344
complement the βca1 βca2 βca4 phenotype (Figure 9). Seventy percent (14/20) of
345
ProA9:βCA1.4/βca1 βca2 βca4 plants (Figure 3F), 60.0% (12/20) of ProA9:βCA1.4T54A/βca1
346
βca2 βca4 plants (Figure 9A and 9D) and 63.3% (19/30) of ProA9:βCA1.4T69A/βca1 βca2 βca4
347
plants (Figure 9A and 9E) produced seeds. Similar to the wild type (Figure 9G and 9M), all of
348
these complemented plants showed viable pollen grains (Figure 3M; Figure 9J and 9K) and
349
normal tapetal cell differentiation (Figure 4E; Figure 9P and 9Q). However, all examined
350
ProA9:βCA1.4T35A/βca1 βca2 βca4 (13 total) and ProA9:βCA1.4S189A/βca1 βca2 βca4 (15 total)
351
plants showed the same defects in seed production (Figure 9B, 9C and 9F), pollen viability
352
(Figure 9H, 9I and 9L) and tapetal cell differentiation (Figure 9N, 9O and 9R) as those of βca1
353
βca2 βca4 plants. The expression levels of these mutated βCA1.4 transgenes were similar to that
354
of the βCA1.4 transgene in the βca1 βca2 βca4 background (Supplemental Figure 12). Therefore,
355
our results suggest that Thr35 and Ser189 are critical for the function of βCA1.4 in anther cell
356
differentiation.
357
Our in vitro studies showed that the enzyme activity of βCA1.4S189D was significantly
358
increased compared to wild type in both the absence and presence of EMS1 (Figure 7E). Thus,
359
we further studied the effect of phosphorylation of βCA1.4 on tapetal cell differentiation by
360
generating ProA9:βCA1.4S189D plants (Figure 10). Forty percent (32/80) of ProA9:βCA1.4S189D
361
plants showed reduced fertility (indicated by shorter siliques) compared with wild type (Figure
362
10A and 10B), as well as a reduced number of viable pollen grains (Figure 10D and 10E).
363
Moreover, 15.0% (12/80) of plants were completely sterile (Figure 10C) and did not produce
364
pollen grains (Figure 10F). Like the Pro4x35S-βCA1:βCA1 lines, anther sections from
365
completely sterile plants revealed an increased number of tapetal cells at stage 7 (Figure 10G to
366
10I). The expression levels of βCA1.4S189D and βCA1.4 in the transgenic plants were similar
367
(Supplemental Figure 13). Thus, our results suggest that phosphorylation of βCA1 by EMS1 is
368
essential for tapetal cell differentiation.
369 370 371
βCAs Regulate the pH of Tapetal Cells
372
13
373
To examine the possible mechanism by which βCAs control anther cell differentiation, we
374
investigated the pH of tapetal cells using the Carboxy SNARF®-1 pH indicator (Zhang et al.,
375
2001; Leshem et al., 2006; Sano et al., 2009). Tapetal cells had a lower pH than epidermal cells
376
in wild-type anthers (Figure 11). In βca1 βca2 βca4 anthers, the pH of epidermal cells was not
377
significantly altered; however, the pH of tapetal cells was significantly reduced compared to wild
378
type (Figure 11). Our findings suggest that βCAs might control tapetal cell differentiation by
379
regulating the pH of tapetal cells.
380
Taken together, our results provide several lines of evidence supporting the notion that
381
β-carbonic anhydrases serve as the direct downstream signaling molecules of the receptor-like
382
kinase EMS1 to control anther cell differentiation in Arabidopsis.
383 384 385
DISCUSSION
386 387
In this study, we identified βCAs as direct downstream players of the LRR-RLK EMS1 to
388
control anther cell differentiation in Arabidopsis. LRR-RLKs, with 223 members in Arabidopsis
389
(McCarty and Chory, 2000; Shiu et al., 2004; Torii, 2004; Zhao, 2009; Ma et al., 2016; Li et al.,
390
2017), are involved in various growth and developmental processes as well as defense responses,
391
including the regulation of shoot and root meristem sizes (Clark et al., 1997; Bommert et al.,
392
2013; Shinohara et al., 2016), cell fate determination and patterning (Canales et al., 2002; Zhao
393
et al., 2002; Kwak et al., 2005; Shpak et al., 2005; Jia et al., 2008), steroid hormone signaling (Li
394
and Chory, 1997; Li et al., 2002; Nam and Li, 2002; Hothorn et al., 2011; She et al., 2011;
395
Santiago et al., 2013), vascular patterning (Clay and Nelson, 2002; Fisher and Turner, 2007),
396
organ size and shape regulation (Torii et al., 1996; Xu et al., 2008), organ abscission (Jinn et al.,
397
2000; Leslie et al., 2010; Kumpf et al., 2013), pollen tube reception (Takeuchi and Higashiyama,
398
2016; Wang et al., 2016b), defense responses (Song et al., 1995; Gomez-Gomez and Boller,
399
2000; Lee et al., 2011; Sun et al., 2013; Halter et al., 2014; Zorzatto et al., 2015), plant
400
transpiration (Masle et al., 2005), nodulation (Endre et al., 2002; Searle et al., 2003), and
401
nitrogen acquisition (Tabata et al., 2014). To date, only a few direct downstream signaling
402
molecules have been identified for an increasing number of functionally investigated
403
LRR-RLKs. 14
Brassinosteroid (BR) signaling is one of the best-understood LRR-RLK-linked signal
404 405
transduction pathways in plants. Unlike animals, which use nuclear receptors to perceive steroids,
406
plants employ plasma membrane LRR-RLKs, the BRI1 receptor and its four redundant
407
co-receptors, SERK1-4, for BR signaling (Li et al., 2002; Nam and Li, 2002; Wang et al., 2008;
408
Hothorn et al., 2011; She et al., 2011; Gou et al., 2012; Santiago et al., 2013; Sun et al., 2013). In
409
the absence of BRs, BRI1 exists as an inhibitory homodimer via associating with the anchor
410
protein BKI1 in the cytoplasm (Wang and Chory, 2006). When BRs are present, the binding of
411
BRs to the extracellular LRR domain of BRI1 causes the release of BKI1 and consequent
412
heterodimerization between BRI1 and one of the SERK co-receptors (Santiago et al., 2013). The
413
activated BRI1 then activates BSK (BR signaling kinase), the BSU1 phosphatase, and the BIN2
414
GSK3-like kinase cascade (Mora-Garcia et al., 2004; Vert and Chory, 2006; Tang et al., 2008;
415
Kim et al., 2011; Sreeramulu et al., 2013). Finally, the BIN2-regulated BES1/BZR1/HAT1
416
transcription factors transmit BR signaling to the nucleus via activating the expression of BR
417
target genes (Yin et al., 2002; He et al., 2005; Vert and Chory, 2006; Oh et al., 2012; Zhang et al.,
418
2013).
419
Perception of the oligo-peptide ligand flg22 by FLAGELLIN SENSITIVE2 (FLS2) triggers
420
the association of FLS2 with its co-receptor BRI1-ASSOCIATED RECEPTOR KINASE1
421
(BAK1), followed by the phosphorylation of BOTRYTIS-INDUCED KINASE1 (BIK1)
422
(Gomez-Gomez and Boller, 2000; Chinchilla et al., 2007; Lu et al., 2010; Sun et al., 2013). The
423
phosphorylated BIK1 phosphorylates and activates the NADPH oxidase RBOHD, leading to the
424
production of reactive oxygen species (ROS) (Kadota et al., 2014). BAK1 also phosphorylates
425
two closely related E3 ubiquitin ligases, PUB12 and PUB13, which polyubiquitinate FLS2 and
426
promote flagellin-induced degradation of FLS2 (Lu et al., 2011). The ERECTA (ER/ERL)
427
family and TOO MANY MOUTHS (TMM) perceive the secreted cysteine-rich peptides
428
EPIDERMAL PATTERNING FACTOR (EPF)/EPF-LIKE to control the proper density and
429
spacing of stomata (Shpak et al., 2005; Sugano et al., 2010; Lee et al., 2015). ER/ERL signaling
430
activates the downstream YODA (YDA), MKK4/5 and MPK3/6 MAPK cascade, resulting in the
431
inhibition of the SPEECHLESS (SPCH)/SCREAM (SCRM) transcriptional module (Lampard et
432
al., 2008; Lee et al., 2012; Torii, 2015). Moreover, ER/ERL and MPK3/MPK6 function in anther
433
lobe formation and anther cell differentiation (Hord et al., 2008). In the presence of the
434
CLAVATA3 (CLV3) ligand, CLV1 acts together with CORYNE (CRN)/CLV2, possibly via a 15
435
MAPK cascade, to inhibit the expression of WUSCHEL (WUS), WUSCHEL-RELATED
436
HOMEOBOX (WOX), and HAIRY MERISTEM (HAM) (Clark et al., 1997; Fletcher et al., 1999;
437
Schoof et al., 2000; Muller et al., 2008; Ogawa et al., 2008; Zhou et al., 2015). In the current
438
study, we established βCAs as downstream effectors of LRR-RLK, shedding light on
439
LRR-RLK-linked signal transduction pathways. EMS1 directly targets an enzyme rather than the
440
MAPK cascade, which might be required for rapid responses to environmental and
441
developmental cues. Moreover, our results indicate that phosphorylation mediated by a
442
receptor-like kinase is essential for the functional regulation of βCA.
443
Although decades of research have implicated CAs in various metabolic, cellular and
444
physiological processes in both plants and animals, the role of CAs in plant cell signaling
445
remains a mystery. CAs regulate the levels of CO2, protons and bicarbonate (HCO3‾), which are
446
important for carboxylation/decarboxylation, pH regulation, inorganic carbon transport, ion
447
transport, as well as water and electrolyte balance. In animals, CAs, which are found in the
448
plasma membrane, cytosol and mitochondria, are essential for gluconeogenesis, lipogenesis,
449
ureagenesis, and tumor formation. For example, the membrane-localized CA IX, an established
450
marker for tumor hypoxia and aggressive cancers, plays a critical role in tumor cell proliferation
451
by regulating pH (Alterio et al., 2009; Swietach et al., 2009). Besides their photosynthetic
452
functions (Hatch and Burnell, 1990; Badger and Price, 1994; Smith and Griffiths, 2000; Studer et
453
al., 2014; Clement et al., 2016; Ho et al., 2016; Ludwig, 2016), recent studies found that βCAs
454
are critical for CO2-induced stomatal development and movement (Hu et al., 2010; Engineer et
455
al., 2014; Hu et al., 2015; Matrosova et al., 2015; Wang et al., 2016a) as well as for amino acid
456
biosynthesis (DiMario et al., 2016).
457
Arabidopsis has six βCAs, each with a variable number of splice variants (Ferreira et al.,
458
2008; Wang et al., 2014). Our studies identified the βCA1 splice variant βCA1.4. Protein
459
localization studies show that βCA1 and βCA5 are found in plastids, including chloroplasts,
460
while βCA2 and βCA3 are localized to the cytosol, βCA4 is found at the plasma membrane, and
461
βCA6 is localized in mitochondria (Fabre et al., 2007). Furthermore, βCA1 is found in the
462
vicinity of the plasma membrane and chloroplasts (Hu et al., 2010) or near the plasma membrane
463
and cytoplasm when the first 65 amino acids (chloroplast signal peptide) are removed (Hu et al.,
464
2015). We found that βCA1.3 was localized in chloroplasts and at the plasma membrane,
465
whereas βCA1.4 was localized at the plasma membrane and in the cytoplasm. The multiplicity of 16
466
βCA isoforms and their diverse localizations suggest that βCAs might possess additional
467
functions.
468
CA activity in animals is regulated by phosphorylation. Phosphorylation stimulated by
469
3’-5’-cyclic adenosine monophosphate (cAMP) increases the activity of CAs from rat gastric
470
tissue (Bersimbaev et al., 1975) and rat astroglial cell cultures (Church et al., 1980). When
471
phosphorylated by protein kinase A (PKA) (Narumi and Miyamoto, 1974) and protein kinase G
472
(PKG) (Carrie and Gilmour, 2016), the activity of CA is enhanced in bovine erythrocytes and
473
rainbow trout gill, respectively. Human CA IX is a tumor-associated transmembrane carbonic
474
anhydrase. Phosphorylation on Thr443 is required for the function of CA IX in hypoxic tumor
475
cells (Ditte et al., 2011). In the single-cell algae, Chlamydomonas reinhardtii, Cah3, an
476
intracellular α-CA, is localized in the thylakoid lumen and its activity is also regulated by
477
phosphorylation (Blanco-Rivero et al., 2012). In this study, we found that phosphorylation by the
478
receptor-like kinase increases the activity of CAs in flowering plants. We also identified four
479
phosphorylation sites (Thr35, Thr54, Thr69, and Ser189) in βCA1.
480
Both the phosphorylation-blocking mutation T35A and the phosphomimic mutation
481
T35D in βCA1.4 caused the loss of enzyme activity, even after EMS1 treatment. T54A, T69A or
482
S189A mutation did not significantly alter βCA1 activity, but the enhancement of activity by
483
phosphorylation was significantly affected by these mutations. In particular, the activity of
484
βCA1.4T189A remained unchanged without or with EMS1 treatment; however, the S189D
485
mutation resulted in a significant increase in βCA1.4 activity. In addition, EMS1 treatment
486
further enhanced the activity of βCA1.4S189D, suggesting that phosphorylation of S189 is critical
487
for the regulation of βCA1 activity. It would be worthwhile to investigate how phosphorylation
488
of these residues affects βCA activity in the future.
489
Our loss- and gain-of-function studies of βCAs showed that βCAs are required for
490
tapetal differentiation. Tapetal cells produce elaioplasts (tapetum-specific plastids) and
491
tapetosomes (an ER-derived organelle rich in triacylglycerols and oleosins) (Dickinson, 1973;
492
Wu et al., 1997; Hsieh and Huang, 2007). Plastids in tapetal cells are essential for pollen wall
493
and pollen coat formation (Owen and Makaroff, 1995; Pacini, 1997; Clement and Pacini, 2001).
494
Lipids are the main precursors for components of pollen exine, such as sporopollenin. In tapetal
495
cells, pro-plastids undergo division during early tapetum development and subsequently develop
496
into non-green plastids (elaioplast) that are involved in the biosynthesis of tapetal lipids as well 17
497
as starch accumulation and/or mobilization (Dickinson, 1973; Pacini et al., 1992; Weber, 1992;
498
Clement et al., 1998; Wu et al., 1999; Clement and Pacini, 2001). In Brassicaceae species
499
including Arabidopsis, fully differentiated tapetal cells accumulate elaioplasts and tapetosomes.
500
In the male sterile1 mutant, tapetal cells produce greatly reduced numbers of elaioplasts and
501
tapetosomes (Ito et al., 2007; Yang et al., 2007). Mutations in the Arabidopsis MS2 and rice
502
DEFECTIVE POLLEN WALL (DPW) genes, which encode plastid-localized fatty acid
503
reductases, result in abnormal tapetum and pollen development (Aarts et al., 1997; Shi et al.,
504
2011). Disruption of phosphoenolpyruvate/phosphate translocator1 (PPT1) and the
505
plastid-localized enolase1 (ENO1) affect sporopollenin formation (Prabhakar et al., 2010). We
506
found that elaioplast and tapetosomes production was reduced when the function of βCAs was
507
disrupted. In animals, the importance of CAs increases in pathological states. Hypoxia-induced
508
CA IX facilitates cancer cell survival and proliferation by combating the high rate of glycolytic
509
metabolism to keep up with the increased energy demand for ATP and biosynthetic precursors
510
(Parks et al., 2013). Like tumor cells, tapetal cells might require high βCA activity to maintain
511
their highly active metabolic state. HCO3‾ is important for lipid formation. Based on our results,
512
the phosphorylation of βCA1 by EMS1 significantly enhances its activity. The highly active
513
βCAs might be required for tapetum development via affecting the formation of elaioplasts and
514
tapetosomes.
515
It is also possible that βCAs regulate tapetal cell pH, which might be important for tapetal
516
cell differentiation and the maintenance of tapetal function. The regulation of extracellular (pHe)
517
and intracellular pH (pHi) is critical for cell division, differentiation, and survival. In animals,
518
CAs play a key role in buffering cellular pH via regulating HCO3‾ and H+ concentrations (Alterio
519
et al., 2009; Chiche et al., 2009; Swietach et al., 2009; Swietach et al., 2010; Parks et al., 2011;
520
Benej et al., 2014). In plants, in addition to H+ pumps, such as P-type H+-adenosine
521
triphosphatase, vacuolar H+-ATPase, and H+-pyrophosphatase (Li et al., 2005), EMS1-regulated
522
βCAs might be especially important for moderating pH in tapetal cells, because they are highly
523
active in metabolism. In fact, our data revealed that the pH of epidermal cells and tapetal cells
524
differed in wild-type anthers. Moreover, loss-of-function of βCAs caused a significant decrease
525
in tapetal cell pH. Auxin signaling is highly active in the tapetum (Aloni et al., 2006; Cecchetti et
526
al., 2017), suggesting that auxin might be important for tapetal cell differentiation. Auxin
527
represses chloroplast and amyloplast development (Miyazawa et al., 1999; Kobayashi et al., 18
528
2012). Therefore, auxin might regulate tapetal cell differentiation and function via affecting the
529
formation of elaioplasts (Sakata et al.; Miyazawa et al., 1999; Cecchetti et al., 2008; Kobayashi
530
et al., 2012). Furthermore, auxin is essential for pollen development and filament elongation
531
(Sakata et al.; Cecchetti et al., 2008). The H+ gradient maintained by βCAs might be important
532
for auxin transport during anther development. Normal tapetum development is required for
533
sexual reproduction and high yield in plants under both normal and stress conditions (Smith and
534
Zhao, 2016). Future research should be focused on investigating the molecular mechanisms by
535
which βCAs control tapetal cell differentiation and pollen development.
536 537 538
METHODS
539 540
Plant Materials and Growth Condition
541 542
The Arabidopsis thaliana Landsberg erecta (Ler) and Columbia (Col-0) ecotypes were used in
543
this study. Plants were grown in Metro-Mix 360 (Sun-Gro Horticulture) in growth chambers
544
(Philips PLUS T8 High Output 8-foot cool white fluorescent lamps and 100μmol m−2 s−1
545
photon density) under a 16-hr light/8-hr dark photoperiod at 22°C and 50% humidity.
546 547
Phylogenetic Analysis
548 549
Alignment of βCA1 to βCA6 protein sequences was carried out with MUSCLE, followed by
550
manual adjustment. Phylogenetic analysis was performed by PhyML using the Maximum
551
Likelihood method with default parameters (Dereeper et al., 2008). TreeDyn was used to display
552
the phylogenetic tree.
553 554 555
Yeast Two-Hybrid (Y2H) Screening
556 557
The ProQuest Two-Hybrid system with Gateway technology (Invitrogen) was employed to
558
identify EMS1 interacting proteins. Briefly, the EMS1 kinase domain (852-1192), which was 19
559
cloned into the pDEST 32 vector, was used as the bait. To enrich the potential EMS1-interacting
560
proteins, mRNA was prepared from young buds containing stage-5/6 anthers in the Landsberg
561
erecta (Ler) background. A cDNA library was then constructed using the SuperScript™ Plasmid
562
System with the Gateway® Technology. According to the manufacturer’s manual,
563
protein-protein interactions were assayed on the synthetic drop-out medium minus Leu, Trp and
564
His, as well as containing 25 mM 3-Amino-1,2,4-Triazole (3AT), using the yeast strain MaV203.
565 566 567
Generation of Constructs and Transgenic Plants
568 569
All DNAs and cDNAs were amplified using Phusion High-Fidelity DNA Polymerase (New
570
England Biolabs). To test interactions between EMS1 and βCAs by Y2H, βCA1.4, βCA2.2,
571
βCA3 and βCA4.3 were cloned into the pENTR/D-TOPO vector (Invitrogen, Cat#: K240020),
572
respectively, followed by introduction into the pDEST22 vector using Gateway LR recombinase
573
II enzyme mix (Invitrogen, Cat#: 11791100). pSAT vectors (Tzfira et al., 2005) were used for
574
the subcellular localization study and the BiFC assay in the protoplast system. To generate the
575
BiFC constructs, cEYFP (C-terminus of EYFP; pSAT1-cEYFP-C1-B) was fused to the
576
full-length EMS1 and BRI1, and the nEYFP (N-terminal of EYFP; pSAT4-nEYFP-N1) was
577
fused to βCA1.4, βCA2.2, βCA3, βCA4.1, βCA5.2, βCA6.1 and αCA1.1. To generate constructs
578
for the Förster Resonance Energy Transfer (FRET) assay, the βCA1.4 was fused to
579
pSAT6-EYFP-N1 to produce βCA1-EYFP. The EYFP in pSAT6-EYFP-N1 was replaced by
580
CFP to generate pSAT6-CFP-N1. Full-length EMS1 was then inserted into pSAT6-CFP-N1 to
581
generate EMS1-CFP. To investigate subcellular localization, βCA1, βCA1.3 and βCA1.4 were
582
cloned into pSAT6--EYFP-N1. βCA1.4T35A, βCA1.4T54A, βCA1.4T69A, βCA1.4S189A, βCA1.4T35D,
583
βCA1.4T54D, βCA1.4T69D and βCA1.4S189D were generated by overlapping PCR and cloned into
584
the pENTR/D-TOPO vector. βCA1.4T35A and βCA1.4S189A were cloned into pSAT6--EYFP-N1
585
for subcellular localization analysis. To test the dimerization of βCA1.4 and its mutated versions,
586
βCA1.4, βCA1.4T35A and βCA1.4S189A were fused to nEYFP and cEYFP, respectively.
587
For the co-immunoprecipitation assay using protoplasts, βCA1.4-Flag was
588
PCR-amplified and inserted into pSAT6-EYFP-N1 after removing the EYFP tag to generate
589
pSAT6- βCA1.4-Flag-N1. 20
590
For βCA protein localization studies in planta, genomic DNA fragments including
591
promoters, introns and exons were cloned into the pENTR/D-TOPO vector. The βCA1 genomic
592
fragment without the last intron and exon was used to construct ProβCA1:βCA1-GFP. Genomic
593
DNA fragments with all exons were used to construct ProβCA2:βCA2-GFP,
594
ProβCA3:βCA3-GFP and ProβCA4:βCA4-GFP, respectively.
595
For the complementation experiments, 4884-bp, 5865-bp and 5207-bp genomic DNA
596
fragments of βCA1, βCA2 and βCA4 were cloned into the pENTR/D-TOPO vector, respectively.
597
The A9 (At5g07230) gene promoter (Feng and Dickinson, 2010) was amplified and cloned into
598
the pENTR/D-TOPO vector. The βCA1.4 cDNA was inserted after the A9 promoter to produce
599
ProA9:βCA1.4. The same method was used to produce the ProA9:βCA1.4T35A,
600
ProA9:βCA1.4T54A, ProA9:βCA1.4T69A, ProA9:βCA1.4S189A, ProA9:βCA1.4T35D and
601
ProA9:βCA1.4S189D constructs. AimrβCA1-4 for generating artificial microRNAs targeting βCA1
602
to βCA4 transcripts (βCA1-4) was designed as described previously (Schwab et al., 2006). The
603
aimrβCA1-4 fragment was cloned into the pENTR/D-TOPO vector to generate
604
Pro35S:aimrβCA1-4 and ProA9:aimrβCA1- 4, respectively.
605
For produce the Pro4x35S-βCA1:βCA1 construct, the CaMV 35S enhancer was amplified
606
from the pSK1015 vector (Weigel et al., 2000) and cloned into the pENTR/D-TOPO vector. A
607
4884-bp fragment of genomic DNA, including a 2-kb βCA1 promoter region, was inserted after
608
the CaMV 35S enhancer.
609
For the phosphorylation assays, a 950-bp cDNA fragment of EMS1 was amplified and
610
cloned into the pGEX-4T-2 vector (GE Healthcare, Cat#: 27-1542-01) to generate EMS1-KD
611
(kinase domain)-GST. βCA1.4 was amplified and cloned into the pET28a vector (EMD Millipore,
612
Cat#: 69864-3) to produce the βCA1.4-His protein. βCA1.4, βCA2.2, βCA4.1, βCA1.4T35A,
613
βCA1.4T54A, βCA1.4T69A, βCA1.4S189A, βCA1.4T35D, βCA1.4T54D, βCA1.4T69D and βCA1.4S189D were
614
PCR-amplified and cloned into the pGEX-4T-2 vector to produce GST fusion proteins,
615
respectively.
616
All constructs generated with pENTR/D-TOPO vectors were introduced into Gateway
617
Binary vectors (pGBWs or pEARLEYs) using Gateway LR recombinase II enzyme mix
618
(Invitrogen). Detailed information for all constructs and primers is shown in Supplemental Data
619
Set 1 and 2. The resulting constructs were transformed into Agrobacterium strain GV3101. Plant
620
transformation was performed using the floral dip method (Clough and Bent, 1998). 21
621
Transformants were screened on 50 µg/mL of kanamycin and 25 µg/mL of hygromycin or 1%
622
Basta.
623 624 625
Protoplast Transfection and BiFC assay
626 627
For high-quality plasmid preparation, the PureYieldTM Plasmid Midiprep System (Promega, Cat#:
628
A2495) was used to isolate plasmids. A pair of constructs was used to co-transfect Arabidopsis
629
protoplasts prepared from 4-week-old leaves (Yoo et al., 2007). The Pro35S:EYFP construct was
630
used to monitor transfection efficiency. At least three replicates were performed for each assay.
631
Samples were observed after 16 hours under a Leica TCS SP2 laser scanning confocal
632
microscope using a 63×/1.4 water immersion objective lens (Huang et al., 2016c).
633 634 635
Co-immunoprecipitation
636 637
The Co-immunoprecipitation assay was performed as previously described (Jia et al., 2008; Li et
638
al., 2017). Young buds were harvested from wild-type, ProEMS1:EMS1-4xcMyc,
639
ProβCA1:βCA1-Flag and ProEMS1:EMS1-4xcMyc ProβCA1:βCA1-Flag plants. βCA1.4-Flag
640
was also transiently expressed in protoplasts from Pro35S:EMS1-cMyc and Pro35S:EMS1-cMyc
641
Pro35S:TPD1sp-GFP-∆TPD1 leaves. Non-transfected wild-type and Pro35S:EMS1-cMyc
642
protoplasts were used as controls. Membrane proteins were extracted and incubated with 30 µl of
643
Myc-Trap® beads (Chromotek, Cat# yta-10) overnight at 4°C. After three washes, the eluted
644
immunoprecipitates were used for protein gel blot analysis to detect EMS1-4xcMyc (anti-cMyc,
645
Roche; Cat. #: 11667149001; 1:500 dilution), TPD1sp-GFP-∆TPD1 (anti-GFP, Torrey Pines
646
Biolabs) and βCA1-Flag (anti-Flag, GeneScript; Cat#: A01428-100; 1:1000 dilution),
647
respectively.
648 649 650
FRET Assay
651
22
652
For the FRET assay, EMS1-CFP (donor) and βCA1.4-EYFP (acceptor) were co-transfected into
653
Arabidopsis protoplasts. Control reactions included single-plasmid transfections and mock
654
transfections (without plasmid). Acquisition of spectrally resolved fluorescence images, spectral
655
unmixing of donor, acceptor, and autofluorescence signals, as well as determination of FRET
656
efficiency were performed as previously described (Raicu and Singh, 2013; King et al., 2016;
657
Stoneman et al., 2016). Briefly, an optical micro-spectroscope (OptiMiS) was used to capture
658
images with linear unmixing of CFP and EYFP from cells co-expressing EMS1-CFP and
659
βCA1-EYFP. OptiMiS is featured with the line-scan excitation powered by a femtosecond laser
660
emitting near-IR light pulses with tunable emission wavelength between 780 nm and 1040 nm
661
(Raicu et al., 2009; Biener et al., 2014). The spectral images were used to generate
662
two-dimensional fluorescence maps of donor and acceptor (Raicu and Singh, 2013). The CFP,
663
EYFP and auto-fluorescence spectra were acquired from cells expressing EMS1-CFP only,
664
βCA1-EYFP only, and that were devoid of any fluorescent protein, respectively. The unmixing
665
protocol has been optimized for cells presenting substantial autofluorescence emission with the
666
final goal to determine donor and acceptor emission images free from any autofluorescence
667
(Mannan et al., 2013). Since the laser wavelength of 860 nm (selected as optimal for donor
668
excitation) also caused slight excitation of the acceptor, we did not attempt to compute the FRET
669
efficiency for every image pixel from the donor and acceptor fluorescence maps (Raicu and
670
Singh, 2013), or else the FRET efficiency map would have been contaminated by such spurious
671
excitation. Instead, we used a method introduced recently, which employed two excitation
672
wavelengths (860 nm and 960 nm) to account for the direct excitation of acceptor and to
673
compute the average FRET efficiency for regions of interest (at the plasma membrane) (King et
674
al., 2016; Stoneman et al., 2016). In order to reduce noise and avoid singularities at pixels with
675
zero intensity values, a signal-to-noise threshold of 1.5 SD was applied to FRET efficiency
676
calculations.
677 678 679
Pollen Staining and Semi-thin Sectioning of Anthers
680 681
To identify viable pollen grains, mature anthers prior to dehiscence were collected for Alexander
682
staining of pollen (Zhao et al., 2002; Xin et al., 2017). To study anther cell differentiation, 23
683
semi-thin sectioning was carried out essentially as previously described (Zhao et al., 2002;
684
Huang et al., 2016b). For semi-thin sectioning, dissected young buds were fixed in 2.5% of
685
(vol/vol) glutaraldehyde in 0.1 M HEPES (N-2-Hydroxyethyl pipera-zine-N-2-ethanesulfonic
686
acid) buffer (pH 7.2) and 0.02% of TritonX-100 overnight at 4°C. Samples were washed twice in
687
0.1 M HEPES buffer (pH 7.2) and post-fixed in 1% of OsO4 overnight. The samples were
688
dehydrated in a graded acetone series (10% of increments), infiltrated, and embedded in Spurr’s
689
resin. Semi-thin (1 µm) sections were made using an Ultracut E ultramicrotome (Reichert-Jung)
690
and stained with 0.25% of Toluidine Blue O. The number of tapetal cells was counted under a
691
microscope from central sections of upper two lobes. Sections from 20 individual buds were
692
examined.
693 694 695
Phosphorylation Assay and Identification of βCA1 Phosphorylation Sites
696 697
EMS1-KD-GST, βCA2.2-GST and βCA4.1-GST were expressed in BL21DE3 according to the
698
manufacturer’s instructions (GE Healthcare) and protein extraction was performed using the
699
Pierce™ Glutathione Agarose (Thermo Scientific™, Cat#: 16100). . The expression and
700
purification of βCA1.4-His protein were performed as previously described (Idrees et al., 2016).
701
βCA1.4-His protein was extracted by the HIS-Select® Nickel Affinity Gel (Sigma, Cat#: P6611)
702
and further purified by Hi Trap DEAE FF column (GE healthcare, Cat#: 17-5055-01). For the
703
phosphorylation reaction, 1 µg of EMS1-KD-GST and 5 µg of βCA1.4-His were incubated with
704
[γ-32P] ATP in the kinase buffer (HEPES at pH7.4, 10 mM MgCl2, 10 mM MnCl2, 1 mM DTT)
705
for 30 min at 25 oC (Zhao et al., 2002; Li et al., 2017). Proteins were separated by 10% of
706
SDS-PAGE and the gel was then analyzed by autoradiography. Additionally, 1 µg of
707
EMS1-KD-GST and 5 µg of βCA2.2-GST and βCA4.1-GST were incubated in kinase buffer for
708
30 min at 25 oC. Proteins were separated by 10% SDS-PAGE and protein phosphorylation was
709
detected with anti-phospho-(Ser/Thr) antibody (Abcam, Cat#: ab17464, 1:400 dilution).
710
For mass spectrometer analysis, the purified recombinant βCA1.4-His, βCA2.2-GST and
711
βCA4.1-GST proteins were trans-phosphorylated by EMS1-KD-GST in vitro in the presence of
712
ATP. After SDS-PAGE, protein bands were excised from the gels, followed by in-gel digestion
713
using trypsin. Samples were analyzed with an Agilent 1100 series LC/MSD Trap SL coupled to 24
714
the LTQ Orbitrap XL (Thermo Scientific) mass spectrometer (University of Wisconsin-Madison
715
Biotechnology Center).
716 717 718
Carbonic Anhydrase Activity Assay
719 720
Phosphorylation of βCA1.4, βCA1.4T35A, βCA1.4T54A, βCA1.4T69A, βCA1.4S189A, βCA1.4T35D,
721
βCA1.4T54D, βCA1.4T69D and βCA1.4S189D was performed as described above. Four-week-old
722
leaves and young buds from wild-type and Pro35S:TPD1 Pro35S:EMS1 transgenic plants
723
(Huang et al., 2016d) were harvested. The CA activity assay was performed as previously
724
described (Wilbur and Anderson, 1948; Hu et al., 2010). CA activity was defined as WA unit/mg
725
protein. Protein concentration was determined by the Bradford method.
726 727 728
RT-PCR, qRT-PCR, and RNA in situ Hybridization
729 730
Total RNA was isolated from various plant tissues/organs using an RNeasy Plant Mini Kit
731
(Qiagen, Cat#: 74904). RNA quantification, reverse transcription, PCR, qRT-PCR (DNA Engine
732
Opticon 2 system), and data analysis were performed as described previously (Liu et al., 2009;
733
Huang et al., 2016a). Expression of βCA1.1, βCA1.2, βCA1.3 and βCA1.4 was determined by
734
amplifying each full-length CDS. RNA in situ hybridization was performed using anthers from
735
wild-type, βca1 βca2 βca4, Pro35S:amirβCA1-4, ProA9:amirβCA1-4, ProA9:βCA1.4/βca1 βca2
736
βca4 and Pro4x35S-βCA1:βCA1 plants (Zhao et al., 2002; Liu et al., 2010). An SP6/T7 DIG
737
RNA Labeling Kit (Roche, Cat#: 11175025910) was used to generate A9 sense and antisense
738
probes. Primers for PCR, qRT-PCR and in situ hybridization are listed in Supplemental Data Set
739
2.
740 741 742
Microscopy
743 744
Images of pollen staining and semi-thin sections were photographed under an Olympus BX51 25
745
microscope equipped with an Olympus DP 70 digital camera (Jia et al., 2008; Huang et al.,
746
2016d). For confocal microscopy analysis, samples were observed under a Leica TCS SP2 laser
747
scanning confocal microscope. Samples were mounted in water and observed with a 20× lens or
748
a 63×/1.4 water immersion objective lens. Protoplasts were observed using a 63×/1.4 water
749
immersion objective lens. A 488-nm laser was used to excite GFP, EYFP, and chlorophyll. The
750
emission was captured using PMTs set at 505–530 nm, 500–550 nm, and 644–719 nm,
751
respectively.
752 753 754
pH Measurement
755 756
Anthers were dissected from flower buds and stained with 20 µm of SNARF-1-AM (Invitrogen,
757
Cat#: c1271) in MES/KCl buffer (5mM KCl, 50 µM CaCl2, 10mM MES buffered to pH 6.15
758
with KOH) for 30 minutes (Zhang et al., 2001). The anthers were washed three times with
759
SNARF-1-AM-free buffer. Confocal imaging was performed on a Leica SP2 confocal laser
760
scanning microscope (Leica Microsystems). The excitation was set at 488 nm. The emission was
761
set at 540-590 nm for channel 1 and 610-700 nm for channel 2 (Sano et al., 2009). Images were
762
analyzed by ImageJ. The intensity ratio of channel 1/channel 2 was converted to pH according to
763
the calibration graph (Zhang et al., 2001; Leshem et al., 2006).
764 765
Accession Numbers
766 767
Sequence data from this article can be found in The Arabidopsis Information Resource
768
(http://www.arabidopsis.org/) under the following accession numbers: EMS1 (At5G07280),
769
TPD1 (AT4G24972), BRI1 (AT4G39400), βCA1 (At3g01500), βCA2 (AT5G14740), βCA3
770
(AT1G23730), βCA4 (AT1G70410), βCA5 (AT4G33580), βCA6 (AT1G58180), αCA1
771
(AT3G52720), A9 (At5G07230), and ATA7 (At4g28395).
772 773 774
Supplemental Data
775
26
776
Supplemental Figure 1. Alternative splice variants of βCA1.
777
Supplemental Figure 2. EMS1 interacts with four βCA1 variants, possibly at different levels.
778
Supplemental Figure 3. Phylogenetic analysis of βCAs in Arabidopsis.
779
Supplemental Figure 4. Expression of EMS1-cEYFP and βCAs-nEYFP in a BiFC experiment.
780
Supplemental Figure 5. TPD1 is not required for EMS1 and βCA1 interaction.
781
Supplemental Figure 6. Expression of four βCA1 splice variants in different tissues.
782
Supplemental Figure 7. Identification of the βca2 mutant.
783
Supplemental Figure 8. βCA2 and βCA4 complement the βca1 βca2 βca4 mutant phenotype.
784
Supplemental Figure 9. Overexpression of βCA1 affects male fertility.
785
Supplemental Figure 10. βCA2.2 and βCA4.1 are phosphorylated by EMS1 in vitro.
786
Supplemental Figure 11. Phosphorylation sites identified in βCA1.4.
787
Supplemental Figure 12. Expression of wild-type and mutated βCA1.4 transgenes in
788
complementation experiments.
789
Supplemental Figure 13. Expression of βCA1.4 and βCA1.4S189D in transgenic plants.
790
Supplemental Table 1. Calculation of FRET efficiency using a two-excitation wavelength
791
scheme.
792
Supplemental Table 2. Phosphorylation sites in βCA1.4 identified by LC/MS/MS.
793
Supplemental Data Set 1. Constructs generated in this study.
794
Supplemental Data Set 2. Primers used in this study.
795
Supplemental File 1. Alignment used to produce the phylogenetic tree shown in Supplemental
796
Figure 3.
797 798
ACKNOWLEDGMENTS
799 800
We thank T. Schuck, J. Gonnering and P. Engevold for plant care, P. He, and J. Sheen for the
801
protoplast transfection, H. Owen for confocal microscopy, the ABRC for EMS1, TPD1, A9,
802
βCA1, βCA2, βCA3, βCA4, βCA5, βCA6, BRI1, αCA1 BAC clones and cDNAs, T. Nakagawa for
803
the pGBW vectors, C. Pikaard for the pEARLEY vectors, J. Haseloff for the pBIN Gal4–
804
mGFP5er vector, D. Weigel for the pRS300 and pSK1015 vectors, S. Gelvin for the pSAT
805
vectors, and S. Ishiguro for the FIB1a-GFP and GRP17-GFP vectors. This work was supported
806
by the National Science Foundation grants IOS-0721192 and IOS-1322796 to D.Z., the Research 27
807
Growth Initiative (RGI) at the University of Wisconsin-Milwaukee to D.Z. and V.R., and the
808
UW-Madison/UW-Milwaukee Intercampus Research Incentive Grants Program to D.Z. D.Z.
809
also gratefully acknowledges the support of the Shaw Scientist Award from the Greater
810
Milwaukee Foundation, the Bradley Catalyst Award from the UWM Research Foundation, and
811
the CAS/SAFEA International Partnership Program for Creative Research Teams. The optical
812
microspectroscopy imaging facility used for this study was developed with financial support
813
from the National Science Foundation, Major Research Instrumentation Program (Grant:
814
PHY-1126386 awarded to V.R.). The FRET assays were partially funded by the National
815
Science Foundation, Physics of Living Systems Program (Grant: PHY-1058470 awarded to
816
V.R.).
817 818 819
Competing Finacial Interests
820 821
Dr. Raicu is a co-founder of Aurora Spectral Technologies LLC (AST), which commercializes
822
the OptiMiS technology used in this work. All the other authors declare no competing financial
823
interests.
824 825 826
AUTHOR CONTRIBUTIONS
827 828
J.H. and D.Z. designed the research. J.H., Z.L., G.B., E.X., S.M., N.E., C.Z., H.K. and V.R.
829
performed research. All authors analyzed data. J.H. and D.Z. wrote the article with contributions
830
from the other authors.
831 832 833 834
FIGURE LEGENDS
835 836 837
Figure 1. EMS1 interacts with βCAs. 28
838 839
(A) Yeast two-hybrid (Y2H) assay using the EMS1 kinase domain (EMS1-KD) fused to BD
840
(DNA binding domain, BD-EMS1-KD) and six AD (DNA activation domain)-βCA fusions.
841
Yeast cells were grown on synthetic drop-out medium lacking Leu, Trp, and His with 25 mM of
842
3-Amino-1,2,4-Triazole (3AT).
843
(B) to (I) Bimolecular fluorescence complementation (BiFC) assay using Arabidopsis mesophyll
844
protoplasts. Confocal images showing that EMS1 interacts with βCA1.4 (B), βCA2.2 (C), and
845
βCA4.1 (E) at the plasma membrane; EMS1 does not interact with βCA3 (D), βCA5.2 (F),
846
βCA6.1 (G), or αCA1.1 (H, control); and no interaction between BRI1 and βCA1.4 (I, control).
847
(J) to (L) Förster Resonance Energy Transfer (FRET) assay showing fluorescence images of the
848
donor (EMS1-CFP) obtained only following excitation at 860 nm (J) and the acceptor
849
(βCA1.4-EYFP) only following excitation at 860 nm (K) and 960 nm (L). Signals were obtained
850
by spectral unmixing of images in cells co-expressing EMS1-CFP and βCA1.4-EYFP.
851
(M) Membrane proteins were extracted from young buds of the wild type (Lane 1),
852
ProEMS1:EMS1-4xcMyc (Lane 2), ProβCA1:βCA1-Flag (Lane 3) and ProEMS1:EMS1-4xcMyc
853
ProβCA1:βCA1-Flag transgenic plants (Lane 4). After immunoprecipitation by anti-cMyc
854
antibody, βCA1-Flag was detected from proteins of ProEMS1:EMS1-4xcMyc
855
ProβCA1:βCA1-Flag plants.
856 857
Figure 2. Expression analyses of βCAs in anthers.
858 859
(A) RT-PCR showing the expression of four splice variants of βCA1 in wild-type (WT) young
860
buds as well as in WT and ems1 anthers.
861
(B) RT-PCR showing the expression of βCA2, βCA3, and βCA4 in WT and ems1 anthers. PCR
862
products represent total transcripts of βCA2, βCA3, and βCA4. The ACTIN2 gene was used as an
863
internal control.
864
(C) to (H) Confocal images showing the localization of βCA1 protein in ProβCA1:βCA1-GFP
865
anthers. Green: GFP signal, red: autofluorescence from chloroplasts, S: anther stage. βCA1 was
866
detected at low levels in the epidermis at stage 4 (C) and at high levels in tapetal cells at stage 5
867
(D) and stage 6 (E), which are two key stages for tapetal cell differentiation. βCA1 levels in
868
tapetal cells gradually decreased at stage 7 (F) and stage 8 (G). No βCA1 was observed in 29
869
stage-10 anthers (H). Bars = 50 µm.
870
(I) to (K) βCA1 was detected at the plasma membrane and in the cytoplasm of tapetal cells in
871
stage-5 anthers (I). (J) FM4-64 stained image of (I). (K) Merged image of (I) and (J). Insets
872
showing tapetal cell at high magnification (Arrowheads: plasma membrane). Bars = 10 µm.
873
(L) to (N) Confocal images showing relatively weak GFP signals in both tapetal cells and
874
microsporocytes in ProβCA2:βCA2-GFP (L) and ProβCA4:βCA4-GFP (N) anthers at stage 5,
875
but no GFP signals in ProβCA3:βCA3-GFP anthers (M). E: epidermis, M: microsporocyte, and
876
T: tapetal cell. Bars = 50 µm.
877 878 879
Figure 3. Downregulation of βCAs impairs male fertility.
880 881
(A) to (G) Wild type (A), small and sterile βca1 βca2 βca4 (B), Pro35S:amirβCA1-4 (C) and
882
ProA9:βCA1.3/βca1 βca2 βca4 (G) plants; normal growth but sterile ProA9:amirβCA1-4 plants
883
(D), normal growth and fertile ProβCA1:βCA1/βca1 βca2 βca4 plants (E), and fertile but small
884
ProA9:βCA1.4/βca1 βca2 βca4 plants (F) (All six weeks old). Arrows: fertile siliques,
885
Arrowheads: sterile siliques; Bars = 2 cm.
886
(H) to (N) Alexander staining of pollen in mature anthers showing viable pollen grains in
887
wild-type (H), ProβCA1:βCA1/βca1 βca2 βca4 (L) and ProA9:βCA1.4/βca1 βca2 βca4 (M)
888
anthers, but no viable pollen grains in βca1 βca2 βca4 (I), Pro35S:amirβCA1-4 (J),
889
ProA9:amirβCA1-4 (K) and ProA9:βCA1.3/βca1 βca2 βca4 (N) anthers. Bars = 50 µm.
890 891 892
Figure 4. Downregulation of βCAs causes abnormal tapetal cell differentiation.
893 894
(A) to (E) Semi-thin sections of stage-6 anther lobes showing four somatic cell layers, including
895
epidermis (E), endothecium (En), the middle layer (ML) and tapetum (T), as well as reproductive
896
microsporocytes (M) in the centers of wild-type (A) and ProA9:βCA1.4/βca1 βca2 βca4 (E)
897
anthers, but vacuolated tapetal-like (TL) cells in βca1 βca2 βca4 (B), Pro35S:amirβCA1-4 (C),
898
and ProA9:amirβCA1-4 (D) anthers.
899
(F) to (J) Semi-thin sections of stage-7 anther lobes showing normal tapetum and tetrads (Tds) in 30
900
wild-type (F) and ProA9:βCA1.4/βca1 βca2 βca4 (J) anthers, but highly vacuolated tapetal-like
901
cells and degenerating microsporocytes (DM, no tetrads formed) in βca1 βca2 βca4 (G),
902
Pro35S:amirβCA1-4 (H), and ProA9:amirβCA1-4 (I) anthers.
903
(K) to (O) Semi-thin sections of stage-9 anther lobes showing a thin tapetum and normal
904
microspores (Ms) in wild-type (K) and ProA9:βCA1.4/βca1 βca2 βca4 (O) anthers, but
905
degenerated tapetal-like cells and microsporocytes (anther lobes are empty) in βca1 βca2 βca4
906
(L), Pro35S:amirβCA1-4 (M), and ProA9:amirβCA1-4 (N) anthers.
907
(P) to (T) At stage 6, in situ hybridization results showing that the tapetal cell marker gene A9
908
was expressed strongly in a monolayer of tapetal cells in wild-type (P) and ProA9:βCA1.4/βca1
909
βca2 βca4 (T) anthers, but weakly expressed in vacuolated tapetal-like cells of βca1 βca2 βca4
910
(Q), Pro35S:amirβCA1-4 (R) and ProA9:amirβCA1-4 (S) anthers.
911
(U) and (V) Confocal images showing that the expression of the elaioplast marker gene
912
FIB1a-GFP in tapetal cells of wild-type anthers (U) was higher than that in tapetal-like cells of
913
ProA9:amirβCA1-4 anthers (V) at stage 8. Arrows indicate GFP signals in tapetal cells. Bars =
914
10 µm.
915
(W) and (X) Confocal images showing that expression of the tapetosome marker gene
916
GRP17-GFP was higher in tapetal cells of wild-type anthers (W) than in tapetal-like cells of
917
ProA9:amirβCA1-4 anthers (X) at stage 10. Arrows indicate GFP signals in tapetal cells. Bars =
918
10 µm.
919
(Y) qRT-PCR results showing the expression of A9 in wild-type, βca1 βca2 βca4,
920
Pro35S:amirβCA1-4, and ProA9:amirβCA1-4 anthers. Numbers indicate independent lines of
921
Pro35S:amirβCA1-4 and ProA9:amirβCA1-4. Three independent transgenic lines were used for
922
each transgenic plant and three samples were collected from each transgenic line. Bars indicate
923
SD. Asterisks indicate significant difference (P