Genes and gene ontologies common to airflow obstruction and emphysema in the lungs of patients with COPD

Chronic obstructive pulmonary disease (COPD) is a major public health problem with increasing prevalence worldwide. The primary aim of this study was to identify genes and gene ontologies associated with COPD severity. Gene expression profiling was performed on total RNA extracted from lung tissue of 18 former smokers with COPD. Class comparison analysis on mild (n = 9, FEV₁ 80–110% predicted) and moderate (n = 9, FEV₁ 50–60% predicted) COPD patients identified 46 differentially expressed genes (p<0.01), of which 14 genes were technically confirmed by quantitative real-time-PCR. Biological replication in an independent test set of 58 lung samples confirmed the altered expression of ten genes with increasing COPD severity, with eight of these genes (NNMT, THBS1, HLA-DPB1, IGHD, ETS2, ELF1, PTGDS and CYRBD1) being differentially expressed by greater than 1.8 fold between mild and moderate COPD, identifying these as candidate determinants of COPD severity. These genes belonged to ontologies potentially implicated in COPD including angiogenesis, cell migration, proliferation and apoptosis. Our secondary aim was to identify gene ontologies common to airway obstruction, indicated by impaired FEV₁ and KCO. Using gene ontology enrichment analysis we have identified relevant biological and molecular processes including regulation of cell-matrix adhesion, leukocyte activation, cell and substrate adhesion, cell adhesion, angiogenesis, cell activation that are enriched among genes involved in airflow obstruction. Exploring the functional significance of these genes and their gene ontologies will provide clues to molecular changes involved in severity of COPD, which could be developed as targets for therapy or biomarkers for early diagnosis.     

Ectopic pregnancy as a model to identify endometrial genes and signaling pathways important in decidualization and regulated by local trophoblast

The endometrium in early pregnancy undergoes decidualization and functional changes induced by local trophoblast, which are not fully understood. We hypothesized that endometrium from tubal ectopic pregnancy (EP) could be interrogated to identify novel genes and pathways involved in these processes. Gestation-matched endometrium was collected from women with EP (n = 11) and intrauterine pregnancies (IUP) (n = 13). RNA was extracted from the tissue. In addition, tissues were prepared for histological analysis for degree of decidualization. We compared a) the samples from EP that were decidualized (n = 6) with non-decidualized samples (n = 5), and b) the decidualized EP (n = 6) with decidualization-matched IUP (n = 6) samples using an Affymetrix gene array platform, with Ingenuity Pathway Analysis, combined with quantitative RT-PCR. Expression of PRL and IGFBP1 was used to confirm the degree of decidualization in each group. There were no differences in PRL or IGFBP1 expression in the decidualization-matched samples but a marked reduction (P<0.001) in the non-decidualized samples. Decidualization was associated with increased expression of 428 genes including SCARA5 (181-fold), DKK1 (71-fold) and PROK1 (32-fold), and decreased expression of 230 genes including MMP-7 (35-fold) and SFRP4 (21-fold). The top canonical pathways associated with these differentially expressed genes were Natural Killer Cell and Wnt/b-Catenin signaling. Local trophoblast was associated with much less alteration of endometrial gene expression with an increase in 56 genes, including CSH1 (8-fold), and a reduction in 29 genes including CRISP3 (8-fold). The top associated canonical pathway was Antigen Presentation. The study of endometrium from tubal EP may promote novel insights into genes involved in decidualization and those influenced by factors from neighboring trophoblast. This has afforded unique information not highlighted by previous studies and adds to our understanding of the endometrium in early pregnancy.     

Genetic association study identifies HSPB7 as a risk gene for idiopathic dilated cardiomyopathy

