CHASE domain-containing receptors play an essential role in the cytokinin response of the moss Physcomitrella patens

Cytokinin and ALOG proteins regulate pluripotent stem cell identity in the moss Physcomitrium patens

The shoot apical meristem (SAM) contains pluripotent stem cells that produce all the aerial parts of the plant. Stem cells undergo asymmetric cell divisions to self-renew and to produce differentiating cells. Our research focused on unraveling the mechanisms governing the specification of these two distinct cell fates following the stem cell division. For this purpose, we used the model organism Physcomitrium patens, which features a singular pluripotent stem cell known as the gametophore apical cell. We show that the activity of cytokinins, critical stem cell regulators, is restricted to the gametophore apical cell due to the specific localization of PpLOG, the enzyme responsible for cytokinin activation. In turn, PpTAW, which promotes differentiating cell identity of the merophyte, is excluded from the gametophore apical cell by the action of cytokinins. We propose a cytokinin-based model for the establishment of asymmetry in the pluripotent stem cell division.

Cytokinin: From autoclaved DNA to two-component signaling

Since its first identification in the 1950s as a regulator of cell division, cytokinin has been linked to many physiological processes in plants, spanning growth and development and various responses to the environment. Studies from the last two and one-half decades have revealed the pathways underlying the biosynthesis and metabolism of cytokinin and have elucidated the mechanisms of its perception and signaling, which reflects an ancient signaling system evolved from two-component elements in bacteria. Mutants in the genes encoding elements involved in these processes have helped refine our understanding of cytokinin functions in plants. Further, recent advances have provided insight into the mechanisms of intracellular and long-distance cytokinin transport and the identification of several proteins that operate downstream of cytokinin signaling. Here, we review these processes through a historical lens, providing an overview of cytokinin metabolism, transport, signaling, and functions in higher plants.

The cytokinin receptor OHK4/OsHK4 regulates inflorescence architecture in rice via an IDEAL PLANT ARCHITECTURE1/WEALTHY FARMER'S PANICLE-mediated positive feedback circuit

Inflorescence architecture is important for rice (Oryza sativa) grain yield. The phytohormone cytokinin (CK) has been shown to regulate rice inflorescence development; however, the underlying mechanism mediated by CK perception is still unclear. Employing a forward genetic approach, we isolated an inactive variant of the CK receptor OHK4/OsHK4 gene named panicle length1, which shows decreased panicle size due to reduced inflorescence meristem (IM) activity. A 2-amino acid deletion in the long α-helix stalk of the sensory module of OHK4 impairs the homodimerization and ligand-binding capacity of the receptor, even though the residues do not touch the ligand-binding domain or the dimerization interface. This deletion impairs CK signaling that occurs through the type-B response regulator OsRR21, which acts downstream of OHK4 in controlling inflorescence size. Meanwhile, we found that IDEAL PLANT ARCHITECTURE1(IPA1)/WEALTHY FARMER'S PANICLE (WFP), encoding a positive regulator of IM development, acts downstream of CK signaling and is directly activated by OsRR21. Additionally, we revealed that IPA1/WFP directly binds to the OHK4 promoter and upregulates its expression through interactions with 2 TCP transcription factors, forming a positive feedback circuit. Altogether, we identified the OHK4-OsRR21-IPA1 regulatory module, providing important insights into the role of CK signaling in regulating rice inflorescence architecture.

The ratio of auxin to cytokinin controls leaf development and meristem initiation in Physcomitrium patens

Crosstalk between auxin and cytokinin contributes to widespread developmental processes, including root and shoot meristem maintenance, phyllotaxy, and vascular patterning. However, our understanding of crosstalk between these hormones is limited primarily to angiosperms. The moss Physcomitrium patens (formerly Physcomitrella patens) is a powerful system for studying plant hormone function. Auxin and cytokinin play similar roles in regulating moss gametophore (shoot) architecture, to those in flowering plant shoots. However, auxin-cytokinin crosstalk is poorly understood in moss. Here we find that the ratio of auxin to cytokinin is an important determinant of development in P. patens, especially during leaf development and branch stem cell initiation. Addition of high levels of auxin to P. patens gametophores blocks leaf outgrowth. However, simultaneous addition of high levels of both auxin and cytokinin partially restores leaf outgrowth, suggesting that the ratio of these hormones is the predominant factor. Likewise, during branch initiation and outgrowth, chemical inhibition of auxin synthesis phenocopies cytokinin application. Finally, cytokinin-insensitive mutants resemble plants with altered auxin signaling and are hypersensitive to auxin. In summary, our results suggest that the ratio between auxin and cytokinin signaling is the basis for developmental decisions in the moss gametophore.

