Paired guide RNA CRISPR-Cas9 screening for protein-coding genes and lncRNAs involved in transdifferentiation of human B-cells to macrophages.

CRISPR-Cas9 screening libraries have arisen as a powerful tool to identify protein-coding (pc) and non-coding genes playing a role along different processes. In particular, the usage of a nuclease active Cas9 coupled to a single gRNA has proven to efficiently impair the expression of pc-genes by generating deleterious frameshifts. Here, we first demonstrate that targeting the same gene simultaneously with two guide RNAs (paired guide RNAs, pgRNAs) synergistically enhances the capacity of the CRISPR-Cas9 system to knock out pc-genes. We next design a library to target, in parallel, pc-genes and lncRNAs known to change expression during the transdifferentiation from pre-B cells to macrophages. We show that this system is able to identify known players in this process, and also predicts 26 potential novel ones, of which we select four (two pc-genes and two lncRNAs) for deeper characterization. Our results suggest that in the case of the candidate lncRNAs, their impact in transdifferentiation may be actually mediated by enhancer regions at the targeted loci, rather than by the lncRNA transcripts themselves. The CRISPR-Cas9 coupled to a pgRNAs system is, therefore, a suitable tool to simultaneously target pc-genes and lncRNAs for genomic perturbation assays.

MRI phenotypes with high neurodegeneration are associated with peripheral blood B-cell changes.

Little is known about the mechanisms leading to neurodegeneration in multiple sclerosis (MS) and the role of peripheral blood cells in this neurodegenerative component. We aimed to correlate brain radiological phenotypes defined by high and low neurodegeneration with gene expression profiling of peripheral blood mononuclear cells (PBMC) from MS patients. Magnetic resonance imaging (MRI) scans from 64 patients with relapsing-remitting MS (RRMS) were classified into radiological phenotypes characterized by low (N = 27) and high (N = 37) neurodegeneration according to the number of contrast-enhancing lesions, the relative volume of non-enhancing black holes on T1-weighted images, and the brain parenchymal fraction. Gene expression profiling was determined in PBMC using microarrays, and validation of selected genes was performed by polymerase chain reaction (PCR). B-cell immunophenotyping was conducted by flow cytometry. Microarray analysis revealed the B-cell specific genes FCRL1, FCRL2, FCRL5 (Fc receptor-like 1, 2 and 5 respectively), and CD22 as the top differentially expressed genes between patients with high and low neurodegeneration. Levels for these genes were significantly down-regulated in PBMC from patients with MRI phenotypes characterized by high neurodegeneration and microarray findings were validated by PCR. In patients with high neurodegeneration, immunophenotyping showed a significant increase in the expression of the B-cell activation markers CD80 in naïve B cells (CD45+/CD19+/CD27-/IgD+), unswitched memory B cells (CD45+/CD19+/CD27+/IgD+), and switched memory B cells (CD45+/CD19+/CD27+/IgD-), and CD86 in naïve and switched memory B cells. These results suggest that RRMS patients with radiological phenotypes showing high neurodegeneration have changes in B cells characterized by down-regulation of B-cell-specific genes and increased activation status.

NLRP3 inflammasome is associated with the response to IFN-ß in patients with multiple sclerosis.

