An ordinal logistic regression model with misclassification of the outcome variable and categorical covariate

Ordinal outcomes are frequently employed in diagnosis and clinical trials. Clinical trials of Alzheimer's disease (AD) treatments are a case in point using the status of mild, moderate or severe disease as outcome measures. As in many other outcome oriented studies, the disease status may be misclassified. This study estimates the extent of misclassification in an ordinal outcome such as disease status. Also, this study estimates the extent of misclassification of a predictor variable such as genotype status. An ordinal logistic regression model is commonly used to model the relationship between disease status, the effect of treatment, and other predictive factors. A simulation study was done. First, data based on a set of hypothetical parameters and hypothetical rates of misclassification was created. Next, the maximum likelihood method was employed to generate likelihood equations accounting for misclassification. The Nelder-Mead Simplex method was used to solve for the misclassification and model parameters. Finally, this method was applied to an AD dataset to detect the amount of misclassification present. The estimates of the ordinal regression model parameters were close to the hypothetical parameters. β1 was hypothesized at 0.50 and the mean estimate was 0.488, β2 was hypothesized at 0.04 and the mean of the estimates was 0.04. Although the estimates for the rates of misclassification of X1 were not as close as β1 and β2, they validate this method. X 1 0-1 misclassification was hypothesized as 2.98% and the mean of the simulated estimates was 1.54% and, in the best case, the misclassification of k from high to medium was hypothesized at 4.87% and had a sample mean of 3.62%. In the AD dataset, the estimate for the odds ratio of X 1 of having both copies of the APOE 4 allele changed from an estimate of 1.377 to an estimate 1.418, demonstrating that the estimates of the odds ratio changed when the analysis includes adjustment for misclassification. ^

Multivariate methods for correlated effects sizes

Current statistical methods for estimation of parametric effect sizes from a series of experiments are generally restricted to univariate comparisons of standardized mean differences between two treatments. Multivariate methods are presented for the case in which effect size is a vector of standardized multivariate mean differences and the number of treatment groups is two or more. The proposed methods employ a vector of independent sample means for each response variable that leads to a covariance structure which depends only on correlations among the $p$ responses on each subject. Using weighted least squares theory and the assumption that the observations are from normally distributed populations, multivariate hypotheses analogous to common hypotheses used for testing effect sizes were formulated and tested for treatment effects which are correlated through a common control group, through multiple response variables observed on each subject, or both conditions.^ The asymptotic multivariate distribution for correlated effect sizes is obtained by extending univariate methods for estimating effect sizes which are correlated through common control groups. The joint distribution of vectors of effect sizes (from $p$ responses on each subject) from one treatment and one control group and from several treatment groups sharing a common control group are derived. Methods are given for estimation of linear combinations of effect sizes when certain homogeneity conditions are met, and for estimation of vectors of effect sizes and confidence intervals from $p$ responses on each subject. Computational illustrations are provided using data from studies of effects of electric field exposure on small laboratory animals. ^

Assessment of occupational exposures to organic chemicals in a plastics manufacturing facility

The objective of this dissertation was to design and implement strategies for assessment of exposures to organic chemicals used in the production of a styrene-butadiene polymer at the Texas Plastics Company (TPC). Linear statistical retrospective exposure models, univariate and multivariate, were developed based on the validation of historical industrial hygiene monitoring data collected by industrial hygienists at TPC, and additional current industrial hygiene monitoring data collected for the purposes of this study. The current monitoring data served several purposes. First, it provided information on current exposure data, in the form of unbiased estimates of mean exposure to organic chemicals for each job title included. Second, it provided information on homogeneity of exposure within each job title, through the use of a carefully designed sampling scheme which addressed variability of exposure both between and within job titles. Third, it permitted the investigation of how well current exposure data can serve as an evaluation tool for retrospective exposure estimation. Finally, this dissertation investigated the simultaneous evaluation of exposure to several chemicals, as well as the use of values below detection limits in a multivariate linear statistical model of exposures. ^

The exchange of carbon, nitrogen, and phosphorus in dwarf and fringe mangroves of the oligotrophic southern Everglades

