Search

Abstract A scalable platform for cell typing in the glioma microenvironment can improve tumor subtyping and immune landscape detection as successful immunotherapy strategies continue to be sought and evaluated. DNA methylation (DNAm) biomarkers for molecular classification of tumor subtypes have been developed for clinical use. However, tools that predict the cellular landscape of the tumor are not well-defined or readily available. We developed the Glioma Immune Microenvironment Composition Calculator (GIMiCC), an approach for deconvolution of cell types in gliomas using DNAm data. Using data from 17 isolated cell types, we describe the derivation of the deconvolution libraries in the biological context of selected genomic regions and validate deconvolution results using independent datasets. We utilize GIMiCC to illustrate that DNAm-based estimates of immune composition are clinically relevant and scalable for potential clinical implementation. In addition, we utilize GIMiCC to identify composition-independent DNAm alterations that are associated with high immune infiltration. Our future work aims to optimize GIMiCC and advance the clinical evaluation of glioma.

Abstract Background Over the past decade, DNA methylation (DNAm)-based deconvolution methods that leverage cell-specific DNAm markers of immune cell types have been developed to provide accurate estimates of the proportions of leukocytes in peripheral blood. Immune cell phenotyping using DNAm markers, termed immunomethylomics or methylation cytometry, offers a solution for determining the body’s immune cell landscape that does not require fresh blood and is scalable to large sample sizes. Despite significant advances in DNAm-based deconvolution, references at the population level are needed for clinical and research interpretation of these additional immune layers. Here we aim to provide some references for immune populations in a group of multi-ethnic post-menopausal American women. Results We applied DNAm-based deconvolution to a large sample of post-menopausal women enrolled in the Women’s Health Initiative (baseline, N = 58) or the ancillary Long Life Study (WHI-LLS, N = 1237) to determine the reference ranges of 58 immune parameters, including proportions and absolute counts for 19 leukocyte subsets and 20 derived cell ratios. Participants were 50–94 years old at the time of blood draw, and N = 898 (69.3%) self-identified as White. Using linear regression models, we observed significant associations between age at blood draw and absolute counts and proportions of naïve B, memory CD4+, naïve CD4+, naïve CD8+, memory CD8+ memory, neutrophils, and natural killer cells. We also assessed the same immune profiles in a subset of paired longitudinal samples collected 14–18 years apart across N = 52 participants. Our results demonstrate high inter-individual variability in rates of change of leukocyte subsets over this time. And, when conducting paired t tests to test the difference in counts and proportions between the baseline visit and LLS visit, there were significant changes in naïve B, memory CD4+, naïve CD4+, naïve CD8+, memory CD8+ cells and neutrophils, similar to the results seen when analyzing the association with age in the entire cohort. Conclusions Here, we show that derived cell counts largely reflect the immune profile associated with proportions and that these novel methods replicate the known immune profiles associated with age. Further, we demonstrate the value this methylation cytometry approach can add as a potential application in epidemiological studies.

Numerous methods and approaches have been developed for generating time-to-event data from the Cox Proportional Hazards (CPH) model; however, they often require specification of a parametric distribution for the baseline hazard even though the CPH model itself makes no assumptions on the distribution of the baseline hazards. In line with the semi-parametric nature of the CPH model, a recently proposed method called the Flexible Hazards Method generates time-to-event data from a CPH model using a non-parametric baseline hazard function. While the initial results of this method are promising, it has not yet been comprehensively assessed with increasing covariates or against data generated under parametric baseline hazards. To fill this gap, we conducted a comprehensive study to benchmark the performance of the Flexible Hazards Method for generating data from a CPH model against parametric methods. Our results showed that with a single covariate and large enough assumed maximum time, the bias in the Flexible Hazards Method is 0.02 (with respect to the log hazard ratio) with a 95% confidence interval having coverage of 84.4%. This bias increases to 0.054 when there are 10 covariates under the same settings and the coverage of the 95% confidence interval decreases to 46.7%. In this paper, we explain the plausible reasons for this observed increase in bias and decrease in coverage as the number of covariates are increased, both empirically and theoretically, and provide readers and potential users of this method with some suggestions on how to best address these issues. In summary, the Flexible Hazards Method performs well when there are few covariates and the user wishes to simulate data from a non-parametric baseline hazard.