Evidence exists for a potential modulation of inflammasome activity by interferon beta. Here, we investigated the roles of inflammasomes [absent in melanoma 2 (AIM2); NLR family, CARD domain containing 4 (NLRC4); NLR family, pyrin domain containing 1 and 3 (NLRP1 and NLRP3)] and related cytokines (IL1B, IL10, IL18) in the response to interferon beta in patients with relapsing-remitting multiple sclerosis. Ninety-seven patients treated with interferon beta were classified into responders and non-responders according to clinical criteria after 24 months and clinical-radiological criteria after 12 months of treatment. Messenger RNA expression levels of inflammasomes and cytokines were determined by real-time polymerase chain reaction in peripheral blood mononuclear cells collected before treatment with interferon beta. In a subgroup of patients, NLRP3 and IL1B expression was also determined after 3 months (n = 32) and 12 months (n = 20) of interferon beta treatment. A polymorphism located in the NLRP3 gene, rs35829419, was genotyped in 789 multiple sclerosis patients treated with interferon beta. Baseline mRNA expression levels for NLRP3 and IL1B were increased in peripheral blood mononuclear cells from non-responders compared to responders classified according to clinical criteria after 24 months (P = 0.02 and P = 0.001, respectively). No significant differences were observed for other inflammasomes and related cytokines. Differences in NLRP3 and IL1B expression remained significant following a clinical-radiological classification after 12 months (P = 0.007 and P = 0.02, respectively). After treatment with interferon beta, NLRP3 and IL1B expression was increased in responders but unchanged in non-responders. A trend for association was observed between rs35829419 and interferon beta response (pM-H = 0.08). These results point to a role of the NLRP3 inflammasome and its related cytokine IL1B in the response to interferon beta in patients with relapsing-remitting multiple sclerosis.

Chitinase 3-like 1: prognostic biomarker in clinically isolated syndromes.

Chitinase 3-like 1 (CHI3L1) has been proposed as a biomarker associated with the conversion to clinically definite multiple sclerosis in patients with clinically isolated syndromes, based on the finding of increased cerebrospinal fluid CHI3L1 levels in clinically isolated syndrome patients who later converted to multiple sclerosis compared to those who remained as clinically isolated syndrome. Here, we aimed to validate CHI3L1 as a prognostic biomarker in a large cohort of patients with clinically isolated syndrome. This is a longitudinal cohort study of clinically isolated syndrome patients with clinical, magnetic resonance imaging, and cerebrospinal fluid data prospectively acquired. A total of 813 cerebrospinal fluid samples from patients with clinically isolated syndrome were recruited from 15 European multiple sclerosis centres. Cerebrospinal fluid CHI3L1 levels were measured by enzyme-linked immunosorbent assay. Multivariable Cox regression models were used to investigate the association between cerebrospinal fluid CHI3L1 levels and time to conversion to multiple sclerosis and time to reach Expanded Disability Status Scale 3.0. CHI3L1 levels were higher in patients who converted to clinically definite multiple sclerosis compared to patients who continued as clinically isolated syndrome (P = 8.1 × 10(-11)). In the Cox regression analysis, CHI3L1 levels were a risk factor for conversion to multiple sclerosis (hazard ratio = 1.7; P = 1.1 × 10(-5) using Poser criteria; hazard ratio = 1.6; P = 3.7 × 10(-6) for McDonald criteria) independent of other covariates such as brain magnetic resonance imaging abnormalities and presence of cerebrospinal fluid oligoclonal bands, and were the only significant independent risk factor associated with the development of disability (hazard ratio = 3.8; P = 2.5 × 10(-8)). High CHI3L1 levels were associated with shorter time to multiple sclerosis (P = 3.2 × 10(-9) using Poser criteria; P = 5.6 × 10(-11) for McDonald criteria) and more rapid development of disability (P = 1.8 × 10(-10)). These findings validate cerebrospinal fluid CHI3L1 as a biomarker associated with the conversion to multiple sclerosis and development of disability and reinforce the prognostic role of CHI3L1 in patients with clinically isolated syndrome. We propose that determining cerebrospinal fluid chitinase 3-like 1 levels at the time of a clinically isolated syndrome event will help identify those patients with worse disease prognosis.

Common genetic variants associated with urinary phthalate levels in children: A genome-wide study.

