Immunogenomic analysis of insect bite hypersensitivity in a model horse population.

Equine insect bite hypersensitivity (IBH) is a seasonal IgE-mediated dermatosis caused by bites of insects of the genus Culicoides. A familial predisposition for the disease has been shown but, except for the MHC, the genes involved have not been identified so far. An immunogenomic analysis of IBH was performed in a model population of Old Kladruby horses, all living in the same environment. Clinical signs of IBH were used as phenotypic manifestation of IBH. Furthermore, total serum IgE levels were determined in the sera of these horses and used as an independent phenotypic marker for the immunogenetic analysis. Single nucleotide polymorphisms (SNPs) in candidate immunity-related genes were used for association analyses. Genotypes composed of two to five genes encoding interferon gamma -IFNG, transforming growth factor beta 1 -TGFB1, Janus kinase 2 -JAK2, thymic stromal lymphopoietin -TSLP, and involucrin -IVL were associated with IBH, indicating a role of the genes in the pathogenesis of IBH. These findings were supported by analysis of gene expression in skin biopsies of 15 affected and 15 unaffected horses. Two markers associated with IBH, IFNG and TGFB1, showed differences in mRNA expression in skin biopsies from IBH-affected and non-affected horses (p<0.05). Expression of the gene coding for the CD14 receptor molecule -CD14 was different in skin biopsies at p<0.06. When total IgE levels were treated as binary traits, genotypes of IGHE, ELA-DRA, and IL10/b were associated with this trait. When treated as a continuous trait, total IgE levels were associated with genes IGHE, FCER1A, IL4, IL4R, IL10, IL1RA, and JAK2. This first report on non-MHC genes associated with IBH in horses is thus supported by differences in expression of genes known to play a role in allergy and immunity.

The additional facet of immunoscore: immunoprofiling as a possible predictive tool for cancer treatment

Recent investigations of the tumor microenvironment have shown that many tumors are infiltrated by inflammatory and lymphocytic cells. Increasing evidence suggests that the number, type and location of these tumor-infiltrating lymphocytes in primary tumors has prognostic value, and this has led to the development of an 'immunoscore. As well as providing useful prognostic information, the immunoscore concept also has the potential to help predict response to treatment, thereby improving decision- making with regard to choice of therapy. This predictive aspect of the tumor microenvironment forms the basis for the concept of immunoprofiling, which can be described as 'using an individual's immune system signature (or profile) to predict that patient's response to therapy' The immunoprofile of an individual can be genetically determined or tumor-induced (and therefore dynamic). Ipilimumab is the first in a series of immunomodulating antibodies and has been shown to be associated with improved overall survival in patients with advanced melanoma. Other immunotherapies in development include anti-programmed death 1 protein (nivolumab), anti-PD-ligand 1, anti-CD137 (urelumab), and anti-OX40. Biomarkers that can be used as predictive factors for these treatments have not yet been clinically validated. However, there is already evidence that the tumor microenvironment can have a predictive role, with clinical activity of ipilimumab related to high baseline expression of the immune-related genes FoxP3 and indoleamine 2,3-dioxygenase and an increase in tumor-infiltrating lymphocytes. These biomarkers could represent the first potential proposal for an immunoprofiling panel in patients for whom anti-CTLA-4 therapy is being considered, although prospective data are required. In conclusion, the evaluation of systemic and local immunological biomarkers could offer useful prognostic information and facilitate clinical decision making. The challenge will be to identify the individual immunoprofile of each patient and the consequent choice of optimal therapy or combination of therapies to be used.

The Pyruvate decarboxylase1 gene of Arabidopsis is required during anoxia but not other environmental stresses

Ethanolic fermentation is classically associated with flooding tolerance when plant cells switch from respiration to anaerobic fermentation. However, recent studies have suggested that fermentation also has important functions in the presence of oxygen, mainly in germinating pollen and during abiotic stress. Pyruvate decarboxylase (PDC), which catalyzes the first step in this pathway, is thought to be the main regulatory enzyme. Here, we characterize the PDC gene family in Arabidopsis. PDC is encoded by four closely related genes. By using real-time quantitative polymerase chain reaction, we determined the expression levels of each individual gene in different tissues, under normal growth conditions, and when the plants were subjected to anoxia or other environmental stress conditions. We show that PDC1 is the only gene induced under oxygen limitation among the PDC1 gene family and that a pdc1 null mutant is comprised in anoxia tolerance but not other environmental stresses. We also characterize the expression of the aldehyde dehydrogenase (ALDH) gene family. None of the three genes is induced by anoxia but ALDH2B7 reacts strongly to ABA application and dehydration, suggesting that ALDH may play a role in aerobic detoxification of acetaldehyde. We discuss the possible role of ethanolic fermentation as a robust back-up energy production pathway under adverse conditions when mitochondrial function is disturbed.

