Sex effects on DNA methylation affect discovery in epigenome-wide association study of schizophrenia.

Sex differences in the epidemiology and clinical characteristics of schizophrenia are well-known; however, the molecular mechanisms underlying these differences remain unclear. Further, the potential advantages of sex-stratified meta-analyses of epigenome-wide association studies (EWAS) of schizophrenia have not been investigated. Here, we performed sex-stratified EWAS meta-analyses to investigate whether sex stratification improves discovery, and to identify differentially methylated regions (DMRs) in schizophrenia. Peripheral blood-derived DNA methylation data from 1519 cases of schizophrenia (male n = 989, female n = 530) and 1723 controls (male n = 997, female n = 726) from three publicly available datasets, and the TOP cohort were meta-analyzed to compare sex-specific, sex-stratified, and sex-adjusted EWAS. The predictive power of each model was assessed by polymethylation score (PMS). The number of schizophrenia-associated differentially methylated positions identified was higher for the sex-stratified model than for the sex-adjusted one. We identified 20 schizophrenia-associated DMRs in the sex-stratified analysis. PMS from sex-stratified analysis outperformed that from sex-adjusted analysis in predicting schizophrenia. Notably, PMSs from the sex-stratified and female-only analyses, but not those from sex-adjusted or the male-only analyses, significantly predicted schizophrenia in males. The findings suggest that sex-stratified EWAS meta-analyses improve the identification of schizophrenia-associated epigenetic changes and highlight an interaction between sex and schizophrenia status on DNA methylation. Sex-specific DNA methylation may have potential implications for precision psychiatry and the development of stratified treatments for schizophrenia.

Neurobiological mechanisms of ECT and TMS treatment in depression: study protocol of a multimodal magnetic resonance investigation.

BACKGROUND: Noninvasive neurostimulation treatments are increasingly being used to treat major depression, which is a common cause of disability worldwide. While electroconvulsive therapy (ECT) and transcranial magnetic stimulation (TMS) are both effective in treating depressive episodes, their mechanisms of action are, however, not completely understood. ECT is given under general anesthesia, where an electrical pulse is administered through electrodes placed on the patient's head to trigger a seizure. ECT is used for the most severe cases of depression and is usually not prescribed before other options have failed. With TMS, brain stimulation is achieved through rapidly changing magnetic fields that induce electric currents underneath a ferromagnetic coil. Its efficacy in depressive episodes has been well documented. This project aims to identify the neurobiological underpinnings of both the effects and side effects of the neurostimulation techniques ECT and TMS. METHODS: The study will utilize a pre-post case control longitudinal design. The sample will consist of 150 subjects: 100 patients (bipolar and major depressive disorder) who are treated with either ECT (N = 50) or TMS (N = 50) and matched healthy controls (N = 50) not receiving any treatment. All participants will undergo multimodal magnetic resonance imaging (MRI) as well as neuropsychological and clinical assessments at multiple time points before, during and after treatment. Arterial spin labeling MRI at baseline will be used to test whether brain perfusion can predict outcomes. Signs of brain disruption, potentiation and rewiring will be explored with resting-state functional MRI, magnetic resonance spectroscopy and multishell diffusion weighted imaging (DWI). Clinical outcome will be measured by clinician assessed and patient reported outcome measures. Memory-related side effects will be investigated, and specific tests of spatial navigation to test hippocampal function will be administered both before and after treatment. Blood samples will be stored in a biobank for future analyses. The observation time is 6 months. Data will be explored in light of the recently proposed disrupt, potentiate and rewire (DPR) hypothesis. DISCUSSION: The study will contribute data and novel analyses important for our understanding of neurostimulation as well as for the development of enhanced and more personalized treatment. TRIAL REGISTRATION: ClinicalTrials.gov Identifier: NCT05135897.

Mental health symptoms during the first months of the COVID-19 outbreak in Norway: A cross-sectional survey study.

