Recurrent MSC E116K mutations in ALK-negative anaplastic large cell lymphoma.

Anaplastic large cell lymphomas (ALCLs) represent a relatively common group of T-cell non-Hodgkin lymphomas (T-NHLs) that are unified by similar pathologic features but demonstrate marked genetic heterogeneity. ALCLs are broadly classified as being anaplastic lymphoma kinase (ALK)+ or ALK-, based on the presence or absence of ALK rearrangements. Exome sequencing of 62 T-NHLs identified a previously unreported recurrent mutation in the musculin gene, MSC E116K, exclusively in ALK- ALCLs. Additional sequencing for a total of 238 T-NHLs confirmed the specificity of MSC E116K for ALK- ALCL and further demonstrated that 14 of 15 mutated cases (93%) had coexisting DUSP22 rearrangements. Musculin is a basic helix-loop-helix (bHLH) transcription factor that heterodimerizes with other bHLH proteins to regulate lymphocyte development. The E116K mutation localized to the DNA binding domain of musculin and permitted formation of musculin-bHLH heterodimers but prevented their binding to authentic target sequence. Functional analysis showed MSCE116K acted in a dominant-negative fashion, reversing wild-type musculin-induced repression of MYC and cell cycle inhibition. Chromatin immunoprecipitation-sequencing and transcriptome analysis identified the cell cycle regulatory gene E2F2 as a direct transcriptional target of musculin. MSCE116K reversed E2F2-induced cell cycle arrest and promoted expression of the CD30-IRF4-MYC axis, whereas its expression was reciprocally induced by binding of IRF4 to the MSC promoter. Finally, ALCL cells expressing MSC E116K were preferentially targeted by the BET inhibitor JQ1. These findings identify a novel recurrent MSC mutation as a key driver of the CD30-IRF4-MYC axis and cell cycle progression in a unique subset of ALCLs.

Next-Generation Sequencing of CYP2C19 in Stent Thrombosis: Implications for Clopidogrel Pharmacogenomics.

PURPOSE: Describe CYP2C19 sequencing results in the largest series of clopidogrel-treated cases with stent thrombosis (ST), the closest clinical phenotype to clopidogrel resistance. Evaluate the impact of CYP2C19 genetic variation detected by next-generation sequencing (NGS) with comprehensive annotation and functional studies. METHODS: Seventy ST cases on clopidogrel identified from the PLATO trial (n = 58) and Mayo Clinic biorepository (n = 12) were matched 1:1 with controls for age, race, sex, diabetes mellitus, presentation, and stent type. NGS was performed to cover the entire CYP2C19 gene. Assessment of exonic variants involved measuring in vitro protein expression levels. Intronic variants were evaluated for potential splicing motif variations. RESULTS: Poor metabolizers (n = 4) and rare CYP2C19*8, CYP2C19*15, and CYP2C19*11 alleles were identified only in ST cases. CYP2C19*17 heterozygote carriers were observed more frequently in cases (n = 29) than controls (n = 18). Functional studies of CYP2C19 exonic variants (n = 11) revealed 3 cases and only 1 control carrying a deleterious variant as determined by in vitro protein expression studies. Greater intronic variation unique to ST cases (n = 169) compared with controls (n = 84) was observed with predictions revealing 13 allele candidates that may lead to a potential disruption of splicing and a loss-of-function effect of CYP2C19 in ST cases. CONCLUSION: NGS detected CYP2C19 poor metabolizers and paradoxically greater number of so-called rapid metabolizers in ST cases. Rare deleterious exonic variation occurs in 4%, and potentially disruptive intronic alleles occur in 16% of ST cases. Additional studies are required to evaluate the role of these variants in platelet aggregation and clopidogrel metabolism.

Integrated mate-pair and RNA sequencing identifies novel, targetable gene fusions in peripheral T-cell lymphoma.

