Genome-Wide Analysis of BBX Gene Family in Three Medicago Species Provides Insights into Expression Patterns under Hormonal and Salt Stresses.

BBX protein is a class of zinc finger transcription factors that have B-box domains at the N-terminus, and some of these proteins contain a CCT domain at the C-terminus. It plays an important role in plant growth, development, and metabolism. However, the expression pattern of BBX genes in alfalfa under hormonal and salt stresses is still unclear. In this study, we identified a total of 125 BBX gene family members by the available Medicago reference genome in diploid alfalfa (Medicago sativa spp. Caerulea), a model plant (M. truncatula), and tetraploid alfalfa (M. sativa), and divided these members into five subfamilies. We found that the conserved motifs of BBXs of the same subfamily reveal similarities. We analyzed the collinearity relationship and duplication mode of these BBX genes and found that the expression pattern of BBX genes is specific in different tissues. Analysis of the available transcriptome data suggests that some members of the BBX gene family are involved in multiple abiotic stress responses, and the highly expressed genes are often clustered together. Furthermore, we identified different expression patterns of some BBX genes under salt, ethylene, salt and ethylene, salicylic acid, and salt and salicylic acid treatments, verified by qRT-PCR, and analyzed the subcellular localization of MsBBX2, MsBBX17, and MsBBX32 using transient expression in tobacco. The results showed that BBX genes were localized in the nucleus. This study systematically analyzed the BBX gene family in Medicago plants, which provides a basis for the study of BBX gene family tolerance to abiotic stresses.

Identification, evolution, and expression of GDSL-type Esterase/Lipase (GELP) gene family in three cotton species: a bioinformatic analysis.

BACKGROUND: GDSL esterase/lipases (GELPs) play important roles in plant growth, development, and response to biotic and abiotic stresses. Presently, an extensive and in-depth analysis of GELP family genes in cotton is still not clear enough, which greatly limits the further understanding of cotton GELP function and regulatory mechanism. RESULTS: A total of 389 GELP family genes were identified in three cotton species of Gossypium hirsutum (193), G. arboreum (97), and G. raimondii (99). These GELPs could be classified into three groups and eight subgroups, with the GELPs in same group to have similar gene structures and conserved motifs. Evolutionary event analysis showed that the GELP family genes tend to be diversified at the spatial dimension and certain conservative at the time dimension, with a trend of potential continuous expansion in the future. The orthologous or paralogous GELPs among different genomes/subgenomes indicated the inheritance from genome-wide duplication during polyploidization, and the paralogous GELPs were derived from chromosomal segment duplication or tandem replication. GELP genes in the A/D subgenome underwent at least three large-scale replication events in the evolutionary process during the period of 0.6-3.2 MYA, with two large-scale evolutionary events between 0.6-1.8 MYA that were associated with tetraploidization, and the large-scale duplication between 2.6-9.1 MYA that occurred during diploidization. The cotton GELPs indicated diverse expression patterns in tissue development, ovule and fiber growth, and in response to biotic and abiotic stresses, combining the existing cis-elements in the promoter regions, suggesting the GELPs involvements of functions to be diversification and of the mechanisms to be a hormone-mediated manner. CONCLUSIONS: Our results provide a systematic and comprehensive understanding the function and regulatory mechanism of cotton GELP family, and offer an effective reference for in-depth genetic improvement utilization of cotton GELPs.

Genome-wide analysis of Hsp40 and Hsp70 gene family in four cotton species provides insights into their involvement in response to Verticillium dahliae and abiotic stress.

