Identification of expression profiles of tapping panel dryness (TPD) associated genes from the latex of rubber tree (Hevea brasiliensis Muell. Arg.)

Bark transcriptome analyses reveals molecular mechanisms involved in tapping panel dryness occurrence and development in rubber tree (Hevea brasiliensis)

Tapping panel dryness (TPD) has become the mostimportant limiting factor for increasing natural rubber yield, whereas illuminating the molecular mechanisms underlying TPD is the prerequisite for solving the problem of TPD. However, molecular mechanisms underlying TPD are largely unknown. In this study, healthy and different stages of TPD-affected rubber trees were utilized to analyze TPD for the first time. We found that the changing tendencies of key latex physiological parameters were closely related to TPD occurrence and development. To reveal the molecular mechanisms underlying TPD, we sequenced and compared bark transcriptomes among healthy rubber tree, and TPD-affected ones at initial and advanced stages. In total, 8607 genes were identified as TPD-related genes in contrast to healthy rubber tree. According to gene expression profiles, the five samples were divided into three groups including healthy rubber tree, and TPD-affected rubber tree in the initial and advanced stages, which was consistent with the stages of TPD occurrence and development. Interestingly, only asmall proportionof the TPD-related genes were constantly down- or up-regulated with TPD occurrence and development. The TPD-related genes in KEGG pathways significantly enriched were closely associated with protein metabolism, cell division and differentiation, PCD, stress responses, terpene biosynthesis, and various metabolism processes. Moreover, overexpression of HbAPX2 identified as a TPD-related gene enhanced oxidative stress tolerance in S. cerevisiae. The typical symptoms of TPD, partial or complete dry zone (no latex flow) on tapping panel, might attribute to lower IPP available for rubber biosynthesis, and downregulation of the genes in post-IPP steps of rubber biosynthesis and the genes involved in latex flow. Our results not only provide new insights into molecular mechanisms underlying TPD occurrence and development but also contribute to developing effective measures to control TPD in rubber trees.

Construction and analysis of the tapping panel dryness-related lncRNA/circRNA-miRNA-mRNA ceRNA network in latex of Hevea brasiliensis

Tapping panel dryness (TPD) results in a severe reduction in latex yield in Hevea brasiliensis. However, the molecular regulatory mechanisms of TPD occurrence are still largely unclear. In this study, whole-transcriptome sequencing was carried out on latex from TPD and healthy trees. In total, 7078 long noncoding RNAs (lncRNAs), 3077 circular RNAs (circRNAs), 4956 miRNAs, and 25041 mRNAs were identified in latex, among which 435 lncRNAs, 68 circRNAs, 320 miRNAs, and 1574 mRNAs were differentially expressed in the latex of TPD trees. GO and KEGG analyses indicated that plant hormone signal transduction, MAPK signaling pathway, and ubiquitin-mediated proteolysis were the key pathways associated with TPD onset. Phytohormone profiling revealed significant changes in the contents of 28 hormonal compounds, among which ACC, ABA, IAA, GA, and JA contents were increased, while SA content was reduced in TPD latex, suggesting that hormone homeostasis is disrupted in TPD trees. Furthermore, we constructed a TPD-related competitive endogenous RNA (ceRNA) regulatory network of lncRNA/circRNA-miRNA-mRNA with 561 edges and 434 nodes (188 lncRNAs, 5 circRNAs, 191 miRNAs, and 50 mRNAs) and identified two hub lncRNAs (MSTRG.11908.1 and MSTRG.8791.1) and four hub miRNAs (hbr-miR156, miR156-x, miRf10477-y, and novel-m0452-3p). Notably, the lncRNA-miR156/157-SPL module containing three hubs probably plays a crucial role in TPD onset. The expression of network hubs and the lncRNA-miR156/157-SPL module were further validated by qRT-PCR. Our results reveal the TPD-associated ceRNA regulatory network of lncRNA/circRNA-miRNA-mRNA in latex and lay a foundation for further investigation of molecular regulatory mechanisms for TPD onset in H. brasiliensis.

Genome-wide identification of lncRNAs, miRNAs, mRNAs and their regulatory networks involved in tapping panel dryness in rubber tree (Hevea brasiliensis)

Noncoding RNAs (ncRNAs) play pivotal roles in various biological processes in plants. However, the role of ncRNAs in tapping panel dryness (TPD) of rubber tree (Hevea brasiliensis Muell. Arg.) is largely unknown. Here, the whole transcriptome analyses of bark tissues from healthy and TPD trees were performed to identify differentially expressed long ncRNAs (DELs), microRNAs/miRNAs (DEMs), genes (DEGs) and their regulatory networks involved in TPD. A total of 263 DELs, 174 DEMs and 1574 DEGs were identified in the bark of TPD tree compared with that of healthy tree. Kyoto Encyclopedia of Genes and Genomes analysis revealed that most of the DEGs and targets of DELs and DEMs were mainly enriched in metabolic pathways, biosynthesis of secondary metabolites and plant hormone signal transduction. Additionally, the majority of DEGs and DELs related to rubber biosynthesis were downregulated in TPD trees. Furthermore, 98 DEGs and 44 DELs were targeted by 54 DEMs, 190 DEGs were identified as putative targets of 56 DELs, and 2 and 44 DELs were predicted as precursors and endogenous target mimics of 2 and 6 DEMs, respectively. Based on these, the DEL-DEM-DEG regulatory network involved in TPD was constructed, and 13 hub DELs, 3 hub DEMs and 2 hub DEGs were identified. The results provide novel insights into the regulatory roles of ncRNAs underlying TPD and lay a foundation for future functional characterization of long ncRNAs, miRNAs and genes involved in TPD in rubber tree.

Transcriptome analysis of Pará rubber tree (H. brasiliensis) seedlings under ethylene stimulation

BACKGROUND

Natural rubber (cis-1,4-polyioprene, NR) is an indispensable industrial raw material obtained from the Pará rubber tree (H. brasiliensis). Natural rubber cannot be replaced by synthetic rubber compounds because of the superior resilience, elasticity, abrasion resistance, efficient heat dispersion, and impact resistance of NR. In NR production, latex is harvested by periodical tapping of the trunk bark. Ethylene enhances and prolongs latex flow and latex regeneration. Ethephon, which is an ethylene-releasing compound, applied to the trunk before tapping usually results in a 1.5- to 2-fold increase in latex yield. However, intense mechanical damage to bark tissues by excessive tapping and/or over-stimulation with ethephon induces severe oxidative stress in laticifer cells, which often causes tapping panel dryness (TPD) syndrome. To enhance NR production without causing TPD, an improved understanding of the molecular mechanism of the ethylene response in the Pará rubber tree is required. Therefore, we investigated gene expression in response to ethephon treatment using Pará rubber tree seedlings as a model system.

