Bruce W. Eckloff

Human catechol O-methyltransferase genetic variation: gene resequencing and functional characterization of variant allozymes

Catechol O-methyltransferase (COMT) plays an important role in the metabolism of catecholamines, catecholestrogens and catechol drugs. A common COMT G472A genetic polymorphism (Val108/158Met) that was identified previously is associated with decreased levels of enzyme activity and has been implicated as a possible risk factor for neuropsychiatric disease. We set out to ‘resequence’ the human COMT gene using DNA samples from 60 African-American and 60 Caucasian-American subjects. A total of 23 single nucleotide polymorphisms (SNPs), including a novel nonsynonymous cSNP present only in DNA from African-American subjects, and one insertion/deletion were observed. The wild type (WT) and two variant allozymes, Thr52 and Met108, were transiently expressed in COS-1 and HEK293 cells. There was no significant change in level of COMT activity for the Thr52 variant allozyme, but there was a 40% decrease in the level of activity in cells transfected with the Met108 construct. Apparent Km values of the WT and variant allozymes for the two reaction cosubstrates differed slightly, but significantly, for 3,4-dihydroxybenzoic acid but not for S-adenosyl-L-methionine. The Met108 allozyme displayed a 70–90% decrease in immunoreactive protein when compared with WT, but there was no significant change in the level of immunoreactive protein for Thr52. A significant decrease in the level of immunoreactive protein was also observed in hepatic biopsy samples from patients homozygous for the allele encoding Met108. These observations represent steps toward an understanding of molecular genetic mechanisms responsible for variation in COMT level and/or properties, variation that may contribute to the pathophysiology of neuropsychiatric disease.

The genomic landscape of small intestine neuroendocrine tumors

Small intestine neuroendocrine tumors (SI-NETs) are the most common malignancy of the small bowel. Several clinical trials target PI3K/Akt/mTOR signaling; however, it is unknown whether these or other genes are genetically altered in these tumors. To address the underlying genetics, we analyzed 48 SI-NETs by massively parallel exome sequencing. We detected an average of 0.1 somatic single nucleotide variants (SNVs) per 106 nucleotides (range, 0-0.59), mostly transitions (C>T and A>G), which suggests that SI-NETs are stable cancers. 197 protein-altering somatic SNVs affected a preponderance of cancer genes, including FGFR2, MEN1, HOOK3, EZH2, MLF1, CARD11, VHL, NONO, and SMAD1. Integrative analysis of SNVs and somatic copy number variations identified recurrently altered mechanisms of carcinogenesis: chromatin remodeling, DNA damage, apoptosis, RAS signaling, and axon guidance. Candidate therapeutically relevant alterations were found in 35 patients, including SRC, SMAD family genes, AURKA, EGFR, HSP90, and PDGFR. Mutually exclusive amplification of AKT1 or AKT2 was the most common event in the 16 patients with alterations of PI3K/Akt/mTOR signaling. We conclude that sequencing-based analysis may provide provisional grouping of SI-NETs by therapeutic targets or deregulated pathways.

Gemcitabine Pharmacogenomics: Cytidine Deaminase and Deoxycytidylate Deaminase Gene Resequencing and Functional Genomics

