The information bottleneck method for genome-wide association studies

In population studies, most current methods focus on identifying one outcome-related SNP at a time by testing for differences of genotype frequencies between disease and healthy groups or among different population groups. However, testing a great number of SNPs simultaneously has a problem of multiple testing and will give false-positive results. Although, this problem can be effectively dealt with through several approaches such as Bonferroni correction, permutation testing and false discovery rates, patterns of the joint effects by several genes, each with weak effect, might not be able to be determined. With the availability of high-throughput genotyping technology, searching for multiple scattered SNPs over the whole genome and modeling their joint effect on the target variable has become possible. Exhaustive search of all SNP subsets is computationally infeasible for millions of SNPs in a genome-wide study. Several effective feature selection methods combined with classification functions have been proposed to search for an optimal SNP subset among big data sets where the number of feature SNPs far exceeds the number of observations. ^ In this study, we take two steps to achieve the goal. First we selected 1000 SNPs through an effective filter method and then we performed a feature selection wrapped around a classifier to identify an optimal SNP subset for predicting disease. And also we developed a novel classification method-sequential information bottleneck method wrapped inside different search algorithms to identify an optimal subset of SNPs for classifying the outcome variable. This new method was compared with the classical linear discriminant analysis in terms of classification performance. Finally, we performed chi-square test to look at the relationship between each SNP and disease from another point of view. ^ In general, our results show that filtering features using harmononic mean of sensitivity and specificity(HMSS) through linear discriminant analysis (LDA) is better than using LDA training accuracy or mutual information in our study. Our results also demonstrate that exhaustive search of a small subset with one SNP, two SNPs or 3 SNP subset based on best 100 composite 2-SNPs can find an optimal subset and further inclusion of more SNPs through heuristic algorithm doesn't always increase the performance of SNP subsets. Although sequential forward floating selection can be applied to prevent from the nesting effect of forward selection, it does not always out-perform the latter due to overfitting from observing more complex subset states. ^ Our results also indicate that HMSS as a criterion to evaluate the classification ability of a function can be used in imbalanced data without modifying the original dataset as against classification accuracy. Our four studies suggest that Sequential Information Bottleneck(sIB), a new unsupervised technique, can be adopted to predict the outcome and its ability to detect the target status is superior to the traditional LDA in the study. ^ From our results we can see that the best test probability-HMSS for predicting CVD, stroke,CAD and psoriasis through sIB is 0.59406, 0.641815, 0.645315 and 0.678658, respectively. In terms of group prediction accuracy, the highest test accuracy of sIB for diagnosing a normal status among controls can reach 0.708999, 0.863216, 0.639918 and 0.850275 respectively in the four studies if the test accuracy among cases is required to be not less than 0.4. On the other hand, the highest test accuracy of sIB for diagnosing a disease among cases can reach 0.748644, 0.789916, 0.705701 and 0.749436 respectively in the four studies if the test accuracy among controls is required to be at least 0.4. ^ A further genome-wide association study through Chi square test shows that there are no significant SNPs detected at the cut-off level 9.09451E-08 in the Framingham heart study of CVD. Study results in WTCCC can only detect two significant SNPs that are associated with CAD. In the genome-wide study of psoriasis most of top 20 SNP markers with impressive classification accuracy are also significantly associated with the disease through chi-square test at the cut-off value 1.11E-07. ^ Although our classification methods can achieve high accuracy in the study, complete descriptions of those classification results(95% confidence interval or statistical test of differences) require more cost-effective methods or efficient computing system, both of which can't be accomplished currently in our genome-wide study. We should also note that the purpose of this study is to identify subsets of SNPs with high prediction ability and those SNPs with good discriminant power are not necessary to be causal markers for the disease.^

Psychosocial measures reported by parents in studies of skin cancer prevention