BACKGROUND: The COVID-19 pandemic has led to major social and economic changes that could impact public mental health. The main aim of the current study was to investigate mental health in Norway during the COVID-19 outbreak (since the first confirmed case on 26 February 2020). METHODS: The results are from the first wave of the data collection (1 April-2 June 2020), which took place during the outbreak along with its initial restrictions. A total of 19,372 (11,883 students) people participated in a cross-sectional web-based survey. RESULTS: A total of 21.8% scored above the cut-off for depression and 23.7% for anxiety. Severity of symptoms was associated with the accumulation of risk factors, such as possible/confirmed infection for oneself or one's family, female/other sex, students, having mental health problems, increased use of tobacco, increased use of alcohol, less exercise, losing one's job, suffering economic impact and lower education. CONCLUSIONS: COVID-19 could have a negative association with public mental health, especially for certain risk groups. Future data-collection waves will provide further insight into the development of symptoms following the pandemic.

SBMLmod: a Python-based web application and web service for efficient data integration and model simulation.

BACKGROUND: Systems Biology Markup Language (SBML) is the standard model representation and description language in systems biology. Enriching and analysing systems biology models by integrating the multitude of available data, increases the predictive power of these models. This may be a daunting task, which commonly requires bioinformatic competence and scripting. RESULTS: We present SBMLmod, a Python-based web application and service, that automates integration of high throughput data into SBML models. Subsequent steady state analysis is readily accessible via the web service COPASIWS. We illustrate the utility of SBMLmod by integrating gene expression data from different healthy tissues as well as from a cancer dataset into a previously published model of mammalian tryptophan metabolism. CONCLUSION: SBMLmod is a user-friendly platform for model modification and simulation. The web application is available at http://sbmlmod.uit.no , whereas the WSDL definition file for the web service is accessible via http://sbmlmod.uit.no/SBMLmod.wsdl . Furthermore, the entire package can be downloaded from https://github.com/MolecularBioinformatics/sbml-mod-ws . We envision that SBMLmod will make automated model modification and simulation available to a broader research community.

Interactively illustrating polymerization using three-level model fusion.

BACKGROUND: Research in cell biology is steadily contributing new knowledge about many aspects of physiological processes, both with respect to the involved molecular structures as well as their related function. Illustrations of the spatio-temporal development of such processes are not only used in biomedical education, but also can serve scientists as an additional platform for in-silico experiments. RESULTS: In this paper, we contribute a new, three-level modeling approach to illustrate physiological processes from the class of polymerization at different time scales. We integrate physical and empirical modeling, according to which approach best suits the different involved levels of detail, and we additionally enable a form of interactive steering, while the process is illustrated. We demonstrate the suitability of our approach in the context of several polymerization processes and report from a first evaluation with domain experts. CONCLUSION: We conclude that our approach provides a new, hybrid modeling approach for illustrating the process of emergence in physiology, embedded in a densely filled environment. Our approach of a complementary fusion of three systems combines the strong points from the different modeling approaches and is capable to bridge different spatial and temporal scales.

Model of tryptophan metabolism, readily scalable using tissue-specific gene expression data.

Tryptophan is utilized in various metabolic routes including protein synthesis, serotonin, and melatonin synthesis and the kynurenine pathway. Perturbations in these pathways have been associated with neurodegenerative diseases and cancer. Here we present a comprehensive kinetic model of the complex network of human tryptophan metabolism based upon existing kinetic data for all enzymatic conversions and transporters. By integrating tissue-specific expression data, modeling tryptophan metabolism in liver and brain returned intermediate metabolite concentrations in the physiological range. Sensitivity and metabolic control analyses identified expected key enzymes to govern fluxes in the branches of the network. Combining tissue-specific models revealed a considerable impact of the kynurenine pathway in liver on the concentrations of neuroactive derivatives in the brain. Moreover, using expression data from a cancer study predicted metabolite changes that resembled the experimental observations. We conclude that the combination of the kinetic model with expression data represents a powerful diagnostic tool to predict alterations in tryptophan metabolism. The model is readily scalable to include more tissues, thereby enabling assessment of organismal tryptophan metabolism in health and disease.

