Survivin prevents the polycomb repressor complex 2 from methylating histone 3 lysine 27.

This study investigates the role of survivin in epigenetic control of gene transcription through interaction with the polycomb repressive complex 2 (PRC2). PRC2 is responsible for silencing gene expression by trimethylating lysine 27 on histone 3. We observed differential expression of PRC2 subunits in CD4+ T cells with varying levels of survivin expression, and ChIP-seq results indicated that survivin colocalizes with PRC2 along DNA. Inhibition of survivin resulted in a significant increase in H3K27 trimethylation, implying that survivin prevents PRC2 from functioning. Peptide microarray showed that survivin interacts with peptides from PRC2 subunits, and machine learning revealed that amino acid composition contains relevant information for predicting survivin interaction. NMR and BLI experiments supported the interaction of survivin with PRC2 subunit EZH2. Finally, protein-protein docking revealed that the survivin-EZH2 interaction interface overlaps with catalytic residues of EZH2, potentially inhibiting its H3K27 methylation activity. These findings suggest that survivin inhibits PRC2 function.

Survivin promotes a glycolytic switch in CD4+ T cells by suppressing the transcription of PFKFB3 in rheumatoid arthritis.

In this study, we explore the role of nuclear survivin in maintaining the effector phenotype of IFNγ-producing T cells acting through the transcriptional control of glucose utilization. High expression of survivin in CD4+T cells was associated with IFNγ-dependent phenotype and anaerobic glycolysis. Transcriptome of CD4+ cells and sequencing of survivin-bound chromatin showed that nuclear survivin had a genome-wide and motif-specific binding to regulatory regions of the genes controlling cell metabolism. Survivin coprecipitates with transcription factors IRF1 and SMAD3, which repressed the transcription of the metabolic check-point enzyme phosphofructokinase 2 gene PFKFB3 and promoted anaerobic glycolysis. Combining transcriptome analyses of CD4+ cells and functional studies in glucose metabolism, we demonstrated that the inhibition of survivin reverted PFKFB3 production, inhibited glucose uptake, and reduces interferon effects in CD4+ cells. These results present a survivin-dependent mechanism in coordinating the metabolic adaptation of CD4+T cells and propose an attractive strategy to counteract IFNγ-dependent inflammation in autoimmunity.

Cohesin-Mediated Chromatin Interactions and Autoimmunity.

Proper physiological functioning of any cell type requires ordered chromatin organization. In this context, cohesin complex performs important functions preventing premature separation of sister chromatids after DNA replication. In partnership with CCCTC-binding factor, it ensures insulator activity to organize enhancers and promoters within regulatory chromatin. Homozygous mutations and dysfunction of individual cohesin proteins are embryonically lethal in humans and mice, which limits in vivo research work to embryonic stem cells and progenitors. Conditional alleles of cohesin complex proteins have been generated to investigate their functional roles in greater detail at later developmental stages. Thus, genome regulation enabled by action of cohesin proteins is potentially crucial in lineage cell development, including immune homeostasis. In this review, we provide current knowledge on the role of cohesin complex in leukocyte maturation and adaptive immunity. Conditional knockout and shRNA-mediated inhibition of individual cohesin proteins in mice demonstrated their importance in haematopoiesis, adipogenesis and inflammation. Notably, these effects occur rather through changes in transcriptional gene regulation than through expected cell cycle defects. This positions cohesin at the crossroad of immune pathways including NF-kB, IL-6, and IFNγ signaling. Cohesin proteins emerged as vital regulators at early developmental stages of thymocytes and B cells and after antigen challenge. Human genome-wide association studies are remarkably concordant with these findings and present associations between cohesin and rheumatoid arthritis, multiple sclerosis and HLA-B27 related chronic inflammatory conditions. Furthermore, bioinformatic prediction based on protein-protein interactions reveal a tight connection between the cohesin complex and immune relevant processes supporting the notion that cohesin will unearth new clues in regulation of autoimmunity.

Glioblastoma gene network reconstruction and ontology analysis by online bioinformatics tools.

