A general hazard model for lifetime data in the presence of cure rate

Historically, the cure rate model has been used for modeling time-to-event data within which a significant proportion of patients are assumed to be cured of illnesses, including breast cancer, non-Hodgkin lymphoma, leukemia, prostate cancer, melanoma, and head and neck cancer. Perhaps the most popular type of cure rate model is the mixture model introduced by Berkson and Gage [1]. In this model, it is assumed that a certain proportion of the patients are cured, in the sense that they do not present the event of interest during a long period of time and can found to be immune to the cause of failure under study. In this paper, we propose a general hazard model which accommodates comprehensive families of cure rate models as particular cases, including the model proposed by Berkson and Gage. The maximum-likelihood-estimation procedure is discussed. A simulation study analyzes the coverage probabilities of the asymptotic confidence intervals for the parameters. A real data set on children exposed to HIV by vertical transmission illustrates the methodology.

Prevalence of trachoma in a population of the upper Rio Negro Basin and risk factors for active disease

Purpose: To determine the prevalence of trachoma in Sao Gabriel da Cachoeira (SGC), the only urban community of the upper Rio Negro Basin of the Amazon state in Brazil, near the Colombian border, and to investigate the risk factors associated with the active forms of the disease. Methods: A total of 1702 people (440 children up to 9 years and 1069 adults aged 15 years and above) were examined. The sample was selected from a probabilistic household sampling procedure based on census data and a previous study of trachoma prevalence in Sao Gabriel da Cachoeira. A two-stage probabilistic household cluster sample was drawn. Household units were randomly selected within each cluster. A variety of socioeconomic and hygiene variables were studied in order to determine the risk factors for active trachoma in a household. Results: The total prevalence of trachoma was 8.9%. Prevalence of active trachoma (TF and/or TI) in children aged 1-9 years was 11.1% and trachomatous trichiasis in adults aged 15 years and above was 0.19%. Trachomatous scarring reached a peak of 22.4% for subjects between 50 to 60 years of age. Corneal opacity occurred in subjects aged 50 years and older with a prevalence of 2.0%. No sex effect was found on the overall prevalence of trachoma in SGC. Risk factors associated with active trachoma were mainly related to poor socioeconomic indicators. Conclusions: Despite the ubiquitous presence of water, the analysis of the risk factors associated with the active forms of the disease supports the idea that a low personal standard of hygiene and not water availability per se, is the key factor associated with trachoma.

Detection of Auditory Cortex Activity by fMRI Using a Dependent Component Analysis

Functional MRI (fMRI) data often have low signal-to-noise-ratio (SNR) and are contaminated by strong interference from other physiological sources. A promising tool for extracting signals, even under low SNR conditions, is blind source separation (BSS), or independent component analysis (ICA). BSS is based on the assumption that the detected signals are a mixture of a number of independent source signals that are linearly combined via an unknown mixing matrix. BSS seeks to determine the mixing matrix to recover the source signals based on principles of statistical independence. In most cases, extraction of all sources is unnecessary; instead, a priori information can be applied to extract only the signal of interest. Herein we propose an algorithm based on a variation of ICA, called Dependent Component Analysis (DCA), where the signal of interest is extracted using a time delay obtained from an autocorrelation analysis. We applied such method to inspect functional Magnetic Resonance Imaging (fMRI) data, aiming to find the hemodynamic response that follows neuronal activation from an auditory stimulation, in human subjects. The method localized a significant signal modulation in cortical regions corresponding to the primary auditory cortex. The results obtained by DCA were also compared to those of the General Linear Model (GLM), which is the most widely used method to analyze fMRI datasets.

Development and agronomic performance of common bean lines simultaneously resistant to anthracnose, angular leaf spot and rust

The common bean is affected by several pathogens that can cause severe yield losses. Here we report the introgression of resistance genes to anthracnose, angular leaf spot and rust in the `carioca-type` bean cultivar `Ruda`. Initially, four backcross (BC) lines were obtained using `TO`, `AB 136`, `Ouro Negro` and `AND 277` as donor parents. Molecular fingerprinting was used to select the lines genetically closer to the recurrent parent. The relative genetic distances between `Ruda` and the BC lines varied between 0.0% and 1.99%. The BC lines were intercrossed and molecular markers linked to the resistance genes were used to identify the plants containing the genes of interest. These plants were selfed to obtain the F(2), F(3) and F(4) plants which were selected based on the presence of the molecular markers mentioned and resistance was confirmed in the F(4) generation by inoculation. Four F(4:7) pyramid lines with all the resistance genes showed resistance spectra equivalent to those of their respective donor parents. Yield tests showed that these lines are as productive as the best `carioca-type` cultivars.

