# CHASE-Containing Histidine Kinase Receptors in Apple Tree: From a Common Receptor Structure to Divergent Cytokinin Binding Properties and Specific Functions

Dimitri Daudu, Elsa Allion, Franziska Liesecke, Nicolas Papon, Vincent Courdavault, Thomas Dugé de Bernonville, Céline Mélin, Audrey Oudin,

Marc Clastre, Arnaud Lanoue, et al.

► To cite this version:

Dimitri Daudu, Elsa Allion, Franziska Liesecke, Nicolas Papon, Vincent Courdavault, et al. CHASE-Containing Histidine Kinase Receptors in Apple Tree: From a Common Receptor Structure to Divergent Cytokinin Binding Properties and Specific Functions. Frontiers in Plant Science, 2017, 8, pp.First online. 10.3389/fpls.2017.01614. hal-01608088

HAL Id: hal-01608088

https://hal.science/hal-01608088

Submitted on 25 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
<figure><img src="image_2.png" /></figure>

Check for
updates

# CHASE-Containing Histidine Kinase Receptors in Apple Tree: From a Common Receptor Structure to Divergent Cytokinin Binding Properties and Specific Functions

Dimitri Daudu¹, Elsa Allion¹, Franziska Liesecke¹, Nicolas Papon², Vincent Courdavault¹, Thomas Dugé de Bernonville¹, Céline Mélin¹, Audrey Oudin¹, Marc Clastre¹, Arnaud Lanoue¹, Martine Courtois¹, Olivier Pichon¹, David Giron³, Sabine Carpin⁴, Nathalie Giglioli-Guivarc'h¹, Joël Crèche¹, Sébastien Besseau¹ and Gaëlle Glévarec¹*

¹ EA 2106 Biomolécules et Biotechnologies Végétales, Université François-Rabelais, Tours, France, ² EA 3142 Groupe d'Etude des Interactions Hôte-Pathogène, Université Angers, Angers, France, ³ UMR 7261 Institut de Recherche sur la Biologie de l'Insecte, Centre National de la Recherche Scientifique (CNRS), Université François-Rabelais, Tours, France, ⁴ EA 1207 Laboratoire de Biologie des Ligneux et des Grandes Cultures, Université d'Orléans, Orléans, France

Cytokinin signaling is a key regulatory pathway of many aspects in plant development and environmental stresses. Herein, we initiated the identification and functional characterization of the five CHASE-containing histidine kinases (CHK) in the economically important *Malus domestica* species. These cytokinin receptors named MdCHK2, MdCHK3a/MdCHK3b, and MdCHK4a/MdCHK4b by homology with *Arabidopsis* AHK clearly displayed three distinct profiles. The three groups exhibited architectural variations, especially in the N-terminal part including the cytokinin sensing domain. Using a yeast complementation assay, we showed that MdCHK2 perceives a broad spectrum of cytokinins with a substantial sensitivity whereas both MdCHK4 homologs exhibit a narrow spectrum. Both MdCHK3 homologs perceived some cytokinins but surprisingly they exhibited a basal constitutive activity. Interaction studies revealed that MdCHK2, MdCHK4a, and MdCHK4b homodimerized whereas MdCHK3a and MdCHK3b did not. Finally, qPCR analysis and bioinformatics approach pointed out contrasted expression patterns among the three MdCHK groups as well as distinct sets of co-expressed genes. Our study characterized for the first time the five cytokinin receptors in apple tree and provided a framework for their further functional studies.

**Keywords:** CHASE-containing histidine kinase, cytokinin, yeast complementation assay, RNAseq data, protein-protein interaction, *Malus domestica*

## INTRODUCTION

Cytokinins are essential adenine-derived plant hormones, gathering more than 40 structures substituted at the N⁶-position by an isoprenoid or aromatic chain (Spíchal, 2012; Osugi and Sakakibara, 2015). They are involved in numerous physiological processes such as cell division, delayed senescence, vascular tissue development, root architecture and light responses (Sakakibara, 2006; Kieber and Schaller, 2014; Zürcher and Müller, 2016). Cytokinins also play roles in the

## OPEN ACCESS

**Edited by:**
Ján A. Miernyk,
Agricultural Research Service (USDA),
United States

**Reviewed by:**
Jan Hejatko,
Masaryk University, Czechia
Dong Xu,
University of Missouri, United States

*Correspondence:
Gaëlle Glévarec
gaëlle.glevarec@univ-tours.fr

**Specialty section:**
This article was submitted to
Plant Physiology,
a section of the journal
Frontiers in Plant Science

Received: 16 June 2017
Accepted: 04 September 2017
Published: 20 September 2017

**Citation:**
Daudu D, Allion E, Liesecke F, Papon N, Courdavault V, Dugé de Bernonville T, Mélin C, Oudin A, Clastre M, Lanoue A, Courtois M, Pichon O, Giron D, Carpin S, Giglioli-Guivarc'h N, Crèche J, Besseau S and Glévarec G (2017) CHASE-Containing Histidine Kinase Receptors in Apple Tree: From a Common Receptor Structure to Divergent Cytokinin Binding Properties and Specific Functions. *Front. Plant Sci.* 8:1614. doi: 10.3389/fpls.2017.01614
interaction with both biotic and abiotic factors (Frugier et al., 2008; Giron and Glévarec, 2014; Naseem et al., 2014; Zwack and Rashotte, 2015).

Plant cytokinin perception is mediated by CHASE domain-
containing histidine kinase receptors (CHK) as first actors
of cytokinin signaling (Inoue et al., 2001). These receptors
display a complex multidomain structure with a N-terminal part
including at least two hydrophobic membrane-spanning domains
(TM) that border an extracytosolic sensing domain referred
to as CHASE (Cyclase/Histidine kinase Associated Sensory
Extracellular) (Anantharaman and Aravind, 2001; Mougel and
Zhulin, 2001) as well as a cytoplasmic C-terminal part containing
a catalytic histidine kinase domain (HK) and both receiver
and pseudo-receiver domains (REC and REC-like, respectively)
(Ueguchi et al., 2001). The HK domain is composed of an
HK dimerization and phosphoacceptor domain (HisKA) and
an HK catalytic domain called the HK-like ATPase domain
(HATPase). The cytokinin perception by the CHASE domain
leads to the autophosphorylation of a conserved histidine within
the HK domain. The phosphate residue is then transferred
to the REC domain on a conserved aspartate residue (Inoue
et al., 2001). Although the pseudo-receiver domain of CHK is
structurally similar to the REC domain, its functionality has
not been yet elucidated (Ueguchi et al., 2001; Lomin et al.,
2012). Subsequently, the signal is transferred by phosphorelay
to Response Regulators (RR) through histidine-containing
phosphotransfer shuttle proteins (HPt). While type-B RR5 (RRB)
are transcription factors that play a positive role in mediating
cytokinin-regulated gene expression, type-A RR5 (RRA) act as
negative regulators of cytokinin responses (To et al., 2004; Mason
et al., 2005; Ginis et al., 2012). In addition, Cytokinin Response
Factors (CRF) interact directly with HPts and were reported to
influence a subset of cytokinin responses (Cutcliffe et al., 2011;
Raines et al., 2016).

Cytokinin receptors were shown to localize mainly to the
endoplasmic reticulum (ER) both in *Arabidopsis thaliana* and
*Zea mays* (Caesar et al., 2011; Lomin et al., 2011; Wulfetange
et al., 2011). They are supposed to interact with each other,
forming potential homo- and hetero-dimers probably enabling
the *trans*-phosphorylation of the HK domain following cytokinin
perception (Dortay et al., 2006; Hothorn et al., 2011). However,
the signal transmission process across the membrane remains
unknown. The ligand-binding properties of the cytokinin
receptors have been investigated mostly using heterologous
assay systems through their expression in *Escherichia coli* or
*Saccharomyces cerevisiae* cells (Inoue et al., 2001; Romanov
et al., 2006; Stolz et al., 2011; von Schwartzenberg et al., 2016).
More recently, a plant assay system has been developed to
overcome the problem of alien membrane environment and the
difficulty to express some membrane receptors in bacteria or yeast
(Lomin et al., 2015). Overall, the cytokinin receptors differ in
their preference toward cytokinin forms (Yonekura-Sakakibara
et al., 2004; Lomin et al., 2011; Kuderová et al., 2015) but their
functional and specific properties as well as the structural changes
caused by cytokinin binding remain to be elucidated.

While it is established that the CHK receptors operate mostly
in a redundant fashion, the extensive studies of *Arabidopsis*

mutants have attributed some specific roles to single receptors.
Among others, AHK4 is the main regulator of primary root
growth and vascular morphogenesis whereas AHK2 and AHK3
are commonly involved in chlorophyll retention during leaf
senescence (Kim et al., 2006; Riefler et al., 2006). CHKs are
also involved in response to environmental changes (Zwack and
Rashotte, 2015). The three AHKs are also known to function
as negative regulators in osmotic stress responses (Tran et al.,
2007; Kumar and Verslues, 2015). AHK2 and AHK3 play an
additional negative regulatory role in cold stress (Jeon et al.,
2016) and ensure a protective function during light stress
(Cortleven et al., 2014). Moreover, cytokinin receptors also take
part in a large range of responses to biotic interactions. In
legume plants, cytokinin receptors regulate nodule formation
(Tirichine et al., 2007; Held et al., 2014; Boivin et al., 2016). In
*Arabidopsis*, the success of the pathogens *Rhodococcus fascians*
and *Hyaloperonospora arabidopsidis* depends on some AHKs
(Petry et al., 2009; Argueso et al., 2012). Finally, NaCHK2
and NaCHK3 modulate herbivory-induced defense signaling and
defenses in *Nicotiana attenuata* (Schäfer et al., 2015). If the
knowledge on the cytokinin receptors is increasingly important,
their study in various plant models is necessary for a complete
understanding of their biological functions.

