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.

Cloning and Expression Analysis of HbPR-1b and HbPR-3 in Hevea brasiliensis During Inoculation with Rigidoporus microporus

BACKGROUND AND OBJECTIVE

Pathogenesis-related (PR) proteins are dramatically accumulated after pathogen infection. Innate defense response through increasing PR-proteins is important for rubber rootstock selection that is tolerant to the white root disease caused by Rigidoporus microporus. This study was aimed to investigate the expression levels of PR-1 and PR-3 genes in tolerant (PB5/51) and susceptible (BPM24 and RRIM600) rubber clones after R. microporus infection.

MATERIALS AND METHODS

The mRNA of HbPR-1b and HbPR-3 was isolated and characterized from rubber leaves. Gene expression levels of HbPR-1b and HbPR-3 were compared among three rubber clones (PB5/51, BPM24 and RRIM600) after R. microporus infection at 0, 12, 24, 48, 72 and 96 h using quantitative real-time PCR. The relative transcript abundances between inoculated and control plants were compared using the means of gene expression between time points and by Tukey's HSD test. A probability value (p<0.05) was used to determine the significance of difference between time points.

RESULTS

The open reading frame of HbPR-1b is 492 bp with deduced 163 amino acid residues and the phylogenetic analysis showed it shared significant evolutionary history and clustering into group I of PR-protein. Moreover, the partial HbPR-3 was isolated with 390 bp. Gene expression levels of HbPR-1b and HbPR-3 showed marked differences in both transcripts depending on the rubber clones. Two genes demonstrated up-regulation of both tolerance and susceptibility in response to attack by R. microporus. The highest expression levels were found in seedlings of PB5/51 after inoculation. In RRIM600, low expression levels of HbPR-1b and HbPR-3 were initially observed but gradually increased at 24 h post inoculation. The transcription profile of HbPR-1b was stable expression in BPM24.

CONCLUSION

The results demonstrated that the level ofHbPR-1b and HbPR-3 transcription can distinguish between tolerant and susceptible clones. The candidate defense genes to the white root disease were observed in PB5/51 seedlings, particularly HbPR-1b.

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.

Concentrations, Stability, and Isolation of the Furan Fatty Acid 9-(3-Methyl-5-pentylfuran-2-yl)-nonanoic Acid from Disposable Latex Gloves

Because of their antioxidant properties, furan fatty acids (furan-FAs) are valuable minor compounds with a widespread occurrence in all living matter. Unfortunately, pure standards are not readily available, because they usually contribute only 1% to the lipid fraction. A known exception of this is the milky fluid of Hevea brasiliensis, commonly known as latex, in which the furan-FA 9-(3-methyl-5-pentylfuran-2-yl)-nonanoic acid (9M5) contributes about 90% to the triacylglycerides. In this study, we investigated the content of 9M5 in 30 different disposable latex gloves, which ranged from 0.7 to 8.2 mg/g of glove. The light degradability of 9M5 in latex gloves was investigated, and different amounts of 9M5 in disposable latex gloves were attributed to varying exposure time to light. Additionally, over 100 mg of the methyl or ethyl ester of 9M5 (purity of >98%) could be extracted from disposable latex gloves, employing cold extraction and silver ion chromatography. With this method, standards for the quantitation of furan-FAs are obtained easily and rapidly in all laboratories.

Comparative Proteomics of Rubber Latex Revealed Multiple Protein Species of REF/SRPP Family Respond Diversely to Ethylene Stimulation among Different Rubber Tree Clones

