Multifaceted immune dysregulation characterizes individuals at-risk for rheumatoid arthritis

Molecular markers of autoimmunity, such as antibodies to citrullinated protein antigens (ACPA), are detectable prior to inflammatory arthritis (IA) in rheumatoid arthritis (RA) and may define a state that is 'at-risk' for future RA. Here we present a cross-sectional comparative analysis among three groups that include ACPA positive individuals without IA (At-Risk), ACPA negative individuals and individuals with early, ACPA positive clinical RA (Early RA). Differential methylation analysis among the groups identifies non-specific dysregulation in peripheral B, memory and naïve T cells in At-Risk participants, with more specific immunological pathway abnormalities in Early RA. Tetramer studies show increased abundance of T cells recognizing citrullinated (cit) epitopes in At-Risk participants, including expansion of T cells reactive to citrullinated cartilage intermediate layer protein I (cit-CILP); these T cells have Th1, Th17, and T stem cell memory-like phenotypes. Antibody-antigen array analyses show that antibodies targeting cit-clusterin, cit-fibrinogen and cit-histone H4 are elevated in At-Risk and Early RA participants, with the highest levels of antibodies detected in those with Early RA. These findings indicate that an ACPA positive at-risk state is associated with multifaceted immune dysregulation that may represent a potential opportunity for targeted intervention.

Hippocampal RNA sequencing in mice selectively bred for high and low activity

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.

Genome-wide analysis of oxylipins and oxylipin profiles in a pediatric population

Background

Oxylipins are inflammatory biomarkers derived from omega-3 and-6 fatty acids implicated in inflammatory diseases but have not been studied in a genome-wide association study (GWAS). The aim of this study was to identify genetic loci associated with oxylipins and oxylipin profiles to identify biologic pathways and therapeutic targets for oxylipins.

Methods

We conducted a GWAS of plasma oxylipins in 316 participants in the Diabetes Autoimmunity Study in the Young (DAISY). DNA samples were genotyped using the TEDDY-T1D Exome array, and additional variants were imputed using the Trans-Omics for Precision Medicine (TOPMed) multi-ancestry reference panel. Principal components analysis of 36 plasma oxylipins was used to capture oxylipin profiles. PC1 represented linoleic acid (LA)- and alpha-linolenic acid (ALA)-related oxylipins, and PC2 represented arachidonic acid (ARA)-related oxylipins. Oxylipin PC1, PC2, and the top five loading oxylipins from each PC were used as outcomes in the GWAS (genome-wide significance: p < 5×10-8).

Results

The SNP rs143070873 was associated with (p < 5×10-8) the LA-related oxylipin 9-HODE, and rs6444933 (downstream of CLDN11) was associated with the LA-related oxylipin 13 S-HODE. A locus between MIR1302-7 and LOC100131146, rs10118380 and an intronic variant in TRPM3 were associated with the ARA-related oxylipin 11-HETE. These loci are involved in inflammatory signaling cascades and interact with PLA2, an initial step to oxylipin biosynthesis.

Conclusion

Genetic loci involved in inflammation and oxylipin metabolism are associated with oxylipin levels.

Propensity scores as a novel method to guide sample allocation and minimize batch effects during the design of high throughput experiments

BACKGROUND

We developed a novel approach to minimize batch effects when assigning samples to batches. Our algorithm selects a batch allocation, among all possible ways of assigning samples to batches, that minimizes differences in average propensity score between batches. This strategy was compared to randomization and stratified randomization in a case-control study (30 per group) with a covariate (case vs control, represented as β1, set to be null) and two biologically relevant confounding variables (age, represented as β2, and hemoglobin A1c (HbA1c), represented as β3). Gene expression values were obtained from a publicly available dataset of expression data obtained from pancreas islet cells. Batch effects were simulated as twice the median biological variation across the gene expression dataset and were added to the publicly available dataset to simulate a batch effect condition. Bias was calculated as the absolute difference between observed betas under the batch allocation strategies and the true beta (no batch effects). Bias was also evaluated after adjustment for batch effects using ComBat as well as a linear regression model. In order to understand performance of our optimal allocation strategy under the alternative hypothesis, we also evaluated bias at a single gene associated with both age and HbA1c levels in the 'true' dataset (CAPN13 gene).