The origin and early evolution of cytokinin signaling

Angiosperms, especially Arabidopsis and rice, have long been at the center of plant research. However, technological advances in sequencing have led to a dramatic increase in genome and transcriptome data availability across land plants and, more recently, among green algae. These data allowed for an in-depth study of the evolution of different protein families - including those involved in the metabolism and signaling of phytohormones. While most early studies on phytohormone evolution were phylogenetic, those studies have started to be complemented by genetic and biochemical studies in recent years. Examples of such functional analyses focused on ethylene, jasmonic acid, abscisic acid, and auxin. These data have been summarized recently. In this review, we will focus on the progress in our understanding of cytokinin biology. We will use these data to synthesize key points about the evolution of cytokinin metabolism and signaling, which might apply to the evolution of other phytohormones as well.

Molecular Regulation and Evolution of Cytokinin Signaling in Plant Abiotic Stresses

The sustainable production of crops faces increasing challenges from global climate change and human activities, which leads to increasing instances of many abiotic stressors to plants. Among the abiotic stressors, drought, salinity and excessive levels of toxic metals cause reductions in global agricultural productivity and serious health risks for humans. Cytokinins (CKs) are key phytohormones functioning in both normal development and stress responses in plants. Here, we summarize the molecular mechanisms on the biosynthesis, metabolism, transport and signaling transduction pathways of CKs. CKs act as negative regulators of both root system architecture plasticity and root sodium exclusion in response to salt stress. The functions of CKs in mineral-toxicity tolerance and their detoxification in plants are reviewed. Comparative genomic analyses were performed to trace the origin, evolution and diversification of the critical regulatory networks linking CK signaling and abiotic stress. We found that the production of CKs and their derivatives, pathways of signal transduction and drought-response root growth regulation are evolutionarily conserved in land plants. In addition, the mechanisms of CK-mediated sodium exclusion under salt stress are suggested for further investigations. In summary, we propose that the manipulation of CK levels and their signaling pathways is important for plant abiotic stress and is, therefore, a potential strategy for meeting the increasing demand for global food production under changing climatic conditions.

How was apical growth regulated in the ancestral land plant? Insights from the development of non-seed plants

Land plant life cycles are separated into distinct haploid gametophyte and diploid sporophyte stages. Indeterminate apical growth evolved independently in bryophyte (moss, liverwort, and hornwort) and fern gametophytes, and tracheophyte (vascular plant) sporophytes. The extent to which apical growth in tracheophytes co-opted conserved gametophytic gene networks, or exploited ancestral sporophytic networks, is a long-standing question in plant evolution. The recent phylogenetic confirmation of bryophytes and tracheophytes as sister groups has led to a reassessment of the nature of the ancestral land plant. Here, we review developmental genetic studies of apical regulators and speculate on their likely evolutionary history.

Exogenous 6-benzylaminopurine inhibits tip growth and cytokinesis via regulating actin dynamics in the moss Physcomitrium patens

Exogenous BAP but not 2iP disrupts actin structures and induces tip-growth retardation and cytokinesis failure in the moss Physcomitrium patens. Synthetic cytokinins have been widely used to address hormonal responses during plant development. However, exogenous cytokinins can cause a variety of cellular effects. A detailed characterization of such effects has not been well studied. Here, using Physcomitrium patens as a model, we show that the aromatic cytokinin 6-benzylaminopurine (BAP) inhibits tip growth at concentrations above 0.2 µM. At higher concentrations (0.6-1 µM), BAP can additionally block mitotic entry and induce cytokinesis defects and cell death. These effects are associated with altered actin dynamics and structures. By contrast, 2-isopentenyladenine (2iP) does not cause marked defects at various concentrations up to 10 µM, while t-zeatin (tZ) can moderately inhibit moss growth. Our results provide mechanistic insight into the inhibitory effects of BAP on cell growth and cell division and call for attention to the use of synthetic cytokinins for bioassays.

