Dynamic Remodeling of the Stressed Heart: Role of Protein Degradation Pathways

The heart is a remarkable organ. In order to maintain its function, it remodels in response to a variety of environmental stresses, including pressure overload, volume overload, mechanical or pharmacological unloading and hormonal or metabolic disturbances. All these responses are linked to the inherent capacity of the heart to rebuild itself. Particularly, cardiac pressure overload activates signaling pathways of both protein synthesis and degradation. While much is known about regulators of protein synthesis, little is known about regulators of protein degradation in hypertrophy. The ubiquitin-proteasome system (UPS) selectively degrades unused and abnormal intracellular proteins. I speculated that the UPS may play an important role in both qualitative and quantitative changes in the composition of heart muscle during hypertrophic remodeling. My study hypothesized that cardiac remodeling in response to hypertrophic stimuli is a dynamic process that requires activation of highly regulated mechanisms of protein degradation as much as it requires protein synthesis. My first aim was to adopt a model of left ventricular hypertrophy and determine its gene expression and structural changes. Male Sprague-Dawley rats were submitted to ascending aortic banding and sacrificed at 7 and 14 days after surgery. Sham operated animals served as controls. Effective aortic banding was confirmed by hemodynamic assessment by Doppler flow measurements in vivo. Banded rats showed a four-fold increase in peak stenotic jet velocities. Histomorphometric analysis revealed a significant increase in myocyte size as well as fibrosis in the banded animals. Transcript analysis showed that banded animals had reverted to the fetal gene program. My second aim was to assess if the UPS is increased and transcriptionally regulated in hypertrophic left ventricular remodeling. Protein extracts from the left ventricles of the banded and control animals were used to perform an in vitro peptidase assay to assess the overall catalytic activity of the UPS. The results showed no difference between hypertrophied and control animals. Transcript analysis revealed decreases in transcript levels of candidate UPS genes in the hypertrophied hearts at 7 days post-banding but not at 14 days. However, protein expression analysis showed no difference at either time point compared to controls. These findings indicate that elements of the UPS are downregulated in the early phase of hypertrophic remodeling and normalizes in a later phase. The results provide evidence in support of a dynamic transcriptional regulation of a major pathway of intracellular protein degradation in the heart. The discrepancy between transcript levels on the one hand and protein levels on the other hand supports post-transcriptional regulation of the UPS pathway in the hypertrophied heart. The exact mechanisms and the functional consequences remain to be elucidated.

POU domain transcription factor Brn3b is critical for the normal development and survival of retinal ganglion cells

The POU domain transcription factor Brn3b/POU4F2 plays a critical role regulating gene expression in mouse retinal ganglion cells (RGCs). Previous investigations have shown that Brn3b is not required for initial cell fate specification or migration; however, it is essential for normal RGC differentiation. In contrast to wild type axons, the mutant neurites were phenotypically different: shorter, rougher, disorganized, and poorly fasciculated. Wild type axons stained intensely with axon specific marker tau-1, while mutant projections were weakly stained and the mutant projections showed strong labeling with dendrite specific marker MAP2. Brn-3b mutant axonal projections contained more microtubules and fewer neurofilaments, a dendritic characteristic, than the wild type. The mutant neurites also exhibited significantly weaker staining of neurofilament low-molecular-weight (NF-L) in the axon when compared to the wild type, and NF-L accumulation in the neuron cell body. The absence of Brn-3b results in an inability to form normal axons and enhanced apoptosis in RGCs, suggesting that Brn-3b may control a set of genes involved in axon formation. ^ Brn3b contains several distinct sequence motifs: a glycine/serine rich region, two histidine rich regions, and a fifteen amino acid conserved sequence shared by all Brn3 family members in the N-terminus and a POU specific and POU homeodomain in the C-terminus. Brn3b activates a Luciferase reporter over 25 fold in cell culture when binding to native brn3 binding sites upstream of a minimal promoter. When fused to the Gal4 DNA Binding domain (DBD) and driven by either a strong (CMV) or weaker (pAHD) promoter, the N-terminal of Brn3b is capable of similar activation when binding to Gal4 UAS sites, indicating a presumptive activator of transcription. Both full length Brn3b or the C-terminus fused to the Gal4DBD and driven by pCMV repressed a Luciferase reporter downstream of UAS binding sites. Lower levels of expression of the fusion protein driven by pADH resulted in an alleviation of repression. This repression appears to be a limitation of this system of transcriptional analysis and a potential pitfall in conventional pCMV based transfection assays. ^

