Phenotyping Hepatic Immune-Related Adverse Events in the Setting of Immune Checkpoint Inhibitor Therapy.

PURPOSE: We present and validate a rule-based algorithm for the detection of moderate to severe liver-related immune-related adverse events (irAEs) in a real-world patient cohort. The algorithm can be applied to studies of irAEs in large data sets. METHODS: We developed a set of criteria to define hepatic irAEs. The criteria include: the temporality of elevated laboratory measurements in the first 2-14 weeks of immune checkpoint inhibitor (ICI) treatment, steroid intervention within 2 weeks of the onset of elevated laboratory measurements, and intervention with a duration of at least 2 weeks. These criteria are based on the kinetics of patients who experienced moderate to severe hepatotoxicity (Common Terminology Criteria for Adverse Events grades 2-4). We applied these criteria to a retrospective cohort of 682 patients diagnosed with hepatocellular carcinoma and treated with ICI. All patients were required to have baseline laboratory measurements before and after the initiation of ICI. RESULTS: A set of 63 equally sampled patients were reviewed by two blinded, clinical adjudicators. Disagreements were reviewed and consensus was taken to be the ground truth. Of these, 25 patients with irAEs were identified, 16 were determined to be hepatic irAEs, 36 patients were nonadverse events, and two patients were of indeterminant status. Reviewers agreed in 44 of 63 patients, including 19 patients with irAEs (0.70 concordance, Fleiss' kappa: 0.43). By comparison, the algorithm achieved a sensitivity and specificity of identifying hepatic irAEs of 0.63 and 0.81, respectively, with a test efficiency (percent correctly classified) of 0.78 and outcome-weighted F1 score of 0.74. CONCLUSION: The algorithm achieves greater concordance with the ground truth than either individual clinical adjudicator for the detection of irAEs.

Baseline correlates of frailty and its association with survival in United States veterans with acute myeloid leukemia.

Frailty is an important construct to measure in acute myeloid leukemia (AML). We used the Veterans Affairs Frailty Index (VA-FI) - calculated using readily available data within the VA's electronic health records - to measure frailty in U.S. veterans with AML. Of the 1166 newly diagnosed and treated veterans with AML between 2012 and 2022, 722 (62%) veterans with AML were classified as frail (VA-FI > 0.2). At a median follow-up of 252.5 days, moderate-severely frail veterans had significantly worse survival than mildly frail, and non-frail veterans (median survival 179 vs. 306 vs. 417 days, p < .001). Increasing VA-FI severity was associated with higher mortality. A model with VA-FI in addition to the European LeukemiaNet (ELN) risk classification and other covariates statistically outperformed a model containing the ELN risk and other covariates alone (p < .001). These findings support the VA-FI as a tool to expand frailty measurement in research and clinical practice for informing prognosis in veterans with AML.

Machine learning-based natural language processing to extract PD-L1 expression levels from clinical notes.

Introduction: PD-L1 expression is used to determine oncology patients' response to and eligibility for immunologic treatments; however, PD-L1 expression status often only exists in unstructured clinical notes, limiting ability to use it in population-level studies. Methods: We developed and evaluated a machine learning based natural language processing (NLP) tool to extract PD-L1 expression values from the nationwide Veterans Affairs electronic health record system. Results: The model demonstrated strong evaluation performance across multiple levels of label granularity. Mean precision of the overall PD-L1 positive label was 0.859 (sd, 0.039), recall 0.994 (sd, 0.013), and F1 0.921 (0.024). When a numeric PD-L1 value was identified, the mean absolute error of the value was 0.537 on a scale of 0 to 100. Conclusion: We presented an accurate NLP method for deriving PD-L1 status from clinical notes. By reducing the time and manual effort needed to review medical records, our work will enable future population-level studies in cancer immunotherapy.

Prognostic value of the veterans affairs frailty index in older patients with non-small cell lung cancer.

BACKGROUND: Older patients with non-small cell lung cancer (NSCLC) are a heterogeneous population with varying degrees of frailty. An electronic frailty index such as the Veterans Affairs Frailty Index (VA-FI) can potentially help identify vulnerable patients at high risk of poor outcomes. METHODS: NSCLC patients ≥65 years old and diagnosed in 2002-2017 were identified using the VA Central Cancer Registry. The VA-FI was calculated using administrative codes from VA electronic health records data linked with Medicare and Medicaid data. We assessed associations between the VA-FI and times to mortality, hospitalization, and emergency room (ER) visit following diagnosis by Kaplan-Meier analysis and multivariable stratified Cox models. We also evaluated the change in discrimination and calibration of reference prognostic models after adding VA-FI. RESULTS: We identified a cohort of 42,204 older NSCLC VA patients, in which 55.5% were classified as frail (VA-FI >0.2). After adjustment, there was a strong association between VA-FI and the risk of mortality (HR = 1.23 for an increase of four deficits or, equivalently, an increase of 0.129 on VA-FI, p < 0.001), hospitalization (HR = 1.16 for four deficits, p < 0.001), and ER visit (HR = 1.18 for four deficits, p < 0.001). Adding VA-FI to baseline prognostic models led to statistically significant improvements in time-dependent area under curves and did not have a strong impact on calibration. CONCLUSION: Older NSCLC patients with higher VA-FI have significantly elevated risks of mortality, hospitalizations, and ER visits following diagnosis. An electronic frailty index can serve as an accessible tool to identify patients with vulnerabilities to inform clinical care and research.

Association of COVID-19 Vaccination With SARS-CoV-2 Infection in Patients With Cancer: A US Nationwide Veterans Affairs Study.

Importance: Patients with cancer are at increased risk for severe COVID-19, but it is unknown whether SARS-CoV-2 vaccination is effective for them. Objective: To determine the association between SARS-CoV-2 vaccination and SARS-CoV-2 infections among a population of Veterans Affairs (VA) patients with cancer. Design, Setting, and Participants: Retrospective, multicenter, nationwide cohort study of SARS-CoV-2 vaccination and infection among patients in the VA health care system from December 15, 2020, to May 4, 2021. All adults with solid tumors or hematologic cancer who received systemic cancer-directed therapy from August 15, 2010, to May 4, 2021, and were alive and without a documented SARS-CoV-2 positive result as of December 15, 2020, were eligible for inclusion. Each day between December 15, 2020, and May 4, 2021, newly vaccinated patients were matched 1:1 with unvaccinated or not yet vaccinated controls based on age, race and ethnicity, VA facility, rurality of home address, cancer type, and treatment type/timing. Exposures: Receipt of a SARS-CoV-2 vaccine. Main Outcomes and Measures: The primary outcome was documented SARS-CoV-2 infection. A proxy for vaccine effectiveness was defined as 1 minus the risk ratio of SARS-CoV-2 infection for vaccinated individuals compared with unvaccinated controls. Results: A total of 184 485 patients met eligibility criteria, and 113 796 were vaccinated. Of these, 29 152 vaccinated patients (median [IQR] age, 74.1 [70.2-79.3] years; 95% were men; 71% were non-Hispanic White individuals) were matched 1:1 to unvaccinated or not yet vaccinated controls. As of a median 47 days of follow-up, 436 SARS-CoV-2 infections were detected in the matched cohort (161 infections in vaccinated patients vs 275 in unvaccinated patients). There were 17 COVID-19-related deaths in the vaccinated group vs 27 COVID-19-related deaths in the unvaccinated group. Overall vaccine effectiveness in the matched cohort was 58% (95% CI, 39% to 72%) starting 14 days after the second dose. Patients who received chemotherapy within 3 months prior to the first vaccination dose were estimated to have a vaccine effectiveness of 57% (95% CI, -23% to 90%) starting 14 days after the second dose vs 76% (95% CI, 50% to 91%) for those receiving endocrine therapy and 85% (95% CI, 29% to 100%) for those who had not received systemic therapy for at least 6 months prior. Conclusions and Relevance: In this cohort study, COVID-19 vaccination was associated with lower SARS-CoV-2 infection rates in patients with cancer. Some immunosuppressed subgroups may remain at early risk for COVID-19 despite vaccination, and consideration should be given to additional risk reduction strategies, such as serologic testing for vaccine response and a third vaccine dose to optimize outcomes.

Multi-scale Pan-cancer Integrative Analyses Identify the STAT3-VSIR Axis as a Key Immunosuppressive Mechanism in Head and Neck Cancer.

PURPOSE: VSIR is a novel immune checkpoint protein whose expression on tumor cells across cancers remains largely uncharacterized. Here we purposed to decode the pan-cancer biologic and clinical significance of VSIR overexpression in the tumor compartment. EXPERIMENTAL DESIGN: We performed multi-omics integrative analyses of 9,735 tumor samples to identify cancers with non-leukocytic expression of VSIR (VSIR High), followed by association with overall survival and immune cell infiltration levels. Orthogonal assessments of VSIR protein expression and lymphocytic infiltration were performed using quantitative immunofluorescence (QIF). RESULTS: Integrative modeling identified a subset of cancer types as being enriched for VSIR High tumors. VSIR High tumors were associated with significantly poorer overall survival in immunogenic ovarian serous adenocarcinoma (SA) and oral cavity squamous cell carcinoma (SCC). QIF assessments in an independent validation cohort confirmed overexpression of VSIR as being associated with poorer overall survival within immunogenic oral cavity SCC. VSIR overexpression was associated with lower CD4 helper T-cell infiltration in both ovarian SA and oral cavity SCC, but did not impact CD8 T-cell infiltration. VSIR overexpressing tumors in both cancer types exhibited significantly higher STAT3 signaling activity. Pharmacologic inhibition of STAT3 signaling resulted in dose-dependent reduction of VSIR expression in ovarian SA and oral cavity SCC cells. CONCLUSIONS: The STAT3-VSIR axis is a potentially significant immunomodulatory mechanism in oral cavity and ovarian cancers, whose activation is associated with poorer survival and an immune microenvironment marked by decreased CD4 helper T-cell activity. The role of VSIR as a tumor-intrinsic modulator of resistance to immunotherapy warrants further exploration.

Prevalence and Outcome of COVID-19 Infection in Cancer Patients: A National Veterans Affairs Study.

BACKGROUND: Emerging data suggest variability in susceptibility and outcome to coronavirus disease 2019 (COVID-19) infection. Identifying risk factors associated with infection and outcomes in cancer patients is necessary to develop healthcare recommendations. METHODS: We analyzed electronic health records of the US Veterans Affairs Healthcare System and assessed the prevalence of COVID-19 infection in cancer patients. We evaluated the proportion of cancer patients tested for COVID-19 who were positive, as well as outcome attributable to COVID-19, and stratified by clinical characteristics including demographics, comorbidities, cancer treatment, and cancer type. All statistical tests are 2-sided. RESULTS: Of 22 914 cancer patients tested for COVID-19, 1794 (7.8%) were positive. The prevalence of COVID-19 was similar across age. Higher prevalence was observed in African American (15.0%) compared with White (5.5%; P < .001) and in patients with hematologic malignancy compared with those with solid tumors (10.9% vs 7.8%; P < .001). Conversely, prevalence was lower in current smokers and patients who recently received cancer therapy (<6 months). The COVID-19-attributable mortality was 10.9%. Higher attributable mortality rates were observed in older patients, those with higher Charlson comorbidity score, and in certain cancer types. Recent (<6 months) or past treatment did not influence attributable mortality. Importantly, African American patients had 3.5-fold higher COVID-19-attributable hospitalization; however, they had similar attributable mortality as White patients. CONCLUSION: Preexistence of cancer affects both susceptibility to COVID-19 infection and eventual outcome. The overall COVID-19-attributable mortality in cancer patients is affected by age, comorbidity, and specific cancer types; however, race or recent treatment including immunotherapy do not impact outcome.

Dual HDAC and PI3K Inhibitor CUDC-907 Downregulates MYC and Suppresses Growth of MYC-dependent Cancers.

Upregulation of MYC is a common driver event in human cancers, and some tumors depend on MYC to maintain transcriptional programs that promote cell growth and proliferation. Preclinical studies have suggested that individually targeting upstream regulators of MYC, such as histone deacetylases (HDAC) and phosphoinositide 3-kinases (PI3K), can reduce MYC protein levels and suppress the growth of MYC-driven cancers. Synergy between HDAC and PI3K inhibition in inducing cancer cell death has also been reported, but the involvement of MYC regulation is unclear. In this study, we demonstrated that HDAC and PI3K inhibition synergistically downregulates MYC protein levels and induces apoptosis in "double-hit" (DH) diffuse large B-cell lymphoma (DLBCL) cells. Furthermore, CUDC-907, a small-molecule dual-acting inhibitor of both class I and II HDACs and class I PI3Ks, effectively suppresses the growth and survival of MYC-altered or MYC-dependent cancer cells, such as DH DLBCL and BRD-NUT fusion-positive NUT midline carcinoma (NMC) cells, and MYC protein downregulation is an early event induced by CUDC-907 treatment. Consistently, the antitumor activity of CUDC-907 against multiple MYC-driven cancer types was also demonstrated in animal models, including DLBCL and NMC xenograft models, Myc transgenic tumor syngeneic models, and MYC-amplified solid tumor patient-derived xenograft (PDX) models. Our findings suggest that dual function HDAC and PI3K inhibitor CUDC-907 is an effective agent targeting MYC and thus may be developed as potential therapy for MYC-dependent cancers. Mol Cancer Ther; 16(2); 285-99. ©2016 AACR.

Estrogen Drives Cellular Transformation and Mutagenesis in Cells Expressing the Breast Cancer-Associated R438W DNA Polymerase Lambda Protein.

Repair of DNA damage is critical for maintaining the genomic integrity of cells. DNA polymerase lambda (POLL/Pol λ) is suggested to function in base excision repair (BER) and nonhomologous end-joining (NHEJ), and is likely to play a role in damage tolerance at the replication fork. Here, using next-generation sequencing, it was discovered that the POLL rs3730477 single-nucleotide polymorphism (SNP) encoding R438W Pol λ was significantly enriched in the germlines of breast cancer patients. Expression of R438W Pol λ in human breast epithelial cells induces cellular transformation and chromosomal aberrations. The role of estrogen was assessed as it is commonly used in hormone replacement therapies and is a known breast cancer risk factor. Interestingly, the combination of estrogen treatment and the expression of the R438W Pol λ SNP drastically accelerated the rate of transformation. Estrogen exposure produces 8-oxoguanine lesions that persist in cells expressing R438W Pol λ compared with wild-type (WT) Pol λ-expressing cells. Unlike WT Pol λ, which performs error-free bypass of 8-oxoguanine lesions, expression of R438W Pol λ leads to an increase in mutagenesis and replicative stress in cells treated with estrogen. Together, these data suggest that individuals who carry the rs3730477 POLL germline variant have an increased risk of estrogen-associated breast cancer. IMPLICATIONS: The Pol λ R438W mutation can serve as a biomarker to predict cancer risk and implicates that treatment with estrogen in individuals with this mutation may further increase their risk of breast cancer. Mol Cancer Res; 14(11); 1068-77. ©2016 AACR.

Phase 1b study of the mammalian target of rapamycin inhibitor sirolimus in combination with nanoparticle albumin-bound paclitaxel in patients with advanced solid tumors.

BACKGROUND: The optimal weekly oral dose of sirolimus and intravenous nanoparticle albumin-bound paclitaxel (nab-paclitaxel) were evaluated. METHODS: A phase 1b study was performed to evaluate escalating doses of oral sirolimus (5-60 mg) on days 2, 9, and 16 with intravenous nab-paclitaxel (100 mg/m(2) ) on days 1, 8, and 15 in a 28-day cycle. A run-in treatment of nab-paclitaxel (day -14) and sirolimus (day -7) was administered for pharmacokinetic and pharmacodynamic assessments. Clinical trial endpoints included dose-limiting toxicities (DLTs), maximum tolerated doses, and response rates. Pharmacodynamics included immunohistochemistry for phosphatase and tensin homolog, mammalian target of rapamycin (mTOR), AKT, phosphorylated AKT, S6K1, and phosphorylated S6K1; exploratory gene expression analysis; and [(18) F]fludeoxyglucose (FDG) positron emission tomography. RESULTS: Twenty-three patients with advanced solid tumors were treated. Fifteen patients had prior taxane therapy. Twenty-two patients were evaluable for responses. One patient had a complete response, and 5 patients had a partial response (3 confirmed). DLTs were seen in 1 patient each at 10 (grade 3 dyspnea/hypoxia) and 40 mg (grade 4 leukopenia/neutropenia) and in 2 patients at 60 mg (grade 3 fatigue and grade 4 pericardial effusion). Patients with higher expression of posttreatment AKT and a greater decline in FDG activity were more likely to have a treatment response or stable disease. CONCLUSIONS: Sirolimus showed an acceptable safety profile at a weekly dose of 40 mg with weekly intravenous nab-paclitaxel at 100 mg/m(2) on days 1, 8, and 15 every 28 days. The posttreatment AKT score and changes in FDG activity may have roles as early predictors of responses to mTOR inhibitors.

Molecular phenotypes in triple negative breast cancer from African American patients suggest targets for therapy.

Triple negative breast cancer (TNBC) is characterized by high proliferation, poor differentiation and a poor prognosis due to high rates of recurrence. Despite lower overall incidence African American (AA) patients suffer from higher breast cancer mortality in part due to the higher proportion of TNBC cases among AA patients compared to European Americans (EA). It was recently shown that the clinical heterogeneity of TNBC is reflected by distinct transcriptional programs with distinct drug response profiles in preclinical models. In this study, gene expression profiling and immunohistochemistry were used to elucidate potential differences between TNBC tumors of EA and AA patients on a molecular level. In a retrospective cohort of 136 TNBC patients, a major transcriptional signature of proliferation was found to be significantly upregulated in samples of AA ethnicity. Furthermore, transcriptional profiles of AA tumors showed differential activation of insulin-like growth factor 1 (IGF1) and a signature of BRCA1 deficiency in this cohort. Using signatures derived from the meta-analysis of TNBC gene expression carried out by Lehmann et al., tumors from AA patients were more likely of basal-like subtypes whereas transcriptional features of many EA samples corresponded to mesenchymal-like or luminal androgen receptor driven subtypes. These results were validated in The Cancer Genome Atlas mRNA and protein expression data, again showing enrichment of a basal-like phenotype in AA tumors and mesenchymal subtypes in EA tumors. In addition, increased expression of VEGF-activated genes together with elevated microvessel area determined by the AQUA method suggest that AA patients exhibit higher tumor vascularization. This study confirms the existence of distinct transcriptional programs in triple negative breast cancer in two separate cohorts and that these programs differ by racial group. Differences in TNBC subtypes and levels of tumor angiogenesis in AA versus EA patients suggest that targeted therapy choices should be considered in the context of race.

PAM50 proliferation score as a predictor of weekly paclitaxel benefit in breast cancer.

To identify a group of patients who might benefit from the addition of weekly paclitaxel to conventional anthracycline-containing chemotherapy as adjuvant therapy of node-positive operable breast cancer. The predictive value of PAM50 subtypes and the 11-gene proliferation score contained within the PAM50 assay were evaluated in 820 patients from the GEICAM/9906 randomized phase III trial comparing adjuvant FEC to FEC followed by weekly paclitaxel (FEC-P). Multivariable Cox regression analyses of the secondary endpoint of overall survival (OS) were performed to determine the significance of the interaction between treatment and the (1) PAM50 subtypes, (2) PAM50 proliferation score, and (3) clinical and pathological variables. Similar OS analyses were performed in 222 patients treated with weekly paclitaxel versus paclitaxel every 3 weeks in the CALGB/9342 and 9840 metastatic clinical trials. In GEICAM/9906, with a median follow up of 8.7 years, OS of the FEC-P arm was significantly superior compared to the FEC arm (unadjusted HR = 0.693, p = 0.013). A benefit from paclitaxel was only observed in the group of patients with a low PAM50 proliferation score (unadjusted HR = 0.23, p < 0.001; and interaction test, p = 0.006). No significant interactions between treatment and the PAM50 subtypes or the various clinical-pathological variables, including Ki-67 and histologic grade, were identified. Finally, similar OS results were obtained in the CALGB data set, although the interaction test did not reach statistical significance (p = 0.109). The PAM50 proliferation score identifies a subset of patients with a low proliferation status that may derive a larger benefit from weekly paclitaxel.

Hypoxia-induced protein CAIX is associated with somatic loss of BRCA1 protein and pathway activity in triple negative breast cancer.

The purpose of this study is to explore the relationship between tumor hypoxia assessed by CA IX protein expression and loss of BRCA1 function in triple negative breast cancer (TNBC). Protein expression of CA IX and BRCA1 was evaluated by AQUA™ technology on two breast cancer cohorts: an unselected cohort of 637 breast cancer patients and a TNBC cohort of 120 patients. Transcriptional profiling was performed on FFPE samples from the TNBC cohort to evaluate a gene expression signature associated with BRCA1 mutation (van't Veer et al., Nature 415(6871):530-536, 2002). CA IX is expressed in 7 % of the unselected breast cancer cohort and in 25 % of the TNBCs and is significantly associated with the triple negative phenotype. CA IX protein expression and BRCA1 protein expression are inversely correlated in both cohorts. Patients expressing high levels of CA IX show significantly worse overall survival (p = 0.02). Importantly, high CA IX protein expression occurs in patients who show the BRCA1 mutant signature and low levels of BRCA1 protein. These data suggest that elevated CA IX protein in TNBC is associated with a BRCA1 mutant signature and loss of BRCA1 function. CA IX may be a useful biomarker to identify triple negative patients with defective homologous recombination, who might benefit from PARP inhibitor therapy.

Loss of the retinoblastoma binding protein 2 (RBP2) histone demethylase suppresses tumorigenesis in mice lacking Rb1 or Men1.

Aberrations in epigenetic processes, such as histone methylation, can cause cancer. Retinoblastoma binding protein 2 (RBP2; also called JARID1A or KDM5A) can demethylate tri- and dimethylated lysine 4 in histone H3, which are epigenetic marks for transcriptionally active chromatin, whereas the multiple endocrine neoplasia type 1 (MEN1) tumor suppressor promotes H3K4 methylation. Previous studies suggested that inhibition of RBP2 contributed to tumor suppression by the retinoblastoma protein (pRB). Here, we show that genetic ablation of Rbp2 decreases tumor formation and prolongs survival in Rb1(+/-) mice and Men1-defective mice. These studies link RBP2 histone demethylase activity to tumorigenesis and nominate RBP2 as a potential target for cancer therapy.

Modeling the clonal heterogeneity of stem cells.

Recent experimental studies suggest that tissue stem cell pools are composed of functionally diverse clones. Metapopulation models in ecology concentrate on collections of populations and their role in stabilizing coexistence and maintaining selected genetic or epigenetic variation. Such models are characterized by expansion and extinction of spatially distributed populations. We develop a mathematical framework derived from the multispecies metapopulation model of Tilman et al (1994) to study the dynamics of heterogeneous stem cell metapopulations. In addition to normal stem cells, the model can be applied to cancer cell populations and their response to treatment. In our model disturbances may lead to expansion or contraction of cells with distinct properties, reflecting proliferation, apoptosis, and clonal competition. We first present closed-form expressions for the basic model which defines clonal dynamics in the presence of exogenous global disturbances. We then extend the model to include disturbances which are periodic and which may affect clones differently. Within the model framework, we propose a method to devise an optimal strategy of treatments to regulate expansion, contraction, or mutual maintenance of cells with specific properties.

Downstream targets of HOXB4 in a cell line model of primitive hematopoietic progenitor cells.

Enforced expression of the homeobox transcription factor HOXB4 has been shown to enhance hematopoietic stem cell self-renewal and expansion ex vivo and in vivo. To investigate the downstream targets of HOXB4 in hematopoietic progenitor cells, HOXB4 was constitutively overexpressed in the primitive hematopoietic progenitor cell line EML. Two genome-wide analytical techniques were used: RNA expression profiling using microarrays and chromatin immunoprecipitation (ChIP)-chip. RNA expression profiling revealed that 465 gene transcripts were differentially expressed in KLS (c-Kit(+), Lin(-), Sca-1(+))-EML cells that overexpressed HOXB4 (KLS-EML-HOXB4) compared with control KLS-EML cells that were transduced with vector alone. In particular, erythroid-specific gene transcripts were observed to be highly down-regulated in KLS-EML-HOXB4 cells. ChIP-chip analysis revealed that the promoter region for 1910 genes, such as CD34, Sox4, and B220, were occupied by HOXB4 in KLS-EML-HOXB4 cells. Side-by-side comparison of the ChIP-chip and RNA expression profiling datasets provided correlative information and identified Gp49a and Laptm4b as candidate "stemness-related" genes. Both genes were highly ranked in both dataset lists and have been previously shown to be preferentially expressed in hematopoietic stem cells and down-regulated in mature hematopoietic cells, thus making them attractive candidates for future functional studies in hematopoietic cells.

Distinct human papillomavirus type 16 methylomes in cervical cells at different stages of premalignancy.

Human papillomavirus (HPV) gene expression is dramatically altered during cervical carcinogenesis. Because dysregulated genes frequently show abnormal patterns of DNA methylation, we hypothesized that comprehensive mapping of the HPV methylomes in cervical samples at different stages of progression would reveal patterns of clinical significance. To test this hypothesis, thirteen HPV16-positive samples were obtained from women undergoing routine cervical cancer screening. Complete methylation data were obtained for 98.7% of the HPV16 CpGs in all samples by bisulfite-sequencing. Most HPV16 CpGs were unmethylated or methylated in only one sample. The other CpGs were methylated at levels ranging from 11% to 100% of the HPV16 copies per sample. The results showed three major patterns and two variants of one pattern. The patterns showed minimal or no methylation (A), low level methylation in the E1 and E6 genes (B), and high level methylation at many CpGs in the E5/L2/L1 region (C). Generally, pattern A was associated with negative cytology, pattern B with low-grade lesions, and pattern C with high-grade lesions. The severity of the cervical lesions was then ranked by the HPV16 DNA methylation patterns and, independently, by the pathologic diagnoses. Statistical analysis of the two rating methods showed highly significant agreement. In conclusion, analysis of the HPV16 DNA methylomes in clinical samples of cervical cells led to the identification of distinct methylation patterns which, after validation in larger studies, could have potential utility as biomarkers of neoplastic cervical progression.

Quantitative assessment of tissue biomarkers and construction of a model to predict outcome in breast cancer using multiple imputation.

Missing data pose one of the greatest challenges in the rigorous evaluation of biomarkers. The limited availability of specimens with complete clinical annotation and quality biomaterial often leads to underpowered studies. Tissue microarray studies, for example, may be further handicapped by the loss of data points because of unevaluable staining, core loss, or the lack of tumor in the histospot. This paper presents a novel approach to these common problems in the context of a tissue protein biomarker analysis in a cohort of patients with breast cancer. Our analysis develops techniques based on multiple imputation to address the missing value problem. We first select markers using a training cohort, identifying a small subset of protein expression levels that are most useful in predicting patient survival. The best model is obtained by including both protein markers (including COX6C, GATA3, NAT1, and ESR1) and lymph node status. The use of either lymph node status or the four protein expression levels provides similar improvements in goodness-of-fit, with both significantly better than a baseline clinical model. Using the same multiple imputation strategy, we then validate the results out-of-sample on a larger independent cohort. Our approach of integrating multiple imputation with each stage of the analysis serves as an example that may be replicated or adapted in future studies with missing values.

Global variation in CYP2C8-CYP2C9 functional haplotypes.

We have studied the global frequency distributions of 10 single nucleotide polymorphisms (SNPs) across 132 kb of CYP2C8 and CYP2C9 in approximately 2500 individuals representing 45 populations. Five of the SNPs were in noncoding sequences; the other five involved the more common missense variants (four in CYP2C8, one in CYP2C9) that change amino acids in the gene products. One haplotype containing two CYP2C8 coding variants and one CYP2C9 coding variant reaches an average frequency of 10% in Europe; a set of haplotypes with a different CYP2C8 coding variant reaches 17% in Africa. In both cases these haplotypes are found in other regions of the world at <1%. This considerable geographic variation in haplotype frequencies impacts the interpretation of CYP2C8/CYP2C9 association studies, and has pharmacogenomic implications for drug interactions.