Peripheral T-cell lymphomas (PTCLs) represent a heterogeneous group of T-cell malignancies that generally demonstrate aggressive clinical behavior, often are refractory to standard therapy, and remain significantly understudied. The most common World Health Organization subtype is PTCL, not otherwise specified (NOS), essentially a "wastebasket" category because of inadequate understanding to assign cases to a more specific diagnostic entity. Identification of novel fusion genes has contributed significantly to improving the classification, biologic understanding, and therapeutic targeting of PTCLs. Here, we integrated mate-pair DNA and RNA next-generation sequencing to identify chromosomal rearrangements encoding expressed fusion transcripts in PTCL, NOS. Two of 11 cases had novel fusions involving VAV1, encoding a truncated form of the VAV1 guanine nucleotide exchange factor important in T-cell receptor signaling. Fluorescence in situ hybridization studies identified VAV1 rearrangements in 10 of 148 PTCLs (7%). These were observed exclusively in PTCL, NOS (11%) and anaplastic large cell lymphoma (11%). In vitro, ectopic expression of a VAV1 fusion promoted cell growth and migration in a RAC1-dependent manner. This growth was inhibited by azathioprine, a clinically available RAC1 inhibitor. We also identified novel kinase gene fusions, ITK-FER and IKZF2-ERBB4, as candidate therapeutic targets that show similarities to known recurrent oncogenic ITK-SYK fusions and ERBB4 transcript variants in PTCLs, respectively. Additional novel and potentially clinically relevant fusions also were discovered. Together, these findings identify VAV1 fusions as recurrent and targetable events in PTCLs and highlight the potential for clinical sequencing to guide individualized therapy approaches for this group of aggressive malignancies.

Homozygous/Compound Heterozygous Triadin Mutations Associated With Autosomal-Recessive Long-QT Syndrome and Pediatric Sudden Cardiac Arrest: Elucidation of the Triadin Knockout Syndrome.

BACKGROUND: Long-QT syndrome (LQTS) may result in syncope, seizures, or sudden cardiac arrest. Although 16 LQTS-susceptibility genes have been discovered, 20% to 25% of LQTS remains genetically elusive. METHODS AND RESULTS: We performed whole-exome sequencing child-parent trio analysis followed by recessive and sporadic inheritance modeling and disease-network candidate analysis gene ranking to identify a novel underlying genetic mechanism for LQTS. Subsequent mutational analysis of the candidate gene was performed with polymerase chain reaction, denaturing high-performance liquid chromatography, and DNA sequencing on a cohort of 33 additional unrelated patients with genetically elusive LQTS. After whole-exome sequencing and variant filtration, a homozygous p.D18fs*13 TRDN-encoded triadin frameshift mutation was discovered in a 10-year-old female patient with LQTS with a QTc of 500 milliseconds who experienced recurrent exertion-induced syncope/cardiac arrest beginning at 1 year of age. Subsequent mutational analysis of TRDN revealed either homozygous or compound heterozygous frameshift mutations in 4 of 33 unrelated cases of LQTS (12%). All 5 TRDN-null patients displayed extensive T-wave inversions in precordial leads V1 through V4, with either persistent or transient QT prolongation and severe disease expression of exercise-induced cardiac arrest in early childhood (≤3 years of age) and required aggressive therapy. The overall yield of TRDN mutations was significantly greater in patients ≤10 years of age (5 of 10, 50%) compared with older patients (0 of 24, 0%; P=0.0009). CONCLUSIONS: We identified TRDN as a novel underlying genetic basis for recessively inherited LQTS. All TRDN-null patients had strikingly similar phenotypes. Given the recurrent nature of potential lethal arrhythmias, patients fitting this phenotypic profile should undergo cardiac TRDN genetic testing.

Genetic alterations associated with progression from pancreatic intraepithelial neoplasia to invasive pancreatic tumor.