Dilated cardiomyopathy (DCM) is a structural heart disease with strong genetic background. Monogenic forms of DCM are observed in families with mutations located mostly in genes encoding structural and sarcomeric proteins. However, strong evidence suggests that genetic factors also affect the susceptibility to idiopathic DCM. To identify risk alleles for non-familial forms of DCM, we carried out a case-control association study, genotyping 664 DCM cases and 1,874 population-based healthy controls from Germany using a 50K human cardiovascular disease bead chip covering more than 2,000 genes pre-selected for cardiovascular relevance. After quality control, 30,920 single nucleotide polymorphisms (SNP) were tested for association with the disease by logistic regression adjusted for gender, and results were genomic-control corrected. The analysis revealed a significant association between a SNP in HSPB7 gene (rs1739843, minor allele frequency 39%) and idiopathic DCM (p = 1.06×10⁻⁶, OR = 0.67 [95% CI 0.57–0.79] for the minor allele T). Three more SNPs showed p < 2.21×10⁻⁵. De novo genotyping of these four SNPs was done in three independent case-control studies of idiopathic DCM. Association between SNP rs1739843 and DCM was significant in all replication samples: Germany (n = 564, n = 981 controls, p = 2.07×10⁻³, OR = 0.79 [95% CI 0.67–0.92]), France 1 (n = 433 cases, n = 395 controls, p = 3.73×10⁻³, OR = 0.74 [95% CI 0.60–0.91]), and France 2 (n = 249 cases, n = 380 controls, p = 2.26×10⁻⁴, OR = 0.63 [95% CI 0.50–0.81]). The combined analysis of all four studies including a total of n = 1,910 cases and n = 3,630 controls showed highly significant evidence for association between rs1739843 and idiopathic DCM (p = 5.28×10⁻¹³, OR = 0.72 [95% CI 0.65–0.78]). None of the other three SNPs showed significant results in the replication stage.This finding of the HSPB7 gene from a genetic search for idiopathic DCM using a large SNP panel underscores the influence of common polymorphisms on DCM susceptibility.     

Comparative genotypic and phenotypic characterisation of methicillin-resistant Staphylococcus aureus ST398 isolated from animals and humans

The high prevalence of methicillin-resistant Staphylococcus aureus (MRSA) ST398 among pigs in certain European countries and North America and its occurrence in other animal species raises a question concerning the molecular mechanisms mediating the success of this lineage. In this study a panel of S. aureus strains belonging to sequence type (ST) 5 (n = 4), ST8 (n = 5), ST15 (n = 5), ST22 (n = 8), clonal complex (CC) 30 (n = 8), CC97 (n = 8), CC130 (n = 4), CC151 (n = 4) and ST398 (n = 18) were screened by DNA microarray and PCR for the carriage of virulence and antimicrobial resistance genes. Isolates belonging to the same sequence type/clonal complex (ST/CC) were found to share similar virulence gene profiles. The ST398 lineage displayed the lowest content of virulence genes, which consisted mainly of genes detected among the majority or all of the analysed lineages. All MRSA ST398 isolates lacked accessory virulence genes that were detected in other ST/CC. In contrast to virulence genotype, the antimicrobial resistance genes profiles varied between isolates belonging to the same ST/CC and profile similarities could be observed for isolates from different lineages. MRSA ST398 isolates in particular displayed significant diversity and high content of antimicrobial resistance genes. This was comparable with certain MRSA belonging to other sequence types particularly the equine MRSA ST8. The apparent lack of significant virulence genes among MRSA ST398 strains, demonstrates that the lineage features a unique genetic background but no ST398-specific virulence markers could be identified.     

The Guinea-Bissau family of Mycobacterium tuberculosis complex revisited

The Guinea-Bissau family of strains is a unique group of the Mycobacterium tuberculosis complex that, although genotypically closely related, phenotypically demonstrates considerable heterogeneity. We have investigated 414 M. tuberculosis complex strains collected in Guinea-Bissau between 1989 and 2008 in order to further characterize the Guinea-Bissau family of strains. To determine the strain lineages present in the study sample, binary outcomes of spoligotyping were compared with spoligotypes existing in the international database SITVIT2. The major circulating M. tuberculosis clades ranked in the following order: AFRI (n = 195, 47.10%), Latin-American-Mediterranean (LAM) (n = 75, 18.12%), ill-defined T clade (n = 53, 12.8%), Haarlem (n = 37, 8.85%), East-African-Indian (EAI) (n = 25, 6.04%), Unknown (n = 12, 2.87%), Beijing (n = 7, 1.68%), X clade (n = 4, 0.96%), Manu (n = 4, 0.97%), CAS (n = 2, 0.48%). Two strains of the LAM clade isolated in 2007 belonged to the Cameroon family (SIT61). All AFRI isolates except one belonged to the Guinea-Bissau family, i.e. they have an AFRI_1 spoligotype pattern, they have a distinct RFLP pattern with low numbers of IS6110 insertions, and they lack the regions of difference RD7, RD8, RD9 and RD10, RD701 and RD702. This profile classifies the Guinea-Bissau family, irrespective of phenotypic biovar, as part of the M. africanum West African 2 lineage, or the AFRI_1 sublineage according to the spoligtyping nomenclature. Guinea-Bissau family strains display a variation of biochemical traits classically used to differentiate M. tuberculosis from M. bovis. Yet, the differential expression of these biochemical traits was not related to any genes so far investigated (narGHJI and pncA). Guinea-Bissau has the highest prevalence of M. africanum recorded in the African continent, and the Guinea-Bissau family shows a high phylogeographical specificity for Western Africa, with Guinea-Bissau being the epicenter. Trends over time however indicate that this family of strains is waning in most parts of Western Africa, including Guinea-Bissau (p = 0.048).     