A Comparison of Bayesian and Sampling Theory Inferences in a Probit Model

HE PROBIT MODEL IS A POPULAR DEVICE for explaining binary choice decisions in econometrics. It has been used to describe choices such as labor force participation, travel mode, home ownership, and type of education. These and many more examples can be found in papers by Amemiya (1981) and Maddala (1983). Given the contribution of economics towards explaining such choices, and given the nature of data that are collected, prior information on the relationship between a choice probability and several explanatory variables frequently exists. Bayesian inference is a convenient vehicle for including such prior information. Given the increasing popularity of Bayesian inference it is useful to ask whether inferences from a probit model are sensitive to a choice between Bayesian and sampling theory techniques. Of interest is the sensitivity of inference on coefficients, probabilities, and elasticities. We consider these issues in a model designed to explain choice between fixed and variable interest rate mortgages. Two Bayesian priors are employed: a uniform prior on the coefficients, designed to be noninformative for the coefficients, and an inequality restricted prior on the signs of the coefficients. We often know, a priori, whether increasing the value of a particular explanatory variable will have a positive or negative effect on a choice probability. This knowledge can be captured by using a prior probability density function (pdf) that is truncated to be positive or negative. Thus, three sets of results are compared:those from maximum likelihood (ML) estimation, those from Bayesian estimation with an unrestricted uniform prior on the coefficients, and those from Bayesian estimation with a uniform prior truncated to accommodate inequality restrictions on the coefficients.

Evaluation of short dietary questions from the 1995 National Nutrition Survey

Six of the short dietary questions used in the 1995 National Nutrition Survey (see box below) were evaluated for relative validity both directly and indirectly and for consistency, by documenting the differences in mean intakes of foods and nutrients as measured on the 24-hour recall, between groups with different responses to the short questions. 1. Including snacks, how many times do you usually have something to eat in a day including evenings? 2. How many days per week do you usually have something to eat for breakfast? 3. In the last 12 months, were there any times that you ran out of food and couldn’t afford to buy more? 4. What type of milk do you usually consume? 5. How many serves of vegetables do you usually eat each day? (a serve = 1/2 cup cooked vegetables or 1 cup of salad vegetables) 6. How many serves of fruit do you usually eat each day? (a serve = 1 medium piece or 2 small pieces of fruit or 1 cup of diced pieces) These comparisons were made for males and females overall and for population sub-groups of interest including: age, socio-economic disadvantage, region of residence, country of birth, and BMI category. Several limitations to this evaluation of the short questions, as discussed in the report, need to be kept in mind including: · The method for comparison available (24-hour recall) was not ideal (gold standard); as it measures yesterday’s intake. This limitation was overcome by examining only mean differences between groups of respondents, since mean intake for a group can provide a reasonable approximation for ‘usual’ intake. · The need to define and identify, post-hoc, from the 24-hour recall the number of eating occasions, and occasions identified by the respondents as breakfast. · Predetermined response categories for some of the questions effectively limited the number of categories available for evaluation. · Other foods and nutrients, not selected for this evaluation, may have an indirect relationship with the question, and might have shown stronger and more consistent responses. · The number of responses in some categories of the short questions eg for food security may have been too small to detect significant differences between population sub-groups. · No information was available to examine the validity of these questions for detecting differences over time (establishing trends) in food habits and indicators of selected nutrient intakes. By contrast, the strength of this evaluation was its very large sample size, (atypical of most validation studies of dietary assessment) and thus, the opportunity to investigate question performance in a range of broad population sub-groups compared with a well-conducted, quantified survey of intakes. The results of the evaluation are summarised below for each of the questions and specific recommendations for future testing, modifications and use provided for each question. The report concludes with some general recommendations for the further development and evaluation of short dietary questions.

Lanczos subspace filter diagonalization: Homogeneous recursive filtering and a low-storage method for the calculation of matrix elements

