A computerized statistical framework for coalescent analysis

Coalescent theory represents the most significant progress in theoretical population genetics in the past three decades. The coalescent theory states that all genes or alleles in a given population are ultimately inherited from a single ancestor shared by all members of the population, known as the most recent common ancestor. It is now widely recognized as a cornerstone for rigorous statistical analyses of molecular data from population [1]. The scientists have developed a large number of coalescent models and methods[2,3,4,5,6], which are not only applied in coalescent analysis and process, but also in today’s population genetics and genome studies, even public health. The thesis aims at completing a statistical framework based on computers for coalescent analysis. This framework provides a large number of coalescent models and statistic methods to assist students and researchers in coalescent analysis, whose results are presented in various formats as texts, graphics and printed pages. In particular, it also supports to create new coalescent models and statistical methods. ^

A cohort study of occupational asbestos-exposure related neoplasms in Texas Gulf Coast area

A cohort study was conducted in Texas and Louisiana Gulf Coast area on individual workers who have been exposed to asbestos for 15 years or more. Most of these workers were employed in petrochemical industries. Of the 15,742 subjects initially selected for the cohort study, 3,258 had positive chest X-ray findings believed to be related to prolonged asbestos exposure. These subjects were further investigated. Their work out included detailed medical and occupational history, laboratory tests and spirometry. One thousand eight-hundred and three cases with positive chest X-ray findings whose data files were considered complete at the end of May 1986 were analyzed and their findings included in this report.^ The prevalence of lung cancer and cancer of the following sights: skin, stomach, oropharyngeal, pancreas and kidneys were significantly increased when compared to data from Connecticut Tumor Registry. The prevalence of other chronic conditions such as hypertension, emphysema, heart disease and peptic ulcer was also significantly high when compared to data for the U.S. and general population furnished by the National Center for Health Statistics (NCHS). In most instances the occurrence of cancer and the chronic ailment previously mentioned appeared to follow 15-25 years of exposure to asbestos. ^

Effects of asbestos and silica on superoxide anion production in the guinea pig alveolar macrophage

Asbestos and silica are important industrial hazards. Exposure to these dusts can result in pulmonary fibrosis and, in the case of asbestos, cancer. Although the hazards of asbestos and silica exposure have long been known, the pathogenesis of dust-related disease is not well understood. Both silica and asbestos are thought to alter the function of the alveolar macrophage, but the nature of the biochemical alteration is unknown. Therefore, this study examined the effect of asbestos and silica on the activation pathway of the guinea pig alveolar macrophage. Activation of macrophages by physiological agents results in stimulation of phospholipase C causing phosphatidyl inositol turnover and intracellular calcium mobilization. Phosphatidyl inositol turnover produces diacylglycerol which activates protein kinase C causing superoxide anion production.^ Chrysotile stimulated alveolar macrophages to produce superoxide anion. This stimulation proceeded via phospholipase C, since chrysotile stimulated phosphatidyl inositol turnover and intracellular calcium mobilization. The possible involvement of a coupling protein was evaluated by pretreating cells with pertussis toxin. Pertussis toxin pretreatment partially inhibited chrysotile stimulation, suggesting that chrysotile activates a coupling protein in an non-classical manner. Potential binding sites for chrysotile stimulation were examined using a series of nine lectins. Chrysotile-stimulated superoxide anion production was blocked by pretreatment with lectins which bound to N-acetylglucosamine, but not by lectins which bound to mannose, fucose, or N-acetylgalactosamine. In addition, incubation with the N-acetylglucosamine polymer, chitin, inhibited chrysotile-stimulated superoxide anion production, suggesting that chrysotile stimulated superoxide anion production by binding to N-acetylglucosamine residues.^ On the other hand, silica did not stimulate superoxide anion production. The effect of silica on agonist stimulation of this pathway was examined using two stimulants of superoxide anion production, N-formyl-nle-leu-phe (FNLP, which stimulates through phospholipase C) and phorbol-12,13-dibutyrate (which directly activates protein kinase C). Sublethal doses of silica inhibited FNLP-stimulated superoxide anion production, but did not affect phorbol-12,13-dibutyrate-stimulated superoxide anion production, suggesting that the site of inhibition precedes protein kinase C. This inhibition was not due to cell membrane damage, since cell permeability to calcium-45 and rubidium-86 was not increased. It is concluded that chrysotile binds to N-acetylglucosamine residues on macrophage surface glycoproteins to stimulate the physiological pathway resulting in superoxide anion production. In contrast, silica does not stimulate superoxide anion production, but it did inhibit FNLP-stimulated superoxide anion production. ^

