- Software
- Open Access

# catmap: Case-control And TDT Meta-Analysis Package

- Kristin K Nicodemus
^{1, 2, 3}Email author

**9**:130

https://doi.org/10.1186/1471-2105-9-130

© Nicodemus; licensee BioMed Central Ltd. 2008

**Received:**30 July 2007**Accepted:**28 February 2008**Published:**28 February 2008

## Abstract

### Background

Risk for complex disease is thought to be controlled by multiple genetic risk factors, each with small individual effects. Meta-analyses of several independent studies may be helpful to increase the ability to detect association when effect sizes are modest. Although many software options are available for meta-analysis of genetic case-control data, no currently available software implements the method described by Kazeem and Farrall (2005), which combines data from independent family-based and case-control studies.

### Results

I introduce the package catmap for the R statistical computing environment that implements fixed- and random-effects pooled estimates for case-control and transmission disequilibrium methods, allowing for the use of genetic association data across study types. In addition, catmap may be used to create forest and funnel plots and to perform sensitivity analysis and cumulative meta-analysis. catmap is available from the Comprehensive R Archive Network http://www.r-project.org.

### Conclusion

catmap allows researchers to synthesize data to assess evidence for association in studies of genetic polymorphisms, facilitating the use of pooled data analyses which may increase power to detect moderate genetic associations.

## Keywords

- Transmission Disequilibrium Test
- Portable Document Format
- Multiple Independent Study
- Single Genetic Effect
- Multiple Genetic Risk

## Background

Two study designs are commonly employed in genetic association studies: a case-control and a family-based approach. The case-control design compares frequencies of alleles carried among cases with a disease and among controls that are free of disease. The family-based design compares the frequency of alleles transmitted to an affected offspring by their parents with alleles carried by the parents but not passed to the offspring; this type of statistical analysis is often called a transmission disequilibrium test (TDT). Complex diseases are likely to be under the influence of several genetic risk factors; therefore, the contribution of a single gene to risk of disease is expected to be modest. Meta-analysis allows for the pooling of independent studies that examine similar hypotheses; for example, that a particular allele at a SNP is associated with disease status, and thus may improve power to detect moderate effect sizes. Although a few methods have been developed to combine family-based data with data from unrelated controls and/or unrelated cases obtained by the same study [1–10], most of these methods are not specifically designed to pool results from multiple independent studies and thus do not have a built-in test for heterogeneity of effect, which may be used to determine whether pooling of independent samples is reasonable. Exceptions to this statement include the method implemented in Genie[7, 8], which conducts a genotype-based test that allows for the pooling of multiple study types across independent samples (meta-analysis), relying on empirical p-values which can be computationally intensive and SCOUT[5], which allows for the pooling of family-based data with data from unrelated cases and/or controls and performs a formal test of whether differing data types can be sensibly pooled. The likelihood-based method implemented in SCOUT[5] assumes that pooled samples are from the same source population with the same mating-type frequencies, and thus is inappropriate to use in a meta-analysis setting where independent samples have been drawn from differing populations; however, the use of additional programs such as the R package meta [10] subsequent to using SCOUT may be performed to pool estimates. A recent paper outlined a method for the combination of odds ratios (ORs) from independent case-control and TDT-based studies using a fixed-effect approach [6], using an allele-based model. Although several programs exist to conduct meta-analyses of case-control genetic data, no software exists that implements the Kazeem & Farrall (2005)[6] method to conduct case-control and TDT meta-analyses of independent samples. I have implemented this model [6], plus extended the method to the random-effects model of DerSimonian and Laird [11] in the freely-available R [12] package catmap.

## Implementation

*w*

_{ i }= 1/

*σ*

_{i}

^{2}), which gives increased weight to larger studies (because their variances are generally smaller). The pooled estimate of the logarithm of the odds ratios using fixed effects methods (OR

_{FE}) is given by[6]:

*w*

_{ i }are the study-specific weights as described above. In the presence of heterogeneity of genetic effects random-effects estimates of the pooled OR and variance may be more appropriate. Random-effects models do not assume that a common genetic effect exists across all samples. Following DerSimonian and Laird [11, 13], the random-effects model implemented in catmap includes an estimate of the between study variance,

*τ*

^{2}, as:

*χ*

^{2}statistic (see [6] for review), K is the number of studies and the

*w*

_{i}are the study weights calculated using fixed-effects estimates as described above. If

*τ*

^{2}≤ 0,

*τ*

^{2}is set to 0 and random-effects estimates are identical to fixed-effects estimates. If

*τ*

^{2}> 0 the random-effects weight for study

*i*is

*r*

_{ i }= 1/(w

_{i}

^{-1}+

*τ*

^{2}) and the random-effects pooled estimate of the logarithm of the odds ratio (OR

_{RE}) can be written as:

and the standard error of the random-effects pooled ln(OR_{RE}) is simply 1/(Σ*r*_{
i
})^{1/2}.

## Results and Discussion

Table 1

name | study | t | nt | caserisk | controlrisk | casenotrisk | controlnotrisk |
---|---|---|---|---|---|---|---|

Yu,2004 | 1 | 78 | 21 | NA | NA | NA | NA |

Wu,2006 | 2 | 0 | 0 | 45 | 56 | 120 | 354 |

**t**and

**nt**columns should contain counts and the additional 4 columns may contain 0 or NA; for case-control studies, the columns

**t**and

**nt**should contain 0 or NA and the remaining 4 columns should contain the allele count data for cases and controls. catmap allows for the meta-analysis of all case-control or all family-based studies and for the combination of both family-based and case-control studies. In situations where both family-based and case-control data exist for the same set of cases only one type of study design (either family-based or case-control) should be used as input to catmap. catmap calculates fixed- and random-effects estimates, a