The normal range of the leukocyte count: Implications for surveillance of occupational and environmental exposure and the practice of preventive medicine

Leukopenia, the leukocyte count, and prognosis of disease are interrelated; a systematic search of the literature was undertaken to ascertain the strength of the evidence. One hundred seventy-one studies were found from 1953 onward pertaining to the predictive capabilities of the leukocyte count. Of those studies, 42 met inclusion criteria. An estimated range of 2,200cells/μL to 7,000cells/μL was determined as that which indicates good prognosis in disease and indicates the least amount of risk to an individual overall. Tables of the evidence are included indicating the disparate populations examined and the possible degree of association. ^

The effect of genotype-by-environment interaction on longitudinal triglyceride profiles in the Bogalusa Heart Study

Triglyceride levels are a component of plasma lipids that are thought to be an important risk factor for coronary heart disease and are influenced by genetic and environmental factors, such as single nucleotide polymorphisms (SNPs), alcohol intake, and smoking. This study used longitudinal data from the Bogalusa Heart Study, a biracial community-based survey of cardiovascular disease risk factors. A sample of 1191 individuals, 4 to 38 years of age, was measured multiple times from 1973 to 2000. The study sample consisted of 730 white and 461 African American participants. Individual growth models were developed in order to assess gene-environment interactions affecting plasma triglycerides over time. After testing for inclusion of significant covariates and interactions, final models, each accounting for the effects of a different SNP, were assessed for fit and normality. After adjustment for all other covariates and interactions, LIPC -514C/T was found to interact with age3, age2, and age and a non-significant interaction of CETP -971G/A genotype with smoking status was found (p = 0.0812). Ever-smokers had higher triglyceride levels than never smokers, but persons heterozygous at this locus, about half of both races, had higher triglyceride levels after smoking cessation compared to current smokers. Since tobacco products increase free fatty acids circulating in the bloodstream, smoking cessation programs have the potential to ultimately reduce triglyceride levels for many persons. However, due to the effect of smoking cessation on the triglyceride levels of CETP -971G/A heterozygotes, the need for smoking prevention programs is also demonstrated. Both smoking cessation and prevention programs would have a great public health impact on minimizing triglyceride levels and ultimately reducing heart disease. ^

A comparison of lipoprotein measurements from the Dallas Heart Study

Conventional cholesterol markers in clinical practice today may systematically underestimate the true atherosclerotic risk of populations with high prevalence of metabolic perturbations. It has been suggested that atherogenic risk indexes that measure the concentration of atherogenic particle concentration rather then cholesterol may improve the recognition of atherogenic risk in a clinical setting. Particle concentration is strongly correlated with cholesterol markers, but only a fair concordance with cholesterol has been seen in male populations with low prevalence of metabolic perturbations. Little is known about the concordance of particle concentration and cholesterol markers in multiethnic populations with high prevalence of metabolic perturbations including both men and women. Furthermore, no study has looked at atherosclerosis while exploring the concordance of particle concentration and cholesterol. NMR total atherogenic particle concentration (LipoScience, Inc.), Non-HDL-C, and coronary CT were performed on 3054 subjects ages 30-65 from the Dallas Heart Study, a multi-ethnic probability-based population study. Patients were stratified into four groups: subjects with a low Non-HDL-C and low particle concentration (n = 929), subjects with high Non-HDL-C and low particle concentration (n = 88), subjects with low Non-HDL-C and high particle concentration, and subjects with high Non-HDL-C and high particle concentration (n = 950). When discordance was defined as two quintiles or more of disagreement, discordant groups were relatively small (n= 389, 12.6% of population). There was no statistically significant difference in prevalence of coronary calcification for the group with high Non-HDL-C and low particle concentration compared to the group with low Non-HDL-C and low particle concentration. The discordant group with low Non-HDL-C and low particle concentration, which included 88 subjects, had the highest prevalence of coronary calcification out of the four groups. Out of the 3054 subjects tested in this study, 88 subjects were considered to be part of the discordant group with low Non-HDL-C and a high particle concentration. Although this group is relatively small and comprise approximately 3% of the total population, they did have the highest prevalence of coronary calcification.^

