Novel Insights into the Molecular Genetic Background of Colorectal Tumors

Many of the genes predisposing to highly penetrant colorectal cancer (CRC) syndromes, including hereditary non-polyposis colorectal cancer (MLH1, MSH2, MSH6, PMS2), familial adenomatous polyposis (APC), Peutz-Jeghers syndrome (LKB1), juvenile polyposis (SMAD4, BMPR1A), MYH-associated polyposis (MYH), and Cowden syndrome (PTEN) have already been discovered. Identification of these genes has allowed a more precise classification of the hereditary CRC syndromes and provided a means for predictive genetic testing and surveillance. Some of the genes are also involved in sporadic cancer forms, and therefore the investigation of the rare CRC syndromes has been a breakthrough for general cancer research. Despite the accumulating knowledge on hereditary cancer syndromes, a significant number of familial CRCs remain molecularly unexplained after genetic testing, reflecting the possibility of other predisposing genes or existence of novel syndromes. Moreover, genetic variants conferring low-penetrance risk are still largely unknown. In this study, we examined the role of some new high- and low-penetrance alleles on CRC predisposition. We identified disease causing MYH mutations in a subset (9%) of patients with APC and AXIN2 mutation negative adenomatous polyposis. Due to differences in the pattern of inheritance and clinical manifestation, screening for mutations in MYH is beneficial in view of genetic counselling and surveillance. A novel functionally deficient MYH founder mutation A459D was identified in the Finnish population, and this finding had immediate clinical implications for genetic counselling of at risk families. Many patients with hamartomatous polyposis remain without molecular diagnosis due to atypical phenotypes. We therefore sought to classify 49 patients with unexplained hamartomatous or hyperplastic/mixed polyposis by extensive molecular analyses of PTEN, LKB1, BMPR1A, SMAD4, ENG, BRAF, MYH, and BHD along with revision of polyp histology. Mutations were identified in 11/49 (22%) of the patients. In 6 cases the molecular diagnosis was re-classified guiding surveillance and decisions for prophylactic surgery. Re-evaluation of polyp histology with subsequent more accurate selection of candidate gene analyses is beneficial and can be recommended for patients with unexplained polyposis. Furthermore, germline mutations in ENG underlying juvenile polyposis were described for the first time, characterizing a possible novel genetically defined form of hereditary CRC. Association analyses on two putative low-penetrance alleles, NOD2 3020insC and MDM2 SNP309 were performed in a population-based series of 1042 Finnish CRC patients and in cancer-free controls. In contrast to previous results, NOD2 3020insC did not associate with CRC or age at disease onset in the Finnish population. These data suggest that NOD2 3020insC alone might not be sufficient for CRC predisposition. MDM2 SNP309 was as common in the CRC cohort as in the healthy controls. Interesting trends, however, were observed, which after correction for multiple testing did not reach statistical significance. SNP309 was more common in female CRC patients and a trend towards an earlier age at disease onset was observed in women with SNP309. Subsequent studies have supported this observation and SNP309 could affect gender- or hormone-related tumorigenesis. Finally, a large-scale unbiased effort was designed to characterize the complete mutatome of CRC with microsatellite instability (MSI). Using an approach combining expression microarray and genome database searches, we were able to identify putative MSI target genes. Further characterization of one of the genes suggested that it might play a role also in microsatellite stable CRC and Peutz-Jeghers syndrome pathogenesis.

Finding genes for economically important traits: Brahman cattle puberty.

Age at puberty is an important component of reproductive performance in beef cattle production systems. Brahman cattle are typically late-pubertal relative to Bos taurus cattle and so it is of economic relevance to select for early age at puberty. To assist selection and elucidate the genes underlying puberty, we performed a genome-wide association study (GWAS) using the BovineSNP50 chip (similar to 54 000 polymorphisms) in Brahman bulls (n = 1105) and heifers (n = 843) and where the heifers were previously analysed in a different study. In a new attempt to generate unbiased estimates of single-nucleotide polymorphism (SNP) effects and proportion of variance explained by each SNP, the available data were halved on the basis of year and month of birth into a calibration and validation set. The traits that defined age at puberty were, in heifers, the age at which the first corpus luteum was detected (AGECL, h(2) = 0.56 +/- 0.11) and in bulls, the age at a scrotal circumference of 26 cm (AGE26, h(2) = 0.78 +/- 0.10). At puberty, heifers were on average older (751 +/- 142 days) than bulls (555 +/- 101 days), but AGECL and AGE26 were genetically correlated (r = 0.20 +/- 0.10). There were 134 SNPs associated with AGECL and 146 SNPs associated with AGE26 (P < 0.0001). From these SNPs, 32 (similar to 22%) were associated (P < 0.0001) with both traits. These top 32 SNPs were all located on Chromosome BTA 14, between 21.95 Mb and 28.4 Mb. These results suggest that the genes located in that region of BTA 14 play a role in pubertal development in Brahman cattle. There are many annotated genes underlying this region of BTA 14 and these are the subject of current research. Further, we identified a region on Chromosome X where markers were associated (P < 1.00E-8) with AGE26, but not with AGECL. Information about specific genes and markers add value to our understanding of puberty and potentially contribute to genomic selection. Therefore, identifying these genes contributing to genetic variation in AGECL and AGE26 can assist with the selection for early onset of puberty.