Drug susceptibility in Leishmania isolates following miltefosine treatment in cases of visceral leishmaniasis and post kala-azar dermal leishmaniasis

Background

With widespread resistance to antimonials in Visceral Leishmaniasis (VL) in the Indian subcontinent, Miltefosine (MIL) has been introduced as the first line therapy. Surveillance of MIL susceptibility in natural populations of Leishmania donovani is vital to preserve it and support the VL elimination program.

Methodology and Principal Findings

We measured in vitro susceptibility towards MIL and paromomycin (PMM) in L. donovani isolated from VL and PKDL, pre- and post-treatment cases, using an amastigote-macrophage model. MIL susceptibility of post-treatment isolates from cured VL cases (n = 13, mean IC₅₀±SD = 2.43±1.44 µM), was comparable (p>0.05) whereas that from relapses (n = 3, mean IC₅₀ = 4.72±1.99 µM) was significantly higher (p = 0.04) to that of the pre-treatment group (n = 6, mean IC₅₀ = 1.86±0.75 µM). In PKDL, post-treatment isolates (n = 3, mean IC₅₀ = 16.13±2.64 µM) exhibited significantly lower susceptibility (p = 0.03) than pre-treatment isolates (n = 5, mean IC₅₀ = 8.63±0.94 µM). Overall, PKDL isolates (n = 8, mean IC₅₀ = 11.45±4.19 µM) exhibited significantly higher tolerance (p<0.0001) to MIL than VL isolates (n = 22, mean IC₅₀ = 2.58±1.58 µM). Point mutations in the miltefosine transporter (LdMT) and its beta subunit (LdRos3) genes previously reported in parasites with experimentally induced MIL resistance were not present in the clinical isolates. Further, the mRNA expression profile of these genes was comparable in the pre- and post-treatment isolates. Parasite isolates from VL and PKDL cases were uniformly susceptible to PMM with respective mean IC₅₀ = 7.05±2.24 µM and 6.18±1.51 µM.

Conclusion

The in vitro susceptibility of VL isolates remained unchanged at the end of MIL treatment; however, isolates from relapsed VL and PKDL cases had lower susceptibility than the pre-treatment isolates. PKDL isolates were more tolerant towards MIL in comparison with VL isolates. All parasite isolates were uniformly susceptible to PMM. Mutations in the LdMT and LdRos3 genes as well as changes in the expression of these genes previously correlated with experimental resistance to MIL could not be verified for the field isolates.     

Frequent Amplification of CENPF, GMNN and CDK13 Genes in Hepatocellular Carcinomas

Genomic changes frequently occur in cancer cells during tumorigenesis from normal cells. Using the Illumina Human NS-12 single-nucleotide polymorphism (SNP) chip to screen for gene copy number changes in primary hepatocellular carcinomas (HCCs), we initially detected amplification of 35 genes from four genomic regions (1q21–41, 6p21.2–24.1, 7p13 and 8q13–23). By integrated screening of these genes for both DNA copy number and gene expression in HCC and colorectal cancer, we selected CENPF (centromere protein F/mitosin), GMNN (geminin, DNA replication inhibitor), CDK13 (cyclin-dependent kinase 13), and FAM82B (family with sequence similarity 82, member B) as common cancer genes. Each gene exhibited an amplification frequency of ∼30% (range, 20–50%) in primary HCC (n = 57) and colorectal cancer (n = 12), as well as in a panel of human cancer cell lines (n = 70). Clonogenic and invasion assays of NIH3T3 cells transfected with each of the four amplified genes showed that CENPF, GMNN, and CDK13 were highly oncogenic whereas FAM82B was not. Interestingly, the oncogenic activity of these genes (excluding FAM82B) was highly correlated with gene-copy numbers in tumor samples (correlation coefficient, r>0.423), indicating that amplifications of CENPF, GMNN, and CDK13 genes are tightly linked and coincident in tumors. Furthermore, we confirmed that CDK13 gene copy number was significantly associated with clinical onset age in patients with HCC (P = 0.0037). Taken together, our results suggest that coincidently amplified CDK13, GMNN, and CENPF genes can play a role as common cancer-driver genes in human cancers.     

