BackgroundEarly diagnosis of rheumatoid arthritis (RA) is crucial to providing effective therapy and often hampered by unspecific clinical manifestations. Elevated levels of extracellular circulating DNA (cirDNA) in patients with autoimmune disease were found to be associated with etiopathogenesis. To our knowledge, this is the first study to investigate the putative diagnostic use of cirDNA in RA and its association with disease activity.MethodsBlood samples were taken from 63 healthy subjects (HS) and 74 patients with RA. cirDNA was extracted from plasma and cell surface-bound cirDNA fractions (csbDNA). cirDNA concentration was measured by quantitative real-time polymerase chain reaction. Rheumatoid factor was analyzed by immunonephelometry, whereas C-reactive protein and anticitrullinated protein/peptide antibodies (ACPA) were detected by enzyme-linked immunosorbent assay.ResultsPlasma cirDNA was significantly elevated in patients with RA compared with HS (12.0 versus 8.4 ng/ml, p < 0.01). In contrast, nuclear csbDNA (n-csbDNA) was significantly decreased (24.0 versus 50.8 ng/ml, p < 0.01), whereas mitochondrial csbDNA (m-csbDNA) was elevated (1.44 × 106 copies/ml versus 0.58 × 106 copies/ml, p < 0.05) in RA. The combination of csbDNA (mitochondrial + nuclear) with ACPA reveals the best positive/negative likelihood ratios (LRs) for the discrimination RA from HS (LR+ 61.00, LR− 0.03) in contrast to ACPA (LR+ 9.00, LR− 0.19) or csbDNA (LR+ 8.00, LR− 0.18) alone.ConclusionsNuclear and mitochondrial cirDNA levels in plasma and on the surface of blood cells are modulated in RA. Combination of cirDNA values with ACPA can improve the serological diagnosis of RA.Electronic supplementary materialThe online version of this article (doi:10.1186/s13075-017-1295-z) contains supplementary material, which is available to authorized users.

Chromatin immunoprecipitation followed by deep sequencing (ChIP-Seq) became a method of choice to locate DNA segments bound by different regulatory proteins. ChIP-Seq produces extremely valuable information to study transcriptional regulation. The wet-lab workflow is often supported by downstream computational analysis including construction of models of nucleotide sequences of transcription factor binding sites in DNA, which can be used to detect binding sites in ChIP-Seq data at a single base pair resolution. The most popular TFBS model is represented by positional weight matrix (PWM) with statistically independent positional weights of nucleotides in different columns; such PWMs are constructed from a gapless multiple local alignment of sequences containing experimentally identified TFBSs. Modern high-throughput techniques, including ChIP-Seq, provide enough data for careful training of advanced models containing more parameters than PWM. Yet, many suggested multiparametric models often provide only incremental improvement of TFBS recognition quality comparing to traditional PWMs trained on ChIP-Seq data. We present a novel computational tool, diChIPMunk, that constructs TFBS models as optimal dinucleotide PWMs, thus accounting for correlations between nucleotides neighboring in input sequences. diChIPMunk utilizes many advantages of ChIPMunk, its ancestor algorithm, accounting for ChIP-Seq base coverage profiles ("peak shape") and using the effective subsampling-based core procedure which allows processing of large datasets. We demonstrate that diPWMs constructed by diChIPMunk outperform traditional PWMs constructed by ChIPMunk from the same ChIP-Seq data. Software website: http://autosome.ru/dichipmunk/

No abstract

A vast amount of SNPs derived from genome-wide association studies are represented by non-coding ones, therefore exacerbating the need for effective identification of regulatory SNPs (rSNPs) among them. However, this task remains challenging since the regulatory part of the human genome is annotated much poorly as opposed to coding regions. Here we describe an approach aggregating the whole set of ENCODE ChIP-seq data in order to search for rSNPs, and provide the experimental evidence of its efficiency. Its algorithm is based on the assumption that the enrichment of a genomic region with transcription factor binding loci (ChIP-seq peaks) indicates its regulatory function, and thereby SNPs located in this region are more likely to influence transcription regulation. To ensure that the approach preferably selects functionally meaningful SNPs, we performed enrichment analysis of several human SNP datasets associated with phenotypic manifestations. It was shown that all samples are significantly enriched with SNPs falling into the regions of multiple ChIP-seq peaks as compared with the randomly selected SNPs. For experimental verification, 40 SNPs falling into overlapping regions of at least 7 TF binding loci were selected from OMIM. The effect of SNPs on the binding of the DNA fragments containing them to the nuclear proteins from four human cell lines (HepG2, HeLaS3, HCT-116, and K562) has been tested by EMSA. A radical change in the binding pattern has been observed for 29 SNPs, besides, 6 more SNPs also demonstrated less pronounced changes. Taken together, the results demonstrate the effective way to search for potential rSNPs with the aid of ChIP-seq data provided by ENCODE project.

Fragments of rRNA, mitochondrial transcripts, microRNAs, fragments of scRNAs, snRNA and snoRNA, fragments of several mRNAs as well as the set of newly discovered transcripts were found to be permanent representatives of human blood plasma RNAs. Advanced mapping allowed to identify circulating herpes virus and enterobacterial transcripts. Documented profile of circulating RNA of healthy individuals provides basis for development of new approaches in research and diagnosis of human pathology.

