Multi-Omics Analysis Revealed the rSNPs Potentially Involved in T2DM Pathogenic Mechanism and Metformin Response.

The goal of our study was to identify and assess the functionally significant SNPs with potentially important roles in the development of type 2 diabetes mellitus (T2DM) and/or their effect on individual response to antihyperglycemic medication with metformin. We applied a bioinformatics approach to identify the regulatory SNPs (rSNPs) associated with allele-asymmetric binding and expression events in our paired ChIP-seq and RNA-seq data for peripheral blood mononuclear cells (PBMCs) of nine healthy individuals. The rSNP outcomes were analyzed using public data from the GWAS (Genome-Wide Association Studies) and Genotype-Tissue Expression (GTEx). The differentially expressed genes (DEGs) between healthy and T2DM individuals (GSE221521), including metformin responders and non-responders (GSE153315), were searched for in GEO RNA-seq data. The DEGs harboring rSNPs were analyzed using the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). We identified 14,796 rSNPs in the promoters of 5132 genes of human PBMCs. We found 4280 rSNPs to associate with both phenotypic traits (GWAS) and expression quantitative trait loci (eQTLs) from GTEx. Between T2DM patients and controls, 3810 rSNPs were detected in the promoters of 1284 DEGs. Based on the protein-protein interaction (PPI) network, we identified 31 upregulated hub genes, including the genes involved in inflammation, obesity, and insulin resistance. The top-ranked 10 enriched KEGG pathways for these hubs included insulin, AMPK, and FoxO signaling pathways. Between metformin responders and non-responders, 367 rSNPs were found in the promoters of 131 DEGs. Genes encoding transcription factors and transcription regulators were the most widely represented group and many were shown to be involved in the T2DM pathogenesis. We have formed a list of human rSNPs that add functional interpretation to the T2DM-association signals identified in GWAS. The results suggest candidate causal regulatory variants for T2DM, with strong enrichment in the pathways related to glucose metabolism, inflammation, and the effects of metformin.

Exploring the Genetic Predisposition to Epigenetic Changes in Alzheimer's Disease.

Alzheimer's disease (AD) is a prevalent type of dementia in elderly populations with a significant genetic component. The accumulating evidence suggests that AD involves a reconfiguration of the epigenetic landscape, including DNA methylation, post-translational modification of histone proteins, and chromatin remodeling. Along with environmental factors, individual specific genetic features play a considerable role in the formation of epigenetic architecture. In this study, we attempt to identify the non-coding regulatory SNPs (rSNPs) able to affect the epigenetic mechanisms in AD. To this end, the multi-omics approach is used. The GEO (Gene Expression Omnibus) available data (GSE153875) for AD patients and controls are integrated to reveal the rSNPs that display allele-specific features in both ChIP-seq profiles of four histone modifications and RNA-seq. Furthermore, we analyze the presence of rSNPs in the promoters of genes reported to be differentially expressed between AD and the normal brain (AD-related genes) and involved in epigenetic regulation according to the EpiFactors database. We also searched for the rSNPs in the promoters of the genes coding for transcription regulators of the identified AD-related genes. These regulators were selected based on the corresponding ChIP-seq peaks (ENCODE) in the promoter regions of these genes. Finally, we formed a panel of rSNPs localized to the promoters of genes that contribute to the epigenetic landscape in AD and, thus, to the genetic predisposition for this disease.

A Panel of rSNPs Demonstrating Allelic Asymmetry in Both ChIP-seq and RNA-seq Data and the Search for Their Phenotypic Outcomes through Analysis of DEGs.

