Letter to the Editor: on the stability and internal consistency of component-wise sparse mixture regression-based clustering.

Understanding the relationship between molecular markers and a phenotype of interest is often obfuscated by patient-level heterogeneity. To address this challenge, Chang et al. recently published a novel method called Component-wise Sparse Mixture Regression (CSMR), a regression-based clustering method that promises to detect heterogeneous relationships between molecular markers and a phenotype of interest under high-dimensional settings. In this Letter to the Editor, we raise awareness to several issues concerning the assessment of CSMR in Chang et al., particularly its assessment in settings where the number of features, P, exceeds the study sample size, N, and advocate for additional metrics/approaches when assessing the performance of regression-based clustering methodologies.

optima: an open-source R package for the Tapestri platform for integrative single cell multiomics data analysis.

SUMMARY: The Tapestri platform offers DNA and protein analysis at the single-cell level. Integrating both types of data is beneficial for studying multiple cell populations in heterogeneous microenvironments, such as tumor tissues. Here, we present optima, an R package for the processing and analysis of data generated from the Tapestri platform. This package provides streamlined functionality for raw data filtering, integration, normalization, transformation, and visualization. Insights gained from the optima package help users to identify unique cell populations and uncover surface protein expression patterns. The results generated by optima help researchers elucidate dynamic changes at the single-cell level in heterogeneous microenvironments. AVAILABILITY AND IMPLEMENTATION: This package is available in Github: https://github.com/rachelgriffard/optima.

Improving the accuracy and internal consistency of regression-based clustering of high-dimensional datasets.

Component-wise Sparse Mixture Regression (CSMR) is a recently proposed regression-based clustering method that shows promise in detecting heterogeneous relationships between molecular markers and a continuous phenotype of interest. However, CSMR can yield inconsistent results when applied to high-dimensional molecular data, which we hypothesize is in part due to inherent limitations associated with the feature selection method used in the CSMR algorithm. To assess this hypothesis, we explored whether substituting different regularized regression methods (i.e. Lasso, Elastic Net, Smoothly Clipped Absolute Deviation (SCAD), Minmax Convex Penalty (MCP), and Adaptive-Lasso) within the CSMR framework can improve the clustering accuracy and internal consistency (IC) of CSMR in high-dimensional settings. We calculated the true positive rate (TPR), true negative rate (TNR), IC and clustering accuracy of our proposed modifications, benchmarked against the existing CSMR algorithm, using an extensive set of simulation studies and real biological datasets. Our results demonstrated that substituting Adaptive-Lasso within the existing feature selection method used in CSMR led to significantly improved IC and clustering accuracy, with strong performance even in high-dimensional scenarios. In conclusion, our modifications of the CSMR method resulted in improved clustering performance and may thus serve as viable alternatives for the regression-based clustering of high-dimensional datasets.

Antibodies to oral pathobionts and colon cancer risk in the CLUE I cohort study.

Periodontitis has been associated with an increased risk for gastrointestinal cancers. The objective of our study was to investigate the association of antibodies to oral bacteria and the risk of colon cancer in a cohort setting. Using the CLUE I cohort, a prospective cohort initiated in 1974 in Washington County, Maryland, we conducted a nested case-control study to examine the association of levels of IgG antibodies to 11 oral bacterial species (13 total strains) with risk of colon cancer diagnosed a median of 16 years later (range: 1-26 years). Antibody response was measured using checkerboard immunoblotting assays. We included 200 colon cancer cases and 200 controls matched on age, sex, cigarette smoking status, time of blood draw and pipe or cigar smoking status. Controls were selected using incidence density sampling. Conditional logistic regression models were used to assess the association between antibody levels and colon cancer risk. In the overall analysis, we observed significant inverse associations for 6 of the 13 antibodies measured (P-trends <.05) and one positive association for antibody levels to Aggregatibacter actinomycetemcomitans (ATCC 29523; P-trend = .04). While we cannot rule out a role for periodontal disease in colon cancer risk, findings from our study suggest that a strong adaptive immune response may be associated with a lower risk of colon cancer. More studies will need to examine whether the positive associations we observed with antibodies to A. actinomycetemcomitans reflect a true causal association for this bacterium.

Assessing Breast Cancer Molecular Subtypes Using Extracellular Vesicles' mRNA.