RESULTS

After ethephon treatment, 3270 genes showed significant differences in expression compared with the mock treatment. Genes associated with carotenoids, flavonoids, and abscisic acid biosynthesis were significantly upregulated by ethephon treatment, which might contribute to an increase in latex flow. Genes associated with secondary cell wall formation were downregulated, which might be because of the reduced sugar supply. Given that sucrose is an important molecule for NR production, a trade-off may arise between NR production and cell wall formation for plant growth and for wound healing at the tapping panel.

CONCLUSIONS

Dynamic changes in gene expression occur specifically in response to ethephon treatment. Certain genes identified may potentially contribute to latex production or TPD suppression. These data provide valuable information to understand the mechanism of ethylene stimulation, and will contribute to improved management practices and/or molecular breeding to attain higher yields of latex from Pará rubber trees.

Post-transcriptional regulation of several biological processes involved in latex production in Hevea brasiliensis

Background

Small RNAs modulate plant gene expression at both the transcriptional and post-transcriptional level, mostly through the induction of either targeted DNA methylation or transcript cleavage, respectively. Small RNA networks are involved in specific plant developmental processes, in signaling pathways triggered by various abiotic stresses and in interactions between the plant and viral and non-viral pathogens. They are also involved in silencing maintenance of transposable elements and endogenous viral elements. Alteration in small RNA production in response to various environmental stresses can affect all the above-mentioned processes. In rubber trees, changes observed in small RNA populations in response to trees affected by tapping panel dryness, in comparison to healthy ones, suggest a shift from a transcriptional to a post-transcriptional regulatory pathway. This is the first attempt to characterise small RNAs involved in post-transcriptional silencing and their target transcripts in Hevea.

Methods

Genes producing microRNAs (MIR genes) and loci producing trans-activated small interfering RNA (ta-siRNA) were identified in the clone PB 260 re-sequenced genome. Degradome libraries were constructed with a pool of total RNA from six different Hevea tissues in stressed and non-stressed plants. The analysis of cleaved RNA data, associated with genomics and transcriptomics data, led to the identification of transcripts that are affected by 20–22 nt small RNA-mediated post-transcriptional regulation. A detailed analysis was carried out on gene families related to latex production and in response to growth regulators.

Results

Compared to other tissues, latex cells had a higher proportion of transcript cleavage activity mediated by miRNAs and ta-siRNAs. Post-transcriptional regulation was also observed at each step of the natural rubber biosynthesis pathway. Among the genes involved in the miRNA biogenesis pathway, our analyses showed that all of them are expressed in latex. Using phylogenetic analyses, we show that both the Argonaute and Dicer-like gene families recently underwent expansion. Overall, our study underlines the fact that important biological pathways, including hormonal signalling and rubber biosynthesis, are subject to post-transcriptional silencing in laticifers.

Complete genome sequence of a novel capillovirus infecting Hevea brasiliensis in China

Tapping panel dryness (TPD) is a complex disorder that causes partial or complete cessation of latex drainage upon tapping of rubber trees (Hevea brasiliensis). In this work, we determined the complete genome sequences of a novel virus identified in a rubber tree with TPD syndrome in China. The genome of the virus consists of 6811 nt and possesses two overlapping open reading frames (ORF1 and ORF2), encoding a polyprotein and a movement protein, respectively. The polyprotein shares 37% amino acid sequence identity with cherry virus A (CVA, ARQ83874.1) over 99% coverage. The genome architecture is similar to that of members of the genus Capillovirus (family Betaflexiviridae). Phylogenetic analysis of the replicase proteins showed that the virus clustered together with members of the genus Capillovirus. The new virus is tentatively called "rubber tree virus 1" (RTV1). RTV1 is the first virus reported to infect rubber trees. This work lays a foundation for research into finding the potential causal agent of TPD in Hevea brasiliensis.

Identification of the targets of HbEIN3/EILs in genomic wide in Hevea brasiliensis

EIN3/EILs are key regulators in ET signaling pathway. In this work, 4 members of EIN3/EILs of Hevea brasiliensis (HbEIN3/EILs) showed interaction with two F box proteins, HbEBF1 and HbEBF2. HbEIN3 located in nucleus and exhibited strong transcriptional activity. HbEIN3 was induced by ET treatment in C-serum, but not in B-serum of latex. HbEIN3/EILs bound to G-box cis-element. To globally search the potential targets of HbEIN3/EILs, genomic sequences of H. brasiliensis was re-annotated and an HCES (Hevea Cis-Elements Scanning) program was developed ( www.h-brasiliensis.com ). HCES scanning results showed that ET- and JA- responsive cis-elements distribute overlapping in gene promoters. 3146 genes containing G-box in promoters are potential targets of HbEIN3, including 41 genes involved in biosynthesis and drainage of latex, of which 7 rate-limiting genes of latex production were regulated by both ET and JA, suggesting that ET and JA signaling pathways coordinated the latex biosynthesis and drainage in H. brasiliensis. Abbreviations: ABRE: ABA responsive elements; bHLH: basic helix-loop-helix; COG: Orthologous Groups; DRE: dehydration response element; ERE: ethylene responsive element; ET: Ethylene; GO: Gene Ontology; HCES: Hevea Cis-Elements Scanning; JA: jasmonates; JRE: Jasmonate-responsive element; KEGG: Kyoto Encyclopedia of Genes and Genomes; NR: non-redundant database; PLACE: Plant Cis-acting Regulatory DNA Elements; qRT-PCR: quantitative real-time RT-PCR.

Involvement of HbMC1-mediated cell death in tapping panel dryness of rubber tree (Hevea brasiliensis)