Previous works on *M. domestica* reported a large accumulation
of cytokinins in the leaves infected by the insect *Phyllonorycter
blancardella*. This increase is responsible for the preservation of
nutrient green tissues when leaves are otherwise turning yellow
(Giron et al., 2007; Kaiser et al., 2010; Zhang et al., 2016). Based
on the involvement of cytokinins in this plant-biotic interaction,
we initiated the study of cytokinin signaling in apple tree with
a special focus on cytokinin receptors. Indeed, considering that
apple tree is one of the most cultivated fruit-tree with a continual
worldwide production increase, a greater knowledge of cytokinin
signaling pathway in this species could provide new opportunities
for agronomical and economical purposes. This study discloses
an overall and complete characterization of the five *M. domestica*
CHASE Histidine Kinases (MdCHKs).

MATERIALS AND METHODS

In Silico Sequence Analysis and
Receptor Identification

To identify MdCHK receptors, BLAST searches were performed
against the Genome Database for Rosaceae (GDR; Jung et al.,
2014) using *A. thaliana* cytokinin receptor sequences as queries
(AHK2, AHK3, and AHK4). Five sequences were identified based
on genome analysis and corresponding cDNA were amplified
from various plant organs using specific primers (Supplementary
Table S1). Sequences were registered in Genbank as *MdCHK2*
(KM114879), *MdCHK3a* (KM114880), *MdCHK3b* (KM114881),
*MdCHK4a* (KM114883) and *MdCHK4b* (KM114882).

*MdCHKs* gene organization has been visualized using the
FancyGene program (Rambaldi and Ciccarelli, 2009). Phylogeny
analyses were performed on conserved domains and local
similarities among proteins sequences. To this aim, multiple
protein sequence alignments were done using the COBALT tool
(Papadopoulos and Agarwala, 2007) and sequences were curated with Gblocks prior the construction of a bootstrap neighbor joining tree. Protein domain predictions were acquired using the SMART (Letunic et al., 2015) and PROSITE (Sigrist et al., 2002) programs, and transmembrane regions were identified with TMHMM (Krogh et al., 2001) and TMPred tools (Hofmann and Stoffel, 1993). Visualization of the transmembrane helices has been performed with a helical wheel drawing program<fnref n="1" />.

## Yeast Complementation Assay

The *S. cerevisiae* strain YIL147C, deficient in SLN1 receptor (*MAT*a/α, *ura3*, *leu2*, *his3*, *can1*Δ::*LEU2-MFA1*pro-HIS3/CAN1, *sln1*Δ::*KanMX/SLN1*) was used in complementation assays. Full-length coding sequences of MdCHKs were amplified and cloned into the yeast expression vector pYES2 under the control of the *GAL1*gene promoter using the NotI restriction site (for primers, see Supplementary Table S1). The *S. Cerevisiae* strain was transformed as follows. Cells were grown in 200 mL YPD liquid medium (150 rpm, 28°C) to 0.4–0.6 OD₆₀₀, harvested by centrifugation (3000 g, 10 min) and resuspended in 0.1 M lithium acetate, 10 mM Tris-HCl, pH 7.5, 1 mM EDTA, 10 mM DTT. After 1 h incubation at 30°C, cells were washed twice with ice-cold 1 M sorbitol and resuspended in 1–5 mL ice-cold 1 M sorbitol. Plasmid DNA (0.5–1 µg DNA) was added to a 180 µL cell suspension and transferred to a 0.2 cm gap width electroporation cuvette. Electroporation was performed using a Bio-Rad gene pulser with an electric pulse of 2.5 kV, 25 µF and 200 Ω. Cells were immediately washed out from the cuvettes after the electroporation with YPD, plated on selective CSM-URA medium and incubated at 28°C for 3–5 days. Fresh colonies were then grown 3 h in liquid YPD at 30°C and meiosis was induced by pouring the suspension cell on ACK medium (10 g/L potassium acetate, 2.5 g/L yeast extract, 20 g/L agar). These plates were incubated 1 week at 20°C and haploids were finally selected on MMAS medium plates [20 g/L galactose, 7 g/L YNB (WA), 0.6 g/L DOB-LEU-HIS-ARG-URA, 0.06 g/L L-canavanine, 0.2 g/L G418, 20 g/L agar] supplemented with 10 µM *trans*-zeatin at 28°C for 3–5 days. Suspensions of transformants were then spotted onto dropout media containing or not 10 µM of *trans*-zeatin with 2% galactose and grown for 48 h at 28°C. For specificity and sensitivity assays, complemented-yeast growth was carried out in liquid YCGal (7 g/L YNB, 0.8 g/L CSM-URA, 20 g/L galactose) supplemented with various types and concentrations of cytokinins for 48 h at 28°C. Cell growth was measured at 630 nm (BioHit Reader BP800).

## Chemicals

Pure standards of isopentenyladenosine 5'-monophosphate, *trans*-zeatin riboside 5'-monophosphate, *cis*-zeatin riboside 5'-monophosphate, isopentenyladenine, *trans*-zeatin, *cis*-zeatin, dihydrozeatin, isopentenyladenosine, *trans*-zeatin riboside, *cis*-zeatin riboside, dihydrozeatin riboside, dihydrozeatin

riboside 5'-monophosphate, *trans*-zeatin O-glucoside, *trans*-zeatin O-glucoside riboside, *trans*-zeatin N7-glucoside, 2-methylthio-isopentyladenine, 2-methylthio-*trans*-zeatin, 2-methylthio-*cis*-zeatin, 2-methylthio-isopentyladenosine, 2-methylthio-*trans*-zeatin riboside, 2-methylthio-*cis*-zeatin riboside were purchased from Olchemim (Olomouc, Czechia).

## RNA Isolation and Gene Expression Analysis

Extraction of total RNA from *M. domestica* organs was performed using the NucleoSpin RNA extraction kit (Marcherey-Nagel), with improved lysis step (McKenzie et al., 1997). First-strand cDNA were synthesized from 1 µg of total RNA using the iScript cDNA Synthesis Kit (Bio-Rad). Quantitative real-time PCR measurements were carried out in triplicate using SsoAdvanced Universal SYBR Green (Bio-Rad) in a 15 µL final volume containing 6 µL diluted template cDNA and specific primers (0.5 µM) (Supplementary Table S1). Amplification was performed on a CFX96 Touch real-time PCR system (Bio-Rad) with the following conditions: 95°C for 7 min and 40 cycles at 95°C for 10 s and 60°C for 40 s. Amplification was followed by a melt curve analysis. Absolute quantification of transcript copy number was assessed with calibration curves. Transcript levels were then normalized with EF1α.

## Subcellular Localization Experiments

Subcellular localization of MdCHK receptors were studied in *Catharanthus roseus* C20D cells transiently transformed using plasmid-coated particles bombardment as described in Guirimand et al. (2009). The full length MdCHK sequences were amplified and cloned into the *SpeI* restriction site of pSCA-YFP plasmid (for primers, see Supplementary Table S1), in frame with the 5' extremity of the YFP coding sequence. The endoplasmic reticulum (ER) cyan fluorescent protein (CFP) marker (Guirimand et al., 2010) was used in co-transformation assays.

Dynamic localization of MdCHK receptors was also studied in yeast *S. cerevisiae* strain WT303 (*MAT*a/α, *leu2*, *trp1*, *ura3*, *ade2*, *his3*) transformed by pESC-LEU plasmids (Foureau et al., 2016) containing *MdCHK* sequences, except for MdCHK2, cloned in pYES2 and transformed in *sln1*Δ yeast strain to bypass the sequence toxicity in microorganisms. The CYP450 T16H2 sequence was clone in the pESC-TRP plasmid in fusion with the 5'end of the CFP sequence and used as an ER marker (Besseau et al., 2013). *MdCHK* sequences were cloned under galactose inducible promoter and fused at the 5' end with the YFP sequence. Transformed colonies were cultivated on selective plates (CSM-LEU or CSM-URA, supplemented by 2% glucose, respectively, for pESC-LEU and pYES2) at 30°C for 48 h and then transferred in inducing liquid media (CSM-LEU or CSM-URA, supplemented by 2% galactose, respectively, for pESC-LEU and pYES2) with or without iP (5 µM) for additional overnight culture at 28°C.

An Olympus BX51 epifluorescence microscope equipped with the Olympus DP71 digital camera and Cell\*D imaging software (Soft Imaging System Olympus) was used for image capture and

<footnote n="1">¹http://www-nmr.cabm.rutgers.edu/bioinformatics/Proteomic_tools/Helical_ wheel/</footnote>
merging false-colored images of both C20D cells and *S. cerevisiae* colonies expressing YFP.

## BiFC Interaction Assays

Bimolecular Fluorescent Complementation (BiFC) experiments were conducted using SPYNE and SPYCE plasmids (Waadt and Kudla, 2008). MdCHK sequences were amplified and cloned into the *SpeI* restriction site (for primers, see Supplementary Table S1), in frame with the 5' extremity of a truncated YFP coding sequence. Transient transformation of *C. roseus* cells by particle bombardment and YFP imaging were performed according to Guirimand et al. (2009) with adaptation for BiFC assays (Guirimand et al., 2010). Interactions were tested in triplicates using three independent plasmid clones.

## RNA-seq Data Analysis