Purpose: Gemcitabine is a nucleoside analogue with activity against solid tumors. Gemcitabine metabolic inactivation is catalyzed by cytidine deaminase (CDA) or, after phosphorylation, by deoxycytidylate deaminase (DCTD). We set out to study the pharmacogenomics of CDA and DCTD. Experimental Design: The genes encoding CDA and DCTD were resequenced using DNA from 60 African American and 60 Caucasian American subjects. Expression constructs were created for nonsynonymous coding single nucleotide polymorphisms (cSNP) and reporter gene constructs were created for 5′-flanking region polymorphisms. Functional genomic studies were then conducted after the transfection of mammalian cells. Results:CDA resequencing revealed 17 polymorphisms, including one common nonsynonymous cSNP, 79 A>C (Lys27Gln). Recombinant Gln27 CDA had 66 ± 5.1% (mean ± SE) of the wild-type (WT) activity for gemcitabine but without a significant decrease in level of immunoreactive protein. The apparent Km (397 ± 40 μmol/L) for the Gln27 allozyme was significantly higher than that for the WT (289 ± 20 μmol/L; P < 0.025). CDA 5′-flanking region reporter gene studies showed significant differences among 5′-flanking region haplotypes in their ability to drive transcription. There were 29 SNPs in DCTD, including one nonsynonymous cSNP, 172 A>G (Asn58Asp), in Caucasian American DNA. Recombinant Asp58 DCTD had 11 ± 1.4% of WT activity for gemcitabine monophosphate with a significantly elevated level of immunoreactive protein. No DCTD polymorphisms were observed in the initial 500 bp of the 5′-flanking region. Conclusions: These results suggest that pharmacogenomic variation in the deamination of gemcitabine and its monophosphate might contribute to variation in therapeutic response to this antineoplastic agent.

Genome-wide analysis reveals recurrent structural abnormalities of TP63 and other p53-related genes in peripheral T-cell lymphomas

Peripheral T-cell lymphomas (PTCLs) are aggressive malignancies of mature T lymphocytes with 5-year overall survival rates of only ∼ 35%. Improvement in outcomes has been stymied by poor understanding of the genetics and molecular pathogenesis of PTCL, with a resulting paucity of molecular targets for therapy. We developed bioinformatic tools to identify chromosomal rearrangements using genome-wide, next-generation sequencing analysis of mate-pair DNA libraries and applied these tools to 16 PTCL patient tissue samples and 6 PTCL cell lines. Thirteen recurrent abnormalities were identified, of which 5 involved p53-related genes (TP53, TP63, CDKN2A, WWOX, and ANKRD11). Among these abnormalities were novel TP63 rearrangements encoding fusion proteins homologous to ΔNp63, a dominant-negative p63 isoform that inhibits the p53 pathway. TP63 rearrangements were seen in 11 (5.8%) of 190 PTCLs and were associated with inferior overall survival; they also were detected in 2 (1.2%) of 164 diffuse large B-cell lymphomas. As TP53 mutations are rare in PTCL compared with other malignancies, our findings suggest that a constellation of alternate genetic abnormalities may contribute to disruption of p53-associated tumor suppressor function in PTCL.

Human Arsenic Methyltransferase (AS3MT) Pharmacogenetics GENE RESEQUENCING AND FUNCTIONAL GENOMICS STUDIES

Arsenic contaminates ground water worldwide. Methylation is an important reaction in the biotransformation of arsenic. We set out to study the pharmacogenetics of human arsenic methyltransferase (AS3MT, previously CYT19). After cloning the human AS3MT cDNA, we annotated the human gene and resequenced its 5′-flanking region, exons, and splice junctions using 60 DNA samples from African-American (AA) and 60 samples from Caucasian-American (CA) subjects. We observed 26 single nucleotide polymorphisms (SNPs), including 3 non-synonymous cSNPs, as well as a variable number of tandem repeats in exon 1 within an area encoding the cDNA 5′-untranslated region. The nonsynonymous cSNPs included T860C (M287T) with frequencies of 10.8 and 10% in AA and CA subjects, respectively, as well as C517T (A173W) in one AA and C917T (T306I) in one CA sample. Haplotype analysis showed that Ile306 was linked to Thr287, so this double variant allozyme was also studied functionally. After expression in COS-1 cells and correction for transfection efficiency, the Trp173 allozyme displayed 31%, Thr287 350%, Ile306 4.8%, and Thr287/Ile306 6.2% of the activity of the wild type (WT) allozyme, with 20, 190, 4.4, and 7.9% of the level of WT immunoreactive protein, respectively. Apparent Km values for S-adenosyl-l-methionine were 4.6, 3.1, and 11 μm for WT, Trp173, and Thr287 allozymes, with Km values for sodium arsenite with the same allozymes of 11.8, 8.9, and 4.5μm. The Ile306 and Thr287/Ile306 allozymes expressed too little activity for inclusion in the substrate kinetic studies. Expression of reporter gene constructs for the 5′-flanking region and the variable number of tandem repeats in the 5′-untranslated region demonstrated cell line-dependent variation in reporter gene expression, with shorter repeats associated with increased transcription in HepG2 cells. These results raise the possibility that inherited variation in AS3MT may contribute to variation in arsenic metabolism and, perhaps, arsenic-dependent carcinogenesis in humans.