Tapping panel dryness (TPD) causes a significant reduction in the latex yield of rubber tree (Hevea brasiliensis Muell. Arg.). It is reported that TPD is a typical programmed cell death (PCD) process. Although PCD plays a vital role in TPD occurrence, there is a lack of detailed and systematic study. Metacaspases are key regulators of diverse PCD in plants. Based on our previous result that HbMC1 was associated with TPD, we further elucidate the roles of HbMC1 on rubber tree TPD in this study. HbMC1 was up-regulated by TPD-inducing factors including wounding, ethephon and H2O2. Moreover, the expression level of HbMC1 was increased along with TPD severity in rubber tree, suggesting a positive correlation between HbMC1 expression and TPD severity. To investigate its biological function, HbMC1 was overexpressed in yeast (Saccharomyces cerevisiae) and tobacco (Nicotiana benthamiana). Transgenic yeast and tobacco overexpressing HbMC1 showed growth retardation compared with controls under H2O2-induced oxidative stress. In addition, overexpression of HbMC1 in yeast and tobacco reduced cell survival after high-concentration H2O2 treatment and resulted in enhanced H2O2-induced leaf cell death, respectively. A total of 11 proteins, rbcL, TM9SF2-like, COX3, ATP9, DRP, HbREF/Hevb1, MSSP2-like, SRC2, GATL8, CIPK14-like and STK, were identified and confirmed to interact with HbMC1 by yeast two-hybrid screening and co-transformation in yeast. The 11 proteins mentioned above are associated with many biological processes, including rubber biosynthesis, stress response, autophagy, carbohydrate metabolism, signal transduction, etc. Taken together, our results suggest that HbMC1-mediated PCD plays an important role in rubber tree TPD, and the identified HbMC1-interacting proteins provide valuable information for further understanding the molecular mechanism of HbMC1-mediated TPD in rubber tree.

Identification and Analysis of a CPYC-Type Glutaredoxin Associated with Stress Response in Rubber Trees

Glutaredoxins (GRXs) are a class of small oxidoreductases which modulate various biological processes in plants. Here, we isolated a GRX gene from the rubber tree (Hevea brasiliensis Müll. Arg.), named as HbSRGRX1, which encoded 107 amino acid residues with a CPYC active site. Phylogenetic analysis displayed that HbSRGRX1 was more correlated with GRXs from Manihot esculenta Crantz. and Ricinus communis L. HbSRGRX1 was localized in the nuclei of tobacco cells, and its transcripts were preferentially expressed in male flowers and in the high-yield variety Reyan 7-33-97 with strong resistance against cold. The expression levels of HbSRGRX1 significantly decreased in tapping panel dryness (TPD) trees. Furthermore, HbSRGRX1 was regulated by wounding, hydrogen peroxide (H2O2), and multiple hormones. Altogether, these results suggest important roles of HbSRGRX1 in plant development and defense response to TPD and multiple stresses.

Over-expression of 3-hydroxy-3- methylglutaryl-coenzyme A reductase 1 (hmgr1) gene under super-promoter for enhanced latex biosynthesis in rubber tree (Hevea brasiliensis Muell. Arg.)

Natural rubber (cis-1, 4-polyisoprene) is being produced from bark laticifer cells of Hevea brasiliensis and the popular high latex yielding Indian rubber clones are easily prone to onset of tapping panel dryness syndrome (TPD) which is considered as a physiological syndrome affecting latex production either partially or completely. This report describes an efficient protocol for development of transgenic rubber plants by over-expression of 3-hydroxy 3-methylglutaryl Co-enzyme A reductase 1 (hmgr1) gene which is considered as rate limiting factor for latex biosynthesis via Agrobacterium-mediated transformation. The pBIB plasmid vector containing hmgr1 gene cloned under the control of a super-promoter was used for genetic transformation using embryogenic callus. Putatively transgenic cell lines were obtained on selection medium and produced plantlets with 44% regeneration efficiency. Transgene integration was confirmed by PCR amplification of 1.8 kb hmgr1 and 0.6 kb hpt genes from all putatively transformed callus lines as well as transgenic plants. Southern blot analysis showed the stable integration and presence of transgene in the transgenic plants. Over expression of hmgr1 transgene was determined by Northern blot hybridization, semi-quantitative PCR and real-time PCR (qRT-PCR) analysis. Accumulation of hmgr1 mRNA transcripts was more abundant in transgenic plants than control. Increased level of photosynthetic pigments, protein contents and HMGR enzyme activity was also noticed in transgenic plants over control. Interestingly, the latex yield was significantly enhanced in all transgenic plants compared to the control. The qRT-PCR results exhibit that the hmgr1 mRNA transcript levels was 160-fold more abundance in transgenic plants over untransformed control. These results altogether suggest that there is a positive correlation between latex yield and accumulation of mRNA transcripts level as well as HMGR enzyme activity in transgenic rubber plants. It is presumed that there is a possibility for enhanced level of latex biosynthesis in transgenic plants as the level of mRNA transcripts and HMGR enzyme activity is directly correlated with latex yield in rubber tree. Further, the present results clearly suggest that the quantification of HMGR enzyme activity in young seedlings will be highly beneficial for early selection of high latex yielding plants in rubber breeding programs.

Transcriptome analysis in Hevea brasiliensis latex revealed changes in hormone signalling pathways during ethephon stimulation and consequent Tapping Panel Dryness

Tapping Panel Dryness (TPD) affects latex production in Hevea brasiliensis. This physiological syndrome involves the agglutination of rubber particles, which leads to partial or complete cessation of latex flow. Latex harvesting consists in tapping soft bark. Ethephon can be applied to stimulate latex flow and its regeneration in laticifers. Several studies have reported transcriptome changes in bark tissues. This study is the first report on deep RNA sequencing of latex to compare the effect of ethephon stimulation and TPD severity. Trees were carefully selected for paired-end sequencing using an Illumina HiSeq 2000. In all, 43 to 60 million reads were sequenced for each treatment in three biological replicates (slight TPD trees without ethephon stimulation, and slight and severe TPD trees with ethephon treatment). Differentially expressed genes were identified and annotated, giving 8,111 and 728 in response to ethephon in slight TPD trees and in ethephon-induced severe TPD trees, respectively. A biological network of responses to ethephon and TPD highlighted the major influence of metabolic processes and the response to stimulus, especially wounding and jasmonate depression in TPD-affected trees induced by ethephon stimulation.

Transcriptome-Wide Identification and Characterization of MYB Transcription Factor Genes in the Laticifer Cells of Hevea brasiliensis