Available RNAseq data for *M. domestica* was downloaded from NCBI via SRA toolkit2.6.2. The recovered SRA files were transformed in fastq format with the "fastq-dump" command from SRA toolkit. The files were cleaned with Trimmomatic 0.36 with default parameters and using provided adapter sequences for TruSeq2 and TruSeq3. The transcription quantification was performed with Salmon 0.6.1 using the Variational Bayesian EM algorithm and biase correction. TPM (transcripts per million) values from the resulting quant.sf files were combined under R 3.3.0 in an expression matrix containing 95,232 predicted genes (*Malus domestica* v3.0) × 250 experimental conditions. Using the expression matrix, Pearson Correlation Coefficients (PCC) and further Highest Reciprocal Ranks (HRR) computation were performed using a homemade program written in C [HRR (gene A, gene B)] = max [rank (gene A, gene B), rank (gene B, gene A)] to establish the co-expressed genes lists for MdCHK and GO enrichment tests. For each MdCHK highly co-expressed genes, i.e., genes with a HRR ≤ 500 were selected. The procedure was repeated with publicly available *A. thaliana* RNAseq data. We similarly prepared an expression matrix containing 33,604 transcripts (*Arabidopsis* TAIR v10 genome annotation) and 1,676 samples. Co-expressed genes lists were obtained after calculating PCC and ranking them with HRR. For each AtCHK (AT5G35750.1, AT1G27320.1, AT2G01830.2), gene pairs having an HRR < 500 were considered to be significantly co-expressed. Orthology between *Arabidopsis* and apple tree was obtained from Plaza 3.0 (Van Bel et al., 2012). The functions represented by coexpressed genes of each MdCHK were analyzed with the Gene Ontology classification. A BlastX was performed on the *M. domestica* genome v3.0 to recover correspondent protein sequences and Pfam domains were identified using Hmmer. The functional annotations of the *M. domestica* genome v3.0 were generated using Trinotate on the previous data. In order to determine potential functional enrichment for every target gene, enrichment of GO terms was tested by comparing effectives to a hypergeometric distribution (p-value cut-off = 0.001) using the R "phyper" function. To compare redundancies in the five co-expressed genes lists, a Venn diagram was drawn using the venneuler package 1.1 for R.

# RESULTS

## Identification of Five *Malus domestica* CHASE-Containing Histidine Kinases (MdCHKs)

Based on the CHASE domain of the three *A. thaliana* CHKs, apple tree genome was browsed to identify putative CHK sequences in *M. domestica*. Five candidates were identified and named MdCHK2 (locus tag MDP0000258078), MdCHK3a (locus tag MDP0000310800), MdCHK3b (locus tag MDP0000155347), MdCHK4a (locus tag MDP0000151825) and MdCHK4b (locus tag MDP0000242242) according to their distribution within the three classical groups homologous to AHK2, AHK3, and AHK4, as shown by phylogenetic analysis (Figure 1). Genomic sequences revealed that MdCHK2 possesses 13 exons and is located on chromosome 9 (Supplementary Figure S1). Both MdCHK3a and MdCHK3b, respectively, located on chromosomes 16 and 13, display 10 exons with a similar organization regarding

<figure><img src="image_3.png" /><figcaption>**FIGURE 1** | Phylogenetic analysis of *Malus domestica* cytokinin receptors. Apple tree CHASE Histidine Kinase receptors MdCHK2 (KM114879), MdCHK3a (KM114880), MdCHK3b (KM114881), MdCHK4a (KM114883), and MdCHK4b (KM114882) were compared to cytokinin receptors characterized in angiosperms. The tree was constructed by neighbor-joining distance analysis on conserved protein sequence domains. Brassica napus: BnCHK1 (KF621029), BnCHK2 (KF621030), BnCHK3 (KF621031), BnCHK4 (KF621032), BnCHK5 (KF621033); *Arabidopsis thaliana*: AHK2 (At5g35750), AHK3 (At1g27320), AHK4 (At2g01830); *Oryza sativa*: OsHK5 (Os02g50480), OsHK2 (Os10g21810), OsHK3 (Os01g69920), OsHK4 (Os03g50860); *Vitis vinifera*: VvCyt3 (CAO42401), VvCyt1 (GSVIVT01030058001), VvCyt2 (CAO66151); *Zea mays*: ZmHK1 (NP_001104859), ZmHK2 (NP_001104866), ZMHK3a (AB102957), ZmHK1a2 (NP_001105857), ZmHK1b1 (NP_001105858), ZmHK1b2 (NP_001105913).</figcaption></figure>
<figure><img src="image_4.png" /><figcaption>**FIGURE 2** | Domain structure of the five MdCHK receptors. MdCHK receptors display a classical CHK organization with a CHASE sensing domain, transmembrane domains (TM), a histidine kinase (HK) domain composed of a HK dimerization and phosphoacceptor domain (HisKA) and a HK catalytic domain (HATPase), a receiver (REC)-like domain and a REC domain. The conserved residues of each domain are indicated (arrows).</figcaption></figure>

intron positions and intron/exon sizes. The same observation was made for *MdCHK4a* and *MdCHK4b* respectively located on chromosomes 13 and 10 and containing 11 exons. Such similarity may reflect the gene duplication events leading to the couples of CHKs homologous *MdCHK3a/MdCHK3b* and *MdCHK4a/MdCHK4b* which share respectively 94.67 and 95.67% nucleotide identity. The full-length cDNAs of the five MdCHKs were cloned and deposited at NCBI under the GenBank accession numbers KM114879 to KM114883. They contained large open-reading frames ranging from 3027 to 3612 bp encoding proteins of 1008 to 1203 amino acids.

The computational analysis of the protein sequences of MdCHKs revealed the presence of the four basic conserved domains called Cyclases/Histidine kinases Associated SEnsory (CHASE), Histidine Kinase (HK), Receiver (REC), and Receiver-like (REC-like) (Figure 2 and Supplementary Figure S2). While these modular MdCHKs share a similar multidomain architecture in their cytoplasmic C-terminal part (HK, REC-like and REC domains), they differ in the TM domain topology of the N-terminus part. Indeed, the MdCHK4a and MdCHK4b have two predicted α-helices transmembrane domains TM1 and TM2 (Supplementary Figure S3) bordering the CHASE domain (Figure 2). MdCHK3a and MdCHK3b possess equivalent TM1/CHASE/TM2 organization with an additional predicted transmembrane helice (TM3) in N-terminus. Finally, the N-terminus of MdCHK2 includes a fourth predicted transmembrane domain

(TM4). Consequently, in addition of the extracytoplasmic loop containing the sensing CHASE domain, the TM3 and TM4 domains border a supplemental extracytoplasmic loop (140 aa) absent in MdCHK3a/b and MdCHK4a/b structures (Figure 2). Regarding the cytoplasmic part, each MdCHK possesses an HK domain containing the conserved phosphorylatable histidine residue as well as the C-terminal REC domain including the conserved phospho-accepting aspartate residue (Figure 2 and Supplementary Figure S2). In addition, the five receptors contain a second receiver domain located between the HK domain and the REC domain called pseudo-receiver domain or REC-like (Figure 2). Interestingly, the putative phospho-accepting aspartate residue in the REC-like domain of MdCHK2 is conserved, suggesting that this domain may be functional in terms of phosphorelay reaction, whereas the corresponding residues in MdCHK3a/MdCHK3b and MdCHK4a/MdCHK4b are substituted with glutamate (Supplementary Figure S2).

## MdCHKs Function as Cytokinin Receptors in a Yeast Complementation Assay

To assess the function of the five MdCHKs as cytokinin receptors, we exploited the *S. cerevisiae sln1Δ* deletion mutant strain which carries a lethal mutation in the *SLN1* gene encoding its unique osmosensing histidine kinase (Inoue et al., 2001).
<table><thead><tr><th rowspan="2"></th><th rowspan="2"></th><th colspan="5">B</th></tr><tr><th>MdCHK2</th><th>MdCHK3a</th><th>MdCHK3b</th><th>MdCHK4a</th><th>MdCHK4b</th></tr></thead><tbody><tr><td rowspan="5">A</td><td rowspan="5">sln1Δ</td><td>- tZ</td><td>- tZ</td><td>- tZ</td><td>- tZ</td><td>- tZ</td></tr><tr><td>MdCHK2</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>MdCHK4a</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>MdCHK4b</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>MdCHK3a</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="5">B</td><td rowspan="5"></td><td>Free-base</td><td>iP</td><td>0,1 nM</td><td>500 nM</td><td>250 nM</td><td>50 nM</td><td>50 nM</td></tr><tr><td></td><td>IZ</td><td>1 nM</td><td>100 nM</td><td>500 nM</td><td>500 nM</td><td>500 nM</td></tr><tr><td></td><td>cZ</td><td>100 nM</td><td>*</td><td>10 μM</td><td>*</td><td>10 μM</td></tr><tr><td></td><td>DHZ</td><td>100 nM</td><td>5 μM</td><td>5 μM</td><td>*</td><td>10 μM</td></tr><tr><td></td><td>Free-sides</td><td>iPA</td><td>100 nM</td><td>*</td><td>*</td><td>*</td><td>*</td></tr><tr><td rowspan="3">A</td><td rowspan="3"></td><td rowspan="3"></td><td>iZR</td><td>500 nM</td><td>*</td><td>*</td><td>*</td><td>*</td></tr><tr><td>cZR</td><td>1 μM</td><td>*</td><td>*</td><td>*</td><td>*</td></tr><tr><td>DHZR</td><td>10 μM</td><td>*</td><td>*</td><td>*</td><td>*</td></tr><tr><td rowspan="4">B</td><td rowspan="4"></td><td rowspan="4"></td><td>iPRMP</td><td>1 μM</td><td>*</td><td>*</td><td>*</td><td>*</td></tr><tr><td>iZRMP</td><td>5 μM</td><td>*</td><td>*</td><td>*</td><td>*</td></tr><tr><td>cZRMP</td><td>1 μM</td><td>*</td><td>*</td><td>*</td><td>*</td></tr><tr><td>DHZRMP</td><td>10 μM</td><td>*</td><td>*</td><td>*</td><td>*</td></tr><tr><td rowspan="4">B</td><td rowspan="4"></td><td rowspan="4"></td><td>Gluco</td><td>tZOG</td><td>100 nM</td><td>*</td><td>*</td><td>*</td><td>*</td></tr><tr><td>tZOGR</td><td>5 μM</td><td>*</td><td>*</td><td>*</td><td>*</td></tr><tr><td>tZ7G</td><td>*</td><td>*</td><td>*</td><td>*</td><td>*</td></tr><tr><td>2MeSCK</td><td>2MeSiP</td><td>0,1 nM</td><td>5 μM</td><td>500 nM</td><td>5 nM</td><td>5 nM</td></tr><tr><td rowspan="5">A</td><td rowspan="5"></td><td rowspan="5"></td><td>2MeStZ</td><td>5 nM</td><td>5 μM</td><td>5 μM</td><td>500 nM</td><td>250 nM</td></tr><tr><td>2MeScZ</td><td>25 nM</td><td>*</td><td>*</td><td>5 μM</td><td>5 μM</td></tr><tr><td>2MeSiPA</td><td>25 nM</td><td>*</td><td>*</td><td>5 μM</td><td>1 μM</td></tr><tr><td>2MeStZR</td><td>1 μM</td><td>*</td><td>*</td><td>*</td><td>*</td></tr><tr><td>2MeScZR</td><td>1 μM</td><td>*</td><td>*</td><td>*</td><td>*</td></tr></tbody></table>

