Current status of pharmacogenomics testing for anti-tumor drug therapies : approaches to non-melanoma skin cancer

Skin cancer is one of the most commonly occurring cancer types, with substantial social, physical, and financial burdens on both individuals and societies. Although the role of UV light in initiating skin cancer development has been well characterized, genetic studies continue to show that predisposing factors can influence an individual's susceptibility to skin cancer and response to treatment. In the future, it is hoped that genetic profiles, comprising a number of genetic markers collectively involved in skin cancer susceptibility and response to treatment or prognosis, will aid in more accurately informing practitioners' choices of treatment. Individualized treatment based on these profiles has the potential to increase the efficacy of treatments, saving both time and money for the patient by avoiding the need for extensive or repeated treatment. Increased treatment responses may in turn prevent recurrence of skin cancers, reducing the burden of this disease on society. Currently existing pharmacogenomic tests, such as those that assess variation in the metabolism of the anticancer drug fluorouracil, have the potential to reduce the toxic effects of anti-tumor drugs used in the treatment of non-melanoma skin cancer (NMSC) by determining individualized appropriate dosage. If the savings generated by reducing adverse events negate the costs of developing these tests, pharmacogenomic testing may increasingly inform personalized NMSC treatment.

Locked nucleic acid (LNA) single nucleotide polymorphism (SNP) genotype analysis and validation using real-time PCR

With an increased emphasis on genotyping of single nucleotide polymorphisms (SNPs) in disease association studies, the genotyping platform of choice is constantly evolving. In addition, the development of more specific SNP assays and appropriate genotype validation applications is becoming increasingly critical to elucidate ambiguous genotypes. In this study, we have used SNP specific Locked Nucleic Acid (LNA) hybridization probes on a real-time PCR platform to genotype an association cohort and propose three criteria to address ambiguous genotypes. Based on the kinetic properties of PCR amplification, the three criteria address PCR amplification efficiency, the net fluorescent difference between maximal and minimal fluorescent signals and the beginning of the exponential growth phase of the reaction. Initially observed SNP allelic discrimination curves were confirmed by DNA sequencing (n = 50) and application of our three genotype criteria corroborated both sequencing and observed real-time PCR results. In addition, the tested Caucasian association cohort was in Hardy-Weinberg equilibrium and observed allele frequencies were very similar to two independently tested Caucasian association cohorts for the same tested SNP. We present here a novel approach to effectively determine ambiguous genotypes generated from a real-time PCR platform. Application of our three novel criteria provides an easy to use semi-automated genotype confirmation protocol.

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.

Saliva-derived DNA performs well in large-scale, high-density single-nucleotide polymorphism microarray studies

As of June 2009, 361 genome-wide association studies (GWAS) had been referenced by the HuGE database. GWAS require DNA from many thousands of individuals, relying on suitable DNA collections. We recently performed a multiple sclerosis (MS) GWAS where a substantial component of the cases (24%) had DNA derived from saliva. Genotyping was done on the Illumina genotyping platform using the Infinium Hap370CNV DUO microarray. Additionally, we genotyped 10 individuals in duplicate using both saliva- and blood-derived DNA. The performance of blood- versus saliva-derived DNA was compared using genotyping call rate, which reflects both the quantity and quality of genotyping per sample and the “GCScore,” an Illumina genotyping quality score, which is a measure of DNA quality. We also compared genotype calls and GCScores for the 10 sample pairs. Call rates were assessed for each sample individually. For the GWAS samples, we compared data according to source of DNA and center of origin. We observed high concordance in genotyping quality and quantity between the paired samples and minimal loss of quality and quantity of DNA in the saliva samples in the large GWAS sample, with the blood samples showing greater variation between centers of origin. This large data set highlights the usefulness of saliva DNA for genotyping, especially in high-density single-nucleotide polymorphism microarray studies such as GWAS.

Lipoxygenase metabolism : roles in tumor progression and survival

