Cost and power analysis of linkage disequilibrium methods for locating genes influencing quantitative traits in humans

Linkage disequilibrium methods can be used to find genes influencing quantitative trait variation in humans. Linkage disequilibrium methods can require smaller sample sizes than linkage equilibrium methods, such as the variance component approach to find loci with a specific effect size. The increase in power is at the expense of requiring more markers to be typed to scan the entire genome. This thesis compares different linkage disequilibrium methods to determine which factors influence the power to detect disequilibrium. The costs of disequilibrium and equilibrium tests were compared to determine whether the savings in phenotyping costs when using disequilibrium methods outweigh the additional genotyping costs.^ Nine linkage disequilibrium tests were examined by simulation. Five tests involve selecting isolated unrelated individuals while four involved the selection of parent child trios (TDT). All nine tests were found to be able to identify disequilibrium with the correct significance level in Hardy-Weinberg populations. Increasing linked genetic variance and trait allele frequency were found to increase the power to detect disequilibrium, while increasing the number of generations and distance between marker and trait loci decreased the power to detect disequilibrium. Discordant sampling was used for several of the tests. It was found that the more stringent the sampling, the greater the power to detect disequilibrium in a sample of given size. The power to detect disequilibrium was not affected by the presence of polygenic effects.^ When the trait locus had more than two trait alleles, the power of the tests maximized to less than one. For the simulation methods used here, when there were more than two-trait alleles there was a probability equal to 1-heterozygosity of the marker locus that both trait alleles were in disequilibrium with the same marker allele, resulting in the marker being uninformative for disequilibrium.^ The five tests using isolated unrelated individuals were found to have excess error rates when there was disequilibrium due to population admixture. Increased error rates also resulted from increased unlinked major gene effects, discordant trait allele frequency, and increased disequilibrium. Polygenic effects did not affect the error rates. The TDT, Transmission Disequilibrium Test, based tests were not liable to any increase in error rates.^ For all sample ascertainment costs, for recent mutations ($<$100 generations) linkage disequilibrium tests were less expensive than the variance component test to carry out. Candidate gene scans saved even more money. The use of recently admixed populations also decreased the cost of performing a linkage disequilibrium test. ^

Effect of nonnormal random components on level and power in group-randomized trials

The use of group-randomized trials is particularly widespread in the evaluation of health care, educational, and screening strategies. Group-randomized trials represent a subset of a larger class of designs often labeled nested, hierarchical, or multilevel and are characterized by the randomization of intact social units or groups, rather than individuals. The application of random effects models to group-randomized trials requires the specification of fixed and random components of the model. The underlying assumption is usually that these random components are normally distributed. This research is intended to determine if the Type I error rate and power are affected when the assumption of normality for the random component representing the group effect is violated. ^ In this study, simulated data are used to examine the Type I error rate, power, bias and mean squared error of the estimates of the fixed effect and the observed intraclass correlation coefficient (ICC) when the random component representing the group effect possess distributions with non-normal characteristics, such as heavy tails or severe skewness. The simulated data are generated with various characteristics (e.g. number of schools per condition, number of students per school, and several within school ICCs) observed in most small, school-based, group-randomized trials. The analysis is carried out using SAS PROC MIXED, Version 6.12, with random effects specified in a random statement and restricted maximum likelihood (REML) estimation specified. The results from the non-normally distributed data are compared to the results obtained from the analysis of data with similar design characteristics but normally distributed random effects. ^ The results suggest that the violation of the normality assumption for the group component by a skewed or heavy-tailed distribution does not appear to influence the estimation of the fixed effect, Type I error, and power. Negative biases were detected when estimating the sample ICC and dramatically increased in magnitude as the true ICC increased. These biases were not as pronounced when the true ICC was within the range observed in most group-randomized trials (i.e. 0.00 to 0.05). The normally distributed group effect also resulted in bias ICC estimates when the true ICC was greater than 0.05. However, this may be a result of higher correlation within the data. ^

AN ESTIMATION OF IN UTERO DRUG EXPOSURE BY DETERMINANTS OF FREQUENCY AND DISTRIBUTION OF MATERNAL DRUG CONSUMPTION

Data derived from 1,194 gravidas presenting at the observation unit of a city/county hospital between October 11, 1979 through December 7, 1979 were evaluated with respect to the proportion ingesting drugs during pregnancy. The mean age of the mother at the time of the interview was 22.0 years; 43.0 percent were Black; 34.0 percent Latin-American, 21.0 percent White and 2.0 percent other; mean gravida was 2.5 pregnancies; mean parity was 1.0; and mean number of previous abortions was 0.34. Completed interview data was available for 1,119 gravida, corresponding urinalyses for 997 subjects. Ninety and one-tenth percent (90.1 percent) of the subjects reported ingestion of one or more drug preparation(s) (prescription, OTC, or substances used for recreational purposes) during pregnancy with a range of 0 to 11 substances and a mean of 2.7. Dietary supplements (vitamins and minerals) were most frequently reported followed by non-narcotic analgesics. Seventy-six and one tenth percent (76.1 percent) of the population reported consumption of prescription medication, 42.5 percent reported consumption of over-the-counter medications, 45.7 percent reported consumption of a substance for recreational purposes and 4.3 percent reported illicit consumption of a substance. For selected substances, no measurable difference was found between obtaining the information from the interview method or from a urinalysis assay. ^

Evidence for the role of agonist binding frequency in receptor -mediated activation of adenylate cyclase

$\beta$-adrenergic receptor-mediated activation of adenylate cyclase exhibits an agonist-specific separation between the dose/response curve (characterized by the EC$\sb{50}$) and the dose/binding curve (characterized by the K$\sb{\rm d}$). Cyclase activity can be near-maximal when receptor occupancy is quite low (EC$\sb{50}$ $\ll$ K$\sb{\rm d}$). This separation between the binding and response curves can be explained by the assumption that the rate of cyclase activation is proportional to the concentration of agonist-bound receptors, since the receptor is mobile and can activate more than one cyclase (the Collision Coupling Model of Tolkovsky and Levitzki). Here it is established that agonist binding frequency plays an additional role in adenylate cyclase activation in S49 murine lymphoma cells. Using epinephrine (EC$\sb{50}$ = 10 nM, K$\sb{\rm d}$ = 2 $\mu$M), the rate of cyclase activation decreased by 80% when a small (1.5%) receptor occupancy was restricted (by addition of the antagonist propranolol) to a small number (1.5%) of receptors rather than being proportionally distributed among the cell's entire population of receptors. Thus adenylate cyclase activity is not proportional to receptor occupancy in all circumstances. Collisions between receptor and cyclase pairs apparently occur a number of times in rapid sequence (an encounter); the high binding frequency of epinephrine ensures that discontiguous regions of the cell surface experience some period of agonist-bound receptor activity per small unit time minimizing "wasted" collisions between activated cyclase and bound receptor within an encounter. A contribution of agonist binding frequency to activation is thus possible when: (1) the mean lifetime of the agonist-receptor complex is shorter than the mean encounter time, and (2) the absolute efficiency (intrinsic ability to promote cyclase activation per collision) of the agonist-receptor complex is high. These conclusions are supported by experiments using agonists of different efficiencies and binding frequencies. These results are formalized in the Encounter Coupling Model of adenylate cyclase activation, which takes into explicit account the agonist binding frequency, agonist affinity for the $\beta$-adrenergic receptor, agonist efficiency, encounter frequency and the encounter time between receptor and cyclase. ^