INTRODUCTION: Phthalates, or dieters of phthalic acid, are a ubiquitous type of plasticizer used in a variety of common consumer and industrial products. They act as endocrine disruptors and are associated with increased risk for several diseases. Once in the body, phthalates are metabolized through partially known mechanisms, involving phase I and phase II enzymes. OBJECTIVE: In this study we aimed to identify common single nucleotide polymorphisms (SNPs) and copy number variants (CNVs) associated with the metabolism of phthalate compounds in children through genome-wide association studies (GWAS). METHODS: The study used data from 1,044 children with European ancestry from the Human Early Life Exposome (HELIX) cohort. Ten phthalate metabolites were assessed in a two-void pooled urine collected at the mean age of 8 years. Six ratios between secondary and primary phthalate metabolites were calculated. Genome-wide genotyping was done with the Infinium Global Screening Array (GSA) and imputation with the Haplotype Reference Consortium (HRC) panel. PennCNV was used to estimate copy number variants (CNVs) and CNVRanger to identify consensus regions. GWAS of SNPs and CNVs were conducted using PLINK and SNPassoc, respectively. Subsequently, functional annotation of suggestive SNPs (p-value < 1E-05) was done with the FUMA web-tool. RESULTS: We identified four genome-wide significant (p-value < 5E-08) loci at chromosome (chr) 3 (FECHP1 for oxo-MiNP_oh-MiNP ratio), chr6 (SLC17A1 for MECPP_MEHHP ratio), chr9 (RAPGEF1 for MBzP), and chr10 (CYP2C9 for MECPP_MEHHP ratio). Moreover, 115 additional loci were found at suggestive significance (p-value < 1E-05). Two CNVs located at chr11 (MRGPRX1 for oh-MiNP and SLC35F2 for MEP) were also identified. Functional annotation pointed to genes involved in phase I and phase II detoxification, molecular transfer across membranes, and renal excretion. CONCLUSION: Through genome-wide screenings we identified known and novel loci implicated in phthalate metabolism in children. Genes annotated to these loci participate in detoxification, transmembrane transfer, and renal excretion.

An encyclopedia of enhancer-gene regulatory interactions in the human genome.

Identifying transcriptional enhancers and their target genes is essential for understanding gene regulation and the impact of human genetic variation on disease1-6. Here we create and evaluate a resource of >13 million enhancer-gene regulatory interactions across 352 cell types and tissues, by integrating predictive models, measurements of chromatin state and 3D contacts, and largescale genetic perturbations generated by the ENCODE Consortium7. We first create a systematic benchmarking pipeline to compare predictive models, assembling a dataset of 10,411 elementgene pairs measured in CRISPR perturbation experiments, >30,000 fine-mapped eQTLs, and 569 fine-mapped GWAS variants linked to a likely causal gene. Using this framework, we develop a new predictive model, ENCODE-rE2G, that achieves state-of-the-art performance across multiple prediction tasks, demonstrating a strategy involving iterative perturbations and supervised machine learning to build increasingly accurate predictive models of enhancer regulation. Using the ENCODE-rE2G model, we build an encyclopedia of enhancer-gene regulatory interactions in the human genome, which reveals global properties of enhancer networks, identifies differences in the functions of genes that have more or less complex regulatory landscapes, and improves analyses to link noncoding variants to target genes and cell types for common, complex diseases. By interpreting the model, we find evidence that, beyond enhancer activity and 3D enhancer-promoter contacts, additional features guide enhancerpromoter communication including promoter class and enhancer-enhancer synergy. Altogether, these genome-wide maps of enhancer-gene regulatory interactions, benchmarking software, predictive models, and insights about enhancer function provide a valuable resource for future studies of gene regulation and human genetics.

Treatment with MOG-DNA vaccines induces CD4+CD25+FoxP3+ regulatory T cells and up-regulates genes with neuroprotective functions in experimental autoimmune encephalomyelitis.

