Volume 1 Supplement 1

## Genetic Analysis Workshop 15: Gene Expression Analysis and Approaches to Detecting Multiple Functional Loci

# A method dealing with a large number of correlated traits in a linkage genome scan

- Tao Feng
^{1}, - Shuanglin Zhang
^{1, 2}and - Qiuying Sha
^{1}Email author

**1(Suppl 1)**:S84

https://doi.org/10.1186/1753-6561-1-S1-S84

© Feng et al; licensee BioMed Central Ltd. 2007

**Published: **18 December 2007

## Abstract

We propose a method to perform linkage genome scans for many correlated traits in the Genetic Analysis Workshop 15 (GAW15) data. The proposed method has two steps: first, we use a clustering method to find the tight clusters of the traits and use the first principal component (PC) of the traits in each cluster to represent the cluster; second, we perform a linkage scan for each cluster by using the representative trait of the cluster. The results of applying the method to the GAW15 Problem 1 data indicate that most of the traits in the same cluster have the same regulators, and the representative trait measure, the first PC, can explain a large part of the total variation of all the traits in each cluster. Furthermore, considering one cluster of traits at a time may yield more linkage signals than considering traits individually.

## Background

From yeast to humans, the expression level of many genes shows familial aggregation and a simple segregation pattern in [1–3], suggesting an inherited contribution. Morley et al. [3] used microarray to measure the baseline expression level of genes in immortalized B cell from members of 14 Centre d'Etude du Polymorphism Humain (CEPH) pedigrees. They mapped the quantitative trait loci (QTL) of many quantitative traits (gene expression phenotypes) to chromosomal locations using genome scans. The data set of Morley et al. was made available to us by Genetic Analysis Workshop 15 (GAW15). This data set provides genome scan data for 3554 quantitative traits, or expression phenotypes of 3554 genes.

Morley et al. carried out a linkage genome scan for each of the 3554 traits individually and found that 142 and 984 of the 3554 traits have at least one significant linkage signal when the point-wise *p*-value thresholds 4.3 × 10^{-7} and 3.7 × 10^{-5} were used, respectively. Considering the correlated nature of the 3554 traits, we proposed a method to find clusters of genes whose expressions appear to be highly correlated and to localize the genetic determinants of each of the clusters. The proposed method has two steps: 1) use a clustering method to find tight clusters of the traits, and in each cluster, use the first principal component (PC) of the traits to represent this cluster; and 2) carry out a linkage scan using the first PC trait for each cluster and consider the linkage evidence of the first PC trait as the linkage evidence of all the traits in the cluster. Upon applying the method to the GAW15 Problem 1, we found that: 1) the traits in the same cluster have linkage signals at the same chromosomal regions; 2) the first PC can explain more than 53% of the total variation of all the traits in each cluster, which indicates that it is reasonable to represent a tight cluster of traits by using the first PC; and 3) this clustering of highly correlated traits may show more linkage signals than considering the traits individually.

## Methods

### Tight clustering

Tight clustering, first proposed by Tseng and Wong [4], is a method that produces tight and stable clusters without forcing all points into clusters. Tseng and Wong used a resample-based algorithm. Briefly, in each resample, the *K*-mean clustering algorithm is used to group the sampled genes into clusters. A similarity measure between two genes is defined as the number of resamples in which the two genes are in the same cluster. Then, this similarity measure is used to define a tight cluster. Here, to find a cluster of highly correlated genes, we propose a new tight clustering algorithm, which is computationally simpler than Tseng and Wong's algorithm.

Suppose there are *n* individuals and each individual has *m* traits (expressions of *m* genes). Let *x*_{
ij
}denote the trait value of the *j*^{th} trait of the *i*^{th} individual, and *x*_{
j
}= (*x*_{1j},..., *x*_{
nj
}) denote the values of the *j*^{th} trait of all individuals. A set of traits is said to be a tight cluster if the average correlation coefficient between traits in this set is larger than a given threshold *α*, and the number of traits in this set is between *m*_{1} and *m*_{2}. Briefly, for given parameters *α*_{0}, *m*_{1} and *m*_{2}, our "tight clustering" algorithm involves the following steps: 1) for a given *α* (≥ *α*_{0}), find the largest tight cluster. If a tight cluster for the given *α* cannot be found, reduce the value of *α*. If a tight cluster still cannot be found when *α* = *α*_{0}, stop the algorithm. 2) Remove the traits in the tight cluster from the data set and repeat Step 1.