The metabolism of arachidonic acid through lipoxygenase pathways leads to the generation of various biologically active eicosanoids. The expression of these enzymes vary throughout the progression of various cancers, and thereby they have been shown to regulate aspects of tumor development. Substantial evidence supports a functional role for lipoxygenase-catalyzed arachidonic and linoleic acid metabolism in cancer development. Pharmacologic and natural inhibitors of lipoxygenases have been shown to suppress carcinogenesis and tumor growth in a number of experimental models. Signaling of hydro[peroxy]fatty acids following arachidonic or linoleic acid metabolism potentially effect diverse biological phenomenon regulating processes such as cell growth, cell survival, angiogenesis, cell invasion, metastatic potential and immunomodulation. However, the effects of distinct LOX isoforms differ considerably with respect to their effects on both the individual mechanisms described and the tumor being examined. 5-LOX and platelet type 12-LOX are generally considered pro-carcinogenic, with the role of 15-LOX-1 remaining controversial, while 15-LOX-2 suppresses carcinogenesis. In this review, we focus on the molecular mechanisms regulated by LOX metabolism in some of the major cancers. We discuss the effects of LOXs on tumor cell proliferation, their roles in cell cycle control and cell death induction, effects on angiogenesis, migration and the immune response, as well as the signal transduction pathways involved in these processes. Understanding the molecular mechanisms underlying the anti-tumor effect of specific, or general, LOX inhibitors may lead to the design of biologically and pharmacologically targeted therapeutic strategies inhibiting LOX isoforms and/or their biologically active metabolites, that may ultimately prove useful in the treatment of cancer, either alone or in combination with conventional therapies. © 2007 Springer Science+Business Media, LLC.

Facile mutant identification via a single parental backcross method and application of whole genome sequencing based mapping pipelines

Forward genetic screens have identified numerous genes involved in development and metabolism, and remain a cornerstone of biological research. However, to locate a causal mutation, the practice of crossing to a polymorphic background to generate a mapping population can be problematic if the mutant phenotype is difficult to recognize in the hybrid F2 progeny, or dependent on parental specific traits. Here in a screen for leaf hyponasty mutants, we have performed a single backcross of an Ethane Methyl Sulphonate (EMS) generated hyponastic mutant to its parent. Whole genome deep sequencing of a bulked homozygous F2 population and analysis via the Next Generation EMS mutation mapping pipeline (NGM) unambiguously determined the causal mutation to be a single nucleotide polymorphisim (SNP) residing in HASTY, a previously characterized gene involved in microRNA biogenesis. We have evaluated the feasibility of this backcross approach using three additional SNP mapping pipelines; SHOREmap, the GATK pipeline, and the samtools pipeline. Although there was variance in the identification of EMS SNPs, all returned the same outcome in clearly identifying the causal mutation in HASTY. The simplicity of performing a single parental backcross and genome sequencing a small pool of segregating mutants has great promise for identifying mutations that may be difficult to map using conventional approaches.

Description of plant tRNA-derived RNA fragments (tRFs) associated with argonaute and identification of their putative targets

tRNA-derived RNA fragments (tRFs) are 19mer small RNAs that associate with Argonaute (AGO) proteins in humans. However, in plants, it is unknown if tRFs bind with AGO proteins. Here, using public deep sequencing libraries of immunoprecipitated Argonaute proteins (AGO-IP) and bioinformatics approaches, we identified the Arabidopsis thaliana AGO-IP tRFs. Moreover, using three degradome deep sequencing libraries, we identified four putative tRF targets. The expression pattern of tRFs, based on deep sequencing data, was also analyzed under abiotic and biotic stresses. The results obtained here represent a useful starting point for future studies on tRFs in plants. © 2013 Loss-Morais et al.; licensee BioMed Central Ltd.

Hairpin RNAs derived from RNA polymerase II and polymerase III promoter-directed transgenes are processed differently in plants

RNA polymerase III (Pol III) as well as Pol II (35S) promoters are able to drive hairpin RNA (hpRNA) expression and induce target gene silencing in plants. siRNAs of 21 nt are the predominant species in a 35S Pol II line, whereas 24- and/or 22-nucleotide (nt) siRNAs are produced by a Pol III line. The 35S line accumulated the loop of the hpRNA, in contrast to full-length hpRNA in the Pol III line. These suggest that Pol II and Pol III-transcribed hpRNAs are processed by different pathways. One Pol III transgene produced only 24-nt siRNAs but silenced the target gene efficiently, indicating that the 24-nt siRNAs can direct mRNA degradation; specific cleavage was confirmed by 59 rapid amplification of cDNA ends (RACE). Both Pol II- and Pol III-directed hpRNA transgenes induced cytosine methylation in the target DNA. The extent of methylation is not correlated with the level of 21-nt siRNAs, suggesting that they are not effective inducers of DNA methylation. The promoter of a U6 transgene was significantly methylated, whereas the promoter of the endogenous U6 gene was almost free of cytosine methylation, suggesting that endogenous sequences are more resistant to de novo DNA methylation than are transgene constructs. Published by Cold Spring Harbor Laboratory Press. Copyright © 2008 RNA Society.

Regulation of dormancy in barley by blue light and after-ripening : effects on abscisic acid and gibberellin metabolism

