Fat intake interacts with polymorphisms of Caspase9, FasLigand and PPARgamma apoptotic genes in modulating Crohn's disease activity

Background & aims: Crohn’s disease (CD) is a multifactorial disease where resistance to apoptosis is one major defect. Also, dietary fat intake has been shown to modulate disease activity. We aimed to explore the interaction between four single nucleotide polymorphisms (SNPs) in apoptotic genes and dietary fat intake in modulating disease activity in CD patients. Methods: Polymerase Chain Reaction (PCR) and Restriction Fragment Length Polymorphism (RFLP) techniques were used to analyze Caspase9þ93C/T, FasLigand-843C/T, Peroxisome Proliferator-Activated Receptor gammaþ161C/T and Peroxisome Proliferator-Activated Receptor gamma Pro12Ala SNPs in 99 patients with CD and 116 healthy controls. Interactions between SNPs and fat intake in modulating disease activity were analyzed using regression analysis. Results: None of the polymorphisms analyzed influenced disease susceptibility and/or activity, but a high intake of total, saturated and monounsaturated fats and a higher ratio of n-6/n-3 polyunsaturated fatty acids (PUFA), was associated with a more active phenotype (p < 0.05). We observed that the detrimental effect of a high intake of total and trans fat was more marked in wild type carriers of the Caspase9þ93C/T polymorphism [O.R (95%CI) 4.64 (1.27e16.89) and O.R (95%CI) 4.84 (1.34e17.50)]. In the Peroxisome Proliferator-Activated Receptor gamma Pro12Ala SNP, we also observed that a high intake of saturated and monounsaturated fat was associated to a more active disease in wild type carriers [OR (95%CI) 4.21 (1.33e13.26) and 4.37 (1.52e12.51)]. Finally, a high intake of n-6 PUFA was associated with a more active disease in wild type carriers for the FasLigand-843C/T polymorphism [O.R (95%CI) 5.15 (1.07e24.74)]. Conclusions: To our knowledge, this is the first study to disclose a synergism between fat intake and SNPs in apoptotic genes in modulating disease activity in CD patients.

Relationship of PON1 192 and 55 gene polymorphisms to calcific valvular aortic stenosis

Introduction and Objectives - Paraoxonases may exert anti-atherogenic action by reducing lipid peroxidation. Previous studies examined associations between polymorphisms in the paraoxonase 1 (PON1) gene and development of coronary artery disease (CAD), with inconsistent results. Given the similarities in clinical and pathophysiological risk factors of CAD and calcific aortic valve stenosis (CAVS), we postulated a link between PON1 alleles and CAVS progression. Methods - We investigated the association between PON1 55 and 192 single nucleotide polymorphisms (SNPs), their enzyme activity, and CAVS progression assessed by aortic valve area and transvalvular peak velocity in 67 consecutive patients with moderate CAVS and 251 healthy controls. Results - PON1 paraoxonase activity was higher in CAVS patients (P<0.001). The PON1 genotype Q192R SNP (P=0.03) and variant allele (R192) (P=0.01) frequencies differed between CAVS patients and controls. Significant association existed between PON1 enzyme activity, phenotypic effects of PON1 192 genotype polymorphisms, and CAVS progression, but not between PON1 55 and high-density lipoprotein (P=0.44) or low-density lipoprotein cholesterol (P=0.12), between 192 genotype and high-density lipoprotein (P=0.24) or low-density lipoprotein cholesterol (P=0.52). Conclusion - The PON1 genotype Q192R SNP has an important effect on CAVS disease progression. This study helps outline a genotype-phenotype relationship for PON1 in this unique population.

Clinical and genetic factors predicting response to therapy in patients with Crohn’s disease