Introduction: Cotton is an important economic crop to provide natural fibers as raw materials to textile industry, and is significantly affected by biotic and abiotic stress during the whole growth stage, in which Verticillium wilt (VW) caused by Verticillium dahliae is one of the most destructive disease to lead to a significant yield reduction. Heat shock proteins (Hsps) are important molecular chaperones, and play crucial roles in plant growth, development, resistance to biotic and abiotic stress. Hsp40 and Hsp70 are two key Hsps in cell chaperone network, however, the function and regulatory mechanism of Hsp40 and Hsp70 members in VW resistance and abiotic stress in cotton are largely unknown. Methods and Results: Herein, a systematic and comprehensive analysis of Hsp40s and Hsp70s in four cotton species of Gossypium arboretum, G. raimondii, G. hirsutum, and G. barbadense were performed. A total of 291 Hsp40s and 171 Hsp70s identified in four Gossypium species. Sequence analysis revealed that all Hsp40 proteins contained J domain that provides the binding sites to Hsp70. Protein-protein interaction prediction analysis displayed that GhHsp40-55 might interact with GhHsp70-2 and GhHsp70-13, suggesting their potential function as protein complex. Promoter cis-acting element analysis demonstrated that multiple cis-elements related to disease and stress response consists in GhHsp40 and GhHsp70 promoters. Further expression analysis showed that eight GhHsp40s (Hsp40-2,4,8,11,20,23,53,55) and seven GhHsp70s (Hsp70-2,3,6,8,13,19,22) were up-regulated after V. dahliae infection. In addition, five GhHsp40s (Hsp40-2,8,11,53,55) and four GhHsp70s (Hsp70-3,6,8,13) were up-regulated after salt treatment, six GhHsp40s (Hsp40-4,11,20,23) and three GhHsp70s (Hsp70-2,8,19) were up-regulated after drought treatment, four GhHsp40s (Hsp40-2,11,20,23) and four GhHsp70s (Hsp70-3,6,19,22) were up-regulated after temperature treatment, suggesting these Hsps have possible important function in the process of abiotic stress response. Discussion: Our results lay a foundation for understanding the function of Hsp40 and Hsp70 in the resistance against V. dahliae and abiotic stress, and elucidating the regulatory mechanism of the protein complex, evolution and molecular mechanism under stress.

Multicopper oxidases GbAO and GbSKS are involved in the Verticillium dahliae resistance in Gossypium barbadense.

Ascorbate oxidase (AO) and skewed5 (SKU5)-similar (SKS) proteins belong to the multicopper oxidase (MCO) family and play important roles in plants in response to environmental stress via modulation of oxidoreduction homeostasis. Currently, reports on the response of Gossypium barbadense MCO to Verticillium wilt (VW) caused by Verticillium dahliae are still limited. Herein, RNA sequencing of two G. barbadense cultivars of VW-resistant XH21 and VW-susceptible XH7 under V. dahliae treatment, combined with physiological and genetic analysis, was performed to analyze the function and mechanism of multicopper oxidases GbAO and GbSKS involved in V. dahliae resistance. The identified differentially expressed genes are mainly involved in the regulation of oxidoreduction reaction, and extracellular components and signaling. Interestingly, ascorbate oxidase family members were discovered as the most significantly upregulated genes after V. dahliae treatment, including GbAO3A/D, GbSKS3A/D, and GbSKS16A/D. H2O2 and Asc contents, especially reductive Asc in both XH21 and XH7, were shown to be increased. Silenced expression of respective GbAO3A/D, GbSKS3A/D, and GbSKS16A/D in virus-induced gene silencing (VIGS) cotton plants significantly decreased the resistance to V. dahliae, coupled with the reduced contents of pectin and lignin. Our results indicate that AO might be involved in cotton VW resistance via the regulation of cell wall components.

Reconstruction of karyotypic evolution in Saccharum spontaneum species by comparative oligo-FISH mapping.

BACKGROUND: Karyotype dynamics driven by chromosomal rearrangements has long been considered as a fundamental question in the evolutionary genetics. Saccharum spontaneum, the most primitive and complex species in the genus Saccharum, has reportedly undergone at least two major chromosomal rearrangements, however, its karyotypic evolution remains unclear. RESULTS: In this study, four representative accessions, i.e., hypothetical diploid sugarcane ancestor (sorghum, x = 10), Sa. spontaneum Np-X (x = 10, tetraploid), 2012-46 (x = 9, hexaploid) and AP85-441 (x = 8, tetraploid), were selected for karyotype evolution studies. A set of oligonucleotide (oligo)-based barcode probes was developed based on the sorghum genome, which allowed universal identification of all chromosomes from sorghum and Sa. spontaneum. By comparative FISH assays, we reconstructed the karyotype evolutionary history and discovered that although chromosomal rearrangements resulted in greater variation in relative lengths of some chromosomes, all chromosomes maintained a conserved metacentric structure. Additionally, we found that the barcode oligo probe was not applicable for chromosome identification in both Sa. robustum and Sa. officinarum species, suggesting that sorghum is more distantly related to Sa. robustum and Sa. officinarum compared with Sa. spontaneum species. CONCLUSIONS: Our study demonstrated that the barcode oligo-FISH is an efficient tool for chromosome identification and karyotyping research, and expanded our understanding of the karyotypic and chromosomal evolution in the genus Saccharum.

Proteomics analysis of the effects for different salt ions in leaves of true halophyte Sesuvium portulacastrum.