White light strongly promotes dormancy in freshly harvested cereal grains, whereas dark and after-ripening have the opposite effect. We have analyzed the interaction of light and after-ripening on abscisic acid (ABA) and gibberellin (GA) metabolism genes and dormancy in barley (Hordeum vulgare ‘Betzes’). Analysis of gene expression in imbibed barley grains shows that different ABA metabolism genes are targeted by white light and after-ripening. Of the genes examined, white light promotes the expression of an ABA biosynthetic gene, HvNCED1, in embryos. Consistent with this result, enzyme-linked immunosorbent assays show that dormant grains imbibed under white light have higher embryo ABA content than grains imbibed in the dark. After-ripening has no effect on expression of ABA biosynthesis genes, but promotes expression of an ABA catabolism gene (HvABA8′OH1), a GA biosynthetic gene (HvGA3ox2), and a GA catabolic gene (HvGA2ox3) following imbibition. Blue light mimics the effects of white light on germination, ABA levels, and expression of GA and ABA metabolism genes. Red and far-red light have no effect on germination, ABA levels, or HvNCED1. RNA interference experiments in transgenic barley plants support a role of HvABA8′OH1 in dormancy release. Reduced HvABA8′OH1 expression in transgenic HvABA8′OH1 RNAi grains results in higher levels of ABA and increased dormancy compared to nontransgenic grains.

Efficient silencing of endogenous microRNAs using artificial microRNAs in Arabidopsis thaliana

We report here that the expression of endogenous microRNAs (miRNAs) can be efficiently silenced in Arabidopsis thaliana (Arabidopsis) using artificial miRNA (amiRNA) technology. We demonstrate that an amiRNA designed to target a mature miRNA directs silencing against all miRNA family members, whereas an amiRNA designed to target the stem-loop region of a miRNA precursor transcript directs silencing against only the individual family member targeted. Furthermore, our results indicate that amiRNAs targeting both the mature miRNA and stem-loop sequence direct RNA silencing through cleavage of the miRNA precursor transcript, which presumably occurs in the nucleus of a plant cell during the initial stages of miRNA biogenesis. This suggests that small RNA (sRNA)-guided RNA cleavage in plants occurs not only in the cytoplasm, but also in the nucleus. Many plant miRNA gene families have been identified via sequencing and bioinformatic analysis, but, to date, only a small tranche of these have been functionally characterized due to a lack of effective forward or reverse genetic tools. Our findings therefore provide a new and powerful reverse-genetic tool for the analysis of miRNA function in plants. © The Author 2010. Published by the Molecular Plant Shanghai Editorial Office in association with Oxford University Press on behalf of CSPP and IPPE, SIBS, CAS.

Expression of Caenorhabditis elegans RNA-directed RNA polymerase in transgenic Drosophila melanogaster does not affect morphological development

Drosophila melanogaster, along with all insects and the vertebrates, lacks an RdRp gene. We created transgenic strains of Drosophila melanogaster in which the rrf-1 or ego-1 RdRp genes from C. elegans were placed under the control of the yeast GAL4 upstream activation sequence. Activation of the gene was performed by crossing these lines to flies carrying the GAL4 transgene under the control of various Drosophila enhancers. RT-PCR confirmed the successful expression of each RdRp gene. The resulting phenotypes indicated that introduction of the RdRp genes had no effect on D. melanogaster morphological development. © 2010 Springer Science+Business Media B.V.

The complete nucleotide sequence of subterranean clover mottle virus

The complete nucleotide sequence of Subterranean clover mottle virus (SCMoV) genomic RNA has been determined. The SCMoV genome is 4,258 nucleotides in length. It shares most nucleotide and amino acid sequence identity with the genome of Lucerne transient streak virus (LTSV). SCMoV RNA encodes four overlapping open reading frames and has a genome organisation similar to that of Cocksfoot mottle virus (CfMV). ORF1 and ORF4 are predicted to encode single proteins. ORF2 is predicted to encode two proteins that are derived from a -1 translational frameshift between two overlapping reading frames (ORF2a and ORF2b). A search of amino acid databases did not find a significant match for ORF1 and the function of this protein remains unclear. ORF2a contains a motif typical of chymotrypsin-like serine proteases and ORF2b has motifs characteristically present in positive-stranded RNA-dependent RNA polymerases. ORF4 is likely to be expressed from a subgenomic RNA and encodes the viral coat protein. The ORF2a/ORF2b overlapping gene expression strategy used by SCMoV and CfMV is similar to that of the poleroviruses and differ from that of other published sobemoviruses. These results suggest that the sobemoviruses could now be divided into two distinct subgroups based on those that express the RNA-dependent RNA polymerase from a single, in-frame polyprotein, and those that express it via a -1 translational frameshifting mechanism.