RESULTS

Pre-batch correction, under the null hypothesis (β1), maximum absolute bias and root mean square (RMS) of maximum absolute bias, were minimized using the optimal allocation strategy. Under the alternative hypothesis (β2 and β3 for the CAPN13 gene), maximum absolute bias and RMS of maximum absolute bias were also consistently lower using the optimal allocation strategy. ComBat and the regression batch adjustment methods performed well as the bias estimates moved towards the true values in all conditions under both the null and alternative hypotheses. Although the differences between methods were less pronounced following batch correction, estimates of bias (average and RMS) were consistently lower using the optimal allocation strategy under both the null and alternative hypotheses.

CONCLUSIONS

Our algorithm provides an extremely flexible and effective method for assigning samples to batches by exploiting knowledge of covariates prior to sample allocation.

An Oxylipin-Related Nutrient Pattern and Risk of Type 1 Diabetes in the Diabetes Autoimmunity Study in the Young (DAISY)

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.

COVID-19 Mortality in the Colorado Center for Personalized Medicine Biobank

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.

Evaluation and characterization of expression quantitative trait analysis methods in the Hybrid Rat Diversity Panel

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 four plausible approaches towards modeling quantitative traits in the HRDP and quantify their operating characteristics. In particular, we investigate three approaches based on genome-wide mixed model analysis, and one approach based on ordinary least squares linear regression. 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. In simulation of eQTL data in the HRDP, we find that a mixed model approach that leverages leave-one-chromosome-out kinship estimation attains the highest power while controlling type I error.

Changes in the Coexpression of Innate Immunity Genes During Persistent Islet Autoimmunity Are Associated With Progression of Islet Autoimmunity: Diabetes Autoimmunity Study in the Young (DAISY)

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.

Loss of Frizzled 9 in Lung Cells Alters Epithelial Phenotype and Promotes Premalignant Lesion Development

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.

Relationship Between a Vitamin D Genetic Risk Score and Autoantibodies Among First-Degree Relatives of Probands With Rheumatoid Arthritis and Systemic Lupus Erythematosus

Objective

Higher 25-hydroxyvitamin D (25(OH)D) levels have been associated with reduced risk for autoimmune diseases and are influenced by vitamin D metabolism genes. We estimated genetically-determined vitamin D levels by calculating a genetic risk score (GRS) and investigated whether the vitamin D GRS was associated with the presence of autoantibodies related to rheumatoid arthritis (RA) and systemic lupus erythematosus (SLE) in those at increased risk for developing RA and SLE, respectively.

Methods

In this cross-sectional study, we selected autoantibody positive (aAb+) and autoantibody negative (aAb-) individuals from the Studies of the Etiologies of Rheumatoid Arthritis (SERA), a cohort study of first-degree relatives (FDRs) of individuals with RA (189 RA aAb+, 181 RA aAb-), and the Lupus Family Registry and Repository (LFRR), a cohort study of FDRs of individuals with SLE (157 SLE aAb+, 185 SLE aAb-). Five SNPs known to be associated with serum 25(OH)D levels were analyzed individually as well as in a GRS: rs4588 (GC), rs12785878 (NADSYN1), rs10741657 (CYP2R1), rs6538691 (AMDHD1), and rs8018720 (SEC23A).

Results