Divergent epigenetic responses to perinatal asphyxia in severe mental disorders.

Epigenetic modifications influenced by environmental exposures are molecular sources of phenotypic heterogeneity found in schizophrenia and bipolar disorder and may contribute to shared etiopathogenetic mechanisms of these two disorders. Newborns who experienced perinatal asphyxia have suffered reduced oxygen delivery to the brain around the time of birth, which increases the risk of later psychiatric diagnosis. This study aimed to investigate DNA methylation in blood cells for associations with a history of perinatal asphyxia, a neurologically harmful condition occurring within the biological environment of birth. We utilized prospective data from the Medical Birth Registry of Norway to identify incidents of perinatal asphyxia in 643 individuals with schizophrenia or bipolar disorder and 676 healthy controls. We performed an epigenome wide association study to distinguish differentially methylated positions associated with perinatal asphyxia. We found an interaction between methylation and exposure to perinatal asphyxia on case-control status, wherein having a history of perinatal asphyxia was associated with an increase of methylation in healthy controls and a decrease of methylation in patients on 4 regions of DNA important for brain development and function. The differentially methylated regions were observed in genes involved in oligodendrocyte survival and axonal myelination and functional recovery (LINGO3); assembly, maturation and maintenance of the brain (BLCAP;NNAT and NANOS2) and axonal transport processes and neural plasticity (SLC2A14). These findings are consistent with the notion that an opposite epigenetic response to perinatal asphyxia, in patients compared with controls, may contribute to molecular mechanisms of risk for schizophrenia and bipolar disorder.

Differences in white blood cell proportions between schizophrenia cases and controls are influenced by medication and variations in time of day.

Cases with schizophrenia (SCZ) and healthy controls show differences in white blood cell (WBC) counts and blood inflammation markers. Here, we investigate whether time of blood draw and treatment with psychiatric medications are related to differences in estimated WBC proportions between SCZ cases and controls. DNA methylation data from whole blood was used to estimate proportions of six subtypes of WBCs in SCZ patients (n = 333) and healthy controls (n = 396). We tested the association of case-control status with estimated cell-type proportions and the neutrophil-to-lymphocyte ratio (NLR) in 4 models: with/without adjusting for time of blood draw, and then compared results from blood samples drawn during a 12-h (07:00-19:00) or 7-h (07:00-14:00) period. We also investigated WBC proportions in a subgroup of medication-free patients (n = 51). Neutrophil proportions were significantly higher in SCZ cases (mean=54.1%) vs. controls (mean=51.1%; p = <0.001), and CD8+T lymphocyte proportions were lower in SCZ cases (mean=12.1%) vs. controls (mean=13.2%; p = 0.001). The effect sizes in the 12-h sample (07:00-19:00) showed a significant difference between SCZ vs. controls for neutrophils, CD4+T, CD8+T, and B-cells, which remained significant after adjusting for time of blood draw. In the samples matched for time of blood draw during 07.00-14.00, we also observed an association with neutrophils, CD4+T, CD8+T, and B-cells that was unaffected by further adjustment for time of blood draw. In the medication-free patients, we observed differences that remained significant in neutrophils (p = 0.01) and CD4+T (p = 0.01) after adjusting for time of day. The association of SCZ with NLR was significant in all models (range: p < 0.001 to p = 0.03) in both medicated and unmedicated patients. In conclusion, controlling for pharmacological treatment and circadian cycling of WBC is necessary for unbiased estimates in case-control studies. Nevertheless, the association of WBC with SCZ remains, even after adjusting for the time of day.

An epigenetic association analysis of childhood trauma in psychosis reveals possible overlap with methylation changes associated with PTSD.