For a given data set and the values of the parameters *α*_{0}, *m*_{1} and *m*_{2}, it is a computational challenge to search the whole state space to find the largest tight cluster. We propose to use a backwards elimination approach to find an approximate solution. Let *ρ*_{
j
}denote the average correlation coefficient between the *j*^{th} trait and all the other traits. The backward elimination has the following three steps: 1) rank the traits by *ρ*_{
j
}from the largest to the smallest and delete a proportion (5% in our implementation) of traits with the smallest values of *ρ*_{
j
}. 2) Recalculate *ρ*_{
j
}in the new data set (5% of traits were deleted) and repeat Step 1. 3) Repeat Steps 1 and 2 until the average correlation coefficient is larger than *α* and the number of traits is between *m*_{1} and *m*_{2}.

In our implementation, we used *m*_{1} = 20, *m*_{2} = 200, and *α*_{0} = 0.5 (see Discussion section for the effect of different parameter values). We first chose *α* = 0.8, and then reduced *α* to 0.7, 0.6, and 0.5. We then found the first PC of each cluster, which, because traits within are highly correlated, should explain a large part of the total variation.

### Linkage genome scan

We used the program *Merlin-regress* [5] to perform the linkage genome scan for each quantitative trait. *Merlin-regress* implements a method proposed by Sham et al. [6], which is based on a regression of estimated identity-by-descent (IBD) sharing between relative pairs on the squared sums and squared differences of trait values of the relative pairs. This program allows quick computation and has similar power of variance-component approaches.

*Merlin-regress*produces a

*p*-value for each trait-marker pair. A linkage signal is declared if the

*p*-value is less than a cut-off value. We propose to obtain a cut-off value of individual

*p*-values by controlling the false-discovery rate (FDR). Suppose we have

*M*trait-marker pairs and denote the ordered

*p*-values by

*p*

_{(1)},...,

*p*

_{(M)}. To control the FDR at 5%, the cut-off value of the individual

*p*-value will be as described by Benjamini and Yekutieli [7]:

## Results

We applied the proposed method to the data set given by the GAW15 Problem 1. The data set contained members of 14 CEPH Utah families. Each individual had 3554 traits and genotypes at 2882 SNPs across the genome. We deleted some SNPs with genotypes that showed Mendelian inconsistency. There were 2761 SNPs left for further analysis.

The 18 clusters found by the tight clustering method

Cluster ID | # of traits | Average correlation |
---|---|---|

1 | 21 | 0.8 |

2 | 43 | 0.7 |

3 | 32 | 0.7 |

4 | 27 | 0.7 |

5 | 28 | 0.7 |

6 | 35 | 0.6 |

7 | 35 | 0.6 |

8 | 47 | 0.6 |

9 | 33 | 0.6 |

10 | 27 | 0.6 |

11 | 64 | 0.5 |

12 | 42 | 0.5 |

13 | 84 | 0.5 |

14 | 59 | 0.5 |

15 | 44 | 0.5 |

16 | 41 | 0.5 |

17 | 37 | 0.5 |

18 | 44 | 0.5 |

We carried out a linkage scan for each of the 18 first PC traits by using the program *Merlin-regress*. For comparison, we also performed a linkage scan for each of the 743 traits in the 18 clusters. For a linkage genome scan, the program *Merlin-regress* gave *p*-values of the 2761 SNPs. We used the Eq. (1) to find the cut-off value *β*_{0} of the individual *p*-values to control the FDR at 5%. For the 18 PC traits, we got *M* = 18 × 2761 trait-marker pairs and the cut-off value *β*_{0} = 1.67 × 10^{-5}. For the 743 individual traits in the 18 clusters, there were *M* = 743 × 2761 trait-marker pairs, and the cut-off value was *β*_{0} = 5.34 × 10^{-8}.

*p*-value less than the cut-off value) are considered to be regulators. The number of regulators of each of the 18 PC traits and the number of traits with at least one regulator in each cluster are summarized in Table 2. Table 2 shows that 6 out of 18 (33%) PC traits have at least one regulator, while 16 out of 743 (2.2%) individual traits have at least one regulator. We conclude that considering one cluster at a time may show more linkage signals than considering traits individually. Figure 2 gives

