Search

Longitudinal changes in gene expression during islet autoimmunity (IA) may provide insight into biological processes that explain progression to type 1 diabetes (T1D). We identified individuals from Diabetes Autoimmunity Study in the Young (DAISY) who developed IA, autoantibodies present on two or more visits. Illumina’s NovaSeq 6000 was used to quantify gene expression in whole blood. With linear mixed models we tested for changes in expression after IA that differed across individuals who progressed to T1D (progressors) (n = 25), reverted to an autoantibody-negative stage (reverters) (n = 47), or maintained IA positivity but did not develop T1D (maintainers) (n = 66). Weighted gene coexpression network analysis was used to identify coexpression modules. Gene Ontology pathway analysis of the top 150 differentially expressed genes (nominal P < 0.01) identified significantly enriched pathways including leukocyte activation involved in immune response, innate immune response, and regulation of immune response. We identified a module of 14 coexpressed genes with roles in the innate immunity. The hub gene, LTF, is known to have immunomodulatory properties. Another gene within the module, CAMP, is potentially relevant based on its role in promoting β-cell survival in a murine model. Overall, results provide evidence of alterations in expression of innate immune genes prior to onset of T1D.

Methods, Supp Figure Legends, Supp Table Legend

Elevated tumor expression of the cell surface GPI-linked CD24 protein signals poor patient prognosis in many tumor types. However, some cancer cells selected to be negative for surface CD24 (surCD24−) still retain aggressive phenotypes in vitro and in vivo. Here, we resolve this apparent paradox with the discovery of biologically active, nuclear CD24 (nucCD24) and finding that its levels are unchanged in surCD24− cells. Using the complementary techniques of biochemical cellular fractionation and immunofluorescence, we demonstrate a signal for CD24 in the nucleus in cells from various histologic types of cancer. Nuclear-specific expression of CD24 (NLS-CD24) increased anchorage-independent growth in vitro and tumor formation in vivo. Immunohistochemistry of patient tumor samples revealed the presence of nucCD24, whose signal intensity correlated positively with the presence of metastatic disease. Analysis of gene expression between cells expressing CD24 and NLS-CD24 revealed a unique nucCD24 transcriptional signature. The median score derived from this signature was able to stratify overall survival in four patient datasets from bladder cancer and five patient datasets from colorectal cancer. Patients with high scores (more nucCD24-like) had reduced survival. These findings define a novel and functionally important intracellular location of CD24; they explain why surCD24− cells can remain aggressive, and they highlight the need to consider nucCD24 in both fundamental research and therapeutic development. Cancer Res; 77(18); 4858–67. ©2017 AACR.

Additional cell lines and GFP-CD24 confirm effects of glutaraldehyde to fixation protocol and the nuclear CD24 signal.

List of 304 genes which were found to have statistically significant differential expression between UMUC13 cells expressing NLS-CD24 and Wt-CD24.

Examples of masking used to quantify total cellular and nuclear CD24.

Additional cancer patient datasets showing the nucCD24 gene signature correlates negatively with survival

Isolated surCD24- cells maintain loss of surCD24.

Inhibiting nuclear export with Leptomycin B leads to nuclear accumulation of CD24 protein.

Over 6.37 million people have died from COVID-19 worldwide, but factors influencing COVID-19-related mortality remain understudied. We aimed to describe and identify risk factors for COVID-19 mortality in the Colorado Center for Personalized Medicine (CCPM) Biobank using integrated data sources, including Electronic Health Records (EHRs). We calculated cause-specific mortality and case-fatality rates for COVID-19 and common pre-existing health conditions defined by diagnostic phecodes and encounters in EHRs. We performed multivariable logistic regression analyses of the association between each pre-existing condition and COVID-19 mortality. Of the 155,859 Biobank participants enrolled as of July 2022, 20,797 had been diagnosed with COVID-19. Of 5334 Biobank participants who had died, 190 were attributed to COVID-19. The case-fatality rate was 0.91% and the COVID-19 mortality rate was 122 per 100,000 persons. The odds of dying from COVID-19 were significantly increased among older men, and those with 14 of the 61 pre-existing conditions tested, including hypertensive chronic kidney disease (OR: 10.14, 95% CI: 5.48, 19.16) and type 2 diabetes with renal manifestations (OR: 5.59, 95% CI: 3.42, 8.97). Male patients who are older and have pre-existing kidney diseases may be at higher risk for death from COVID-19 and may require special care.

High and Low Activity strains of mice were bidirectionally selected for differences in open-field activity (DeFries et al., 1978, Behavior Genetics, 8: 3-13) and subsequently inbred to use as a genetic model for studying anxiety-like behaviors (Booher et al., 2021, Genes, Brain and Behavior, 20: e12730). Hippocampal RNA-sequencing of the High and Low Activity mice identified 3901 differentially expressed protein-coding genes, with both sex-dependent and sex-independent effects. Functional enrichment analysis (PANTHER) highlighted 15 gene ontology terms, which allowed us to create a narrow list of 264 top candidate genes. Of the top candidate genes, 46 encoded four Complexes (I, II, IV and V) and two electron carriers (cytochrome c and ubiquinone) of the mitochondrial oxidative phosphorylation process. The most striking results were in the female high anxiety, Low Activity mice, where 39/46 genes relating to oxidative phosphorylation were upregulated. In addition, comparison of our top candidate genes with two previously curated High and Low Activity gene lists highlight 24 overlapping genes, where Ndufa13, which encodes the supernumerary subunit A13 of complex I, was the only gene to be included in all three lists. Mitochondrial dysfunction has recently been implicated as both a cause and effect of anxiety-related disorders and thus should be further explored as a possible novel pharmaceutical treatment for anxiety disorders.

