BioMed Research International
Volume 2016 (2016), Article ID 5164347, 7 pages
http://dx.doi.org/10.1155/2016/5164347

Research Article

Detecting Susceptibility to Breast Cancer with SNP-SNP Interaction Using BPSOHS and Emotional Neural Networks

Xiao Wang, Qinke Peng, and Yue Fan

Systems Engineering Institute and School of Electronic and Information Engineering, Xi’an Jiaotong University, Xi’an, Shaanxi 710049, China

Received 23 March 2016; Revised 18 April 2016; Accepted 20 April 2016

Academic Editor: Gang Liu

Copyright © 2016 Xiao Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Studies for the association between diseases and informative single nucleotide polymorphisms (SNPs) have received great attention. However, most of them just use the whole set of useful SNPs and fail to consider the SNP-SNP interactions, while these interactions have already been proven in biology experiments. In this paper, we use a binary particle swarm optimization with hierarchical structure (BPSOHS) algorithm to improve the effective of PSO for the identification of the SNP-SNP interactions. Furthermore, in order to use these SNP interactions in the susceptibility analysis, we propose an emotional neural network (ENN) to treat SNP interactions as emotional tendency. Different from the normal architecture, just as the emotional brain, this architecture provides a specific path to treat the emotional value, by which the SNP interactions can be considered more quickly and directly. The ENN helps us use the prior knowledge about the SNP interactions and other influence factors together. Finally, the experimental results prove that the proposed BPSOHS_ENN algorithm can detect the informative SNP-SNP interaction and predict the breast cancer risk with a much higher accuracy than existing methods.

1. Introduction

Breast cancer is a major cause of death among women. Some genes (e.g., BRCA1 and BRCA2) have already been known as the cause of breast cancer [1]. However, only 5% case of breast cancer has these mutations; in this way, these symbols failed to be used for most of women. Fortunately, besides, on these rare mutations, increasing evidence shows that the risk of breast cancer can be measured by the SNPs, which is one of the most common types of mutations for human being [2–4]. Furthermore, with the development of SNP microarray and genome-wide associations studies (GWAS), the research about breast cancer and SNPs becomes more popular.

For complex disease like breast cancer, the effect of an individual SNP is small. Researchers generally focus on the joint genetic effect of SNP combinations which may increase susceptibility to the cancer. However, most of studies just care about the whole set of disease related SNPs and treat all the SNPs in the dataset equally [5]. In this way, they fail to consider the small scale interaction among SNPs. On the other hand, increasing evidence proves that SNP-SNP interactions exist [6–8]. To change this situation, a method to use the SNP-SNP interaction in breast cancer susceptibility identification is in need.

With the development of microarray, one of the most important challenges to detect SNP-SNP interactions is the complex combination of data with increasing SNPs number. In order to improve the effectiveness to identify these interactions, researchers try to use many different algorithms, for example, [9] uses the Genetic Algorithm (GA), [10] uses the Particle Swarm Optimization (PSO) algorithm, and [11] uses the Polymorphism Interaction Analysis (PIA). These methods have the ability to detect SNP interactions in high dimensional dataset; however, due to the random generator initial values and optimizing process, they generally need lots of iterations and are easily trapped into the local optima. So an improved algorithm for solving this puzzle is essential. Here, we propose BPSOHS to improve the performance of identifying the SNP interactions. Inspired by the decision-making process in the real society, we divide the particle swarm as “leader” part and “follower” part, which have different evolutionary strategies. Our former work proved that this novel method is faster than other swarm intelligent algorithms and much easier to converge to the globally optimal solution [12]. So, by this improved PSO, after encoding and matching, we can improve the performance of identifying the SNP interactions related to the cancer.

The defined SNP-SNP interaction is an important tendency factor about the breast cancer [13]. However, there is still no specific method to use this kind of information into the cancer susceptibility. ENN is a novel method inspired by the emotional process of the human brain, which is usually used in the face recognition [14, 15]. In this paper, by designing the emotional value based on the SNP interactions, we explore the range of ENN and use it for susceptibility analysis. So, in this novel ENN, the SNP interaction features can be regarded as the prior tendency of classification and influence the result of classification, just like what emotions do in the process of people making choices.

