A Pseudolikelihood Approach for Simultaneous Analysis of Array Comparative Genomic Hybridizations (aCGH)

DNA sequence copy number has been shown to be associated with cancer development and progression. Array-based Comparative Genomic Hybridization (aCGH) is a recent development that seeks to identify the copy number ratio at large numbers of markers across the genome. Due to experimental and biological variations across chromosomes and across hybridizations, current methods are limited to analyses of single chromosomes. We propose a more powerful approach that borrows strength across chromosomes and across hybridizations. We assume a Gaussian mixture model, with a hidden Markov dependence structure, and with random effects to allow for intertumoral variation, as well as intratumoral clonal variation. For ease of computation, we base estimation on a pseudolikelihood function. The method produces quantitative assessments of the likelihood of genetic alterations at each clone, along with a graphical display for simple visual interpretation. We assess the characteristics of the method through simulation studies and through analysis of a brain tumor aCGH data set. We show that the pseudolikelihood approach is superior to existing methods both in detecting small regions of copy number alteration and in accurately classifying regions of change when intratumoral clonal variation is present.

RANDOM EFFECTS MODELS IN A META-ANALYSIS OF THE ACCURACY OF DIAGNOSTIC TESTS WITHIN A GOLD STANDARD IN THE PRESENCE OF MISSING DATA

In evaluating the accuracy of diagnosis tests, it is common to apply two imperfect tests jointly or sequentially to a study population. In a recent meta-analysis of the accuracy of microsatellite instability testing (MSI) and traditional mutation analysis (MUT) in predicting germline mutations of the mismatch repair (MMR) genes, a Bayesian approach (Chen, Watson, and Parmigiani 2005) was proposed to handle missing data resulting from partial testing and the lack of a gold standard. In this paper, we demonstrate an improved estimation of the sensitivities and specificities of MSI and MUT by using a nonlinear mixed model and a Bayesian hierarchical model, both of which account for the heterogeneity across studies through study-specific random effects. The methods can be used to estimate the accuracy of two imperfect diagnostic tests in other meta-analyses when the prevalence of disease, the sensitivities and/or the specificities of diagnostic tests are heterogeneous among studies. Furthermore, simulation studies have demonstrated the importance of carefully selecting appropriate random effects on the estimation of diagnostic accuracy measurements in this scenario.

Biological safety concepts of genetically modified live bacterial vaccines

Live vaccines possess the advantage of having access to induce cell-mediated and antibody-mediated immunity; thus in certain cases they are able to prevent infection, and not only disease. Furthermore, live vaccines, particularly bacterial live vaccines, are relatively cheap to produce and easy to apply. Hence they are suitable to immunize large communities or herds. The induction of both cell-mediated immunity as well as antibody-mediated immunity, which is particularly beneficial in inducing mucosal immune responses, is obtained by the vaccine-strain's ability to colonize and multiply in the host without causing disease. For this reason, live vaccines require attenuation of virulence of the bacterium to which immunity must be induced. Traditionally attenuation was achieved simply by multiple passages of the microorganism on growth medium, in animals, eggs or cell cultures or by chemical or physical mutagenesis, which resulted in random mutations that lead to attenuation. In contrast, novel molecular methods enable the development of genetically modified organisms (GMOs) targeted to specific genes that are particularly suited to induce attenuation or to reduce undesirable effects in the tissue in which the vaccine strains can multiply and survive. Since live vaccine strains (attenuated by natural selection or genetic engineering) are potentially released into the environment by the vaccinees, safety issues concerning the medical as well as environmental aspects must be considered. These involve (i) changes in cell, tissue and host tropism, (ii) virulence of the carrier through the incorporation of foreign genes, (iii) reversion to virulence by acquisition of complementation genes, (iv) exchange of genetic information with other vaccine or wild-type strains of the carrier organism and (v) spread of undesired genes such as antibiotic resistance genes. Before live vaccines are applied, the safety issues must be thoroughly evaluated case-by-case. Safety assessment includes knowledge of the precise function and genetic location of the genes to be mutated, their genetic stability, potential reversion mechanisms, possible recombination events with dormant genes, gene transfer to other organisms as well as gene acquisition from other organisms by phage transduction, transposition or plasmid transfer and cis- or trans-complementation. For this, GMOs that are constructed with modern techniques of genetic engineering display a significant advantage over random mutagenesis derived live organisms. The selection of suitable GMO candidate strains can be made under in vitro conditions using basic knowledge on molecular mechanisms of pathogenicity of the corresponding bacterial species rather than by in vivo testing of large numbers of random mutants. This leads to a more targeted safety testing on volunteers and to a reduction in the use of animal experimentation.