Rubber elongation factor (REF) and small rubber particle protein (SRPP) are two key factors for natural rubber biosynthesis. To further understand the roles of these proteins in rubber formation, six different genes for latex abundant REF or SRPP proteins, including REF_138,175,258 and SRPP_117,204,243, were characterized from Hevea brasiliensis Reyan (RY) 7-33-97. Sequence analysis showed that REFs have a variable and long N-terminal, whereas SRPPs have a variable and long C-terminal beyond the REF domain, and REF₂₅₈ has a β subunit of ATPase in its N-terminal. Through two-dimensional electrophoresis (2-DE), each REF/SRPP protein was separated into multiple protein spots on 2-DE gels, indicating they have multiple protein species. The abundance of REF/SRPP proteins was compared between ethylene and control treatments or among rubber tree clones with different levels of latex productivity by analyzing 2-DE gels. The total abundance of each REF/SRPP protein decreased or changed a little upon ethylene stimulation, whereas the abundance of multiple protein species of the same REF/SRPP changed diversely. Among the three rubber tree clones, the abundance of the protein species also differed significantly. Especially, two protein species of REF₁₇₅ or REF₂₅₈ were ethylene-responsive only in the high latex productivity clone RY 8-79 instead of in RY 7-33-97 and PR 107. Some individual protein species were positively related to ethylene stimulation and latex productivity. These results suggested that the specific protein species could be more important than others for rubber production and post-translational modifications might play important roles in rubber biosynthesis.

Cropping system diversification for food production in Mindanao rubber plantations: a rice cultivar mixture and rice intercropped with mungbean

Including food production in non-food systems, such as rubber plantations and biofuel or bioenergy crops, may contribute to household food security. We evaluated the potential for planting rice, mungbean, rice cultivar mixtures, and rice intercropped with mungbean in young rubber plantations in experiments in the Arakan Valley of Mindanao in the Philippines. Rice mixtures consisted of two- or three-row strips of cultivar Dinorado, a cultivar with higher value but lower yield, and high-yielding cultivar UPL Ri-5. Rice and mungbean intercropping treatments consisted of different combinations of two- or three-row strips of rice and mungbean. We used generalized linear mixed models to evaluate the yield of each crop alone and in the mixture or intercropping treatments. We also evaluated a land equivalent ratio for yield, along with weed biomass (where Ageratum conyzoides was particularly abundant), the severity of disease caused by Magnaporthe oryzae and Cochliobolus miyabeanus, and rice bug (Leptocorisa acuta) abundance. We analyzed the yield ranking of each cropping system across site-year combinations to determine mean relative performance and yield stability. When weighted by their relative economic value, UPL Ri-5 had the highest mean performance, but with decreasing performance in low-yielding environments. A rice and mungbean intercropping system had the second highest performance, tied with high-value Dinorado but without decreasing relative performance in low-yielding environments. Rice and mungbean intercropped with rubber have been adopted by farmers in the Arakan Valley.

Identification and subcellular localization analysis of two rubber elongation factor isoforms on Hevea brasiliensis rubber particles

Rubber elongation factor (REF) is the most abundant protein found on the rubber particles or latex from Hevea brasiliensis (the Para rubber tree) and is considered to play important roles in natural rubber (cis-polyisoprene) biosynthesis. 16 BAC (benzyldimethyl-n-hexadecylammonium chloride)/SDS-PAGE separations and mass spectrometric identification had revealed that two REF isoforms shared similar amino acid sequences and common C-terminal sequences. In this study, the gene sequences encoding these two REF isoforms (one is 23.6 kDa in size with 222 amino acid residues and the other is 27.3 kDa in size with 258 amino acid residues) were obtained. Their proteins were relatively enriched by sequential extraction of the rubber particle proteins and separated by 16 BAC/SDS-PAGE. The localization of these isoforms on the surfaces of rubber particles was further verified by western blotting and immunogold electron microscopy, which demonstrated that these two REF isoforms are mainly located on the surfaces of larger rubber particles and that they bind more tightly to rubber particles than the most abundant REF and SRPP (small rubber particle protein).

Function of Hevea brasiliensis NAC1 in dehydration-induced laticifer differentiation and latex biosynthesis

MAIN CONCLUSIONS