Extracellular vesicles (EVs) carry RNA cargo that is believed to be associated with the cell-of-origin and thus have the potential to serve as a minimally invasive liquid biopsy marker for supplying molecular information to guide treatment decisions (i.e., precision medicine). We report the affinity isolation of EV subpopulations with monoclonal antibodies attached to the surface of a microfluidic chip that is made from a plastic to allow for high-scale production. The EV microfluidic affinity purification (EV-MAP) chip was used for the isolation of EVs sourced from two-orthogonal cell types and was demonstrated for its utility in a proof-of-concept application to provide molecular subtyping information for breast cancer patients. The orthogonal selection process better recapitulated the epithelial tumor microenvironment by isolating two subpopulations of EVs: EVEpCAM (epithelial cell adhesion molecule, epithelial origin) and EVFAPα (fibroblast activation protein α, mesenchymal origin). The EV-MAP provided recovery >80% with a specificity of 99 ± 1% based on exosomal mRNA (exo-mRNA) and real time-droplet digital polymerase chain reaction results. When selected from the plasma of healthy donors and breast cancer patients, EVs did not differ in size or total RNA mass for both markers. On average, 0.5 mL of plasma from breast cancer patients yielded ∼2.25 ng of total RNA for both EVEpCAM and EVFAPα, while in the case of cancer-free individuals, it yielded 0.8 and 1.25 ng of total RNA from EVEpCAM and EVFAPα, respectively. To assess the potential of these two EV subpopulations to provide molecular information for prognostication, we performed the PAM50 test (Prosigna) on exo-mRNA harvested from each EV subpopulation. Results suggested that EVEpCAM and EVFAPα exo-mRNA profiling using subsets of the PAM50 genes and a novel algorithm (i.e., exo-PAM50) generated 100% concordance with the tumor tissue.

Placental microRNAs relate to early childhood growth trajectories.

BACKGROUND: Poor placental function is a common cause of intrauterine growth restriction, which in turn is associated with increased risks of adverse health outcomes. Our prior work suggests that birthweight and childhood obesity-associated genetic variants functionally impact placental function and that placental microRNA are associated with birthweight. To address the influence of the placenta beyond birth, we assessed the relationship between placental microRNAs and early childhood growth. METHODS: Using the SITAR package, we generated two parameters that describe individual weight trajectories of children (0-5 years) in the New Hampshire Birth Cohort Study (NHBCS, n = 238). Using negative binomial generalized linear models, we identified placental microRNAs that relate to growth parameters (FDR < 0.1), while accounting for sex, gestational age at birth, and maternal parity. RESULTS: Genes targeted by the six growth trajectory-associated microRNAs are enriched (FDR < 0.05) in growth factor signaling (TGF/beta: miR-876; EGF/R: miR-155, Let-7c; FGF/R: miR-155; IGF/R: Let-7c, miR-155), calmodulin signaling (miR-216a), and NOTCH signaling (miR-629). CONCLUSIONS: Growth-trajectory microRNAs target pathways affecting placental proliferation, differentiation and function. Our results suggest a role for microRNAs in regulating placental cellular dynamics and supports the Developmental Origins of Health and Disease hypothesis that fetal environment can have impacts beyond birth. IMPACT: We found that growth trajectory associated placenta microRNAs target genes involved in signaling pathways central to the formation, maintenance and function of placenta; suggesting that placental cellular dynamics remain critical to infant growth to term and are under the control of microRNAs. Our results contribute to the existing body of research suggesting that the placenta plays a key role in programming health in the offspring. This is the first study to relate molecular patterns in placenta, specifically microRNAs, to early childhood growth trajectory.

Human placental microRNAs dysregulated by cadmium exposure predict neurobehavioral outcomes at birth.

BACKGROUND: Prenatal cadmium (Cd) exposure has been implicated in both placental toxicity and adverse neurobehavioral outcomes. Placental microRNAs (miRNAs) may function to developmentally program adverse pregnancy and newborn health outcomes in response to gestational Cd exposure. METHODS: In a subset of the Rhode Island Child Health Study (RICHS, n = 115) and the New Hampshire Birth Cohort Study (NHBCS, = 281), we used small RNA sequencing and trace metal analysis to identify Cd-associated expression of placental miRNAs using negative binomial generalized linear models. We predicted mRNAs targeted by Cd-associated miRNAs and relate them to neurobehavioral outcomes at birth through the integration of transcriptomic data and summary scores from the NICU Network Neurobehavioral Scale (NNNS). RESULTS: Placental Cd concentrations are significantly associated with the expression level of five placental miRNAs in NHBCS, with similar effect sizes in RICHS. These miRNA target genes overrepresented in nervous system development, and their expression is correlated with NNNS metrics suggestive of atypical neurobehavioral outcomes at birth. CONCLUSIONS: Gestational Cd exposure is associated with the expression of placental miRNAs. Predicted targets of these miRNAs are involved in nervous system development and may also regulate placental physiology, allowing their dysregulation to modify developmental programming of early life health outcomes. IMPACT: This research aims to address the poor understanding of the molecular mechanisms governing adverse pregnancy and newborn health outcomes in response to Gestational cadmium (Cd) exposure. Our results outline a robust relationship between Cd-associated placental microRNA expression and NICU Network Neurobehavioral Scales (NNNS) at birth indicative of atypical neurobehavior. This study utilized healthy mother-infant cohorts to describe the role of Cd-associated dysregulation of placental microRNAs as a potential mechanism by which adverse neurobehavioral outcomes are developmentally programmed.