BACKGROUND & AIMS: Increasing grade of pancreatic intraepithelial neoplasia (PanIN) has been associated with progression to pancreatic ductal adenocarcinoma (PDAC). However, the mechanisms that control progression from PanINs to PDAC are not well understood. We investigated the genetic alterations involved in this process. METHODS: Genomic DNA samples from laser-capture microdissected PDACs and adjacent PanIN2 and PanIN3 lesions from 10 patients with pancreatic cancer were analyzed by exome sequencing. RESULTS: Similar numbers of somatic mutations were identified in PanINs and tumors, but the mutational load varied greatly among cases. Ten of the 15 isolated PanINs shared more than 50% of somatic mutations with associated tumors. Mutations common to tumors and clonally related PanIN2 and PanIN3 lesions were identified as genes that could promote carcinogenesis. KRAS and TP53 frequently were altered in PanINs and tumors, but few other recurrently modified genes were detected. Mutations in DNA damage response genes were prevalent in all samples. Genes that encode proteins involved in gap junctions, the actin cytoskeleton, the mitogen-activated protein kinase signaling pathway, axon guidance, and cell-cycle regulation were among the earliest targets of mutagenesis in PanINs that progressed to PDAC. CONCLUSIONS: Early stage PanIN2 lesions appear to contain many of the somatic gene alterations required for PDAC development.

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.

Novel TRAF1-ALK fusion identified by deep RNA sequencing of anaplastic large cell lymphoma.

Chromosomal translocations leading to expression of abnormal fusion proteins play a major role in the pathogenesis of various hematologic malignancies. The recent development of high-throughput, "deep" sequencing has allowed discovery of novel translocations leading to a rapid increase in understanding these diseases. Translocations involving the anaplastic lymphoma kinase (ALK) gene leading to ALK fusion proteins originally were discovered in anaplastic large cell lymphomas (ALCLs). Among ALCLs, NPM1-ALK fusions are most common and lead to nuclear localization of the fusion protein. Here, we present a 50-year-old male with ALCL demonstrating cytoplasmic ALK immunoreactivity only, suggesting the presence of a non-NPM1 fusion partner. We performed deep RNA sequencing of tumor tissue from this patient and identified a novel transcript fusing Exon 6 of TRAF1 to Exon 20 of ALK. The TRAF1-ALK fusion transcript was confirmed at the mRNA level by Sanger sequencing and the fusion protein was visualized by Western blot. The discovery of this TRAF1-ALK fusion expands the diversity of known ALK fusion partners and highlights the power of deep sequencing for fusion transcript discovery. © 2013 Wiley Periodicals, Inc.

FKBP5 genetic variation: association with selective serotonin reuptake inhibitor treatment outcomes in major depressive disorder.

OBJECTIVES: FKBP51 (51 kDa immunophilin) acts as a modulator of the glucocorticoid receptor and a negative regulator of the Akt pathway. Genetic variation in FKBP5 plays a role in antidepressant response. The aim of this study was to comprehensively assess the role of genetic variation in FKBP5, identified by both Sanger and Next Generation DNA resequencing, as well as genome-wide single nucleotide polymorphisms (SNPs) associated with FKBP5 expression in the response to the selective serotonin reuptake inhibitor (SSRI) treatment of major depressive disorder. METHODS: We identified 657 SNPs in FKBP5 by Next Generation sequencing of 96 DNA samples from white patients, and 149 SNPs were selected for the genotyping together with 235 SNPs that were trans-associated with variation in FKBP5 expression in lymphoblastoid cells. A total of 529 DNA samples from the Mayo Clinic PGRN-SSRI Pharmacogenomic trial for which genome-wide SNPs had already been obtained were genotyped for these 384 SNPs, and associations with treatment outcomes were determined. The most significant SNPs were genotyped using 96 DNA samples from white non-Hispanic patients of the NIMH-supported Sequenced Treatment Alternatives to Relieve Depression (STAR*D) study to attempt replication, followed by functional genomic studies. RESULTS: Genotype-phenotype association analysis indicated that rs352428 was associated with both 8-week treatment response in the Mayo study (odds ratio=0.49; P=0.003) and 6-week response in the STAR*D replication study (odds ratio=0.74; P=0.05). The electrophoresis mobility shift assay and the reporter gene assay confirmed the possible role of this SNP in transcription regulation. CONCLUSION: This comprehensive FKBP5 sequence study provides insight into the role of common genetic polymorphisms that might influence SSRI treatment outcomes in major depressive disorder patients.

