Genetic and functional association of FAM5C with myocardial infarction
- Jessica J Connelly1,
- Svati H Shah1, 2,
- Jennifer F Doss1,
- Shera Gadson1,
- Sarah Nelson1,
- David R Crosslin1,
- A Brent Hale1,
- Xuemei Lou1,
- Ty Wang1,
- Carol Haynes1,
- David Seo3,
- David C Crossman4,
- Vincent Mooser5,
- Christopher B Granger2,
- Christopher JH Jones6,
- William E Kraus1, 2,
- Elizabeth R Hauser1 and
- Simon G Gregory1Email author
© Connelly et al; licensee BioMed Central Ltd. 2008
Received: 24 January 2008
Accepted: 22 April 2008
Published: 22 April 2008
We previously identified a 40 Mb region of linkage on chromosome 1q in our early onset coronary artery disease (CAD) genome-wide linkage scan (GENECARD) with modest evidence for linkage (n = 420, LOD 0.95). When the data are stratified by acute coronary syndrome (ACS), this modest maximum in the overall group became a well-defined LOD peak (maximum LOD of 2.17, D1S1589/D1S518). This peak overlaps a recently identified inflammatory biomarker (MCP-1) linkage region from the Framingham Heart Study (maximum LOD of 4.27, D1S1589) and a region of linkage to metabolic syndrome from the IRAS study (maximum LOD of 2.59, D1S1589/D1S518). The overlap of genetic screens in independent data sets provides evidence for the existence of a gene or genes for CAD in this region.
A peak-wide association screen (457 SNPs) was conducted of a region 1 LOD score down from the peak marker (168–198 Mb) in a linkage peak for acute coronary syndrome (ACS) on chromosome 1, within a family-based early onset coronary artery disease (CAD) sample (GENECARD).
Polymorphisms were identified within the 'family with sequence similarity 5, member C' gene (FAM5C) that show genetic linkage to and are associated with myocardial infarction (MI) in GENECARD. The association was confirmed in an independent CAD case-control sample (CATHGEN) and strong association with MI was identified with single nucleotide polymorphisms (SNPs) in the 3' end of FAM5C. FAM5C genotypes were also correlated with expression of the gene in human aorta. Expression levels of FAM5C decreased with increasing passage of proliferating aortic smooth muscle cells (SMC) suggesting a role for this molecule in smooth muscle cell proliferation and senescence.
These data implicate FAM5C alleles in the risk of myocardial infarction and suggest further functional studies of FAM5C are required to identify the gene's contribution to atherosclerosis.
Coronary artery disease (CAD), and its extreme manifestation of myocardial infarction (MI), is the leading cause of death in the United States  and, concomitant with the epidemic of obesity and diabetes is rapidly becoming the leading cause of death in many developing countries . The genetic predilection of CAD has been well-established; family history has been repeatedly shown to be a independent risk factor even after adjustment for shared environmental factors [3, 4]. The heritability of CAD is particularly strong in early-onset forms of the disease (< 50 years of age) where the relative risk of developing early-onset CAD in the first-degree sibling is between 3.8 and 12.1, depending on the age-of-onset in the proband .
Myocardial infarction, a serious manifestation of CAD, has also been the subject of genetic analysis. In fact, the only major successes in mapping genes for common, complex cardiovascular disease have focused on the MI phenotype, suggesting that this extreme (but common) manifestation of CAD may be most powerful for disease gene identification. For example, investigators in Iceland carried out a genome wide linkage scan for MI in an Icelandic cohort . Using fine-mapping techniques, the arachidonate 5-lipoxygenase-activating protein (ALOX5AP) gene was identified as the major susceptibility gene for MI on chromosome 13q. Although initial results implicated a single gene, expanded analysis of other leukotriene pathway genes identified additional CAD/MI susceptibility genes . Another study of MI revealed evidence for linkage to chromosome 1p34  and later identified apolipoprotein E receptor 2 (LRP8) as an underlying susceptibility gene for familial and premature CAD and MI . Most recently, four CAD/MI genome wide association studies (GWAS) [10–13] concomitantly identified genetic association of a region on chromosome 9p21.3 adjacent to the cyclin-dependent kinase inhibitor 2A (CDKN2A) and cyclin-dependent kinase inhibitor 2B (CDKN2B) genes within independent MI and CAD populations. A study by Broadbent et al. validated these findings and identified the antisense noncoding RNA in the INK4 locus (ANRIL) as the potential disease-causing gene in this region . Despite these successes, these risk variants in combination do not account for the entire underlying genetic component of MI risk, suggesting further studies to identify additional genetic risk factors are necessary.
We previously identified a 40 megabase region of linkage on chromosome 1q25-31 in 420 families from our early onset CAD genome-wide linkage scan (GENECARD) , showing modest evidence for linkage (n = 420 families, LOD 0.95). A priori stratification by sibling pairs concordant for acute coronary syndrome (ACS; an extreme manifestation of CAD primarily composed of MI) revealed that this modest maximum in the overall group became a well-defined linkage peak (n = 228 families, LOD 2.17, at peak microsatellite marker D1S1589/D1S518). This peak coincides exactly with a quantitative trait locus (QTL) for inflammatory biomarker (monocyte chemoattractant protein 1, MCP-1) linkage in the Framingham Heart Study (LOD 4.27, at peak marker D1S1589) , as well as a linkage peak to metabolic syndrome (MetS) from the IRAS family study (LOD 2.59, at peak marker D1S1589/D1S518) . MCP-1 is a potent inflammatory molecule that promotes monocyte/macrophage accumulation in atherosclerotic plaques, which can lead to plaque instability through chronic inflammation and smooth muscle cell proliferation . Metabolic syndrome is a known strong risk factor for development of CAD and MI [19–21]. The overlap of genetic screens in independent data sets provides strong evidence for the existence of a susceptibility gene or genes for CAD in this region.
We conducted a peak-wide association screen of 457 SNPs across the region one LOD score down from the peak marker (168–198 Mb) on chromosome 1q and identified polymorphisms within the 'family with sequence similarity 5, member C' (FAM5C) gene that show genetic linkage to and are associated with MI. We confirmed association within an independent CAD case-control dataset (CATHGEN) [22, 23] and have identified strong association between SNPs within the 3' end of FAM5C and myocardial infarction. We show that FAM5C expression is correlated with FAM5C genotype in human aortas, suggesting that genetically heritable levels of this molecule in the cells that make up the aorta may play a role in atherosclerosis. Importantly, FAM5C is known to promote proliferation, migration, and invasion of pituitary tumors, a phenotype relevant to the cellular changes of smooth muscle cells that are associated with the formation and vulnerability of an atherosclerotic plaque [25, 26]. We identified FAM5C expression in proliferating aortic smooth muscle cells (AoSMCs) and show that levels of this molecule decrease with increasing passage, suggesting that FAM5C levels may play a similar role in SMCs as in pituitary cell tumors. We hypothesize that polymorphisms within this gene may alter FAM5C levels in smooth muscle cells, enhancing the atherosclerotic smooth muscle cell phenotype, thus leading to plaque instability and changes in MCP-1 levels.
Early-onset CAD family-based sample (GENECARD)
Clinical characteristics of samples used in this study
(n = 422)
(n = 228)
(n = 368)
(n = 289)
(n = 88)
Age of onset (SD)
69.0 (7.0) exam
Family history of CAD
Body Mass Index (SD)
Data not available
Sys/Dias BP (SD)
Total Cholesterol (SD)
MI case-control sample (CATHGEN)
CATHGEN participants were recruited sequentially through the cardiac catheterization laboratories at Duke University Hospital (Durham, North Carolina, United States) with approval from the Duke Institutional Review Board. All participants undergoing catheterization were offered participation in the study and signed informed consent. Medical history and clinical data were collected and stored in the Duke Information System for Cardiovascular Care database maintained at the Duke Clinical Research Institute . Controls and cases were chosen on the basis of extent of CAD as measured by the CAD index (CADi). CADi is a numerical summary of coronary angiographic data that incorporates the extent and anatomical distribution of coronary disease . CADi has been shown to be a better predictor of clinical outcome than extent of CAD . Affected status was determined by the presence of significant CAD defined as a CADi ≥ 32 , as we have used for previous studies [22, 23]. For patients older than 55 years of age, a higher CADi threshold (CADi ≥ 74) was used to adjust for the higher baseline extent of CAD in this group. Given the linkage to MI in the GENECARD study, we also selected a subgroup of Caucasian individuals who experienced MI from the overall CATHGEN group, similar to the GENECARD ACS families. MI was defined as documentation of MI in the medical history. Controls were defined as ≥60 years of age, with no CAD as demonstrated by coronary angiography and no documented history of cerebrovascular or peripheral vascular disease, myocardial infarction, transplant, or interventional or surgical coronary revascularization procedures . A comparison of clinical characteristics between GENECARD probands and CATHGEN MI cases and unaffected controls is presented in Table 1. As a replication data set, the CATHGEN sample provides 80% power to detect effect sizes between 1.25 and 1.45, for allele frequencies ranging from 0.1 to 0.5. The CATHGEN MI cases were also stratified into a young affected group (CAD AOO ≤ 55), which provides a consistent comparison for the GENECARD family study.
Human donor aorta samples and expression
Human aorta samples were collected from heart transplant donors as previously described . DNA and RNA were isolated from each sample. Eighty-eight samples in total were analyzed representing fifty-eight unique samples, as 30 individuals had more than one sample. As harvested tissue was obtained from deceased heart donors, the clinical data associated with these aortas is very limited and consists of age, sex, and race (Table 1). Genotyping was performed using Applied Biosystems Taqman allelic discrimination assays and expression profiling was performed using Affymetrix GeneChip U95Av2 (Affymetrix, Santa Clara, CA). Expression signal intensity values were log2 transformed and normalized using quantile normalization. We analyzed cis effects of FAM5C variants on FAM5C expression using the Affymetrix tag 34442_at, representing the 3' end of FAM5C. Because some subjects had multiple samples while others did not, we treated each individual and sample separately. The expression level of the tag was modeled using multiple linear regression including age, sex, race and additive genotype. To account for repeated measures (i.e., multiple sections per subjects) and as a validation, a mixed model as implemented in the SAS PROC MIXED procedure was utilized (SAS Institute Inc., Cary, NC).
Peak-wide SNP selection
The approximate 30 Mb region of linkage on chromosome 1q contains 138 genes (NCBI Build 36 in Ensembl). We used an iterative fine mapping approach that involved selecting SNPs with a minor allele frequency of at least 5% in the Caucasian population spaced at 100 kilobases (kb) . These SNPs were intercalated with at least one previously identified and confirmed (validated) SNP to represent positional candidate genes from NCBI Build 35 (120 genes). When possible, coding SNPs, SNPs that altered a splice site or that were contained within putative regulatory regions of multi-species sequence conservation (chimp, mouse, rat, dog, chicken, Fugu, zebrafish) were selected for genotyping. If the gene was greater than 100 kb in size, multiple SNPs were selected using the same criteria. Additionally, SNPs were selected from regions that do not contain annotated genes at an average SNP to SNP distance of 100 kb, while taking into account regions of multi-species sequence conservation.
Using these criteria, 231 intragenic and 226 intergenic SNPs were selected for genotype analysis within the chromosome 1q linkage region with an average spacing of 66 kb. Twenty-seven genes contained validated coding SNPs with allele frequencies >5%; twenty-two genes contained synonymous SNPs with allele frequencies >5%; eleven SNPs are found within untranslated regions (5' or 3' UTR); and one hundred and seventy-one intronic SNPs were selected, the reference allele is conserved within fifty-five of these (Additional files, Table 1).
A total of 457 SNPs were genotyped at two genotyping facilities. 384 SNPs were genotyped at the University of Washington through a grant from SeattleSNPs (see Availability and requirements section for URL). Genotyping in the GENECARD sample was performed by SeattleSNPs using the Illumina BeadStation 500 G SNP genotyping system (Illumina, San Diego, California, United States). Each Sentrix Array generates 384 genotypes for 96 individuals; within each individual array experiment, four samples purchased from Coriell were included as quality control samples. The Coriell samples were of Centre d'Etude du Polymorphisme Humain (CEPH) (Caucasian) origin, Asian and Hispanic ethnicity and were used for gender checking and to assess sample quality. In addition, one sample per plate was duplicated for internal quality control. These samples were used to identify possible sample plating errors and genotype calling inconsistencies. Twenty-two SNPs did not cluster well and were not called. The call rate for the remaining 362 SNPs was 99.8%. One replicate error was identified in SNP rs2370025, all other SNPs yielded greater than 99.999% concordance among the replicates.
Ninety-five additional inter- and intragenic SNPs were assayed at the Center for Human Genetics using Taqman allelic discrimination assays. A total of 15 quality control samples – composed of six reference genotype controls in duplicate, two CEPH pedigree individuals, and one no-template sample – were included in each quadrant of the 384-well plate. SNPs that showed mismatches on quality-control samples were reviewed by an independent genotyping supervisor for potential genotyping errors. All SNPs examined were successfully genotyped for 95% or more of the individuals in the study. Error rate estimates for SNPs meeting the quality control benchmarks were determined to be less than 0.2%.
All SNPs were tested for deviations (p ≤ 0.0001) from Hardy-Weinberg equilibrium (HWE) in the affected and unaffected race-stratified groups. Two SNPs deviated from HWE in the peak-wide association study, rs1503122 (cases) and rs438781 (cases and controls). Neither of these SNPs was statistically significant in our association analysis. Interestingly, rs438781 lies within a reported copy number polymorphism, which is most likely the cause of the disequilibrium. In CATHGEN, one SNP, rs9427746, significantly deviated from HWE (p < 0.0001) in the Caucasian control group. This deviation was neither due to mixed racial backgrounds as all of the samples were from Caucasian origin, nor was it due to assay dropout as all of our controls are randomly interspersed with our cases in the 384 well plates. This deviation is either a significant departure representing an association or due to random chance. Linkage disequilibrium between pairs of SNPs was assessed using the Graphical Overview of Linkage Disequilibrium package  and displayed using Haploview .
In the initial screen, two-point linkage using a dominant and recessive model was performed using Fastlink [37, 38] and Homog  and family-based association was tested using the Association in the Presence of Linkage (APL) test . We analyzed the GENECARD families together regardless of race because these families formed the basis of the original linkage evidence. We chose to use the APL test for family based association because this test incorporates data from affected sibling pairs with available parental data and unaffected siblings in the analyses, effectively using all available information in the GENECARD families. The APL software infers missing parental genotypes, appropriately accounts for the non-independence of affected siblings and calculates a robust estimate of the variance. APL results from markers with variance estimates of less than five are viewed as less reliable . A single SNP in TROVE2 was significant in our original screen with a variance less than five and did not replicate.
In addition to APL, family-based association of FAM5C SNPs was performed using the Pedigree Disequilibrium Test (PDT) [42, 43]. We employed this additional test at the gene level in order to obtain the most information about each polymorphism. The PDT test uses related trios and discordant sibships from extended pedigrees and like other family-based association tests is robust to population stratification. In order to determine if replicated SNPs could account for the MI linkage signal we employed Linkage and Association Modeling in Pedigrees (LAMP). LAMP estimates the degree of linkage disequilibrium (LD) between a candidate SNP and the putative disease locus through joint modeling of linkage and association [44, 45].
Allelic association in the CATHGEN MI cases and controls was evaluated using multivariable logistic regression modeling adjusted for sex, and known CAD risk factors (history of hypertension, history of diabetes mellitus, body mass index, history of dyslipidemia, and smoking history) as covariates. These adjustments could hypothetically allow us to control for competing genetic pathways that are independent risk factors for CAD, thereby allowing us to detect a separate CAD genetic effect. SAS 9.1 (SAS Institute, Cary, North Carolina, United States) was used for statistical analysis.
LAMP Analysis in GENECARD ACS families and replication in CATHGEN MI cases and controls
LINKAGE TEST (df = 3)
ASSOCIATION TEST (df = 1)
OTHER LINKED VARIANTS (df = 2)
CATHGEN MI Logistic regression†
CATHGEN Minor Allele Frequency
Allele OR (95% CI)*
Aortic endothelial and smooth muscle cell culture
Cryopreserved human aortic smooth muscle (AoSMC) and human aortic endothelial cells (HAEC) were obtained from Lonza (Walkersville, MD) at passage 3 and cultured using SGM-2 or EGM-2 media (Lonza) until the indicated passage. Cells were grown in T-75 SoLo Flasks Nunclon™Δ (NUNC, Rochester, NY) according to manufacturer's instructions (see Availability and requirements section for URL).
Reverse transcriptase PCR (RT-PCR)
RNA from 20 human tissues (Human Total RNA Master Panel II) was obtained from Clontech (Mountain View, CA). Lot number 6020107 contained RNA from the following tissue types: Fetal Brain, Spinal Cord, Skeletal Muscle, Lung, Trachea, Fetal Liver, Liver, Bone Marrow, Kidney, Thymus, Thyroid, Adrenal Gland, Salivary Gland, Uterus, Testes, Prostate, and Placenta. RNA from cultured aortic smooth muscle and endothelial cells was isolated using a Ribopure Kit from Applied Biosystems (Foster City, CA). cDNA was synthesized using the protocol for cDNA synthesis in the Illumina Total Prep RNA Amplification kit from Applied Biosystems (Foster City, CA). Primers corresponding to exon 6 and exon 7 of FAM5C were used for PCR from cDNA (5'-ACAACAGTGACTTTGAGGAATCAGA-3' and 5'-GCGCTGAAAATTAGAATCCATTG-3'). Primers corresponding to exons 8 and 9 of GAPDH were used as a template control (5'-CTCCTCCACCTTTGACGCTG-3', and 5'-AGGGGAGATTCAGTGTGGTG-3'). 20 ng of cDNA template was used for each reaction and Platinum Taq polymerase (Invitrogen, Carlsbad, CA) was used for all PCR reactions. PCR reaction conditions were as follows: Step 1: 95°C for 5 minutes, Step 2: 30 cycles (FAM5C) or 20 cycles (GAPDH) of 95°C for 1 minute followed by 62°C for 1 minute, and 72°C for 1 min, and Step 3, 72°C 10 min.
FAM5C real time RT-PCR
Applied Biosystems Taqman Gene expression assays were used to perform quantitative (real time) RT-PCR (FAM5C, Hs00982332_m1 and GAPDH, Hs02758991_g1). The following reaction components were used for each probe: 2 uL cDNA, 5 ul Custom TaqMan SNP Master Mix (Applied Biosystems, Foster City, CA), 0.5 ul of assay, and 2.5 ul water. Reactions were performed in a single 384 well plate in triplicate using an ABI PRISM® 7900 HT Sequence Detection System. PCR reaction conditions were as follows: Step 1: 50°C for 2 minutes, Step 2: 95°C for 10 minutes, Step 3: 40 cycles of 95°C for 15 seconds followed by 60°C for 1 minute. Expression relative to GAPDH was calculated using 2ΔCt  and levels were normalized to heart expression.
Peak-wide genotyping and analysis in GENECARD ACS linkage Peak
457 SNPs, at an average spacing of 66 kilobases, were genotyped in the one LOD down support interval (169–199 Mb) of the GENECARD chromosome 1q ACS linkage peak (see Additional files, Table 1). Family-based association performed in the ACS families using APL identified sixteen SNPs that were significant at p ≤ 0.05 (Figure 1, Additional files, Table 3), three of which also displayed linkage: rs10912660 in the ankyrin repeat domain 45 gene (ANKRD45), rs12092963 in the ral guanine nucleotide dissociation gene (RGL1), and rs1891586 in FAM5C (Figure 1, bold). The clinical characteristics comparing the ACS probands in GENECARD to the MI cases in CATHGEN show the suitability of this group for replication (Table 1). Three of the sixteen SNPs replicated in the CATHGEN MI cohort at p ≤ 0.05 (Figure 1, underlined). Two SNPs, rs2239816 and rs1891586, reside within genes, chromosome 1 open reading frame 9 protein (C1orf9) and FAM5C, and the third within an intergenic region (rs1324713).
A SNP found within the first intron of the tumor necrosis factor superfamily 4 gene (TNFSF4), rs3850641, residing within the chromosome 1q linkage peak, has previously shown to be associated with MI in two independent populations , but although it was linked, was not associated in the GENECARD samples (Additional files, Table 4). LAMP analysis determined that this SNP does not account for the linkage observed in these families (data not shown).
tagSNP genotyping and analysis
We chose haplotype tagging SNPs in C1orf9, FAM5C and the intergenic region 100 kilobases upstream and downstream of rs1324713. None of the additional genotyped tagSNPs in C1orf9 or the intergenic region displayed both association and linkage (Additional files, Table 4). However, the original C1orf9 SNP identified in the screen, rs2239816, is significantly associated in both GENECARD ACS (APL, p = 0.05) and CATHGEN MI samples (allele p = 0.04) and results in an amino acid change at position 11 (Pro11Ser). This SNP does not, however, explain the linkage seen in our ACS family subset.
Sequencing FAM5C conserved regions in intron 7
We noted blocks of multi-species sequence conservation in the last intron of FAM5C and hypothesized that the causal SNP, in LD with rs10920501, may lie in one of these regions. The allele frequency of rs10920501 is approximately 16% in the Caucasian population. We rationalized that a causal SNP that was in LD with rs10920501 would be at a similar allele frequency. Given Kruglyak and Nickerson's estimates  we performed sequencing across 32 individuals that should allow for a SNP detection rate between 99% and 99.9% at a minor allele frequency of 10%. The seven conserved regions (see Additional files, Figure 1, arrows) were PCR amplified from 16 Caucasian cases and 16 Caucasian controls using standard conditions and sequenced. No novel SNPs were identified (data not shown).
FAM5C allele specific gene expression
FAM5C is expressed in proliferating human aortic smooth muscle cells and expression changes with passage
We noted that the expression levels of FAM5C appeared to vary with AoSMC passage (compare p.5 and p.8 AoSMC RT-PCR Figure 5 and Additional files, Figure 2) and therefore hypothesized that FAM5C levels may change with proliferative capacity and cellular senescence. Sub-culturing of primary cells, such as AoSMC, results in a progressive reduction of proliferation until the cells reach senescence, described as Hayflick's limit . We serially cultured AoSMC derived from a 29 year old Caucasian male donor until senescence (passage 20) indicated by cessation of proliferation and an enlarged flattened cell morphology. We measured the level of FAM5C expression in triplicate during each passage (passage 8 RNA was depleted) using real time RT-PCR. We found that FAM5C levels decrease as cells age in culture (Figure 5B).
Discussion and Conclusion
We herein report the results of fine-mapping a linkage peak for MI and identifying a novel MI susceptibility gene, FAM5C. We show that polymorphisms in FAM5C are associated with MI in two independent samples and that the SNP, which partially accounts for the linkage in our GENECARD ACS peak, is also significantly associated with changes in FAM5C gene expression in the human aorta. To begin to address the functional role of FAM5C in MI, we observed that FAM5C is expressed in the human aorta and that its transcript levels decrease with increasing passage of AoSMC, suggesting that the level of gene expression may play a role in proliferation and senescence of this cell type.
FAM5C was originally identified in the mouse brain as a gene that is induced by bone morphogenic protein and retinoic acid signaling (BRINP3 gene)  and has recently been shown to be over-expressed in human pituitary tumors . In the aforementioned study, the authors show that FAM5C is localized to the mitochondria and that over-expression of this molecule leads to increased proliferation, migration, and invasion of non-tumorogenic pituitary cells (the opposite being true for knockdown of FAM5C). Through complex signaling cascades, mitochondria have the ability to activate multiple pathways that modulate both cell proliferation and, conversely, promote cell arrest and programmed cell death (reviewed in [52, 53]). FAM5C localization to the mitochondria and its putative role in regulating cell proliferation and migration provide an intriguing hypothesis for the role of FAM5C in smooth muscle cells and the formation and vulnerability of the atherosclerotic plaque. SMCs move from a contractile state to a proliferative, migratory state in the presence of endothelial cell dysfunction at the initiation phase of plaque formation. They migrate into the region of plaque formation and form the fibrous cap that covers the plaque which is prone to rupture (reviewed in ). FAM5C levels could play a role in the initiation of smooth muscle cell proliferation and migration and/or in the disintegration of the smooth muscle cell in the fibrous cap. We are currently conducting experiments to determine if FAM5C levels have effects on smooth muscle cell proliferation, migration, apoptosis and senescence.
While our analysis was primarily aimed at an MI phenotype, our ACS family sample included a limited number of individuals with unstable angina. As a result the power to detect a distinct MI effect may be somewhat lower in our ACS family data set. Nevertheless, we observed consistent and replicated results in the MI case-control set and in the ACS family data set, suggesting that our results are robust to potential phenotypic heterogeneity. In addition, two recent WGA studies from the Framingham Heart Study have identified significant association between atrial fibrillation (AF) and FAM5C (rs1604355), and also for left atrial size (LAS) and FAM5C (rs1935881) [54, 55]. Both of these SNPs are in low LD (r2 < 0.1) with rs10920501; however, a SNP within the linkage disequilibrium bin containing rs10920501 was not represented on the SNP chip used in these studies, leaving open the possibility that rs10920501 may also play a role in these phenotypes. Although atrial fibrillation and left atrial size are very different phenotypes from MI, these data, together with the data presented in this paper, suggest that further genetic and functional studies are warranted to determine the role of FAM5C and variation within this gene in diverse diseases in cardiovascular tissues.
Three separate studies identified linkage to the same region of chromosome 1, the Framingham Heart Study , the IRAS family study , and the GENECARD study . Through careful fine mapping of a linkage peak for an extreme CAD phenotype, MI, we have defined a single gene that can partially account for the linkage in this peak. The relationship between MCP-1 levels, metabolic syndrome, MI, and FAM5C still remains to be defined. Inflammation lies at the heart of these pathophysiologies. We know that inflammation is an important part of metabolic syndrome , individuals with metabolic syndrome are at higher coronary disease risk , and metabolic syndrome increases the risk of an event (MI, revascularization, or cardiac death) in individuals with a family history of MI . In addition, vascular SMC pathophysiology is related to inflammation which could explain the MCP-1 QTL linkage. Further work is clearly indicated in this area.
One of the more difficult aspects of genetic association analysis is how to appropriately correct for multiple comparisons in evaluating the statistical significance of any given result. Possible approaches range from the most conservative Bonferroni correction, to estimation of false discovery rate , to weighted corrections of combined data, to no correction at all. We are aware that the p-values presented in this paper were not corrected for multiple comparisons. Our enthusiasm for FAM5C as an MI susceptibility gene is based on the observed replication of association results from multiple independent data sets and different approaches. The line of reasoning has several parts. First, rs10920501 represents a significant marker that can partially account for the linkage signal in our ACS families, and the statistical significance of the observation replicates in an independent case-control sample. Second, the probability of observing three or more significant results out of sixteen in the independent case-control sample under the null hypothesis that there is no association (p ≤ 0.05) is 0.043. Third, we have observed that the level of the transcript in the aorta is associated with genetic variation this SNP (or a SNP LD with it). Last, we observed that levels of FAM5C change with passage in AoSMC, implying that genetic variation in the gene might be involved in AoSMC cellular processes involved in the formation of unstable atherosclerotic plaque leading to MI. Further work to elucidate the function of FAM5C in MI, including the MI-associated polymorphisms and its downstream targets, is necessary in order to begin to understand the role this mitochondrial molecule plays in this disease.
Availability and requirements
We are appreciative of the subject volunteers' participation in the GENECARD and CATHGEN studies. We would also like to acknowledge the essential contributions of the following individuals to making this publication possible: Elaine Dowdy, the GENECARD Investigators Network, the CATHGEN Steering Committee Members, Margaret Pericak-Vance, Jeffery M. Vance, Pascal Goldschmidt-Clermont, Neil Freedman, Beth Sutton and the staff at the Center for Human Genetics for their innumerable contributions to this manuscript. Genotyping provided through the SeattleSNPs PGA Early Investigator Program (HL66682). This work was supported by NIH grants K99HL089412 (Connelly), HL073389 (Hauser), and HL73042 (Kraus).
- Rosamond W, Flegal K, Friday G, Furie K, Go A, Greenlund K, Haase N, Ho M, Howard V, Kissela B, Kittner S, Lloyd-Jones D, McDermott M, Meigs J, Moy C, Nichol G, O'Donnell CJ, Roger V, Rumsfeld J, Sorlie P, Steinberger J, Thom T, Wasserthiel-Smoller S, Hong Y, for the American Heart Association Statistics Committee and Stroke Statistics Subcommittee: Heart Disease and Stroke Statistics--2007 Update: A Report From the American Heart Association Statistics Committee and Stroke Statistics Subcommittee. Circulation. 2007, 115: e69-171. 10.1161/CIRCULATIONAHA.106.179918.View ArticlePubMedGoogle Scholar
- Organization WH: The World Health Report 2002 - Reducing Risks, Promoting Healthy Life. 2002Google Scholar
- Shea S, Ottman R, Gabrieli C, Stein Z, Nichols A: Family history as an independent risk factor for coronary artery disease. J Am Coll Cardiol. 1984, 4: 793-801.View ArticlePubMedGoogle Scholar
- Ciruzzi M, Schargrodsky H, Rozlosnik J, Pramparo P, Delmonte H, Rudich V, Piskorz D, Negri E, Soifer S, La Vecchia C: Frequency of family history of acute myocardial infarction in patients with acute myocardial infarction. Argentine FRICAS (Factores de Riesgo Coronario en America del Sur) Investigators. Am J Cardiol. 1997, 80: 122-127. 10.1016/S0002-9149(97)00304-4.View ArticlePubMedGoogle Scholar
- Rissanen AM: Familial occurrence of coronary heart disease: effect of age at diagnosis. Am J Cardiol. 1979, 44: 60-66. 10.1016/0002-9149(79)90251-0.View ArticlePubMedGoogle Scholar
- Helgadottir A, Manolescu A, Thorleifsson G, Gretarsdottir S, Jonsdottir H, Thorsteinsdottir U, Samani NJ, Gudmundsson G, Grant SF, Thorgeirsson G, Sveinbjornsdottir S, Valdimarsson EM, Matthiasson SE, Johannsson H, Gudmundsdottir O, Gurney ME, Sainz J, Thorhallsdottir M, Andresdottir M, Frigge ML, Topol EJ, Kong A, Gudnason V, Hakonarson H, Gulcher JR, Stefansson K: The gene encoding 5-lipoxygenase activating protein confers risk of myocardial infarction and stroke. Nat Genet. 2004, 36: 233-239. 10.1038/ng1311.View ArticlePubMedGoogle Scholar
- Helgadottir A, Manolescu A, Helgason A, Thorleifsson G, Thorsteinsdottir U, Gudbjartsson DF, Gretarsdottir S, Magnusson KP, Gudmundsson G, Hicks A, Jonsson T, Grant SF, Sainz J, O'Brien SJ, Sveinbjornsdottir S, Valdimarsson EM, Matthiasson SE, Levey AI, Abramson JL, Reilly MP, Vaccarino V, Wolfe ML, Gudnason V, Quyyumi AA, Topol EJ, Rader DJ, Thorgeirsson G, Gulcher JR, Hakonarson H, Kong A, Stefansson K: A variant of the gene encoding leukotriene A4 hydrolase confers ethnicity-specific risk of myocardial infarction 1. Nat Genet. 2006, 38: 68-74. 10.1038/ng1692.View ArticlePubMedGoogle Scholar
- Wang Q, Rao S, Shen GQ, Li L, Moliterno DJ, Newby LK, Rogers WJ, Cannata R, Zirzow E, Elston RC, Topol EJ: Premature Myocardial Infarction Novel Susceptibility Locus on Chromosome 1P34-36 Identified by Genomewide Linkage Analysis. Am J Hum Genet. 2004, 74: 262-271. 10.1086/381560.View ArticlePubMedPubMed CentralGoogle Scholar
- Shen GQ, Li L, Girelli D, Seidelmann SB, Rao S, Fan C, Park JE, Xi Q, Li J, Hu Y, Olivieri O, Marchant K, Barnard J, Corrocher R, Elston R, Cassano J, Henderson S, Hazen SL, Plow EF, Topol EJ, Wang QK: An LRP8 Variant Is Associated with Familial and Premature Coronary Artery Disease and Myocardial Infarction. Am J Hum Genet. 2007, 81: 780-791. 10.1086/521581.View ArticlePubMedPubMed CentralGoogle Scholar
- Helgadottir A, Thorleifsson G, Manolescu A, Gretarsdottir S, Blondal T, Jonasdottir G, Jonasdottir A, Sigurdsson A, Baker A, Palsson A, Masson G, Gudbjartsson D, Magnusson KP, Andersen K, Levey AI, Backman VM, Matthiasdottir S, Jonsdottir T, Palsson ST, H E, Gunnarsdottir S, Gylfason A, Vaccarino V, Hooper W, Reilly MP, CB G, Austin HA, Rader DJ, SH S, AA A, Gulcher J, Thorgeirsson G, Thorsteinsdottir U, Kong A, Stefansson K: A common variant on chromosome 9p21 affects the risk of myocardial infarction. Science. 2007.Google Scholar
- McPherson R, Pertsemlidis A, Kavaslar N, Stewart A, Roberts R, Cox DR, Hinds DA, Pennacchio LA, Tybjaerg-Hansen A, Folsom AR, Boerwinkle E, Hobbs HH, Cohen JC: A Common Allele on Chromosome 9 Associated with Coronary Heart Disease. Science. 2007Google Scholar
- Samani NJ, Erdmann J, Hall AS, Hengstenberg C, Mangino M, Mayer B, Dixon RJ, Meitinger T, Braund P, Wichmann HE, Barrett JH, Konig IR, Stevens SE, Szymczak S, Tregouet DA, Iles MM, Pahlke F, Pollard H, Lieb W, Cambien F, Fischer M, Ouwehand W, Blankenberg S, Balmforth AJ, Baessler A, Ball SG, Strom TM, Braenne I, Gieger C, Deloukas P, Tobin MD, Ziegler A, Thompson JR, Schunkert H: Genomewide association analysis of coronary artery disease. N Engl J Med. 2007, 357: 443-453. 10.1056/NEJMoa072366.View ArticlePubMedPubMed CentralGoogle Scholar
- Consortium WTCC: Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls. Nature. 2007, 447: 661-678. 10.1038/nature05911.View ArticleGoogle Scholar
- Broadbent HM, Peden JF, Lorkowski S, Goel A, Ongen H, Green F, Clarke R, Collins R, Franzosi MG, Tognoni G, Seedorf U, Rust S, Eriksson P, Hamsten A, Farrall M, Watkins H: Susceptibility to coronary artery disease and diabetes is encoded by distinct, tightly linked SNPs in the ANRIL locus on chromosome 9p. Hum Mol Genet. 2008, 17: 806-814. 10.1093/hmg/ddm352.View ArticlePubMedGoogle Scholar
- Hauser ER, Crossman DC, Granger CB, Haines JL, Jones CJ, Mooser V, McAdam B, Winkelmann BR, Wiseman AH, Muhlestein JB, Bartel AG, Dennis CA, Dowdy E, Estabrooks S, Eggleston K, Francis S, Roche K, Clevenger PW, Huang L, Pedersen B, Shah S, Schmidt S, Haynes C, West S, Asper D, Booze M, Sharma S, Sundseth S, Middleton L, Roses AD, Hauser MA, Vance JM, Pericak-Vance MA, Kraus WE: A genomewide scan for early-onset coronary artery disease in 438 families: the GENECARD Study. Am J Hum Genet. 2004, 75: 436-447. 10.1086/423900.View ArticlePubMedPubMed CentralGoogle Scholar
- Dupuis J, Larson MG, Vasan RS, Massaro JM, Wilson PW, Lipinska I, Corey D, Vita JA, Keaney JF, Benjamin EJ: Genome scan of systemic biomarkers of vascular inflammation in the Framingham Heart Study: evidence for susceptibility loci on 1q. Atherosclerosis. 2005, 182: 307-314. 10.1016/j.atherosclerosis.2005.02.015.View ArticlePubMedGoogle Scholar
- Langefeld CD, Wagenknecht LE, Rotter JI, Williams AH, Hokanson JE, Saad MF, Bowden DW, Haffner S, Norris JM, Rich SS, Mitchell BD: Linkage of the metabolic syndrome to 1q23-q31 in Hispanic families: the Insulin Resistance Atherosclerosis Study Family Study. Diabetes. 2004, 53: 1170-1174. 10.2337/diabetes.53.4.1170.View ArticlePubMedGoogle Scholar
- Porreca E, DiFebbo C, Reale M, Castellani ML, Baccante G, Barbacane R, Conti P, Cuccurullo F, Poggi A: Monocyte chemotactic protein 1 (MCP-1) is a mitogen for cultured rat vascular smooth muscle cells. Journal of Vascular Research. 1997, 34: 58-65.PubMedGoogle Scholar
- Ford ES: The metabolic syndrome and mortality from cardiovascular disease and all-causes: findings from the National Health and Nutrition Examination Survey II Mortality Study. Atherosclerosis. 2004, 173: 309-314. 10.1016/j.atherosclerosis.2003.12.022.View ArticlePubMedGoogle Scholar
- Pyorala M, Miettinen H, Halonen P, Laakso M, Pyorala K: Insulin resistance syndrome predicts the risk of coronary heart disease and stroke in healthy middle-aged men - The 22-year follow-up results of the Helsinki Policemen Study. Arteriosclerosis Thrombosis and Vascular Biology. 2000, 20: 538-544.View ArticleGoogle Scholar
- Fan AZ: Metabolic syndrome and progression of atherosclerosis among middle-aged US adults. Journal of Atherosclerosis and Thrombosis. 2006, 13: 46-54.View ArticlePubMedGoogle Scholar
- Connelly JJ, Wang T, Cox JE, Haynes C, Wang L, Shah SH, Crosslin DR, Hale AB, Nelson S, Crossman DC, Granger CB, Haines JL, Jones CJ, Vance JM, Goldschmidt-Clermont PJ, Kraus WE, Hauser ER, Gregory SG: GATA2 Is Associated with Familial Early-Onset Coronary Artery Disease. PLoS Genet. 2006, 2:Google Scholar
- Wang L, Hauser ER, Shah SH, Pericak-Vance MA, Haynes C, Crosslin D, Harris M, Nelson S, Hale AB, Granger CB, Haines JL, Jones CJ, Crossman D, Seo D, Gregory SG, Kraus WE, Goldschmidt-Clermont PJ, Vance JM: Peakwide mapping on chromosome 3q13 identifies the kalirin gene as a novel candidate gene for coronary artery disease. Am J Hum Genet. 2007, 80: 650-663. 10.1086/512981.View ArticlePubMedPubMed CentralGoogle Scholar
- Shorts-Cary L, Xu M, Ertel J, Kleinschmidt-DeMasters BK, Lillehei K, Matsuoka I, Nielsen-Preiss S, Wierman ME: Bone morphogenetic protein and retinoic acid-inducible neural specific protein-3 is expressed in gonadotrope cell pituitary adenomas and induces proliferation, migration, and invasion. Endocrinology. 2007, 148: 967-975. 10.1210/en.2006-0905.View ArticlePubMedGoogle Scholar
- Watkins H, Farrall M: Genetic susceptibility to coronary artery disease: from promise to progress. Nat Rev Genet. 2006, 7: 163-173. 10.1038/nrg1805.View ArticlePubMedGoogle Scholar
- Clarke MCH, Figg N, Maguire JJ, Davenport AP, Goddard M, Littlewood TD, Bennett MR: Apoptosis of vascular smooth muscle cells induces features of plaque vulnerability in atherosclerosis. Nature Medicine. 2006, 12: 1075-1080. 10.1038/nm1459.View ArticlePubMedGoogle Scholar
- Hauser ER, Mooser V, Crossman DC, Haines JL, Jones CH, Winkelmann BR, Schmidt S, Scott WK, Roses AD, Pericak-Vance MA, Granger CB, Kraus WE: Design of the genetics of early onset cardiovascular disease (GENECARD) study. Am Heart J. 2003, 145: 602-613. 10.1067/mhj.2003.13.View ArticlePubMedGoogle Scholar
- Fortin DF, Califf RM, Pryor DB, Mark DB: The way of the future redux. Am J Cardiol. 1995, 76: 1177-1182. 10.1016/S0002-9149(99)80331-2.View ArticlePubMedGoogle Scholar
- Smith LR, Harrell FE, Rankin JS, Califf RM, Pryor DB, Muhlbaier LH, Lee KL, Mark DB, Jones RH, Oldham HN, Glower DD, Reves JG, Sabiston DC: Determinants of Early Versus Late Cardiac Death in Patients Undergoing Coronary-Artery Bypass Graft-Surgery. Circulation. 1991, 84: 245-253.Google Scholar
- Kong DF, Shaw LK, Harrell FE, Muhlbaier LH, Lee KL, Califf RM, Jones RH: Predicting survival from the coronary arteriogram: an experience-based statistical index of coronary artery disease severity. Journal of the American College of Cardiology. 2002, 39(Suppl A): 327A-10.1016/S0735-1097(02)81469-3.View ArticleGoogle Scholar
- Felker GM, Shaw LK, O'Connor CM: A standardized definition of ischemic cardiomyopathy for use in clinical research. J Am Coll Cardiol. 2002, 39: 210-218. 10.1016/S0735-1097(01)01738-7.View ArticlePubMedGoogle Scholar
- Luo AK, Jefferson BK, Garcia MJ, Ginsburg GS, Topol EJ: Challenges in the phenotypic characterisation of patients in genetic studies of coronary artery disease. Journal of Medical Genetics. 2007, 44: 161-165. 10.1136/jmg.2006.045732.View ArticlePubMedGoogle Scholar
- Seo D, Wang T, Dressman H, Herderick EE, Iversen ES, Dong C, Vata K, Milano CA, Rigat F, Pittman J, Nevins JR, West M, Goldschmidt-Clermont PJ: Gene Expression Phenotypes of Atherosclerosis. Arterioscler Thromb Vasc Biol. 2004, 24: 1922-1927. 10.1161/01.ATV.0000141358.65242.1f.View ArticlePubMedGoogle Scholar
- Oliveira S, Li YJ, Noureddine M, Zuchner S, Qin X, Pericak-Vance MA, Vance JM: Identification of risk and age-at-onset genes on chromosome 1p in Parkinson disease. Am J Hum Genet. 2005, 77: 252-264. 10.1086/432588.View ArticlePubMedPubMed CentralGoogle Scholar
- Abecasis GR, Cookson WO: GOLD--graphical overview of linkage disequilibrium. BioInformatics. 2000, 16: 182-183. 10.1093/bioinformatics/16.2.182.View ArticlePubMedGoogle Scholar
- Barrett JC, Fry B, Maller J, Daly MJ: Haploview: analysis and visualization of LD and haplotype maps. BioInformatics. 2005, 21: 263-265. 10.1093/bioinformatics/bth457.View ArticlePubMedGoogle Scholar
- Cottingham RW, Idury RM, Schaffer AA: Faster sequential genetic linkage computations. Am J Hum Genet. 1993, 53: 252-263.PubMedPubMed CentralGoogle Scholar
- Schaffer AA, Gupta SK, Shriram K, Cottingham RW: Avoiding recomputation in linkage analysis. Hum Hered. 1994, 44: 225-237.View ArticlePubMedGoogle Scholar
- Ott J: Analysis of Human Genetic Linkage. Third Edition. 1999, Baltimore, Maryland, The Johns Hopkins University Press, 3Google Scholar
- Martin ER, Bass MP, Hauser ER, Kaplan NL: Accounting for linkage in family-based tests of association with missing parental genotypes. Am J Hum Genet. 2003, 73: 1016-1026. 10.1086/378779.View ArticlePubMedPubMed CentralGoogle Scholar
- Chung RH, Hauser ER, Martin ER: The APL Test: Extension to General Nuclear Families and Haplotypes and Examination of Its Robustness. Hum Hered. 2006, 61: 189-199. 10.1159/000094774.View ArticlePubMedGoogle Scholar
- Martin ER, Monks SA, Warren LL, Kaplan NL: A test for linkage and association in general pedigrees: the pedigree disequilibrium test. Am J Hum Genet. 2000, 67: 146-154. 10.1086/302957.View ArticlePubMedPubMed CentralGoogle Scholar
- Martin ER, Bass MP, Gilbert JR, Pericak-Vance MA, Hauser ER: Genotype-based association test for general pedigrees: the genotype-PDT. Genet Epidemiol. 2003, 25: 203-213. 10.1002/gepi.10258.View ArticlePubMedGoogle Scholar
- Li M, Boehnke M, Abecasis GR: Joint modeling of linkage and association: identifying SNPs responsible for a linkage signal. Am J Hum Genet. 2005, 76: 934-949. 10.1086/430277.View ArticlePubMedPubMed CentralGoogle Scholar
- Li M, Boehnke M, Abecasis GR: Efficient study designs for test of genetic association using sibship data and unrelated cases and controls. Am J Hum Genet. 2006, 78: 778-792. 10.1086/503711.View ArticlePubMedPubMed CentralGoogle Scholar
- Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.View ArticlePubMedGoogle Scholar
- Wang X, Ria M, Kelmenson PM, Eriksson P, Higgins DC, Samnegard A, Petros C, Rollins J, Bennet AM, Wiman B, De Faire U, Wennberg C, Olsson PG, Ishii N, Sugamura K, Hamsten A, Forsman-Semb K, Lagercrantz J, Paigen B: Positional identification of TNFSF4, encoding OX40 ligand, as a gene that influences atherosclerosis susceptibility. Nat Genet. 2005, 37: 365-372. 10.1038/ng1524.View ArticlePubMedGoogle Scholar
- de Bakker PI, Yelensky R, Pe'er I, Gabriel SB, Daly MJ, Altshuler D: Efficiency and power in genetic association studies. Nat Genet. 2005, 37: 1217-1223. 10.1038/ng1669.View ArticlePubMedGoogle Scholar
- Kruglyak L, Nickerson DA: Variation is the spice of life. Nat Genet. 2001, 27: 234-236. 10.1038/85776.View ArticlePubMedGoogle Scholar
- Hayflick L, Moorhead PS: Serial Cultivation of Human Diploid Cell Strains. Experimental Cell Research. 1961, 25: 585-&. 10.1016/0014-4827(61)90192-6.View ArticlePubMedGoogle Scholar
- Kawano H, Nakatani T, Mori T, Ueno S, Fukaya M, Abe A, Kobayashi M, Toda F, Watanabe M, Matsuoka I: Identification and characterization of novel developmentally regulated neural-specific proteins, BRINP family. Molecular Brain Research. 2004, 125: 60-75. 10.1016/j.molbrainres.2004.04.001.View ArticlePubMedGoogle Scholar
- McBride HM, Neuspiel M, Wasiak S: Mitochondria: More than just a powerhouse. Current Biology. 2006, 16: R551-R560. 10.1016/j.cub.2006.06.054.View ArticlePubMedGoogle Scholar
- Rustin P: Mitochondria, from cell death to proliferation. Nature Genetics. 2002, 30: 352-353. 10.1038/ng0402-352.View ArticlePubMedGoogle Scholar
- Larson MG, Atwood LD, Benjamin EJ, Cupples LA, D'Agostino RB, Fox CS, Govindaraju DR, Guo CY, Heard-Costa NL, Hwang SJ, Murabito JM, Newton-Cheh C, O'Donnell CJ, Seshadri S, Vasan RS, Wang TJ, Wolf PA, Levy D: Framingham Heart Study 100K project: genome-wide associations for cardiovascular disease outcomes. Bmc Medical Genetics. 2007, 8:Google Scholar
- Vasan RS, Larson MG, Aragam J, Wang TJ, Mitchell GF, Kathiresan S, Newton-Cheh C, Vita JA, Keyes MJ, O'Donnell CJ, Levy D, Benjamin EJ: Genome-wide association of echocardiographic dimensions, brachial artery endothelial function and treadmill exercise responses in the Framingham Heart Study. Bmc Medical Genetics. 2007, 8:Google Scholar
- Gustafson B, Hammarstedt A, Andersson CX, Smith U: Inflamed adipose tissue - A culprit underlying the metabolic syndrome and atherosclerosis. Arteriosclerosis Thrombosis and Vascular Biology. 2007, 27: 2276-2283. 10.1161/ATVBAHA.107.147835.View ArticleGoogle Scholar
- Espinola-Klein C, Rupprecht HJ, Bickel C, Post F, Genth-Zotz S, Lackner K, Munzel T, Blankenberg S: Impact of metabolic syndrome on atherosclerotic burden and cardiovascular prognosis. American Journal of Cardiology. 2007, 99: 1623-1628. 10.1016/j.amjcard.2007.01.049.View ArticlePubMedGoogle Scholar
- Reinhard W, Holmer SR, Fischer M, Gloeckner C, Hubauer U, Baessler A, Mayer B, Schunkert H, Riegger GA, Hengstenberg C: Association of the metabolic syndrome with early coronary disease in families with frequent myocardial infarction. Am J Cardiol. 2006, 97: 964-967. 10.1016/j.amjcard.2005.10.063.View ArticlePubMedGoogle Scholar
- Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc B. 1995, 57: 289-300.Google Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2350/9/33/prepub
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/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.