NDN and CD1A are novel prognostic methylation markers in patients with head and neck squamous carcinomas

Genome-wide DNA methylation profiling of HPV-negative leukoplakia and gingivobuccal complex cancers

BACKGROUND

Gingivobuccal complex oral squamous cell carcinoma (GBC-OSCC) is an aggressive malignancy with high mortality often preceded by premalignant lesions, including leukoplakia. Previous studies have reported genomic drivers in OSCC, but much remains to be elucidated about DNA methylation patterns across different stages of oral carcinogenesis.

RESULTS

There is a serious lack of biomarkers and clinical application of biomarkers for early detection and prognosis of gingivobuccal complex cancers. Hence, in search of novel biomarkers, we measured genome-wide DNA methylation in 22 normal oral tissues, 22 leukoplakia, and 74 GBC-OSCC tissue samples. Both leukoplakia and GBC-OSCC had distinct methylation profiles as compared to normal oral tissue samples. Aberrant DNA methylation increases during the different stages of oral carcinogenesis, from premalignant lesions to carcinoma. We identified 846 and 5111 differentially methylated promoters in leukoplakia and GBC-OSCC, respectively, with a sizable fraction shared between the two sets. Further, we identified potential biomarkers from integrative analysis in gingivobuccal complex cancers and validated them in an independent cohort. Integration of genome, epigenome, and transcriptome data revealed candidate genes with gene expression synergistically regulated by copy number and DNA methylation changes. Regularised Cox regression identified 32 genes associated with patient survival. In an independent set of samples, we validated eight genes (FAT1, GLDC, HOXB13, CST7, CYB5A, MLLT11, GHR, LY75) from the integrative analysis and 30 genes from previously published reports. Bisulfite pyrosequencing validated GLDC (P = 0.036), HOXB13 (P < 0.0001) promoter hypermethylation, and FAT1 (P < 0.0001) hypomethylation in GBC-OSCC compared to normal controls.

CONCLUSIONS

Our findings identified methylation signatures associated with leukoplakia and gingivobuccal complex cancers. The integrative analysis in GBC-OSCC identified putative biomarkers that enhance existing knowledge of oral carcinogenesis and may potentially help in risk stratification and prognosis of GBC-OSCC.

The Epigenesis of Salivary Glands Carcinoma: From Field Cancerization to Carcinogenesis

Salivary gland carcinomas (SGCs) are a diverse collection of malignant tumors with marked differences in biological activity, clinical presentation and microscopic appearance. Although the etiology is varied, secondary radiation, oncogenic viruses as well as chromosomal rearrangements have all been linked to the formation of SGCs. Epigenetic modifications may also contribute to the genesis and progression of SGCs. Epigenetic modifications are any heritable changes in gene expression that are not caused by changes in DNA sequence. It is now widely accepted that epigenetics plays an important role in SGCs development. A basic epigenetic process that has been linked to a variety of pathological as well as physiological conditions including cancer formation, is DNA methylation. Transcriptional repression is caused by CpG islands hypermethylation at gene promoters, whereas hypomethylation causes overexpression of a gene. Epigenetic changes in SGCs have been identified, and they have been linked to the genesis, progression as well as prognosis of these neoplasms. Thus, we conduct a thorough evaluation of the currently known evidence on the involvement of epigenetic processes in SGCs.

Prognostic and therapeutic prediction by screening signature combinations from transcriptome-methylome interactions in oral squamous cell carcinoma

DNA methylation pattern in oral squamous cell carcinoma (OSCC) remains poorly described. This study aimed to perform a genome-wide integrated analysis of the transcriptome and methylome and assess the efficacy of their prognostic signature model in patients with OSCC. We analyzed transcriptome and methylome data from 391 OSCC samples and 41 adjacent normal samples. A total of 8074 differentially expressed genes (DEGs) and 10,084 differentially expressed CpGs (DMCpGs) were identified. Then 241 DEGs with DMCpGs were identified. According to the prognostic analysis, the prognostic signature of methylation-related differentially expressed genes (mrDEGPS) was established. mrDEGPS consisted of seven prognostic methylation-related genes, including ESRRG, CCNA1, SLC20A1, COL6A6, FCGBP, CDKN2A, and ZNF43. mrDEGPS was a significant stratification factor of survival (P < 0.00001) irrespective of the clinical stage. The immune effector components, including B cells, CD4⁺ T cells, and CD8⁺ T cells, were decreased in the tumor environment of patients with high mrDEGPS. Immune checkpoint expressions, including CTLA-4, PD-1, LAG3, LGALS9, HAVCR2, and TIGHT, were comprehensively elevated (P < 0.001). The estimated half-maximal inhibitory concentration difference between low- and high-risk patients was inconsistent among chemotherapeutic drugs. In conclusion, the transcriptome–methylome interaction pattern in OSCC is complex. mrDEGPS can predict patient survival and responses to immunotherapy and chemotherapy and facilitate clinical decision-making in patients with OSCC.

