The calcium-dependent protein kinase (CDPK) and CDPK-related kinase gene families in Hevea brasiliensis-comparison with five other plant species in structure, evolution, and expression

The rubber tree kinome: Genome-wide characterization and insights into coexpression patterns associated with abiotic stress responses

The protein kinase (PK) superfamily constitutes one of the largest and most conserved protein families in eukaryotic genomes, comprising core components of signaling pathways in cell regulation. Despite its remarkable relevance, only a few kinase families have been studied in Hevea brasiliensis. A comprehensive characterization and global expression analysis of the PK superfamily, however, is currently lacking. In this study, with the aim of providing novel inferences about the mechanisms associated with the stress response developed by PKs and retained throughout evolution, we identified and characterized the entire set of PKs, also known as the kinome, present in the Hevea genome. Different RNA-sequencing datasets were employed to identify tissue-specific expression patterns and potential correspondences between different rubber tree genotypes. In addition, coexpression networks under several abiotic stress conditions, such as cold, drought and latex overexploitation, were employed to elucidate associations between families and tissues/stresses. A total of 1,809 PK genes were identified using the current reference genome assembly at the scaffold level, and 1,379 PK genes were identified using the latest chromosome-level assembly and combined into a single set of 2,842 PKs. These proteins were further classified into 20 different groups and 122 families, exhibiting high compositional similarities among family members and with two phylogenetically close species Manihot esculenta and Ricinus communis. Through the joint investigation of tandemly duplicated kinases, transposable elements, gene expression patterns, and coexpression events, we provided insights into the understanding of the cell regulation mechanisms in response to several conditions, which can often lead to a significant reduction in rubber yield.

Functional Characterization of the Cystine-Rich-Receptor-like Kinases (CRKs) and Their Expression Response to Sclerotinia sclerotiorum and Abiotic Stresses in Brassica napus

Cysteine-rich receptor-like kinases (CRKs) are transmembrane proteins that bind to the calcium ion to regulate stress-signaling and plant development-related pathways, as indicated by several pieces of evidence. However, the CRK gene family hasn’t been inadequately examined in Brassica napus. In our study, 27 members of the CRK gene family were identified in Brassica napus, which are categorized into three phylogenetic groups and display synteny relationship to the Arabidopsis thaliana orthologs. All the CRK genes contain highly conserved N-terminal PKINASE domain; however, the distribution of motifs and gene structure were variable conserved. The functional divergence analysis between BnaCRK groups indicates a shift in evolutionary rate after duplication events, demonstrating that BnaCRKs might direct a specific function. RNA-Seq datasets and quantitative real-time PCR (qRT-PCR) exhibit the complex expression profile of the BnaCRKs in plant tissues under multiple stresses. Nevertheless, BnaA06CRK6-1 and BnaA08CRK8 from group B were perceived to play a predominant role in the Brassica napus stress signaling pathway in response to drought, salinity, and Sclerotinia sclerotiorum infection. Insights gained from this study improve our knowledge about the Brassica napus CRK gene family and provide a basis for enhancing the quality of rapeseed.

Transcription level differences in Taxus wallichiana var. mairei elicited by Ce3+, Ce4+ and methyl jasmonate

Taxol is a precious and effective anticancer drug. Cerium and methyl jasmonate (MJ) have been shown to increase the yield of taxol in taxus cells. However, the mechanisms of cerium-mediated and MJ-mediated taxol biosynthesis remain unknown. RNA-Seq was applied to study the overall regulation mechanism of cerium and MJ on taxol biosynthesis and analyze the differences among T. mairei cells elicited by Ce³⁺, Ce⁴⁺ and MJ on transcriptional level . Using sequence homology, 179 unigenes were identified as taxol synthesis genes. Under the condition of 100 μM MJ, taxol synthesis genes were up-regulated. Notably, taxol synthesis genes were down-regulated expression at 1 mM Ce³⁺ and 1 mM Ce⁴⁺. Differential expression genes involved in some related functions were analyzed, such as MAPK signaling pathway and plant-pathogen interaction. Sequence alignment and phylogenetic analysis of nine differentially expressed WRKYs in our data were carried out.

The role of CDPKs in plant development, nutrient and stress signaling