The Hybrid Rat Diversity Panel (HRDP) is a stable and well-characterized set of more than 90 inbred rat strains that can be leveraged for systems genetics approaches to understanding the genetic and genomic variation associated with complex disease. The HRDP exhibits substantial between-strain diversity while retaining substantial within-strain isogenicity, allowing for the precise mapping of genetic variation associated with complex phenotypes and providing statistical power to identify associated variants. In order to robustly identify associated genetic variants, it is important to account for the population structure induced by inbreeding. To this end, we investigate the performance of various plausible approaches towards modeling quantitative traits in the HRDP and quantify their operating characteristics. Towards facilitating study planning and design, we conduct extensive simulations to investigate the power of genetic association analyses in the HRDP, and characterize the impressive attained power.

Omics studies frequently use samples collected during cohort studies. Conditioning on sample availability can cause selection bias if sample availability is nonrandom. Inverse probability weighting (IPW) is purported to reduce this bias. We evaluated IPW in an epigenome-wide analysis testing the association between DNA methylation (261,435 probes) and age in healthy adolescent subjects (n = 114). We simulated age and sex to be correlated with sample selection and then evaluated four conditions: complete population/no selection bias (all subjects), naïve selection bias (no adjustment), and IPW selection bias (selection bias with IPW adjustment). Assuming the complete population condition represented the "truth," we compared each condition to the complete population condition. Bias or difference in associations between age and methylation was reduced in the IPW condition versus the naïve condition. However, genomic inflation and type 1 error were higher in the IPW condition relative to the naïve condition. Postadjustment using bacon, type 1 error and inflation were similar across all conditions. Power was higher under the IPW condition compared with the naïve condition before and after inflation adjustment. IPW methods can reduce bias in genome-wide analyses. Genomic inflation is a potential concern that can be minimized using methods that adjust for inflation.

The androgen receptor (AR) is important in the development of both experimental and human bladder cancer. However, the role of AR in bladder cancer growth and progression is less clear, with literature indicating that more advanced stage and grade disease are associated with reduced AR expression. To determine the mechanisms underlying these relationships, we profiled AR-expressing human bladder cancer cells by AR chromatin immunoprecipitation sequencing and complementary transcriptomic approaches in response to in vitro stimulation by the synthetic androgen R1881. In vivo functional genomics consisting of pooled shRNA or pooled open reading frame libraries was employed to evaluate 97 genes that recapitulate the direction of expression associated with androgen stimulation. Interestingly, we identified CD44, the receptor for hyaluronic acid, a potent biomarker and driver of progressive disease in multiple tumor types, as significantly associated with androgen stimulation. CRISPR-based mutagenesis of androgen response elements associated with CD44 identified a novel silencer element leading to the direct transcriptional repression of CD44 expression. In human patients with bladder cancer, tumor AR and CD44 mRNA and protein expression were inversely correlated, suggesting a clinically relevant AR–CD44 axis. Collectively, our work describes a novel mechanism partly explaining the inverse relationship between AR and bladder cancer tumor progression and suggests that AR and CD44 expression may be useful for prognostication and therapeutic selection in primary bladder cancer. Significance: This study describes novel AREs that suppress CD44 and an expected inverse correlation of AR-CD44 expression observed in human bladder tumors.