Sesuvium portulacastrum is a true halophyte and shows an optimal development under moderate salinity with large amounts of salt ions in its leaves. However, the specific proteins in response to salt ions are remained unknown. In this study, comparative physiological and proteomic analyses of different leaves subject to NaCl, KCl, NaNO3 and KNO3 were performed. Chlorophyll content was decreased under the above four kinds of salt treatments. Starch and soluble sugar contents changed differently under different salt treatments. A total of 53 differentially accumulated proteins (DAPs) were identified by mass spectrometry. Among them, 13, 25, 26 and 25 DAPs were identified after exposure to KCl, NaCl, KNO3, and NaNO3, respectively. These DAPs belong to 47 unique genes, and 37 of them are involved in protein-protein interactions. These DAPs displayed different expression patterns after treating with different salt ions. Functional annotation revealed they are mainly involved in photosynthesis, carbohydrate and energy metabolism, lipid metabolism, and biosynthesis of secondary metabolites. Genes and proteins showed different expression profiles under different salt treatments. Enzyme activity analysis indicated P-ATPase was induced by KCl, NaCl and NaNO3, V-ATPase was induced by KCl and NaCl, whereas V-PPase activity was significantly increased after application of KNO3, but sharply inhibited by NaCl. These results might deepen our understanding of responsive mechanisms in the leaves of S. portulacastrum upon different salt ions.

Selection of the reference genes for quantitative gene expression by RT-qPCR in the desert plant Stipagrostis pennata.

The desert pioneer plant Stipagrostis pennata plays an important role in sand fixation, wind prevention, and desert ecosystem recovery. An absence of reference genes greatly limits investigations into the regulatory mechanism by which S. pennata adapts to adverse desert environments at the molecular and genetic levels. In this study, eight candidate reference genes were identified from rhizosheath development transcriptome data from S. pennata, and their expression stability in the rhizosheaths at different development stages, in a variety of plant tissues, and under drought stress was evaluated using four procedures, including geNorm, NormFinder, BestKeeper, and RefFinder. The results showed that GAPDH and elF were the most stable reference genes under drought stress and in rhizosheath development, and ARP6 and ALDH were relatively stable in all plant tissues. In addition, elF was the most suitable reference gene for all treatments. Analysis of the consistency between the reverse transcription-quantitative PCR (RT-qPCR) and RNA sequencing data showed that the identified elF and GAPDH reference genes were stable during rhizosheath development. These results provide reliable reference genes for assuring the accuracy of RT-qPCR and offer a foundation for further investigations into the genetic responses of S. pennata to abiotic stress.

Comparison of Morphological Characteristics and Determination of Different Patterns for Rubber Particles in Dandelion and Different Rubber Grass Varieties.

Russian dandelion Taraxacum kok-saghyz (TKS) is one promising alternative crop for natural rubber production. However, it is easily confused with other dandelions. In this study, we performed a systematical comparison of the morphological characteristics for different TKS varieties and common dandelion Taraxacum officinale (TO). Our results demonstrated that several obvious differences in morphology can be found between TKS and TO. TO leaf is a pinnate shape, its margin is heavily jagged and its base is cuneate, but TKS leaf is more cuneate and its leaf margin is nearly smooth and round. There are obvious differences for the outer bracts of TO and TKS flower buds. TKS bracts are oblanceolate, apex obtuse, margin smooth and sinuate, and its outer layer of flower buds and faceplate involucre sepal is buckled inward to form a certain angle. TKS is self-incompatible, and its seeds are spindle-shaped achene and show upright plumpness. A large amount of laticifer cells and rubber particles can be detected from many TKS tissues, and dry roots of TKS contain high contents of natural rubber. Laticifer cells and rubber particles can only be examined in the vein, stem, and roots of TKS. Our statical results also revealed that the numbers of laticifer cells and rubber particles have a positive relationship with the rubber content in TKS roots. These morphological features can help us to easily distinguish TKS from common dandelion and approximately estimate the rubber content in the roots of different TKS varieties for TKS breeding in future.

Identification and Characterization of Glycoproteins and Their Responsive Patterns upon Ethylene Stimulation in the Rubber Latex.