The second messenger calcium (Ca2+) is a ubiquitous intracellular signaling molecule found in eukaryotic cells. In plants, the multigene family of calcium-dependent protein kinases (CDPKs) plays an important role in regulating plant growth, development, and stress tolerance. CDPKs sense changes in intracellular Ca2+ concentration and translate them into phosphorylation events that initiate downstream signaling processes. Several functional and expression studies on different CDPKs and their encoding genes have confirmed their multifunctional role in stress. Here, we provide an overview of the signal transduction mechanisms and functional roles of CDPKs. This review includes details on the regulation of secondary metabolites, nutrient uptake, regulation of flower development, hormonal regulation, and biotic and abiotic stress responses.

Identification and Expression Analysis of Calcium-Dependent Protein Kinases Gene Family in Potato Under Drought Stress

Calcium-dependent protein kinases (CDPKs) are a class of serine/threonine protein kinases encoded by several gene families that play key roles in stress response and plant growth and development. In this study, the BLAST method was used to search for protein sequences of the potato Calcium-dependent protein kinase gene family. The chromosome location, phylogeny, gene structures, gene duplication, cis-acting elements, protein-protein interaction, and expression profiles were analyzed. Twenty-five CDPK genes in the potato genome were identified based on RNA-seq data and were clustered into four groups (I-IV) based on their structural features and phylogenetic analysis. The result showed the composition of the promoter region of the StCDPKs gene, including light-responsive elements such as Box4, hormone-responsive elements such as ABRE, and stress-responsive elements such as MBS. Four pairs of segmental duplications were found in StCDPKs genes and the Ka/Ks ratios were below 1, indicating a purifying selection of the genes. The protein-protein interaction network revealed defense-related proteins such as; respiratory burst oxidase homologs (RBOHs) interacting with potato CDPKs. Transcript abundance was measured via RT-PCR between the two cultivars and their relative expression of CDPK genes was analyzed after 15, 20, and 25 days of drought. There were varied expression patterns of StCDPK3/13/21 and 23, between the two potato cultivars under mannitol induced-drought conditions. Correlation analysis showed that StCDPK21/22 and StCDPK3 may be the major differentially expressed genes involved in the regulation of malondialdehyde (MDA) and proline content in response to drought stress, opening a new research direction for genetic improvement of drought resistance in potato.

Identification and expression assay of calcium-dependent protein kinase family genes in Hevea brasiliensis and determination of HbCDPK5 functions in disease resistance

Calcium (Ca2+) signaling is one of the earliest factors to coordinate plant adaptive responses. As direct sensors and activators of Ca2+ signals, calcium-dependent protein kinases (CDPKs) were reported to be widely involved in regulating different biotic and abiotic stress stimuli. In this study, 32 Hevea brasiliensis CDPK (HbCDPK) genes were predicted and classified into four subgroups. Among them, the full-length coding sequences of 28 HbCDPK genes were confirmed by RT-PCR and verified by sequencing. Putative cis-elements assay in the promoters of HbCDPKs showed that most of the HbCDPK genes contained gibberellic acid-responsive element (GARE), abscisic acid-responsive element (ABRE), salicylic acid-responsive element (SARE), defense and stress responsive element (TC-rich repeats) and low-temperature response element (LTR), which could be activated by different biotic and abiotic stresses. Real-time PCR analysis indicated that 28 HbCDPK genes respond to infection of pathogenic fungi and a variety of phytohormones. Subcellular localization was observed with most HbCDPKs located in cell membrane, cytoplasm or organelles. Some HbCDPKs were confirmed to cause reactive oxygen species (ROS) production and accumulation in rubber tree mesophyll protoplast directly. HbCDPK5 was strongly induced by the inoculation with Colletotrichum gloeosporioides and was chosen for further analysis. HbCDPK5 localized to the cell membrane and cytoplasm, and obviously regulated the accumulation of ROS in rubber tree mesophyll protoplast. Overexpression of HbCDPK5 in Arabidopsis enhanced the resistance to Botrytis cinerea. These results indicate that rubber tree CDPK genes play important roles in plant disease resistance.

Molecular Characterization Reveals the Involvement of Calcium Dependent Protein Kinases in Abiotic Stress Signaling and Development in Chickpea (Cicer arietinum)