In the end, we propose a pipeline for the selection and usage of SNP-SNP interactions. At the beginning, the BPSOHS algorithm is used to identify the informative SNP combinations related to the cancer. Then, we transform these new features as an additional vector of each sample. Finally, by treating the vector as tendency emotional value in the ENN, the novel network can particularly consider the SNP-SNP interaction in a more directly and quickly way. According to the output of the networks, our method can measure the breast cancer risk for each sample. The case-control study in 10000 people suggests that our pipeline can detect the useful SNP interactions and effectively consider them in susceptibility analysis to breast cancer.

2. Material and Methods

2.1. Dataset Preparation

The datasets we use in this paper were obtained from the breast cancer case-control study in [16]. This dataset has 5000 controls and 5000 cases, with the SNPs selected by the biological research.

In order to be convenient for the further discussion, we assume that the dataset has samples and each sample has SNPs. The SNP value can be represented as , where 1 represents the major homozygous sites, 2 represents the minor homozygous sites, and 3 represents the heterozygous sites. Then, denote an individual sample and represent the whole dataset.

2.2. Particle Swarm Optimization

PSO is a swarm algorithm developed by Kennedy and Eberhart [17]. For basic PSO, each solution is corresponding to the position of a particle, and the velocity and direction of each particle can be adjusted according to both its own parameter and the best particles in the swarm. The idea of this process is particles “flying” through a multidimensional possible search space and looking for the global optimum. By both individual memory and global memory, a particle calculates and moves to its next position. The basic elements of PSO are described as follows:

where and are the velocity and position of the particle in bit, , are the random numbers within , and , are the acceleration constants. is the best position of the particle in bit, and is the bit of best position of the best particle in the swarm. The velocity is within to make sure that the particle is flying in the range of possible solution space.

2.3. Binary Particle Swarm Optimization with Hierarchical Structure

The basic PSO is just designed for the continuous problems. In order to make it more widely used, in 1997, Kennedy and Eberhart propose BPSO (binary particle swarm optimization) to use in the discrete problems [18]. Different from the basic conventional PSO, in BPSO, each bit of a particle can just move in a state space to 0 or 1 as in the following function:

In this function, denote the probability of choosing 1 and is a random real number within .

In basic PSO and BPSO algorithm, the particles are treated equally. However, this process works much differently from the real society. The researchers about sociology point out that the leaders in the group seem to have the stronger say in making decision. Most of people prefer to follow their leaders, so the particles seem to also have status in the swarm [19]. Inspired by this idea, we propose BPSOHS [12] with two kinds of particles: “leader” particles and “follower” particles, and they can be regarded as the “leader” and “followers” in the society.

At the beginning, there are particles defined as “leader” in the swarm, while others are “follower.” According to the fitness function, the “leader” and “follower” can switch in these two statuses, to make sure that the better particles are “leaders.” Thus, at the iteration, the followers can walk toward the leaders based on the following formula:

where , indicates the “leader,” and indicates the “follower.” And is a parameter to limit the followers’ speed. According to these formulas, when we update the position of particle, we need consider both its own position and the decisions of “leader.”

2.4. Encoding Schemes and Performance Measurement

In BPSOHS, the SNPs in particle can be represented by two bits: . and can be 0 and 1, so a SNP can be represented as follows:

In this way, the combination of can represent the four different states of the SNPs. The particle can be described by

For example, if , this particle represents that we choose the 2nd and 4th SNP in genotypes 3 and 1. Thus, the phenotype of samples can be described by the particles. In this BPSOHS process, the value of fitness function presents the importance of particle. In the susceptibility analysis, we employ the SNP-SNP interactions to influence the tendency of the classification; thus, the ratio difference between the cases and controls is useful. Referring to [9], we use the following formula to measure the importance of SNP interactions:

where denote the elements number in a dataset, while is the dataset subset in the breast cancer case group with specific interaction , and is the dataset subset in the control group with specific interaction .