HbNAC1 is a transcription factor in rubber plants whose expression is induced by dehydration, leading to latex biosynthesis. Laticifer is a special tissue in Hevea brasiliensis where natural rubber is biosynthesized and accumulated. In young stems of epicormic shoots, the differentiation of secondary laticifers can be induced by wounding, which can be prevented when the wounding site is wrapped. Using this system, differentially expressed genes were screened by suppression subtractive hybridization (SSH) and macroarray analyses. This led to the identification of several dehydration-related genes that could be involved in laticifer differentiation and/or latex biosynthesis, including a NAC transcription factor (termed as HbNAC1). Tissue sections confirmed that local tissue dehydration was a key signal for laticifer differentiation. HbNAC1 was localized at the nucleus and showed strong transcriptional activity in yeast, suggesting that HbNAC1 is a transcription factor. Furthermore, HbNAC1 was found to bind to the cis-element CACG in the promoter region of the gene encoding the small rubber particle protein (SRPP). Transgenic experiments also confirmed that HbNAC1 interacted with the SRPP promoter when co-expressed, and enhanced expression of the reporter gene β-glucuronidase occurred in planta. In addition, overexpression of HbNAC1 in tobacco plants conferred drought tolerance. Together, the data suggest that HbNAC1 might be involved in dehydration-induced laticifer differentiation and latex biosynthesis.

Effect of HbDHN1 and HbDHN2 Genes on Abiotic Stress Responses in Arabidopsis

Dehydrin is a type of late embryogenesis abundant (LEA) protein. The dehydrin genes, HbDHN1 and HbDHN2, in Hevea brasiliensis were previously found to be induced at the wounding site of epicormic shoots, with local tissue dehydration identified as the key signal for laticifer differentiation. However, the exact role of the HbDHNs remains unknown. In this study, HbDHN1 and HbDHN2 expression was examined under multiple abiotic stresses; namely, cold, salt, drought, wounding, abscisic acid (ABA), ethylene (ET), and jasmonic acid (JA) treatment. Although, both HbDHNs were defined as SK2-type dehydrin, they showed different cellular localizations. Overexpression of the HbDHNs in Arabidopsis thaliana further revealed a significant increase in tolerance to salt, drought and osmotic stresses. Increased accumulation of proline and a reduction in electrolyte leakage were also observed under salt and drought stress, and a higher water content was indicated under osmotic stress. The transgenic plants also showed higher activity levels of ascorbate peroxidase (APX), superoxide dismutase (SOD) and catalase, and accumulated less hydrogen peroxide (H₂O₂) and superoxide ([Formula: see text]). Given that reactive oxygen species (ROS) are thought to be a key signal for laticifer differentiation, these findings suggest that HbDHNs act as ROS scavengers, directly or indirectly affecting laticifer differentiation. Both HbDHNs therefore influence physiological processes, improving plant tolerance to multiple abiotic stresses.

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

Calcium‐dependent protein kinases (CDPKs or CPKs) play important roles in various physiological processes of plants, including growth and development, stress responses and hormone signaling. Although the CDPK gene family has been characterized in several model plants, little is known about this gene family in Hevea brasiliensis (the Para rubber tree). Here, we characterize the entire H. brasiliensis CDPK and CDPK‐related kinase (CRK) gene families comprising 30 CDPK genes (HbCPK1 to 30) and nine CRK genes (HbCRK1 to 9). Structure and phylogeny analyses of these CDPK and CRK genes demonstrate evolutionary conservation in these gene families across H. brasiliensis and other plant species. The expression of HbCPK and HbCRK genes was investigated via Solexa sequencing in a range of experimental conditions (different tissues, phases of leaf development, ethylene treatment, and various abiotic stresses). The results suggest that HbCPK and HbCRK genes are important components in growth, development, and stress responses of H. brasiliensis. Parallel studies on the CDPK and CRK gene families were also extended to five other plant species (Arabidopsis thaliana, Oryza sativa, Populus trichocarpa, Manihot esculenta, and Ricinus communis). The CDPK and CRK genes from different plant species that exhibit similar expression patterns tend to cluster together, suggesting a coevolution of gene structure and expression behavior in higher plants. The results serve as a foundation to further functional studies of these gene families in H. brasiliensis as well as in the whole plant kingdom.

Identification and functional characterization of HbOsmotin from Hevea brasiliensis