An evaluation of sleep hypnotic and stimulant effect on human performance using a computerized psychological test battery

This study evaluated the administration-time-dependent effects of a stimulant (Dexedrine 5-mg), a sleep-inducer (Halcion 0.25-mg) and placebo (control) on human performance. The investigation was conducted on 12 diurnally active (0700-2300) male adults (23-38 yrs) using a double-blind, randomized sixway-crossover three-treatment, two-timepoint (0830 vs 2030) design. Performance tests were conducted hourly during sleepless 13-hour studies using a computer generated, controlled and scored multi-task cognitive performance assessment battery (PAB) developed at the Walter Reed Army Institute of Research. Specific tests were Simple and Choice Reaction Time, Serial Addition/Subtraction, Spatial Orientation, Logical Reasoning, Time Estimation, Response Timing and the Stanford Sleepiness Scale. The major index of performance was "Throughput", a combined measure of speed and accuracy.^ For the Placebo condition, Single and Group Cosinor Analysis documented circadian rhythms in cognitive performance for the majority of tests, both for individuals and for the group. Performance was best around 1830-2030 and most variable around 0530-0700 when sleepiness was greatest (0300).^ Morning Dexedrine dosing marginally enhanced performance an average of 3% with reference to the corresponding in time control level. Dexedrine AM also increased alertness by 10% over the AM control. Dexedrine PM failed to improve performance with reference to the corresponding PM control baseline. With regard to AM and PM Dexedrine administrations, AM performance was 6% better with subjects 25% more alert.^ Morning Halcion administration caused a 7% performance decrement and 16% increase in sleepiness and a 13% decrement and 10% increase in sleepiness when administered in the evening compared to corresponding in time control data. Performance was 9% worse and sleepiness 24% greater after evening versus morning Halcion administration.^ These results suggest that for evening Halcion dosing, the overnight sleep deprivation occurring in coincidence with the nadir in performance due to circadian rhythmicity together with the CNS depressant effects combine to produce performance degradation. For Dexedrine, morning administration resulted in only marginal performance enhancement; Dexedrine in the evening was less effective, suggesting the 5-mg dose level may be too low to counteract the partial sleep deprivation and nocturnal nadir in performance. ^

Cytogenetic assessment of occupational exposure to antineoplastic agents

This study has evaluated the utility of measuring effects of low level occupational exposure of nursing personnel to antineoplastic agents. The effect measured in this study is chromosomal damage in peripheral lymphocytes (chromosomal breakage and micronuclei frequency).^ Using nursing personnel in three exposure classifications (low, moderate and high) and breast cancer patients before and after treatment with antineoplastic agents, a weak but statistically significant association was found between exposure and chromosomal damage. Of special interest was the finding that consistent glove usage was negatively associated with increased chromosomal damage.^ The study also demonstrated a statistically significant association between the two measures of chromosomal damage: chromosomal breakage and micronuclei frequency. This suggests that the micronucleus method is a useful test for studying cytogenetic effects in lymphocytes. ^

Genetic epidemiology of childhood brain tumors

A rare familial cancer syndrome involving childhood brain tumors (CBT), breast cancer, sarcomas and an array of other tumors has been described (Li and Fraumeni 1969, 1975, 1982, 1987). A survey of CBT identified through the Connnecticut Tumor Registry in 1984 revealed a high frequency of CBT, leukemia and other childhood cancer in siblings of CBT patients (Farwell and Flannery, 1984). Other syndromes such as neurofibromatosis and nevoid basal cell carcinoma syndrome have also been associated with CBT; however, no systematic family studies have been conducted to determine the extent to which cancer aggregates in family members of CBT patients. This family study was designed to determine the frequency of cancer aggregation overall or at specific sites, to determine the frequency of known or potentially hereditary syndromes in families of CBT patients, and to determine a genetic model to characterize familial cancer syndromes and to identify specific kindreds to which such a model(s) might apply. This study includes 244 confirmed CBT patients referred to the University of Texas M. D. Anderson Cancer Center between the years 1944 and 1983, diagnosed under the age of 15 years and resident in the U.S. or Canada. Family histories were obtained on the proband's first (parents, siblings and offspring) and second degree (proband's aunts, uncles and grandparents) relatives following sequential sampling scheme rules. To determine if cancer aggregates in families, we compared the cancer experience in the population to that expected in the general population using Connecticut Tumor Registry calendar year, age, race and sex-specific rates. The standardized incidence ratio (SIR) for cancer overall was 0.91 (41 observed (O) and 44.94 expected (E); 95% Confidence Interval (CI) = 0.65-1.24). We observed a significant excess of colon cancer among the proband's first degree relatives (O/E = 5/1.64; 95% CI = 1.01-7.65), in particular those under age 45 year. Segregation analysis showed evidence for multifactorial inheritance in the small percentage (N = 5) of the families. ^