In the majority of colorectal cancer (CRC) cases, the genetic basis of predisposition remains unexplained. The goal of the study was to assess the regulatory SNPs (rSNPs) in the human genome and to reveal СRC drivers based on the available chromatin immunoprecipitation sequencing (ChIP-Seq, ChIA-PET) and transcriptional profiling (RNA-Seq) data. We combined positional (locations within genome regulatory elements) and functional (associated with allele-specific binding and expression) criteria followed by an analysis using genome-wide association studies (GWAS) and minor allele frequency (MAF) datasets. DeSeq2 analysis through 70 CRC patients reinforced the regulatory potential. rSNPs (1,476) that were associated with significant (P < 0.01) allele-specific events resulting in thirty that exhibited a link with CRC according to the MAF and 27, with a risk of malignancy in general according to GWAS. Selected rSNPs may modify the expression of genes for tumor suppressors and the regulators of signaling pathways, including noncoding RNAs. However, the rSNPs from the most represented group affect the expression of genes related to splicing. Our findings strongly suggest that the identified variants might contribute to CRC susceptibility, which indicates that aberrant splicing is one of the key mechanisms for unraveling disease etiopathogenesis and provides useful inputs for interpreting how genotypic variation corresponds to phenotypic outcome.

Chronic stress is a risk factor for major depression. Social defeat stress is a well-validated murine model of depression. However, little is known about the gene activity dynamics during the development of a depression-like state. We analyzed the effects of social defeat stress of varying duration (10 and 30 days) on the behavioral patterns and prefrontal-cortex transcriptome of C57BL/6 mice. The 10-day exposure to social defeat stress resulted in a high level of social avoidance with no signs of depression-associated behavior. Most animals exposed to 30 days of social defeat stress demonstrated clear hallmarks of depression, including a higher level of social avoidance, increased immobility in the forced swimming test, and anhedonic behavior. The monitoring of transcriptome changes revealed widespread alterations in gene expression on the 10th day. Surprisingly, the expression of only a few genes were affected by the 30th day of stress, apparently due to a reversal of the majority of the early stress-induced changes to the original basal state. Moreover, we have found that glucocorticoid-sensitive genes are clearly stimulated targets on the 10th day of stress, but these genes stop responding to the elevated corticosterone level by the 30th day of stress. The majority of genes altered by the 30-day stress were downregulated, with the most relevant ones participating in chromatin modifications and neuroplasticity (e.g., guanine nucleotide exchange factors of the Rho-family of GTPases). Very different molecular responses occur during short-term and long-term social stress in mice. The early-stress response is associated with social avoidance and with upregulation and downregulation of many genes, including those related to signal transduction and cell adhesion pathways. Downregulation of a few genes, in particular, genes for histone-modifying methyltransferases, is a signature response to prolonged stress that induces symptoms of depression. Altogether, our data show that the development of depression under social stress conditions is correlated with suppression of the overactive molecular response to induced stress, involving gene regulatory resistance to glucocorticoid molecules, potentially via a chromatin remodeling mechanism.

BackgroundEtiology of complex disorders, such as cataract and neurodegenerative diseases including age-related macular degeneration (AMD), remains poorly understood due to the paucity of animal models, fully replicating the human disease. Previously, two quantitative trait loci (QTLs) associated with early cataract, AMD-like retinopathy, and some behavioral aberrations in senescence-accelerated OXYS rats were uncovered on chromosome 1 in a cross between OXYS and WAG rats. To confirm the findings, we generated interval-specific congenic strains, WAG/OXYS-1.1 and WAG/OXYS-1.2, carrying OXYS-derived loci of chromosome 1 in the WAG strain. Both congenic strains displayed early cataract and retinopathy but differed clinically from OXYS rats. Here we applied a high-throughput RNA sequencing (RNA-Seq) strategy to facilitate nomination of the candidate genes and functional pathways that may be responsible for these differences and can contribute to the development of the senescence-accelerated phenotype of OXYS rats.ResultsFirst, the size and map position of QTL-derived congenic segments were determined by comparative analysis of coding single-nucleotide polymorphisms (SNPs), which were identified for OXYS, WAG, and congenic retinal RNAs after sequencing. The transferred locus was not what we expected in WAG/OXYS-1.1 rats. In rat retina, 15442 genes were expressed. Coherent sets of differentially expressed genes were identified when we compared RNA-Seq retinal profiles of 20-day-old WAG/OXYS-1.1, WAG/OXYS-1.2, and OXYS rats. The genes most different in the average expression level between the congenic strains included those generally associated with the Wnt, integrin, and TGF-β signaling pathways, widely involved in neurodegenerative processes. Several candidate genes (including Arhgap33, Cebpg, Gtf3c1, Snurf, Tnfaip3, Yme1l1, Cbs, Car9 and Fn1) were found to be either polymorphic in the congenic loci or differentially expressed between the strains. These genes may contribute to the development of cataract and retinopathy.ConclusionsThis study is the first RNA-Seq analysis of the rat retinal transcriptome generated with 40 mln sequencing read depth. The integration of QTL and transcriptomic analyses in our study forms the basis of future research into the relationship between the candidate genes within the congenic regions and specific changes in the retinal transcriptome as possible causal mechanisms that underlie age-associated disorders.