Bayesian estimation of two-part joint models for a longitudinal semicontinuous biomarker and a terminal event with INLA: Interests for cancer clinical trial evaluation

Two-part joint models for a longitudinal semicontinuous biomarker and a terminal event have been recently introduced based on frequentist estimation. The biomarker distribution is decomposed into a probability of positive value and the expected value among positive values. Shared random effects can represent the association structure between the biomarker and the terminal event. The computational burden increases compared to standard joint models with a single regression model for the biomarker. In this context, the frequentist estimation implemented in the R package frailtypack can be challenging for complex models (i.e., a large number of parameters and dimension of the random effects). As an alternative, we propose a Bayesian estimation of two-part joint models based on the Integrated Nested Laplace Approximation (INLA) algorithm to alleviate the computational burden and fit more complex models. Our simulation studies confirm that INLA provides accurate approximation of posterior estimates and to reduced computation time and variability of estimates compared to frailtypack in the situations considered. We contrast the Bayesian and frequentist approaches in the analysis of two randomized cancer clinical trials (GERCOR and PRIME studies), where INLA has a reduced variability for the association between the biomarker and the risk of event. Moreover, the Bayesian approach was able to characterize subgroups of patients associated with different responses to treatment in the PRIME study. Our study suggests that the Bayesian approach using the INLA algorithm enables to fit complex joint models that might be of interest in a wide range of clinical applications.     

Contribution of the different planktonic microbial assemblages to ETS activity in the Ligurian frontal area: north-west Mediterranean Sea

Spatial and size distribution of micro-organisms and their ETSactivity has been investigated in Ligurian Sea surface watersalong the Nice-Calvi transect across frontal areas from 18 to37 km offshore (TOMOFRONT 1 and 2 cruises, April 1988 and April-May1989 respectively). Aplastidic and plastidic nanoflagellatesand aplastidic picoflagellates were present in numbers closeto 0.25 x 10⁴ cells ml^–1, whereas plastidic picoflagellatesaccounted for about half this number. Correlations have beenevidenced between plastidic and aplastidic micro-organisms withinthe same size group, suggesting that they belong to a well-definedecosystem. The highest correlation between total ETS activityand abundance of the considered size groups was observed fornanoflagellates (r = 0.94, n = 22, and r = 0.90, n = 22 foraplastidic and plastidic cells respectively). The importanceof the role of nanoflagellates in surface waters, with respectto the overall ETS activity, was supported by results from sizefractionation which assigned to the 3–10 µm sizerange a 73.3% contribution to overall ETS activity. Resultsemphasize analysing global ETS activity of natural samples inorder to derive relationships between the different populationspresent in the sampled water. It is suggested that couplingflow cytometry to the ETS approach should be very helpful inthat respect.     

A Quantitative Technique to Compare and Classify Humpback Whale (Megaptera novaeangliae) Sounds

In an attempt to minimize observer bias, numerical taxonomy methods were used to describe and classify humpback whale sounds. The spectrograms (N = 1255) were digitized into a 16 × 21 binary matrix. The rows were 16 frequencies selected on a logarithmic scale (0.12–8 kHz). The columns were 21 time samples taken every 0.1 s. Each point of the matrix was coded 1 if it lay over part of the sound. Other binary variables were added to code for relative intensity within a sound, frequency modulation and amplitude modulation. The sounds were then compared using the Jaccard similarity coefficient for binary data, and classified with average linkage cluster analysis. This technique produced 115 clusters, which were compared with my aural and visual impressions of the sounds. I agreed with most major categories identified by cluster analysis, but many small clusters had to be fused to other categories. This was partially due to the technique used, and to the complexity of the repertoire under study. Improvements are proposed to further reduce observer bias in classification of sounds, and thus make studies of animal communication performed by different researchers or on different species more easily comparable.     

Characterization of anther-expressed genes encoding a major class of extracellular oleosin-like proteins in the pollen coat of Brassicaceae†

A large, heterogeneous, highly expressed gene family encoding oleosin-like proteins is described in the Brassicaceae. íeven related cDNA sequences were isolated from Brassica napus anther mRNA using RACE-PCR and compared with other recently described anther-specific oleosin-like genes from B. napus. The expression patterns of four representative members of this diverse gene family were analyzed by Northern blotting and in situ hybridization. In all cases, the genes were expressed specifically in the tapetum of 3–5 mm B. napus buds, which contained microspores at the late-vacuolate and bicellular stages of development. The predicted protein products are ordered into subclasses, each of which has a characteristic C-terminal domain, containing different amino acid motifs or repeated residues. Tryphine (pollen coat) fractions from mature B. napus pollen were found to be particularly enriched in polypeptides of apparent molecular weights 32–38 kDa, plus numerous less abundant polypeptides of less than 15 kDa. The N-terminal 15–20 residues of three of these polypeptides (12, 32 and 38 kDa) were found by microsequencing to be identical to parts of the predicted amino acid sequences of three of the tapetal-expressed oleosin-like genes. This indicates the possibility of post-translational modification of these proteins resulting in a cleavage of the primary translation products in order to generate the mature tryphine polypeptides. These data imply that a large and diverse group of oleosin-like proteins is synthesized in the tapeturn of B. napus anthers and that following tapetal degradation, these proteins, possibly in modified form, then relocate to the developing microspores where they eventually constitute some of the major components of the extracellular tryphine of mature pollen grains. These proteins share a conserved 70 amino acid residue hydrophobic domain and are related structurally to the seed-specific intracellular oleosins, although their biological function may be different.     

