Using the ratio of means as the effect size measure in combining results of microarray experiments

<p>Abstract</p> <p>Background</p> <p>Development of efficient analytic methodologies for combining microarray results is a major challenge in gene expression analysis. The widely used effect size models are thought to provide an efficient modeling framework for this purpose, where the measures of association for each study and each gene are combined, weighted by the standard errors. A significant disadvantage of this strategy is that the quality of different data sets may be highly variable, but this information is usually neglected during the integration. Moreover, it is widely known that the estimated standard deviations are probably unstable in the commonly used effect size measures (such as standardized mean difference) when sample sizes in each group are small.</p> <p>Results</p> <p>We propose a re-parameterization of the traditional mean difference based effect measure by using the log ratio of means as an effect size measure for each gene in each study. The estimated effect sizes for all studies were then combined under two modeling frameworks: the quality-unweighted random effects models and the quality-weighted random effects models. We defined the quality measure as a function of the detection p-value, which indicates whether a transcript is reliably detected or not on the Affymetrix gene chip. The new effect size measure is evaluated and compared under the quality-weighted and quality-unweighted data integration frameworks using simulated data sets, and also in several data sets of prostate cancer patients and controls. We focus on identifying differentially expressed biomarkers for prediction of cancer outcomes.</p> <p>Conclusion</p> <p>Our results show that the proposed effect size measure (log ratio of means) has better power to identify differentially expressed genes, and that the detected genes have better performance in predicting cancer outcomes than the commonly used effect size measure, the standardized mean difference (SMD), under both quality-weighted and quality-unweighted data integration frameworks. The new effect size measure and the quality-weighted microarray data integration framework provide efficient ways to combine microarray results.</p

Linkage and association analysis in pedigrees from different populations

Using the Genetic Analysis Workshop 14 simulated datasets we carried out nonparametric linkage analyses and applied a log-linear method for analysis of case-parent-triad data with stratification on parental mating type. We proposed and applied a random effect modelling approach to explore the impact of population heterogeneity on tests of association between genetic markers and disease status. The estimated genetic effect may appear to be strongly significant in one population but nonsignificant in another population, leading to confusion about interpretation. However, when results are interpreted in the light of a random effects model, both studies may be making similar statements about a genetic effect that varies depending on environment and background

Pathway analysis for genetic association studies: to do, or not to do? That is the question

In Genetic Analysis Workshop 18 data, we used a 3-stage approach to explore the benefits of pathway analysis in improving a model to predict 2 diastolic blood pressure phenotypes as a function of genetic variation. At stage 1, gene-based tests of association in family data of approximately 800 individuals found over 600 genes associated at p<0.05 for each phenotype. At stage 2, networks and enriched pathways were estimated with Cytoscape for genes from stage 1, separately for the 2 phenotypes, then examining network overlap. This overlap identified 4 enriched pathways, and 3 of these pathways appear to interact, and are likely candidates for playing a role in hypertension. At stage 3, using 157 maximally unrelated individuals, partial least squares regression was used to find associations between diastolic blood pressure and single-nucleotide polymorphisms in genes highlighted by the pathway analyses. However, we saw no improvement in the adjusted cross-validated R(2). Although our pathway-motivated regressions did not improve prediction of diastolic blood pressure, merging gene networks did identify several plausible pathways for hypertension

Identifying cis- and trans-acting single-nucleotide polymorphisms controlling lymphocyte gene expression in humans

Assuming multiple loci play a role in regulating the expression level of a single phenotype, we propose a new approach to identify cis- and trans-acting loci that regulate gene expression. Using the Problem 1 data set made available for Genetic Analysis Workshop 15 (GAW15), we identified many expression phenotypes that have significant evidence of association and linkage to one or more chromosomal regions. In particular, six of ten phenotypes that we found to be regulated by cis- and trans-acting loci were also mapped by a previous analysis of these data in which a total of 27 phenotypes were identified with expression levels regulated by cis-acting determinants. However, in general, the p-values associated with these regulators identified in our study were larger than in their studies, since we had also identified other factors regulating expression. In fact, we found that most of the gene expression phenotypes are influenced by at least one trans-acting locus. Our study also shows that much of the observable heritability in the phenotypes could be explained by simple single-nucleotide polymorphism associations; residual heritability was reduced and the remaining heritability may represent complex regulation systems with interactions or noise

A genome scan for parent-of-origin linkage effects in alcoholism

BACKGROUND: Alcoholism is a complex disease in which genomic imprinting may play an important role in its susceptibility. OBJECTIVE: To conduct a genome-wide search for loci that may have strong parent-of-origin linkage effects in alcoholism; to compare the linkage results between the microsatellites and the two single-nucleotide polymorphism (SNP) platforms. METHODS: Nonparametric linkage analyses were performed using ALLEGRO with the three sets of markers provided by the Genetic Analysis Workshop 14 for the Collaborative Study on the Genetics of Alcoholism Problem 1 data. Both sex-averaged and sex-specific genetic maps were used. We also provided a valid statistical test to determine whether the parental allele sharing differed significantly. RESULTS: Significant maternal linkage effects (paternal imprinting) were observed on chromosome 12 using either the microsatellite markers or the two SNP panels. The two SNP sets did not improve the linkage signals compared to the results from the microsatellite markers on chromosome 12. Possible paternal linkage effects (maternal imprinting) on chromosome 7 and maternal linkage effects (paternal imprinting) on chromosome 10 were found using the two SNP panels. CONCLUSION: For diseases which may have parent-of-origin effects, linkage analysis looking at parental sharing separately may reduce locus heterogeneity and increase the ability to identify that which can not be identified with usual linkage analysis