By using this fitness function and BPSOHS, we get some SNP-SNP interactions in significant association with breast cancer; the details of some of them are shown in Table 1.

Table 1: SNP interactions about breast cancer.

SNP (gene) (chromosome/position)
SNP type
Case number
Control number
Difference
Odds ratio
95% CI
value

rs3020314 (ESR1) (6/152270672)
rs500760 (PGR) (11/100909991)
CC-
AA
1230
1404
174
0.836
0.7644–0.9135
0.001

rs3020314 (ESR1) (6/152270672)
rs2017591 (STS) (X/7158114)
CC-
TC
984
1152
168
0.818
0.7436–0.9008
0.003

rs3020314 (ESR1) (6/152270672)
rs2077647 (ESR1) (6/152129077)
CT-
GG
1055
1213
158
0.835
0.7602–0.9170
0.008

rs500760 (PGR) (11/100909991)
rs2017591 (STS) (X/7158114)
AG-
CC
1326
1476
150
0.862
0.7896–0.9404
0.001

rs2077647 (ESR1) (6/152129077)
rs2017591 (STS) (X/7158114)
AG-
TC
1117
1263
146
0.851
0.7762–0.9333
0.040

rs3020314 (ESR1) (6/152270672)
rs660149 (PGR) (11/100934314)
rs11571171 (PGR) (11/100974887)
CT-
CC-
TT
602
509
93
1.208
1.0657–1.3687
0.003

rs3020314 (ESR1) (6/152270672)
rs500760 (PGR) (11/100909991)
rs2017591 (STS) (X/7158114)
CC-
AA-
TC
571
699
128
0.793
0.7048–0.8929
0.012

rs6269 (COMT) (22/19949952)
rs2175898 (ESR1) (6/152196952)
rs660149 (PGR) (11/100934314)
AA-
AG-
CC
1035
1158
123
0.866
0.7877–0.9522
0.003

rs3020314 (ESR1) (6/152270672)
rs1543404 (ESR1) (6/152428838)
rs2747652 (ESR1) (6/152437016)
rs9340799 (ESR1) (6/152163381)
rs1709182 (ESR1) (6/152175357)
rs9478249 (ESR1) (6/152194431)
rs660149 (PGR) (11/100934314)
rs11571171 (PGR) (11/100974887)
rs858518 (SHBG) (17/7533025)
rs858524 (SHBG) (17/7511287)
rs2017591 (STS) (X/7158114)
CT-
TC-
CT-
AG-
TT-
TG-
CC-
TT-
TC-
AG-
TT
7
1
6
7.01
0.8800–151.68
0.034

In addition to these interactions in 2 or 3 interactions, we still find some interactions with more SNPs. For example, the combination, rs3020314 (CT)-rs1543404 (TC)-rs2747652 (CT)-rs9340799 (AG)-rs1709182 (TT)-rs9478249 (TG)-rs660149 (CC)-rs11571171 (TT)-rs858518 (TC)-rs858524 (AG)-rs2017591 (TT), appears 7 times in the case group and only once in the control group.

2.5. Emotional Neural Networks Based on SNP-SNP Interactions

Neural networks are popular in the field of bioinformatics and disease susceptibility analysis. However, the basic fully connected feed forward neural networks just treat all the input in the same level. Although the powers will adjust by the importance of input data, this architecture still fails to use the inside information about data. Emotional neural networks are a novel method to deal with this problem. According to the related works [20], researches about brain characterize that there are a short path for the emotional signal in brain, which help the feedback become more directly and quickly. The simple modeling of the emotional brain is shown in Figure 1.

Figure 1: Simple modeling of emotional brain.

Referring to the works of Lotfi and Akbarzadeh-T [14], we can use the neural network in Figure 2 to simulate this architecture. In the long path from the thalamus, sensory cortex OFC to the amygdala is the path for the general input. And the short path from the thalamus directly to the amygdala is the path for the tendency emotional input.

Figure 2: Emotional neural networks for susceptibility analysis to breast cancer.