Analysis of human V H gene repertoire expression in peripheral CD19+ B cells

Using CD19 B-cell selection and polymerase chain reaction-amplified cDNA libraries, we analyzed the peripheral immunoglobulin heavy chain variable repertoire of three healthy adult donors. Here we report that most of the CD19+ circulating B cells expressed unmutated V _H-D-J_H rearrangements. By specific V _H family hybridization, we show that V _H gene family utilization in the periphery roughly corresponds to the complexity of these families in the germline and appears to be relatively constant among the analyzed subjects. However, sequence data of clones picked at random from one IgM cDNA library reveals that in spite of this random utilization, the V _H gene expression in naive circulating B cells is highly biased towards the expression of a limited set of V _H genes. As previously reported by others, this restricted mechanism is also found for the D and J _H segments.The nucleotide sequence data reported in this paper have been submitted to the GenBank/EMBL nucleotide sequence database and have been assigned the accession numbers Z47213-Z47243 and Z47349     

Zolpidem Displays Heterogeneity in Its Binding to the Nonhuman Primate Benzodiazepine Receptor In Vivo

Abstract: The distinctive pharmacological activity of zolpidem in rats compared with classical benzodiazepines has been related to its differential affinity for benzodiazepine receptor (BZR) subtypes. By contrast, in nonhuman primates the pharmacological activity of zolpidem was found to be quite similar to that of classical BZR agonists. In an attempt to explain this discrepancy, we examined the ability of zolpidem to differentiate BZR subtypes in vivo in primate brain using positron emission tomography. The BZRs were specifically labeled with [¹¹C]flumazenil. Radiotracer displacement by zolpidem was monophasic in cerebellum and neocortex, with in vivo Hill coefficients close to 1. Conversely, displacement of [¹¹C]flumazenil was biphasic in hippocampus, amygdala, septum, insula, striatum, and pons, with Hill coefficients significantly smaller than 1, suggesting two different binding sites for zolpidem. In these cerebral regions, the half-maximal inhibitory doses for the high-affinity binding site were similar to those found in cerebellum and neocortex and ～100-fold higher for the low-affinity binding site. The low-affinity binding site accounted for <32% of the specific [¹¹C]-flumazenil binding. Such zolpidem binding characteristics contrast with those reported for rodents, where three different binding sites were found. Species differences in binding characteristics may explain why zolpidem has a distinctive pharmacological activity in rodents, whereas its pharmacological activity in primates is quite similar to that of classical BZR agonists, except for the absence of severe effects on memory functions, which may be due to the lack of substantial zolpidem affinity for a distinct BZR subtype in cerebral structures belonging to the limbic system.     

The DnaJ domain of polyomavirus large T antigen is required to regulate Rb family tumor suppressor function.

Tumor suppressors of the retinoblastoma susceptibility gene family regulate cell growth and differentiation. Polyomavirus large T antigens (large T) bind Rb family members and block their function. Mutations of large T sequences conserved with the DnaJ family affect large T binding to a cellular DnaK, heat shock protein 70. The same mutations abolish large T activation of E2F-containing promoters and Rb binding-dependent large T activation of cell cycle progression. Cotransfection of a cellular DnaJ domain blocks wild-type large T action, showing that the connection between the chaperone system and tumor suppressors is direct. Although they are inactive in assays dependent on Rb family binding, mutants in the J region retain the ability to associate with pRb, p107, and p130. This suggests that binding of Rb family members by large T is not sufficient for their inactivation and that a functional J domain is required as well. This work connects the DnaJ and DnaK molecular chaperones to regulation of tumor suppressors by polyomavirus large T.     

Modeling of adaptations to physical training by using a recursive least squares algorithm

Busso, Thierry, Christian Denis, Régis Bonnefoy,André Geyssant, and Jean-René Lacour. Modeling ofadaptations to physical training by using a recursive least squaresalgorithm. J. Appl. Physiol. 82(5):1685-1693, 1997.The present study assesses the usefulnessof a systems model with time-varying parameters for describing theresponses of physical performance to training. Data for two subjectswho undertook a 14-wk training on a cycle ergometer were used to testthe proposed model, and the results were compared with a model withtime-invariant parameters. Two 4-wk periods of intensive training wereseparated by a 2-wk period of reduced training and followed by a 4-wkperiod of reduced training. The systems input ascribed to the trainingdoses was made up of interval exercises and computed in arbitraryunits. The systems output was evaluated one to five times per week byusing the endurance time at a constant workload. The time-invariantparameters were fitted from actual performances by using the leastsquares method. The time-varying parameters were fitted by using arecursive least squares algorithm. The coefficients of determinationr² were 0.875 and0.879 for the two subjects using the time-varying model, higher thanthe values of 0.682 and 0.666, respectively, obtained with thetime-invariant model. The variations over time in the model parametersresulting from the expected reduction in the residuals appearedgenerally to account for changes in responses to training. Such a modelwould be useful for investigating the underlying mechanisms ofadaptation and fatigue.