The association of single nucleotide polymorphisms with endometrial cancer risk and prognosis

Endometrial cancer is one of the most common female diseases in developed nations and is the most commonly diagnosed gynaecological cancer in Australia. The disease is commonly classified by histology: endometrioid or non-endometrioid endometrial cancer. While non-endometrioid endometrial cancers are accepted to be high-grade, aggressive cancers, endometrioid cancers (comprising 80% of all endometrial cancers diagnosed) generally carry a favourable patient prognosis. However, endometrioid endometrial cancer patients endure significant morbidity due to surgery and radiotherapy used for disease treatment, and patients with recurrent disease have a 5-year survival rate of less than 50%. Genetic analysis of women with endometrial cancer could uncover novel markers associated with disease risk and/or prognosis, which could then be used to identify women at high risk and for the use of specialised treatments. Proteases are widely accepted to play an important role in the development and progression of cancer. This PhD project hypothesised that SNPs from two protease gene families, the matrix metalloproteases (MMPs, including their tissue inhibitors, TIMPs) and the tissue kallikrein-related peptidases (KLKs) would be associated with endometrial cancer susceptibility and/or prognosis. In the first part of this study, optimisation of the genotyping techniques was performed. Results from previously published endometrial cancer genetic association studies were attempted to be validated in a large, multicentre replication set (maximum cases n = 2,888, controls n = 4,483, 3 studies). The rs11224561 progesterone receptor SNP (PGR, A/G) was observed to be associated with increased endometrial cancer risk (per A allele OR 1.31, 95% CI 1.12-1.53; p-trend = 0.001), a result which was initially reported among a Chinese sample set. Previously reported associations for the remaining 8 SNPs investigated for this section of the PhD study were not confirmed, thereby reinforcing the importance of validation of genetic association studies. To examine the effect of SNPs from the MMP and KLK families on endometrial cancer risk, we selected the most significantly associated MMP and KLK SNPs from genome-wide association study analysis (GWAS) to be genotyped in the GWAS replication set (cases n = 4,725, controls n = 9,803, 13 studies). The significance of the MMP24 rs932562 SNP was unchanged after incorporation of the stage 2 samples (Stage 1 per allele OR 1.18, p = 0.002; Combined Stage 1 and 2 OR 1.09, p = 0.002). The rs10426 SNP, located 3' to KLK10 was predicted by bioinformatic analysis to effect miRNA binding. This SNP was observed in the GWAS stage 1 result to exhibit a recessive effect on endometrial cancer risk, a result which was not validated in the stage 2 sample set (Stage 1 OR 1.44, p = 0.007; Combined Stage 1 and 2 OR 1.14, p = 0.08). Investigation of the regions imputed surrounding the MMP, TIMP and KLK genes did not reveal any significant targets for further analysis. Analysis of the case data from the endometrial cancer GWAS to identify genetic variation associated with cancer grade did not reveal SNPs from the MMP, TIMP or KLK genes to be statistically significant. However, the representation of SNPs from the MMP, TIMP and KLK families by the GWAS genotyping platform used in this PhD project was examined and observed to be very low, with the genetic variation of four genes (MMP23A, MMP23B, MMP28 and TIMP1) not captured at all by this technique. This suggests that comprehensive candidate gene association studies will be required to assess the role of SNPs from these genes with endometrial cancer risk and prognosis. Meta-analysis of gene expression microarray datasets curated as part of this PhD study identified a number of MMP, TIMP and KLK genes to display differential expression by endometrial cancer status (MMP2, MMP10, MMP11, MMP13, MMP19, MMP25 and KLK1) and histology (MMP2, MMP11, MMP12, MMP26, MMP28, TIMP2, TIMP3, KLK6, KLK7, KLK11 and KLK12). In light of these findings these genes should be prioritised for future targeted genetic association studies. Two SNPs located 43.5 Mb apart on chromosome 15 were observed from the GWAS analysis to be associated with increased endometrial cancer grade, results that were validated in silico in two independent datasets. One of these SNPs, rs8035725 is located in the 5' untranslated region of a MYC promoter binding protein DENND4A (Stage 1 OR 1.15, p = 9.85 x 10P -5 P, combined Stage 1 and in silico validation OR 1.13, p = 5.24 x 10P -6 P). This SNP has previously been reported to alter the expression of PTPLAD1, a gene involved in the synthesis of very long fatty acid chains and in the Rac1 signaling pathway. Meta-analysis of gene expression microarray data found PTPLAD1 to display increased expression in the aggressive non-endometrioid histology compared with endometrioid endometrial cancer, suggesting that the causal SNP underlying the observed genetic association may influence expression of this gene. Neither rs8035725 nor significant SNPs identified by imputation were predicted bioinformatically to affect transcription factor binding sites, indicating that further studies are required to assess their potential effect on other regulatory elements. The other grade- associated SNP, rs6606792, is located upstream of an inferred pseudogene, ELMO2P1 (Stage 1 OR 1.12, p = 5 x 10P -5 P; combined Stage 1 and in silico validation OR 1.09, p = 3.56 x 10P -5 P). Imputation of the ±1 Mb region surrounding this SNP revealed a cluster of significantly associated variants which are predicted to abolish various transcription factor binding sites, and would be expected to decrease gene expression. ELMO2P1 was not included on the microarray platforms collected for this PhD, and so its expression could not be investigated. However, the high sequence homology of ELMO2P1 with ELMO2, a gene important to cell motility, indicates that ELMO2 could be the parent gene for ELMO2P1 and as such, ELMO2P1 could function to regulate the expression of ELMO2. Increased expression of ELMO2 was seen to be associated with increasing endometrial cancer grade, as well as with aggressive endometrial cancer histological subtypes by microarray meta-analysis. Thus, it is hypothesised that SNPs in linkage disequilibrium with rs6606792 decrease the transcription of ELMO2P1, reducing the regulatory effect of ELMO2P1 on ELMO2 expression. Consequently, ELMO2 expression is increased, cell motility is enhanced leading to an aggressive endometrial cancer phenotype. In summary, these findings have identified several areas of research for further study. The results presented in this thesis provide evidence that a SNP in PGR is associated with risk of developing endometrial cancer. This PhD study also reports two independent loci on chromosome 15 to be associated with increased endometrial cancer grade, and furthermore, genes associated with these SNPs to be differentially expressed according in aggressive subtypes and/or by grade. The studies reported in this thesis support the need for comprehensive SNP association studies on prioritised MMP, TIMP and KLK genes in large sample sets. Until these studies are performed, the role of MMP, TIMP and KLK genetic variation remains unclear. Overall, this PhD study has contributed to the understanding of genetic variation involvement in endometrial cancer susceptibility and prognosis. Importantly, the genetic regions highlighted in this study could lead to the identification of novel gene targets to better understand the biology of endometrial cancer and also aid in the development of therapeutics directed at treating this disease.