Identification of Gene Mutations in Autosomal Dominant Polycystic Kidney Disease through Targeted Resequencing

Mutations in two large multi-exon genes, PKD1 and PKD2, cause autosomal dominant polycystic kidney disease (ADPKD). The duplication of PKD1 exons 1-32 as six pseudogenes on chromosome 16, the high level of allelic heterogeneity, and the cost of Sanger sequencing complicate mutation analysis, which can aid diagnostics of ADPKD. We developed and validated a strategy to analyze both the PKD1 and PKD2 genes using next-generation sequencing by pooling long-range PCR amplicons and multiplexing bar-coded libraries. We used this approach to characterize a cohort of 230 patients with ADPKD. This process detected definitely and likely pathogenic variants in 115 (63%) of 183 patients with typical ADPKD. In addition, we identified atypical mutations, a gene conversion, and one missed mutation resulting from allele dropout, and we characterized the pattern of deep intronic variation for both genes. In summary, this strategy involving next-generation sequencing is a model for future genetic characterization of large ADPKD populations.

Human sulfotransferase SULT2A1 pharmacogenetics: genotype-to-phenotype studies

SULT2A1 catalyzes the sulfate conjugation of dehydroepiandrosterone (DHEA) as well as other steroids. As a step toward pharmacogenetic studies, we have ‘resequenced’ SULT2A1 using 60 DNA samples from African-American and 60 samples from Caucasian-American subjects. All exons, splice junctions and approximately 370 bp located 5′ of the site of transcription initiation were sequenced. We observed 15 single nucleotide polymorphisms (SNPs), including three non-synonymous coding SNPs (cSNPs) that were present only in DNA from African-American subjects. Linkage analysis revealed that two of the nonsynonymous cSNPs were tightly linked. Expression constructs were created for all nonsynonymous cSNPs observed, including a ‘double variant’ construct that included the two linked cSNPs, and those constructs were expressed in COS-1 cells. SULT2A1 activity was significantly decreased for three of the four variant allozymes. Western blot analysis demonstrated that decreased levels of immunoreactive protein appeared to be the major mechanism responsible for decreases in activity, although apparent Km values also varied among the recombinant allozymes. In addition, the most common of the nonsynonymous cSNPs disrupted the portion of SULT2A1 involved with dimerization, and this variant allozyme behaved as a monomer rather than a dimer during gel filtration chromatography. These observations indicate that common genetic polymorphisms for SULT2A1 can result in reductions in levels of both activity and enzyme protein. They also raise the possibility of ethnic-specific pharmacogenetic variation in SULT2A1-catalyzed sulfation of both endogenous and exogenous substrates for this phase II drug-metabolizing enzyme.

Human estrogen sulfotransferase (SULT1E1) pharmacogenomics: gene resequencing and functional genomics

