CAMCR: Computer assisted mixture model analysis for capture-recapture count data

Population size estimation with discrete or nonparametric mixture models is considered, and reliable ways of construction of the nonparametric mixture model estimator are reviewed and set into perspective. Construction of the maximum likelihood estimator of the mixing distribution is done for any number of components up to the global nonparametric maximum likelihood bound using the EM algorithm. In addition, the estimators of Chao and Zelterman are considered with some generalisations of Zelterman’s estimator. All computations are done with CAMCR, a special software developed for population size estimation with mixture models. Several examples and data sets are discussed and the estimators illustrated. Problems using the mixture model-based estimators are highlighted.

A capture-recapture approach for screening using two diagnostic tests with availability of disease status for the test positives only

The article considers screening human populations with two screening tests. If any of the two tests is positive, then full evaluation of the disease status is undertaken; however, if both diagnostic tests are negative, then disease status remains unknown. This procedure leads to a data constellation in which, for each disease status, the 2 x 2 table associated with the two diagnostic tests used in screening has exactly one empty, unknown cell. To estimate the unobserved cell counts, previous approaches assume independence of the two diagnostic tests and use specific models, including the special mixture model of Walter or unconstrained capture-recapture estimates. Often, as is also demonstrated in this article by means of a simple test, the independence of the two screening tests is not supported by the data. Two new estimators are suggested that allow associations of the screening test, although the form of association must be assumed to be homogeneous over disease status. These estimators are modifications of the simple capture-recapture estimator and easy to construct. The estimators are investigated for several screening studies with fully evaluated disease status in which the superior behavior of the new estimators compared to the previous conventional ones can be shown. Finally, the performance of the new estimators is compared with maximum likelihood estimators, which are more difficult to obtain in these models. The results indicate the loss of efficiency as minor.

Gene flow, risk assessment and the environmental release of transgenic plants

The release of genetically modified plants is governed by regulations that aim to provide an assessment of potential impact on the environment. One of the most important components of this risk assessment is an evaluation of the probability of gene flow. In this review, we provide an overview of the current literature on gene flow from transgenic plants, providing a framework of issues for those considering the release of a transgenic plant into the environment. For some plants gene flow from transgenic crops is well documented, and this information is discussed in detail in this review. Mechanisms of gene flow vary from plant species to plant species and range from the possibility of asexual propagation, short- or long-distance pollen dispersal mediated by insects or wind and seed dispersal. Volunteer populations of transgenic plants may occur where seed is inadvertently spread during harvest or commercial distribution. If there are wild populations related to the transgenic crop then hybridization and eventually introgression in the wild may occur, as it has for herbicide resistant transgenic oilseed rape (Brassica napus). Tools to measure the amount of gene flow, experimental data measuring the distance of pollen dispersal, and experiments measuring hybridization and seed survivability are discussed in this review. The various methods that have been proposed to prevent gene flow from genetically modified plants are also described. The current "transgenic traits'! in the major crops confer resistance to herbicides and certain insects. Such traits could confer a selective advantage (an increase in fitness) in wild plant populations in some circumstances, were gene flow to occur. However, there is ample evidence that gene flow from crops to related wild species occurred before the development of transgenic crops and this should be taken into account in the risk assessment process.

Application of one-list capture-recapture models to scrapie surveillance data in Great Britain

In this paper, we apply one-list capture-recapture models to estimate the number of scrapie-affected holdings in Great Britain. We applied this technique to the Compulsory Scrapie Flocks Scheme dataset where cases from all the surveillance sources monitoring the presence of scrapie in Great Britain, the abattoir survey, the fallen stock survey and the statutory reporting of clinical cases, are gathered. Consequently, the estimates of prevalence obtained from this scheme should be comprehensive and cover all the different presentations of the disease captured individually by the surveillance sources. Two estimators were applied under the one-list approach: the Zelterman estimator and Chao's lower bound estimator. Our results could only inform with confidence the scrapie-affected holding population with clinical disease; this moved around the figure of 350 holdings in Great Britain for the period under study, April 2005-April 2006. Our models allowed the stratification by surveillance source and the input of covariate information, holding size and country of origin. None of the covariates appear to inform the model significantly. Crown Copyright (C) 2008 Published by Elsevier B.V. All rights reserved.