pwrBRIDGE: a user-friendly web application for power and sample size estimation in batch-confounded microarray studies with dependent samples.

Batch effect Reduction of mIcroarray data with Dependent samples usinG Empirical Bayes (BRIDGE) is a recently developed statistical method to address the issue of batch effect correction in batch-confounded microarray studies with dependent samples. The key component of the BRIDGE methodology is the use of samples run as technical replicates in two or more batches, "bridging samples", to inform batch effect correction/attenuation. While previously published results indicate a relationship between the number of bridging samples, M, and the statistical power of downstream statistical testing on the batch-corrected data, there is of yet no formal statistical framework or user-friendly software, for estimating M to achieve a specific statistical power for hypothesis tests conducted on the batch-corrected data. To fill this gap, we developed pwrBRIDGE, a simulation-based approach to estimate the bridging sample size, M, in batch-confounded longitudinal microarray studies. To illustrate the use of pwrBRIDGE, we consider a hypothetical, longitudinal batch-confounded study whose goal is to identify Alzheimer's disease (AD) progression-associated genes from amnestic mild cognitive impairment (aMCI) to AD in human blood after a 5-year follow-up. pwrBRIDGE helps researchers design and plan batch-confounded microarray studies with dependent samples to avoid over- or under-powered studies.

Epigenome-wide association study using peripheral blood leukocytes identifies genomic regions associated with periodontal disease and edentulism in the Atherosclerosis Risk in Communities study.

AIM: To investigate individual susceptibility to periodontitis by conducting an epigenome-wide association study using peripheral blood. MATERIALS AND METHODS: We included 1077 African American and 457 European American participants of the Atherosclerosis Risk in Communities (ARIC) study who had completed a dental examination or reported being edentulous at Visit 4 and had available data on DNA methylation from Visit 2 or 3. DNA methylation levels were compared by periodontal disease severity and edentulism through discovery analyses and subsequent testing of individual CpGs. RESULTS: Our discovery analysis replicated findings from a previous study reporting a region in gene ZFP57 (6p22.1) that was significantly hypomethylated in severe periodontal disease compared with no/mild periodontal disease in European American participants. Higher methylation levels in a separate region in an unknown gene (located in Chr10: 743,992-744,958) was associated with significantly higher odds of edentulism compared with no/mild periodontal disease in African American participants. In subsequent CpG testing, four CpGs in a region previously associated with periodontitis located within HOXA4 were significantly hypermethylated in severe periodontal disease compared with no/mild periodontal disease in African American participants (odds ratio per 1 SD increase in methylation level: cg11015251: 1.28 (1.02, 1.61); cg14359292: 1.24 (1.01, 1.54); cg07317062: 1.30 (1.05, 1.61); cg08657492: 1.25 (1.01, 1.55)). CONCLUSIONS: Our study highlights epigenetic variations in ZPF57 and HOXA4 that are significantly and reproducibly associated with periodontitis. Future studies should evaluate gene regulatory mechanisms in the candidate regions of these loci.

HiTIMED: hierarchical tumor immune microenvironment epigenetic deconvolution for accurate cell type resolution in the tumor microenvironment using tumor-type-specific DNA methylation data.

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.

Synergistic anti-proliferative activity of JQ1 and GSK2801 in triple-negative breast cancer.

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.

Batch effect reduction of microarray data with dependent samples using an empirical Bayes approach (BRIDGE).