JQ1 and GSK2801 are bromo domain inhibitors (BDI) known to exhibit enhanced anti-cancer activity when combined with other agents. However, the underlying molecular mechanisms behind such enhanced activity remain unclear. We used network-pharmacology approaches to understand the shared molecular mechanisms behind the enhanced activity of JQ1 and GSK2801 when used together to treat breast cancer (BC). The gene targets of JQ1 and GSK2801 were intersected with known BC-targets and their putative targets against BC were derived. The key genes were explored through gene-ontology-enrichment, Protein-Protein-Interaction (PPI) networking, survival analysis, and molecular modeling simulations. The genes, CTSB, MAPK14, MET, PSEN2 and STAT3, were found to be common targets for both drugs. In total, 49 biological processes, five molecular functions and 61 metabolic pathways were similarly enriched for JQ1 and GSK2801 BC targets among which several terms are related to cancer: IL-17, TNF and JAK-STAT signaling pathways. Survival analyses revealed that all five putative synergistic targets are significantly associated with survival in BC (log-rank p

Abstract Background Our randomized controlled clinical trial will explore the potential of bazedoxifene plus conjugated estrogen to modulate breast tissue-based risk biomarkers as a surrogate for breast cancer risk reduction. This paper investigates the statistical design features of the trial and the rationale for the final choice of its design. Group sequential designs are a popular design approach to allow a trial to stop early for success or futility, potentially saving time and money over a fixed trial design. While Bayesian adaptive designs enjoy the same properties as group sequential designs, they have the added benefit of using prior information as well as inferential interpretation conditional on the data. Whether a frequentist or Bayesian trial, most adaptive designs have interim analyses that allow for early stopping, typically utilizing only the primary endpoint. A drawback to this approach is that the study may not have enough data for adequate comparisons of a single, key secondary endpoint. This can happen, for example, if the secondary endpoint has a smaller effect than the primary endpoint. Methods In this paper, we investigate a trial design called two-endpoint adaptive, which stops early only if a criterion is met for primary and secondary endpoints. The approach focuses the final analysis on the primary endpoint but ensures adequate data for the secondary analysis. Our study has two arms with a primary (change in mammographic fibroglandular volume) and secondary endpoint (change in mammary tissue Ki-67). Results We present operating characteristics including power, trial duration, and type I error rate and discuss the value and risks of modeling Bayesian group sequential designs with primary and secondary endpoints, comparing against alternative designs. The results indicate that the two-endpoint adaptive design has better operating characteristics than competing designs if one is concerned about having adequate information for a key secondary endpoint. Discussion Our approach balances trial speed and the need for information on the single, key secondary endpoint.

Abstract Background Cellular compositions of solid tumor microenvironments are heterogeneous, varying across patients and tumor types. High-resolution profiling of the tumor microenvironment cell composition is crucial to understanding its biological and clinical implications. Previously, tumor microenvironment gene expression and DNA methylation-based deconvolution approaches have been shown to deconvolve major cell types. However, existing methods lack accuracy and specificity to tumor type and include limited identification of individual cell types. Results We employed a novel tumor-type-specific hierarchical model using DNA methylation data to deconvolve the tumor microenvironment with high resolution, accuracy, and specificity. The deconvolution algorithm is named HiTIMED. Seventeen cell types from three major tumor microenvironment components can be profiled (tumor, immune, angiogenic) by HiTIMED, and it provides tumor-type-specific models for twenty carcinoma types. We demonstrate the prognostic significance of cell types that other tumor microenvironment deconvolution methods do not capture. Conclusion We developed HiTIMED, a DNA methylation-based algorithm, to estimate cell proportions in the tumor microenvironment with high resolution and accuracy. HiTIMED deconvolution is amenable to archival biospecimens providing high-resolution profiles enabling to study of clinical and biological implications of variation and composition of the tumor microenvironment.

Selenium is an important micronutrient for foetal development. MicroRNAs play an important role in the function of the placenta, in communication between the placenta and maternal systems, and their expression can be altered through environmental and nutritional cues. To investigate the associations between placental selenium concentration and microRNA expression in the placenta, our observational study included 393 mother-child pairs from the New Hampshire Birth Cohort Study (NHBCS) and the Rhode Island Child Health Study (RICHS). Placental selenium concentrations were quantified using inductively coupled plasma mass spectrometry, and microRNA transcripts were measured using RNA-seq. We fit negative binomial additive models for assessing the association between selenium and microRNAs. We used the microRNA Data Integration Portal (mirDIP) to predict the target mRNAs of the differentially expressed microRNAs and verified the relationships between miRNA and mRNA targets in a subset of samples using existing whole transcriptome data (N = 199). We identified a non-monotonic association between selenium concentration and the expression of miR-216a-5p/miR-217-5p cluster (effective degrees of freedom, EDF = 2.44 and 2.08; FDR = 3.08 × 10−5) in placenta. Thirty putative target mRNAs of miR-216a-5p and/or miR-217-5p were identified computationally and empirically and were enriched in selenium metabolic pathways (driven by selenoprotein coding genes, TXNRD2 and SELENON). Our findings suggest that selenium influences placental microRNA expression. Further, miR-216a-5p and its putative target mRNAs could be the potential mechanistic targets of the health effect of selenium.

Glucocorticoids, such as dexamethasone, are used as anti-inflammatory and immunosuppressive drugs, however patients may exhibit resistance or side effects. Here the authors propose that a dexamethasone related neutrophil-specific DNA methylation index can be used as a marker of glucocorticoid exposure and response and as a prognostic factor in brain tumor survival.

IntroductionThe human brain comprises heterogeneous cell types whose composition can be altered with physiological and pathological conditions. New approaches to discern the diversity and distribution of brain cells associated with neurological conditions would significantly advance the study of brain-related pathophysiology and neuroscience. Unlike single-nuclei approaches, DNA methylation-based deconvolution does not require special sample handling or processing, is cost-effective, and easily scales to large study designs. Existing DNA methylation-based methods for brain cell deconvolution are limited in the number of cell types deconvolvedMethodsUsing DNA methylation profiles of the top cell-type-specific differentially methylated CpGs, we employed a hierarchical modeling approach to deconvolve GABAergic neurons, glutamatergic neurons, astrocytes, microglial cells, oligodendrocytes, endothelial cells, and stromal cells.ResultsWe demonstrate the utility of our method by applying it to data on normal tissues from various brain regions and in aging and diseased tissues, including Alzheimer’s disease, autism, Huntington’s disease, epilepsy, and schizophrenia.DiscussionWe expect that the ability to determine the cellular composition in the brain using only DNA from bulk samples will accelerate understanding brain cell type composition and cell-type-specific epigenetic states in normal and diseased brain tissues.

Abstract Background Triple-negative breast cancer (TNBC) constitutes 10–20% of breast cancers and is challenging to treat due to a lack of effective targeted therapies. Previous studies in TNBC cell lines showed in vitro growth inhibition when JQ1 or GSK2801 were administered alone, and enhanced activity when co-administered. Given their respective mechanisms of actions, we hypothesized the combinatorial effect could be due to the target genes affected. Hence the target genes were characterized for their expression in the TNBC cell lines to prove the combinatorial effect of JQ1 and GSK2801. Methods RNASeq data sets of TNBC cell lines (MDA-MB-231, HCC-1806 and SUM-159) were analyzed to identify the differentially expressed genes in single and combined treatments. The topmost downregulated genes were characterized for their downregulated expression in the TNBC cell lines treated with JQ1 and GSK2801 under different dose concentrations and combinations. The optimal lethal doses were determined by cytotoxicity assays. The inhibitory activity of the drugs was further characterized by molecular modelling studies. Results Global expression profiling of TNBC cell lines using RNASeq revealed different expression patterns when JQ1 and GSK2801 were co-administered. Functional enrichment analyses identified several metabolic pathways (i.e., systemic lupus erythematosus, PI3K-Akt, TNF, JAK-STAT, IL-17, MAPK, Rap1 and signaling pathways) enriched with upregulated and downregulated genes when combined JQ1 and GSK2801 treatment was administered. RNASeq identified downregulation of PTPRC, MUC19, RNA5-8S5, KCNB1, RMRP, KISS1 and TAGLN (validated by RT-qPCR) and upregulation of GPR146, SCARA5, HIST2H4A, CDRT4, AQP3, MSH5-SAPCD1, SENP3-EIF4A1, CTAGE4 and RNASEK-C17orf49 when cells received both drugs. In addition to differential gene regulation, molecular modelling predicted binding of JQ1 and GSK2801 with PTPRC, MUC19, KCNB1, TAGLN and KISS1 proteins, adding another mechanism by which JQ1 and GSK2801 could elicit changes in metabolism and proliferation. Conclusion JQ1-GSK2801 synergistically inhibits proliferation and results in selective gene regulation. Besides suggesting that combinatorial use could be useful therapeutics for the treatment of TNBC, the findings provide a glimpse into potential mechanisms of action for this combination therapy approach.

Background: DNA methylation markers have been associated with lung cancer risk and may identify aetiologically relevant genomic regions, or alternatively, be markers of disease risk factors or biological processes associated with disease development. Methods: In a nested case–control study, we measured blood leukocyte DNA methylation levels in pre-diagnostic samples collected from 430 participants (208 cases; 222 controls) in the 1989 CLUE II cohort. We compared DNA methylation levels with case/control status to identify novel genomic regions, both single CpG sites and differentially methylated regions (DMRs), while controlling for known DNA methylation changes associated with smoking using a previously described pack-years-based smoking methylation score. Stratification analyses were conducted over time from blood draw to diagnosis, histology, and smoking status. Results: We identified 16 single CpG sites and 40 DMRs significantly associated with lung cancer risk (q

The mammalian kidney has extensive repair capacity; however, identifying adult renal stem cells has proven elusive. We applied an epigenetic marker of foetal cell origin (FCO) in diverse human tissues as a probe for developmental cell persistence, finding a 5.4-fold greater FCO proportion in kidney. Normal kidney FCO proportions averaged 49% with extensive interindividual variation. FCO proportions were significantly negatively correlated with immune-related gene expression and positively correlated with genes expressed in the renal medulla, including those involved in renal organogenesis (e.g., FGF2, PAX8, and HOXB7). FCO associated genes also mapped to medullary nephron segments in mouse and rat, suggesting evolutionary conservation of this cellular compartment. Renal cancer patients whose tumours contained non-zero FCO scores survived longer. The kidney appears unique in possessing substantial foetal epigenetic features. Further study of FCO-related gene methylation may elucidate regenerative regulatory programmes in tissues without apparent discrete stem cell compartments.

Deconvolution algorithms facilitate studying cell type-specific changes using bulk data from complex tissues. Here, the authors present a deconvolution method that predicts DNA methylation levels in 12 leukocyte subtypes using human microarray data and apply it to various examples.

DNA methylation (DNAm) age may reflect age-related variations in biological changes and abnormalities related to ageing. DNAm age acceleration measures have been associated with a number of cancers, but to our knowledge, have not been examined in relation to pancreatic cancer risk or survival. DNAm levels in leukocytes of prediagnostic blood samples of 393 pancreatic cancer cases and 431 matched controls, pooled from three large prospective cohort studies, were used to estimate DNAm age, epigenetic age acceleration (AA), and intrinsic epigenetic age acceleration (IEAA) metrics. Logistic regression and Cox proportional hazard regression models were used to examine the relationship between the various AA and IEAA metrics and pancreatic cancer risk and survival, respectively. The results showed that pancreatic cancer risk was significantly increased across all IEAA metrics, ranging from 83% to 95% increased risk when comparing the third and highest quartiles to the lowest quartile of IEAA. Consistent with these findings, the results from multivariate spline regression analyses showed non-linear relationships between all three IEAA metrics and pancreatic cancer risk with apparent threshold effect including two turning points at minimal and at maximal risks, respectively. There is no evidence of a significant association between pancreatic cancer survival and any of the epigenetic AA or IEAA metrics. Our results indicate DNAm age acceleration, measured in blood prior to cancer diagnosis, is associated with an increased risk of pancreatic cancer in a complex nonlinear, dose–response manner. Epigenetic IEAA metrics may be a useful addition to current methods for pancreatic cancer risk prediction.

Background: It is important to identify when two exposures impact a molecular marker (e.g., a gene’s expression) in similar ways, for example, to learn that a new drug has a similar effect to an existing drug. Currently, statistically robust approaches for making comparisons of equivalence of effect sizes obtained from two independently run treatment vs. control comparisons have not been developed.Results: Here, we propose two approaches for evaluating the question of equivalence between effect sizes of two independent studies: a bootstrap test of the Equivalent Change Index (ECI), which we previously developed, and performing Two One-Sided t-Tests (TOST) on the difference in log-fold changes directly. The ECI of a gene is computed by taking the ratio of the effect size estimates obtained from the two different studies, weighted by the maximum of the two p-values and giving it a sign indicating if the effects are in the same or opposite directions, whereas TOST is a test of whether the difference in log-fold changes lies outside a region of equivalence. We used a series of simulation studies to compare the two tests on the basis of sensitivity, specificity, balanced accuracy, and F1-score. We found that TOST is not efficient for identifying equivalently changed gene expression values (F1-score = 0) because it is too conservative, while the ECI bootstrap test shows good performance (F1-score = 0.95). Furthermore, applying the ECI bootstrap test and TOST to publicly available microarray expression data from pancreatic cancer showed that, while TOST was not able to identify any equivalently or inversely changed genes, the ECI bootstrap test identified genes associated with pancreatic cancer. Additionally, when investigating publicly available RNAseq data of smoking vs. vaping, no equivalently changed genes were identified by TOST, but ECI bootstrap test identified genes associated with smoking.Conclusion: A bootstrap test of the ECI is a promising new statistical approach for determining if two diverse studies show similarity in the differential expression of genes and can help to identify genes which are similarly influenced by a specific treatment or exposure. The R package for the ECI bootstrap test is available at https://github.com/Hecate08/ECIbootstrap.

The Accelerate Cancer Education (ACE) summer research program at The University of Kansas Cancer Center (KUCC) is a six-week, cancer-focused, summer research experience for high school students from historically marginalized populations in the Kansas City metropolitan area. Cancer affects all populations and continues to be the second leading cause of death in the United States, and a large number of disparities impact racial and ethnic minorities, including increased cancer incidence and mortality. Critically, strategies to bolster diversity, equity, inclusion, and accessibility are needed to address persistent cancer disparities. The ACE program offers an educational opportunity for a population of students who otherwise would not have easy access onto a medical center campus to make connections with cancer physicians and researchers and provides a vital response to the need for a more diverse and expansive oncology workforce. Students grow their technical, social, and professional skills and develop self-efficacy and long-lasting connections that help them matriculate and persist through post-secondary education. Developed in 2018, the ACE program has trained 37 high school junior and senior students. This article describes the need for and how we successfully developed and implemented the ACE program.

Reference-based deconvolution methods use reference libraries of cell-specific DNA methylation (DNAm) measurements as a means toward deconvoluting cell proportions in heterogeneous biospecimens (e.g., whole-blood). As the accuracy of such methods depends highly on the CpG loci comprising the reference library, recent research efforts have focused on the selection of libraries to optimize deconvolution accuracy. While existing approaches for library selection work extremely well, the best performing approaches require a training data set consisting of both DNAm profiles over a heterogeneous cell population and gold-standard measurements of cell composition (e.g., flow cytometry) in the same samples. Here, we present a framework for reference library selection without a training dataset (RESET) and benchmark it against the Legacy method (minfi:pickCompProbes), where libraries are constructed based on a pre-specified number of cell-specific differentially methylated loci (DML). RESET uses a modified version of the Dispersion Separability Criteria (DSC) for comparing different libraries and has four main steps: 1) identify a candidate set of cell-specific DMLs, 2) randomly sample DMLs from the candidate set, 3) compute the Modified DSC of the selected DMLs, and 4) update the selection probabilities of DMLs based on their contribution to the Modified DSC. Steps 2–4 are repeated many times and the library with the largest Modified DSC is selected for subsequent reference-based deconvolution. We evaluated RESET using several publicly available datasets consisting of whole-blood DNAm measurements with corresponding measurements of cell composition. We computed the RMSE and R2 between the predicted cell proportions and their measured values. RESET outperformed the Legacy approach in selecting libraries that improve the accuracy of deconvolution estimates. Additionally, reference libraries constructed using RESET resulted in cellular composition estimates that explained more variation in DNAm as compared to the Legacy approach when evaluated in the context of epigenome-wide association studies (EWAS) of several publicly available data sets. This finding has implications for the statistical power of EWAS. RESET combats potential challenges associated with existing approaches for reference library assembly and thus, may serve as a viable strategy for library construction in the absence of a training data set.

