Population genetic structure and demography analysis in palinurus elephas populations from atlantic and mediterranean waters

Report for the scientific sojourn at the University of Reading, United Kingdom, from January until May 2008. The main objectives have been firstly to infer population structure and parameters in demographic models using a total of 13 microsatellite loci for genotyping approximately 30 individuals per population in 10 Palinurus elephas populations both from Mediterranean and Atlantic waters. Secondly, developing statistical methods to identify discrepant loci, possibly under selection and implement those methods using the R software environment. It is important to consider that the calculation of the probability distribution of the demographic and mutational parameters for a full genetic data set is numerically difficult for complex demographic history (Stephens 2003). The Approximate Bayesian Computation (ABC), based on summary statistics to infer posterior distributions of variable parameters without explicit likelihood calculations, can surmount this difficulty. This would allow to gather information on different demographic prior values (i.e. effective population sizes, migration rate, microsatellite mutation rate, mutational processes) and assay the sensitivity of inferences to demographic priors by assuming different priors.

Life tables and reproductive parameters of Lutzomyia spinicrassa (Diptera: Psychodidae) under laboratory conditions

Lutzomyia spinicrassa is a vector of Leishmania braziliensis in Colombia. This sand fly has a broad geographical distribution in Colombia and Venezuela and it is found mainly in coffee plantations. Baseline biological growth data of L. spinicrassa were obtained under experimental laboratory conditions. The development time from egg to adult ranged from 59 to 121 days, with 12.74 weeks in average. Based on cohorts of 100 females, horizontal life table was constructed. The following predictive parameters were obtained: net rate of reproduction (8.4 females per cohort female), generation time (12.74 weeks), intrinsic rate of population increase (0.17), and finite rate of population increment (1.18). The reproductive value for each class age of the cohort females was calculated. Vertical life tables were elaborated and mortality was described for the generation obtained of the field cohort. In addition, for two successive generations, additive variance and heritability for fecundity were estimated.

Evolution of life history traits in a natural population of the shrew Crocidura russula