Batch-effects present challenges in the analysis of high-throughput molecular data and are particularly problematic in longitudinal studies when interest lies in identifying genes/features whose expression changes over time, but time is confounded with batch. While many methods to correct for batch-effects exist, most assume independence across samples; an assumption that is unlikely to hold in longitudinal microarray studies. We propose Batch effect Reduction of mIcroarray data with Dependent samples usinGEmpirical Bayes (BRIDGE), a three-step parametric empirical Bayes approach that leverages technical replicate samples profiled at multiple timepoints/batches, so-called "bridge samples", to inform batch-effect reduction/attenuation in longitudinal microarray studies. Extensive simulation studies and an analysis of a real biological data set were conducted to benchmark the performance of BRIDGE against both ComBat and longitudinalComBat. Our results demonstrate that while all methods perform well in facilitating accurate estimates of time effects, BRIDGE outperforms both ComBat and longitudinal ComBat in the removal of batch-effects in data sets with bridging samples, and perhaps as a result, was observed to have improved statistical power for detecting genes with a time effect. BRIDGE demonstrated competitive performance in batch effect reduction of confounded longitudinal microarray studies, both in simulated and a real data sets, and may serve as a useful preprocessing method for researchers conducting longitudinal microarray studies that include bridging samples.

Low variability in the underlying cellular landscape adversely affects the performance of interaction-based approaches for conducting cell-specific analyses of DNA methylation in bulk samples.

Statistical methods that allow for cell type specific DNA methylation (DNAm) analyses based on bulk-tissue methylation data have great potential to improve our understanding of human disease and have created unprecedented opportunities for new insights using the wealth of publicly available bulk-tissue methylation data. These methodologies involve incorporating interaction terms formed between the phenotypes/exposures of interest and proportions of the cell types underlying the bulk-tissue sample used for DNAm profiling. Despite growing interest in such "interaction-based" methods, there has been no comprehensive assessment how variability in the cellular landscape across study samples affects their performance. To answer this question, we used numerous publicly available whole-blood DNAm data sets along with extensive simulation studies and evaluated the performance of interaction-based approaches in detecting cell-specific methylation effects. Our results show that low cell proportion variability results in large estimation error and low statistical power for detecting cell-specific effects of DNAm. Further, we identified that many studies targeting methylation profiling in whole-blood may be at risk to be underpowered due to low variability in the cellular landscape across study samples. Finally, we discuss guidelines for researchers seeking to conduct studies utilizing interaction-based approaches to help ensure that their studies are adequately powered.

A ketogenic diet differentially affects neuron and astrocyte transcription.

Ketogenic diets (KDs) alter brain metabolism. Multiple mechanisms may account for their effects, and different brain regions may variably respond. Here, we considered how a KD affects brain neuron and astrocyte transcription. We placed male C57Bl6/N mice on either a 3-month KD or chow diet, generated enriched neuron and astrocyte fractions, and used RNA-Seq to assess transcription. Neurons from KD-treated mice generally showed transcriptional pathway activation while their astrocytes showed a mix of transcriptional pathway suppression and activation. The KD especially affected pathways implicated in mitochondrial and endoplasmic reticulum function, insulin signaling, and inflammation. An unbiased analysis of KD-associated expression changes strongly implicated transcriptional pathways altered in AD, which prompted us to explore in more detail the potential molecular relevance of a KD to AD. Our results indicate a KD differently affects neurons and astrocytes, and provide unbiased evidence that KD-induced brain effects are potentially relevant to neurodegenerative diseases such as AD.

Tracing human stem cell lineage during development using DNA methylation.

Stem cell maturation is a fundamental, yet poorly understood aspect of human development. We devised a DNA methylation signature deeply reminiscent of embryonic stem cells (a fetal cell origin signature, FCO) to interrogate the evolving character of multiple human tissues. The cell fraction displaying this FCO signature was highly dependent upon developmental stage (fetal versus adult), and in leukocytes, it described a dynamic transition during the first 5 yr of life. Significant individual variation in the FCO signature of leukocytes was evident at birth, in childhood, and throughout adult life. The genes characterizing the signature included transcription factors and proteins intimately involved in embryonic development. We defined and applied a DNA methylation signature common among human fetal hematopoietic progenitor cells and have shown that this signature traces the lineage of cells and informs the study of stem cell heterogeneity in humans under homeostatic conditions.

Pre-diagnosis neutrophil-to-lymphocyte ratio and mortality in individuals who develop lung cancer.