Re-evaluation of pulmonary titanium dioxide nanoparticle distribution using the "relative deposition index": Evidence for clearance through microvasculature

ABSTRACT: BACKGROUND: Translocation of nanoparticles (NP) from the pulmonary airways into other pulmonary compartments or the systemic circulation is controversially discussed in the literature. In a previous study it was shown that titanium dioxide (TiO2) NP were "distributed in four lung compartments (air-filled spaces, epithelium/endothelium, connective tissue, capillary lumen) in correlation with compartment size". It was concluded that particles can move freely between these tissue compartments. To analyze whether the distribution of TiO2 NP in the lungs is really random or shows a preferential targeting we applied a newly developed method for comparing NP distributions. METHODS: Rat lungs exposed to an aerosol containing TiO2 NP were prepared for light and electron microscopy at 1 h and at 24 h after exposure. Numbers of TiO2 NP associated with each compartment were counted using energy filtering transmission electron microscopy. Compartment size was estimated by unbiased stereology from systematically sampled light micrographs. Numbers of particles were related to compartment size using a relative deposition index and chi-squared analysis. RESULTS: Nanoparticle distribution within the four compartments was not random at 1 h or at 24 h after exposure. At 1 h the connective tissue was the preferential target of the particles. At 24 h the NP were preferentially located in the capillary lumen. CONCLUSION: We conclude that TiO2 NP do not move freely between pulmonary tissue compartments, although they can pass from one compartment to another with relative ease. The residence time of NP in each tissue compartment of the respiratory system depends on the compartment and the time after exposure. It is suggested that a small fraction of TiO2 NP are rapidly transported from the airway lumen to the connective tissue and subsequently released into the systemic circulation.

Conversion of CT numbers into tissue parameters for Monte Carlo dose calculations: a multi-centre study

The conversion of computed tomography (CT) numbers into material composition and mass density data influences the accuracy of patient dose calculations in Monte Carlo treatment planning (MCTP). The aim of our work was to develop a CT conversion scheme by performing a stoichiometric CT calibration. Fourteen dosimetrically equivalent tissue subsets (bins), of which ten bone bins, were created. After validating the proposed CT conversion scheme on phantoms, it was compared to a conventional five bin scheme with only one bone bin. This resulted in dose distributions D(14) and D(5) for nine clinical patient cases in a European multi-centre study. The observed local relative differences in dose to medium were mostly smaller than 5%. The dose-volume histograms of both targets and organs at risk were comparable, although within bony structures D(14) was found to be slightly but systematically higher than D(5). Converting dose to medium to dose to water (D(14) to D(14wat) and D(5) to D(5wat)) resulted in larger local differences as D(5wat) became up to 10% higher than D(14wat). In conclusion, multiple bone bins need to be introduced when Monte Carlo (MC) calculations of patient dose distributions are converted to dose to water.

Confidence Bands for Convex Median Curves using Sign-Tests

Suppose that one observes pairs (x1,Y1), (x2,Y2), ..., (xn,Yn), where x1 < x2 < ... < xn are fixed numbers while Y1, Y2, ..., Yn are independent random variables with unknown distributions. The only assumption is that Median(Yi) = f(xi) for some unknown convex or concave function f. We present a confidence band for this regression function f using suitable multiscale sign tests. While the exact computation of this band seems to require O(n4) steps, good approximations can be obtained in O(n2) steps. In addition the confidence band is shown to have desirable asymptotic properties as the sample size n tends to infinity.

Comparison of amplification methods for transcriptomic analyses of low abundance prokaryotic RNA sources

Microarrays have established as instrumental for bacterial detection, identification, and genotyping as well as for transcriptomic studies. For gene expression analyses using limited numbers of bacteria (derived from in vivo or ex vivo origin, for example), RNA amplification is often required prior to labeling and hybridization onto microarrays. Evaluation of the fidelity of the amplification methods is crucial for the robustness and reproducibility of microarray results. We report here the first utilization of random primers and the highly processive Phi29 phage polymerase to amplify material for transcription profiling analyses. We compared two commercial amplification methods (GenomiPhi and MessageAmp kits) with direct reverse-transcription as the reference method, focusing on the robustness of mRNA quantification using either microarrays or quantitative RT-PCR. Both amplification methods using either poly-A tailing followed by in vitro transcription, or direct strand displacement polymerase, showed appreciable linearity. Strand displacement technique was particularly affordable compared to in vitro transcription-based (IVT) amplification methods and consisted in a single tube reaction leading to high amplification yields. Real-time measurements using low-, medium-, and highly expressed genes revealed that this simple method provided linear amplification with equivalent results in terms of relative messenger abundance as those obtained by conventional direct reverse-transcription.