Water management has altered both the natural timing and volume of freshwater delivered to Everglades National Park. This is especially true for Taylor Slough and the C-111 basin, as hypersaline events in Florida Bay have been linked to reduced freshwater flow in this area. In light of recent efforts to restore historical flows to the eastern Everglades, an understanding of the impact of this hydrologic shift is needed in order to predict the trajectory of restoration. I conducted a study to assess the importance of season, water chemistry, and hydrologic conditions on the exchange of nutrients in dwarf and fringe mangrove wetlands along Taylor Slough. I also performed mangrove leaf decomposition studies to determine the contribution of biotic and abiotic processes to mass loss, the effect of salinity and season on degradation rates, and the importance of this litter component as a rapid source of nutrients. ^ Dwarf mangrove wetlands consistently imported total nutrients (C, N, and P) and released NO2− +NO3 −, with enhanced release during the dry season. Ammonium flux shifted from uptake to release over the study period. Dissolved phosphate activity was difficult to discern in either wetland, as concentrations were often below detection limits. Fluxes of dissolved inorganic nitrogen in the fringe wetland were positively related to DIN concentrations. The opposite was found for total nitrogen in the fringe wetland. A dynamic budget revealed a net annual export of TN to Florida Bay that was highest during the wet season. Simulated increases and decreases in freshwater flow yielded reduced exports of TN to Florida Bay as a result of changes in subsystem and water flux characteristics. Finally, abiotic processes yielded substantial nutrient and mass losses from senesced leaves with little influence of salinity. Dwarf mangrove leaf litter appeared to be a considerable source of nutrients to the water column of this highly oligotrophic wetland. To summarize, nutrient dynamics at the subsystem level were sensitive to short-term changes in hydrologic and seasonal conditions. These findings suggest that increased freshwater flow has the potential to lead to long-term, system-level changes that may reach as far as eastern Florida Bay. ^

Ecosystem responses to hydrologic change and mechanisms of nitrogen sequestration in seasonally flooded tree islands of the southern Everglades

In the Florida Everglades, tree islands are conspicuous heterogeneous elements in a complex wetland landscape. I investigated the effects of increased freshwater flow in southern Everglades seasonally flooded tree islands, and characterized biogeochemical interactions among tree islands and the marsh landscape matrix, specifically examining hydrologic flows of nitrogen (N), and landscape N sequestration capacity. I utilized ecological trajectories of key ecosystem variables to differentiate effects of increased sheetflow and hydroperiod. I utilized stable isotope analyses and nutrient content of tree island ecosystem components to test the hypothesis that key processes in tree island nitrogen cycling would favor ecosystem N sequestration. I combined estimates of tree island ecosystem N standing stocks and fluxes, soil and litter N transformation rates, and hydrologic inputs of N to quantify the net sequestration of N by a seasonally flooded tree island. ^ Results show that increased freshwater flow to seasonally flooded tree islands promoted ecosystem oligotrophy, whereas reduced flows allowed some plant species to cycle P less efficiently. As oligotrophy is a defining parameter of Everglades wetlands, and likely promotes belowground production and peat development, reintroducing freshwater flow from an upstream canal had a favorable effect on ecosystem dynamics of tree islands in the study area. Important factors influencing the stable isotopic composition of nitrogen and carbon were: (1) a contribution to soil N by soil invertebrates, animal excrement, and microbes, (2) a possible NO3 source from an upstream canal and an "open" ecosystem N cycle, and (3) greater availability of phosphorus in tree islands relative to the marsh landscape, suggesting that tree island N cycling favors N sequestration. Hydrologic sources of N were dominated by surface water loads of NO3- and NH 4+, and an important soil N transformation promoting the net loss of surface water DIN was nitrate immobilization associated with soils and surficial leaf litter. The net inorganic N sequestration capacity of a seasonally flooded tree island was 50 g yr-1 m -2. Thus, tree islands likely have an important function in landscape sequestration of inorganic N, and may reduce significant anthropogenic N loads to downstream coastal systems. ^

Genetic predisposition to prostate cancer