Calcium-dependent protein kinases (CDPKs) are a major group of calcium (Ca²⁺⁾ sensors in plants. CDPKs play a dual function of "Ca²⁺ sensor and responder." These sensors decode the "Ca²⁺ signatures" generated in response to adverse growth conditions such as drought, salinity, and cold and developmental processes. However, knowledge of the CDPK family in the legume crop chickpea is missing. Here, we have identified a total of 22 CDPK genes in the chickpea genome. The phylogenetic analysis of the chickpea CDPK family with other plants revealed their evolutionary conservation. Protein homology modeling described the three-dimensional structure of chickpea CDPKs. Defined arrangements of α-helix, β-strands, and transmembrane-helix represent important structures like kinase domain, inhibitory junction domain, N and C-lobes of EF-hand motifs. Subcellular localization analysis revealed that CaCDPK proteins are localized mainly at the cytoplasm and in the nucleus. Most of the CaCDPK promoters had abiotic stress and development-related cis-regulatory elements, suggesting the functional role of CaCDPKs in abiotic stress and development-related signaling. RNA sequencing (RNA-seq) expression analysis indicated the role of the CaCDPK family in various developmental stages, including vegetative, reproductive development, senescence stages, and during seed stages of early embryogenesis, late embryogenesis, mid and late seed maturity. The real-time quantitative PCR (qRT-PCR) analysis revealed that several CaCDPK genes are specifically as well as commonly induced by drought, salt, and Abscisic acid (ABA). Overall, these findings indicate that the CDPK family is probably involved in abiotic stress responses and development in chickpeas. This study provides crucial information on the CDPK family that will be utilized in generating abiotic stress-tolerant and high-yielding chickpea varieties.

The Calcium Sensor Calcineurin B-Like Proteins -Calcineurin B-Like Interacting Protein Kinases Is Involved in Leaf Development and Stress Responses Related to Latex Flow in Hevea brasiliensis

Latex flow in Hevea brasiliensis (the Para rubber tree), the sole commercial source of natural rubber (cis-1,4-polyisoprene, NR), renders it uniquely suited for the study of plant stress responses. Calcineurin B-like interacting protein kinases (CIPK) serving as calcium-sensor protein kinases react with calcineurin B-like proteins (CBL) to play crucial roles in hormone signaling transduction and response to abiotic stress in plant developmental processes. However, little is known about their functions in Hevea. In this study, a total of twelve CBL (HbCBL) and thirty CIPK (HbCIPK) genes were identified from the Hevea genome. Structure and phylogenetic analysis assigned these CIPKs to five groups and CBLs to four groups, and mapped onto fourteen of the eighteen Hevea chromosomes. RNA-seq and qPCR analysis showed that the expressions of HbCBL and HbCIPK genes varied in the seven Hevea tissues examined, i.e., latex (cytoplasm of rubber-producing laticifers), bark, leaf, root, seed, female flower, and male flower. The expressions of two HbCBL and sixteen HbCIPK genes showed upward trends during leaf development. Following ethylene yield stimulation and the latex tapping treatment, both practices invoking stress, the expression levels of most latex-expressed genes were significantly altered. Yeast two-hybrid test revealed interactions for multiple combinations of HbCBLs and HbCIPKs with substantial gene expression in latex or other Hevea tissues. However, all the HbCBL-HbCIPK complexes examined did not recruit HbSOS1 or AtSOS1 to form functional salt tolerance SOS pathway in yeast cells. Taken together, the results suggested a role of the Hevea CBL-CIPK network as a point of convergence for several different signaling pathways in growth, development, and stress responses in relation to latex production.

Plant protein-coding gene families: Their origin and evolution

Steady advances in genome sequencing methods have provided valuable insights into the evolutionary processes of several gene families in plants. At the core of plant biodiversity is an extensive genetic diversity with functional divergence and expansion of genes across gene families, representing unique phenomena. The evolution of gene families underpins the evolutionary history and development of plants and is the subject of this review. We discuss the implications of the molecular evolution of gene families in plants, as well as the potential contributions, challenges, and strategies associated with investigating phenotypic alterations to explain the origin of plants and their tolerance to environmental stresses.

Molecular Evolution of Calcium Signaling and Transport in Plant Adaptation to Abiotic Stress

Adaptation to unfavorable abiotic stresses is one of the key processes in the evolution of plants. Calcium (Ca2+) signaling is characterized by the spatiotemporal pattern of Ca2+ distribution and the activities of multi-domain proteins in integrating environmental stimuli and cellular responses, which are crucial early events in abiotic stress responses in plants. However, a comprehensive summary and explanation for evolutionary and functional synergies in Ca2+ signaling remains elusive in green plants. We review mechanisms of Ca2+ membrane transporters and intracellular Ca2+ sensors with evolutionary imprinting and structural clues. These may provide molecular and bioinformatics insights for the functional analysis of some non-model species in the evolutionarily important green plant lineages. We summarize the chronological order, spatial location, and characteristics of Ca2+ functional proteins. Furthermore, we highlight the integral functions of calcium-signaling components in various nodes of the Ca2+ signaling pathway through conserved or variant evolutionary processes. These ultimately bridge the Ca2+ cascade reactions into regulatory networks, particularly in the hormonal signaling pathways. In summary, this review provides new perspectives towards a better understanding of the evolution, interaction and integration of Ca2+ signaling components in green plants, which is likely to benefit future research in agriculture, evolutionary biology, ecology and the environment.

High-density SNP map facilitates fine mapping of QTLs and candidate genes discovery for Aspergillus flavus resistance in peanut (Arachis hypogaea)

Two novel resistant QTLs mapped and candidate genes identified for Aspergillus flavus resistance in cultivated peanut using SLAF-seq. Aflatoxin contamination in peanuts caused by Aspergillus flavus is a serious food safety issue for human health around the world. Host plant resistance to fungal infection and reduction in aflatoxin are crucial for mitigating this problem. Identification of the resistance-linked markers can be used in marker-assisted breeding for varietal development. Here we report construction of two high-density genetic linkage maps with 1975 SNP loci and 5022 SNP loci, respectively. Two consistent quantitative trait loci (QTL) were identified as qRAF-3-1 and qRAF-14-1, which located on chromosomes A03 and B04, respectively. QTL qRAF-3-1 was mapped within 1.67 cM and had more than 19% phenotypic variance explained (PVE), while qRAF-14-1 was located within 1.34 cM with 5.15% PVE. While comparing with the reference genome, the mapped QTLs, qRAF-3-1 and qRAF-14-1, were located within a physical distance of 1.44 Megabase pair (Mbp) and 2.22 Mbp, harboring 67 and 137 genes, respectively. Among the identified candidate genes, six genes with the same function were found within both QTLs regions. In addition, putative disease resistance RPP13-like protein 1 (RPP13), lipoxygenase (Lox), WRKY transcription factor (WRKY) and cytochrome P450 71B34 genes were also identified. Using microarray analysis, genes responded to A. flavus infection included coding for RPP13, pentatricopeptide repeat-containing-like protein, and Lox which may be possible candidate genes for resistance to A. flavus. The QTLs and candidate genes will further facilitate marker development and validation of genes for deployment in the molecular breeding programs against A. flavus in peanuts.

Genome-wide identification and expression analysis of the phosphatase 2A family in rubber tree (Hevea brasiliensis)

The protein phosphatase 2As (PP2As) play a key role in manipulating protein phosphorylation. Although a number of proteins in the latex of laticifers are phosphorylated during latex regeneration in rubber tree, information about the PP2A family is limited. In the present study, 36 members of the HbPP2A family were genome-wide identified. They were clustered into five subgroups: the subgroup HbPP2AA (4), HbPP2AB' (14), HbPP2AB'' (6), HbPP2AB55 (4), and HbPP2AC (8). The members within the same subgroup shared highly conserved gene structures and protein motifs. Most of HbPP2As possessed ethylene- and wounding-responsive cis-acting elements. The transcripts of 29 genes could be detected in latex by using published high-throughput sequencing data. Of the 29 genes, seventeen genes were significantly down-regulated while HbPP2AA1-1 and HbPP2AB55α/Bα-1were up-regulated by tapping. Of the 17 genes, 14 genes were further significantly down-regulated by ethrel application. The down-regulated expression of a large number of HbPP2As may attribute to the enhanced phosphorylation of the proteins in latex from the tapped trees and the trees treated with ethrel application.

Evolutionary Analysis of Calcium-Dependent Protein Kinase in Five Asteraceae Species

Calcium-dependent protein kinase (CPK) is crucial in Ca²⁺ signal transduction, and is a large gene family in plants. In our previous work, we reported Hevea brasiliensis CPKs were important for natural rubber biosynthesis. However, this CPK gene family in other rubber producing plants has not been investigated. Here, we report the CPKs in five representative Asteraceae species, including three rubber-producing and two non-rubber species. A total of 34, 34, 40, 34 and 30 CPKs were identified from Taraxacum koksaghyz, Lactuca sativa, Helianthus annuus, Chrysanthemum nankingense and Cynara cardunculus, respectively. All CPKs were classified into four individual groups (group I to IV). In addition, 10 TkCPKs, 11 LsCPKs, 20 HaCPKs, 13 CnCPKs and 7 CcCPKs duplicated paralogs were identified. Further evolutionary analysis showed that, compared to other subfamilies, the group III had been expanded in the Asteraceae species, especially in the rubber-producing species. Meanwhile, the CPKs in group III from Asteraceae species tend to expand with low calcium binding capacity. This study provides a systematical evolutionary investigation of the CPKs in five representative Asteraceae species, suggesting that the sub-family specific expansion of CPKs might be related to natural rubber producing.

Insights on Calcium-Dependent Protein Kinases (CPKs) Signaling for Abiotic Stress Tolerance in Plants

Abiotic stresses are the major limiting factors influencing the growth and productivity of plants species. To combat these stresses, plants can modify numerous physiological, biochemical, and molecular processes through cellular and subcellular signaling pathways. Calcium-dependent protein kinases (CDPKs or CPKs) are the unique and key calcium-binding proteins, which act as a sensor for the increase and decrease in the calcium (Ca) concentrations. These Ca flux signals are decrypted and interpreted into the phosphorylation events, which are crucial for signal transduction processes. Several functional and expression studies of different CPKs and their encoding genes validated their versatile role for abiotic stress tolerance in plants. CPKs are indispensable for modulating abiotic stress tolerance through activation and regulation of several genes, transcription factors, enzymes, and ion channels. CPKs have been involved in supporting plant adaptation under drought, salinity, and heat and cold stress environments. Diverse functions of plant CPKs have been reported against various abiotic stresses in numerous research studies. In this review, we have described the evaluated functions of plant CPKs against various abiotic stresses and their role in stress response signaling pathways.

Calcium-mediation of jasmonate biosynthesis and signaling in plants

Jasmonates (JAs) play vital roles in regulating a range of plant growth and development processes including seed germination, seedling development, reproduction, formation and development of storage organs, and senescence. JAs are also involved in the regulation of plant responses to environmental stimuli. The biosynthesis of JAs takes place in three different subcellular compartments, namely, the chloroplast, peroxisome, and cytoplasm. JAs activate the expression of JA-responsive genes by degrading jasmonate zinc-finger-inflorescence meristem (Zim) domain (JAZ) repressors via the E3 ubiquitin-ligase Skp/Cullin/F-box protein CORONATINE INSENSITIVE1 (COI1) complex (SCF^COI1) by using 26S proteasome. Calcium, reactive oxygen species (ROS), mitogen-activated protein kinase (MAPK), and nitric oxide (NO) are involved in the regulation of the biosynthesis and signaling of JAs in plants. Among these signaling molecules, calcium is one of the most important within plant cells. In plants, intracellular calcium levels change in response to JAs, resulting in calcium signatures with temporal and spatial features. Calcium channels are involved in the generation of calcium signatures. Calcium sensors, including calmodulins (CaMs), CaM-like proteins (CMLs), calcineurin B-like proteins (CBLs), and calcium-dependent protein kinases (CDPKs), can act to regulate the biosynthesis and signaling of JAs.

Calcium-Dependent Protein Kinase Genes in Glycyrrhiza Uralensis Appear to be Involved in Promoting the Biosynthesis of Glycyrrhizic Acid and Flavonoids under Salt Stress

As calcium signal sensors, calcium-dependent protein kinases (CPKs) play vital roles in stimulating the production of secondary metabolites to participate in plant development and response to environmental stress. However, investigations of the Glycyrrhiza uralensis CPK family genes and their multiple functions are rarely reported. In this study, a total of 23 GuCPK genes in G. uralensis were identified, and their phylogenetic relationships, evolutionary characteristics, gene structure, motif distribution, and promoter cis-acting elements were analyzed. Ten GuCPKs showed root-specific preferential expressions, and GuCPKs indicated different expression patterns under treatments of CaCl2 and NaCl. In addition, under 2.5 mM of CaCl2 and 30 mM of NaCl treatments, the diverse, induced expression of GuCPKs and significant accumulations of glycyrrhizic acid and flavonoids suggested the possible important function of GuCPKs in regulating the production of glycyrrhizic acid and flavonoids. Our results provide a genome-wide characterization of CPK family genes in G. uralensis, and serve as a foundation for understanding the potential function and regulatory mechanism of GuCPKs in promoting the biosynthesis of glycyrrhizic acid and flavonoids under salt stress.

Functional Analysis of the Arabidopsis thaliana CDPK-Related Kinase Family: AtCRK1 Regulates Responses to Continuous Light

The Calcium-Dependent Protein Kinase (CDPK)-Related Kinase family (CRKs) consists of eight members in Arabidopsis. Recently, AtCRK5 was shown to play a direct role in the regulation of root gravitropic response involving polar auxin transport (PAT). However, limited information is available about the function of the other AtCRK genes. Here, we report a comparative analysis of the Arabidopsis CRK genes, including transcription regulation, intracellular localization, and biological function. AtCRK transcripts were detectable in all organs tested and a considerable variation in transcript levels was detected among them. Most AtCRK proteins localized at the plasma membrane as revealed by microscopic analysis of 35S::cCRK-GFP (Green Fluorescence Protein) expressing plants or protoplasts. Interestingly, 35S::cCRK1-GFP and 35S::cCRK7-GFP had a dual localization pattern which was associated with plasma membrane and endomembrane structures, as well. Analysis of T-DNA insertion mutants revealed that AtCRK genes are important for root growth and control of gravitropic responses in roots and hypocotyls. While Atcrk mutants were indistinguishable from wild type plants in short days, Atcrk1-1 mutant had serious growth defects under continuous illumination. Semi-dwarf phenotype of Atcrk1-1 was accompanied with chlorophyll depletion, disturbed photosynthesis, accumulation of singlet oxygen, and enhanced cell death in photosynthetic tissues. AtCRK1 is therefore important to maintain cellular homeostasis during continuous illumination.

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 Identification and Functional Analyses of the CRK Gene Family in Cotton Reveals GbCRK18 Confers Verticillium Wilt Resistance in Gossypium barbadense

Cysteine-rich receptor-like kinases (CRKs) are a large subfamily of plant receptor-like kinases that play a critical role in disease resistance in plants. However, knowledge about the CRK gene family in cotton and its function against Verticillium wilt (VW), a destructive disease caused by Verticillium dahliae that significantly reduces cotton yields is lacking. In this study, we identified a total of 30 typical CRKs in a Gossypium barbadense genome (GbCRKs). Eleven of these (>30%) are located on the A06 and D06 chromosomes, and 18 consisted of 9 paralogous pairs encoded in the A and D subgenomes. Phylogenetic analysis showed that the GbCRKs could be classified into four broad groups, the expansion of which has probably been driven by tandem duplication. Gene expression profiling of the GbCRKs in resistant and susceptible cotton cultivars revealed that a phylogenetic cluster of nine of the GbCRK genes were up-regulated in response to V. dahliae infection. Virus-induced gene silencing of each of these nine GbCRKs independently revealed that the silencing of GbCRK18 was sufficient to compromise VW resistance in G. barbadense. GbCRK18 expression could be induced by V. dahliae infection or jasmonic acid, and displayed plasma membrane localization. Therefore, our expression analyses indicated that the CRK gene family is differentially regulated in response to Verticillium infection, while gene silencing experiments revealed that GbCRK18 in particular confers VW resistance in G. barbadense.

The SWEET gene family in Hevea brasiliensis - its evolution and expression compared with four other plant species

SWEET proteins play an indispensable role as a sugar efflux transporter in plant development and stress responses. The SWEET genes have previously been characterized in several plants. Here, we present a comprehensive analysis of this gene family in the rubber tree, Hevea brasiliensis. There are 36 members of the SWEET gene family in this species, making it one of the largest families in plant genomes sequenced so far. Structure and phylogeny analyses of these genes in Hevea and in other species demonstrated broad evolutionary conservation. RNA‐seq analyses revealed that SWEET2, 16, and 17 might represent the main evolutionary direction of SWEET genes in plants. Our results in Hevea suggested the involvement of HbSWEET1a, 2e, 2f, and 3b in phloem loading, HbSWEET10a and 16b in laticifer sugar transport, and HbSWEET9a in nectary‐specific sugar transport. Parallel studies of RNA‐seq analyses extended to three other plant species (Manihot esculenta, Populus trichocarpa, and Arabidopsis thaliana) produced findings which implicated MeSWEET10a, 3a, and 15b in M. esculenta storage root development, and the involvement of PtSWEET16b and PtSWEET16d in P. trichocarpa xylem development. RT‐qPCR results further revealed that HbSWEET10a, 16b, and 1a play important roles in phloem sugar transport. The results from this study provide a foundation not only for further investigation into the functionality of the SWEET gene family in Hevea, especially in its sugar transport for latex production, but also for related studies of this gene family in the plant kingdom.