Genetic analysis of chromosomal regions 2q33, 7q32 and 19q13 in multiple sclerosis susceptibility

Multiple sclerosis (MS) is an immune-mediated demyelinating disorder of the central nervous system (CNS) affecting 0.1-0.2% of Northern European descent population. MS is considered to be a multifactorial disease, both environment and genetics play a role in its pathogenesis. Despite several decades of intense research, the etiological and pathogenic mechanisms underlying MS remain still largely unknown and no curative treatment exists. The genetic architecture underlying MS is complex with multiple genes involved. The strongest and the best characterized predisposing genetic factors for MS are located, as in other immune-mediated diseases, in the major histocompatibility complex (MHC) on chromosome 6. In humans MHC is called human leukocyte antigen (HLA). Alleles of the HLA locus have been found to associate strongly with MS and remained for many years the only consistently replicable genetic associations. However, recently other genes located outside the MHC region have been proposed as strong candidates for susceptibility to MS in several studies. In this thesis a new genetic locus located on chromosome 7q32, interferon regulatory factor 5 (IRF5), was identified in the susceptibility to MS. In particular, we found that common variation of the gene was associated with the disease in three different populations, Spanish, Swedish and Finnish. We also suggested a possible functional role for one of the risk alleles with impact on the expression of the IRF5 locus. Previous studies have pointed out a possible role played by chromosome 2q33 in the susceptibility to MS and other autoimmune disorders. The work described here also investigated the involvement of this chromosomal region in MS predisposition. After the detection of genetic association with 2q33 (article-1), we extended our analysis through fine-scale single nucleotide polymorphism (SNP) mapping to define further the contribution of this genomic area to disease pathogenesis (article-4). We found a trend (p=0.04) for association to MS with an intronic SNP located in the inducible T-cell co-stimulator (ICOS) gene, an important player in the co-stimulatory pathway of the immune system. Expression analysis of ICOS revealed a novel, previously uncharacterized, alternatively spliced isoform, lacking the extracellular domain that is needed for ligand binding. The stability of the newly-identified transcript variant and its subcellular localization were analyzed. These studies indicated that the novel isoform is stable and shows different subcellular localization as compared to full-length ICOS. The novel isoform might have a regulatory function, but further studies are required to elucidate its function. Chromosome 19q13 has been previously suggested as one of the genomic areas involved in MS predisposition. In several populations, suggestive linkage signals between MS predisposition and 19q13 have been obtained. Here, we analysed the role of allelic variation in 19q13 by family based association analysis in 782 MS families collected from Finland. In this dataset, we were not able to detect any statistically significant associations, although several previously suggested markers were included to the analysis. Replication of the previous findings on the basis of linkage disequilibrium between marker allele and disease/risk allele appears notoriously difficult because of limitations such as allelic heterogeneity. Re-sequencing based approaches may be required for elucidating the role of chromosome 19q13 with MS. This thesis has resulted in the identification of a new MS susceptibility locus (IRF5) previously associated with other inflammatory or autoimmune disorders, such as SLE. IRF5 is one of the mediators of interferons biological function. In addition to providing new insight in the possible pathogenetic pathway of the disease, this finding suggests that there might be common mechanisms between different immune-mediated disorders. Furthermore the work presented here has uncovered a novel isoform of ICOS, which may play a role in regulatory mechanisms of ICOS, an important mediator of lymphocyte activation. Further work is required to uncover its functions and possible involvement of the ICOS locus in MS susceptibility.

Gene expression profile of epithelial cells and mesenchymal cells derived from limbal explant culture

Purpose: Limbal stem cell deficiency is a challenging clinical problem and the current treatment involves replenishing the depleted limbal stem cell (LSC) pool by either limbal tissue transplantation or use of cultivated limbal epithelial cells (LEC). Our experience of cultivating the LEC on denuded human amniotic membrane using a feeder cell free method, led to identification of mesenchymal cells of limbus (MC-L), which showed phenotypic resemblance to bone marrow derived mesenchymal stem cells (MSC-BM). To understand the transcriptional profile of these cells, microarray experiments were carried out.Methods: RNA was isolated from cultured LEC, MC-L and MSC-BM and microarray experiments were carried out by using Agilent chip (4x44 k). The microarray data was validated by using Realtime and semiquntitative reverse transcription polymerase chain reaction. Results: The microarray analysis revealed specific gene signature of LEC and MC-L, and also their complementary role related to cytokine and growth factor profile, thus supporting the nurturing roles of the MC-L. We have also observed similar and differential gene expression between MC-L and MSC-BM.Conclusions: This study represents the first extensive gene expression analysis of limbal explant culture derived epithelial and mesenchymal cells and as such reveals new insight into the biology, ontogeny, and in vivo function of these cells.