**FIGURE 3** | sln1Δ mutated yeast complementation assays with MdCHK. (A) Cytokinin-dependant growth phenotype of mutant yeast sln1Δ transformed with pYES2-MdCHK. Suspensions of transformed yeasts were spotted onto plates containing or not tZ (10 µM). (B) Cytokinin binding specificity and sensitivity profiles of MdCHKs. MdCHK-complemented yeasts were grown in liquid culture medium containing different cytokinin types with different concentrations for 48 h. Minimum concentration inducing yeast growth is indicated for each cytokinin tested. iP, isopentenyladenine; tZ, trans-zeatin; cZ, cis-zeatin; DHZ, dihydrozeatin; iPA, isopentenyladenosine; tZR, trans-zeatin riboside; cZR, cis-zeatin riboside; DHZR, dihydrozeatin riboside; iPRMP, isopentenyladenosine 5'-monophosphate; tZRMP, trans-zeatin riboside 5'-monophosphate; cZRMP, cis-zeatin riboside 5'-monophosphate; DHZ, dihydrozeatin riboside 5'-monophosphate; tZOG, trans-zeatin O-glucoside; tZOGR; trans-zeatin O-glucoside riboside; tZ7G, trans-zeatin N7-glucoside; 2MeSiP, 2-methylthio-isopentenyladenine; 2MeStZ, 2-methylthio-trans-zeatin; 2MeScZ, 2-methylthio-cis-zeatin; 2MeSiPA, 2-methylthio-isopentenyladenosine; 2MeStZR, 2-methylthio-trans-zeatin riboside; 2MeScZR, 2-methylthio-cis-zeatin riboside. *, no growth observed. Basal growth of MdCHK3a- and MdCHK3b-complemented yeast is also indicated (●).

The sln1Δ yeast mutants carrying pYES2-MdCHK2, pYES2-MdCHK4a or pYES2-MdCHK4b were lethal (Figure 3A). However, the addition of trans-zeatin (tZ) in culture medium allowed recovering a normal growth (Figure 3A). By depending on the presence and perception of cytokinin to complement sln1 mutation, the three recombinant yeast strains clearly demonstrated that MdCHK2, MdCHK4a, and MdCHK4b act as cytokinin receptors in this heterologous system. Concerning the recombinant strains carrying pYES2-MdCHK3a and pYES2-MdCHK3b, they displayed an original phenotype since they exhibited a basal growth in absence of cytokinin, especially for MdCHK3b (Figure 3A and Supplementary Figure S4). However, the addition of tZ clearly induced the yeast growth pointing out that both MdCHK3a and MdCHK3b sense cytokinins (Figure 3A and Supplementary Figure S4). The basal constitutive activity of these two cytokinin receptors raised the question of their putative additional sensing function.

## MdCHK Receptors Show Different Binding Specificities toward Cytokinin

We further evaluated the substrate specificity and sensitivity of M. domestica cytokinin receptors. In this way, we measured the growth of the yeast cells in presence of various cytokinin-types at different concentrations including free-bases, ribosides, glucosides, and methylthio cytokinins. The specificity as well as the sensitivity (the minimal cytokinin concentration that induced yeast growth) were reported on Figure 3B and Supplementary Figure S4. Three distinct profiles corresponding to MdCHK2, MdCHK3, and MdCHK4 groups were observed.

First, MdCHK2 clearly perceived a wide range of cytokinin forms since each cytokinin type activated the receptor except tZ7G. Furthermore, this receptor presented a remarkable higher sensitivity than other MdCHKs (until 0,1 nM for iP and 2MeSiP forms) and was the only one to be activated by some ribosides and glucosides cytokinin-types. Secondly, the strictly cytokinin-dependent MdCHK4a and MdCHK4b were activated by the free-bases iP and tZ and the methylthio-forms 2MeSiP, 2MeStZ, 2MeScZ, and 2MeSiPA (Supplementary Figure S4). High concentrations (10 µM) of cZ and DHZ were also effective on MdCHK4b. Nevertheless, the cytokinin-sensitivity of both MdCHK4a and MdCHK4b was obviously lower than MdCHK2. Finally, even if MdCHK3a and MdCHK3b showed a basal growth in absence of cytokinin, we were able to detect a significant difference of growth in presence of iP, tZ, DHZ, 2MeSiP, and 2MeStZ (Supplementary Figure S4).

## MdCHKs Exhibit Distinct Expression Patterns

To examine the gene expression of MdCHKs, RT-qPCR was carried out using distinct plant organs including roots, stems, leaves, flower buds and flowers. Transcripts of the five MdCHKs were detected in all the tested organs, but with distinct expression pattern. Thus, MdCHK2 reached higher expression level in leaves and stems (Figure 4). MdCHK4a and MdCHK4b displayed similar expression profiles with high expression level in stems whereas gene expression was hardly detected in flowers (Figure 4). MdCHK3a and MdCHK3b disclosed differential pattern of expression. While MdCHK3a was mainly
Normalized transcript copy numbers

<figure><img src="image_5.png" /><figcaption>**FIGURE 4** | MdCHKs display various expression patterns in apple tree organs. Total transcript numbers of *MdCHK2*, *MdCHK3a*, *MdCHK3b*, *MdCHK4a*, and *MdCHK4b* were determined by real-time RT-PCR analyses performed on total RNA extracted from various *M. domestica* organs. R, roots; S, stems; L, leaves; FB, flower buds; F, flowers. Transcript abundance of each gene was normalized against *EF1α* reference gene. Data have been analyzed by ANOVA (p < 0.05), after running Shapiro-Wilk and Bartlett tests, followed by HSD Tukey test. Error bars were calculated from triplicates.</figcaption></figure>

expressed in roots, *MdCHK3b* showed its highest expression in flowers (Figure 4). Finally, substantial differences in the overall expression level of each *MdCHKs* were also observed. Interestingly, *MdCHK2* displayed the higher expression level whilst *MdCHK3b/MdCHK4a* and *MdCHK3a/MdCHK4b* retained a 10- and a 100-fold lower expression, respectively.

To complete qPCR analysis, we used the available RNAseq data to generate an expression matrix of apple tree genes (Supplementary Table S2). The best co-expressed genes with each MdCHK through our expression matrix were investigated and compared (Supplementary Table S3). A very weak degree of overlap among lists of genes co-expressed with each MdCHK was found. MdCHK2 shared up to 18 genes with other cytokinin receptors. MdCHK3 and MdCHK4 groups did not have common co-expressed genes. MdCHK3a/MdCHK3b homologs as well as MdCHK4a/MdCHK4b homologs shared respectively 88 and 26 genes (Figure 5). This low overlapping of co-associated genes might support the limited functional redundancy of MdCHKs. Each MdCHK co-expressed genes list was compared to the list established for their *Arabidopsis* ortholog in order to highlight potential shared genes. We

found a relatively weak overlap between functions associated to either CHK (Supplementary Figure S5A and Table S4). For example, only 19 genes were similarly co-expressed between MdCHK2 (1.8%) and AtCHK2 (8.7%). While such a weak overlap could be due to the initial datasets used to calculate correlations which differ in size and experiments, conserved co-expressed genes may be good candidates for a further investigation of the cytokinin pathway. In addition, we found very small overlaps between co-expressed gene lists of AtCHKs, as observed for MdCHKs, reinforcing a potential specificity in CHK functions (Supplementary Figure S5B). Enrichment tests of Gene Ontology (GO) terms performed on each list of co-expressed genes also gave an overview of possible specific physiological processes associated with each receptor. For example, “Embryo development ending in seed dormancy” and “Response to cadmium ion” were exclusively enriched for MdCHK2 whereas “Plant-type secondary cell wall biogenesis” and “Regulation of growth” were specifically enriched for MdCHK4a and MdCHK3b, respectively (Figure 6 and Supplementary Figure S6).

