Evolution of developmental gene regulatory networks in echinoids

Developmental gene regulatory networks (dGRNs) are assemblages of regulatory genes that direct embryonic development of animal body plans and their morpho-logical structures. dGRNs exhibit recursively-wired circuitry that is encoded in the genome and executed during development. Alteration to the regulatory architecture of dGRNs causes variation in developmental programs both during the development of an individual organism and during the evolution of an individual lineage. The ex-planatory power of these networks is best exemplified by the global dGRN directing early development of the euechinoid sea urchin Strongylocentrotus purpuratus. This network consists of numerous regulatory genes engaging in hundreds of genomic regulatory transactions that collectively direct the delineation of early embryonic domains and the specification of cell lineages. Research on closely-related euechi-noid sea urchins, e.g. Lytechinus variegatus and Paracentrotus lividus, has revealed marked conservation of dGRN architecture in echinoid development, suggesting little appreciable alteration has occurred since their divergence in evolution at least 90 million years ago (mya).

We sought to test whether this observation extends to all sea urchins (echinoids) and undertook a systematic analysis of over 50 regulatory genes in the cidaroid sea urchin Eucidaris tribuloides, surveing their regulatory activity and function in a sea urchin that diverged from euechinoid sea urchins at least 268 mya. Our results revealed extensive alterations have occurred to all levels of echinoid dGRN archi-tecture since the cidaroid-euechinoid divergence. Alterations to mesodermal sub-circuits were particularly striking, including functional di˙erences in specification of non-skeletogenic mesenchyme (NSM), skeletogenic mesenchyme (SM), and en-domesodermal segregation. Specification of endomesodermal embryonic domains revealed that, while their underlying network circuitry had clearly diverged, regu-latory states established in pregastrular embryos of these two groups are strikingly similar. Analyses of E. tribuloides specification leading to the estab-lishment of dorsal-ventral (aboral-oral) larval polarity indicated that regulation of regulatory genes expressed in mesodermal embryonic domains had incurred significantly more alterations than those expressed in endodermal and ectodermal domains. Taken together, this study highlights the ability of dGRN architecture to buffer extensive alterations in the evolution and early development of echinoids and adds further support to the notion that alterations can occur at all levels of dGRN architecture and all stages of embryonic development.

Targeted gene modification in Fragaria vesca mediated by CRISPR/Cas9 system

