PubMed Health⌕ Search

Biomedical subjects

L P Zhao

Publications and source records attributed to L P Zhao.

At least 19 recordsLinked to original sources

A statistically driven approach for image segmentation and signal extraction in cDNA microarrays.

The increasing use of cDNA microarrays necessitates the development of methods for extracting quality data. Here, we set forth hurdles to overcome in image analysis of microarrays. We emphasize the importance of objective data extraction methods resulting in reliable signal estimates. Based on statistical principles, we describe a method for automated grid alignment, spot detection, background estimation, flagging, and signal extraction. A software application that we call SignalViewer has been implemented for this method. We identify areas where we improved upon current methods used for array image analysis at each step in the process. Finally, we give examples to illustrate the performance of our algorithms on raw data.

Algorithms↗

SignalViewer: analyzing microarray images.

Microarray technology is now routinely used to monitor genome-wide expression profiles. However, current microarray imaging and analysis packages typically require manual intervention and assumptions on alignments. Unfortunately, limitations and assumptions are typically undisclosed and methods are not published. To facilitate exploration of image data, we developed SignalViewer. This paper presents a description of the application.

Artifacts↗

Statistical modeling of large microarray data sets to identify stimulus-response profiles.

A statistical modeling approach is proposed for use in searching large microarray data sets for genes that have a transcriptional response to a stimulus. The approach is unrestricted with respect to the timing, magnitude or duration of the response, or the overall abundance of the transcript. The statistical model makes an accommodation for systematic heterogeneity in expression levels. Corresponding data analyses provide gene-specific information, and the approach provides a means for evaluating the statistical significance of such information. To illustrate this strategy we have derived a model to depict the profile expected for a periodically transcribed gene and used it to look for budding yeast transcripts that adhere to this profile. Using objective criteria, this method identifies 81% of the known periodic transcripts and 1,088 genes, which show significant periodicity in at least one of the three data sets analyzed. However, only one-quarter of these genes show significant oscillations in at least two data sets and can be classified as periodic with high confidence. The method provides estimates of the mean activation and deactivation times, induced and basal expression levels, and statistical measures of the precision of these estimates for each periodic transcript.

CDC28 Protein Kinase, S cerevisiae↗

Incorporating age at onset information into the transmission/disequilibrium test.

In this paper, we proposed two approaches for incorporating the age at onset into the transmission/disequilibrium test (TDT). Trios (affected offspring and their parents) were extracted from the first four replicate data sets of the general population type. Focusing on chromosome 6 where MG6 and MG7 reside, we compared the usual TDT with the newly proposed tests in terms of gene localization.

Adult↗

An efficient and robust statistical modeling approach to discover differentially expressed genes using genomic expression profiles.

We have developed a statistical regression modeling approach to discover genes that are differentially expressed between two predefined sample groups in DNA microarray experiments. Our model is based on well-defined assumptions, uses rigorous and well-characterized statistical measures, and accounts for the heterogeneity and genomic complexity of the data. In contrast to cluster analysis, which attempts to define groups of genes and/or samples that share common overall expression profiles, our modeling approach uses known sample group membership to focus on expression profiles of individual genes in a sensitive and robust manner. Further, this approach can be used to test statistical hypotheses about gene expression. To demonstrate this methodology, we compared the expression profiles of 11 acute myeloid leukemia (AML) and 27 acute lymphoblastic leukemia (ALL) samples from a previous study (Golub et al. 1999) and found 141 genes differentially expressed between AML and ALL with a 1% significance at the genomic level. Using this modeling approach to compare different sample groups within the AML samples, we identified a group of genes whose expression profiles correlated with that of thrombopoietin and found that genes whose expression associated with AML treatment outcome lie in recurrent chromosomal locations. Our results are compared with those obtained using t-tests or Wilcoxon rank sum statistics.

Acute Disease↗

Assessing changes in ages at onset over successive generation: an application to breast cancer.