Natural rubber is an important industrial material, which is obtained from the only commercially cultivated rubber tree, Hevea brasiliensis. In rubber latex production, ethylene has been extensively used as a stimulant. Recent research showed that post-translational modifications (PTMs) of latex proteins, such as phosphorylation, glycosylation and ubiquitination, are crucial in natural rubber biosynthesis. In this study, comparative proteomics was performed to identify the glycosylated proteins in rubber latex treated with ethylene for different days. Combined with Pro-Q Glycoprotein gel staining and mass spectrometry techniques, we provided the first visual profiling of glycoproteomics of rubber latex and finally identified 144 glycosylated protein species, including 65 differentially accumulated proteins (DAPs) after treating with ethylene for three and/or five days. Gene Ontology (GO) functional annotation showed that these ethylene-responsive glycoproteins are mainly involved in cell parts, membrane components and metabolism. Pathway analysis demonstrated that these glycosylated rubber latex proteins are mainly involved in carbohydrate metabolism, energy metabolism, degradation function and cellular processes in rubber latex metabolism. Protein-protein interaction analysis revealed that these DAPs are mainly centered on acetyl-CoA acetyltransferase and hydroxymethylglutaryl-CoA synthase (HMGS) in the mevalonate pathway for natural rubber biosynthesis. In our glycoproteomics, three protein isoforms of HMGS2 were identified from rubber latex, and only one HMGS2 isoform was sharply increased in rubber latex by ethylene treatment for five days. Furthermore, the HbHMGS2 gene was over-expressed in a model rubber-producing grass Taraxacum Kok-saghyz and rubber content in the roots of transgenic rubber grass was significantly increased over that in the wild type plant, indicating HMGS2 is the key component for natural rubber production.

Proteomic Landscape Has Revealed Small Rubber Particles Are Crucial Rubber Biosynthetic Machines for Ethylene-Stimulation in Natural Rubber Production.

Rubber particles are a specific organelle for natural rubber biosynthesis (NRB) and storage. Ethylene can significantly improve rubber latex production by increasing the generation of small rubber particles (SRPs), regulating protein accumulation, and activating many enzyme activities. We conducted a quantitative proteomics study of different SRPs upon ethylene stimulation by differential in-gel electrophoresis (DIGE) and using isobaric tags for relative and absolute quantification (iTRAQ) methods. In DIGE, 79 differentially accumulated proteins (DAPs) were determined as ethylene responsive proteins. Our results show that the abundance of many NRB-related proteins has been sharply induced upon ethylene stimulation. Among them, 23 proteins were identified as rubber elongation factor (REF) and small rubber particle protein (SRPP) family members, including 16 REF and 7 SRPP isoforms. Then, 138 unique phosphorylated peptides, containing 129 phosphorylated amino acids from the 64 REF/SRPP family members, were identified, and most serine and threonine were phosphorylated. Furthermore, we identified 226 DAPs from more than 2000 SRP proteins by iTRAQ. Integrative analysis revealed that almost all NRB-related proteins can be detected in SRPs, and many proteins are positively responsive to ethylene stimulation. These results indicate that ethylene may stimulate latex production by regulating the accumulation of some key proteins. The phosphorylation modification of REF and SRPP isoforms might be crucial for NRB, and SRP may act as a complex natural rubber biosynthetic machine.

Proteomic Landscape of the Mature Roots in a Rubber-Producing Grass Taraxacum Kok-saghyz.

The rubber grass Taraxacum kok-saghyz (TKS) contains large amounts of natural rubber (cis-1,4-polyisoprene) in its enlarged roots and it is an alternative crop source of natural rubber. Natural rubber biosynthesis (NRB) and storage in the mature roots of TKS is a cascade process involving many genes, proteins and their cofactors. The TKS genome has just been annotated and many NRB-related genes have been determined. However, there is limited knowledge about the protein regulation mechanism for NRB in TKS roots. We identified 371 protein species from the mature roots of TKS by combining two-dimensional gel electrophoresis (2-DE) and mass spectrometry (MS). Meanwhile, a large-scale shotgun analysis of proteins in TKS roots at the enlargement stage was performed, and 3545 individual proteins were determined. Subsequently, all identified proteins from 2-DE gel and shotgun MS in TKS roots were subject to gene ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses and most proteins were involved in carbon metabolic process with catalytic activity in membrane-bounded organelles, followed by proteins with binding ability, transportation and phenylpropanoid biosynthesis activities. Fifty-eight NRB-related proteins, including eight small rubber particle protein (SRPP) and two rubber elongation factor(REF) members, were identified from the TKS roots, and these proteins were involved in both mevalonate acid (MVA) and methylerythritol phosphate (MEP) pathways. To our best knowledge, it is the first high-resolution draft proteome map of the mature TKS roots. Our proteomics of TKS roots revealed both MVA and MEP pathways are important for NRB, and SRPP might be more important than REF for NRB in TKS roots. These findings would not only deepen our understanding of the TKS root proteome, but also provide new evidence on the roles of these NRB-related proteins in the mature TKS roots.