Racial differences in heart failure with preserved ejection fraction in the ambulatory heart failure population

Racial differences in heart failure with preserved ejection fraction (HFpEF) have rarely been studied in an ambulatory, financially "equal access" cohort, although the majority of such patients are treated as outpatients. ^ Retrospective data was collected from 2,526 patients (2,240 Whites, 286 African American) with HFpEF treated at 153 VA clinics, as part of the VA External Peer Review Program (EPRP) between October 2000 and September 2002. Kaplan Meier curves (stratified by race) were created for time to first heart failure (HF) hospitalization, all cause hospitalization and death and Cox proportional multivariate regression models were constructed to evaluate the effect of race on these outcomes. ^ African American patients were younger (67.7 ± 11.3 vs. 71.2 ± 9.8 years; p < 0.001), had lower prevalence of atrial fibrillation (24.5 % vs. 37%; p <0.001), chronic obstructive pulmonary disease (23.4 % vs. 36.9%, p <0.001), but had higher blood pressure (systolic blood pressure > 120 mm Hg 77.6% vs. 67.8%; p < 0.01), glomerular filtration rate (67.9 ± 31.0 vs. 61.6 ± 22.6 mL/min/1.73 m2; p < 0.001), anemia (56.6% vs. 41.7%; p <0.001) as compared to whites. African Americans were found to have higher risk adjusted rate of HF hospitalization (HR 1.52, 95% CI 1.1 - 2.11; p = 0.01), with no difference in risk-adjusted all cause hospitalization (p = 0.80) and death (p= 0.21). ^ In a financially "equal access" setting of the VA, among ambulatory patients with HFpEF, African Americans have similar rates of mortality and all cause hospitalization but have an increased risk of HF hospitalizations compared to whites.^

A conditional Mdm2 allele shows the importance of Mdm2 in heart development

Over 50% of sporadic tumors in humans have a p53 mutation highlighting its importance as a tumor suppressor. Considering additional mutations in other genes involved in p53 pathways, every tumor probably has mutant p53 or impaired p53-mediated functions. In response to a variety of cellular and genotoxic stresses, p53, mainly through its transcriptional activity, induces pathways involved in apoptosis and growth arrest. In these circumstances and under normal situations, p53 must be tightly regulated. Mdm2 is an important regulator of p53. Mdm2 inhibits p53 function by binding and blocking its transactivation domain. In addition, Mdm2 helps target p53 for degradation through its E3 ligase activity. Mdm2 null mice are embryonic lethal due to apoptosis in the blastocysts. However, a p53 null background rescues this lethality demonstrating the importance of the p53-Mdm2 interaction, particularly during development. The lethality of the Mdm2 null mouse prior to implantation limits the ability to investigate the role of Mdm2 in regulating p53 in a temporal and tissue specific manner. Does p53 need to be regulated in all tissues throughout the life of a mouse? Does Mdm2 always have to regulate it? To address these questions, we created a conditional Mdm2 allele. The conditional allele, Mdm2FM, in the presence of Cre recombinase results in the deletion of exons 5 and 6 of Mdm2 (most of the p53 binding domain) and represents a null allele. ^ The Mdm2FM allele was crossed with a heart muscle specific Cre expressing mouse (α-myosin heavy chain promoter driven Cre) to ask whether Mdm2 acts as a negative regulator of p53 in the heart. The heart is the most prominent organ early in embryogenesis and is shaped by cell death and proliferation. p53 does not appear to be active in the heart in response to some types of stress, so it remained to be determined if it has to be regulated in normal heart development. Loss of Mdm2 in the heart results in heart defects as early as E9.5. Loss of Mdm2 results in stabilized p53 and apoptosis. This apoptosis leads to a thinning of the myocardial wall particularly in the ventricles and abnormal ventricular structure. Eventually the abnormal heart fails resulting in lethality by E13.5. The embryonic lethality is rescued in a p53 null background. Thus, Mdm2 is important in regulating p53 in the development of the heart. ^

Citric acid cycle flux and pressure-volume work in the isolated working rat heart