*p*-value curves of PC traits and individual traits in cluster 5. The figure indicates that the linkage peaks of all individual traits and the PC trait of the cluster appear to be very close to each other. It indicates that the traits in one cluster have similar regulators, and we can find these regulators by using the PC trait. Furthermore, we can see from Table 2 that in 16 out of 18 clusters, if the PC trait has regulators, there will be some individual traits with regulators, and vice versa. This indicates that the PC trait is a reasonable measure to represent all the traits in the corresponding cluster.

Details of the 18 clusters.

Cluster ID | No. regulators for the PC trait | No. traits | No. RE traits |
---|---|---|---|

1 | 1 | 21 | 0 |

2 | 0 | 43 | 0 |

3 | 2 | 32 | 6 |

4 | 0 | 27 | 0 |

5 | 1 | 28 | 0 |

6 | 0 | 35 | 0 |

7 | 0 | 35 | 0 |

8 | 2 | 47 | 5 |

9 | 0 | 33 | 0 |

10 | 2 | 27 | 1 |

11 | 0 | 64 | 0 |

12 | 0 | 42 | 0 |

13 | 0 | 84 | 2 |

14 | 1 | 59 | 2 |

15 | 0 | 44 | 0 |

16 | 0 | 41 | 0 |

17 | 0 | 47 | 0 |

18 | 0 | 44 | 0 |

## Discussion

In order to analyze a large number of correlated traits given in the GAW15 Problem 1, we have developed a method that first clusters the traits into tight clusters and then finds a representative measure for each cluster. We then performed a linkage scan for the representative trait. The results of analyzing the GAW15 Problem 1 data set indicate that the traits in the same cluster have the same regulators, and that the first PC in one cluster is a good representative measure of this cluster. The results also show that performing a genome scan for one cluster at a time may show more linkage signals than that for genome scans of individual traits.

One remaining problem is how to choose the values of parameters *α*_{0}, *m*_{1}, and *m*_{2} in the tight clustering algorithm. We have applied the algorithm to the GAW15 data with different parameter values. Our results indicate that when *m*_{2} varies from 100 to 1000, there is almost no effect on the cluster assignment and linkage analyses. When *m*_{1} is small, say 2, there will be more small clusters and more linkage signals. Because two biologically unrelated genes may have highly correlated traits by chance, by using a small *m*_{1}, the tight clustering algorithm may cluster biologically unrelated genes into one cluster. We suggest choosing *m*_{1} between 10 and 20, and *m*_{2} between 100 and 400. Our independent simulation studies indicate that if *α*_{0} is too small, the FDR cannot be controlled, and the first PC trait will not be a good representative of the cluster, because the first PC can only explain a small part of the total variation. Thus we suggest using *α*_{0} ≥ 0.5.

## Declarations

### Acknowledgements

This work was supported by National Institute of Health (NIH) grants R01 GM069940, R03 HG003613, R03 AG024491, and R01 HG003054.

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.

## Authors’ Affiliations

## References

- 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
- Cheung VG, Spielman RS, Ewens KG, Weber TM, Morley M, Burdick JT: Mapping determinants of human gene expression by regional and whole genome association. Nature. 2005, 437: 1365-1369. 10.1038/nature04244.View ArticlePubMed CentralPubMedGoogle Scholar
- Morley M, Molony CM, Weber T, 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
- Tseng GC, Wong WH: Tight clustering: a resampling-based approach for identifying stable and tight patterns in data. Biometrics. 2005, 61: 10-16. 10.1111/j.0006-341X.2005.031032.x.View ArticlePubMedGoogle Scholar
- Wigginton JE, Abecasis GR: PEDSTATS: descriptive statistics, graphics and quality assessment for gene mapping data. Bioinformatics. 2005, 21: 3445-3447. 10.1093/bioinformatics/bti529.View ArticlePubMedGoogle 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
- Benjamini Y, Yekutieli D: The control of the false discovery rate in multiple testing under dependency. Ann Stat. 2001, 29: 1165-1188. 10.1214/aos/1013699998.View ArticleGoogle Scholar

## Copyright

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.