Patients with a severe mental disorder report significantly higher levels of childhood trauma (CT) than healthy individuals. Studies have suggested that CT may affect brain plasticity through epigenetic mechanisms and contribute to developing various psychiatric disorders. We performed a blood-based epigenome-wide association study using the Childhood Trauma Questionnaire-short form in 602 patients with a current severe mental illness, investigating DNA methylation association separately for five trauma subtypes and the total trauma score. The median trauma score was set as the predefined cutoff for determining whether the trauma was present or not. Additionally, we compared our genome-wide results with methylation probes annotated to candidate genes previously associated with CT. Of the patients, 83.2% reported CT above the cutoff in one or more trauma subtypes, and emotional neglect was the trauma subtype most frequently reported. We identified one significant differently methylated position associated with the gene TANGO6 for physical neglect. Seventeen differentially methylated regions (DMRs) were associated with different trauma categories. Several of these DMRs were annotated to genes previously associated with neuropsychiatric disorders such as post-traumatic stress disorder and cognitive impairments. Our results support a biomolecular association between CT and severe mental disorders. Genes that were previously identified as differentially methylated in CT-exposed subjects with and without psychosis did not show methylation differences in our analysis. We discuss this inconsistency, the relevance of our findings, and the limitations of our study.

Cannabis Use Is Associated With Increased Levels of Soluble gp130 in Schizophrenia but Not in Bipolar Disorder.

The complex effects of plant cannabinoids on human physiology is not yet fully understood, but include a wide spectrum of effects on immune modulation. The immune system and its inflammatory effector pathways are recently emerging as possible causative factors in psychotic disorders. The present study aimed to investigate whether self-administered Cannabis use was associated with changes in circulating immune and neuroendocrine markers in schizophrenia (SCZ) and bipolar disorder (BD) patients. A screening of 13 plasma markers reflecting different inflammatory pathways was performed in SCZ (n = 401) and BD patients (n = 242) after subdividing each group into Cannabis user and non-user subgroups. We found that i) soluble gp130 (sgp130) concentrations were significantly elevated among Cannabis users in the SCZ group (p = 0.002) after multiple testing correction, but not in BD. ii) Nominally significant differences were observed in the levels of IL-1RA (p = 0.0059), YKL40 (p = 0.0069), CatS (p = 0.013), sTNFR1 (p = 0.031), and BDNF (p = 0.020), where these factors exhibited higher plasma levels in Cannabis user SCZ patients than in non-users. iii) These differences in systemic levels were not reflected by altered mRNA expression of genes encoding sgp130, IL-1RA, YKL40, CatS, sTNFR1, and BDNF in whole blood. Our results show that Cannabis self-administration is associated with markedly higher sgp130 levels in SCZ, but not in BD, and that this phenomenon is independent of the modulation of peripheral immune cells. These findings warrant further investigation into the potential IL-6 trans-signaling modulatory, anti-inflammatory, neuroimmune, and biobehavioral-cognitive effects of Cannabis use in SCZ.

Genome-wide transcription analysis of histidine-related cataract in Atlantic salmon (Salmo salar L).

PURPOSE: Elevated levels of dietary histidine have previously been shown to prevent or mitigate cataract formation in farmed Atlantic salmon (Salmo salar L). The aim of this study was to shed light on the mechanisms by which histidine acts. Applying microarray analysis to the lens transcriptome, we screened for differentially expressed genes in search for a model explaining cataract development in Atlantic salmon and possible markers for early cataract diagnosis. METHODS: Adult Atlantic salmon (1.7 kg) were fed three standard commercial salmon diets only differing in the histidine content (9, 13, and 17 g histidine/kg diet) for four months. Individual cataract scores for both eyes were assessed by slit-lamp biomicroscopy. Lens N-acetyl histidine contents were measured by high performance liquid chromatography (HPLC). Total RNA extracted from whole lenses was analyzed using the GRASP 16K salmonid microarray. The microarray data were analyzed using J-Express Pro 2.7 and validated by quantitative real-time polymerase chain reaction (qRT-PCR). RESULTS: Fish developed cataracts with different severity in response to dietary histidine levels. Lens N-acetyl histidine contents reflected the dietary histidine levels and were negatively correlated to cataract scores. Significance analysis of microarrays (SAM) revealed 248 significantly up-regulated transcripts and 266 significantly down-regulated transcripts in fish that were fed a low level of histidine compared to fish fed a higher histidine level. Among the differentially expressed transcripts were metallothionein A and B as well as transcripts involved in lipid metabolism, carbohydrate metabolism, regulation of ion homeostasis, and protein degradation. Hierarchical clustering and correspondence analysis plot confirmed differences in gene expression between the feeding groups. The differentially expressed genes could be categorized as "early" and "late" responsive according to their expression pattern relative to progression in cataract formation. CONCLUSIONS: Dietary histidine regimes affected cataract formation and lens gene expression in adult Atlantic salmon. Regulated transcripts selected from the results of this genome-wide transcription analysis might be used as possible biological markers for cataract development in Atlantic salmon.