A genome scan for parent-of-origin linkage effects in alcoholism

BACKGROUND: Alcoholism is a complex disease in which genomic imprinting may play an important role in its susceptibility. OBJECTIVE: To conduct a genome-wide search for loci that may have strong parent-of-origin linkage effects in alcoholism; to compare the linkage results between the microsatellites and the two single-nucleotide polymorphism (SNP) platforms. METHODS: Nonparametric linkage analyses were performed using ALLEGRO with the three sets of markers provided by the Genetic Analysis Workshop 14 for the Collaborative Study on the Genetics of Alcoholism Problem 1 data. Both sex-averaged and sex-specific genetic maps were used. We also provided a valid statistical test to determine whether the parental allele sharing differed significantly. RESULTS: Significant maternal linkage effects (paternal imprinting) were observed on chromosome 12 using either the microsatellite markers or the two SNP panels. The two SNP sets did not improve the linkage signals compared to the results from the microsatellite markers on chromosome 12. Possible paternal linkage effects (maternal imprinting) on chromosome 7 and maternal linkage effects (paternal imprinting) on chromosome 10 were found using the two SNP panels. CONCLUSION: For diseases which may have parent-of-origin effects, linkage analysis looking at parental sharing separately may reduce locus heterogeneity and increase the ability to identify that which can not be identified with usual linkage analysis

Recursive partitioning models for linkage in COGA data

We have developed a recursive-partitioning (RP) algorithm for identifying phenotype and covariate groupings that interact with the evidence for linkage. This data-mining approach for detecting gene × environment interactions uses genotype and covariate data on affected relative pairs to find evidence for linkage heterogeneity across covariate-defined subgroups. We adapted a likelihood-ratio based test of linkage parameterized with relative risks to a recursive partitioning framework, including a cross-validation based deviance measurement for choosing optimal tree size and a bootstrap sampling procedure for choosing robust tree structure. ALDX2 category 5 individuals were considered affected, categories 1 and 3 unaffected, and all others unknown. We sampled non-overlapping affected relative pairs from each family; therefore, we used 144 affected pairs in the RP model. Twenty pair-level covariates were defined from smoking status, maximum drinks, ethnicity, sex, and age at onset. Using the all-pairs score in GENEHUNTER, the nonparametric linkage tests showed no regions with suggestive linkage evidence. However, using the RP model, several suggestive regions were found on chromosomes 2, 4, 6, 14, and 20, with detection of associated covariates such as sex and age at onset

Chromosome-breakage genomic instability and chromothripsis in breast cancer

Background: Chromosomal breakage followed by faulty DNA repair leads to gene amplifications and deletions in cancers. However, the mere assessment of the extent of genomic changes, amplifications and deletions may reduce the complexity of genomic data observed by array comparative genomic hybridization (array CGH). We present here a novel approach to array CGH data analysis, which focuses on putative breakpoints responsible for rearrangements within the genome. Results: We performed array comparative genomic hybridization in 29 primary tumors from high risk patients with breast cancer. The specimens were flow sorted according to ploidy to increase tumor cell purity prior to array CGH. We describe the number of chromosomal breaks as well as the patterns of breaks on individual chromosomes in each tumor. There were differences in chromosomal breakage patterns between the 3 clinical subtypes of breast cancers, although the highest density of breaks occurred at chromosome 17 in all subtypes, suggesting a particular proclivity of this chromosome for breaks. We also observed chromothripsis affecting various chromosomes in 41% of high risk breast cancers. Conclusions: Our results provide a new insight into the genomic complexity of breast cancer. Genomic instability dependent on chromosomal breakage events is not stochastic, targeting some chromosomes clearly more than others. We report a much higher percentage of chromothripsis than described previously in other cancers and this suggests that massive genomic rearrangements occurring in a single catastrophic event may shape many breast cancer genomes

How old is this mutation? - a study of three Ashkenazi Jewish founder mutations

Abstract Background Several founder mutations leading to increased risk of cancer among Ashkenazi Jewish individuals have been identified, and some estimates of the age of the mutations have been published. A variety of different methods have been used previously to estimate the age of the mutations. Here three datasets containing genotype information near known founder mutations are reanalyzed in order to compare three approaches for estimating the age of a mutation. The methods are: (a) the single marker method used by Risch et al., (1995); (b) the intra-allelic coalescent model known as DMLE, and (c) the Goldgar method proposed in Neuhausen et al. (1996), and modified slightly by our group. The three mutations analyzed were MSH2*1906 G->C, APC*I1307K, and BRCA2*6174delT. Results All methods depend on accurate estimates of inter-marker recombination rates. The modified Goldgar method allows for marker mutation as well as recombination, but requires prior estimates of the possible haplotypes carrying the mutation for each individual. It does not incorporate population growth rates. The DMLE method simultaneously estimates the haplotypes with the mutation age, and builds in the population growth rate. The single marker estimates, however, are more sensitive to the recombination rates and are unstable. Mutation age estimates based on DMLE are 16.8 generations for MSH2 (95% credible interval (13, 23)), 106 generations for I1037K (86-129), and 90 generations for 6174delT (71-114). Conclusions For recent founder mutations where marker mutations are unlikely to have occurred, both DMLE and the Goldgar method can give good results. Caution is necessary for older mutations, especially if the effective population size may have remained small for a long period of time