Prostate cancer (PC) is a significant economic and health burden in the U.S. and Europe but its causes are largely unknown. The most significant risk factors (after gender) are age and family history of the disease. A gene with high penetrance but low frequency on chromosome 1q, HPC 1, has been suggested to cause a proportion of the familial aggregation of PC but other more common genes, conferring less risk, are also thought to contribute to disease predisposition. We have pursued a strategy to study both types of genetic risk in PC. To identify high penetrance genes, affected men from thirteen families have been genotyped for genetic linkage analysis at six microsatellite markers spanning 45 cM of 1q24-25. Both LOD score and non-parametric statistics provide no significant support for HPC1 in this genomic region, although 3 of the families did combine to produce a LOD score of 0.9. These families will be included in a genome wide search for other PC predisposition genes as part of a multinational collaboration.^ For study of common genetic factors in PC development, leukocyte DNA samples from an unselected series of 55 patients and 67 controls have been examined for genetic differences in two other candidate genes, the androgen receptor gene, hAR, at Xq11-12, and the vitamin D receptor gene, hVDR, at 12q12-14. hAR was typed for two trinucleotide repeat length polymorphisms, (CAG)$\rm\sb{n}$ and (GGC)$\rm\sb{n},$ encoding polyglutamine and polyglycine tracts, respectively, which have been implicated in PC susceptibility. These data, combined with similarly processed patients and controls from the U.K. show no consistent association of allele length with PC risk. A novel finding, however, has been a significant association between the number of GGC repeats and the length of time between diagnosis and relapse in stage T1-T4 Caucasian patients irrespective of therapy and age of the patient. Of 49 patients who relapsed out of 108 entering the study, those with 16 or fewer GGC repeats had an average relapse-free-period of 101 (+/$-$7.7) months while for those with more than 16 repeats the period averaged 48 (+/$-$2.9) months, a difference of 2.1 fold or 4.4 years.^ The second gene, hVDR, was genotyped at two polymorphisms, a synonymous C/T substitution in exon 9 identified by differential TaqI enzymatic digestion and a variable length polyA tract in the 3$\sp\prime$ UTR. Although these polymorphisms are in strong linkage disequilibrium only the polyA region showed a possible association with PC risk. Men homozygous for alleles with fewer than 18 A's had an increased risk (OR = 3.0, p = 0.0578) compared to controls. This result is opposite to the findings of others and may either indicate off-setting random errors which together balance out to no significant overall effect or reflect more complex genetic and/or environmental associations.^ Overall, this research suggests that single gene familial predisposition may be less prominent in PC than in other cancers and that the characteristics of PC pathology may be useful in identifying the effects of common genetic factors. ^

Studying phenotypic evolution using multivariate quantitative genetics

Quantitative genetics provides a powerful framework for studying phenotypic evolution and the evolution of adaptive genetic variation. Central to the approach is G, the matrix of additive genetic variances and covariances. G summarizes the genetic basis of the traits and can be used to predict the phenotypic response to multivariate selection or to drift. Recent analytical and computational advances have improved both the power and the accessibility of the necessary multivariate statistics. It is now possible to study the relationships between G and other evolutionary parameters, such as those describing the mutational input, the shape and orientation of the adaptive landscape, and the phenotypic divergence among populations. At the same time, we are moving towards a greater understanding of how the genetic variation summarized by G evolves. Computer simulations of the evolution of G, innovations in matrix comparison methods, and rapid development of powerful molecular genetic tools have all opened the way for dissecting the interaction between allelic variation and evolutionary process. Here I discuss some current uses of G, problems with the application of these approaches, and identify avenues for future research.

Comparison of the Salmonella/Microsome Microsuspension Assay with the new Microplate Fluctuation Protocol for Testing the Mutagenicity of Environmental Samples

The objective of this study was to compare the responses of the Salmonella/microsome microsuspension assay with the new microplate fluctuation protocol (MPF) for the evaluation of the mutagenic activity of environmental samples. Organic extracts of total particulate atmospheric air samples, surface waters, and effluents were tested in dose-response experiments. The assays were performed with strain TA98 in the absence and presence of S9 mix. Both protocols produced similar results, despite the fact that the maximum score of the MPF is limited to 48 wells, whereas in the regular plate assay it is possible to count up to 1,500 colonies using an automatic counter. Similar sensitivities based on the lowest dose that resulted in a positive response were obtained for both assays. The MPF procedure is less laborious (e.g., all-liquid format, use of multi-channel pipettors) and allows for automation of the pipetting and dispensing steps, thus, reducing time of the analysis which is particularly important in environmental quality monitoring programs or in effect-directed analysis. The results show that the MPF procedure is a promising tool to test environmental samples for mutagenic activity. Environ. Mol. Mutagen. 51:31-38, 2010. (C) 2009 Wiley-Liss, Inc.