Division site determination during asymmetric cell division in plants

During development, both animals and plants exploit asymmetric cell division (ACD) to increase tissue complexity, a process that usually generates cells dissimilar in size, morphology, and fate. Plants lack the key regulators that control ACD in animals. Instead, plants have evolved two unique cytoskeletal structures to tackle this problem: the preprophase band (PPB) and phragmoplast. The assembly of the PPB and phragmoplast and their contributions to division plane orientation have been extensively studied. However, how the division plane is positioned off the cell center during asymmetric division is poorly understood. Over the past 20 years, emerging evidence points to a critical role for polarly localized membrane proteins in this process. Although many of these proteins are species- or cell type specific, and the molecular mechanism underlying division asymmetry is not fully understood, common features such as morphological changes in cells, cytoskeletal dynamics, and nuclear positioning have been observed. In this review, we provide updates on polarity establishment and nuclear positioning during ACD in plants. Together with previous findings about symmetrically dividing cells and the emerging roles of developmental cues, we aim to offer evolutionary insight into a common framework for asymmetric division-site determination and highlight directions for future work.

Cytokinin-CLAVATA cross-talk is an ancient mechanism regulating shoot meristem homeostasis in land plants

SignificancePlants grow from their tips. The gametophore (shoot-like organ) tip of the moss Physcomitrium patens is a single cell that performs the same functions as those of multicellular flowering plants, producing the cells that make leaves and regenerating new stem cells to maintain the shoot tip. Several pathways, including CLAVATA and cytokinin hormonal signaling, regulate stem cell abundance in flowering plants and in mosses, although the mechanisms whereby these pathways regulate stem cell abundance and their conservation between these plant lineages is poorly understood. Using moss, we investigated how PpCLAVATA and cytokinin signaling interact. Overall, we found evidence that PpCLAVATA and cytokinin signaling interact similarly in moss and flowering plants, despite their distinct anatomies, life cycles, and evolutionary distance.

Red Light and Glucose Enhance Cytokinin-Mediated Bud Initial Formation in Physcomitrium patens

Growth and development of Physcomitrium patens is endogenously regulated by phytohormones such as auxin and cytokinin. Auxin induces the transition of chloronema to caulonema. This transition is also regulated by additional factors such as quantity and quality of light, carbon supply, and other phytohormones such as strigolactones and precursors of gibberrelic acid. On the other hand, cytokinins induce the formation of bud initials following caulonema differentiation. However, the influence of external factors such as light or nutrient supply on cytokinin-mediated bud initial formation has not been demonstrated in Physcomitrium patens. This study deals with the effect of light quality and nutrient supply on cytokinin-mediated bud initial formation. Bud initial formation has been observed in wild type plants in different light conditions such as white, red, and blue light in response to exogenously supplied cytokinin as well as glucose. In addition, budding assay has been demonstrated in the cry1a mutant of Physcomitrium in different light conditions. The results indicate that carbon supply and red light enhance the cytokinin response, while blue light inhibits this process in Physcomitrium.

Cytokinin Perception in Ancient Plants beyond Angiospermae

Cytokinins (CKs) control many plant developmental processes and responses to environmental cues. Although the CK signaling is well understood, we are only beginning to decipher its evolution. Here, we investigated the CK perception apparatus in early-divergent plant species such as bryophyte Physcomitrium patens, lycophyte Selaginella moellendorffii, and gymnosperm Picea abies. Of the eight CHASE-domain containing histidine kinases (CHKs) examined, two CHKs, PpCHK3 and PpCHK4, did not bind CKs. All other CHK receptors showed high-affinity CK binding (K_D of nM range), with a strong preference for isopentenyladenine over other CK nucleobases in the moss and for trans-zeatin over cis-zeatin in the gymnosperm. The pH dependences of CK binding for these six CHKs showed a wide range, which may indicate different subcellular localization of these receptors at either the plasma- or endoplasmic reticulum membrane. Thus, the properties of the whole CK perception apparatuses in early-divergent lineages were demonstrated. Data show that during land plant evolution there was a diversification of the ligand specificity of various CHKs, in particular, the rise in preference for trans-zeatin over cis-zeatin, which indicates a steadily increasing specialization of receptors to various CKs. Finally, this distinct preference of individual receptors to different CK versions culminated in vascular plants, especially angiosperms.