A multiplicative model for standardizing vital rates without using external standards

Standardization is a common method for adjusting confounding factors when comparing two or more exposure category to assess excess risk. Arbitrary choice of standard population in standardization introduces selection bias due to healthy worker effect. Small sample in specific groups also poses problems in estimating relative risk and the statistical significance is problematic. As an alternative, statistical models were proposed to overcome such limitations and find adjusted rates. In this dissertation, a multiplicative model is considered to address the issues related to standardized index namely: Standardized Mortality Ratio (SMR) and Comparative Mortality Factor (CMF). The model provides an alternative to conventional standardized technique. Maximum likelihood estimates of parameters of the model are used to construct an index similar to the SMR for estimating relative risk of exposure groups under comparison. Parametric Bootstrap resampling method is used to evaluate the goodness of fit of the model, behavior of estimated parameters and variability in relative risk on generated sample. The model provides an alternative to both direct and indirect standardization method. ^

Mitral valve prolapse and its relationship with stroke

The purpose of this study was to elucidate the relationship between mitral valve prolapse and stroke. A population-based historical cohort investigation was conducted among residents of Olmsted County, Minnesota who had an initial echocardiographic diagnosis of mitral valve prolapse from 1975 through 1989. This cohort (N = 1085) was followed for stroke outcomes using the resources of an operational medical record linkage system. There was an overall two-fold increase in the incidence of stroke among individuals with mitral valve prolapse relative to a standard population (standardized morbidity ratio = 2.12, 95% confidence limits = 1.33-3.21). When the data were partitioned by duration of follow-up from the diagnosis of mitral valve prolapse, or by the calendar years at echocardiographic diagnosis, respectively, the association between mitral valve prolapse and stroke was not modified. Mitral valve prolapse subjects 85 years and older were at highest increased risk of developing strokes relative to the general population (standardized morbidity ratio = 5.47, 95% confidence limits = 2.20-11.24). Coronary heart disease, atrial fibrillation, diabetes mellitus and hypertension, were unlikely to have confounded the association between mitral valve prolapse and stroke.^ The cumulative risk of first stroke among individuals initially diagnosed with mitral valve prolapse age 15 to 64 years, given survival to 15.2 years of follow-up, was 4.0%. The cumulative risk of first stroke among individuals initially diagnosed with mitral valve prolapse age 65 to 74 years, given survival to 11.2 years of follow-up, was 13.2%. The cumulative risk of first stroke among individuals initially diagnosed with mitral valve prolapse age 75 years and older, given survival to 6.7 years of follow-up, was 30.6%.^ Among individuals with mitral valve prolapse, age, diabetes, and atrial fibrillation were associated with an increased risk of stroke. Atrial fibrillation was associated with a four-fold rate of stroke and diabetes associated with a seven-fold rate of stroke.^ Findings from this research support the hypothesis that mitral valvular heart prolapse is linked with a stroke sequela. ^

Human DNA ligase and DNA polymerase as molecular targets for heavy metals and anticancer drugs

DNA ligase and DNA polymerase play important roles in DNA replication, repair, and recombination. Frequencies of spontaneous and chemical- and physical-induced mutations are correlated to the fidelity of DNA replication. This dissertation elucidates the mechanisms of the DNA ligation reaction by DNA ligases and demonstrates that human DNA ligase I and DNA polymerase $\alpha$ are the molecular targets for two metal ions, Zn$\sp{2+}$ and Cd$\sp{2+},$ and an anticancer drug, F-ara-ATP.^ Human DNA ligases were purified to homogeneity and their AMP binding domains were mapped. Although their AMP-binding domains are similar, there could be difference between the two ligases in their DNA binding domains.^ The formation of the AMP-DNA intermediate and the successive ligation reaction by human DNA ligases were analyzed. Both reactions showed their substrate specificity for ligases I and II, required Mg2+, and were inhibited by ATP.^ A protein inhibitor from HeLa cells and specific for human DNA ligase I but not ligase II and T4 ligase was discovered. It reversibly inhibited DNA ligation activity but not the AMP-binding activity due to the formation of a reversible ligase I-inhibitor complex.^ F-ara-ATP inhibited human DNA ligase I activity by competing with ATP for the AMP-binding site of DNA ligase I, forming a ligase I-F-ara-AMP complex, as well as when it was incorporated at 3$\sp\prime$-terminus of DNA nick by DNA polymerase $\alpha.$^ All steps of the DNA ligation reaction were inhibited by Zn$\sp{2+}$ and Cd$\sp{2+}$ in a concentration-dependent manner. Both ions did not show the ability to change the fidelity of DNA ligation reaction catalyzed by human DNA ligase I. However, Zn$\sp{2+}$ and Cd$\sp{2+}$ showed their contradictory effects on the fidelity of the reaction by human DNA polymerase $\alpha.$ Zn$\sp{2+}$ decreased the frequency of misinsertion but less affected that of mispair extension. On the contrary, Cd$\sp{2+}$ increased the frequencies of both misinsertion and mispair extension at very low concentration. Our data provided strong evidence in the molecular mechanisms for the mutagenicity of zinc and cadmium, and were comparable with the results previously reported. ^