Both cohorts had similar demographic characteristics, with significantly older and a higher proportion of males in the aAb+ FDRs. The vitamin D GRS was inversely associated with RA aAb+ (OR = 0.85, 95% CI = 0.74-0.99), suggesting a possible protective factor for RA aAb positivity in FDRs of RA probands. The vitamin D GRS was not associated with SLE aAb+ in the LFRR (OR = 1.09, 95% CI = 0.94-1.27). The SEC23A SNP was associated with RA aAb+ in SERA (OR = 0.65, 95% CI = 0.43-0.99); this SNP was not associated with SLE aAb+ in LFRR (OR = 1.41, 95% CI = 0.90 - 2.19).

Conclusion

Genes associated with vitamin D levels may play a protective role in the development of RA aAbs in FDRs of RA probands, perhaps through affecting lifelong vitamin D status. The GRS and the SEC23A SNP may be of interest for future investigation in pre-clinical RA. In contrast, these results do not support a similar association in SLE FDRs, suggesting other mechanisms involved in the relationship between vitamin D and SLE aAbs not assessed in this study.

An effective processing pipeline for harmonizing DNA methylation data from Illumina's 450K and EPIC platforms for epidemiological studies

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.

A Mediation Approach to Discovering Causal Relationships between the Metabolome and DNA Methylation in Type 1 Diabetes

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 <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.

The oxylipin profile is associated with development of type 1 diabetes: the Diabetes Autoimmunity Study in the Young (DAISY)

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.

Phospholipid Levels at Seroconversion Are Associated With Resolution of Persistent Islet Autoimmunity: The Diabetes Autoimmunity Study in the Young

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.

Inverse probability weighting is an effective method to address selection bias during the analysis of high dimensional data

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.

HDAC Inhibition Reverses Preexisting Diastolic Dysfunction and Blocks Covert Extracellular Matrix Remodeling

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.

Metabolomics-related nutrient patterns at seroconversion and risk of progression to type 1 diabetes

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.

Dietary factors during pregnancy and atopic outcomes in childhood: A systematic review from the European Academy of Allergy and Clinical Immunology

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.

DNA methylation near the INS gene is associated with INS genetic variation (rs689) and type 1 diabetes in the Diabetes Autoimmunity Study in the Young

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 < .05 and validation P < .05). Probes were associated with the GSDMB, C1QTNF6, IL27, and INS genes. The cg03366382 (OR: 1.9, meta-P = .0495), cg21574853 (OR: 2.5, meta-P = .0232), and cg25336198 (odds ratio: 6.6, meta-P = .0081) probes were significantly associated with T1D. The three probes were located upstream from the INS transcription start site.

CONCLUSIONS

We confirmed an association between DNA methylation and rs689 that has been identified in related studies. Measurements in our study preceded the onset of T1D suggesting methylation may have a role in the relationship between INS variation and T1D development.

Longitudinal DNA methylation differences precede type 1 diabetes

DNA methylation may be involved in development of type 1 diabetes (T1D), but previous epigenome-wide association studies were conducted among cases with clinically diagnosed diabetes. Using multiple pre-disease peripheral blood samples on the Illumina 450 K and EPIC platforms, we investigated longitudinal methylation differences between 87 T1D cases and 87 controls from the prospective Diabetes Autoimmunity Study in the Young (DAISY) cohort. Change in methylation with age differed between cases and controls in 10 regions. Average longitudinal methylation differed between cases and controls at two genomic positions and 28 regions. Some methylation differences were detectable and consistent as early as birth, including before and after the onset of preclinical islet autoimmunity. Results map to transcription factors, other protein coding genes, and non-coding regions of the genome with regulatory potential. The identification of methylation differences that predate islet autoimmunity and clinical diagnosis may suggest a role for epigenetics in T1D pathogenesis; however, functional validation is warranted.

Multi-Omic Approaches Identify Metabolic and Autophagy Regulators Important in Ovarian Cancer Dissemination