Latex in the laticiferous cell network of Hevea brasiliensis tree is composed of cytoplasm that synthesizes natural rubber. Ethylene stimulation of the tree bark enhances latex production partly by prolonging the duration of latex flow during the tapping process. Here, we identified an osmotin-like cDNA sequence (HbOsmotin) from H. brasiliensis that belongs to the pathogenesis-related 5 (PR-5) gene family. The HbOsmotin protein is present in the lutoids of latex in H. brasiliensis, whereas in onion epidermal cells, this protein is predominantly distributed around the cell wall, suggesting that it may be secreted from the cytoplasm. We investigated the effects of exogenous ethylene on HbOsmotin transcription and protein accumulation in rubber latex, and further determined the protein function after osmotic stress in Arabidopsis. In regularly tapped trees, HbOsmotin expression was drastically inhibited in rubber latex after tapping, although the expression was subsequently recovered by ethylene stimulation. However, in virgin plants that had never been tapped, exogenous ethylene application slightly decreased HbOsmotin expression. HbOsmotin overexpression in Arabidopsis showed that HbOsmotin reduced the osmotic stress tolerance of the plant, which likely occurred by raising the water potential. These data indicated that HbOsmotin may contribute to osmotic regulation in laticiferous cells.

Characterization of the interaction between Oidium heveae and Arabidopsis thaliana

Oidium heveae, an obligate biotrophic pathogen of rubber trees (Hevea brasiliensis), causes significant yield losses of rubber worldwide. However, the molecular mechanisms underlying the interplay between O. heveae and rubber trees remain largely unknown. In this study, we isolated an O. heveae strain, named HN1106, from cultivated H. brasiliensis in Hainan, China. We found that O. heveae HN1106 triggers the hypersensitive response in a manner that depends on the effector-triggered immunity proteins EDS1 (Enhanced Disease Susceptibility 1) and PAD4 (Phytoalexin Deficient 4) and on salicylic acid (SA) in the model plant Arabidopsis thaliana. However, SA-independent resistance also appears to limit O. heveae infection of Arabidopsis, because the pathogen does not produce conidiospores on npr1 (nonexpressor of pr1), sid2 (SA induction deficient 2) and NahG plants, which show disruptions in SA signalling. Furthermore, we found that the callose synthase PMR4 (Powdery Mildew Resistant 4) prevents O. heveae HN1106 penetration into leaves in the early stages of infection. To elucidate the potential mechanism of resistance of Arabidopsis to O. heveae HN1106, we inoculated 47 different Arabidopsis accessions with the pathogen, and analysed the plant disease symptoms and O. heveae HN1106 hyphal growth and conidiospore formation on the leaves. We found that the accession Lag2-2 showed significant susceptibility to O. heveae HN1106. Overall, this study provides a basis for future research aimed at combatting powdery mildew caused by O. heveae in rubber trees.

Diversity of root-associated arbuscular mycorrhizal fungal communities in a rubber tree plantation chronosequence in Northeast Thailand

Rubber tree (Hevea brasiliensis) is of major economic importance in Southeast Asia and for small land holders in Thailand in particular. Due to the high value of latex, plantations are expanding into unsuitable areas, such as the northeast province of Thailand where soil fertility is very low and therefore appropriate management practices are of primary importance. Arbuscular mycorrhizal fungi (AMF) contribute to plant growth through a range of mechanisms and could play a key role in a more sustainable management of the rubber plantations. We described the diversity of AMF associated with rubber tree roots in Northeast Thailand in relation to tree age and soil parameters along a chronosequence of rubber tree plantations. Cassava fields were included for comparison. Rubber tree and cassava roots harbored high diversity of AMF (111 Virtual Taxa, VT), including 20 novel VT. AMF VT richness per sample was consistently high (per site mean 16 to 21 VT per sample) along the chronosequence and was not related to soil properties. The composition of AMF communities differed between cassava and rubber tree plantations and was influenced by soil texture and nutrient content (sand, K, P, Ca). AMF community composition gradually shifted with the age of the trees. Our results suggest that the high diversity of AMF in this region is potentially significant for maintaining high functionality of AMF communities.

Identification and reconstitution of the rubber biosynthetic machinery on rubber particles from Hevea brasiliensis