Abstract Life history traits encompass all the decisions concerning fitness an individual is faced with during his life. The study of these traits is crucial to understand the factors shaping the biology of living organisms. Up until now, most of the information on the evolution of life history traits comes from laboratory studies. While these studies are interesting to test the effect of specific parameters, their conclusions are difficult to extrapolate to natural populations. Investigating the evolution of life history traits in natural populations is of great interest. This may be tricky because it requires information on reproduction, survival and morphology of individuals. Mark-recapture methods allow most of this information to be obtained. However, when direct observations of a species are not possible due to its ecology, indirect methods must be used to infer lifetime reproductive success. In this case, molecular markers are particularly helpful in assessing the genetic relationships between individuals and allow the construction of a pedigree. This thesis focuses on a natural population of a small insectivorous mammal, the greater white-toothed shrew, Crocidura russula. Because of its hidden lifestyle, the two complementary techniques mentioned above were combined to gather information on this population. The data were used to explore diverse aspects of evolutionary biology. We demonstrated that the high genetic variance displayed by the species was not maintained by its mating system because this shrew was less monogamous than previously thought. The large genetic diversity was most likely promoted by gene flow from the neighborhood. Dispersal was thus a central topic in this thesis. We showed that dispersal was not driven by inbreeding avoidance. In addition, we did not find any inbreeding depression in the population. Dispersal was promoted by a high number of vacant territories in the population for both sexes, meaning that territory acquisition played an important role in driving dispersal. Moreover, dispersal propensity was shown to have a genetic basis and, once achieved, to have no effect on individual fitness. Body mass was found to be a life history trait strongly influenced by sexual and viability selection in both sexes. Larger individuals had higher access to reproduction through territory acquisition and defense than lighter ones. By contrast, intermediate size individuals were favored by viability selection presumably because of ecological constraints and metabolic costs. Finally, we demonstrated that the majority of the life history traits in our shrew population has the potential to evolve because they maintained substantial amounts of additive genetic variance. Nonetheless, life history traits had no significant heritability due to their high level of nonadditive or environmental variance. Résumé Les traits d'histoire de vie comprennent toutes les décisions auxquelles un individu est confronté au cours de sa vie et qui concernent sa valeur adaptative. L'étude de ces traits est cruciale pour comprendre les facteurs qui façonnent la biologie des êtres vivants. Jusqu'à ce jour, la majorité des informations sur l'évolution des traits d'histoire de vie provient d'études réalisées en laboratoire. Alors que ces études sont intéressantes pour tester l'effet de paramètres spécifiques, leurs conclusions sont difficilement extrapolables aux populations naturelles. Il est particulièrement intéressant d'étudier l'évolution des traits d'histoire de vie dans des populations naturelles. Toutefois, ces études peuvent se révéler difficiles parce qu'elles requièrent des informations sur la reproduction, la survie et la morphologie des individus. Des méthodes de marquage-recapture permettent d'obtenir ces informations. Cependant, lorsque l'écologie de l'espèce rend les obervations directes impossibles, des méthodes indirectes doivent être utilisées pour obtenir le succès reproducteur des individus. Dans ce cas, les marqueurs moléculaires sont particulièrement utiles pour évaluer les relations génétiques entre individus et permettre la construction d'un pedigree. Cette thèse porte sur une population naturelle d'un petit mammifère insectivore, la musaraigne musette, Crocidura russula. Parce que cette espèce présente un mode de vie souterrain, les deux techniques complémentaires mentionnées ci-dessus ont été combinées pour acquérir les informations nécessaires. Les données ont été utilisées pour explorer divers aspects de biologie evolutive. Nous avons montré que la grande quantité de variance génétique trouvée chez cette espèce n'est pas maintenue par son système d'appariement. Celle-ci s'est en effet avérée être moins monogame que ce qui était admis jusqu'ici. Sa grande diversité génétique est plutôt entretenue par le flux de gènes provenant du voisinage. La dispersion a donc été un sujet phare dans cette thèse. Nous avons montré qu'elle n'est pas provoquée par un évitement de la consanguinité et nous n'avons pas trouvé de dépression de consanguité dans notre population. L'acquisition d'un territoire joue par contre un rôle important dans la dispersion. En outre, la dispersion possède une base génétique chez cette espèce. De plus, une fois qu'ils ont dispersé, les individus n'ont pas une valeur adaptative differente d'individus philopatriques. Le poids s'est avéré être un trait d'histoire de vie fortement influencé par la sélection sexuelle et de viabilité chez les deux sexes. Les gros individus ont accès à la reproduction parce qu'ils acquièrent et défendent un territoire plus facilement que les plus légers. Au contraire, les individus de taille intermédiaire sont favorisés par la sélection de viabilité, certainement à cause de contraintes écologiques et de coûts métaboliques. Finalement, nous avons montré que la majorité des traits d'histoire de vie dans notre population a le potentiel d'évoluer parce qu'elle maintient des quantités considérables de variance génétique additive. Néanmoins, l'héritabilité de ces traits d'histoire de vie n'est pas significative à cause de la grande quantité de variance non-additive ou environmentale associée à ces traits.

Population pharmacokinetics of imatinib in CML and GIST patients under long-term treatment

Imatinib (Glivec®) has transformed the treatment and short-term prognosis of chronic myeloid leukaemia (CML) and gastrointestinal stromal tumour (GIST). However, the treatment must be taken indefinitely and is not devoid of inconvenience and toxicity. Moreover, resistance or escape from disease control occurs in a significant number of patients. Imatinib is a substrate of the cytochromes P450 CYP3A4/5 and of the multidrug transporter P glycoprotein (product of the MDR1 gene), and is also bound to the alpha1-acid glycoprotein (AAG) in plasma. Considering the large inter-individual differences in the expression and function of those systems, the disposition and clinical activity of imatinib can be expected to vary widely among patients, calling for dosage individualisation. The aim of this exploratory study was to determine the average pharmacokinetic parameters characterizing the disposition of imatinib in the target population, to assess their inter-individual variability, and to identify influential factors affecting them. A total of 321 plasma concentrations were measured in 59 patients receiving Glivec® at diverse dosage regimens, using a validated chromatographic method developed for this study. The results were analysed by non-linear mixed effect modelling (NONMEM). A one-compartment model with first-order absorption described the data appropriately, with an average apparent clearance of 12.4 l/h, a volume of distribution of 268 l and an absorption constant of 0.47 h-1. The clearance was affected by body weight, age and sex. No influences of interacting drugs were found. DNA samples were used for pharmacogenetic explorations. The MDR1 polymorphism 3435C>T and the AAG phenotype appears to modulate the disposition of imatinib. Large inter-individual variability (CV %) remained unexplained by the demographic covariates considered, both on clearance (40%) and distribution volume (71%). Together with intra-patient variability (34%), this translates into an 8-fold width of the 90%-prediction interval of plasma concentrations expected under a fixed dosing regimen. This is a strong argument to further investigate the possible usefulness of a therapeutic drug monitoring programme for imatinib. It may help in individualising the dosing regimen before overt disease progression or observation of treatment toxicity, thus improving both the long-term therapeutic effectiveness and tolerability of this drug.