Epigenetic Regulations of Perineural Invasion in Head and Neck Squamous Cell Carcinoma

Carcinomas of the oral cavity and oropharynx belong among the ten most common malignancies in the human population. The prognosis of head and neck squamous cell carcinoma (HNSCC) is determined by the degree of invasiveness of the primary tumor and by the extent of metastatic spread into regional and distant lymph nodes. Moreover, the level of the perineural invasion itself associates with tumor localization, invasion's extent, and the presence of nodal metastases. Here, we summarize the current knowledge about different aspects of epigenetic changes, which can be associated with HNSCC while focusing on perineural invasion (PNI). We review epigenetic modifications of the genes involved in the PNI process in HNSCC from the omics perspective and specific epigenetic modifications in OSCC or other neurotropic cancers associated with perineural invasion. Moreover, we summarize DNA methylation status of tumor-suppressor genes, methylation and demethylation enzymes and histone post-translational modifications associated with PNI. The influence of other epigenetic factors on the HNSCC incidence and perineural invasion such as tobacco, alcohol and oral microbiome is overviewed and HPV infection is discussed as an epigenetic factor associated with OSCC and related perineural invasion. Understanding epigenetic regulations of axon growth that lead to tumorous spread or uncovering the molecular control of axon interaction with cancer tissue can help to discover new therapeutic targets for these tumors.

Necdin, one of the important pathway proteins in the regulation of osteosarcoma progression by microRNA-200c

MicroRNA-200c (miR-200c) generally acts as a tumor suppressor in multiple cancer types and a promising therapeutic target in tumorigenesis. However, only a few studies have explained the role of miR-200c in the development of osteosarcoma (OS). In this study, we investigated the role of miR-200c in OS progression and identified the regulatory pathway protein NDN involved in inhibiting the occurrence and development of OS. Firstly, we found that miR-200c is downregulated in OS cells and tissues. As well, in vitro and in vivo experiments showed that upregulating miR-200c inhibits the proliferation, invasion, metastasis of Saos-2 cells, promotes the apoptosis of Saos-2 cells and suppresses tumor growth in mice, indicating miR-200c plays a major role in regulating the OS progression. Furthermore, bioinformatics analysis showed that an anti-tumor protein, necdin (NDN), might be a potential target by miR-200c. To verify this hypothesis, we measured the expression level of NDN in OS cells and tissues and found NDN is downregulated, suggesting NDN is functional in OS progression. Moreover, we found that the expression levels of NDN and miR-200c in in vivo and in vitro experiments were positively correlated. However, the results of dual-luciferase reporter gene experiment showed miR-200c does not directly act on the 3ʹ untranslated region (UTR) of NDN gene, indicating that NDN might be an important pathway protein which regulates OS progression in the presence of miR-200c. Therefore, miR-200c/NDN could be potential targets for developing effective treatment against OS.

Novel immune-related signature for risk stratification and prognosis in prostatic adenocarcinoma

