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.

Human cytomegalovirus alters immune cell profile with potential implications for patient survival in head and neck cancer.

Cytomegalovirus (CMV) is a highly prevalent human herpes virus that exerts a strong influence on immune repertoire which may influence cancer risk. We have tested whether CMV immunoglobulin G (IgG) serostatus is associated with immune cell proportions (n = 132 population controls), human papillomavirus (HPV) co-infection and head and neck cancer risk (n = 184 cancer cases and 188 controls) and patient survival. CMV status was not associated with the proportion of Natural Killer cells, B cells or the neutrophil-to-lymphocyte ratio. However, CD8+ T cells increased with increasing categories of IgG titers (P =1.7 × 10-10), and titers were inversely associated with the CD4:CD8 ratio (P = 5.6 × 10-5). Despite these differences in T cell proportions, CMV was not associated with HPV16 co-infection. CMV seropositivity was similar in cases (52%) and controls (47%) and was not associated with patient survival (hazard ratio [HR] 1.14, 95% confidence interval [CI]: 0.70 to 1.86). However, those patients with the highest titers had the worst survival (HR 1.91, 95% CI: 1.13 to 3.23). Tumor-based data from The Cancer Genome Atlas demonstrated that the presence of CMV transcripts was associated with worse patient survival (HR 1.79, 95% CI: 0.96 to 2.78). These findings confirm that a history of CMV infection alters T cell proportions, but this does not translate to HPV16 co-infection or head and neck cancer risk. Our data suggest that high titers and active CMV virus in the tumor environment may confer worse survival.

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.

Patients with type 2 diabetes and elevated fibrosis-4 are under-referred to hepatology and have unrecognized hepatic decompensation.

BACKGROUND AND AIM: The American Association for the Study of Liver Diseases recommends a high index of suspicion for nonalcoholic steatohepatitis and advanced fibrosis in patients with type 2 diabetes (T2D) and an elevated fibrosis-4 index (FIB-4). We investigated the referral pattern of patients with T2D and FIB4 > 3.25 to the hepatology clinic and evaluated the clinical benefits to the patient. METHODS: We included patients aged 18-80 years with T2D and a FIB4 score >3.25 who had visited the internal medicine, family medicine, endocrinology clinic from 01/01/2014-5/31/2019. The first time point of high-risk FIB-4 was identified as the baseline for time-to-event analysis. The patients were classified based on whether they had visited the hepatology clinic (referred vs not referred). RESULTS: Of the 2174 patients, 290 (13.3%) were referred to the hepatology clinic, and 1884 (86.7%) were not referred. In multivariate analyses, the referred patients had a lower overall mortality risk (Hazard Ratio: 0.57; 95% CI: 0.38-87). Notably, the referred patients had the same rate of biochemical decompensation, as measured by progression to MELD ≥ 14, but a substantially higher rate of diagnosis in cirrhosis (27, 19-38) and cirrhosis complications, including ascites (2.9, 2.0-4.1), hepatic encephalopathy (99, 13-742), and liver cancer (14, 5-38). CONCLUSIONS: We found that patients with T2D and high-risk FIB4 are associated with better overall survival after referral to a hepatology clinic. We speculate that the survival difference is due to the increased recognition of cirrhosis and cirrhosis complications in the referred populations.

Antifungal drug price increases in the United States, 2000-2019.

BACKGROUND: Antifungal drugs treat a variety of conditions, ranging from localised dermatologic disease to life-threatening systemic infections. Some common antifungal drugs experienced large price increases in recent years, however, factors contributing to these price increases are poorly understood. We sought to examine trends in antifungal drug prices and determine underlying drivers of price changes. METHODS: Antifungal drug products in the United States were identified using the Food and Drug Administration (FDA) Label database. For each product, we determined the wholesale acquisition cost per unit over time between 2000 and 2019, adjusting for inflation, and examined variables that could impact price: route of administration, number of FDA indications, the quantity of professional guideline recommendations, use for prophylaxis, number of FDA-approved manufacturers, and whether it was compounded. Price trajectories were clustered into four groups: (1) stable, 2) moderate, (3) high, and (4) extreme price increases. RESULTS: Of 139 identified drug products, one outlier was removed due to exorbitant price increases. Cluster 1 (n = 31) demonstrated the most stable prices with a 25% mean price increase. Clusters 2 (n = 97), 3 (n = 7), and 4 (n = 3) demonstrated moderate, high, and extreme price increases with 52%, 318%, and 900% mean price increases, respectively. Atypical routes of administration and compounding were over-represented in clusters 3 and 4. There was no correlation between the number of manufacturers and price changes. CONCLUSIONS: Antifungal drugs exhibited large, inflation-adjusted price increases. Atypical routes of administration and compounding were over-represented within clusters exhibiting extraordinary price increases. Our data support policies aiming to curb large price increases for medically important drugs.

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.