MYB transcription factors hold vital roles in the regulation of plant secondary metabolic pathways. Laticifers in rubber trees (Hevea brasiliensis) are of primary importance in natural rubber production because natural rubber is formed and stored within these structures. To understand the role of MYB transcription factors in the specialized cells, we identified 44 MYB genes (named HblMYB1 to HblMYB44) by using our previously obtained transcriptome database of rubber tree laticifer cells and the public rubber tree genome database. Expression profiles showed that five MYB genes were highly expressed in the laticifers. HblMYB19 and HblMYB44 were selected for further study. HblMYB19 and HblMYB44 bound the promoters of HbFDPS1, HbSRPP, and HRT1 in yeast. Furthermore, the transient overexpression of HblMYB19 and HblMYB44 in tobacco plants significantly increased the activity of the promoters of HbFDPS1, HbSRPP, and HRT1. Basing on this information, we proposed that HblMYB19 and HblMYB44 are the regulators of HbFDPS1, HbSRPP, and HRT1, which are involved in the biosynthesis pathway of natural rubber.

Survey of the rubber tree genome reveals a high number of cysteine protease-encoding genes homologous to Arabidopsis SAG12

Arabidopsis thaliana SAG12, a senescence-specific gene encoding a cysteine protease, is widely used as a molecular marker for the study of leaf senescence. To date, its potential orthologues have been isolated from several plant species such as Brassica napus and Nicotiana tabacum. However, little information is available in rubber tree (Hevea brasiliensis), a rubber-producing plant of the Euphorbiaceae family. This study presents the identification of SAG12-like genes from the rubber tree genome. Results showed that an unexpected high number of 17 rubber orthologues with a single intron were found, contrasting the single copy with two introns in Arabidopsis. The gene expansion was also observed in another two Euphorbiaceae plants, castor bean (Ricinus communis) and physic nut (Jatropha curcas), both of which contain 8 orthologues. In accordance with no occurrence of recent whole-genome duplication (WGD) events, most duplicates in castor and physic nut were resulted from tandem duplications. In contrast, the duplicated HbSAG12H genes were derived from tandem duplications as well as the recent WGD. Expression analysis showed that most HbSAG12H genes were lowly expressed in examined tissues except for root and male flower. Furthermore, HbSAG12H1 exhibits a strictly senescence-associated expression pattern in rubber tree leaves, and thus can be used as a marker gene for the study of senescence mechanism in Hevea.

Reactive oxygen species in Hevea brasiliensis latex and relevance to Tapping Panel Dryness

Environmental stress can lead to oxidative stress resulting from an increase in reactive oxygen species (ROS) and involves redox adjustments. Natural rubber is synthesized in laticifers, which is a non-photosynthetic tissue particularly prone to oxidative stress. This paper reviews the current state of knowledge on the ROS production and ROS-scavenging systems in laticifers. These regulations have been the subject of intense research into a physiological syndrome, called Tapping Panel Dryness (TPD), affecting latex production in Hevea brasiliensis. In order to prevent TPD occurrence, monitoring thiol content appeared to be a crucial factor of latex diagnosis. Thiols, ascorbate and γ-tocotrienol are the major antioxidants in latex. They are involved in membrane protection from ROS and likely have an effect on the quality of raw rubber. Some transcription factors might play a role in the redox regulatory network in Hevea, in particular ethylene response factors, which have been the most intensively studied given the role of ethylene on rubber production. Current challenges for rubber research and development with regard to redox systems will involve improving antioxidant capacity using natural genetic variability.

Roles of the Translationally Controlled Tumor Protein (TCTP) in Plant Development

The Translationally Controlled Tumor Protein (TCTP) is a conserved protein which expression was associated with several biochemical and cellular functions. Loss-of-function mutants are lethal both in animals and in plants, making the identification of its exact role difficult. Recent data using the model plant Arabidopsis thaliana provided the first viable adult knockout for TCTP and helped addressing the biological role of TCTP during organ development and the functional conservation between plants and animals. This chapter summarizes our up to date knowledge about the role of TCTP in plants and discuss about conserved functions and mechanisms between plants and animals.

Genome-wide identification and expression analysis of the metacaspase gene family in Hevea brasiliensis

Metacaspases, a family of cysteine proteases, have been suggested to play important roles in programmed cell death (PCD) during plant development and stress responses. To date, no systematic characterization of this gene family has been reported in rubber tree (Hevea brasiliensis). In the present study, nine metacaspase genes, designated as HbMC1 to HbMC9, were identified from whole-genome sequence of rubber tree. Multiple sequence alignment and phylogenetic analyses suggested that these genes were divided into two types: type I (HbMC1-HBMC7) and type II (HbMC8 and HbMC9). Gene structure analysis demonstrated that type I and type II HbMCs separately contained four and two introns, indicating the conserved exon-intron organization of HbMCs. Quantitative real-time PCR analysis revealed that HbMCs showed distinct expression patterns in different tissues, suggesting the functional diversity of HbMCs in various tissues during development. Most of the HbMCs were regulated by drought, cold, and salt stress, implying their possible functions in regulating abiotic stress-induced cell death. Of the nine HbMCs, HbMC1, HbMC2, HbMC5, and HbMC8 displayed a significantly higher relative transcript accumulation in barks of tapping panel dryness (TPD) trees compared with healthy trees. In addition, the four genes were up-regulated by ethephon (ET) and methyl jasmonate (MeJA), indicating their potential involvement in TPD resulting from ET- or JA-induced PCD. In summary, this work provides valuable information for further functional characterization of HbMC genes in rubber tree.

Transcriptome analyses reveal molecular mechanism underlying tapping panel dryness of rubber tree (Hevea brasiliensis)

Tapping panel dryness (TPD) is a serious threat to natural rubber yields from rubber trees, but the molecular mechanisms underlying TPD remain poorly understood. To identify TPD-related genes and reveal these molecular mechanisms, we sequenced and compared the transcriptomes of bark between healthy and TPD trees. In total, 57,760 assembled genes were obtained and analyzed in details. In contrast to healthy rubber trees, 5652 and 2485 genes were up- or downregulated, respectively, in TPD trees. The TPD-related genes were significantly enriched in eight GO terms and five KEGG pathways and were closely associated with ROS metabolism, programmed cell death and rubber biosynthesis. Our results suggest that rubber tree TPD is a complex process involving many genes. The observed lower rubber yield from TPD trees might result from lower isopentenyl diphosphate (IPP) available for rubber biosynthesis and from downregulation of the genes in post-IPP steps of rubber biosynthesis pathway. Our results not only extend our understanding of the complex molecular events involved in TPD but also will be useful for developing effective measures to control TPD of rubber trees.

Involvement of Ethylene in the Latex Metabolism and Tapping Panel Dryness of Hevea brasiliensis