Comparison against 186 canid whole genome sequences reveals survival strategies of an ancient clonally transmissible canine tumor.

Canine transmissible venereal tumor (CTVT) is a parasitic cancer clone that has propagated for thousands of years via sexual transfer of malignant cells. Little is understood about the mechanisms that converted an ancient tumor into the world's oldest known continuously propagating somatic cell lineage. We created the largest existing catalog of canine genome-wide variation and compared it against two CTVT genome sequences, thereby separating alleles derived from the founder's genome from somatic drivers of clonal transmissibility. We show that CTVT has undergone continuous adaptation to its transmissible allograft niche, with overlapping mutations at every step of immunosurveillance, particularly self-antigen presentation and apoptosis. We also identified chronologically early somatic mutations in oncogenesis- and immune-related genes that may represent key initiators of clonal transmissibility. Thus, we provide the first insights into the specific genomic aberrations that underlie CTVT's dogged perseverance in canids around the world.

The transporter GAT1 plays an important role in GABA-mediated carbon-nitrogen interactions in Arabidopsis

Glutamate derived γ-aminobutyric acid (GABA) is synthetized in the cytosol prior to delivery to the mitochondria where it is catabolized via the TCA cycle. GABA accumulates under various environmental conditions, but an increasing number of studies show its involvement at the crossroad between C and N metabolism. To assess the role of GABA in modulating cellular metabolism, we exposed seedlings of A. thaliana GABA transporter gat1 mutant to full nutrition medium and media deficient in C and N combined with feeding of different concentrations (0.5 and 1 mM) of exogenous GABA. GC-MS based metabolite profiling showed an expected effect of medium composition on the seedlings metabolism of mutant and wild type alike. That being said, a significant interaction between GAT1 deficiency and medium composition was determined with respect to magnitude of change in relative amino acid levels. The effect of exogenous GABA treatment on metabolism was contingent on both the medium and the genotype, leading for instance to a drop in asparagine under full nutrition and low C conditions and glucose under all tested media, but not to changes in GABA content. We additionally assessed the effect of GAT1 deficiency on the expression of glutamate metabolism related genes and genes involved in abiotic stress responses. These results suggest a role for GAT1 in GABA-mediated metabolic alterations in the context of the C-N equilibrium of plant cells.

Signal signature of aboveground-induced resistance upon belowground herbivory in maize

Plants activate local and systemic defence mechanisms upon exposure to stress. This innate immune response is partially regulated by plant hormones, and involves the accumulation of defensive metabolites. Although local defence reactions to herbivores are well studied, less is known about the impact of root herbivory on shoot defence. Here, we examined the effects of belowground infestation by the western corn rootworm Diabrotica virgifera virgifera on aboveground resistance in maize. Belowground herbivory by D. v. virgifera induced aboveground resistance against the generalist herbivore Spodoptera littoralis, and the necrotrophic pathogen Setosphaeria turcica. Furthermore, D. v. virgifera increased shoot levels of 2,4-dihydroxy-7-methoxy-1,4-benzoxazin-3-one (DIMBOA), and primed the induction of chlorogenic acid upon subsequent infestation by S. littoralis. To gain insight into the signalling network behind this below- and aboveground defence interaction, we compiled a set of 32 defence-related genes, which can be used as transcriptional marker systems to detect activities of different hormone-response pathways. Belowground attack by D. v. virgifera triggered an ABA-inducible transcription pattern in the shoot. The quantification of defence hormones showed a local increase in the production of oxylipins after root and shoot infestation by D. v. virgifera and S. littoralis, respectively. On the other hand, ABA accumulated locally and systemically upon belowground attack by D. v. virgifera. Furthermore, D. v. virgifera reduced the aboveground water content, whereas the removal of similar quantities of root biomass had no effect. Our study shows that root herbivory by D. v. virgifera specifically alters the aboveground defence status of a maize, and suggests that ABA plays a role in the signalling network mediating this interaction.