A substantial proportion of prostatic adenocarcinoma (PRAD) patients experience biochemical failure (BCF) after radical prostatectomy (RP). The immune microenvironment plays a vital role in carcinogenesis and the development of PRAD. This study aimed to identify a novel immune-related gene (IRG)-based signature for risk stratification and prognosis of BCF in PRAD. Weighted gene coexpression network analysis was carried out to identify a BCF-related module in a discovery cohort of patients who underwent RP at the Massachusetts General Hospital. The median follow-up time was 70.32 months. Random forest and multivariate stepwise Cox regression analyses were used to identify an IRG-based signature from the specific module. Risk plot analyses, Kaplan-Meier curves, receiver operating characteristic curves, univariate and multivariate Cox regression analyses, stratified analysis, and Harrell's concordance index were used to assess the prognostic value and predictive accuracy of the IRG-based signature in the internal discovery cohort; The Cancer Genome Atlas database was used as a validation cohort. Tumor immune estimation resource database analysis and CIBERSORT algorithm were used to assess the immunophenotype of PRAD. A novel IRG-based signature was identified from the specific module. Five IRGs (BUB1B, NDN, NID1, COL4A6, and FLRT2) were verified as components of the risk signature. The IRG-based signature showed good prognostic value and predictive accuracy in both the discovery and validation cohorts. Infiltrations of various immune cells were significantly different between low-risk and high-risk groups in PRAD. We identified a novel IRG-based signature that could function as an index for assessing tumor immune status and risk stratification in PRAD.

Portrait of DNA methylated genes predictive of poor prognosis in head and neck cancer and the implication for targeted therapy

In addition to chronic infection with human papilloma virus (HPV) and exposure to environmental carcinogens, genetic and epigenetic factors act as major risk factors for head and neck cancer (HNC) development and progression. Here, we conducted a systematic review in order to assess whether DNA hypermethylated genes are predictive of high risk of developing HNC and/or impact on survival and outcomes in non-HPV/non-tobacco/non-alcohol associated HNC. We identified 85 studies covering 32,187 subjects where the relationship between DNA methylation, risk factors and survival outcomes were addressed. Changes in DNA hypermethylation were identified for 120 genes. Interactome analysis revealed enrichment in complex regulatory pathways that coordinate cell cycle progression (CCNA1, SFN, ATM, GADD45A, CDK2NA, TP53, RB1 and RASSF1). However, not all these genes showed significant statistical association with alcohol consumption, tobacco and/or HPV infection in the multivariate analysis. Genes with the most robust HNC risk association included TIMP3, DCC, DAPK, CDH1, CCNA1, MGMT, P16, MINT31, CD44, RARβ. From these candidates, we further validated CD44 at translational level in an independent cohort of 100 patients with tongue cancer followed-up beyond 10 years. CD44 expression was associated with high-risk of tumor recurrence and metastasis (P = 0.01) in HPV-cases. In summary, genes regulated by methylation play a modulatory function in HNC susceptibility and it represent a critical therapeutic target to manage patients with advanced disease.

DNA methylation biomarkers in peripheral blood of patients with head and neck squamous cell carcinomas. A systematic review

BACKGROUND

Head and neck squamous cell carcinomas (HNSCC) are often diagnosed in advanced stages. In search of new diagnostic tools, focus has shifted towards the biological properties of the HNSCC, and the number of different biomarkers under investigation is rapidly growing.

OBJECTIVES

The objective was to review the current literature regarding aberrantly methylated DNA found in peripheral blood plasma or serum in patients with HNSCC and to evaluate the diagnostic accuracy of these changes.

METHODS

The inclusion criteria were clinical studies involving patients with verified HNSCC that reported findings of aberrantly methylated DNA in peripheral blood serum or plasma. We systematically searched PubMed, OVID Embase and Cochrane Library. In addition to the search, we performed forward and backward chaining in references and Web of Science. The protocol was registered in PROSPERO: CRD42019135406. Two authors independently extracted data. The quality and the risk of bias of the included studies were assessed by the QUADAS-2 tool.

RESULTS

A total of 1,743 studies were found eligible for screening, while ultimately seven studies were included. All studies were found to have methodological weaknesses, mainly concerning patient selection bias. The best individual marker of HNSCC was Septin 9 in plasma with a sensitivity of 57% and a specificity of 95%.

CONCLUSIONS

None of the aberrantly methylated genes found in the retrieved studies are applicable as single diagnostic markers for HNSCC and the best gene-panels still lack diagnostic accuracy. Future studies may benefit from newer sequencing techniques but validation studies with well-designed cohorts are also needed in the process of developing epigenetic based diagnostic tests for HNSCC.

Epigenetic Modifications in Head and Neck Cancer