Peripheral Blood Gene Expression as a Novel Genomic Biomarker in Complicated Sarcoidosis

Sarcoidosis, a systemic granulomatous syndrome invariably affecting the lung, typically spontaneously remits but in ∼20% of cases progresses with severe lung dysfunction or cardiac and neurologic involvement (complicated sarcoidosis). Unfortunately, current biomarkers fail to distinguish patients with remitting (uncomplicated) sarcoidosis from other fibrotic lung disorders, and fail to identify individuals at risk for complicated sarcoidosis. We utilized genome-wide peripheral blood gene expression analysis to identify a 20-gene sarcoidosis biomarker signature distinguishing sarcoidosis (n = 39) from healthy controls (n = 35, 86% classification accuracy) and which served as a molecular signature for complicated sarcoidosis (n = 17). As aberrancies in T cell receptor (TCR) signaling, JAK-STAT (JS) signaling, and cytokine-cytokine receptor (CCR) signaling are implicated in sarcoidosis pathogenesis, a 31-gene signature comprised of T cell signaling pathway genes associated with sarcoidosis (TCR/JS/CCR) was compared to the unbiased 20-gene biomarker signature but proved inferior in prediction accuracy in distinguishing complicated from uncomplicated sarcoidosis. Additional validation strategies included significant association of single nucleotide polymorphisms (SNPs) in signature genes with sarcoidosis susceptibility and severity (unbiased signature genes - CX3CR1, FKBP1A, NOG, RBM12B, SENS3, TSHZ2; T cell/JAK-STAT pathway genes such as AKT3, CBLB, DLG1, IFNG, IL2RA, IL7R, ITK, JUN, MALT1, NFATC2, PLCG1, SPRED1). In summary, this validated peripheral blood molecular gene signature appears to be a valuable biomarker in identifying cases with sarcoidoisis and predicting risk for complicated sarcoidosis.     

MicroRNA expression aberration as potential peripheral blood biomarkers for schizophrenia

Since brain tissue is not readily accessible, a new focus in search of biomarkers for schizophrenia is blood-based expression profiling of non-protein coding genes such as microRNAs (miRNAs), which regulate gene expression by inhibiting the translation of messenger RNAs. This study aimed to identify potential miRNA signature for schizophrenia by comparing genome-wide miRNA expression profiles in patients with schizophrenia vs. healthy controls. A genome-wide miRNA expression profiling was performed using a Taqman array of 365 human miRNAs in the mononuclear leukocytes of a learning set of 30 cases and 30 controls. The discriminating performance of potential biomarkers was validated in an independent testing set of 60 cases and 30 controls. The expression levels of the miRNA signature were then evaluated for their correlation with the patients'' clinical symptoms, neurocognitive performances, and neurophysiological functions. A seven-miRNA signature (hsa-miR-34a, miR-449a, miR-564, miR-432, miR-548d, miR-572 and miR-652) was derived from a supervised classification with internal cross-validation, with an area under the curve (AUC) of receiver operating characteristics of 93%. The putative signature was then validated in the testing set, with an AUC of 85%. Among these miRNAs, miR-34a was differentially expressed between cases and controls in both the learning (P = 0.005) and the testing set (P = 0.002). These miRNAs were differentially correlated with patients'' negative symptoms, neurocognitive performance scores, and event-related potentials. The results indicated that the mononuclear leukocyte-based miRNA profiling is a feasible way to identify biomarkers for schizophrenia, and the seven-miRNA signature warrants further investigation.     

The expression of VEGF-A is down regulated in peripheral blood mononuclear cells of patients with secondary progressive multiple sclerosis

Background

Most patients with relapsing-remitting multiple sclerosis (RRMS) eventually enter a secondary progressive (SPMS) phase, characterized by increasing neurological disability. The mechanisms underlying transition to SPMS are unknown and effective treatments and biomarkers are lacking. Vascular endothelial growth factor-A (VEGF-A) is an angiogenic factor with neuroprotective effects that has been associated with neurodegenerative diseases. SPMS has a prominent neurodegenerative facet and we investigated a possible role for VEGF-A during transition from RRMS to SPMS.