In vitro analysis of osteogenic inhibition in non-union spinal fusion caused by nucleus pulposus cells.

Discectomy and spinal fusion is the gold standard for spinal surgery to relieve pain. However, fusion can be hindered for yet unknown reasons that lead to non-fusions with pseudo-arthrosis. Clinical observations indicate that presence of residual intervertebral disc (IVD) tissue might hinder the ossification. We hypothesize that BMP-antagonists are constantly secreted by IVD cells and potentially prevent the ossification process. Furthermore, L51P, the engineered BMP2 variant, stimulates osseo-induction of bone marrow-derived mesenchymal stem cells (MSC) by antagonizing BMP-inhibitors. Human MSCs, primary nucleus pulposus (NPC) and annulus pulposus cells (AFC) were isolated and expanded in monolayer cultures up to passage 3. IVD cells were seeded in 1.2% alginate beads (4Mio/mL) and separated by culture inserts from MSCs. MSCs were kept in 1:control medium, 2:osteogenic medium±alginate beads, 3:osteogenic medium+NPC (±L51P) and 4:osteogenic medium+AFC (±L51P) for 21 days. Relative gene expression of bone-related genes, alkaline phosphatase assay and histological staining were performed. Osteogenesis of MSCs was hindered as shown by reduced alizarin red staining in the presence of NPC. No such inhibition was observed if co-cultured with alginate only or in the presence of AFC. The results were confirmed on the RNA and protein level. Addition of L51Pto the co- cultures, however, induced mineralization of MSCs in presence of NPC. We demonstrated that NPC secrete BMP-antagonists that prevent osteogenesis of MSCs and L51P can antagonize BMP-antagonists and induce bone formation.

Investigation of bone inhibition in non-union spinal fusion caused by intervertebral disc cells in vitro

Introduction: Discectomy and spinal fusion is the gold standard for spinal surgery to relieve pain. However, fusion can be hindered for yet unknown reasons that lead to non-fusions with pseudo-arthrose. It is hence appealing to develop biomaterials that can enhance bone formation. Clinical observations indicate that presence of residual intervertebral disc (IVD) tissue might hinder the ossification. We hypothesize that BMP-antagonists are constantly secreted by IVD cells and potentially prevent the ossification process. Furthermore, L51P, the engineered BMP2 variant, stimulates osteoinduction of bone marrow-derived mesenchymal stem cells (MSC) by antagonizing BMP-inhibitors. Methods: Human MSCs, primary nucleus pulposus (NPC) and annulus pulposus cells (AFC) were isolated and expanded in monolayer cultures up to passage 3. IVD cells were seeded in 1.2% alginate beads (4Mio/mL) and separated by culture inserts from MSCs in a co-culture set-up. MSCs were kept in 1:control medium, 2:osteogenic medium+alginate control, 3:osteogenic medium+NPC (±L51P) and 4:osteogenic medium+AFC (±L51P) for 21 days. Relative gene expression of bone-related genes, Alkaline Phosphatase (ALP) assay and histological staining were performed. Results: Osteogenesis of MSCs was hindered as shown by reduced alizarin red staining in the presence of NPC. No such inhibition was observed if co-cultured with alginate only or in the presence of AFC. The results were confirmed on the RNA and protein level. Addition of L51P to the co-cultures induced mineralization of MSCs, however a reduced ALP was observed. Conclusion: We demonstrated that NPC secrete BMP-antagonists that prevent osteogenesis of MSCs and L51P can antagonize BMP-antagonists and induce bone formation.

Osteogenic differentiation of bone marrow stromal cells is hindered by the presence of intervertebral disc cells.

