In Utero Exposure to Antiretroviral Drugs: Effect on Birth Weight and Growth Among HIV-exposed Uninfected Children in Brazil.

BACKGROUND There are concerns about the effects of in utero exposure to antiretroviral drugs (ARVs) on the development of HIV-exposed but uninfected (HEU) children. The aim of this study was to evaluate whether in utero exposure to ARVs is associated with lower birth weight/height and reduced growth during the first 2 years of life. METHODS This cohort study was conducted among HEU infants born between 1996 and 2010 in Tertiary children's hospital in Rio de Janeiro, Brazil. Weight was measured by mechanical scale, and height was measured by measuring board. Z-scores for weight-for-age (WAZ), length-for-age (LAZ) and weight-for-length were calculated. We modeled trajectories by mixed-effects models and adjusted for mother's age, CD4 cell count, viral load, year of birth and family income. RESULTS A total of 588 HEU infants were included of whom 155 (26%) were not exposed to ARVs, 114 (19%) were exposed early (first trimester) and 319 (54%) later. WAZ were lower among infants exposed early compared with infants exposed later: adjusted differences were -0.52 (95% confidence interval [CI]: -0.99 to -0.04, P = 0.02) at birth and -0.22 (95% CI: -0.47 to 0.04, P = 0.10) during follow-up. LAZ were lower during follow-up: -0.35 (95% CI: -0.63 to -0.08, P = 0.01). There were no differences in weight-for-length scores. Z-scores of infants exposed late during pregnancy were similar to unexposed infants. CONCLUSIONS In HEU children, early exposure to ARVs was associated with lower WAZ at birth and lower LAZ up to 2 years of life. Growth of HEU children needs to be monitored closely.

Assessing homogeneity in workplace exposure under stationary and non-stationary scenarios

In most epidemiological studies, historical monitoring data are scant and must be pooled to identify occupational groups with homogeneous exposures. Homogeneity of exposure is generally assessed in a group of workers who share a common job title or work in a common area. While published results suggest that the degree of homogeneity varies widely across job groups, less is known whether such variation differs across industrial sectors, classes of contaminants, or in the methods used to group workers. Relying upon a compilation of results presented in the literature, patterns of homogeneity among nearly 500 occupational groups of workers were evaluated on the basis of type of industry and agent. Additionally, effects of the characteristics of the sampling strategy on estimated indicators of homogeneity of exposure were assessed. ^ Exposure profiles for occupational groups of workers have typically been assessed under the assumption of stationarity, i.e., the mean exposure level and variance of the distribution that describes the underlying population of exposures are constant over time. Yet, the literature has shown that occupational exposures have declined in the last decades. This renders traditional methods for the description of exposure profiles inadequate. Thus, work was needed to develop appropriate methods to assess homogeneity for groups of workers whose exposures have changed over time. A study was carried out applying mixed effects models with a term for temporal trend to appropriately describe exposure profiles of groups of workers in the nickel-producing industry over a 20-year period. Using a sub-set of groups of nickel-exposed workers, another study was conducted to develop and apply a framework to evaluate the assumption of stationarity of the variances in the presence of systematic changes in exposure levels over time. ^

Bayesian joint modeling of longitudinal and survival data

The joint modeling of longitudinal and survival data is a new approach to many applications such as HIV, cancer vaccine trials and quality of life studies. There are recent developments of the methodologies with respect to each of the components of the joint model as well as statistical processes that link them together. Among these, second order polynomial random effect models and linear mixed effects models are the most commonly used for the longitudinal trajectory function. In this study, we first relax the parametric constraints for polynomial random effect models by using Dirichlet process priors, then three longitudinal markers rather than only one marker are considered in one joint model. Second, we use a linear mixed effect model for the longitudinal process in a joint model analyzing the three markers. In this research these methods were applied to the Primary Biliary Cirrhosis sequential data, which were collected from a clinical trial of primary biliary cirrhosis (PBC) of the liver. This trial was conducted between 1974 and 1984 at the Mayo Clinic. The effects of three longitudinal markers (1) Total Serum Bilirubin, (2) Serum Albumin and (3) Serum Glutamic-Oxaloacetic transaminase (SGOT) on patients' survival were investigated. Proportion of treatment effect will also be studied using the proposed joint modeling approaches. ^ Based on the results, we conclude that the proposed modeling approaches yield better fit to the data and give less biased parameter estimates for these trajectory functions than previous methods. Model fit is also improved after considering three longitudinal markers instead of one marker only. The results from analysis of proportion of treatment effects from these joint models indicate same conclusion as that from the final model of Fleming and Harrington (1991), which is Bilirubin and Albumin together has stronger impact in predicting patients' survival and as a surrogate endpoints for treatment. ^