Ethephon, an ethylene releaser, is used to stimulate latex production in Hevea brasiliensis. Ethylene induces many functions in latex cells including the production of reactive oxygen species (ROS). The accumulation of ROS is responsible for the coagulation of rubber particles in latex cells, resulting in the partial or complete stoppage of latex flow. This study set out to assess biochemical and histological changes as well as changes in gene expression in latex and phloem tissues from trees grown under various harvesting systems. The Tapping Panel Dryness (TPD) susceptibility of Hevea clones was found to be related to some biochemical parameters, such as low sucrose and high inorganic phosphorus contents. A high tapping frequency and ethephon stimulation induced early TPD occurrence in a high latex metabolism clone and late occurrence in a low latex metabolism clone. TPD-affected trees had smaller number of laticifer vessels compared to healthy trees, suggesting a modification of cambial activity. The differential transcript abundance was observed for twenty-seven candidate genes related to TPD occurrence in latex and phloem tissues for ROS-scavenging, ethylene biosynthesis and signalling genes. The predicted function for some Ethylene Response Factor genes suggested that these candidate genes should play an important role in regulating susceptibility to TPD.

Transcriptome sequencing and analysis of rubber tree (Hevea brasiliensis Muell.) to discover putative genes associated with tapping panel dryness (TPD)

Background

Tapping panel dryness (TPD) involves in the partial or complete cessation of latex flow thus seriously affect latex production in rubber tree (Hevea brasiliensis). Numerous studies have been conducted to define its origin and nature, but the molecular nature and mechanism of TPD occurrence remains unknown. This study is committed to de novo sequencing and comparative analysis of the transcriptomes of healthy (H) and TPD-affected (T) rubber trees to identify the genes and pathways related to the TPD.

Results

Total raw reads of 34,632,012 and 35,913,020 bp were obtained from H and T library, respectively using Illumina Hiseq 2000 sequencing technology. De novo assemblies yielded 141,456 and 169,285 contigs, and 96,070 and 112,243 unigenes from H and T library, respectively.

Among 73597 genes, 22577 genes were identified as differential expressed genes between H and T library via comparative transcript profiling. A majority of genes involved in natural rubber biosynthesis and jasmonate synthesis with most potential relevance in TPD occurrence were found to be differentially expressed.

Conclusions

In TPD-affected trees, the expression of most genes related to the latex biosynthesis and jasmonate synthesis was severely inhibited and is probably the direct cause of the TPD. These new de novo transcriptome data sets provide a significant resource for the discovery of genes related to TPD and improve our understanding of the occurrence and maintainace of TPD.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1562-9) contains supplementary material, which is available to authorized users.

Molecular cloning and expression of a novel MYB transcription factor gene in rubber tree

MYB family proteins regulate a variety of cellular processes in plants. Tapping panel dryness (TPD) in rubber tree (Hevea brasiliensis Muell. Arg.) affects latex biosynthesis and causes serious losses to rubber producers. In this study, a novel SANT/MYB transcription factor gene down-regulated in TPD rubber tree, named as HbSM1, was isolated from rubber tree. The complete HbSM1 open reading frame (ORF) was 948 bp in length. The deduced HbSM1 protein is 315 amino acids. HbSM1 belonged to 1RMYB subfamily with a single SANT domain. Sequence alignment revealed that HbSM1 had high homology with MYB members from Ricinus communis and Manihot esculenta, with 72 and 78 % identity, respectively. Moreover, HbSM1 shared 56 % identity with Glycine max GmMYB176. Phylogenetic analysis revealed that HbSM1, GmMYB176, rice OsMYBS2, and OsMYBS3 fell into the same cluster with 93 % bootstrap support value. Comparing expression among different tissues demonstrated that HbSM1 was ubiquitously expressed in all tissues, but it appeared to be preferentially expressed in leaf and latex. Furthermore, HbSM1 transcripts were significantly induced by various phytohormones (including gibberellic acid, ethephon, methyl jasmonate, salicylic acid, and abscisic acid) and wounding treatments. These results suggested that HbSM1 might play multiple roles in plant development via different phytohormones signaling pathways.

The small RNA profile in latex from Hevea brasiliensis trees is affected by tapping panel dryness

Natural rubber is harvested by tapping Hevea brasiliensis (Willd. ex A. Juss.) Müll. Arg. Harvesting stress can lead to tapping panel dryness (TPD). MicroRNAs (miRNAs) are induced by abiotic stress and regulate gene expression by targeting the cleavage or translational inhibition of target messenger RNAs. This study set out to sequence miRNAs expressed in latex cells and to identify TPD-related putative targets. Deep sequencing of small RNAs was carried out on latex from trees affected by TPD using Solexa technology. The most abundant small RNA class size was 21 nucleotides for TPD trees compared with 24 nucleotides in healthy trees. By combining the LeARN pipeline, data from the Plant MicroRNA database and Hevea EST sequences, we identified 19 additional conserved and four putative species-specific miRNA families not found in previous studies on rubber. The relative transcript abundance of the Hbpre-MIR159b gene increased with TPD. This study revealed a small RNA-specific signature of TPD-affected trees. Both RNA degradation and a shift in miRNA biogenesis are suggested to explain the general decline in small RNAs and, particularly, in miRNAs.

The function of Rad6 gene in Hevea brasiliensis extends beyond DNA repair

The Rad6 gene of Saccharomyces cerevisiae encodes an ubiquitin-conjugating enzyme (E2) which is required for DNA repair, damage-induced mutagenesis, sporulation, etc. In this study, one Rad6 homolog, designated HbRad6, was cloned in rubber tree (Hevea brasiliensis). The putative protein sequence of HbRad6 contains 152 amino acids, a conserved UBC domain, and a conserved active-site cysteine in the UBC domain, which is required for E2 enzymes catalytic activity. HbRad6 shared high similarity with Rad6 from other species. It shared the highest similarity with rice OsRad6 and Arabidopsis thaliana AtUBC2 with 96.05% identical residues, and 63.16% sequence identity with yeast Rad6 (excluding the acidic tail). Comparing expression among different Hevea tissues demonstrated that HbRad6 was ubiquitously expressed in all tissues, but it revealed a preferential expression in the latex. Furthermore, HbRad6 expression was markedly induced by DNA-damaging agent H2O2, the latex stimulator ethephon (ET), and methyl jasmonate (MeJA), while NaCl and wounding treatments had relatively minor effect upon its expression. Genetic complementation experiment revealed that HbRad6 had minor effects on the complementation of the UV sensitivity of yeast rad6 null mutant, indicating that the Hevea Rad6 protein may partially suppress the UV sensitivity of the yeast rad6 mutant. These results suggested that HbRad6 was a multifunction gene involved in DNA damage repair, hormones and stress responses in rubber tree.