In the related works about ENN architecture [21, 22], the input of the sensory cortex and amygdala is the same, which is different from the real process of our brain. Generally, in the human emotional brain, the emotional value is the a priori tendency information [23, 24]. So, in this paper, based on the BPSOHS and the SNP-SNP interaction detection, we try to use this kind of a priori information in susceptibility analysis. As it is shown in Table 1, the specific phenotype has different emergence probability in case and control class. So these SNP-SNP interactions can influence the cancer risk of the samples.

Assume that we use SNP-SNP interactions as the emotional value; then, the tendency vector is as follows: . According to (7), we can get the final sequence, , as follows:

where is the odds ratio of the interaction, which can measure its importance. Inspired by the ENN architecture focus on the image recognition [25], we design the output of the thalamus as , and the final output of the neural networks is calculated by

where is the weight for the amygdala, while , , is the related weight of , and , , is the related weight about the emotional value . is the weight for the OFC.

In this architecture, these parameters are updated as in the following formula:

where the learning rates are and . The decay rate is . is a binary input about reinforce, and is the reward value calculated according to the BELPIC model [26]:

where

In the training process, the value is used when the goal is presented, while the value is used for other situations. In the test process, the final output is calculated according to (8). By this process, we treat the SNP interactions as the emotional value and use these tendency features to reflect the final output of the neural network. By this architecture, we get the score to measure the risk for breast cancer, and the details will be shown in the next section.

3. Results and Discussion

3.1. Evaluation Function

To measure the performance of the SNP-SNP interactions and the novel neural networks, we need an effective evaluation function. Referring to other works about classification, in this paper, we use Sn (sensitivity), Sp (specificity), Acc (accuracy), RR (risk ratio), and OR (odds ratio); they are defined as follows:

where TP (True Positive) is the number of positive samples that are predicted as positive, TN (True Negative) is the number of negative samples that are predicted as negative, FP (False Positive) is the number of negative samples that are predicted as positive, and FN (False Negative) is the number of positive samples that are predicted as negative.

3.2. Experimental Results and Comparison

In order to prove that the SNP-SNP interaction feature is useful for the susceptibility analysis, we use 5-fold cross-validation to evaluate the performance of our method. We compare the performance of some popular basic neural networks (BP (Back Propagation), RBF (Radial Basis Function), and PNN (Pattern Recognition Neural networks)) with or without the novel information. The details of the experiment are shown in Table 2.

Table 2: The performance of ENN compared with basic neural networks.

Method
Sn
Sp
Acc
RR
OR

BP
0.5486
0.4864
0.5175
1.0681
1.1505

RBF
0.5394
0.5152
0.5273
1.1126
1.2445

PNN
0.5448
0.5066
0.5257
1.1508
1.2808

BPSOHS_BP
0.6462
0.5606
0.6034
1.4706
2.3302

BPSOHS_RBF
0.6290
0.5826
0.6058
1.5069
2.3664

BPSOHS_PNN
0.6406
0.5888
0.6147
1.5579
2.5523

BPSOHS_ENN
0.7268
0.7106
0.7187
2.5114
6.5322

Figure 3 shows that these methods use the SNP-SNP interaction features to achieve better performance than those methods which do not use the small scale information. By considering this tendency information in the specific architecture, our BPSOHS_ENN gets the best performance compared with these basic methods.

Figure 3: Bar graph of the performance of different methods.

To prove that our method is suitable for this problem, we also use our method to compare with some newly published papers. Majority of these researches just give useful SNP interaction, while only part of samples has these interactions. So, to make sure of the fairness of the comparison, we filter these results based on interactions that only appeared in less than 10% of samples, and in Table 3 we compare our result with the best of the remaining ones.

Table 3: The performance of ENN compared with published methods.

Method
Sn
Sp
Acc
RR
OR

IGA [9]
—
—
0.5351
1.08
1.17

IPSO [16]
—
—
0.49
0.82
0.79

IBBA [27]
0.132
0.940
0.619
—
2.384

BPSOHS_ENN
0.7268
0.7106
0.7187
2.5114
6.5322