We develop a new iterative filter diagonalization (FD) scheme based on Lanczos subspaces and demonstrate its application to the calculation of bound-state and resonance eigenvalues. The new scheme combines the Lanczos three-term vector recursion for the generation of a tridiagonal representation of the Hamiltonian with a three-term scalar recursion to generate filtered states within the Lanczos representation. Eigenstates in the energy windows of interest can then be obtained by solving a small generalized eigenvalue problem in the subspace spanned by the filtered states. The scalar filtering recursion is based on the homogeneous eigenvalue equation of the tridiagonal representation of the Hamiltonian, and is simpler and more efficient than our previous quasi-minimum-residual filter diagonalization (QMRFD) scheme (H. G. Yu and S. C. Smith, Chem. Phys. Lett., 1998, 283, 69), which was based on solving for the action of the Green operator via an inhomogeneous equation. A low-storage method for the construction of Hamiltonian and overlap matrix elements in the filtered-basis representation is devised, in which contributions to the matrix elements are computed simultaneously as the recursion proceeds, allowing coefficients of the filtered states to be discarded once their contribution has been evaluated. Application to the HO2 system shows that the new scheme is highly efficient and can generate eigenvalues with the same numerical accuracy as the basic Lanczos algorithm.

Quantitative PCR-ELAHA for the determination of retroviral vector transduction efficiency

Current methods to detect transduction efficiency during the routine use of integrating retroviral vectors in gene therapy applications may require the use of radioactivity and usually rely upon subjective determination of the results. We have developed two competitive quantitative assays that use an enzyme-linked, amplicon hybridization assay (ELAHA) to detect the products of PCR-amplified regions of transgene from cells transduced with Moloney murine leukemia virus vectors. The quantitative assays (PCR-ELAHA) proved to be simple, rapid, and sensitive, avoiding the need for Southern hybridization, complex histochemical stains, or often subjective and time-consuming tissue culture and immunofluorescence assays. The PCR-ELAHA systems can rapidly detect proviral DNA from any retroviral vector carrying the common selective and marker genes neomycin phosphotransferase and green fluorescent protein, and the methods described are equally applicable to other sequences of interest, providing a cheaper alternative to the evolving real-time PCR methods. The results revealed the number of copies of retrovector provirus present per stably transduced cell using vectors containing either one or both qPCR targets.

Development of disabled, replication-defective gene transfer vectors from the Jembrana disease virus, a new infectious agent of cattle

Jembrana disease virus (JDV) is a newly isolated and characterised bovine lentivirus. It causes an acute disease in Ball cattle (Bos javanicus). which can be readily transmitted to susceptible cattle with 17% mortality. There is as yet no treatment or preventive vaccine. We have developed a gene transfer vector system based on JDV that has three components. The first of the components is a bicistronic transfer vector plasmid that was constructed to contain cis-sequences from the JDV genome, including 5 '- and 3 ' -long terminal repeats (LTRs), 0.4 kb of truncated gag and 1.1 kb of 3 ' -env, a multiple cloning site to accommodate the gene(s) of interest for transfer, and an internal ribosome entry site plus the neomycin phosphotransferase (Neo) gene cassette for antibiotic selection. The second element is a packaging plasmid that contains trans-sequences. including gag, pol. vif, tar and rev: but without the env and packaging signals. The third is a plasmid encoding the G glycoprotein of vesicular stomatitis virus (VSV-G) to supply the vector an envelope for pseudotyping. Cotransfection of 293T cells with these three plasmid components produced VSV-G pseudotyped. disabled, replication defective, bicistronic JDV vectors encoding the green fluorescent protein (EGFP) and the Neo resistance selection maker simultaneously with a titre range of (0.4-1.2) x 10(6) CFU/ml. Transduction of several replicating primary and transformed cells from cattle, primate and human sources and importantly growth-arrested cells with the JDV vectors showed high efficiency of EGFP gene transfer at 35-75%, which was stable and the expression of EGFP was long term. Furthermore, these JDV vectors were designed to suit the inclusion and expression of genes corresponding to JDV specific proteins, such as gag or env, for the development of vaccines for Jembrana disease. This strategy should also be applicable to other bovine diseases as wall. The design and construction of the JDV vector system should facilitate the study of the lentivirology and pathogenesis of the diseases associated with JDV or other bovine virus infections. To our knowledge, this is the first such vector system developed from a cattle virus. (C) 2001 Elsevier Science B.V. All rights reserved.