Methodology/Principal Findings

VEGF-A mRNA expression in peripheral blood mononuclear (PBMC) and cerebrospinal fluid (CSF) cells from RRMS (n = 128), SPMS (n = 55) and controls (n = 116) were analyzed using real time PCR. We demonstrate reduced expression of VEGF-A mRNA in MS CSF cells compared to controls (p<0.001) irrespective of disease course and expression levels are restored by natalizumab treatment(p<0.001). VEGF-A was primarily expressed in monocytes and our CSF findings in part may be explained by effects on relative monocyte proportions. However, VEGF-A mRNA expression was also down regulated in the peripheral compartment of SPMS (p<0.001), despite unchanged monocyte counts, demonstrating a particular phenotype differentiating SPMS from RRMS and controls. A possible association of allelic variability in the VEGF-A gene to risk of MS was also studied by genotyping for six single nucleotide polymorphisms (SNPs) in MS (n = 1114) and controls (n = 1234), which, however, did not demonstrate any significant association between VEGF-A alleles and risk of MS.

Conclusions/Significance

Expression of VEGF-A in CSF cells is reduced in MS patients compared to controls irrespective of disease course. In addition, SPMS patients display reduced VEGF-A mRNA expression in PBMC, which distinguish them from RRMS and controls. This indicates a possible role for VEGF-A in the mechanisms regulating transition to SPMS. Decreased levels of PBMC VEGF-A mRNA expression should be further evaluated as a biomarker for SPMS.     

A novel survival-based tissue microarray of pancreatic cancer validates MUC1 and mesothelin as biomarkers

Background

One–fifth of patients with seemingly ‘curable’ pancreatic ductal adenocarcinoma (PDA) experience an early recurrence and death, receiving no definable benefit from a major operation. Some patients with advanced stage tumors are deemed ‘unresectable’ by conventional staging criteria (e.g. liver metastasis), yet progress slowly. Effective biomarkers that stratify PDA based on biologic behavior are needed. To help researchers sort through the maze of biomarker data, a compendium of ∼2500 published candidate biomarkers in PDA was compiled (PLoS Med, 2009. 6(4) p. e1000046).

Methods and Findings

Building on this compendium, we constructed a survival tissue microarray (termed s-TMA) comprised of short-term (cancer-specific death <12 months, n = 58) and long-term survivors (>30 months, n = 79) who underwent resection for PDA (total, n = 137). The s-TMA functions as a biological filter to identify bona fide prognostic markers associated with survival group extremes (at least 18 months separate survival groups). Based on a stringent selection process, 13 putative PDA biomarkers were identified from the public biomarker repository. Candidates were tested against the s-TMA by immunohistochemistry to identify the best markers of tumor biology. In a multivariate model, MUC1 (odds ratio, OR = 28.95, 3+ vs. negative expression, p = 0.004) and MSLN (OR = 12.47, 3+ vs. negative expression, p = 0.01) were highly predictive of early cancer-specific death. By comparison, pathologic factors (size, lymph node metastases, resection margin status, and grade) had ORs below three, and none reached statistical significance. ROC curves were used to compare the four pathologic prognostic features (ROC area = 0.70) to three univariate molecular predictors (MUC1, MSLN, MUC2) of survival group (ROC area = 0.80, p = 0.07).

Conclusions

MUC1 and MSLN were superior to pathologic features and other putative biomarkers as predicting survival group. Molecular assays comparing cancers from short and long survivors are an effective strategy to screen biomarkers and prioritize candidate cancer genes for diagnostic and therapeutic studies.     

Novel loci for adiponectin levels and their influence on type 2 diabetes and metabolic traits: a multi-ethnic meta-analysis of 45,891 individuals