An improved protein extraction method applied to cotton leaves is compatible with 2-DE and LC-MS.

BACKGROUND: Two-dimensional electrophoresis (2-DE) and liquid chromatography-tandem mass spectrometry (LC-MS/MS) are widely used in plant proteomics research. However, these two techniques cannot be simultaneously satisfied by traditional protein extraction methods when investigate cotton leaf proteome. RESULTS: Here, we evaluated the efficiency of three different protein extraction methods for 2-DE and LC-MS/MS analyses of total proteins obtained from cotton leaves. The protein yield of the borax/PVPP/phenol (BPP) method (0.14%) was significantly lower than the yields of the trichloroacetic acid/acetone (TCA) precipitation method (1.42%) and optimized TCA combined with BPP (TCA-B) method (0.47%). The BPP method was failed to get a clear 2-DE electrophoretogram. Fifty pairs of protein spots were randomly selected from the 2-DE gels of TCA- and TCA-B-extracted proteins for identification by MALDI TOF/TOF, and the results of 42 pairs were consistent. High-throughput proteomic analysis showed that 6339, 9282 and 9697 unique proteins were identified from the total cotton leaf proteins extracted by the TCA, BPP and TCA-B methods, respectively. Gene Ontology (GO) analysis revealed that the proteins specifically identified by TCA method were primarily distributed in the plasma membrane, while BPP and TCA-B methods specific proteins distributed in the cytosol, indicating the sub-cellular preference of different protein extraction methods. Further, ATP-dependent zinc metalloprotease FTSH 8 could be observed in the 2-DE gels of TCA and TCA-B methods, and could only be detected in the LC-MS/MS results of the BPP and TCA-B methods, showing that TCA-B method might be the optimized choice for both 2-DE and LC-MS/MS. CONCLUSION: Our data provided an improved TCA-B method for protein extraction that is compatible with 2-DE and LC-MS/MS for cotton leaves and similar plant tissues which is rich in polysaccharides and polyphenols.

Subcellular proteome profiles of different latex fractions revealed washed solutions from rubber particles contain crucial enzymes for natural rubber biosynthesis.

Rubber particle (RP) is a specific organelle for natural rubber biosynthesis (NRB) and storage in rubber tree Hevea brasiliensis. NRB is processed by RP membrane-localized proteins, which were traditionally purified by repeated washing. However, we noticed many proteins in the discarded washing solutions (WS) from RP. Here, we compared the proteome profiles of WS, C-serum (CS) and RP by 2-DE, and identified 233 abundant proteins from WS by mass spectrometry. Many spots on 2-DE gels were identified as different protein species. We further performed shotgun analysis of CS, WS and RP and identified 1837, 1799 and 1020 unique proteins, respectively. Together with 2-DE, we finally identified 1825 proteins from WS, 246 were WS-specific. These WS-specific proteins were annotated in Gene Ontology, indicating most abundant pathways are organic substance metabolic process, protein degradation, primary metabolic process, and energy metabolism. Protein-protein interaction analysis revealed these WS-specific proteins are mainly involved in ribosomal metabolism, proteasome system, vacuolar protein sorting and endocytosis. Label free and Western blotting revealed many WS-specific proteins and protein complexes are crucial for NRB initiation. These findings not only deepen our understanding of WS proteome, but also provide new evidences on the roles of RP membrane proteins in NRB. SIGNIFICANCE: Natural rubber is stored in rubber particle from the rubber tree. Rubber particles were traditionally purified by repeated washing, but many proteins were identified from the washing solutions (WS). We obtained the first visualization proteome profiles with 1825 proteins from WS, including 246 WS-specific ones. These WS proteins contain almost all enzymes for polyisoprene initiation and may play important roles in rubber biosynthesis.

Calcium-Dependent Protein Kinase Family Genes Involved in Ethylene-Induced Natural Rubber Production in Different Hevea brasiliensis Cultivars.