In Table 3, IBBA method uses the dataset derived by 7 SNPs from CXCL12-related genes, and IGA and IPSO use the same dataset as this paper. The tests on the different dataset also prove the effectiveness of our method. By using the small scale feature and other features in the same neural network architecture, the process of the susceptibility analysis becomes much easier. The tendency of SNP combination can give a more reasonable initial value for the classifier. So consider this a priori information in a fast and direct way to help improve the performance of the neural network.

The experimental results prove that our method has the power to consider these useful SNP-SNP interactions together, which is much better than these single barcodes. What is more, Table 3 shows that our method is competitive compared with the other published methods. Compared with them, our method has significant advantages in identifying the case and control, which is useful for us to predict and prevent the cancer.

4. Conclusions

The susceptibility analysis to disease is normally based on a group of SNPs. However, the small scale relationship in the group is rarely being mentioned. In this paper, we use the BPSOHS to pick out the neglected small scale SNP-SNPs. In order to consider the relationship about the SNP interaction, we propose specific partially connected architecture neural networks. By simulating the process of human brain deal with the emotional value, the related SNPs have the chance to work together with other normal features to calculate possibility of the samples suffering from breast cancer. According to the cancer-related SNP interactions and output of the novel neural network, we can measure the cancer risk of the samples, which is useful for us to prevent the possible cancer.

Competing Interests

The authors declare that there are no competing interests regarding the publication of this paper.

Acknowledgments

Thanks are due to Hsueh-Wei Chang, Ph.D., for providing the dataset in this paper. This work was jointly supported by the National Natural Science Foundation of China (Grants nos. 61173111 and 60774086) and the Ph.D. Programs Foundation of Ministry of Education of China (Grant no. 20090201110027).

References

A. M. Dunning, C. S. Healey, P. D. P. Pharoah, M. D. Teare, B. A. J. Ponder, and D. F. Easton, “A systematic review of genetic polymorphisms and breast cancer risk,” Cancer Epidemiology Biomarkers & Prevention, vol. 8, no. 10, pp. 843–854, 1999.

R. Fagerholm, M. K. Schmidt, S. Khan et al., “The SNP rs6500843 in 16p13.3 is associated with survival specifically among chemotherapy-treated breast cancer patients,” Oncotarget, vol. 6, no. 10, pp. 7390–7407, 2015.

C. C. Ho, K.-S. Mun, and R. Naidu, “SNP array technology: an array of hope in breast cancer research,” Malaysian Journal of Pathology, vol. 35, no. 1, pp. 33–43, 2013.

T. Wang, B. Huang, R. Guo et al., “A let-7b binding site SNP in the 3′-UTR of the Bcl-xL gene enhances resistance to 5-fluorouracil and doxorubicin in breast cancer cells,” Oncology Letters, vol. 9, no. 4, pp. 1907–1911, 2015.

B. Wei, Q. K. Peng, Q. W. Zhang, and C. Y. Li, “Identification of a combination of SNPs associated with Graves' disease using swarm intelligence,” Science China Life Sciences, vol. 54, no. 2, pp. 139–145, 2011.

J. Shang, Y. Sun, S. Li, J.-X. Liu, C.-H. Zheng, and J. Zhang, “An improved opposition-based learning particle swarm optimization for the detection of SNP-SNP interactions,” BioMed Research International, vol. 2015, Article ID 524821, 12 pages, 2015.

J. R. Groom, J. Richmond, and T. T. Murooka, “CXCR3 chemokine receptor-ligand interactions in the lymph node optimize CD4⁺ T helper 1 cell differentiation,” Immunity, vol. 37, no. 6, pp. 1091–1103, 2012.

P.-Y. Hsu, H.-K. Hsu, X. Lan et al., “Amplification of distant estrogen response elements deregulates target genes associated with tamoxifen resistance in breast cancer,” Cancer Cell, vol. 24, no. 2, pp. 197–212, 2013.

C. H. Yang, Y.-D. Lin, L.-Y. Chuang, and H.-W. Chang, “Evaluation of breast cancer susceptibility using improved genetic algorithms to generate genotype SNP barcodes,” Transactions on Computational Biology Bioinformatics, vol. 10, no. 2, pp. 361–371, 2013.