Fundamental mechanisms of the stem cell regulation in land plants: lesson from shoot apical cells in bryophytes

This review compares the molecular mechanisms of stem cell control in the shoot apical meristems of mosses and angiosperms and reveals the conserved features and evolution of plant stem cells. The establishment and maintenance of pluripotent stem cells in the shoot apical meristem (SAM) are key developmental processes in land plants including the most basal, bryophytes. Bryophytes, such as Physcomitrium (Physcomitrella) patens and Marchantia polymorpha, are emerging as attractive model species to study the conserved features and evolutionary processes in the mechanisms controlling stem cells. Recent studies using these model bryophyte species have started to uncover the similarities and differences in stem cell regulation between bryophytes and angiosperms. In this review, we summarize findings on stem cell function and its regulation focusing on different aspects including hormonal, genetic, and epigenetic control. Stem cell regulation through auxin, cytokinin, CLAVATA3/EMBRYO SURROUNDING REGION-RELATED (CLE) signaling and chromatin modification by Polycomb Repressive Complex 2 (PRC2) and PRC1 is well conserved. Several transcription factors crucial for SAM regulation in angiosperms are not involved in the regulation of the SAM in mosses, but similarities also exist. These findings provide insights into the evolutionary trajectory of the SAM and the fundamental mechanisms involved in stem cell regulation that are conserved across land plants.

Phytohormone biosynthesis and signaling pathways of mosses

Most known phytohormones regulate moss development. We present a comprehensive view of the synthesis and signaling pathways for the most investigated of these compounds in mosses, focusing on the model Physcomitrium patens. The last 50 years of research have shown that most of the known phytohormones are synthesized by the model moss Physcomitrium patens (formerly Physcomitrella patens) and regulate its development, in interaction with responses to biotic and abiotic stresses. Biosynthesis and signaling pathways are best described in P. patens for the three classical hormones auxins, cytokinins and abscisic acid. Furthermore, their roles in almost all steps of development, from early filament growth to gametophore development and sexual reproduction, have been the focus of much research effort over the years. Evidence of hormonal roles exist for ethylene and for CLE signaling peptides, as well as for salicylic acid, although their possible effects on development remain unclear. Production of brassinosteroids by P. patens is still debated, and modes of action for these compounds are even less known. Gibberellin biosynthesis and signaling may have been lost in P. patens, while gibberellin precursors such as ent-kaurene derivatives could be used as signals in a yet to discover pathway. As for jasmonic acid, it is not used per se as a hormone in P. patens, but its precursor OPDA appears to play a corresponding role in defense against abiotic stress. We have tried to gather a comprehensive view of the biosynthesis and signaling pathways for all these compounds in mosses, without forgetting strigolactones, the last class of plant hormones to be reported. Study of the strigolactone response in P. patens points to a novel signaling compound, the KAI2-ligand, which was likely employed as a hormone prior to land plant emergence.

Red light-regulated interaction of Per-Arnt-Sim histidine kinases with partner histidine-containing phosphotransfer proteins in Physcomitrium patens

