Volume 15 Supplement 12

Selected articles from the IX International Conference on the Bioinformatics of Genome Regulation and Structure\Systems Biology (BGRS\SB-2014): Genomics

Open Access

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

  • Elena E Korbolina1, 2Email author,
  • Nikita I Ershov1,
  • Leonid O Bryzgalov1, 2 and
  • Natalia G Kolosova1, 2
BMC Genomics201415(Suppl 12):S3

https://doi.org/10.1186/1471-2164-15-S12-S3

Published: 19 December 2014

Abstract

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.

Keywords

Senescence-accelerated OXYS rats Age-related macular degeneration Quantitative trait locus Congenic strain RNA-Seq Non-synonymous SNP Retinal transcriptome

Background

The last years have seen great progress in understanding the pathophysiology of complex age-related diseases such as cataract and age-related macular degeneration (AMD): two of the major leading causes of visual impairment and blindness of elderly people in industrial countries. The molecular pathways underlying their onset and progression have yet to be described. Early stages of these diseases cannot be studied in humans, and existing animal models have severe limitations in recapitulating the disease progression [13].

We have previously shown that OXYS rats spontaneously develop a phenotype similar to human aging-associated disorders including retinopathy and cataract, with clinical, morphological, and molecular features similar to human AMD and senile cataract [48]. As we reported recently, neurodegenerative processes observed in OXYS rats are similar to those seen in Alzheimer's disease [911]. This animal model is successfully used to study the pathways and molecular alterations implicated in the development and progression of these disorders as well as to test new therapeutic interventions [1116]. Nevertheless, the genetic architecture of this senescence-accelerated phenotype of OXYS rats remains poorly understood.

Analysis of quantitative trait loci (QTLs) is an unbiased genetic approach to studies of susceptibility genes and molecular pathways involved in diseases of complex etiology with a strong genetic component. The usual paradigm is to produce segregating populations derived from "affected" and "control" strains and to search for linkage of a complex trait to genetic markers using sophisticated statistical techniques. The methods of QTL analysis have long been used successfully in the studies of hypertension [1719], neurodegeneration [20, 21] and modifiers of pathological ocular phenotypes including retinal degeneration in mouse models [2225]. On the other hand, some difficulties arise with the identification of QTL alleles at the level of causative genes [26]. Despite being confirmed and refined in location to a small interval spanning several megabases, a single QTL still harbors a number of genes (or does not generally contain annotated genes) and can act through a confluence of multiple interactions and molecular mechanisms.

Further testing for expression differences is one of the approaches aimed at identifying the genes underlying phenotypic differences [27]. In systems biology, identification of the underlying molecular pathways can be facilitated by high-throughput analyses, particularly at the proteomic and transcriptomic levels. RNA-Seq analysis, based on the next-generation sequencing technology, provides a far more precise quantification of transcripts than do other methods and allows for detection of novel transcriptomic features, e.g., novel exons and alternative splicing variants. The RNA-Seq technology has successfully been used for exploration of complex traits and, particularly, for identifying the genes related to retinal development and diseases [2832].

Our previous QTL studies were focused on chromosome 1 and led to identification of two QTLs associated with early cataract, retinopathy similar to human AMD, and some behavioral aberrations in senescence-accelerated OXYS rats in a cross between OXYS and WAG rats. One QTL (named QTL1) was mapped to the medial region (100.6-188.0 Mbp) and the other QTL (named QTL2) - to the distal region (188.0-250.4 Mbp) of chromosome 1 [33]. Control WAG rats have no appreciable signs of cataract and retinopathy (unpublished data). Two interval-specific congenic strains, carrying OXYS-derived loci of chromosome 1 on the WAG strain, were generated to confirm our findings: WAG/OXYS-1.1 and WAG/OXYS-1.2. Both congenic strains develop early cataract and retinopathy but are clinically different from OXYS rats [33]. According to the results of a preliminary histological study (unpublished), the degenerative changes in the retina of the congenic rats occur in the context of mass migration of mononuclear phagocytes, e.g., into the ganglionic layer. This finding can serve as evidence of the development of an inflammatory process in the retinas of the congenic rats that is typical for patients with AMD, but not for OXYS rats.

The purpose of the present study was to identify candidate genes likely to be involved in the development of early cataract and retinopathy in rats of congenic and OXYS strains by analyzing the RNA-Seq data obtained from retinal RNA samples.

Results

A previous study [33, 34] showed that early cataract and AMD-like retinopathy in senescence-accelerated OXYS rats are associated with two QTLs located on chromosome 1, with the cumulative genomic length of 149.8 Mbp. We generated two interval-specific congenic strains, WAG/OXYS-1.1 and WAG/OXYS-1.2, carrying OXYS-derived chromosome 1 loci in the WAG genetic background. To map QTLs to the congenic loci and to facilitate nomination of the candidate genes responsible for the OXYS phenotype, we examined the retinal transcriptome of 20-day-old WAG/OXYS-1.1, WAG/OXYS-1.2, and OXYS strains using RNA-Seq. One sample of the WAG strain (three 3-month-old animals pooled) was processed, sequenced, and the reads were aligned to the rat genome (Rnor_5.0 assembly) to obtain a characteristic set of single-nucleotide polymorphisms (SNPs) for the second parental strain.

Identifying the congenic regions on chromosome 1

The comparative analysis of SNPs specific to the OXYS, WAG, and congenic strains allowed us to identify the regions of chromosome 1 that came from the OXYS strain and were introgressed into the WAG/OXYS-1.1 and WAG/OXYS-1.2 constructed strains (Figure 1).
Figure 1

The distribution of SNPs specific to the parental strains, on chromosome 1 of the congenic rats: WAG/OXYS-1.1 (panel A) and WAG/OXYS-1.2 (panel C). Panel B shows a plot of distribution of retinopathy linkage probability (LOD score) among the F2 progeny derived from OXYS and WAG strains at age 3 months (using microsatellite markers on chromosome 1). The dotted vertical line at LOD = 2.0 is the threshold for suggestive statistical significance. Red bars correspond to individual SNPs, specific to WAG rats, blue - to OXYS rats. COUNT - a number of coding SNPs for 1 Mbp. The OXYS-derived loci of chromosome 1 transferred to the recipient WAG genome during the construction of the corresponding congenic strains are shown as green (expected) and blue (the congenic segment) inserts, to the left of the linkage map. The blue dotted lines are designed to compare the certain congenic loci between the panels. Note, that the transferred locus was not what we expected in WAG/OXYS-1.1 rats.

We found that the OXYS segments introgressed into the WAG/OXYS-1.1 strain corresponded approximately to the positions of 8.9 × 107-9.7 × 107 bp plus 1.04 × 108-1.05 × 108 bp of chromosome 1. The OXYS segment introgressed into the WAG/OXYS-1.2 strain corresponded approximately to the positions of 1.78 × 108-2.1 × 108 bp plus 2.34 × 108-2.75 × 108 bp of chromosome 1. We found no additional transferred extended loci on other chromosomes in each case (Figure 2). So, the QTL-derived congenic segments were confined to the cumulative genomic segment of approximately 81.1 Mbp. There were 195 and 1013 ENSEMBL IDs within the WAG/OXYS-1.1 and WAG/OXYS-1.2 congenic segments, respectively, according to the results of data mining with the Encode project [35, 36].
Figure 2

The distribution of SNPs, specific to the parental OXYS and WAG strains, on all chromosomes of the congenic rats: WAG/OXYS-1.1 (panel A) and WAG/OXYS-1.2 (panel B). The vertical axis shows counts-the number of coding SNPs for 10 Mbp; the horizontal axis shows the length (in Mbp) of all chromosomes in genome, arranged in numerical order, where the chromosome boundaries are indicated by the red vertical dotted lines. The loci derived from OXYS parental strain are located on chromosome 1 of the congenic strains. Note, that no additional transferred extended loci are found on other chromosomes in each case.

The global transcriptome

Of 26405 genes in the reference genome rno5, our reads were uniquely mapped to 15442 genes (Additional file 1), with at least 10 counts on average. Differential expression of genes among the congenic WAG/OXYS-1.1 and WAG/OXYS-1.2 and parental OXYS strains was assessed using the DESeq package [37], resulting in the lists of 81 and 231 differentially expressed (DE) genes identified in the WAG/OXYS-1.1 versus OXYS and WAG/OXYS-1.2 versus OXYS comparisons, respectively (at padj<0.1 and the cutoff of |log2FC|≥1). We identified only 3 ENSEMBL IDs for genes significantly differentially expressed between the WAG/OXYS-1.1 and WAG/OXYS-1.2 congenic strains (padj<0.1): ENSRNOG00000009373, ENSRNOG00000042565, and ENSRNOG00000001273. Thus, in the WAG/OXYS-1.1 versus WAG/OXYS-1.2 comparison the strategy was to use Gene Ontology (GO) enrichment tools to test for functional categories that are enriched with the genes that are most different in the average expression level. Figure 3 shows the Venn diagram visualizing the overlapping results between the sets of DAVID IDs for DE genes found in the WAG/OXYS-1.1 versus OXYS, WAG/OXYS-1.2 versus OXYS, and in WAG/OXYS-1.1 versus WAG/OXYS-1.2 (speculative observations) comparisons. The DE genes were subjected to various further analyses.
Figure 3

The Venn diagram visualizing the overlapping results between the sets of ENSEMBL IDs for differentially expressed (DE) genes found in the WAG/OXYS-1.1 versus OXYS, WAG/OXYS-1.2 versus OXYS, and in WAG/OXYS-1.1 versus WAG/OXYS-1.2 (speculative observations) comparisons. DESeq used only uniquely mapped to the reference genome reads at the cutoff of |log2FC |≥1.0; padj<0.1 for the WAG/OXYS-1.1 versus OXYS and WAG/OXYS-1.2 versus OXYS DE genes. In the WAG/OXYS-1.1 versus WAG/OXYS-1.2 comparison, the genes most different in the average expression level (pvalue < 0.05; |log2FC |≥1.0) were tested. The number of genes in each group is shown in the corresponding field.

Pathway analysis

To identify pathways and biological functions associated (more than by chance) with the DE genes, we carried out gene annotation enrichment analysis using 2 online bioinformatic tools: DAVID [38] and WebGestalt [39, 40]. The most representative GO terms are shown in the corresponding list of Additional file 2. A complete set of DAVID analyses of the entire lists of DE genes is presented in Additional file 2. Pathway analysis was also performed for separate up- and down-regulated DE groups (Additional file 3).

In the WAG/OXYS-1.1 versus OXYS comparison, DAVID revealed significant enrichment in GO terms relevant to the Srb-mediator complex, peroxisomes, and regulation of neurotransmitter release. Significant enrichment in GO terms, including the MHC class I protein complex, antigen processing and presentation, wound healing, a response to wounding, negative regulation of natural killer cell-mediated cytotoxicity, and Ca2+binding was obtained in other than the WAG/OXYS-1.2 versus OXYS DE gene set (Additional file 2). Some GO terms obtained using these two DE gene sets were relevant. Aside from GO terms related to extracellular matrix components, the analysis of both WAG/OXYS-1.1 versus OXYS and WAG/OXYS-1.2 versus OXYS DE sets of genes revealed categories related to fatty acid metabolic processes, the inflammatory response, and DNA ligation during DNA repair.

DAVID revealed enrichment in GO terms relevant to neuronal differentiation, blood vessel development, melanin biosynthetic processes, cell-cell signaling, membrane rafts, wide pore channel activity, regulation of transmission of nerve impulses, and transcription factor activity in a set of genes most variable in expression between the WAG/OXYS-1.1 and WAG/OXYS-1.2 congenic strains. The enriched KEGG pathway categories unique for this comparison were the Wnt signaling pathway, TGF-β signaling pathway, cadherin signaling pathway, melanogenesis, and Alzheimer's disease-presenilin pathway, where some named categories shared common genes from the five: Bace2, Wnt16, Wnt6, Cdh1, and Cdh3, members of the set of notably variable in average expression genes (Additional file 2, 3).

Cross-referencing of QTL expression differences

The overlap between a gene's nomination via its presence within a QTL-derived congenic segment and its differential expression was analyzed first by cross-referencing of the lists of QTL genes with 81 and 231 DE genes identified in the WAG/OXYS-1.1 versus OXYS and WAG/OXYS-1.2 versus OXYS comparisons, respectively (Additional file 1). Thereby we identified 4 DE genes located only within the second congenic locus (introgressed into the WAG/OXYS-1.2 strain): Zp2, Ifit1, Col17a1, and Snurf.

Polymorphisms

An allelic effect could be brought about at a different level of gene expression and as a result of an altered coding sequence. Further, the comparative analysis of non-synonymous coding SNPs (revealed for OXYS and congenic strains but not for WAG strain) was conducted with the focus only on the genes located within the congenic loci. Analysis using the online Variant Effect Predictor tool [41] revealed that non-synonymous substitution of amino acid residues in the protein products of 6 genes located within the congenic loci in WAG/OXYS-1.1 rats and of 13 genes located within the congenic loci in WAG/OXYS-1.2 rats may affect the function of the proteins (Table 1). The possible impact of an amino acid substitution on the structure and function of the protein was predicted using the SIFT algorithm, scale-invariant feature transform [42]. These 19 new candidate genes were not among the genes differentially expressed between the congenic carrier and OXYS strains (Additional file 1).
Table 1