The effect of life events on blood pressure and body mass index in children

In Conroe, Texas, 492 students ages 5 to 15 participated in a screening examination for cardiovascular risk factor study. Among 492 students, 141 elementary and junior high students participated in the present sub-study to investigate the effect of the number of recent life events on blood pressure and on body mass index. Using the elementary and junior high school Coddington scales, life events occurring in the past 12 months were measured for students ages 9 to 14 years, no significant differences in life events were observed by age and sex. The number of life events was not related to blood pressure but was positively correlated to body mass index in children and adolescents. ^

Modification of mutagen sensitivity by chemopreventive agents {\it in vitro\/}

The clinical application of chemopreventive agents is expected to prevent the appearance of cancer by arresting carcinogenesis or reversing it in the precancerous stages. The hypothesis of the present investigations was that chemopreventive agents (retinoids and antioxidant vitamins) may counteract the clastogenic effects of bleomycin in vitro in both lymphoblastoid cell lines and primary lymphocyte cultures and that a similar phenomenon can be detected in lymphocytes from individuals treated with 13-cis-retinoic acid. The efficacy of 13-cis-retinoic acid, n-(4-hydroxyphenyl)-retinamide, ascorbic acid, n-acetyl-l-cysteine, alpha-tocopherol, and alpha-tocopherol-acid succinate was tested against bleomycin-induced chromosomal breakage.^ The results provided direct evidence of the concentration-related protective effects of these agents against bleomycin-induced clastogenicity in cultures of human lymphoblastoid cell lines in vitro. Similar anticlastogenic protection was demonstrated with 13-cis-retinoic acid, ascorbic acid, n-acetyl-l-cysteine, and alpha-tocopherol-acid succinate in primary lymphocyte cultures in vitro. The in vitro anticlastogenic effect of 13-cis-retinoic acid was also demonstrated in lymphocyte cultures from peripheral blood samples from patients treated with this retinoid.^ An important consideration is that the concentrations used in the present investigations are comparable to those achieved in clinical situations.^ The in vitro anticlastogenic effect of these retinoids and antioxidants may constitute an important element of their chemopreventive properties. The results corroborate the hypothesis that these compounds may be effective in clinical chemoprevention trials. The bleomycin-assay may also be used as a short-term test to evaluate the antimutagenic effects of various agents. ^

Modeling the survival of T-cell lymphocytes using compound branching processes

Radiotherapy has been a method of choice in cancer treatment for a number of years. Mathematical modeling is an important tool in studying the survival behavior of any cell as well as its radiosensitivity. One particular cell under investigation is the normal T-cell, the radiosensitivity of which may be indicative to the patient's tolerance to radiation doses.^ The model derived is a compound branching process with a random initial population of T-cells that is assumed to have compound distribution. T-cells in any generation are assumed to double or die at random lengths of time. This population is assumed to undergo a random number of generations within a period of time. The model is then used to obtain an estimate for the survival probability of T-cells for the data under investigation. This estimate is derived iteratively by applying the likelihood principle. Further assessment of the validity of the model is performed by simulating a number of subjects under this model.^ This study shows that there is a great deal of variation in T-cells survival from one individual to another. These variations can be observed under normal conditions as well as under radiotherapy. The findings are in agreement with a recent study and show that genetic diversity plays a role in determining the survival of T-cells. ^

Fate and effects of drilling wastes in a shallow estuarine mesocosm