Natural rubber (NR) is stored in latex as rubber particles (RPs), rubber molecules surrounded by a lipid monolayer. Rubber transferase (RTase), the enzyme responsible for NR biosynthesis, is believed to be a member of the cis-prenyltransferase (cPT) family. However, none of the recombinant cPTs have shown RTase activity independently. We show that HRT1, a cPT from Heveabrasiliensis, exhibits distinct RTase activity in vitro only when it is introduced on detergent-washed HeveaRPs (WRPs) by a cell-free translation-coupled system. Using this system, a heterologous cPT from Lactucasativa also exhibited RTase activity, indicating proper introduction of cPT on RP is the key to reconstitute active RTase. RP proteomics and interaction network analyses revealed the formation of the protein complex consisting of HRT1, rubber elongation factor (REF) and HRT1-REF BRIDGING PROTEIN. The RTase activity enhancement observed for the complex assembled on WRPs indicates the HRT1-containing complex functions as the NR biosynthetic machinery.

DOI:http://dx.doi.org/10.7554/eLife.19022.001

Identification of differentially expressed genes and signalling pathways in bark of Hevea brasiliensis seedlings associated with secondary laticifer differentiation using gene expression microarray

The natural rubber of Para rubber tree, Hevea brasiliensis, is the main crop involved in industrial rubber production due to its superior quality. The Hevea bark is commercially exploited to obtain latex, which is produced from the articulated secondary laticifer. The laticifer is well defined in the aspect of morphology; however, only some genes associated with its development have been reported. We successfully induced secondary laticifer in the jasmonic acid (JA)-treated and linolenic acid (LA)-treated Hevea bark but secondary laticifer is not observed in the ethephon (ET)-treated and untreated Hevea bark. In this study, we analysed 27,195 gene models using NimbleGen microarrays based on the Hevea draft genome. 491 filtered differentially expressed (FDE) transcripts that are common to both JA- and LA-treated bark samples but not ET-treated bark samples were identified. In the Eukaryotic Orthologous Group (KOG) analysis, 491 FDE transcripts belong to different functional categories that reflect the diverse processes and pathways involved in laticifer differentiation. In the Kyoto Encyclopedia of Genes and Genomes (KEGG) and KOG analysis, the profile of the FDE transcripts suggest that JA- and LA-treated bark samples have a sufficient molecular basis for secondary laticifer differentiation, especially regarding secondary metabolites metabolism. FDE genes in this category are from the cytochrome (CYP) P450 family, ATP-binding cassette (ABC) transporter family, short-chain dehydrogenase/reductase (SDR) family, or cinnamyl alcohol dehydrogenase (CAD) family. The data includes many genes involved in cell division, cell wall synthesis, and cell differentiation. The most abundant transcript in FDE list was SDR65C, reflecting its importance in laticifer differentiation. Using the Basic Local Alignment Search Tool (BLAST) as part of annotation and functional prediction, several characterised as well as uncharacterized transcription factors and genes were found in the dataset. Hence, the further characterization of these genes is necessary to unveil their role in laticifer differentiation. This study provides a platform for the further characterization and identification of the key genes involved in secondary laticifer differentiation.

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.

Molecular evolution and functional characterisation of haplotypes of an important rubber biosynthesis gene in Hevea brasiliensis

Hydroxy-methylglutaryl coenzyme-A synthase (HMGS) is a rate-limiting enzyme in the cytoplasmic isoprenoid biosynthesis pathway leading to natural rubber production in Hevea brasiliensis (rubber). Analysis of the structural variants of this gene is imperative to understand their functional significance in rubber biosynthesis so that they can be properly utilised for ongoing crop improvement programmes in Hevea. We report here allele richness and diversity of the HMGS gene in selected popular rubber clones. Haplotypes consisting of single nucleotide polymorphisms (SNPs) from the coding and non-coding regions with a high degree of heterozygosity were identified. Segregation and linkage disequilibrium analysis confirmed that recombination is the major contributor to the generation of allelic diversity, rather than point mutations. The evolutionarily conserved nature of some SNPs was identified by comparative DNA sequence analysis of HMGS orthologues from diverse taxa, demonstrating the molecular evolution of rubber biosynthesis genes in general. In silico three-dimensional structural studies highlighting the structural positioning of non-synonymous SNPs from different HMGS haplotypes revealed that the ligand-binding site on the enzyme remains impervious to the reported sequence variations. In contrast, gene expression results indicated the possibility of association between specific haplotypes and HMGS expression in Hevea clones, which may have a downstream impact up to the level of rubber production. Moreover, haplotype diversity of the HMGS gene and its putative association with gene expression can be the basis for further genetic association studies in rubber. Furthermore, the data also show the role of SNPs in the evolution of candidate genes coding for functional traits in plants.