Genes (within the congenic loci) containing non-synonymous coding SNPs with potential relevance to the senescence-accelerated phenotype in OXYS rats.

Gene

symbol

Official gene name

Allele

CDS position

Protein position

AA change

Codon

change

SIFT prediction (score)

Co-located variation

Congenic locus in WAG/OXYS-1.1 strain

Hps5

Hermansky-Pudlak syndrome 5

C

667

223

R/G

Aga/Gga

tolerated (0.46)

-

Hps5

Hermansky-Pudlak syndrome 5

T

188

63

R/K

aGg/aAg

tolerated (0.25)

-

Cic

capicua transcriptional repressor

A

1517

506

C/Y

tGc/tAc

tolerated (0.08)

-

Arhgap33

RhoGTPase activating protein 33

T

3466

1156

D/N

Gat/Aat

deleterious (0)

-

Lin37

lin-37 homolog (Caenorhabditis elegans)

A

547

183

P/S

Cca/Tca

tolerated (0.18)

-

Cebpg

CCAAT/enhancer binding protein (C/EBP), γ

T

410

137

S/N

aGc/aAc

tolerated (0.45)

-

Nudt19

nudix (nucleoside diphosphate linked moiety X)-type motif 19

C

1015

339

S/G

Agc/Ggc

tolerated (0.42)

-

Congenic locus in WAG/OXYS-1.2 strain

Ccp110

centriolar coiled coil protein 110 kDa

A

868

290

D/N

Gat/Aat

tolerated (0.17)

rs197588212

Gtf3c1

general transcription factor IIIC, polypeptide 1, α

A

6379

2127

R/C

Cgt/Tgt

-

rs8156681

Gtf3c1

general transcription factor IIIC, polypeptide 1, α

C

5692

1898

T/A

Aca/Gca

tolerated (0.75)

rs198509438

Gtf3c1

general transcription factor IIIC, polypeptide 1, α

C

5021

1674

N/S

aAt/aGt

tolerated (0.74)

rs106585674

Rabep2

rabaptin, RAB GTPase binding effector protein 2

C

1007

336

I/T

aTc/aCc

tolerated (1)

rs8153744

Aldoa

aldolase A, fructose-bisphosphate

C

493

165

M/V

Atg/Gtg

tolerated (0.58)

rs8160964

Inpp5f

inositol polyphosphate-5-phosphatase F

C

3075

1025

L/F

ttG/ttC

tolerated (0.44)

-

Psat1

phosphoserine aminotransferase 1

T

316

106

A/T

Gct/Act

tolerated (0.61)

-

Tjp2

tight junction protein 2

T

767

256

R/H

cGc/cAc

tolerated (0.59)

rs198995028

Ak3

adenylate kinase 3

T

368

123

L/Q

cTa/cAa

tolerated (0.66)

rs197195051

Ric1

RAB6A GEF complex partner 1

G

229

77

N/D

Aac/Gac

tolerated (1)

rs197794754

RGD1311595

similar to KIAA2026 protein

C

2613

871

I/M

atA/atG

tolerated (1)

rs198237726

Sfxn3

sideroflexin 3

A

927

309

N/K

aaC/aaA

tolerated (0.92)

rs8163692

Nolc1

nucleolar and coiled-body phosphoprotein 1

A

421

141

E/K

Gag/Aag

tolerated (0.61)

rs8165446

Pdcd11

programmed cell death 11

C

3005

1002

I/T

aTc/aCc

tolerated (1)

-

The SNP quality integral parameter threshold (in Samtools/Bcftools) was Qual> 180, and the consequence according to the Variant Effect Predictor tool, VEP [41] was "missense". Only 1 SNP, located in the Arhgap33 gene, was predicted to significantly affect protein structure (marked as "deleterious" by VEP). Position = position in base pairs (for CDS sequence) or in amino acid residues (for protein); AA = amino acid residue; colocated variation is presented according to the ENCODE data.

Discussion

To identify candidate genes that are potentially involved in the development of early cataract and retinopathy in OXYS and derived congenic rats, we integrated transcriptomic data (obtained from retinas of 2 parental inbred: senescence-accelerated OXYS and WAG strains; and from retinas of 2 derivative congenic strains) with the results of previous QTL mapping of cataract and AMD-like retinopathy on chromosome 1 [33, 34]. The analyses converged on a set of genes, including transcription factors (Table 1, Additional file 1), as the catching candidates (located within any of QTL-derived congenic regions or representing DE candidates), likely to be associated with the senescence-accelerated phenotype in rats.

Congenic loci

In our previous study [33] the identification of genetic loci responsible for the trait of senescence-accelerated phenotype in OXYS rats was performed by QTL-analysis on hybrid F2 population aged 3-4 months and bred by a reciprocal cross between OXYS and WAG rats. Chromosome 1 showed significant associations between retinopathy and two loci: a) the locus marked by microsatellite markers D1Rat30 - D1Rat219 (named QTL1); and b) the locus marked by microsatellite markers D1Rat219 - D1Rat81 (named QTL2) (see Figure 1, panel B). Early cataract development was found significally associated with the locus in the region of the microsatellite markers D1Rat219 - D1Rat81 (QTL2). To determine the effects of these two QTLs, we generated two interval-specific congenic strains by transferring chromosome 1 regions from OXYS on WAG background with the transfer controlled by limiting microsatellite markers. 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 eye diseases in OXYS rats.

In the present study we determined the size and map position of QTL-derived congenic segments by comparative analysis of coding single-nucleotide polymorphisms (SNPs), which were identified for OXYS, WAG, and congenic retinal RNAs by the means of RNA-Seq. We found that a certain congenic locus, introgressed into the WAG/OXYS-1.1 strain, unexpectedly was not the locus in the vicinity of microsatellite markers D1Rat30 and D1Rat219 (Figure 1, panel B), but rather the locus marked by D1Rat224 - D1Rat30, which yielded LOD of 1.38 for linkage with cataract development in OXYS rats according to our unpublished data. Figuring out the reason for this discrepancy is beyond the objectives of the present study because WAG/OXYS-1.1 rats carry the OXYS-derived locus of chromosome 1 and develop cataract and retinopathy. Thus, the influence of genes located in the transferred locus on the disease development in OXYS and congenic WAG/OXYS-1.1 rats can be considered confirmed [43]. A possible explanation is that the congenic strain resulting from selecting only 1 locus of interest will have flanking donor chromosomal segments equal on average to 100/N cM on each side [43, 44], where N is the number of backcrosses and 1 cM is equivalent to 2 × 106 bp of DNA in mice (and presumably in rats). After 8 backcrosses during the production of congenic strains in our study [33], the flanking DNA is on average ~12.5 cM on each side of the selected marker or 25 cM in size.

The locus, transferred from the OXYS to WAG/OXYS-1.2 strain was the expected one when introgressing the locus marked by D1Rat30 and D1Rat219 microsatellites. On other chromosomes, we did not find any extended OXYS-derived loci that were additionally transferred to any congenic strain (see Figure 2).

Here we report that the QTL-derived segments introgressed from OXYS rats are confined to the cumulative genomic length of ~81.1 Mbp and distinctly differentiate the WAG/OXYS-1.1 strain from WAG/OXYS-1.2 strain.

The transcriptome

Because it has been suggested that variation in gene expression is an important contributor to the genetic architecture of complex traits [27], unraveling the changes in transcriptomic signatures should provide further understanding of the susceptible genetic loci that are involved in the pathogenesis of age-associated disorders, including AMD [32, 45, 46]. Recent genetic discoveries demonstrate that the pathogenic pathways in AMD are related to those in Alzheimer's disease among other disorders [47, 48], thereby supporting the notion that AMD is an ocular manifestation of systemic disease processes [49], with a role for inflammation and immune processes in AMD pathogenesis [46, 50, 51].

We reported previously that AMD-like retinopathy is accompanied by downregulation of immune response genes and by changes in mRNA levels of a number of genes related to mitochondrial function in the OXYS retina at ages 3 and 18 months [52]. Impairment of the appropriate inflammatory response in OXYS rats was also observed in other tissues in a model of collagen-induced arthritis (our unpublished data, manuscript in preparation). In addition, we previously reported reduced strength of a delayed hypersensitivity reaction and a decline of T-cell-mediated immunity in OXYS rats [53]. One possible cause of the observed immune imbalance can be attributed to the accelerated thymus involution in OXYS rats [13]. Accordingly, we reported that the development of AMD-like retinopathy in OXYS rats is associated with an imbalance in immune and inflammatory responses [52] and carefully hypothesized that the impairment of immune function in OXYS rats may be associated with systemic neurodegenerative processes. The results of the present study support this hypothesis. When it comes over the GO terms, the imbalance of immune and inflammatory responses is observed at age 20 days between OXYS and both derivative congenic strains. Thus, the GO terms related to negative regulation of leukocyte-mediated cytotoxicity and positive regulation of T-cell-mediated cytotoxicity are downregulated in the WAG/OXYS-1.2 retina compared to OXYS retina; however, the GO term related to MHC class I protein complex is upregulated in this comparison. As for genes from the "Antigen processing and presentation" GO category, these are among up- and down-regulated sets between WAG/OXYS-1.2 and OXYS rats (Additional file 3). The results of pathway analysis are also suggestive of differences in DNA and lipid metabolism between the congenic and OXYS rats.

Candidate genes in the congenic regions

Integration of RNA-Seq data from congenic rats with the results of QTL mapping reduces the number of positional candidates from the genes located throughout two QTL regions [9, 34] to 19 candidate genes within the defined congenic loci. These genes contain non-synonymous coding polymorphisms that were revealed in the cDNA sequence of OXYS and derived congenic rats, but were not detected in the cDNA sequence of WAG rats (Table 1). None of 19 genes were defined as differentially expressed between the congenic carrier and parental OXYS strain (Additional file 1). The identified polymorphisms likely have functional consequences according to web-enabled Variant Effect Predictor Tool, thus may contribute to the senescence-accelerated phenotype of OXYS rats and to the differences in cataract and retinopathy development between congenic and OXYS strains. The consequences for only one of the coding non-synonymous SNPs described here, located in the Arhgap33 gene within the WAG/OXYS-1.1 congenic locus, are expected to significantly affect protein structure (labeled "deleterious" by the Variant Effect Predictor tool). Arhgap33 codes the member of the Rho family of GTPases, which can trigger signal transduction events in the regulation of actin cytoskeleton [54]. Members of this family modulate endocytic trafficking and α-secretase cleavage of the amyloid precursor protein (APP) [55], including interaction with BACE1 [56].

The label for the other missense SNPs in the congenic loci is "tolerated". Nevertheless, both congenic loci contain positional candidate genes that are of interest because of their potential contribution to the disease-associated pathways.

Nudt19, nudix (nucleoside diphosphate-linked moiety X)-type motif 19, is a potential biomarker of early stages of Alzheimer's disease [57] and, along with multiple Nudix family hydrolases, may function in mRNA decapping and maintenance of mRNA stability [58]. The CCAAT/enhancer-binding proteins (C/EBPs) are basic leucine zipper transcription factors that perform important functions in regulation of cell growth and differentiation, the inflammatory response, and hematopoiesis [5962]. Biological functions of the protein product of Cebpg, CCAAT/enhancer-binding protein (C/EBP) γ, are not well understood, but a senescence-suppressing activity of C/EBPγ, which requires its ability to heterodimerize with C/EBPβ, was reported recently [63]. CEBPG is described as a modulator of C/EBP activity [64, 65], a component of transcriptional regulation of B-cell development [62], and as a transcription factor primarily responsible for regulating transcription of key antioxidant and DNA repair genes in normal bronchial epithelial cells [66]. The protein product of Lin37, the lin-37 DREAM MuvB core complex component, is required for regulation of chondrocyte proliferation [67]. Another positional candidate gene within the WAG/OXYS-1.1 congenic locus is Cic, an ortholog of a Drosophila gene that is predominantly expressed during mouse granule cell development [68]. This gene is implicated in c-erbB (Egfr) signaling and is involved in abnormal migration of retinal pigment epithelium (RPE) [69]. As demonstrated on Timp3−/−mice, inhibition of EGFR signaling most likely reduces pathological retinal neovascularization in a mouse model of oxygen-induced retinopathy [70]. The EGFR signaling pathway is considered a possible link between extracellular signaling and regeneration-associated gene expression in the injured retina [71].

In the WAG/OXYS-1.2 congenic locus, we identified 3 non-synonymous coding SNPs within the Gtf3c1 gene, which is the general transcription factor IIIC, polypeptide 1α, coding a subunit of the general transcription factor IIIC [72]. The latter may regulate rearrangement of the nuclear architecture, allowing for coordinated expression of activity-dependent neuronal genes [73].

