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.

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.

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.

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.

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.

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.

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.

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.

Intergenic, gene terminal, and intragenic CpG islands in the human genome.

BACKGROUND: Recently, it has been discovered that the human genome contains many transcription start sites for non-coding RNA. Regulatory regions related to transcription of this non-coding RNAs are poorly studied. Some of these regulatory regions may be associated with CpG islands located far from transcription start-sites of any protein coding gene. The human genome contains many such CpG islands; however, until now their properties were not systematically studied. RESULTS: We studied CpG islands located in different regions of the human genome using methods of bioinformatics and comparative genomics. We have observed that CpG islands have a preference to overlap with exons, including exons located far from transcription start site, but usually extend well into introns. Synonymous substitution rate of CpG-containing codons becomes substantially reduced in regions where CpG islands overlap with protein-coding exons, even if they are located far downstream from transcription start site. CAGE tag analysis displayed frequent transcription start sites in all CpG islands, including those found far from transcription start sites of protein coding genes. Computational prediction and analysis of published ChIP-chip data revealed that CpG islands contain an increased number of sites recognized by Sp1 protein. CpG islands containing more CAGE tags usually also contain more Sp1 binding sites. This is especially relevant for CpG islands located in 3' gene regions. Various examples of transcription, confirmed by mRNAs or ESTs, but with no evidence of protein coding genes, were found in CAGE-enriched CpG islands located far from transcription start site of any known protein coding gene. CONCLUSIONS: CpG islands located far from transcription start sites of protein coding genes have transcription initiation activity and display Sp1 binding properties. In exons, overlapping with these islands, the synonymous substitution rate of CpG containing codons is decreased. This suggests that these CpG islands are involved in transcription initiation, possibly of some non-coding RNAs.

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.

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.

Activation of the hTERT expression in squamous cell cervical carcinoma is not associated with gene amplification.

The hTERT gene encodes the telomerase catalytic subunit that plays a key role in cancer cell immortalization. Earlier, hTERT amplification was detected in squamous cell cervical carcinomas (SCC), however possible relations between elevated hTERT mRNA level and gene amplification was not studied. Here, we compared the hTERT expression and copy number in the same tumors by quantitative real-time PCR. The hTERT DNA copy number was virtually unchanged in all 33 studied tumors, when compared to normal tissues. This result was confirmed using two reference genes beta-actin and beta-D-glucuronidase. Nevertheless, the activation of hTERT expression was found in 80% of cases (37/46, p<0.001). There was no correlation between the degree of mRNA increase and the tumor size and/or presence of metastases. No hTERT gene expression was observed in 20% of cases (9/46), while the control GADPH expression was unchanged. The detected elevation of the hTERT mRNA level was found using primers specific to functionally active full-length isoform of mRNA. Similar results were obtained with SCC cell lines carrying human papilloma virus (HPV) genomes. We conclude that frequent activation of hTERT expression in SCC is not associated with gene amplification.

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.

Common and specific amino acid residues in the prokaryotic polypeptide release factors RF1 and RF2: possible functional implications.

Termination of protein synthesis is promoted in ribosomes by proper stop codon discrimination by class 1 polypeptide release factors (RFs). A large set of prokaryotic RFs differing in stop codon specificity, RF1 for UAG and UAA, and RF2 for UGA and UAA, was analyzed by means of a recently developed computational method allowing identification of the specificity-determining positions (SDPs) in families composed of proteins with similar but not identical function. Fifteen SDPs were identified within the RF1/2 superdomain II/IV known to be implicated in stop codon decoding. Three of these SDPs had particularly high scores. Five residues invariant for RF1 and RF2 [invariant amino acid residues (IRs)] were spatially clustered with the highest-scoring SDPs that in turn were located in two zones within the SDP/IR area. Zone 1 (domain II) included PxT and SPF motifs identified earlier by others as 'discriminator tripeptides'. We suggest that IRs in this zone take part in the recognition of U, the first base of all stop codons. Zone 2 (domain IV) possessed two SDPs with the highest scores not identified earlier. Presumably, they also take part in stop codon binding and discrimination. Elucidation of potential functional role(s) of the newly identified SDP/IR zones requires further experiments.

Invariant amino acids essential for decoding function of polypeptide release factor eRF1.

In eukaryotic ribosome, the N domain of polypeptide release factor eRF1 is involved in decoding stop signals in mRNAs. However, structure of the decoding site remains obscure. Here, we specifically altered the stop codon recognition pattern of human eRF1 by point mutagenesis of the invariant Glu55 and Tyr125 residues in the N domain. The 3D structure of generated eRF1 mutants was not destabilized as demonstrated by calorimetric measurements and calculated free energy perturbations. In mutants, the UAG response was most profoundly and selectively affected. Surprisingly, Glu55Arg mutant completely retained its release activity. Substitution of the aromatic ring in position 125 reduced response toward all stop codons. This result demonstrates the critical importance of Tyr125 for maintenance of the intact structure of the eRF1 decoding site. The results also suggest that Tyr125 is implicated in recognition of the 3d stop codon position and probably forms an H-bond with Glu55. The data point to a pivotal role played by the YxCxxxF motif (positions 125-131) in purine discrimination of the stop codons. We speculate that eRF1 decoding site is formed by a 3D network of amino acids side chains.