Decreased age at onset in successive generations has been observed for a number of diseases. Two nonparametric matched and unmatched test statistics are proposed, taking into account not only current age or age at death for unaffected individuals and age at disease onset for affected individuals, but also possible correlations among family members. Both are asymptotically normal with readily estimated variances from the data. A simulation study is conducted to compare the proposed tests with the commonly used paired t-test and log-rank test. It has been shown that the proposed test statistics yield valid conclusions in assessing genetic anticipation under all situations considered. However, the paired t-test is valid only when the censoring distributions are comparable between two generations, whereas the log-rank test is valid when the correlation among family members is weak. As expected, the matched test is most powerful when the data are heterogeneous, and the unmatched and the log-rank tests are most powerful when the data are homogeneous and the correlation is weak. Lastly, a population-based family study of breast cancer conducted at the Fred Hutchinson Cancer Research Center is used for illustration of the proposed and the log-rank tests. The preliminary analysis suggests that there appears a decreased age at onset over the successive generations in breast cancer.

Adolescent↗

An efficient, robust and unified method for mapping complex traits (III): combined linkage/linkage-disequilibrium analysis.

Extending the method for linkage analysis [Zhao et al., 1998a: Am. J. Med. Genet. 77:366-383; 1998b: Am. J. Med. Genet. 79:49-61], this article describes a method for the linkage-disequilibrium analysis, and for combining linkage and linkage-disequilibrium analyses. As highly dense markers are increasingly used in genome scans, one or more markers are not only linked with the disease genes if they exist, but also likely in linkage-disequilibrium with those putative genes. Hence, linkage-disequilibrium analysis potentially offers additional information about positions of putative disease genes. Combining both linkage and linkage-disequilibrium signals, this approach is able to improve positional signals. As before, the proposed method is a model-based approach, but semiparametric via the estimating equation technique. Under the assumptions of penetrance and allele frequency, this method efficiently estimates recombination fractions for linkage analysis and odds ratios for linkage-disequilibrium analysis. As described in two previous papers, this method is relatively more robust than the lod score methods, since it requires weaker assumption than conditional independence. While the estimated recombination fractions are used for inference as part of linkage analysis, the estimated odds ratios are used for linkage-disequilibrium inference and combined linkage, and linkage-disequilibrium parameters can be used to test combined linkage/linkage-disequilibrium analysis. This approach has been implemented, named gSCAN, and its compiled version is available for trial on request via the web site (http:/lynx.fhcrc.org/qge). We applied this new approach to affected sib-pair data collected for the genome scan to localize type 1 diabetes genes. Under an assumed autosomal dominant gene model, the linkage analysis confirms an earlier suggestion of one major gene around D6S281. Interestingly, the linkage-disequilibrium analysis suggests several additional signals around D6S250, GATA30, D6S311, D6S441, D6S442, D6S415, D6S411, D6S305, and a290xh9. The linkage analysis, on the other hand, suggests a signal around D6S281, while providing supporting evidence for several other marker loci. However, the combined analysis did not provide strong support for any of the findings, implying that linkage and linkage-disequilibrium findings are not consistent.

Alleles↗

Mapping alcoholism genes using linkage/linkage disequilibrium analysis.

Using a recently developed semiparametric method for combined linkage/linkage-disequilibrium analysis, we analyzed the Collaborative Study on the Genetics of Alcoholism data subset developed for Genetic Analysis Workshop 11 (GAW11). This semiparametric approach estimates recombination fractions for linkage, marker log odds ratios for linkage-disequilibrium, their product for combined linkage/linkage-disequilibrium, and corresponding z-scores. We used two outcomes: alcohol dependence and "alcoholism-free" and a genome-wide significance level of 4.1 (which corresponds to a genome-wide lod score of 3.6). For the alcohol dependence outcome, we observed significant linkage signals at D1S1588-D1S1631, D1S547, D2S399, D2S425, D4S2361, D7S1796, and D7S1824. We also found significant linkage-disequilibrium signals at D1S547 and D7S1795. For the "alcoholism-free" outcome, we found significant linkage signals at D4S2457, D41651 (both flank ADH3), D11S2359, and D16S47 and significant linkage-disequilibrium signals at D4S2361, FABP2, D11S2359, D19S431 and D19S47-D19S198-D19S601.

Alcoholism↗

A genome-wide scan for a simulated data set using two newly developed methods.

A genome-wide scan of a simulated data set for fictitious disease genes was conducted using both semiparametric and nonparametric methods. The semiparametric model-based method, which tests for linkage/linkage disequilibrium separately and together, correctly identified all three underlying disease loci along with two false positives through the linkage analysis. However, the nonparametric model-free method which tests combined linkage/linkage disequilibrium, failed to yield any results due to the lack of linkage disequilibrium information in the data.

Chromosome Mapping↗

On the assessment of statistical significance in disease-gene discovery.