Population pharmacokinetics of mefloquine, piperaquine and artemether-lumefantrine in Cambodian and Tanzanian malaria patients.

BACKGROUND: Inter-individual variability in plasma concentration-time profiles might contribute to differences in anti-malarial treatment response. This study investigated the pharmacokinetics of three different forms of artemisinin combination therapy (ACT) in Tanzania and Cambodia to quantify and identify potential sources of variability. METHODS: Drug concentrations were measured in 143 patients in Tanzania (artemether, dihydroartemisinin, lumefantrine and desbutyl-lumefantrine), and in 63 (artesunate, dihydroartemisinin and mefloquine) and 60 (dihydroartemisinin and piperaquine) patients in Cambodia. Inter- and intra-individual variabilities in the pharmacokinetic parameters were assessed and the contribution of demographic and other covariates was quantified using a nonlinear mixed-effects modelling approach (NONMEM®). RESULTS: A one-compartment model with first-order absorption from the gastrointestinal tract fitted the data for all drugs except piperaquine (two-compartment). Inter-individual variability in concentration exposure was about 40% and 12% for mefloquine. From all the covariates tested, only body weight (for all antimalarials) and concomitant treatment (for artemether only) showed a significant influence on these drugs' pharmacokinetic profiles. Artesunate and dihydroartemisinin could not be studied in the Cambodian patients due to insufficient data-points. Modeled lumefantrine kinetics showed that the target day 7 concentrations may not be achieved in a substantial proportion of patients. CONCLUSION: The marked variability in the disposition of different forms of ACT remained largely unexplained by the available covariates. Dosing on body weight appears justified. The concomitance of unregulated drug use (residual levels found on admission) and sub-optimal exposure (variability) could generate low plasma levels that contribute to selecting for drug-resistant parasites.

Nutritional status using the Mini Nutritional Assessment questionnaire and its relationship with bone quality in a population of institutionalized elderly women.