Cryptic vicariance in Gulf of California fishes parallels vicariant patterns found in Baja, California mammals and reptiles

Comparisons across multiple taxa can often clarify the histories of biogeographic regions. In particular, historic barriers to movement should affect multiple species and, thus, result in a pattern of concordant intraspecific genetic divisions among species. A striking example of such comparative phylogeography is the recent observation that populations of many small mammals and reptiles living on the Baja, California peninsula have a large genetic break between northern and southern peninsular populations. In the present study, I demonstrate that five species of near-shore fishes living on the Baja coastline of the Gulf of California share this genetic pattern. The simplest explanation for this concordant genetic division within both terrestrial and marine vertebrates is that the Baja peninsula was fragmented by a Plio-Pleistocene marine seaway and that this seaway posed a substantial barrier to movement for near-shore fishes. The genetic divisions within Gulf of California fishes also coincide with recognized biogeographic regions based on fish community composition and several environmental factors. It is likely that adaptation to regional environments and present-day oceanographic circulation limits gene exchange between biogeographic regions and helps maintain evidence of past vicariance.

Genetic covariance between indices of body condition and immunocompetence in a passerine bird

Background: Condition-dependence is a ubiquitous feature of animal life histories and has important implications for both natural and sexual selection. Mate choice, for instance, is typically based on condition-dependent signals. Theory predicts that one reason why condition-dependent signals may be special is that they allow females to scan for genes that confer high parasite resistance. Such explanations require a genetic link between immunocompetence and body condition, but existing evidence is limited to phenotypic associations. It remains unknown, therefore, whether females selecting males with good body condition simply obtain a healthy mate, or if they acquire genes for their offspring that confer high immunocompetence. Results: Here we use a cross-foster experimental design to partition the phenotypic covariance in indices of body condition and immunocompetence into genetic, maternal and environmental effects in a passerine bird, the zebra finch Taeniopygia guttata. We show that there is significant positive additive genetic covariance between an index of body condition and an index of cell-mediated immune response. In this case, genetic variance in the index of immune response explained 56% of the additive genetic variance in the index of body condition. Conclusion: Our results suggest that, in the context of sexual selection, females that assess males on the basis of condition-dependent signals may gain genes that confer high immunocompetence for their offspring. More generally, a genetic correlation between indices of body condition and imuunocompetence supports the hypothesis that parasite resistance may be an important target of natural selection. Additional work is now required to test whether genetic covariance exists among other aspects of both condition and immunocompetence.

Natural selection and quantitative genetics of life-history traits in Western women: A twin study

Whether contemporary human populations are still evolving as a result of natural selection has been hotly debated. For natural selection to cause evolutionary change in a trait, variation in the trait must be correlated with fitness and be genetically heritable and there must be no genetic constraints to evolution. These conditions have rarely been tested in human populations. In this study, data from a large twin cohort were used to assess whether selection Will cause a change among women in contemporary Western population for three life-history traits: age at menarche, age at first reproduction, and age at menopause. We control for temporal variation in fecundity (the baby boom phenomenon) and differences between women in educational background and religious affiliation. University-educated women have 35% lower fitness than those with less than seven years education, and Roman Catholic women have about 20% higher fitness than those of other religions. Although these differences were significant, education and religion only accounted for 2% and 1% of variance in fitness, respectively. Using structural equation modeling, we reveal significant genetic influences for all three life-history traits, with heritability estimates of 0.50, 0.23, and 0.45, respectively. However, strong genetic covariation with reproductive fitness could only be demonstrated for age at first reproduction, with much weaker covariation for age at menopause and no significant covariation for age at menarche. Selection may, therefore, lead to the evolution of earlier age at first reproduction in this population. We also estimate substantial heritable variation in fitness itself, with approximately 39% of the variance attributable to additive genetic effects, the remainder consisting of unique environmental effects and small effects from education and religion. We discuss mechanisms that could be maintaining such a high heritability for fitness. Most likely is that selection is now acting on different traits from which it did in pre-industrial human populations.