Comparison of survival rates of captive-reared and wild-bred Mauritius kestrels (falco punctatus) in a re-introduced population

Re-introduction is a technique widely used in the conservation of threatened bird species. With advances in aviculture the use of captive-produced individuals as the release stock is becoming more commonplace, and ideally, survival of captive-produced, released individuals should be no different from their wild-bred counterparts. During the late 1980s the Critically Endangered Mauritius kestrel (Falco punctatus) was successfully re-introduced into the Bambous mountain range, Mauritius, some 30 years after its local extinction. Between 1987 and 2001 the developing population was closely monitored enabling us to construct re-sighting histories for 88 released and 284 wild-bred kestrels. We used age-structured models in the survival analysis software program MARK to determine if an individual's origin influenced its subsequent survival. Our analysis indicated no compelling evidence for reduced survival among juvenile captive-reared and released individuals, relative to their wild-bred counterparts, across the majority of cohorts and only limited evidence of a cohort-specific effect. This study illustrates that despite the lack of a formal experimental approach it is still feasible to conduct an assessment of re-introduction outcomes and techniques. (C) 2003 Elsevier Ltd. All rights reserved.

Declining survival rates in a reintroduced population of the Mauritius kestrel: evidence for non-linear density dependence and environmental stochasticity

1. We studied a reintroduced population of the formerly critically endangered Mauritius kestrel Falco punctatus Temmink from its inception in 1987 until 2002, by which time the population had attained carrying capacity for the study area. Post-1994 the population received minimal management other than the provision of nestboxes. 2. We analysed data collected on survival (1987-2002) using program MARK to explore the influence of density-dependent and independent processes on survival over the course of the population's development. 3.We found evidence for non-linear, threshold density dependence in juvenile survival rates. Juvenile survival was also strongly influenced by climate, with the temporal distribution of rainfall during the cyclone season being the most influential climatic variable. Adult survival remained constant throughout. 4. Our most parsimonious capture-mark-recapture statistical model, which was constrained by density and climate, explained 75.4% of the temporal variation exhibited in juvenile survival rates over the course of the population's development. 5. This study is an example of how data collected as part of a threatened species recovery programme can be used to explore the role and functional form of natural population regulatory processes. With the improvements in conservation management techniques and the resulting success stories, formerly threatened species offer unique opportunities to further our understanding of the fundamental principles of population ecology.

Fragrance release from the surface of branched poly(amide)s

Enzymes are powerful tools in organic synthesis that are able to catalyse a wide variety of selective chemical transformations under mild and environmentally friendly conditions. Enzymes such as the lipases have also found applications in the synthesis and degradation of polymeric materials. However, the use of these natural catalysts in the synthesis and the post-synthetic modification of dendrimers and hyperbranched molecules is an application of chemistry yet to be explored extensively. In this study the use of two hydrolytic enzymes, a lipase from Candida cylindracea and a cutinase from Fusarium solani pisii, were investigated in the selective cleavage of ester groups situated on the peripheral layer of two families of branched polyamides. These branched polyamides were conjugated to simple fragrances citronellol and L-menthol via ester linkages. Hydrolysis of the ester linkage between the fragrances and the branched polyamide support was carried out in aqueous buffered systems at slightly basic pH values under the optimum operative conditions for the enzymes used. These preliminary qualitative investigations revealed that partial cleavage of the ester functionalities from the branched polyamide support had occurred. However, the ability of the enzymes to interact with the substrates decreased considerably as the branching density, the rigidity of the structure and the bulkiness of the polyamide-fragrance conjugates increased.

Thermoselective phase property of silica tethered with fluorinated chains for controlled 'release' and 'capture' of catalytic fluorous tin species