Malnutrition, a risk factor for osteoporotic fractures, is frequent in elderly people and, is underdiagnosed and undertreated. There are only few studies on the nutritional status of elderly people in Europe. The Mini Nutritional Assessment (MNA) is a non invasive and validated questionnaire to evaluate nutritional status in elderly people, classified in three groups: 1 degree score < 17: malnourished, 2 degrees score >17 and < 24: at risk of malnutrition, 3 degrees score >24: well-nourished, with a maximum of 30 points. Quantitative ultrasound of bone (QUS) is a method for assessing quality of bone which can be easily performed in nursing homes. Therefore, these two tests allowed to study the relationships between nutritional status and ultrasonic parameters of bone in 78 institutionalized women aged 86 +/- 6 years, living in 11 nursing homes around Lausanne (Switzerland). All were assessed by the MNA, had a measurement of the tricipital skin fold and of the grip strength. Functional status was evaluated by the scale "Activity of Daily Living" (ADL), and serum albumin level was measured when permitted. All had QUS of the calcaneus (with an Achilles, GE Lunar). The measured parameters are the Broadband Ultrasound Attenuation (BUA), attenuation of a band of ultrasonic frequencies through the medium, expressed in dB/MHz, and the Speed of Sound (SOS), speed of the ultrasounds through the medium, expressed in m/s. A third parameter, the stiffness index (SI), expressed as a percentage of the values obtained by the manufacturer in a young population and derived from BUA and SOS, was calculated automatically : SI = (0.67xBUA) + (0.28xSOS) - 420, expressed in percent compared to a young adult population (%YA). Fifteen percent of the women were undernourished and 58% were at risk of malnutrition. As expected, compared with the well-nourished minority, undernourished subjects had significant lower body mass index (BMI), tricipital skin fold (TSF), ADL score and albumin level (p < 0,01). The subjects "at risk of malnutrition" had significant lower BMI, ADL score (p < 0.01), tricipital skin fold and serum albumin (p < 0.05). Ultrasound parameters were low independently of the nutritional status. MNA score correlated significantly with tricipital skin fold (r = 0.508, p < 0.01), ADL (r = 0.538, p < 0.01) and albumin serum level (r = 0.409, p = 0.01). There was a trend for a correlation between the MNA and the ultrasound parameter BUA (r = 0.207, p = 0.07), whereas no correlation was found with SOS and SI. A multivariate analysis showed that tricipital skin fold and ADL explained 61% of the variance of the MNA. In conclusion, using simple and non invasive methods, this study showed that malnutrition and osteoporosis are frequent in institutionalized elderly persons in our country, and the ultrasound parameters are influenced by many others factors in addition to nutrition, especially at this age and in elderly residents of nursing homes.

Electrocardiographic findings in a middle-aged African population in the Seychelles islands.

This study describes major electrocardiogram (ECG) measurements and diagnoses in a population of African individuals; most reference data have been collected in Caucasian populations and evidence exists for interethnic differences in ECG findings. This study was conducted in the Seychelles islands (Indian Ocean) and included 709 black individuals (343 men and 366 women) aged 25 to 64 years randomly selected from the general population. Resting ECG were recorded by using a validated ECG unit equipped with a measurement and interpretation software (Cardiovit AT-6, Schiller, Switzerland). The epidemiology of 14 basic ECG measurements, 6 composite criteria for left ventricular hypertrophy and 19 specific ECG diagnoses including abnormal rhythms, conduction abnormalities, repolarization abnormalities, and myocardial infarction were examined. Substantial gender and age differences were found for several ECG parameters. Moreover, tracings recorded in African individuals of the Seychelles differed from those collected similarly in Caucasian populations in many respects. For instance, heart rate was approximately 5 beats per minute lower in the African individuals than in selected Caucasian populations, prevalence of first degree atrio-ventricular block was especially high (4.8%), and the average Sokolow-Lyon voltage was markedly higher in African individuals of the Seychelles compared with black and white Americans. The integrated interpretation software detected "old myocardial infarction" in 3.8% of men and 0% of women and "old myocardial infarction possible" in 6.1% and 3%, respectively. Cardiac infarction injury scores are also provided. In conclusion, the study provides reference values for ECG findings in a specific population of people of African descent and stresses the need to systematically consider gender, age, and ethnicity when interpreting ECG tracings in individuals.

Relationship among epidemiological parameters of six childhood infections in a non-immunized Brazilian community

Epidemiological parameters, such as age-dependent force of infection and average age at infection () were estimated for rubella, varicella, rotavirus A, respiratory syncytial virus, hepatitis A and parvovirus B19 infections for a non-immunized Brazilian community, using the same sera samples. The for the aforementioned diseases were 8.45 years (yr) [95% CI: (7.23, 9.48) yr], 3.90 yr [95% CI: (3.51, 4.28) yr], 1.03 yr [95% CI: (0.96, 1.09) yr], 1.58 yr [95% CI: (1.39, 1.79) yr], 7.17 yr [95% CI: (6.48, 7.80) yr] and 7.43 yr [95% CI: (5.68, 9.59) yr], respectively. The differences between average ages could be explained by factors such as differences in the effectiveness of the protection conferred to newborns by maternally derived antibodies, competition between virus species and age-dependent host susceptibility. Our seroprevalence data may illustrate a case of the above-mentioned mechanisms working together within the same population.

Population pharmacokinetics and effects of efavirenz in patients with human immunodeficiency virus infection.