High-grade serous ovarian cancers (HGSOCs) arise from exfoliation of transformed cells from the fallopian tube, indicating that survival in suspension, and potentially escape from anoikis, is required for dissemination. We report here the results of a multi-omic study to identify drivers of anoikis escape, including transcriptomic analysis, global non-targeted metabolomics, and a genome-wide CRISPR/Cas9 knockout (GeCKO) screen of HGSOC cells cultured in adherent and suspension settings. Our combined approach identified known pathways, including NOTCH signaling, as well as novel regulators of anoikis escape. Newly identified genes include effectors of fatty acid metabolism, ACADVL and ECHDC2, and an autophagy regulator, ULK1. Knockdown of these genes significantly inhibited suspension growth of HGSOC cells, and the metabolic profile confirmed the role of fatty acid metabolism in survival in suspension. Integration of our datasets identified an anoikis-escape gene signature that predicts overall survival in many carcinomas.

Unsupervised, Statistically Based Systems Biology Approach for Unraveling the Genetics of Complex Traits: A Demonstration with Ethanol Metabolism

BACKGROUND

A statistical pipeline was developed and used for determining candidate genes and candidate gene coexpression networks involved in 2 alcohol (i.e., ethanol [EtOH]) metabolism phenotypes, namely alcohol clearance and acetate area under the curve in a recombinant inbred (RI) (HXB/BXH) rat panel. The approach was also used to provide an indication of how EtOH metabolism can impact the normal function of the identified networks.

METHODS

RNA was extracted from alcohol-naïve liver tissue of 30 strains of HXB/BXH RI rats. The reconstructed transcripts were quantitated, and data were used to construct gene coexpression modules and networks. A separate group of rats, comprising the same 30 strains, were injected with EtOH (2 g/kg) for measurement of blood EtOH and acetate levels. These data were used for quantitative trait loci (QTL) analysis of the rate of EtOH disappearance and circulating acetate levels. The analysis pipeline required calculation of the module eigengene values, the correction of these values with EtOH metabolism rates and acetate levels across the rat strains, and the determination of the eigengene QTLs. For a module to be considered a candidate for determining phenotype, the module eigengene values had to have significant correlation with the strain phenotypic values and the module eigengene QTLs had to overlap the phenotypic QTLs.

RESULTS

Of the 658 transcript coexpression modules generated from liver RNA sequencing data, a single module satisfied all criteria for being a candidate for determining the alcohol clearance trait. This module contained 2 alcohol dehydrogenase genes, including the gene whose product was previously shown to be responsible for the majority of alcohol elimination in the rat. This module was also the only module identified as a candidate for influencing circulating acetate levels. This module was also linked to the process of generation and utilization of retinoic acid as related to the autonomous immune response.

CONCLUSIONS

We propose that our analytical pipeline can successfully identify genetic regions and transcripts which predispose a particular phenotype and our analysis provides functional context for coexpression module components.

Voluntary exposure to a toxin: the genetic influence on ethanol consumption

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.

Nuclear CD24 Drives Tumor Growth and Is Predictive of Poor Patient Prognosis

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.

The sequenced rat brain transcriptome--its use in identifying networks predisposing alcohol consumption

A quantitative genetic approach, which involves correlation of transcriptional networks with the phenotype in a recombinant inbred (RI) population and in selectively bred lines of rats, and determination of coinciding quantitative trait loci for gene expression and the trait of interest, has been applied in the present study. In this analysis, a novel approach was used that combined DNA-Seq data, data from brain exon array analysis of HXB/BXH RI rat strains and six pairs of rat lines selectively bred for high and low alcohol preference, and RNA-Seq data (including rat brain transcriptome reconstruction) to quantify transcript expression levels, generate co-expression modules and identify biological functions that contribute to the predisposition of consuming varying amounts of alcohol. A gene co-expression module was identified in the RI rat strains that contained both annotated and unannotated transcripts expressed in the brain, and was associated with alcohol consumption in the RI panel. This module was found to be enriched with differentially expressed genes from the selected lines of rats. The candidate genes within the module and differentially expressed genes between high and low drinking selected lines were associated with glia (microglia and astrocytes) and could be categorized as being related to immune function, energy metabolism and calcium homeostasis, as well as glial-neuronal communication. The results of the present study show that there are multiple combinations of genetic factors that can produce the same phenotypic outcome. Although no single gene accounts for predisposition to a particular level of alcohol consumption in every animal model, coordinated differential expression of subsets of genes in the identified pathways produce similar phenotypic outcomes.