Discovery of recurrent t(6;7)(p25.3;q32.3) translocations in ALK-negative anaplastic large cell lymphomas by massively parallel genomic sequencing.

The genetics of peripheral T-cell lymphomas are poorly understood. The most well-characterized abnormalities are translocations involving ALK, occurring in approximately half of anaplastic large cell lymphomas (ALCLs). To gain insight into the genetics of ALCLs lacking ALK translocations, we combined mate-pair DNA library construction, massively parallel ("Next Generation") sequencing, and a novel bioinformatic algorithm. We identified a balanced translocation disrupting the DUSP22 phosphatase gene on 6p25.3 and adjoining the FRA7H fragile site on 7q32.3 in a systemic ALK-negative ALCL. Using fluorescence in situ hybridization, we demonstrated that the t(6;7)(p25.3;q32.3) was recurrent in ALK-negative ALCLs. Furthermore, t(6;7)(p25.3;q32.3) was associated with down-regulation of DUSP22 and up-regulation of MIR29 microRNAs on 7q32.3. These findings represent the first recurrent translocation reported in ALK-negative ALCL and highlight the utility of massively parallel genomic sequencing to discover novel translocations in lymphoma and other cancers.

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.

Relationship between BCL2 mutations and follicular lymphoma outcome in the chemoimmunotherapy era.

How to identify follicular lymphoma (FL) patients with low disease burden but high risk for early progression is unclear. Building on a prior study demonstrating the early transformation of FLs with high variant allele frequency (VAF) BCL2 mutations at activation-induced cytidine deaminase (AICDA) sites, we examined 11 AICDA mutational targets, including BCL2, BCL6, PAX5, PIM1, RHOH, SOCS, and MYC, in 199 newly diagnosed grade 1 and 2 FLs. BCL2 mutations with VAF ≥20% occurred in 52% of cases. Among 97 FL patients who did not initially receive rituximab-containing therapy, nonsynonymous BCL2 mutations at VAF ≥20% were associated with increased transformation risk (HR 3.01, 95% CI 1.04-8.78, p = 0.043) and a trend toward shorter event-free survival (EFS, median 20 months with mutations versus 54 months without, p = 0.052). Other sequenced genes were less frequently mutated and did not increase the prognostic value of the panel. Across the entire population, nonsynonymous BCL2 mutations at VAF ≥20% were associated with decreased EFS (HR 1.55, 95% CI 1.02-2.35, p = 0.043 after correction for FLIPI and treatment) and decreased overall survival after median 14-year follow-up (HR 1.82, 95% CI 1.05-3.17, p = 0.034). Thus, high VAF nonsynonymous BCL2 mutations remain prognostic even in the chemoimmunotherapy era.

Serine hydroxymethyltransferase 1 and 2: gene sequence variation and functional genomic characterization.

Serine hydroxymethyltransferase (SHMT) catalyzes the transfer of a ß-carbon from serine to tetrahydrofolate to form glycine and 5,10-methylene-tetrahydrofolate. This reaction plays an important role in neurotransmitter synthesis and metabolism. We set out to resequence SHMT1 and SHMT2, followed by functional genomic studies. We identified 87 and 60 polymorphisms in SHMT1 and SHMT2, respectively. We observed no significant functional effect of the 13 non-synonymous single-nucleotide polymorphism (SNPs) in these genes, either on catalytic activity or protein quantity. We imputed additional variants across the two genes using '1000 Genomes' data, and identified 14 variants that were significantly associated (p<1.0E-10) with SHMT1 messenger RNA expression in lymphoblastoid cell lines. Many of these SNPs were also significantly correlated with basal SHMT1 protein expression in 268 human liver biopsy samples. Reporter gene assays suggested that the SHMT1 promoter SNP, rs669340, contributed to this variation. Finally, SHMT1 and SHMT2 expression were significantly correlated with those of other Folate and Methionine Cycle genes at both the messenger RNA and protein levels. These experiments represent a comprehensive study of SHMT1 and SHMT2 gene sequence variation and its functional implications. In addition, we obtained preliminary indications that these genes may be co-regulated with other Folate and Methionine Cycle genes.