Candidate gene association studies : a comprehensive guide to useful in silico tools

The candidate gene approach has been a pioneer in the field of genetic epidemiology, identifying risk alleles and their association with clinical traits. With the advent of rapidly changing technology, there has been an explosion of in silico tools available to researchers, giving them fast, efficient resources and reliable strategies important to find casual gene variants for candidate or genome wide association studies (GWAS). In this review, following a description of candidate gene prioritisation, we summarise the approaches to single nucleotide polymorphism (SNP) prioritisation and discuss the tools available to assess functional relevance of the risk variant with consideration to its genomic location. The strategy and the tools discussed are applicable to any study investigating genetic risk factors associated with a particular disease. Some of the tools are also applicable for the functional validation of variants relevant to the era of GWAS and next generation sequencing (NGS).

Expression Profiles of Mitochondrial Genes in the Frontal Cortex and the Caudate Nucleus of Developing Humans and Mice Selectively Bred for High and Low Fear

A growing body of evidence suggests that mitochondrial function may be important in brain development and psychiatric disorders. However, detailed expression profiles of those genes in human brain development and fear-related behavior remain unclear. Using microarray data available from the public domain and the Gene Ontology analysis, we identified the genes and the functional categories associated with chronological age in the prefrontal cortex (PFC) and the caudate nucleus (CN) of psychiatrically normal humans ranging in age from birth to 50 years. Among those, we found that a substantial number of genes in the PFC (115) and the CN (117) are associated with the GO term: mitochondrion (FDR qv <0.05). A greater number of the genes in the PFC (91%) than the genes in the CN (62%) showed a linear increase in expression during postnatal development. Using quantitative PCR, we validated the developmental expression pattern of four genes including monoamine oxidase B (MAOB), NADH dehydrogenase flavoprotein (NDUFV1), mitochondrial uncoupling protein 5 (SLC25A14) and tubulin beta-3 chain (TUBB3). In mice, overall developmental expression pattern of MAOB, SLC25A14 and TUBB3 in the PFC were comparable to the pattern observed in humans (p<0.05). However, mice selectively bred for high fear did not exhibit normal developmental changes of MAOB and TUBB3. These findings suggest that the genes associated with mitochondrial function in the PFC play a significant role in brain development and fear-related behavior.