Understanding disparities in emergency department visits for asthma

Asthma is the most common chronic disorder in childhood, affecting an estimated 6.2 million children under 18 years (1). The purpose of this study was to look at individual- and community-level characteristics simultaneously to examine and explain the factors that contribute to the use of emergency department services by children 18 years old or less and to determine if there was an association between air quality and ED visits in the same population, from 2005-2007 in Houston/Harris County. Data were collected from the Houston Safety Net Hospital Emergency Department Use Study and the 2000 US Census. Bivariate and multivariate logistic regression models and mixed effects models were used to analyze data that was collected during the study period.^ There were 704,902 ED visits made by children 18 and younger, who were living in Houston from January 1, 2005 to December 31, 2007. Of those, 19,098 had a primary discharge diagnosis of asthma. Asthma ED visits varied by season, with proportions of ED visits for asthma highest from September-December. African-American children were 2.6 (95% CI, 2.43-2.66) times more likely to have an ED visit for asthma compared to White children. Poverty, single parent headed households, and younger age all a greater likelihood of having gone to the ED for asthma treatment. Compared to Whites living in lightly-monitored pollution areas, African-Americans and Hispanics living in heavily monitored areas were 1.15 (95% CI, 1.04-1.28) times more likely to have an ED visit for asthma.^ Race and poverty seem to account for a large portion of the disparities in ED use found among children. This was true even after accounting for multiple individual- and community-level variables. These results suggest that racial disparities in asthma continue to pose risks for African American children, and they point to the need for additional research into potential explanations and remedies. Programs to reduce inappropriate ED use must be sensitive to an array of complex socioeconomic issues within minority and income populations. ^

The role of cardiorespiratory fitness and weight status on 24-hr ambulatory blood pressure among adolescents

BACKGROUND: This observational research study investigated the association of cardiorespiratory fitness and weight status with repeated measures of 24-hr ambulatory blood pressure (24-hr ABP). Little is known about these associations and few data exist examining the interaction between cardiorespiratory fitness and weight status and the contributions of each on 24-hr ABP in youth. ^ METHODS: This research study used secondary analysis data from the "Adolescent Blood Pressure and Anger: Ethnic Differences" study. This current study sample included 374 African-American, Anglo-American, and Mexican-American adolescents 11-16 years of age. Mixed-effects models were used for testing the relationship between weight status and cardiorespiratory fitness and repeated measures of ambulatory blood pressure over 24 hours (24-hr ABP). Weight status was categorized into "normal weight" (BMI<85th percentile), "overweight" (85th≤BMI<95th), and "obese" (BMI≥95th). Cardiorespiratory fitness, determined by heart rate recovery (HRR), was defined as the difference between heart rate at peak exercise and heart rate at two minutes post-exercise, as measured by a height-adjusted step test and stratified into two groups: low and high fitness, using a median split. Ambulatory blood pressure (ABP) was monitored for a 24-hr period on a school day using the Spacelabs ambulatory monitor (Model 90207). Blood pressure and heart rate were recorded at 30 minute intervals throughout the day of recording and at 60 minute intervals during sleep. ^ RESULTS: No significant associations were found between weight status and mean 24-hr systolic blood pressure (SBP) or mean arterial pressure (MAP). A significant and inverse association between weight status and mean 24-hr diastolic blood pressure (DBP) was revealed. Cardiorespiratory fitness was significantly and inversely associated with mean 24-hr ABP. High fitness adolescents had significantly lower mean 24-hr SPB, DBP, and MAP measurements than low fitness adolescents. Compared to low fitness adolescents, high fitness adolescents had 1.90 mmHg, 1.16 mmHg, and 1.68 mmHg lower mean 24-hr SBP, DBP, and MAP, respectively. Additionally, high fitness appeared to afford protection from higher mean 24-hr SBP and MAP, irrespective of weight status. Among normal weight adolescents, low fitness resulted in higher mean 24-hr SBP and MAP, compared to their fit counterparts. Among adolescents categorized as high fitness, increasing weight status did not appear to result in higher mean 24-hr SBP or MAP. Cardiorespiratory fitness, rather than weight status, appeared to be a more dominant predictor of mean 24-hr SBP and MAP. ^ CONCLUSIONS: To our knowledge, this research is the first study to investigate the independent and combined contributions of cardiorespiratory fitness and weight status on 24-hr ABP, all objectively measured. The results of this study may potentially guide and inform future research. It appears that early cardiovascular disease (CVD) prevention should focus on improving cardiorespiratory fitness levels among all adolescents, particularly those adolescents least fit, regardless of their weight status, while obesity prevention efforts continue.^