It has been demonstrated previously that the mammalian heart cannot sustain physiologic levels of pressure-volume work if ketone bodies are the only substrates for respiration. In order to determine the metabolic derangement responsible for contractile failure in hearts utilizing ketone bodies, rat hearts were prefused at a near-physiologic workload in a working heart apparatus with acetoacetate and competing or alternate substrates including glucose, lactate, pyruvate, propionate, leucine, isoleucine, valine and acetate. While the pressure-volume work for hearts utilizing glucose was stable for 60 minutes of perfusion, performance fell by 30 minutes for hearts oxidizing acetoacetate as the sole substrate. The tissue content of 2-oxoglutarate and its transamination product, glutamate, were elevated in hearts utilizing acetoacetate while succinyl-CoA was decreased suggesting impaired flux through the citric acid cycle at the level of 2-oxoglutarate dehydrogenase. Further studies indicated that the inhibition of 2-oxoglutarate dehydrogenase developed prior to the onset of contractile failure and that the inhibition of the enzyme may be related to sequestration of the required cofactor, coenzyme A, as the thioesters acetoacetyl-CoA and acetyl-CoA. The contractile failure was not observed when glucose, lactate, pyruvate, propionate, valine or isoleucine were present together with acetoacetate, but the addition of acetate or leucine to acetoacetate did not improve performance indicating that improved performance is not mediated through the provision of additional acetyl-CoA. Furthermore, addition of competing substrates that improved function did not relieve the inhibition of 2-oxoglutarate dehydrogenase and actually resulted in the further accumulation of citric acid cycle intermediates "upstream" of 2-oxoglutarate dehydrogenase (2-oxoglutarate, glutamate, citrate and malate). Studies with (1-$\sp{14}$C) pyruvate indicate that the utilization of ketone bodies is associated with activation of NADP$\sp+$dependent malic enzyme and enrichment of the C4 pool of the citric acid cycle. The results suggest that contractile failure induced by ketone bodies in rat heart results from inhibition of 2-oxoglutarate dehydrogenase and that reversal of contractile failure is dissociated from relief of the inhibition, but rather is due to the entry of carbon units into the citric acid cycle as compounds other than acetyl-CoA. This mechanism of enrichment (anaplerosis) provides oxaloacetate for condensation with acetyl-CoA derived from ketone bodies allowing continued energy production by sustaining flux through a span of the citric acid cycle up to the point of inhibition at 2-oxoglutarate dehydrogenase for energy production thereby producing the reducing equivalents necessary to sustain oxidative phosphorylation. ^

'Your Heart Can't Wait'---A worksite program: Increasing the effectiveness of worksite coordinators to implement a community-wide health education program

A worksite health education program called “Your Heart Can't Wait,” was designed by the American Heart Association Gulf Coast Area (AHA). The objectives were to educate individuals about the signs and symptoms of heart attacks and the actions they should take to improve heart attack victims' chances for survival. AHA volunteers agreed to serve as mentors for this program. ^ A study was designed to determine if worksite coordinators who had the assistance of experienced AHA volunteers had higher rates of program adoption and implementation than worksite coordinators without assistance. Ninety-seven companies participated in the study. Twelve AHA volunteers were randomly assigned to work with forty-three of the worksite coordinators. Mentor/mentee contact forms were used to assess the mentoring process during the course of the study. Program adoption forms were used to measure rates of program adoption and follow-up questionnaires were used to measure rates of program implementation after the study was completed. The twelve mentors were interviewed to provide information for improving future mentoring efforts. ^ Thirty-eight companies completed program adoption forms and fifty-one companies reported using YHCW program components. For the most part, the volunteer mentors did not spend a significant amount of time contacting or working with their assigned worksite coordinators. As a result, the planned analysis comparing the implemented programs between worksite coordinators with and without assistance could not be completed. ^ Additional analyses were performed comparing the implemented programs based upon whether the companies had existing health education/health promotion programs and whether the worksite coordinators had experience using AHA Heart At Work program components. ^ Recommendations based on the mentor interviews were made to improve the success of volunteer assistance programs in the future. ^

The roles of Smads and Pitx2 in cardiac development