The Hevea brasiliensis XIP aquaporin subfamily: genomic, structural and functional characterizations with relevance to intensive latex harvesting

X-Intrinsic Proteins (XIP) were recently identified in a narrow range of plants as a full clade within the aquaporins. These channels reportedly facilitate the transport of a wide range of hydrophobic solutes. The functional roles of XIP in planta remain poorly identified. In this study, we found three XIP genes (HbXIP1;1, HbXIP2;1 and HbXIP3;1) in the Hevea brasiliensis genome. Comprehensive bioinformatics, biochemical and structural analyses were used to acquire a better understanding of this AQP subfamily. Phylogenetic analysis revealed that HbXIPs clustered into two major groups, each distributed in a specific lineage of the order Malpighiales. Tissue-specific expression profiles showed that only HbXIP2;1 was expressed in all the vegetative tissues tested (leaves, stem, bark, xylem and latex), suggesting that HbXIP2;1 could take part in a wide range of cellular processes. This is particularly relevant to the rubber-producing laticiferous system, where this isoform was found to be up-regulated during tapping and ethylene treatments. Furthermore, the XIP transcriptional pattern is significantly correlated to latex production level. Structural comparison with SoPIP2;1 from Spinacia oleracea species provides new insights into the possible role of structural checkpoints by which HbXIP2;1 ensures glycerol transfer across the membrane. From these results, we discuss the physiological involvement of glycerol and HbXIP2;1 in water homeostasis and carbon stream of challenged laticifers. The characterization of HbXIP2;1 during rubber tree tapping lends new insights into molecular and physiological response processes of laticifer metabolism in the context of latex exploitation.

HbIDI, SlIDI and EcIDI: A comparative study of isopentenyl diphosphate isomerase activity and structure

In this study, we cloned, expressed and purified the isopentenyl diphosphate isomerases (IDIs) from two plants, Hevea brasiliensis and Solanum lycopersicum, and compared them to the already well characterized Escherichia coli IDI. Phylogenetic analysis showed high homology between the three enzymes. Their catalytic activity was investigated in vitro with recombinant purified enzymes and in vivo by complementation colorimetric tests. The three enzymes displayed consistent activities both in vitro and in vivo. In term of structure, studied by ATR-FTIR and molecular modeling, it is clear that both plant enzymes are more related to their human homologue than to E. coli IDI. But it is assumed that EcIDI represent the minimalistic part of the catalytic core, as both plant enzymes present a supplementary sequence forming an extra α-helice surrounding the catalytic site that could facilitate the biocatalysis. New potential biotechnological applications may be envisaged.

Micromorphological characterization and label-free quantitation of small rubber particle protein in natural rubber latex

Commercial natural rubber is traditionally supplied by Hevea brasiliensis, but now there is a big energy problem because of the limited resource and increasing demand. Intensive study of key rubber-related substances is urgently needed for further research of in vitro biosynthesis of natural rubber. Natural rubber is biosynthesized on the surface of rubber particles. A membrane protein called small rubber particle protein (SRPP) is a key protein associated closely with rubber biosynthesis; however, SRPP in different plants has been only qualitatively studied, and there are no quantitative reports so far. In this work, H. brasiliensis was chosen as a model plant. The microscopic distribution of SRPP on the rubber particles during the washing process was investigated by transmission electron microscopy-immunogold labeling. A label-free surface plasmon resonance (SPR) immunosensor was developed to quantify SRPP in H. brasiliensis for the first time. The immunosensor was then used to rapidly detect and analyze SRPP in dandelions and prickly lettuce latex samples. The label-free SPR immunosensor can be a desirable tool for rapid quantitation of the membrane protein SRPP, with excellent assay efficiency, high sensitivity, and high specificity. The method lays the foundation for further study of the functional relationship between SRPP and natural rubber content.

Species Turnover through Time: Colonization and Extinction Dynamics across Metacommunities