Head and neck squamous cell carcinoma (HNSCC) is the sixth most common human malignancy in the world, with high mortality and poor prognosis for patients. Among the risk factors are tobacco and alcohol intake, human papilloma virus, and also genetic and epigenetic modifications. Many studies show that epigenetic events play an important role in HNSCC development and progression, including DNA methylation, chromatin remodeling, histone posttranslational covalent modifications, and effects of non-coding RNA. Epigenetic modifications may influence silencing of tumor suppressor genes by promoter hypermethylation, regulate transcription by microRNAs and changes in chromatin structure, or induce genome instability through hypomethylation. Moreover, getting to better understand aberrant patterns of methylation may provide biomarkers for early detection and diagnosis, while knowledge about target genes of microRNAs may improve the therapy of HNSCC and extend overall survival. The aim of this review is to present recent studies which demonstrate the role of epigenetic regulation in the development of HNSCC.

Soy Isoflavone Supplementation Increases Long Interspersed Nucleotide Element-1 (LINE-1) Methylation in Head and Neck Squamous Cell Carcinoma

AIM

Soy isoflavones have been suggested as epigenetic modulating agents with effects that could be important in carcinogenesis. Hypomethylation of LINE-1 has been associated with head and neck squamous cell carcinoma (HNSCC) development from oral premalignant lesions and with poor prognosis. To determine if neoadjuvant soy isoflavone supplementation could modulate LINE-1 methylation in HNSCC, we undertook a clinical trial.

METHODS

Thirty-nine patients received 2-3 weeks of soy isoflavone supplements (300 mg/day) orally prior to surgery. Methylation of LINE-1, and 6 other genes was measured by pyrosequencing in biopsy, resection, and whole blood (WB) specimens. Changes in methylation were tested using paired t tests and ANOVA. Median follow up was 45 months.

RESULTS

LINE-1 methylation increased significantly after soy isoflavone (P < 0.005). Amount of change correlated positively with days of isoflavone taken (P = 0.04). Similar changes were not seen in corresponding WB samples. No significant changes in tumor or blood methylation levels were seen in the other candidate genes.

CONCLUSION

This is the first demonstration of in vivo increases in tissue-specific global methylation associated with soy isoflavone intake in patients with HNSCC. Prior associations of LINE-1 hypomethylation with genetic instability, carcinogenesis, and prognosis suggest that soy isoflavones maybe potential chemopreventive agents in HNSCC.

IRX-2 natural cytokine biologic for immunotherapy in patients with head and neck cancers

Head and neck squamous cell carcinoma (HNSCC) is an immunosuppressive malignancy characterized by tumor-driven immune-system abnormalities that contribute to disease progression. For patients with surgically resectable HNSCC, treatment is often curative surgery followed by irradiation or chemoradiation in high-risk settings to reduce the risk of recurrence. Poor survival and considerable morbidity of current treatments suggest the need for new therapeutic modalities that can improve outcomes. Defects in antitumor immunity of HNSCC patients include suppressed dendritic cell (DC) maturation, deficient antigen-presenting cell function, compromised natural killer (NK)-cell cytotoxicity, increased apoptosis of activated T lymphocytes, and impaired immune-cell migration to tumor sites. Strategies for relieving immunosuppression and restoring antitumor immune functions could benefit HNSCC patients. IRX-2 is a primary cell-derived biologic consisting of physiologic levels of T-helper type 1 cytokines produced by stimulating peripheral blood mononuclear cells of normal donors with phytohemagglutinin. The primary active components in IRX-2 are IL2, IL1β, IFNγ, and TNFα. In vitro, IRX-2 acts on multiple immune-system cell types, including DCs, T cells, and NK cells, to overcome tumor-mediated immunosuppression. In clinical settings, IRX-2 is administered as part of a 21-day neoadjuvant regimen, which includes additional pharmacologic agents (low-dose cyclophosphamide, indomethacin, and zinc) to promote anticancer immunoresponses. In a Phase IIA trial in 27 patients with surgically resectable, previously untreated HNSCC, neoadjuvant IRX-2 increased infiltration of T cells, B cells, and DCs into tumors and was associated with radiological reductions in tumor size. Event-free survival was 64% at 2 years, and overall 5-year survival was 65%. Follow-up and data analysis are under way in the multicenter, randomized, Phase IIB INSPIRE trial evaluating the IRX-2 regimen as a stand-alone therapy for activating the immune system to recognize and attack tumors.

Individualized outcome prognostication for patients with laryngeal cancer

BACKGROUND

Accurate prognostication is essential to the optimal management of laryngeal cancer. Predictive models have been developed to calculate the risk of oncologic outcomes, but extensive external validation of accuracy and reliability is necessary before implementing them into clinical practice.