Resolution of natural groups using iterative assignment tests: an example from two species of Australian native rats (Rattus)

Sympatric individuals of Rattus fuscipes and Rattus leucopus, two Australian native rats from the tropical wet forests of north Queensland, are difficult to distinguish morphologically and are often confused in the field. When we started a study on fine-scale movements of these species, using microsatellite markers, we found that the species as identified in the field did not form coherent genetic groups. In this study, we examined the potential of an iterative process of genetic assignment to separate specimens from distinct (e.g. species, populations) natural groups. Five loci with extensive overlap in allele distributions between species were used for the iterative process. Samples were randomly distributed into two starting groups of equal size and then subjected to the test. At each iteration, misassigned samples switched groups, and the output groups from a given round of assignment formed the input groups for the next round. All samples were assigned correctly on the 10th iteration, in which two genetic groups were clearly separated. Mitochondrial DNA sequences were obtained from samples from each genetic group identified by assignment, together with those of museum voucher specimens, to assess which species corresponded to which genetic group. The iterative procedure was also used to resolve groups within species, adequately separating the genetically identified R. leucopus from our two sampling sites. These results show that the iterative assignment process can correctly differentiate samples into their appropriate natural groups when diagnostic genetic markers are not available, which allowed us to resolve accurately the two R. leucopus and R. fuscipes species. Our approach provides an analytical tool that may be applicable to a broad variety of situations where genetic groups need to be resolved.

Genetics of melanoma predisposition

Predisposition to melanoma is genetically heterogeneous. Two high penetrance susceptibility genes, CDKN2A and CDK4, have so far been identified and mapping is ongoing to localize and identify others. With the advent of a catalogue of millions of potential DNA polymorphisms, attention is now also being focused on identification of genes that confer a more modest contribution to melanoma risk, such as those encoding proteins involved in pigmentation, DNA repair, cell growth and differentiation or detoxification of metabolites. One such pigmentation gene, MC1R, has not only been found to be a low penetrance melanoma gene but has also been shown to act as a genetic modifier of melanoma risk in individuals carrying CDKN2A mutations. Most recently, an environmental agent, ultraviolet radiation, has also been established as a modifier of melanoma risk in CDKN2A mutation carriers. Hence, melanoma is turning out to be an excellent paradigm for studying gene-gene and gene-environment interactions.

Zebrafish embryos as in vivo model for toxicity evaluation: screening of nanoparticle formulations DODAB:MO and DODAC:MO

Dissertação de mestrado em Molecular Genetics

Assessing the quality of alternative energy sources: Energy Return On the Investment (EROI), the metabolic pattern of societies and energy statistics

This paper presents an initial challenge to tackle the every so "tricky" points encountered when dealing with energy accounting, and thereafter illustrates how such a system of accounting can be used when assessing for the metabolic changes in societies. The paper is divided in four main sections. The first three, present a general discussion on the main issues encountered when conducting energy analyses. The last section, subsequently, combines this heuristic approach to the actual formalization of it, in quantitative terms, for the analysis of possible energy scenarios. Section one covers the broader issue of how to account for the relevant categories used when accounting for Joules of energy; emphasizing on the clear distinction between Primary Energy Sources (PES) (which are the physical exploited entities that are used to derive useable energy forms (energy carriers)) and Energy Carriers (EC) (the actual useful energy that is transmitted for the appropriate end uses within a society). Section two sheds light on the concept of Energy Return on Investment (EROI). Here, it is emphasized that, there must already be a certain amount of energy carriers available to be able to extract/exploit Primary Energy Sources to thereafter generate a net supply of energy carriers. It is pointed out that this current trend of intense energy supply has only been possible to the great use and dependence on fossil energy. Section three follows up on the discussion of EROI, indicating that a single numeric indicator such as an output/input ratio is not sufficient in assessing for the performance of energetic systems. Rather an integrated approach that incorporates (i) how big the net supply of Joules of EC can be, given an amount of extracted PES (the external constraints); (ii) how much EC needs to be invested to extract an amount of PES; and (iii) the power level that it takes for both processes to succeed, is underlined. Section four, ultimately, puts the theoretical concepts at play, assessing for how the metabolic performances of societies can be accounted for within this analytical framework.