DATABASE

The datasets supporting the results of the present study are available at http://phenogen.ucdenver.edu.

Genome-wide association discoveries of alcohol dependence

OBJECTIVE

To report the genome-wide significant and/or replicable risk variants for alcohol dependence and explore their potential biological functions.

METHODS

We searched in PubMed for all genome-wide association studies (GWASs) of alcohol dependence. The following three types of the results were extracted: genome-wide significant associations in an individual sample, the combined samples, or the meta-analysis (p < 5 × 10(-8) ); top-ranked associations in an individual sample (p < 10(-5) ) that were nominally replicated in other samples (p < .05); and nominally replicable associations across at least three independent GWAS samples (p < .05). These results were meta-analyzed. cis-eQTLs in human, RNA expression in rat and mouse brains and bioinformatics properties of all of these risk variants were analyzed.

RESULTS

The variants located within the alcohol dehydrogenase (ADH) cluster were significantly associated with alcohol dependence at the genome-wide level (p < 5 × 10(-8) ) in at least one sample. Some associations with the ADH cluster were replicable across six independent GWAS samples. The variants located within or near SERINC2, KIAA0040, MREG-PECR or PKNOX2 were significantly associated with alcohol dependence at the genome-wide level (p < 5 × 10(-8) ) in meta-analysis or combined samples, and these associations were replicable across at least one sample. The associations with the variants within NRD1, GPD1L-CMTM8 or MAP3K9-PCNX were suggestive (5 × 10(-8)  < p < 10(-5) ) in some samples, and nominally replicable in other samples. The associations with the variants at HTR7 and OPA3 were nominally replicable across at least three independent GWAS samples (10(-5)  < p < .05). Some risk variants at the ADH cluster, SERINC2, KIAA0040, NRD1, and HTR7 had potential biological functions.

CONCLUSION

The most robust risk locus was the ADH cluster. SERINC2, KIAA0040, NRD1, and HTR7 were also likely to play important roles in alcohol dependence. PKNOX2, MREG, PECR, GPD1L, CMTM8, MAP3K9, PCNX, and OPA3 might play less important roles in risk for alcohol dependence based on the function analysis. This conclusion will significantly contribute to the post-GWAS follow-up studies on alcohol dependence.

Influence of sex on genetic regulation of "drinking in the dark" alcohol consumption

The ILSXISS (LXS) recombinant inbred (RI) panel of mice is a valuable resource for genetic mapping studies of complex traits, due to its genetic diversity and large number of strains. Male and female mice from this panel were used to investigate genetic influences on alcohol consumption in the "drinking in the dark" (DID) model. Male mice (38 strains) and female mice (36 strains) were given access to 20% ethanol during the early phase of their circadian dark cycle for four consecutive days. The first principal component of alcohol consumption measures on days 2, 3, and 4 was used as a phenotype (DID phenotype) to calculate QTLs, using a SNP marker set for the LXS RI panel. Five QTLs were identified, three of which included a significant genotype by sex interaction, i.e., a significant genotype effect in males and not females. To investigate candidate genes associated with the DID phenotype, data from brain microarray analysis (Affymetrix Mouse Exon 1.0 ST Arrays) of male LXS RI strains were combined with RNA-Seq data (mouse brain transcriptome reconstruction) from the parental ILS and ISS strains in order to identify expressed mouse brain transcripts. Candidate genes were determined based on common eQTL and DID phenotype QTL regions and correlation of transcript expression levels with the DID phenotype. The resulting candidate genes (in particular, Arntl/Bmal1) focused attention on the influence of circadian regulation on the variation in the DID phenotype in this population of mice.