Glioblastoma is the most aggressive type of brain tumors resistant to a number of antitumor drugs. The problem of therapy and drug treatment course is complicated by extremely high heterogeneity in the benign cell populations, the random arrangement of tumor cells, and polymorphism of their nuclei. The pathogenesis of gliomas needs to be studied using modern cellular technologies, genome- and transcriptome-wide technologies of high-throughput sequencing, analysis of gene expression on microarrays, and methods of modern bioinformatics to find new therapy targets. Functional annotation of genes related to the disease could be retrieved based on genetic databases and cross-validated by integrating complementary experimental data. Gene network reconstruction for a set of genes (proteins) proved to be effective approach to study mechanisms underlying disease progression. We used online bioinformatics tools for annotation of gene list for glioma, reconstruction of gene network and comparative analysis of gene ontology categories. The available tools and the databases for glioblastoma gene analysis are discussed together with the recent progress in this field.

Prognostic Significance of BIRC5/Survivin in Breast Cancer: Results from Three Independent Cohorts.

Breast cancer (BC) histological and molecular classifications significantly improved the treatment strategy and prognosis. Inhibitor of apoptosis BIRC5/survivin is often overexpressed in cancers, however, indications of its importance in BC are inconsistent. We integrate BIRC5 protein and mRNA measures with clinical associates and long-term outcome in three independent cohorts Protein levels of BIRC5 were measured in primary lysates of 845 patients of the West Swedish BC cohort (VGR-BC) and linked to 5- and 27-years survival. The results were externally validated in transcriptomic data from METABRIC and SCAN-B cohorts. Survival analysis showed that high levels of BIRC5 were consistently associated with a poor probability of 5-year overall survival. High BIRC5 in VGR-BC contributed negatively to the disease-specific survival at 5 and 27 years. Subsets with different status by ER (estrogen receptor) expression and presence of nodal metastasis supported independent association of high BIRC5 with poor prognosis in all cohorts. In METABRIC and SCAN-B cohorts, high levels of BIRC5 mRNA were associated with the basal-like and luminal B molecular BC subtypes and with increasing histologic grade. BIRC5 is a sensitive survival marker that acts independent of ER and nodal status, and its levels need to be considered when making treatment decisions.

Antibodies against Phosphorylcholine and Malondialdehyde during the First Two Years of Life.

Abs against phosphorylcholine (anti-PC) and Abs against malondialdehyde (anti-MDA) may be protective in chronic inflammation, like atherosclerosis and cardiovascular disease. It is not known how they develop early in life. Ab titers were measured using ELISA in healthy women (n = 105; born into life study) and their children. Plasma samples were collected from the mothers before conception and from the children at birth as well as at 1 and 2 y after birth. Extracted Abs were compared using a proteomics de novo sequencing approach. It was observed that children were born with very low levels of IgM anti-PC, whereas IgM anti-MDA was present at birth. Both IgM anti-PC and anti-MDA increased during the first 2 y of life, but IgM anti-PC in contrast to IgM anti-MDA was still significantly lower than in the mothers. IgG anti-PC decreased after 1 y but reached similar levels as mothers' after 2 y, whereas IgG anti-MDA reached similar levels as mothers' already after 1 y. Proteomics peptide sequencing analysis indicated large peptide sequence variation without specific clone expression during the early stage of life compared with the adult stage for which specific peptide sequences dominated. IgM anti-PC levels develop much slower than anti-MDA and are still relatively low at 2 y. We hypothesize that anti-PC is developed by a combination of preprogramming and exposure to the external world, in which infectious agents may play a role. For anti-MDA, preprogramming is likely to play a major role and at an earlier stage than for anti-PC.

IgM antibodies against malondialdehyde and phosphorylcholine in different systemic rheumatic diseases.

