Accelerated epigenetic age is associated with whole-brain functional connectivity and impaired cognitive performance in older adults

While chronological age is a strong predictor for health-related risk factors, it is an incomplete metric that fails to fully characterize the unique aging process of individuals with different genetic makeup, neurodevelopment, and environmental experiences. Recent advances in epigenomic array technologies have made it possible to generate DNA methylation-based biomarkers of biological aging, which may be useful in predicting a myriad of cognitive abilities and functional brain network organization across older individuals. It is currently unclear which cognitive domains are negatively correlated with epigenetic age above and beyond chronological age, and it is unknown if functional brain organization is an important mechanism for explaining these associations. In this study, individuals with accelerated epigenetic age (i.e. AgeAccelGrim) performed worse on tasks that spanned a wide variety of cognitive faculties including both fluid and crystallized intelligence (N = 103, average age = 68.98 years, 73 females, 30 males). Additionally, fMRI connectome-based predictive models suggested a mediating mechanism of functional connectivity on epigenetic age acceleration-cognition associations primarily in medial temporal lobe and limbic structures. This research highlights the important role of epigenetic aging processes on the development and maintenance of healthy cognitive capacities and function of the aging brain.

Spatial transcriptomics reveals segregation of tumor cell states in glioblastoma and marked immunosuppression within the perinecrotic niche

Glioblastoma (GBM) remains an untreatable malignant tumor with poor patient outcomes, characterized by palisading necrosis and microvascular proliferation. While single-cell technology made it possible to characterize different lineage of glioma cells into neural progenitor-like (NPC-like), oligodendrocyte-progenitor-like (OPC-like), astrocyte-like (AC-like) and mesenchymal like (MES-like) states, it does not capture the spatial localization of these tumor cell states. Spatial transcriptomics empowers the study of the spatial organization of different cell types and tumor cell states and allows for the selection of regions of interest to investigate region-specific and cell-type-specific pathways. Here, we obtained paired 10x Chromium single-nuclei RNA-sequencing (snRNA-seq) and 10x Visium spatial transcriptomics data from three GBM patients to interrogate the GBM microenvironment. Integration of the snRNA-seq and spatial transcriptomics data reveals patterns of segregation of tumor cell states. For instance, OPC-like tumor and NPC-like tumor significantly segregate in two of the three samples. Our differentially expressed gene and pathway analyses uncovered significant pathways in functionally relevant niches. Specifically, perinecrotic regions were more immunosuppressive than the endogenous GBM microenvironment, and perivascular regions were more pro-inflammatory. Our gradient analysis suggests that OPC-like tumor cells tend to reside in areas closer to the tumor vasculature compared to tumor necrosis, which may reflect increased oxygen requirements for OPC-like cells. In summary, we characterized the localization of cell types and tumor cell states, the gene expression patterns, and pathways in different niches within the GBM microenvironment. Our results provide further evidence of the segregation of tumor cell states and highlight the immunosuppressive nature of the necrotic and perinecrotic niches in GBM.

Subventricular zone stem cell niche injury is associated with intestinal perforation in preterm infants and predicts future motor impairment

Brain injury is highly associated with preterm birth. Complications of prematurity, including spontaneous or necrotizing enterocolitis (NEC)-associated intestinal perforations, are linked to lifelong neurologic impairment, yet the mechanisms are poorly understood. Early diagnosis of preterm brain injuries remains a significant challenge. Here, we identified subventricular zone echogenicity (SVE) on cranial ultrasound in preterm infants following intestinal perforations. The development of SVE was significantly associated with motor impairment at 2 years. SVE was replicated in a neonatal mouse model of intestinal perforation. Examination of the murine echogenic subventricular zone (SVZ) revealed NLRP3-inflammasome assembly in multiciliated FoxJ1+ ependymal cells and a loss of the ependymal border in this postnatal stem cell niche. These data suggest a mechanism of preterm brain injury localized to the SVZ that has not been adequately considered. Ultrasound detection of SVE may serve as an early biomarker for neurodevelopmental impairment after inflammatory disease in preterm infants.

PI3K/mTOR is a therapeutically targetable genetic dependency in diffuse intrinsic pontine glioma

Diffuse midline glioma (DMG), including tumors diagnosed in the brainstem (diffuse intrinsic pontine glioma; DIPG), are uniformly fatal brain tumors that lack effective treatment. Analysis of CRISPR/Cas9 loss-of-function gene deletion screens identified PIK3CA and MTOR as targetable molecular dependencies across patient derived models of DIPG, highlighting the therapeutic potential of the blood-brain barrier-penetrant PI3K/Akt/mTOR inhibitor, paxalisib. At the human-equivalent maximum tolerated dose, mice treated with paxalisib experienced systemic glucose feedback and increased insulin levels commensurate with patients using PI3K inhibitors. To exploit genetic dependence and overcome resistance while maintaining compliance and therapeutic benefit, we combined paxalisib with the antihyperglycemic drug metformin. Metformin restored glucose homeostasis and decreased phosphorylation of the insulin receptor in vivo, a common mechanism of PI3K-inhibitor resistance, extending survival of orthotopic models. DIPG models treated with paxalisib increased calcium-activated PKC signaling. The brain penetrant PKC inhibitor enzastaurin, in combination with paxalisib, synergistically extended the survival of multiple orthotopic patient-derived and immunocompetent syngeneic allograft models; benefits potentiated in combination with metformin and standard-of-care radiotherapy. Therapeutic adaptation was assessed using spatial transcriptomics and ATAC-Seq, identifying changes in myelination and tumor immune microenvironment crosstalk. Collectively, this study has identified what we believe to be a clinically relevant DIPG therapeutic combinational strategy.

Transcriptome Analysis Identifies An ASD-Like Phenotype In Oligodendrocytes And Microglia From C58/J Amygdala That Is Dependent On Sex and Sociability

Background

Autism Spectrum Disorder (ASD) is a group of neurodevelopmental disorders with higher incidence in males and is characterized by atypical verbal/nonverbal communication, restricted interests that can be accompanied by repetitive behavior, and disturbances in social behavior. This study investigated brain mechanisms that contribute to sociability deficits and sex differences in an ASD animal model.

Methods

Sociability was measured in C58/J and C57BL/6J mice using the 3-chamber social choice test. Bulk RNA-Seq and snRNA-Seq identified transcriptional changes in C58/J and C57BL/6J amygdala within which DMRseq was used to measure differentially methylated regions in amygdala.

Results

C58/J mice displayed divergent social strata in the 3-chamber test. Transcriptional and pathway signatures revealed immune-related biological processes differ between C58/J and C57BL/6J amygdala. Hypermethylated and hypomethylated genes were identified in C58/J versus C57BL/6J amygdala. snRNA-Seq data in C58/J amygdala identified differential transcriptional signatures within oligodendrocytes and microglia characterized by increased ASD risk gene expression and predicted impaired myelination that was dependent on sex and sociability. RNA velocity, gene regulatory network, and cell communication analysis showed diminished oligodendrocyte/microglia differentiation. Findings were verified using bulk RNA-Seq and demonstrated oxytocin's beneficial effects on myelin gene expression.

Limitations

Our findings are significant. However, limitations can be noted. The cellular mechanisms linking reduced oligodendrocyte differentiation and reduced myelination to an ASD phenotype in C58/J mice need further investigation. Additional snRNA-Seq and spatial studies would determine if effects in oligodendrocytes/microglia are unique to amygdala or if this occurs in other brain regions. Oxytocin's effects need further examination to understand its potential as an ASD therapeutic.

Conclusions

Our work demonstrates the C58/J mouse model's utility in evaluating the influence of sex and sociability on the transcriptome in concomitant brain regions involved in ASD. Our single-nucleus transcriptome analysis elucidates potential pathological roles of oligodendrocytes and microglia in ASD. This investigation provides details regarding regulatory features disrupted in these cell types, including transcriptional gene dysregulation, aberrant cell differentiation, altered gene regulatory networks, and changes to key pathways that promote microglia/oligodendrocyte differentiation. Our studies provide insight into interactions between genetic risk and epigenetic processes associated with divergent affiliative behavior and lack of positive sociability.

UBE2N is essential for maintenance of skin homeostasis and suppression of inflammation

UBE2N, a Lys63-ubiquitin conjugating enzyme, plays critical roles in embryogenesis and immune system development and function. However, its roles in adult epithelial tissue homeostasis and pathogenesis are unclear. We generated conditional mouse models that deleted Ube2n in skin cells in a temporally and spatially controlled manner. We found that Ube2n-knockout (KO) in the adult skin keratinocytes induced a range of inflammatory skin defects characteristic of psoriatic and actinic keratosis. These included eczematous inflammation, epidermal and dermal thickening, parakeratosis, and increased immune cell infiltration, as well as signs of edema and blistering. Single cell transcriptomic analyses and RT-qPCR showed that Ube2n KO keratinocytes expressed elevated myeloid cell chemo-attractants such as Cxcl1 and Cxcl2 and decreased the homeostatic T lymphocyte chemo-attractant, Ccl27a. Consistently, the infiltrating immune cells of Ube2n-KO skin were predominantly myeloid-derived cells including neutrophils and M1-like macrophages that were highly inflammatory, as indicated by expression of Il1β and Il24. Pharmacological blockade of the IL-1 receptor associated kinases (IRAK1/4) alleviated eczema, epidermal and dermal thickening, and immune infiltration of the Ube2n mutant skin. Together, these findings highlight a key role of keratinocyte-UBE2N in maintenance of epidermal homeostasis and skin immunity and identify IRAK1/4 as potential therapeutic target for inflammatory skin disorders.

Ataxia-telangiectasia mutated ( Atm ) disruption sensitizes spatially-directed H3.3K27M/TP53 diffuse midline gliomas to radiation therapy