S.-J. Wu, L.-Y. Chuang, Y.-D. Lin et al., “Particle swarm optimization algorithm for analyzing SNP-SNP interaction of renin-angiotensin system genes against hypertension,” Molecular Biology Reports, vol. 40, no. 7, pp. 4227–4233, 2013.

J. E. Goodman, L. E. Mechanic, B. T. Luke, S. Ambs, S. Chanock, and C. C. Harris, “Exploring SNP-SNP interactions and colon cancer risk using polymorphism interaction analysis,” International Journal of Cancer, vol. 118, no. 7, pp. 1790–1797, 2006.

W. Bin, P. Qinke, J. Zhao, and X. Chen, “A binary particle swarm optimization algorithm inspired by multi-level organizational learning behavior,” European Journal of Operational Research, vol. 219, no. 2, pp. 224–233, 2012.

V. U. Onay, L. Briollais, J. A. Knight et al., “SNP-SNP interactions in breast cancer susceptibility,” BMC Cancer, vol. 6, no. 1, article 114, 2006.

E. Lotfi and M.-R. Akbarzadeh-T, “Practical emotional neural networks,” Neural Networks the Official Journal of the International Neural Network Society, vol. 59, no. 11, pp. 61–72, 2014.

A. Khashman, “Modeling cognitive and emotional processes: a novel neural network architecture,” Neural Networks, vol. 23, no. 10, pp. 1155–1163, 2010.

L.-Y. Chuang, Y.-D. Lin, H.-W. Chang, and C.-H. Yang, “An improved PSO algorithm for generating protective SNP barcodes in breast cancer,” PLoS ONE, vol. 7, no. 5, Article ID e37018, 2012.

J. Kennedy and R. Eberhart, “Particle swarm optimization,” in Proceedings of the IEEE International Conference on Neural Networks, vol. 4, pp. 1942–1948, 1995.

J. Kennedy and R. C. Eberhart, “A discrete binary version of the particle swarm algorithm,” in Proceedings of the IEEE International Conference on Systems, Man, and Cybernetics, Computational Cybernetics and Simulation, vol. 5, pp. 4104–4108, Orlando, Fla, USA, October 1997.

S. D. Dionne and P. J. Dionne, “Levels-based leadership and hierarchical group decision optimization: a simulation,” Leadership Quarterly, vol. 19, no. 2, pp. 212–234, 2008.

E. Lotfi and A. Keshavarz, “Gene expression microarray classification using PCA-BEL,” Computers in Biology & Medicine, vol. 54, pp. 180–187, 2014.

E. Lotfi and M.-R. Akbarzadeh-T, “A winner-take-all approach to emotional neural networks with universal approximation property,” Information Sciences, vol. 346-347, pp. 369–388, 2016.

E. Lotfi and M.-R. Akbarzadeh-T, “Brain emotional learning-based pattern recognizer,” Cybernetics & Systems, vol. 44, no. 5, pp. 402–421, 2013.

A. Khashman, “A modified backpropagation learning algorithm with added emotional coefficients,” IEEE Transactions on Neural Networks, vol. 19, no. 11, pp. 1896–1909, 2008.

A. Khashman, “Credit risk evaluation using neural networks: emotional versus conventional models,” Applied Soft Computing, vol. 11, no. 8, pp. 5477–5484, 2011.

E. Lotfi and M.-R. Akbarzadeh-T, “Adaptive brain emotional decayed learning for online prediction of geomagnetic activity indices,” Neurocomputing, vol. 126, no. 3, pp. 188–196, 2014.

E. Lotfi, S. Setayeshi, and S. Taimory, “A neural basis computational model of emotional brain for online visual object recognition,” Applied Artificial Intelligence, vol. 28, no. 8, pp. 814–834, 2014.

L.-Y. Chuang, H.-W. Chang, M.-C. Lin, and C.-H. Yang, “Improved branch and bound algorithm for detecting SNP-SNP interactions in breast cancer,” Journal of Clinical Bioinformatics, vol. 3, no. 1, article 4, 2013.