Is the alcohol deprivation effect genetically mediated? Studies with HXB/BXH recombinant inbred rat strains

BACKGROUND

Two features of alcohol addiction that have been widely studied in animal models are relapse drinking following periods of alcohol abstinence and the escalation of alcohol consumption after chronic continuous or intermittent alcohol exposure. The genetic contribution to these phenotypes has not been systematically investigated.

METHODS

HXB/BXH recombinant inbred (RI) rat strains were given access to alcohol sequentially as follows: alcohol (10%) as the only fluid for 1 week; alcohol (10%) and water in a 2-bottle choice paradigm for 7 weeks ("pre-alcohol deprivation effect [ADE] alcohol consumption"); 2 weeks of access to water only (alcohol deprivation); and 2 weeks of reaccess to 10% alcohol and water ("post-ADE alcohol consumption"). The periods of deprivation and reaccess to alcohol were repeated 3 times. The ADE was defined as the amount of alcohol consumed in the first 24 hours after deprivation minus the average daily amount of alcohol consumed in the week prior to deprivation. Heritability of the phenotypes was determined by analysis of variance, and quantitative trait loci (QTLs) were identified.

RESULTS

All strains showed increased alcohol consumption, compared to the predeprivation period, in the first 24 hours after each deprivation (ADE). Broad-sense heritability of the ADEs was low (ADE1, 9.1%; ADE2, 26.2%; ADE3, 16.3%). Alcohol consumption levels were relatively stable over weeks 2 to 7. Post-ADE alcohol consumption levels consistently increased in some strains and were decreased or unchanged in others. Heritability of pre- and post-ADE alcohol consumption was high and increased over time (week 2, 38.5%; week 7, 51.1%; week 11, 56.8%; week 15, 63.3%). QTLs for pre- and post-ADE alcohol consumption were similar, but the strength of the QTL association with the phenotype decreased over time.

CONCLUSIONS

In the HXB/BXH RI rat strains, genotypic variance does not account for a large proportion of phenotypic variance in the ADE phenotype (low heritability), suggesting a role of environmental factors. In contrast, a large proportion of the variance across the RI strains in pre- and post-ADE alcohol consumption is due to genetically determined variance (high heritability).

Rat mammary extracellular matrix composition and response to ibuprofen treatment during postpartum involution by differential GeLC-MS/MS analysis

Breast cancer patients diagnosed within five years following pregnancy have increased metastasis and decreased survival. A hallmark of postpartum biology that may contribute to this poor prognosis is mammary gland involution, involving massive epithelial cell death and dramatic stromal remodeling. Previous studies show pro-tumorigenic properties of extracellular matrix (ECM) isolated from rodent mammary glands undergoing postpartum involution. More recent work demonstrates systemic ibuprofen treatment during involution decreases its tumor-promotional nature. Utilizing a proteomics approach, we identified relative differences in the composition of mammary ECM isolated from nulliparous rats and those undergoing postpartum involution, with and without ibuprofen treatment. GeLC-MS/MS experiments resulted in 20327 peptide identifications that mapped to 884 proteins with a <0.02% false discovery rate. Label-free quantification yielded several ECM differences between nulliparous and involuting glands related to collagen-fiber organization, cell motility and attachment, and cytokine regulation. Increases in known pro-tumorigenic ECM proteins osteopontin, tenascin-C, and laminin-α1 and pro-inflammatory proteins STAT3 and CD68 further identify candidate mediators of breast cancer progression specific to the involution window. With postpartum ibuprofen treatment, decreases in tenascin-C and three laminin chains were revealed. Our data suggest novel ECM mediators of breast cancer progression and demonstrate a protective influence of ibuprofen on mammary ECM composition.