Currently, the detection of the allele asymmetry of gene expression from RNA-seq data or the transcription factor binding from ChIP-seq data is one of the approaches used to identify the functional genetic variants that can affect gene expression (regulatory SNPs or rSNPs). In this study, we searched for rSNPs using the data for human pulmonary arterial endothelial cells (PAECs) available from the Sequence Read Archive (SRA). Allele-asymmetric binding and expression events are analyzed in paired ChIP-seq data for H3K4me3 mark and RNA-seq data obtained for 19 individuals. Two statistical approaches, weighted z-scores and predicted probabilities, were used to improve the efficiency of finding rSNPs. In total, we identified 14,266 rSNPs associated with both allele-specific binding and expression. Among them, 645 rSNPs were associated with GWAS phenotypes; 4746 rSNPs were reported as eQTLs by GTEx, and 11,536 rSNPs were located in 374 candidate transcription factor binding motifs. Additionally, we searched for the rSNPs associated with gene expression using an SRA RNA-seq dataset for 281 clinically annotated human postmortem brain samples and detected eQTLs for 2505 rSNPs. Based on these results, we conducted Gene Ontology (GO), Disease Ontology (DO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses and constructed the protein-protein interaction networks to represent the top-ranked biological processes with a possible contribution to the phenotypic outcome.

Novel approach to functional SNPs discovery from genome-wide data reveals promising variants for colon cancer risk.

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.

Novel functional variants at the GWAS-implicated loci might confer risk to major depressive disorder, bipolar affective disorder and schizophrenia.

BACKGROUND: A challenge of understanding the mechanisms underlying cognition including neurodevelopmental and neuropsychiatric disorders is mainly given by the potential severity of cognitive disorders for the quality of life and their prevalence. However, the field has been focused predominantly on protein coding variation until recently. Given the importance of tightly controlled gene expression for normal brain function, the goal of the study was to assess the functional variation including non-coding variation in human genome that is likely to play an important role in cognitive functions. To this end, we organized and utilized available genome-wide datasets from genomic, transcriptomic and association studies into a comprehensive data corpus. We focused on genomic regions that are enriched in regulatory activity-overlapping transcriptional factor binding regions and repurpose our data collection especially for identification of the regulatory SNPs (rSNPs) that showed associations both with allele-specific binding and allele-specific expression. We matched these rSNPs to the nearby and distant targeted genes and then selected the variants that could implicate the etiology of cognitive disorders according to Genome-Wide Association Studies (GWAS). Next, we use DeSeq 2.0 package to test the differences in the expression of the certain targeted genes between the controls and the patients that were diagnosed bipolar affective disorder and schizophrenia. Finally, we assess the potential biological role for identified drivers of cognition using DAVID and GeneMANIA. RESULTS: As a result, we selected fourteen regulatory SNPs locating within the loci, implicated from GWAS for cognitive disorders with six of the variants unreported previously. Grouping of the targeted genes according to biological functions revealed the involvement of processes such as 'posttranscriptional regulation of gene expression', 'neuron differentiation', 'neuron projection development', 'regulation of cell cycle process' and 'protein catabolic processes'. We identified four rSNP-targeted genes that showed differential expression between patient and control groups depending on brain region: NRAS-in schizophrenia cohort, CDC25B, DDX21 and NUCKS1-in bipolar disorder cohort. CONCLUSIONS: Overall, our findings are likely to provide the keys for unraveling the mechanisms that underlie cognitive functions including major depressive disorder, bipolar disorder and schizophrenia etiopathogenesis.

Genes of susceptibility to early neurodegenerative changes in the rat retina and brain: analysis by means of congenic strains.

BACKGROUND: There has been considerable interest in discovery of the genetic architecture of complex traits, particularly age-related neurodegenerative disorders. To predict disease risk and to understand its genetic basis in humans, it is necessary to study animal models. Our previous research on the accelerated-senescence OXYS strain has revealed two quantitative trait loci (QTLs) on rat chromosome 1 that are associated with early cataract and/or retinopathy as well as with behavioral abnormalities. Each locus was partially mapped within the introgressed segments in a certain congenic strain: WAG/OXYS-1.1 or WAG/OXYS-1.2. Retinal transcriptome profiling of 20-day-old congenic and OXYS rats by high-throughput RNA sequencing uncovered relevant candidate genes and pathways. Nonetheless, the question remained open whether the same genetic components simultaneously have effects on various manifestations of the accelerated-senescence phenotype in OXYS rats. The present study was designed to analyze the genes of susceptibility to early neurodegenerative processes taking place in the OXYS rat retina and brain and to assess their potential functional clustering. The study was based on the findings from recent publications (including mapping of quantitative trait loci) and on comparative phenotyping of congenic rat strains. RESULTS: The backcrossing of Wistar Albino Glaxo (WAG) and OXYS strains to generate the congenics resulted in two congenic strains with high susceptibility to cataract and retinopathy but with no obvious signs of Alzheimer's disease-like brain pathology that are specific for OXYS rats. Thus, the genes of susceptibility to brain neurodegeneration were not introgressed into the congenic strains or there is a strong effect of the genetic background on the disease phenotype. Moreover, the progression of retinopathy with age was relatively less severe in the WAG background compared to the OXYS background. A comparative analysis of previously defined QTLs and congenic segments led to identification of candidate genes with a suspected effect on brain neurodegeneration including the genes showing differential expression in the congenic strains. CONCLUSION: Overall, our findings suggest that the cause of the cataract and the cause of retinopathy phenotypes in OXYS rats may be genetically linked to each other within the introgressed segments in the WAG/OXYS-1.1 and/or WAG/OXYS-1.2 congenic strains.

Application of quantitative trait locus mapping and transcriptomics to studies of the senescence-accelerated phenotype in rats.

BACKGROUND: Etiology 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. RESULTS: First, 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. CONCLUSIONS: This 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.

Association of AMD-like retinopathy development with an Alzheimer's disease metabolic pathway in OXYS rats.

The main cause of vision loss in older individuals is age-related macular degeneration (AMD)--a complex multifactorial disease, whose etiology and pathogenesis are not completely understood. This is due to the impossibility of investigating the early stages of AMD and paucity of biological models. The senescence-accelerated OXYS rats develop retinopathy with clinical and morphological manifestations similar to AMD. But the genetic determinants of its development are not known. Previously we identified quantitative trait loci (QTLs) associated with the development of cataract, retinopathy, and behavioral signs in OXYS rat. In this study, we used bioinformatic analysis to show the enrichment of QTL region with genes associated with neurodegeneration, including a pathway of Alzheimer's disease. For selected list of candidate genes we designed oligonucleotide DNA chips. Using them we found small but significant changes in expression of several genes in OXYS retina compared to disease-free Wistar rats. Among the genes with altered expression were Picalm and Apba2, known to be participants in the processing of the beta-amyloid (Ab). Measurement of Ab 1-42 in the retina showed that its level increases with age in rats, and at advanced stages of retinopathy in OXYS rats, its expression becomes significantly higher than that of disease-free Wistar rats. Based on functional annotation of QTL, microarray, and ELISA results we suggest that accumulation of Ab may have a role in the pathogenesis of retinopathy in OXYS rats.

Identification of functional networks associated with cell death in the retina of OXYS rats during the development of retinopathy.

Age-related macular degeneration (AMD) is a major cause of blindness in developed countries, and the molecular pathogenesis of early events in AMD is poorly understood. Senescence-accelerated OXYS rats develop AMD-like retinopathy. The aim of this study was to explore the differences in retinal gene expression between OXYS and Wistar (control) rats at age 20 d and to identify the pathways of retinal cell death involved in the OXYS retinopathy initiation and progression. Retinal mRNA profiles of 20-day-old OXYS and Wistar rats were generated at the sequencing read depth 40 mln, in triplicate, using Illumina GAIIx. A terminal deoxynucleotidyl transferase-mediated deoxyuridine triphosphate nick end labeling (TUNEL) assay was performed to measure the apoptosis level. GeneMANIA was used to construct interaction networks for differentially expressed (DE) apoptosis-related genes at ages 20 d and 3 and 18 months. Functional analysis was suggestive of a developmental process, signal transduction, and cell differentiation as the most enriched biological processes among 245 DE genes at age 20 d An increased level of apoptosis was observed in OXYS rats at age 20 d but not at advanced stages. We identified functional clusters in the constructed interaction networks and possible hub genes (Rasa1, cFLAR, Birc3, Cdk1, Hspa1b, Erbb3, and Ntf3). We also demonstrated the significance of the extrinsic apoptotic pathway at preclinical, early, and advanced stages of retinopathy development. Besides the cell death signaling pathways, immune system-related processes and lipid-metabolic processes showed overrepresentation in the clusters of all networks. These characteristics of the expression profile of the genes functionally associated with apoptosis may contribute to the pathogenesis of AMD-like retinopathy in senescence-accelerated OXYS rats.

Amyloid accumulation is a late event in sporadic Alzheimer's disease-like pathology in nontransgenic rats.

The amyloid cascade hypothesis posits that deposition of the amyloid ß (Aß) peptide in the brain is a key event in the initiation of Alzheimer's disease (AD). Nonetheless, it now seems increasingly unlikely that amyloid toxicity is the cause of sporadic AD, which leads to cognitive decline. Here, using accelerated-senescence nontransgenic OXYS rats, we confirmed that aggregation of Aß is a later event in AD-like pathology. We showed that an age-dependent increase in the levels of Aß1₋42 and extracellular Aß deposits in the brain of OXYS rats occur later than do synaptic losses, neuronal cell death, mitochondrial structural abnormalities, and hyperphosphorylation of the tau protein. We identified the variants of the genes that are strongly associated with the risk of either late-onset or early-onset AD, including App, Apoe4, Bace1, Psen1, Psen2, and Picalm. We found that in OXYS rats nonsynonymous SNPs were located only in the genes Casp3 and Sorl1. Thus, we present proof that OXYS rats may be a model of sporadic AD. It is possible that multiple age-associated pathological processes may precede the toxic amyloid accumulation, which in turn triggers the final stage of the sporadic form of AD and becomes a hallmark event of the disease.

Senescence-accelerated OXYS rats: a model of age-related cognitive decline with relevance to abnormalities in Alzheimer disease.

Senescence-accelerated OXYS rats are an experimental model of accelerated aging that was established from Wistar stock via selection for susceptibility to cataractogenic effects of a galactose-rich diet and via subsequent inbreeding of highly susceptible rats. Currently, we have the 102nd generation of OXYS rats with spontaneously developing cataract and accelerated senescence syndrome, which means early development of a phenotype similar to human geriatric disorders, including accelerated brain aging. In recent years, our group found strong evidence that OXYS rats are a promising model for studies of the mechanisms of brain aging and neurodegenerative processes similar to those seen in Alzheimer disease (AD). The manifestation of behavioral alterations and learning and memory deficits develop since the fourth week of age, i.e., simultaneously with first signs of neurodegeneration detectable on magnetic resonance imaging and under a light microscope. In addition, impaired long-term potentiation has been demonstrated in OXYS rats by the age of 3 months. With age, neurodegenerative changes in the brain of OXYS rats become amplified. We have shown that this deterioration happens against the background of overproduction of amyloid precursor protein (AßPP), accumulation of ß-amyloid (Aß), and hyperphosphorylation of the tau protein in the hippocampus and cortex. The development of AMD-like retinopathy in OXYS rats is also accompanied by increased accumulation of Aß in the retina. These published data suggest that the OXYS strain may serve as a spontaneous rat model of AD-like pathology and could help to decipher the pathogenesis of AD.

Rat retinal transcriptome: effects of aging and AMD-like retinopathy.

Pathogenesis of age-related macular degeneration (AMD), the leading cause of vision loss in the elderly, remains poorly understood due to the paucity of animal models that fully replicate the human disease. Recently, we showed that senescence-accelerated OXYS rats develop a retinopathy similar to human AMD. To identify alterations in response to normal aging and progression of AMD-like retinopathy, we compared gene expression profiles of retina from 3- and 18-mo-old OXYS and control Wistar rats by means of high-throughput RNA sequencing (RNA-Seq). We identified 160 and 146 age-regulated genes in Wistar and OXYS retinas, respectively. The majority of them are related to the immune system and extracellular matrix turnover. Only 24 age-regulated genes were common for the two strains, suggestive of different rates and mechanisms of aging. Over 600 genes showed significant differences in expression between the two strains. These genes are involved in disease-associated pathways such as immune response, inflammation, apoptosis, Ca ( 2+) homeostasis and oxidative stress. The altered expression for selected genes was confirmed by qRT-PCR analysis. To our knowledge, this study represents the first analysis of retinal transcriptome from young and old rats with biologic replicates generated by RNA-Seq technology. We can conclude that the development of AMD-like retinopathy in OXYS rats is associated with an imbalance in immune and inflammatory responses. Aging alters the expression profile of numerous genes in the retina, and the genetic background of OXYS rats has a profound impact on the development of AMD-like retinopathy.

Quantitative trait loci on chromosome 1 for cataract and AMD-like retinopathy in senescence-accelerated OXYS rats.

Age-related macular degeneration (AMD) and cataract are common age-related diseases in humans. Previously we showed that senescence-accelerated OXYS rats develop retinopathy and cataract, which are comparable to human AMD and senile cataract. Here we focused on the identification of quantitative trait loci (QTLs), which affect early-onset cataract and retinopathy in OXYS rats, using F2 hybrids bred by a reciprocal cross (OXYS×WAG and WAG×OXYS). Chromosome 1 showed significant associations between retinopathy and loci in the regions of markers D1Rat30 and D1Rat219 (QTL1) as well as D1Rat219 and D1Rat81 (QTL2); and between early cataract development with the locus in the region of the markers D1Rat219 and D1Rat81 (QTL2). To determine the effects of these QTLs, we generated two congenic strains by transferring chromosome 1 regions from OXYS into WAG background. Both congenic strains (named WAG/OXYS-1.1 and WAG/OXYS-1.2, respectively) display early cataract and retinopathy development. Thus, we confirmed that genes located in the analyzed regions of chromosome 1 are associated with the development of these diseases in OXYS rats.