Expression pattern of the cannabinoid receptor genes in the frontal cortex of mood disorder patients and mice selectively bred for high and low fear

Although the endocannabinoid system (ECS) has been implicated in brain development and various psychiatric disorders, precise mechanisms of the ECS on mood and anxiety disorders remain unclear. Here, we have investigated developmental and disease-related expression pattern of the cannabinoid receptor 1 (CB1) and the cannabinoid receptor 2 (CB2) genes in the dorsolateral prefrontal cortex (PFC) of humans. Using mice selectively bred for high and low fear, we further investigated potential association between fear memory and the cannabinoid receptor expression in the brain. The CB1, not the CB2, mRNA levels in the PFC gradually decrease during postnatal development ranging in age from birth to 50 years (r 2 > 0.6 & adj. p < 0.05). The CB1 levels in the PFC of major depression patients were higher when compared to the age-matched controls (adj. p < 0.05). In mice, the CB1, not the CB2, levels in the PFC were positively correlated with freezing behavior in classical fear conditioning (p < 0.05). These results suggest that the CB1 in the PFC may play a significant role in regulating mood and anxiety symptoms. Our study demonstrates the advantage of utilizing data from postmortem brain tissue and a mouse model of fear to enhance our understanding of the role of the cannabinoid receptors in mood and anxiety disorders

Gemcitabine reactivates epigenetically silenced genes and functions as a DNA methyltransferase inhibitor

Gemcitabine is indicated in combination with cisplatin as first-line therapy for solid tumours including non-small cell lung cancer (NSCLC), bladder cancer and mesothelioma. Gemcitabine is an analogue of pyrimidine cytosine and functions as an anti-metabolite. Structurally, however, gemcitabine has similarities to 5-aza-2-deoxycytidine (decitabine/Dacogen®), a DNA methyltransferase inhibitor (DNMTi). NSCLC, mesothelioma and prostate cancer cell lines were treated with decitabine and gemcitabine. Reactivation of epigenetically silenced genes was examined by RT-PCR/qPCR. DNA methyltransferase activity in nuclear extracts and recombinant proteins was measured using a DNA methyltransferase assay, and alterations in DNA methylation status were examined using methylation-specific PCR (MS-PCR) and pyrosequencing. We observe a reactivation of several epigenetically silenced genes including GSTP1, IGFBP3 and RASSF1A. Gemcitabine functionally inhibited DNA methyltransferase activity in both nuclear extracts and recombinant proteins. Gemcitabine dramatically destabilised DNMT1 protein. However, DNA CpG methylation was for the most part unaffected by gemcitabine. In conclusion, gemcitabine both inhibits and destabilises DNA methyltransferases and reactivates epigenetically silenced genes having activity equivalent to decitabine at concentrations significantly lower than those achieved in the treatment of patients with solid tumours. This property may contribute to the anticancer activity of gemcitabine.

De Novo Transcriptome Sequence Assembly and Analysis of RNA Silencing Genes of Nicotiana benthamiana