Natural rubber latex production can be improved by ethylene stimulation in the rubber tree (Hevea brasiliensis). However, the expression levels of most functional proteins for natural rubber biosynthesis are not induced after ethylene application, indicating that post-translational modifications, especially protein phosphorylation, may play important roles in ethylene signaling in Hevea. Here, we performed a comprehensive investigation on evolution, ethylene-induced expression and protein-protein interaction of calcium-dependent protein kinases (CPKs), an important serine/threonine protein kinase family, in Hevea. Nine duplication events were determined in the 30 identified HbCPK genes. Expression profiling of HbCPKs in three rubber tree cultivars with low, medium and high ethylene sensitivity showed that HbCPK6, 17, 20, 22, 24, 28 and 30 are induced by ethylene in at least one cultivar. Evolution rate analysis suggested accelerated evolution rates in two paralogue pairs, HbCPK9/18 and HbCPK19/20. Analysis of proteomic data for rubber latex after ethylene treatment showed that seven HbCPK proteins could be detected, including six ethylene-induced ones. Protein-protein interaction analysis of the 493 different abundant proteins revealed that protein kinases, especially calcium-dependent protein kinases, possess most key nodes of the interaction network, indicating that protein kinase and protein phosphorylation play important roles in ethylene signaling in latex of Hevea. In summary, our data revealed the expression patterns of HbCPK family members and functional divergence of two HbCPK paralogue pairs, as well as the potential important roles of HbCPKs in ethylene-induced rubber production improvement in Hevea.

Genome-wide investigation and expression profiling of APX gene family in Gossypium hirsutum provide new insights in redox homeostasis maintenance during different fiber development stages.

Ascorbate peroxidase (APX) is a member of heme-containing peroxidases which catalyze the H2O2-dependent oxidation of a wide range of substrates in plants and animals. As is known, H2O2 acts as a signaling molecule in the regulation of fiber development. Our previous work reported that ascorbate peroxidase 1 (GhAPX1) was important for cotton fiber elongation. However, knowledge about APX gene family members and their evolutionary and functional characteristics in cotton is limited. Here, we report 26 GhAPX genes by genome-wide investigation of tetraploid cotton Gossypium hirsutum. Phylogenetic and gene structure analyses classified these APX members into five clades and syntenic analysis suggested two duplication events. Expression profiling of the 26 APXs revealed that ten members are expressed in cotton fibers. Notably, GhAPX10A, GhAPX10D, GhAPX12A, and GhAPX12D showed high expression levels in 30-day fiber, while GhAPX1A/D, GhAPX3A/D, and GhAPX6A/D showed very low expression levels. The enzyme activity and H2O2 content assays revealed that cotton fiber kept high enzyme activity and the lowest H2O2 level in 30-day fibers, indicating that other than GhAPX1, the newly reported APX members are responsible for the reactive oxygen species homeostasis in the cotton fiber maturation stages. Expression profiling of ten fiber-expressed APXs after phytohormone treatments revealed their regulation patterns by different stimuli, suggesting that GhAPX1, GhAPX12A, and GhAPX12D are responsible to most phytohormone treatments. Our data provided evolutionary and functional information of GhAPX gene family members and revealed that different members are responsible to redox homeostasis during different cotton fiber development stages.

Expression Profiling of Mitogen-Activated Protein Kinase Genes Reveals Their Evolutionary and Functional Diversity in Different Rubber Tree (Hevea brasiliensis) Cultivars.

Rubber tree (Heveabrasiliensis) is the only commercially cultivated plant for producing natural rubber, one of the most essential industrial raw materials. Knowledge of the evolutionary and functional characteristics of kinases in H. brasiliensis is limited because of the long growth period and lack of well annotated genome information. Here, we reported mitogen-activated protein kinases in H.brasiliensis (HbMPKs) by manually checking and correcting the rubber tree genome. Of the 20 identified HbMPKs, four members were validated by proteomic data. Protein motif and phylogenetic analyses classified these members into four known groups comprising Thr-Glu-Tyr (TEY) and Thr-Asp-Tyr (TDY) domains, respectively. Evolutionary and syntenic analyses suggested four duplication events: HbMPK3/HbMPK6, HbMPK8/HbMPK9/HbMPK15, HbMPK10/HbMPK12 and HbMPK11/HbMPK16/HbMPK19. Expression profiling of the identified HbMPKs in roots, stems, leaves and latex obtained from three cultivars with different latex yield ability revealed tissue- and variety-expression specificity of HbMPK paralogues. Gene expression patterns under osmotic, oxidative, salt and cold stresses, combined with cis-element distribution analyses, indicated different regulation patterns of HbMPK paralogues. Further, Ka/Ks and Tajima analyses suggested an accelerated evolutionary rate in paralogues HbMPK10/12. These results revealed HbMPKs have diverse functions in natural rubber biosynthesis, and highlighted the potential possibility of using MPKs to improve stress tolerance in future rubber tree breeding.