Dive depth frequency of 23 adult male southern elephant seals (Mirounga leonina) from Marion Island and King George Island

Access to different environments may lead to inter-population behavioural changes within a species that allow populations to exploit their immediate environments. Elephant seals from Marion Island (MI) and King George Island (KGI) (Isla 25 de Mayo) forage in different oceanic environments and evidently employ different foraging strategies. This study elucidates some of the factors influencing the diving behaviour of male southern elephant seals from these populations tracked between 1999 and 2002. Mixed-effects models were used to determine the influence of bathymetry, population of origin, body length (as a proxy for size) and individual variation on the diving behaviour of adult male elephant seals from the two populations. Males from KGI and MI showed differences in all dive parameters. MI males dived deeper and longer (median: 652.0 m and 34.00 min) than KGI males (median: 359.1 m and 25.50 min). KGI males appeared to forage both benthically and pelagically while MI males in this study rarely reached depths close to the seafloor and appeared to forage pelagically. Model outputs indicate that males from the two populations showed substantial differences in their dive depths, even when foraging in areas of similar water depth. Whereas dive depths were not significantly influenced by the size of the animals, size played a significant role in dive durations, though this was also influenced by the population that elephant seals originated from. This study provides some support for inter-population differences in dive behaviour of male southern elephant seals.

(Table 2) Carbon and nitrogen isotopic ratios in blood sampled from arctic foxes (Vulpes lagopus) on Bylot Island between 2003-2008

Inter-individual variation in diet within generalist animal populations is thought to be a widespread phenomenon but its potential causes are poorly known. Inter-individual variation can be amplified by the availability and use of allochthonous resources, i.e., resources coming from spatially distinct ecosystems. Using a wild population of arctic fox as a study model, we tested hypotheses that could explain variation in both population and individual isotopic niches, used here as proxy for the trophic niche. The arctic fox is an opportunistic forager, dwelling in terrestrial and marine environments characterized by strong spatial (arctic-nesting birds) and temporal (cyclic lemmings) fluctuations in resource abundance. First, we tested the hypothesis that generalist foraging habits, in association with temporal variation in prey accessibility, should induce temporal changes in isotopic niche width and diet. Second, we investigated whether within-population variation in the isotopic niche could be explained by individual characteristics (sex and breeding status) and environmental factors (spatiotemporal variation in prey availability). We addressed these questions using isotopic analysis and Bayesian mixing models in conjunction with linear mixed-effects models. We found that: i) arctic fox populations can simultaneously undergo short-term (i.e., within a few months) reduction in both isotopic niche width and inter-individual variability in isotopic ratios, ii) individual isotopic ratios were higher and more representative of a marine-based diet for non-breeding than breeding foxes early in spring, and iii) lemming population cycles did not appear to directly influence the diet of individual foxes after taking their breeding status into account. However, lemming abundance was correlated to proportion of breeding foxes, and could thus indirectly affect the diet at the population scale.