Background: Nicotiana benthamiana has been widely used for transient gene expression assays and as a model plant in the study of plant-microbe interactions, lipid engineering and RNA silencing pathways. Assembling the sequence of its transcriptome provides information that, in conjunction with the genome sequence, will facilitate gaining insight into the plant's capacity for high-level transient transgene expression, generation of mobile gene silencing signals, and hyper-susceptibility to viral infection. Methodology/Results: RNA-seq libraries from 9 different tissues were deep sequenced and assembled, de novo, into a representation of the transcriptome. The assembly, of16GB of sequence, yielded 237,340 contigs, clustering into 119,014 transcripts (unigenes). Between 80 and 85% of reads from all tissues could be mapped back to the full transcriptome. Approximately 63% of the unigenes exhibited a match to the Solgenomics tomato predicted proteins database. Approximately 94% of the Solgenomics N. benthamiana unigene set (16,024 sequences) matched our unigene set (119,014 sequences). Using homology searches we identified 31 homologues that are involved in RNAi-associated pathways in Arabidopsis thaliana, and show that they possess the domains characteristic of these proteins. Of these genes, the RNA dependent RNA polymerase gene, Rdr1, is transcribed but has a 72 nt insertion in exon1 that would cause premature termination of translation. Dicer-like 3 (DCL3) appears to lack both the DEAD helicase motif and second dsRNA binding motif, and DCL2 and AGO4b have unexpectedly high levels of transcription. Conclusions: The assembled and annotated representation of the transcriptome and list of RNAi-associated sequences are accessible at www.benthgenome.com alongside a draft genome assembly. These genomic resources will be very useful for further study of the developmental, metabolic and defense pathways of N. benthamiana and in understanding the mechanisms behind the features which have made it such a well-used model plant. © 2013 Nakasugi et al.

Production of transgenic rice with rice ragged stunt virus synthetic resistance genes

Rice ragged stunt virus (RRSV) is an important pathogen of rice affecting its cultivation in South and South East Asia. An approach based on pathogen derived resistance (PDR) was used to produce RRSV resistant rice cultivars. Sequences from the coding region of RRSV genome segments 7 and 10 (non-structural genes), and 5, 8 and 9 (structural genes) were placed in sense or antisense orientation behind the plant expression promoters CaMV35S, RolC, Ubil, Actl and RBTV. Rice cultivars Taipei 309 and Chinsurah Boro II were transformed by biolistic and/or Agrobacterium-mediated delivery of one or more of these PDR gene constructs. A large number of transgenic lines were produced from calli derived from mature or immature embryos, co-bombarded with the marker gene hph encoding hygromycin resistance and RRSV PDR genes or co-cultivated with strains having the binary vector containing these two genes. Both Mendelian and non-Mendelian segregations were observed in transgenic progeny, especially with transgenic lines produced by biolistics. Preliminary tests conducted in China on selected transgenic lines indicate that plants with RRSV segment 5 antisense PDR gene confer RRSV resistance.

Agrobacterium tumefaciens-mediated transformation of an elite Australian barley cultivar with virus resistance and reporter genes

Efficient transformation of barley cv. Schooner was achieved using Agrobacterium delivery, hygromycin or bialaphos selection and embryogenic callus. Using this system, transgenic plants were generated that contained either the green fluorescent protein gene, or transgenes derived from barley yellow dwarf (BYDV) and cereal yellow dwarf (CYDV) viruses. Many of these plants contained 1-3 transgene copies that were inherited in a simple Mendelian manner. Some plants containing BYDV and/or CYDV derived transgenes showed reduced virus symptoms and rates of viral replication when challenged with the appropriate virus. The ability to transform Schooner is a significant advance for the Australian barley industry, as this elite malting variety is, and has for the last 15 years been, the most widely grown barley variety in eastern Australia.

A rapid and simple method of assaying plants transformed with hygromycin or PPT resistance genes

A very simple leaf assay is described that rapidly and reliably identifies transgenic plants expressing the hygromycin resistance gene, hph or the phosphinothricin resistance gene, bar. Leaf tips were cut from plants propagated either in the glasshouse or in tissue culture and the cut surface embedded in solid medium containing the appropriate selective agent. Non-transgenic barley or rice leaf tips had noticeable symptoms of either bleaching or necrosis after three days on the medium and were completely bleached or necrotic after one week. Transgenic leaf tips remained green and healthy over this period. This gave unambiguous discrimination between transgenic and non-transgenic plants. The leaf assay was also effective for dicot plants tested (tobacco and peas).

Comparison of promoters and selectable marker genes for use in Indica rice transformation

