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. ^

ISOELECTRIC FOCUSING STUDIES OF STEROID BINDING PROTEINS

Steroid binding proteins are an obvious choice in the search for genetic factors in plasma that might predispose to upper body obesity, a risk factor for non-insulin dependent diabetes and cardiovascular disease. The two steroid binding proteins studied by isoelectric focusing were sex hormone binding globulin (SHBG), the transport protein for sex hormones and corticosteroid binding globulin (CBG), the transport protein for corticosteroids. Auto-radiography and immunoblotting on polyacrylamide gels were used to detect polymorphism in SHBG. Immunoblotting on agarose gels was used to visualize corticosteroid binding globulin. SHBG showed similar structural variation in American Caucasians, American Blacks and Canadian Indians. Two alleles (1, 2) were hypothesized with highly polymorphic frequencies in all three ethnic groups. CBG was not found to be polymorphic, but two variants were found in Caucasian male twins and in a Black individual. The finding of a good assay and a polymorphic system for SHBG are the first steps for additional studies into disease associations. ^

RELATION OF TYPE OF REPORTED ILLNESSES TO PATTERN OF USE OF THE STUDENT HEALTH CENTER BY THE INTERNATIONAL STUDENTS AT TEXAS SOUTHERN UNIVERSITY, HOUSTON, TEXAS

The state of knowledge on the relation of stress factors, health problems and health service utilization among university students is limited. Special problems of stress exist for the international students due to their having to adjust to a new environment. It is this latter problem area that provides the focus for this study. Recognizing there are special stress factors affecting the international students, it is first necessary to see if the problems of cultural adaptation affect them to any greater degree than American students attending the same university.^ To make the comparison, the study identified a number of health problems of both American and international students and related their frequency to the use of the Student Health Center. The expectation was that there would be an association between the number of health problems and the number of life change events experienced by these students and between the number of health problems and stresses from social factors. It was also expected that the number of health problems would decline with the amount of social support.^ The population chosen were students newly enrolled in Texas Southern University, Houston, Texas in the Fall Semester of 1979. Two groups were selected at random: 126 international and 126 American students. The survey instrument was a self-administered questionnaire. The response rate was 90% (114) for the international and 94% (118) for the American students.^ Data analyses consisted of both descriptive and inferential statistics. Chi-squares and correlation coefficients were the statistics used in comparing the international students and the American students.^ There was a weak association between the number of health problems and the number of life change events, as reported by both the international and the American students. The study failed to show any statistically significant association between the number of stress from social factors and the number of health problems. It also failed to show an association between the number of health problems and the amount of social support. These findings applied to both the international and the American students.^ One unexpected finding was that certain health problems were reported by more American than international students. There were: cough, diarrhea, and trouble in sleeping. Another finding was that those students with health insurance had a higher level of utilization of the Health Center than those without health insurance. More international than American students utilized the Student Health Center.^ In comparing the women students, there was no statistical significant difference in their reported fertility related health problems.^ The investigator recommends that in follow-up studies, instead of grouping all international students together, that they be divided by major nationalities represented in the student body; that is, Iranians, Nigerians and others. ^

DEVELOPMENT AND APPLICATION OF SHORT-TERM MUTAGENICITY AND CYTOTOXICITY BIOASSAYS FOR INTEGRATED HEALTH ASSESSMENT STUDIES OF COAL FLY ASH EMISSIONS

Two respirable coal fly ash samples ((LESSTHEQ) 3(mu)m), one from a pressurized fluidized-bed combustion miniplant and one from a conventional combustion power plant, were investigated for physical properties, chemical composition and biological activity. Electron microscopy illustrated irregularity in fluidized-bed combustion fly ash and sphericity in conventional combustion fly ash. Elemental analysis of these samples showed differences in trace elements. Both fly ash samples were toxic in rabbit alveolar macrophage and Chinese hamster ovary cell systems in vitro. The macrophages were more sensitive to toxicity of fly ash than the ovary cells. For measuring the cytotoxicity of fly ash, the most sensitive parameters were adenosine triphosphate in the alveolar macrophage system and viability index in the hamster ovary system. Intact fluidized-bed combustion fly-ash particles showed mutagenicity only in strains TA98 and TA1538 without metabolic activation in the Ames Salmonella assay. No mutagenicity was detected in bioassay of conventional combustion fly ash particles. Solvent extraction yielded more mass from fluidized-bed combustion fly ash than from conventional combustion fly ash. The extracts of fluidized-bed combustion fly ash showed higher mutagenic activity than conventional combustion fly ash. These samples contained direct-acting, frameshift mutagens.^ Fly ash samples collected from the same fluidized-bed source by cyclones, a fabric filter, and a electrostatic precipitator at various temperatures were compared for particle size, toxicity, and mutagenicity. Results demonstrated that the biological activity of coal fly ash were affected by the collection site, device, and temperature.^ Coal fly ash vapor-coated with 1-nitropyrene was developed as a model system to study the bioavailability and recovery of nitroaromatic compounds in fly ash. The effects of vapor deposition on toxicity and mutagenicity of fly ash were examined. The nitropyrene coating did not significantly alter the ash's cytotoxicity. Nitropyrene was bioavailable in the biological media, and a significant percentage was not recovered after the coated fly ash was cultured with alveolar macrophages. 1-Nitropyrene loss increased as the number of macrophages was increased, suggesting that the macrophages are capable of metabolizing or binding 1-nitropyrene present in coal fly ash. ^

REGRESSION ANALYSIS FOR STANDARDIZED MORTALITY RATIOS WITH APPLICATION TO OCCUPATIONAL FOLLOW-UP STUDIES

Traditional comparison of standardized mortality ratios (SMRs) can be misleading if the age-specific mortality ratios are not homogeneous. For this reason, a regression model has been developed which incorporates the mortality ratio as a function of age. This model is then applied to mortality data from an occupational cohort study. The nature of the occupational data necessitates the investigation of mortality ratios which increase with age. These occupational data are used primarily to illustrate and develop the statistical methodology.^ The age-specific mortality ratio (MR) for the covariates of interest can be written as MR(,ij...m) = ((mu)(,ij...m)/(theta)(,ij...m)) = r(.)exp (Z('')(,ij...m)(beta)) where (mu)(,ij...m) and (theta)(,ij...m) denote the force of mortality in the study and chosen standard populations in the ij...m('th) stratum, respectively, r is the intercept, Z(,ij...m) is the vector of covariables associated with the i('th) age interval, and (beta) is a vector of regression coefficients associated with these covariables. A Newton-Raphson iterative procedure has been used for determining the maximum likelihood estimates of the regression coefficients.^ This model provides a statistical method for a logical and easily interpretable explanation of an occupational cohort mortality experience. Since it gives a reasonable fit to the mortality data, it can also be concluded that the model is fairly realistic. The traditional statistical method for the analysis of occupational cohort mortality data is to present a summary index such as the SMR under the assumption of constant (homogeneous) age-specific mortality ratios. Since the mortality ratios for occupational groups usually increase with age, the homogeneity assumption of the age-specific mortality ratios is often untenable. The traditional method of comparing SMRs under the homogeneity assumption is a special case of this model, without age as a covariate.^ This model also provides a statistical technique to evaluate the relative risk between two SMRs or a dose-response relationship among several SMRs. The model presented has application in the medical, demographic and epidemiologic areas. The methods developed in this thesis are suitable for future analyses of mortality or morbidity data when the age-specific mortality/morbidity experience is a function of age or when there is an interaction effect between confounding variables needs to be evaluated. ^