Gemcitabine metabolic pathway genetic polymorphisms and response in patients with non-small cell lung cancer.

BACKGROUND AND OBJECTIVE: Gemcitabine is widely used to treat non-small cell lung cancer (NSCLC). The aim of this study was to assess the pharmacogenomic effects of the entire gemcitabine metabolic pathway, we genotyped single nucleotide polymorphisms (SNPs) within the 17 pathway genes using DNA samples from patients with NSCLC treated with gemcitabine to determine the effect of genetic variants within gemcitabine pathway genes on overall survival (OS) of patients with NSCLC after treatment of gemcitabine. METHODS: Eight of the 17 pathway genes were resequenced with DNA samples from Coriell lymphoblastoid cell lines (LCLs) using Sanger sequencing for all exons, exon-intron junctions, and 5'-, 3'-UTRs. A total of 107 tagging SNPs were selected on the basis of the resequencing data for the eight genes and on HapMap data for the remaining nine genes, followed by successful genotyping of 394 NSCLC patient DNA samples. Association of SNPs/haplotypes with OS was performed using the Cox regression model, followed by functional studies performed with LCLs and NSCLC cell lines. RESULTS: Five SNPs in four genes (CDA, NT5C2, RRM1, and SLC29A1) showed associations with OS of those patients with NSCLC, as well as nine haplotypes in four genes (RRM1, RRM2, SLC28A3, and SLC29A1) with a P value of less than 0.05. Genotype imputation using the LCLs was performed for a region of 200 kb surrounding those SNPs, followed by association studies with gemcitabine cytotoxicity. Functional studies demonstrated that downregulation of SLC29A1, NT5C2, and RRM1 in NSCLC cell lines altered cell susceptibility to gemcitabine. CONCLUSION: These studies help in identifying biomarkers to predict gemcitabine response in NSCLC, a step toward the individualized chemotherapy of lung cancer.

Human liver methionine cycle: MAT1A and GNMT gene resequencing, functional genomics, and hepatic genotype-phenotype correlation.

The "methionine cycle" plays a critical role in the regulation of concentrations of (S)-adenosylmethionine (AdoMet), the major biological methyl donor. We set out to study sequence variation in genes encoding the enzyme that synthesizes AdoMet in liver, methionine adenosyltransferase 1A (MAT1A) and the major hepatic AdoMet using enzyme, glycine N-methyltransferase (GNMT), as well as functional implications of that variation. We resequenced MAT1A and GNMT using DNA from 288 subjects of three ethnicities, followed by functional genomic and genotype-phenotype correlation studies performed with 268 hepatic biopsy samples. We identified 44 and 42 polymorphisms in MAT1A and GNMT, respectively. Quantitative Western blot analyses for the human liver samples showed large individual variation in MAT1A and GNMT protein expression. Genotype-phenotype correlation identified two genotyped single-nucleotide polymorphisms (SNPs), reference SNP (rs) 9471976 (corrected p = 3.9 × 10(-10)) and rs11752813 (corrected p = 1.8 × 10(-5)), and 42 imputed SNPs surrounding GNMT that were significantly associated with hepatic GNMT protein levels (corrected p values < 0.01). Reporter gene studies showed that variant alleles for both genotyped SNPs resulted in decreased transcriptional activity. Correlation analyses among hepatic protein levels for methionine cycle enzymes showed significant correlations between GNMT and MAT1A (p = 1.5 × 10(-3)) and between GNMT and betaine homocysteine methyltransferase (p = 1.6 × 10(-7)). Our discovery of SNPs that are highly associated with hepatic GNMT protein expression as well as the "coordinate regulation" of methionine cycle enzyme protein levels provide novel insight into the regulation of this important human liver biochemical pathway.

RNA-Seq Reveals Differences in Expressed Tumor Mutation Burden in Colorectal and Endometrial Cancers with and without Defective DNA-Mismatch Repair.