Abstract Background A low level of methylation at cg05575921 in the aryl-hydrocarbon receptor repressor (AHRR) gene is robustly associated with smoking, and some studies have observed associations between cg05575921 methylation and increased lung cancer risk and mortality. To prospectively examine whether decreased methylation at cg05575921 may identify high risk subpopulations for lung cancer screening among heavy smokers, and mortality in cases, we evaluated associations between cg05575921 methylation and lung cancer risk and mortality, by histotype, in heavy smokers. Methods The β-Carotene and Retinol Efficacy Trial (CARET) included enrollees ages 45–69 with ≥ 20 pack-year smoking histories and/or occupational asbestos exposure. A subset of CARET participants had cg05575921 methylation available from HumanMethylationEPIC assays of blood collected on average 4.3 years prior to lung cancer diagnosis in cases. Cg05575921 methylation β-values were treated continuously for a 10% methylation decrease and as quintiles, where quintile 1 (Q1, referent) represents high methylation and Q5, low methylation. We used conditional logistic regression models to examine lung cancer risk overall and by histotype in a nested case-control study including 316 lung cancer cases (diagnosed through 2005) and 316 lung cancer-free controls matched on age (±5 years), sex, race/ethnicity, enrollment year, current/former smoking, asbestos exposure, and follow-up time. Mortality analyses included 372 lung cancer cases diagnosed between 1985 and 2013 with available methylation data. We used Cox proportional hazards models to examine mortality overall and by histotype. Results Decreased cg05575921 methylation was strongly associated with smoking, even in our population of heavy smokers. We did not observe associations between decreased pre-diagnosis cg05575921 methylation and increased lung cancer risk, overall or by histotype. We observed linear increasing trends for lung cancer-specific mortality across decreasing cg05575921 methylation quintiles for adenocarcinoma and small cell carcinoma (P-trends = 0.01 and 0.04, respectively). Conclusions In our study of heavy smokers, decreased cg05575921 methylation was strongly associated with smoking but not increased lung cancer risk. The observed association between cg05575921 methylation and increased mortality in adenocarcinoma and small cell histotypes requires further examination. Our results do not support using decreased cg05575921 methylation as a biomarker for lung cancer screening risk stratification.