The temperature dependent mixing of organic and fluorous phases is one of the key principals of fluorous biphasic systems (FBS). Given the high cost of the perfluorous solvents and their impacts to the environment, it is apparent that elimination of these solvents in bulk quantity in the FBS is advantageous. We report for the first time, the surface coverage of silica with a fluorous solvent like material that traps (at ambient temperatures) and releases (at elevated temperatures) a fluorous tin bromide in organic solvent. Here, we demonstrate the catalytic utilisation of this species for the hydrocyclisation of 6-bromo-1-hexene with NaBH4. (C) 2002 Elsevier Science B.V. All rights reserved.

G protein ss gamma subunits mediate presynaptic inhibition of transmitter release from rat superior cervical ganglion neurones in culture

The activation of presynaptic G protein-coupled receptors (GPCRs) is widely reported to inhibit transmitter release; however, the lack of accessibility of many presynaptic terminals has limited direct analysis of signalling mediators. We studied GPCR-mediated inhibition of fast cholinergic transmission between superior cervical ganglion neurones (SCGNs) in culture. The adrenoceptor agonist noradrenaline (NA) caused a dose-related reduction in evoked excitatory postsynaptic potentials (EPSPs). NA-induced EPSP decrease was accompanied by effects on the presynaptic action potential (AP), reducing AP duration and amplitude of the after-hyperpolarization (AHP), without affecting the pre- and postsynaptic membrane potential. All effects of NA were blocked by yohimbine and synaptic transmission was reduced by clonidine, consistent with an action at presynaptic alpha 2-adrenoceptors. NA-induced inhibition of transmission was sensitive to pre-incubation of SCGNs with pertussis toxin (PTX), implicating the involvement of G alpha(i)/(o)beta y subunits. Expression of G alpha transducin, an agent which sequesters G protein beta gamma (G beta y) subunits, in the presynaptic neurone caused a time-dependent attenuation of NA-induced inhibition. Injection of purified G beta gamma subunits into the presynaptic neurone inhibited transmission, and also reduced the AHP amplitude. Furthermore, NA-induced inhibition was occluded by pre-injection of G beta gamma subunits. The Ca2+ channel blocker Cd2+ mimicked NA effects on transmitter release. Cd2+, NA and G beta gamma subunits also inhibited somatic Ca2+ current. In contrast to effects on AP-evoked transmitter release, NA had no clear action on AP-independent EPSPs induced by hypertonic solutions. These results demonstrate that G beta gamma subunits functionally mediate inhibition of transmitter release by alpha 2-adrenoceptors and represent important regulators of synaptic transmission at mammalian presynaptic terminals.

The relationship between membrane damage, release of protein and loss of viability in Escherichia coli exposed to high hydrostatic pressure

The aim of this work was to examine a possible association between resistance of two Escherichia coli strains to high hydrostatic pressure and the susceptibility of their cell membranes to pressure-induced damage. Cells were exposed to pressures between 100 and 700 MPa at room temperature (~20C) in phosphate-buffered-saline. In the more pressure-sensitive strain E. coli 8164, loss of viability occurred at pressures between 100 MPa and 300 MPa and coincided with irreversible loss of membrane integrity as indicated by uptake of propidium iodide (PI) and leakage of protein of molecular mass between 9 and 78 kDa from the cells. Protein release increased to a maximum at 400 MPa then decreased, possibly due to intracellular aggregation at the higher pressures. In the pressure-resistant strain E. coli J1, PI was taken up during pressure treatment but not after decompression indicating that cells were able to reseal their membranes. Loss of viability in strain J1 coincided with the transient loss of membrane integrity between approximately 200 MPa and 600 MPa. In E. coli J1 leakage of protein occurred before loss of viability and the released protein was of low molecular mass, between 8 and 11 kDa and may have been of periplasmic origin. In these two strains differences in pressure resistance appeared to be related to differences in the ability of their membranes to withstand disruption by pressure. However it appears that transient loss of membrane integrity during pressure can lead to cell death irrespective of whether cells can reseal their membranes afterwards.