Aim - To identify clinical and/or genetic predictors of response to several therapies in Crohn’s disease (CD) patients. Methods - We included 242 patients with CD (133 females) aged (mean ± standard deviation) 39 ± 12 years and a disease duration of 12 ± 8 years. The single-nucleotide polymorphisms (SNPs) studied were ABCB1 C3435T and G2677T/A, IL23R G1142A, C2370A, and G9T, CASP9 C93T, Fas G670A and LgC844T, and ATG16L1 A898G. Genotyping was performed with real-time PCR with Taqman probes. Results - Older patients responded better to 5-aminosalicylic acid (5-ASA) and to azathioprine (OR 1.07, p = 0.003 and OR 1.03, p = 0.01, respectively) while younger ones responded better to biologicals (OR 0.95, p = 0.06). Previous surgery negatively influenced response to 5-ASA compounds (OR 0.25, p = 0.05), but favoured response to azathioprine (OR 2.1, p = 0.04). In respect to genetic predictors, we observed that heterozygotes for ATGL16L1 SNP had a significantly higher chance of responding to corticosteroids (OR 2.51, p = 0.04), while homozygotes for Casp9 C93T SNP had a lower chance of responding both to corticosteroids and to azathioprine (OR 0.23, p = 0.03 and OR 0.08, p = 0.02,). TT carriers of ABCB1 C3435T SNP had a higher chance of responding to azathioprine (OR 2.38, p = 0.01), while carriers of ABCB1 G2677T/A SNP, as well as responding better to azathioprine (OR 1.89, p = 0.07), had a lower chance of responding to biologicals (OR 0.31, p = 0.07), which became significant after adjusting for gender (OR 0.75, p = 0.005). Conclusions - In the present study, we were able to identify a number of clinical and genetic predictors of response to several therapies which may become of potential utility in clinical practice. These are preliminary results that need to be replicated in future pharmacogenomic studies.

Common genetic polymorphisms in the ABCB1 gene are associated with risk of major depressive disorder in male Portuguese individuals

Major depressive disorder (MDD) is a highly prevalent disorder, which has been associated with an abnormal response of the hypothalamus–pituitary–adrenal (HPA) axis. Reports have argued that an abnormal HPA axis response can be due to an altered P-Glycoprotein (P-GP) function. This argument suggests that genetic polymorphisms in ABCB1 may have an effect on the HPA axis activity; however, it is still not clear if this influences the risk of MDD. Our study aims to evaluate the effect of ABCB1 C1236T, G2677TA and C3435T genetic polymorphisms on MDD risk in a subset of Portuguese patients. DNA samples from 80 MDD patients and 160 control subjects were genotyped using TaqMan SNP Genotyping assays. A significant protection for MDD males carrying the T allele was observed (C1236T: odds ratio (OR) = 0.360, 95% confidence interval [CI]: [0.140– 0.950], p = 0.022; C3435T: OR= 0.306, 95% CI: [0.096–0.980], p = 0.042; and G2677TA: OR= 0.300, 95% CI: [0.100– 0.870], p = 0.013). Male Portuguese individuals carrying the 1236T/2677T/3435T haplotype had nearly 70% less risk of developing MDD (OR = 0.313, 95% CI: [0.118–0.832], p = 0.016, FDR p = 0.032). No significant differences were observed regarding the overall subjects. Our results suggest that genetic variability of the ABCB1 is associated with MDD development in male Portuguese patients. To the best of our knowledge, this is the first report in Caucasian samples to analyze the effect of these ABCB1 genetic polymorphisms on MDD risk.

Deregulated expression of selected histone methylases and demethylases in prostate carcinoma

Prostate cancer (PCa), a leading cause of cancer-related morbidity and mortality, arises through the acquisition of genetic and epigenetic alterations. Deregulation of histone methyltransferases (HMTs) or demethylases (HDMs) has been associated with PCa development and progression. However, the precise influence of altered HMTs or HDMs expression and respective histone marks in PCa onset and progression remains largely unknown. To clarify the role of HMTs and HDMs in prostate carcinogenesis, expression levels of 37 HMTs and 20 HDMs were assessed in normal prostate and PCa tissue samples by RT-qPCR. SMYD3, SUV39H2, PRMT6, KDM5A, and KDM6A were upregulated, whereas KMT2A-E (MLL1-5) and KDM4B were downregulated in PCa, compared with normal prostate tissues. Remarkably, PRMT6 was the histone modifier that best discriminated normal from tumorous tissue samples. Interestingly, EZH2 and SMYD3 expression levels significantly correlated with less differentiated and more aggressive tumors. Remarkably, SMYD3 expression levels were of independent prognostic value for the prediction of disease-specific survival of PCa patients with clinically localized disease submitted to radical prostatectomy. We concluded that expression profiling of HMTs and HDMs, especially SMYD3, might be of clinical usefulness for the assessment of PCa patients and assist in pre-therapeutic decision-making.

An unsupervised approach to feature discretization and selection