In addition, 4 candidate genes without coding polymorphisms, identified in the present study, were located in the second congenic locus and showed differential expression (Additional file 1) between the WAG/OXYS-1.2 carrier and parental OXYS strain: Zp2, Ifit1, Col17a1, and Snurf. Note, that this fact suggests a diverse transcriptional regulation which occurs against different genetic background.

The ZP251-149 sperm-binding domain is necessary for human and mouse gamete recognition and penetration through the zona pellucida, the extracellular glycocalyx that surrounds an oocyte. GO annotations related to this gene include "acrosin binding" and "coreceptor activity" [74]. The Col17a1gene encodes the α chain of type XVII collagen, a transmembrane protein that is a structural component of hemidesmosomes, whose mutations result in a blistering disorder non-Herlitz junctional epidermolysis bullosa (JEB-nH) [75]. Ifit1 encodes a protein containing tetratrico peptide repeats that was originally identified as a protein induced upon treatment with interferon; it belongs to the "cytokine signaling in the immune system" category according to REACTOME pathway database [76]. This protein may inhibit viral replication and translation initiation [77, 78].

Snurf, the SNRPN (small nuclear ribonucleoprotein N) gene, was found strongly upregulated in WAG/OXYS-1.2 retina when compared to OXYS (log2FC = -3,01324). This is the best characterized of genes associated with Prader-Willi syndrome (PWS) and Angelman syndrome (AS): neurobehavioral disorders that are caused by a deficiency in imprinted-gene expression from paternal and maternal chromosome 15q11-q13, respectively. The data from human and mouse models point to SNURF's being a protein that is produced along with SmN from a bicistronic transcript (this phenomenon is rarely observed in eukaryotes), with each cistron potentially playing a role in the disease [79, 80]. The highly conserved activator sequence in Snrpn intron 1 demonstrates developmental dynamic changes of DNA methylation and has a methylation-sensitive enhancer activity in rodents; thereby it may control tissue-specific expression of transcripts [81]. Note, that the ubiquitously expressed transcription factor 1, Sp1, interacts with the SNURF protein, according to the Biological General Repository for Interaction Datasets, BioGRID [82]. Sp1 is reported to play a crucial role at early stages of hematopoietic specification [83] and to be involved in photoreceptor-specific transcription either directly or through competition [84] as well as in the regulation of expression of the vascular endothelial growth factor receptor 3 (VEGFR-3) [85].

Candidate DE genes

Recently, we showed that senescence-accelerated OXYS rats are a promising model for studies of the mechanisms of the neurodegenerative processes similar to those seen in Alzheimer's disease [911]. In these rats, behavioral alterations and learning and memory deficits develop by the age of 3 months, accompanied by mitochondrial structural abnormalities, such as formation of cristae-free regions inside mitochondria by the age of 4 months. Moreover, we established a relationship between the retinopathy development in OXYS rats with alterations in the Alzheimer's disease metabolic pathway observed in the retina at the ages of 3 and 18 months [9].

In the present study, we uncovered candidate genes that are already differentially expressed in the retina of congenic rats at the age of 20 days - when compared with each other and with OXYS rats--and are associated with the Alzheimer's disease pathway or with maintenance of mitochondrion structure and function. Some of these genes are mentioned below.

The comparison between WAG/OXYS-1.1 and WAG/OXYS-1.2 strains

Inbred WAG rats used in our study are of normal phenotype. So this comparison was intended to give the candidate DE genes which are associated with the different OXYS-derived loci of chromosome 1 transferred to congenic strains and potentially drive early cataract and retinopathy onset and progression in congenic rats. Simultaneously, we hypothesized that the same DE candidate genes can contribute to the development of the senescence-accelerated phenotype of OXYS rats.

As mentioned above, we identified only 3 ENSEMBL IDs for genes significantly differentially expressed between the WAG/OXYS-1.1 and WAG/OXYS-1.2 congenic strains (padj<0.1): ENSRNOG00000009373, ENSRNOG00000042565, and ENSRNOG00000001273. The last is related to Psmg3, proteasome (prosome, macropain) assembly chaperone 3, which is 1 of 14 distinct gene products essential for biogenesis of mammalian 20S proteasomes [86]. 20S proteasome is a catalytic core of the 26S proteasome, a central enzyme in the degradation of ubiquitin-conjugated proteins. This protein-degradation system might be a potential therapeutic strategy for AD and other tauopathies. The hypothesis that failure of proteasomal and non-proteasomal proteolytic clearance mechanisms leads to tau accumulation and progressive neurofibrillary degeneration in AD has been supported in a number of papers [8789].

However, we could make some speculative assumptions regarding the results of the present study, based on the set of genes that are most variable in expression between the WAG/OXYS-1.1 and WAG/OXYS-1.2 congenic strains (with the cutoff of pvalue < 0.05; |log2FC|≥1). The list of most variable in expression genes from this comparison includes Bace2, β-site APP-cleaving enzyme 2, a paralog of Bace1 (a promising drug target in AD). We can hypothesize that such a difference is due to the differences in molecular networks, underlying the two congenic loci against the WAG background. Although normal functions of both enzymes, Bace1 and Bace2, are still unclear [90], it has been shown that astrocytes possess β-secretase activity and produce Aβ because of the activity of BACE2 [91]. BACE2 variants have an effect on the age of onset of dementia in people with Down syndrome [92]. Among other genes most variable in expression in this comparison are Wnt16 and Wnt6 (upregulated in the WAG/OXYS-1.1 group), the members of wingless-type MMTV integration site family and of the Wnt signaling pathway. The Wnt proteins perform a key function during various stages of retinal development and disease, e.g., they mediate angiogenesis and choroidal neovascularization, a severe complication of AMD [93]. Wnt proteins have also been implicated in AD [94, 95]. The other gene upregulated in the WAG/OXYS-1.1 strain is the Dkk1gene, the product of which influences eye development from a defined developmental time point on, and is essential for lens separation from the surface ectoderm [96]. The function of Dkk1 has not been analyzed in detail in terms of eye development, but most likely the resulting protein acts in a dose-dependent manner, as a competent antagonist of canonical Wnt signaling [97, 98].

NFATC4 (a nuclear factor of activated T-cells, cytoplasmic, calcineurin-dependent 4) is a transcriptional factor involved in the control of the flow of genetic information and in modulation of diverse cellular activities [99]. Nfatc4 exhibits a higher-than-average expression level in the WAG/OXYS-1.1 retina compared to WAG/OXYS-1.2. Accumulating evidence has demonstrated that NFATc4 exerts a proapoptotic effect in multiple diseases, and might participate in retinal ganglion cell apoptosis [100]. The proteoglycan decorin, encoded by the Dcn gene from the group, may also play a role in the differentiation of retinal ganglion cells [101] and is involved in hypoxic retinal damage, significantly reduce in expression in oxygen-induced retinopathy [102, 103].

The comparisons of retinal transcriptome profiles of the congenic and OXYS strains

These comparisons were intended to give the DE genes associated with phenotypic differences between parental OXYS and derivative congenic strains. Candidates from this set are potentially involved in molecular pathways underlying cataract and retinopathy progression, as well as in the development of other features of the senescence-accelerated phenotype in OXYS rats.

Candidate DE genes unique in the WAG/OXYS-1.1 versus OXYS comparison

Ephx2, soluble epoxide hydrolase - a key enzyme in the metabolism and termination of action of epoxyeicosatrienoic acids, is significantly downregulated in WAG/OXYS-1.1 rats when compared to OXYS. Transgenic expression of Ephx2 in the endothelium is known to result in impaired endothelium-dependent vasodilation in the cerebral circulation [104].

Cabp1, a Ca2+-binding protein (CaBP), and the calcium sensor protein caldendrin [105] are alternatively spliced variants of a subfamily of CaBPs with high homology to calmodulin. Caldendrin is abundantly expressed in neurons and is thought to perform an important function at different stages of synaptodendritic Ca2+ signaling [106, 107]. Thus, caldendrin immunoreactivity is displayed by subpopulations of most retinal cell classes, with the exception of glial cells [108]. In the mouse brain, CaBP1/caldendrin is localized both pre- and post-synaptically, where it may regulate Ca2+ signaling and excitability in selected groups of excitatory and inhibitory neurons [109].

The product of Elovl1 gene ortholog, ELOVL fatty acid elongase 1, participates in the metabolism of very long-chain fatty acids (VLCFAs) and a variety of inherited diseases, such as ichthyosis, macular degeneration, myopathy, mental retardation, and demyelination, which are caused by mutations in the genes encoding VLCFA-metabolizing enzymes [110, 111].

Cd14, encoding the CD14 molecule, is significantly downregulated in WAG/OXYS-1.1 rats compared to OXYS strain. It is shown that fibrillar amyloid β (Aβ) can directly interact with Toll-like receptor 2 (TLR2), TLR4, and CD14 to induce microglial Aβ phagocytosis at the initial stages and neuroinflammatory responses at the late stages of Alzheimer's disease [112]. In particular, AD is characterized by the formation of insoluble deposits of Aβ associated with a robust microglia-mediated inflammatory response within the parenchyma of the brain. The involvement of CD14and TLRs in microglial activation was also reported, and it was suggested, that CD14 is a critical regulator of the microglial inflammatory response that acts to modulate Aβ deposition [113, 114].

Is notable, that Mrps10, MRPS10 mitochondrial ribosomal protein S10, from this set of DE genes (upregulated in WAG/OXYS-1.1 rats when compared to OXYS) shows a tendency to be involved in regulation of the lifespan according to reverse transcriptase quantitative PCR (RT-qPCR) analysis using 2 genetic backgrounds [115].

Candidate DE genes unique in the WAG/OXYS-1.2 versus OXYS comparison

The plasma membrane in eukaryotic cells contains microdomains that are enriched in certain glycosphingolipids, gangliosides, and sterols (such as cholesterol) and form membrane/lipid rafts that play a role in several cellular functions or characteristics including polarity of endothelial and epithelial cells, cell migration, mechano transduction, lymphocyte activation, neuronal growth, and signaling in a variety of disease settings [116]. Sdc1, syndecan 1, downregulated in the WAG/OXYS-1.2 strain when compared to OXYS rats, is a highly conserved multifunctional receptor [117]. Syndecan 1 was previously shown to undergo raft-dependent endocytosis upon clustering, regulating atherogenic postprandial remnant lipoproteins and molecules implicated in AD, during uptake of the biologically and medically important ligands, such as HIV-1 [118].

Tnfaip3, tumor necrosis factor α-induced protein 3 from this group, restricts and terminates inflammatory responses through modulation of the ubiquitination status of central components in the NF-κB [119], IRF3 (interferon-regulatory factor 3), and apoptosis signaling cascades. Tnfaip3 expression is necessary for prevention of chronic inflammation and autoimmune pathology according to studies on mice with full or conditional gene deletion [120]. Although Tnfaip3 was ruled out as the candidate gene for the chromosome 10 QTL for age-related retinal degeneration in mice [23], this gene was shown to downregulate adhesion markers, chemokine production, and adventitial angiogenesis, all of which are required for macrophage trafficking to sites of vascular injury in rats [121]. Recent experimental evidence shows that A20/TNFAIP3 is a susceptibility locus for human inflammatory and autoimmune disorders as well as a crucial gatekeeper preserving tissue homeostasis [122].

A human homolog of Yme1l1 has a similarity to mitochondrial conserved ATP-dependent proteinases (termed AAA proteases) that are embedded in the mitochondrial inner membrane, with the catalytic domain facing the mitochondrial intermembrane space [123]. This protein is required for resistance to apoptosis, cristae morphogenesis, and cell proliferation [124]. Nonetheless, how YME1L regulates mammalian mitochondrial function is still not clear.

Hypoxia-inducible factors (HIFs) mediate tissue-specific CBS (cystathionine β-synthase) expression and may augment cerebral vasodilation as an adaptive response to chronic hypoxia [125]. Moreover, ischemia/hypoxia increases accumulation of the CBS proteins in mitochondria [126], being a part of a possible general mechanism of oxygen-sensitive regulation of mitochondrial functions by linking the oxygenation level to accumulation/degradation of mitochondrial heme proteins. The Cbs gene is upregulated in the OXYS retina when compared to WAG/OXYS-1.2 rats.

Carbonic anhydrase IX (CA IX)-deficient (Car9−/−) mice exhibit abnormal locomotor activity, poor performance on a memory test, and vacuolar degenerative changes in the brain; these alterations seem to be age-dependent [127]. Our data show that Car9 expression in OXYS rats is downregulated in comparison with WAG/OXYS-1.2 rats.

The Elf3 codes transcription factor which is highly expressed in the retinal pigment epithelium (RPE) and can regulate important ocular genes, such as TIMP3, in vitro [128]. Variants of Tap2 have been shown to influence the antigenicity of MHC class I molecules by altering the MHC class I ligandome and to lead to reduced negative selection of CD8SP cells in rats [129].