Background. Because it is important to minimize children's sun exposure to reduce skin cancer risk, much of the extensive skin cancer prevention literature consists of studies of children's sun protection, sun avoidance and ultraviolet radiation (UVR) exposure. Little attention has been focused on the measurement of psychosocial constructs in these studies. Identification of the psychosocial correlates or determinants of children's skin cancer risk or risk-reduction behavior is critical to more fully understand and predict behavior. Furthermore, psychosocial variables may be influenced by interventions to reduce risk. Thus, it is important to examine the psychosocial measures used in studies of children's skin cancer prevention. Information on the validity and reliability of psychosocial measures may increase confidence in study findings based on these measures. In particular, self-efficacy and barriers are key constructs in several major theoretical frameworks and parental measures have been associated with children's sun protection. However, there is conceptual overlap of self-efficacy and barriers measures and little is known about the psychometric properties of these measures.^ Study Aims and Methods. The overall goal of this dissertation was to examine the measurement of psychosocial constructs relevant to children's skin cancer prevention. Because children depend primarily on their parents for skin cancer prevention, measures of parents' psychosocial constructs are the focus. Study 1 was a systematic review of parental psychosocial measures used in studies of children's sun protection, sun avoidance and UVR exposure. The specific aims of Study 1 were to (1) describe psychosocial measures reported by parents, including available information on the psychometric properties of these measures and their use in analyses and (2) provide recommendations for the development, refinement and standardized reporting of measures. ^ Study 2 examined the psychometric properties of measures of parental self-efficacy and barriers regarding children's sun protection. Melanoma patients (N=205) who were parents of children ≤ 12 years of age completed a telephone interview that included self-efficacy and barriers measures specific to sunscreen, clothing, shade and limiting time outdoors. The specific aims of Study 2 were to (1) use a confirmatory factor analytic approach to examine the factorial validity of parental self-efficacy and barriers measures, (2) examine the convergent and discriminant validity of behavior-specific measures of self-efficacy and barriers and (3) assess the reliability of item and scale measures.^ Results. In Study 1, a search of standard databases yielded 48 eligible studies. Most studies assessed only one or two psychosocial constructs. Knowledge was measured most frequently. There was little discussion of measure source, development, theoretical background or psychometric properties, besides internal consistency reliability. There was conceptual overlap of some measures. In Study 2, confirmatory factor analytic findings supported the factorial validity of the self-efficacy and barriers measures. When all eight self-efficacy and barriers measures were included in the same model, a modified eight-factor model adequately fit the data, providing preliminary evidence that the measures are distinct. Measure associations supported the convergent validity of all measures and the discriminant validity of most measures. The self-efficacy and barriers measures were reliable.^ Conclusions. Recommendations based on the literature review include developing and refining psychosocial measures based on theory. Describing a measure's theoretical basis and psychometric properties would facilitate critical evaluation. Standardized reporting of source, development, theory, construct, items and analytic role would facilitate comparison of findings, continual refinement and future applications of measures. In the validation study, self-efficacy and barriers measures were examined in a sample of parents with a personal history of melanoma. Findings suggested that these measures are valid and reliable for use in studies of children's sun protection. There was preliminary evidence that these measures are distinct but additional study is needed. ^

Interim analysis of clinical trials: Simulation studies of fractional Brownian motion

Interim clinical trial monitoring procedures were motivated by ethical and economic considerations. Classical Brownian motion (Bm) techniques for statistical monitoring of clinical trials were widely used. Conditional power argument and α-spending function based boundary crossing probabilities are popular statistical hypothesis testing procedures under the assumption of Brownian motion. However, it is not rare that the assumptions of Brownian motion are only partially met for trial data. Therefore, I used a more generalized form of stochastic process, called fractional Brownian motion (fBm), to model the test statistics. Fractional Brownian motion does not hold Markov property and future observations depend not only on the present observations but also on the past ones. In this dissertation, we simulated a wide range of fBm data, e.g., H = 0.5 (that is, classical Bm) vs. 0.5< H <1, with treatment effects vs. without treatment effects. Then the performance of conditional power and boundary-crossing based interim analyses were compared by assuming that the data follow Bm or fBm. Our simulation study suggested that the conditional power or boundaries under fBm assumptions are generally higher than those under Bm assumptions when H > 0.5 and also matches better with the empirical results. ^

