Comparison of two population pharmacokinetic programs, NONMEM and P-PHARM, for tacrolimus

Objectives: To compare the population modelling programs NONMEM and P-PHARM during investigation of the pharmacokinetics of tacrolimus in paediatric liver-transplant recipients. Methods: Population pharmacokinetic analysis was performed using NONMEM and P-PHARM on retrospective data from 35 paediatric liver-transplant patients receiving tacrolimus therapy. The same data were presented to both programs. Maximum likelihood estimates were sought for apparent clearance (CL/F) and apparent volume of distribution (V/F). Covariates screened for influence on these parameters were weight, age, gender, post-operative day, days of tacrolimus therapy, transplant type, biliary reconstructive procedure, liver function tests, creatinine clearance, haematocrit, corticosteroid dose, and potential interacting drugs. Results: A satisfactory model was developed in both programs with a single categorical covariate - transplant type - providing stable parameter estimates and small, normally distributed (weighted) residuals. In NONMEM, the continuous covariates - age and liver function tests - improved modelling further. Mean parameter estimates were CL/F (whole liver) = 16.3 1/h, CL/F (cut-down liver) = 8.5 1/h and V/F = 565 1 in NONMEM, and CL/F = 8.3 1/h and V/F = 155 1 in P-PHARM. Individual Bayesian parameter estimates were CL/F (whole liver) = 17.9 +/- 8.8 1/h, CL/F (cutdown liver) = 11.6 +/- 18.8 1/h and V/F = 712 792 1 in NONMEM, and CL/F (whole liver) = 12.8 +/- 3.5 1/h, CL/F (cut-down liver) = 8.2 +/- 3.4 1/h and V/F = 221 1641 in P-PHARM. Marked interindividual kinetic variability (38-108%) and residual random error (approximately 3 ng/ml) were observed. P-PHARM was more user friendly and readily provided informative graphical presentation of results. NONMEM allowed a wider choice of errors for statistical modelling and coped better with complex covariate data sets. Conclusion: Results from parametric modelling programs can vary due to different algorithms employed to estimate parameters, alternative methods of covariate analysis and variations and limitations in the software itself.

Sensitivity and specificity of pooled faecal culture and serology as flock-screening tests for detection of ovine paratuberculosis in Australia