The effectiveness of different promoters for use in Indica rice transformation was compared. Plasmids encoding the Escherichia coli uidA (gus) gene under the control of CaMV 35S, Emu, Act1 or Ubi1 promoters were delivered into cell suspension cultures by particle bombardment. Transient gene expression, 48 h after delivery, was greatest from plasmids utilising the constitutive promoters, Act1 and Ubi1. Gene expression in stably transformed tissue was examined by bombarding embryogenic Indica rice calli with a pUbil-gus plasmid and a plasmid containing either the selectable marker gene, hph, which confers hygromycin resistance, or bar, which confers resistance to the herbicide phosphinothricin (BASTA) each under the control of the CaMV 35S, Emu, Act1 or the Ubi1 promoters. The bombarded calli were placed on the appropriate selection media and stained for GUS activity at 1 day, 3 weeks and 5 weeks after shooting. Callus bombarded with the pUbi1-hph or the pEmu-hph constructs gave a dramatic increase in the size of the GUS staining areas with time. No such increase in the size of GUS staining areas was observed in calli co-bombarded with pUbi1-gus and any of the bar containing constructs. Co-bombardment of calli with either the pEmu-hph or pUbi1-hph construct and a virus minor coat protein (cp) gene construct resulted in many fertile transgenic Indica rice plants, containing one to eight copies of both the hph and cp genes. These genes were stably inherited by the T 1 generation.

Bcl-2 genes and growth factors in the pathology of ischaemic acute renal failure

For the past decade, an attempt has been made by many research groups to define the roles of the growing number of Bcl-2 gene family proteins in the apoptotic process. The Bcl-2 family consists of pro-apoptotic (or cell death) and anti-apoptotic (or cell survival) genes and it is the balance in expression between these gene lineages that may determine the death or survival of a cell. The majority of studies have analysed the role/s of the Bcl-2 genes in cancer development. Equally important is their role in normal tissue development, homeostasis and non-cancer disease states. Bcl-2 is crucial for normal development in the kidney, with a deficiency in Bcl-2 producing such malformation that renal failure and death result. As a corollary, its role in renal disease states in the adult has been sought. Ischaemia is one of the most common causes of both acute and chronic renal failure. The section of the kidney that is most susceptible to ischaemic damage is the outer zone of the outer medulla. Within this zone the proximal tubules are most sensitive and often die by necrosis or desquamate. In the distal nephron, apoptosis is the more common form of cell death. Recent results from our laboratory have indicated that ischaemia-induced acute renal failure is associated with up-regulation of two anti-apoptotic Bcl-2 proteins (Bcl-2 and Bcl-XL) in the damaged distal tubule and occasional up-regulation of Bax in the proximal tubule. The distal tubule is a known reservoir for several growth factors important to renal growth and repair, such as insulin-like growth factor-1 (IGF-1) and epidermal growth factor (EGF). One of the likely possibilities for the anti-cell death action of the Bcl-2 genes is that the protected distal cells may be able to produce growth factors that have a further reparative or protective role via an autocrine mechanism in the distal segment and a paracrine mechanism in the proximal cells. Both EGF and IGF-1 are also up-regulated in the surviving distal tubules and are detected in the surviving proximal tubules, where these growth factors are not usually synthesized. As a result, we have been using in vitro methods to test: (i) the relative sensitivities of renal distal and proximal epithelial cell populations to injury caused by mechanisms known to act in ischaemia-reperfusion; (ii) whether a Bcl-2 anti-apoptotic mechanism acts in these cells; and (iii) whether an autocrine and/or paracrine growth factor mechanism is initiated. The following review discusses the background to these studies as well as some of our preliminary results.

Genetic analysis of GRIA2 and GRIA4 genes in migraine

Background Migraine is a brain disorder affecting ∼12% of the Caucasian population. Genes involved in neurological, vascular, and hormonal pathways have all been implicated in predisposing individuals to developing migraine. The migraineur presents with disabling head pain and varying symptoms of nausea, emesis, photophobia, phonophobia, and occasionally visual sensory disturbances. Biochemical and genetic studies have demonstrated dysfunction of neurotransmitters: serotonin, dopamine, and glutamate in migraine susceptibility. Glutamate mediates the transmission of excitatory signals in the mammalian central nervous system that affect normal brain function including cognition, memory and learning. The aim of this study was to investigate polymorphisms in the GRIA2 and GRIA4 genes, which encode subunits of the ionotropic AMPA receptor for association in an Australian Caucasian population. Methods Genotypes for each polymorphism were determined using high resolution melt analysis and the RFLP method. Results Statistical analysis showed no association between migraine and the GRIA2 and GRIA4 polymorphisms investigated. Conclusions Although the results of this study showed no significant association between the tested GRIA gene variants and migraine in our Australian Caucasian population further investigation of other components of the glutamatergic system may help to elucidate if there is a relationship between glutamatergic dysfunction and migraine.