BACKGROUND Clinical observations indicate that the presence of nucleus pulposus (NP) tissue during spinal fusion hinders the rate of disc ossification. While the underlying mechanism remains unknown, this observation could be due to incomplete removal of NP cells (NPCs) that secrete factors preventing disc calcification, such as bone morphogenetic protein (BMP) antagonists including noggin and members of the DAN (differential screening selected gene aberrative in neuroblastoma) family. METHODS Monolayer human bone marrow-derived mesenchymal stem cells (MSCs) were cocultured withNPCs and annulus fibrosus cells (AFCs) embedded in alginate for 21 days. At the end of coculture, MSCs were stained for mineral deposition by alizarin red, and relative expression of bone-related genes [Runt-related transcription factor 2, (RUNX2), Osteopontin (OPN), and Alkaline phosphatase (ALP)] and ALP activity were analyzed. Relative expression of three BMP antagonists, chordin (CHRD), gremlin (GREM1), and noggin (NOG), was determined in primary human NPCs and AFCs. These cells were also stained for Gremlin and Noggin by immunocytochemistry. RESULTS Alizarin red staining showed that MSC osteogenesis in monolayer cultures was inhibited by coculture with NPCs or AFCs. ALP activity and RT-PCR analyses confirmed these results and demonstrated inhibition of osteogenesis of MSC in the presence of disc cells. NOG was significantly up-regulated in MSCs after coculture. Relative gene expression of intervertebral disc (IVD) cells showed higher expression of GREM1 in NPCs than in AFCs. CONCLUSIONS We show that primary IVD cells inhibit osteogenesis of MSCs. BMP inhibitors NOG, GREM1 and CHRD were expressed in IVD cells. GREM1 appears to be differentially expressed in NPCs and AFCs. Our results have implications for the design and development of treatments for non-union in spinal fusion.

Thyroid disruption in zebrafish (Danio rerio) larvae: Different molecular response patterns lead to impaired eye development and visual functions

The vertebrate thyroid system is important for multiple developmental processes, including eye development. Thus, its environmentally induced disruption may impact important fitness-related parameters like visual capacities and behaviour. The present study investigated the relation between molecular effects of thyroid disruption and morphological and physiological changes of eye development in zebrafish (Danio rerio). Two test compounds representing different molecular modes of thyroid disruption were used: propylthiouracil (PTU), which is an enzyme-inhibitor of thyroid hormone synthesis, and tetrabromobisphenol A (TBBPA), which interacts with the thyroid hormone receptors. Both chemicals significantly altered transcript levels of thyroid system-related genes (TRα, TRβ, TPO, TSH, DIO1, DIO2 and DIO3) in a compound-specific way. Despite these different molecular response patterns, both treatments resulted in similar pathological alterations of the eyes such as reduced size, RPE cell diameter and pigmentation, which were concentration-dependent. The morphological changes translated into impaired visual performance of the larvae: the optokinetic response was significantly and concentration-dependently decreased in both treatments, together with a significant increase of light preference of PTU-treated larvae. In addition, swimming activity was impacted. This study provides first evidence that different modes of molecular action of the thyroid disruptors can be associated with uniform apical responses. Furthermore, this study is the first to show that pathological eye development, as it can be induced by exposure to thyroid disruptors, indeed translates into impaired visual capacities of zebrafish early life stages.

Gene expression profile in newborn rat lungs after two days of recovery of mechanical ventilation.

BACKGROUND Preterm infants having immature lungs often require respiratory support, potentially leading to bronchopulmonary dysplasia (BPD). Conventional BPD rodent models based on mechanical ventilation (MV) present outcome measured at the end of the ventilation period. A reversible intubation and ventilation model in newborn rats recently allowed discovering that different sets of genes modified their expression related to time after MV. In a newborn rat model, the expression profile 48 h after MV was analyzed with gene arrays to detect potentially interesting candidates with an impact on BPD development. METHODS Rat pups were injected P4-5 with 2 mg/kg lipopolysaccharide (LPS). One day later, MV with 21 or 60% oxygen was applied during 6 h. Animals were sacrified 48 h after end of ventilation. Affymetrix gene arrays assessed the total gene expression profile in lung tissue. RESULTS In fully treated animals (LPS + MV + 60% O(2)) vs. controls, 271 genes changed expression significantly. All modified genes could be classified in six pathways: tissue remodeling/wound repair, immune system and inflammatory response, hematopoiesis, vasodilatation, and oxidative stress. Major alterations were found in the MMP and complement system. CONCLUSION MMPs and complement factors play a central role in several of the pathways identified and may represent interesting targets for BPD treatment/prevention.Bronchopulmonary dysplasia (BPD) is a chronic lung disease occurring in ~30% of preterm infants born less than 30 wk of gestation (1). Its main risk factors include lung immaturity due to preterm delivery, mechanical ventilation (MV), oxygen toxicity, chorioamnionitis, and sepsis. The main feature is an arrest of alveolar and capillary formation (2). Models trying to decipher genes involved in the pathophysiology of BPD are mainly based on MV and oxygen application to young mammals with immature lungs of different species (3). In newborn rodent models, analyses of lung structure and gene and protein expression are performed for practical reasons directly at the end of MV (4,5,6). However, later appearing changes of gene expression might also have an impact on lung development and the evolution towards BPD and cannot be discovered by such models. Recently, we developed a newborn rat model of MV using an atraumatic (orotracheal) intubation technique that allows the weaning of the newborn animal off anesthesia and MV, the extubation to spontaneous breathing, and therefore allows the evaluation of effects of MV after a ventilation-free period of recovery (7). Indeed, applying this concept of atraumatic intubation by direct laryngoscopy, we recently were able to show significant differences between gene expression changes appearing directly after MV compared to those measured after a ventilation-free interval of 48 h. Immediately after MV, inflammation-related genes showed a transitory modified expression, while another set of more structurally related genes changed their expression only after a delay of 2 d (7). Lung structure, analyzed by conventional 2D histology and also by 3D reconstruction using synchrotron x-ray tomographic microscopy revealed, 48 h after end of MV, a reduced complexity of lung architecture compared to the nonventilated rat lungs, similar to the typical findings in BPD. To extend these observations about late gene expression modifications, we performed with a similar model a full gene expression profile of lung tissue 48 h after the end of MV with either room air or 60% oxygen. Essentially, we measured changes in the expression of genes related to the MMPs and complement system which played a role in many of the six identified mostly affected pathways.