A high-density association screen of 155 ion transport genes for involvement with common migraine

The clinical overlap between monogenic Familial Hemiplegic Migraine (FHM) and common migraine subtypes, and the fact that all three FHM genes are involved in the transport of ions, suggest that ion transport genes may underlie susceptibility to common forms of migraine. To test this leading hypothesis, we examined common variation in 155 ion transport genes using 5257 single nucleotide polymorphisms (SNPs) in a Finnish sample of 841 unrelated migraine with aura cases and 884 unrelated non-migraine controls. The top signals were then tested for replication in four independent migraine case-control samples from the Netherlands, Germany and Australia, totalling 2835 unrelated migraine cases and 2740 unrelated controls. SNPs within 12 genes (KCNB2, KCNQ3, CLIC5, ATP2C2, CACNA1E, CACNB2, KCNE2, KCNK12, KCNK2, KCNS3, SCN5A and SCN9A) with promising nominal association (0.00041 < P < 0.005) in the Finnish sample were selected for replication. Although no variant remained significant after adjusting for multiple testing nor produced consistent evidence for association across all cohorts, a significant epistatic interaction between KCNB2 SNP rs1431656 (chromosome 8q13.3) and CACNB2 SNP rs7076100 (chromosome 10p12.33) (pointwise P = 0.00002; global P = 0.02) was observed in the Finnish case-control sample. We conclude that common variants of moderate effect size in ion transport genes do not play a major role in susceptibility to common migraine within these European populations, although there is some evidence for epistatic interaction between potassium and calcium channel genes, KCNB2 and CACNB2. Multiple rare variants or trans-regulatory elements of these genes are not ruled out.

Expression Analysis of Host-Induced Genes and Proteins of Potato Pathogen Pectobacterium atrosepticum

Pectobacterium atrosepticum on Gram-negatiivinen bakteeri, joka aiheuttaa perunan tyvi- ja märkämätää. P. atrosepticum bakteerin optimilämpötila on melko alhainen ja se on yleinen lauhkeilla alueilla. Tyvimätä leviää pääasiassa siemenperunan välityksellä ja siksi se on ongelma erityisesti siemenperunan tuotannossa. P. atrosepticum kannan SCRI1043 genomi on julkaistu ja sitä tutkitaan malliorganismina märkä- ja tyvimädän taudinaiheuttamisen ymmärtämiseksi. Tämä opportunistinen taudinaiheuttaja voi elää isäntäkasvissa kuukausia piilevänä, aiheuttamatta näkyviä oireita. Suotuisissa olosuhteissa bakteerit alkavat jakautua ja tuottaa kasvin kudoksia hajottavia entsyymejä. Mädäntyvä kasvimassa tarjoaa ravinteita bakteerien kasvuun ja mahdollistaa isäntäkasvin asuttamisen. Soluseiniä hajottavien entsyymien merkitys taudinaiheuttamisessa on hyvin tunnettu, mutta oireettomasta jaksosta ja taudin alkuvaiheista tiedätään vain vähän. Bakteerin genomi sisältää monia toksiineja, adhesiineja, hemolysiineja ja muita proteiineja, joilla saattaa olla merkitys taudinaiheuttamisessa. Tässä työssä käytettiin proteomiikkaa ja mikrosiruanalysiä P. atrosepticum bakteerin erittyvien proteiinien ja geeniekspression tutkimiseen. Proteiinit, jotka eritetään ulos bakteerista, toimivat todennäköisesti taudinaiheuttamisessa, koska ne ovat suorassa kontaktissa isäntäkasvin kanssa. Analyysit suoritettiin olosuhteissa, jotka muistuttavat kasvin soluvälitilaa: matala pH, vähän ravinteita ja matala lämpötila. Isäntäkasvin läsnäolon vaikutusta proteiinien tuottoon ja geeniekspressioon tutkittiin lisäämällä perunauutetta kasvatusalustaan. Tutkimuksessa tunnistettiin P. atrosepticum bakteerin monia jo tunnettuja ja mahdollisesti taudinaiheuttamiseen liittyviä proteiineja. Perunauute lisäsi hiljattain tunnistetun, proteiinien eritysreittiä (tyyppi VI sekreetio, T6SS) koodaavien geenien ilmentymistä. Lisäksi bakteerin havaittiin erittävän useita T6SS:n liittyviä proteiineja kasvualustaan, johon oli lisätty perunauutetta. T6SS:n merkitys bakteereille on vielä epäselvä ja sen vaikutuksesta taudinaiheuttamiseen on julkaistu ristiriitaisia tuloksia. Märkä- ja tyvimädän ymmärtäminen molekulaarisella tasolla luo pohjan tautien kontrollointiin tähtäävään soveltavaan tutkimukseen. Tämä tutkimus lisää tietoa kasvi-patogeeni- interaktiosta ja sitä voidaan tulevaisuudessa käyttää hyväksi esimerkiksi diagnostiikassa, resistenttien perunalajikkeiden jalostuksessa tai viljely- ja varastointiolosuhteiden parantamisessa.