## MdCHKs Mainly Localize to the Endoplasmic Reticulum in Plant and Show a Dynamic Behavior in Response to Cytokinin in Yeast

We investigated the subcellular distribution of the MdCHKs using C-terminal YFP tagging to ensure the correct anchoring of the transmembrane domains. MdCHK-YFP constructs were transiently expressed in C. roseus cells that constitute a reliable model for studying protein subcellular localization (Foureau et al., 2016). In transiently transformed cells, the fusion proteins displayed a fluorescence signal located in

<figure><img src="image_6.png" /><figcaption>**FIGURE 5** | Venn diagram of MdCHK co-expressed gene lists. Best co-expressed gene lists (HRR < 500) were separately established for each MdCHK and compared (see Materials and Methods).</figcaption></figure>
<figure><img src="image_7.png" /></figure>

**FIGURE 6** | Functional annotation of best co-expressed genes (HRR < 500) with Gene Ontology (GO) terms. Only GO terms from << Biological process >> represented by 5 or more genes were kept. P-values correspond to GO term enrichment tests which were performed by comparison to hypergeometrical distribution.

the endoplasmic reticulum (ER) network throughout the cell as well as in the perinuclear space (Figures 7A,E,I,M,Q). The signal perfectly co-localized (Figures 7C,G,K,O,S) with the specific ER-CFP marker (Figures 7B,F,J,N,R), confirming that the five MdCHKs are located in the ER of plant cells. Noteworthy, cytokinin addition in plant cell medium did not alter the ER localization of the MdCHK-YFP fusions (data not shown).

Since plant cells may produce their own pool of cytokinins preventing the study influence of exogenous cytokinins on the localization of the MdCHKs, we therefore investigated the subcellular distribution of MdCHKs in the yeast *S. cerevisiae*,

by using YFP fusion proteins. In absence of cytokinins, a punctate fluorescence pattern was observed for the five MdCHKs (Figures 8G1,I1,K1,M1,O1), which accumulated in the ER forming structures comparable as organized smooth ER (Snapp et al., 2003) that is described to result from protein interactions. Upon cytokinin treatment, the ER localization of the five MdCHKs did not change. But interestingly, a reorganization of the fluorescent pattern was observed reflecting the decrease or disappearance of aggregate structures (Figures 8H1,J1,L1,N1,P1). Indeed, MdCHK2 and MdCHK4b displayed a strong perinuclear localization (Figures 8H1,P1). Concerning MdCHK3a, MdCHK3b, and MdCHK4a, the
<figure><img src="image_8.png" /><figcaption>**FIGURE 7** | MdCHK-YFP fusion proteins are localized at the Endoplasmic Reticulum (ER) membrane. *Catharanthus roseus* cells were transiently co-transformed with plasmids expressing MdCHK-YFP (*A*,*E*,*I*,*M*,*Q*) and endoplasmic reticulum-CFP marker (*B*,*F*,*J*,*N*,*R*). Co-localization of the two fluorescence signals appeared on the merged image (*C*,*G*,*K*,*O*,*S*). The morphology is observed with differential interference contrast (DIC; *D*,*H*,*L*,*P*,*T*). Scale bars: 10 µm.</figcaption></figure>

fluorescence signal appeared in a discontinuous pattern as well as in the perinuclear space (Figures 8J1,L1,N1). In order to ensure that cytokinins themselves had no impact on the architecture of the ER, we used T16H2-CFP construct as a specific ER marker (Besseau et al., 2013). Upon cytokinin treatment, no redistribution of fluorescence signal was observed, reinforcing the plausibility of a specific reorganization of MdCHKs in response to cytokinin signal (Figures 8E1,F1).

## Homo- and Heterodimerization Are Common Features of MdCHK2 and MdCHK4 in Contrast to the MdCHK3 Pair That Only Displays Specific Heterodimerization Characteristics

Cytokinin receptors were previously proposed to interact each other to enable the *trans*-phosphorylation of the HK domain after cytokinin perception (Dortay et al., 2006;
<figure><img src="image_9.png" /></figure>

Caesar et al., 2011; Hothorn et al., 2011; Wulfetange et al., 2011). Considering the multiple possibilities of interactions between the five MdCHKs, homo- and hetero-dimerization were investigated by BiFC assays in planta. The full coding

sequences of MdCHKs were cloned upstream of the coding sequence of the two split-YFP fragments (YFPN and YFPc) to generate the MdCHK-YFPN and MdCHK-YFPc fusion proteins. BiFC analysis revealed that MdCHK2, MdCHK4a, and MdCHK4b were able to form homodimers within the ER network (Figures 9A,S,Y) whereas no BiFC complex reconstitution was observed when testing the MdCHK3a and MdCHK3b homodimers (Figures 9G,M). Moreover, no signal was detected with the MdCHK3a and MdCHK3b heterodimer combination (Figures 9H,L). Additionally, a cytokinin application did not result in the formation of a fluorescent signal within the three configurations (data not shown). Nevertheless, MdCHK3a and MdCHK3b were able to heterodimerize with MdCHK2, MdCHK4a and MdCHK4b within the ER (Figures 9B,C,F,I-K,N,O,Q,R,V,W). In addition, MdCHK2 and both MdCHK4 homologs shape heterodimers with each other (Figures 9D,E,P,U,T,X).

DISCUSSION

Five CHASE domain-containing Histidine Kinases belonging to the three classical groups were identified in *M. domestica*: *MdCHK2*, the homolog of *Arabidopsis AHK2* and the two pairs *MdCHK3a/MdCHK3b* and *MdCHK4a/MdCHK4b*, homologous to *AHK3* and *AHK4*, respectively. These two pairs might be the direct consequence of apple genome-wide duplication (Velasco et al., 2010) as it was found in another hardwood tree *Populus trichocarpa* (Nieminen et al., 2008).

To confirm the functionality of MdCHKs as cytokinin receptors, we conducted a cytokinin-responsive assay based on the use of the *sln1Δ S. cerevisiae* strain mutant (Maeda et al., 1994). We thus showed that MdCHK2, MdCHK4a, and MdCHK4b restore the viability of the *sln1Δ* mutant in a strictly cytokinin-dependent fashion providing convincing evidences of their cytokinin receptor function. By contrast, both MdCHK3a and MdCHK3b, they presented an unexpected profile since they conferred a basal growth to the *sln1Δ* mutant in absence of cytokinin. However, the enhancement of the yeast growth in presence of cytokinin confirmed their cytokinin receptor function. It is important to emphasize that the basal constitutive activity confers an originality for MdCHK3a and MdCHK3b compared to MdCHK2, MdCHK4a, and MdCHK4b. To our knowledge, a constitutive activity for a CHK protein has never been reported before.

The extensive exploitation of our five MdCHK-complemented yeast strains revealed that the receptors differed greatly in their cytokinin specificity and sensitivity. Three distinct specificity and sensitivity profiles clearly emerged. The most remarkable result comes from MdCHK2 which perceives an unprecedented range of cytokinins including nucleotide-type precursors as well as free-base-, riboside-, O-glucoside- and methylthio-forms, with a substantial sensitivity for iP, tZ, 2MeSiP, and 2MeStZ (Figure 3B and Supplementary Figure S4). Until now, due to its presumed toxicity in bacteria, few studies examined the ligand-binding properties of the full-length *Arabidopsis* AHK2 and its homologs in other plant species. Recently, an *E. coli* cytokinin-binding assay
<figure><img src="image_10.png" /><figcaption>**FIGURE 9** | MdCHKs interact to form homodimers and heterodimers in BiFC assays. Cells of *C. roseus* were transiently co-transformed using the plasmids expressing the different MdCHK proteins fused with the YFP<sup>N</sup> and YFP<sup>C</sup> split in C-terminal. Homodimerizations (*A*,*G*,*M*,*S*,*Y*) and heterodimerizations (*B*–*F*,*H*–*L*,*N*–*R*,*T*–*X*) have been tested, with both YFP<sup>N</sup> and YFP<sup>C</sup> combinations. Three independent plasmid clones were used to test interactions. Scale bars: 10 µm.</figcaption></figure>

revealed that the full-length BnCHK1 and BnCHK3, two AHK2 homologs in *Brassica napus*, showed high affinity for *tZ*, iP, and *tZR* (Kuderová et al., 2015). A tobacco membrane assay also revealed that *tZ* and iP strongly interacted with AHK2, whereas their conjugated forms did not, suggesting that free bases were the sole biologically active cytokinin compounds (Lomin et al., 2015). Even if our experiments used a heterologous system in which the yeast membrane environment can potentially differ from those of plant, we highlighted not only the receptor ability to bind hormones but also their activation through cytokinin perception activating the phosphorelay in yeast. Thus, we can assume that our assay reflects the biological activity of cytokinin riboside and nucleotide forms on the receptor. Moreover, artefactual cytokinin activation or conversion occurrence in yeast can be omitted, since MdCHK4a and MdCHK4b are not activated in presence of the nucleotide-, riboside- or O-glucoside types. Therefore, the broad cytokinin spectrum of MdCHK2 raises the question of

its possible central role in *M. domestica*. Such hypothesis is also reinforced by the high expression level of MdCHK2 compared to other MdCHKs (Figure 4). Concerning MdCHK4a and MdCHK4b, they perceived a restricted spectrum of cytokinins such as free-base forms and some of the methylthiolated-forms with a lower sensitivity than MdCHK2 supporting previous works obtained with *Arabidopsis* AHK4 and AHK2 (Stolz et al., 2011). Concerning MdCHK3a and MdCHK3b, they perceived some cytokinin free-bases and methylthiolated forms. However, regarding their constitutive activity, we must consider that we could have under estimated their real cytokinin binding effectiveness and that it might not really reflect their complete capacity to perceive the diverse structures of cytokinins. The original activities of MdCHK3a and MdCHK3b were definitely interesting and need to be further investigated. In particular, did the monomeric or heterodimeric forms of MdCHK3a/b influence their cytokinin perception? In any case, our results
clearly supported previous works reporting that cytokinin-
binding properties of AHK3 differed from those of AHK2 and
AHK4 (Spíchal et al., 2004; Romanov et al., 2006; Stolz et al., 2011;
Heyl et al., 2012).