PEGylated amyloid peptide nanocontainer delivery and release system

A micellar nanocontainer delivery and release system is designed on the basis of a peptide-polymer conjugate. The hybrid molecules self-assemble into micelles comprising a modified amyloid peptide core surrounded by a PEG corona. The modified amyloid peptide previously studied in our group forms helical ribbons based on a beta-sheet motif and contains beta-amino acids that are excluded from the beta-sheet structure, thus being potentially useful as fibrillization inhibitors. In the model peptide-PEG hybrid system studied, enzymatic degradation using alpha-chymotrypsin leads to selective cleavage close to the PEG-peptide linkage, break up of the micelles, and release of peptides in unassociated form. The release of monomeric peptide is useful because aggregation of the released peptide into beta-sheet amyloid fibrils is not observed. This concept has considerable potential in the targeted delivery of peptides for therapeutic applications.

Multicore-enabling the MPJ express messaging library

With the transition to multicore processors almost complete, the parallel processing community is seeking efficient ways to port legacy message passing applications on shared memory and multicore processors. MPJ Express is our reference implementation of Message Passing Interface (MPI)-like bindings for the Java language. Starting with the current release, the MPJ Express software can be configured in two modes: the multicore and the cluster mode. In the multicore mode, parallel Java applications execute on shared memory or multicore processors. In the cluster mode, Java applications parallelized using MPJ Express can be executed on distributed memory platforms like compute clusters and clouds. The multicore device has been implemented using Java threads in order to satisfy two main design goals of portability and performance. We also discuss the challenges of integrating the multicore device in the MPJ Express software. This turned out to be a challenging task because the parallel application executes in a single JVM in the multicore mode. On the contrary in the cluster mode, the parallel user application executes in multiple JVMs. Due to these inherent architectural differences between the two modes, the MPJ Express runtime is modified to ensure correct semantics of the parallel program. Towards the end, we compare performance of MPJ Express (multicore mode) with other C and Java message passing libraries---including mpiJava, MPJ/Ibis, MPICH2, MPJ Express (cluster mode)---on shared memory and multicore processors. We found out that MPJ Express performs signicantly better in the multicore mode than in the cluster mode. Not only this but the MPJ Express software also performs better in comparison to other Java messaging libraries including mpiJava and MPJ/Ibis when used in the multicore mode on shared memory or multicore processors. We also demonstrate effectiveness of the MPJ Express multicore device in Gadget-2, which is a massively parallel astrophysics N-body siimulation code.

Constraint release in entangled binary blends of linear polymers: a molecular dynamics study

We present extensive molecular dynamics simulations of the dynamics of diluted long probe chains entangled with a matrix of shorter chains. The chain lengths of both components are above the entanglement strand length, and the ratio of their lengths is varied over a wide range to cover the crossover from the chain reptation regime to tube Rouse motion regime of the long probe chains. Reducing the matrix chain length results in a faster decay of the dynamic structure factor of the probe chains, in good agreement with recent neutron spin echo experiments. The diffusion of the long chains, measured by the mean square displacements of the monomers and the centers of mass of the chains, demonstrates a systematic speed-up relative to the pure reptation behavior expected for monodisperse melts of sufficiently long polymers. On the other hand, the diffusion of the matrix chains is only weakly perturbed by the diluted long probe chains. The simulation results are qualitatively consistent with the theoretical predictions based on constraint release Rouse model, but a detailed comparison reveals the existence of a broad distribution of the disentanglement rates, which is partly confirmed by an analysis of the packing and diffusion of the matrix chains in the tube region of the probe chains. A coarse-grained simulation model based on the tube Rouse motion model with incorporation of the probability distribution of the tube segment jump rates is developed and shows results qualitatively consistent with the fine scale molecular dynamics simulations. However, we observe a breakdown in the tube Rouse model when the short chain length is decreased to around N-S = 80, which is roughly 3.5 times the entanglement spacing N-e(P) = 23. The location of this transition may be sensitive to the chain bending potential used in our simulations.