The flock-level sensitivity of pooled faecal culture and serological testing using AGID for the detection of ovine Johne's disease-infected flocks were estimated using non-gold-standard methods. The two tests were compared in an extensive field trial in 296 flocks in New South Wales during 1998. In each flock, a sample of sheep was selected and tested for ovine Johne's disease using both the AGID and pooled faecal culture. The flock-specificity of pooled faecal culture also was estimated from results of surveillance and market-assurance testing in New South Wales. The overall flock-sensitivity of pooled faecal culture was 92% (95% CI: 82.4 and 97.4%) compared to 61% (50.5 and 70.9%) for serology (assuming that both tests were 100% specific). In low-prevalence flocks (estimated prevalence

A new species of Phyllurus (Lacertilia : Gekkonidae) and a revised phylogeny and key for the Australian leaf-tailed geckos

Phyllurus gulbaru, sp. nov., is a highly distinct species of leaf-tailed gecko restricted to rocky rainforest of Pattersons Gorge, north-west of Townsville. The possession of a cylindrical, non-depressed, tapering original and regenerated tail separates P. gulbaru from all congeners except P. caudiannulatus. From this species P. gulbaru is separated by having a partially divided, as opposed to fully divided, rostral scale. Furthermore, the very small spinose body tubercles of P. gulbaru are in marked contrast to the large spinose body scales of P. caudiannulatus. An analysis of 729 bp of mitochondrial 12S rRNA and cytochrome b genes reveals P. gulbaru to be a deeply divergent lineage with closer affinities to mid-east Queensland congeners than the geographically neighbouring P. amnicola on Mt Elliot. In conservation terms, P. gulbaru is clearly at risk. Field surveys of Pattersons Gorge and the adjacent ranges indicate that this species is restricted to a very small area of highly fragmented habitat, of which only a small proportion receives a degree of protection in State forest. Further, there is ongoing, unchecked destruction of dry rainforest habitat by fire. Under current IUCN criteria, P. gulbaru warrants an Endangered ( B1, 2) listing.

Bridging the gap between different statistical approaches: An integrated framework for modelling

This paper proposes a template for modelling complex datasets that integrates traditional statistical modelling approaches with more recent advances in statistics and modelling through an exploratory framework. Our approach builds on the well-known and long standing traditional idea of 'good practice in statistics' by establishing a comprehensive framework for modelling that focuses on exploration, prediction, interpretation and reliability assessment, a relatively new idea that allows individual assessment of predictions. The integrated framework we present comprises two stages. The first involves the use of exploratory methods to help visually understand the data and identify a parsimonious set of explanatory variables. The second encompasses a two step modelling process, where the use of non-parametric methods such as decision trees and generalized additive models are promoted to identify important variables and their modelling relationship with the response before a final predictive model is considered. We focus on fitting the predictive model using parametric, non-parametric and Bayesian approaches. This paper is motivated by a medical problem where interest focuses on developing a risk stratification system for morbidity of 1,710 cardiac patients given a suite of demographic, clinical and preoperative variables. Although the methods we use are applied specifically to this case study, these methods can be applied across any field, irrespective of the type of response.

Phylogeny and classification of the Digenea (Platyhelminthes : Trematoda)

Complete small subunit ribosomal RNA gene (ssrDNA) and partial (D1-D3) large subunit ribosomal RNA gene (lsrDNA) sequences were used to estimate the phylogeny of the Digenea via maximum parsimony and Bayesian inference. Here we contribute 80 new ssrDNA and 124 new lsrDNA sequences. Fully complementary data sets of the two genes were assembled from newly generated and previously published sequences and comprised 163 digenean taxa representing 77 nominal families and seven aspidogastrean outgroup taxa representing three families. Analyses were conducted on the genes independently as well as combined and separate analyses including only the higher plagiorchiidan taxa were performed using a reduced-taxon alignment including additional characters that could not be otherwise unambiguously aligned. The combined data analyses yielded the most strongly supported results and differences between the two methods of analysis were primarily in their degree of resolution. The Bayesian analysis including all taxa and characters, and incorporating a model of nucleotide substitution (general-time-reversible with among-site rate heterogeneity), was considered the best estimate of the phylogeny and was used to evaluate their classification and evolution. In broad terms, the Digenea forms a dichotomy that is split between a lineage leading to the Brachylaimoidea, Diplostomoidea and Schistosomatoidea (collectively the Diplostomida nomen novum (nom. nov.)) and the remainder of the Digenea (the Plagiorchiida), in which the Bivesiculata nom. nov. and Transversotremata nom. nov. form the two most basal lineages, followed by the Hemiurata. The remainder of the Plagiorchiida forms a large number of independent lineages leading to the crown clade Xiphidiata nom. nov. that comprises the Allocreadioidea, Gorgoderoidea, Microphalloidea and Plagiorchioidea, which are united by the presence of a penetrating stylet in their cercariae. Although a majority of families and to a lesser degree, superfamilies are supported as currently defined, the traditional divisions of the Echinostomida, Plagiorchiida and Strigeida were found to comprise non-natural assemblages. Therefore, the membership of established higher taxa are emended, new taxa erected and a revised, phylogenetically based classification proposed and discussed in light of ontogeny, morphology and taxonomic history. (C) 2003 Australian Society for Parasitology Inc. Published by Elsevier Science Ltd. All rights reserved.

A molecular phylogeny of the Australian skink genera Eulamprus, Gnypetoscincus and Nangura

Skinks from the genera Eulamprus, Gnypetoscincus and Nangura are a prominent component of the reptile fauna of the mesic forests of the east coast of Australia and have been the subject of numerous ecological studies. Highly conserved morphology and the retention of ancestral traits have limited our understanding of the relationships within and among these genera beyond an initial identification of species groups within Eulamprus. To address this deficit and to explore the relationships between Eulamprus and the monotypic genera Nangura and Gnypetoscincus, sections of two mitochondrial genes (ND4 and 16S rRNA) were sequenced and subjected to Bayesian phylogenetic analysis. This phylogenetic analysis supports recognition of the three species groups proposed for Eulamprus (murrayi, quoyii and tenuis) and indicates that this genus is paraphyletic, with Gnypetoscincus and Nangura being proximal to basal lineages of the tenuis group. To resolve these and broader problems of paraphyly, we suggest that each of the species groups from 'Eulamprus' should be recognised as a distinct genus. The phylogenetically and ecologically distinct water skinks of the quoyii group would be retained within Eulamprus and the diverse species of the tenuis group allocated to Concinnia. We suggest placing the monophyletic murrayi group, endemic to the rainforests of central eastern Australia, in a new genus ( yet to be formally described). The sequencing data also revealed the existence of a genetically divergent but morphologically cryptic lineage within E. murrayi and substantial diversity within E. quoyii. There is evidence for two major habitat shifts from rainforest towards drier habitats, one leading to the quoyii group and the second defining a clade of three species within the tenuis complex. These ecological transitions may represent adaptations to general drying across eastern Australia during the late Miocene - Pliocene. Each of the major areas of east coast tropical or subtropical rainforest contains multiple phylogenetically diverse endemic species, reflecting the long-term persistence and high conservation value of wet forest habitats in each area.

Gender differences in CF phenotype associated with low bone mineral density

Adolescents and adults with CF have lower bone mineral density (BMD) than normal, but its relationship with phenotype is not well understood. Point FEV1% predicted (FEV) and rate of change of FEV are biased estimates of disease severity, because progressively older subjects represent a selected survivor population, with females at greater risk of death than males. To investigate the relationship between BMD and phenotype we used an index (predicted age at death) derived from Bayesian estimates of slope and intercept of FEV, age at last measurement and survival status. Predictive equations for the index were derived from 97 subjects (78 survivors) from the RCH CF clinic, and applied to a group of 102 comparable subjects who had BMD measured, classified as having‘mild’ ()75th), ‘moderate’ (25– 75th), or ‘severe’ (-25th centile) phenotype. Total body (TB) and lumbar spine (LS) BMD z-scores (Z) were compared, adjustingfor gender effects, using 2-way ANOVA. Annual mean change in FEV segregated, as expected, according to phenotype, ‘severe’ (ns25), ‘moderate’ (ns51) and ‘mild’ (ns25) y3.01(y3.73 to y2.30)%, y0.85(y1.36 to y0.35)%, 2.70(1.92 to 3.46)%, respectively, with no gender difference. LS and TB BMDZ were different in each phenotype (P-s 0.002), LS BMDZ for ‘severe’, ‘moderate’ and ‘mild’ y1.63(CI: y2.07 to y 1.19), y0.86(CI: y1.17 to y0.55), y0.06(CI: y0.54 to 0.41). Males had lower LS BMDZ than females overall (y1.22 (CI: y1.54 to y0.91) vs. y0.48(CI: y 0.84 to y0.12) Ps0.002). In the ‘severe’ group, males had lower TB BMDZ and LS BMDZ (PF0.002). Low BMD is associated with ‘moderate’ and ‘severe’ phenotypes, with relative preservation in females in the ‘severe’ group. Female biology (reproductive fitness) might promote resistance to bone resorption at a critical level of BMD loss.