Multi-step phosphorelay (MSP) is a broadly distributed signaling system in organisms. In MSP, histidine kinases (HKs) receive various environmental signals and transmit them by autophosphorylation followed by phosphotransfer to partner histidine-containing phosphotransfer proteins (HPts). Previously, we reported that Per-Arnt-Sim (PAS) domain-containing HK1 (PHK1) and PHK2 of the moss Physcomitrium (Physcomitrella) patens repressed red light-induced protonema branching, a critical step in the moss life cycle. In plants, PHK homolog-encoding genes are conserved only in early-diverging lineages such as bryophytes and lycophytes. PHKs-mediated signaling machineries attract attention especially from an evolutionary viewpoint, but they remain uninvestigated. Here, we studied the P. patens PHKs focusing on their subcellular patterns of localization and interaction with HPts. Yeast two-hybrid analysis, a localization assay with a green fluorescent protein, and a bimolecular fluorescence complementation analysis together showed that PHKs are localized and interact with partner HPts mostly in the nucleus, as unprecedented features for plant HKs. Additionally, red light triggered the interactions between PHKs and HPts in the cytoplasm, and light co-repressed the expression of PHK1 and PHK2 as well as genes encoding their partner HPts. Our results emphasize the uniqueness of PHKs-mediated signaling machineries, and functional implications of this uniqueness are discussed.

There and back again: An evolutionary perspective on long-distance coordination of plant growth and development

Vascular plants, unlike bryophytes, have a strong root-shoot dichotomy in which the tissue systems are mutually interdependent; roots are completely dependent on shoots for photosynthetic sugars, and shoots are completely dependent on roots for water and mineral nutrients. Long-distance communication between shoot and root is therefore critical for the growth, development and survival of vascular plants, especially with regard to variable environmental conditions. However, this long-distance signalling does not appear an ancestral feature of land plants, and has likely arisen in vascular plants to service the radical alterations in body-plan seen in this taxon. In this review, we examine the defined hormonal root-to-shoot and shoot-to-root signalling pathways that coordinate the growth of vascular plants, with a particular view to understanding how these pathways may have evolved. We highlight the completely divergent roles of isopentenyl-adenine and trans-zeatin cytokinin species in long-distance signalling, and ask whether cytokinin can really be considered as a single class of hormones in the light of recent research. We also discuss the puzzlingly sparse evidence for auxin as a shoot-to-root signal, the evolutionary re-purposing of strigolactones and gibberellins as hormonal signals, and speculate on the possible role of sugars as long-distance signals. We conclude by discussing the 'design principles' of long-distance signalling in vascular plants.

The 2D to 3D growth transition in the moss Physcomitrella patens

The colonization of land by plants coincided with and was most likely facilitated by the evolution of 3-dimensional (3D) growth. 3D growth is a pivotal feature of all land plants, but most develop in a way that precludes genetic investigation. In the moss Physcomitrella patens, 3D growth (gametophores) is preceded by an extended 2-dimensional (2D) growth phase (protonemata) that can be propagated indefinitely. Studies using P. patens have thus elucidated some of the molecular mechanisms underlying 3D growth regulation. This review summarizes the known molecular mechanisms underlying both the formation of gametophore initial cells and the development of the 3D growth in gametophores.

Chiral events in developing gametophores of Physcomitrella patens and other moss species are driven by an unknown, universal direction-sensing mechanism

PREMISE OF THE STUDY

We used the model species Physcomitrella patens to examine chirality in moss gametophores. Chirality is manifested in the direction of consecutive apical cell divisions, cell plate configurations, and deviations of leaf connecting lines from the vertical course. However, the frequencies of chiral configurations of all these processes as well as their mutual dependence-especially in the case of gametophore branching-are not known. Other moss species were checked to determine the universality of our findings.

METHODS

The gametophore structure of Physcomitrella patens grown in the laboratory under controlled conditions was investigated using light microscopy and compared with that of other moss species collected from their natural stands.

KEY RESULTS

In all investigated moss species, the tetrahedral apical cell exhibits either clockwise or counterclockwise consecutive divisions, and selection of this directionality in the primary axis is random. It is, however, related to cell plate configuration. If the plate is skewed, leaf-producing segments arising from the apical cell cleavage exhibit circumferential rotation. Three parallel lines connecting the leaves deviate from a vertical course, but always in the same direction as that of leaf initiation; thus, the angular distance between consecutive leaves increases to >120°. Lateral branches are exclusively antidromous.

CONCLUSIONS