The heart is the first organ to form in vertebrates during embryogenesis, and its circulatory function is essential to embryonic survival. Cardiac morphogenesis comprises a complex series of interactions involving cells from several embryonic origins. These cell-cell interactions are regulated temporally and spatially by programs of inductive signaling events, including BMP signaling transduced by Smads and left-right asymmetry signaling mediated by Pitx2. Disruptions of BMP signaling and left-right asymmetry signaling result in abnormal cardiac morphogenesis that causes congenital heart disease in humans. In this study, conventional and conditional gene targeting approaches were employed to dissect the functions of Smad8 and Smad1, intracellular BMP signaling transducers, and Pitx2, a direct target of left-right signaling, in cardiac development. We generated the Smad8mt mutant allele and the Smad8lacZ knock-in allele. Smad8 homozygous mutant mice were viable and fertile without obvious abnormalities. The Smad8lacZ knock-in allele showed that Smad8 was expressed in the myocardium of cardiac outflow tract and atrioventricular cushions. We did not find defects in these Smad8-expressing cardiac regions in Smad8mt/mt and Smad8lacZ/lacZ mutants, indicating that Smad8 is dispensable for cardiac development. Conditional knockout of Smad1 using the Nkx2.5Cre allele in cardiac mesoderm resulted in partial inactivation of Smad1 in the myocardium and complete deletion of Smad1 in the epicardium, and caused ventricular hypoplasia featured with a thinner compact zone, suggesting that Smad1 signaling in the epicardium is required for myocardial morphogenesis in ventricles. Previous data have shown that Pitx2 null mutants exhibit defects in the cardiac outflow tract, a region populated with cells from the cardiac mesoderm and the cardiac neural crest. We found that the cardiac neural crest normally populated into the outflow tract in Pitx2 null mutant. Moreover, specific deletion of Pitx2 in the neural crest resulted in normal heart formation. Deletion of Pitx2 in the cardiac mesoderm caused defective outflow tract, revealing that the function of Pitx2 in the cardiac outflow tract resides in splanchnic and branchial arch mesoderm, and is independent of cardiac neural crest cells. ^

Identification and characterization of new genes involved in Drosophila heart development

Heart development is a crucial and conserved process that is related to the major type of human birth defects. Dorsal vessel, the Drosophila heart, has been regarded as an insightful system to identify new genes and study gene functions involved in heart development. Using heart-specific GFP transgenes, I did a genetic screen for cardiogenic genes on Drosophila chromosome II. Drosophila mutants that carry chromosome II deficiencies were tested for their phenotypes of heart development. Based on the screen results, chromosome regions containing genes required for heart development were identified. Fly strains with single gene mutations located within the defined deficiency regions were tested further. Seven genes have been identified to be involved in heart development. ^ The LIM homeodomain transcription factor gene tailup (tup) was further studied for its function in heart development. Based on this study, tup is expressed in cardioblasts and pericardial cells of the heart tube, as well as in associated lymph glands and alary muscles. In depth analysis of tup mutant phenotypes demonstrated tup is required for normal development of both heart and lymph glands. Tup was shown to bind to two DNA recognition sequences in the dorsal vessel enhancer of the Hand bHLH transcription factor gene, with one site proven essential for the expression of Hand in lymph glands, pericardial cells, and Svp/Doc cardioblasts. Together, these studies demonstrate that Tup is a critical new transcription factor in dorsal vessel morphogenesis and lymph gland formation, and strongly suggest Tup is a direct regulator of the expression of Hand in these developmental processes. ^

The differential staurosporine mediated G1 arrest in normal versus tumor cells is dependent on the retinoblastoma and Chk1 proteins

The ability to regulate cell cycle progression is one of the differences that separates normal from tumor cells. A protein, which is frequently mutated or deleted in a majority of tumor cells, is the retinoblastoma protein (pRb). Previously, we reported that normal cells, which have a wild-type Rb pathway, can be reversibly arrested in the G1 phase of the cell cycle by staurosporine (ST), while tumor cells were unaffected by this treatment. As a result, ST may be used to protect normal cells against the toxic affects of chemotherapy. Here we set out to determine the mechanism(s) by which ST can mediate a reversible G1 arrest in pRb positive cells. To this end, we used an isogenic cell model system of normal human mammary epithelial cells (HMEC) with either intact pRb+ (p53-) or p53+ (pRb-) treated with ST. Our results show that pRb+ cells treated with low concentrations of ST, arrested in the G1 phase of the cell cycle; however, in pRb - cells there was no response. This was verified as a true G 1 arrest in pRb+ cells by two different methods for monitoring cell cycle kinetics and in two additional model systems for Rb (i.e. pRb -/- mouse embryo fibroblasts, and downregulation of RB with siRNA). Our results indicated that ST-mediated G1 arrest required pRb, which in turn initiated a cascade of events leading to inhibition of CDK4 and CDK2 activities and up-regulation of p21 protein. Further assessment of this pathway revealed the novel finding that Chk1 expression and activity were required for the Rb-dependent, ST-mediated G1 arrest. In fact, overexpression of Chk1 facilitated recovery from ST-mediated G1 arrest, an effect only observed in RB+ cells. Collectively, our data suggest pRb is able to cooperate with Chk1 to mediate a G1 arrest in pRb+ cells, but not in pRb- cells. The elucidation of this pathway can help identify novel agents that can be used to protect cancer patients against the debilitating affects of chemotherapy, by targeting only the normal proliferating cells in the body that are otherwise destroyed. ^