Water temperature data from southern elephant seal from Marion Island with links to datasets

As the Antarctic Circumpolar Current crosses the South-West Indian Ocean Ridge, it creates an extensive eddy field characterised by high sea level anomaly variability. We investigated the diving behaviour of female southern elephant seals from Marion Island during their post-moult migrations in relation to this eddy field in order to determine its role in the animals' at-sea dispersal. Most seals dived within the region significantly more often than predicted by chance, and these dives were generally shallower and shorter than dives outside the eddy field. Mixed effects models estimated reductions of 44.33 ± 3.00 m (maximum depth) and 6.37 ± 0.10 min (dive duration) as a result of diving within the region, along with low between-seal variability (maximum depth: 5.5 % and dive duration: 8.4 %). U-shaped dives increased in frequency inside the eddy field, whereas W-shaped dives with multiple vertical movements decreased. Results suggest that Marion Island's adult female elephant seals' dives are characterised by lowered cost-of-transport when they encounter the eddy field during the start and end of their post-moult migrations. This might result from changes in buoyancy associated with varying body condition upon leaving and returning to the island. Our results do not suggest that the eddy field is a vital foraging ground for Marion Island's southern elephant seals. However, because seals preferentially travel through this area and likely forage opportunistically while minimising transport costs, we hypothesise that climate-mediated changes in the nature or position of this region may alter the seals' at-sea dispersal patterns.

Prospective longitudinal study of subcortical brain volumes in individuals at high familial risk of mood disorders with or without subsequent onset of depression

Subcortical volumetric brain abnormalities have been observed in mood disorders. However, it is unknown whether these reflect adverse effects predisposing to mood disorders or emerge at illness onset. Magnetic resonance imaging was conducted at baseline and after two years in 111 initially unaffected young adults at increased risk of mood disorders because of a close family history of bipolar disorder and 93 healthy controls (HC). During the follow-up, 20 high-risk subjects developed major depressive disorder (HR-MDD), with the others remaining well (HR-well). Volumes of the lateral ventricles, caudate, putamen, pallidum, thalamus, hippocampus and amygdala were extracted for each hemisphere. Using linear mixed-effects models, differences and longitudinal changes in subcortical volumes were investigated between groups (HC, HR-MDD, HR-well). There were no significant differences for any subcortical volume between groups controlling for multiple testing. Additionally, no significant differences emerged between groups over time. Our results indicate that volumetric subcortical brain abnormalities of these regions using the current method appear not to form familial trait markers for vulnerability to mood disorders in close relatives of bipolar disorder patients over the two-year time period studied. Moreover, they do not appear to reduce in response to illness onset at least for the time period studied.

Heat stress, heat strain, and productivity in Washington State tree fruit harvesters