The M(r) 43K major capsid protein of rice ragged stunt oryzavirus is a post-translationally processed product of a M(r) 67 348 polypeptide encoded by genome segment 8

The nucleotide sequence of DNA complementary to rice ragged stunt oryzavirus (RRSV) genome segment 8 (S8) of an isolate from Thailand was determined. RRSV S8 is 1 914 bp in size and contains a single large open reading frame (ORF) spanning nucleotides 23 to 1 810 which is capable of encoding a protein of M(r) 67 348. The N-terminal amino acid sequence of a ~43K virion polypeptide matched to that inferred for an internal region of the S8 coding sequence. These data suggest that the 43K protein is encoded by S8 and is derived by a proteolytic cleavage. Predicted polypeptide sizes from this possible cleavage of S8 protein are 26K and 42K. Polyclonal antibodies raised against a maltose binding protein (MBP)-S8 fusion polypeptide (expressed in Escherichia coli) recognised four RRSV particle associated polypeptides of M(r) 67K, 46K, 43K and 26K and all except the 26K polypeptide were also highly immunoreactive to polyclonal antibodies raised against purified RRSV particles. Cleavage of the MBP-S8 fusion polypeptide with protease Factor X produced the expected 40K MBP and two polypeptides of apparent M(r) 46K and 26K. Antibodies to purified RRSV particles reacted strongly with the intact fusion protein and the 46K cleavage product but weakly to the 26K product. Furthermore, in vitro transcription and translation of the S8 coding region revealed a post-translational self cleavage of the 67K polypeptide to 46K and 26K products. These data indicate that S8 encodes a structural polypeptide, the majority of which is auto- catalytically cleaved to 26K and 46K proteins. The data also suggest that the 26K protein is the self cleaving protease and that the 46K product is further processed or undergoes stable conformational changes to a ~43K major capsid protein.

Nucleotide sequence of coat protein gene for GPV isolate of barley yellow dwarf virus and construction of expression plasmid for plant

GPV is a Chinese serotype isolate of barley yellow dwarf virus (BYDV) that has no reaction with antiserum of MAV, PAV, SGV, RPV and RMV The sequence of the coat protein (CP) of GPV isolate of BYDV was identified and its amino acid sequence was deduced. The coding region for the putative GPV CP is 603 bases nucleotides and encodes a Mr 22 218 (22 ku) protein. The same as MAV, PAV and RPV, GPV contained a second ORF within the coat protein coding region. This protein of 17 024 Mr (17 ku) is thought to correspond to the Virion protein genome linked (Vpg). Sequence comparisons of the CP coding region between the GPV isolate of BYDV and other isolates of BYDV have been done. The nucleotide and amino acid sequence homology of GPV has a greater identity to the sequence of RPV than those of PAV and MAV. The GPV CP sequence stored 83.7% of nucleotide similarity and 77.5% of deduced amino acid similarity, whereas that of the PAV and MAV shared 56.9%, 53.2% and 44.1%, 43.8% respectively. According to BYDV-GPV CP sequence, two primers were designed. The cDNA of CP was produced by RT-PCR. Full-length cDNA of CP was inserted into plasmid to construct expression plasmids named pPPI1, pPPI2 and pPPI5 based on different promoters. The recombinant plasmids were identified by using α-32P-dATP labelled CP probe, α-32P-ATP labelled GPV RNA probe and sequencing to confirm real GPV CP gene cDNA in plasmids.

Nucleotide sequences of an Australian and a Canadian isolate of potato leafroll luteovirus and their relationships with two European isolates

The genomes of an Australian and a Canadian isolate of potato leafroll virus have been cloned and sequenced. The sequences of both isolates are similar (about 93%), but the Canadian isolate (PLRV-C) is more closely related (about 98% identity) to a Scottish (PLRV-S) and a Dutch isolate (PLRV-N) than to the Australian isolate (PLRV-A). The 5'-terminal 18 nucleotide residues of PLRV-C, PLRV-A, PLRV-N and beet western yellows virus have 17 residues in common. In contrast, PLRV-S shows no obvious similarity in this region. PLRV-A and PLRV-C genomic sequences have localized regions of marked diversity, in particular a 600 nucleotide residue sequence in the polymerase gene. These data provide a world-wide perspective on the molecular biology of PLRV strains and their comparison with other luteoviruses and related RNA plant viruses suggests that there are two major subgroups in the plant luteoviruses.