Diffuse midline gliomas (DMGs) are lethal brain tumors characterized by p53-inactivating mutations and oncohistone H3.3K27M mutations that rewire the cellular response to genotoxic stress, which presents therapeutic opportunities. We used RCAS/tv-a retroviruses and Cre recombinase to inactivate p53 and induce K27M in the native H3f3a allele in a lineage- and spatially-directed manner, yielding primary mouse DMGs. Genetic or pharmacologic disruption of the DNA damage response kinase Ataxia-telangiectasia mutated (ATM) enhanced the efficacy of focal brain irradiation, extending mouse survival. This finding suggests that targeting ATM will enhance the efficacy of radiation therapy for p53-mutant DMG but not p53-wildtype DMG. We used spatial in situ transcriptomics and an allelic series of primary murine DMG models with different p53 mutations to identify transactivation-independent p53 activity as a key mediator of such radiosensitivity. These studies deeply profile a genetically faithful and versatile model of a lethal brain tumor to identify resistance mechanisms for a therapeutic strategy currently in clinical trials.

20-αHydroxycholesterol, an oxysterol in human breast milk, reverses mouse neonatal white matter injury through Gli-dependent oligodendrogenesis

White matter injuries (WMIs) are the leading cause of neurologic impairment in infants born premature. There are no treatment options available. The most common forms of WMIs in infants occur prior to the onset of normal myelination, making its pathophysiology distinctive, thus requiring a tailored approach to treatment. Neonates present a unique opportunity to repair WMIs due to a transient abundance of neural stem/progenitor cells (NSPCs) present in the germinal matrix with oligodendrogenic potential. We identified an endogenous oxysterol, 20-αHydroxycholesterol (20HC), in human maternal breast milk that induces oligodendrogenesis through a sonic hedgehog (shh), Gli-dependent mechanism. Following WMI in neonatal mice, injection of 20HC induced subventricular zone-derived oligodendrogenesis and improved myelination in the periventricular white matter, resulting in improved motor outcomes. Targeting the oligodendrogenic potential of postnatal NSPCs in neonates with WMIs may be further developed into a novel approach to mitigate this devastating complication of preterm birth.

The RNA helicase DDX39B activates FOXP3 RNA splicing to control T regulatory cell fate

Genes associated with increased susceptibility to multiple sclerosis (MS) have been identified, but their functions are incompletely understood. One of these genes codes for the RNA helicase DExD/H-Box Polypeptide 39B (DDX39B), which shows genetic and functional epistasis with interleukin-7 receptor-α gene (IL7R) in MS-risk. Based on evolutionary and functional arguments, we postulated that DDX39B enhances immune tolerance thereby decreasing MS risk. Consistent with such a role we show that DDX39B controls the expression of many MS susceptibility genes and important immune-related genes. Among these we identified Forkhead Box P3 (FOXP3), which codes for the master transcriptional factor in CD4+/CD25+ T regulatory cells. DDX39B knockdown led to loss of immune-regulatory and gain of immune-effector expression signatures. Splicing of FOXP3 introns, which belong to a previously unrecognized type of introns with C-rich polypyrimidine tracts, was exquisitely sensitive to DDX39B levels. Given the importance of FOXP3 in autoimmunity, this work cements DDX39B as an important guardian of immune tolerance.

Adolescent peer struggles predict accelerated epigenetic aging in midlife

This study examined struggles to establish autonomy and relatedness with peers in adolescence and early adulthood as predictors of advanced epigenetic aging assessed at age 30. Participants (N = 154; 67 male and 87 female) were observed repeatedly, along with close friends and romantic partners, from ages 13 through 29. Observed difficulty establishing close friendships characterized by mutual autonomy and relatedness from ages 13 to 18, an interview-assessed attachment state of mind lacking autonomy and valuing of attachment at 24, and self-reported difficulties in social integration across adolescence and adulthood were all linked to greater epigenetic age at 30, after accounting for chronological age, gender, race, and income. Analyses assessing the unique and combined effects of these factors, along with lifetime history of cigarette smoking, indicated that each of these factors, except for adult social integration, contributed uniquely to explaining epigenetic age acceleration. Results are interpreted as evidence that the adolescent preoccupation with peer relationships may be highly functional given the relevance of such relationships to long-term physical outcomes.

Ganglioglioma deep transcriptomics reveals primitive neuroectoderm neural precursor-like population

Gangliogliomas are brain tumors composed of neuron-like and macroglia-like components that occur in children and young adults. Gangliogliomas are often characterized by a rare population of immature astrocyte-appearing cells expressing CD34, a marker expressed in the neuroectoderm (neural precursor cells) during embryogenesis. New insights are needed to refine tumor classification and to identify therapeutic approaches. We evaluated five gangliogliomas with single nucleus RNA-seq, cellular indexing of transcriptomes and epitopes by sequencing, and/or spatially-resolved RNA-seq. We uncovered a population of CD34+ neoplastic cells with mixed neuroectodermal, immature astrocyte, and neuronal markers. Gene regulatory network interrogation in these neuroectoderm-like cells revealed control of transcriptional programming by TCF7L2/MEIS1-PAX6 and SOX2, similar to that found during neuroectodermal/neural development. Developmental trajectory analyses place neuroectoderm-like tumor cells as precursor cells that give rise to neuron-like and macroglia-like neoplastic cells. Spatially-resolved transcriptomics revealed a neuroectoderm-like tumor cell niche with relative lack of vascular and immune cells. We used these high resolution results to deconvolute clinically-annotated transcriptomic data, confirming that CD34+ cell-associated gene programs associate with gangliogliomas compared to other glial brain tumors. Together, these deep transcriptomic approaches characterized a ganglioglioma cellular hierarchy-confirming CD34+ neuroectoderm-like tumor precursor cells, controlling transcription programs, cell signaling, and associated immune cell states. These findings may guide tumor classification, diagnosis, prognostication, and therapeutic investigations.

JAK-STAT activation contributes to cytotoxic T cell-mediated basal cell death in human chronic lung allograft dysfunction

Chronic lung allograft dysfunction (CLAD) is the leading cause of death in lung transplant recipients. CLAD is characterized clinically by a persistent decline in pulmonary function and histologically by the development of airway-centered fibrosis known as bronchiolitis obliterans. There are no approved therapies to treat CLAD, and the mechanisms underlying its development remain poorly understood. We performed single-cell RNA-Seq and spatial transcriptomic analysis of explanted tissues from human lung recipients with CLAD, and we performed independent validation studies to identify an important role of Janus kinase-signal transducer and activator of transcription (JAK-STAT) signaling in airway epithelial cells that contributes to airway-specific alloimmune injury. Specifically, we established that activation of JAK-STAT signaling leads to upregulation of major histocompatibility complex 1 (MHC-I) in airway basal cells, an important airway epithelial progenitor population, which leads to cytotoxic T cell-mediated basal cell death. This study provides mechanistic insight into the cell-to-cell interactions driving airway-centric alloimmune injury in CLAD, suggesting a potentially novel therapeutic strategy for CLAD prevention or treatment.

Genetic and epigenetic signatures associated with plasma oxytocin levels in children and adolescents with autism spectrum disorder

Oxytocin (OT), the brain's most abundant neuropeptide, plays an important role in social salience and motivation. Clinical trials of the efficacy of OT in autism spectrum disorder (ASD) have reported mixed results due in part to ASD's complex etiology. We investigated whether genetic and epigenetic variation contribute to variable endogenous OT levels that modulate sensitivity to OT therapy. To carry out this analysis, we integrated genome-wide profiles of DNA-methylation, transcriptional activity, and genetic variation with plasma OT levels in 290 participants with ASD enrolled in a randomized controlled trial of OT. Our analysis identified genetic variants with novel association with plasma OT, several of which reside in known ASD risk genes. We also show subtle but statistically significant association of plasma OT levels with peripheral transcriptional activity and DNA-methylation profiles across several annotated gene sets. These findings broaden our understanding of the effects of the peripheral oxytocin system and provide novel genetic candidates for future studies to decode the complex etiology of ASD and its interaction with OT signaling and OT-based interventions. LAY SUMMARY: Oxytocin (OT) is an abundant chemical produced by neurons that plays an important role in social interaction and motivation. We investigated whether genetic and epigenetic factors contribute to variable OT levels in the blood. To this, we integrated genetic, gene expression, and non-DNA regulated (epigenetic) signatures with blood OT levels in 290 participants with autism enrolled in an OT clinical trial. We identified genetic association with plasma OT, several of which reside in known autism risk genes. We also show statistically significant association of plasma OT levels with gene expression and epigenetic across several gene pathways. These findings broaden our understanding of the factors that influence OT levels in the blood for future studies to decode the complex presentation of autism and its interaction with OT and OT-based treatment.

Aging and obesity prime the methylome and transcriptome of adipose stem cells for disease and dysfunction

The epigenome of stem cells occupies a critical interface between genes and environment, serving to regulate expression through modification by intrinsic and extrinsic factors. We hypothesized that aging and obesity, which represent major risk factors for a variety of diseases, synergistically modify the epigenome of adult adipose stem cells (ASCs). Using integrated RNA- and targeted bisulfite-sequencing in murine ASCs from lean and obese mice at 5- and 12-months of age, we identified global DNA hypomethylation with either aging or obesity, and a synergistic effect of aging combined with obesity. The transcriptome of ASCs in lean mice was relatively stable to the effects of age, but this was not true in obese mice. Functional pathway analyses identified a subset of genes with critical roles in progenitors and in diseases of obesity and aging. Specifically, Mapt, Nr3c2, App, and Ctnnb1 emerged as potential hypomethylated upstream regulators in both aging and obesity (AL vs. YL and AO vs. YO), and App, Ctnnb1, Hipk2, Id2, and Tp53 exhibited additional effects of aging in obese animals. Furthermore, Foxo3 and Ccnd1 were potential hypermethylated upstream regulators of healthy aging (AL vs. YL), and of the effects of obesity in young animals (YO vs. YL), suggesting that these factors could play a role in accelerated aging with obesity. Finally, we identified candidate driver genes that appeared recurrently in all analyses and comparisons undertaken. Further mechanistic studies are needed to validate the roles of these genes capable of priming ASCs for dysfunction in aging- and obesity-associated pathologies.