Abstract Background In silico functional genomics have become a driving force in the way we interpret and use gene expression data, enabling researchers to understand which biological pathways are likely to be affected by the treatments or conditions being studied. There are many approaches to functional genomics, but a number of popular methods determine if a set of modified genes has a higher than expected overlap with genes known to function as part of a pathway (functional enrichment testing). Recently, researchers have started to apply such analyses in a new way: to ask if the data they are collecting show similar disruptions to biological functions compared to reference data. Examples include studying whether similar pathways are perturbed in smokers vs. users of e-cigarettes, or whether a new mouse model of schizophrenia is justified, based on its similarity in cytokine expression to a previously published model. However, there is a dearth of robust statistical methods for testing hypotheses related to these questions and most researchers resort to ad hoc approaches. The goal of this work is to develop a statistical approach to identifying gene pathways that are equivalently (or inversely) changed across two experimental conditions. Results We developed Equivalent Change Enrichment Analysis (ECEA). This is a new type of gene enrichment analysis based on a statistic that we call the equivalent change index (ECI). An ECI of 1 represents a gene that was over or under-expressed (compared to control) to the same degree across two experiments. Using this statistic, we present an approach to identifying pathways that are changed in similar or opposing ways across experiments. We compare our approach to current methods on simulated data and show that ECEA is able to recover pathways exhibiting such changes even when they exhibit complex patterns of regulation, which other approaches are unable to do. On biological data, our approach recovered pathways that appear directly connected to the condition being studied. Conclusions ECEA provides a new way to perform gene enrichment analysis that allows researchers to compare their data to existing datasets and determine if a treatment will cause similar or opposing genomic perturbations.