Padrões espaço-temporais da mortalidade por câncer de mama em mulheres no Estado do Espírito Santo - um estudo ecológico

O câncer de mama é a principal neoplasia maligna que acomete o sexo feminino no Brasil. O câncer de mama é hoje uma doença de extrema importância para a saúde pública nacional, motivando ampla discussão em torno das medidas que promova o seu diagnóstico precoce, a redução em sua morbidade e mortalidade. A presente pesquisa possui três objetivos, cujos resultados encontram-se organizados em artigos. O primeiro objetivo buscou analisar a completude dos dados do Sistema de Informação de Mortalidade sobre os óbitos por câncer de mama em mulheres no Espírito Santo, Sudeste e Brasil (1998 a 2007). Realizou-se um estudo descritivo analítico baseado em dados secundários, onde foi analisado o número absoluto e percentual de não preenchimento das variáveis nas declarações de óbitos. Adotou-se escore para avaliar os graus de não completude. Os resultados para as variáveis sexo e idade foram excelentes tanto para o Espírito Santo, Sudeste e Brasil. O preenchimento das variáveis raça/cor, grau de escolaridade e estado civil apresentam problemas no Espírito Santo. Enquanto no Sudeste e Brasil as variáveis raça/cor e escolaridade têm tendência decrescente para a não completude, no Espírito Santo a tendência se mantém estável. Para a variável estado civil, a não completude tem tendência crescente no Estado do Espírito Santo. O segundo objetivo foi analisar a evolução das taxas de mortalidade por câncer de mama, em mulheres no Espírito Santo no período de 1980 a 2007. Estudo de série temporal, cujos dados sobre óbitos foram obtidos do Sistema de Informação de Mortalidade e as estimativas populacionais segundo idade e anos-calendário, do Instituto Brasileiro Geografia e Estatística. Os coeficientes específicos 9 de mortalidade, segundo faixa etária, foram calculados anualmente. A análise de tendência foi realizada por meio da padronização das taxas de mortalidade pelo método direto, em que a população do senso IBGE-2000, foi considerada padrão. No período de estudo, ocorreram 2.736 óbitos por câncer de mama. O coeficiente de mortalidade neste período variou de 3,41 a 10,99 por 100.000 mulheres. Os resultados indicam que há tendência de mortalidade por câncer de mama ao longo da série (p=0,001 com crescimento de 75,42%). Todas as faixas etárias a partir de 30 anos apresentaram tendência de crescimento da mortalidade estatisticamente significante (p=0,001). Os percentuais de crescimento foram aumentando, segundo as idades mais avançadas, sendo 48,4% na faixa de 40 a 49 anos, chegando a 92,3%, na faixa de 80 anos e mais. O terceiro objetivo foi realizar a análise espacial dos óbitos em mulheres por câncer de mama no estado do Espírito Santo, nos anos de 2003 a 2007, com análise das correlações espaciais dessa mortalidade e componentes do município. O cenário foi o Estado do Espírito Santo, composto por 78 municípios. Para análise dos dados, utilizou-se a abordagem bayesiana (métodos EBest Global e EBest Local) para correção de taxas epidemiológicas. Calculou-se o índice I de Moran, para dependência espacial em nível global e a estatística Moran Local. As maiores taxas estão concentradas em 19 municípios pertencentes às Microrregiões: Metropolitana (Fundão, Vitória, Vila Velha, Viana, Cariacica e Guarapari), Metrópole Expandida Sul (Anchieta, Alfredo Chaves), Pólo Cachoeiro (Vargem Alta, Rio Novo do Sul, Mimoso do Sul, Cachoeiro de Itapemirim, Castelo, Jerônimo Monteiro, Bom Jesus do Norte, Apiacá e Muqui) e Caparaó (Alegre e São José do Calçado). Os resultados da Estimação Bayesiana (Índice de Moran) dos óbitos por câncer de mama em mulheres ocorridos no estado do Espírito Santo, segundo os dados brutos e 10 ajustados indicam a existência de correlação espacial significativa para o mapa Local (I = 0,573; p = 0,001) e Global (I = 0,118; p = 0,039). Os dados brutos não apresentam correlação espacial (I = 0,075; p = 0,142).