Evolution and detection of cyanobacterial hepatotoxin synthetase genes

Mass occurrences (blooms) of cyanobacteria are common in aquatic environments worldwide. These blooms are often toxic, due to the presence of hepatotoxins or neurotoxins. The most common cyanobacterial toxins are hepatotoxins: microcystins and nodularins. In freshwaters, the main producers of microcystins are Microcystis, Anabaena, and Planktothrix. Nodularins are produced by strains of Nodularia spumigena in brackish waters. Toxic and nontoxic strains of cyanobacteria co-occur and cannot be differentiated by conventional microscopy. Molecular biological methods based on microcystin and nodularin synthetase genes enable detection of potentially hepatotoxic cyanobacteria. In the present study, molecular detection methods for hepatotoxin-producing cyanobacteria were developed, based on microcystin synthetase gene E (mcyE) and the orthologous nodularin synthetase gene F (ndaF) sequences. General primers were designed to amplify the mcyE/ndaF gene region from microcystin-producing Anabaena, Microcystis, Planktothrix, and Nostoc, and nodularin-producing Nodularia strains. The sequences were used for phylogenetic analyses to study how cyanobacterial mcy genes have evolved. The results showed that mcy genes and microcystin are very old and were already present in the ancestor of many modern cyanobacterial genera. The results also suggested that the sporadic distribution of biosynthetic genes in modern cyanobacteria is caused by repeated gene losses in the more derived lineages of cyanobacteria and not by horizontal gene transfer. Phylogenetic analysis also proposed that nda genes evolved from mcy genes. The frequency and composition of the microcystin producers in 70 lakes in Finland were studied by conventional polymerase chain reaction (PCR). Potential microcystin producers were detected in 84% of the lakes, using general mcyE primers, and in 91% of the lakes with the three genus-specific mcyE primers. Potential microcystin-producing Microcystis were detected in 70%, Planktothrix in 63%, and Anabaena in 37% of the lakes. The presence and co-occurrence of potential microcystin producers were more frequent in eutrophic lakes, where the total phosphorus concentration was high. The PCR results could also be associated with various environmental factors by correlation and regression analyses. In these analyses, the total nitrogen concentration and pH were both associated with the presence of multiple microcystin-producing genera and partly explained the probability of occurrence of mcyE genes. In general, the results showed that higher nutrient concentrations increased the occurrence of potential microcystin producers and the risk for toxic bloom formation. Genus-specific probe pairs for microcystin-producing Anabaena, Microcystis, Planktothrix, and Nostoc, and nodularin-producing Nodularia were designed to be used in a DNA-chip assay. The DNA-chip can be used to simultaneously detect all these potential microcystin/nodularin producers in environmental water samples. The probe pairs detected the mcyE/ndaF genes specifically and sensitively when tested with cyanobacterial strains. In addition, potential microcystin/nodularin producers were identified in lake and Baltic Sea samples by the DNA-chip almost as sensitively as by quantitative real-time PCR (qPCR), which was used to validate the DNA-chip results. Further improvement of the DNA-chip assay was achieved by optimization of the PCR, the first step in the assay. Analysis of the mcy and nda gene clusters from various hepatotoxin-producing cyanobacteria was rewarding; it revealed that the genes were ancient. In addition, new methods detecting all the main producers of hepatotoxins could be developed. Interestingly, potential microcystin-producing cyanobacterial strains of Microcystis, Planktothrix, and Anabaena, co-occurred especially in eutrophic and hypertrophic lakes. Protecting waters from eutrophication and restoration of lakes may thus decrease the prevalence of toxic cyanobacteria and the frequency of toxic blooms.