Abstract Background Exposure to the herbicide Agent Orange during the Vietnam War was widespread and is associated with numerous adverse health outcomes. A continuing concern of veterans is the possibility that exposure to the dioxin-containing herbicide might induce adverse reproductive outcomes. We sought to assess whether exposure to Agent Orange in Vietnam was associated with changes in DNA methylation in sperm in a subset of Vietnam veterans who participated in the Air Force Health Study (AFHS). Methods We studied 37 members of the AFHS chosen to have no, low, medium or high exposure to Agent Orange, based upon serum dioxin levels obtained during a series of examinations. DNA from stored semen was extracted and DNA methylation assessed on the Illumina 450 K platform. Results Initial epigenome-wide analysis returned no loci that survived control for false discovery. However, the TEAD3 gene had four different CpG sites that showed loss of DNA methylation associated with dioxin exposure. Analysis assessing regional DNA methylation changes revealed 36 gene regions, including the region of the imprinted gene H19 to have altered DNA methylation associated with high exposure compared to the low exposure group. Additional comparison of our data with sperm DNA methylation data from Russian boys exposed to dioxin found an additional 5 loci that were altered in both studies and exhibited a consistent direction of association. Conclusions Studying a small number of sperm samples from veterans enrolled in the AFHS, we did not find evidence of significant epigenome-wide alterations associated with exposure to Agent Orange. However, additional analysis showed that the H19 gene region is altered in the sperm of Agent Orange-exposed Ranch Hand veterans. Our study also replicated several findings of a prior study of dioxin-exposed Russian boys. These results provide additional candidate loci for further investigation and may have implications for the reproductive health of dioxin-exposed individuals.