Estrogens are used as drugs and estrogen exposure is a risk factor for hormone‐dependent diseases such as breast cancer. Sulfate conjugation is an important pathway for estrogen metabolism. The sulfotransferase (SULT) enzyme SULT1E1 has the lowest Km values for estrogens and catecholestrogens of the 10 known human SULT isoforms. We previously cloned and characterized the human SULT1E1 cDNA and gene as steps toward pharmacogenetic studies. In the present experiments, we set out to determine whether common, functionally significant genetic polymorphisms might exist for SULT1E1. As a first step, we ‘resequenced’ the eight SULT1E1 exons and exon–intron splice junctions as well as portions of the 5′‐flanking region using DNA from 60 African‐American and 60 Caucasian‐American subjects. In all, 23 polymorphisms, 22 single nucleotide polymorphisms (SNPs) and one insertion deletion were observed. There were three nonsynonymous coding SNPs (cSNPs) that altered the following encoded amino acids: Asp22Tyr, Ala32Val and Pro253His. Among these, 12 pairs of SNPs were tightly linked. In addition, 12 unambiguous SULT1E1 haplotypes were identified, including six that were common to both populations studied. Transient expression in COS‐1 cells of constructs containing the three nonsynonymous cSNPs showed significant decreases in SULT1E1 activity for the Tyr22 and Val32 allozymes, with corresponding decreases in levels of immunoreactive protein. There were no changes in levels of either activity or immunoreactive protein for the His253 allozyme. Apparent Km values of the Val32 allozyme for the two cosubstrates for the reaction, 17β‐estradiol and 3′‐phosphoadenosine 5′‐phosphosulfate, were not significantly different from those of the wild‐type enzyme, but there was a two‐ to three‐fold increase in Km values for the His253 allozyme and a greater than five‐fold increase for the Tyr22 allozyme. These observations raise the possibility that genetically determined variation in SULT1E1‐catalyzed estrogen sulfation might contribute to the pathophysiology of estrogen‐dependent diseases as well as variation in the biotransformation of exogenously administered estrogens.

GLUTATHIONE S-TRANSFERASE OMEGA 1 AND OMEGA 2 PHARMACOGENOMICS

Glutathione S-transferase omega 1 and omega 2 (GSTO1 and GSTO2) catalyze monomethyl arsenate reduction, the rate-limiting reaction in arsenic biotransformation. As a step toward pharmacogenomic studies of these phase II enzymes, we resequenced human GSTO1 and GSTO2 using DNA samples from four ethnic groups. We identified 31 and 66 polymorphisms in GSTO1 and GSTO2, respectively, with four nonsynonymous-coding single nucleotide polymorphisms (cSNPs) in each gene. There were striking variations among ethnic groups in polymorphism frequencies and types. Expression constructs were created for all eight nonsynonymous cSNPs, as well as a deletion of codon 155 in GSTO1, and those constructs were used to transfect COS-1 cells. Quantitative Western blot analysis, after correction for transfection efficiency, showed a reduction in protein level of greater than 50% for the GSTO1 Tyr32 variant allozyme compared with wild type (WT), whereas levels for the Asp140, Lys208, Val236, and codon 155 deletion variant constructs were similar to that of the WT. For GSTO2, the Tyr130 and Ile158 variant allozymes showed 50 and 84% reductions in levels of expression, respectively, compared with WT, whereas the Ile41 and Asp142 allozymes displayed levels similar to that of WT GSTO2. Rabbit reticulocyte lysate degradation studies showed that the GSTO1 Tyr32 and the GSTO2 Tyr130, Ile158, and Asp142/Ile158 variant allozymes were degraded more rapidly than were their respective WT allozymes. These observations raise the possibility of functionally significant pharmacogenomic variation in the expression and function of GSTO1 and GSTO2.

Multi-Platform Analysis of MicroRNA Expression Measurements in RNA from Fresh Frozen and FFPE Tissues

MicroRNAs play a role in regulating diverse biological processes and have considerable utility as molecular markers for diagnosis and monitoring of human disease. Several technologies are available commercially for measuring microRNA expression. However, cross-platform comparisons do not necessarily correlate well, making it difficult to determine which platform most closely represents the true microRNA expression level in a tissue. To address this issue, we have analyzed RNA derived from cell lines, as well as fresh frozen and formalin-fixed paraffin embedded tissues, using Affymetrix, Agilent, and Illumina microRNA arrays, NanoString counting, and Illumina Next Generation Sequencing. We compared the performance within- and between the different platforms, and then verified these results with those of quantitative PCR data. Our results demonstrate that the within-platform reproducibility for each method is consistently high and although the gene expression profiles from each platform show unique traits, comparison of genes that were commonly detectable showed that detection of microRNA transcripts was similar across multiple platforms.