*χ*

^{2}test for heterogeneity and associated p-values and must be used to create a catmap object for use in downstream functions. To create a forest plot using either the pooled fixed- or random-effects OR and 95% CI the function catmap.forest is used. catmap.sense calculates leave-one-out sensitivity analyses using either fixed- or random-effects estimates and creates a Portable Document Format (PDF) plot of the estimate of the pooled OR and 95% CIs with the removed study listed on the y-axis (Figure 1). Because the first report of a genetic association between a polymorphism and disease status often is larger than those reported by follow-up studies [14, 15], the function catmap.cumulative may be used to calculate cumulative ORs and associated CIs and to create a PDF plot of the pooled estimates, adding one study at a time, using either fixed- or random-effects analyses. A graphic in PDF format to assess publication bias can be created using catmap.funnel.

## Conclusion

Pooling studies of genetic polymorphisms data via meta-analysis may improve the ability of researchers to detect modest association signals. Since case-control studies are not robust to spurious effects caused by confounding via population stratification, before combining results from case-control studies researchers should assess evidence for substructure in individual studies and also examine the similarity of effect sizes and allele frequencies across study type and population, along with the formal test of heterogeneity of genetic effect implemented in catmap. catmap implements a random-effects estimate of the pooled odds ratio, which assumes a Normal distribution of genetic effects across studies instead of a single genetic effect, which may be preferred when combining studies across sub-groups or when the heterogeneity *χ*^{2} indicates significant between-study differences in effect sizes, although when evidence for heterogeneity exists the pooling of data may not be indicated. Since two types of study design (case-control and family-based) are common in statistical genetics, methods and freely-available software that can combine estimates from both designs should prove useful to applied researchers.

## Availability and Requirements

**Project name**: catmap

**Project home page**: Contributed package at http://www.r-project.org.

**Operating systems**: MS Windows, Linux, Mac

**Programming language**: R

**Other requirements**: R 2.4 or higher

**License**: GPL

**Any restrictions to use by non-academics**: none

## Declarations

## Authors’ Affiliations

## References

- Schaid DJ, Rowland C: Use of parents, sibs, and unrelated controls for detection of associations between genetic markers and disease.
*Am J Hum Genet*1998, 63(5):1492–506. 10.1086/302094PubMed CentralView ArticlePubMedGoogle Scholar - Becker T, Knapp M: Maximum-likelihood estimation of haplotype frequencies in nuclear families.
*Genet Epidemiol*2004, 27: 21–32. 10.1002/gepi.10323View ArticlePubMedGoogle Scholar - Nagelkerke NJ, Hoebee B, Teunis P, Kimman TG: Combining the transmission disequilibrium test and case-control methodology using generalized logistic regression.
*Eur J Hum Genet*2004, 12: 964–970. 10.1038/sj.ejhg.5201255View ArticlePubMedGoogle Scholar - Becker T, Knapp M: Impact of missing genotype data on Monte-Carlo simulation based haplotype analysis.
*Hum Hered*2005, 59: 185–189. 10.1159/000086696View ArticlePubMedGoogle Scholar - Epstein MP, Veal CD, Trembath RC, Barker JNWN, Li C, Satten GA: Genetic association analysis using data from triads and unrelated subjects.
*Am J Hum Genet*2005, 76: 592–608. 10.1086/429225PubMed CentralView ArticlePubMedGoogle Scholar - Kazeem GR, Farrall M: Integrating case-control and TDT studies.
*Ann Hum Genet*2005, 69: 329–335. 10.1046/J.1469-1809.2005.00156.xView ArticlePubMedGoogle Scholar - Allen-Brady K, Wong J, Camp NJ: PedGenie: an analysis approach for genetic association testing in extended pedigrees and genealogies of arbitrary size.
*BMC Bioinformatics*2006, 7: 209. 10.1186/1471-2105-7-209PubMed CentralView ArticlePubMedGoogle Scholar - Curtin K, Wong J, Allen-Brady K, Camp NJ: PedGenie: meta genetic association testing in mixed family and case-control designs.
*BMC Bioinformatics*2007, 8: 448. 10.1186/1471-2105-8-448PubMed CentralView ArticlePubMedGoogle Scholar - Dudbridge F: UNPHASED user guide. In Technical Report 2006/5. MRC Biostatistics Unit, Cambridge UK;Google Scholar
- Schwarzer G: The meta Package. [R package version 0.8–2]Google Scholar
- DerSimonian R, Laird N: Meta-analysis in clinical trials.
*Control Clin Trials*1986, 7: 177–188. 10.1016/0197-2456(86)90046-2View ArticlePubMedGoogle Scholar - R Development Core Team:
*R: A Language and Environment for Statistical Computing.*R Foundation for Statistical Computing, Vienna, Austria; 2007.Google Scholar - Egger M, Davey Smith G, Altman DG:
*Systematic reviews in health care: Meta-analysis in context*. London: BMJ Publishing Group; 2001.View ArticleGoogle Scholar - Ioannidis JP, Ntzani EE, Trikalinos TA, Contopoulos-Ioannidis DG: Replication validity of genetic association studies.
*Nat Genet*2001, 29: 306–309. 10.1038/ng749View ArticlePubMedGoogle Scholar - Trikalinos TA, Ntzani EE, Contopoulos-Ioannidis DG, Ioannidis JP: Establishment of genetic associations for complex diseases is independent of early study findings.
*Eur J Hum Genet*2004, 12: 762–769. 10.1038/sj.ejhg.5201227View ArticlePubMedGoogle 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.