Fluorescence Microscopy Imaging Denoising with Log-Euclidean Priors and Photobleaching Compensation

Fluorescent protein microscopy imaging is nowadays one of the most important tools in biomedical research. However, the resulting images present a low signal to noise ratio and a time intensity decay due to the photobleaching effect. This phenomenon is a consequence of the decreasing on the radiation emission efficiency of the tagging protein. This occurs because the fluorophore permanently loses its ability to fluoresce, due to photochemical reactions induced by the incident light. The Poisson multiplicative noise that corrupts these images, in addition with its quality degradation due to photobleaching, make long time biological observation processes very difficult. In this paper a denoising algorithm for Poisson data, where the photobleaching effect is explicitly taken into account, is described. The algorithm is designed in a Bayesian framework where the data fidelity term models the Poisson noise generation process as well as the exponential intensity decay caused by the photobleaching. The prior term is conceived with Gibbs priors and log-Euclidean potential functions, suitable to cope with the positivity constrained nature of the parameters to be estimated. Monte Carlo tests with synthetic data are presented to characterize the performance of the algorithm. One example with real data is included to illustrate its application.

Elements of a methodology to assess the alignment of core-values in collaborative networks

Collaborative networks are typically formed by heterogeneous and autonomous entities, and thus it is natural that each member has its own set of core-values. Since these values somehow drive the behaviour of the involved entities, the ability to quickly identify partners with compatible or common core-values represents an important element for the success of collaborative networks. However, tools to assess or measure the level of alignment of core-values are lacking. Since the concept of 'alignment' in this context is still ill-defined and shows a multifaceted nature, three perspectives are discussed. The first one uses a causal maps approach in order to capture, structure, and represent the influence relationships among core-values. This representation provides the basis to measure the alignment in terms of the structural similarity and influence among value systems. The second perspective considers the compatibility and incompatibility among core-values in order to define the alignment level. Under this perspective we propose a fuzzy inference system to estimate the alignment level, since this approach allows dealing with variables that are vaguely defined, and whose inter-relationships are difficult to define. Another advantage provided by this method is the possibility to incorporate expert human judgment in the definition of the alignment level. The last perspective uses a belief Bayesian network method, and was selected in order to assess the alignment level based on members' past behaviour. An example of application is presented where the details of each method are discussed.

