Modeling the interactions of Alzheimer-related genes from the whole brain microarray data and diffusion tensor images of human brain
BMC Bioinformatics volume 13, Article number: S10 (2012)
In recent years the genome-wide microarray-based gene expression profiles and diffusion tensor images (DTI) in human brain have been made available with accompanying anatomic and histology data. The challenge is to integrate various types of data to investigate the interactions of genes that are associated with specific neurological disorder.
In this study, we analyzed the whole brain microarray data and the physical connectivity of the hippocampus with other brain regions to identify the genes related to Alzheimer's disease and their interactions with proteins. We generated a physical connectivity map of the left and right hippocampuses with 12 other brain regions and identified 33 Alzheimer-related genes that interact with many proteins. These genes are highly linked to the development of Alzheimer's disease.
In Alzheimer's brain both brain regions and inter-regional communications through the white matter are often hampered. So far the connectivity of regions in Alzheimer's brain has been studied mostly at the functional level using functional MRI (fMRI). Analyzing the inter-regional fiber connectivity without tracking crossing-fiber regions often provides coarse and inaccurate results. A few deep brain fibers were analyzed but the inter-regional fiber connectivity was not analyzed in their studies. The inter-regional fiber connectivity analysis can provide comprehensive and measurable degradation of fiber tracts in AD patients' brains, but is not easy to perform. We tracked crossing-fiber regions and identified genes with high expression levels in the fiber pathways of the hippocampus. The interactions of the genes with other proteins can provide comprehensive and measurable degradation of fiber tracts in Alzheimer brains. To the best of our knowledge, this is the first attempt to integrate the whole brain microarray data with DTI data to identify specific genes and their interactions.
Alzheimer's disease (AD) is one of the most devastating diseases for people in advanced age between 40 to 60 years. It affects parts of the brain that control thought, memory, and sometimes language. The most affected part is the memory which is primarily controlled by the hippocampus. Thus, a comprehensive study of the connectivity of the hippocampus with other brain regions may reveal a new pattern specific to AD.
So far the connectivity of regions in Alzheimer's brain has been studied mostly at the functional level using functional MRI (fMRI) [1, 2]. It is until recently that the importance of the physical connectivity of regions in Alzheimer's brain has become evident [3–6]. The cause of AD is now believed to be the excessive storage of Amyloid Beta (Aβ) plaque  into the white fiber tracts. So, the study of relationship between fiber tracts and Amyloid Beta is very important. For the past few years diffusion tensor imaging (DTI) has received much attention from brain studies since it can capture the fine fiber orientations of brain white matter. Water molecules diffuse in the direction of the fiber more rapidly due to the myelination property of the tracts, and DTI records the diffusion. DTI including information of diffusion can provide inter-regional fiber pathway.
Each of the DTI data and the microarray data for a human brain takes much time and effort to generate. Thus, most DTI data was generated independently from brain microarray data so far. Recently both DTI and gene expression data were generated from two human brains that have no known neuropsychiatric or neuropathological history . Thus, it is timely to track inter-region physical connectivity pathways from the DTI data and find Alzheimer-related genes expressed in the pathways to identify the interactions of the genes. Analyzing genes expressed in specific fiber pathways will help identify the interactions of genes in neurons, which in turn will provide a valuable resource to research and treatment of neurological disorders such as Alzheimer's disease.
In this study, we tracked inter-region fiber tracts using a new method for detecting crossing fibers in voxels and streamlined the fibers by probabilistic tractography  to generate an inter-region physical connectivity map of the hippocampus. The primary focus of this study is to investigate the connectivity of the left and right hippocampus with other brain regions. The secondary focus of this study is to find the interactions of Alzheimer-associated genes that are expressed in the fiber pathways of hippocampuses. The rest of this paper presents our method of tracking the fiber pathways and finding genes highly expressed in the fiber pathways and major results of our study.
Both the brain microarray data and the DTI data for a 24 year old male were obtained from the Allen Brain Atlas (ABA, http://human.brain-map.org/). In the microarray, there are 58,692 probes to detect for human brain genes. We selected 473 genes that are known to be relevant with the Alzheimer disease. The diffusion tensor image (DTI) for the same brain was obtained from ABA.
To identify protein-protein interactions (PPIs) that involve brain genes, we extracted 39,194 PPIs from the Human Protein Reference Database (HPRD). 473 Alzheimer's disease genes are associated with 96 human proteins in HPRD. 81 out of the 96 proteins have 1,506 interactions.
Human brain fiber tractography
The left part of Figure 1 shows the process of performing tractography from the DTI data of human brain. First, we converted the DTI data from .hdr and .img formats to the FSL compatible NiFTi format (*.nii) using the LONI-Inspector (http://www.loni.ucla.edu/Software/LONI-Inspector) program. After conversion of data format, we used the FMRIB Software Library (FSL) for processing DTI. FSL provides a complete set of tools for the processing and analysis of DTI .
3D NiFTi files converted by FSL were merged into a single 4D NiFTi image. The 4D NiFTi image was corrected to remove distortions caused by eddy currents and simple head motion. From corrected 4D image, we extracted standalone 3D image with no diffusion. This 3D image was used as a reference after process. For the non-linear registration with the standard template, we discarded the skull from the reference 3D image using the Brain Extraction Tool (BET) . After discarding the skull we obtained a pure brain without the skull.
Recently the technique of detecting crossing fiber regions that appear multiple principal diffusion direction at each voxel has been widely used for the analysis of DTI. We estimated direction of diffusion using the Bayesian Estimation of Diffusion Parameters Obtained using Sampling Techniques (BEDPOSTX)  to detect two crossing fibers at each voxel. BEDPOSTX made a multi fiber diffusion model using crossing fiber regions.
After detecting crossing fibers, we used the FMRIB's Non-linear Image Registration Tool (FNIRT) for registration. FNIRT generated a warp field and its corresponding inversed warp field for a tractography in the standard space. We selected the MNI152-T1 1mm template for the registration and the probabilistic tractography. The region of interest (ROI), which was used as the seed and target in fiber tracking, was defined in this standard template.
The hippocampus is known to be closely related to learning and memory. The early symptom of Alzheimer's disease is memory impairment, so we chose the hippocampus for the seed region of probabilistic tractography. To identify seed and target ROIs, we used the FMRIB's Integrated Registration and Segmentation Tool (FIRST)  and extracted 15 ROIs including the hippocampus (left and right hippocampus, left and right nucleus accumbens, left and right amygdala, left and right caudate nucleus, left and right globus pallidus, left and right putamen, left and right thalamus, and brainstem) from the standard template (Figure 2). Tractography was then performed in the standard space from every voxel of the seed ROI for other all target ROIs except the brainstem. We traced the pathway by hemisphere (left to left and right to right). Thus we obtained 2 (left and right) × 1 (seed) × 6 (target) = 12 pathways.
Finding samples in tractography
Each of the 946 samples of the human brain microarray data contains 58,692 probes. To identify the genes expressed in fiber pathways, we examined the positions of the voxels in the 946 samples, and selected the samples that are placed in the pathways. The selected samples contain 473 genes that are related to Alzheimer's disease, and 81 genes out of the 473 genes have protein-protein interaction data available at HPRD.
The human brain microarray data has the expression values of the genes as the log2 intensity. For each probe we computed the z-score by the following equation.
where x is the log2 intensity, μ is the mean, and σ is the standard deviation.
Since multiple samples can be placed in a pathway, a gene can have multiple z-scores. In this case, we calculated the average of the z-scores of gene expression levels in each pathway.
where n is the number of samples of a gene in a pathway.
Results and discussion
We traced 12 fiber pathways from the left and right hippocampuses to 6 target regions in each hemisphere. No fiber pathways were found from the hippocampuses to the nucleus accumbens. Table 1 shows 10 fiber pathways of the hippocampuses. Figures 3 and 4 show the pathways visualized by FSL. Figure 3 shows the pathways (red) from the hippocampus (yellow) to the amygdala (sky-blue), the caudate nucleus (yellow-green) and the globus pallidus (orange). Figure 4 shows the pathways (red) from the hippocampus to the putamen (blue) and the thalamus (green). Figure 3A shows the left and right pathways from the hippocampus to the amygdala on the MNI152-T1 1mm template. Likewise, Figures 3B, C, 4A and 4B show the pathways from the left and right hippocampus to target regions.
Seven samples were found in the pathway from the right hippocampus to the right amygdala. Figure 5 shows the z-scores of the expression of the amyloid precursor protein (APP) in each of the 7 samples. APP is known to induce the amyloid beta , which is directly related to Alzheimer's disease. The figure shows different expression levels of the samples. The AVG in Figure 5 represents the average z-scores of the 7 samples. Among the 473 genes that are related to Alzheimer's disease, 81 genes have protein-protein interaction data available at HPRD. Figure 6 shows the average z-scores of the expression levels of the 81 Alzheimer-related genes. G50 in Figure 6, for example, represents the Amyloid beta A4 precursor protein-binding family A member 1 (APBA1), which directly interacts with APP [13, 14].
Among the 81 Alzheimer-related genes, we selected 33 genes that have more than 10 protein-protein interactions (Table 2). Figure 7 shows average z-scores of the expression levels of the 33 genes in the left and right fiber pathways. As shown in Figure 7, all the genes have different expression levels in the left and right fiber pathways of the hippocampus. The 33 genes also have different expression levels from the remaining 48 genes (Figure 8). Figure 9 shows the protein-protein interaction network of the 33 Alzheimer-related genes (yellow nodes) that have more than 10 protein-protein interactions.
The protein-protein interaction network in Figure 10 includes four types of proteins. Yellow nodes (AKT1, CASP3, CYCS, GRIA1, GRIN1, and GSK3B) represent the proteins related to the amyloid precursor protein (APP). Finding CASP3 and GRIA1 as APP-related proteins agrees with a recent study result that activation of caspase-3 (CASP3) via mitochondria is required for the long-term depression and AMPA (GRIA1 in the network) receptor internalization in hippocampal neurons . Red nodes (AR, PRKCZ, CTNNB1, PSEN1, CASP8, and CASP9) represent the proteins that interact with three APP-related proteins (KT1, CASP3, CYCS, GRIA1, GRIN1, and GSK3B). Green nodes represent the proteins that interact with two proteins among APP, AKT1, CASP3, CYCS, GRIA1, GRIN1, and GSK3B. The remaining nodes represent the proteins that interact with the APP-related proteins only. The microarray data was obtained from a male brain, and the node AR in the network represents the androgen receptor which is critical for the development and maintenance of the male sexual phenotype. The average z-scores of the proteins shown in red (AR, PRKCZ, CTNNB1, PSEN1, CASP8, CASP9) are available in Additional file 1.
In this study, we tracked inter-region fibers tracts using a new method for detecting crossing fibers in voxels and streamlined the fibers by probabilistic tractography to generate an inter-region physical connectivity map starting from the hippocampus. The primary result of this study is to investigate the connectivity of the left and right hippocampus with other brain regions. We tracked fiber pathways from the left and right hippocampus to 12 other regions of interest (ROI) of brain. The secondary result of this study is to investigate the microarray gene expression of the fiber pathways and to find protein-protein interactions of the genes, which are related to Alzheimer disease.
Although the microarray data DTI data for a single normal control was analyzed in this study, we plan to investigate more cases in the future study when more data is available. Although preliminary, the genes and their interactions found in our study will be useful to design biochemical experiments for further investigation.
Golby A, Silverberg G, Race E, Gabrieli S, O'Shea J, Knierim K, Stebbins G, Gabrieli J: Memory encoding in Alzheimer's disease: an fMRI study of explicit and implicit memory. Brain 2005, 128(Pt 4):773–787.
Petrella JR, Wang L, Krishnan S, Slavin MJ, Prince SE, Tran TT, Doraiswamy PM: Cortical deactivation in mild cognitive impairment: high-field-strength functional MR imaging. Radiology 2007, 245(1):224–235. 10.1148/radiol.2451061847
Liu Y, Spulber G, Lehtimaki KK, Kononen M, Hallikainen I, Grohn H, Kivipelto M, Hallikainen M, Vanninen R, Soininen H: Diffusion tensor imaging and Tract-Based Spatial Statistics in Alzheimer's disease and mild cognitive impairment. Neurobiol Aging 2011, 32(9):1558–1571. 10.1016/j.neurobiolaging.2009.10.006
Serra L, Cercignani M, Lenzi D, Perri R, Fadda L, Caltagirone C, Macaluso E, Bozzali M: Grey and white matter changes at different stages of Alzheimer's disease. J Alzheimers Dis 2010, 19(1):147–159.
Zhou Y, Dougherty JH Jr, Hubner KF, Bai B, Cannon RL, Hutson RK: Abnormal connectivity in the posterior cingulate and hippocampus in early Alzheimer's disease and mild cognitive impairment. Alzheimers Dement 2008, 4(4):265–270. 10.1016/j.jalz.2008.04.006
Qiu A, Oishi K, Miller MI, Lyketsos CG, Mori S, Albert M: Surface-based analysis on shape and fractional anisotropy of white matter tracts in Alzheimer's disease. PLoS One 2010, 5(3):e9811. 10.1371/journal.pone.0009811
Shankar GM, Li S, Mehta TH, Garcia-Munoz A, Shepardson NE, Smith I, Brett FM, Farrell MA, Rowan MJ, Lemere CA, Regan CM, Walsh DM, Sabatini BL, Selkoe DJ: Amyloid-beta protein dimers isolated directly from Alzheimer's brains impair synaptic plasticity and memory. Nat Med 2008, 14(8):837–842. 10.1038/nm1782
Jones AR, Overly CC, Sunkin SM: The Allen Brain Atlas: 5 years and beyond. Nat Rev Neurosci 2009, 10(11):821–828. 10.1038/nrn2722
Behrens TE, Berg HJ, Jbabdi S, Rushworth MF, Woolrich MW: Probabilistic diffusion tractography with multiple fibre orientations: what can we gain? Neuroimage 2007, 34(1):144–155. 10.1016/j.neuroimage.2006.09.018
Smith SM, Jenkinson M, Woolrich MW, Beckmann CF, Behrens TE, Johansen-Berg H, Bannister PR, De Luca M, Drobnjak I, Flitney DE, Niazy RK, Saunders J, Vickers J, Zhang Y, De Stefano N, Brady JM, Matthews PM: Advances in functional and structural MR image analysis and implementation as FSL. Neuroimage 2004, 23(Suppl 1):S208-S219.
Smith SM: Fast robust automated brain extraction. Hum Brain Mapp 2002, 17(3):143–155. 10.1002/hbm.10062
Patenaude B, Smith SM, Kennedy DN, Jenkinson M: A Bayesian model of shape and appearance for subcortical brain segmentation. Neuroimage 2011, 56(3):907–922. 10.1016/j.neuroimage.2011.02.046
Biederer T, Cao X, Sudhof TC, Liu X: Regulation of APP-dependent transcription complexes by Mint/X11s: differential functions of Mint isoforms. J Neurosci 2002, 22(17):7340–7351.
Borg JP, Ooi J, Levy E, Margolis B: The phosphotyrosine interaction domains of X11 and FE65 bind to distinct sites on the YENPTY motif of amyloid precursor protein. Mol Cell Biol 1996, 16(11):6229–6241.
Li Z, Jo J, Jia JM, Lo SC, Whitcomb DJ, Jiao S, Cho K, Sheng M: Caspase-3 activation via mitochondria is required for long-term depression and AMPA receptor internalization. Cell 2010, 141(5):859–871. 10.1016/j.cell.2010.03.053
This work was supported by the Key Research Institute Program (2011-0018394) and in part by the Basic Science Research Program (2011-0003766) through the National Research Foundation of Korea (NRF) funded by the Ministry of Education, Science and Technology.
This article has been published as part of BMC Bioinformatics Volume 13 Supplement 7, 2012: Advanced intelligent computing theories and their applications in bioinformatics. Proceedings of the 2011 International Conference on Intelligent Computing (ICIC 2011). The full contents of the supplement are available online at http://www.biomedcentral.com/bmcbioinformatics/supplements/13/S7.
The authors declare that they have no competing interests.
BP analyzed the human brain microarray data and prepared the first draft of the manuscript. WL analyzed DTI data to track the fiber pathway and prepared the manuscript together. KH supervised the work and rewrote the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
About this article
Cite this article
Park, B., Lee, W. & Han, K. Modeling the interactions of Alzheimer-related genes from the whole brain microarray data and diffusion tensor images of human brain. BMC Bioinformatics 13 (Suppl 7), S10 (2012). https://doi.org/10.1186/1471-2105-13-S7-S10