Abstract Background Differentiated cells that arise from stem cells in early development contain DNA methylation features that provide a memory trace of their fetal cell origin (FCO). The FCO signature was developed to estimate the proportion of cells in a mixture of cell types that are of fetal origin and are reminiscent of embryonic stem cell lineage. Here we implemented the FCO signature estimation method to compare the fraction of cells with the FCO signature in tumor tissues and their corresponding nontumor normal tissues. Methods We applied our FCO algorithm to discovery data sets obtained from The Cancer Genome Atlas (TCGA) and replication data sets obtained from the Gene Expression Omnibus (GEO) data repository. Wilcoxon rank sum tests, linear regression models with adjustments for potential confounders and non-parametric randomization-based tests were used to test the association of FCO proportion between tumor tissues and nontumor normal tissues. P-values of

Abstract Background When designing an epigenome-wide association study (EWAS) to investigate the relationship between DNA methylation (DNAm) and some exposure(s) or phenotype(s), it is critically important to assess the sample size needed to detect a hypothesized difference with adequate statistical power. However, the complex and nuanced nature of DNAm data makes direct assessment of statistical power challenging. To circumvent these challenges and to address the outstanding need for a user-friendly interface for EWAS power evaluation, we have developed pwrEWAS. Results The current implementation of pwrEWAS accommodates power estimation for two-group comparisons of DNAm (e.g. case vs control, exposed vs non-exposed, etc.), where methylation assessment is carried out using the Illumina Human Methylation BeadChip technology. Power is calculated using a semi-parametric simulation-based approach in which DNAm data is randomly generated from beta-distributions using CpG-specific means and variances estimated from one of several different existing DNAm data sets, chosen to cover the most common tissue-types used in EWAS. In addition to specifying the tissue type to be used for DNAm profiling, users are required to specify the sample size, number of differentially methylated CpGs, effect size(s) (Δ β ), target false discovery rate (FDR) and the number of simulated data sets, and have the option of selecting from several different statistical methods to perform differential methylation analyses. pwrEWAS reports the marginal power, marginal type I error rate, marginal FDR, and false discovery cost (FDC). Here, we demonstrate how pwrEWAS can be applied in practice using a hypothetical EWAS. In addition, we report its computational efficiency across a variety of user settings. Conclusion Both under- and overpowered studies unnecessarily deplete resources and even risk failure of a study. With pwrEWAS, we provide a user-friendly tool to help researchers circumvent these risks and to assist in the design and planning of EWAS. Availability The web interface is written in the R statistical programming language using Shiny (RStudio Inc., 2016) and is available at https://biostats-shinyr.kumc.edu/pwrEWAS/. The R package for pwrEWAS is publicly available at GitHub (https://github.com/stefangraw/pwrEWAS).

Abstract Background Unsupervised clustering represents one of the most widely applied methods in analysis of high-throughput ‘omics data. A variety of unsupervised model-based or parametric clustering methods and non-parametric clustering methods have been proposed for RNA-seq count data, most of which perform well for large samples, e.g. N ≥ 500. A common issue when analyzing limited samples of RNA-seq count data is that the data follows an over-dispersed distribution, and thus a Negative Binomial likelihood model is often used. Thus, we have developed a Negative Binomial model-based (NBMB) clustering approach for application to RNA-seq studies. Results We have developed a Negative Binomial Model-Based (NBMB) method to cluster samples using a stochastic version of the expectation-maximization algorithm. A simulation study involving various scenarios was completed to compare the performance of NBMB to Gaussian model-based or Gaussian mixture modeling (GMM). NBMB was also applied for the clustering of two RNA-seq studies; type 2 diabetes study (N = 96) and TCGA study of ovarian cancer (N = 295). Simulation results showed that NBMB outperforms GMM applied with different transformations in majority of scenarios with limited sample size. Additionally, we found that NBMB outperformed GMM for small clusters distance regardless of sample size. Increasing total number of genes with fixed proportion of differentially expressed genes does not change the outperformance of NBMB, but improves the overall performance of GMM. Analysis of type 2 diabetes and ovarian cancer tumor data with NBMB found good agreement with the reported disease subtypes and the gene expression patterns. This method is available in an R package on CRAN named NB.MClust. Conclusion Use of Negative Binomial model based clustering is advisable when clustering over dispersed RNA-seq count data.