IgM antibodies against phosphorylcholine (anti-PC) and malondialdehyde (anti-MDA) may have protective properties in cardiovascular and rheumatic diseases. We here compare these antibodies in systemic rheumatic conditions and study their properties. Anti-PC and anti-MDA was measured using ELISA in patients with SLE (374), RA (354), Mixed connective tissue disease (MCTD, 77), Systemic sclerosis (SSc, 331), Sjögren's syndrome (SjS, 324), primary antiphospholipid syndrome (PAPs, 65), undifferentiated connective tissue disease (UCTD, 118) and 515 matched healthy controls (HC). Cardiovascular score (CV) was broadly defined based on clinical disease symptoms. Anti-PC and anti-MDA peptide/protein characterization were compared using a proteomics de novo sequencing approach. anti-MDA and anti-PC were extracted from total IgM. The proportion of Treg cells was determined by flow cytometry. The maximal difference between cases and controls was shown for MCTD: significantly lower IgM Anti-PC but not anti-MDA among patients (median 49.3RU/ml vs 70.4 in healthy controls, p(t-test) = 0.0037). IgM low levels were more prevalent in MCTD, SLE, SjS, SSc and UCTD. IgM anti-PC variable region profiles were different from and more homologous than anti-MDA. Anti-PC but not anti-MDA were significantly negatively correlated with CV in the whole patient group. In contrast to IgM anti-PC, anti-MDA did not promote polarization of Tregs. Taken together, Anti-PC is decreased in MCTD and also in SLE, SjS and SSc but not in other studied diseases. Anti-PC may thus differentiate between these. In contrast, anti-MDA did not show these differences between diseases studied. Anti-PC level is negatively correlated with CV in the patient group cohort. In contrast to anti-PC, anti-MDA did not promote Treg polarization. These findings could have both diagnostic and therapeutic implications, one possibility being active or passive immunization with PC in some rheumatic conditions.

An update on the genetics of systemic lupus erythematosus.

PURPOSE OF REVIEW: The aim of this study is to update on the most recent findings on the genetics of systemic lupus erythematosus. RECENT FINDINGS: Our overview focuses particularly on results from expression quantitative trait loci, exome sequencing, and rare variants and their impact on disease. SUMMARY: Systemic lupus erythematosus is a systemic autoimmune disease for which a significant number of susceptibility genes have been identified. Several genome-wide association studies were recently published in different populations that provide a better picture of the molecular mechanisms. It is becoming clear that the genetic architecture of lupus is quite well established but more information is required on the role of rare variants.

Trans-Ethnic Mapping of BANK1 Identifies Two Independent SLE-Risk Linkage Groups Enriched for Co-Transcriptional Splicing Marks.

BANK1 is a susceptibility gene for several systemic autoimmune diseases in several populations. Using the genome-wide association study (GWAS) data from Europeans (EUR) and African Americans (AA), we performed an extensive fine mapping of ankyrin repeats 1 (BANK1). To increase the SNP density, we used imputation followed by univariate and conditional analysis, combined with a haplotypic and expression quantitative trait locus (eQTL) analysis. The data from Europeans showed that the associated region was restricted to a minimal and dependent set of SNPs covering introns two and three, and exon two. In AA, the signal found in the Europeans was split into two independent effects. All of the major risk associated SNPs were eQTLs, and the risks were associated with an increased BANK1 gene expression. Functional annotation analysis revealed the enrichment of repressive B cell epigenomic marks (EZH2 and H3K27me3) and a strong enrichment of splice junctions. Furthermore, one eQTL located in intron two, rs13106926, was found within the binding site for RUNX3, a transcriptional activator. These results connect the local genome topography, chromatin structure, and the regulatory landscape of BANK1 with co-transcriptional splicing of exon two. Our data defines a minimal set of risk associated eQTLs predicted to be involved in the expression of BANK1 modulated through epigenetic regulation and splicing. These findings allow us to suggest that the increased expression of BANK1 will have an impact on B-cell mediated disease pathways.

Whole Exome Sequencing of Patients from Multicase Families with Systemic Lupus Erythematosus Identifies Multiple Rare Variants.