These different properties between MdCHK receptors raised
the question of their specialized functions in apple tree.
MdCHKs clearly showed an organ-specific gene expression
pattern (Figure 4) and the analysis of gene co-expression profiles
with each MdCHK unequivocally shed light on the weak degree
of overlap among lists of co-expressed genes that strongly
underlines the distinct roles of the five receptors in physiological
processes (Figure 5 and Supplementary Figure S6). Besides these
distinct patterns, we can point out that both pairs of homologs
MdCHK3 and MdCHK4 might acquire distinct functions after
duplication.

From a structural point of view, the five MdCHKs possess
the conserved CHASE, HK, REC-like and REC domains also
found in *Arabidopsis*, maize and rapeseed (Ueguchi et al.,
2001; Yonekura-Sakakibara et al., 2004; Kuderová et al.,
2015). While the architecture of cytoplasmic C-terminal part
is similar within MdCHK, the topology of their N-terminus
differs in the number of predicted transmembrane domains
which surround the CHASE sensing domain (Figure 2).
MdCHK4a/MdCHK4b and MdCHK3a/MdCHK3b possess
two and three transmembrane domains, respectively, whereas
MdCHK2 exhibits a fourth transmembrane helix that forms a
unique additional extracytoplasmic loop. This variability might
reflect the specific sensing activities of MdCHK receptors. It
is well established that structural variations in the CHASE
domain result in different ligand specificities of *Arabidopsis*
receptors (Romanov et al., 2006; Heyl et al., 2007; Stolz et al.,
2011). Nevertheless, the organization of the surrounding
environment of the CHASE domain might also be important for
receptor functioning (Steklov et al., 2013). Indeed, the CHASE
flanking regions including transmembrane helices are assumed
to play a substantial role in localization and intramolecular
signaling (Steklov et al., 2013). Directed mutagenesis on
a transmembrane helix highlighted its importance for the
AHK4 receptor activation (Miwa et al., 2007). Moreover,
experiments on the CHASE domain of AHK4 compared with
the full-length receptor revealed significant differences in the
binding affinity, highlighting the importance of the CHASE
environment in cytokinin perception (Stolz et al., 2011).
Thus, it cannot be excluded that architectural variations in
the N-terminal part of the MdCHKs somehow influence their
distinct properties.

Regarding the cytoplasmic C-terminal part, only MdCHK2
harbors a phospho-accepting aspartate in the receiver-like
domain (Supplementary Figure S2). Its presence is a common
feature shared with AHK2 and the four AHK2 homologs
in *B. napus* BnCHK1, BnCHK2, BnCHK3, and BnCHK4
(Ueguchi et al., 2001; Kuderová et al., 2015). However,
the function of this receiver-like domain has not been yet
elucidated. As reported above, MdCHK2 perceives a broad
cytokinin spectrum with a substantial sensitivity compared
to other MdCHKs, thus it would be interesting to further
investigate if the second phosphorylatable aspartate contribute

to the properties of this receptor. Indeed, if we consider
that the MdCHK2 receiver-like domain is functional, it might
optimize the phosphorelay reaction in addition to the receiver
domain or might guide specific interactions with downstream
HPts.

As previously described in other species, the five MdCHKs
are located at the ER membrane and the perinuclear space.
Moreover, we revealed for the first time, a dynamic redistribution
of cytokinin receptors in response to cytokinin application.
The development of an approach in *S. cerevisiae* allowed us
to overcome the use of plant cells, which probably produce
their own cytokinins and prevent studying cytokinin influence.
More precisely, we showed that MdCHKs relocalized through
the ER network, especially by getting closer to the nucleus
(Figure 8). This result clearly supports the current concept of
phosphotransfer enhancement through a perinuclear localization
which overcomes intracellular distance and optimizes the
signal transduction (Caesar et al., 2011; Wulfetange et al.,
2011).

We also report herein a complete analysis of the full-
length cytokinin receptors interactions *in planta*. Until now,
only homodimerization of full-length AHK2 was demonstrated
*in planta* (Wulfetange et al., 2011) and a partial study
in yeast two-hybrid system based on full-length receptors
showed the AHK3/AHK4 interaction as well as the formation
of AHK3 homo-oligomers (Caesar et al., 2011). Here, we
examined the homo- and the hetero-dimerization of the
MdCHKs *in planta* since histidine kinases are supposed
to act as dimers. Not only MdCHK2, MdCHK4a, and
MdCHK4b homodimerize, but they form heterodimers with
each other. Surprisingly, MdCHK3a and MdCHK3b do not
form homodimers. This feature might explain their singularity
in cytokinin perception in our yeast system compared to
MdCHK2, MdCHK4a, and MdCHK4b. Furthermore, MdCHK3a
and MdCHK3b did not heterodimerize with each other, but
exclusively heterodimerize with MdCHK2, MdCHK4a, and
MdCHK4b. This particularity needs to be deeply addressed
for the complete understanding of MdCHK3 functioning.
To date, our overview clearly emphasized the complexity of
cytokinin perception in *M. domestica* since MdCHKs are able
to form not only homodimers, but also heterodimers as well
as monomers. Regarding the heterodimers in plant cells, they
provided a new layer of intricacy since most of the histidine
kinases form homodimers in order to autophosphorylate (Capra
and Laub, 2012). Nevertheless, their physiological relevance
*in planta* needs to be further determined. It cannot be
excluded that CHK heterodimers might operate in cytokinin
perception contributing to specify the cytokinin signaling
pathways in order to regulate distinct physiological processes.
For instance, the low gene expression of the MdCHK3
and MdCHK4 pairs compared to MdCHK2 suggested a
MdCHK2 dimerization ratio in favor of homodimerization. Such
homodimers could ensure signaling for the major physiological
processes associated to cytokinins while MdCHK2 heterodimers
would be associated to more discrete functions. In this
way, the putative functions of MdCHK3 predicted through
gene correlation analysis (Figure 6) would be assumed by
heterodimers with MdCHK2 or MdCHK4 since the MdCHK3 pair does not homo- or heterodimerize and is potentially not able to active the phosphorelay.

## CONCLUSION

This work provided a framework for further functional studies of cytokinin receptors in apple tree. In particular, it will be greatly interesting to focus on their involvement in response to the pathogens of apple tree. Furthermore, a structural approach would also contribute to gain insights into the key aspects of the mechanisms by which MdCHKs are differentially activated by cytokinin signal.

## AUTHOR CONTRIBUTIONS

DD, EA, CM, and GG conducted experiments. FL and TDdB achieved bioinformatics analyses. NP, VC, AO, AL, MCI, OP, and SB participated in the design of the study and interpretation. DG, SC, NG-G, MCo, JC, and SB assisted in the supervision of this work. GG conceived, supervised and coordinated the work. DD and GG wrote the first draft of the manuscript, to which all authors contributed.

## REFERENCES

Anantharaman, V., and Aravind, L. (2001). The CHASE domain: a predicted ligand-binding module in plant cytokinin receptors and other eukaryotic and bacterial receptors. *Trends Biochem. Sci.* 26, 579–582. doi: 10.1016/S0968-0004(01)01968-5

Argueso, C. T., Ferreira, F. J., Epple, P., To, J. P. C., Hutchison, C. E., Schaller, G. E., et al. (2012). Two-component elements mediate interactions between cytokinin and salicylic acid in plant immunity. *PLOS Genet.* 8:e1002448. doi: 10.1371/journal.pgen.1002448

Besseau, S., Kellner, F., Lanoue, A., Thamm, A. M. K., Salim, V., Schneider, B., et al. (2013). A pair of tabersonine 16-hydroxylases initiates the synthesis of vindoline in an organ-dependent manner in *Catharanthus roseus*. *Plant Physiol.* 163, 1792–1803. doi: 10.1104/pp.113.222828

Boivin, S., Kazmierczak, T., Brault, M., Wen, J., Gamas, P., Mysore, K. S., et al. (2016). Different cytokinin histidine kinase receptors regulate nodule initiation as well as later nodule developmental stages in *Medicago truncatula*. *Plant Cell Environ.* 39, 2198–2209. doi: 10.1111/pce. 12779

Caesar, K., Thamm, A. M. K., Witthoft, J., Elgass, K., Huppenberger, P., Grefen, C., et al. (2011). Evidence for the localization of the *Arabidopsis* cytokinin receptors AHK3 and AHK4 in the endoplasmic reticulum. *J. Exp. Bot.* 62, 5571–5580. doi: 10.1093/jxb/err238

Capra, E. J., and Laub, M. T. (2012). Evolution of two-component signal transduction systems. *Annu. Rev. Microbiol.* 66, 325–347. doi: 10.1146/ annurev-micro-092611-150039

Cortleven, A., Nitschke, S., Klaumunzer, M., AbdElgawad, H., Asard, H., Grimm, B., et al. (2014). A novel protective function for cytokinin in the light stress response is mediated by the *Arabidopsis* histidine kinase2 and *Arabidopsis* histidine kinase3 receptors. *Plant Physiol.* 164, 1470–1483. doi: 10.1104/pp.113. 224667