Linkage of Paget disease of bone to a novel region on human chromosome 18q23

Paget disease of bone (PDB) is characterized by increased osteoclast activity and localized abnormal bone remodeling. PDB has a significant genetic component, with evidence of linkage to chromosomes 6p21.3 (PDB1) and 18q21-22 (PDB2) in some pedigrees. There is evidence of genetic heterogeneity, with other pedigrees showing negative linkage to these regions. TNFRSF11A, a gene that is essential for osteoclast formation and that encodes receptor activator of nuclear factor-kappa B (RANK), has been mapped to the PDB2 region. TNFRSF11A mutations that segregate in pedigrees with either familial expansile osteolysis or familial PDB have been identified; however, linkage studies and mutation screening have excluded the involvement of RANK in the majority of patients with PDB. We have excluded linkage, both to PDB1 and to PDB2, in a large multigenerational pedigree with multiple family members affected by PDB. We have conducted a genomewide scan of this pedigree, followed by fine mapping and multipoint analysis in regions of interest. The peak two-point LOD scores from the genomewide scan were 2.75, at D7S507, and 1.76, at D18S70. Multipoint and haplotype analysis of markers flanking D7S507 did not support linkage to this region. Haplotype analysis of markers flanking D18S70 demonstrated a haplotype segregating with PDB in a large subpedigree. This subpedigree had a significantly lower age at diagnosis than the rest of the pedigree (51.2 +/- 8.5 vs. 64.2 +/- 9.7 years; P = .0012). Linkage analysis of this subpedigree demonstrated a peak two-point LOD score of 4.23, at marker D18S1390 (theta = 0), and a peak multipoint LOD score of 4.71, at marker D18S70. Our data are consistent with genetic heterogeneity within the pedigree and indicate that 18q23 harbors a novel susceptibility gene for PDB.

An inverse design of an open, head/neck RF coil for MRI

Radio-frequency (RF) coils are a necessary component of magnetic resonance imaging (MRI) systems. When used in transmit operation, they act to generate a homogeneous RF magnetic field within a volume of interest and when in receive operation, they act to receive the nuclear magnetic resonance signal from the RF-excited specimen. This paper outlines a procedure for the design of open RF coils using the time-harmonic inverse method. This method entails the calculation of an ideal current density on a multipaned planar surface that would generate a specified magnetic field within the volume of interest. Because of the averaging effect of the regularization technique in the matrix solution, the specified magnetic field is shaped within an iterative procedure until the generated magnetic field matches the desired magnetic field. The stream-function technique is used to ascertain conductor positions and a method of moments package is then used to finalize the design. An open head/neck coil was designed to operate in a clinical 2T MRI system and the presented results prove the efficacy of this design methodology.

The anammox case - A new experimental manifesto for microbiological eco-physiology

The anaerobic ammonium oxidation process is a new process for ammonia removal from wastewater. It is also a new microbial physiology that was previously believed to be impossible. The identification of Candidatus Brocadia anammoxidans and its relatives as the responsible bacteria was only possible with the development of a new experimental approach. That approach is the focus of this paper. The approach is a modernisation of the Winogradsky/Beyerinck strategy of selective enrichment and is based on the introduction of the molecular toolbox and modern bioreactor engineering to microbial ecology. It consists of five steps: (1) postulation of an ecological niche based on thermodynamic considerations and macro-ecological field data; (2) engineering of this niche into a laboratory bioreactor for enrichment culture; (3) black-box physiological characterisation of the enrichment culture as a whole; (4) phylogenetic characterisation of the enriched community using molecular tools; (5) physical separation of the dominant members of the enrichment culture using gradient centrifugation and the identification of the species of interest in accordance with Koch's postulates; (6) verification of the in situ importance of these species in the actual ecosystems. The power of this approach is illustrated with a case study: the identification of the planctomycetes responsible for anaerobic ammonium oxidation. We argue that this was impossible using molecular ecology or conventional 'cultivation based techniques' alone. We suggest that the approach might also be used for the microbiological study of many interesting microbes such as anaerobic methane oxidisers.