In an effort to identify rare alleles associated with SLE, we have performed whole exome sequencing of the most distantly related affected individuals from two large Icelandic multicase SLE families followed by Ta targeted genotyping of additional relatives. We identified multiple rare likely pathogenic variants in nineteen genes co-segregating with the disease through multiple generations. Gene co-expression and protein-protein interaction analysis identified a network of highly connected genes comprising several loci previously implicated in autoimmune diseases. These genes were significantly enriched for immune system development, lymphocyte activation, DNA repair, and V(D)J gene recombination GO-categories. Furthermore, we found evidence of aggregate association and enrichment of rare variants at the FAM71E1/EMC10 locus in an independent set of 4,254 European SLE-cases and 4,349 controls. Our study presents evidence supporting that multiple rare likely pathogenic variants, in newly identified genes involved in known disease pathogenic pathways, segregate with SLE at the familial and population level.

Prediction of amino acid positions specific for functional groups in a protein family based on local sequence similarity.

The exchange of single amino acid residue in protein can substantially affect the specificity of molecular recognition. Many protein families can be divided into the groups based on specificity to recognized ligands. Prediction of group-discriminating residues within the certain family is extremely necessary for theoretical studies, enzyme engineering, drug design, and so on. The most existing methods use the multiple sequence alignment. They have the limitations in prediction accuracy due to the family sequence divergence and ligand-based grouping. We developed a new method SPrOS (Specificity Projection On Sequence) for estimating the specificity of residues to user-defined groups. SPrOS compares the sequence segments from the test protein and training proteins. Contrary to other segment-comparison approaches extracting the string motifs, SPrOS calculates the scores for single positions by the similarity of their surroundings. The method was evaluated on the simulated sequences and real protein families. The high-prediction accuracy was achieved for simulated sequences, in which SPrOS detected specific positions not predicted with the alignment-based method. For bacterial transcription factors (LacI/GalR) clearly divided into functional groups, the predicted specific residues corresponded to the published experimental data. In a more complicated case of protein kinases classified by inhibitor specificity, the positions predicted with high significance were located in ligand-binding areas. As the ligand specificity is not necessary coincided with phylogeny, evolutionary-coupled mutations could disturb the detection of ligand-specific residues. Excluding proximate homologs of the test protein kinase from the training set, we improved the prediction of the ligand-specific residues. The SPrOS is available at http://www.way2drug.com/spros/

PXK locus in systemic lupus erythematosus: fine mapping and functional analysis reveals novel susceptibility gene ABHD6.

OBJECTIVES: To perform fine mapping of the PXK locus associated with systemic lupus erythematosus (SLE) and study functional effects that lead to susceptibility to the disease. METHODS: Linkage disequilibrium (LD) mapping was conducted by using 1251 SNPs (single nucleotide polymorphism) covering a 862 kb genomic region on 3p14.3 comprising the PXK locus in 1467 SLE patients and 2377 controls of European origin. Tag SNPs and genotypes imputed with IMPUTE2 were tested for association by using SNPTEST and PLINK. The expression QTLs data included three independent datasets for lymphoblastoid cells of European donors: HapMap3, MuTHER and the cross-platform eQTL catalogue. Correlation analysis of eQTLs was performed using Vassarstats. Alternative splicing for the PXK gene was analysed on mRNA from PBMCs. RESULTS: Fine mapping revealed long-range LD (>200 kb) extended over the ABHD6, RPP14, PXK, and PDHB genes on 3p14.3. The highly correlated variants tagged an SLE-associated haplotype that was less frequent in the patients compared with the controls (OR=0.89, p=0.00684). A robust correlation between the association with SLE and enhanced expression of ABHD6 gene was revealed, while neither expression, nor splicing alterations associated with SLE susceptibility were detected for PXK. The SNP allele frequencies as well as eQTL pattern analysed in the CEU and CHB HapMap3 populations indicate that the SLE association and the effect on ABHD6 expression are specific to Europeans. CONCLUSIONS: These results confirm the genetic association of the locus 3p14.3 with SLE in Europeans and point to the ABHD6 and not PXK, as the major susceptibility gene in the region. We suggest a pathogenic mechanism mediated by the upregulation of ABHD6 in individuals carrying the SLE-risk variants.