Cutcliffe, J. W., Hellmann, E., Heyl, A., and Rashotte, A. M. (2011). CRFs form protein-protein interactions with each other and with members of the cytokinin signalling pathway in *Arabidopsis* via the CRF domain. *J. Exp. Bot.* 62, 4995–5002. doi: 10.1093/jxb/err199

## FUNDING

This study was supported by the Région Centre-Val de Loire, France (SiSCyLi grant). Doctoral Fellow attributed to DD was jointly funded by the Région Centre-Val de Loire, France and the Ministère de l'Enseignement Supérieur et de la Recherche, France.

## ACKNOWLEDGMENTS

We thank Marie-Antoinette Marquet, Evelyne Danos, and Emeline Marais (EA2106 Biomolécules et Biotechnologies Végétales) for their help in maintaining cell cultures and Emilien Foureau for his technical help in subcellular localization. We also thank François Héricourt (LBLGC) for discussions concerning Histidine-Kinase receptors. We would also like to acknowledge the Fédération CaSciModOT (CCSC, Orléans, France) for accessing the Région Centre-Val de Loire computing grid.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fpls.2017.01614/ full#supplementary-material

Dortay, H., Mehnert, N., Bürkle, L., Schmülling, T., and Heyl, A. (2006). Analysis of protein interactions within the cytokinin-signaling pathway of *Arabidopsis thaliana*. *FEBS J.* 273, 4631–4644. doi: 10.1111/j.1742-4658.2006.05467.x

Foureau, E., Carqueijeiro, I., de Bernonville, T. D., Melin, C., Lafontaine, F., Besseau, S., et al. (2016). Prequels to synthetic biology: from candidate gene identification and validation to enzyme subcellular localization in plant and yeast cells. *Methods Enzymol.* 516, 167–206. doi: 10.1016/bs.mie.2016.02.013

Frugier, F., Kosuta, S., Murray, J. D., Crespi, M., and Szczyglowski, K. (2008). Cytokinin: secret agent of symbiosis. *Trends Plant Sci.* 13, 115–120. doi: 10.1016/j.tplants.2008.01.003

Ginis, O., Oudin, A., Guirimand, G., Chebbi, M., Courdavault, V., Glévarec, G., et al. (2012). A type-B response regulator drives the expression of the hydroxymethylbutenyl diphosphate synthase gene in periwinkle. *J. Plant Physiol.* 169, 1571–1574. doi: 10.1016/j.jplph.2012.07.008

Giron, D., and Glévarec, G. (2014). Cytokinin-induced phenotypes in plant-insect interactions: learning from the bacterial world. *J. Chem. Ecol.* 40, 826–835. doi: 10.1007/s10886-014-0466-5

Giron, D., Kaiser, W., Imbault, N., and Casas, J. (2007). Cytokinin-mediated leaf manipulation by a leafminer caterpillar. *Biol. Lett.* 3, 340–343. doi: 10.1098/rsbl. 2007.0051

Guirimand, G., Burlat, V., Oudin, A., Lanoue, A., St-Pierre, B., and Courdavault, V. (2009). Optimization of the transient transformation of *Catharanthus roseus* cells by particle bombardment and its application to the subcellular localization of hydroxymethylbutenyl 4-diphosphate synthase and geraniol 10-hydroxylase. *Plant Cell Rep.* 28, 1215–1234. doi: 10.1007/s00299-009- 0722-2

Guirimand, G., Courdavault, V., Lanoue, A., Mahroug, S., Guihur, A., Blanc, N., et al. (2010). Strictosidine activation in Apocynaceae: towards a 'nuclear time bomb'? *BMC Plant Biol.* 10:182. doi: 10.1186/1471-2229-10-182

Held, M., Hou, H., Miri, M., Huynh, C., Ross, L., Hossain, M. S., et al. (2014). *Lotus japonicus* cytokinin receptors work partially redundantly to mediate nodule formation. *Plant Cell* 26, 678–694. doi: 10.1105/tpc.113.119362

Heyl, A., Riefler, M., Romanov, G. A., and Schmülling, T. (2012). Properties, functions and evolution of cytokinin receptors. *Eur. J. Cell Biol.* 91, 246–256. doi: 10.1016/j.ejcb.2011.02.009
Heyl, A., Wulfetange, K., Pils, B., Nielsen, N., Romanov, G. A., and Schmülling, T. (2007). Evolutionary proteomics identifies amino acids essential for ligand-binding of the cytokinin receptor CHASE domain. *BMC Evol. Biol.* **7**:62. doi: **10.1186/1471-2148-7-62**

Hofmann, K., and Stoffel, W. (1993). TM base-A database of membrane spanning proteins segments. *Biol. Chem. Hoppe Seyler* **374**, 166.

Hothorn, M., Dabi, T., and Chory, J. (2011). Structural basis for cytokinin recognition by *Arabidopsis thaliana* histidine kinase 4. *Nat. Chem. Biol.* **7**, 766–768. doi: **10.1038/nchembio.667**

Inoue, T., Higuchi, M., Hashimoto, Y., Seki, M., Kobayashi, M., Kato, T., et al. (2001). Identification of CRE1 as a cytokinin receptor from *Arabidopsis*. *Nature* **409**, 1060–1063. doi: **10.1038/35059117**

Jeon, J., Cho, C., Lee, M. R., Van Binh, N., and Kim, J. (2016). *CYTOKININ RESPONSE FACTOR 2 (CRF2) and CRF3* regulate lateral root development in response to cold stress in *Arabidopsis*. *Plant Cell* **28**, 1828–1843. doi: **10.1105/tpc.15.00909**

Jung, S., Ficklin, S. P., Lee, T., Cheng, C. H., Blenda, A., Zheng, P., et al. (2014). The genome database for rosaceae (GDR): year 10 update. *Nucleic Acids Res.* **42**, 1237–1244. doi: **10.1093/nar/gkt1012**

Kaiser, W., Huguet, E., Casas, J., Commin, C., and Giron, D. (2010). Plant green-island phenotype induced by leaf-miners is mediated by bacterial symbionts. *Proc. Biol. Sci.* **277**, 2311–2319. doi: **10.1098/rspb.2010.0214**

Kieber, J. J., and Schaller, G. E. (2014). Cytokinins. *Arabidopsis Book* **12**:e0168. doi: **10.1199/tab.0168**

Kim, H. J., Ryu, H., Hong, S. H., Woo, H. R., Lim, P. O., Lee, I. C., et al. (2006). Cytokinin-mediated control of leaf longevity by AHK3 through phosphorylation of ARR2 in *Arabidopsis*. *Proc. Natl. Acad. Sci. U.S.A.* **103**, 814–819. doi: **10.1073/pnas.0505150103**

Krogh, A., Larsson, B., von Heijne, G., and Sonnhammer, E. L. (2001). Predicting transmembrane protein topology with a hidden markov model: application to complete genomes. *J. Mol. Biol.* **305**, 567–580. doi: **10.1006/jmb.2000.4315**

Kumar, M. N., and Verslues, P. E. (2015). Stress physiology functions of the *Arabidopsis* histidine kinase cytokinin receptors. *Physiol. Plant*. **154**, 369–380. doi: **10.1111/ppl.12290**

Kuderová, A., Gallová, L., Kuricová, K., Nejedlá, E., Čurdová, A., Micenková, L., et al. (2015). Identification of AHK2- and AHK3-like cytokinin receptors in *Brassica napus* reveals two subfamilies of AHK2 orthologues. *J. Exp. Bot.* **66**, 339–353. doi: **10.1093/jxb/eru422**

Letunic, I., Doerks, T., and Bork, P. (2015). SMART: recent updates, new developments and status in 2015. *Nucleic Acids Res.* **43**, 257–260. doi: **10.1093/nar/gku949**

Lomin, S. N., Krivosheev, D. M., Steklov, M. Y., Arkhipov, D. V., Osolodkin, D. I., Schmülling, T., et al. (2015). Plant membrane assays with cytokinin receptors underpin the unique role of free cytokinin bases as biologically active ligands. *J. Exp. Bot.* **66**, 1851–1863. doi: **10.1093/jxb/eru522**

Lomin, S. N., Krivosheev, D. M., Steklov, M. Y., Osolodkin, D. I., and Romanov, G. A. (2012). Receptor properties and features of cytokinin signaling. *Acta Naturae* **43**, 31–45.

Lomin, S. N., Yonekura-Sakakibara, K., Romanov, G. A., and Sakakibara, H. (2011). Ligand-binding properties and subcellular localization of maize cytokinin receptors. *J. Exp. Bot.* **62**, 5149–5159. doi: **10.1093/jxb/err220**

Maeda, T., Wurgler-Murphy, S. M., and Saito, H. (1994). A two-component system that regulates an osmosensing MAP kinase cascade in yeast. *Nature* **369**, 242–245. doi: **10.1038/369242a0**

Mason, M. G., Mathews, D. E., Argyros, D. A., Maxwell, B. B., Kieber, J. J., Alonso, J. M., et al. (2005). Multiple type-B response regulators mediate cytokinin signal transduction in *Arabidopsis*. *Plant Cell* **17**, 3007–3018. doi: **10.1105/tpc.105.035451**

McKenzie, D. J., McLean, M. A., Mukerji, S., and Green, M. (1997). Improved RNA extraction from woody plants for the detection of viral pathogens by reverse transcription-polymerase chain reaction. *Plant Dis.* **81**, 222–226. doi: **10.1094/PDIS.1997.81.2.222**

Miwa, K., Ishikawa, K., Terada, K., Yamada, H., Suzuki, T., Yamashino, T., et al. (2007). Identification of amino acid substitutions that render the *Arabidopsis* cytokinin receptor histidine kinase AHK4 constitutively active. *Plant Cell Physiol.* **48**, 1809–1814. doi: **10.1093/pcp/pcm145**