Pharmacokinetics and transcriptional effects of the anti-salmon lice drug emamectin benzoate in Atlantic salmon (Salmo salar L.).

BACKGROUND: Emamectin benzoate (EB) is a dominating pharmaceutical drug used for the treatment and control of infections by sea lice (Lepeophtheirus salmonis) on Atlantic salmon (Salmo salar L). Fish with an initial mean weight of 132 g were experimentally medicated by a standard seven-day EB treatment, and the concentrations of drug in liver, muscle and skin were examined. To investigate how EB affects Atlantic salmon transcription in liver, tissues were assessed by microarray and qPCR at 7, 14 and 35 days after the initiation of medication. RESULTS: The pharmacokinetic examination revealed highest EB concentrations in all three tissues at day 14, seven days after the end of the medication period. Only modest effects were seen on the transcriptional levels in liver, with small fold-change alterations in transcription throughout the experimental period. Gene set enrichment analysis (GSEA) indicated that EB treatment induced oxidative stress at day 7 and inflammation at day 14. The qPCR examinations showed that medication by EB significantly increased the transcription of both HSP70 and glutathione-S-transferase (GST) in liver during a period of 35 days, compared to un-treated fish, possibly via activation of enzymes involved in phase II conjugation of metabolism in the liver. CONCLUSION: This study has shown that a standard seven-day EB treatment has only a modest effect on the transcription of genes in liver of Atlantic salmon. Based on GSEA, the medication seems to have produced a temporary oxidative stress response that might have affected protein stability and folding, followed by a secondary inflammatory response.

Global Gene Expression Response in Peripheral Blood Cells of Petroleum Workers Exposed to Sub-Ppm Benzene Levels.

Altered gene expression in pathways relevant to leukaemogenesis, as well as reduced levels of circulating lymphocytes, have been reported in workers that were exposed to benzene concentrations below 1 ppm. In this study, we analysed whole blood global gene expression patterns in a worker cohort with altered levels of T cells and immunoglobulins IgM and IgA at three time points; pre-shift, post-shift (after three days), and post-recovery (12 hours later). Eight benzene exposed tank workers performing maintenance work in crude oil cargo tanks with a mean benzene exposure of 0.3 ppm (range 0.1⁻0.5 ppm) and five referents considered to be unexposed were examined by gene expression arrays. By using our data as independent validation, we reanalysed selected genes that were reported to be altered from previous studies of workers being exposed to sub-ppm benzene levels Four out of six genes previously proposed as marker genes in chronically exposed workers separated benzene exposed workers from unexposed referents (CLEC5, ACSL1, PRG2, IFNB1). Even better separation of benzene exposed workers and referents was observed for short-term exposure for genes in the Jak-STAT pathway, particularly elevated expression of IL6 and reduced expression of IL19.

Effect of substrate competition in kinetic models of metabolic networks.

Substrate competition can be found in many types of biological processes, ranging from gene expression to signal transduction and metabolic pathways. Although several experimental and in silico studies have shown the impact of substrate competition on these processes, it is still often neglected, especially in modelling approaches. Using toy models that exemplify different metabolic pathway scenarios, we show that substrate competition can influence the dynamics and the steady state concentrations of a metabolic pathway. We have additionally derived rate laws for substrate competition in reversible reactions and summarise existing rate laws for substrate competition in irreversible reactions.