Macromolecular interaction studies on FF1-3 of human CA150 and STAT1-importin alpha5 using biophysical and biochemical methods

Macromolecular interactions, such as protein-protein interactions and protein-DNA interactions, play important roles in executing biological functions in cells. However the complexity of such interactions often makes it very challenging to elucidate the structural details of these subjects. In this thesis, two different research strategies were applied on two different two macromolecular systems: X-ray crystallography on three tandem FF domains of transcription regulator CA150 and electron microscopy on STAT1-importin α5 complex. The results from these studies provide novel insights into the function-structure relationships of transcription coupled RNA splicing mediated by CA150 and the nuclear import process of the JAK-STAT signaling pathway. ^ The first project aimed at the protein-protein interaction module FF domain, which often occurs as tandem repeats. Crystallographic structure of the first three FF domains of human CA150 was determined to 2.7 Å resolution. This is the only crystal structure of an FF domain and the only structure on tandem FF domains to date. It revealed a striking connectivity between an FF domain and the next. Peptide binding assay with the potential binding ligand of FF domains was performed using fluorescence polarization. Furthermore, for the first time, FF domains were found to potentially interact with DNA. DNA binding assays were also performed and the results were supportive to this newly proposed functionality of an FF domain. ^ The second project aimed at understanding the molecular mechanism of the nuclear import process of transcription factor STAT1. The first structural model of pSTAT1-importin α5 complex in solution was built from the images of negative staining electron microscopy. Two STAT1 molecules were observed to interact with one molecule of importin α5 in an asymmetric manner. This seems to imply that STAT1 interacts with importin α5 with a novel mechanism that is different from canonical importin α-cargo interactions. Further in vitro binding assays were performed to obtain more details on the pSTAT1-importin α5 interaction. ^

In vivo longitudinal studies of spinal cord injury and vascular endothelial growth factor treatment

Approximately 12,000 new cases of spinal cord injury (SCI) are added each year to the estimated 259,000 Americans living with SCI. The majority of these patients return to society, their lives forever changed by permanent loss of sensory and motor function. While there are no FDA approved drugs for the treatment of SCI or a universally accepted standard therapy, the current though controversial treatment includes the delivery of high dosages of the corticosteroid methyliprednisolone sodium succinate, surgical interventions to stabilize the spinal column, and physical rehabilitation. It is therefore critically important to fully understand the pathology of injury and determine novel courses and rationally-based therapies for SCI. ^ Vascular endothelial growth factor (VEGF) is an attractive target for treating central nervous system (CNS) injury and disease because it has been shown to influence angiogenesis and neuroprotection. Preliminary studies have indicated that increased vasculature may be associated with functional recovery; therefore exogenous delivery of a pro-angiogenic growth factor such as VEGF may improve neurobehavioral outcome. In addition, VEGF may provide protection from secondary injury and result in increased survival and axonal sprouting. ^ In these studies, SCI rats received acute intraspinal injections of VEGF, the antibody to VEGF, or vehicle control. The effect of these various agents was investigated using longitudinalmulti-modal magnetic resonance imaging (MRI), neuro- and sensory behavioral assays, and end point immunohistochemistry. We found that rats that received VEGF after SCI had increased tissue sparing and improved white matter integrity at the earlier time points as shown by advanced magnetic resonance imaging (MRI) techniques. However, these favorable effects of VEGF were not maintained, suggesting that additional treatments with VEGF at multiple time points may be more beneficial, Histological examinations revealed that VEGF treatment may result in increased oligodendrogenesis and therefore may eventually lead to remyelination and improved functional outcome. ^ On the neurobehavioral studies, treatments with VEGF and Anti-VEGF did not significantly affect performance on tests of open-field locomotion, grid walk, inclined plane, or rearing. However, VEGF treatment resulted in significantly increased incidence of chronic neuropathic pain. This phenomenon could possibly be attributed to the fact that VEGF treatment may promote axonal sprouting and also results in tissue sparing, thereby providing a substrate for the growth of new axons. New connections made by these sprouting axons may involve components of pathways involved in the transmission of pain and therefore result in increased pain in those animals. ^