Tumor mutation burden (TMB) is an emerging biomarker of immunotherapy response. RNA sequencing in FFPE tissue samples was used for determining TMB in microsatellite-stable (MSS) and microsatellite instability-high (MSI-H) tumors in patients with colorectal or endometrial cancer. Tissue from tumors and paired normal tissue from 46 MSI-H and 12 MSS cases were included. Of the MSI-H tumors, 29 had defective DNA mismatch-repair mutations, and 17 had MLH1 promoter hypermethylation. TMB was measured using the expressed somatic nucleotide variants (eTMB). A method of accurate measurement of eTMB was developed that removes FFPE-derived artifacts by leveraging mutation signatures. There was a significant difference in the median eTMB values observed between MSI-H and MSS cases: 27.3 versus 6.7 mutations/megabase (mut/Mb) (P = 3.5 × 10-9). Among tumors with defective DNA-mismatch repair, those with mismatch-repair mutations had a significantly higher median eTMB than those with hypermethylation: 28.1 versus 17.5 mut/Mb (P = 0.037). Multivariate analysis showed that MSI status, tumor type (endometrial or colorectal), and age were significantly associated with eTMB. Additionally, using whole-exome sequencing in a subset of these patients, it was determined that DNA TMB correlated well with eTMB (Spearman correlation coefficient, 0.83). These results demonstrate that RNA sequencing can be used for measuring eTMB in FFPE tumor specimens.

Human S-adenosylhomocysteine hydrolase: common gene sequence variation and functional genomic characterization.

S-Adenosylhomocysteine hydrolase (AHCY) is the only mammalian enzyme known to catalyze the hydrolysis of S-adenosylhomocysteine. We have used a genotype-to-phenotype strategy to study this important enzyme by resequencing AHCY in 240 DNA samples from four ethnic groups. Thirty-nine polymorphisms were identified - 28 of which were novel. Functional genomic studies for wild type AHCY and the three variant allozymes identified showed that two variant allozymes had slight, but significant decreases in enzyme activity, but with no significant differences in levels of immunoreactive protein. Luciferase reporter gene assays for common 5'-flanking region haplotypes revealed that one haplotype with a frequency of approximately 2% in Caucasian-American subjects displayed a decreased ability to drive transcription. The variant nucleotide at 5'-flanking region single nucleotide polymorphism (SNP) (-34) in that haplotype altered the DNA-protein binding pattern during electrophoresis mobility shift assay. Finally, an AHCY genotype-phenotype association study for expression in lymphoblastoid cells identified four SNPs that were associated with decreased expression. For the IVS6 (intervening sequence 6, i.e., intron 6) G56 > C SNP among those four, electrophoresis mobility shift assay showed that a C > G nucleotide change resulted in an additional shifted band. These results represent a step toward understanding the functional consequences of common genetic variation in AHCY for the regulation of neurotransmitter, drug and macromolecule methylation.

Proteasome beta subunit pharmacogenomics: gene resequencing and functional genomics.

PURPOSE: The proteasome is a multisubunit cellular organelle that functions as a nonlysosomal threonine protease. Proteasomes play a critical role in the degradation of proteins, regulating a variety of cellular processes, and they are also the target for antineoplastic proteasome inhibitors. Genetic variation in proteasome subunits could influence both proteasome function and response to drug therapy. EXPERIMENTAL DESIGN: We resequenced genes encoding the three active proteasome beta subunits using 240 DNA samples from four ethnic groups and the beta 5 subunit gene in 79 DNA samples from multiple myeloma patients who had been treated with the proteasome inhibitor bortezomib. Resequencing was followed by functional studies of polymorphisms identified in the coding region and 3'-flanking region (3'-FR) of PSMB5, the gene encoding the target for clinically useful proteasome inhibitors. RESULTS: Resequencing of 240 DNA samples identified a series of novel ethnic-specific polymorphisms that are not represented in public databases. The PSMB5 3'-FR 1042 G allele significantly increased transcription during reporter gene studies, observations confirmed by genotype-phenotype correlations between single nucleotide polymorphisms (SNP) in PSMB5 and mRNA expression in the 240 lymphoblastoid cell lines from which the resequenced DNA was obtained. Studies with patient DNA samples identified additional novel PSMB5 polymorphisms, including a SNP and an insertion in the 3'-FR. Reporter-gene studies indicated that these two novel polymorphisms might decrease transcription. CONCLUSIONS: These results show that nonsynonymous coding SNPs in the PSMB5 gene did not show significant effects on proteasome activity, but SNPs did influence transcription. Future studies might focus on regulatory region polymorphisms.