Convex Total Variation Denoising of Poisson Fluorescence Confocal Images With Anisotropic Filtering

Fluorescence confocal microscopy (FCM) is now one of the most important tools in biomedicine research. In fact, it makes it possible to accurately study the dynamic processes occurring inside the cell and its nucleus by following the motion of fluorescent molecules over time. Due to the small amount of acquired radiation and the huge optical and electronics amplification, the FCM images are usually corrupted by a severe type of Poisson noise. This noise may be even more damaging when very low intensity incident radiation is used to avoid phototoxicity. In this paper, a Bayesian algorithm is proposed to remove the Poisson intensity dependent noise corrupting the FCM image sequences. The observations are organized in a 3-D tensor where each plane is one of the images acquired along the time of a cell nucleus using the fluorescence loss in photobleaching (FLIP) technique. The method removes simultaneously the noise by considering different spatial and temporal correlations. This is accomplished by using an anisotropic 3-D filter that may be separately tuned in space and in time dimensions. Tests using synthetic and real data are described and presented to illustrate the application of the algorithm. A comparison with several state-of-the-art algorithms is also presented.

Role of malnutrition and parasite infections in the spatial variation in children’s anaemia risk in Northern Angola

Anaemia has a significant impact on child development and mortality and is a severe public health problem in most countries in sub-Saharan Africa. Nutritional and infectious causes of anaemia are geographically variable and anaemia maps based on information on the major aetiologies of anaemia are important for identifying communities most in need and the relative contribution of major causes. We investigated the consistency between ecological and individual-level approaches to anaemia mapping, by building spatial anaemia models for children aged ≤15 years using different modeling approaches. We aimed to a) quantify the role of malnutrition, malaria, Schistosoma haematobium and soil-transmitted helminths (STH) for anaemia endemicity in children aged ≤15 years and b) develop a high resolution predictive risk map of anaemia for the municipality of Dande in Northern Angola. We used parasitological survey data on children aged ≤15 years to build Bayesian geostatistical models of malaria (PfPR≤15), S. haematobium, Ascaris lumbricoides and Trichuris trichiura and predict small-scale spatial variation in these infections. The predictions and their associated uncertainty were used as inputs for a model of anemia prevalence to predict small-scale spatial variation of anaemia. Stunting, PfPR≤15, and S. haematobium infections were significantly associated with anaemia risk. An estimated 12.5%, 15.6%, and 9.8%, of anaemia cases could be averted by treating malnutrition, malaria, S. haematobium, respectively. Spatial clusters of high risk of anaemia (>86%) were identified. Using an individual-level approach to anaemia mapping at a small spatial scale, we found that anaemia in children aged ≤15 years is highly heterogeneous and that malnutrition and parasitic infections are important contributors to the spatial variation in anemia risk. The results presented in this study can help inform the integration of the current provincial malaria control program with ancillary micronutrient supplementation and control of neglected tropical diseases, such as urogenital schistosomiasis and STH infection.