Gametophore chiral configuration appears to be useful in resolving problems of moss modular growth and branching. Morphological and anatomical evidence strongly suggests that an unknown direction-sensing mechanism controls the development of moss axial organs. We propose that leaves are responsible for a horizontal gradient of sugar signals that develops along the gametophore circumference, thus influencing branching-unit chirality.

RNA-Binding Protein RBP-P Is Required for Glutelin and Prolamine mRNA Localization in Rice Endosperm Cells

In developing rice (Oryza sativa) endosperm, mRNAs of the major storage proteins, glutelin and prolamine, are transported and anchored to distinct subdomains of the cortical endoplasmic reticulum. RNA binding protein RBP-P binds to both glutelin and prolamine mRNAs, suggesting a role in some aspect of their RNA metabolism. Here, we show that rice lines expressing mutant RBP-P mislocalize both glutelin and prolamine mRNAs. Different mutant RBP-P proteins exhibited varying degrees of reduced RNA binding and/or protein-protein interaction properties, which may account for the mislocalization of storage protein RNAs. In addition, partial loss of RBP-P function conferred a broad phenotypic variation ranging from dwarfism, chlorophyll deficiency, and sterility to late flowering and low spikelet fertility. Transcriptome analysis highlighted the essential role of RBP-P in regulating storage protein genes and several essential biological processes during grain development. Overall, our data demonstrate the significant roles of RBP-P in glutelin and prolamine mRNA localization and in the regulation of genes important for plant growth and development through its RNA binding activity and cooperative regulation with interacting proteins.

Light-regulated PAS-containing histidine kinases delay gametophore formation in the moss Physcomitrella patens

Two-component systems (TCSs) are signal transduction mechanisms for responding to various environmental stimuli. In angiosperms, TCSs involved in phytohormone signaling have been intensively studied, whereas there are only a few reports on TCSs in basal land plants. The moss Physcomitrella patens possesses several histidine kinases (HKs) that are lacking in seed plant genomes. Here, we studied two of these unique HKs, PAS-histidine kinase 1 (PHK1) and its paralog PHK2, both of which have PAS (Per-Arnt-Sim) domains, which are known to show versatile functions such as sensing light or molecular oxygen. We found homologs of PHK1 and PHK2 only in early diverged clades such as bryophytes and lycophytes, but not in seed plants. The PAS sequences of PHK1 and PHK2 are more similar to a subset of bacterial PAS sequences than to any angiosperm PAS sequences. Gene disruption lines that lack either PHK1 or PHK2 or both formed gametophores earlier than the wild-type, and consistently, more caulonema side branches were induced in response to light in the disruption lines. Therefore, PHK1 and PHK2 delay the timing of gametophore development, probably by suppressing light-induced caulonema branching. This study provides new insights into the evolution of TCSs in plants.

The effects of repeated whole genome duplication events on the evolution of cytokinin signaling pathway

BACKGROUND

It is thought that after whole-genome duplications (WGDs), a large fraction of the duplicated gene copies is lost over time while few duplicates are retained. Which factors promote survival or death of a duplicate remains unclear and the underlying mechanisms are poorly understood. According to the model of gene dosage balance, genes encoding interacting proteins are predicted to be preferentially co-retained after WGDs. Among these are genes encoding proteins involved in complexes or in signal transduction.

RESULTS

We have investigated the way that repeated WGDs during land plant evolution have affected cytokinin signaling to study patterns of gene duplicability and co-retention in this important signal transduction pathway. Through the integration of phylogenetic analyses with comparisons of genome collinearity, we have found that signal input mediated by cytokinin receptors proved to be highly conserved over long evolutionary time-scales, with receptors showing predominantly gene loss after repeated WGDs. However, the downstream elements, e,g. response regulators, were mainly retained after WGDs and thereby formed gene families in most plant lineages.

CONCLUSIONS

Gene dosage balance between the interacting components indicated by co-retention after WGDs seems to play a minor role in the evolution of cytokinin signaling pathway. Overall, core genes of cytokinin signaling show a highly heterogeneous pattern of gene retention after WGD, reflecting complex relationships between the various factors that shape the long-term fate of a duplicated gene.

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

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.