BACKGROUND: DNA vaccines represent promising therapeutic strategies in autoimmune disorders such as multiple sclerosis (MS). However, the precise mechanisms by which DNA vaccines induce immune regulation remain largely unknown. Here, we aimed to expand previous knowledge existing on the mechanisms of action of DNA vaccines in the animal model of MS, experimental autoimmune encephalomyelitis (EAE), by treating EAE mice with a DNA vaccine encoding the myelin oligodendrocyte glycoprotein (MOG), and exploring the therapeutic effects on the disease-induced inflammatory and neurodegenerative changes. METHODS: EAE was induced in C57BL6/J mice by immunization with MOG35₋55 peptide. Mice were intramuscularly treated with a MOG-DNA vaccine or vehicle in prophylactic and therapeutic approaches. Histological studies were performed in central nervous system (CNS) tissue. Cytokine production and regulatory T cell (Treg) quantification were achieved by flow cytometry. Gene expression patterns were determined using microarrays, and the main findings were validated by real-time PCR. RESULTS: MOG-DNA treatment reduced the clinical and histopathological signs of EAE when administered in both prophylactic and therapeutic settings. Suppression of clinical EAE was associated with dampening of antigen (Ag)-specific proinflammatory Th1 and Th17 immune responses and, interestingly, expansion of Treg in the periphery and upregulation in the CNS of genes encoding neurotrophic factors and proteins involved in remyelination. CONCLUSIONS: These results suggest for the first time that the beneficial effects of DNA vaccines in EAE are not limited to anti-inflammatory mechanisms, and DNA vaccines may also exert positive effects through hitherto unknown neuroprotective mechanisms.

Tyrosine kinase 2 variant influences T lymphocyte polarization and multiple sclerosis susceptibility.

The tyrosine kinase 2 variant rs34536443 has been established as a genetic risk factor for multiple sclerosis in a variety of populations. However, the functional effect of this variant on disease pathogenesis remains unclear. This study replicated the genetic association of tyrosine kinase 2 with multiple sclerosis in a cohort of 1366 French patients and 1802 controls. Furthermore, we assessed the functional consequences of this polymorphism on human T lymphocytes by comparing the reactivity and cytokine profile of T lymphocytes isolated from individuals expressing the protective TYK2(GC) genotype with the disease-associated TYK2(GG) genotype. Our results demonstrate that the protective C allele infers decreased tyrosine kinase 2 activity, and this reduction of activity is associated with a shift in the cytokine profile favouring the secretion of Th2 cytokines. These findings suggest that the rs34536443 variant effect on multiple sclerosis susceptibility might be mediated by deviating T lymphocyte differentiation toward a Th2 phenotype. This impact of tyrosine kinase 2 on effector differentiation is likely to be of wider importance because other autoimmune diseases also have been associated with polymorphisms within tyrosine kinase 2. The modulation of tyrosine kinase 2 activity might therefore represent a new therapeutic approach for the treatment of autoimmune diseases.

PLANdbAffy: probe-level annotation database for Affymetrix expression microarrays.