Hepatic in vitro toxicity assessment of PBDE congeners BDE47, BDE153 and BDE154 in Atlantic salmon (Salmo salar L.).

The brominated flame retardant congeners BDE47, BDE153 and BDE154 are among the congeners accumulating to the highest degree in fish. In order to gain knowledge about the toxicological effects of PBDEs in fish, microarray-based transcriptomic and 2D-DIGE/MALDI-TOF/TOF proteomic approaches were used to screen for effects in primary Atlantic salmon hepatocytes exposed to these congeners alone or in combination (PBDE-MIX). A small set of stress related transcripts and proteins were differentially expressed in the PBDE exposed hepatocytes. The PBDE-MIX, and BDE153 to a lesser degree, seems to have induced metabolic disturbances by affecting several pathways related to glucose homeostasis. Further, effects on cell cycle control and proliferation signal pathways in PBDE-MIX-exposed hepatocytes clearly suggest that the PBDE exposure affected cell proliferation processes. CYP1A was 7.41- and 7.37-fold up-regulated in hepatocytes exposed to BDE47 and PBDE-MIX, respectively, and was the only biotransformation pathway affected by the PBDE exposure. The factorial design and PLS regression analyses of the effect of the PBDE-MIX indicated that BDE47 contributed the most to the observed CYP1A response, suggesting that this congener should be incorporated in the toxic equivalent (TEQ) concept in future risk assessment of dioxin-like chemicals. Additionally, a significant up-regulation of the ER-responsive genes VTG and ZP3 was observed in cells exposed to BDE47 and PBDE-MIX. Further analyses suggested that BDE47 and BDE154 have an estrogenic effect in male fish. The data also suggested an antagonistic interaction between BDE153 and BDE154. In conclusion, this study shows that PBDEs can affect several biological systems in Atlantic salmon cells, and demonstrates the need for more studies on the simultaneous exposure to chemical mixtures to identify combined effects of chemicals.

Modulation of transcriptional and inflammatory responses in murine macrophages by the Mycobacterium tuberculosis mammalian cell entry (Mce) 1 complex.

The outcome of many infections depends on the initial interactions between agent and host. Aiming at elucidating the effect of the M. tuberculosis Mce1 protein complex on host transcriptional and immunological responses to infection with M. tuberculosis, RNA from murine macrophages at 15, 30, 60 min, 4 and 10 hrs post-infection with M. tuberculosis H37Rv or Δ-mce1 H37Rv was analyzed by whole-genome microarrays and RT-QPCR. Immunological responses were measured using a 23-plex cytokine assay. Compared to uninfected controls, 524 versus 64 genes were up-regulated by 15 min post H37Rv- and Δ-mce1 H37Rv-infection, respectively. By 15 min post-H37Rv infection, a decline of 17 cytokines combined with up-regulation of Ccl24 (26.5-fold), Clec4a2 (23.2-fold) and Pparγ (10.5-fold) indicated an anti-inflammatory response initiated by IL-13. Down-regulation of Il13ra1 combined with up-regulation of Il12b (30.2-fold), suggested switch to a pro-inflammatory response by 4 hrs post H37Rv-infection. Whereas no significant change in cytokine concentration or transcription was observed during the first hour post Δ-mce1 H37Rv-infection, a significant decline of IL-1b, IL-9, IL-13, Eotaxin and GM-CSF combined with increased transcription of Il12b (25.1-fold) and Inb1 (17.9-fold) by 4 hrs, indicated a pro-inflammatory response. The balance between pro-and anti-inflammatory responses during the early stages of infection may have significant bearing on outcome.

Switch from stress response to homeobox transcription factors in adipose tissue after profound fat loss.