Many learning problems require handling high dimensional datasets with a relatively small number of instances. Learning algorithms are thus confronted with the curse of dimensionality, and need to address it in order to be effective. Examples of these types of data include the bag-of-words representation in text classification problems and gene expression data for tumor detection/classification. Usually, among the high number of features characterizing the instances, many may be irrelevant (or even detrimental) for the learning tasks. It is thus clear that there is a need for adequate techniques for feature representation, reduction, and selection, to improve both the classification accuracy and the memory requirements. In this paper, we propose combined unsupervised feature discretization and feature selection techniques, suitable for medium and high-dimensional datasets. The experimental results on several standard datasets, with both sparse and dense features, show the efficiency of the proposed techniques as well as improvements over previous related techniques.

DNA nanoprobes for molecular detection

Dissertação apresentada para obtenção do Grau de Doutor em Engenharia Biológica – especialidade Engenharia Genética, pela Universidade Nova de Lisboa, Faculdade de Ciências e Tecnologia

Susceptibilidade genética para cancro da mama

RESUMO: O cancro da mama é a patologia oncológica mais frequente nas mulheres sendo o responsável pela maior taxa de mortalidade por cancro no sexo feminino. Contudo, as causas inerentes a esta patologia permanecem por esclarecer. Nos últimos anos tem-se verificado que o risco para patologia neoplásica depende de factores ambientais e genéticos, estando estes últimos associados à variabilidade genética inter-individual. Polimorfismos genéticos em genes envolvidos no metabolismo de hormonas sexuais, de cancerígenos ambientais e na reparação da lesão genética, são potenciais candidatos a estarem associados à susceptibilidade individual para esta patologia. Assim, neste trabalho desenvolveram-se estudos de associação caso-controlo na população Portuguesa, com vista a avaliar-se o papel atribuído aos polimorfismos na susceptibilidade para cancro da mama. Foram seleccionados polimorfismos em genes envolvidos em diferentes vias mecanicistas: destoxificação de cancerígenos, metabolismo de estrogénios, reparação por excisão de bases, reparação por excisão de nucleótidos, reparação mismatch e reparação por recombinação homóloga. Os resultados obtidos revelaram associação entre os seguintes polimorfismos e a susceptibilidade individual para cancro da mama: os dois SNPs estudados no gene XRCC1 (Arg194Trp e Arg399Gln) e o SNP no gene XRCC3 (Thr241Met) após estratificação pelo status menopausico. Mediante estratificação por status de amamentação os SNPs identificados nos genes MnSOD (Val16Ala) e XRCC2 (Arg118His); um SNP no gene MLH3 (Leu844Pro), e por fim como resultado de interacção gene-gene as interacções descritas por MSH3 Ala1045Thr/MSH6 Gly39Glu e MSH4 Ala97Thr/MLH3 Leu844Pro. Os resultados obtidos e apresentados na presente dissertação, revelam que o estudo de polimorfismos pode representar um papel determinante na etiologia do cancro da mama. No entanto, mais estudos envolvendo estes mesmos polimorfismos em populações casuisticamente superiores serão uma mais-valia nos estudos de associação para esta neoplasia. Adicionalmente, a utilização da metodologia de Pools de DNA, poderá ser uma ferramenta útil na pré-selecção dos polimorfismos mais relevantes a estudar, na medida em que permite estimar a frequência alélica de cada SNP numa determinada população.-----------------------------------ABSTRACT: Breast cancer is the most common form of cancer among women, being the responsible for the highest mortality rate from cancer among the female sex. However, the main causes related to this pathology remain unclear. The risk of neoplasic disease has been connected with genetic and environmental factors. In fact, genes and the environment share the stage for most, if not all, common non-familial cancers, and are related to individual susceptibility. Genetic polymorphisms identified in genes encoding enzymes involved in estrogen metabolism, xenobiotics and DNA repair pathways are believed to be candidates for associations with breast cancer. Therefore, it was our intention to develop case-control studies among the Portuguese population, in order to evaluate the potential role of several genetic polymorphisms in breast cancer susceptibility. We selected polymorphisms in genes involved in different pathways: carcinogenic detoxification, estrogen metabolism, base excision repair, nucleotide excision repair, mismatch repair and double strand break repair by homologous recombination. The results obtained revealed potential associations between some polymorphisms studied and individual susceptibility to breast cancer. Regarding this fact, our results suggest the potential involvement of two XRCC1 gene polymorphisms (Arg194Trp and Arg399Gln) and XRCC3 gene polymorphism (Thr241Met) after stratification to menopausal status and after stratification to breastfeeding status an association of MnSOD gene polymorphism (Val16Ala) and XRCC2 (Arg188His) with the disease. The SNP identified in MLH3 gene (Leu844Pro), and the interaction gene-gene described by MSH3 Ala1045Thr/MSH6 Gly39Glu and MSH4 Ala97Thr/MLH3 Leu844Pro were also related to breast cancer susceptibility. The results shown in the present dissertation have revealed the potential role of polymorphisms in breast cancer etiology. However, further studies will be needed with larger populations to confirm these results. Additionally, the use of DNA pools methodology, as a pre-selection tool, could allow the identification of the most relevant polymorphisms to be studied, estimating the allelic frequency of each SNPs in different populations.