Non-random DNA fragmentation in next-generation sequencing.

Next Generation Sequencing (NGS) technology is based on cutting DNA into small fragments, and their massive parallel sequencing. The multiple overlapping segments termed "reads" are assembled into a contiguous sequence. To reduce sequencing errors, every genome region should be sequenced several dozen times. This sequencing approach is based on the assumption that genomic DNA breaks are random and sequence-independent. However, previously we showed that for the sonicated restriction DNA fragments the rates of double-stranded breaks depend on the nucleotide sequence. In this work we analyzed genomic reads from NGS data and discovered that fragmentation methods based on the action of the hydrodynamic forces on DNA, produce similar bias. Consideration of this non-random DNA fragmentation may allow one to unravel what factors and to what extent influence the non-uniform coverage of various genomic regions.

Altered expression of multiple genes involved in retinoic acid biosynthesis in human colorectal cancer.

All-trans-retinoic acid (atRA), the oxidized form of vitamin A (retinol), regulates a wide variety of biological processes, such as cell proliferation and differentiation. Multiple alcohol, retinol and retinaldehyde dehydrogenases (ADHs, RDHs, RALDHs) as well as aldo-keto reductases (AKRs) catalyze atRA production. The reduced atRA biosynthesis has been observed in several human tumors, including colorectal cancer. However, subsets of atRA-synthesizing enzymes have not been determined in colorectal tumors. We investigated the expression patterns of genes involved in atRA biosynthesis in normal human colorectal tissues, primary carcinomas and cancer cell lines by RT-PCR. These genes were identified using transcriptomic data analysis (expressed sequence tags, RNA-sequencing, microarrays). Our results indicate that each step of the atRA biosynthesis pathway is dysregulated in colorectal cancer. Frequent and significant decreases in the mRNA levels of the ADH1B, ADH1C, RDHL, RDH5 and AKR1B10 genes were observed in a majority of colorectal carcinomas. The expression levels of the RALDH1 gene were reduced, and the expression levels of the cytochrome CYP26A1 gene increased. The human colon cancer cell lines showed a similar pattern of changes in the mRNA levels of these genes. A dramatic reduction in the expression of genes encoding the predominant retinol-oxidizing enzymes could impair atRA production. The most abundant of these genes, ADH1B and ADH1C, display decreased expression during progression from adenoma to early and more advanced stage of colorectal carcinomas. The diminished atRA biosynthesis may lead to alteration of cell growth and differentiation in the colon and rectum, thus contributing to the progression of colorectal cancer.

Mutation in ribosomal protein S5 leads to spectinomycin resistance in Neisseria gonorrhoeae.

Spectinomycin remains a useful reserve option for therapy of gonorrhea. The emergence of multidrug-resistant Neisseria gonorrhoeae strains with decreased susceptibility to cefixime and to ceftriaxone makes it the only medicine still effective for treatment of gonorrhea infection in analogous cases. However, adoption of spectinomycin as a routinely used drug of choice was soon followed by reports of spectinomycin resistance. The main molecular mechanism of spectinomycin resistance in N. gonorrhoeae was C1192T substitution in 16S rRNA genes. Here we reported a Thr-24→Pro mutation in ribosomal protein S5 (RPS5) found in spectinomycin resistant clinical N. gonorrhoeae strain, which carried no changes in 16S rRNA. In a series of experiments, the transfer of rpsE gene allele encoding the mutant RPS5 to the recipient N. gonorrhoeae strains was analyzed. The relatively high rate of transformation [ca. 10(-5) colony-forming units (CFUs)] indicates the possibility of spread of spectinonycin resistance within gonococcal population due to the horizontal gene transfer (HGT).

Prediction of amino acid residues participated in substrate recognition by cytochrome P450 subfamilies with broad substrate specificity.