Breast Cancer Metastasis Suppressor 1 (BRMS1) expression is associated with longer patient survival in multiple cancer types. Understanding BRMS1 functionality will provide insights into both mechanism of action and will enhance potential therapeutic development. In this study, we confirmed that the C-terminus of BRMS1 is critical for metastasis suppression and hypothesized that critical protein interactions in this region would explain its function. Phosphorylation status at S237 regulates BRMS1 protein interactions related to a variety of biological processes, phenotypes [cell cycle (e.g., CDKN2A), DNA repair (e.g., BRCA1)], and metastasis [(e.g., TCF2 and POLE2)]. Presence of S237 also directly decreased MDA-MB-231 breast carcinoma migration in vitro and metastases in vivo. The results add significantly to our understanding of how BRMS1 interactions with Sin3/HDAC complexes regulate metastasis and expand insights into BRMS1’s molecular role, as they demonstrate BRMS1 C-terminus involvement in distinct protein-protein interactions.

Background: Oral microbiota is believed to play important roles in systemic diseases, including cancer. Methods: We collected oral samples (tongue, buccal, supragingival, and saliva) and pancreatic tissue or intestinal samples from 52 subjects, and characterized 16S rRNA genes using high-throughput DNA sequencing. Results: Bray–Curtis plot showed clear separations between bacterial communities in the oral cavity and those in intestinal and pancreatic tissue samples. PERMANOVA tests indicated that bacterial communities from buccal samples were similar to supragingival and saliva samples, and pancreatic duct samples were similar to pancreatic tumor samples, but all other samples were significantly different from each other. A total of 73 unique Amplicon Sequence Variants (ASVs) were shared between oral and pancreatic or intestinal samples. Only four ASVs showed significant concordance, and two specific bacterial species (Gemella morbillorum and Fusobacterium nucleatum subsp. vincentii) showed consistent presence or absence patterns between oral and intestinal or pancreatic samples, after adjusting for within-subject correlation and disease status. Lastly, microbial co-abundance analyses showed distinct strain-level cluster patterns among microbiome members in buccal, saliva, duodenum, jejunum, and pancreatic tumor samples. Conclusions: Our findings indicate that oral, intestinal, and pancreatic bacterial microbiomes overlap but exhibit distinct co-abundance patterns in patients with pancreatic cancer and other gastrointestinal diseases.

1. Abstract: Early phase clinical trials are the first step in testing new medications and therapeutics developed by clinical and biomedical investigators. These trials aim to find a safe dose of a newly targeted drug (phase I) or find out more about the side effects and early signals of treatment efficacy (phase II). In a research institute, many biomedical investigators in oncology are encouraged to initiate such trials early in their careers as part of developing their research portfolio. These investigator-initiated trials (IITs) are funded internally by the University of Kansas Cancer Center or partially funded by pharmaceutical companies. As financial, administrative, and practical considerations play an essential role in the successful completion of IITs, it is imperative to efficiently allocate resources to plan, design, and execute these studies within the allotted time. This manuscript describes monitoring tools and processes to improve the efficiency, cost-effectivness, and reliability of IITs. The contributions of this team to processes such as: participant recruitment, feasibility analysis, clinical trial design, accrual monitoring, data management, interim analysis support, and final analysis and reporting are described in detail. This manuscript elucidates how, through the aid of technology and dedicated personnel support, the efficiency of IIT-related processes can be improved. Early results of these initiatives look promising, and the Biostatistics and Informatics team intends to continue fostering innovative methodologies to enhance cancer research by improving the efficiency of IITs.