Circulating levels of adiponectin, a hormone produced predominantly by adipocytes, are highly heritable and are inversely associated with type 2 diabetes mellitus (T2D) and other metabolic traits. We conducted a meta-analysis of genome-wide association studies in 39,883 individuals of European ancestry to identify genes associated with metabolic disease. We identified 8 novel loci associated with adiponectin levels and confirmed 2 previously reported loci (P = 4.5×10⁻⁸–1.2×10⁻⁴³). Using a novel method to combine data across ethnicities (N = 4,232 African Americans, N = 1,776 Asians, and N = 29,347 Europeans), we identified two additional novel loci. Expression analyses of 436 human adipocyte samples revealed that mRNA levels of 18 genes at candidate regions were associated with adiponectin concentrations after accounting for multiple testing (p<3×10⁻⁴). We next developed a multi-SNP genotypic risk score to test the association of adiponectin decreasing risk alleles on metabolic traits and diseases using consortia-level meta-analytic data. This risk score was associated with increased risk of T2D (p = 4.3×10⁻³, n = 22,044), increased triglycerides (p = 2.6×10⁻¹⁴, n = 93,440), increased waist-to-hip ratio (p = 1.8×10⁻⁵, n = 77,167), increased glucose two hours post oral glucose tolerance testing (p = 4.4×10⁻³, n = 15,234), increased fasting insulin (p = 0.015, n = 48,238), but with lower in HDL-cholesterol concentrations (p = 4.5×10⁻¹³, n = 96,748) and decreased BMI (p = 1.4×10⁻⁴, n = 121,335). These findings identify novel genetic determinants of adiponectin levels, which, taken together, influence risk of T2D and markers of insulin resistance.     

Clinical, biological and genetic analysis of anorchia in 26 boys

Background

Anorchia is defined as the absence of testes in a 46,XY individual with a male phenotype. The cause is unknown.

Methods

We evaluated the clinical and biological presentation, and family histories of 26 boys with anorchia, and sequenced their SRY, NR5A1, INSL3, MAMLD1 genes and the T222P variant for LGR8.

Results

No patient had any associated congenital anomaly. At birth, testes were palpable bilaterally or unilaterally in 13 cases and not in 7; one patient presented with bilateral testicular torsion immediately after birth. The basal plasma concentrations of anti-Müllerian hormone (AMH, n = 15), inhibin B (n = 7) and testosterone (n = 19) were very low or undetectable in all the patients evaluated, as were the increases in testosterone after human chorionic gonadotropin (hCG, n = 12). The basal plasma concentrations of follicle stimulating hormone (FSH) were increased in 20/25, as was that of luteinising hormone in 10/22 cases. Family members of 7/26 cases had histories of primary ovarian failure in the mother (n = 2), or sister 46,XX, together with fetal malformations of the only boy with microphallus and secondary foot edema (n = 1), secondary infertility in the father (n = 2), or cryptorchidism in first cousins (n = 2). The sequences of all the genes studied were normal.

Conclusion

Undetectable plasma concentrations of AMH and inhibin B and an elevated plasma FSH, together with 46,XY complement are sufficient for diagnosis of anorchia. The hCG test is unnecessary. NR5A1 and other genes implicated in gonadal development and testicle descent were not mutated, which suggests that other genes involved in these developments contribute to the phenotypes.     

Interaction analysis between HLA-DRB1 shared epitope alleles and MHC class II transactivator CIITA gene with regard to risk of rheumatoid arthritis

HLA-DRB1 shared epitope (SE) alleles are the strongest genetic determinants for autoantibody positive rheumatoid arthritis (RA). One of the key regulators in expression of HLA class II receptors is MHC class II transactivator (CIITA). A variant of the CIITA gene has been found to associate with inflammatory diseases.We wanted to explore whether the risk variant rs3087456 in the CIITA gene interacts with the HLA-DRB1 SE alleles regarding the risk of developing RA. We tested this hypothesis in a case-control study with 11767 individuals from four European Caucasian populations (6649 RA cases and 5118 controls).We found no significant additive interaction for risk alleles among Swedish Caucasians with RA (n = 3869, attributable proportion due to interaction (AP) = 0.2, 95%CI: −0.2–0.5) or when stratifying for anti-citrullinated protein antibodies (ACPA) presence (ACPA positive disease: n = 2945, AP = 0.3, 95%CI: −0.05–0.6, ACPA negative: n = 2268, AP = −0.2, 95%CI: −1.0–0.6). We further found no significant interaction between the main subgroups of SE alleles (DRB1*01, DRB1*04 or DRB1*10) and CIITA. Similar analysis of three independent RA cohorts from British, Dutch and Norwegian populations also indicated an absence of significant interaction between genetic variants in CIITA and SE alleles with regard to RA risk.Our data suggest that risk from the CIITA locus is independent of the major risk for RA from HLA-DRB1 SE alleles, given that no significant interaction between rs3087456 and SE alleles was observed. Since a biological link between products of these genes is evident, the genetic contribution from CIITA and class II antigens in the autoimmune process may involve additional unidentified factors.