Single-cell genome-wide association reveals that a nonsynonymous variant in ERAP1 confers increased susceptibility to influenza virus

During pandemics, individuals exhibit differences in risk and clinical outcomes. Here, we developed single-cell high-throughput human in vitro susceptibility testing (scHi-HOST), a method for rapidly identifying genetic variants that confer resistance and susceptibility. We applied this method to influenza A virus (IAV), the cause of four pandemics since the start of the 20th century. scHi-HOST leverages single-cell RNA sequencing (scRNA-seq) to simultaneously assign genetic identity to cells in mixed infections of cell lines of European, African, and Asian origin, reveal associated genetic variants for viral burden, and identify expression quantitative trait loci. Integration of scHi-HOST with human challenge and experimental validation demonstrated that a missense variant in endoplasmic reticulum aminopeptidase 1 (ERAP1; rs27895) increased IAV burden in cells and human volunteers. rs27895 exhibits population differentiation, likely contributing to greater permissivity of cells from African populations to IAV. scHi-HOST is a broadly applicable method and resource for decoding infectious-disease genetics.

Abstract 646: Liquid biopsy transcriptomics identify pathways associated with poor outcomes and immune phenotypes in men with mCRPC

Androgen receptor signaling inhibitors (ARSi) and taxanes are mainstays for patients with metastatic castration-resistant prostate cancer (mCRPC). However, patient response is heterogeneous, and the molecular underpinnings of treatment resistance are not well elucidated.

To identify clinically meaningful mechanisms of treatment resistance, we performed transcriptome analysis of circulating tumor cells (CTCs) isolated from mCRPC patients enrolled in two independent prospective clinical trials: PROPHECY, a clinical study of patients (n=118) treated with abiraterone or enzalutamide followed by docetaxel; and TAXYNERGY where patients were randomized to docetaxel or cabazitaxel treatment. CTCs were obtained at baseline (before treatment), on treatment and at progression and their comprehensive transcriptomic analysis was correlated with clinical outcomes. To uncover potential involvement of the circulating immune macroenvironment (CIME) in treatment resistance, we performed transcriptomic analysis of matching peripheral blood mononuclear cells (PBMCs) using an established, rigorous, blood-derived transcriptional modular framework.

In PROPHECY, CTC RNA-seq identified that RB loss concurrently with enhanced E2F signaling networks were associated with intrinsic ARSi resistance. Using single sample GSEA (ssGSEA) score, we identified that the RB/E2F common signature at baseline was associated with short PFS (median PFS=6.5 months) and OS (median OS=24.5 months) (hazard ratio (HR) = 3.5; 95% CI 1.5-8.2) in men with mCRPC. We further developed a BRCA-loss transcriptional signature, and validated it in the SU2C mCRPC patient cohort, expanding the identification of patients with BRCA-loss phenotypes beyond genomic loss. Applying this signature to PROPHECY baseline samples, we showed that men with high BRCA-loss scores experienced shorter OS (HR=2.42; 95% CI=1-5.9). Through the comparison of CTC transcriptomic profiles at progression with baseline, we identified an inflammatory response signature in CTCs which was significantly associated with acquired ARSi resistance. Transcriptomic PBMC analysis further identified enrichment of inflammasome gene signatures at progression, with concurrent downregulation of CD8+ T and NK cells.

Furthermore, preliminary data from both clinical trials, showed a significant upregulation of TGF-β1 and corresponding TGFβ-Receptor signaling pathway in CTCs from patients at progression following taxane treatment, suggesting a role for TGFβ pathway in clinical response to taxane chemotherapy.

Taken together, these data demonstrate that liquid biopsy transcriptomics of both tumor cells and immune cells can identify molecular pathways associated with treatment resistance paving the way for treatment optimization and the development of novel precision therapies in patients with mCRPC.