STUDIES ON CHOLINERGIC FUNCTION AND MEMORY IN MICE (CHOLINOTOXINS, WORKING, REFERENCE MEMORY, RADIAL ARM MAZE)

Considerable evidence suggests that central cholinergic neurons participate in either acquisition, storage or retrieval of information. Experiments were designed to evaluate information processing in mice following either reversible or irreversible impairment in central cholinergic activity. The cholinergic receptor antagonists, atropine and methylatropine were used to reversibly inhibit cholinergic transmission. Irreversible impairment in central cholinergic function was achieved by central administration of the cholinergic-specific neurotoxins, N-ethyl-choline aziridinium (ECA) and N-ethyl-acetylcholine aziridinium (EACA).^ ECA and EACA appear to act by irreversible inhibition of high affinity choline uptake (proposed rate-limiting step in acetylcholine synthesis). Intraventricular administration of ECA or EACA produced persistent reduction in hippocampal choline acetyltransferase activity. Other neuronal systems and brain regions showed no evidence of toxicity.^ Mice treated with either ECA or EACA showed behavioral deficits associated with cholinergic dysfunction. Passive avoidance behavior was significantly impaired by cholinotoxin treatment. Radial arm maze performance was also significantly impaired in cholinotoxin-treated animals. Deficits in radial arm maze performance were transient, however, such that rapid and apparent complete behavioral recovery was seen during retention testing. The centrally active cholinergic receptor antagonist atropine also caused significant impairment in radial arm maze behavior, while equivalent doses of methylatropine were without effect.^ The relative effects of cholinotoxin and receptor antagonist treatment on short-term (working) memory and long-term (reference) memory in radial arm maze behavior were examined. Maze rotation studies indicated that there were at least two different response strategies which could result in accurate maze performance. One strategy involved the use of response algorithms and was considered to be a function of reference memory. Another strategy appeared to be primarily dependent on spatial working memory. However, all behavioral paradigms with multiple trails have reference memory requirements (i.e. information useful over all trials). Performance was similarly affected following either cholinotoxin or anticholinergic treatment, regardless of the response strategy utilized. In addition, rates of behavioral recovery following cholinotoxin treatment were similar between response groups. It was concluded that both cholinotoxin and anticholinergic treatment primarily resulted in impaired reference memory processes. ^

Functional data analysis approaches for genotype-phenotype association studies from next-generation sequencing

Next-generation DNA sequencing platforms can effectively detect the entire spectrum of genomic variation and is emerging to be a major tool for systematic exploration of the universe of variants and interactions in the entire genome. However, the data produced by next-generation sequencing technologies will suffer from three basic problems: sequence errors, assembly errors, and missing data. Current statistical methods for genetic analysis are well suited for detecting the association of common variants, but are less suitable to rare variants. This raises great challenge for sequence-based genetic studies of complex diseases.^ This research dissertation utilized genome continuum model as a general principle, and stochastic calculus and functional data analysis as tools for developing novel and powerful statistical methods for next generation of association studies of both qualitative and quantitative traits in the context of sequencing data, which finally lead to shifting the paradigm of association analysis from the current locus-by-locus analysis to collectively analyzing genome regions.^ In this project, the functional principal component (FPC) methods coupled with high-dimensional data reduction techniques will be used to develop novel and powerful methods for testing the associations of the entire spectrum of genetic variation within a segment of genome or a gene regardless of whether the variants are common or rare.^ The classical quantitative genetics suffer from high type I error rates and low power for rare variants. To overcome these limitations for resequencing data, this project used functional linear models with scalar response to develop statistics for identifying quantitative trait loci (QTLs) for both common and rare variants. To illustrate their applications, the functional linear models were applied to five quantitative traits in Framingham heart studies. ^ This project proposed a novel concept of gene-gene co-association in which a gene or a genomic region is taken as a unit of association analysis and used stochastic calculus to develop a unified framework for testing the association of multiple genes or genomic regions for both common and rare alleles. The proposed methods were applied to gene-gene co-association analysis of psoriasis in two independent GWAS datasets which led to discovery of networks significantly associated with psoriasis.^