Candidate DE genes common for the WAG/OXYS-1.1 versus OXYS and WAG/OXYS-1.2 versus OXYS comparisons

Mediator (MED) complex is a multiprotein assembly playing a key role in eukaryotic transcription. Alterations of MED function are implicated in human diseases, e.g., alternative splicing of any MED gene may contribute to neovascularization [130]. Although neovascularization characterizes the end-stage of "dry" and "wet" forms of AMD, MED8 was found upregulated and MED20 - downregulated in both congenic strains compared to OXYS rats at the age of 20 days.

Two DE genes from this group are associated with maintenance of mitochondrion function: Ehhadh and Acadm. Ehhadh, L-bifunctional enzyme, is indispensable for the production of medium-chain dicarboxylic acids, thus providing an explanation for the coordinated induction of mitochondrial and peroxisomal oxidative pathways during fasting [131]. A synonymous polymorphic variation of ACADM, acyl-CoA dehydrogenase, C-4 to C-12 straight chain (exon 11) affects splicing efficiency and may alter fatty acid oxidation [132].

Candidate DE genes in the Alzheimer's disease network

The DE genes found in WAG/OXYS-1.1 versus OXYS and WAG/OXYS-1.2 versus OXYS comparisons were also compared with human genes and their first-order interacting partners of the Alzheimer's disease network annotated in the NetAge database, which is used for the study of aging, longevity and age-related diseases [133]. Among several genes, we have got a potentially relevant candidate: fibronectin, Fn1, in the overlap. This gene encodes a glycoprotein, a member of a common signaling network which comprises over 40% of all proteins with multiple interactions in the human interactome [134]. Fibronectin binds to a complement protein C1q [135, 136], which is induced in the brain in response to a variety of neuronal injuries, including AD, and blocks fibrillar amyloid-β (fAβ) neurotoxicity [137]. C1q is also well-known risk factor for retinal diseases related to oxidative stress, inflammation, and the complement cascade, e.g. is associated with complement alternative pathway of RPE cell death in AMD [138, 139].

Summarizing our results, the QTL-derived loci of chromosome 1 in both congenic strains contain strong polymorphic candidates. Moreover, promising DE candidates are associated with each OXYS-derived locus on WAG/OXYS-1.1 and WAG/OXYS-1.2 genetic background. It appears that the trait of accelerated senescence under study is truly polygenic, with more than one gene contributing to OXYS-derived loci in each congenic strain even when some of these genes turn out to be not causative. As both congenic strains displayed early cataract and retinopathy but differed clinically from OXYS rats, it can be argued that the senescence-accelerated phenotype in OXYS rats is associated not just with the congenic loci under study. Thorough validation studies are necessary to address this question.

In general, the studies of variation in transcriptome snapshots in various tissues under different influences in congenic and parental strains can be used to explore and test molecular networks underlying variation in pathological phenotypes. It is possible that the systemic factors, such as environmental [140142], have an influence on cataract and retinopathy development and contribute to differential expression of genes between strains. For some complex diseases, viral, bacterial, or environmental causes may ultimately prove to be of greater importance than disease-associated genes [141, 143145]. So, the combination of genetic factors with known environmental determinants indicates the highly complex nature of disease, but at the same time also offers insights into risk prediction [146]. Finally, some of the genes described here as candidates might be involved in the characteristic features of the senescence-accelerated phenotype that are not manifested during eye disease progression (e.g., behavioral variation observed in standard tests). Some of candidate genes may become expressed with age during development of a pathological phenotype. Thorough validation studies are needed to address this question, yet we can speculate that secondary changes in gene expression patterns between congenic and parental OXYS strains are triggered in the network associated with QTL-causing genes. In the present study, genes encoding transcriptional regulators that are polymorphic or differentially expressed in association with the congenic loci are especially good candidates.