Standard Affymetrix technology evaluates gene expression by measuring the intensity of mRNA hybridization with a panel of the 25-mer oligonucleotide probes, and summarizing the probe signal intensities by a robust average method. However, in many cases, signal intensity of the probe does not correlate with gene expression. This could be due to the hybridization of the probe to a transcript of another gene, mapping of the probe to an intron, alternative splicing, single nucleotide polymorphisms and other reasons. We have developed a database, PLANdbAffy (available at http://affymetrix2.bioinf.fbb.msu.ru), that contains the results of the alignment of probe sequences from five Affymetrix expression microarrays to the human genome. We have determined the probes matching the transcript-coding regions in the correct orientation. For each such probe alignment region, we determined the mRNA and EST sequences that contain the probe sequence. In the textual part of the database interface we summarize the data on the sequences that cover the probe alignment region and SNPs that are located inside it. The graphical part of our database interface is implemented as custom tracks to the UCSC genome browser that allows one to utilize all the data that are offered by UCSC browser.

Positive selection in alternatively spliced exons of human genes.

Alternative splicing is a well-recognized mechanism of accelerated genome evolution. We have studied single-nucleotide polymorphisms and human-chimpanzee divergence in the exons of 6672 alternatively spliced human genes, with the aim of understanding the forces driving the evolution of alternatively spliced sequences. Here, we show that alternatively spliced exons and exon fragments (alternative exons) from minor isoforms experience lower selective pressure at the amino acid level, accompanied by selection against synonymous sequence variation. The results of the McDonald-Kreitman test suggest that alternatively spliced exons, unlike exons constitutively included in the mRNA, are also subject to positive selection, with up to 27% of amino acids fixed by positive selection.

NEAT1 Long Isoform Is Highly Expressed in Chronic Lymphocytic Leukemia Irrespectively of Cytogenetic Groups or Clinical Outcome.

The biological role and therapeutic potential of long non-coding RNAs (lncRNAs) in chronic lymphocytic leukemia (CLL) are still open questions. Herein, we investigated the significance of the lncRNA NEAT1 in CLL. We examined NEAT1 expression in 310 newly diagnosed Binet A patients, in normal CD19+ B-cells, and other types of B-cell malignancies. Although global NEAT1 expression level was not statistically different in CLL cells compared to normal B cells, the median ratio of NEAT1_2 long isoform and global NEAT1 expression in CLL samples was significantly higher than in other groups. NEAT1_2 was more expressed in patients carrying mutated IGHV genes. Concerning cytogenetic aberrations, NEAT1_2 expression in CLL with trisomy 12 was lower with respect to patients without alterations. Although global NEAT1 expression appeared not to be associated with clinical outcome, patients with the lowest NEAT1_2 expression displayed the shortest time to first treatment; however, a multivariate regression analysis showed that the NEAT1_2 risk model was not independent from other known prognostic factors, particularly the IGHV mutational status. Overall, our data prompt future studies to investigate whether the increased amount of the long NEAT1_2 isoform detected in CLL cells may have a specific role in the pathology of the disease.

Rodent-specific alternative exons are more frequent in rapidly evolving genes and in paralogs.

BACKGROUND: Alternative splicing is an important mechanism for generating functional and evolutionary diversity of proteins in eukaryotes. Here, we studied the frequency and functionality of recently gained, rodent-specific alternative exons. RESULTS: We projected the data about alternative splicing of mouse genes to the rat, human, and dog genomes, and identified exons conserved in the rat genome, but missing in more distant genomes. We estimated the frequency of rodent-specific exons while controlling for possible residual conservation of spurious exons. The frequency of rodent-specific exons is higher among predominantly skipped exons and exons disrupting the reading frame. Separation of all genes by the rate of sequence evolution and by gene families has demonstrated that rodent-specific cassette exons are more frequent in rapidly evolving genes and in rodent-specific paralogs. CONCLUSION: Thus we demonstrated that recently gained exons tend to occur in fast-evolving genes, and their inclusion rate tends to be lower than that of older exons. This agrees with the theory that gain of alternative exons is one of the major mechanisms of gene evolution.

The effects of death and post-mortem cold ischemia on human tissue transcriptomes.

Post-mortem tissues samples are a key resource for investigating patterns of gene expression. However, the processes triggered by death and the post-mortem interval (PMI) can significantly alter physiologically normal RNA levels. We investigate the impact of PMI on gene expression using data from multiple tissues of post-mortem donors obtained from the GTEx project. We find that many genes change expression over relatively short PMIs in a tissue-specific manner, but this potentially confounding effect in a biological analysis can be minimized by taking into account appropriate covariates. By comparing ante- and post-mortem blood samples, we identify the cascade of transcriptional events triggered by death of the organism. These events do not appear to simply reflect stochastic variation resulting from mRNA degradation, but active and ongoing regulation of transcription. Finally, we develop a model to predict the time since death from the analysis of the transcriptome of a few readily accessible tissues.

Conserved and species-specific alternative splicing in mammalian genomes.

BACKGROUND: Alternative splicing has been shown to be one of the major evolutionary mechanisms for protein diversification and proteome expansion, since a considerable fraction of alternative splicing events appears to be species- or lineage-specific. However, most studies were restricted to the analysis of cassette exons in pairs of genomes and did not analyze functionality of the alternative variants. RESULTS: We analyzed conservation of human alternative splice sites and cassette exons in the mouse and dog genomes. Alternative exons, especially minor-isofom ones, were shown to be less conserved than constitutive exons. Frame-shifting alternatives in the protein-coding regions are less conserved than frame-preserving ones. Similarly, the conservation of alternative sites is highest for evenly used alternatives, and higher when the distance between the sites is divisible by three. The rate of alternative-exon and site loss in mouse is slightly higher than in dog, consistent with faster evolution of the former. The evolutionary dynamics of alternative sites was shown to be consistent with the model of random activation of cryptic sites. CONCLUSION: Consistent with other studies, our results show that minor cassette exons are less conserved than major-alternative and constitutive exons. However, our study provides evidence that this is caused not only by exon birth, but also lineage-specific loss of alternative exons and sites, and it depends on exon functionality.

Overlapping alternative donor splice sites in the human genome.

Over 50% of donor splice sites in the human genome have a potential alternative donor site at a distance of three to six nucleotides. Conservation of these potential sites is determined by the consensus requirements and by its exonic or intronic location. Several hundred pairs of overlapping sites are confirmed to be alternatively spliced as both sites in a pair are supported by a protein, by a full-length mRNA, or by expressed sequence tags (ESTs) from at least two independent clone libraries. Overlapping sites may clash with consensus requirements. Pairs with a site shift of four nucleotides are the most abundant, despite the frameshift in the protein-coding region that they introduce. The site usage in pairs is usually uneven, and the major site is more frequently conserved in other mammalian genomes. Overlapping alternative donor sites and acceptor sites may have different functional roles: alternative splicing of overlapping acceptor sites leads mainly to microvariations in protein sequences; whereas alternative donor sites often lead to frameshifts and thus either yield major differences in the protein sequence and structure, or generate nonsense-mediated decay-inducing mRNA isoforms likely involved in regulated unproductive splicing pathways.

Fast rate of evolution in alternatively spliced coding regions of mammalian genes.

BACKGROUND: At least half of mammalian genes are alternatively spliced. Alternative isoforms are often genome-specific and it has been suggested that alternative splicing is one of the major mechanisms for generating protein diversity in the course of evolution. Another way of looking at alternative splicing is to consider sequence evolution of constitutive and alternative regions of protein-coding genes. Indeed, it turns out that constitutive and alternative regions evolve in different ways. RESULTS: A set of 3029 orthologous pairs of human and mouse alternatively spliced genes was considered. The rate of nonsynonymous substitutions (dN), the rate of synonymous substitutions (dS), and their ratio (omega = dN/dS) appear to be significantly higher in alternatively spliced coding regions compared to constitutive regions. When N-terminal, internal and C-terminal alternatives are analysed separately, C-terminal alternatives appear to make the main contribution to the observed difference. The effects become even more pronounced in a subset of fast evolving genes. CONCLUSION: These results provide evidence of weaker purifying selection and/or stronger positive selection in alternative regions and thus one more confirmation of accelerated evolution in alternative regions. This study corroborates the theory that alternative splicing serves as a testing ground for molecular evolution.

No statistical support for correlation between the positions of protein interaction sites and alternatively spliced regions.

BACKGROUND: Alternative splicing is an efficient mechanism for increasing the variety of functions fulfilled by proteins in a living cell. It has been previously demonstrated that alternatively spliced regions often comprise functionally important and conserved sequence motifs. The objective of this work was to test the hypothesis that alternative splicing is correlated with contact regions of protein-protein interactions. RESULTS: Protein sequence spans involved in contacts with an interaction partner were delineated from atomic structures of transient interaction complexes and juxtaposed with the location of alternatively spliced regions detected by comparative genome analysis and spliced alignment. The total of 42 alternatively spliced isoforms were identified in 21 amino acid chains involved in biomolecular interactions. Using this limited dataset and a variety of sophisticated counting procedures we were not able to establish a statistically significant correlation between the positions of protein interaction sites and alternatively spliced regions. CONCLUSIONS: This finding contradicts a naïve hypothesis that alternatively spliced regions would correlate with points of contact. One possible explanation for that could be that all alternative splicing events change the spatial structure of the interacting domain to a sufficient degree to preclude interaction. This is indirectly supported by the observed lack of difference in the behaviour of relatively short regions affected by alternative splicing and cases when large portions of proteins are removed. More structural data on complexes of interacting proteins, including structures of alternative isoforms, are needed to test this conjecture.

Baseline gene expression signatures in monocytes from multiple sclerosis patients treated with interferon-beta.

BACKGROUND: A relatively large proportion of relapsing-remitting multiple sclerosis (RRMS) patients do not respond to interferon-beta (IFNb) treatment. In previous studies with peripheral blood mononuclear cells (PBMC), we identified a subgroup of IFNb non-responders that was characterized by a baseline over-expression of type I IFN inducible genes. Additional mechanistic experiments carried out in IFNb non-responders suggested a selective alteration of the type I IFN signaling pathway in the population of blood monocytes. Here, we aimed (i) to investigate whether the type I IFN signaling pathway is up-regulated in isolated monocytes from IFNb non-responders at baseline; and (ii) to search for additional biological pathways in this cell population that may be implicated in the response to IFNb treatment. METHODS: Twenty RRMS patients classified according to their clinical response to IFNb treatment and 10 healthy controls were included in the study. Monocytes were purified from PBMC obtained before treatment by cell sorting and the gene expression profiling was determined with oligonucleotide microarrays. RESULTS AND DISCUSSION: Purified monocytes from IFNb non-responders were characterized by an over-expression of type I IFN responsive genes, which confirms the type I IFN signature in monocytes suggested from previous studies. Other relevant signaling pathways that were up-regulated in IFNb non-responders were related with the mitochondrial function and processes such as protein synthesis and antigen presentation, and together with the type I IFN signaling pathway, may also be playing roles in the response to IFNb.

Search for specific biomarkers of IFNß bioactivity in patients with multiple sclerosis.

Myxovirus A (MxA), a protein encoded by the MX1 gene with antiviral activity, has proven to be a sensitive measure of IFNß bioactivity in multiple sclerosis (MS). However, the use of MxA as a biomarker of IFNß bioactivity has been criticized for the lack of evidence of its role on disease pathogenesis and the clinical response to IFNß. Here, we aimed to identify specific biomarkers of IFNß bioactivity in order to compare their gene expression induction by type I IFNs with the MxA, and to investigate their potential role in MS pathogenesis. Gene expression microarrays were performed in PBMC from MS patients who developed neutralizing antibodies (NAB) to IFNß at 12 and/or 24 months of treatment and patients who remained NAB negative. Nine genes followed patterns in gene expression over time similar to the MX1, which was considered the gold standard gene, and were selected for further experiments: IFI6, IFI27, IFI44L, IFIT1, HERC5, LY6E, RSAD2, SIGLEC1, and USP18. In vitro experiments in PBMC from healthy controls revealed specific induction of selected biomarkers by IFNß but not IFNγ, and several markers, in particular USP18 and HERC5, were shown to be significantly induced at lower IFNß concentrations and more selective than the MX1 as biomarkers of IFNß bioactivity. In addition, USP18 expression was deficient in MS patients compared with healthy controls (p = 0.0004). We propose specific biomarkers that may be considered in addition to the MxA to evaluate IFNß bioactivity, and to further explore their implication in MS pathogenesis.

Probe-Level Universal Search (PLUS) algorithm for gender differentiation in affymetrix datasets.

Affymetrix microarrays measure gene expression based on the intensity of hybridization of a panel of oligonucleotide probes (probe set) with mRNA. The signals from all probes within a probe set are converted into a single measure that represents the expression value of a gene. This step diminishes the number of independently measured parameters and eliminates from consideration individual "good-working" probes. We propose a new feature selection algorithm (Probe Level Universal Search or PLUS algorithm) for probe-level analysis of gene expression datasets. The algorithm evaluates the intensities of perfect-match Affymetrix probes individually and selects probes that allow one to distinguish two given classes of samples. The algorithm was used to differentiate the samples according to their gender ("gender differentiation"). The universal gender differentiating set of 3' Gene Affymetrix microarray probes was selected; the set consists of 38 probes from XIST gene of X-chromosome and 17 probes from five Y-chromosome genes: RPS4Y1, EIF1A, DDX3Y, JARID1D and USP9Y. The selection procedure based on the probes selected by PLUS algorithm differentiates the sex chromosome karyotype of the sample, reveals samples with incorrect gender labels and samples from patients with hereditary syndromes or cancer-associated chromosome abnormalities.