Studies of biological rhythms in responses to influenza infection and vaccination

Three studies examined seasonal or circadian variations in selected responses to influenza infection or vaccination. The first, a seroepidemiologic study, evaluated temporal patterns of antibody titers to influenza A/Texas. Human umbilical cord bloods were sampled over a two-year period when the virus was not present in the community. No endogenous seasonal pattern was detected. The second study included three experiments on circadian rhythms in mice. Neither susceptibility nor protection from inactivated or attenuated vaccine varied significantly according to time of administration. A slight effect, however, was suggested with inactivated vaccine. Three human vaccine trials comprised the third study. Outcome variables included rise in antibody titer, final antibody titer, incidence of adverse reactions, and protection from community infection. Patterns in antibody response and protection variables were inconsistent, and generally not clinically significant. Local reactions to inactivated vaccine were more frequent if injections were received in the afternoon as compared to morning. This was true to adults that had been previously vaccinated. ^

Genetic studies of the human plasma orosomucoid (alpha-1 acid glycoprotein)

Orosomucoid (ORM) or alpha-1 acid glycoprotein is an acute phase protein of human plasma whose function is suggested to be the competitive inhibition of cellular recognition by infective agents. Isoelectric focusing (IEF) and immunoblotting have been combined and optimum conditions have been determined for reliable classification of different ORM phenotypes. Addition of 6 M urea in an IEF gel revealed additional microheterogeneity in the ORM system which has not been previously reported. 1,667 individuals from different native ethnic groups of North and South America, Africa and New Guinea have been screened to determine the distribution of ORM alleles. Two common alleles, ORM1*1 and ORM1*2 have been observed and their frequencies were determined. Genetically independent variation consistent with expression of the ORM2 locus was observed in American and African blacks but was not observed in other sampled populations. The population allele frequencies for this new locus were 0.958, 0.025, 0.006, 0.011, for alleles ORM2*1, ORM2*2, ORM2*3, ORM2*4, respectively. Family studies confirm the autosomal codominant inheritance of the phenotypes observed at both ORM loci. ^

Design of mixed-longitudinal studies

Cross-sectional designs, longitudinal designs in which a single cohort is followed over time, and mixed-longitudinal designs in which several cohorts are followed for a shorter period are compared by their precision, potential for bias due to age, time and cohort effects, and feasibility. Mixed longitudinal studies have two advantages over longitudinal studies: isolation of time and age effects and shorter completion time. Though the advantages of mixed-longitudinal studies are clear, choosing an optimal design is difficult, especially given the number of possible combinations of the number of cohorts and number of overlapping intervals between cohorts. The purpose of this paper is to determine the optimal design for detecting differences in group growth rates.^ The type of mixed-longitudinal study appropriate for modeling both individual and group growth rates is called a "multiple-longitudinal" design. A multiple-longitudinal study typically requires uniform or simultaneous entry of subjects, who are each observed till the end of the study.^ While recommendations for designing pure-longitudinal studies have been made by Schlesselman (1973b), Lefant (1990) and Helms (1991), design recommendations for multiple-longitudinal studies have never been published. It is shown that by using power analyses to determine the minimum number of occasions per cohort and minimum number of overlapping occasions between cohorts, in conjunction with a cost model, an optimal multiple-longitudinal design can be determined. An example of systolic blood pressure values for cohorts of males and cohorts of females, ages 8 to 18 years, is given. ^

The use of the computerized and revised vital record in conducting etiologic studies on neural tube defects