Longitudinal changes in gene expression during islet autoimmunity (IA) may provide insight into biological processes that explain progression to type 1 diabetes (T1D). We identified individuals from DAISY who developed IA, autoantibodies present on two or more visits. Illumina’s NovaSEQ 600TM was used to quantify gene expression in whole blood. Linear mixed models tested for changes in expression after IA that differed across individuals who progressed to T1D (progressors, n=25), reverted to an autoantibody negative stage (reverters, n=47), or maintained IA positivity but did not develop T1D (maintainers, n=66). Weighted gene co-expression network analysis was used to identify co-expression modules. Gene Ontology pathway analysis of the top 150 differentially expressed genes (nominal pLTF, is known to have immunomodulatory properties. Another gene within the module, CAMP, is potentially relevant based on its role in promoting beta cell survival in a murine model. Overall, results provide evidence of alterations in expression of innate immune genes prior to onset of T1D.

Lower blood 25 (OH) D concentration is associated with increased risk of islet autoimmunity (IA) in some cohorts; it is not known if 25 (OH) D metabolites are important after the onset of IA. We investigated novel vitamin D metabolites in development of IA and progression from IA to T1D in the Diabetes Autoimmunity Study in the Young (DAISY) . IA was defined as two or more consecutive visits positive for ≥ 1 autoantibody (IAA, IA-2A, GAD, ZnT8) . In a nested case-control study, we examined plasma vitamin D metabolites at a pre-seroconversion visit in 114 IA cases and 116 controls. We then followed 144 IA-positive children, of whom 46 developed T1D during a mean follow-up of 8.4 years. Plasma was quantified for 25 (OH) D3, 3-epi-25 (OH) D3, 25 (OH) D2, 24,25 (OH) 2D3, and 1,25 (OH) 2D3, using targeted LC-MS/MS. Logistic regression was used to examine pre-seroconversion metabolites and IA, and Cox regression was used to examine metabolites at seroconversion and progression to T1D. Both models adjusted for age, ethnicity, T1D family history and HLA-DR3/4. Higher pre-seroconversion level of 24,25 (OH) 2D3, an inactive catabolite of 25 (OH) D3, was associated with increased odds of IA (OR: 1.35; 95%CI: 1.03,1.78 per SD of 24,25 (OH) 2D3; p=0.031) . The association between 3-epi-25 (OH) D3 at seroconversion and progression from IA to T1D was modified by age at seroconversion. Higher 3-epi-25 (OH) D3 level was more protective against progression to T1D among children who seroconverted at an older age (HR: 0.47; 95%CI: 0.26, 0.85 for children at the 75th percentile of age at seroconversion (ie years)) compared to children who seroconverted at a younger age (HR: 0.88; 95%CI: 0.63, 1.24 for children at the 25th percentile (ie 3 years)) . Neither IA nor progression to T1D was associated with 25 (OH) D3, 25 (OH) D2 or 1,25 (OH) 2D3.Increased 24,25 (OH) 2D3 may play a role in IA development. The age-related heterogeneity of the relationship between 3-epi-25 (OH) D3 and T1D progression requires further investigation. Disclosure S.Kim: None. K.Kechris: None. O.Fiehn: None. M.Rewers: Consultant; Janssen Research & Development, LLC, Medscape, Provention Bio, Inc., Research Support; Dexcom, Inc., JDRF, Roche Diagnostics USA. J.M.Norris: None. N.L.Zwick: None. P.M.Carry: None. B.C.Defelice: None. R.K.Johnson: None. F.Dong: None. T.Buckner: None. L.A.Vanderlinden: None. B.I.Frohnert: Advisory Panel; Provention Bio, Inc. Funding National Institutes of Health (R01-DK104351, R01-DK32493, RI142483)

The R7 and R8 photoreceptor cells of the Drosophila compound eye mediate color vision. Throughout the majority of the eye, these cells occur in two principal types of ommatidia. Approximately 35% of ommatidia are of the pale type and express Rh3 in R7 cells and Rh5 in R8 cells. The remaining 65% are of the yellow type and express Rh4 in R7 cells and Rh6 in R8 cells. The specification of an R8 cell in a pale or yellow ommatidium depends on the fate of the adjacent R7 cell. However, pale and yellow R7 cells are specified by a stochastic process that requires the genes spineless, tango and klumpfuss. To identify additional genes involved in this process we performed genetic screens using a collection of 480 P{EP} transposon insertion strains. We identified genes in gain of function and loss of function screens that significantly altered the percentage of Rh3 expressing R7 cells (Rh3%) from wild-type. 36 strains resulted in altered Rh3% in the gain of function screen where the P{EP} insertion strains were crossed to a sevEP-GAL4 driver line. 53 strains resulted in altered Rh3% in the heterozygous loss of function screen. 4 strains showed effects that differed between the two screens, suggesting that the effect found in the gain of function screen was either larger than, or potentially masked by, the P{EP} insertion alone. Analyses of homozygotes validated many of the candidates identified. These results suggest that R7 cell fate specification is sensitive to perturbations in mRNA transcription, splicing and localization, growth inhibition, post-translational protein modification, cleavage and secretion, hedgehog signaling, ubiquitin protease activity, GTPase activation, actin and cytoskeletal regulation, and Ser/Thr kinase activity, among other diverse signaling and cell biological processes.

Objective Our aim was to elucidate the role of diet in type 1 diabetes (T1D) by examining combinations of nutrient intake in the progression from islet autoimmunity (IA) to T1D. Methods We measured 2457 metabolites and dietary intake at the time of seroconversion in 132 IA-positive children in the prospective Diabetes Autoimmunity Study in the Young. IA was defined as the first of two consecutive visits positive for at least one autoantibody (insulin, GAD, IA-2, or ZnT8). By December 2018, 40 children progressed to T1D. Intakes of 38 nutrients were estimated from semiquantitative food frequency questionnaires. We tested the association of each metabolite with progression to T1D using multivariable Cox regression. Nutrient patterns that best explained variation in candidate metabolites were identified using reduced rank regression (RRR), and their association with progression to T1D was tested using Cox regression adjusting for age at seroconversion and high-risk HLA genotype. Results In stepwise selection, 22 nutrients significantly predicted at least two of the 13 most significant metabolites associated with progression to T1D, and were included in RRR. A nutrient pattern corresponding to intake lower in linoleic acid, niacin, and riboflavin, and higher in total sugars, explained 18% of metabolite variability. Children scoring higher on this metabolite-related nutrient pattern at seroconversion had increased risk for progressing to T1D (HR = 3.17, 95%CI = 1.42-7.05). Conclusions Combinations of nutrient intake reflecting candidate metabolites are associated with increased risk of T1D, and may help focus dietary prevention efforts.

Background Oxylipins are formed from oxidation of omega-6 (n6) and omega-3 (n3) fatty acids (FAs). Evidence for inflammatory effects comes mostly from adults. Methods Oxylipins from n6 FA (27 n6-oxylipins) and n3 FA (12 n3-oxylipins) were measured through ultra-high-performance liquid chromatography-mass spectrometry (LC-MS/MS) in plasma from 111 children at risk of type 1 diabetes (age 1-17 years) studied longitudinally. Oxylipin precursor FAs (arachidonic acid, linoleic acid, alpha-linolenic acid, docosahexaenoic acid, eicosapentaenoic acid) were measured in red blood cell (RBC) membrane and plasma. Precursor FAs dietary intake was measured through food frequency questionnaire and environmental tobacco smoke (ETS) through questionnaires. Linear mixed models were used to test oxylipins with predictors. Results Age associated with 15 n6- and 6 n3-oxylipins; race/ethnicity associated with 3 n6- and 1 n3-oxylipins; sex associated with 2 n6-oxylipins. ETS associated with lipoxin-A4. Oxylipins associated with precursor FAs in plasma more often than RBC. RBC levels and dietary intake of precursor FAs more consistently associated with n3-oxylipins than with n6-oxylipins. Conclusions In healthy children, oxylipin levels change with age. Oxylipins associated with precursor FAs more often in plasma than RBC or diet, suggesting that inflammatory regulation leading to FA release into plasma may also be a determinant of oxylipin generation. Impact This is the first study to examine predictors of oxylipins in healthy children at risk of type 1 diabetes.In healthy children at risk of type 1 diabetes, many oxylipins change with age, and most oxylipins do not differ by sex or race/ethnicity.Environmental tobacco smoke exposure was associated with the presence of lipoxin A4.Omega-6- and omega-3-related oxylipin levels were consistently associated with their respective precursor fatty acid levels measured in the plasma.Proportionally more omega-3 compared to omega-6 oxylipins were associated with dietary intake and red blood cell membrane levels of the respective precursor fatty acid.

OBJECTIVE Mechanisms underlying the role of non-human leukocyte antigen (HLA) genetic risk variants in type 1 diabetes (T1D) are poorly understood. We aimed to test the association between methylation and non-HLA genetic risk. METHODS We conducted a methylation quantitative trait loci (mQTL) analysis in a nested case-control study from the Dietary Autoimmunity Study in the Young. Controls (n = 83) were frequency-matched to T1D cases (n = 83) based on age, race/ethnicity, and sample availability. We evaluated 13 non-HLA genetic markers known be associated with T1D. Genome-wide methylation profiling was performed on peripheral blood samples collected prior to T1D using the Illumina 450 K (discovery set) and infinium methylation EPIC beadchip (EPIC validation) platforms. Linear regression models, adjusting for age and sex, were used to test to each single nucleotide polymorphism (SNP) -probe combination. Logistic regression models were used to test the association between T1D and methylation levels among probes with a significant mQTL. A meta-analysis was used to combine odds ratios from the two platforms. RESULTS We identified 10 SNP-methylation probe pairs (false discovery rate (FDR) adjusted P

Oxylipins, pro-inflammatory and pro-resolving lipid mediators, are associated with the risk of type 1 diabetes (T1D) and may be influenced by diet. This study aimed to develop a nutrient pattern related to oxylipin profiles and test their associations with the risk of T1D among youth. The nutrient patterns were developed with a reduced rank regression in a nested case-control study (n = 335) within the Diabetes Autoimmunity Study in the Young (DAISY), a longitudinal cohort of children at risk of T1D. The oxylipin profiles (adjusted for genetic predictors) were the response variables. The nutrient patterns were tested in the case-control study (n = 69 T1D cases, 69 controls), then validated in the DAISY cohort using a joint Cox proportional hazards model (n = 1933, including 81 T1D cases). The first nutrient pattern (NP1) was characterized by low beta cryptoxanthin, flavanone, vitamin C, total sugars and iron, and high lycopene, anthocyanidins, linoleic acid and sodium. After adjusting for T1D family history, the HLA genotype, sex and race/ethnicity, NP1 was associated with a lower risk of T1D in the nested case-control study (OR: 0.44, p = 0.0126). NP1 was not associated with the risk of T1D (HR: 0.54, p-value = 0.1829) in the full DAISY cohort. Future studies are needed to confirm the nested case-control findings and investigate the modifiable factors for oxylipins.

The transmembrane receptor Frizzled 9 (FZD9) is important for fetal neurologic and bone development through both canonical and non-canonical WNT/FZD signaling. In the adult lung, however, Fzd9 helps to maintain a normal epithelium by signaling through peroxisome proliferator activated receptor γ (PPARγ). The effect of FZD9 loss on normal lung epithelial cells and regulators of its expression in the lung are unknown. We knocked down FZD9 in human bronchial epithelial cell (HBEC) lines and found that downstream EMT targets and PPARγ activity are altered. We used a FZD9-/- mouse in the urethane lung adenocarcinoma model and found FZD9-/- adenomas had more proliferation, increased EMT signaling, decreased activation of PPARγ, increased expression of lung cancer associated genes, increased transformed growth, and increased potential for invasive behavior. We identified PPARγ as a transcriptional regulator of FZD9. We also demonstrated that extended cigarette smoke exposure in HBEC leads to decreased FZD9 expression, decreased activation of PPARγ, and increased transformed growth, and found that higher exposure to cigarette smoke in human lungs leads to decreased FZD9 expression. These results provide evidence for the role of FZD9 in lung epithelial maintenance and in smoking related malignant transformation. We identified the first transcriptional regulator of FZD9 in the lung and found FZD9 negative lesions are more dangerous. Loss of FZD9 creates a permissive environment for development of premalignant lung lesions, making it a potential target for intervention.

We assessed associations between infant diet (e.g., breastfeeding and introduction to solid foods) and DNA methylation in infancy and childhood. We measured DNA methylation in peripheral blood collected in infancy (9–15 months of age) in 243 children, and in a subset of 50 children, we also measured methylation in childhood (6–9 years of age) to examine persistence, and at birth (in cord blood) to examine temporality. We performed multivariable linear regression of infant diet on the outcome of methylation using epigenome-wide and candidate site approaches. We identified six novel CpG sites associated with breastfeeding duration using an EWAS approach. One differentially methylated site presented directionally consistent associations with breastfeeding (cg00574958, CPT1A) in infancy and childhood but not at birth. Two differentially methylated sites in infancy (cg19693031, TXNIP, cg23307264, KHSRP) were associated with breastfeeding and were not present at birth, however, these associations did not persist into childhood. Associations between infant diet and methylation in infancy at three sites (cg22369607, AP001525.1, cg2409200, TBCD, cg27173510, PGBD5) were also present at birth, suggesting the influence of exposures other than infant diet. Infant diet exposures are associated with persistent methylation differences in CPT1A, which may be one mechanism behind infant diet’s long-term health effects.

BackgroundCardiac hypertrophy can be considered a maladaptive response which in many cases results in heart failure and sudden death. Genetic predisposition to cardiac hypertrophy is ill-defined, but current research has given credence to a role of non-protein-coding RNAs in the development of cardiac hypertrophy.ResultsWe used microarrays, RNA-Seq and quantitative genetics with the spontaneously hypertensive, SHR (SHR/OlaIpcv) rat strain, the normotensive Brown Norway (BN-Lx/Cub) rat strain, and a recombinant inbred (RI) panel of rats (HXB/BXH) derived from these strains to examine the areas of the genome associated with cardiac hypertrophy. We identified circular (circ) RNAs coded within these areas. Of the 122 differentially expressed circRNAs in left ventricle of SHR and BN-Lx rats, three were transcribed from areas corresponding to the QTLs we identified for cardiac hypertrophy using the HXB/BXH rat panel. We then identified microRNAs which are “sponged” by the circ RNAs and the mRNAs which are destabilized by the microRNAs. 265 miRNAs could be identified as targets for the three circRNAs. We focused on the four miRNAs that were also differentially expressed between SHR and BN-Lx rats. One of the miRNAs (Mir-210-5p) was a target of the differentially expressed circ H2afy and circ H2afy was located within the QTL on Chr 17 (genome wide p-value = 0.011). Mir-210-5p has a binding site on the 3’ UTR of DR6, which is differentially expressed and in turn controls the expression level of NFκB. NFκB is an important component of the oxidative stress response leading to hypertrophy.ConclusionsOur work identified a novel circRNA that may be the mediator of a cascade of events that influence a cardiac hypertrophy phenotype. The circRNA and miRNA which we identified may become useful as markers of the risk for cardiac hypertrophy.

AIMS/HYPOTHESIS: Oxylipins are lipid mediators derived from polyunsaturated fatty acids. Some oxylipins are proinflammatory (e.g. those derived from arachidonic acid [ARA]), others are pro-resolving of inflammation (e.g. those derived from α-linolenic acid [ALA], docosahexaenoic acid [DHA] and eicosapentaenoic acid [EPA]) and others may be both (e.g. those derived from linoleic acid [LA]). The goal of this study was to examine whether oxylipins are associated with incident type 1 diabetes. METHODS: We conducted a nested case–control analysis in the Diabetes Autoimmunity Study in the Young (DAISY), a prospective cohort study of children at risk of type 1 diabetes. Plasma levels of 14 ARA-derived oxylipins, ten LA-derived oxylipins, six ALA-derived oxylipins, four DHA-derived oxylipins and two EPA-related oxylipins were measured by ultra-HPLC–MS/MS at multiple timepoints related to autoantibody seroconversion in 72 type 1 diabetes cases and 71 control participants, which were frequency matched on age at autoantibody seroconversion (of the case), ethnicity and sample availability. Linear mixed models were used to obtain an age-adjusted mean of each oxylipin prior to type 1 diabetes. Age-adjusted mean oxylipins were tested for association with type 1 diabetes using logistic regression, adjusting for the high risk HLA genotype HLA-DR3/4,DQB1*0302. We also performed principal component analysis of the oxylipins and tested principal components (PCs) for association with type 1 diabetes. Finally, to investigate potential critical timepoints, we examined the association of oxylipins measured before and after autoantibody seroconversion (of the cases) using PCs of the oxylipins at those visits. RESULTS: The ARA-related oxylipin 5-HETE was associated with increased type 1 diabetes risk. Five LA-related oxylipins, two ALA-related oxylipins and one DHA-related oxylipin were associated with decreased type 1 diabetes risk. A profile of elevated LA- and ALA-related oxylipins (PC1) was associated with decreased type 1 diabetes risk (OR 0.61; 95% CI 0.40, 0.94). A profile of elevated ARA-related oxylipins (PC2) was associated with increased diabetes risk (OR 1.53; 95% CI 1.03, 2.29). A critical timepoint analysis showed type 1 diabetes was associated with a high ARA-related oxylipin profile at post-autoantibody-seroconversion but not pre-seroconversion. CONCLUSIONS/INTERPRETATION: The protective association of higher LA- and ALA-related oxylipins demonstrates the importance of both inflammation promotion and resolution in type 1 diabetes. Proinflammatory ARA-related oxylipins may play an important role once the autoimmune process has begun.

Epigenetic mechanisms may contribute to idiopathic scoliosis (IS). We identified 8 monozygotic twin pairs with IS, 6 discordant (Cobb angle difference > 10°) and 2 concordant (Cobb angle difference ≤ 2°). Genome-wide methylation in blood was measured with the Infinium HumanMethylation EPIC Beadchip. We tested for differences in methylation and methylation variability between discordant twins and tested the association between methylation and curve severity in all twins. Differentially methylated region (DMR) analyses identified gene promoter regions. Methylation at cg12959265 (chr. 7 DPY19L1) was less variable in cases (false discovery rate (FDR) = 0.0791). We identified four probes (false discovery rate, FDR < 0.10); cg02477677 (chr. 17, RARA gene), cg12922161 (chr. 2 LOC150622 gene), cg08826461 (chr. 2), and cg16382077 (chr. 7) associated with curve severity. We identified 57 DMRs where hyper- or hypo-methylation was consistent across the region and 28 DMRs with a consistent association with curve severity. Among DMRs, 21 were correlated with bone methylation. Prioritization of regions based on methylation concordance in bone identified promoter regions for WNT10A (WNT signaling), NPY (regulator of bone and energy homeostasis), and others predicted to be relevant for bone formation/remodeling. These regions may aid in understanding the complex interplay between genetics, environment, and IS.

Reversion of islet autoimmunity (IA) may point to mechanisms that prevent IA progression. We followed 199 individuals who developed IA during the Diabetes Autoimmunity Study in the Young. Untargeted metabolomics was performed in serum samples following IA. Cox proportional hazards models were used to test whether the metabolites (2,487) predicted IA reversion: two or more consecutive visits negative for all autoantibodies. We conducted a principal components analysis (PCA) of the top metabolites; |hazard ratio (HR) >1.25| and nominal P < 0.01. Phosphatidylcholine (16:0_18:1(9Z)) was the strongest individual metabolite (HR per 1 SD 2.16, false discovery rate (FDR)-adjusted P = 0.0037). Enrichment analysis identified four clusters (FDR P < 0.10) characterized by an overabundance of sphingomyelin (d40:0), phosphatidylcholine (16:0_18:1(9Z)), phosphatidylcholine (30:0), and l-decanoylcarnitine. Overall, 63 metabolites met the criteria for inclusion in the PCA. PC1 (HR 1.4, P < 0.0001), PC2 (HR 0.85, P = 0.0185), and PC4 (HR 1.28, P = 0.0103) were associated with IA reversion. Given the potential influence of diet on the metabolome, we investigated whether nutrients were correlated with PCs. We identified 20 nutrients that were correlated with the PCs (P < 0.05). Total sugar intake was the top nutrient. Overall, we identified an association between phosphatidylcholine, sphingomyelin, and carnitine levels and reversion of IA.

Low molecular weight polycyclic aromatic hydrocarbons (LMW PAHs

Environmental factors including viruses, diet, and the metabolome have been linked with the appearance of islet autoimmunity (IA) that precedes development of type 1 diabetes (T1D). We measured global DNA methylation (DNAm) and untargeted metabolomics prior to IA and at the time of seroconversion to IA in 92 IA cases and 91 controls from the Diabetes Autoimmunity Study in the Young (DAISY). Causal mediation models were used to identify seven DNAm probe-metabolite pairs in which the metabolite measured at IA mediated the protective effect of the DNAm probe measured prior to IA against IA risk. These pairs included five DNAm probes mediated by histidine (a metabolite known to affect T1D risk), one probe (cg01604946) mediated by phostidyl choline p-32:0 or o-32:1, and one probe (cg00390143) mediated by sphingomyelin d34:2. The top 100 DNAm probes were over-represented in six reactome pathways at the FDR &lt, 0.1 level (q = 0.071), including transport of small molecules and inositol phosphate metabolism. While the causal pathways in our mediation models require further investigation to better understand the biological mechanisms, we identified seven methylation sites that may improve our understanding of epigenetic protection against T1D as mediated by the metabolome.

Background: Diastolic dysfunction (DD) is associated with the development of heart failure and contributes to the pathogenesis of other cardiac maladies, including atrial fibrillation. Inhibition of histone deacetylases (HDACs) has been shown to prevent DD by enhancing myofibril relaxation. We addressed the therapeutic potential of HDAC inhibition in a model of established DD with preserved ejection fraction. Methods: Four weeks after uninephrectomy and implantation with deoxycorticosterone acetate pellets, when DD was clearly evident, 1 cohort of mice was administered the clinical-stage HDAC inhibitor ITF2357/Givinostat. Echocardiography, blood pressure measurements, and end point invasive hemodynamic analyses were performed. Myofibril mechanics and intact cardiomyocyte relaxation were assessed ex vivo. Cardiac fibrosis was evaluated by picrosirius red staining and second harmonic generation microscopy of left ventricle (LV) sections, RNA sequencing of LV mRNA, mass spectrometry–based evaluation of decellularized LV biopsies, and atomic force microscopy determination of LV stiffness. Mechanistic studies were performed with primary rat and human cardiac fibroblasts. Results: HDAC inhibition normalized DD without lowering blood pressure in this model of systemic hypertension. In contrast to previous models, myofibril relaxation was unimpaired in uninephrectomy/deoxycorticosterone acetate mice. Furthermore, cardiac fibrosis was not evident in any mouse cohort on the basis of picrosirius red staining or second harmonic generation microscopy. However, mass spectrometry revealed induction in the expression of >100 extracellular matrix proteins in LVs of uninephrectomy/deoxycorticosterone acetate mice, which correlated with profound tissue stiffening based on atomic force microscopy. ITF2357/Givinostat treatment blocked extracellular matrix expansion and LV stiffening. The HDAC inhibitor was subsequently shown to suppress cardiac fibroblast activation, at least in part, by blunting recruitment of the profibrotic chromatin reader protein BRD4 (bromodomain-containing protein 4) to key gene regulatory elements. Conclusions: These findings demonstrate the potential of HDAC inhibition as a therapeutic intervention to reverse existing DD and establish blockade of extracellular matrix remodeling as a second mechanism by which HDAC inhibitors improve ventricular filling. Our data reveal the existence of pathophysiologically relevant covert or hidden cardiac fibrosis that is below the limit of detection of histochemical stains such as picrosirius red, highlighting the need to evaluate fibrosis of the heart using diverse methodologies.

Introduction: Diastolic dysfunction is associated with the development of heart failure with preserved ejection fraction and other cardiac maladies, including atrial fibrillation. Inhibition of histone deacetylase (HDAC) activity has previously been shown to prevent diastolic dysfunction by enhancing myofibril relaxation. Here, we further addressed the therapeutic potential of HDAC inhibition in a model of established diastolic dysfunction with preserved ejection fraction. Methods: Four weeks following uninephrectomy (UNX) and implantation with deoxycorticosterone acetate (DOCA) pellets, when diastolic dysfunction was clearly evident, one cohort of mice was administered the clinical-stage HDAC inhibitor ITF2357/Givinostat. Serial echocardiography, blood pressure measurements, and endpoint invasive hemodynamic analyses were performed, along with the evaluation of pathological cardiac fibrosis and remodeling. Results: HDAC inhibition completely normalized diastolic function without lowering blood pressure in this model of systemic hypertension. Cardiac fibrosis was not evident in any mouse cohorts based on picrosirius red staining or second harmonic generation microscopy. However, mass spectrometry revealed induction in the expression of more than one hundred extracellular matrix (ECM) proteins in LVs of UNX/DOCA mice, which correlated with profound tissue stiffening based on atomic force microscopy. Remarkably, ITF2357/Givinostat treatment entirely blocked ECM expansions and LV stiffening. The HDAC inhibitor was subsequently shown to suppress cardiac fibroblast activation, at least in part, by blunting recruitment of the pro-fibrotic chromatin reader protein, BRD4, to key gene regulatory elements. Conclusions: These findings demonstrate the potential of HDAC inhibition as a therapeutic intervention to reverse existing diastolic dysfunction, and establish blockade of ECM remodeling as a second mechanism by which HDAC inhibitors improve ventricular filling. Additionally, our data reveal the existence of pathophysiologically relevant ‘covert’ cardiac fibrosis that is below the limit of detection of histochemical stains, highlighting the need to evaluate fibrosis of the heart using diverse methodologies.

Objective Illumina BeadChip arrays are commonly used to generate DNA methylation data for large epidemiological studies. Updates in technology over time create challenges for data harmonization within and between studies, many of which obtained data from the older 450K and newer EPIC platforms. The pre-processing pipeline for DNA methylation is not trivial, and influences the downstream analyses. Incorporating different platforms adds a new level of technical variability that has not yet been taken into account by recommended pipelines. Our study evaluated the performance of various tools on different versions of platform data harmonization at each step of pre-processing pipeline, including quality control (QC), normalization, batch effect adjustment, and genomic inflation. We illustrate our novel approach using 450K and EPIC data from the Diabetes Autoimmunity Study in the Young (DAISY) prospective cohort. Results We found normalization and probe filtering had the biggest effect on data harmonization. Employing a meta-analysis was an effective and easily executable method for accounting for platform variability. Correcting for genomic inflation also helped with harmonization. We present guidelines for studies seeking to harmonize data from the 450K and EPIC platforms, which includes the use of technical replicates for evaluating numerous pre-processing steps, and employing a meta-analysis.

Rationale: Allergic diseases are an increasing public health concern, and early life environment is critical to immune development. Maternal diet during pregnancy has been linked to offspring allergy risk. In turn, maternal diet is a potentially modifiable factor, which could be targeted as an allergy prevention strategy. In this systematic review, we focused on non-allergen-specific modifying factors of the maternal diet in pregnancy on allergy outcomes in their offspring. Methods: We undertook a systematic review of studies investigating the association between maternal diet during pregnancy and allergic outcomes (asthma/wheeze, hay fever/allergic rhinitis/seasonal allergies, eczema/atopic dermatitis (AD), food allergies, and allergic sensitization) in offspring. Studies evaluating the effect of food allergen intake were excluded. We searched three bibliographic databases (MEDLINE, EMBASE, and Web of Science) through February 26, 2019. Evidence was critically appraised using modified versions of the Cochrane Collaboration Risk of Bias tool for intervention trials and the National Institute for Clinical Excellence methodological checklist for cohort and case-control studies and meta-analysis performed from RCTs. Results: We identified 95 papers: 17 RCTs and 78 observational (case-control, cross-sectional, and cohort) studies. Observational studies varied in design and dietary intakes and often had contradictory findings. Based on our meta-analysis, RCTs showed that vitamin D supplementation (OR: 0.72; 95% CI: 0.56-0.92) is associated with a reduced risk of wheeze/asthma. A positive trend for omega-3 fatty acids was observed for asthma/wheeze, but this did not reach statistical significance (OR: 0.70; 95% CI: 0.45-1.08). Omega-3 supplementation was also associated with a non-significant decreased risk of allergic rhinitis (OR: 0.76; 95% CI: 0.56-1.04). Neither vitamin D nor omega-3 fatty acids were associated with an altered risk of AD or food allergy. Conclusions: Prenatal supplementation with vitamin D may have beneficial effects for prevention of asthma. Additional nutritional factors seem to be required for modulating the risk of skin and gastrointestinal outcomes. We found no consistent evidence regarding other dietary factors, perhaps due to differences in study design and host features that were not considered. While confirmatory studies are required, there is also a need for performing RCTs beyond single nutrients/foods.

Phthalates are known endocrine disruptors and associated with decreased fecundity, pregnancy loss, and adverse obstetrical outcomes, however the underlying mechanisms remain to be established. Environmental factors can influence gene expression and cell function by modifying epigenetic marks, impacting the developing embryo as well as future generations of offspring. The impact of phthalates on placental gene methylation and expression is largely unknown. We studied the effect of maternal phthalate exposure on the human placental DNA methylome and transcriptome. We determined epigenome-wide DNA methylation marks (Illumina Infinium Human Methylation 850k BeadChip) and gene expression (Agilent whole human genome array) associated with phthalate exposure in first trimester placenta. Integrative genomic analysis of candidate genes was performed to define gene methylation-expression relationships. We identified 39 genes with significantly altered methylation and gene expression in the high phthalate exposure group. Most of these relationships were inversely correlated. This analysis identified epidermal growth factor receptor (EGFR) as a critical candidate gene mediating the effects of phthalates on early placental function. Although additional studies are needed to determine the functional consequences of these changes, our findings are consistent with the model that phthalates impact placental function by modulating the expression of critical placental genes through epigenetic regulation.

The intrauterine period is a critical time wherein developmental exposure can influence risk for chronic disease including childhood obesity. Using umbilical cord-derived mesenchymal stem cells (uMSC) from offspring born to normal-weight and obese mothers, we tested the hypothesis that changes in infant body composition over the first 5 months of life correspond with differences in cellular metabolism and transcriptomic profiles at birth. Higher long-chain acylcarnitine concentrations, lipid transport gene expression, and indicators of oxidative stress in uMSC-adipocytes were related to higher adiposity at 5 months of age. In uMSC-myocytes, lower amino acid concentrations and global differential gene expression for myocyte growth, amino acid biosynthesis, and oxidative stress were related to lower infant percent fat-free mass at 5 months of age, particularly in offspring of obese mothers. This is the first evidence of human infant adipocyte- or myocyte-related alterations in cellular metabolic pathways that correspond with increased adiposity and lower fat-free mass in early infancy. These pathways might reflect the effects of an adverse maternal metabolic environment on the fetal metabolome and genome. Our findings suggest that programmed differences in infant stem cell metabolism correspond with differences in body composition in early life, a known contributor to obesity risk.

Ethyl alcohol is a toxin that, when consumed at high levels, produces organ damage and death. One way to prevent or ameliorate this damage in humans is to reduce the exposure of organs to alcohol by reducing alcohol ingestion. Both the propensity to consume large volumes of alcohol and the susceptibility of human organs to alcohol-induced damage exhibit a strong genetic influence. We have developed an integrative genetic/genomic approach to identify transcriptional networks that predispose complex traits, including propensity for alcohol consumption and propensity for alcohol-induced organ damage. In our approach, the phenotype is assessed in a panel of recombinant inbred (RI) rat strains, and quantitative trait locus (QTL) analysis is performed. Transcriptome data from tissues/organs of naïve RI rat strains are used to identify transcriptional networks using Weighted Gene Coexpression Network Analysis (WGCNA). Correlation of the first principal component of transcriptional coexpression modules with the phenotype across the rat strains, and overlap of QTLs for the phenotype and the QTLs for the coexpression modules (module eigengene QTL) provide the criteria for identification of the functionally related groups of genes that contribute to the phenotype (candidate modules). While we previously identified a brain transcriptional module whose QTL overlapped with a QTL for levels of alcohol consumption in HXB/BXH RI rat strains and 12 selected rat lines, this module did not account for all of the genetic variation in alcohol consumption. Our search for QTL overlap and correlation of coexpression modules with phenotype can, however, be applied to any organ in which the transcriptome has been measured, and this represents a holistic approach in the search for genetic contributors to complex traits. Previous work has implicated liver/brain interactions, particularly involving inflammatory/immune processes, as influencing alcohol consumption levels. We have now analyzed the liver transcriptome of the HXB/BXH RI rat panel in relation to the behavioral trait of alcohol consumption. We used RNA-Seq and microarray data to construct liver transcriptional networks, and identified a liver candidate transcriptional coexpression module that explained 24% of the genetic variance in voluntary alcohol consumption. The transcripts in this module focus attention on liver secretory products that influence inflammatory and immune signaling pathways. We propose that these liver secretory products can interact with brain mechanisms that affect alcohol consumption, and targeting these pathways provides a potential approach to reducing high levels of alcohol intake and also protecting the integrity of the liver and other organs.