The reduced expression of the HADH2 protein causes X-linked mental retardation, choreoathetosis, and abnormal behavior

Recently, we defined a new syndromic form of X-linked mental retardation in a 4-generation family with a unique clinical phenotype characterized by mild mental retardation, choreoathetosis, and abnormal behavior (MRXS10). Linkage analysis in this family revealed a candidate region of 13.4 Mb between markers DXS1201 and DXS991 on Xp11; therefore, mutation analysis was performed by direct sequencing in most of the 135 annotated genes located in the region. The gene (HADH2) encoding L-3-hydroxyacyl-CoA dehydrogenase II displayed a sequence alteration (c.574 C-->A; p.R192R) in all patients and carrier females that was absent in unaffected male family members and could not be found in 2,500 control X chromosomes, including in those of 500 healthy males. The silent C-->A substitution is located in exon 5 and was shown by western blot to reduce the amount of HADH2 protein by 60%-70% in the patient. Quantitative in vivo and in vitro expression studies revealed a ratio of splicing transcript amounts different from those normally seen in controls. Apparently, the reduced expression of the wild-type fragment, which results in the decreased protein expression, rather than the increased amount of aberrant splicing fragments of the HADH2 gene, is pathogenic. Our data therefore strongly suggest that reduced expression of the HADH2 protein causes MRXS10, a phenotype different from that caused by 2-methyl-3-hydroxybutyryl-CoA dehydrogenase deficiency, which is a neurodegenerative disorder caused by missense mutations in this multifunctional protein.

Genomic organisation, activity and distribution analysis of the microbial putrescine oxidase degradation pathway

The catalytic action of putrescine specific amine oxidases acting in tandem with 4-aminobutyraldehyde dehydrogenase is explored as a degradative pathway in Rhodococcus opacus. By limiting the nitrogen source, increased catalytic activity was induced leading to a coordinated response in the oxidative deamination of putrescine to 4-aminobutyraldehyde and subsequent dehydrogenation to 4-aminobutyrate. Isolating the dehydrogenase by ion exchange chromatography and gel filtration revealed that the enzyme acts principally on linear aliphatic aldehydes possessing an amino moiety. Michaelis-Menten kinetic analysis delivered a Michaelis constant (KM=0.014mM) and maximum rate (Vmax=11.2μmol/min/mg) for the conversion of 4-aminobutyraldehyde to 4-aminobutyrate. The dehydrogenase identified by MALDI-TOF mass spectrometric analysis (E value=0.031, 23% coverage) belongs to a functionally related genomic cluster that includes the amine oxidase, suggesting their association in a directed cell response. Key regulatory, stress and transport encoding genes have been identified, along with candidate dehydrogenases and transaminases for the further conversion of 4-aminobutyrate to succinate. Genomic analysis has revealed highly similar metabolic gene clustering among members of Actinobacteria, providing insight into putrescine degradation notably among Micrococcaceae, Rhodococci and Corynebacterium by a pathway that was previously uncharacterised in bacteria.

The role of enoyl reductase genes in phloridzin biosynthesis in apple

Phloridzin is the predominant polyphenol in apple (Malus× domestica Borkh.) where it accumulates to high concentrations in many tissues including the leaves, bark, roots and fruit. Despite its relative abundance in apple the biosynthesis of phloridzin and other related dihydrochalcones remains only partially understood. The key unidentified enzyme in phloridzin biosynthesis is a putative carbon double bond reductase which is thought to act on p-coumaroyl-CoA to produce the dihydro p-coumaroyl-CoA precursor. A functional screen of six apple enoyl reductase-like (ENRL) genes was carried out using transient infiltration into tobacco and gene silencing by RNA interference (RNAi) in order to determine carbon double bond reductase activity and contribution to foliar phloridzin concentrations. The ENRL-3 gene caused a significant increase in phloridzin concentration when infiltrated into tobacco leaves whilst a second protein ENRL-5, with over 98% amino acid sequence similarity to ENRL-3, showed p-coumaroyl-CoA reductase activity in enzyme assays. Finally, an RNAi study showed that reducing the transcript levels of ENRL-3 in transgenic 'Royal Gala' led to a 66% decrease in the concentration of dihydrochalcones in the leaves in the one available silenced line. Overall these results suggest that ENRL-3, and its close homolog ENRL-5, may contribute to the biosynthesis of phloridzin in apple.