Functional studies of insulin-like growth factor binding protein 2 pathway in glioma progression

Comparison of gene expressing profiles between gliomas with different grades revealed frequent overexpression of insulin-like growth factor binding protein 2 (IGFBP2) in glioblastomas (GBM), in which uncontrolled cell proliferation, angiogenesis, invasion and anti-apoptosis are hallmarks. Using the glia-specific gene transfer transgenic mouse and the stable LN229(BP2) GBM cell lines, we found that IGFBP2 by itself cannot transform cells in vitro and in vivo. IGFBP2 had growth inhibitory effects on mouse primary neural progenitors, but overexpression of IGFBP2 had no effect on GBM cells. ^ Although IGFBP2 does not initiate gliomagenesis, using tissue array technology, we observed strong correlation between IGFBP2 overexpression and VEGF up-regulation in human diffuse gliomas. Furthermore, overexpression of IGFBP2 in GBM cells not only enhanced VEGF expression but also increased the malignant potential of U87 MG cells in our angiogenesis xenograft animal model. ^ In parallel to these studies, using established stable SNB19 GBM cells that overexpress IGFBP2, we found that IGFBP2 significantly increased invasion by induction of matrix metalloproteinase-2 (MMP-2) as well as other invasion related genes, providing evidence that IGFBP2 contributes to glioma progression in part by enhancing MMP-2 gene transcription and in turn tumor cell invasion. ^ Finally, we found that primary filial cells infected with an anti-sense IGFBP2 construct have markedly increased sensitivity to γ irradiation and reduced Akt activation. On the other hand, SNB19(BP2) stable lines have consistently increased levels of Akt and NFkB activation, suggesting that one possible mechanism for anti-apoptosic function of IGFBP2 is through the activation of Akt and NFkB. Beside this, what is especially interesting is the finding that Akt protein was cleaved and inactivated during apoptosis by caspases, and IGFBP2 can prevent Akt cleavage, revealing another possible mechanism through it IGFBP2 exhibit strong antiapoptotic effects. Our data showed that IGFBP2 is a specific substrate for caspase-3, raising the possibility that IGFBP2 may inhibit apoptosis by a suicide mechanism. ^ In summary, using cellular, genomics, and molecular approaches, this thesis documented the potential roles of IGFBP2 in glioma progression. Our findings shed light on an important biological aspect of glioma progression and may provide new insights useful for the design of novel mechanism-based therapies for GBM. ^

Reactive oxygen species-dependent mechanisms of N-(4-hydroxyphenyl)retinamide (4HPR)-induced apoptosis in human carcinoma cells