We also identified significant transcript differences with potential relevance to disease development in a number of candidate genes, but we have not found SNPs within the coding regions of these DE genes. It is likely that other factors, e.g., the presence of non-coding SNPs in elements that regulate gene expression (including promoter and enhancer regions and the sites of transcription factor's binding) contribute to the observed transcriptional differences. Their detection requires further integrative studies.

Identification of master regulators of biological processes, including disease susceptibility and mapping of downstream gene networks, remains a big challenge for systems biology. Therefore, the efforts of researchers are focused on integrating results from diverse experiments, especially the data generated using high-throughput methods: genome-wide association (GWA) studies, chromatin immunoprecipitation followed by sequencing (ChIP-Seq), and RNA-Seq [147149]. These approaches provide more complex information on expression levels, differential splicing, allele-specific expression, RNA editing, and fusion transcripts resulting from chromosomal translocations and other mutations [150153], and also are helpful to investigate protein-DNA interactions [154]. As for the rat model, transcriptomic changes across multiple organs or developmental stages are recently beginning to be reported and integrated into open-access databases [149, 155, 156]. To our knowledge, this study represents the first analysis of the rat retinal transcriptome generated with 40 mln sequencing read depth by RNA-Seq technology.

Conclusion

In this study, we carried out a comprehensive analysis of RNA-Seq data obtained for retinas of rats from parental and derivative congenic strains differing in the development of cataract and retinopathy. Furthermore, the genetic background of the trait under study obtained in previously published QTL-analysis was determined more accurately, and we overlaid the transcriptomic data obtained for 20-day-old rat retinas on the results. Thus, we obtained an unbiased list of candidate genes that may affect disease-associated pathways. In summary, our research of the congenic strains may facilitate dissecting of the complex genetic architecture of the senescence-accelerated phenotype in OXYS rats, although there is a need for further characterization of the congenic rats. With an improved understanding of the underlying genetic susceptibility, researchers can hope to identify therapeutic targets to nip AMD in the bud and to prevent its progression and vision loss in humans.

Materials and methods

Animals

All animal procedures were in compliance with the Association for Research in Vision and Ophthalmology statement for the Use of Animals in Ophthalmic and Vision Research and the European Communities Council Directive 86/609/EES. All manipulations of rats were approved by the Institutional Review Board #9 of the Institute of Cytology and Genetics, the Siberian Branch of the Russian Academy of Sciences (SB RAS), according to The Guidelines for Manipulations of Experimental Animals. Twenty-day-old male senescence-accelerated OXYS (n = 3), age-matched male WAG/OXYS-1.2 (n = 3) and WAG/OXYS-1.1 (n = 3) congenic rats, and 3-month-old WAG rats (n = 3) were obtained from the Breeding Experimental Animal Laboratory of the Institute of Cytology and Genetics, SB RAS (Novosibirsk, Russia). The OXYS rat strain was derived from the Wistar strain at the Institute of Cytology and Genetics as described earlier [157] and registered in the Rat Genome Database [158]. At this point, we have the 105th generation of OXYS rats who exhibit spontaneously developing cataract and accelerated senescence syndrome inherited in a linked manner. The congenic rats were originally created as described earlier [33, 34]. At this point, we have the 13th generation of the congenic rats who exhibit spontaneously developing cataract and retinopathy. At age 4 weeks, the pups were weaned, housed in groups of 5 animals per cage (57 × 36 × 20 cm), and kept under standard laboratory conditions (22°C ± 2°C, 60% relative humidity, and 12 h light/12 h dark cycle). The animals were provided with standard rodent feed (PK-120-1, Ltd, Laboratorsnab, Russia) and water ad libitum. All experimental procedures were in compliance with the European Communities Council Directive of 24 November 1986 (86/609/EEC). All efforts were made to minimize the number of animals used and their discomfort.

Tissue preparation

All rats were euthanized by CO2 asphyxiation. After decapitation the chorioretinal complex was excised rapidly on ice, placed in RNALater reagent (cat. #AM7020, Ambion), incubated at 4°C overnight in accordance with the manufacturer's protocol, frozen, and stored at −70°C prior to further processing.

Sampling

The chorioretinal complex samples from three 20-day-old OXYS and three 20-day-old WAG/OXYS-1.2 rats were collected, processed, and sequenced individually, in triplicate. For 20-day-old WAG/OXYS-1.1 rats and 3-month-old WAG rats, 3 chorioretinal complex samples from 3 animals per strain were pooled to reduce the cost of the experiment. The RNA-Seq data obtained for WAG strain were not used for comparative analysis of retinal transcriptional profiles, as WAG rats were of other age group.

RNA isolation

Frozen rat tissues were lysed in the TRIzol Reagent (cat. #15596-018, Invitrogen), and total RNA was isolated according to the manufacturer's protocols. RNA quality and quantity were assessed using Agilent Bioanalyser (Agilent) and the resulting RNA concentrations were determined using a Nanodrop 8000 spectrophotometer. The samples were stored at −70°C until further processing.

Illumina sequencing

More than 40 mln single-end reads 50 bp long were obtained for each sample of retinal RNAs, using Illumina non-stranded sequencing on an Illumina GA IIx instrument at the "Genoanalitika" Lab, Moscow [159] in accordance with standard Illumina protocols (mRNA-Seq Sample Prep Kit, cat. #1004816). Briefly, polyA-tailed mRNA was purified from total RNA using Sera-Mag Magnetic Oligo (dT) beads and then fragmented into small pieces using divalent cations and heating. Using reverse transcriptase and random primers, first- and second-strand cDNAs were synthesized. The cDNA was processed in an end repair reaction with T4 DNA polymerase and Klenow DNA polymerase in order to blunt the termini. An "A" base was then added to the 3′ end of the blunt phosphorylated DNA fragments, and an Illumina adaptor with a single T overhang at its 3′ end was then ligated to the end of the DNA fragment, for hybridization in a single-read flow cell. After that, a size range of cDNA templates was selected, and these fragments were amplified on a cluster station using Single-Read Cluster Generation Kit v2. Sequencing-by-synthesis (SBS) of 50-nucleotide length was performed using SBS v4 reagents on a Genome Analyzer IIx running the SCS2.8 software (Illumina, cat. #FC-940-4001).

Mapping and SNP discovery

After barcode trimming, the sequencing data were assessed for quality using the FastQC software and were mapped to the Rattus norvegicus reference genome assembly RGSC 5.0 (Ensemble release 75) using Bowtie 2 or TopHat v2.0.4. The SNP positions within the aligned reads (compared to the reference genome) were identified using the pileup function in SAMtools (v. 0.1.17) utilities [160]. Using the various filter commands, we predicted SNPs for various positions using the minimal mapping quality (Qual) of 100. These parameters ensured high-quality reliable mapping of the reads, which is important for variant calling. Using custom-made Perl scripts, we converted the VCF files into MySQL tables. Further analysis was based on the comparison of SNPs specifically present in OXYS, WAG/OXYS-1.1, WAG/OXYS-1.2, and WAG rats using MySQL 5.0 query, Oracle Corporation [161]. The free programming language R was used for the final data processing and visualization [162].

Prediction of the SNP phenotypes

The SNPs--present in the specific congenic locus of chromosome 1 in OXYS and congenic rats (WAG/OXYS-1.1 or WAG/OXYS-1.2) but absent in control WAG rats--were considered possibly associated with the senescence-accelerated phenotype. The effect of a variant (amino acid substitution) on protein function was predicted using the Variant Effect Predictor web tool [41]; the consequence type, SIFT score, and prediction were calculated for each variant. SIFT scores of 0-0.05 were classified as "deleterious" and 0.05-1 - as "tolerated."

Gene expression analysis

Sequencing data were preprocessed with Cutadapt [163] tool in order to remove adapters and low-quality sequences. The resulting reads were mapped on Rnor_5.0 reference genome assembly using TopHat2 [164]. The data were then converted into gene count tables using ENSEMBL gene annotaton data. The resulting tables were subjected to the analysis of differential gene expression with DESeq package [37]. The Benjamini-Hochberg correction for multiple testing was applied to the resulting p-values, and the genes with adjusted p-value < 0.1 were selected as differentially expressed for further study. Genes that failed to converge to Generalized Linear Model (GLM) in DE analysis were excluded.

Pathway analysis

A bioinformatic approach was used to determine the biological context of the large amounts of gene expression (RNA-Seq) data. Gene lists from comparative analyses showing significant differences in gene expression were submitted to the free Database for Annotation, Visualization, and Integrated Discovery, DAVID [38]. In addition, functional annotation was conducted using the WEB-based Gene SeT AnaLysis Toolkit, WebGestalt [39, 40]. All rat genes were set as a background. The NetAge database [133] was surfed to compare the rat DE genes from this study with human orthologs and their partners from the annotated Alzheimer's disease network.

Availability of supporting data

The data sets supporting the results of this article are included within the article and its additional files.

List of abbreviations

Mln: 

million

Mb: 

megabase (106 base pairs)

cM: 

Centimorgan

AMD: 

age-related macular degeneration

RPE: 

retinal pigment epithelium

QTL: 

quantitative trait locus

LOD: 

LOD score: log10 likelihood ratio, comparing single-QTL model to the "no QTL anywhere" model

DE: 

differential expression

DEGs: 

DE genes: differentially expressed genes

NGS: 

next generation sequencing

GO: 

gene ontology

KEGG: 

Kyoto Encyclopedia of Genes and Genomes

SNP: 

single-nucleotide polymorphism

Declarations

Acknowledgements

This work was supported by the Ministry of Education and Science of the Russian Federation and partially (when regards to instrumental high-throughput sequencing) by Grants of the Government of the Russian Federation N 2012-220-03-435 and N 14.B25.31.0033. The computations were performed at the Siberian Supercomputer Center SB RAS (Novosibirsk, Russia). We wish to thank the Breeding Experimental Animal Laboratory of the Institute of Cytology and Genetics, SB RAS (Novosibirsk, Russia) for their care and maintenance of the experimental animals used in this study. The study was conducted in the Center for Genetic Resources of Laboratory Animals ICG SB RAS (RFMEFI61914X0005).

Declarations

Publication of this article has been funded by the BGRS\SB-2014 Organizing Committee and the Russian Science Foundation (RSF).

This article has been published as part of BMC Genomics Volume 15 Supplement 12, 2014: Selected articles from the IX International Conference on the Bioinformatics of Genome Regulation and Structure\Systems Biology (BGRS\SB-2014): Genomics. The full contents of the supplement are available online at http://www.biomedcentral.com/bmcgenomics/supplements/15/S12.

Authors’ Affiliations

(1)
Institute of Cytology and Genetics, SB RAS
(2)
Novosibirsk State University

References

  1. Elizabeth Rakoczy P, Yu MJT, Nusinowitz S, Chang B, Heckenlively JR: Mouse models of age-related macular degeneration. Exp Eye Res. 2006, 82: 741-52. 10.1016/j.exer.2005.10.012.PubMedGoogle Scholar
  2. Dubrovsky Y V, Samsa WE, Kondratov R V: Deficiency of circadian protein CLOCK reduces lifespan and increases age-related cataract development in mice. Aging (Albany NY). 2010, 2: 936-44.Google Scholar
  3. Strunnikova N V, Maminishkis A, Barb JJ, Wang F, Zhi C, Sergeev Y, Chen W, Edwards AO, Stambolian D, Abecasis G, Swaroop A, Munson PJ, Miller SS: Transcriptome analysis and molecular signature of human retinal pigment epithelium. Hum Mol Genet. 2010, 19: 2468-86. 10.1093/hmg/ddq129.PubMedPubMed CentralGoogle Scholar
  4. Kolosova NG, Lebedev PA, Aidagulova S V, Morozkova TS: OXYS rats as a model of senile cataract. Bull Exp Biol Med. 2003, 136: 415-9.PubMedGoogle Scholar
  5. Kolosova NG, Lebedev PA, Fursova AZ, Moroskova TS, Gusarevich OG: [Prematurely aging OXYS rats as an animal model of senile cataract in human]. Adv Gerontol = Uspekhi Gerontol/Ross Akad Nauk Gerontol Obs. 2003, 12: 143-8.Google Scholar
  6. Markovets AM, Saprunova VB, Zhdankina AA, Fursova AZ, Bakeeva LE, Kolosova NG: Alterations of retinal pigment epithelium cause AMD-like retinopathy in senescence-accelerated OXYS rats. Aging (Albany NY). 2011, 3: 44-54.Google Scholar
  7. Rumyantseva Y V, Fursova AZ, Fedoseeva LA, Kolosova NG: Changes in physicochemical parameters and alpha-crystallin expression in the lens during cataract development in OXYS rats. 5. Biochem Biokhimiiâ. 2008, 73: 1176-82. 10.1134/S0006297908110023.Google Scholar
  8. Zhdankina AA, Fursova AZ, Logvinov S V, Kolosova NG: Clinical and morphological characteristics of chorioretinal degeneration in early aging OXYS rats. Bull Exp Biol Med. 2008, 146: 455-8. 10.1007/s10517-009-0298-4.PubMedGoogle Scholar
  9. Kozhevnikova OS, Korbolina EE, Stefanova NA, Muraleva NA, Orlov YL, Kolosova NG: Association of AMD-like retinopathy development with an Alzheimer's disease metabolic pathway in OXYS rats. Biogerontology. 2013, 14: 753-62. 10.1007/s10522-013-9439-2.PubMedGoogle Scholar
  10. Stefanova NA, Kozhevnikova OS, Vitovtov AO, Maksimova KY, Logvinov S V, Rudnitskaya EA, Korbolina EE, Muraleva NA, Kolosova NG: Senescence-accelerated OXYS rats: a model of age-related cognitive decline with relevance to abnormalities in Alzheimer disease. Cell Cycle. 2014, 13: 898-909. 10.4161/cc.28255.PubMedPubMed CentralGoogle Scholar
  11. Stefanova NA, Muraleva NA, Skulachev VP, Kolosova NG: Alzheimer's disease-like pathology in senescence-accelerated OXYS rats can be partially retarded with mitochondria-targeted antioxidant SkQ1. J Alzheimers Dis. 2014, 38: 681-94.PubMedGoogle Scholar
  12. Neroev V V, Archipova MM, Bakeeva LE, Fursova AZ, Grigorian EN, Grishanova AY, Iomdina EN, Ivashchenko ZN, Katargina LA, Khoroshilova-Maslova IP, Kilina O V, Kolosova NG, Kopenkin EP, Korshunov SS, Kovaleva NA, Novikova YP, Philippov PP, Pilipenko DI, Robustova O V, Saprunova VB, Senin II, Skulachev M V, Sotnikova LF, Stefanova NA, Tikhomirova NK, Tsapenko I V, Shchipanova AI, Zinovkin RA, Skulachev VP: Mitochondria-targeted plastoquinone derivatives as tools to interrupt execution of the aging program. 4. Age-related eye disease. SkQ1 returns vision to blind animals. Biochem Biokhimiiâ. 2008, 73: 1317-28. 10.1134/S0006297908120043.Google Scholar
  13. Obukhova LA, Skulachev VP, Kolosova NG: Mitochondria-targeted antioxidant SkQ1 inhibits age-dependent involution of the thymus in normal and senescence-prone rats. Aging (Albany NY). 2009, 1: 389-401.Google Scholar
  14. Stefanova NA, Fursova AZ, Sarsenbaev KN, Kolosova NG: Effects of Cistanche deserticola on behavior and signs of cataract and retinopathy in senescence-accelerated OXYS rats. J Ethnopharmacol. 2011, 138: 624-32. 10.1016/j.jep.2011.10.017.PubMedGoogle Scholar
  15. Perepechaeva ML, Grishanova AY, Rudnitskaya EA, Kolosova NG: The Mitochondria-Targeted Antioxidant SkQ1 Downregulates Aryl Hydrocarbon Receptor-Dependent Genes in the Retina of OXYS Rats with AMD-Like Retinopathy. J Ophthalmol. 2014, 2014: 530943-PubMedPubMed CentralGoogle Scholar
  16. Maslov MY, Chernysheva GA, Smol'jakova VI, Aliev OI, Kolosova NG, Plotnikov MB: Hemorheological parameters and their correlations in OXYS rats: A new model of hyperviscosity syndrome. Clin Hemorheol Microcirc. 2014Google Scholar
  17. Simino J, Shi G, Kume R, Schwander K, Province MA, Gu CC, Kardia S, Chakravarti A, Ehret G, Olshen RA, Turner ST, Ho L-T, Zhu X, Jaquish C, Paltoo D, Cooper RS, Weder A, Curb JD, Boerwinkle E, Hunt SC, Rao DC: Five blood pressure loci identified by an updated genome-wide linkage scan: meta-analysis of the Family Blood Pressure Program. Am J Hypertens. 2011, 24: 347-54. 10.1038/ajh.2010.238.PubMedPubMed CentralGoogle Scholar
  18. Chauvet C, Ménard A, Xiao C, Aguila B, Blain M, Roy J, Deng AY: Novel genes as primary triggers for polygenic hypertension. J Hypertens. 2012, 30: 81-6. 10.1097/HJH.0b013e32834dddb1.PubMedGoogle Scholar
  19. Buys ES, Raher MJ, Kirby A, Shahid M, Mohd S, Baron DM, Hayton SR, Tainsh LT, Sips PY, Rauwerdink KM, Yan Q, Tainsh RET, Shakartzi HR, Stevens C, Decaluwé K, Rodrigues-Machado M da G, Malhotra R, Van de Voorde J, Wang T, Brouckaert P, Daly MJ, Bloch KD: Genetic modifiers of hypertension in soluble guanylate cyclase α1-deficient mice. J Clin Invest. 2012, 122: 2316-25. 10.1172/JCI60119.PubMedPubMed CentralGoogle Scholar
  20. Swanberg M, Harnesk K, Ström M, Diez M, Lidman O, Piehl F: Fine mapping of gene regions regulating neurodegeneration. PLoS One. 2009, 4: e5906-10.1371/journal.pone.0005906.PubMedPubMed CentralGoogle Scholar
  21. Jiao Y, Lu L, Williams RW, Smeyne RJ: Genetic dissection of strain dependent paraquat-induced neurodegeneration in the substantia nigra pars compacta. PLoS One. 2012, 7: e29447-10.1371/journal.pone.0029447.PubMedPubMed CentralGoogle Scholar
  22. Won J, Charette JR, Philip VM, Stearns TM, Zhang W, Naggert JK, Krebs MP, Nishina PM: Genetic modifier loci of mouse Mfrp(rd6) identified by quantitative trait locus analysis. Exp Eye Res. 2014, 118: 30-5.PubMedPubMed CentralGoogle Scholar
  23. Ogando DG, Dahlquist KD, Alizadeh M, Kunchithapautham K, Li J, Yu N, LaVail MM, Rohrer B, Vollrath D, Danciger M: Candidate genes for chromosomes 6 and 10 quantitative trait loci for age-related retinal degeneration in mice. Mol Vis. 2010, 16: 1004-18.PubMedPubMed CentralGoogle Scholar
  24. Danciger M, Ogando D, Yang H, Matthes MT, Yu N, Ahern K, Yasumura D, Williams RW, Lavail MM: Genetic modifiers of retinal degeneration in the rd3 mouse. Invest Ophthalmol Vis Sci. 2008, 49: 2863-9. 10.1167/iovs.08-1715.PubMedPubMed CentralGoogle Scholar
  25. Wolf N, Galecki A, Lipman R, Chen S, Smith-Wheelock M, Burke D, Miller R: Quantitative trait locus mapping for age-related cataract severity and synechia prevalence using four-way cross mice. Invest Ophthalmol Vis Sci. 2004, 45: 1922-9. 10.1167/iovs.03-0435.PubMedGoogle Scholar
  26. Flint J, Valdar W, Shifman S, Mott R: Strategies for mapping and cloning quantitative trait genes in rodents. Nat Rev Genet. 2005, 6: 271-86. 10.1038/nrg1576.PubMedGoogle Scholar
  27. Mackay TFC, Stone EA, Ayroles JF: The genetics of quantitative traits: challenges and prospects. Nat Rev Genet. 2009, 10: 565-77. 10.1038/nrg2612.PubMedGoogle Scholar
  28. Ratnapriya R, Swaroop A: Genetic architecture of retinal and macular degenerative diseases: the promise and challenges of next-generation sequencing. Genome Med. 2013, 5: 84-10.1186/gm488.PubMedPubMed CentralGoogle Scholar
  29. Farkas MH, Grant GR, White JA, Sousa ME, Consugar MB, Pierce EA: Transcriptome analyses of the human retina identify unprecedented transcript diversity and 3.5 Mb of novel transcribed sequence via significant alternative splicing and novel genes. BMC Genomics. 2013, 14: 486-10.1186/1471-2164-14-486.PubMedPubMed CentralGoogle Scholar
  30. Gamsiz ED, Ouyang Q, Schmidt M, Nagpal S, Morrow EM: Genome-wide transcriptome analysis in murine neural retina using high-throughput RNA sequencing. Genomics. 2012, 99: 44-51. 10.1016/j.ygeno.2011.09.003.PubMedPubMed CentralGoogle Scholar
  31. Kandpal RP, Rajasimha HK, Brooks MJ, Nellissery J, Wan J, Qian J, Kern TS, Swaroop A: Transcriptome analysis using next generation sequencing reveals molecular signatures of diabetic retinopathy and efficacy of candidate drugs. Mol Vis. 2012, 18: 1123-46.PubMedPubMed CentralGoogle Scholar
  32. Cai L, Lyu YL: Analysis of Retinal Development and Diseases Using RNA-Seq. Cell Dev Biol. 2012, 1:Google Scholar
  33. Korbolina EE, Kozhevnikova OS, Stefanova NA, Kolosova NG: Quantitative trait loci on chromosome 1 for cataract and AMD-like retinopathy in senescence-accelerated OXYS rats. Aging (Albany NY). 2012, 4: 49-59.Google Scholar
  34. Oĭdopova OS, Polygalova NE, Korbolina EE, Kolosova NG: [A search for genetic determinants of premature aging in OXYS rats]. Adv Gerontol = Uspekhi Gerontol/Ross Akad Nauk Gerontol Obs. 2008, 21: 499-500.Google Scholar
  35. A user's guide to the encyclopedia of DNA elements (ENCODE). PLoS Biol. 2011, 9: e1001046-10.1371/journal.pbio.1001046.Google Scholar
  36. The ENCODE (ENCyclopedia Of DNA Elements) Project. Science. 2004, 306: 636-40.Google Scholar
  37. Anders S, Huber W: Differential expression analysis for sequence count data. Genome Biol. 2010, 11: R106-10.1186/gb-2010-11-10-r106.PubMedPubMed CentralGoogle Scholar
  38. Dennis G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA: DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003, 4: P3-10.1186/gb-2003-4-5-p3.PubMedGoogle Scholar
  39. Zhang B, Kirov S, Snoddy J: WebGestalt: an integrated system for exploring gene sets in various biological contexts. Nucleic Acids Res. 2005, 33: (Web Server issue):W741-8. 10.1093/nar/gki475.PubMedGoogle Scholar
  40. Wang J, Duncan D, Shi Z, Zhang B: WEB-based GEne SeT AnaLysis Toolkit (WebGestalt): update 2013. Nucleic Acids Res. 2013, W77-83. 41 Web ServerGoogle Scholar
  41. McLaren W, Pritchard B, Rios D, Chen Y, Flicek P, Cunningham F: Deriving the consequences of genomic variants with the Ensembl API and SNP Effect Predictor. Bioinformatics. 2010, 26: 2069-70. 10.1093/bioinformatics/btq330.PubMedPubMed CentralGoogle Scholar
  42. Ng PC, Henikoff S: Predicting deleterious amino acid substitutions. Genome Res. 2001, 11: 863-74. 10.1101/gr.176601.PubMedPubMed CentralGoogle Scholar
  43. Rapp JP: Genetic analysis of inherited hypertension in the rat. Physiol Rev. 2000, 80: 135-72.PubMedGoogle Scholar
  44. Drake TA, Schadt EE, Lusis AJ: Integrating genetic and gene expression data: application to cardiovascular and metabolic traits in mice. Mamm Genome. 2006, 17: 466-79. 10.1007/s00335-005-0175-z.PubMedPubMed CentralGoogle Scholar
  45. Fritsche LG, Fariss RN, Stambolian D, Abecasis GR, Curcio CA, Swaroop A: Age-Related Macular Degeneration: Genetics and Biology Coming Together. Annu Rev Genomics Hum Genet. 2014Google Scholar
  46. Tuo J, Grob S, Zhang K, Chan CC: Genetics of immunological and inflammatory components in age-related macular degeneration. Ocul Immunol Inflamm. 2012, 20: 27-36. 10.3109/09273948.2011.628432.PubMedPubMed CentralGoogle Scholar
  47. Logue MW, Schu M, Vardarajan BN, Farrell J, Lunetta KL, Jun G, Baldwin CT, Deangelis MM, Farrer LA: A search for age-related macular degeneration risk variants in Alzheimer disease genes and pathways. Neurobiol Aging. 2014, 35: 1510.e7-18. 10.1016/j.neurobiolaging.2013.12.007.Google Scholar
  48. Ohno-Matsui K: Parallel findings in age-related macular degeneration and Alzheimer's disease. Prog Retin Eye Res. 2011, 30: 217-38. 10.1016/j.preteyeres.2011.02.004.PubMedGoogle Scholar
  49. Cheung CMG, Wong TY: Is age-related macular degeneration a manifestation of systemic disease? New prospects for early intervention and treatment. J Intern Med. 2014, 276: 140-53. 10.1111/joim.12227.PubMedGoogle Scholar
  50. Stanton CM, Wright AF: Inflammatory biomarkers for AMD. Adv Exp Med Biol. 2014, 801: 251-7. 10.1007/978-1-4614-3209-8_32.PubMedGoogle Scholar
  51. Camelo S: Potential Sources and Roles of Adaptive Immunity in Age-Related Macular Degeneration: Shall We Rename AMD into Autoimmune Macular Disease?. Autoimmune Dis. 2014, 2014: 532487-PubMedPubMed CentralGoogle Scholar
  52. Kozhevnikova OS, Korbolina EE, Ershov NI, Kolosova NG: Rat retinal transcriptome: effects of aging and AMD-like retinopathy. Cell Cycle. 2013, 12: 1745-61. 10.4161/cc.24825.PubMedPubMed CentralGoogle Scholar
  53. Markova E V, Obukhova LA, Kolosova NG: Activity of cell immune response and open field behavior in Wistar and OXYS rats. Bull Exp Biol Med. 2003, 136: 377-9.PubMedGoogle Scholar
  54. Clark EA, King WG, Brugge JS, Symons M, Hynes RO: Integrin-mediated signals regulated by members of the rho family of GTPases. J Cell Biol. 1998, 142: 573-86. 10.1083/jcb.142.2.573.PubMedPubMed CentralGoogle Scholar
  55. Schöbel S, Neumann S, Hertweck M, Dislich B, Kuhn P-H, Kremmer E, Seed B, Baumeister R, Haass C, Lichtenthaler SF: A novel sorting nexin modulates endocytic trafficking and alpha-secretase cleavage of the amyloid precursor protein. J Biol Chem. 2008, 283: 14257-68. 10.1074/jbc.M801531200.PubMedGoogle Scholar
  56. Zhao Y, Wang Y, Yang J, Wang X, Zhao Y, Zhang X, Zhang Y: Sorting nexin 12 interacts with BACE1 and regulates BACE1-mediated APP processing. Mol Neurodegener. 2012, 7: 30-10.1186/1750-1326-7-30.PubMedPubMed CentralGoogle Scholar
  57. Arisi I, D'Onofrio M, Brandi R, Felsani A, Capsoni S, Drovandi G, Felici G, Weitschek E, Bertolazzi P, Cattaneo A: Gene expression biomarkers in the brain of a mouse model for Alzheimer's disease: mining of microarray data by logic classification and feature selection. J Alzheimers Dis. 2011, 24: 721-38.PubMedGoogle Scholar
  58. Song MG, Bail S, Kiledjian M: Multiple Nudix family proteins possess mRNA decapping activity. RNA. 2013, 19: 390-9. 10.1261/rna.037309.112.PubMedPubMed CentralGoogle Scholar
  59. Nookala AR, Shah A, Noel RJ, Kumar A: HIV-1 Tat-mediated induction of CCL5 in astrocytes involves NF-κB, AP-1, C/EBPα and C/EBPγ transcription factors and JAK, PI3K/Akt and p38 MAPK signaling pathways. PLoS One. 2013, 8: e78855-10.1371/journal.pone.0078855.PubMedPubMed CentralGoogle Scholar
  60. Yan C, Wang X, Cao J, Wu M, Gao H: CCAAT/enhancer-binding protein γ is a critical regulator of IL-1β-induced IL-6 production in alveolar epithelial cells. PLoS One. 2012, 7: e35492-10.1371/journal.pone.0035492.PubMedPubMed CentralGoogle Scholar
  61. Cloutier A, Guindi C, Larivée P, Dubois CM, Amrani A, McDonald PP: Inflammatory cytokine production by human neutrophils involves C/EBP transcription factors. J Immunol. 2009, 182: 563-71. 10.4049/jimmunol.182.1.563.PubMedGoogle Scholar
  62. Yamanaka R, Lekstrom-Himes J, Barlow C, Wynshaw-Boris A, Xanthopoulos KG: CCAAT/enhancer binding proteins are critical components of the transcriptional regulation of hematopoiesis (Review). Int J Mol Med. 1998, 1: 213-21.PubMedGoogle Scholar
  63. Huggins CJ, Malik R, Lee S, Salotti J, Thomas S, Martin N, Quiñones OA, Alvord WG, Olanich ME, Keller JR, Johnson PF: C/EBPγ suppresses senescence and inflammatory gene expression by heterodimerizing with C/EBPβ. Mol Cell Biol. 2013, 33: 3242-58. 10.1128/MCB.01674-12.PubMedPubMed CentralGoogle Scholar
  64. Parkin SE, Baer M, Copeland TD, Schwartz RC, Johnson PF: Regulation of CCAAT/enhancer-binding protein (C/EBP) activator proteins by heterodimerization with C/EBPgamma (Ig/EBP). J Biol Chem. 2002, 277: 23563-72. 10.1074/jbc.M202184200.PubMedGoogle Scholar
  65. Cooper C, Henderson A, Artandi S, Avitahl N, Calame K: Ig/EBP (C/EBP gamma) is a transdominant negative inhibitor of C/EBP family transcriptional activators. Nucleic Acids Res. 1995, 23: 4371-7. 10.1093/nar/23.21.4371.PubMedPubMed CentralGoogle Scholar
  66. Mullins DN, Crawford EL, Khuder SA, Hernandez D-A, Yoon Y, Willey JC: CEBPG transcription factor correlates with antioxidant and DNA repair genes in normal bronchial epithelial cells but not in individuals with bronchogenic carcinoma. BMC Cancer. 2005, 5: 141-10.1186/1471-2407-5-141.PubMedPubMed CentralGoogle Scholar
  67. Forristal C, Henley SA, MacDonald JI, Bush JR, Ort C, Passos DT, Talluri S, Ishak CA, Thwaites MJ, Norley CJ, Litovchick L, DeCaprio JA, DiMattia G, Holdsworth DW, Beier F, Dick FA: Loss of the Mammalian DREAM Complex Deregulates Chondrocyte Proliferation. Mol Cell Biol. 2014, 34: 2221-34. 10.1128/MCB.01523-13.PubMedPubMed CentralGoogle Scholar
  68. Lee CJ, Chan WI, Scotting PJ: CIC, a gene involved in cerebellar development and ErbB signaling, is significantly expressed in medulloblastomas. J Neurooncol. 2005, 73: 101-8. 10.1007/s11060-004-4598-2.PubMedGoogle Scholar
  69. Zhang L, Wang F, Jiang Y, Xu S, Lu F, Wang W, Sun X, Sun X: Migration of retinal pigment epithelial cells is EGFR/PI3K/AKT dependent. Front Biosci (Schol Ed). 2013, 5: 661-71.Google Scholar
  70. Hewing NJ, Weskamp G, Vermaat J, Farage E, Glomski K, Swendeman S, Chan RVP, Chiang MF, Khokha R, Anand-Apte B, Blobel CP: Intravitreal injection of TIMP3 or the EGFR inhibitor erlotinib offers protection from oxygen-induced retinopathy in mice. Invest Ophthalmol Vis Sci. 2013, 54: 864-70. 10.1167/iovs.12-10954.PubMedPubMed CentralGoogle Scholar
  71. Wan J, Ramachandran R, Goldman D: HB-EGF is necessary and sufficient for Müller glia dedifferentiation and retina regeneration. Dev Cell. 2012, 22: 334-47. 10.1016/j.devcel.2011.11.020.PubMedPubMed CentralGoogle Scholar
  72. Lagna G, Kovelman R, Sukegawa J, Roeder RG: Cloning and characterization of an evolutionarily divergent DNA-binding subunit of mammalian TFIIIC. Mol Cell Biol. 1994, 14: 3053-64.PubMedPubMed CentralGoogle Scholar
  73. Crepaldi L, Policarpi C, Coatti A, Sherlock WT, Jongbloets BC, Down TA, Riccio A: Binding of TFIIIC to sine elements controls the relocation of activity-dependent neuronal genes to transcription factories. PLoS Genet. 2013, 9: e1003699-10.1371/journal.pgen.1003699.PubMedPubMed CentralGoogle Scholar
  74. Avella MA, Baibakov B, Dean J: A single domain of the ZP2 zona pellucida protein mediates gamete recognition in mice and humans. J Cell Biol. 2014, 205: 801-9. 10.1083/jcb.201404025.PubMedPubMed CentralGoogle Scholar
  75. Pasmooij AMG, Nijenhuis M, Brander R, Jonkman MF: Natural gene therapy may occur in all patients with generalized non-Herlitz junctional epidermolysis bullosa with COL17A1 mutations. J Invest Dermatol. 2012, 132: 1374-83. 10.1038/jid.2011.477.PubMedGoogle Scholar
  76. Reactome is a free, open-source, curated and peer reviewed pathway database.Google Scholar
  77. Diamond MS: IFIT1: A dual sensor and effector molecule that detects non-2'-O methylated viral RNA and inhibits its translation. Cytokine Growth Factor Rev. 2014Google Scholar
  78. Roby JA, Clarke BD, Khromykh AA: Loop de loop: viral RNA evades IFIT1 targeting. Trends Microbiol. 2014, 22: 171-3. 10.1016/j.tim.2014.02.014.PubMedGoogle Scholar
  79. Bervini S, Herzog H: Mouse models of Prader-Willi Syndrome: a systematic review. Front Neuroendocrinol. 2013, 34: 107-19. 10.1016/j.yfrne.2013.01.002.PubMedGoogle Scholar
  80. Gray TA, Saitoh S, Nicholls RD: An imprinted, mammalian bicistronic transcript encodes two independent proteins. Proc Natl Acad Sci USA. 1999, 96: 5616-21. 10.1073/pnas.96.10.5616.PubMedPubMed CentralGoogle Scholar
  81. Miyazaki K, Mapendano CK, Fuchigami T, Kondo S, Ohta T, Kinoshita A, Tsukamoto K, Yoshiura K, Niikawa N, Kishino T: Developmentally dynamic changes of DNA methylation in the mouse Snurf/Snrpn gene. Gene. 2009, 432: 97-101. 10.1016/j.gene.2008.11.019.PubMedGoogle Scholar
  82. Stark C, Breitkreutz BJ, Reguly T, Boucher L, Breitkreutz A, Tyers M: BioGRID: a general repository for interaction datasets. Nucleic Acids Res. 2006, D535-9. 34 DatabaseGoogle Scholar
  83. Gilmour J, Assi SA, Jaegle U, Kulu D, van de Werken H, Clarke D, Westhead DR, Philipsen S, Bonifer C: A crucial role for the ubiquitously expressed transcription factor Sp1 at early stages of hematopoietic specification. Development. 2014, 141: 2391-401. 10.1242/dev.106054.PubMedPubMed CentralGoogle Scholar
  84. Lerner LE, Peng G-H, Gribanova YE, Chen S, Farber DB: Sp4 is expressed in retinal neurons, activates transcription of photoreceptor-specific genes, and synergizes with Crx. J Biol Chem. 2005, 280: 20642-50. 10.1074/jbc.M500957200.PubMedGoogle Scholar
  85. Hertel J, Hirche C, Wissmann C, Ebert MP, Höcker M: Transcription of the vascular endothelial growth factor receptor-3 (VEGFR3) gene is regulated by the zinc finger proteins Sp1 and Sp3 and is under epigenetic control: transcription of vascular endothelial growth factor receptor 3. Cell Oncol (Dordr). 2014, 37: 131-45.Google Scholar
  86. Hirano Y, Hayashi H, Iemura S-I, Hendil KB, Niwa S-I, Kishimoto T, Kasahara M, Natsume T, Tanaka K, Murata S: Cooperation of multiple chaperones required for the assembly of mammalian 20S proteasomes. Mol Cell. 2006, 24: 977-84. 10.1016/j.molcel.2006.11.015.PubMedGoogle Scholar
  87. Wang Y, Mandelkow E: Degradation of tau protein by autophagy and proteasomal pathways. Biochem Soc Trans. 2012, 40: 644-52. 10.1042/BST20120071.PubMedGoogle Scholar
  88. Chadwick L, Gentle L, Strachan J, Layfield R: Review: unchained maladie - a reassessment of the role of Ubb(+1) -capped polyubiquitin chains in Alzheimer's disease. Neuropathol Appl Neurobiol. 2012, 38: 118-31. 10.1111/j.1365-2990.2011.01236.x.PubMedGoogle Scholar
  89. Lee MJ, Lee JH, Rubinsztein DC: Tau degradation: the ubiquitin-proteasome system versus the autophagy-lysosome system. Prog Neurobiol. 2013, 105: 49-59.PubMedGoogle Scholar
  90. Southan C, Hancock JM: A tale of two drug targets: the evolutionary history of BACE1 and BACE2. Front Genet. 2013, 4: 293-PubMedPubMed CentralGoogle Scholar
  91. Bettegazzi B, Mihailovich M, Di Cesare A, Consonni A, Macco R, Pelizzoni I, Codazzi F, Grohovaz F, Zacchetti D: β-Secretase activity in rat astrocytes: translational block of BACE1 and modulation of BACE2 expression. Eur J Neurosci. 2011, 33: 236-43. 10.1111/j.1460-9568.2010.07482.x.PubMedGoogle Scholar
  92. Mok KY, Jones EL, Hanney M, Harold D, Sims R, Williams J, Ballard C, Hardy J: Polymorphisms in BACE2 may affect the age of onset Alzheimer's dementia in Down syndrome. Neurobiol Aging. 2014, 35: 1513.e1-5. 10.1016/j.neurobiolaging.2013.12.022.Google Scholar
  93. Hu Y, Chen Y, Lin M, Lee K, Mott RA, Ma J: Pathogenic role of the Wnt signaling pathway activation in laser-induced choroidal neovascularization. Invest Ophthalmol Vis Sci. 2013, 54: 141-54. 10.1167/iovs.12-10281.PubMedPubMed CentralGoogle Scholar
  94. Inestrosa NC, Varela-Nallar L: Wnt signaling in the nervous system and in Alzheimer's disease. J Mol Cell Biol. 2014, 6: 64-74. 10.1093/jmcb/mjt051.PubMedGoogle Scholar
  95. Rosso SB, Inestrosa NC: WNT signaling in neuronal maturation and synaptogenesis. Front Cell Neurosci. 2013, 7: 103-PubMedPubMed CentralGoogle Scholar
  96. Lieven O, Rüther U: The Dkk1 dose is critical for eye development. Dev Biol. 2011, 355: 124-37. 10.1016/j.ydbio.2011.04.023.PubMedGoogle Scholar
  97. Scott EL, Brann DW: Estrogen regulation of Dkk1 and Wnt/β-Catenin signaling in neurodegenerative disease. Brain Res. 2013, 1514: 63-74.PubMedPubMed CentralGoogle Scholar
  98. Lyros O, Rafiee P, Nie L, Medda R, Jovanovic N, Schmidt J, Mackinnon A, Venu N, Shaker R: Dickkopf-1, the Wnt antagonist, is induced by acidic pH and mediates epithelial cellular senescence in human reflux esophagitis. Am J Physiol Gastrointest Liver Physiol. 2014, 306: G557-74. 10.1152/ajpgi.00153.2013.PubMedPubMed CentralGoogle Scholar
  99. Li L, Ke K, Tan X, Xu W, Shen J, Zhai T, Xu L, Rui Y, Zheng H, Zhai P, Zhao J, Cao M: Up-regulation of NFATc4 involves in neuronal apoptosis following intracerebral hemorrhage. Cell Mol Neurobiol. 2013, 33: 893-905. 10.1007/s10571-013-9955-2.PubMedGoogle Scholar
  100. Xu Y, Yang L, Yu S, Shu Q, Yang C, Wang J, Xu F, Sang A, Liang X: Spatiotemporal changes in NFATc4 expression of retinal ganglion cells after light-induced damage. J Mol Neurosci. 2014, 53: 69-77. 10.1007/s12031-013-0198-y.PubMedGoogle Scholar
  101. Inatani M, Tanihara H, Honjo M, Hangai M, Kresse H, Honda Y: Expression of proteoglycan decorin in neural retina. Invest Ophthalmol Vis Sci. 1999, 40: 1783-91.PubMedGoogle Scholar
  102. Lee NE, Park YJ, Chung IY, Seo SW, Park JM, Yoo JM, Song JK: Gene expression changes in a rat model of oxygen-induced retinopathy. Korean J Ophthalmol. 2011, 25: 42-7. 10.3341/kjo.2011.25.1.42.PubMedPubMed CentralGoogle Scholar
  103. Park YJ, Kim YH, Choi WS, Chung IY, Yoo JM: Treatment with triamcinolone acetonide prevents decreased retinal levels of decorin in a rat model of oxygen-induced retinopathy. Curr Eye Res. 2010, 35: 657-63. 10.3109/02713681003760143.PubMedGoogle Scholar
  104. Zhang W, Davis CM, Edin ML, Lee CR, Zeldin DC, Alkayed NJ: Role of endothelial soluble epoxide hydrolase in cerebrovascular function and ischemic injury. PLoS One. 2013, 8: e61244-10.1371/journal.pone.0061244.PubMedPubMed CentralGoogle Scholar
  105. Laube G, Seidenbecher CI, Richter K, Dieterich DC, Hoffmann B, Landwehr M, Smalla KH, Winter C, Böckers TM, Wolf G, Gundelfinger ED, Kreutz MR: The neuron-specific Ca2+-binding protein caldendrin: gene structure, splice isoforms, and expression in the rat central nervous system. Mol Cell Neurosci. 2002, 19: 459-75. 10.1006/mcne.2001.1078.PubMedGoogle Scholar
  106. Mikhaylova M, Sharma Y, Reissner C, Nagel F, Aravind P, Rajini B, Smalla K-H, Gundelfinger ED, Kreutz MR: Neuronal Ca2+ signaling via caldendrin and calneurons. Biochim Biophys Acta. 2006, 1763: 1229-37. 10.1016/j.bbamcr.2006.08.047.PubMedGoogle Scholar
  107. Landwehr M, Redecker P, Dieterich DC, Richter K, Böckers TM, Gundelfinger ED, Kreutz MR: Association of Caldendrin splice isoforms with secretory vesicles in neurohypophyseal axons and the pituitary. FEBS Lett. 2003, 547: 189-92. 10.1016/S0014-5793(03)00713-0.PubMedGoogle Scholar
  108. Seidenbecher CI, Reissner C, Kreutz MR: Caldendrins in the inner retina. Adv Exp Med Biol. 2002, 514: 451-63. 10.1007/978-1-4615-0121-3_27.PubMedGoogle Scholar
  109. Kim KY, Scholl ES, Liu X, Shepherd A, Haeseleer F, Lee A: Localization and expression of CaBP1/caldendrin in the mouse brain. Neuroscience. 2014, 268: 33-47.PubMedPubMed CentralGoogle Scholar
  110. Sassa T, Kihara A: Metabolism of very long-chain Fatty acids: genes and pathophysiology. Biomol Ther (Seoul). 2014, 22: 83-92. 10.4062/biomolther.2014.017.Google Scholar
  111. Kihara A: Very long-chain fatty acids: elongation, physiology and related disorders. J Biochem. 2012, 152: 387-95. 10.1093/jb/mvs105.PubMedGoogle Scholar
  112. Gambuzza ME, Sofo V, Salmeri FM, Soraci L, Marinoc S, Bramanti P: Toll-Like Receptors in Alzheimer's Disease: a Therapeutic Perspective. CNS Neurol Disord Drug Targets. 2014Google Scholar
  113. Reed-Geaghan EG, Reed QW, Cramer PE, Landreth GE: Deletion of CD14 attenuates Alzheimer's disease pathology by influencing the brain's inflammatory milieu. J Neurosci. 2010, 30: 15369-73. 10.1523/JNEUROSCI.2637-10.2010.PubMedPubMed CentralGoogle Scholar
  114. Landreth GE, Reed-Geaghan EG: Toll-like receptors in Alzheimer's disease. Curr Top Microbiol Immunol. 2009, 336: 137-53.PubMedPubMed CentralGoogle Scholar
  115. Wit J, Frydenberg J, Sarup P, Loeschcke V: Confirming candidate genes for longevity by RT-qPCR using two different genetic backgrounds and selection methods. J Insect Physiol. 2013, 59: 255-62. 10.1016/j.jinsphys.2013.01.005.PubMedGoogle Scholar
  116. Head BP, Patel HH, Insel PA: Interaction of membrane/lipid rafts with the cytoskeleton: impact on signaling and function: membrane/lipid rafts, mediators of cytoskeletal arrangement and cell signaling. Biochim Biophys Acta. 2014, 1838: 532-45. 10.1016/j.bbamem.2013.07.018.PubMedGoogle Scholar
  117. Chen K, Williams KJ: Molecular mediators for raft-dependent endocytosis of syndecan-1, a highly conserved, multifunctional receptor. J Biol Chem. 2013, 288: 13988-99. 10.1074/jbc.M112.444737.PubMedPubMed CentralGoogle Scholar
  118. Fabelo N, Martín V, Marín R, Moreno D, Ferrer I, Díaz M: Altered lipid composition in cortical lipid rafts occurs at early stages of sporadic Alzheimer's disease and facilitates APP/BACE1 interactions. Neurobiol Aging. 2014, 35: 1801-12. 10.1016/j.neurobiolaging.2014.02.005.PubMedGoogle Scholar
  119. Yamaguchi N, Oyama M, Kozuka-Hata H, Inoue J: Involvement of A20 in the molecular switch that activates the non-canonical NF-κB pathway. Sci Rep. 2013, 3: 2568-PubMedPubMed CentralGoogle Scholar
  120. Vereecke L, Beyaert R, van Loo G: Genetic relationships between A20/TNFAIP3, chronic inflammation and autoimmune disease. Biochem Soc Trans. 2011, 39: 1086-91. 10.1042/BST0391086.PubMedGoogle Scholar
  121. Damrauer SM, Fisher MD, Wada H, Siracuse JJ, da Silva CG, Moon K, Csizmadia E, Maccariello ER, Patel VI, Studer P, Essayagh S, Aird WC, Daniel S, Ferran C: A20 inhibits post-angioplasty restenosis by blocking macrophage trafficking and decreasing adventitial neovascularization. Atherosclerosis. 2010, 211: 404-8. 10.1016/j.atherosclerosis.2010.03.029.PubMedPubMed CentralGoogle Scholar
  122. Catrysse L, Vereecke L, Beyaert R, van Loo G: A20 in inflammation and autoimmunity. Trends Immunol. 2014, 35: 22-31. 10.1016/j.it.2013.10.005.PubMedGoogle Scholar
  123. Ruan Y, Li H, Zhang K, Jian F, Tang J, Song Z: Loss of Yme1L perturbates mitochondrial dynamics. Cell Death Dis. 2013, 4: e896-10.1038/cddis.2013.414.PubMedPubMed CentralGoogle Scholar
  124. Stiburek L, Cesnekova J, Kostkova O, Fornuskova D, Vinsova K, Wenchich L, Houstek J, Zeman J: YME1L controls the accumulation of respiratory chain subunits and is required for apoptotic resistance, cristae morphogenesis, and cell proliferation. Mol Biol Cell. 2012, 23: 1010-23. 10.1091/mbc.E11-08-0674.PubMedPubMed CentralGoogle Scholar
  125. Takano N, Peng Y-J, Kumar GK, Luo W, Hu H, Shimoda LA, Suematsu M, Prabhakar NR, Semenza GL: Hypoxia-inducible factors regulate human and rat cystathionine β-synthase gene expression. Biochem J. 2014, 458: 203-11. 10.1042/BJ20131350.PubMedPubMed CentralGoogle Scholar
  126. Teng H, Wu B, Zhao K, Yang G, Wu L, Wang R: Oxygen-sensitive mitochondrial accumulation of cystathionine β-synthase mediated by Lon protease. Proc Natl Acad Sci USA. 2013, 110: 12679-84. 10.1073/pnas.1308487110.PubMedPubMed CentralGoogle Scholar
  127. Pan P, Parkkila AK, Autio S, Hilvo M, Sormunen R, Pastorekova S, Pastorek J, Haapasalo H, Parkkila S: Brain phenotype of carbonic anhydrase IX-deficient mice. Transgenic Res. 2012, 21: 163-76. 10.1007/s11248-011-9520-z.PubMedGoogle Scholar
  128. Jobling AI, Fang Z, Koleski D, Tymms MJ: Expression of the ETS transcription factor ELF3 in the retinal pigment epithelium. Invest Ophthalmol Vis Sci. 2002, 43: 3530-7.PubMedGoogle Scholar
  129. Tuncel J, Haag S, Yau ACY, Norin U, Baud A, Lönnblom E, Maratou K, Ytterberg AJ, Ekman D, Thordardottir S, Johannesson M, Gillett A, Stridh P, Jagodic M, Olsson T, Fernández-Teruel A, Zubarev RA, Mott R, Aitman TJ, Flint J, Holmdahl R: Natural polymorphisms in Tap2 influence negative selection and CD4:CD8 lineage commitment in the rat. PLoS Genet. 2014, 10: e1004151-10.1371/journal.pgen.1004151.PubMedPubMed CentralGoogle Scholar
  130. Rienzo M, Costa V, Scarpato M, Schiano C, Casamassimi A, Grimaldi V, Ciccodicola A, Napoli C: RNA-Seq for the identification of novel Mediator transcripts in endothelial progenitor cells. Gene. 2014, 547: 98-105. 10.1016/j.gene.2014.06.034.PubMedGoogle Scholar
  131. Houten SM, Denis S, Argmann CA, Jia Y, Ferdinandusse S, Reddy JK, Wanders RJA: Peroxisomal L-bifunctional enzyme (Ehhadh) is essential for the production of medium-chain dicarboxylic acids. J Lipid Res. 2012, 53: 1296-303. 10.1194/jlr.M024463.PubMedPubMed CentralGoogle Scholar
  132. Bruun GH, Doktor TK, Andresen BS: A synonymous polymorphic variation in ACADM exon 11 affects splicing efficiency and may affect fatty acid oxidation. Mol Genet Metab. 110: 122-8.Google Scholar
  133. Tacutu R, Budovsky A, Fraifeld VE: The NetAge database: a compendium of networks for longevity, age-related diseases and associated processes. Biogerontology. 2010, 11: 513-22. 10.1007/s10522-010-9265-8.PubMedGoogle Scholar
  134. Wolfson M, Budovsky A, Tacutu R, Fraifeld V: The signaling hubs at the crossroad of longevity and age-related disease networks. Int J Biochem Cell Biol. 2009, 41: 516-20. 10.1016/j.biocel.2008.08.026.PubMedGoogle Scholar
  135. Bing DH, Almeda S, Isliker H, Lahav J, Hynes RO: Fibronectin binds to the C1q component of complement. Proc Natl Acad Sci USA. 1982, 79: 4198-201. 10.1073/pnas.79.13.4198.PubMedPubMed CentralGoogle Scholar
  136. Moussavi-Harami SF, Annis DS, Ma W, Berry SM, Coughlin EE, Strotman LN, Maurer LM, Westphall MS, Coon JJ, Mosher DF, Beebe DJ: Characterization of molecules binding to the 70K N-terminal region of fibronectin by IFAST purification coupled with mass spectrometry. J Proteome Res. 2013, 12: 3393-404. 10.1021/pr400225p.PubMedGoogle Scholar
  137. Benoit ME, Hernandez MX, Dinh ML, Benavente F, Vasquez O, Tenner AJ: C1q-induced LRP1B and GPR6 proteins expressed early in Alzheimer disease mouse models, are essential for the C1q-mediated protection against amyloid-β neurotoxicity. J Biol Chem. 2013, 288: 654-65. 10.1074/jbc.M112.400168.PubMedPubMed CentralGoogle Scholar
  138. Yang P, Baciu P, Kerrigan BCP, Etheridge M, Sung E, Toimil BA, Berchuck JE, Jaffe GJ: Retinal pigment epithelial cell death by the alternative complement cascade: role of membrane regulatory proteins, calcium, PKC, and oxidative stress. Invest Ophthalmol Vis Sci. 2014, 55: 3012-21. 10.1167/iovs.13-13554.PubMedPubMed CentralGoogle Scholar
  139. Skeie JM, Mahajan VB: Proteomic Landscape of the Human Choroid-Retinal Pigment Epithelial Complex. JAMA Ophthalmol. 2014Google Scholar
  140. Sobrin L, Seddon JM: Nature and nurture- genes and environment- predict onset and progression of macular degeneration. Prog Retin Eye Res. 2014, 40: 1-15.PubMedGoogle Scholar
  141. Seddon JM: Genetic and environmental underpinnings to age-related ocular diseases. Invest Ophthalmol Vis Sci. 2013, 54: ORSF28-30. 10.1167/iovs.13-13234.PubMedPubMed CentralGoogle Scholar
  142. Keilhauer CN, Fritsche LG, Guthoff R, Haubitz I, Weber BH: Age-related macular degeneration and coronary heart disease: evaluation of genetic and environmental associations. Eur J Med Genet. 2013, 56: 72-9. 10.1016/j.ejmg.2012.10.005.PubMedGoogle Scholar
  143. Cavar I, Lovrić S, Vukojević M, Sesar I, Petric-Vicković I, Sesar A: Metabolic risk factors, coping with stress, and psychological well-being in patients with age-related macular degeneration. Acta Clin Croat. 2014, 53: 79-87.PubMedGoogle Scholar
  144. Ersoy L, Ristau T, Hahn M, Karlstetter M, Langmann T, Dröge K, Caramoy A, den Hollander AI, Fauser S: Genetic and environmental risk factors for age-related macular degeneration in persons 90 years and older. Invest Ophthalmol Vis Sci. 2014, 55: 1842-7. 10.1167/iovs.13-13420.PubMedGoogle Scholar
  145. Peña AS: Genetics of inflammatory bowel disease. The candidate gene approach: susceptibility versus disease heterogeneity. Dig Dis. 16: 356-63.Google Scholar
  146. Baird PN, Chakrabarti S: How genetic studies have advanced our understanding of age-related macular degeneration and their impact on patient care: a review. Clin Experiment Ophthalmol. 42: 53-64.Google Scholar
  147. Janky R, Verfaillie A, Imrichová H, Van de Sande B, Standaert L, Christiaens V, Hulselmans G, Herten K, Naval Sanchez M, Potier D, Svetlichnyy D, Kalender Atak Z, Fiers M, Marine J-C, Aerts S: iRegulon: from a gene list to a gene regulatory network using large motif and track collections. PLoS Comput Biol. 2014, 10: e1003731-10.1371/journal.pcbi.1003731.PubMedPubMed CentralGoogle Scholar
  148. Hawkins RD, Hon GC, Ren B: Next-generation genomics: an integrative approach. Nat Rev Genet. 2010, 11: 476-86.PubMedPubMed CentralGoogle Scholar
  149. Qiu J, Kleineidam A, Gouraud S, Tieng Yao S, Greenwood M, Hoe SZ, Hindmarch C, Murphy D: The use of protein-DNA, chromatin immunoprecipitation, and transcriptome arrays to describe transcriptional circuits in the dehydrated male rat hypothalamus. Endocrinology. 2014, en20141448-Google Scholar
  150. Costa V, Angelini C, De Feis I, Ciccodicola A: Uncovering the complexity of transcriptomes with RNA-Seq. J Biomed Biotechnol. 2010, 2010: 853916-PubMedPubMed CentralGoogle Scholar
  151. Mutz KO, Heilkenbrinker A, Lönne M, Walter JG, Stahl F: Transcriptome analysis using next-generation sequencing. Curr Opin Biotechnol. 2013, 24: 22-30. 10.1016/j.copbio.2012.09.004.PubMedGoogle Scholar
  152. Su Z, Li Z, Chen T, Li QZ, Fang H, Ding D, Ge W, Ning B, Hong H, Perkins RG, Tong W, Shi L: Comparing next-generation sequencing and microarray technologies in a toxicological study of the effects of aristolochic acid on rat kidneys. Chem Res Toxicol. 2011, 24: 1486-93. 10.1021/tx200103b.PubMedGoogle Scholar
  153. Su Z, Ning B, Fang H, Hong H, Perkins R, Tong W, Shi L: Next-generation sequencing and its applications in molecular diagnostics. Expert Rev Mol Diagn. 2011, 11: 333-43.PubMedGoogle Scholar
  154. Xing H, Liao W, Mo Y, Zhang MQ: A novel Bayesian change-point algorithm for genome-wide analysis of diverse ChIPseq data types. J Vis Exp. 2012, e4273-Google Scholar
  155. Yu Y, Fuscoe JC, Zhao C, Guo C, Jia M, Qing T, Bannon DI, Lancashire L, Bao W, Du T, Luo H, Su Z, Jones WD, Moland CL, Branham WS, Qian F, Ning B, Li Y, Hong H, Guo L, Mei N, Shi T, Wang KY, Wolfinger RD, Nikolsky Y, Walker SJ, Duerksen-Hughes P, Mason CE, Tong W, Thierry-Mieg J, et al: A rat RNA-Seq transcriptomic BodyMap across 11 organs and 4 developmental stages. Nat Commun. 2014, 5: 3230-PubMedPubMed CentralGoogle Scholar
  156. Cowley AW, Moreno C, Jacob HJ, Peterson CB, Stingo FC, Ahn KW, Liu P, Vannucci M, Laud PW, Reddy P, Lazar J, Evans L, Yang C, Kurth T, Liang M: Characterization of biological pathways associated with a 1.37 Mbp genomic region protective of hypertension in Dahl S rats. Physiol Genomics. 2014, 46: 398-410. 10.1152/physiolgenomics.00179.2013.PubMedPubMed CentralGoogle Scholar
  157. Stefanova NA, Fursova AZ, Kolosova NG: Behavioral effects induced by mitochondria-targeted antioxidant SkQ1 in Wistar and senescence-accelerated OXYS rats. J Alzheimers Dis. 2010, 21: 479-91.PubMedGoogle Scholar
  158. Laulederkind SJF, Hayman GT, Wang S-J, Smith JR, Lowry TF, Nigam R, Petri V, de Pons J, Dwinell MR, Shimoyama M, Munzenmaier DH, Worthey EA, Jacob HJ: The Rat Genome Database 2013--data, tools and users. Brief Bioinform. 2013, 14: 520-6. 10.1093/bib/bbt007.PubMedPubMed CentralGoogle Scholar
  159. Genoanalitika Lab, Moscow. [http://www.genoanalytica.ru/]
  160. Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R: The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009, 25: 2078-9. 10.1093/bioinformatics/btp352.PubMedPubMed CentralGoogle Scholar
  161. Oracle Corporation. [http://www.oracle.com/index.html]
  162. R, a free language and software environment for statistical computing and graphics.Google Scholar
  163. Martin M: Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnetjournal, North Am. 2011, 9:Google Scholar
  164. Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL: TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013, 14: R36-10.1186/gb-2013-14-4-r36.PubMedPubMed CentralGoogle Scholar

Copyright

© Korbolina et al.; licensee BioMed Central Ltd. 2014

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.