Molecular cloning, expression profiles and characterization of a novel translationally controlled tumor protein in rubber tree (Hevea brasiliensis)

The translationally controlled tumor protein (TCTP) is a multi-functioning protein that carries out vital roles in various life processes. In this study, a new TCTP gene, designated as HbTCTP1, was isolated in Hevea brasiliensis. The full-length complementary DNA (cDNA) of HbTCTP1 contained a maximum open reading frame (ORF) of 507base pair (bp) encoding 168 amino acids. The sequence comparison showed that the deduced HbTCTP1 indicated high identities to plant TCTP proteins, and clustered in the dicot cluster of plant TCTPs. Although HbTCTP1 and human TCTP proteins did not parallel in overall sequence similarity, they indicated highly similar 3D structures with a nearly identical spatial organization of α-helices, β-sheets, and coil regions. Real time reverse-transcription PCR (RT-PCR) analyses showed that HbTCTP1 was expressed throughout different tissues and developmental stages of leaves. Besides being related to tapping panel dryness (TPD), the HbTCTP1 transcripts were regulated by various treatments, including drought, low temperature, high salt, ethrel (ET), wounding, H2O2, and methyl jasmonate (Me-JA) treatments. The recombinant HbTCTP1 fusion protein was shown to protect supercoiled plasmid DNA from damages induced by metal-catalyzed generation of reactive oxygen species. The (45)Ca(2+)-overlay assay showed that HbTCTP1 was a calcium-binding protein. Our results are greatly helpful in understanding the molecular characterization and expression profiles of HbTCTP1, and lay the foundation for further analyzing the function of HbTCTP1 in rubber tree.

Over-expression of a cytosolic isoform of the HbCuZnSOD gene in Hevea brasiliensis changes its response to a water deficit

Hevea brasiliensis is the main commercial source of natural rubber. Reactive oxygen species (ROS) scavenging systems are involved in various biotic and abiotic stresses. Genetic engineering was undertaken to study the strengthening of plant defences by antioxidants. To that end, Hevea transgenic plant lines over-expressing a Hevea brasiliensis cytosolic HbCuZnSOD gene were successfully established and regenerated. Over-expression of the HbCuZnSOD gene was not clearly related to an increase in SOD activity in plant leaves. The impact of HbCuZnSOD gene over-expression in somatic embryogenesis and in plant development are presented and discussed. The water deficit tolerance of two HbCuZnSOD over-expressing lines was evaluated. The physiological parameters of transgenic plantlets subjected to a water deficit suggested that plants from line TS4T8An displayed lower stomatal conductance and a higher proline content. Over-expression of the HbCuZnSOD gene and activation of all ROS-scavenging enzymes also suggested that protection against ROS was more efficient in the TS4T8An transgenic line.

Transcriptome and gene expression analysis in cold-acclimated guayule (Parthenium argentatum) rubber-producing tissue

Natural rubber biosynthesis in guayule (Parthenium argentatum Gray) is associated with moderately cold night temperatures. To begin to dissect the molecular events triggered by cold temperatures that govern rubber synthesis induction in guayule, the transcriptome of bark tissue, where rubber is produced, was investigated. A total of 11,748 quality expressed sequence tags (ESTs) were obtained. The vast majority of ESTs encoded proteins that are similar to stress-related proteins, whereas those encoding rubber biosynthesis-related proteins comprised just over one percent of the ESTs. Sequence information derived from the ESTs was used to design primers for quantitative analysis of the expression of genes that encode selected enzymes and proteins with potential impact on rubber biosynthesis in field-grown guayule plants, including 3-hydroxy-3-methylglutaryl-CoA synthase, 3-hydroxy-3-methylglutaryl-CoA reductase, farnesyl pyrophosphate synthase, squalene synthase, small rubber particle protein, allene oxide synthase, and cis-prenyl transferase. Gene expression was studied for field-grown plants during the normal course of seasonal variation in temperature (monthly average maximum 41.7 °C to minimum 0 °C, from November 2005 through March 2007) and rubber transferase enzymatic activity was also evaluated. Levels of gene expression did not correlate with air temperatures nor with rubber transferase activity. Interestingly, a sudden increase in night temperature 10 days before harvest took place in advance of the highest CPT gene expression level.

Development and characterization of intron-flanking EST-PCR markers in rubber tree (Hevea brasiliensis Muell. Arg.)

With a long-term goal of constructing a linkage map enriched with gene-specific markers in rubber tree (Hevea brasiliensis Muell. Arg.), we utilized rubber tree ESTs associated with tapping panel dryness (TPD) to develop intron-flanking PCR markers. After downloading and assembling the rubber tree ESTs associated with TPD, we predicted the exon/exon junction sites (E/E) by aligning rubber tree transcripts with the genomic sequences of castor bean (Ricinus communis L.). Based on the predicted E/E, the primers flanking intron(s) and no intron were designed. Compared with the markers designed by conventional method, the PCR success rate of the markers designed with the predicted E/E increased 28-30%, whereas the polymorphism rate of intron-flanking EST-PCR markers was approximately 3.43-fold increase. Therefore, the intron-flanking marker was more polymorphism-generating efficient than the markers designed by conventional methods. In addition, analyzing the polymorphic information content (PIC) among Hevea germplasm showed that the polymorphism of wild rubber tree accessions was higher than one of cultivated rubber tree clones and Hevea species. This study enriches the categories and numbers of molecular markers in rubber tree, and the markers developed in this research will have a wide application in DNA fingerprinting, marker-assisted selection and genetic mapping in rubber tree. This research also indicates that it is possible to develop intron-flanking EST-PCR markers of rubber tree with castor bean genome as reference sequences, which provides new insights into developing intron-flanking EST-PCR markers for rubber tree or other plant species without genomic information.

De novo assembly and characterization of bark transcriptome using Illumina sequencing and development of EST-SSR markers in rubber tree (Hevea brasiliensis Muell. Arg.)

Background

In rubber tree, bark is one of important agricultural and biological organs. However, the molecular mechanism involved in the bark formation and development in rubber tree remains largely unknown, which is at least partially due to lack of bark transcriptomic and genomic information. Therefore, it is necessary to carried out high-throughput transcriptome sequencing of rubber tree bark to generate enormous transcript sequences for the functional characterization and molecular marker development.