Association genetics in Corymbia citriodora subsp variegata identifies single nucleotide polymorphisms affecting wood growth and cellulosic pulp yield

Wood is an important biological resource which contributes to nutrient and hydrology cycles through ecosystems, and provides structural support at the plant level. Thousands of genes are involved in wood development, yet their effects on phenotype are not well understood. We have exploited the low genomic linkage disequilibrium (LD) and abundant phenotypic variation of forest trees to explore allelic diversity underlying wood traits in an association study. Candidate gene allelic diversity was modelled against quantitative variation to identify SNPs influencing wood properties, growth and disease resistance across three populations of Corymbia citriodora subsp. variegata, a forest tree of eastern Australia. Nine single nucleotide polymorphism (SNP) associations from six genes were identified in a discovery population (833 individuals). Associations were subsequently tested in two smaller populations (130160 individuals), validating our findings in three cases for actin 7 (ACT7) and COP1 interacting protein 7 (CIP7). The results imply a functional role for these genes in mediating wood chemical composition and growth, respectively. A flip in the effect of ACT7 on pulp yield between populations suggests gene by environment interactions are at play. Existing evidence of gene function lends strength to the observed associations, and in the case of CIP7 supports a role in cortical photosynthesis.

Regulation of oxidative-stress responsive genes by arecoline in human keratinocytes

Background and Objective: Arecoline, an arecanut alkaloid present in the saliva of betel quid chewers, has been implicated in the pathogenesis of a variety of inflammatory oral diseases, including oral submucous fibrosis and periodontitis. To understand the molecular b asis of arecoline action in epithelial changes associated with these diseases, we investigated the effects of arecoline on human keratinocytes with respect to cell growth regulation and the expression of stress-responsive genes.Material and Methods:Human keratinocyte cells (of the HaCaT cell line) were treated with arecoline, following which cell viability was assessed using the Trypan Blue dye-exclusion assay, cell growth and proliferation were analyzed using the MTT (3-[4,5-dimethylthiazol-2-yl]-2,5-diphenyltetrazolium bromide) and 5-bromo-2-deoxyuridine incorporation assays, cell cycle arrest and generation of reactive oxygen species were examined using flow cytometry, and gene expression changes were investigated using the reverse transcription-polymerase chain reaction technique. The role of oxidative stress, muscarinic acetylcholine receptor and mitogen-activated protein kinase (MAPK) pathways were studied using specific inhibitors. Western blot analysis was performed to study p38 MAPK activation.Results:Arecoline induced the generation of reactive oxygen species and cell cycle arrest at the G1/G0 phase in HaCaT cells without affecting the expression of p21/Cip1. Arecoline-induced epithelial cell death at higher concentrations was caused by oxidative trauma without eliciting apoptosis. Sublethal concentrations of arecoline upregulated the expression of the following stress-responsive genes: heme oxygenase-1; ferritin light chain; glucose-6-phosphate dehydrogenase; glutamate-cysteine ligase catalytic subunit; and glutathione reductase.Additionally, there was a dose-dependent induction of interleukin-1alfa mRNA by arecoline via oxidative stress and p38 MAPK activation. Conclusion:our data highlight the role of oxidative stress in arecoline-mediated cell death, gene regulation and inflammatory processes in human keratinocytes.

Sorghum dwarfing genes can affect radiation capture and radiation use efficiency

Presence of the dw3 sorghum dwarfing gene had negative effects on grain yield in some genetic backgrounds and environments. In a previous study we showed that this was due to a significant reduction in shoot biomass (mainly via reduced stem mass), which in turn negatively affected grain size. The current study examines whether shoot biomass was reduced via effects of dw3 on traits associated with resource capture, such as leaf area index (LAI), light interception (LI), and canopy extinction coefficient (k) or with resource use efficiency, such as radiation use efficiency (RUE). Three pairs of near-isogenic sorghum lines differing only in the presence or absence of the dwarfing allele dw3 (3-dwarfs vs 2-dwarfs) were grown in large field plots. Biomass accumulation and LI were measured for individual canopy layers to examine canopy characteristics of tall and short types. Similar to the previously reported effects on grain yield, the effects of dw3 on RUE, LI and k varied among genetic backgrounds and environments. Interactions between dw3 and genetic background, but also interactions with environment are likely to have modulated the extent to which RUE, LI, or k contributed to biomass differences between tall and short sorghum. © 2013 .

