Search

The current study presents the analysis of stromal cells obtained from an hyperplastic left-ovary of a Holstein cow. Cultured hyperplastic stromal cells displayed a fibroblast-like morphology and ceased proliferation after the 8th passage. The non-cancerous nature of stromal cells was confirmed by in vitro cell proliferation and migration assays. Negligible amounts of E2 were detected in the spent media of cultured stromal cells, which suggests that stromal cells were non-estradiol synthesizing cells. As revealed in immunofluorescence and gene expression analysis, the hyperplastic stromal cells explicitly expressed vimentin in their cytoskeleton. Upon hematoxylin staining, a highly dense population of stromal cells was observed in the stromal tissue of the hyperplastic ovary. To explore genome-wide alterations, mRNA microarray analysis was performed using Affymetrix Bovine Gene 1.0ST Arrays compared to normal ovarian derived stromal cells. The microarray identified 1396 differentially expressed genes, of which 733 were up- and 663 down-regulated in hyperplastic stromal cells. Importantly, asporin (ASPN) and vascular cell adhesion molecule 1 (VCAM1) were among the highly up-regulated genes. Higher expression of ASPN was also confirmed by immunohistochemistry and RT-qPCR analysis. Ingenuity pathway analysis (IPA) identified about 98 significantly enriched (-log (p value ≥ 1.3) canonical pathways, importantly of which the "Sirutin Signaling Pathway" and "Mitochondrial Dysfunction" were highly activated while "Oxidative phosphorylation" was inhibited. Additionally, higher proportion of hyperplastic stromal cells in the S-phase of cell cycle, could be attributed to higher expression levels of cell proliferation genes such as CCND2 and CDK6., (© 2024. The Author(s).)

Background: Nutrition has not only an impact on the general wellbeing of an animal but can also affect reproductive processes. In cattle, feeding regimes can influence the age of puberty onset and alter gonadal development. We analyzed effects of different milk replacer (MR) feeding regimes during rearing on ovarian physiology with specific emphasis on the numbers as well as gene expression characteristics of granulosa cells (GCs) at the age of puberty onset. Two groups of calves received either 10% or 20% of bodyweight MR per day during their first 8 weeks. After weaning, both groups were fed the same mixed ration ad libitum until slaughter at 8 months., Results: Animals of the 20% feeding group had a significantly higher body weight, but the proportion of animals having a corpus luteum at the time of slaughter was not different between groups, suggesting a similar onset of puberty. Calves of the 10% group showed a constant GC count regardless of the number of follicles (r = 0.23) whereas in the 20% group increasing numbers of GCs were detected with a higher follicle count (r = 0.71). As a first effort to find a possible molecular explanation for this unexpected limitation of GC numbers in the 10% group, we comparatively analyzed GC transcriptomes in both diet groups. The mRNA microarray analysis revealed a total of 557 differentially expressed genes comparing both groups (fold change > |1.5| and p < 0.05). OAS1X, MX2 and OAS1Z were among the top downregulated genes in the 20% vs. the 10% group, whereas top upregulated genes comprised BOLA and XCL1. All of these genes are known to be regulated by interferon. Subsequent signaling pathway analysis revealed the involvement of several immune response mechanisms in accordance with a number of interferons as upstream regulators., Conclusions: The results indicate that the plane of MR feeding in early life has an impact on the number and physiology of GCs later in life. This might influence the overall reproductive life initiated by the onset of puberty in cattle. In addition, the observed alterations in GCs of calves fed less MR might be a consequence of interferon regulated immunological pathways., (© 2023. BioMed Central Ltd., part of Springer Nature.)

Background: The immune response is a crucial factor for mediating the benefit of cardiac cell therapies. Our previous research showed that cardiomyocyte transplantation alters the cardiac immune response and, when combined with short-term pharmacological CCR2 inhibition, resulted in diminished functional benefit. However, the specific role of innate immune cells, especially CCR2 macrophages on the outcome of cardiomyocyte transplantation, is unclear., Methods: We compared the cellular, molecular, and functional outcome following cardiomyocyte transplantation in wildtype and T cell- and B cell-deficient Rag2 del mice. The cardiac inflammatory response was assessed using flow cytometry. Gene expression profile was assessed using single-cell and bulk RNA sequencing. Cardiac function and morphology were determined using magnetic resonance tomography and immunohistochemistry respectively., Results: Compared to wildtype mice, Rag2 del mice show an increased innate immune response at steady state and disparate macrophage response after MI. Subsequent single-cell analyses after MI showed differences in macrophage development and a lower prevalence of CCR2 expressing macrophages. Cardiomyocyte transplantation increased NK cells and monocytes, while reducing CCR2 - MHC-II lo macrophages. Consequently, it led to increased mRNA levels of genes involved in extracellular remodelling, poor graft survival, and no functional improvement. Using machine learning-based feature selection, Mfge8 and Ccl7 were identified as the primary targets underlying these effects in the heart., Conclusions: Our results demonstrate that the improved functional outcome following cardiomyocyte transplantation is dependent on a specific CCR2 macrophage response. This work highlights the need to study the role of the immune response for cardiomyocyte cell therapy for successful clinical translation., (© 2023. BioMed Central Ltd., part of Springer Nature.)

Background: Multiple sclerosis (MS) is a chronic, inflammatory and neurodegenerative disease that leads to irreversible damage to the brain and spinal cord. The goal of so-called "immune reconstitution therapies" (IRTs) is to achieve long-term disease remission by eliminating a pathogenic immune repertoire through intense short-term immune cell depletion. B cells are major targets for effective immunotherapy in MS., Objectives: The aim of this study was to analyze the gene expression pattern of B cells before and during IRT (i.e., before B-cell depletion and after B-cell repopulation) to better understand the therapeutic effects and to identify biomarker candidates of the clinical response to therapy., Methods: B cells were obtained from blood samples of patients with relapsing-remitting MS (n = 50), patients with primary progressive MS (n = 13) as well as healthy controls (n = 28). The patients with relapsing MS received either monthly infusions of natalizumab (n = 29) or a pulsed IRT with alemtuzumab (n = 15) or cladribine (n = 6). B-cell subpopulation frequencies were determined by flow cytometry, and transcriptome profiling was performed using Clariom D arrays. Differentially expressed genes (DEGs) between the patient groups and controls were examined with regard to their functions and interactions. We also tested for differences in gene expression between patients with and without relapse following alemtuzumab administration., Results: Patients treated with alemtuzumab or cladribine showed on average a > 20% lower proportion of memory B cells as compared to before IRT. This was paralleled by profound transcriptome shifts, with > 6000 significant DEGs after adjustment for multiple comparisons. The top DEGs were found to regulate apoptosis, cell adhesion and RNA processing, and the most highly connected nodes in the network of encoded proteins were ESR2, PHB and RC3H1. Higher mRNA levels of BCL2, IL13RA1 and SLC38A11 were seen in patients with relapse despite IRT, though these differences did not pass the false discovery rate correction., Conclusions: We show that B cells circulating in the blood of patients with MS undergoing IRT present a distinct gene expression signature, and we delineated the associated biological processes and gene interactions. Moreover, we identified genes whose expression may be an indicator of relapse risk, but further studies are needed to verify their potential value as biomarkers., (© 2023. The Author(s).)

Background: Metabolic stress, as negative energy balance on one hand or obesity on the other hand can lead to increased levels of free fatty acids in the plasma and follicular fluid of animals and humans. In an earlier study, we showed that increased oleic acid (OA) concentrations affected the function of cultured bovine granulosa cells (GCs). Here, we focus on genome wide effects of increased OA concentrations., Results: Our data showed that 413 genes were affected, of which 197 were down- and 216 up-regulated. Specifically, the expression of FSH-regulated functional key genes, CCND2, LHCGR, INHA and CYP19A1 and 17-β-estradiol (E2) production were reduced by OA treatment, whereas the expression of the fatty acid transporter CD36 was increased and the morphology of the cells was changed due to lipid droplet accumulation. Bioinformatic analysis revealed that associated pathways of the putative upstream regulators "FSH" and "Cg (choriogonadotropin)" were inhibited and activated, respectively. Down-regulated genes are over-represented in GO terms "reproductive structure/system development", "ovulation cycle process", and "(positive) regulation of gonadotropin secretion", whereas up-regulated genes are involved in "circulatory system development", "vasculature development", "angiogenesis" or "extracellular matrix/structure organization"., Conclusions: From these data we conclude that besides inhibiting GC functionality, increased OA levels seemingly promote angiogenesis and tissue remodelling, thus suggestively initiating a premature fulliculo-luteal transition. In vivo this may lead to impeded folliculogenesis and ovulation, and cause sub-fertility.

Objectives: In the bovine placenta, intimate fetomaternal contact is restricted to placentomes. Within the placentomes fetal chorionic villi interdigitate with corresponding maternal caruncular crypts. The trophoblast epithelium covering the chorionic villi consists of 80% uninucleate trophoblast cells (UTCs) and 20% trophoblast giant cells (TGCs). TGCs migrate toward the endometrium and fuse with endometrial cells to form short-lived fetomaternal hybrid cells. Thereby the TGCs transport molecules of fetal origin across the placental barrier into the maternal compartment. The UTC/TGC ratio is constant during pregnancy because UTCs can differentiate into new TGCs to replace spent TGCs. However, our understanding of this differentiation process was sparse. Therefore, we collected the data to study the gene expression profiles in UTCs and TGCs and to identify differently expressed genes between the two trophoblast cell populations. Using Gene Ontology analysis, we wanted to identify biological processes and pathways that play an important role in the differentiation of UTCs into TGCs., Data Description: Bovine placentas were from days 118 to 130 of gestation. We obtained virtually pure UTCs and TGCs using a fluorescence-activated cell sorting (FACS) method. Total RNA was extracted from the UTC and TGC isolates, labeled and hybridized to Affymetrix Bovine Gene 1.0 ST Arrays.

Background: In the bovine placenta, intimate fetomaternal contacts are restricted to discrete placentomes. Here, widely branched fetal chorionic villi interdigitate with corresponding maternal caruncular crypts. The fetal trophoblast epithelium covering the chorionic villi consists of approximately 80% uninucleate trophoblast cells (UTCs) and 20% binuclear trophoblast giant cells (TGCs). The weakly invasive TGCs migrate toward the caruncle epithelium and eventually fuse with individual epithelial cells to form short-lived fetomaternal hybrid cells. In this way, molecules of fetal origin are transported across the placental barrier and released into the maternal compartment. The UTC/TGC ratio in the trophoblast remains almost constant because approximately as many new TGCs are produced from UTCs as are consumed by the fusions. The process of developing TGCs from UTCs was insufficiently understood. Therefore, we aimed to detect differentially expressed genes (DEGs) between UTCs and TGCs and identify molecular functions and biological processes regulated by DEGs., Results: We analyzed gene expression patterns in virtually pure UTC and TGC isolates using gene arrays and detected 3193 DEGs (p < 0.05; fold change values < - 1.5 or > 1.5). Of these DEGs, 1711 (53.6%) were upregulated in TGCs and 1482 (46.4%) downregulated. Gene Ontology (GO) analyses revealed that molecular functions and biological processes regulated by DEGs are related to the extracellular matrix (ECM) and its interactions with cellular receptors, cell migration and signal transduction. Furthermore, there was some evidence that O-glycan biosynthesis in TGCs may produce sialylated short-chain O-glycans (Tn antigen, core 1 O-glycans), while the synthesis of other O-glycan core structures required for the formation of complex (i.e., branched and long-chain) O-glycans appears to be decreased in TGCs., Conclusion: The differentiation of UTCs into TGCs particularly regulates genes that enable trophoblast cells to interact with their environment. Significant differences between UTCs and TGCs in ECM composition indicate reduced anchoring of TGCs in the surrounding matrix, which might contribute to their migration and their weakly invasive interaction with the maternal endometrium. Furthermore, increased expression of sialylated short chain O-glycans by TGCs could facilitate the modulation of maternal immune tolerance.

Background: Previously, we could show that L-lactate affects cultured bovine granulosa cells (GC) in a specific manner driving the cells into an early pre-ovulatory phenotype. Here we studied genome wide effects in L-lactate-treated GC to further elucidate the underlying mechanisms that are responsible for the L-lactate induced transformation. Cultured estrogen producing GC treated either with L-lactate or vehicle control were subjected to mRNA microarray analysis., Results: The analysis revealed 487 differentially expressed clusters, representing 461 annotated genes. Of these, 333 (= 318 genes) were identified as up- and 154 (= 143 genes) as down-regulated. As the top up-regulated genes we detected TXNIP, H19 and AHSG as well as our previously established marker transcripts RGS2 and PTX3. The top down-regulated genes included VNN1, SLC27A2 and GFRA1, but also MYC and the GC marker transcript CYP19A1. Pathway analysis with differentially expressed genes indicated "cAMP-mediated signaling" and "Axon guidance signaling" among the most affected pathways. Furthermore, estradiol, progesterone and Vegf were identified as potential upstream regulators. An effector network analysis by IPA provided first hints that processes of "angiogenesis" and "vascularization", but also "cell movement" appeared to be activated, whereas "organismal death" was predicted to be inhibited., Conclusions: Our data clearly show that L-lactate alters gene expression in cultured bovine GC in a broad, but obviously specific manner. Pathway analysis revealed that the mode of L-lactate action in GC initiates angiogenic processes, but also migratory events like cell movement and axonal guidance signaling, thus supporting the transformation of GC into an early luteal phenotype.

Background: Many genes important for reproductive performance are shared by both sexes. However, fecundity indices are primarily based on female parameters such as litter size. We examined a fertility mouse line (FL2), which has a considerably increased number of offspring and a total litter weight of 180% compared to a randomly bred control line (Ctrl) after more than 170 generations of breeding. In the present study, we investigated whether there might be a parallel evolution in males after more than 40 years of breeding in this outbred mouse model., Results: Males of the fertility mouse line FL2 showed reduced sperm motility performance in a 5 h thermal stress experiment and reduced birth rate in the outbred mouse line. Transcriptional analysis of the FL2 testis showed the differential expression of genes associated with steroid metabolic processes (Cyp1b1, Cyp19a1, Hsd3b6, and Cyp21a1) and female fecundity (Gdf9), accompanied by 150% elevated serum progesterone levels in the FL2 males. Cluster analysis revealed the downregulation of genes of the kallikrein-related peptidases (KLK) cluster located on chromosome 7 in addition to alterations in gene expression with serine peptidase activity, e.g., angiotensinogen (Agt), of the renin-angiotensin system essential for ovulation. Although a majority of functional annotations map to female reproduction and ovulation, these genes are differentially expressed in FL2 testis., Conclusions: These data indicate that selection for primary female traits of increased litter size not only affects sperm characteristics but also manifests as transcriptional alterations of the male side likely with direct long-term consequences for the reproductive performance of the mouse line.

Background: ADAM23 is widely expressed in the embryonic central nervous system and plays an important role in tissue formation., Results: In this study, we showed that ADAM23 contributes to cell survival and is involved in neuronal differentiation during the differentiation of human neural progenitor cells (hNPCs). Upregulation of ADAM23 in hNPCs was found to increase the number of neurons and the length of neurite, while its downregulation decreases them and triggers cell apoptosis. RNA microarray analysis revealed mechanistic insights into genes and pathways that may become involved in multiple cellular processes upon up- or downregulation of ADAM23., Conclusions: Our results suggest that ADAM23 regulates neuronal differentiation by triggering specific signaling pathways during hNPC differentiation.

Background: In previous studies it has been shown that bovine granulosa cells (GC) cultured at a high plating density dramatically change their physiological and molecular characteristics, thus resembling an early stage of luteinization. During the present study, these specific effects on the GC transcriptome were comprehensively analysed to clarify the underlying mechanisms., Methods: GC were cultured in serum free medium with FSH and IGF-1 stimulation at different initial plating density. The estradiol and progesterone production was determined by radioimmunoassays and the gene expression profiles were analysed by mRNA microarray analysis after 9 days. The data were statistically analysed and the abundance of selected, differentially expressed transcripts was re-evaluated by qPCR. Bioinformatic pathway analysis of density affected transcripts was done using Ingenuity Pathway Analysis., Results: The data showed that at high plating density the expression of 1510 annotated genes, represented by 1575 transcript clusters, showed highly altered expression levels. Nearly two-thirds were up- and one third down-regulated. Within the top up-regulated genes VNN2, RGS2 and PTX3 could be identified, as well as HBA or LOXL2. Down-regulated genes included important key genes of folliculogenesis like CYP19A1 and FSHR. Ingenuity pathway analysis identified "AMPK signaling" as well as "cAMP-mediated signaling" as major pathways affected by the alteration of the expression profile. Main putative upstream regulators were TGFB1 and VEGF, thus indicating a connection with cell differentiation and angiogenesis. A detailed cluster analysis revealed one single cluster that was highly associated with the upstream regulator beta-estradiol. Within this cluster key genes of steroid biosynthesis were not included, but instead, other genes importantly involved in follicular development, like OXT and VEGFA as well as the three most down-regulated genes TXNIP, PAG11 and ARRDC4 were identified., Conclusions: From these data we hypothesize that high density conditions induce a stage of differentiation in cultured GC that is similar to early post-LH conditions in vivo. Furthermore we hypothesize that specific cell-cell-interactions led to this differentiation including transformations necessary to promote angiogenesis.

Background: Because inheritance is recognized as playing a role in age at menarche and natural menopause, the development of chemotherapy-induced amenorrhea (CIA) might depend on inherited genetic factors; however, studies that explore such a correlation are few and have received scant attention. Given the importance of this topic we conducted a comprehensive genotype study in young women (≤45 years) with early-stage breast cancer., Methods: Our approach tested the effect of variant polymorphisms in drug metabolism enzymes (DMEs) using a predesigned pharmacogenomics panel (TaqMan® OpenArray®, Life Technologies GmbH, Darmstadt, Germany) in premenopausal women (n = 50). Patients received contemporary chemotherapy; in all cases a cyclophosphamide-based regimen with a dose of at least 500 mg/m(2) for six cycles. CIA was considered to be present in women with no resumption of menstrual bleeding within 12 months after completion of chemotherapy or goserelin., Results: Twenty-six patients (52 %) showed CIA during follow-up whereas 24 women (48 %) remained premenopausal. Of all the DMEs studied, only the SLCO1B1*5 (rs4149056) genotype was associated with the development of CIA (P = 0.017). Of the 26 patients who were homozygous for the T/T allele SLCO1B1*5, 18 (69.2 %) developed CIA compared with 8 (30.8 %) of the 22 patients who were heterozygous (C/T allele). The association of heterozygous SLCO1B1*5 allele (OR 0.038; 95%CI: 0.05-0.92) with a lower risk of developing CIA remained significant in a binary logistic regression analysis that include age, SLCO1B1*5 allele variants, and goserelin therapy., Conclusions: Patient age and SLCO1B1*5 allele variants predict the likelihood of young women with breast cancer developing CIA.

Background: Pooled samples are frequently used in experiments measuring gene expression. In this method, RNA from different individuals sharing the same experimental conditions and explanatory variables is blended and their concentrations are jointly measured. As a matter of principle, individuals are represented in equal shares in each pool. However, some degree of disproportionality may arise from the limits of technical precision. As a consequence a special kind of technical error occurs, which can be modelled by a respective variance component. Previously published theory - allowing for variable pool sizes - has been applied to four microarray gene expression data sets from different species in order to assess the practical relevance of this type of technical error in terms of significance and size of this variance component., Results: The number of transcripts with a significant variance component due to imperfect blending was found to be 4329 (23 %) in mouse data and 7093 (49 %) in honey bees, but only 6 in rats and none whatsoever in human data. These results correspond to a false discovery rate of 5 % in each data set. The number of transcripts found to be differentially expressed between treatments was always higher when the blending error variance was neglected. Simulations clearly indicated overly-optimistic (anti-conservative) test results in terms of false discovery rates whenever this source of variability was not represented in the model., Conclusions: Imperfect equality of shares when blending RNA from different individuals into joint pools of variable size is a source of technical variation with relevance for experimental design, practice at the laboratory bench and data analysis. Its potentially adverse effects, incorrect identification of differentially expressed transcripts and overly-optimistic significance tests, can be fully avoided, however, by the sound application of recently established theory and models for data analysis.

Mutations of mediator subcomplex 12 (MED12) and of high mobility group protein AT-hook 2 (HMGA2) are driver mutations in uterine leiomyomas (UL) that have not been observed to coexist in one tumor and even rarely coexist in different UL tumors of one patient. Here we describe a patient who underwent hysterectomy because of multiple leiomyomas which were studied by cytogenetics, MED12 hotspot sequencing, and copy number variation arrays. Two of the UL tumors had different HMGA2 rearrangements not detected by G-banding. Two UL tumors had deletions of the long arm of chromosome 3, in one case associated with a MED12 mutation. Both deletions lead to the loss of MED12L showing strong similarity with MED12. It remains to be determined if this gene can play a role in leiomyomagenesis independent of MED12. In summary, the patient presented exhibits an unusual coincidence of different driver mutations among her leiomyomas.

Background: Recent findings on genetic changes in uterine leiomyomas suggest these benign tumors being a heterogeneous group of diseases in terms of molecular pathogenesis with those showing karyotype alterations as well as those characterized only by cytogenetically invisible mutations of mediator subcomplex 12 (MED12). Herein, five uterine leiomyomas (UL) with an apparently normal karyotype that lacked MED12-mutations were investigated by copy number variation arrays along with their matching myometrium to search for small genomic imbalances., Results: Of five tumors one showed chromothripsis-like phenomena with numerous gains and losses of small segments mainly clustered to five chromosomal regions i.e. 2p14-2pter, 2q33.1-2q37.3, 5q31.3-5qter,11q14.1-11qter, and 18p11.21-18q2.3. Apparently, these cells had escaped detection by classical cytogenetics. Histologically, the tumor presented as a cellular leiomyoma with extended hyalinization. Of the remaining four tumors, one had a small intragenic deletion of the HMGA2 gene that was lacking in the corresponding myometrium. The other three tumors did not show relevant copy number alterations at all., Conclusions: Overall, the results suggest that leiomyomas with an apparently normal karyotype based on classical cytogenetics and lacking MED12 mutations represent a heterogeneous group of diseases. While the HMGA2 deletion detected in one of the tumors likely represents the driver mutation and, due to its size, has escaped detection by classical cytogenetics, the extended genomic imbalances detected in one of the other cases cannot be overlooked by this method suggesting an inability of the affected cells to divide in vitro. Of particular interest in that case is the occurrence of so-called "chromothripsis" or "firestorms" without involvement of the loci of common chromosomal rearrangements in UL, as e.g. 12q14 ~ 15 and 6p21. While chromothripsis was initially described as a hallmark of malignancy, the etiology and significance of this phenomenon in benign tumors still remain obscure. In uterine smooth muscle tumors, these changes per se do not indicate malignancy.

Background: Network inference of gene expression data is an important challenge in systems biology. Novel algorithms may provide more detailed gene regulatory networks (GRN) for complex, chronic inflammatory diseases such as rheumatoid arthritis (RA), in which activated synovial fibroblasts (SFBs) play a major role. Since the detailed mechanisms underlying this activation are still unclear, simultaneous investigation of multi-stimuli activation of SFBs offers the possibility to elucidate the regulatory effects of multiple mediators and to gain new insights into disease pathogenesis., Methods: A GRN was therefore inferred from RA-SFBs treated with 4 different stimuli (IL-1 β, TNF- α, TGF- β, and PDGF-D). Data from time series microarray experiments (0, 1, 2, 4, 12 h; Affymetrix HG-U133 Plus 2.0) were batch-corrected applying 'ComBat', analyzed for differentially expressed genes over time with 'Limma', and used for the inference of a robust GRN with NetGenerator V2.0, a heuristic ordinary differential equation-based method with soft integration of prior knowledge., Results: Using all genes differentially expressed over time in RA-SFBs for any stimulus, and selecting the genes belonging to the most significant gene ontology (GO) term, i.e., 'cartilage development', a dynamic, robust, moderately complex multi-stimuli GRN was generated with 24 genes and 57 edges in total, 31 of which were gene-to-gene edges. Prior literature-based knowledge derived from Pathway Studio or manual searches was reflected in the final network by 25/57 confirmed edges (44%). The model contained known network motifs crucial for dynamic cellular behavior, e.g., cross-talk among pathways, positive feed-back loops, and positive feed-forward motifs (including suppression of the transcriptional repressor OSR2 by all 4 stimuli., Conclusion: A multi-stimuli GRN highly concordant with literature data was successfully generated by network inference from the gene expression of stimulated RA-SFBs. The GRN showed high reliability, since 10 predicted edges were independently validated by literature findings post network inference. The selected GO term 'cartilage development' contained a number of differentiation markers, growth factors, and transcription factors with potential relevance for RA. Finally, the model provided new insight into the response of RA-SFBs to multiple stimuli implicated in the pathogenesis of RA, in particular to the 'novel' potent growth factor PDGF-D.

Introduction: Discrimination of rheumatoid arthritis (RA) patients from patients with other inflammatory or degenerative joint diseases or healthy individuals purely on the basis of genes differentially expressed in high-throughput data has proven very difficult. Thus, the present study sought to achieve such discrimination by employing a novel unbiased approach using rule-based classifiers., Methods: Three multi-center genome-wide transcriptomic data sets (Affymetrix HG-U133 A/B) from a total of 79 individuals, including 20 healthy controls (control group - CG), as well as 26 osteoarthritis (OA) and 33 RA patients, were used to infer rule-based classifiers to discriminate the disease groups. The rules were ranked with respect to Kiendl's statistical relevance index, and the resulting rule set was optimized by pruning. The rule sets were inferred separately from data of one of three centers and applied to the two remaining centers for validation. All rules from the optimized rule sets of all centers were used to analyze their biological relevance applying the software Pathway Studio., Results: The optimized rule sets for the three centers contained a total of 29, 20, and 8 rules (including 10, 8, and 4 rules for 'RA'), respectively. The mean sensitivity for the prediction of RA based on six center-to-center tests was 96% (range 90% to 100%), that for OA 86% (range 40% to 100%). The mean specificity for RA prediction was 94% (range 80% to 100%), that for OA 96% (range 83.3% to 100%). The average overall accuracy of the three different rule-based classifiers was 91% (range 80% to 100%). Unbiased analyses by Pathway Studio of the gene sets obtained by discrimination of RA from OA and CG with rule-based classifiers resulted in the identification of the pathogenetically and/or therapeutically relevant interferon-gamma and GM-CSF pathways., Conclusion: First-time application of rule-based classifiers for the discrimination of RA resulted in high performance, with means for all assessment parameters close to or higher than 90%. In addition, this unbiased, new approach resulted in the identification not only of pathways known to be critical to RA, but also of novel molecules such as serine/threonine kinase 10.

Background: Epitheloid leiomyoma is a rare subtype of benign smooth muscle tumors., Results: Herein, we present the results of classical cytogenetics, MED12 mutation analysis, and copy number variation array evaluation in one such case. Whereas cytogenetic did not show evidence for clonal chromosome abnormalities and no MED12 mutation in the "fibroid hot spot" region was detected, array hybridization revealed multiple abnormalities. Most noteworthy, almost all chromosomes showed copy-number neutral loss of heterozygosity. As examples of further abnormalities, trisomies of chromosomes 8, 12, 20, and X were noted., Discussion: The data presented suggest a near-haploid karyotype of the tumor as the initial genetic alteration followed by secondary duplications of large parts of the genome. The absence of any clonal karyotypic alterations after performing classical cytogenetics is likely explained by a reduced ability of the tumor cells to proliferate in vitro. However, to the best of our knowledge this is the first report of an uterine leiomyoma showing extended uniparental disomy. It remains to be determined if this is a more common phenomenon in epithelioid leiomyomas or even subsets of "ordinary" leiomyomas.

Background: Glatiramer acetate (GA) is a mixture of synthetic peptides used in the treatment of patients with relapsing-remitting multiple sclerosis (RRMS). The aim of this study was to investigate the effects of GA therapy on the gene expression of monocytes., Methods: Monocytes were isolated from the peripheral blood of eight RRMS patients. The blood was obtained longitudinally before the start of GA therapy as well as after one day, one week, one month and two months. Gene expression was measured at the mRNA level by microarrays., Results: More than 400 genes were identified as up-regulated or down-regulated in the course of therapy, and we analyzed their biological functions and regulatory interactions. Many of those genes are known to regulate lymphocyte activation and proliferation, but only a subset of genes was repeatedly differentially expressed at different time points during treatment., Conclusions: Overall, the observed gene regulatory effects of GA on monocytes were modest and not stable over time. However, our study revealed several genes that are worthy of investigation in future studies on the molecular mechanisms of GA therapy.

Background: A subset of patients with multiple sclerosis (MS) shows an increased endogenous IFN-like activity before initiation of IFN-beta treatment. The molecular basis of this phenomenon and its relevance to predict individual therapy outcomes are not yet fully understood. We studied the expression patterns of these patients, the prognostic value of an elevated IFN-like activity, and the gene regulatory effects of exogenously administered IFN-beta., Methods: Microarray gene expression profiling was performed for 61 MS patients using peripheral blood mononuclear cells obtained before and after 1 month of IFN-beta therapy. Expression levels of genes involved in pathways either inducing or being activated by IFN-beta were compared between patients with high (MX1(high) cohort) and low (MX1(low) cohort) endogenous IFN-like activity. Patients were followed for 5 years and relapses as well as progression on the expanded disability status scale (EDSS) were documented., Results: Before the start of therapy, 11 patients presented elevated mRNA levels of IFN-stimulated genes indicative of a relatively high endogenous IFN-like activity (MX1(high)). In these patients, pathogen receptors (for example, TLR7, RIG-I and IFIH1) and transcription factors were also expressed more strongly, which could be attributed to an overactivity of IFN-stimulated gene factor 3 (ISGF3, a complex formed by STAT1, STAT2 and IFN regulatory factor 9). After 1 month of IFN-beta therapy, the expression of many pathway genes was significantly induced in MX1(low) patients, but remained unaltered in MX1(high) patients. During follow-up, relapse rate and changes in EDSS were comparable between both patient groups, with differences seen between different types of IFN-beta drug application., Conclusions: Therapeutic IFN-beta induces the transcription of several genes involved in IFN-related pathways. In a subgroup of MS patients, the expression of these genes is already increased before therapy initiation, possibly driven by an overexpression of ISGF3. Patients with high and low endogenous IFN-like activity showed similar clinical long-term courses of disease. Different results were obtained for different IFN-beta drug preparations, and this merits further investigation.

Background: Batch effects due to sample preparation or array variation (type, charge, and/or platform) may influence the results of microarray experiments and thus mask and/or confound true biological differences. Of the published approaches for batch correction, the algorithm "Combating Batch Effects When Combining Batches of Gene Expression Microarray Data" (ComBat) appears to be most suitable for small sample sizes and multiple batches., Methods: Synovial fibroblasts (SFB; purity > 98%) were obtained from rheumatoid arthritis (RA) and osteoarthritis (OA) patients (n = 6 each) and stimulated with TNF-α or TGF-β1 for 0, 1, 2, 4, or 12 hours. Gene expression was analyzed using Affymetrix Human Genome U133 Plus 2.0 chips, an alternative chip definition file, and normalization by Robust Multi-Array Analysis (RMA). Data were batch-corrected for different acquiry dates using ComBat and the efficacy of the correction was validated using hierarchical clustering., Results: In contrast to the hierarchical clustering dendrogram before batch correction, in which RA and OA patients clustered randomly, batch correction led to a clear separation of RA and OA. Strikingly, this applied not only to the 0 hour time point (i.e., before stimulation with TNF-α/TGF-β1), but also to all time points following stimulation except for the late 12 hour time point. Batch-corrected data then allowed the identification of differentially expressed genes discriminating between RA and OA. Batch correction only marginally modified the original data, as demonstrated by preservation of the main Gene Ontology (GO) categories of interest, and by minimally changed mean expression levels (maximal change 4.087%) or variances for all genes of interest. Eight genes from the GO category "extracellular matrix structural constituent" (5 different collagens, biglycan, and tubulointerstitial nephritis antigen-like 1) were differentially expressed between RA and OA (RA > OA), both constitutively at time point 0, and at all time points following stimulation with either TNF-α or TGF-β1., Conclusion: Batch correction appears to be an extremely valuable tool to eliminate non-biological batch effects, and allows the identification of genes discriminating between different joint diseases. RA-SFB show an upregulated expression of extracellular matrix components, both constitutively following isolation from the synovial membrane and upon stimulation with disease-relevant cytokines or growth factors, suggesting an "imprinted" alteration of their phenotype.

Background: Udder infections with environmental pathogens like Escherichia coli are a serious problem for the dairy industry. Reduction of incidence and severity of mastitis is desirable and mild priming of the immune system either through vaccination or with low doses of immune stimulants such as lipopolysaccharide LPS was previously found to dampen detrimental effects of a subsequent infection. Monocytes/macrophages are known to develop tolerance towards the endotoxin LPS (endotoxin tolerance, ET) as adaptation strategy to prevent exuberant inflammation.We have recently observed that infusion of 1 μg of LPS into the quarter of an udder effectively protected for several days against an experimentally elicited mastitis. We have modelled this process in primary cultures of mammary epithelial cells (MEC) from the cow. MEC are by far the most abundant cells in the healthy udder coming into contact with invading pathogens and little is known about their role in establishing ET., Results: We primed primary MEC cultures for 12 h with LPS (100 ng/ml) and stimulated three cultures either 12 h or 42 h later with 107/ml particles of heat inactivated E. coli bacteria for six hours. Priming-related alterations in the global transcriptome of those cells were quantified with Affymetrix microarrays. LPS priming alone caused differential expression of 40 genes and mediated significantly different response to a subsequent E. coli challenge of 226 genes. Expression of 38 genes was enhanced while that of 188 was decreased. Higher expressed were anti-microbial factors (β-defensin LAP, SLPI), cell and tissue protecting factors (DAF, MUC1, TGM1, TGM3) as well as mediators of the sentinel function of MEC (CCL5, CXCL8). Dampened was the expression of potentially harmful pro-inflammatory master cytokines (IL1B, IL6, TNF-α) and immune effectors (NOS2, matrix metalloproteases). Functional network analysis highlighted the reduced expression of IL1B and of IRF7 as key to this modulation., Conclusion: LPS-primed MEC are fitter to repel pathogens and better protected against misguided attacks of the immune response. Attenuated is the exuberant expression of factors potentially promoting immunopathological processes. MEC therefore recapitulate many aspects of ET known so far from professional immune cells.

ReNcell VM is an immortalized human neural progenitor cell line with the ability to differentiate in vitro into astrocytes and neurons, in which the Wnt/β-catenin pathway is known to be involved. However, little is known about kinetic changes of this pathway in human neural progenitor cell differentiation. In the present study, we provide a quantitative profile of Wnt/β-catenin pathway dynamics showing its spatio-temporal regulation during ReNcell VM cell differentiation. We show first that T-cell factor dependent transcription can be activated by stabilized β-catenin. Furthermore, endogenous Wnt ligands, pathway receptors and signaling molecules are temporally controlled, demonstrating changes related to differentiation stages. During the first three hours of differentiation the signaling molecules LRP6, Dvl2 and β-catenin are spatio-temporally regulated between distinct cellular compartments. From 24 h onward, components of the Wnt/β-catenin pathway are strongly activated and regulated as shown by mRNA up-regulation of Wnt ligands (Wnt5a and Wnt7a), receptors including Frizzled-2, -3, -6, -7, and -9, and co-receptors, and target genes including Axin2. This detailed temporal profile of the Wnt/β-catenin pathway is a first step to understand, control and to orientate, in vitro, human neural progenitor cell differentiation.

Background: Expansion of multi-C2H2 domain zinc finger (ZNF) genes, including the Krüppel-associated box (KRAB) subfamily, paralleled the evolution of tetrapodes, particularly in mammalian lineages. Advances in their cataloging and characterization suggest that the functions of the KRAB-ZNF gene family contributed to mammalian speciation., Results: Here, we characterized the human 8q24.3 ZNF cluster on the genomic, the phylogenetic, the structural and the transcriptome level. Six (ZNF7, ZNF34, ZNF250, ZNF251, ZNF252, ZNF517) of the seven locus members contain exons encoding KRAB domains, one (ZNF16) does not. They form a paralog group in which the encoded KRAB and ZNF protein domains generally share more similarities with each other than with other members of the human ZNF superfamily. The closest relatives with respect to their DNA-binding domain were ZNF7 and ZNF251. The analysis of orthologs in therian mammalian species revealed strong conservation and purifying selection of the KRAB-A and zinc finger domains. These findings underscore structural/functional constraints during evolution. Gene losses in the murine lineage (ZNF16, ZNF34, ZNF252, ZNF517) and potential protein truncations in primates (ZNF252) illustrate ongoing speciation processes. Tissue expression profiling by quantitative real-time PCR showed similar but distinct patterns for all tested ZNF genes with the most prominent expression in fetal brain. Based on accompanying expression signatures in twenty-six other human tissues ZNF34 and ZNF250 revealed the closest expression profiles. Together, the 8q24.3 ZNF genes can be assigned to a cerebellum, a testis or a prostate/thyroid subgroup. These results are consistent with potential functions of the ZNF genes in morphogenesis and differentiation. Promoter regions of the seven 8q24.3 ZNF genes display common characteristics like missing TATA-box, CpG island-association and transcription factor binding site (TFBS) modules. Common TFBS modules partly explain the observed expression pattern similarities., Conclusions: The ZNF genes at human 8q24.3 form a relatively old mammalian paralog group conserved in eutherian mammals for at least 130 million years. The members persisted after initial duplications by undergoing subfunctionalizations in their expression patterns and target site recognition. KRAB-ZNF mediated repression of transcription might have shaped organogenesis in mammalian ontogeny.

Background: Due to the rapid data accumulation on pathogenesis and progression of chronic inflammation, there is an increasing demand for approaches to analyse the underlying regulatory networks. For example, rheumatoid arthritis (RA) is a chronic inflammatory disease, characterised by joint destruction and perpetuated by activated synovial fibroblasts (SFB). These abnormally express and/or secrete pro-inflammatory cytokines, collagens causing joint fibrosis, or tissue-degrading enzymes resulting in destruction of the extra-cellular matrix (ECM).We applied three methods to analyse ECM regulation: data discretisation to filter out noise and to reduce complexity, Boolean network construction to implement logic relationships, and formal concept analysis (FCA) for the formation of minimal, but complete rule sets from the data., Results: First, we extracted literature information to develop an interaction network containing 18 genes representing ECM formation and destruction. Subsequently, we constructed an asynchronous Boolean network with biologically plausible time intervals for mRNA and protein production, secretion, and inactivation. Experimental gene expression data was obtained from SFB stimulated by TGFbeta1 or by TNFalpha and discretised thereafter. The Boolean functions of the initial network were improved iteratively by the comparison of the simulation runs to the experimental data and by exploitation of expert knowledge. This resulted in adapted networks for both cytokine stimulation conditions. The simulations were further analysed by the attribute exploration algorithm of FCA, integrating the observed time series in a fine-tuned and automated manner. The resulting temporal rules yielded new contributions to controversially discussed aspects of fibroblast biology (e.g., considerable expression of TNF and MMP9 by fibroblasts stimulation) and corroborated previously known facts (e.g., co-expression of collagens and MMPs after TNFalpha stimulation), but also revealed some discrepancies to literature knowledge (e.g., MMP1 expression in the absence of FOS)., Conclusion: The newly developed method successfully and iteratively integrated expert knowledge at different steps, resulting in a promising solution for the in-depth understanding of regulatory pathways in disease dynamics. The knowledge base containing all the temporal rules may be queried to predict the functional consequences of observed or hypothetical gene expression disturbances. Furthermore, new hypotheses about gene relations were derived which await further experimental validation.

We examined the repertoire and extent of inflammation dependent gene regulation in a bovine mammary epithelial cell (MEC) model, to better understand the contribution of the MEC in the immune defence of the udder. We challenged primary cultures of MEC from cows with heat inactivated Escherichia coli pathogens and used Affymetrix DNA-microarrays to profile challenge related alterations in their transcriptome. Compared to acute mastitis, the most prominently activated genes comprise those encoding chemokines, interleukins, beta-defensins, serum amyloid A and haptoglobin. Hence, the MEC exert sentinel as well as effector functions of innate immune defence. E. coli stimulated a larger fraction of genes (30%) in the MEC belonging to the functional category Inflammatory Response than we recorded with the same microarrays during acute mastitis in the udder (17%). This observation underscores the exquisite immune capacity of MEC. To more closely examine the adequacy of immunological regulation in MEC, we compared the inflammation dependent regulation of factors contributing to the complement system between the udder versus the MEC. In the MEC we observed only up regulation of several complement factor-encoding genes. Mastitis, in contrast, in the udder strongly down regulates such genes encoding factors contributing to both, the classical pathway of complement activation and the Membrane Attack Complex, while the expression of factors contributing to the alternative pathway may be enhanced. This functionally polarized regulation of the complex complement pathway is not reflected in the MEC models.

Background: In a recent study it was demonstrated that a largely increased ovulation number is responsible for high prolificacy in two mouse lines selected for fertility performance. The objective of the present study was to identify genes that are involved in increasing the ovulation number in one of these lines, FL1. For differential expression profiling, ovaries of FL1 and of a non-selected control line, DUKsi, both lines derived from the same genetic pool, were analyzed with microarray analysis and quantitative polymerase chain reaction (qPCR). Ovaries from 30 animals of each line were collected at the metestrous stage, combined to 6 pools each, and processed for microarray analysis., Results: The actual number of ova shed in FL1 exceeded that of the DUKsi control line more than twofold (26.6 vs. 12.9). 148 differentially expressed ovarian transcripts could be identified, 74 of them up- and 74 down-regulated. Of these, 47 significantly mapped to specific Gene Ontology (GO) terms representing different biological processes as steroid metabolism, folliculogenesis, immune response, intracellular signal transduction (particularly of the G protein signaling cascade), regulation of transcription and translation, cell cycle and others. qPCR was used to re-evaluate selected transcripts and to estimate inter-individual variation of expression levels. These data significantly correlated with microarray data in 12 out of 15 selected transcripts but revealed partly large variations of expression levels between individuals., Conclusion: (1) The abundance of numerous ovarian transcripts was significantly different in FL1 compared to the non-selected control line DUKsi thus suggesting that at least some of the respective genes and corresponding biological processes are involved in improving reproductive traits, particularly by increasing the number of ovulation. (2) Selective qPCR re-evaluation largely confirmed the microarray data and in addition demonstrated that sample pooling can be beneficial to find out group-specific expression profiles despite of large inter-individual variation. (3) The present data will substantially help ongoing genetic association studies to identify candidate genes and causative mutations responsible for increased fertility performance in mice.

Introduction: About 30% of rheumatoid arthritis patients fail to respond adequately to TNFalpha-blocking therapy. There is a medical and socioeconomic need to identify molecular markers for an early prediction of responders and nonresponders., Methods: RNA was extracted from peripheral blood mononuclear cells of 19 rheumatoid arthritis patients before the first application of the TNFalpha blocker etanercept as well as after 72 hours. Clinical response was assessed over 3 months using the 28-joint-count Disease Activity Score and X-ray scans. Supervised learning methods were applied to Affymetrix Human Genome U133 microarray data analysis to determine highly selective discriminatory gene pairs or triplets with prognostic relevance for the clinical outcome evinced by a decline of the 28-joint-count Disease Activity Score by 1.2., Results: Early downregulation of expression levels secondary to TNFalpha neutralization was associated with good clinical responses, as shown by a decline in overall disease activity 3 months after the start of treatment. Informative gene sets include genes (for example, NFKBIA, CCL4, IL8, IL1B, TNFAIP3, PDE4B, PPP1R15A and ADM) involved in different pathways and cellular processes such as TNFalpha signalling via NFkappaB, NFkappaB-independent signalling via cAMP, and the regulation of cellular and oxidative stress response. Pairs and triplets within these genes were found to have a high prognostic value, reflected by prediction accuracies of over 89% for seven selected gene pairs and of 95% for 10 specific gene triplets., Conclusion: Our data underline that early gene expression profiling is instrumental in identifying candidate biomarkers to predict therapeutic outcomes of anti-TNFalpha treatment regimes.

Introduction: Rheumatoid arthritis (RA) is a chronic inflammatory and destructive joint disease characterized by overexpression of pro-inflammatory/pro-destructive genes and other activating genes (for example, proto-oncogenes) in the synovial membrane (SM). The gene expression in disease is often characterized by significant inter-individual variances via specific synchronization/desynchronization of gene expression. To elucidate the contribution of the variance to the pathogenesis of disease, expression variances were tested in SM samples of RA patients, osteoarthritis (OA) patients, and normal controls (NCs)., Method: Analysis of gene expression in RA, OA, and NC samples was carried out using Affymetrix U133A/B oligonucleotide arrays, and the results were validated by real-time reverse transcription-polymerase chain reaction. For the comparison between RA and NC, 568 genes with significantly different variances in the two groups (P

Genome-wide gene expression was comparatively investigated in early-passage rheumatoid arthritis (RA) and osteoarthritis (OA) synovial fibroblasts (SFBs; n = 6 each) using oligonucleotide microarrays; mRNA/protein data were validated by quantitative PCR (qPCR) and western blotting and immunohistochemistry, respectively. Gene set enrichment analysis (GSEA) of the microarray data suggested constitutive upregulation of components of the transforming growth factor (TGF)-beta pathway in RA SFBs, with 2 hits in the top 30 regulated pathways. The growth factor TGF-beta1, its receptor TGFBR1, the TGF-beta binding proteins LTBP1/2, the TGF-beta-releasing thrombospondin 1 (THBS1), the negative effector SkiL, and the smad-associated molecule SARA were upregulated in RA SFBs compared to OA SFBs, whereas TGF-beta2 was downregulated. Upregulation of TGF-beta1 and THBS1 mRNA (both positively correlated with clinical markers of disease activity/severity) and downregulation of TGF-beta2 mRNA in RA SFBs were confirmed by qPCR. TGFBR1 mRNA (only numerically upregulated in RA SFBs) and SkiL mRNA were not differentially expressed. At the protein level, TGF-beta1 showed a slightly higher expression, and the signal-transducing TGFBR1 and the TGF-beta-activating THBS1 a significantly higher expression in RA SFBs than in OA SFBs. Consistent with the upregulated TGF-beta pathway in RA SFBs, stimulation with TGF-beta1 resulted in a significantly enhanced expression of matrix-metalloproteinase (MMP)-11 mRNA and protein in RA SFBs, but not in OA SFBs. In conclusion, RA SFBs show broad, constitutive alterations of the TGF-beta pathway. The abundance of TGF-beta, in conjunction with an augmented mRNA and/or protein expression of TGF-beta-releasing THBS1 and TGFBR1, suggests a pathogenetic role of TGF-beta-induced effects on SFBs in RA, for example, the augmentation of MMP-mediated matrix degradation/remodeling.

Quantitative traits such as complex diseases are controlled by many small-effect genes that are difficult to identify. Here we present a novel strategy to identify the candidate genes for small-effect quantitative trait loci (QTL) in collagen induced arthritis (CIA) using global genome and transcriptome approaches. First, we performed genome linkage analysis in F2 progeny of the CIA susceptible and resistant strains to search for small-effect QTL. Second, we detected gene expression patterns of both strains during CIA. The candidate genes were identified using three criteria: they are located in a genomic region linked to CIA; they are disease-specific differentially expressed during CIA; and they are strain-specific differentially expressed regarding the two parental strains. Eight small-effect QTL controlling CIA severity were identified. Of 22,000 screened genes, 117 were both strain-specific and disease-specific differentially expressed during CIA. Of these 117 genes, 21 were located inside the support intervals of the 8 small-effect QTL and thus were considered as candidate genes.

Background: Biological tissues consist of various cell types that differentially contribute to physiological and pathophysiological processes. Determining and analyzing cell type-specific gene expression under diverse conditions is therefore a central aim of biomedical research. The present study compares gene expression profiles in whole tissues and isolated cell fractions purified from these tissues in patients with rheumatoid arthritis and osteoarthritis., Results: The expression profiles of the whole tissues were compared to computationally reconstituted expression profiles that combine the expression profiles of the isolated cell fractions (macrophages, fibroblasts, and non-adherent cells) according to their relative mRNA proportions in the tissue. The mRNA proportions were determined by trimmed robust regression using only the most robustly-expressed genes (1/3 to 1/2 of all measured genes), i.e. those showing the most similar expression in tissue and isolated cell fractions. The relative mRNA proportions were determined using several different chip evaluation methods, among which the MAS 5.0 signal algorithm appeared to be most robust. The computed mRNA proportions agreed well with the cell proportions determined by immunohistochemistry except for a minor number of outliers. Genes that were either regulated (i.e. differentially-expressed in tissue and isolated cell fractions) or robustly-expressed in all patients were identified using different test statistics., Conclusion: Robust Computational Reconstitution uses an intermediate number of robustly-expressed genes to estimate the relative mRNA proportions. This avoids both the exclusive dependence on the robust expression of individual, highly cell type-specific marker genes and the bias towards an equal distribution upon inclusion of all genes for computation.

A better understanding of the initial mechanisms that lead to arthritic disease could facilitate development of improved therapeutic strategies. We characterized the synovial microcirculation of knee joints in susceptible mouse strains undergoing intradermal immunization with bovine collagen II in complete Freund's adjuvant to induce arthritis (i.e. collagen-induced arthritis [CIA]). Susceptible DBA1/J and collagen II T-cell receptor transgenic mice were compared with CIA-resistant FVB/NJ mice. Before onset of clinical symptoms of arthritis, in vivo fluorescence microscopy of knee joints revealed marked leucocyte activation and interaction with the endothelial lining of synovial microvessels. This initial inflammatory cell response correlated with the gene expression profile at this disease stage. The majority of the 655 differentially expressed genes belonged to classes of genes that are involved in cell movement and structure, cell cycle and signal transduction, as well as transcription, protein synthesis and metabolism. However, 24 adhesion molecules and chemokine/cytokine genes were identified, some of which are known to contribute to arthritis (e.g. CD44 and neutrophil cytosolic factor 1) and some of which are novel in this respect (e.g. CC chemokine ligand-27 and IL-13 receptor alpha1). Online in vivo data on synovial tissue microcirculation, together with gene expression profiling, emphasize the potential role played by early inflammatory events in the development of arthritis.

Collagen-induced arthritis (CIA), an approved animal model for rheumatoid arthritis, is thought to be a T cell-dependent disease. There is evidence that CD8+ T cells are a major subset controlling the pathogenesis of CIA. They probably contribute to certain features of disease, namely tissue destruction and synovial hyperplasia. In this study we examined the role of perforin (pfp), a key molecule of the cytotoxic death pathway that is expressed mainly in CD8+ T cells, for the pathogenesis of CIA. We generated DBA/1J mice suffering from mutations of the pfp molecule, DBA/1J-pfp-/-, and studied their susceptibility to arthritis. As a result, pfp-deficient mice showed a reduced incidence (DBA/1J-pfp+/+, 64%; DBA/1J-pfp-/-, 54%), a slightly delayed onset (onset of disease: DBA/1J-pfp+/+, 53 +/- 3.6; DBA/1J-pfp-/-, 59 +/- 4.9 (mean +/- SEM), and milder form of the disease (maximum disease score: DBA/1J-pfp+/+, 7.3 +/- 1.1; DBA/1J-pfp-/-, 3.4 +/- 1.4 (mean +/- SEM); P < 0.05). Concomitantly, peripheral T cell proliferation in response to the specific antigen bovine collagen II was increased in pfp-/- mice compared with pfp+/+ mice, arguing for an impaired killing of autoreactive T cells caused by pfp deficiency. Thus, pfp-mediated cytotoxicity is involved in the initiation of tissue damage in arthritis, but pfp-independent cytotoxic death pathways might also contribute to CIA.

Arthritis susceptibility genes were sought by analysis of differential gene expression between pristane-induced arthritis (PIA)-susceptible DA rats and PIA-resistant E3 rats. Inguinal lymph nodes of naïve animals and animals 8 days after pristane injection were analyzed for differential gene expression. mRNA expression was investigated by microarray and real-time PCR, and protein expression was analyzed by flow cytometry or ELISA. Twelve genes were significantly differentially expressed when analyzed by at least two independent methods, and an additional five genes showed a strong a tendency toward differential expression. In naïve DA rats IgE, the bone marrow stromal cell antigen 1 (Bst1) and the MHC class II beta-chain (MhcII) were expressed at a higher level, and the immunoglobulin kappa chain (Igkappa) was expressed at a lower level. In pristane-treated DA rats the MHC class II beta-chain, gelatinase B (Mmp9) and the protein tyrosine phosphatase CL100 (Ptpn16) were expressed at a higher level, whereas immunoglobulins, the CD28 molecule (Cd28), the mast cell specific protease 1 (Mcpt1), the carboxylesterase precursor (Ces2), K-cadherin (Cdh6), cyclin G1 (Ccng1), DNA polymerase IV (Primase) and the tumour associated glycoprotein E4 (Tage) were expressed at a lower level. Finally, the differentially expressed mRNA was confirmed with protein expression for some of the genes. In conclusion, the results show that animal models are well suited for reproducible microarray analysis of candidate genes for arthritis. All genes have functions that are potentially important for arthritis, and nine of the genes are located within genomic regions previously associated with autoimmune disease.

To improve our knowledge on the pathophysiology of rheumatoid arthritis (RA), we investigated gene expression patterns in synovial tissue from RA and osteoarthritis (OA) patients. DNA oligonucleotide microarray analysis was employed to identify differentially expressed genes in synovial tissue from pathologically classified tissue samples from RA (n = 20) and OA patients (n = 10). From 7131 gene sets displayed on the microarray chip, 101 genes were found to be upregulated and 300 genes to be downregulated in RA as compared with OA. Semiquantitative reverse-transcription polymerase chain reaction, Western blotting and immunohistochemistry were used to validate microarray expression levels. These experiments revealed that Cys-X-Cys receptor (CXCR)1, CXCR2 and CXCR3 mRNAs, as well as Cys-X-Cys ligand (CXCL)9 (monokine induced by IFN-gamma) and CXCL10 (IFN-gamma inducible protein 10) mRNAs, were significantly upregulated in RA as compared with OA disease. Elevated protein levels in RA synovial tissue were detected for CXCR1 and CXCR3 by Western blotting. Using immunohistochemistry, CXCR3 protein was found to be preferentially expressed on mast cells within synovial tissue from RA patients. These findings suggest that substantial expression of CXCR3 protein on mast cells within synovial tissue from RA patients plays a significant role in the pathophysiology of RA, accompanied by elevated levels of the chemokines CXCL9 and CXCL10. Mature mast cells are likely to contribute to and sustain the inflamed state in arthritic lesions (e.g. by production of inflammatory mediators such as histamine, proteinases, arachidonic acid metabolites and cytokines). Thus, the mast cell could become a potential target in therapeutic intervention.