Results

In this study, more than 30 million sequencing reads were generated using Illumina paired-end sequencing technology. In total, 22,756 unigenes with an average length of 485 bp were obtained with de novo assembly. The similarity search indicated that 16,520 and 12,558 unigenes showed significant similarities to known proteins from NCBI non-redundant and Swissprot protein databases, respectively. Among these annotated unigenes, 6,867 and 5,559 unigenes were separately assigned to Gene Ontology (GO) and Clusters of Orthologous Group (COG). When 22,756 unigenes searched against the Kyoto Encyclopedia of Genes and Genomes Pathway (KEGG) database, 12,097 unigenes were assigned to 5 main categories including 123 KEGG pathways. Among the main KEGG categories, metabolism was the biggest category (9,043, 74.75%), suggesting the active metabolic processes in rubber tree bark. In addition, a total of 39,257 EST-SSRs were identified from 22,756 unigenes, and the characterizations of EST-SSRs were further analyzed in rubber tree. 110 potential marker sites were randomly selected to validate the assembly quality and develop EST-SSR markers. Among 13 Hevea germplasms, PCR success rate and polymorphism rate of 110 markers were separately 96.36% and 55.45% in this study.

Conclusion

By assembling and analyzing de novo transcriptome sequencing data, we reported the comprehensive functional characterization of rubber tree bark. This research generated a substantial fraction of rubber tree transcriptome sequences, which were very useful resources for gene annotation and discovery, molecular markers development, genome assembly and annotation, and microarrays development in rubber tree. The EST-SSR markers identified and developed in this study will facilitate marker-assisted selection breeding in rubber tree. Moreover, this study also supported that transcriptome analysis based on Illumina paired-end sequencing is a powerful tool for transcriptome characterization and molecular marker development in non-model species, especially those with large and complex genomes.

HbMyb1, a Myb transcription factor from Hevea brasiliensis, suppresses stress induced cell death in transgenic tobacco

Tapping panel dryness (TPD) is a complex physiological syndrome found widely in rubber tree (Hevea brasiliensis) plantations that causes severe yield loss in natural rubber-producing countries. In an earlier study, we confirmed that there is a negative correlation between HbMyb1 expression and TPD severity. To further investigate the function of HbMyb1 in TPD, HbMyb1 was over-expressed in tobacco controlled by a CaMV 35S promoter. In transgenic plants expressing HbMyb1, cell death induced by UV-B irradiation, paraquat and the hypersensitive reaction to necrotrophic fungal infection (Botrytis cinerea) was suppressed with a close correlation between HbMyb1 protein levels and the extent of suppression. In addition the nuclear condensation and degradation were observed in laticifer cells of TPD trees, while the nucleus of laticifer cells of healthy trees was morphologically normal. On the basis of the results described above, we propose that HbMyb1 maybe suppress stress induced cell death in rubber trees.

Gene expression pattern in response to wounding, methyl jasmonate and ethylene in the bark of Hevea brasiliensis

Natural rubber production in Hevea brasiliensis is determined by both tapping and ethephon frequencies. It is affected by a complex physiological disorder called tapping panel dryness. This syndrome is likely to be induced by environmental and latex harvesting stresses. Defence responses, including rubber biosynthesis, are dramatically mediated by wounding, jasmonate and ethylene (ET), among other factors. Using real-time RT-PCR, the effects of wounding, methyl jasmonate (MeJA) and ET on the relative transcript abundance of a set of 25 genes involved in their signalling and metabolic pathways were studied in the bark of 3-month-old epicormic shoots. Temporal regulation was found for 9 out of 25 genes. Wounding treatment regulated the transcript abundance of 10 genes. Wounding-specific regulation was noted for the HbMAPK, HbBTF3b, HbCAS1, HbLTPP and HbPLD genes. MeJA treatment regulated the transcript abundance of nine genes. Of these, the HbMYB, HbCAS2, HbCIPK and HbChi genes were shown to be specifically MeJA inducible. ET response was accompanied by regulation of the transcript abundance of eight genes, and six genes, HbETR2, HbEIN2, HbEIN3, HbCaM, HbPIP1 and HbQM, were specifically regulated by ET treatment. Additionally, the transcript level of the HbGP and HbACR genes was enhanced by all three treatments simultaneously. Overall, a large number of genes were found to be regulated 4 h after the treatments were applied. This study nevertheless revealed some jasmonic acid-independent wound signalling pathways in H. brasiliensis, provided a general characterization of signalling pathways and will serve as a new base from which to launch advanced studies of the network of pathways operating in H. brasiliensis.

Molecular characterization of a thioredoxin h gene (HbTRX1) from Hevea brasiliensis showing differential expression in latex between self-rooting juvenile clones and donor clones

The cDNA code of thioredoxin h, designated as HbTRX1, was isolated from Hevea brasiliensis by rapid amplification of cDNA ends. HbTRX1 contained a 542-bp open reading frame encoding 123 amino acids. The deduced HbTRX1 protein showing high identity to thioredoxin h of other plant species was predicted to possess the conserved catalytic site WCXPC. Semiquantitative reverse transcription-polymerase chain reaction analysis revealed that HbTRX1 was constitutively expressed in all tested tissues. HbTRX1 transcripts accumulated at relatively low levels in the flower, somatic embryo, and leaves, while HbTRX1 transcripts accumulated at relatively high levels in the callus and latex. The HbTRX1 transcript was expressed at different levels, with higher levels in self-rooting juvenile clones than in their donor clones. HbTRX1 was expressed in Escherichia coli, and its activity was demonstrated using the dithiothreitol-dependent insulin assay. This work provides a basis for studying the biological function of thioredoxin h in rubber tree.

Identification and characterization of genes associated with tapping panel dryness from Hevea brasiliensis latex using suppression subtractive hybridization

BACKGROUND

Tapping panel dryness (TPD) is one of the most serious threats to natural rubber production. Although a great deal of effort has been made to study TPD in rubber tree, the molecular mechanisms underlying TPD remain poorly understood. Identification and systematical analyses of the genes associated with TPD are the prerequisites for elucidating the molecular mechanisms involved in TPD. The present study is undertaken to generate information about the genes related to TPD in rubber tree.

RESULTS

