Joint study of genetic regulators for expression traits related to breast cancer
© Zheng et al; licensee BioMed Central Ltd. 2007
Published: 18 December 2007
The mRNA expression levels of genes have been shown to have discriminating power for the classification of breast cancer. Studying the heritability of gene expression levels on breast cancer related transcripts can lead to the identification of shared common regulators and inter-regulation patterns, which would be important for dissecting the etiology of breast cancer.
We applied multilocus association genome-wide scans to 18 breast cancer related transcripts and combined the results with traditional linkage scans. Regulatory hotspots for these transcripts were identified and some inter-regulation patterns were observed. We also derived evidence on interacting genetic regulatory loci shared by a number of these transcripts.
In this paper, by restricting to a set of related genes, we were able to employ a more detailed multilocus approach that evaluates both marginal and interaction association signals at each single-nucleotide polymorphism. Interesting inter-regulation patterns and significant overlaps of genetic regulators between transcripts were observed. Interaction association results returned more expression quantitative trait locus hotspots that are significant.
Breast cancer (MIM 114480) is a common and genetically heterogeneous human disorder . Many studies have shown that gene expressions possess discriminating power for the diagnosis of breast cancer (e.g., van't Veer et al. ). Therefore, studying the genetic regulators of breast cancer related gene expression transcripts would shed light on the genetic mechanism of this disorder. Currently, the search for genetic regulators of expression traits in humans has been underway through linkage or association scans [3–5]. Morley et al.  measured the expression levels of 8500 transcripts using the Affymetrix Human Focus Arrays on 194 individuals in 14 Centre d'Etude du Polymorphisme Humain (CEPH) families. 3554 of these transcripts were identified to have greater between-subject variation than within-subject variation and were then used for linkage analysis. For the linkage analysis, genotypes of these CEPH individuals on 2882 single-nucleotide polymorphisms (SNPs) across the genome were obtained from The SNP Consortium http://snp.cshl.org/linkage%20maps/.
In this paper, we study transcripts related to breast cancer susceptibility genes. Among the 3554 selected transcripts in , we identified 18 transcripts that are related to seven genes listed in the overview of breast cancer from OMIM (the Online Mendelian Inheritance in Man). Treating the expression levels of these genes as hereditable traits, we used multilocus association genome scans, combined with linkage scans, to study their regulation patterns and identify joint features and inter-relations. Gene expressions are likely to be complex traits that are regulated by multiple genetic factors. Therefore, we used an extension of our previous multilocus methods [6–8] to extract more association information.
where S i is the trait's rank sum on the n i individuals with genotype i, and n i (n + 1)/2 is the expected value of S i under the null hypothesis that these SNPs are not associated with the phenotype.
qGTD captures the differences between the observed rank sums and those under the null hypothesis. The magnitude of qGTD scores reflects the level of association with the phenotype: the greater the value, the stronger the association . Unassociated SNPs add dimensions to the multilocus genotypes and lower the value of qGTD. Therefore, a greedy screening algorithm is used to screen out SNPs that do not contribute to increase the value of qGTD and retain a cluster of SNPs that contribute important information to the score. As discussed in [6–8], such a screening is not informative for a large number of SNPs simultaneously due to sparseness in high dimensions. A random subspace strategy is then employed, where the greedy algorithm is repeated on a large number of random SNP subsets. SNPs are then ranked by the numbers of times (return frequencies) that they are retained by the screening algorithm, which measure the overall importance of individual SNPs.
To evaluate the importance of the SNPs in gene × gene interactions, we further filtered the retained SNP clusters from the qGTD screening by their qGTD scores and only selected the top 1000 distinctive clusters with the highest qGTD values. Using these 1000 clusters, we computed the qGTD return frequencies for each SNP. As discussed previously, higher value of qGTD indicates stronger joint effects from the SNPs on the quantitative phenotype. SNPs that present more frequently in clusters with higher qGTD values play a more critical role in gene × gene interactions that decide the variation of the phenotype.
In this paper, we apply the above association scan (repeated on 5 million random subsets) using the 56 unrelated grandparents in the 14 CEPH families. For each selected expression trait, we selected the top 30 overall important SNPs with the highest return frequencies, and the top 30 important interaction SNPs with the highest qGTD return frequencies, which give us a comparable number of identified loci as that by the linkage scans.
Linkage analysis was done on all 194 members of 14 CEPH families using the pedigree analysis package MERLIN . The command pedwipe was first used to remove unlikely genotypes in the pedigree data. Regression-based linkage analysis for quantitative traits proposed by Sham et al.  was applied to all 18 expression traits with estimated mean, variance and heritability. The original data only contain physical map. In our analysis, we used linkage map provided by Sung et al. .
Clustering of transcripts based on identified regulatory loci
where m ij is the number of shared interacting loci pairs between transcripts i and j, m i and m j are the total numbers of interacting loci pairs for i and j, respectively, and m is the total number of loci pairs on the genome. We also clustered the transcripts based on their gene expression values with dissimilarity being one minus the correlation.
Results and discussion
Combined association and linkage scans
Genome scan results for transcripts are arranged in the rows of Figure 1. Several interesting patterns are observed. First, association signals frequently cluster with linkage signals. Actually, similar to Roeder et al. , we may use linkage signals to control for false LD signals. However, since we only used single-locus analysis in the linkage scans, some loci may fail to have linkage signals if they are in interactions when deciding the traits' variation and thus have lower marginal signals. Second, both linkage and association signals show overlaps between transcripts that cannot be explained by chance, which is discussed in the next section. The performance of the association scan would also be greatly improved if denser SNP data were available.
Transcriptional hotspots may be observed even when the genotypic data are not linked to or associated with the phenotypes, due to reasons such as highly correlated phenotypes [17, 18]. We calculated the correlation coefficients between the 18 transcripts studied and discovered that the highest value was 0.65. According to simulations, it is less likely to have false hotspots for traits with correlations at this level . Therefore, what we have observed in this paper is more likely due to true regulatory activities.
cis-Acting versus trans-acting
From Figure 1, not much evidence on cis-acting regulators has been identified except for some association signals and for the locus of BRCA1, which may require further investigations. Most interestingly, TP53BP1, a p53 binding protein, is shown to have a genetic regulator near the locus of BRCA1, a breast cancer susceptibility gene. It had been previously discovered by Rauch et al.  that 53BP1 binds to a promoter region of BRCA1 (MIM 113705).
Clustering of transcripts
In this paper, we carried out a detailed joint study on 18 breast cancer related transcripts, using both linkage and association scans. Interesting inter-regulation patterns and significant overlaps of genetic regulators between transcripts were observed. Quantitative backward genotype-trait association (qBGTA), as a nonparametric multilocus association approach, studies eQTL without assuming a trait model while considering interactions . Using qBGTA, we evaluated both marginal and interaction association signals at each SNPs locus and results on interaction association returned more significant eQTL hotspots (Figure 2).
This research is supported in part by NIH grant R01 GM070789.
This article has been published as part of BMC Proceedings Volume 1 Supplement 1, 2007: Genetic Analysis Workshop 15: Gene Expression Analysis and Approaches to Detecting Multiple Functional Loci. The full contents of the supplement are available online at http://www.biomedcentral.com/1753-6561/1?issue=S1.
- Nathanson KI, Weber BL: Breast Cancer. The Genetic Basis of Common Diseases. Edited by: King RA, Rotter JI, Motulsky AG. 2002, New York: Oxford University Press, 670-699. 2Google Scholar
- van't Veer LJ, Dai H, van de Vijver MJ, He YD, Hart AA, Mao M, Peterse HL, van der Kooy K, Marton MJ, Witteveen AT, Schreiber GJ, Kerkhoven RM, Roberts C, Linsley PS, Bernards R, Friend SH: Gene expression profiling predicts clinical outcome of breast cancer. Nature. 2002, 415: 530-536. 10.1038/415530a.View ArticleGoogle Scholar
- Cheung VG, Conlin LK, Weber TM, Arcaro M, Jen KY, Morley M, Spielman RS: Natural variation in human gene expression assessed in lymphoblastoid cells. Nat Genet. 2003, 33: 422-425. 10.1038/ng1094.View ArticlePubMedGoogle Scholar
- Morley M, Molony CM, Weber TM, Devlin JL, Ewens KG, Spielman RS, Cheung VG: Genetic analysis of genome-wide variation in human gene expression. Nature. 2004, 430: 743-747. 10.1038/nature02797.View ArticlePubMed CentralPubMedGoogle Scholar
- Cheung VG, Spielman RS, Ewens KG, Weber TM, Morley M, Burdick JT: Mapping determinants of human gene expression by regional and genome-wide association. Nature. 2005, 437: 1365-1369. 10.1038/nature04244.View ArticlePubMed CentralPubMedGoogle Scholar
- Lo SH, Zheng T: Backward haplotype transmission association (BHTA) algorithm – a fast multiple-marker screening method. Hum Hered. 2002, 53: 197-215. 10.1159/000066194.View ArticlePubMedGoogle Scholar
- Lo SH, Zheng T: A demonstration and findings of a statistical approach through reanalysis of inflammatory bowel disease data. Proc Natl Acad Sci USA. 2004, 101: 10386-10391. 10.1073/pnas.0403662101.View ArticlePubMed CentralPubMedGoogle Scholar
- Zheng T, Wang H, Lo SH: Backward genotype-trait association (BGTA) based dissection of complex traits in case-control designs. Hum Hered. 2006, 62: 196-212. 10.1159/000096995.View ArticlePubMed CentralPubMedGoogle Scholar
- Scheet P, Stephens M: A fast and flexible statistical model for large-scale population genotype data: applications to inferring missing genotypes and haplotypic phase. Am J Hum Genet. 2006, 78: 629-644. 10.1086/502802.View ArticlePubMed CentralPubMedGoogle Scholar
- Wang H: A new approach of detecting influential SNPs for complex phenotypes with genotype data. PhD dissertation. 2005, Columbia University, Department of Statistics, [http://digitalcommons.libraries.columbia.edu/dissertations/AAI3183007/]Google Scholar
- Abecasis GR, Cherny SS, Cookson WO, Cardon LR: Merlin-rapid analysis of dense genetic maps using sparse gene flow trees. Nat Genet. 2006, 30: 97-101. 10.1038/ng786.View ArticleGoogle Scholar
- Sham PC, Purcell S, Cherny SS, Abecasis GR: Powerful regression-based quantitative-trait linkage analysis of general pedigrees. Am J Hum Genet. 2002, 71: 238-253. 10.1086/341560.View ArticlePubMed CentralPubMedGoogle Scholar
- Sung YJ, Di Y, Fu AQ, Rothstein JH, Sieh W, Tong L, Thompson EA, Wijsman EM: Comparison of multipoint linkage analyses for quantitative traits: parametric LOD scores, variance components LOD scores, and Bayes factors in the CEPH data. BMC Proc. 2007, 1 (Suppl 1): S93-View ArticlePubMed CentralPubMedGoogle Scholar
- Gordon AD: Classification. 1999, Boca Raton, FL: CRC Press, 2Google Scholar
- Roeder K, Bacanu SA, Wasserman L, Devlin B: Using linkage genome scans to improve power of association in genome scans. Am J Hum Genet. 2006, 78: 243-252. 10.1086/500026.View ArticlePubMed CentralPubMedGoogle Scholar
- Rauch T, Zhong X, Pfeifer GP, Xu X: 53BP1 is a positive regulator of the BRCA1 promoter. Cell Cycle. 2005, 4: 1078-1083.View ArticlePubMedGoogle Scholar
- Perez-Enciso M: In silico study of transcriptome genetic variation in outbred populations. Genetics. 2004, 166: 547-554. 10.1534/genetics.166.1.547.View ArticlePubMed CentralPubMedGoogle Scholar
- Wang S, Zheng T, Wang Y: Transcription activity hot spot, is it real or an artifact?. BMC Proc. 2007, 1 (Suppl 1): S94-View ArticlePubMed CentralPubMedGoogle Scholar
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.