PURPOSE: The neutrophil-to-lymphocyte ratio (NLR) is a marker of systemic inflammation that has been reported to be associated with survival after chronic disease diagnoses, including lung cancer. We hypothesized that the inflammatory profile reflected by pre-diagnosis NLR, rather than the well-studied pre-treatment NLR at diagnosis, may be associated with increased mortality after lung cancer is diagnosed in high-risk heavy smokers. METHODS: We examined associations between pre-diagnosis methylation-derived NLR (mdNLR) and lung cancer-specific and all-cause mortality in 279 non-small lung cancer (NSCLC) and 81 small cell lung cancer (SCLC) cases from the ß-Carotene and Retinol Efficacy Trial (CARET). Cox proportional hazards models were adjusted for age, sex, smoking status, pack years, and time between blood draw and diagnosis, and stratified by stage of disease. Models were run separately by histotype. RESULTS: Among SCLC cases, those with pre-diagnosis mdNLR in the highest quartile had 2.5-fold increased mortality compared to those in the lowest quartile. For each unit increase in pre-diagnosis mdNLR, we observed 22-23% increased mortality (SCLC-specific hazard ratio [HR] = 1.23, 95% confidence interval [CI]: 1.02, 1.48; all-cause HR = 1.22, 95% CI 1.01, 1.46). SCLC associations were strongest for current smokers at blood draw (Interaction Ps = 0.03). Increasing mdNLR was not associated with mortality among NSCLC overall, nor within adenocarcinoma (N = 148) or squamous cell carcinoma (N = 115) case groups. CONCLUSION: Our findings suggest that increased mdNLR, representing a systemic inflammatory profile on average 4.5 years before a SCLC diagnosis, may be associated with mortality in heavy smokers who go on to develop SCLC but not NSCLC.

Leveraging cell-specific differentially methylated regions to identify leukocyte infiltration in adipose tissue.

Obesity is understood to be an inflammatory condition characterized in part by changes in resident immune cell populations in adipose tissue. However, much of this knowledge has been obtained through experimental animal models. Epigenetic mechanisms, such as DNA methylation may be useful tools for characterizing the changes in immune cell populations in human subjects. In this study, we introduce a simple and intuitive method for assessing cellular infiltration by blood into other heterogeneous, admixed tissues such as adipose tissue, and apply this approach in a large human cohort study. Associations between higher leukocyte infiltration, measured by evaluating a distance measure between the methylation signatures of leukocytes and adipose tissue, and increasing body mass index (BMI) or android fat mass (AFM) were identified and validated in independent replication samples for CD4 (pBMI = 0.009, pAFM = 0.020), monocytes (pBMI = 0.001, pAFM = 4.3 × 10-4 ), and dendritic cells (pBMI = 0.571, pAFM = 0.012). Patterns of depletion with increasing adiposity were observed for plasma B (pBMI = 0.430, pAFM = 0.004) and immature B (pBMI = 0.022, pAFM = 0.042) cells. CD4, dendritic, monocytes, immature B, and plasma B cells may be important agents in the inflammatory process. Finally, the method used to assess leukocyte infiltration in this study is straightforwardly extended to other cell types and tissues in which infiltration might be of interest.

O-GlcNAc homeostasis contributes to cell fate decisions during hematopoiesis.

The addition of a single ß-d-GlcNAc sugar (O-GlcNAc) by O-GlcNAc-transferase (OGT) and O-GlcNAc removal by O-GlcNAcase (OGA) maintain homeostatic O-GlcNAc levels on cellular proteins. Changes in protein O-GlcNAcylation regulate cellular differentiation and cell fate decisions, but how these changes affect erythropoiesis, an essential process in blood cell formation, remains unclear. Here, we investigated the role of O-GlcNAcylation in erythropoiesis by using G1E-ER4 cells, which carry the erythroid-specific transcription factor GATA-binding protein 1 (GATA-1) fused to the estrogen receptor (GATA-1-ER) and therefore undergo erythropoiesis after ß-estradiol (E2) addition. We observed that during G1E-ER4 differentiation, overall O-GlcNAc levels decrease, and physical interactions of GATA-1 with both OGT and OGA increase. RNA-Seq-based transcriptome analysis of G1E-ER4 cells differentiated in the presence of the OGA inhibitor Thiamet-G (TMG) revealed changes in expression of 433 GATA-1 target genes. ChIP results indicated that the TMG treatment decreases the occupancy of GATA-1, OGT, and OGA at the GATA-binding site of the lysosomal protein transmembrane 5 (Laptm5) gene promoter. TMG also reduced the expression of genes involved in differentiation of NB4 and HL60 human myeloid leukemia cells, suggesting that O-GlcNAcylation is involved in the regulation of hematopoietic differentiation. Sustained treatment of G1E-ER4 cells with TMG before differentiation reduced hemoglobin-positive cells and increased stem/progenitor cell surface markers. Our results show that alterations in O-GlcNAcylation disrupt transcriptional programs controlling erythropoietic lineage commitment, suggesting a role for O-GlcNAcylation in regulating hematopoietic cell fate.

Equivalent change enrichment analysis: assessing equivalent and inverse change in biological pathways between diverse experiments.

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.

AHRR methylation in heavy smokers: associations with smoking, lung cancer risk, and lung cancer mortality.

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.