Gene Expression Profiling of Preovulatory Follicle in the Buffalo Cow: Effects of Increased IGF-I Concentration on Periovulatory Events

The preovulatory follicle in response to gonadotropin surge undergoes dramatic biochemical, and morphological changes orchestrated by expression changes in hundreds of genes. Employing well characterized bovine preovulatory follicle model, granulosa cells (GCs) and follicle wall were collected from the preovulatory follicle before, 1, 10 and 22 h post peak LH surge. Microarray analysis performed on GCs revealed that 450 and 111 genes were differentially expressed at 1 and 22 h post peak LH surge, respectively. For validation, qPCR and immunocytochemistry analyses were carried out for some of the differentially expressed genes. Expression analysis of many of these genes showed distinct expression patterns in GCs and the follicle wall. To study molecular functions and genetic networks, microarray data was analyzed using Ingenuity Pathway Analysis which revealed majority of the differentially expressed genes to cluster within processes like steroidogenesis, cell survival and cell differentiation. In the ovarian follicle, IGF-I is established to be an important regulator of the above mentioned molecular functions. Thus, further experiments were conducted to verify the effects of increased intrafollicular IGF-I levels on the expression of genes associated with the above mentioned processes. For this purpose, buffalo cows were administered with exogenous bGH to transiently increase circulating and intrafollicular concentrations of IGF-I. The results indicated that increased intrafollicular concentrations of IGF-I caused changes in expression of genes associated with steroidogenesis (StAR, SRF) and apoptosis (BCL-2, FKHR, PAWR). These results taken together suggest that onset of gonadotropin surge triggers activation of various biological pathways and that the effects of growth factors and peptides on gonadotropin actions could be examined during preovulatory follicle development.

Mycobacterium tuberculosis Response Regulators, DevR and NarL, Interact in Vivo and Co-regulate Gene Expression during Aerobic Nitrate Metabolism

Mycobacterium tuberculosis genes Rv0844c/Rv0845 encoding the NarL response regulator and NarS histidine kinase are hypothesized to constitute a two-component system involved in the regulation of nitrate metabolism. However, there is no experimental evidence to support this. In this study, we established M. tuberculosis NarL/NarS as a functional two-component system and identified His(241) and Asp(61) as conserved phosphorylation sites in NarS and NarL, respectively. Transcriptional profiling between M. tuberculosis H37Rv and Delta narL mutant strain during exponential growth in broth cultures with or without nitrate defined an similar to 30-gene NarL regulon that exhibited significant overlap with DevR-regulated genes, thereby implicating a role for the DevR response regulator in the regulation of nitrate metabolism. Notably, expression analysis of a subset of genes common to NarL and DevR regulons in M. tuberculosis Delta devR, Delta devS Delta dosT, and Delta narL mutant strains revealed that in response to nitrite produced during aerobic nitrate metabolism, the DevRS/DosT regulatory system plays a primary role that is augmented by NarL. Specifically, NarL itself was unable to bind to the narK2, acg, and Rv3130c promoters in phosphorylated or unphosphorylated form; however, its interaction with DevR similar to P resulted in cooperative binding, thereby enabling co-regulation of these genes. These findings support the role of physiologically derived nitrite as a metabolic signal in mycobacteria. We propose NarL-DevR binding, possibly as a heterodimer, as a novel mechanism for co-regulation of gene expression by the DevRS/DosT and NarL/NarS regulatory systems.

Regulation of Growth, Cell Shape, Cell Division, and Gene Expression by Second Messengers (p)ppGpp and Cyclic Di-GMP in Mycobacterium smegmatis