Cytochromes P450 comprise a large superfamily and several of their isoforms play a crucial role in metabolism of xenobiotics, including drugs. Although these enzymes demonstrate broad and cross-substrate specificity, different cytochrome P450 subfamilies exhibit certain selectivity for some types of substrates. Analysis of amino acid residues of the active sites of six cytochrome subfamilies (CYP1А, CYP2А, CYP2С, CYP2D, CYP2E and CYP3А) enables to define subfamily-specific patterns that consist of four residues. These residues are located on the periphery of the active sites of these cytochromes. We suggest that they can form a primary binding site at the entrance to the active site, defining cytochrome substrate recognition.

Impact of recombination on polymorphism of genes encoding Kunitz-type protease inhibitors in the genus Solanum.

BACKGROUND: The group of Kunitz-type protease inhibitors (KPI) from potato is encoded by a polymorphic family of multiple allelic and non-allelic genes. The previous explanations of the KPI variability were based on the hypothesis of random mutagenesis as a key factor of KPI polymorphism. RESULTS: KPI-A genes from the genomes of Solanum tuberosum cv. Istrinskii and the wild species Solanum palustre were amplified by PCR with subsequent cloning in plasmids. True KPI sequences were derived from comparison of the cloned copies. "Hot spots" of recombination in KPI genes were independently identified by DnaSP 4.0 and TOPALi v2.5 software. The KPI-A sequence from potato cv. Istrinskii was found to be 100% identical to the gene from Solanum nigrum. This fact illustrates a high degree of similarity of KPI genes in the genus Solanum. Pairwise comparison of KPI A and B genes unambiguously showed a non-uniform extent of polymorphism at different nt positions. Moreover, the occurrence of substitutions was not random along the strand. Taken together, these facts contradict the traditional hypothesis of random mutagenesis as a principal source of KPI gene polymorphism. The experimentally found mosaic structure of KPI genes in both plants studied is consistent with the hypothesis suggesting recombination of ancestral genes. The same mechanism was proposed earlier for other resistance-conferring genes in the nightshade family (Solanaceae). Based on the data obtained, we searched for potential motifs of site-specific binding with plant DNA recombinases. During this work, we analyzed the sequencing data reported by the Potato Genome Sequencing Consortium (PGSC), 2011 and found considerable inconsistence of their data concerning the number, location, and orientation of KPI genes of groups A and B. CONCLUSIONS: The key role of recombination rather than random point mutagenesis in KPI polymorphism was demonstrated for the first time.

Molecular surveillance of clinical Neisseria gonorrhoeae isolates in Russia.

The choice of adequate methods for epidemiological purposes remains a challenging problem in Neisseria gonorrhoeae molecular monitoring. In this study, the collection of geographically unrelated gonococci (n = 103) isolated in Russian clinics was comparably tested by (i) a traditional serotyping scheme, (ii) por typing, (iii) Neisseria gonorrhoeae multiantigen sequence typing (NG-MAST), and (iv) multilocus sequence typing (MLST). It is shown that, according to sequencing data, a third of the strains carried new porB1 alleles, as well as tbpB ones, and more than half of the samples had new sequence types (STs) as determined by NG-MAST or MLST. The discriminatory power for each typing method was calculated by using the Hunter-Gaston discriminatory index, D. Commonly, modern nucleic acid-based typing methods (por typing, NG-MAST, and MLST) appeared to be more efficient than the classical serotyping scheme. While the traditional serotyping gave a D value of 0.82, the por typing, NG-MAST, and MLST approaches yielded D values of 0.97, 0.98, and 0.91, respectively. Each typing technique revealed the distribution of gonococci slightly correlated with their geographical sources. However, only the MLST method STs were highly associated with certain phenotypes. Although ST1594, ST1892, and ST6720 were typical for susceptible gonococci, ST1901 and ST6716 were undoubtedly associated with a multidrug-resistant phenotype. We conclude that every tested nucleic acid-based typing method is suitable for N. gonorrhoeae molecular surveillance. However, the MLST method seems to serve large-scale epidemiological purposes, whereas the NG-MAST and por typing approaches are more appropriate for the investigation of local outbreaks.