The purpose of this study was to evaluate the adequacy of computerized vital records in Texas for conducting etiologic studies on neural tube defects (NTDs), using the revised and expanded National Centers for Health Statistics vital record forms introduced in Texas in 1989.^ Cases of NTDs (anencephaly and spina bifida) among Harris County (Houston) residents were identified from the computerized birth and death records for 1989-1991. The validity of the system was then measured against cases ascertained independently through medical records and death certificates. The computerized system performed poorly in its identification of NTDs, particularly for anencephaly, where the false positive rate was 80% with little or no improvement over the 3-year period. For both NTDs the sensitivity and predictive value positive of the tapes were somewhat higher for Hispanic than non-Hispanic mothers.^ Case control studies were conducted utilizing the tape set and the independently verified data set, using controls selected from the live birth tapes. Findings varied widely between the data sets. For example, the anencephaly odds ratio for Hispanic mothers (vs. non-Hispanic) was 1.91 (CI = 1.38-2.65) for the tape file, but 3.18 (CI = 1.81-5.58) for verified records. The odds ratio for diabetes was elevated for the tape set (OR = 3.33, CI = 1.67-6.66) but not for verified cases (OR = 1.09, CI = 0.24-4.96), among whom few mothers were diabetic. It was concluded that computerized tapes should not be solely relied on for NTD studies.^ Using the verified cases, Hispanic mother was associated with spina bifida, and Hispanic mother, teen mother, and previous pregnancy terminations were associated with anencephaly. Mother's birthplace, education, parity, and diabetes were not significant for either NTD.^ Stratified analyses revealed several notable examples of statistical interaction. For anencephaly, strong interaction was observed between Hispanic origin and trimester of first prenatal care.^ The prevalence was 3.8 per 10,000 live births for anencephaly and 2.0 for spina bifida (5.8 per 10,000 births for the combined categories). ^

Planning and implementation efficacy as a determinant of health risk communication outcomes: Case studies of the U.S. Public Health Service (PHS)

This study was conducted under the auspices of the Subcommittee on Risk Communication and Education of the Committee to Coordinate Environmental Health and Related Programs (CCEHRP) to determine how Public Health Service (PHS) agencies are communicating information about health risk, what factors contributed to effective communication efforts, and what specific principles, strategies, and practices best promote more effective health risk communication outcomes.^ Member agencies of the Subcommittee submitted examples of health risk communication activities or decisions they perceived to be effective and some examples of cases they thought had not been as effective as desired. Of the 10 case studies received, 7 were submitted as examples of effective health risk communication, and 3, as examples of less effective communication.^ Information contained in the 10 case studies describing the respective agencies' health risk communication strategies and practices was compared with EPA's Seven Cardinal Rules of Risk Communication, since similar rules were not found in any PHS agency. EPA's rules are: (1) Accept and involve the public as a legitimate partner. (2) Plan carefully and evaluate your efforts. (3) Listen to the public's specific concerns. (4) Be honest, frank, and open. (5) Coordinate and collaborate with other credible sources. (6) Meet the needs of the media. (7) Speak clearly and with compassion.^ On the basis of case studies analysis, the Subcommittee, in their attempts to design and implement effective health risk communication campaigns, identified a number of areas for improvement among the agencies. First, PHS agencies should consider developing a focus specific to health risk communication (i.e., office or specialty resource). Second, create a set of generally accepted practices and guidelines for effective implementation and evaluation of PHS health risk communication activities and products. Third, organize interagency initiatives aimed at increasing awareness and visibility of health risk communication issues and trends within and between PHS agencies.^ PHS agencies identified some specific implementation strategies the CCEHRP might consider pursuing to address the major recommendations. Implementation strategies common to PHS agencies emerged in the following five areas: (1) program development, (2) building partnerships, (3) developing training, (4) expanding information technologies, and (5) conducting research and evaluation. ^

Studies of the expression of DNA repair related genes