Genome editing is becoming an important biotechnological tool for gene function analysis and crop improvement, being the CRISPR-Cas9 (Clustered Regularly Interspaced Short Palindromic Repeat-CRISPR associated protein 9) system the most widely used. The natural CRISPR/Cas9 system has been reduced to two components: a single-guide RNA (sgRNA) for target recognition via RNA-DNA base pairing, which is commonly expressed using a promoter for small-RNAs (U6 promoter), and the Cas9 endonuclease for DNA cleavage (1). To validate the CRISPR/Cas9 system in strawberry plants, we designed two sgRNAs directed against the floral homeotic gene APETALA3 (sgRNA-AP3#1 and sgRNA-AP3#2). This gene was selected because ap3 mutations induce clear developmental phenotypes in which petals and stamens are missing or partially converted to sepals and carpels respectively (2). In this work, we used two different U6 promoters to drive the sgRNA-AP3s expression: AtU6-26 from Arabidopsis (4), and a U6 promoter from Fragaria vesca (FvU6) (this work). We also tested two different coding sequences of Cas9: a human- (hSpCas9) (3) and a plant-codon optimized (pSpCas9) (this work). Transient expression experiments using both CRISPR/Cas9 systems (AtU6-26:sgRNA-AP3#1_35S:hSpCas9_AtU6-26:sgRNA-AP3#2 and FvU6:sgRNA-AP3#1_35S:pSpCas9_FvU6:sgRNA-AP3#2) were performed infiltrating Agrobacterium tumefaciens into F. vesca fruits. PCR amplification and sequencing analyses across the target sites showed a deletion of 188-189 bp corresponding to the region comprised between the two cutting sites of Cas9, confirming that the CRISPR/Cas9 system is functional in F. vesca. Remarkably, the two systems showed different mutagenic efficiency that could be related to differences in expression of the U6 promoters as well as differences in the Cas9 transcripts stability and translation. Stable transformants for both F. vesca (2n) and Fragaria X anannassa (8n) are currently being established to test whether is possible to obtain heritable homozygous mutants derived from CRISPR/Cas9 strategies in strawberry. Thus, our work offers a promising tool for genome editing and gene functional analysis in strawberry. This tool might represent a more efficient alternative to the sometimes inefficient RNAi silencing methods commonly used in this species.

Arrangement of gene pairs, retrotransposon insertions, and regulation of gene expression in plants

Plant genomes are extremely complex. Myriad factors contribute to their evolution and organization, as well as to the expression and regulation of individual genes. Here we present investigations into several such factors and their influence on genome structure and gene expression: the arrangement of pairs of physically adjacent genes, retrotransposons closely associated with genes, and the effect of retrotransposons on gene pair evolution. All sequenced plant genomes contain a significant fraction of retrotransposons, including that of rice. We investigated the effects of retrotransposons within rice genes and within a 1 kb putative promoter region upstream of each gene. We found that approximately one-sixth of all rice genes are closely associated with retrotransposons. Insertions within a gene’s promoter region tend to block gene expression, while retrotransposons within genes promote the existence of alternative splicing forms. We also identified several other trends in retrotransposon insertion and its effects on gene expression. Several studies have previously noted a connection among genes between physical proximity and correlated expression profiles. To determine the degree to which this correlation depends on an exact physical arrangement, we studied the expression and interspecies conservation of convergent and divergent gene pairs in rice, Arabidopsis, and Populus trichocarpa. Correlated expression among gene pairs was quite common in all three species, yet conserved arrangement was rare. However, conservation of gene pair arrangement was significantly more common among pairs with strongly correlated expression levels. In order to uncover additional properties of gene pair conservation and rearrangement, we performed a comparative analysis of convergent, divergent, and tandem gene pairs in rice, sorghum, maize, and Brachypodium. We noted considerable differences between gene pair types and species. We also constructed a putative evolutionary history for each pair, which led to several interesting discoveries. To further elucidate the causes of gene pair conservation and rearrangement, we identified retrotransposon insertions in and near rice gene pairs. Retrotransposon-associated pairs are less likely to be conserved, although there are significant differences in the possible effect of different types and locations of retrotransposon insertions. The three types of gene pair also varied in their susceptibility to retrotransposon-associated evolutionary changes.

Understanding ten-eleven translocation-2 in hematological and nervous systems

I proposed the study of two distinct aspects of Ten-Eleven Translocation 2 (TET2) protein for understanding specific functions in different body systems. ^ In Part I, I characterized the molecular mechanisms of Tet2 in the hematological system. As the second member of Ten-Eleven Translocation protein family, TET2 is frequently mutated in leukemic patients. Previous studies have shown that the TET2 mutations frequently occur in 20% myelodysplastic syndrome/myeloproliferative neoplasm (MDS/MPN), 10% T-cell lymphoma leukemia and 2% B-cell lymphoma leukemia. Genetic mouse models also display distinct phenotypes of various types of hematological malignancies. I performed 5-hydroxymethylcytosine (5hmC) chromatin immunoprecipitation sequencing (ChIP-Seq) and RNA sequencing (RNA-Seq) of hematopoietic stem/progenitor cells to determine whether the deletion of Tet2 can affect the abundance of 5hmC at myeloid, T-cell and B-cell specific gene transcription start sites, which ultimately result in various hematological malignancies. Subsequent Exome sequencing (Exome-Seq) showed that disease-specific genes are mutated in different types of tumors, which suggests that TET2 may protect the genome from being mutated. The direct interaction between TET2 and Mutator S Homolog 6 (MSH6) protein suggests TET2 is involved in DNA mismatch repair. Finally, in vivo mismatch repair studies show that the loss of Tet2 causes a mutator phenotype. Taken together, my data indicate that TET2 binds to MSH6 to protect genome integrity. ^ In Part II, I intended to better understand the role of Tet2 in the nervous system. 5-hydroxymethylcytosine regulates epigenetic modification during neurodevelopment and aging. Thus, Tet2 may play a critical role in regulating adult neurogenesis. To examine the physiological significance of Tet2 in the nervous system, I first showed that the deletion of Tet2 reduces the 5hmC levels in neural stem cells. Mice lacking Tet2 show abnormal hippocampal neurogenesis along with 5hmC alternations at different gene promoters and corresponding gene expression downregulation. Through the luciferase reporter assay, two neural factors Neurogenic differentiation 1 (NeuroD1) and Glial fibrillary acidic protein (Gfap) were down-regulated in Tet2 knockout cells. My results suggest that Tet2 regulates neural stem/progenitor cell proliferation and differentiation in adult brain.^

Gene expression profile suggests that pigs (Sus scrofa) are susceptible to Anaplasma phagocytophilum but control infection

BACKGROUND Anaplasma phagocytophilum infects a wide variety of hosts and causes granulocytic anaplasmosis in humans, horses and dogs and tick-borne fever in ruminants. Infection with A. phagocytophilum results in the modification of host gene expression and immune response. The objective of this research was to characterize gene expression in pigs (Sus scrofa) naturally and experimentally infected with A. phagocytophilum trying to identify mechanisms that help to explain low infection prevalence in this species. RESULTS For gene expression analysis in naturally infected pigs, microarray hybridization was used. The expression of differentially expressed immune response genes was analyzed by real-time RT-PCR in naturally and experimentally infected pigs. Results suggested that A. phagocytophilum infection affected cytoskeleton rearrangement and increased both innate and adaptive immune responses by up regulation of interleukin 1 receptor accessory protein-like 1 (IL1RAPL1), T-cell receptor alpha chain (TCR-alpha), thrombospondin 4 (TSP-4) and Gap junction protein alpha 1 (GJA1) genes. Higher serum levels of IL-1 beta, IL-8 and TNF-alpha in infected pigs when compared to controls supported data obtained at the mRNA level. CONCLUSIONS These results suggested that pigs are susceptible to A. phagocytophilum but control infection, particularly through activation of innate immune responses, phagocytosis and autophagy. This fact may account for the low infection prevalence detected in pigs in some regions and thus their low or no impact as a reservoir host for this pathogen. These results advanced our understanding of the molecular mechanisms at the host-pathogen interface and suggested a role for newly reported genes in the protection of pigs against A. phagocytophilum.

Identification of a novel a-L-arabinofuranosidase gene associated with mealiness in apple.

In order to investigate the genetic bases of the physiological syndrome mealiness that causes abnormal fruit softening and juice loss in apples, an integrative approach was devised, consisting of sensory, instrumental, biochemical, genetic, and genomic methods. High levels of activity of a-L-arabinofuranosidase (a-AFase), a hydrolase acting on the pectic component of the cell walls, were found in individuals exhibiting the mealiness phenotype in a segregating population. The expression levels of the previously uncharacterized apple AF gene MdAF3 are higher in fruits from plants consistently showing mealiness symptons and high a-AFase activity. The transcription of MdAF3 is differentially regulated in distinct genomic contexts and appears to be independent of ethylene. Thus, it is likely to be controlled by endogenous developmental mechanisms associated with fruit ripening. The use of integrative approaches has allowed the identification of a novel contributor to the mealiness phenotype in apple and it has been possible to overcome the problems posed by the unavailability of near-isogenic lines to dissect the genetic bases of a complex physiological trait in woody perennial species.

Ruolo di Glutaredossina-2 nella trasformazione tumorale oncocitaria

Gli oncocitomi sono tumori epiteliali caratterizzati da un accumulo di mitocondri strutturalmente e funzionalmente compromessi, a prognosi generalmente benigna. Le cause genetiche della trasformazione oncocitaria sono tuttora sconosciute; pertanto, lo studio di oncocitomi in contesti familiari sindromici è utile nella ricerca dei determinanti genetici predisponenti il fenotipo. Diversi membri di una famiglia affetta da sindrome dell’iperparatiroidismo con tumore della mandibola (HPT-JT), dovuta ad un'ampia delezione in CDC73, hanno mostrato recidiva di tumori paratiroidei oncocitari. Il sequenziamento dell’esoma ha escluso mutazioni private della famiglia; all'interno della delezione ereditata, tuttavia, sono stati individuati elementi regolatori del gene glutaredossina 2 (GLRX2), codificante un'isoforma mitocondriale deputata alla deglutationilazione proteica reversibile -modificazione modulante l’attività di numerosi target- il cui ruolo nel cancro non è noto. La proteina è risultata assente in tutti i tumori e dimezzata nei tessuti sani dei soggetti. Per indagare se la sua assenza alteri la deglutationilazione proteica predisponendo al fenotipo oncocitario, sono stati generati modelli cellulari TPC1 e HCT116 GLRX2 KO in cui sono stati riscontrati un ridotto tasso proliferativo ed un'alterata glutationilazione proteica, particolarmente in seguito a stress ossidativo. Un esperimento pilota in vivo ha mostrato cellule KO oncocitoidi, con mitocondri morfologicamente alterati, suggerendo che l’alterazione redox innescata dall’assenza di GLRX2 possa indurre una disfunzione metabolica mitocondriale tale da mimare quelle osservate negli oncocitomi. L’analisi proteomica ha individuato diversi target di glutationilazione nei campioni KO identificando proteine del ciclo di Krebs e della catena respiratoria mitocondriale. In particolare, una marcata glutationilazione del complesso della piruvato deidrogenasi (PDHc) è stata correlata ad una ridotta sintesi di ATP dipendente da piruvato. Considerando l'importanza dello stress ossidativo nella fisiopatologia del cancro ed il ruolo del glutatione nella risposta antiossidante, GLRX2 rappresenta un potenziale candidato nella regolazione del metabolismo ossidativo nelle cellule tumorali esposte allo stress e nella modulazione del fenotipo tumorale.

Graphene And Carbon Nanotube Nanocomposite For Gene Transfection.

Graphene and carbon nanotube nanocomposite (GCN) was synthesised and applied in gene transfection of pIRES plasmid conjugated with green fluorescent protein (GFP) in NIH-3T3 and NG97 cell lines. The tips of the multi-walled carbon nanotubes (MWCNTs) were exfoliated by oxygen plasma etching, which is also known to attach oxygen content groups on the MWCNT surfaces, changing their hydrophobicity. The nanocomposite was characterised by high resolution scanning electron microscopy; energy-dispersive X-ray, Fourier transform infrared and Raman spectroscopies, as well as zeta potential and particle size analyses using dynamic light scattering. BET adsorption isotherms showed the GCN to have an effective surface area of 38.5m(2)/g. The GCN and pIRES plasmid conjugated with the GFP gene, forming π-stacking when dispersed in water by magnetic stirring, resulting in a helical wrap. The measured zeta potential confirmed that the plasmid was connected to the nanocomposite. The NIH-3T3 and NG97 cell lines could phagocytize this wrap. The gene transfection was characterised by fluorescent protein produced in the cells and pictured by fluorescent microscopy. Before application, we studied GCN cell viability in NIH-3T3 and NG97 line cells using both MTT and Neutral Red uptake assays. Our results suggest that GCN has moderate stability behaviour as colloid solution and has great potential as a gene carrier agent in non-viral based therapy, with low cytotoxicity and good transfection efficiency.

Carbon Nanoparticles For Gene Transfection In Eukaryotic Cell Lines.

For the first time, oxygen terminated cellulose carbon nanoparticles (CCN) was synthesised and applied in gene transfection of pIRES plasmid. The CCN was prepared from catalytic of polyaniline by chemical vapour deposition techniques. This plasmid contains one gene that encodes the green fluorescent protein (GFP) in eukaryotic cells, making them fluorescent. This new nanomaterial and pIRES plasmid formed π-stacking when dispersed in water by magnetic stirring. The frequencies shift in zeta potential confirmed the plasmid strongly connects to the nanomaterial. In vitro tests found that this conjugation was phagocytised by NG97, NIH-3T3 and A549 cell lines making them fluorescent, which was visualised by fluorescent microscopy. Before the transfection test, we studied CCN in cell viability. Both MTT and Neutral Red uptake tests were carried out using NG97, NIH-3T3 and A549 cell lines. Further, we use metabolomics to verify if small amounts of nanomaterial would be enough to cause some cellular damage in NG97 cells. We showed two mechanisms of action by CCN-DNA complex, producing an exogenous protein by the transfected cell and metabolomic changes that contributed by better understanding of glioblastoma, being the major finding of this work. Our results suggested that this nanomaterial has great potential as a gene carrier agent in non-viral based therapy, with low cytotoxicity, good transfection efficiency, and low cell damage in small amounts of nanomaterials in metabolomic tests.

Desquamation Is A Novel Phenomenon For Collective Prostate Epithelial Cell Deletion After Castration.

The mechanism underlying castration-induced prostate regression, which is a classical physiological concept translated into the therapeutic treatment of advanced prostate cancer, involves epithelial cell apoptosis. In searching for events and mechanisms contributing to prostate regression in response to androgen modulation, we have frequently observed the collective deletion of epithelial cells. This work was undertaken to characterize this phenomenon hereafter named desquamation and to verify its presence after 17β-estradiol (E2) administration. Electron microscopy revealed that the desquamating cells had preserved cell-cell junctions and collapsed nuclear contents. The TUNEL reaction was negative for these cells, which were also negative for cleaved caspases-8, -9, -3 and nuclear apoptosis-inducing factor. Detailed analyses revealed that the condensed chromatin was first affected detaching from the nuclear lamina, which was observable after lamin A immunohistochemistry, suggesting the lack of lamin A degradation. A search in animals treated with supraphysiological E2 employed as an alternative anti-androgen treatment revealed no desquamation. The combined treatment (Cas + E2 group) caused changes particular to each treatment, including desquamation. In conclusion, desquamation appeared as a novel phenomenon contributing to collective prostate epithelial cell deletion, distinct from the classical castration-induced apoptosis and particular to the androgen deprivation resulting from surgical castration, and should be considered as part of the mechanisms promoting organ regression.

Identification Of Target Genes Using Gene Expression Profile Of Granulocytes From Patients With Chronic Myeloid Leukemia Treated With Tyrosine Kinase Inhibitors.

Differential gene expression analysis by suppression subtractive hybridization with correlation to the metabolic pathways involved in chronic myeloid leukemia (CML) may provide a new insight into the pathogenesis of CML. Among the overexpressed genes found in CML at diagnosis are SEPT5, RUNX1, MIER1, KPNA6 and FLT3, while PAN3, TOB1 and ITCH were decreased when compared to healthy volunteers. Some genes were identified and involved in CML for the first time, including TOB1, which showed a low expression in patients with CML during tyrosine kinase inhibitor treatment with no complete cytogenetic response. In agreement, reduced expression of TOB1 was also observed in resistant patients with CML compared to responsive patients. This might be related to the deregulation of apoptosis and the signaling pathway leading to resistance. Most of the identified genes were related to the regulation of nuclear factor κB (NF-κB), AKT, interferon and interleukin-4 (IL-4) in healthy cells. The results of this study combined with literature data show specific gene pathways that might be explored as markers to assess the evolution and prognosis of CML as well as identify new therapeutic targets.

Iis--integrated Interactome System: A Web-based Platform For The Annotation, Analysis And Visualization Of Protein-metabolite-gene-drug Interactions By Integrating A Variety Of Data Sources And Tools.

High-throughput screening of physical, genetic and chemical-genetic interactions brings important perspectives in the Systems Biology field, as the analysis of these interactions provides new insights into protein/gene function, cellular metabolic variations and the validation of therapeutic targets and drug design. However, such analysis depends on a pipeline connecting different tools that can automatically integrate data from diverse sources and result in a more comprehensive dataset that can be properly interpreted. We describe here the Integrated Interactome System (IIS), an integrative platform with a web-based interface for the annotation, analysis and visualization of the interaction profiles of proteins/genes, metabolites and drugs of interest. IIS works in four connected modules: (i) Submission module, which receives raw data derived from Sanger sequencing (e.g. two-hybrid system); (ii) Search module, which enables the user to search for the processed reads to be assembled into contigs/singlets, or for lists of proteins/genes, metabolites and drugs of interest, and add them to the project; (iii) Annotation module, which assigns annotations from several databases for the contigs/singlets or lists of proteins/genes, generating tables with automatic annotation that can be manually curated; and (iv) Interactome module, which maps the contigs/singlets or the uploaded lists to entries in our integrated database, building networks that gather novel identified interactions, protein and metabolite expression/concentration levels, subcellular localization and computed topological metrics, GO biological processes and KEGG pathways enrichment. This module generates a XGMML file that can be imported into Cytoscape or be visualized directly on the web. We have developed IIS by the integration of diverse databases following the need of appropriate tools for a systematic analysis of physical, genetic and chemical-genetic interactions. IIS was validated with yeast two-hybrid, proteomics and metabolomics datasets, but it is also extendable to other datasets. IIS is freely available online at: http://www.lge.ibi.unicamp.br/lnbio/IIS/.

Exposure To Bisphenol-a During Pregnancy Partially Mimics The Effects Of A High-fat Diet Altering Glucose Homeostasis And Gene Expression In Adult Male Mice.

Bisphenol-A (BPA) is one of the most widespread EDCs used as a base compound in the manufacture of polycarbonate plastics. The aim of our research has been to study how the exposure to BPA during pregnancy affects weight, glucose homeostasis, pancreatic β-cell function and gene expression in the major peripheral organs that control energy flux: white adipose tissue (WAT), the liver and skeletal muscle, in male offspring 17 and 28 weeks old. Pregnant mice were treated with a subcutaneous injection of 10 µg/kg/day of BPA or a vehicle from day 9 to 16 of pregnancy. One month old offspring were divided into four different groups: vehicle treated mice that ate a normal chow diet (Control group); BPA treated mice that also ate a normal chow diet (BPA); vehicle treated animals that had a high fat diet (HFD) and BPA treated animals that were fed HFD (HFD-BPA). The BPA group started to gain weight at 18 weeks old and caught up to the HFD group before week 28. The BPA group as well as the HFD and HFD-BPA ones presented fasting hyperglycemia, glucose intolerance and high levels of non-esterified fatty acids (NEFA) in plasma compared with the Control one. Glucose stimulated insulin release was disrupted, particularly in the HFD-BPA group. In WAT, the mRNA expression of the genes involved in fatty acid metabolism, Srebpc1, Pparα and Cpt1β was decreased by BPA to the same extent as with the HFD treatment. BPA treatment upregulated Pparγ and Prkaa1 genes in the liver; yet it diminished the expression of Cd36. Hepatic triglyceride levels were increased in all groups compared to control. In conclusion, male offspring from BPA-treated mothers presented symptoms of diabesity. This term refers to a form of diabetes which typically develops in later life and is associated with obesity.

Evolution Of The Deaths Registry System In Brazil: Associations With Changes In The Mortality Profile, Under-registration Of Death Counts, And Ill-defined Causes Of Death.

This paper examines the spatial pattern of ill-defined causes of death across Brazilian regions, and its relationship with the evolution of completeness of the deaths registry and changes in the mortality age profile. We make use of the Brazilian Health Informatics Department mortality database and population censuses from 1980 to 2010. We applied demographic methods to evaluate the quality of mortality data for 137 small areas and correct for under-registration of death counts when necessary. The second part of the analysis uses linear regression models to investigate the relationship between, on the one hand, changes in death counts coverage and age profile of mortality, and on the other, changes in the reporting of ill-defined causes of death. The completeness of death counts coverage increases from about 80% in 1980-1991 to over 95% in 2000-2010 at the same time the percentage of ill-defined causes of deaths reduced about 53% in the country. The analysis suggests that the government's efforts to improve data quality are proving successful, and they will allow for a better understanding of the dynamics of health and the mortality transition.

Leaf-, Panel- And Latex-expressed Sequenced Tags From The Rubber Tree (hevea Brasiliensis) Under Cold-stressed And Suboptimal Growing Conditions: The Development Of Gene-targeted Functional Markers For Stress Response.

Hevea brasiliensis is a native species of the Amazon Basin of South America and the primary source of natural rubber worldwide. Due to the occurrence of South American Leaf Blight disease in this area, rubber plantations have been extended to suboptimal regions. Rubber tree breeding is time-consuming and expensive, but molecular markers can serve as a tool for early evaluation, thus reducing time and costs. In this work, we constructed six different cDNA libraries with the aim of developing gene-targeted molecular markers for the rubber tree. A total of 8,263 reads were assembled, generating 5,025 unigenes that were analyzed; 912 expressed sequence tags (ESTs) represented new transcripts, and two sequences were highly up-regulated by cold stress. These unigenes were scanned for microsatellite (SSR) regions and single nucleotide polymorphisms (SNPs). In total, 169 novel EST-SSR markers were developed; 138 loci were polymorphic in the rubber tree, and 98 % presented transferability to six other Hevea species. Locus duplication was observed in H. brasiliensis and other species. Additionally, 43 SNP markers in 13 sequences that showed similarity to proteins involved in stress response, latex biosynthesis and developmental processes were characterized. cDNA libraries are a rich source of SSR and SNP markers and enable the identification of new transcripts. The new markers developed here will be a valuable resource for linkage mapping, QTL identification and other studies in the rubber tree and can also be used to evaluate the genetic variability of other Hevea species, which are valuable assets in rubber tree breeding.