Assessing DNA methylation profiles in human blood has become a major focus of epidemiologic inquiry. Understanding variability in CpG-specific DNA methylation over moderate periods of time is a critical first step in identifying CpG sites that are candidates for DNA methylation-based etiologic, diagnostic and prognostic predictors of pathogenesis. Using the Illumina MethylationEPIC [850K] BeadArray, DNA methylation was profiled in paired whole blood samples collected approximately 1 year apart from 35 healthy women enrolled in the Nurses Study II cohort. The median intraclass correlation coefficient (ICC) across all CpG loci was 0.19 [Interquartile Range (IQR) 0.00–0.50]; 74.8% of ICCs were in the low range (0–0.5), 16.9% in the mid-range of ICCs (0.5–0.8), and 8.3% in the high-range of ICCs (0.8–1). ICCs were similar for CpG probes on the 450K Illumina array (median 0.17) and the new probes added to the 850K array (median 0.21). ICCs for CpG loci on the sex chromosomes and known metastable epialleles were high (median 0.71, 0.97, respectively), and ICCs among methylation quantitative trait loci (mQTL) CpGs were significantly higher as compared to non-mQTL CpGs (median 0.73, 0.16, respectively, P

Abstract Genome-wide methylation arrays are powerful tools for assessing cell composition of complex mixtures. We compare three approaches to select reference libraries for deconvoluting neutrophil, monocyte, B-lymphocyte, natural killer, and CD4+ and CD8+ T-cell fractions based on blood-derived DNA methylation signatures assayed using the Illumina HumanMethylationEPIC array. The IDOL algorithm identifies a library of 450 CpGs, resulting in an average R2 = 99.2 across cell types when applied to EPIC methylation data collected on artificial mixtures constructed from the above cell types. Of the 450 CpGs, 69% are unique to EPIC. This library has the potential to reduce unintended technical differences across array platforms.

Genetic and epigenetic alterations are required for carcinogenesis and the mutation burden across tumor types has been investigated. Here, we investigate epigenetic alterations with a novel measure of global DNA methylation dysregulation, the methylation dysregulation index (MDI), across 14 cancer types in The Cancer Genome Atlas (TCGA) database. DNA methylation data—obtained using Illumina HumanMethylation450 BeadChip—was accessed from TCGA. We calculated the MDI in 14 tumor types (n = 5,592 tumors), using adjacent normal tissues (n = 701) from each tumor site. Copy number alteration, and mutation burden were retrieved from cBioportal (n = 5,152). We tested the relation of subject MDI across tumors and with age, gender, tumor stage, estimated tumor purity, and copy number alterations for both overall MDI and genomic-context-specific MDI. We also investigated the top most dysregulated loci shared across tumor types. There was a broad range of extent in methylation dysregulation across tumor types (P < 2.2E-16). However, a consistent pattern of methylation dysregulation stratified by genomic context was observed across tumor types where the highest dysregulation occurred at non-CpG island regions. Considering other summary measures of somatic alteration, MDI was correlated with copy number alterations but not with mutation burden. Using the top dysregulated CpG sites in common across tumors, 4 classes of cancer types were observed, and the functional consequences of these alterations to gene expression were confirmed. This work identified the global DNA methylation dysregulation patterns across 14 cancer types showing a higher impact for the non-CpG island areas. The most dysregulated loci across cancer types identified common clusters across cancer types that may have implications for future treatment and prevention measures.

Abstract Background Sex-specific factors play a major role in human health and disease, including responses to environmental stresses such as toxicant exposure. Increasing evidence suggests that such sex differences also exist during fetal development. In a previous report using the resources of the New Hampshire Birth Cohort Study (NHBCS), we found that low-to-moderate in utero exposure to arsenic, a highly toxic and widespread pollutant, was associated with altered expression of several key developmental genes in the fetal portion of the placenta. These associations were sex-dependent, suggesting that in utero arsenic exposure differentially impacts male and female fetuses. In the present study, we investigated the molecular basis for these sex-specific responses to arsenic. Methods Using NanoString technology, we further analyzed the fetal placenta samples from the NHBCS for the expression of genes encoding arsenic transporters and metabolic enzymes. Multivariable linear regression analysis was used to examine their relationship with arsenic exposure and with key developmental genes, after stratification by fetal sex. Results We found that maternal arsenic exposure was strongly associated with expression of the AQP9 gene, encoding an aquaglyceroporin transporter, in female but not male fetal placenta. Moreover, AQP9 expression associated with that of a subset of female-specific arsenic-responsive genes. Conclusions Our results suggest that AQP9 is upregulated in response to arsenic exposure in female, but not male, fetal placenta. Based on these results and prior studies, increased AQP9 expression may lead to increased arsenic transport in the female fetal placenta, which in turn may alter the expression patterns of key developmental genes that we have previously shown to be associated with arsenic exposure. Thus, this study suggests that AQP9 may play a role in the sex-specific effects of in utero arsenic exposure.

Coronavirus disease (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is one of the pressing contemporary public health challenges. Investigations into the genomic structure of SARS-CoV-2 may inform ongoing vaccine development efforts and/or provide insights into vaccine efficacy to fight against COVID-19. Evolutionary analysis of 540 genomes spanning 20 different countries/territories was conducted and revealed an increase in the genomic divergence across successive generations. The ancestor of the phylogeny was found to be the isolate from the 2019/2020 Wuhan outbreak. Its transmission was outlined across 20 countries/territories as per genomic similarity. Our results demonstrate faster evolving variations in the genomic structure of SARS-CoV-2 when compared to the isolates from early stages of the pandemic. Genomic alterations were predominantly located and mapped onto the reported vaccine candidates of structural genes, which are the main targets for vaccine candidates. S protein showed 34, N protein 25, E protein 2, and M protein 3 amino acid variations in 246 genomes among 540. Among identified mutations, 23 in S protein, 1 in E, 2 from M, and 7 from N protein were mapped with the reported vaccine candidates explaining the possible implications on universal vaccines. Hence, potential target regions for vaccines would be ideally chosen from the structural regions of the genome that lack high variation. The increasing variations in the genome of SARS-CoV-2 together with our observations in structural genes have important implications for the efficacy of a successful universal vaccine against SARS-CoV-2.

Although considerable evidence suggests that in utero arsenic exposure affects children's health, these data are mainly from areas of the world where groundwater arsenic levels far exceed the World Health Organization limit of 10 μg/L. We, and others, have found that more common levels of in utero arsenic exposure may also impact children's health. However, the underlying molecular mechanisms are poorly understood. To address this issue, we analyzed the expression of key developmental genes in fetal placenta in a birth cohort of women using unregulated water supplies in a US region with elevated groundwater arsenic. We identified several genes whose expression associated with maternal arsenic exposure in a fetal sex-specific manner. In particular, expression of the HEDGEHOG pathway component, GLI3, in female placentae was both negatively associated with arsenic exposure and positively associated with infant birth weight. This suggests that modulation of GLI3 in the fetal placenta, and perhaps in other fetal tissues, contributes to arsenic's detrimental effects on fetal growth. We showed previously that arsenic-exposed NIH3T3 cells have reduced GLI3 repressor protein. Together, these studies identify GLI3 as a key signaling node that is affected by arsenic, mediating a subset of its effects on developmental signaling and fetal health.