Citation Format: Jiaren Zhang, Bob Zimmermann, Giuseppe Galletti, Susan Halabi, Ada Gjyrezi, Qian Yang, Santosh Gupta, Akanksha Verma, Andrea Sboner, Monika Anand, Daniel J. George, Simon G. Gregory, Prerna Mahtani, Seunghee Hong, Virginia Pascual, Clio P. Mavragani, Emmanuel S. Antonarakis, David M. Nanus, Scott T. Tagawa, Olivier Elemento, Andrew J. Armstrong, Paraskevi Giannakakou. Liquid biopsy transcriptomics identify pathways associated with poor outcomes and immune phenotypes in men with mCRPC [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 646.

Single Cell RNA-Seq Analysis of Human Red Cells

Human red blood cells (RBCs), or erythrocytes, are the most abundant blood cells responsible for gas exchange. RBC diseases affect hundreds of millions of people and impose enormous financial and personal burdens. One well-recognized, but poorly understood feature of RBC populations within the same individual are their phenotypic heterogeneity. The granular characterization of phenotypic RBC variation in normative and disease states may allow us to identify the genetic determinants of red cell diseases and reveal novel therapeutic approaches for their treatment. Previously, we discovered diverse RNA transcripts in RBCs that has allowed us to dissect the phenotypic heterogeneity and malaria resistance of sickle red cells. However, these analyses failed to capture the heterogeneity found in RBC sub-populations. To overcome this limitation, we have performed single cell RNA-Seq to analyze the transcriptional heterogeneity of RBCs from three adult healthy donors which have been stored in the blood bank conditions and assayed at day 1 and day 15. The expression pattern clearly separated RBCs into seven distinct clusters that include one RBC cluster that expresses HBG2 and a small population of RBCs that express fetal hemoglobin (HbF) that we annotated as F cells. Almost all HBG2-expessing cells also express HBB, suggesting bi-allelic expression in single RBC from the HBG2/HBB loci, and we annotated another cluster as reticulocytes based on canonical gene expression. Additional RBC clusters were also annotated based on the enriched expression of NIX, ACVR2B and HEMGN, previously shown to be involved in erythropoiesis. Finally, we found the storage of RBC was associated with an increase in the ACVR2B and F-cell clusters. Collectively, these data indicate the power of single RBC RNA-Seq to capture and discover known and unexpected heterogeneity of RBC population.

Sex-dimorphic gene effects on survival outcomes in people with coronary artery disease

Background:

Ischemic coronary heart disease (IHD) is the leading cause of death worldwide. Genetic variation is presumed to be a major factor underlying sex differences for IHD events, including mortality. The purpose of this study was to identify sex-specific candidate genes associated with all-cause mortality among people diagnosed with coronary artery disease (CAD).

Methods:

We performed a sex-stratified, exploratory genome-wide association (GWAS) screen using existing data from CAD-diagnosed males (n = 510) and females (n = 174) who reported European ancestry from the Duke Catheterization Genetics biorepository. Extant genotype data for 785,945 autosomal SNPs generated with the Human Omni1-Quad BeadChip (Illumina, CA, USA) were analyzed using an additive inheritance model. We estimated instantaneous risk of all-cause mortality by genotype groups across the 11-year follow-up using Cox multivariate regression, covarying for age and genomic ancestry.

Results:

The top GWAS hits associated with all-cause mortality among people with CAD included 8 SNPs among males and 15 among females (p = 1 × 10⁻⁶ or 10⁻⁷), adjusted for covariates. Cross-sex comparisons revealed distinct candidate genes. Biologically relevant candidates included rs9932462 (EMP2/TEKT5) and rs2835913 (KCNJ6) among males and rs7217169 (RAP1GAP2), rs8021816 (PRKD1), rs8133010 (PDE9A), and rs12145981 (LPGAT1) among females.

Conclusions:

We report 20 sex-specific candidate genes having suggestive association with all-cause mortality among CAD-diagnosed subjects. Findings demonstrate proof of principle for identifying sex-associated genetic factors that may help explain differential mortality risk in people with CAD. Replication and meta-analyses in larger studies with more diverse samples will strengthen future work in this area.

Lifetime marijuana use and epigenetic age acceleration: A 17-year prospective examination

AIMS

This study was designed to assess links between lifetime levels of marijuana use and accelerated epigenetic aging.

DESIGN

Prospective longitudinal study, following participants annually from age 13 to age 30.

SETTING AND PARTICIPANTS

A community sample of 154 participants recruited from a small city in the Southeastern United States.

MEASUREMENTS

Participants completed annual assessments of marijuana use from age 13 to age 29 and provided blood samples that yielded two indices of epigenetic aging (DNAmGrimAge and DunedinPoAm) at age 30. Additional covariates examined included history of cigarette smoking, anxiety and depressive symptoms, childhood illness, gender, adolescent-era family income, and racial/ethnic minority status.

FINDINGS

Lifetime marijuana use predicted accelerated epigenetic aging, with effects remaining even after covarying cell counts, demographic factors and chronological age (β's = 0.32 & 0.27, p's < 0.001, 95% CI's = 0.21-0.43 & 0.16-0.39 for DNAmGrimAge and DunedinPoAm, respectively). Predictions remained after accounting for cigarette smoking (β's = 0.25 & 0.21, respectively, p's < 0.001, 95% CI's = 0.14-0.37 & 0.09-0.32 for DNAmGrimAge and DunedinPoAm, respectively). A dose-response effect was observed and there was also evidence that effects were dependent upon recency of use. Effects of marijuana use appeared to be fully mediated by hypomethylation of a site linked to effects of hydrocarbon inhalation (cg05575921).

CONCLUSIONS

Marijuana use predicted epigenetic changes linked to accelerated aging, with evidence suggesting that effects may be primarily due to hydrocarbon inhalation among marijuana smokers. Further research is warranted to explore mechanisms underlying this linkage.

Intranasal Oxytocin in Children and Adolescents with Autism Spectrum Disorder

BACKGROUND

Experimental studies and small clinical trials have suggested that treatment with intranasal oxytocin may reduce social impairment in persons with autism spectrum disorder. Oxytocin has been administered in clinical practice to many children with autism spectrum disorder.

METHODS

We conducted a 24-week, placebo-controlled phase 2 trial of intranasal oxytocin therapy in children and adolescents 3 to 17 years of age with autism spectrum disorder. Participants were randomly assigned in a 1:1 ratio, with stratification according to age and verbal fluency, to receive oxytocin or placebo, administered intranasally, with a total target dose of 48 international units daily. The primary outcome was the least-squares mean change from baseline on the Aberrant Behavior Checklist modified Social Withdrawal subscale (ABC-mSW), which includes 13 items (scores range from 0 to 39, with higher scores indicating less social interaction). Secondary outcomes included two additional measures of social function and an abbreviated measure of IQ.

RESULTS

Of the 355 children and adolescents who underwent screening, 290 were enrolled. A total of 146 participants were assigned to the oxytocin group and 144 to the placebo group; 139 and 138 participants, respectively, completed both the baseline and at least one postbaseline ABC-mSW assessments and were included in the modified intention-to-treat analyses. The least-squares mean change from baseline in the ABC-mSW score (primary outcome) was -3.7 in the oxytocin group and -3.5 in the placebo group (least-squares mean difference, -0.2; 95% confidence interval, -1.5 to 1.0; P = 0.61). Secondary outcomes generally did not differ between the trial groups. The incidence and severity of adverse events were similar in the two groups.

CONCLUSIONS

This placebo-controlled trial of intranasal oxytocin therapy in children and adolescents with autism spectrum disorder showed no significant between-group differences in the least-squares mean change from baseline on measures of social or cognitive functioning over a period of 24 weeks. (Funded by the National Institute of Child Health and Human Development; SOARS-B ClinicalTrials.gov number, NCT01944046.).

Single-Cell RNA Sequencing Reveals Cellular and Transcriptional Changes Associated With M1 Macrophage Polarization in Hidradenitis Suppurativa

Hidradenitis suppurativa (HS) is a chronic inflammatory skin disease characterized by recurrent abscesses, nodules, and sinus tracts in areas of high hair follicle and sweat gland density. These sinus tracts can present with purulent drainage and scar formation. Dysregulation of multiple immune pathways drives the complexity of HS pathogenesis and may account for the heterogeneity of treatment response in HS patients. Using transcriptomic approaches, including single-cell sequencing and protein analysis, we here characterize the innate inflammatory landscape of HS lesions. We identified a shared upregulation of genes involved in interferon (IFN) and antimicrobial defense signaling through transcriptomic overlap analysis of differentially expressed genes (DEGs) in datasets from HS skin, diabetic foot ulcers (DFUs), and the inflammatory stage of normal healing wounds. Overlap analysis between HS- and DFU-specific DEGs revealed an enrichment of gene signatures associated with monocyte/macrophage functions. Single-cell RNA sequencing further revealed monocytes/macrophages with polarization toward a pro-inflammatory M1-like phenotype and increased effector function, including antiviral immunity, phagocytosis, respiratory burst, and antibody-dependent cellular cytotoxicity. Specifically, we identified the STAT1/IFN-signaling axis and the associated IFN-stimulated genes as central players in monocyte/macrophage dysregulation. Our data indicate that monocytes/macrophages are a potential pivotal player in HS pathogenesis and their pathways may serve as therapeutic targets and biomarkers in HS treatment.

Circulating Tumor Cell Chromosomal Instability and Neuroendocrine Phenotype by Immunomorphology and Poor Outcomes in Men with mCRPC Treated with Abiraterone or Enzalutamide

PURPOSE

While the detection of AR-V7 in circulating tumor cells (CTC) is associated with resistance to abiraterone or enzalutamide in men with metastatic castration-resistant prostate cancer (mCRPC), it only accounts for a minority of this resistance. Neuroendocrine (NE) differentiation or chromosomal instability (CIN) may be additional mechanisms that mediate resistance.

EXPERIMENTAL DESIGN

PROPHECY was a multicenter prospective study of men with high-risk mCRPC starting abiraterone or enzalutamide. A secondary objective was to assess Epic CTC CIN and NE phenotypes before abiraterone or enzalutamide and at progression. The proportional hazards (PH) model was used to investigate the prognostic importance of CIN and NE in predicting progression-free survival and overall survival (OS) adjusting for CTC number (CellSearch), AR-V7, prior therapy, and clinical risk score. The PH model was utilized to validate this association of NE with OS in an external dataset of patients treated similarly at Memorial Sloan Kettering Cancer Center (MSKCC; New York, NY).

RESULTS

We enrolled 118 men with mCRPC starting on abiraterone or enzalutamide; 107 were evaluable on the Epic platform. Of these, 36.4% and 8.4% were CIN positive and NE positive, respectively. CIN and NE were independently associated with worse OS [HR, 2.2; 95% confidence interval (CI), 1.2-4.0 and HR 3.8; 95% CI, 1.2-12.3, respectively] when treated with abiraterone/enzalutamide. The prognostic significance of NE positivity for worse OS was confirmed in the MSKCC dataset (n = 173; HR, 5.7; 95% CI, 2.6-12.7).

CONCLUSIONS

A high CIN and NE CTC phenotype is independently associated with worse survival in men with mCRPC treated with abiraterone/enzalutamide, warranting further prospective controlled predictive studies to inform treatment decisions.

A non-canonical type 2 immune response coordinates tuberculous granuloma formation and epithelialization

The central pathogen-immune interface in tuberculosis is the granuloma, a complex host immune structure that dictates infection trajectory and physiology. Granuloma macrophages undergo a dramatic transition in which entire epithelial modules are induced and define granuloma architecture. In tuberculosis, relatively little is known about the host signals that trigger this transition. Using the zebrafish-Mycobacterium marinum model, we identify the basis of granuloma macrophage transformation. Single-cell RNA-sequencing analysis of zebrafish granulomas and analysis of Mycobacterium tuberculosis-infected macaques reveal that, even in the presence of robust type 1 immune responses, countervailing type 2 signals associate with macrophage epithelialization. We find that type 2 immune signaling, mediated via stat6, is absolutely required for epithelialization and granuloma formation. In mixed chimeras, stat6 acts cell autonomously within macrophages, where it is required for epithelioid transformation and incorporation into necrotic granulomas. These findings establish the signaling pathway that produces the hallmark structure of mycobacterial infection.

Circulating Tumor Cell Genomic Evolution and Hormone Therapy Outcomes in Men with Metastatic Castration-Resistant Prostate Cancer

Men with circulating tumor cell (CTC) AR-V7-positive metastatic castration-resistant prostate cancer (mCRPC) have worse outcomes when treated with enzalutamide/abiraterone. However, most men lack CTC AR-V7 detection, and additional predictive biomarkers are needed. We conducted a retrospective secondary analysis of the prospective PROPHECY trial (NCT02269982) of men with mCRPC undergoing treatment with enzalutamide/abiraterone, analyzing pooled CTC and germline DNA for whole-genome copy-number alterations (CNA) in 73 samples from 48 men over time along with pooled CTC and germline whole-exome sequencing on 22 paired samples before and following progression on androgen receptor (AR) inhibitor therapy to identify somatic genomic alterations associated with acquired resistance. We observed broad interpatient and longitudinal CTC genomic heterogeneity from AR-V7-negative men with mCRPC, including common gains of KDM6A, MYCN, and AR, and loss of ZFHX3, BRCA1, and PTEN. Men who had progression-free survival of ≤3 months despite enzalutamide/abiraterone treatment were more likely to have baseline CTC genomic loss of CHD1, PTEN, PHLPP1, and ZFHX3 and gains of BRCA2, KDM5D, MYCN, and SPARC. After progression on abiraterone/enzalutamide, we observed clonal evolution of CTCs harboring TP53 mutations and gain of ATM, KDM6A, and MYC, and loss of NCOR1, PTEN, RB1, and RUNX2. CTC genomic findings were independently confirmed in a separate cohort of mCRPC men who progressed despite prior treatment with abiraterone/enzalutamide (NCT02204943). IMPLICATIONS: We identified common and reproducible genomic alterations in CTCs from AR-V7-negative mCRPC men associated with poor outcomes during enzalutamide/abiraterone treatment, including CNAs in genes linked to lineage plasticity and epigenetic signaling, DNA repair, AR, TP53/RB1, PTEN, and WNT pathways.

U2AF2 binds IL7R exon 6 ectopically and represses its inclusion

Interleukin 7 receptor α-chain is crucial for the development and maintenance of T cells and is genetically associated with autoimmune disorders including multiple sclerosis (MS), a demyelinating disease of the CNS. Exon 6 of IL7R encodes for the transmembrane domain of the receptor and is regulated by alternative splicing: inclusion or skipping of IL7R exon 6 results in membrane-bound or soluble IL7R isoforms, respectively. We previously identified a SNP (rs6897932) in IL7R exon 6, strongly associated with MS risk and showed that the risk allele (C) increases skipping of the exon, resulting in elevated levels of sIL7R. This has important pathological consequences as elevated levels of sIL7R has been shown to exacerbate the disease in the experimental autoimmune encephalomyelitis mouse model of MS. Understanding the regulation of exon 6 splicing provides important mechanistic insights into the pathogenesis of MS. Here we report two mechanisms by which IL7R exon 6 is controlled. First, a competition between PTBP1 and U2AF2 at the polypyrimidine tract (PPT) of intron 5, and second, an unexpected U2AF2-mediated assembly of spicing factors in the exon. We noted the presence of a branchpoint sequence (BPS) (TACTAAT or TACTAAC) within exon 6, which is stronger with the C allele. We also noted that the BPS is followed by a PPT and conjectured that silencing could be mediated by the binding of U2AF2 to that tract. In support of this model, we show that evolutionary conservation of the exonic PPT correlates well with the degree of alternative splicing of exon 6 in two non-human primate species and that U2AF2 binding to this PPT recruits U2 snRNP components to the exon. These observations provide the first explanation for the stronger silencing of IL7R exon 6 with the disease associated C allele at rs6897932.

Genetic, epigenetic, and environmental factors controlling oxytocin receptor gene expression

BACKGROUND

The neuropeptide oxytocin regulates mammalian social behavior. Disruptions in oxytocin signaling are a feature of many psychopathologies. One commonly studied biomarker for oxytocin involvement in psychiatric diseases is DNA methylation at the oxytocin receptor gene (OXTR). Such studies focus on DNA methylation in two regions of OXTR, exon 3 and a region termed MT2 which overlaps exon 1 and intron 1. However, the relative contribution of exon 3 and MT2 in regulating OXTR gene expression in the brain is currently unknown.

RESULTS

Here, we use the prairie vole as a translational animal model to investigate genetic, epigenetic, and environmental factors affecting Oxtr gene expression in a region of the brain that has been shown to drive Oxtr related behavior in the vole, the nucleus accumbens. We show that the genetic structure of Oxtr in prairie voles resembles human OXTR. We then studied the effects of early life experience on DNA methylation in two regions of a CpG island surrounding the Oxtr promoter: MT2 and exon 3. We show that early nurture in the form of parental care results in DNA hypomethylation of Oxtr in both MT2 and exon 3, but only DNA methylation in MT2 is associated with Oxtr gene expression. Network analyses indicate that CpG sites in the 3' portion of MT2 are most highly associated with Oxtr gene expression. We also identify two novel SNPs in exon 3 of Oxtr in prairie voles and a novel alternative transcript originating from the third intron of the gene. Expression of the novel alternative transcript is associated with genotype at SNP KLW2.

CONCLUSIONS

These results identify putative regulatory features of Oxtr in prairie voles which inform future studies examining OXTR in human social behaviors and disorders. These studies indicate that in prairie voles, DNA methylation in MT2, particularly in the 3' portion, is more predictive of Oxtr gene expression than DNA methylation in exon 3. Similarly, in human temporal cortex, we find that DNA methylation in the 3' portion of MT2 is associated with OXTR expression. Together, these results suggest that among the CpG sites studied, DNA methylation of MT2 may be the most reliable indicator of OXTR gene expression. We also identify novel features of prairie vole Oxtr, including SNPs and an alternative transcript, which further develop the prairie vole as a translational model for studies of OXTR.

Single-cell RNA-seq reveals transcriptomic heterogeneity mediated by host-pathogen dynamics in lymphoblastoid cell lines

Lymphoblastoid cell lines (LCLs) are generated by transforming primary B cells with Epstein–Barr virus (EBV) and are used extensively as model systems in viral oncology, immunology, and human genetics research. In this study, we characterized single-cell transcriptomic profiles of five LCLs and present a simple discrete-time simulation to explore the influence of stochasticity on LCL clonal evolution. Single-cell RNA sequencing (scRNA-seq) revealed substantial phenotypic heterogeneity within and across LCLs with respect to immunoglobulin isotype; virus-modulated host pathways involved in survival, activation, and differentiation; viral replication state; and oxidative stress. This heterogeneity is likely attributable to intrinsic variance in primary B cells and host–pathogen dynamics. Stochastic simulations demonstrate that initial primary cell heterogeneity, random sampling, time in culture, and even mild differences in phenotype-specific fitness can contribute substantially to dynamic diversity in populations of nominally clonal cells.

Rationale, design, and methods of the Autism Centers of Excellence (ACE) network Study of Oxytocin in Autism to improve Reciprocal Social Behaviors (SOARS-B)

OBJECTIVE

To describe the rationale, design, and methods of the Autism Centers of Excellence (ACE) network Study of Oxytocin in Autism to improve Reciprocal Social Behaviors (SOARS-B).

METHOD

This phase 2 clinical trial was designed to evaluate the use of intranasal oxytocin treatment to improve social difficulties in individuals with autism spectrum disorder (ASD). In total, 290 participants ages 3 to 17 years with a DSM-5 diagnosis of ASD were enrolled to receive 24 weeks of treatment with either oxytocin or a matched placebo at one of seven collaborating sites. Participants were subsequently treated with open-label oxytocin for 24 additional weeks. Post-treatment assessments were done approximately 4 weeks after treatment discontinuation. Plasma oxytocin and oxytocin receptor gene (OXTR) methylation level were measured at baseline, and week 8, 24 and 36 to explore potential relationships between these biomarkers and treatment response.

RESULTS

This report describes the rationale, design, and methods of the SOARS-B clinical trial.

CONCLUSIONS

There is a tremendous unmet need for safe and effective pharmacological treatment options that target the core symptoms of ASD. Several studies support the hypothesis that intranasal oxytocin could improve social orienting and the salience of social rewards in ASD, thereby enhancing reciprocal social behaviors. However, due to conflicting results from a number of pilot studies on the prosocial effects of exogenous oxytocin, this hypothesis remains controversial and inconclusive. SOARS-B is the best powered study to date to address this hypothesis and promises to improve our understanding of the safety and efficacy of intranasal oxytocin in the treatment of social deficits in children with ASD.

Author Correction: Meteorin-like facilitates skeletal muscle repair through a Stat3/IGF-1 mechanism

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Circulating MicroRNA Profiling in Non-ST Elevated Coronary Artery Syndrome Highlights Genomic Associations with Serial Platelet Reactivity Measurements

Changes in platelet physiology are associated with simultaneous changes in microRNA concentrations, suggesting a role for microRNA in platelet regulation. Here we investigated potential associations between microRNA and platelet reactivity (PR), a marker of platelet function, in two cohorts following a non-ST elevation acute coronary syndrome (NSTE-ACS) event. First, non-targeted microRNA concentrations and PR were compared in a case (N = 77) control (N = 76) cohort within the larger TRILOGY-ACS trial. MicroRNA significant in this analysis plus CVD-associated microRNAs from the literature were then quantified by targeted rt-PCR in the complete TRILOGY-ACS cohort (N = 878) and compared with matched PR samples. Finally, microRNA significant in the non-targeted & targeted analyses were verified in an independent post NSTE-ACS cohort (N = 96). From the non-targeted analysis, 14 microRNAs were associated with PR (Fold Change: 0.91–1.27, p-value: 0.004–0.05). From the targeted analysis, five microRNAs were associated with PR (Beta: −0.09–0.22, p-value: 0.004–0.05). Of the 19 significant microRNAs, three, miR-15b-5p, miR-93 and miR-126, were consistently associated with PR in the TRILOGY-ACS and independent Singapore post-ACS cohorts, suggesting the measurement of circulating microRNA concentrations may report on dynamic changes in platelet biology following a cardiovascular ischemic event.

Meteorin-like facilitates skeletal muscle repair through a Stat3/IGF-1 mechanism

The immune system plays a multifunctional role throughout the regenerative process, regulating both pro-/anti-inflammatory phases and progenitor cell function. In the present study, we identify the myokine/cytokine Meteorin-like (Metrnl) as a critical regulator of muscle regeneration. Mice genetically lacking Metrnl have impaired muscle regeneration associated with a reduction in immune cell infiltration and an inability to transition towards an anti-inflammatory phenotype. Isochronic parabiosis, joining wild-type and whole-body Metrnl knock-out (KO) mice, returns Metrnl expression in the injured muscle and improves muscle repair, providing supportive evidence for Metrnl secretion from infiltrating immune cells. Macrophage-specific Metrnl KO mice are also deficient in muscle repair. During muscle regeneration, Metrnl works, in part, through Stat3 activation in macrophages, resulting in differentiation to an anti-inflammatory phenotype. With regard to myogenesis, Metrnl induces macrophage-dependent insulin-like growth factor 1 production, which has a direct effect on primary muscle satellite cell proliferation. Perturbations in this pathway inhibit efficacy of Metrnl in the regenerative process. Together, these studies identify Metrnl as an important regulator of muscle regeneration and a potential therapeutic target to enhance tissue repair.

Single-Cell RNA Sequencing Identifies Yes-Associated Protein 1-Dependent Hepatic Mesothelial Progenitors in Fibrolamellar Carcinoma

Fibrolamellar carcinoma (FLC) is characterized by in-frame fusion of DnaJ heat shock protein family (Hsp40) member B1 (DNAJB1) with protein kinase cAMP-activated catalytic subunit α (PRKACA) and by dense desmoplasia. Surgery is the only effective treatment because mechanisms supporting tumor survival are unknown. We used single-cell RNA sequencing to characterize a patient-derived FLC xenograft model and identify therapeutic targets. Human FLC cells segregated into four discrete clusters that all expressed the oncogene Yes-associated protein 1 (YAP1). The two communities most enriched with cells coexpressing FLC markers [CD68, A-kinase anchoring protein 12 (AKAP12), cytokeratin 7, epithelial cell adhesion molecule (EPCAM), and carbamoyl palmitate synthase-1] also had the most cells expressing YAP1 and its proproliferative target genes (AREG and CCND1), suggesting these were proliferative FLC cell clusters. The other two clusters were enriched with cells expressing profibrotic YAP1 target genes, ACTA2, ELN, and COL1A1, indicating these were fibrogenic FLC cells. All clusters expressed the YAP1 target gene and mesothelial progenitor marker mesothelin, and many mesothelin-positive cells coexpressed albumin. Trajectory analysis predicted that the four FLC communities were derived from a single cell type transitioning among phenotypic states. After establishing a novel FLC cell line that harbored the DNAJB1-PRKACA fusion, YAP1 was inhibited, which significantly reduced expression of known YAP1 target genes as well as cell growth and migration. Thus, both FLC epithelial and stromal cells appear to arise from DNAJB1-PRKACA fusion in a YAP1-dependent liver mesothelial progenitor, identifying YAP1 as a target for FLC therapy.

Discordant and heterogeneous clinically relevant genomic alterations in circulating tumor cells vs plasma DNA from men with metastatic castration resistant prostate cancer

Circulating tumor cell (CTC) and cell-free (cf) DNA-based genomic alterations are increasingly being used for clinical decision-making in oncology. However, the concordance and discordance between paired CTC and cfDNA genomic profiles remain largely unknown. We performed comparative genomic hybridization (CGH) on CTCs and cfDNA, and low-pass whole genome sequencing (lpWGS) on cfDNA to characterize genomic alterations (CNA) and tumor content in two independent prospective studies of 93 men with mCRPC treated with enzalutamide/abiraterone, or radium-223. Comprehensive analysis of 69 patient CTCs and 72 cfDNA samples from 93 men with mCRPC, including 64 paired samples, identified common concordant gains in FOXA1, AR, and MYC, and losses in BRCA1, PTEN, and RB1 between CTCs and cfDNA. Concordant PTEN loss and discordant BRCA2 gain were associated with significantly worse outcomes in Epic AR-V7 negative men with mCRPC treated with abiraterone/enzalutamide. We identified and externally validated CTC-specific genomic alternations that were discordant in paired cfDNA, even in samples with high tumor content. These CTC/cfDNA-discordant regions included key genomic regulators of lineage plasticity, osteomimicry, and cellular differentiation, including MYCN gain in CTCs (31%) that was rarely detected in cfDNA. CTC MYCN gain was associated with poor clinical outcomes in AR-V7 negative men and small cell transformation. In conclusion, we demonstrated concordance of multiple genomic alterations across CTC and cfDNA platforms; however, some genomic alterations displayed substantial discordance between CTC DNA and cfDNA despite the use of identical copy number analysis methods, suggesting tumor heterogeneity and divergent evolution associated with poor clinical outcomes.

Epigenome-Wide Association Study for All-Cause Mortality in a Cardiovascular Cohort Identifies Differential Methylation in Castor Zinc Finger 1 (CASZ1)

Background DNA methylation is implicated in many chronic diseases and may contribute to mortality. Therefore, we conducted an epigenome-wide association study (EWAS) for all-cause mortality with whole-transcriptome data in a cardiovascular cohort (CATHGEN [Catheterization Genetics]). Methods and Results Cases were participants with mortality≥7 days postcatheterization whereas controls were alive with≥2 years of follow-up. The Illumina Human Methylation 450K and EPIC arrays (Illumina, San Diego, CA) were used for the discovery and validation sets, respectively. A linear model approach with empirical Bayes estimators adjusted for confounders was used to assess difference in methylation (Δβ). In the discovery set (55 cases, 49 controls), 25 629 (6.5%) probes were differently methylated (P<0.05). In the validation set (108 cases, 108 controls), 3 probes were differentially methylated with a false discovery rate-adjusted P<0.10: cg08215811 (SLC4A9; log2 fold change=-0.14); cg17845532 (MATK; fold change=-0.26); and cg17944110 (castor zinc finger 1 [CASZ1]; FC=0.26; P<0.0001; false discovery rate-adjusted P=0.046-0.080). Meta-analysis identified 6 probes (false discovery rate-adjusted P<0.05): the 3 above, cg20428720 (intergenic), cg17647904 (NCOR2), and cg23198793 (CAPN3). Messenger RNA expression of 2 MATK isoforms was lower in cases (fold change=-0.24 [P=0.007] and fold change=-0.61 [P=0.009]). The CASZ1, NCOR2, and CAPN3 transcripts did not show differential expression (P>0.05); the SLC4A9 transcript did not pass quality control. The cg17944110 probe is located within a potential regulatory element; expression of predicted targets (using GeneHancer) of the regulatory element, UBIAD1 (P=0.01) and CLSTN1 (P=0.03), were lower in cases. Conclusions We identified 6 novel methylation sites associated with all-cause mortality. Methylation in CASZ1 may serve as a regulatory element associated with mortality in cardiovascular patients. Larger studies are necessary to confirm these observations.

Associations of osteopontin and NT-proBNP with circulating miRNA levels in acute coronary syndrome

The genomic regulatory networks underlying the pathogenesis of non-ST-segment elevation acute coronary syndrome (NSTE-ACS) are incompletely understood. As intermediate traits, protein biomarkers report on underlying disease severity and prognosis in NSTE-ACS. We hypothesized that integration of dense microRNA (miRNA) profiling with biomarker measurements would highlight potential regulatory pathways that underlie the relationships between prognostic biomarkers, miRNAs, and cardiovascular phenotypes. We performed miRNA sequencing using whole blood from 186 patients from the TRILOGY-ACS trial. Seven circulating prognostic biomarkers were measured: NH₂-terminal pro-B-type natriuretic peptide (NT-proBNP), high-sensitivity C-reactive protein, osteopontin (OPN), myeloperoxidase, growth differentiation factor 15, monocyte chemoattractant protein, and neopterin. We tested miRNAs for association with each biomarker with generalized linear models and controlled the false discovery rate at 0.05. Ten miRNAs, including known cardiac-related miRNAs 25-3p and 423-3p, were associated with NT-proBNP levels (min. P = 7.5 × 10⁻⁴) and 48 miRNAs, including cardiac-related miRNAs 378a-3p, 20b-5p and 320a, -b, and -d, were associated with OPN levels (min. P = 1.6 × 10⁻⁶). NT-proBNP and OPN were also associated with time to cardiovascular death, myocardial infarction (MI), or stroke in the sample. By integrating large-scale miRNA profiling with circulating biomarkers as intermediate traits, we identified associations of known cardiac-related and novel miRNAs with two prognostic biomarkers and identified potential genomic networks regulating these biomarkers. These results, highlighting plausible biological pathways connecting miRNAs with biomarkers and outcomes, may inform future studies seeking to delineate genomic pathways underlying NSTE-ACS outcomes.

Pharmacodynamic study of radium-223 in men with bone metastatic castration resistant prostate cancer

BACKGROUND

Radium-223 is a targeted alpha-particle therapy that improves survival in men with metastatic castration resistant prostate cancer (mCRPC), particularly in men with elevated serum levels of bone alkaline phosphatase (B-ALP). We hypothesized that osteomimicry, a form of epithelial plasticity leading to an osteoblastic phenotype, may contribute to intralesional deposition of radium-223 and subsequent irradiation of the tumor microenvironment.

METHODS

We conducted a pharmacodynamic study (NCT02204943) of radium-223 in men with bone mCRPC. Prior to and three and six months after radium-223 treatment initiation, we collected CTCs and metastatic biopsies for phenotypic characterization and CTC genomic analysis. The primary objective was to describe the impact of radium-223 on the prevalence of CTC B-ALP over time. We measured radium-223 decay products in tumor and surrounding normal bone during treatment. We validated genomic findings in a separate independent study of men with bone metastatic mCRPC (n = 45) and publicly accessible data of metastatic CRPC tissues.

RESULTS

We enrolled 20 men with symptomatic bone predominant mCRPC and treated with radium-223. We observed greater radium-223 radioactivity levels in metastatic bone tumor containing biopsies compared with adjacent normal bone. We found evidence of persistent Cellsearch CTCs and B-ALP (+) CTCs in the majority of men over time during radium-223 therapy despite serum B-ALP normalization. We identified genomic gains in osteoblast mimicry genes including gains of ALPL, osteopontin, SPARC, OB-cadherin and loss of RUNX2, and validated genomic alterations or increased expression at the DNA and RNA level in an independent cohort of 45 men with bone-metastatic CRPC and in 150 metastatic biopsies from men with mCRPC.

CONCLUSIONS

Osteomimicry may contribute in part to the uptake of radium-223 within bone metastases and may thereby enhance the therapeutic benefit of this bone targeting radiotherapy.

Whole Exome Sequencing of Cell-Free DNA for Early Lung Cancer: A Pilot Study to Differentiate Benign From Malignant CT-Detected Pulmonary Lesions

Introduction: Indeterminate pulmonary lesions (IPL) detected by CT pose a significant clinical challenge, frequently necessitating long-term surveillance or biopsy for diagnosis. In this pilot investigation, we performed whole exome sequencing (WES) of plasma cell free (cfDNA) and matched germline DNA in patients with CT-detected pulmonary lesions to determine the feasibility of somatic cfDNA mutations to differentiate benign from malignant pulmonary nodules. Methods: 33 patients with a CT-detected pulmonary lesions were retrospectively enrolled (n = 16 with a benign nodule, n = 17 with a malignant nodule). Following isolation and amplification of plasma cfDNA and matched peripheral blood mononuclear cells (PBMC) from patient blood samples, WES of cfDNA and PBMC DNA was performed. After genomic alignment and filtering, we looked for lung-cancer associated driver mutations and next identified high-confidence somatic variants in both groups. Results: Somatic cfDNA mutations were observed in both groups, with the cancer group demonstrating more variants than the benign group (1083 ± 476 versus 553 ± 519, p < 0.0046). By selecting variants present in >2 cancer patients and not the benign group, we accurately identified 82% (14/17) of cancer patients. Conclusions: This study suggests a potential role for cfDNA for the early identification of lung cancer in patients with CT-detected pulmonary lesions. Importantly, a substantial number of somatic variants in healthy patients with benign pulmonary nodules were also found. Such "benign" variants, while largely unexplored to date, have widespread relevance to all liquid biopsies if cfDNA is to be used accurately for cancer detection.

Prospective Multicenter Validation of Androgen Receptor Splice Variant 7 and Hormone Therapy Resistance in High-Risk Castration-Resistant Prostate Cancer: The PROPHECY Study

PURPOSE

Androgen receptor splice variant 7 (AR-V7) results in a truncated receptor, which leads to ligand-independent constitutive activation that is not inhibited by anti-androgen therapies, including abiraterone or enzalutamide. Given that previous reports suggested that circulating tumor cell (CTC) AR-V7 detection is a poor prognostic indicator for the clinical efficacy of secondary hormone therapies, we conducted a prospective multicenter validation study.

PATIENTS AND METHODS

PROPHECY (ClinicalTrials.gov identifier: NCT02269982) is a multicenter, prospective-blinded study of men with high-risk mCRPC starting abiraterone acetate or enzalutamide treatment. The primary objective was to validate the prognostic significance of baseline CTC AR-V7 on the basis of radiographic or clinical progression free-survival (PFS) by using the Johns Hopkins University modified-AdnaTest CTC AR-V7 mRNA assay and the Epic Sciences CTC nuclear-specific AR-V7 protein assay. Overall survival (OS) and prostate-specific antigen responses were secondary end points.

RESULTS

We enrolled 118 men with mCRPC who were starting abiraterone or enzalutamide treatment. AR-V7 detection by both the Johns Hopkins and Epic AR-V7 assays was independently associated with shorter PFS (hazard ratio, 1.9 [95% CI, 1.1 to 3.3; P = .032] and 2.4 [95% CI, 1.1 to 5.1; P = .020], respectively) and OS (hazard ratio, 4.2 [95% CI, 2.1 to 8.5] and 3.5 [95% CI, 1.6 to 8.1], respectively) after adjusting for CTC number and clinical prognostic factors. Men with AR-V7–positive mCRPC had fewer confirmed prostate-specific antigen responses (0% to 11%) or soft tissue responses (0% to 6%). The observed percentage agreement between the two AR-V7 assays was 82%.

CONCLUSION

Detection of AR-V7 in CTCs by two blood-based assays is independently associated with shorter PFS and OS with abiraterone or enzalutamide, and such men with mCRPC should be offered alternative treatments.

Early nurture epigenetically tunes the oxytocin receptor

Mammalian sociality is regulated in part by the neuropeptide oxytocin. In prairie voles, subtle variation in early life experience changes oxytocin receptor-mediated social behaviors. We report that low levels of early care in voles leads to de novo DNA methylation at specific regulatory sites in the oxytocin receptor gene (Oxtr), impacting gene expression and protein distribution in the nucleus accumbens. DNA methylation state of the blood predicts expression in the brain indicating the utility of the blood as a biomarker for the transcription state of the brain. These experience-sensitive CpG sites are conserved in humans, are related to gene expression in the brain, and have been associated with psychiatric disorders and individual differences in neural response to social stimuli. These results identify a mechanism by which early care regulates later displays of typical prairie vole social behavior and suggest the potential for nurture driven epigenetic tuning of OXTR in humans.

Associations Between Residential Proximity to Traffic and Vascular Disease in a Cardiac Catheterization Cohort

OBJECTIVE

Exposure to mobile source emissions is nearly ubiquitous in developed nations and is associated with multiple adverse health outcomes. There is an ongoing need to understand the specificity of traffic exposure associations with vascular outcomes, particularly in individuals with cardiovascular disease.

APPROACH AND RESULTS

We performed a cross-sectional study using 2124 individuals residing in North Carolina, United States, who received a cardiac catheterization at the Duke University Medical Center. Traffic-related exposure was assessed via 2 metrics: (1) the distance between the primary residence and the nearest major roadway; and (2) location of the primary residence in regions defined based on local traffic patterns. We examined 4 cardiovascular disease outcomes: hypertension, peripheral arterial disease, the number of diseased coronary vessels, and recent myocardial infarction. Statistical models were adjusted for race, sex, smoking, type 2 diabetes mellitus, body mass index, hyperlipidemia, and home value. Results are expressed in terms of the odds ratio (OR). A 23% decrease in residential distance to major roadways was associated with higher prevalence of peripheral arterial disease (OR=1.29; 95% confidence interval, 1.08-1.55) and hypertension (OR=1.15; 95% confidence interval, 1.01-1.31). Associations with peripheral arterial disease were strongest in men (OR=1.42; 95% confidence interval, 1.17-1.74) while associations with hypertension were strongest in women (OR=1.21; 95% confidence interval, 0.99-1.49). Neither myocardial infarction nor the number of diseased coronary vessels were associated with traffic exposure.

CONCLUSIONS

Traffic-related exposure is associated with peripheral arterial disease and hypertension while no associations are observed for 2 coronary-specific vascular outcomes.

Human Epistatic Interaction Controls IL7R Splicing and Increases Multiple Sclerosis Risk

Multiple sclerosis (MS) is an autoimmune disorder where T cells attack neurons in the central nervous system (CNS) leading to demyelination and neurological deficits. A driver of increased MS risk is the soluble form of the interleukin-7 receptor alpha chain gene (sIL7R) produced by alternative splicing of IL7R exon 6. Here, we identified the RNA helicase DDX39B as a potent activator of this exon and consequently a repressor of sIL7R, and we found strong genetic association of DDX39B with MS risk. Indeed, we showed that a genetic variant in the 5' UTR of DDX39B reduces translation of DDX39B mRNAs and increases MS risk. Importantly, this DDX39B variant showed strong genetic and functional epistasis with allelic variants in IL7R exon 6. This study establishes the occurrence of biological epistasis in humans and provides mechanistic insight into the regulation of IL7R exon 6 splicing and its impact on MS risk.

Genome-wide association study identifies three novel loci in Fuchs endothelial corneal dystrophy

The structure of the cornea is vital to its transparency, and dystrophies that disrupt corneal organization are highly heritable. To understand the genetic aetiology of Fuchs endothelial corneal dystrophy (FECD), the most prevalent corneal disorder requiring transplantation, we conducted a genome-wide association study (GWAS) on 1,404 FECD cases and 2,564 controls of European ancestry, followed by replication and meta-analysis, for a total of 2,075 cases and 3,342 controls. We identify three novel loci meeting genome-wide significance (P<5 × 10-8): KANK4 rs79742895, LAMC1 rs3768617 and LINC00970/ATP1B1 rs1200114. We also observe an overwhelming effect of the established TCF4 locus. Interestingly, we detect differential sex-specific association at LAMC1, with greater risk in women, and TCF4, with greater risk in men. Combining GWAS results with biological evidence we expand the knowledge of common FECD loci from one to four, and provide a deeper understanding of the underlying pathogenic basis of FECD.

A genome-wide trans-ethnic interaction study links the PIGR-FCAMR locus to coronary atherosclerosis via interactions between genetic variants and residential exposure to traffic

Air pollution is a worldwide contributor to cardiovascular disease mortality and morbidity. Traffic-related air pollution is a widespread environmental exposure and is associated with multiple cardiovascular outcomes such as coronary atherosclerosis, peripheral arterial disease, and myocardial infarction. Despite the recognition of the importance of both genetic and environmental exposures to the pathogenesis of cardiovascular disease, studies of how these two contributors operate jointly are rare. We performed a genome-wide interaction study (GWIS) to examine gene-traffic exposure interactions associated with coronary atherosclerosis. Using race-stratified cohorts of 538 African-Americans (AA) and 1562 European-Americans (EA) from a cardiac catheterization cohort (CATHGEN), we identify gene-by-traffic exposure interactions associated with the number of significantly diseased coronary vessels as a measure of chronic atherosclerosis. We found five suggestive (P<1x10^-5) interactions in the AA GWIS, of which two (rs1856746 and rs2791713) replicated in the EA cohort (P < 0.05). Both SNPs are in the PIGR-FCAMR locus and are eQTLs in lymphocytes. The protein products of both PIGR and FCAMR are implicated in inflammatory processes. In the EA GWIS, there were three suggestive interactions; none of these replicated in the AA GWIS. All three were intergenic; the most significant interaction was in a regulatory region associated with SAMSN1, a gene previously associated with atherosclerosis and B cell activation. In conclusion, we have uncovered several novel genes associated with coronary atherosclerosis in individuals chronically exposed to increased ambient concentrations of traffic air pollution. These genes point towards inflammatory pathways that may modify the effects of air pollution on cardiovascular disease risk.

Human centromere repositioning within euchromatin after partial chromosome deletion

Centromeres are defined by a specialized chromatin organization that includes nucleosomes that contain the centromeric histone variant centromere protein A (CENP-A) instead of canonical histone H3. Studies in various organisms have shown that centromeric chromatin (i.e., CENP-A chromatin or centrochromatin) exhibits plasticity, in that it can assemble on different types of DNA sequences. However, once established on a chromosome, the centromere is maintained at the same position. In humans, this location is the highly homogeneous repetitive DNA alpha satellite. Mislocalization of centromeric chromatin to atypical locations can lead to genome instability, indicating that restriction of centromeres to a distinct genomic position is important for cell and organism viability. Here, we describe a rearrangement of Homo sapiens chromosome 17 (HSA17) that has placed alpha satellite DNA next to euchromatin. We show that on this mutant chromosome, CENP-A chromatin has spread from the alpha satellite into the short arm of HSA17, establishing a ∼700 kb hybrid centromeric domain that spans both repetitive and unique sequences and changes the expression of at least one gene over which it spreads. Our results illustrate the plasticity of human centromeric chromatin and suggest that heterochromatin normally constrains CENP-A chromatin onto alpha satellite DNA. This work highlights that chromosome rearrangements, particularly those that remove the pericentromere, create opportunities for centromeric nucleosomes to move into non-traditional genomic locations, potentially changing the surrounding chromatin environment and altering gene expression.

Case-Only Survival Analysis Reveals Unique Effects of Genotype, Sex, and Coronary Disease Severity on Survivorship

Survival bias may unduly impact genetic association with complex diseases; gene-specific survival effects may further complicate such investigations. Coronary artery disease (CAD) is a complex phenotype for which little is understood about gene-specific survival effects; yet, such information can offer insight into refining genetic associations, improving replications, and can provide candidate genes for both mortality risk and improved survivorship in CAD. Building on our previous work, the purpose of this current study was to: evaluate LSAMP SNP-specific hazards for all-cause mortality post-catheterization in a larger cohort of our CAD cases; and, perform additional replication in an independent dataset. We examined two LSAMP SNPs—rs1462845 and rs6788787—using CAD case-only Cox proportional hazards regression for additive genetic effects, censored on time-to-all-cause mortality or last follow-up among Caucasian subjects from the Catheterization Genetics Study (CATHGEN; n = 2,224) and the Intermountain Heart Collaborative Study (IMHC; n = 3,008). Only after controlling for age, sex, body mass index, histories of smoking, type 2 diabetes, hyperlipidemia and hypertension (HR = 1.11, 95%CI = 1.01–1.22, p = 0.032), rs1462845 conferred significantly increased hazards of all-cause mortality among CAD cases. Even after controlling for multiple covariates, but in only the primary cohort, rs6788787 conferred significantly improved survival (HR = 0.80, 95% CI = 0.69–0.92, p = 0.002). Post-hoc analyses further stratifying by sex and disease severity revealed replicated effects for rs1462845: even after adjusting for aforementioned covariates and coronary interventional procedures, males with severe burden of CAD had significantly amplified hazards of death with the minor variant of rs1462845 in both cohorts (HR = 1.29, 95% CI = 1.08–1.55, p = 0.00456; replication HR = 1.25, 95% CI = 1.05–1.49, p = 0.013). Kaplan-Meier curves revealed unique cohort-specific genotype effects on survival. Additional analyses demonstrated that the homozygous risk genotype (‘A/A’) fully explained the increased hazard in both cohorts. None of the post-hoc analyses in control subjects were significant for any model. This suggests that genetic effects of rs1462845 on survival are unique to CAD presence. This represents formal, replicated evidence of genetic contribution of rs1462845 to increased risk for all-cause mortality; the contribution is unique to CAD case status and specific to males with severe burden of CAD.

Interaction Between the FOXO1A-209 Genotype and Tea Drinking Is Significantly Associated with Reduced Mortality at Advanced Ages

On the basis of the genotypic/phenotypic data from Chinese Longitudinal Healthy Longevity Survey (CLHLS) and Cox proportional hazard model, the present study demonstrates that interactions between carrying FOXO1A-209 genotypes and tea drinking are significantly associated with lower risk of mortality at advanced ages. Such a significant association is replicated in two independent Han Chinese CLHLS cohorts (p = 0.028-0.048 in the discovery and replication cohorts, and p = 0.003-0.016 in the combined dataset). We found the associations between tea drinking and reduced mortality are much stronger among carriers of the FOXO1A-209 genotype compared to non-carriers, and drinking tea is associated with a reversal of the negative effects of carrying FOXO1A-209 minor alleles, that is, from a substantially increased mortality risk to substantially reduced mortality risk at advanced ages. The impacts are considerably stronger among those who carry two copies of the FOXO1A minor allele than those who carry one copy. On the basis of previously reported experiments on human cell models concerning FOXO1A-by-tea-compounds interactions, we speculate that results in the present study indicate that tea drinking may inhibit FOXO1A-209 gene expression and its biological functions, which reduces the negative impacts of FOXO1A-209 gene on longevity (as reported in the literature) and offers protection against mortality risk at oldest-old ages. Our empirical findings imply that the health outcomes of particular nutritional interventions, including tea drinking, may, in part, depend upon individual genetic profiles, and the research on the effects of nutrigenomics interactions could potentially be useful for rejuvenation therapies in the clinic or associated healthy aging intervention programs.

Metabolomic Quantitative Trait Loci (mQTL) Mapping Implicates the Ubiquitin Proteasome System in Cardiovascular Disease Pathogenesis

Levels of certain circulating short-chain dicarboxylacylcarnitine (SCDA), long-chain dicarboxylacylcarnitine (LCDA) and medium chain acylcarnitine (MCA) metabolites are heritable and predict cardiovascular disease (CVD) events. Little is known about the biological pathways that influence levels of most of these metabolites. Here, we analyzed genetics, epigenetics, and transcriptomics with metabolomics in samples from a large CVD cohort to identify novel genetic markers for CVD and to better understand the role of metabolites in CVD pathogenesis. Using genomewide association in the CATHGEN cohort (N = 1490), we observed associations of several metabolites with genetic loci. Our strongest findings were for SCDA metabolite levels with variants in genes that regulate components of endoplasmic reticulum (ER) stress (USP3, HERC1, STIM1, SEL1L, FBXO25, SUGT1) These findings were validated in a second cohort of CATHGEN subjects (N = 2022, combined p = 8.4x10^-6–2.3x10^-10). Importantly, variants in these genes independently predicted CVD events. Association of genomewide methylation profiles with SCDA metabolites identified two ER stress genes as differentially methylated (BRSK2 and HOOK2). Expression quantitative trait loci (eQTL) pathway analyses driven by gene variants and SCDA metabolites corroborated perturbations in ER stress and highlighted the ubiquitin proteasome system (UPS) arm. Moreover, culture of human kidney cells in the presence of levels of fatty acids found in individuals with cardiometabolic disease, induced accumulation of SCDA metabolites in parallel with increases in the ER stress marker BiP. Thus, our integrative strategy implicates the UPS arm of the ER stress pathway in CVD pathogenesis, and identifies novel genetic loci associated with CVD event risk.

Cardiovascular disease is a strongly heritable trait. Despite application of the latest genomic technologies, the genetic architecture of disease risk remains poorly defined, and mechanisms underlying this susceptibility are incompletely understood. In this study, we performed genome-wide mapping of heart disease-related metabolites measured in the blood as the genetic traits of interest (instead of the disease itself), in a large cohort of 3512 patients at risk of heart disease from the CATHGEN study. Our goal was to discover new cardiovascular disease genes and thereby mechanisms of disease pathogenesis by understanding the genes that regulate levels of these metabolites. These analyses identified novel genetic variants associated with metabolite levels and with cardiovascular disease itself. Importantly, by utilizing an unbiased systems-based approach integrating genetics, gene expression, epigenetics and metabolomics, we uncovered a novel pathway of heart disease pathogenesis, that of endoplasmic reticulum (ER) stress, represented by elevated levels of circulating short-chain dicarboxylacylcarnitine (SCDA) metabolites.

Association of Roadway Proximity with Fasting Plasma Glucose and Metabolic Risk Factors for Cardiovascular Disease in a Cross-Sectional Study of Cardiac Catheterization Patients

BACKGROUND

The relationship between traffic-related air pollution (TRAP) and risk factors for cardiovascular disease needs to be better understood in order to address the adverse impact of air pollution on human health.

OBJECTIVE

We examined associations between roadway proximity and traffic exposure zones, as markers of TRAP exposure, and metabolic biomarkers for cardiovascular disease risk in a cohort of patients undergoing cardiac catheterization.

METHODS

We performed a cross-sectional study of 2,124 individuals residing in North Carolina (USA). Roadway proximity was assessed via distance to primary and secondary roadways, and we used residence in traffic exposure zones (TEZs) as a proxy for TRAP. Two categories of metabolic outcomes were studied: measures associated with glucose control, and measures associated with lipid metabolism. Statistical models were adjusted for race, sex, smoking, body mass index, and socioeconomic status (SES).

RESULTS

An interquartile-range (990 m) decrease in distance to roadways was associated with higher fasting plasma glucose (β = 2.17 mg/dL; 95% CI: -0.24, 4.59), and the association appeared to be limited to women (β = 5.16 mg/dL; 95% CI: 1.48, 8.84 compared with β = 0.14 mg/dL; 95% CI: -3.04, 3.33 in men). Residence in TEZ 5 (high-speed traffic) and TEZ 6 (stop-and-go traffic), the two traffic zones assumed to have the highest levels of TRAP, was positively associated with high-density lipoprotein cholesterol (HDL-C; β = 8.36; 95% CI: -0.15, 16.9 and β = 5.98; 95% CI: -3.96, 15.9, for TEZ 5 and 6, respectively).

CONCLUSION

Proxy measures of TRAP exposure were associated with intermediate metabolic traits associated with cardiovascular disease, including fasting plasma glucose and possibly HDL-C.

Evidence for fumonisin inhibition of ceramide synthase in humans consuming maize-based foods and living in high exposure communities in Guatemala

SCOPE

Fumonisin (FB) occurs in maize and is an inhibitor of ceramide synthase (CerS). We determined the urinary FB1 (UFB1 ) and sphingoid base 1-phosphate levels in blood from women consuming maize in high and low FB exposure communities in Guatemala.

METHODS AND RESULTS

FB1 intake was estimated using the UFB1 . Sphinganine 1-phosphate (Sa 1-P), sphingosine 1-phosphate (So 1-P), and the Sa 1-P/So 1-P ratio were determined in blood spots collected on absorbent paper at the same time as urine collection. In the first study, blood spots and urine were collected every 3 months (March 2011 to February 2012) from women living in low (Chimaltenango and Escuintla) and high (Jutiapa) FB exposure communities (1240 total recruits). The UFB1 , Sa 1-P/So 1-P ratio, and Sa 1-P/mL in blood spots were significantly higher in the high FB1 intake community compared to the low FB1 intake communities. The results were confirmed in a follow-up study (February 2013) involving 299 women living in low (Sacatepéquez) and high (Santa Rosa and Chiquimula) FB exposure communities.

CONCLUSIONS

High levels of FB1 intake are correlated with changes in Sa 1-P and the Sa 1-P/So 1-P ratio in human blood in a manner consistent with FB1 inhibition of CerS.