mRNA genotyping by gold nanoprobes

Dissertação para obtenção do Grau de Mestre em Biotecnologia

Role of the DHH1 Gene in the Regulation of Monocarboxylic Acids Transporters Expression in Saccharomyces cerevisiae

Previous experiments revealed that DHH1, a RNA helicase involved in the regulation of mRNA stability and translation, complemented the phenotype of a Saccharomyces cerevisiae mutant affected in the expression of genes coding for monocarboxylic-acids transporters, JEN1 and ADY2 (Paiva S, Althoff S, Casal M, Leao C. FEMS Microbiol Lett, 1999, 170∶301–306). In wild type cells, JEN1 expression had been shown to be undetectable in the presence of glucose or formic acid, and induced in the presence of lactate. In this work, we show that JEN1 mRNA accumulates in a dhh1 mutant, when formic acid was used as sole carbon source. Dhh1 interacts with the decapping activator Dcp1 and with the deadenylase complex. This led to the hypothesis that JEN1 expression is post-transcriptionally regulated by Dhh1 in formic acid. Analyses of JEN1 mRNAs decay in wild-type and dhh1 mutant strains confirmed this hypothesis. In these conditions, the stabilized JEN1 mRNA was associated to polysomes but no Jen1 protein could be detected, either by measurable lactate carrier activity, Jen1-GFP fluorescence detection or western blots. These results revealed the complexity of the expression regulation of JEN1 in S. cerevisiae and evidenced the importance of DHH1 in this process. Additionally, microarray analyses of dhh1 mutant indicated that Dhh1 plays a large role in metabolic adaptation, suggesting that carbon source changes triggers a complex interplay between transcriptional and post-transcriptional effects.

Epidemiological characterization, antimicrobial resistance and virulence mechanisms of human and animal streptococci

Dissertação para obtenção do Grau de Doutor em Biologia

Ankylosing spondylitis: genomic and functional characterization of candidate genes and their repercussion in clinical practice