BACKGROUND: In obesity, impaired adipose tissue function may promote secondary disease through ectopic lipid accumulation and excess release of adipokines, resulting in systemic low-grade inflammation, insulin resistance and organ dysfunction. However, several of the genes regulating adipose tissue function in obesity are yet to be identified. METHODOLOGY/PRINCIPAL FINDINGS: In order to identify novel candidate genes that may regulate adipose tissue function, we analyzed global gene expression in abdominal subcutaneous adipose tissue before and one year after bariatric surgery (biliopancreatic diversion with duodenal switch, BPD/DS) (n = 16). Adipose tissue from lean healthy individuals was also analyzed (n = 13). Two different microarray platforms (AB 1700 and Illumina) were used to measure the differential gene expression, and the results were further validated by qPCR. Surgery reduced BMI from 53.3 to 33.1 kg/m(2). The majority of differentially expressed genes were down-regulated after profound fat loss, including transcription factors involved in stress response, inflammation, and immune cell function (e.g., FOS, JUN, ETS, C/EBPB, C/EBPD). Interestingly, a distinct set of genes was up-regulated after fat loss, including homeobox transcription factors (IRX3, IRX5, HOXA5, HOXA9, HOXB5, HOXC6, EMX2, PRRX1) and extracellular matrix structural proteins (COL1A1, COL1A2, COL3A1, COL5A1, COL6A3). CONCLUSIONS/SIGNIFICANCE: The data demonstrate a marked switch of transcription factors in adipose tissue after profound fat loss, providing new molecular insight into a dichotomy between stress response and metabolically favorable tissue development. Our findings implicate homeobox transcription factors as important regulators of adipose tissue function.

Analysis of gene-expression data using J-Express.

The J-Express package has been designed to facilitate the analysis of microarray data with an emphasis on efficiency, usability, and comprehensibility. The J-Express system provides a powerful and integrated platform for the analysis of microarray gene expression data. It is platform-independent in that it requires only the availability of a Java virtual machine on the system. The system includes a range of analysis tools and a project management system supporting the organization and documentation of an analysis project. This unit describes the J-Express tool, emphasizing central concepts and principles, and gives examples of how it can be used to explore gene expression data sets.

The influence of the NOD Nss1/Idd5 loci on sialadenitis and gene expression in salivary glands of congenic mice.

The nonobese diabetic (NOD) Nss1 and Idd5 loci have been associated with sialadenitis development in mice. In this study the NOD Nss1 and Idd5 loci were backcrossed onto the healthy control strain B10.Q by using the speed congenic breeding strategy, resulting in three congenic strains: B10.Q.Nss1, B10.Q.Nss1/Idd5 heterozygous and B10.Q.Nss1/Idd5 homozygous. We investigated the effects of the Nss1 and Idd5 loci on sialadenitis and gene expression in NOD congenic mice. One submandibular salivary gland from each mouse was used for histological analysis of sialadenitis, whereas the contralateral salivary gland was used for gene expression profiling with the Applied Biosystems Mouse Genome Survey chip v.1.0. The results were validated using quantitative reverse transcriptase PCR. The NOD Nss1 and Idd5 loci had clear influence on the onset and progression of sialadenitis in congenic mice. Double congenic mice exhibited the most severe phenotype. We successfully identified several genes that are located in the NOD congenic regions to be differentially expressed between the congenic strains and the control strain. Several of these were found to be co-regulated, such as Stat1, complement component C1q genes and Tlr12. Also, a vast contingency of interferon-regulated genes (such as Ltb, Irf7 and Irf8) and cytokine and chemokine genes (such as Ccr7 and Ccl19) were differentially expressed between the congenic strains and the control strain. Over-representation of inflammatory signalling pathways was observed among the differentially expressed genes. We have found that the introgression of the NOD loci Nss1 and Idd5 on a healthy background caused sialadenitis in NOD congenic mouse strains, and we propose that genes within these loci are important factors in the pathogenesis. Furthermore, gene expression profiling has revealed several differentially expressed genes within and outside the NOD loci that are similar to genes found to be differentially expressed in patients with Sjögren's syndrome, and as such are interesting candidates for investigation to enhance our understanding of disease mechanisms and to develop future therapies.