Lower Bounds on the State Complexity of Geometric Goppa Codes

We reinterpret the state space dimension equations for geometric Goppa codes. An easy consequence is that if deg G less than or equal to n-2/2 or deg G greater than or equal to n-2/2 + 2g then the state complexity of C-L(D, G) is equal to the Wolf bound. For deg G is an element of [n-1/2, n-3/2 + 2g], we use Clifford's theorem to give a simple lower bound on the state complexity of C-L(D, G). We then derive two further lower bounds on the state space dimensions of C-L(D, G) in terms of the gonality sequence of F/F-q. (The gonality sequence is known for many of the function fields of interest for defining geometric Goppa codes.) One of the gonality bounds uses previous results on the generalised weight hierarchy of C-L(D, G) and one follows in a straightforward way from first principles; often they are equal. For Hermitian codes both gonality bounds are equal to the DLP lower bound on state space dimensions. We conclude by using these results to calculate the DLP lower bound on state complexity for Hermitian codes.

Production of the baculovirus-expressed dengue virus glycoprotein NS1 can be improved dramatically with optimised regimes for fed-batch cultures and the addition of the insect moulting hormone, 20-Hydroxyecdysone

A perennial problem in recombinant protein expression is low yield of the product of interest. A strategy which has been shown to increase the production of baculovirus-expressed proteins is to utilise fed-batch cultures. One disadvantage of this approach is the time-consuming task of optimising the feeding strategy. Previously, a statistical optimisation routine was applied to develop a feeding strategy that increased the yield of beta-Galactosidase (beta-Gal) by 2.4-fold (Biotechnol. Bioeng, 59 (1998) 178). This involves the single addition of nutrient concentrates (amino acids, lipids. glucose and yeastolate ultrafiltrate) into Sf9 cell cultures grown in SF900II medium. In this study, it is demonstrated that this optimised fed-batch strategy developed for a high-yielding intracellular product beta-Gal could be applied successfully to a relatively low-yielding glycosylated and secreted product such as the dengue virus glycoprotein NS1. Optimised batch infections yielded 4 mug/ml of NS1 at a peak cell density of 4.2 x 10(6) cells/ml. In contrast. optimised fed-batch infections exhibited a 3-fold improvement in yield, with 12 mug ml of NS1 produced at a peak cell density of 11.3 x 10(6) cells/ml. No further improvements in yield were recorded when the feed volumes were doubled and the peak cell density was increased to 23 x 10(6) cells/ml, unless the cultures were stimulated by the addition of 4 mug/ml of 20-Hydroxyecdysone (an insect moulting hormone). In this case, the NS1 yield was increased to 20 mug/ml. which was nearly 5-fold higher than optimised batch cultures. (C) 2002 Elsevier Science B.V. All rights reserved.

Use of heterologously-expressed cytochomre P450 and glutathione transferase enzymes in toxicity assays

Our groups have had a long-term interest in utilizing bacterial systems in the characterization of bioactivation and detoxication reactions catalyzed by cytochrome P450 (P450) and glutathione transferase (GST) enzymes. Bacterial systems remain the first choice for initial screens with new chemicals and have advantages, including high-throughput capability. Most human P450s of interest in toxicology have been readily expressed in Escherichia coli with only minor sequence modification. These enzymes can be readily purified and used in assays of activation of chemicals. Bicistronic systems have been developed in order to provide the auxiliary NADPH-P450 reductase. Alternative systems involve these enzymes expressed together within bacteria. In one approach, a lac selection system is used with E. coli and has been applied to the characterization of inhibitors of P450s 1A2 and 1131, as well as in basic studies involving random mutagenesis. Another approach utilizes induction of the SOS (umu) response in Salmonella typhimurium, and systems have now been developed with human P450s 1A1, 1A2, 1B1, 2C9, 2D6, 2E1, and 3A4, which have been used to report responses from heterocyclic amines. S. typhimurium his reporter systems have also been used with GSTs, first to demonstrate the role of rat GST 5-5 in the activation of dihalomethanes. These systems have been used to compare these GSTs with regard to activation of dihaloalkanes and potential toxicity. (C) 2002 Elsevier Science Ireland Ltd. All rights reserved.