Rhipicephalus (Boophilus) microplus tick in vitro feeding methods for functional (dsRNA) and vaccine candidate (antibody) screening

Rhipicephalus (Boophilus) microplus (Acari: Ixodidae) ticks cause economic losses for cattle industries throughout tropical and subtropical regions of the world estimated at $US2.5 billion annually. Lack of access to efficacious long-lasting vaccination regimes and increases in tick acaricide resistance have led to the investigation of targets for the development of novel tick vaccines and treatments. In vitro tick feeding has been used for many tick species to study the effect of new acaricides on the transmission of tick-borne pathogens. Few studies have reported the use of in vitro feeding for functional genomic studies using RNA interference and/or the effect of specific anti-tick antibodies. In particular, in vitro feeding reports for the cattle tick are limited due to its relatively short hypostome. Previously published methods were further modified to broaden optimal tick sizes/weights, feeding sources including bovine and ovine serum, optimisation of commercially available blood anti-coagulant tubes, and IgG concentrations for effective antibody delivery. Ticks are fed overnight and monitored for ∼5–6 weeks to determine egg output and success of larval emergence using a humidified incubator. Lithium-heparin blood tubes provided the most reliable anti-coagulant for bovine blood feeding compared with commercial citrated (CPDA) and EDTA tubes. Although >30 mg semi-engorged ticks fed more reliably, ticks as small as 15 mg also fed to repletion to lay viable eggs. Ticks which gained less than ∼10 mg during in vitro feeding typically did not lay eggs. One mg/ml IgG from Bm86-vaccinated cattle produced a potent anti-tick effect in vitro (83% efficacy) similar to that observed in vivo. Alternatively, feeding of dsRNA targeting Bm86 did not demonstrate anti-tick effects (11% efficacy) compared with the potent effects of ubiquitin dsRNA. This study optimises R. microplus tick in vitro feeding methods which support the development of cattle tick vaccines and treatments.

Do Phosphine Resistance Genes Influence Movement and Dispersal Under Starvation?

Phosphine resistance alleles might be expected to negatively affect energy demanding activities such as walking and flying, because of the inverse relationship between phosphine resistance and respiration. We used an activity monitoring system to quantify walking of Rhyzopertha dominica (F.) and a flight chamber to estimate their propensity for flight initiation. No significant difference in the duration of walking was observed between the strongly resistant, weakly resistant, and susceptible strains of R. dominica we tested, and females walked significantly more than males regardless of genotype. The walking activity monitor revealed no pattern of movement across the day and no particular time of peak activity despite reports of peak activity of R. dominica and Tribolium castaneum (Herbst) under field conditions during dawn and dusk. Flight initiation was significantly higher for all strains at 28 degrees C and 55% relative humidity than at 25, 30, 32, and 35 degrees C in the first 24 h of placing beetles in the flight chamber. Food deprivation and genotype had no significant effect on flight initiation. Our results suggest that known resistance alleles in R. dominica do not affect insect mobility and should therefore not inhibit the dispersal of resistant insects in the field.

The plasticity of NBS resistance genes in sorghum is driven by multiple evolutionary processes

Background Increased disease resistance is a key target of cereal breeding programs, with disease outbreaks continuing to threaten global food production, particularly in Africa. Of the disease resistance gene families, the nucleotide-binding site plus leucine-rich repeat (NBS-LRR) family is the most prevalent and ancient and is also one of the largest gene families known in plants. The sequence diversity in NBS-encoding genes was explored in sorghum, a critical food staple in Africa, with comparisons to rice and maize and with comparisons to fungal pathogen resistance QTL. Results In sorghum, NBS-encoding genes had significantly higher diversity in comparison to non NBS-encoding genes and were significantly enriched in regions of the genome under purifying and balancing selection, both through domestication and improvement. Ancestral genes, pre-dating species divergence, were more abundant in regions with signatures of selection than in regions not under selection. Sorghum NBS-encoding genes were also significantly enriched in the regions of the genome containing fungal pathogen disease resistance QTL; with the diversity of the NBS-encoding genes influenced by the type of co-locating biotic stress resistance QTL. Conclusions NBS-encoding genes are under strong selection pressure in sorghum, through the contrasting evolutionary processes of purifying and balancing selection. Such contrasting evolutionary processes have impacted ancestral genes more than species-specific genes. Fungal disease resistance hot-spots in the genome, with resistance against multiple pathogens, provides further insight into the mechanisms that cereals use in the “arms race” with rapidly evolving pathogens in addition to providing plant breeders with selection targets for fast-tracking the development of high performing varieties with more durable pathogen resistance.