Mougel, C., and Zhulin, I. B. (2001). CHASE: an extracellular sensing domain common to transmembrane receptors from prokaryotes, lower eukaryotes and plants. *Trends Biochem. Sci.* **26**, 582–584. doi: **10.1016/S0968-0004(01)01969-7**

Naseem, M., Wölfling, M., and Dandekar, T. (2014). Cytokinins for immunity beyond growth, galls and green islands. *Trends Plant Sci.* **19**, 481–484. doi: **10.1016/j.tplants.2014.04.001**

Nieminen, K., Immanen, J., Laxell, M., Kauppinen, L., Tarkowski, P., Dolezal, K., et al. (2008). Cytokinin signaling regulates cambial development in poplar. *Proc. Natl. Acad. Sci. U.S.A.* **105**, 20032–20037. doi: **10.1073/pnas.0805617106**

Osugi, A., and Sakakibara, H. (2015). Q&A: How do plants respond to cytokinins and what is their importance? *BMC Biol.* **13**:102. doi: **10.1186/s12915-015-0214-5**

Papadopoulos, J. S., and Agarwala, R. (2007). COBALT: constraint-based alignment tool for multiple protein sequences. *Bioinformatics* **23**, 1073–1079. doi: **10.1093/bioinformatics/btm076**

Pertry, I., Václavíková, K., Depuydt, S., Galuszka, P., Spíchal, L., Temmerman, W., et al. (2009). Identification of *Rhodococcus fascians* cytokinins and their modus operandi to reshape the plant. *Proc. Natl. Acad. Sci. U.S.A.* **106**, 929–934. doi: **10.1073/pnas.0811683106**

Raines, T., Shanks, C., Cheng, C. Y., McPherson, D., Argueso, C. T., Kim, H. J., et al. (2016). The cytokinin response factors modulate root and shoot growth and promote leaf senescence in *Arabidopsis*. *Plant J.* **85**, 134–147. doi: **10.1111/tpj.13097**

Rambaldi, D., and Ciccarelli, F. D. (2009). FancyGene: dynamic visualization of gene structures and protein domain architectures on genomic loci. *Bioinformatics* **25**, 2281–2282. doi: **10.1093/bioinformatics/btp381**

Riefler, M., Novak, O., Strnad, M., and Schmülling, T. (2006). *Arabidopsis* cytokinin receptor mutants reveal functions in shoot growth, leaf senescence, seed size, germination, root development, and cytokinin metabolism. *Plant Cell* **18**, 40–54. doi: **10.1105/tpc.105.037796**

Romanov, G. A., Lomin, S. N., and Schmülling, T. (2006). Biochemical characteristics and ligand-binding properties of *Arabidopsis* cytokinin receptor AHK3 compared to CRE1/AHK4 as revealed by a direct binding assay. *J. Exp. Bot.* **57**, 4051–4058. doi: **10.1093/jxb/erl179**

Sakakibara, H. (2006). Cytokinins: activity, biosynthesis, and translocation. *Annu. Rev. Plant Biol.* **57**, 431–449. doi: **10.1146/annurev.arplant.57.032905.105231**

Schäfer, M., Meza-Canales, I. D., Brütting, C., Baldwin, I. T., and Meldau, S. (2015). Cytokinin concentrations and CHASE-domain containing His Kinase 2 (NaCHK2)- and NaCHK3-mediated perception modulate herbivory-induced defense signaling and defenses in *Nicotiana attenuata*. *New Phytol.* **207**, 645–658. doi: **10.1111/nph.13404**

Sigrist, C. J., Cerutti, L., Hulo, N., Gattiker, A., Falquet, L., Pagni, M., et al. (2002). PROSITE: a documented database using patterns and profiles as motif descriptors. *Brief. Bioinform.* **3**, 265–274. doi: **10.1093/bib/3.3.265**

Snapp, E. L., Hegde, R. S., Francolini, M., Lombardo, F., Colombo, S., Pedrazzini, E., et al. (2003). Formation of stacked ER cisternae by low affinity protein interactions. *J. Cell Biol.* **163**, 257–269. doi: **10.1083/jcb.200306020**

Spíchal, L. (2012). Cytokinins - recent news and views of evolutionarily old molecules. *Funct. Plant Biol.* **39**, 267–284. doi: **10.1071/FP11276**

Spíchal, L., Rakova, N. Y., Riefler, M., Mizuno, T., Romanov, G. A., Strnad, M., et al. (2004). Two cytokinin receptors of *Arabidopsis thaliana*, CRE1/AHK4 and AHK3, differ in their ligand specificity in a bacterial assay. *Plant Cell Physiol.* **45**, 1299–1305. doi: **10.1093/pcp/pch132**

Steklov, M. Y., Lomin, S. N., Osolodkin, D. I., and Romanov, G. A. (2013). Structural basis for cytokinin receptor signaling: an evolutionary approach. *Plant Cell Rep.* **32**, 781–793. doi: **10.1007/s00299-013-1408-3**

Stolz, A., Riefler, M., Lomin, S. N., Achazi, K., Romanov, G. A., and Schmülling, T. (2011). The specificity of cytokinin signalling in *Arabidopsis thaliana* is mediated by differing ligand affinities and expression profiles of the receptors. *Plant J.* **67**, 157–168. doi: **10.1111/j.1365-313X.2011.04584.x**

Tirichine, L., Sandal, N., Madsen, L. H., Radutoiu, S., Albrektsen, A. S., Sato, S., et al. (2007). A gain-of-function mutation in a cytokinin receptor triggers spontaneous root nodule organogenesis. *Science* **315**, 104–107. doi: **10.1126/science.1132397**

To, J. P., Haberer, G., Ferreira, F. J., Deruère, J., Mason, M. G., Schaller, G. E., et al. (2004). Type-A Arabidopsis response regulators are partially redundant negative regulators of cytokinin signaling. *Plant Cell* **16**, 658–671. doi: **10.1105/tpc.018978**
Tran, L.-S. P., Urao, T., Qin, F., Maruyama, K., Kakimoto, T., Shinozaki, K., et al. (2007). Functional analysis of AHK1/ATHK1 and cytokinin receptor histidine kinases in response to abscisic acid, drought, and salt stress in *Arabidopsis*. *Proc. Natl. Acad. Sci. U.S.A.* 104, 20623–20628. doi: 10.1073/pnas.0706547105

Ueguchi, C., Koizumi, H., Suzuki, T., and Mizuno, T. (2001). Novel family of sensor histidine kinase genes in *Arabidopsis thaliana*. *Plant Cell Physiol.* 42, 231–235. doi: 10.1093/pcp/pce015

Van Bel, M., Proost, S., Wischnitzki, E., Movahedi, S., Scheerlinck, C., Van de Peer, Y., et al. (2012). Dissecting plant genomes with the PLAZA comparative genomics platform. *Plant Physiol* 158, 590–600. doi: 10.1104/pp.111.189514

Velasco, R., Zharkikh, A., Affourtit, J., Dhingra, A., Cestaro, A., Kalyanaraman, A., et al. (2010). The genome of the domesticated apple (Malus × domestica Borkh.). *Nat. Genet.* 42, 833–839. doi: 10.1038/ng.654

von Schwartzenberg, K., Lindner, A. C., Gruhn, N., Šimura, J., Novák, O., Strnad, M., et al. (2016). CHASE domain-containing receptors play an essential role in the cytokinin response of the moss *Physcomitrella patens*. *J. Exp. Bot.* 67, 667–679. doi: 10.1093/jxb/erv479

Waadt, R., and Kudla, J. (2008). In planta visualization of protein interactions using bimolecular fluorescence complementation (BiFC). *CSH Protocol*. 2008:db.rot49953. doi: 10.1101/pdb.prot4995

Wulfetange, K., Lomin, S. N., Romanov, G. A., Stolz, A., Heyl, A., and Schmulling, T. (2011). The cytokinin receptors of *Arabidopsis* are located mainly to the endoplasmic reticulum. *Plant Physiol.* 156, 1808–1818. doi: 10.1104/pp.111.180539

Yonekura-Sakakibara, K., Kojima, M., Yamaya, T., and Sakakibara, H. (2004). Molecular characterization of cytokinin-responsive histidine kinases in maize.

Differential ligand preferences and response to cis-zeatin. *Plant Physiol.* 134, 1654–1661. doi: 10.1104/pp.103.037176

Zhang, H., de Bernonville, T. D., Body, M., Glevarec, G., Reichelt, M., Unsicker, S., et al. (2016). Leaf-mining by Phyllonorycter blancardella reprograms the host-leaf transcriptome to modulate phytohormones associated with nutrient mobilization and plant defense. *J. Insect Physiol.* 84, 114–127. doi: 10.1016/j.jinsphys.2015.06.003

Zürcher, E., and Müller, B. (2016). Cytokinin synthesis, signaling, and function—advances and new insights. *Int. Rev. Cell Mol. Biol.* 324, 1–38. doi: 10.1016/bs.ircmb.2016.01.001

Zwack, P. J., and Rashotte, A. M. (2015). Interactions between cytokinin signalling and abiotic stress responses. *J. Exp. Bot.* 66, 4863–4871. doi: 10.1093/jxb/erv172

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

Copyright © 2017 Daudu, Allion, Liesecke, Papon, Courdavault, Dugé de Bernonville, Mélin, Oudin, Clastre, Lanoue, Courtois, Pichon, Giron, Carpin, Giglioli-Guivarc'h, Crèche, Besseau and Glévarec. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.