One of the major challenges facing genome-scan studies to discover disease genes is the assessment of the genomewide significance. The assessment becomes particularly challenging if the scan involves a large number of markers collected from a relatively small number of meioses. Typically, this assessment has two objectives: to assess genomewide significance under the null hypothesis of no linkage and to evaluate true-positive and false-positive prediction error rates under alternative hypotheses. The distinction between these goals allows one to formulate the problem in the well-established paradigm of statistical hypothesis testing. Within this paradigm, we evaluate the traditional criterion of LOD score 3.0 and a recent suggestion of LOD score 3.6, using the Monte Carlo simulation method. The Monte Carlo experiments show that the type I error varies with the chromosome length, with the number of markers, and also with sample sizes. For a typical setup with 50 informative meioses on 50 markers uniformly distributed on a chromosome of average length (i.e., 150 cM), the use of LOD score 3.0 entails an estimated chromosomewide type I error rate of.00574, leading to a genomewide significance level >.05. In contrast, the corresponding type I error for LOD score 3.6 is.00191, giving a genomewide significance level of slightly <.05. However, with a larger sample size and a shorter chromosome, a LOD score between 3.0 and 3.6 may be preferred, on the basis of proximity to the targeted type I error. In terms of reliability, these two LOD-score criteria appear not to have appreciable differences. These simulation experiments also identified factors that influence power and reliability, shedding light on the design of genome-scan studies.

Chromosome Mapping↗

Integrated designs for gene discovery and characterization.

Recent advances, including near completion of the human genome map, ever improving high-throughput technologies, and successes in discovering chronic disease-related genes, have stimulated the further development of genetic epidemiology. The primary mission of genetic epidemiology is to discover and characterize genes, whether independent of or interactive with environmental factors, that cause human diseases. To accomplish such a mission, genetic epidemiology needs to integrate both genetic and epidemiologic approaches. One of the challenges facing such an integrated approach is the identification of study designs that are efficient for both gene discovery and characterization. Because designs for gene discovery alone and designs for gene characterization alone have been elaborated in the other two panels, the focus of this paper is to describe those designs that may be useful for discovery and characterization jointly, including case-family and case-control-family designs. Examples of integrated designs are described, and studies of breast cancer conducted at the Fred Hutchinson Cancer Research Center are used for illustration. Finally, related analytic issues are also discussed.

Breast Neoplasms↗

An efficient, robust, and unified method for mapping complex traits (II): multipoint linkage analysis.

Extending the method for two-point linkage analysis [Zhao et al., 1998: Am J Med Genet 77:366-383], this paper introduces a semiparametric method for multipoint linkage analysis, expected to gain efficiency by using multiple markers simultaneously. Overcoming the longstanding statistical and computational challenge to the parametric approaches (or lod score methods) for multipoint linkage analysis, this semiparametric approach, based on the estimating equation technique, yields statistically efficient and yet robust estimates and enjoys the computational efficiency in processing multiple markers from large pedigrees. Its computational burden increases linearly with the sizes of pedigrees and with the number of marker loci. To illustrate this semiparametric method, we apply it to marker data gathered for the Breast Cancer Consortium. The result supports the earlier finding of the positive linkage with BRCA1 and has also shown that the multipoint linkage analysis has an improved power. In addition, we have applied this method to analyze genome scanning data that have been used to localize genes responsible for type 1 diabetes. In support of the earlier findings, the genome scanning detects the linkage signals on chromosome 6 but does not support the earlier suggestions of two major genes in that genome segment. Through sensitivity analysis, it appears that the results are robust to misspecification of penetrance and allele frequency.

Automation↗

Embryonic cardiomyocyte hypoplasia and craniofacial defects in G alpha q/G alpha 11-mutant mice.

Heterotrimeric G proteins of the Gq class have been implicated in signaling pathways regulating cardiac growth under physiological and pathological conditions. Knockout mice carrying inactivating mutations in both of the widely expressed G alpha q class genes, G alpha q and G alpha 11, demonstrate that at least two active alleles of these genes are required for extrauterine life. Mice carrying only one intact allele [G alpha q(-/+);G alpha 11(-/-) or G alpha q(-/-);G alpha 11(-/+)] died shortly after birth. These mutants showed a high incidence of cardiac malformation. In addition, G alpha q(-/-);G alpha 11(-/+) newborns suffered from craniofacial defects. Mice lacking both G alpha q and G alpha 11 [G alpha q(-/-);G alpha 11(-/-)] died at embryonic day 11 due to cardiomyocyte hypoplasia. These data demonstrate overlap in G alpha q and G alpha 11 gene functions and indicate that the Gq class of G proteins plays a crucial role in cardiac growth and development.