OBJECTIVE: The reverse transcriptase inhibitor efavirenz is currently used at a fixed dose of 600 mg/d. However, dosage individualization based on plasma concentration monitoring might be indicated. This study aimed to assess the efavirenz pharmacokinetic profile and interpatient versus intrapatient variability in patients who are positive for human immunodeficiency virus, to explore the relationship between drug exposure, efficacy, and central nervous system toxicity and to build up a Bayesian approach for dosage adaptation. METHODS: The population pharmacokinetic analysis was performed by use of NONMEM based on plasma samples from a cohort of unselected patients receiving efavirenz. With the use of a 1-compartment model with first-order absorption, the influence of demographic and clinical characteristics on oral clearance and oral volume of distribution was examined. The average drug exposure during 1 dosing interval was estimated for each patient and correlated with markers of efficacy and toxicity. The population kinetic parameters and the variabilities were integrated into a Bayesian equation for dosage adaptation based on a single plasma sample. RESULTS: Data from 235 patients with a total of 719 efavirenz concentrations were collected. Oral clearance was 9.4 L/h, oral volume of distribution was 252 L, and the absorption rate constant was 0.3 h(-1). Neither the demographic covariates evaluated nor the comedications showed a clinically significant influence on efavirenz pharmacokinetics. A large interpatient variability was found to affect efavirenz relative bioavailability (coefficient of variation, 54.6%), whereas the intrapatient variability was small (coefficient of variation, 26%). An inverse correlation between average drug exposure and viral load and a trend with central nervous system toxicity were detected. This enabled the derivation of a dosing adaptation strategy suitable to bring the average concentration into a therapeutic target from 1000 to 4000 microg/L to optimize viral load suppression and to minimize central nervous system toxicity. CONCLUSIONS: The high interpatient and low intrapatient variability values, as well as the potential relationship with markers of efficacy and toxicity, support the therapeutic drug monitoring of efavirenz. However, further evaluation is needed before individualization of an efavirenz dosage regimen based on routine drug level monitoring should be recommended for optimal patient management.

Population pharmacokinetic-pharmacodynamic modelling of angiotensin receptor blockade in healthy volunteers.

OBJECTIVE: To compare the pharmacokinetic and pharmacodynamic characteristics of angiotensin II receptor antagonists as a therapeutic class. DESIGN: Population pharmacokinetic-pharmacodynamic modelling study. METHODS: The data of 14 phase I studies with 10 different drugs were analysed. A common population pharmacokinetic model (two compartments, mixed zero- and first-order absorption, two metabolite compartments) was applied to the 2685 drug and 900 metabolite concentration measurements. A standard nonlinear mixed effect modelling approach was used to estimate the drug-specific parameters and their variabilities. Similarly, a pharmacodynamic model was applied to the 7360 effect measurements, i.e. the decrease of peak blood pressure response to intravenous angiotensin challenge recorded by finger photoplethysmography. The concentration of drug and metabolite in an effect compartment was assumed to translate into receptor blockade [maximum effect (Emax) model with first-order link]. RESULTS: A general pharmacokinetic-pharmacodynamic (PK-PD) model for angiotensin antagonism in healthy individuals was successfully built up for the 10 drugs studied. Representatives of this class share different pharmacokinetic and pharmacodynamic profiles. Their effects on blood pressure are dose-dependent, but the time course of the effect varies between the drugs. CONCLUSIONS: The characterisation of PK-PD relationships for these drugs gives the opportunity to optimise therapeutic regimens and to suggest dosage adjustments in specific conditions. Such a model can be used to further refine the use of this class of drugs.

Inference of distributional parameters from compositional samples containing nondetects