The alarmone (p)ppGpp regulates transcription, translation, replication, virulence, lipid synthesis, antibiotic sensitivity, biofilm formation, and other functions in bacteria. Signaling nucleotide cyclic di-GMP (c-di-GMP) regulates biofilm formation, motility, virulence, the cell cycle, and other functions. In Mycobacterium smegmatis, both (p) ppGpp and c-di-GMP are synthesized and degraded by bifunctional proteins Rel(Msm) and DcpA, encoded by rel(Msm) and dcpA genes, respectively. We have previously shown that the Delta rel(Msm) and Delta dcpA knockout strains are antibiotic resistant and defective in biofilm formation, show altered cell surface properties, and have reduced levels of glycopeptidolipids and polar lipids in their cell wall (K. R. Gupta, S. Kasetty, and D. Chatterji, Appl Environ Microbiol 81:2571-2578, 2015, http://dx.doi.org/10.1128/AEM.03999-14). In this work, we have explored the phenotypes that are affected by both (p) ppGpp and c-di-GMP in mycobacteria. We have shown that both (p) ppGpp and c-di-GMP are needed to maintain the proper growth rate under stress conditions such as carbon deprivation and cold shock. Scanning electron microscopy showed that low levels of these second messengers result in elongated cells, while high levels reduce the cell length and embed the cells in a biofilm-like matrix. Fluorescence microscopy revealed that the elongated Delta rel(Msm) and Delta dcpA cells are multinucleate, while transmission electron microscopy showed that the elongated cells are multiseptate. Gene expression analysis also showed that genes belonging to functional categories such as virulence, detoxification, lipid metabolism, and cell-wall-related processes were differentially expressed. Our results suggests that both (p) ppGpp and c-di-GMP affect some common phenotypes in M. smegmatis, thus raising a possibility of cross talk between these two second messengers in mycobacteria. IMPORTANCE Our work has expanded the horizon of (p) ppGpp and c-di-GMP signaling in Gram-positive bacteria. We have come across a novel observation that M. smegmatis needs (p) ppGpp and c-di-GMP for cold tolerance. We had previously shown that the Delta rel(Msm) and Delta dcpA strains are defective in biofilm formation. In this work, the overproduction of (p) ppGpp and c-di-GMP encased M. smegmatis in a biofilm-like matrix, which shows that both (p) ppGpp and c-di-GMP are needed for biofilm formation. The regulation of cell length and cell division by (p) ppGpp was known in mycobacteria, but our work shows that c-di-GMP also affects the cell size and cell division in mycobacteria. This is perhaps the first report of c-di-GMP regulating cell division in mycobacteria.

Comparative Expression of Zebrafish Lats1 and Lats2 and Their Implication in Gastrulation Movements

Large tumor suppressor (Lats) is a Ser/Thr kinase, and it presents an important function in tumor suppression. lats was originally identified in Drosophila and recently in mammals. In mammals, it contains two homologues, lats1 and lats2. In the present study, lats1 and lats2 were characterized from zebrafish (Danio rerio), which is the first report of lats in a nonmammalian vertebrate. The primary structure, genomic organization, and phylogenesis of lats from different species were studied, and the results suggest that lats1 is the direct descendant of invertebrate lats, whereas lats2 is formed by genome duplication. In zebrafish, both lats genes are maternally expressed, while they show distinctly different expression profiles during gastrulation. lats1 is almost ubiquitously expressed through development, and lats2 is more prominently expressed in the non-neural ectoderm region of zebrafish gastrula. Most intriguingly, as revealed by cell tracing and gene expression analysis, morpholino-mediated knockdown of either lats1 or lats2 led to obvious defects of cell migration in gastrulation, indicating the functional significance of lats in gastrulation movements. Developmental Dynamics 238:28502859, 2009. (C) 2009 Wiley-Liss, Inc.

Differential expression of three Paralichthys olivaceus Hsp40 genes in responses to virus infection and heat shock

Heat shock proteins (Hsps) are a family of highly conserved cellular proteins present in all organisms, mediating a range of essential housekeeping and cytoprotective functions as well-known molecular chaperons and recently as regulators of the immune response. By subtractive suppression hybridization, three Hsp40 homologues have been identified in the flounder (Paralichthys olivaceus) embryonic cells (FEC) after treatment with UV-inactivated turbot (Scophthalmus maximus L.) rhabdovirus (SMRV), termed PoHsp40A4, PoHsp40B6 and PoHsp40B11, whose encoded proteins all possess the conserved DnaJ domain, a signature motif of the Hsp40 family. Based on different protein structure and phylogenetic analysis, they can be categorized into two subfamilies, PoHsp40A4 for Type I Hsp40, PoHsp40B6 and PoHsp40B11 for Type 11 Hsp40. Further expression analysis revealed two very different types of kinetics in response either to heat shock or to virus infection, with a marked induction for PoHsp4OA4 and a weak one for both PoHsp40B6 and PoHsp40B11. A very distinct tissue distribution of mRNA was also revealed among the three genes, even between PoHsp40B6 and PoHsp40B11. This is the first report on the transcriptional induction of Hsp40 in virally stimulated fish cells, and the differential expressions might reflect their different roles in unstressed and stressed cells. (c) 2005 Elsevier Ltd. All rights reserved.

Molecular characterization and IFN signal pathway analysis of Carassius auratus CaSTAT1 identified from the cultured cells in response to virus infection

Type I interferon (IFN) exerts its pleiotropic effects mainly through the JAK-STAT signaling pathway, which is presently best described in mammals. By subtractive suppression hybridization, two fish signaling factors, JAK1 and STAT1, had been identified in the IFN-induced crucian carp Carassius auratus L. blastulae embryonic (CAB) cells after treatment with UV-inactivated grass carp hemorrhagic virus (GCHV). Further, the full-length cDNA of STAT1, termed CaSTAT1, was obtained. It contains 2926 bp and encodes a protein of 718 aa. CaSTAT1 is most similar to rat STAT1 with 59% identity overall and displays all highly conserved domains that the STAT family possesses. Like human STAT1beta, it lacks the C-terminus acting as transcriptional activation domain in mammals. By contrast, only a single transcript was detected in virus-induced CAB cells. Expression analysis showed that CaSTAT1 could be activated by stimulation of CAB cells with poly I:C, active GCHV, UV-inactivated GCHV or CAB IFN, and displayed diverse expression patterns similar to that of mammalian STATI. Additionally, the expression of an antiviral gene CaMx1 was also induced under the same conditions, and expression difference between CaSTAT1 and CaMx1 was revealed by induction of CAB IFN. These results provide molecular evidence supporting the notion that the fish IFN signaling transduction pathway is similar to that in mammals. Fish IFN exerts its multiple functions, at least antiviral action, through a JAK-STAT pathway. (C) 2004 Elsevier Ltd. All rights reserved.

Altered gene expression and DNA damage in peripheral blood cells from Friedreich's ataxia patients: cellular model of pathology.

The neurodegenerative disease Friedreich's ataxia (FRDA) is the most common autosomal-recessively inherited ataxia and is caused by a GAA triplet repeat expansion in the first intron of the frataxin gene. In this disease, transcription of frataxin, a mitochondrial protein involved in iron homeostasis, is impaired, resulting in a significant reduction in mRNA and protein levels. Global gene expression analysis was performed in peripheral blood samples from FRDA patients as compared to controls, which suggested altered expression patterns pertaining to genotoxic stress. We then confirmed the presence of genotoxic DNA damage by using a gene-specific quantitative PCR assay and discovered an increase in both mitochondrial and nuclear DNA damage in the blood of these patients (p<0.0001, respectively). Additionally, frataxin mRNA levels correlated with age of onset of disease and displayed unique sets of gene alterations involved in immune response, oxidative phosphorylation, and protein synthesis. Many of the key pathways observed by transcription profiling were downregulated, and we believe these data suggest that patients with prolonged frataxin deficiency undergo a systemic survival response to chronic genotoxic stress and consequent DNA damage detectable in blood. In conclusion, our results yield insight into the nature and progression of FRDA, as well as possible therapeutic approaches. Furthermore, the identification of potential biomarkers, including the DNA damage found in peripheral blood, may have predictive value in future clinical trials.

Diagnosis of partial body radiation exposure in mice using peripheral blood gene expression profiles.

In the event of a terrorist-mediated attack in the United States using radiological or improvised nuclear weapons, it is expected that hundreds of thousands of people could be exposed to life-threatening levels of ionizing radiation. We have recently shown that genome-wide expression analysis of the peripheral blood (PB) can generate gene expression profiles that can predict radiation exposure and distinguish the dose level of exposure following total body irradiation (TBI). However, in the event a radiation-mass casualty scenario, many victims will have heterogeneous exposure due to partial shielding and it is unknown whether PB gene expression profiles would be useful in predicting the status of partially irradiated individuals. Here, we identified gene expression profiles in the PB that were characteristic of anterior hemibody-, posterior hemibody- and single limb-irradiation at 0.5 Gy, 2 Gy and 10 Gy in C57Bl6 mice. These PB signatures predicted the radiation status of partially irradiated mice with a high level of accuracy (range 79-100%) compared to non-irradiated mice. Interestingly, PB signatures of partial body irradiation were poorly predictive of radiation status by site of injury (range 16-43%), suggesting that the PB molecular response to partial body irradiation was anatomic site specific. Importantly, PB gene signatures generated from TBI-treated mice failed completely to predict the radiation status of partially irradiated animals or non-irradiated controls. These data demonstrate that partial body irradiation, even to a single limb, generates a characteristic PB signature of radiation injury and thus may necessitate the use of multiple signatures, both partial body and total body, to accurately assess the status of an individual exposed to radiation.

Individual mRNA expression profiles reveal the effects of specific microRNAs

BACKGROUND: MicroRNAs (miRNAs) are oligoribonucleotides with an important role in regulation of gene expression at the level of translation. Despite imperfect target complementarity, they can also significantly reduce mRNA levels. The validity of miRNA target gene predictions is difficult to assess at the protein level. We sought, therefore, to determine whether a general lowering of predicted target gene mRNA expression by endogenous miRNAs was detectable within microarray gene expression profiles. RESULTS: The target gene sets predicted for each miRNA were mapped onto known gene expression data from a range of tissues. Whether considering mean absolute target gene expression, rank sum tests or 'ranked ratios', many miRNAs with significantly reduced target gene expression corresponded to those known to be expressed in the cognate tissue. Expression levels of miRNAs with reduced target mRNA levels were higher than those of miRNAs with no detectable effect on mRNA expression. Analysis of microarray data gathered after artificial perturbation of expression of a specific miRNA confirmed the predicted increase or decrease in influence of the altered miRNA upon mRNA levels. Strongest associations were observed with targets predicted by TargetScan. CONCLUSION: We have demonstrated that the effect of a miRNA on its target mRNAs' levels can be measured within a single gene expression profile. This emphasizes the extent of this mode of regulation in vivo and confirms that many of the predicted miRNA-mRNA interactions are correct. The success of this approach has revealed the vast potential for extracting information about miRNA function from gene expression profiles.

Expression of the phosphonoalanine-degradative gene cluster from Variovorax sp Pal2 is induced by growth on phosphonoalanine and phosphonopyruvate

The phosphonopyruvate hydrolase (PalA) found in Variovorax sp., Pal2, is a novel carbon-phosphorus bond cleavage enzyme, which is expressed even in the presence of high levels of phosphate, thus permitting phosphonopyruvate to be used as the sole carbon and energy source. Analysis of the regions adjacent to the palA gene revealed the presence of the five structural genes that constitute the 2-amino-3-phosphonopropionic acid (phosphonoalanine)-degradative operon. Reverse transcriptase-PCR (RT-PCR) experiments demonstrated that all five genes in the operon are transcribed as a single mRNA and that their transcription is induced by phosphonoalanine or phosphonopyruvate. Transcriptional fusions of the regulatory region of the phosphonoalanine degradative operon with the gfp gene were constructed. Expression analysis indicated that the presence of a LysR-type regulator (encoded by the palR gene) is essential for the transcription of the structural genes of the operon. Similar gene clusters were found in the sequenced genomes of six bacterial species from the Alpha-, Beta- and Gammaproteobacteria, and analysis of metagenomic libraries revealed that sequences related to palA are widely spread in the marine environment.

Association of gene expression with sequential proliferation, differentiation and tumor formation in murine skin

Differential gene expression in two established initiation and promotion skin carcinogenesis models during promotion and tumor formation was determined by microarray technology with the purpose of distinguishing the genes more associated with neoplastic transformation from those linked with proliferation and differentiation. The first model utilized dimethylbenz[a]anthracene initiation and 12-O-tetradecanoylphorbol 13-acetate (TPA) promotion in the FVB/N mouse, and the second TPA promotion of the Tg.Ac mouse, which is endogenously initiated by virtue of an activated Ha-ras transgene. Comparison of gene expression profiles across the two models identified genes whose altered expression was associated with papilloma formation rather than TPA-induced proliferation and differentiation. DMBA suppressed TPA-induced differentiation which allowed identification of those genes associated more specifically with differentiation rather than proliferation. EASE (Expression Analysis Systemic Explorer) indicated a correlation between muscle-associated genes and skin differentiation, whereas genes involved with protein biosynthesis were strongly correlated with proliferation. For verification the altered expression of selected genes were confirmed by RT-PCR; Carbonic anhydrase 2, Thioredoxin 1 and Glutathione S-transferase omega 1 associated with papilloma formation and Enolase 3, Cystatin 6 and Filaggrin associated with TPA-induced proliferation and differentiation. In situ analysis located the papillomas Glutathione S-transferase omega 1 expression to the proliferating areas of the papillomas. Thus we have identified profiles of differential gene expression associated with the tumorigenesis and promotion stages for skin carcinogenesis in the mouse.

Protein expression profiling during chick retinal maturation: a proteomics-based approach

Background: The underlying pathways that drive retinal neurogenesis and synaptogenesis are still relatively poorly understood. Protein expression analysis can provide direct insight into these complex developmental processes. The aim of this study was therefore to employ proteomic analysis to study the developing chick retina throughout embryonic (E) development commencing at day 12 through 13, 17, 19 and post-hatch (P) 1 and 33 days.

Results: 2D proteomic and mass spectrometric analysis detected an average of 1514 spots per gel with 15 spots demonstrating either modulation or constitutive expression identified via MS. Proteins identified included alpha and beta-tubulin, alpha enolase, B-creatine kinase, gamma-actin, platelet-activating factor (PAF), PREDICTED: similar to TGF-beta interacting protein 1, capping protein (actin filament muscle Z line), nucleophosmin 1 (NPM1), dimethylarginine dimethylaminohydrolase, triosphoaphate isomerase, DJ1, stathmin, fatty acid binding protein 7 (FABP7/B-FABP), beta-synuclein and enhancer of rudimentary homologue.

Conclusion: This study builds upon previous proteomic investigations of retinal development and represents the addition of a unique data set to those previously reported. Based on reported bioactivity some of the identified proteins are most likely to be important to normal retinal development in the chick. Continued analysis of the dynamic protein populations present at the early stages and throughout retinal development will increase our understanding of the molecular events underpinning retinogenesis.

New prognostic markers, determined using gene expression analyses, reveal two distinct subtypes of chronic myelomonocytic leukaemia patients

Chronic myelomonocytic leukaemia (CMML) is a heterogeneous haematopoietic disorder characterized by myeloproliferative or myelodysplastic features. At present, the pathogenesis of this malignancy is not completely understood. In this study, we sought to analyse gene expression profiles of CMML in order to characterize new molecular outcome predictors. A learning set of 32 untreated CMML patients at diagnosis was available for TaqMan low-density array gene expression analysis. From 93 selected genes related to cancer and cell cycle, we built a five-gene prognostic index after multiplicity correction. Using this index, we characterized two categories of patients with distinct overall survival (94% vs. 19% for good and poor overall survival, respectively; P = 0.007) and we successfully validated its strength on an independent cohort of 21 CMML patients with Affymetrix gene expression data. We found no specific patterns of association with traditional prognostic stratification parameters in the learning cohort. However, the poor survival group strongly correlated with high-risk treated patients and transformation to acute myeloid leukaemia. We report here a new multigene prognostic index for CMML, independent of the gene expression measurement method, which could be used as a powerful tool to predict clinical outcome and help physicians to evaluate criteria for treatments.