ABSTRACT
DNA methylation (DNAm)-based deconvolution estimates contain relative data, forming a composition, that standard methods (testing directly on cell proportions) are ill-suited to handle. In this study we examined the performance of an alternative method, analysis of compositions of microbiomes (ANCOM), for the analysis of DNAm-based deconvolution estimates. We performed two different simulation studies comparing ANCOM to a standard approach (two sample t-test performed directly on cell proportions) and analyzed a real-world data from the Women's Health Initiative to evaluate the applicability of ANCOM to DNAm-based deconvolution estimates. Our findings indicate that ANCOM can effectively account for the compositional nature of DNAm-based deconvolution estimates. ANCOM adequately controls the false discovery rate while maintaining statistical power comparable to that of standard methods.
DNA methylation (DNAm)-based deconvolution provides highly accurate estimates of the proportion of each cell type in a mixed-cell type biological sample (e.g., whole-blood). These estimates can be used for examining the association between cell type proportions and biological or clinical end points; for example, comparing the estimated neutrophil proportion in whole blood between smokers and non-smokers. Cell proportion data has unique features which present challenges for traditional and widely used statistical methods. In response to this issue, our work presents two simulation studies and a real-world analysis that benchmark the performance of current standard statistical methods against an alternative method called analysis composition of microbes (ANCOM), which was originally developed for the analysis of microbiome data. In our real-world analysis we used DNAm data collected from Women's Health Initiative Long Life Study I and compared the results of each method against a gold-standard that is typically not available for these analyses. In each of our simulation studies, ANCOM was able to detect true differences in cell proportions between the groups being compared but had a much lower rate of false discovery compared with the standard statistical methods. Our real-world analysis demonstrated similar findings. Overall, our study highlights the potential of ANCOM as a powerful and robust method for analyzing DNAm-derived deconvolution estimates when the interest is comparisons of cell type proportions and biological or clinical end points. ANCOM's ability to minimize false discovery while maintaining robust statistical power positions it as a valuable addition to the epigenomic analysis toolkit.
ABSTRACT
Aim: This study addresses the challenge of predicting the response of head and neck squamous cell carcinoma (HNSCC) patients to immunotherapy. Methods: Using DNA methylation cytometry, we analyzed the immune profiles of six HNSCC patients who showed a positive response to immunotherapy over a year without disease progression. Results: There was an initial increase in CD8 T memory cells and natural killer cells during the first four cycles of immunotherapy, which then returned to baseline levels after a year. Baseline CD8 T cell levels were lower in HNSCC immunotherapy responders but became similar to those in healthy subjects after immunotherapy. Conclusion: These findings suggest that monitoring fluctuations in immune profiles could potentially identify biomarkers for immunotherapy response in HNSCC patients.
[Box: see text].
ABSTRACT
Background: Idiopathic interstitial pneumonias (IIPs) such as idiopathic pulmonary fibrosis (IPF) and interstitial pneumonia with autoimmune features (IPAF), present diagnostic and therapeutic challenges due to their heterogeneous nature. This study aimed to identify intrinsic molecular signatures within the lung microenvironment of these IIPs through proteomic analysis of bronchoalveolar lavage fluid (BALF). Methods: Patients with IIP (n=23) underwent comprehensive clinical evaluation including pre-treatment bronchoscopy and were compared to controls without lung disease (n=5). Proteomic profiling of BALF was conducted using label-free quantitative methods. Unsupervised cluster analyses identified protein expression profiles which were then analyzed to predict survival outcomes and investigate associated pathways. Results: Proteomic profiling successfully differentiated IIP from controls. k-means clustering, based on protein expression revealed three distinct IIP clusters, which were not associated with age, smoking history, or baseline pulmonary function. These clusters had unique survival trajectories and provided more accurate survival predictions than the Gender Age Physiology (GAP) index (C-index 0.794 vs. 0.709). The cluster with the worst prognosis featured decreased inflammatory signaling and complement activation, with pathway analysis highlighting altered immune response pathways related to immunoglobulin production and B cell-mediated immunity. Conclusions: The unsupervised clustering of BALF proteomics provided a novel stratification of IIP patients, with potential implications for prognostic and therapeutic targeting. The identified molecular phenotypes underscore the diversity within the IIP classification and the potential importance of personalized treatments for these conditions. Future validation in larger, multi-ethnic cohorts is essential to confirm these findings and to explore their utility in clinical decision-making for patients with IIP.
ABSTRACT
Prior cohort studies assessing cancer risk based on immune cell subtype profiles have predominantly focused on White populations. This limitation obscures vital insights into how cancer risk varies across race. Immune cell subtype proportions were estimated using deconvolution based on leukocyte DNA methylation markers from blood samples collected at baseline on participants without cancer in the Atherosclerosis Risk in Communities (ARIC) Study. Over a mean of 17.5 years of follow-up, 668 incident cancers were diagnosed in 2,467 Black participants. Cox proportional hazards regression was used to examine immune cell subtype proportions and overall cancer incidence and site-specific incidence (lung, breast, and prostate cancers). Higher T regulatory cell proportions were associated with statistically significantly higher lung cancer risk (hazard ratio = 1.22, 95% confidence interval = 1.06-1.41 per percent increase). Increased memory B cell proportions were associated with significantly higher risk of prostate cancer (1.17, 1.04-1.33) and all cancers (1.13, 1.05-1.22). Increased CD8+ naïve cell proportions were associated with significantly lower risk of all cancers in participants ≥55 years (0.91, 0.83-0.98). Other immune cell subtypes did not display statistically significant associations with cancer risk. These results in Black participants align closely with prior findings in largely White populations. Findings from this study could help identify those at high cancer risk and outline risk stratifying to target patients for cancer screening, prevention, and other interventions. Further studies should assess these relationships in other cancer types, better elucidate the interplay of B cells in cancer risk, and identify biomarkers for personalized risk stratification.
ABSTRACT
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.
ABSTRACT
Background: Bladder cancer and therapy responses hinge on immune profiles in the tumor microenvironment (TME) and blood, yet studies linking tumor-infiltrating immune cells to peripheral immune profiles are limited. Methods: DNA methylation cytometry quantified TME and matched peripheral blood immune cell proportions. With tumor immune profile data as the input, subjects were grouped by immune infiltration status and consensus clustering. Results: Immune hot and cold groups had different immune compositions in the TME but not in circulating blood. Two clusters of patients identified with consensus clustering had different immune compositions not only in the TME but also in blood. Conclusion: Detailed immune profiling via methylation cytometry reveals the significance of understanding tumor and systemic immune relationships in cancer patients.
Bladder cancer and treatment outcomes depend on the immune profiles in the tumor and blood. Our study, using DNA methylation cytometry, measured immune cell proportions in both areas. Patients were grouped based on immune status and consensus clustering. Results showed distinct immune compositions in the tumor, but not in blood, for hot and cold groups. Consensus clustering revealed two patient clusters with differing immune compositions in both tumor and blood. This detailed immune profiling highlights the importance of understanding the complex interplay between tumor and systemic immunity in bladder cancer patients.
Subject(s)
Tumor Microenvironment , Urinary Bladder Neoplasms , Humans , Urinary Bladder Neoplasms/genetics , Cluster Analysis , DNA Methylation , Protein Processing, Post-Translational , PrognosisABSTRACT
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.
Subject(s)
Multiomics , Software , Data AnalysisABSTRACT
BACKGROUND: Immune profiles have been associated with bladder cancer outcomes and may have clinical applications for prognosis. However, associations of detailed immune cell subtypes with patient outcomes remain underexplored and may contribute crucial prognostic information for better managing bladder cancer recurrence and survival. METHODS: Bladder cancer case peripheral blood DNA methylation was measured using the Illumina HumanMethylationEPIC array. Extended cell-type deconvolution quantified 12 immune cell-type proportions, including memory, naïve T and B cells, and granulocyte subtypes. DNA methylation clocks determined biological age. Cox proportional hazards models tested associations of immune cell profiles and age acceleration with bladder cancer outcomes. The partDSA algorithm discriminated 10-year overall survival groups from clinical variables and immune cell profiles, and a semi-supervised recursively partitioned mixture model (SS-RPMM) with DNA methylation data was applied to identify a classifier for 10-year overall survival. RESULTS: Higher CD8T memory cell proportions were associated with better overall survival [HR = 0.95, 95% confidence interval (CI) = 0.93-0.98], while higher neutrophil-to-lymphocyte ratio (HR = 1.36, 95% CI = 1.23-1.50), CD8T naïve (HR = 1.21, 95% CI = 1.04-1.41), neutrophil (HR = 1.04, 95% CI = 1.03-1.06) proportions, and age acceleration (HR = 1.06, 95% CI = 1.03-1.08) were associated with worse overall survival in patient with bladder cancer. partDSA and SS-RPMM classified five groups of subjects with significant differences in overall survival. CONCLUSIONS: We identified associations between immune cell subtypes and age acceleration with bladder cancer outcomes. IMPACT: The findings of this study suggest that bladder cancer outcomes are associated with specific methylation-derived immune cell-type proportions and age acceleration, and these factors could be potential prognostic biomarkers.
Subject(s)
Neoplasm Recurrence, Local , Urinary Bladder Neoplasms , Humans , Neoplasm Recurrence, Local/genetics , DNA Methylation , Lymphocytes , Proportional Hazards Models , PrognosisABSTRACT
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.
Subject(s)
Atherosclerosis , Periodontal Diseases , Periodontitis , Humans , Epigenome , Genome-Wide Association Study , Periodontal Diseases/genetics , Atherosclerosis/genetics , Periodontitis/genetics , Leukocytes , GenomicsABSTRACT
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.
Subject(s)
Algorithms , Benchmarking , Cluster Analysis , Computer Simulation , PhenotypeABSTRACT
Introduction: The 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 deconvolved. Methods: Using 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. Results: We 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. Discussion: We 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
DNA methylation (DNAm)-based cell mixture deconvolution (CMD) has become a quintessential part of epigenome-wide association studies where DNAm is profiled in heterogeneous tissue types. Despite being introduced over a decade ago, detection limits, which represent the smallest fraction of a cell type in a mixed biospecimen that can be reliably detected, have yet to be determined in the context of DNAm-based CMD. Moreover, there has been little attention given to approaches for quantifying the uncertainty associated with DNAm-based CMD. Here, analytical frameworks for determining both cell-specific limits of detection and quantification of uncertainty associated with DNAm-based CMD are described. This work may contribute to improved rigor, reproducibility and replicability of epigenome-wide association studies involving CMD.
Subject(s)
DNA Methylation , Epigenesis, Genetic , Humans , Uncertainty , Limit of Detection , Reproducibility of ResultsABSTRACT
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.
Subject(s)
Breast Neoplasms , Extracellular Vesicles , Humans , Female , Breast Neoplasms/genetics , Breast Neoplasms/metabolism , Epithelial Cell Adhesion Molecule/metabolism , RNA, Messenger/genetics , RNA, Messenger/metabolism , Extracellular Vesicles/metabolism , Liquid Biopsy , Tumor MicroenvironmentABSTRACT
BACKGROUND: Oral health is a key indicator of overall health, well-being, and quality of life. Several studies have provided new evidence about the role of oral diseases, specifically periodontitis, in generating risk for various forms of cancers, including lung, colorectal, and pancreatic cancers. METHODS: Incident lung cancer cases (n = 192) and matched controls (n = 192) were selected from participants of the CLUE I and CLUE II cohorts. Archived serum samples collected from participants in 1974 (in CLUE I) were analyzed using immunoblotting for immunoglobulin G (IgG) antibody levels to 13 bacteria of the periodontium. Associations between antibody levels and lung cancer were estimated using conditional logistic regression. RESULTS: Most of the periodontal bacterial antibodies measured were inversely associated with lung cancer risk; of these, 3 were statistically significant (Prevotellaintermedia, Actinomyces naeslundii, and Veillonella parvula). A statistically significant positive association was observed for one of the Porphyromonas gingivalis strains after adjusting for P. intermedia. The sum of the logarithm of antibodies against the 13 measured bacteria was inversely associated with risk of lung cancer when the analysis was restricted to a longer follow-up (31-44 years after blood collection, highest vs lowest quartile: odds ratio = 0.26, 95% confidence interval = 0.08 to 0.84). CONCLUSIONS: Findings from this study highlight the complexity of using serum IgG antibodies to periodontal bacteria to identify associations between oral pathogens and risk of lung cancer. The inverse associations observed for antibodies to periodontal bacteria suggest that these may represent markers of immunity that provide some advantage in reducing the development of lung cancer.
Subject(s)
Lung Neoplasms , Quality of Life , Humans , Immunoglobulin G , Porphyromonas gingivalis , Lung Neoplasms/epidemiology , LungABSTRACT
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.
Subject(s)
DNA Methylation , Postmenopause , Humans , Female , Middle Aged , Aged , Aged, 80 and over , Leukocytes , CD8-Positive T-Lymphocytes , Women's HealthABSTRACT
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.
Subject(s)
Antibodies, Bacterial , Colonic Neoplasms , Humans , Cohort Studies , Case-Control Studies , Prospective Studies , Bacteria , Colonic Neoplasms/epidemiology , Colonic Neoplasms/etiologyABSTRACT
BACKGROUND: Epigenetic age, a robust marker of biological aging, has been associated with obesity, low-grade inflammation and metabolic diseases. However, few studies have examined associations between different epigenetic age measures and risk of lung cancer, despite great interest in finding biomarkers to assist in risk stratification for lung cancer screening. METHODS: A nested case-control study of lung cancer from the CLUE II cohort study was conducted using incidence density sampling with 1:1 matching of controls to lung cancer cases (n = 208 matched pairs). Prediagnostic blood samples were collected in 1989 (CLUE II study baseline) and stored at -70°C. DNA was extracted from buffy coat and DNA methylation levels were measured using Illumina MethylationEPIC BeadChip Arrays. Three epigenetic age acceleration (i.e., biological age is greater than chronological age) measurements (Horvath, Hannum and PhenoAge) were examined in relation to lung cancer risk using conditional logistic regression. RESULTS: We did not observe associations between the three epigenetic age acceleration measurements and risk of lung cancer overall; however, inverse associations for the two Hannum age acceleration measures (intrinsic and extrinsic) were observed in men and among younger participants, but not in women or older participants. We did not observe effect modification by time from blood draw to diagnosis. CONCLUSION: Findings from this study do not support a positive association between three different biological age acceleration measures and risk of lung cancer. Additional studies are needed to address whether epigenetic age is associated with lung cancer in never smokers.
Subject(s)
Epigenesis, Genetic , Lung Neoplasms , Male , Humans , Female , Cohort Studies , Case-Control Studies , Prospective Studies , Early Detection of Cancer , Lung Neoplasms/epidemiology , Lung Neoplasms/genetics , Aging/genetics , DNA MethylationABSTRACT
In the United States, cardiovascular disease is the leading cause of death and the rate of maternal mortality remains among the highest of any industrialized nation. Maternal cardiometabolic health throughout gestation and postpartum is representative of placental health and physiology. Both proper placental functionality and placental microRNA expression are essential to successful pregnancy outcomes, and both are highly sensitive to genetic and environmental sources of variation. Placental pathologies, such as preeclampsia, are associated with maternal cardiovascular health but may also contribute to the developmental programming of chronic disease in offspring. However, the role of more subtle alterations to placental function and microRNA expression in this developmental programming remains poorly understood. We performed small RNA sequencing to investigate microRNA in placentae from the Rhode Island Child Health Study (n = 230). MicroRNA counts were modeled on maternal family history of cardiovascular disease using negative binomial generalized linear models. MicroRNAs were considered to be differentially expressed at a false discovery rate (FDR) less than 0.10. Parallel mRNA sequencing data and bioinformatic target prediction software were then used to identify potential mRNA targets of differentially expressed microRNAs. Nine differentially expressed microRNAs were identified (FDR < 0.1). Bioinformatic target prediction revealed 66 potential mRNA targets of these microRNAs, many of which are implicated in TGFß signaling pathway but also in pathways involving cellular metabolism and immunomodulation. A robust association exists between familial cardiovascular disease and placental microRNA expression which may be implicated in both placental insufficiencies and the developmental programming of chronic disease.
Subject(s)
Cardiovascular Diseases , MicroRNAs , Placenta , Female , Humans , Pregnancy , Cardiovascular Diseases/genetics , Cardiovascular Diseases/metabolism , MicroRNAs/genetics , MicroRNAs/metabolism , Placenta/metabolism , Pregnancy Outcome , RNA, Messenger/metabolismABSTRACT
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.
Subject(s)
MicroRNAs , Placenta , Infant, Newborn , Child , Humans , Pregnancy , Female , Placenta/metabolism , MicroRNAs/genetics , MicroRNAs/metabolism , Cadmium , Cohort Studies , ParturitionABSTRACT
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.