Low concentrations of elements in geochemical analyses have the peculiarity of beingcompositional data and, for a given level of significance, are likely to be beyond thecapabilities of laboratories to distinguish between minute concentrations and completeabsence, thus preventing laboratories from reporting extremely low concentrations of theanalyte. Instead, what is reported is the detection limit, which is the minimumconcentration that conclusively differentiates between presence and absence of theelement. A spatially distributed exhaustive sample is employed in this study to generateunbiased sub-samples, which are further censored to observe the effect that differentdetection limits and sample sizes have on the inference of population distributionsstarting from geochemical analyses having specimens below detection limit (nondetects).The isometric logratio transformation is used to convert the compositional data in thesimplex to samples in real space, thus allowing the practitioner to properly borrow fromthe large source of statistical techniques valid only in real space. The bootstrap method isused to numerically investigate the reliability of inferring several distributionalparameters employing different forms of imputation for the censored data. The casestudy illustrates that, in general, best results are obtained when imputations are madeusing the distribution best fitting the readings above detection limit and exposes theproblems of other more widely used practices. When the sample is spatially correlated, itis necessary to combine the bootstrap with stochastic simulation

Pragmatic approach for interpreting antiretroviral drug concentrations based on a systematic review of population pharmacokinetic studies

Objectives: The study objective was to derive reference pharmacokinetic curves of antiretroviral drugs (ART) based on available population pharmacokinetic (Pop-PK) studies that can be used to optimize therapeutic drug monitoring guided dosage adjustment.¦Methods: A systematic search of Pop-PK studies of 8 ART in adults was performed in PubMed. To simulate reference PK curves, a summary of the PK parameters was obtained for each drug based on meta-analysis approach. Most models used one-compartment model, thus chosen as reference model. Models using bi-exponential disposition were simplified to one-compartment, since the first distribution phase was rapid and not determinant for the description of the terminal elimination phase, mostly relevant for this project. Different absorption were standardized for first-order absorption processes.¦Apparent clearance (CL), apparent volume of distribution of the terminal phase (Vz) and absorption rate constant (ka) and inter-individual variability were pooled into summary mean value, weighted by number of plasma levels; intra-individual variability was weighted by number of individuals in each study.¦Simulations based on summary PK parameters served to construct concentration PK percentiles (NONMEM®).¦Concordance between individual and summary parameters was assessed graphically using Forest-plots. To test robustness, difference in simulated curves based on published and summary parameters was calculated using efavirenz as probe drug.¦Results: CL was readily accessible from all studies. For studies with one-compartment, Vz was central volume of distribution; for two-compartment, Vz was CL/λz. ka was directly used or derived based on the mean absorption time (MAT) for more complicated absorption models, assuming MAT=1/ka.¦The value of CL for each drug was in excellent agreement throughout all Pop-PK models, suggesting that minimal concentration derived from summary models was adequately characterized. The comparison of the concentration vs. time profile for efavirenz between published and summary PK parameters revealed not more than 20% difference. Although our approach appears adequate for estimation of elimination phase, the simplification of absorption phase might lead to small bias shortly after drug intake.¦Conclusions: Simulated reference percentile curves based on such an approach represent a useful tool for interpretating drug concentrations. This Pop-PK meta-analysis approach should be further validated and could be extended to elaborate more sophisticated computerized tool for the Bayesian TDM of ART.

Vertebral fractures in patients with inflammatory bowel disease compared with a healthy population: a prospective case-control study.