Island biogeography and metacommunity theory often use equilibrium assumptions to predict local diversity, yet nonequilibrium dynamics are common in nature. In nonequilibrium communities, local diversity fluctuates through time as the relative importance of colonization and extinction change. Here, we test the prevalence and causes of nonequilibrium dynamics in metacommunities of mites associated with rubber trees distributed over large spatial (>1,000 km) and temporal (>30-60 generations) scales in Brazil. We measured colonization and extinction rates to test species turnover and nonequilibrium dynamics over a growing season. Mite metacommunities exhibited nonequilibrium dynamics for most months of the year, and these dynamics tracked climatic conditions. Monthly shifts in temperature of more than 1°C resulted in nonequilibrium dynamics, as did mean temperatures outside of two critical ranges. Nonequilibrium dynamics were caused by a change in colonization with temperature change and changes in both colonization and extinction with absolute temperature. Species turnover showed different trends; high relative humidity increased both colonization and extinction rates, increasing turnover but not nonequilibrium dynamics. Our study illustrates that testing nonequilibrium dynamics can provide new insights into the drivers of colonization, extinction, and diversity fluctuations in metacommunities.

Molecular cloning and characterization of a novel bi-functional α-amylase/subtilisin inhibitor from Hevea brasiliensis

A novel cDNA encoding a bi-functional α-amylase/subtilisin inhibitor (HbASI) was isolated from rubber (Hevea brasiliensis) leaves cultivar RRIM600. The HbASI had strong homology with the soybean trypsin inhibitor (Kunitz) family of protease inhibitors. Its putative amino acid sequence was similar to that of the α-amylase/subtilisin inhibitor from Ricinus communis (72% identity). Genomic sequencing indicated that the HbASI gene contained no introns. The messenger RNA of HbASI was detected in leaf, hypocotyl and root. The recombinant HbASI expressed extracellularly in Pichia pastoris exhibited inhibitory activity against α-amylase from Aspergillus oryzae, trypsin and subtilisin A. The HbASI gene was induced in the rubber leaves infected with a rubber tree pathogen, Phytophthora palmivora. It was also enhanced by salicylic acid (SA) treatment and mechanical wounding. In addition, the biological activity of the HbASI protein involving in the plant defence responses was also investigated. The HbASI at a concentration of 0.16 mg mL(-1) could inhibit the mycelium growth of P. palmivora. These data suggested that the HbASI protein might play a crucial role in defence against pathogen of rubber trees.

Molecular mechanism of ethylene stimulation of latex yield in rubber tree (Hevea brasiliensis) revealed by de novo sequencing and transcriptome analysis

Background

Rubber tree (Hevea brasiliensis) is an important industrial crop cultivated in tropical areas for natural rubber production. Treatment of the bark of rubber trees with ehephon (an ethylene releaser) has been a routine measure to increase latex yield, but the molecular mechanism behind the stimulation of rubber production by ethylene still remains a puzzle. Deciphering the enigma is of great importance for improvement of rubber tree for high yield.

Results

De novo sequencing and assembly of the bark transciptomes of Hevea brasiliensis induced with ethephon for 8 h (E8) and 24 h (E24) were performed. 51,965,770, 52,303,714 and 53,177,976 high-quality clean reads from E8, E24 and C (control) samples were assembled into 81,335, 80,048 and 80,800 unigenes respectively, with a total of 84,425 unigenes and an average length of 1,101 bp generated. 10,216 and 9,374 differentially expressed genes (DEGs) in E8 and E24 compared with C were respectively detected. The expression of several enzymes in crucial points of regulation in glycolysis were up-regulated and DEGs were not significantly enriched in isopentenyl diphosphate (IPP) biosynthesis pathway. In addition, up-regulated genes of great regulatory importance in carbon fixation (Calvin cycle) were identified.

Conclusions

The rapid acceleration of glycolytic pathway supplying precursors for the biosynthesis of IPP and natural rubber, instead of rubber biosynthesis per se, may be responsible for ethylene stimulation of latex yield in rubber tree. The elevated rate of flux throughout the Calvin cycle may account for some durability of ethylene-induced stimulation. Our finding lays the foundations for molecular diagnostic and genetic engineering for high-yielding improvement of rubber tree.