RESUMO: Introdução: A espondilite anquilosante (EA) é uma doença inflamatória crónica caracterizada pela inflamação das articulações sacroilíacas e da coluna. A anquilose progressiva motiva uma deterioração gradual da função física e da qualidade de vida. O diagnóstico e o tratamento precoces podem contribuir para um melhor prognóstico. Neste contexto, a identificação de biomarcadores, assume-se como sendo muito útil para a prática clínica e representa hoje um grande desafio para a comunidade científica. Objetivos: Este estudo teve como objetivos: 1 - caracterizar a EA em Portugal; 2 - investigar possíveis associações entre genes, MHC e não-MHC, com a suscetibilidade e as características fenotípicas da EA; 3 - identificar genes candidatos associados a EA através da tecnologia de microarray. Material e Métodos: Foram recrutados doentes com EA, de acordo com os critérios modificados de Nova Iorque, nas consultas de Reumatologia dos diferentes hospitais participantes. Colecionaram-se dados demográficos, clínicos e radiológicos e colhidas amostras de sangue periférico. Selecionaram-se de forma aleatória, doentes HLA-B27 positivos, os quais foram tipados em termos de HLA classe I e II por PCR-rSSOP. Os haplótipos HLA estendidos foram estimados pelo algoritmo Expectation Maximization com recurso ao software Arlequin v3.11. As variantes alélicas dos genes IL23R, ERAP1 e ANKH foram estudadas através de ensaios de discriminação alélica TaqMan. A análise de associação foi realizada utilizando testes da Cochrane-Armitage e de regressão linear, tal como implementado pelo PLINK, para variáveis qualitativas e quantitativas, respetivamente. O estudo de expressão génica foi realizado por Illumina HT-12 Whole-Genome Expression BeadChips. Os genes candidatos foram validados usando qPCR-based TaqMan Low Density Arrays (TLDAs). Resultados: Foram incluídos 369 doentes (62,3% do sexo masculino, com idade média de 45,4 ± 13,2 anos, duração média da doença de 11,4 ± 10,5 anos). No momento da avaliação, 49,9% tinham doença axial, 2,4% periférica, 40,9% mista e 7,1% entesopática. A uveíte anterior aguda (33,6%) foi a manifestação extra-articular mais comum. Foram positivos para o HLA-B27, 80,3% dos doentes. Os haplótipo A*02/B*27/Cw*02/DRB1*01/DQB1*05 parece conferir suscetibilidade para a EA, e o A*02/B*27/Cw*01/DRB1*08/DQB1*04 parece conferir proteção em termos de atividade, repercussão funcional e radiológica da doença. Três variantes (2 para IL23R e 1 para ERAP1) mostraram significativa associação com a doença, confirmando a associação destes genes com a EA na população Portuguesa. O mesmo não se verificou com as variantes estudadas do ANKH. Não se verificou associação entre as variantes génicas não-MHC e as manifestações clínicas da EA. Foi identificado um perfil de expressão génica para a EA, tendo sido validados catorze genes - alguns têm um papel bem documentado em termos de inflamação, outros no metabolismo da cartilagem e do osso. Conclusões: Foi estabelecido um perfil demográfico e clínico dos doentes com EA em Portugal. A identificação de variantes génicas e de um perfil de expressão contribuem para uma melhor compreensão da sua fisiopatologia e podem ser úteis para estabelecer modelos com relevância em termos de diagnóstico, prognóstico e orientação terapêutica dos doentes. -----------ABSTRACT: Background: Ankylosing Spondylitis (AS) is a chronic inflammatory disorder characterized by inflammation in the spine and sacroiliac joints leading to progressive joint ankylosis and in progressive deterioration of physical function and quality of life. An early diagnosis and early therapy may contribute to a better prognosis. The identification of biomarkers would be helpful and represents a great challenge for the scientific community. Objectives: The present study had the following aims: 1- to characterize the pattern of AS in Portuguese patients; 2- to investigate MHC and non-MHC gene associations with susceptibility and phenotypic features of AS and; 3- to identify candidate genes associated with AS by means of whole-genome microarray. Material and Methods: AS was defined in accordance to the modified New York criteria and AS cases were recruited from hospital outcares patient clinics. Demographic and clinical data were recorded and blood samples collected. A random group of HLA-B27 positive patients and controls were selected and typed for HLA class I and II by PCR-rSSOP. The extended HLA haplotypes were estimated by Expectation Maximization Algorithm using Arlequin v3.11 software. Genotyping of IL23R, ERAP1 and ANKH allelic variants was carried out with TaqMan allelic discrimination assays. Association analysis was performed using the Cochrane-Armitage and linear regression tests as implemented in PLINK, for dichotomous and quantitative variables, respectively. Gene expression profile was carried out using Illumina HT-12 Whole-Genome Expression BeadChips and candidate genes were validated using qPCR-based TaqMan Low Density Arrays (TLDAs). Results: A total of 369 patients (62.3% male; mean age 45.4±13.2 years; mean disease duration 11.4±10.5 years), were included. Regarding clinical disease pattern, at the time of assessment, 49.9% had axial disease, 2.4% peripheral disease, 40.9% mixed disease and 7.1% isolated enthesopathic disease. Acute anterior uveitis (33.6%) was the most common extra-articular manifestation. 80.3% of AS patients were HLA-B27 positive. The haplotype A*02/B*27/Cw*02/DRB1*01/DQB1*05 seems to confer susceptibility to AS, whereas A*02/B*27/Cw*01/DRB1*08/DQB1*04 seems to provide protection in terms of disease activity, functional and radiological repercussion. Three markers (two for IL23R and one for ERAP1) showed significant single-locus disease associations. Association of these genes with AS in the Portuguese population was confirmed, whereas ANKH markers studied did not show an association with AS. No association was seen between non-MHC genes and clinical manifestations of AS. A gene expression signature for AS was established; among the fourteen validated genes, a number of them have a well-documented inflammatory role or in modulation of cartilage and bone metabolism. Conclusions: A demographic and clinical profile of patients with AS in Portugal was established. Identification of genetic variants of target genes as well as gene expression signatures could provide a better understanding of AS pathophysiology and could be useful to establish models with relevance in terms of susceptibility, prognosis, and potential therapeutic guidance.

Differential Internalization of Amphotericin B - Conjugated Nanoparticles in Human Cells and the Expression of Heat Shock Protein 70