Induced Variations in Brassinosteroid Genes Define Barley Height and Sturdiness, and Expand the Green Revolution Genetic Toolkit

Reduced plant height and culm robustness are quantitative characteristics important for assuring cereal crop yield and quality under adverse weather conditions. A very limited number of short-culm mutant alleles were introduced into commercial crop cultivars during the Green Revolution. We identified phenotypic traits, including sturdy culm, specific for deficiencies in brassinosteroid biosynthesis and signaling in semidwarf mutants of barley (Hordeum vulgare). This set of characteristic traits was explored to perform a phenotypic screen of near-isogenic short-culm mutant lines from the brachytic, breviaristatum, dense spike, erectoides, semibrachytic, semidwarf, and slender dwarf mutant groups. In silico mapping of brassinosteroid-related genes in the barley genome in combination with sequencing of barley mutant lines assigned more than 20 historic mutants to three brassinosteroid-biosynthesis genes (BRASSINOSTEROID-6-OXIDASE, CONSTITUTIVE PHOTOMORPHOGENIC DWARF, and DIMINUTO) and one brassinosteroid-signaling gene (BRASSINOSTEROID-INSENSITIVE1 [HvBRI1]). Analyses of F2 and M2 populations, allelic crosses, and modeling of nonsynonymous amino acid exchanges in protein crystal structures gave a further understanding of the control of barley plant architecture and sturdiness by brassinosteroid-related genes. Alternatives to the widely used but highly temperature-sensitive uzu1.a allele of HvBRI1 represent potential genetic building blocks for breeding strategies with sturdy and climate-tolerant barley cultivars.

Sequence diversity and differential expression of major phenylpropanoid-flavonoid biosynthetic genes among three mango varieties

Background: Mango fruits contain a broad spectrum of phenolic compounds which impart potential health benefits; their biosynthesis is catalysed by enzymes in the phenylpropanoid-flavonoid (PF) pathway. The aim of this study was to reveal the variability in genes involved in the PF pathway in three different mango varieties Mangifera indica L., a member of the family Anacardiaceae: Kensington Pride (KP), Irwin (IW) and Nam Doc Mai (NDM) and to determine associations with gene expression and mango flavonoid profiles. Results: A close evolutionary relationship between mango genes and those from the woody species poplar of the Salicaceae family (Populus trichocarpa) and grape of the Vitaceae family (Vitis vinifera), was revealed through phylogenetic analysis of PF pathway genes. We discovered 145 SNPs in total within coding sequences with an average frequency of one SNP every 316bp. Variety IW had the highest SNP frequency (one SNP every 258bp) while KP and NDM had similar frequencies (one SNP every 369bp and 360bp, respectively). The position in the PF pathway appeared to influence the extent of genetic diversity of the encoded enzymes. The entry point enzymes phenylalanine lyase (PAL), cinnamate 4-mono-oxygenase (C4H) and chalcone synthase (CHS) had low levels of SNP diversity in their coding sequences, whereas anthocyanidin reductase (ANR) showed the highest SNP frequency followed by flavonoid 3'-hydroxylase (F3'H). Quantitative PCR revealed characteristic patterns of gene expression that differed between mango peel and flesh, and between varieties. Conclusions: The combination of mango expressed sequence tags and availability of well-established reference PF biosynthetic genes from other plant species allowed the identification of coding sequences of genes that may lead to the formation of important flavonoid compounds in mango fruits and facilitated characterisation of single nucleotide polymorphisms between varieties. We discovered an association between the extent of sequence variation and position in the pathway for up-stream genes. The high expression of PAL, C4H and CHS genes in mango peel compared to flesh is associated with high amounts of total phenolic contents in peels, which suggest that these genes have an influence on total flavonoid levels in mango fruit peel and flesh. In addition, the particularly high expression levels of ANR in KP and NDM peels compared to IW peel and the significant accumulation of its product epicatechin gallate (ECG) in those extracts reflects the rate-limiting role of ANR on ECG biosynthesis in mango. © 2015 Hoang et al.