The information bottleneck method for genome-wide association studies

In population studies, most current methods focus on identifying one outcome-related SNP at a time by testing for differences of genotype frequencies between disease and healthy groups or among different population groups. However, testing a great number of SNPs simultaneously has a problem of multiple testing and will give false-positive results. Although, this problem can be effectively dealt with through several approaches such as Bonferroni correction, permutation testing and false discovery rates, patterns of the joint effects by several genes, each with weak effect, might not be able to be determined. With the availability of high-throughput genotyping technology, searching for multiple scattered SNPs over the whole genome and modeling their joint effect on the target variable has become possible. Exhaustive search of all SNP subsets is computationally infeasible for millions of SNPs in a genome-wide study. Several effective feature selection methods combined with classification functions have been proposed to search for an optimal SNP subset among big data sets where the number of feature SNPs far exceeds the number of observations. ^ In this study, we take two steps to achieve the goal. First we selected 1000 SNPs through an effective filter method and then we performed a feature selection wrapped around a classifier to identify an optimal SNP subset for predicting disease. And also we developed a novel classification method-sequential information bottleneck method wrapped inside different search algorithms to identify an optimal subset of SNPs for classifying the outcome variable. This new method was compared with the classical linear discriminant analysis in terms of classification performance. Finally, we performed chi-square test to look at the relationship between each SNP and disease from another point of view. ^ In general, our results show that filtering features using harmononic mean of sensitivity and specificity(HMSS) through linear discriminant analysis (LDA) is better than using LDA training accuracy or mutual information in our study. Our results also demonstrate that exhaustive search of a small subset with one SNP, two SNPs or 3 SNP subset based on best 100 composite 2-SNPs can find an optimal subset and further inclusion of more SNPs through heuristic algorithm doesn't always increase the performance of SNP subsets. Although sequential forward floating selection can be applied to prevent from the nesting effect of forward selection, it does not always out-perform the latter due to overfitting from observing more complex subset states. ^ Our results also indicate that HMSS as a criterion to evaluate the classification ability of a function can be used in imbalanced data without modifying the original dataset as against classification accuracy. Our four studies suggest that Sequential Information Bottleneck(sIB), a new unsupervised technique, can be adopted to predict the outcome and its ability to detect the target status is superior to the traditional LDA in the study. ^ From our results we can see that the best test probability-HMSS for predicting CVD, stroke,CAD and psoriasis through sIB is 0.59406, 0.641815, 0.645315 and 0.678658, respectively. In terms of group prediction accuracy, the highest test accuracy of sIB for diagnosing a normal status among controls can reach 0.708999, 0.863216, 0.639918 and 0.850275 respectively in the four studies if the test accuracy among cases is required to be not less than 0.4. On the other hand, the highest test accuracy of sIB for diagnosing a disease among cases can reach 0.748644, 0.789916, 0.705701 and 0.749436 respectively in the four studies if the test accuracy among controls is required to be at least 0.4. ^ A further genome-wide association study through Chi square test shows that there are no significant SNPs detected at the cut-off level 9.09451E-08 in the Framingham heart study of CVD. Study results in WTCCC can only detect two significant SNPs that are associated with CAD. In the genome-wide study of psoriasis most of top 20 SNP markers with impressive classification accuracy are also significantly associated with the disease through chi-square test at the cut-off value 1.11E-07. ^ Although our classification methods can achieve high accuracy in the study, complete descriptions of those classification results(95% confidence interval or statistical test of differences) require more cost-effective methods or efficient computing system, both of which can't be accomplished currently in our genome-wide study. We should also note that the purpose of this study is to identify subsets of SNPs with high prediction ability and those SNPs with good discriminant power are not necessary to be causal markers for the disease.^