METHOD

Four published prognostic calculators that predict 5-year overall survival for patients with laryngeal cancer were evaluated using patient information from a prospective epidemiology study cohort (n = 246; median follow-up, 60 months) with previously untreated, stage I through IVb laryngeal squamous cell carcinoma.

RESULTS

Different calculators yielded substantially different predictions for individual patients. The observed 5-year overall survival was significantly higher than the averaged predicted 5-year overall survival of the 4 calculators (71.9%; 95% confidence interval [CI], 65%-78%] vs 47.7%). Statistical analyses demonstrated the calculators' limited capacity to discriminate outcomes for risk-stratified patients. The area under the receiver operating characteristic curve ranged from 0.68 to 0.72. C-index values were similar for each of the 4 models (range, 0.66-0.68). There was a lower than expected hazard of death for patients who received induction (bioselective) chemotherapy (hazard ratio, 0.46; 95% CI, 0.24-0.88; P = .024) or primary surgical intervention (hazard ratio, 0.43; 95 % CI, 0.21-0.90; P = .024) compared with those who received concurrent chemoradiation.

CONCLUSIONS

Suboptimal reliability and accuracy limit the integration of existing individualized prediction tools into routine clinical decision making. The calculators predicted significantly worse than observed survival among patients who received induction chemotherapy and primary surgery, suggesting a need for updated consideration of modern treatment modalities. Further development of individualized prognostic calculators may improve risk prediction, treatment planning, and counseling for patients with laryngeal cancer. Cancer 2018;124:706-16. © 2017 American Cancer Society.

Whole-Genome DNA Methylation Analyses Revealed Epigenetic Instability in Tumorigenic Human iPS Cell-Derived Neural Stem/Progenitor Cells

Although human induced pluripotent stem cell (hiPSC) derivatives are considered promising cellular resources for regenerative medicine, their tumorigenicity potentially limits their clinical application in hiPSC technologies. We previously demonstrated that oncogenic hiPSC-derived neural stem/progenitor cells (hiPSC-NS/PCs) produced tumor-like tissues that were distinct from teratomas. To gain insight into the mechanisms underlying the regulation of tumorigenicity in hiPSC-NS/PCs, we performed an integrated analysis using the Infinium HumanMethylation450 BeadChip array and the HumanHT-12 v4.0 Expression BeadChip array to compare the comprehensive DNA methylation and gene expression profiles of tumorigenic hiPSC-NS/PCs (253G1-NS/PCs) and non-tumorigenic cells (201B7-NS/PCs). Although the DNA methylation profiles of 253G1-hiPSCs and 201B7-hiPSCs were similar regardless of passage number, the methylation status of the global DNA methylation profiles of 253G1-NS/PCs and 201B7-NS/PCs differed; the genomic regions surrounding the transcriptional start site of the CAT and PSMD5 genes were hypermethylated in 253G1-NS/PCs but not in 201B7-NS/PCs. Interestingly, the aberrant DNA methylation profile was more pronounced in 253G1-NS/PCs that had been passaged more than 15 times. In addition, we identified aberrations in DNA methylation at the RBP1 gene locus; the DNA methylation frequency in RBP1 changed as 253G1-NS/PCs were sequentially passaged. These results indicate that different NS/PC clones have different DNA methylomes and that DNA methylation patterns are unstable as cells are passaged. Therefore, DNA methylation profiles should be included in the criteria used to evaluate the tumorigenicity of hiPSC-NS/PCs in the clinical setting. Stem Cells 2017;35:1316-1327.

Focusing the management of rectal cancer

Rectal cancer treatment has undergone major changes over the last 15 years with a focus on individualized care based around MRI assessment of the relationship of the tumour to the mesorectal fascia, improved surgical techniques and targeted use of pre-operative oncological therapies in patients with locally advanced disease. The recognition that some tumours responded completely to pre-operative chemoradiotherapy, and the selective use of a non-operative policy has led to a quest to further identify those patients and their tumour in whom this approach could be used, irrespective of MRI stage. With no clear patient factors identified, the tumour and its gene expression has become a target for research to identify individual single-nucleotide polymorphisms, which may indicate a response to specific treatment, or not. To date some agents have been identified and trialed, such as cetuximab, with individual tumours being assessed for response allowing directed treatment. The reviewed paper by Sebio and colleagues report a study that links polymorphisms in the DNA repair gene XRCC1 with response to neoadjuvant 5-Fluorouracil treatment in rectal cancer patients. However, genetic heterogeneity alone may not explain the variations of drug response and environmental factors may lead to epigenetic effects and therefore alter responses. Therefore whilst this study demonstrates the impact of different single nucleotide polymorphisms (SNPs), it is only one step forward, but perhaps a step in the right direction.