Electronic supplementary material

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

De novo transcriptomic assembly and profiling of Rigidoporus microporus during saprotrophic growth on rubber wood

BACKGROUND

The basidiomycete Rigidoporus microporus is a fungus that causes the white rot disease of the tropical rubber tree, Hevea brasiliensis, the major source of commercial natural rubber. Besides its lifestyle as a pathogen, the fungus is known to switch to saprotrophic growth on wood with the ability to degrade both lignin and cellulose. There is almost no genomic or transcriptomic information on the saprotrophic abilities of this fungus. In this study, we present the fungal transcriptomic profiles during saprotrophic growth on rubber wood.

RESULTS

A total of 266.6 million RNA-Seq reads were generated from six libraries of the fungus growing either on rubber wood or without wood. De novo assembly produced 34, 518 unigenes with an average length of 2179 bp. Annotation of unigenes using public databases; GenBank, Swiss-Prot, Kyoto Encyclopedia of Genes and Genomes (KEGG), Cluster of Orthologous Groups (COG) and Gene Ontology (GO) produced 25, 880 annotated unigenes. Transcriptomic profiling analysis revealed that the fungus expressed over 300 genes encoding lignocellulolytic enzymes. Among these, 175 genes were up-regulated in rubber wood. These include three members of the glycoside hydrolase family 43, as well as various glycosyl transferases, carbohydrate esterases and polysaccharide lyases. A large number of oxidoreductases which includes nine manganese peroxidases were also significantly up-regulated in rubber wood. Several genes involved in fatty acid metabolism and degradation as well as natural rubber degradation were expressed in the transcriptome. Four genes (acyl-CoA synthetase, enoyl-CoA hydratase, 3-hydroxyacyl-CoA dehydrogenase and acyl-CoA acetyltransferase) potentially involved in rubber latex degradation pathway were also induced. A number of ATP binding cassette (ABC) transporters and hydrophobin genes were significantly expressed in the transcriptome during saprotrophic growth. Some genes related to energy metabolism were also induced.

CONCLUSIONS

The analysed data gives an insight into the activation of lignocellulose breakdown machinery of R. microporus. This study also revealed genes with relevance in antibiotic metabolism (e.g. cephalosporin esterase) as well as those with potential applications in fatty acid degradation. This is the first study on the transcriptomic analysis of R. microporus on rubber wood and should serve as a pioneering resource for future studies of the fungus at the genomic or transcriptomic level.

Identification and Characterization of the Glucose-6-Phosphate Dehydrogenase Gene Family in the Para Rubber Tree, Hevea brasiliensis

As a key enzyme in the pentose phosphate pathway (PPP), glucose-6-phosphate dehydrogenase (G6PDH) provides nicotinamide adenine dinucleotide phosphate (NADPH) and intermediary metabolites for rubber biosynthesis, and plays an important role in plant development and stress responses. In this study, four Hevea brasiliensis (Para rubber tree) G6PDH genes (HbG6PDH1 to 4) were identified and cloned using a genome-wide scanning approach. All four HbG6PDH genes encode functional G6PDH enzymes as shown by heterologous expression in E. coli. Phylogeny analysis and subcellular localization prediction show that HbG6PDH3 is a cytosolic isoform, while the other three genes (HbG6PDH1, 2 and 4) are plastidic isoforms. The subcellular locations of HbG6PDH3 and 4, two latex-abundant isoforms were further verified by transient expression in rice protoplasts. Enzyme activity assay and expression analysis showed HbG6PDH3 and 4 were implicated in PPP during latex regeneration, and to influence rubber production positively in rubber tree. The cytosolic HbG6PDH3 is a predominant isoform in latex, implying a principal role for this isoform in controlling carbon flow and NADPH production in the PPP during latex regeneration. The expression pattern of plastidic HbG6PDH4 correlates well with the degree of tapping panel dryness, a physiological disorder that stops the flow of latex from affected rubber trees. In addition, the four HbG6PDHs responded to temperature and drought stresses in root, bark, and leaves, implicating their roles in maintaining redox balance and defending against oxidative stress.