BACKGROUND A prospective study was performed to compare the prevalence of morphometric vertebral fractures (MVF) between patients with inflammatory bowel disease (IBD) and healthy subjects and to identify predictive factors of fracture. METHODS A total of 107 patients with IBD (53 with Crohn's disease and 54 with ulcerative colitis) and 51 healthy subjects participated in the study. Information about anthropometric parameters, toxins, previous fractures, and parameters related to this disease were evaluated. The index of vertebral deformity, bone mass density (BMD), and biochemical parameters were calculated. RESULTS A total of 72 fractures were detected in 38.32% of patients with IBD, and 10 fractures were detected in 13.73% of healthy subjects; the risk of fracture in patients with IBD was higher than that in control subjects (OR, 4.03; 95% CI, 1.652-9.847; p < 0.002). We found no correlation between fracture and BMD in patients with IBD (lumbar spine, r = -0.103, p = 0.17 and femoral neck, r = -0.138, p = 0.07). Corticosteroid treatment was not associated with prevalent vertebral fractures nor with taking corticosteroids (r = 0.135, p = 0.14) or the duration for which they were taken (r = 0.08, p = 0.38), whereas this relationship was present in the controls (r = -0.365, p = 0.01). In the multivariate analysis, none of the measured parameters were significantly predictive of fracture, only to manifested IBD. Hypovitaminosis D was observed in 55.14% of patients with IBD. CONCLUSIONS The prevalence of morphometric vertebral fractures is higher in patients with IBD than in the healthy population, without association with BMD or corticoid treatment. Simply having IBD was proven to be a predictive factor of fracture. We observed a high incidence of hypovitaminosis D in patients with IBD.

Deficient selenium status of a healthy adult Spanish population.

INTRODUCTION Selenium is an essential micronutrient for human health, being a cofactor for enzymes with antioxidant activity that protect the organism from oxidative damage. An inadequate intake of this mineral has been associated with the onset and progression of chronic diseases such as hypertension, diabetes, coronary diseases, asthma, and cancer. For this reason, knowledge of the plasma and erythrocyte selenium levels of a population makes a relevant contribution to assessment of its nutritional status. OBJECTIVE The objective of the present study was to determine the nutritional status of selenium and risk of selenium deficiency in a healthy adult population in Spain by examining food and nutrient intake and analyzing biochemical parameters related to selenium metabolism, including plasma and erythrocyte levels and selenium-dependent glutathione peroxidase (GPx) enzymatic activity. MATERIAL AND METHODS We studied 84 healthy adults (31 males and 53 females) from the province of Granada, determining their plasma and erythrocyte selenium concentrations and the association of these levels with the enzymatic activity of glutathione peroxidase (GPx) and with life style factors. We also gathered data on their food and nutrient intake and the results of biochemical analyses. Correlations were studied among all of these variables. RESULTS The mean plasma selenium concentration was 76.6 ± 17.3 μg/L (87.3 ± 17.4 μg/L in males, 67.3 ± 10.7 μg/L in females), whereas the mean erythrocyte selenium concentration was 104.6 μg/L (107.9 ± 26.1 μg/L in males and 101.7 ± 21.7 μg/L in females). The nutritional status of selenium was defined by the plasma concentration required to reach maximum GPx activity, establishing 90 μg/L as reference value. According to this criterion, 50% of the men and 53% of the women were selenium deficient. CONCLUSIONS Selenium is subjected to multiple regulation mechanisms. Erythrocyte selenium is a good marker of longer term selenium status, while plasma selenium appears to be a marker of short-term nutritional status. The present findings indicate a positive correlation between plasma selenium concentration and the practice of physical activity. Bioavailability studies are required to establish appropriate reference levels of this mineral for the Spanish population.

Association of HbA1c and cardiovascular and renal disease in an adult Mediterranean population

BACKGROUND: Increasing evidence suggests a mechanistic link between the glycemic environment and renal and cardiovascular events, even below the threshold for diabetes. We aimed to assess the association between HbA1c and chronic kidney disease (CKD) and cardiovascular disease (CVD). METHODS: A cross-sectional study involving a random representative sample of 2270 adults from southern Spain (Malaga) was undertaken. We measured HbA1c, serum creatinine and albuminuria in fasting blood and urine samples. RESULTS: Individuals without diabetes in the upper HbA1c tertile had an unfavorable cardiovascular and renal profile and shared certain clinical characteristics with the patients with diabetes. Overall, a higher HbA1c concentration was strongly associated with CKD or CVD after adjustment for traditional risk factors. The patients with known diabetes had a 2-fold higher odds of CKD or CVD. However, when both parameters were introduced in the same model, the HbA1c concentration was only significantly associated with clinical endpoints (OR: 1.4, 95% CI, 1.1-1.6, P = 0.002). An increase in HbA1c of one percentage point was associated with a 30% to 40% increase in the rate of CKD or CVD. This relationship was apparent in persons with and without known diabetes. ROC curves illustrated that a HbA1c of 37 mmol/mol (5.5%) was the optimal value in terms of sensitivity and specificity for predicting endpoints in this population. CONCLUSION: HbA1c levels were associated with a higher prevalence of CKD and CVD cross-sectionally, regardless of diabetes status. These data support the value of HbA1c as a marker of cardiovascular and renal disease in the general population.