To address concerns expressed about the possible effect of drilling mud discharges on shallow, low-energy estuarine ecosystems, a 12 month study was designed to detect alterations in water quality and sediment geochemistry. Each drilling mud used in the study and sediments from the study site were analyzed in the laboratory for chemical and physical characteristics. Potential water quality impacts were simulated by the EPA-COE elutriation test procedure. Mud toxicity was measured by acute and chronic bioassays with Mysidopsis bahia, Mercenaria mercenaria, and Nereis virens.^ For the field study, a relatively pristine, shallow (1.2 m) estuary (Christmas Bay, TX) without any drilling activity for the last 30 years was chosen for the study site. After a three month baseline study, three stations were selected. Station 1 was an external control. At each treatment station (2, 3), mesocosms were constructed to enclose a 3.5 m$\sp3$ water column. Each treatment station included an internal control site also. Each in situ mesocosm, except the controls, was successively dosed at a mesocosm-specific dose (1:100; 1:1,000; or 1:10,000 v/v) with 4 field collected drilling muds (spud, nondispersed, lightly-treated, and heavily-treated lignosulfonate) in sequential order over 1.5 months. Twenty-four hours after each dose, water exchange was allowed until the next treatment. Station 3 was destroyed by a winter storm. After the last treatment, the enclosures were removed and the remaining sites monitored for 6 months. One additional site was similarly dosed (1:100 v/v) with clean dredged sediment from Christmas Bay for comparison between dredged sediments and drilling muds.^ Results of the analysis of the water samples and field measurements showed that water quality was impacted during the discharges, primarily at the highest dose (1:100 v/v), but that elevated levels of C, Cr (T,F), Cr$\sp{+3}$ (T, F), N, Pb, and Zn returned to ambient levels before the end of the 24 hour exposure period or immediately after water exchange was allowed (Al, Ba(T), Chlorophyll ABC, SS, %T). Barium, from the barite, was used as a geochemical tracer in the sediments to confirm estimated doses by mass balance calculations. Barium reached a maximum of 166x background levels at the high dose mesocosm. Barium levels returned to ambient or only slightly elevated levels at the end of the 6 month monitoring period due to sediment deposition, resuspension, and bioturbation. QA/QC results using blind samples consisting of lab standards and spiked samples for both water and sediment matrices were within acceptable coefficients of variation.^ In order to avoid impacts on water quality and sediment geochemistry in a shallow estuarine ecosystem, this study concluded that a minimal dilution of 1:1,000 (v/v) would be required in addition to existing regulatory constraints. ^

Maximal exercise hemodynamics and risk of mortality in apparently healthy men and women

The use of exercise electrocardiography (ECG) to detect latent coronary heart disease (CHD) is discouraged in apparently healthy populations because of low sensitivity. These recommendations however, are based on the efficacy of evaluation of ischemia (ST segment changes) with little regard for other measures of cardiac function that are available during exertion. The purpose of this investigation was to determine the association of maximal exercise hemodynamic responses with risk of mortality due to all-causes, cardiovascular disease (CVD), and coronary heart disease (CHD) in apparently healthy individuals. Study participants were 20,387 men (mean age = 42.2 years) and 6,234 women (mean age = 41.9 years) patients of a preventive medicine center in Dallas, TX examined between 1971 and 1989. During an average of 8.1 years of follow-up, there were 348 deaths in men and 66 deaths in women. In men, age-adjusted all-cause death rates (per 10,000 person years) across quartiles of maximal systolic blood pressure (SBP) (low to high) were: 18.2, 16.2, 23.8, and 24.6 (p for trend $<$0.001). Corresponding rates for maximal heart rate were: 28.9, 15.9, 18.4, and 15.1 (p trend $<$0.001). After adjustment for confounding variables including age, resting systolic pressure, serum cholesterol and glucose, body mass index, smoking status, physical fitness and family history of CVD, risks (and 95% confidence interval (CI)) of all-cause mortality for quartiles of maximal SBP, relative to the lowest quartile, were: 0.96 (0.70-1.33), 1.36 (1.01-1.85), and 1.37 (0.98-1.92) for quartiles 2-4 respectively. Similar risks for maximal heart rate were: 0.61 (0.44-0.85), 0.69 (0.51-0.93), and 0.60 (0.41-0.87). No associations were noted between maximal exercise rate-pressure product mortality. Similar results were seen for risk of CVD and CHD death. In women, similar trends in age-adjusted all-cause and CVD death rates across maximal SBP and heart rate categories were observed. Sensitivity of the exercise test in predicting mortality was enhanced when ECG results were evaluated together with maximal exercise SBP or heart rate with a concomitant decrease in specificity. Positive predictive values were not improved. The efficacy of the exercise test in predicting mortality in apparently healthy men and women was not enhanced by using maximal exercise hemodynamic responses. These results suggest that an exaggerated systolic blood pressure or an attenuated heart rate response to maximal exercise are risk factors for mortality in apparently healthy individuals. ^

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