To identify the genes related to TPD in rubber tree, forward and reverse cDNA libraries from the latex of healthy and TPD trees were constructed using suppression subtractive hybridization (SSH) method. Among the 1106 clones obtained from the two cDNA libraries, 822 clones showed differential expression in two libraries by reverse Northern blot analyses. Sequence analyses indicated that the 822 clones represented 237 unique genes; and most of them have not been reported to be associated with TPD in rubber tree. The expression patterns of 20 differentially expressed genes were further investigated to validate the SSH data by reverse transcription PCR (RT-PCR) and real-time PCR analysis. According to the Gene Ontology convention, 237 unique genes were classified into 10 functional groups, such as stress/defense response, protein metabolism, transcription and post-transcription, rubber biosynthesis, etc. Among the genes with known function, the genes preferentially expressed were associated with stress/defense response in the reverse library, whereas metabolism and energy in the forward one.

CONCLUSIONS

The genes associated with TPD were identified by SSH method in this research. Systematic analyses of the genes related to TPD suggest that the production and scavenging of reactive oxygen species (ROS), ubiquitin proteasome pathway, programmed cell death and rubber biosynthesis might play important roles in TPD. Therefore, our results not only enrich information about the genes related to TPD, but also provide new insights into understanding the TPD process in rubber tree.

Genes induced in response to mercury-ion-exposure in heavy metal hyperaccumulator Sesbania drummondii

Sesbania drummondii plants have been recognized as a potential mercury (Hg) hyperaccumulator. To identify genes modulated by Hg, two suppressive subtraction hybridization (SSH) cDNA libraries (forward and reverse) were constructed. A total of 348 differentially expressed clones were isolated and 95 of them were identified as Hg responsive. Reverse Northern results showed that 31 clones from forward library were down-regulated and 64 clones from reverse library were up-regulated in Hg-treated plants. Sixty-seven of them showed high homology to genes with known or putative function, and 28 were uncharacterized genes. Two full-length cDNAs coding for a putative metallothionein type 2 protein (SdMT2) and an auxin responsive protein (SdARP) were isolated and characterized. The expression levels of SdMT2 and SdARP increased 3- and 5-fold, respectively. Results suggest that up-regulated expression of SdARP may contribute to the survival of Sesbania plants under mercury stress, whereas SdMT2 is likely to be involved in alleviation of Hg toxicity. The possible correlation between gene expression and heavy metal tolerance of Sesbania plants is discussed.

Molecular cloning and characterization of phosphate (Pi) responsive genes in Gulf ryegrass (Lolium multiflorum L.): a Pi hyperaccumulator

Gulf annual ryegrass has been identified as potential Pi hyperaccumulator, however the molecular mechanism remains largely unknown. A suppression subtractive hybridization (SSH) analysis was used to evaluate the phosphate (Pi) responsive genome expression pattern changes in Gulf annual ryegrass (Lolium multiflorum L.). Differential screening identified 384 putative Pi-starvation induced cDNAs. Bioinformatic analysis revealed that 116 cDNAs are nonredundant unigenes of which 108 exhibited high similarities with Genbank entries. The differential expression patterns of 13 cDNAs, representing diverse functional categories, were confirmed by RNA gel blot analysis. Further, detailed molecular analysis of three genes (LmPAP1, LmIPS1 and LmIDS1) was carried out by cloning and characterization of full-length cDNAs. LmPAP1 is 1,414 bp in length with an open reading frame (ORF) of 1,188 bp capable of encoding an N-terminal signal peptide of 26 amino acids. LmIPS1 gene is a member of TPSI1/Mt4 family that contains 3 short ORFs. The cDNA of LmIDS1 is 346 bp in size including a single ORF of 222 nucleotides that encodes 74 amino acid proteins, exhibiting homology with IDS1 with similarity to type 2 metallothionein like protein. In our preliminary screening of different genotypes of annual ryegrass for hyperaccumulation of Pi in their shoots, Gulf and Urugrary showed significant differences with values of 1.0% and 0.7%, respectively. Since it is logical to assume a plausible correlation that may exist between Pi-accumulation in the shoots and the expression of Pi-responsive genes, the expression of LmPAP1, LmIPS1 and LmIDS1 was evaluated in these two genotypes grown under different Pi regimes. Although there was a significant induction of these genes in both the genotypes grown under Pi-deprived condition, the abundance of LmPAP1 transcripts was relatively higher in the Gulf genotype as compared to that in the Urugrary genotype. A similar trend was observed in qRT-PCR data of other tested genotypes of annual ryegrasses. This suggests the potential role of LmPAP1 in accumulation of Pi in Gulf grass. In addition, Gulf grass genotype revealed higher levels of total P, (33)Pi uptake, and APase activity as compared to Urugrary. Together, these results suggest that the Gulf ryegrass has evolved mechanisms to acquire and hyperaccumulate more Pi under different Pi regimes by activating multiple Pi acquisition and mobilization mechanisms.

Molecular cloning and characterization of phosphorus starvation responsive genes in common bean (Phaseolus vulgaris L.)

Common bean (Phaseolus vulgaris L.) is one of the most important food legumes in the world and its production is limited by low phosphate (Pi) availability in many arable soils. To gain better insight into the molecular mechanisms by which common bean adapts to low Pi availability, we generated a suppression subtractive cDNA library to identify genes involved in P starvation responses. Over 240 putative Pi starvation-responsive genes were identified. The identified clones were sequenced and BLASTx/BLASTn analysis revealed an array of 82 genes showing a high degree of sequence homology to known and unknown proteins in the database. Transcript abundance of seven genes representing different functional categories was examined by Northern blot analysis. Six genes were strongly induced/enhanced under Pi deficiency confirming the results of SSH. Full length cDNAs for three genes, representing PvIDS4-like, PvPS2, and PvPT1 were cloned and characterized. The open reading frame (ORF) of PvIDS4-like encodes a 281-amino acid protein, containing a SPX domain. The ORF of PvPS2 gene encodes a 271-amino acid protein coding for a putative phosphatase. The PvPT1 encodes a 531-amino acid protein exhibiting high homology with high affinity Pi transporters. Expression patterns of these three genes in relation to Pi availability were evaluated with two contrasting genotypes (P-inefficient Dor364 and P-efficient G19833). Both Northern and RT-PCR results showed enhanced accumulation of phosphate transporters and phosphatases in P-efficient genotype, implying that in addition to modified root morphology and architecture, increased P transport and phosphatases activity might contribute to efficient Pi acquisition and translocation in G19833 common bean genotype under limited Pi conditions.