4HPR is a synthetic retinoid that has shown chemopreventive and therapeutic efficacy against premalignant and malignant lesions including oral leukoplakia, ovarian and breast cancer, and neuroblastoma. 4HPR induces apoptosis in various cancer cells and production of reactive oxygen species (ROS) has been suggested as a possible cause underlying these effects. However, the mechanisms governing these effects by 4HPR are not fully elucidated. In this study, we explored the mechanisms of 4HPR-induced ROS increase and apoptosis in human cancer cells. ^ First, we identified genes modulated by 4HPR using oligonucleotide gene expression arrays and found that they fall into specific functional canonical pathways and gene networks using Ingenuity Pathways Analysis®. Further analysis has shown that 4HPR induced up-regulation of Endoplasmic Reticulum (ER)-related genes such as Heat shock proteins 70 and 90 and the transcriptional factor, GADD153. These findings were validated using quantitative real-time PCR. ^ Second, we found that 4HPR induced extensive ER stress evidenced by dilation of the ER and endoribonuclease-mediated splicing and activation of the transcriptional factor, XBP-1. In addition, 4HPR induced the up-regulation of various ER stress-related genes and their protein products, as well as cleavage and activation of the ER specific Caspase-4. Concomitantly with XBP-1 splicing, all of these effects were dependent on ROS generation by 4HPR. Furthermore, chemical inhibition and RNA interference studies revealed a novel pro-apoptotic role for HSP70/A1A in 4HPR-mediated apoptosis. ^ Third, we observed rapid activation of the small GTPase Rac by 4HPR which was upstream of ROS generation. Inhibition of Rac activity or silencing of its expression by RNA interference inhibited ROS generation and apoptosis induction by 4HPR. siRNA targeting PAK1 and expression of a dominant negative Rac, decreased 4HPR-mediated ROS generation, while expression of a constitutive active Rac increased basal and 4HPR-induced ROS generation and PARP cleavage. Furthermore, metastatic cancer cells exhibited higher Rac activation, ROS generation, and cell growth inhibition due to 4HPR exposure compared to their primary cancer cell counterparts. ^ These findings provide novel insights into 4HPR-mediated ROS generation and apoptosis induction and support the use of ROS inducing agents such as 4HPR against metastatic cancer cells. ^

Immune and genetic risk factors in glioma

Despite extensive research, the etiology of adult glioma remains largely unknown. We sought to further explore the role of immune and genetic factors in glioma etiology using data from the Harris County Brain Tumor Study and the first U.S. genome-wide association study of glioma. First, using a case-control study design, we examined the association between adult glioma risk and surrogates of the timing and frequency of common early childhood infections, birth order and sibship size, respectively. We found that each one-unit increase in birth order was associated with a 12% decreased risk of glioma development in adulthood (OR=0.88, 95% CI=0.81-0.96); however, sibship size was not associated with adult glioma risk (OR=0.96, 95% CI=0.91-1.02). Second, we used a multi-strategic approach to explore the relationships between glioma risk, history of asthma/allergies, and 23 functional SNPs in 11 inflammation genes. We found three inflammation gene SNPs to be significantly associated with glioma risk (COX2/PTGS2 rs20417 [OR=1.41]; IL10 rs1800896 [OR=1.57]; and IL13 rs20541 [OR=0.39]). Joint effects analysis of the risk-conferring alleles of these three SNPs revealed a trend of increasing risk with increasing number of adverse alleles among those without asthma/allergies (p<0.0001). Finally, we conducted a case-only study to explore pairwise SNP-SNP interactions in immune-related pathways among a population of 1304 non-Hispanic white glioma cases. After correction for multiple comparisons, we found 279 significant SNP-SNP interactions among polymorphisms of immune-related genes, many of which have not been previously examined. Our results, cumulatively, suggest an important role for immune and genetic factors in glioma etiology and provide several new hypotheses for future studies.^

Network-based genome-wide association study of Crohn's disease

Genome-wide association studies (GWAS) have rapidly become a standard method for disease gene discovery. Many recent GWAS indicate that for most disorders, only a few common variants are implicated and the associated SNPs explain only a small fraction of the genetic risk. The current study incorporated gene network information into gene-based analysis of GWAS data for Crohn's disease (CD). The purpose was to develop statistical models to boost the power of identifying disease-associated genes and gene subnetworks by maximizing the use of existing biological knowledge from multiple sources. The results revealed that Markov random field (MRF) based mixture model incorporating direct neighborhood information from a single gene network is not efficient in identifying CD-related genes based on the GWAS data. The incorporation of solely direct neighborhood information might lead to the low efficiency of these models. Alternative MRF models looking beyond direct neighboring information are necessary to be developed in the future for the purpose of this study.^