This dissertation examines the biological functions and the regulation of expression of DNA ligase I by studying its expression under different conditions.^ The gene expression of DNA ligase I was induced two- to four-fold in S-phase lymphoblastoid cells but was decreased to 15% of control after administration of a DNA damaging agent, 4-nitroquinoline-1-oxide. When cells were induced into differentiation, the expression level of DNA ligase I was decreased to less than 15% of that of the control cells. When the gene of DNA ligase I was examined for tissue specific expression in adult rats, high levels of DNA ligase I mRNA were observed in testis (8-fold), intermediate levels in ovary and brain (4-fold), and low levels were found in intestine, spleen, and liver (1- to 2-fold).^ In confluent cells of normal skin fibroblasts, UV irradiation induced the gene expression of DNA ligase I at 24 and 48 h. The induction of DNA ligase I gene expression requires active p53 protein. Introducing a vector containing the wild type p53 protein in the cells caused an induction of the DNA ligase I protein 24 h after the treatment.^ Our results indicate that, in addition to the regulation by phosphorylation/dephosphorylation, cellular DNA ligase I activity can be regulated at the gene transcription level, and the p53 tumor suppresser is one of the transcription factors for the DNA ligase I gene. Also, our results suggest that DNA ligase I is involved in DNA repair as well as in DNA replication.^ Also, as an early attempt to clone the human homolog of the yeast CDC9 gene which has been shown to be involved in DNA replication, DNA repair, and DNA recombination, we have identified a human gene with mRNA of 1.7 kb. This dissertation studies the gene regulation and the possible biological functions of this new human gene by examining its expression at different stages of the cell cycle, during cell differentiation, and in cellular response to DNA damage.^ The new gene that we recently identified from human cells is highly expressed in brain and reproductive organs (BRE). This BRE gene encodes an mRNA of 1.7-1.9 kb, with an open reading frame of 1,149 bp, and gives rise to a deduced polypeptide of 383 amino acid residues. No extensive homology was found between BRE and sequences from the EMBL-Gene Banks. BRE showed tissue-specific expression in adult rats. The steady state mRNA levels were high in testis (5-6 fold), ovary and brain (3-4 fold) compared to the spleen level, but low in intestine and liver (1-2 fold). The expression of this gene is responsive to DNA damage and/or retinoic acid (RA) treatment. Treatment of fibroblast cells with UV irradiation and 4-nitroquinoline-1-oxide caused more than 90% and 50% decreases in BRE mRNA, respectively. Similar decreases in BRE expression were observed after treatment of the brain glioma cell line U-251 and the promyelocytic cell line HL-60 with retinoic acid. (Abstract shortened by UMI). ^

THE USE OF MULTIVARIATE ANALYSIS OF COVARIANCE IN EVALUATING THE COMPARABILITY OF LABORATORY DETERMINATIONS IN COOPERATIVE STUDIES

The role of clinical chemistry has traditionally been to evaluate acutely ill or hospitalized patients. Traditional statistical methods have serious drawbacks in that they use univariate techniques. To demonstrate alternative methodology, a multivariate analysis of covariance model was developed and applied to the data from the Cooperative Study of Sickle Cell Disease.^ The purpose of developing the model for the laboratory data from the CSSCD was to evaluate the comparability of the results from the different clinics. Several variables were incorporated into the model in order to control for possible differences among the clinics that might confound any real laboratory differences.^ Differences for LDH, alkaline phosphatase and SGOT were identified which will necessitate adjustments by clinic whenever these data are used. In addition, aberrant clinic values for LDH, creatinine and BUN were also identified.^ The use of any statistical technique including multivariate analysis without thoughtful consideration may lead to spurious conclusions that may not be corrected for some time, if ever. However, the advantages of multivariate analysis far outweigh its potential problems. If its use increases as it should, the applicability to the analysis of laboratory data in prospective patient monitoring, quality control programs, and interpretation of data from cooperative studies could well have a major impact on the health and well being of a large number of individuals. ^