Although a variety of nanoparticles (NPs) functionalized with amphotericin B, an antifungal agent widely used in the clinic, have been studied in the last years their cytotoxicity profile remains elusive. Here we show that human endothelial cells take up high amounts of silica nanoparticles (SNPs) conjugated with amphotericin B (AmB) (SNP-AmB) (65.4 12.4 pg of Si per cell) through macropinocytosis while human fibroblasts internalize relatively low amounts (2.3 0.4 pg of Si per cell) because of their low capacity for macropinocytosis. We further show that concentrations of SNP-AmB and SNP up to 400 mg/mL do not substantially affect fibroblasts. In contrast, endothelial cells are sensitive to low concentrations of NPs (above 10 mg/mL), in particular to SNP-AmB. This is because of their capacity to internalize high concentration of NPs and high sensitivity of their membrane to the effects of AmB. Low-moderate concentrations of SNP-AmB (up to 100 mg/mL) induce the production of reactive oxygen species (ROS), LDH release, high expression of pro-inflammatory cytokines and chemokines (IL-8, IL-6, G-CSF, CCL4, IL-1b and CSF2) and high expression of heat shock proteins (HSPs) at gene and protein levels. High concentrations of SNP-AmB (above 100 ug/mL) disturb membrane integrity and kill rapidly human cells(60% after 5 h). This effect is higher in SNP-AmB than in SNP.

PRELIMINARY REPORT ON THE PUTATIVE ASSOCIATION OF IL10 -3575 T/A GENETIC POLYMORPHISM WITH MALARIA SYMPTOMS

Only a small percentage of individuals living in endemic areas develop severe malaria suggesting that host genetic factors may play a key role. This study has determined the frequency of single nucleotide polymorphisms (SNPs) in some pro and anti-inflammatory cytokine gene sequences: IL6 (-174; rs1800795), IL12p40 (+1188; rs3212227), IL4 (+33; rs2070874), IL10 (-3575; rs1800890) and TGFb1 (+869; rs1800470), by means of PCR-RFLP. Blood samples were collected from 104 symptomatic and 37 asymptomatic subjects. Laboratory diagnosis was assessed by the thick blood smear test and nested-PCR. No association was found between IL6 (-174), IL12p40 (+1188), IL4 (+33), IL10 (- 3575), TGFb1 (+869) SNPs and malaria symptoms. However, regarding the IL10 -3575 T/A SNP, there were significantly more AA and AT subjects, carrying the polymorphic allele A, in the symptomatic group (c2 = 4.54, p = 0.01, OR = 0.40 [95% CI - 0.17- 0.94]). When the analysis was performed by allele, the frequency of the polymorphic allele A was also significantly higher in the symptomatic group (c2 = 4.50, p = 0.01, OR = 0.45 [95% CI - 0.21-0.95]). In conclusion, this study has suggested the possibility that the IL10 - 3575 T/A SNP might be associated with the presence and maintenance of malaria symptoms in individuals living in endemic areas. Taking into account that this polymorphism is related to decreased IL10 production, a possible role of this SNP in the pathophysiology of malaria is also suggested, but replication studies with a higher number of patients and evaluation of IL10 levels are needed for confirmation.

Gene Expression Profiling and Association Studies Implicate the Neuregulin Signaling Pathway in Behçet's Disease Susceptibility

Behçet's disease (BD) is a complex disease with genetic and environmental risk factors implicated in its etiology; however, its pathophysiology is poorly understood. To decipher BD's genetic underpinnings, we combined gene expression profiling with pathway analysis and association studies. We compared the gene expression profiles in peripheral blood mononuclear cells (PBMCs) of 15 patients and 14 matched controls using Affymetrix microarrays and found that the neuregulin signaling pathway was over-represented among the differentially expressed genes. The Epiregulin (EREG), Amphiregulin (AREG), and Neuregulin-1 (NRG1) genes of this pathway stand out as they are also among the top differentially expressed genes. Twelve haplotype tagging SNPs at the EREG-AREG locus and 15 SNPs in NRG1 found associated in at least one published BD genome-wide association study were tested for association with BD in a dataset of 976 Iranian patients and 839 controls. We found a novel association with BD for the rs6845297 SNP located downstream of EREG, and replicated three associations at NRG1 (rs4489285, rs383632, and rs1462891). Multifactor dimensionality reduction analysis indicated the existence of epistatic interactions between EREG and NRG1 variants. EREG-AREG and NRG1, which are members of the epidermal growth factor (EGF) family, seem to modulate BD susceptibility through main effects and gene–gene interactions. These association findings support a role for the EGF/ErbB signaling pathway inBD pathogenesis that warrants further investigation and highlight the importance of combining genetic and genomic approaches to dissect the genetic architecture of complex diseases.