DiveRsity: An R package for the estimation and exploration of population genetics parameters and their associated errors

Summary: We present a new R package, diveRsity, for the calculation of various diversity statistics, including common diversity partitioning statistics (?, G) and population differentiation statistics (D, GST ', ? test for population heterogeneity), among others. The package calculates these estimators along with their respective bootstrapped confidence intervals for loci, sample population pairwise and global levels. Various plotting tools are also provided for a visual evaluation of estimated values, allowing users to critically assess the validity and significance of statistical tests from a biological perspective. diveRsity has a set of unique features, which facilitate the use of an informed framework for assessing the validity of the use of traditional F-statistics for the inference of demography, with reference to specific marker types, particularly focusing on highly polymorphic microsatellite loci. However, the package can be readily used for other co-dominant marker types (e.g. allozymes, SNPs). Detailed examples of usage and descriptions of package capabilities are provided. The examples demonstrate useful strategies for the exploration of data and interpretation of results generated by diveRsity. Additional online resources for the package are also described, including a GUI web app version intended for those with more limited experience using R for statistical analysis. © 2013 British Ecological Society.

Space and phase resolved plasma parameters in an industrial dual-frequency capacitively coupled radio-frequency discharge

The dynamics of high energetic electrons (>= 11.7 eV) in a modified industrial confined dual-frequency capacitively coupled RF discharge (Exelan, Lam Research Inc.), operated at 1.937 MHz and 27.118 MHz, is investigated by means of phase resolved optical emission spectroscopy. Operating in a He-O-2. plasma with small rare gas admixtures the emission is measured, with one-dimensional spatial resolution along the discharge axis. Both the low and high frequency RF cycle are resolved. The diagnostic is based on time dependent measurements of the population densities of specifically chosen excited rare gas states. A time dependent model, based on rate equations, describes the dynamics of the population densities of these levels. Based on this model and the comparison of the excitation of various rare gas states, with different excitation thresholds, time and space resolved electron temperature, propagation velocity and qualitative electron density as well as electron energy distribution functions are determined. This information leads to a better understanding of the dual-frequency sheath dynamics and shows, that separate control of ion energy and electron density is limited.

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

WHAT IS ALREADY KNOWN ABOUT THIS SUBJECT center dot 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. center dot 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 center dot 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. center dot 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. center dot 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. 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. 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. 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

Simplification of structured population dynamics in complex ecological communities

Goldstone's idea of slow dynamics resulting from spontaneously broken symmetries is applied to Hubbell's neutral hypothesis of community dynamics, to efficiently simplify stage-structured multi-species models-introducing the quasi-neutral approximation (QNA). Rather than assuming population-dynamical neutrality in the QNA, deviations from ideal neutrality, thought to be small, drive dynamics. The QNA is systematically derived to first and second order in a two-scale singular perturbation expansion. The total reproductive value of species, as computed from the effective life-history parameters resulting from the non-linear interactions with the surrounding community, emerges as the new dynamic variables in this aggregated description. Using a simple stage-structured community-assembly model, the QNA is demonstrated to accurately reproduce population dynamics in large, complex communities. Further, the utility of the QNA in building intuition for management problems is illustrated by estimating the responses of a fish stock to harvesting and variations in fecundity.

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.1lh for total clearance and 285l for volume of distribution, both allometrically modelled for a 70kg adult. Final estimates for absorption rate constant and bioavailability were 1.31h 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.

Potassium canrenoate treatment in paediatric patients: a population pharmacokinetic study using novel dried blood spot sampling

Objective: To characterize the population pharmacokinetics of canrenone following administration of potassium canrenoate (K-canrenoate) in paediatric patients.

Methods: Data were collected prospectively from 37 paediatric patients (median weight 2.9?kg, age range 2 days–0.85 years) who received intravenous K-canrenoate for management of retained fluids, for example in heart failure and chronic lung disease. Dried blood spot (DBS) samples (n?=?213) from these were analysed for canrenone content and the data subjected to pharmacokinetic analysis using nonlinear mixed-effects modelling. Another group of patients (n?=?16) who had 71 matching plasma and DBS samples was analysed separately to compare canrenone pharmacokinetic parameters obtained using the two different matrices.

Results: A one-compartment model best described the DBS data. Significant covariates were weight, postmenstrual age (PMA) and gestational age. The final population models for canrenone clearance (CL/F) and volume of distribution (V/F) in DBS were CL/F (l/h)?=?12.86?×? (WT/70.0)0.75?×?e [0.066?×? (PMA?-?40]) and V/F (l)?=?603.30?×? (WT/70)?×?(GA/40)1.89 where weight is in kilograms. The corresponding values of CL/F and V/F in a patient with a median weight of 2.9?kg are 1.11?l/h and 20.48?l, respectively. Estimated half-life of canrenone based on DBS concentrations was similar to that based on matched plasma concentrations (19.99 and 19.37?h, respectively, in 70?kg patient).

Conclusion: The range of estimated CL/F in DBS for the study population was 0.12–9.62?l/h; hence, bodyweight-based dosage adjustment of K-canrenoate appears necessary. However, a dosing scheme that takes into consideration both weight and age (PMA/gestational age) of paediatric patients seems more appropriate.

Molecular pedigree reconstruction and estimation of evolutionary parameters in a wild Atlantic salmon river system with incomplete sampling: a power analysis

Background: Pedigree reconstruction using genetic analysis provides a useful means to estimate fundamental population biology parameters relating to population demography, trait heritability and individual fitness when combined with other sources of data. However, there remain limitations to pedigree reconstruction in wild populations, particularly in systems where parent-offspring relationships cannot be directly observed, there is incomplete sampling of individuals, or molecular parentage inference relies on low quality DNA from archived material. While much can still be inferred from incomplete or sparse pedigrees, it is crucial to evaluate the quality and power of available genetic information a priori to testing specific biological hypotheses. Here, we used microsatellite markers to reconstruct a multi-generation pedigree of wild Atlantic salmon (Salmo salar L.) using archived scale samples collected with a total trapping system within a river over a 10 year period. Using a simulation-based approach, we determined the optimal microsatellite marker number for accurate parentage assignment, and evaluated the power of the resulting partial pedigree to investigate important evolutionary and quantitative genetic characteristics of salmon in the system.

Results: We show that at least 20 microsatellites (ave. 12 alleles/locus) are required to maximise parentage assignment and to improve the power to estimate reproductive success and heritability in this study system. We also show that 1.5 fold differences can be detected between groups simulated to have differing reproductive success, and that it is possible to detect moderate heritability values for continuous traits (h(2) similar to 0.40) with more than 80% power when using 28 moderately to highly polymorphic markers.

Conclusion: The methodologies and work flow described provide a robust approach for evaluating archived samples for pedigree-based research, even where only a proportion of the total population is sampled. The results demonstrate the feasibility of pedigree-based studies to address challenging ecological and evolutionary questions in free-living populations, where genealogies can be traced only using molecular tools, and that significant increases in pedigree assignment power can be achieved by using higher numbers of markers.

An improved TLBO with elite strategy for parameters identification of PEM fuel cell and solar cell models

Clean and renewable energy generation and supply has drawn much attention worldwide in recent years, the proton exchange membrane (PEM) fuel cells and solar cells are among the most popular technologies. Accurately modeling the PEM fuel cells as well as solar cells is critical in their applications, and this involves the identification and optimization of model parameters. This is however challenging due to the highly nonlinear and complex nature of the models. In particular for PEM fuel cells, the model has to be optimized under different operation conditions, thus making the solution space extremely complex. In this paper, an improved and simplified teaching-learning based optimization algorithm (STLBO) is proposed to identify and optimize parameters for these two types of cell models. This is achieved by introducing an elite strategy to improve the quality of population and a local search is employed to further enhance the performance of the global best solution. To improve the diversity of the local search a chaotic map is also introduced. Compared with the basic TLBO, the structure of the proposed algorithm is much simplified and the searching ability is significantly enhanced. The performance of the proposed STLBO is firstly tested and verified on two low dimension decomposable problems and twelve large scale benchmark functions, then on the parameter identification of PEM fuel cell as well as solar cell models. Intensive experimental simulations show that the proposed STLBO exhibits excellent performance in terms of the accuracy and speed, in comparison with those reported in the literature.

Evaluation of the Retinal Vasculature in Hypertension and Chronic Kidney Disease in an elderly population of Irish nuns

Background: Chronic kidney disease (CKD) and hypertension are global public health problems associated with considerable morbidity, premature mortality and attendant healthcare costs. Previous studies have highlighted that non-invasive examination of the retinal microcirculation can detect microvascular pathology that is associated with systemic disorders of the circulatory system such as hypertension. We examined the associations between retinal vessel caliber (RVC) and fractal dimension (DF), with both hypertension and CKD in elderly Irish nuns.

Methods: Data from 1233 participants in the cross-sectional observational Irish Nun Eye Study (INES) were assessed from digital photographs with a standardized protocol using computer-assisted software. Multivariate regression analyses were used to assess associations with hypertension and CKD, with adjustment for age, body mass index (BMI), refraction, fellow eye RVC, smoking, alcohol consumption, ischemic heart disease (IHD), cerebrovascular accident (CVA), diabetes and medication use.

Results: In total, 1122 (91%) participants (mean age: 76.3 [range: 56-100] years) had gradable retinal images of sufficient quality for blood vessel assessment. Hypertension was significantly associated with a narrower central retinal arteriolar equivalent (CRAE) in a fully adjusted analysis (P = 0.002; effect size= -2.16 μm; 95% confidence intervals [CI]: -3.51, -0.81 μm). No significant associations between other retinal vascular parameters and hypertension or between any retinal vascular parameters and CKD were found.

Conclusions: Individuals with hypertension have significantly narrower retinal arterioles which may afford an earlier opportunity for tailored prevention and treatment options to optimize the structure and function of the microvasculature, providing additional clinical utility. No significant associations between retinal vascular parameters and CKD were detected.