Individualized risk prediction of outcomes for oral cavity cancer patients

BACKGROUND

Optimal management of oral cancer relies upon accurate and individualized risk prediction of relevant clinical outcomes. Individualized prognostic calculators have been developed to guide patient-physician communication and treatment-related decision-making. However it is critical to scrutinize their accuracy prior to integrating into clinical care.

AIM

To compare and evaluate oral cavity cancer prognostic calculators using an independent dataset.

METHODS

Five prognostic calculators incorporating patient and tumor characteristics were identified that evaluated five-year overall survival. A total of 505 patients with previously untreated oral cancer diagnosed between 2003 and 2014 were analyzed. Calculators were applied to each patient to generate individual predicted survival probabilities. Predictions were compared among prognostic tools and with observed outcomes using Kaplan-Meier plots, ROC curves and calibration plots.

RESULTS

Correlation between the five calculators varied from 0.59 to 0.86. There were considerable differences between individual predictions from pairs of calculators, with as many as 64% of patients having predictions that differed by more than 10%. Four of five calculators were well calibrated. For all calculators the predictions were associated with survival outcomes. The area under the ROC curve ranged from 0.65 to 0.71, with C-indices ranging from 0.63 to 0.67. An average of the 5 predictions had slightly better performance than any individual calculator.

CONCLUSION

Five prognostic calculators designed to predict individual outcomes of oral cancer differed significantly in their assessments of risk. Most were well calibrated and had modest discriminatory ability. Given the increasing importance of individualized risk prediction, more robust models are needed.

DNA methylation-based variation between human populations

Several studies have proved that DNA methylation affects regulation of gene expression and development. Epigenome-wide studies have reported variation in methylation patterns between populations, including Caucasians, non-Caucasians (Blacks), Hispanics, Arabs, and numerous populations of the African continent. Not only has DNA methylation differences shown to impact externally visible characteristics, but is also a potential biomarker for underlying racial health disparities between human populations. Ethnicity-related methylation differences set their mark during early embryonic development. Genetic variations, such as single-nucleotide polymorphisms and environmental factors, such as age, dietary folate, socioeconomic status, and smoking, impacts DNA methylation levels, which reciprocally impacts expression of phenotypes. Studies show that it is necessary to address these external influences when attempting to differentiate between populations since the relative impacts of these factors on the human methylome remain uncertain. The present review summarises several reported attempts to establish the contribution of differential DNA methylation to natural human variation, and shows that DNA methylation could represent new opportunities for risk stratification and prevention of several diseases amongst populations world-wide. Variation of methylation patterns between human populations is an exciting prospect which inspires further valuable research to apply the concept in routine medical and forensic casework. However, trans-generational inheritance needs to be quantified to decipher the proportion of variation contributed by DNA methylation. The future holds thorough evaluation of the epigenome to understand quantification, heritability, and the effect of DNA methylation on phenotypes. In addition, methylation profiling of the same ethnic groups across geographical locations will shed light on conserved methylation differences in populations.

Multiple imputation of missing covariates for the Cox proportional hazards cure model

We explore several approaches for imputing partially observed covariates when the outcome of interest is a censored event time and when there is an underlying subset of the population that will never experience the event of interest. We call these subjects 'cured', and we consider the case where the data are modeled using a Cox proportional hazards (CPH) mixture cure model. We study covariate imputation approaches using fully conditional specification. We derive the exact conditional distribution and suggest a sampling scheme for imputing partially observed covariates in the CPH cure model setting. We also propose several approximations to the exact distribution that are simpler and more convenient to use for imputation. A simulation study demonstrates that the proposed imputation approaches outperform existing imputation approaches for survival data without a cure fraction in terms of bias in estimating CPH cure model parameters. We apply our multiple imputation techniques to a study of patients with head and neck cancer. Copyright © 2016 John Wiley & Sons, Ltd.