Alleles↗

Efficient, robust, and unified method for mapping complex traits (I): two-point linkage analysis.

The completion of a preliminary human genome map and development of molecular methods have enabled researchers to assay a large number of polymorphic markers that are evenly spaced along the entire human genome. Among many applications, marker data are valuable for mapping complex traits through linkage or linkage-disequilibrium analysis, the former of which is the focus of this paper, the first in a series on this subject. Formalizing the concept and computation for linkage analysis, Elston and Stewart [1971; Human Heredity 21:523-542] introduced a likelihood function to capture relevant genetic information and a recursive algorithm for computing the likelihood function. However, the computing burden is prohibitive in processing complex pedigrees. Since that fundamental development, improving the computational algorithm and extending the method has been a dynamic area of research. The primary objective of this communication is to introduce a semiparametric method for linkage analysis. It is a particularly suitable approach with desirable properties for mapping complex traits that may be binary, continuous, and partially observed (i.e., censored). It incorporates candidate genes, environmental factors, and their interactions with the putative gene and is expected to be robust and efficient in comparison with likelihood-based methods. The properties of the estimates have been studied in finite samples with a limited simulation study. This method is illustrated with an application to family data contributed to the Breast Cancer Consortium.

Algorithms↗

Mapping of complex traits by single-nucleotide polymorphisms.

Molecular geneticists are developing the third-generation human genome map with single-nucleotide polymorphisms (SNPs), which can be assayed via chip-based microarrays. One use of these SNP markers is the ability to locate loci that may be responsible for complex traits, via linkage/linkage-disequilibrium analysis. In this communication, we describe a semiparametric method for combined linkage/linkage-disequilibrium analysis using SNP markers. Asymptotic results are obtained for the estimated parameters, and the finite-sample properties are evaluated via a simulation study. We also applied this technique to a simulated genome-scan experiment for mapping a complex trait with two major genes. This experiment shows that separate linkage and linkage-disequilibrium analyses correctly detected the signals of both major genes; but the rates of false-positive signals seem high. When linkage and linkage-disequilibrium signals were combined, the analysis yielded much stronger and clearer signals for the presence of two major genes than did two separate analyses.

Biosensing Techniques↗

An efficient protocol for rare mutation genotyping in a large population.

We introduce a method to efficiently detect rare mutations for individual subjects in a large population by pooling samples and retesting subgroups of positive pooled samples. We conducted computer simulations of this method and discovered that it seems efficient for mutation prevalences less than 0.1, regardless of the number of samples. The simulations also indicate that splitting the pooled samples into three to five subgroups at each level is optimal. The expected number of necessary tests and relative efficiency of this method are given, by mutation prevalence and sample size.

Computer Simulation↗

Population-based family study designs: an interdisciplinary research framework for genetic epidemiology.

Most complex traits such as cancer and coronary heart diseases are attributed either to heritable factors or to environmental factors or to both. Dissecting the genetic and environmental etiology of complex traits thus requires an interdisciplinary research strategy. Genetic studies generally involve families and investigate familial aggregations of traits, segregation of major disease genes, and locations of disease genes on the human genome, the latter of which can be identified via linkage analysis. Epidemiologic studies often use population-based case-control studies to establish the role of specific environmental factors. Integrating both objectives, genetic epidemiology is to assess the associations of environmental factors with disease status, to quantify the aggregation of cases within families, to characterize putative disease genes via segregation analysis, and to localize disease genes via linkage analysis with genetic markers. To accomplish these objectives through designed studies, we propose a class of population-based family study designs, which are formed by choosing among sampling designs at three stages. The objectives of sampling at these three stages are 1) combined aggregation and association analysis, 2) combined segregation, aggregation, and association analysis, and 3) combined linkage, segregation, aggregation, and association analysis. These designs form an interdisciplinary research framework for genetic epidemiology. Our preliminary exploration of this framework and related analytic methods indicates that population-based family study designs retain the efficiency of linkage analysis for localizing disease genes without losing the property of being population-based, and they will therefore allow an assessment of a joint contribution of genetic and environmental factors to complex traits.

Case-Control Studies↗