Role of malnutrition and parasite infections in the spatial variation in children’s anaemia risk in northern Angola

Anaemia is known to have an impact on child development and mortality and is a severe public health problem in most countries in sub-Saharan Africa. We investigated the consistency between ecological and individual-level approaches to anaemia mapping by building spatial anaemia models for children aged ≤15 years using different modelling approaches. We aimed to (i) quantify the role of malnutrition, malaria, Schistosoma haematobium and soil-transmitted helminths (STHs) in anaemia endemicity; and (ii) develop a high resolution predictive risk map of anaemia for the municipality of Dande in northern Angola. We used parasitological survey data for children aged ≤15 years to build Bayesian geostatistical models of malaria (PfPR≤15), S. haematobium, Ascaris lumbricoides and Trichuris trichiura and predict small-scale spatial variations in these infections. Malnutrition, PfPR≤15, and S. haematobium infections were significantly associated with anaemia risk. An estimated 12.5%, 15.6% and 9.8% of anaemia cases could be averted by treating malnutrition, malaria and S. haematobium, respectively. Spatial clusters of high risk of anaemia (>86%) were identified. Using an individual-level approach to anaemia mapping at a small spatial scale, we found that anaemia in children aged ≤15 years is highly heterogeneous and that malnutrition and parasitic infections are important contributors to the spatial variation in anaemia risk. The results presented in this study can help inform the integration of the current provincial malaria control programme with ancillary micronutrient supplementation and control of neglected tropical diseases such as urogenital schistosomiasis and STH infections.

Diffuse liver disease classification from ultrasound surface characterization, clinical and laboratorial data

In this work liver contour is semi-automatically segmented and quantified in order to help the identification and diagnosis of diffuse liver disease. The features extracted from the liver contour are jointly used with clinical and laboratorial data in the staging process. The classification results of a support vector machine, a Bayesian and a k-nearest neighbor classifier are compared. A population of 88 patients at five different stages of diffuse liver disease and a leave-one-out cross-validation strategy are used in the classification process. The best results are obtained using the k-nearest neighbor classifier, with an overall accuracy of 80.68%. The good performance of the proposed method shows a reliable indicator that can improve the information in the staging of diffuse liver disease.

Fatty liver characterization and classification by ultrasound

Steatosis, also known as fatty liver, corresponds to an abnormal retention of lipids within the hepatic cells and reflects an impairment of the normal processes of synthesis and elimination of fat. Several causes may lead to this condition, namely obesity, diabetes, or alcoholism. In this paper an automatic classification algorithm is proposed for the diagnosis of the liver steatosis from ultrasound images. The features are selected in order to catch the same characteristics used by the physicians in the diagnosis of the disease based on visual inspection of the ultrasound images. The algorithm, designed in a Bayesian framework, computes two images: i) a despeckled one, containing the anatomic and echogenic information of the liver, and ii) an image containing only the speckle used to compute the textural features. These images are computed from the estimated RF signal generated by the ultrasound probe where the dynamic range compression performed by the equipment is taken into account. A Bayes classifier, trained with data manually classified by expert clinicians and used as ground truth, reaches an overall accuracy of 95% and a 100% of sensitivity. The main novelties of the method are the estimations of the RF and speckle images which make it possible to accurately compute textural features of the liver parenchyma relevant for the diagnosis.