Thesis (Master's)--University of Washington, 2016-06

Environmental Enteropathy and Fibroblast Growth Factor 21 are associated with early childhood growth in Bangladesh

Thesis (Ph.D.)--University of Washington, 2016-06

Toward better outcomes with tacrolimus therapy: Population pharmacokinetics and individualized dosage prediction in adult liver transplantation

Patient outcomes in transplantation would improve if dosing of immunosuppressive agents was individualized. The aim of this study is to develop a population pharmacokinetic model of tacrolimus in adult liver transplant recipients and test this model in individualizing therapy. Population analysis was performed on data from 68 patients. Estimates were sought for apparent clearance (CL/F) and apparent volume of distribution (V/F) using the nonlinear mixed effects model program (NONMEM). Factors screened for influence on these parameters were weight, age, sex, transplant type, biliary reconstructive procedure, postoperative day, days of therapy, liver function test results, creatinine clearance, hematocrit, corticosteroid dose, and interacting drugs. The predictive performance of the developed model was evaluated through Bayesian forecasting in an independent cohort of 36 patients. No linear correlation existed between tacrolimus dosage and trough concentration (r(2) = 0.005). Mean individual Bayesian estimates for CL/F and V/F were 26.5 8.2 (SD) L/hr and 399 +/- 185 L, respectively. CL/F was greater in patients with normal liver function. V/F increased with patient weight. CL/F decreased with increasing hematocrit. Based on the derived model, a 70-kg patient with an aspartate aminotransferase (AST) level less than 70 U/L would require a tacrolimus dose of 4.7 mg twice daily to achieve a steady-state trough concentration of 10 ng/mL. A 50-kg patient with an AST level greater than 70 U/L would require a dose of 2.6 mg. Marked interindividual variability (43% to 93%) and residual random error (3.3 ng/mL) were observed. Predictions made using the final model were reasonably nonbiased (0.56 ng/mL), but imprecise (4.8 ng/mL). Pharmacokinetic information obtained will assist in tacrolimus dosing; however, further investigation into reasons for the pharmacokinetic variability of tacrolimus is required.

Development of a dosage strategy in patients receiving enoxaparin by continuous intravenous infusion using modelling and simulation

Aim To develop an appropriate dosing strategy for continuous intravenous infusions (CII) of enoxaparin by minimizing the percentage of steady-state anti-Xa concentration (C-ss) outside the therapeutic range of 0.5-1.2 IU ml(-1). Methods A nonlinear mixed effects model was developed with NONMEM (R) for 48 adult patients who received CII of enoxaparin with infusion durations that ranged from 8 to 894 h at rates between 100 and 1600 IU h(-1). Three hundred and sixty-three anti-Xa concentration measurements were available from patients who received CII. These were combined with 309 anti-Xa concentrations from 35 patients who received subcutaneous enoxaparin. The effects of age, body size, height, sex, creatinine clearance (CrCL) and patient location [intensive care unit (ICU) or general medical unit] on pharmacokinetic (PK) parameters were evaluated. Monte Carlo simulations were used to (i) evaluate covariate effects on C-ss and (ii) compare the impact of different infusion rates on predicted C-ss. The best dose was selected based on the highest probability that the C-ss achieved would lie within the therapeutic range. Results A two-compartment linear model with additive and proportional residual error for general medical unit patients and only a proportional error for patients in ICU provided the best description of the data. Both CrCL and weight were found to affect significantly clearance and volume of distribution of the central compartment, respectively. Simulations suggested that the best doses for patients in the ICU setting were 50 IU kg(-1) per 12 h (4.2 IU kg(-1) h(-1)) if CrCL < 30 ml min(-1); 60 IU kg(-1) per 12 h (5.0 IU kg(-1) h(-1)) if CrCL was 30-50 ml min(-1); and 70 IU kg(-1) per 12 h (5.8 IU kg(-1) h(-1)) if CrCL > 50 ml min(-1). The best doses for patients in the general medical unit were 60 IU kg(-1) per 12 h (5.0 IU kg(-1) h(-1)) if CrCL < 30 ml min(-1); 70 IU kg(-1) per 12 h (5.8 IU kg(-1) h(-1)) if CrCL was 30-50 ml min(-1); and 100 IU kg(-1) per 12 h (8.3 IU kg(-1) h(-1)) if CrCL > 50 ml min(-1). These best doses were selected based on providing the lowest equal probability of either being above or below the therapeutic range and the highest probability that the C-ss achieved would lie within the therapeutic range. Conclusion The dose of enoxaparin should be individualized to the patients' renal function and weight. There is some evidence to support slightly lower doses of CII enoxaparin in patients in the ICU setting.

Prophylactic ranitidine treatment in critically ill children - a population pharmacokinetic study

Aims - To characterize the population pharmacokinetics of ranitidine in critically ill children and to determine the influence of various clinical and demographic factors on its disposition. Methods - Data were collected prospectively from 78 paediatric patients (n = 248 plasma samples) who received oral or intravenous ranitidine for prophylaxis against stress ulcers, gastrointestinal bleeding or the treatment of gastro-oesophageal reflux. Plasma samples were analysed using high-performance liquid chromatography, and the data were subjected to population pharmacokinetic analysis using nonlinear mixed-effects modelling. Results - A one-compartment model best described the plasma concentration profile, with an exponential structure for interindividual errors and a proportional structure for intra-individual error. After backward stepwise elimination, the final model showed a significant decrease in objective function value (−12.618; P < 0.001) compared with the weight-corrected base model. Final parameter estimates for the population were 32.1 l h−1 for total clearance and 285 l for volume of distribution, both allometrically modelled for a 70 kg adult. Final estimates for absorption rate constant and bioavailability were 1.31 h−1 and 27.5%, respectively. No significant relationship was found between age and weight-corrected ranitidine pharmacokinetic parameters in the final model, with the covariate for cardiac failure or surgery being shown to reduce clearance significantly by a factor of 0.46. Conclusions - Currently, ranitidine dose recommendations are based on children's weights. However, our findings suggest that a dosing scheme that takes into consideration both weight and cardiac failure/surgery would be more appropriate in order to avoid administration of higher or more frequent doses than necessary.

Population pharmacokinetic and pharmacogenetic analysis of 6-mercaptopurine in paediatric patients with acute lymphoblastic leukaemia

WHAT IS ALREADY KNOWN ABOUT THIS SUBJECT • The cytotoxic effects of 6-mercaptopurine (6-MP) were found to be due to drug-derived intracellular metabolites (mainly 6-thioguanine nucleotides and to some extent 6-methylmercaptopurine nucleotides) rather than the drug itself. • Current empirical dosing methods for oral 6-MP result in highly variable drug and metabolite concentrations and hence variability in treatment outcome. WHAT THIS STUDY ADDS • The first population pharmacokinetic model has been developed for 6-MP active metabolites in paediatric patients with acute lymphoblastic leukaemia and the potential demographic and genetically controlled factors that could lead to interpatient pharmacokinetic variability among this population have been assessed. • The model shows a large reduction in interindividual variability of pharmacokinetic parameters when body surface area and thiopurine methyltransferase polymorphism are incorporated into the model as covariates. • The developed model offers a more rational dosing approach for 6-MP than the traditional empirical method (based on body surface area) through combining it with pharmacogenetically guided dosing based on thiopurine methyltransferase genotype. AIMS - To investigate the population pharmacokinetics of 6-mercaptopurine (6-MP) active metabolites in paediatric patients with acute lymphoblastic leukaemia (ALL) and examine the effects of various genetic polymorphisms on the disposition of these metabolites. METHODS - Data were collected prospectively from 19 paediatric patients with ALL (n = 75 samples, 150 concentrations) who received 6-MP maintenance chemotherapy (titrated to a target dose of 75 mg m−2 day−1). All patients were genotyped for polymorphisms in three enzymes involved in 6-MP metabolism. Population pharmacokinetic analysis was performed with the nonlinear mixed effects modelling program (nonmem) to determine the population mean parameter estimate of clearance for the active metabolites. RESULTS - The developed model revealed considerable interindividual variability (IIV) in the clearance of 6-MP active metabolites [6-thioguanine nucleotides (6-TGNs) and 6-methylmercaptopurine nucleotides (6-mMPNs)]. Body surface area explained a significant part of 6-TGNs clearance IIV when incorporated in the model (IIV reduced from 69.9 to 29.3%). The most influential covariate examined, however, was thiopurine methyltransferase (TPMT) genotype, which resulted in the greatest reduction in the model's objective function (P < 0.005) when incorporated as a covariate affecting the fractional metabolic transformation of 6-MP into 6-TGNs. The other genetic covariates tested were not statistically significant and therefore were not included in the final model. CONCLUSIONS - The developed pharmacokinetic model (if successful at external validation) would offer a more rational dosing approach for 6-MP than the traditional empirical method since it combines the current practice of using body surface area in 6-MP dosing with a pharmacogenetically guided dosing based on TPMT genotype.