Human betaine-homocysteine methyltransferase (BHMT) and BHMT2: common gene sequence variation and functional characterization.

Betaine-homocysteine methyltransferase (BHMT) catalyzes the remethylation of homocysteine. BHMT2 encodes a protein 73% identical in amino acid sequence to BHMT, but the function of BHMT2 remains unclear. We set out to identify and functionally characterize common genetic variation in BHMT and BHMT2. Specifically, we sequenced exons, exon-intron splice junctions and the 5'-flanking regions (5'-FRs) of BHMT and BHMT2 using 240 DNA samples from four ethnic groups. Twenty-five single nucleotide polymorphisms (SNPs), including 4 nonsynonymous SNPs, and 39 SNPs, including 4 nonsynonymous, were observed in BHMT and BHMT2, respectively. BHMT wild type (WT) and variant allozymes were expressed in COS-1 cells. Variant allozymes showed no significant differences from WT in levels of enzyme activity or immunoreactive protein, but there were statistically significant differences in apparent K(m) values. Luciferase reporter gene constructs were created for the three most common BHMT 5'-FR haplotypes, and significant variation was observed in the ability of these constructs to drive transcription. Although BHMT2 mRNA has been observed in human liver and kidney, expression of the protein has not been reported. We were unable to express BHMT2 in mammalian cells, and the protein aggregated after bacterial expression. Furthermore, BHMT2 was rapidly degraded in a rabbit reticulocyte lysate, but it could be stabilized by cotransfection of COS-1 cells with BHMT and, after cotransfection, it coprecipitated with BHMT. These studies have defined common genetic variation in BHMT and BHMT2 and functionally characterized BHMT SNPs. They may also help to explain why BHMT2 has not previously been defined functionally.

Gemcitabine pharmacogenomics: deoxycytidine kinase and cytidylate kinase gene resequencing and functional genomics.

Gemcitabine and other cytidine antimetabolites require metabolic activation by phosphorylation. Deoxycytidine kinase (DCK) and cytidine monophosphate kinase (CMPK) catalyze these reactions. We have applied a genotype-to-phenotype strategy to study DCK and CMPK pharmacogenomics. Specifically, we resequenced DCK and CMPK using 240 DNA samples, 60 each from African-American, Caucasian-American, Han Chinese-American, and Mexican-American subjects. We observed 28 DCK polymorphisms and 28 polymorphisms in CMPK, 33 of which were novel. Expression in COS-1 cells showed that variant allozyme enzyme activities ranged from 32 to 105% of the wild type (WT) for DCK and from 78 to 112% of WT for CMPK--with no significant differences in apparent K(m) values for either enzyme except for a DCK Val24/Ser122 double variant allozyme. Relative levels of DCK and CMPK immunoreactive protein in the COS-1 cells paralleled relative levels of enzyme activity and were significantly correlated for DCK (R(p) = 0.89, P = 0.0004) but not for CMPK (R(p) = 0.82, P = 0.095). The results of an analysis of DCK and CMPK structural models were compatible with the observed functional consequences of sequence alterations in variant allozymes. We also confirmed that the CMPK protein expressed in COS-1 cells and in a rabbit reticulocyte lysate was 196 rather than 228 amino acids in length. In summary, we determined common sequence variations in DCK and CMPK and systematically evaluated their functional implications. These gene sequence differences may contribute to variations in the metabolic activation of gemcitabine and other cytidine antimetabolites.