Bootstrapping pairs in Distance-Based Regression

La regressió basada en distàncies és un mètode de predicció que consisteix en dos passos: a partir de les distàncies entre observacions obtenim les variables latents, les quals passen a ser els regressors en un model lineal de mínims quadrats ordinaris. Les distàncies les calculem a partir dels predictors originals fent us d'una funció de dissimilaritats adequada. Donat que, en general, els regressors estan relacionats de manera no lineal amb la resposta, la seva selecció amb el test F usual no és possible. En aquest treball proposem una solució a aquest problema de selecció de predictors definint tests estadístics generalitzats i adaptant un mètode de bootstrap no paramètric per a l'estimació dels p-valors. Incluim un exemple numèric amb dades de l'assegurança d'automòbils.

Detection of live and antibiotic-killed bacteria by quantitative real-time PCR of specific fragments of rRNA.

Assessing bacterial viability by molecular markers might help accelerate the measurement of antibiotic-induced killing. This study investigated whether rRNA could be suitable for this purpose. Cultures of penicillin-susceptible and penicillin-tolerant (Tol1 mutant) Streptococcus gordonii were exposed to mechanistically different penicillin and levofloxacin. Bacterial survival was assessed by viable counts and compared to quantitative real-time PCR amplification of either the 16S rRNA genes or the 16S rRNA, following reverse transcription. Penicillin-susceptible S. gordonii lost > or =4 log(10) CFU/ml of viability over 48 h of penicillin treatment. In comparison, the Tol1 mutant lost < or =1 log(10) CFU/ml. Amplification of a 427-bp fragment of 16S rRNA genes yielded amplicons that increased proportionally to viable counts during bacterial growth but did not decrease during drug-induced killing. In contrast, the same 427-bp fragment amplified from 16S rRNA paralleled both bacterial growth and drug-induced killing. It also differentiated between penicillin-induced killing of the parent and the Tol1 mutant (> or =4 log(10) CFU/ml and < or =1 log(10) CFU/ml, respectively) and detected killing by mechanistically unrelated levofloxacin. Since large fragments of polynucleotides might be degraded faster than smaller fragments, the experiments were repeated by amplifying a 119-bp region internal to the original 427-bp fragment. The amount of 119-bp amplicons increased proportionally to viability during growth but remained stable during drug treatment. Thus, 16S rRNA was a marker of antibiotic-induced killing, but the size of the amplified fragment was critical for differentiation between live and dead bacteria.

Semiquantitative RT-PCR measurement of gene expression in rat tissues including a correction for varying cell size and number

Background: Current methodology of gene expression analysis limits the possibilities of comparison between cells/tissues of organs in which cell size and/or number changes as a consequence of the study (e.g. starvation). A method relating the abundance of specific mRNA copies per cell may allow direct comparison or different organs and/or changing physiological conditions. Methods: With a number of selected genes, we analysed the relationship of the number of bases and the fluorescence recorded at a present level using cDNA standards. A lineal relationship was found between the final number of bases and the length of the transcript. The constants of this equation and those of the relationship between fluorescence and number of bases in cDNA were determined and a general equation linking the length of the transcript and the initial number of copies of mRNA was deduced for a given pre-established fluorescence setting. This allowed the calculation of the concentration of the corresponding mRNAs per g of tissue. The inclusion of tissue RNA and the DNA content per cell, allowed the calculation of the mRNA copies per cell. Results: The application of this procedure to six genes: Arbp, cyclophilin, ChREBP, T4 deiodinase 2, acetyl-CoA carboxylase 1 and IRS-1, in liver and retroperitoneal adipose tissue of food-restricted rats allowed precise measures of their changes irrespective of the shrinking of the tissue, the loss of cells or changes in cell size, factors that deeply complicate the comparison between changing tissue conditions. The percentage results obtained with the present methods were essentially the same obtained with the delta-delta procedure and with individual cDNA standard curve quantitative RT-PCR estimation. Conclusion: The method presented allows the comparison (i.e. as copies of mRNA per cell) between different genes and tissues, establishing the degree of abundance of the different molecular species tested.

Identification by Real-time PCR of 13 mature microRNAs differentially expressed in colorectal cancer and non-tumoral tissues

MicroRNAs (miRNAs) are short non-coding RNA molecules playing regulatory roles by repressing translation or cleaving RNA transcripts. Although the number of verified human miRNA is still expanding, only few have been functionally described. However, emerging evidences suggest the potential involvement of altered regulation of miRNA in pathogenesis of cancers and these genes are thought to function as both tumours suppressor and oncogenes. In our study, we examined by Real-Time PCR the expression of 156 mature miRNA in colorectal cancer. The analysis by several bioinformatics algorithms of colorectal tumours and adjacent non-neoplastic tissues from patients and colorectal cancer cell lines allowed identifying a group of 13 miRNA whose expression is significantly altered in this tumor. The most significantly deregulated miRNA being miR-31, miR-96, miR-133b, miR-135b, miR-145, and miR-183. In addition, the expression level of miR-31 was correlated with the stage of CRC tumor. Our results suggest that miRNA expression profile could have relevance to the biological and clinical behavior of colorectal neoplasia.

Nonspecific PCR amplification by high-fidelity polymerases: implications for next-generation sequencing of AFLP markers.

High-fidelity 'proofreading' polymerases are often used in library construction for next-generation sequencing projects, in an effort to minimize errors in the resulting sequence data. The increased template fidelity of these polymerases can come at the cost of reduced template specificity, and library preparation methods based on the AFLP technique may be particularly susceptible. Here, we compare AFLP profiles generated with standard Taq and two versions of a high-fidelity polymerase. We find that Taq produces fewer and brighter peaks than high-fidelity polymerase, suggesting that Taq performs better at selectively amplifying templates that exactly match the primer sequences. Because the higher accuracy of proofreading polymerases remains important for sequencing applications, we suggest that it may be more effective to use alternative library preparation methods.

Validation of a Low-Cost Human Papillomavirus Genotyping Assay Based on PGMY PCR and Reverse Blotting Hybridization with Reusable Membranes.

The genotyping of human papillomaviruses (HPV) is essential for the surveillance of HPV vaccines. We describe and validate a low-cost PGMY-based PCR assay (PGMY-CHUV) for the genotyping of 31 HPV by reverse blotting hybridization (RBH). Genotype-specific detection limits were 50 to 500 genome equivalents per reaction. RBH was 100% specific and 98.61% sensitive using DNA sequencing as the gold standard (n = 1,024 samples). PGMY-CHUV was compared to the validated and commercially available linear array (Roche) on 200 samples. Both assays identified the same positive (n = 182) and negative samples (n = 18). Seventy-six percent of the positives were fully concordant after restricting the comparison to the 28 genotypes shared by both assays. At the genotypic level, agreement was 83% (285/344 genotype-sample combinations; κ of 0.987 for single infections and 0.853 for multiple infections). Fifty-seven of the 59 discordant cases were associated with multiple infections and with the weakest genotypes within each sample (P < 0.0001). PGMY-CHUV was significantly more sensitive for HPV56 (P = 0.0026) and could unambiguously identify HPV52 in mixed infections. PGMY-CHUV was reproducible on repeat testing (n = 275 samples; 392 genotype-sample combinations; κ of 0.933) involving different reagents lots and different technicians. Discordant results (n = 47) were significantly associated with the weakest genotypes in samples with multiple infections (P < 0.0001). Successful participation in proficiency testing also supported the robustness of this assay. The PGMY-CHUV reagent costs were estimated at $2.40 per sample using the least expensive yet proficient genotyping algorithm that also included quality control. This assay may be used in low-resource laboratories that have sufficient manpower and PCR expertise.

A simple, robust and rapid approach to detect carbapenemases in Gram-negative isolates by MALDI-TOF mass spectrometry: validation with triple quadripole tandem mass spectrometry, microarray and PCR.

Carbapenemases should be accurately and rapidly detected, given their possible epidemiological spread and their impact on treatment options. Here, we developed a simple, easy and rapid matrix-assisted laser desorption ionization-time of flight (MALDI-TOF)-based assay to detect carbapenemases and compared this innovative test with four other diagnostic approaches on 47 clinical isolates. Tandem mass spectrometry (MS-MS) was also used to determine accurately the amount of antibiotic present in the supernatant after 1 h of incubation and both MALDI-TOF and MS-MS approaches exhibited a 100% sensitivity and a 100% specificity. By comparison, molecular genetic techniques (Check-MDR Carba PCR and Check-MDR CT103 microarray) showed a 90.5% sensitivity and a 100% specificity, as two strains of Aeromonas were not detected because their chromosomal carbapenemase is not targeted by probes used in both kits. Altogether, this innovative MALDI-TOF-based approach that uses a stable 10-μg disk of ertapenem was highly efficient in detecting carbapenemase, with a sensitivity higher than that of PCR and microarray.

MGMT methylation analysis of glioblastoma on the Infinium methylation BeadChip identifies two distinct CpG regions associated with gene silencing and outcome, yielding a prediction model for comparisons across datasets, tumor grades, and CIMP-status.

The methylation status of the O(6)-methylguanine-DNA methyltransferase (MGMT) gene is an important predictive biomarker for benefit from alkylating agent therapy in glioblastoma. Recent studies in anaplastic glioma suggest a prognostic value for MGMT methylation. Investigation of pathogenetic and epigenetic features of this intriguingly distinct behavior requires accurate MGMT classification to assess high throughput molecular databases. Promoter methylation-mediated gene silencing is strongly dependent on the location of the methylated CpGs, complicating classification. Using the HumanMethylation450 (HM-450K) BeadChip interrogating 176 CpGs annotated for the MGMT gene, with 14 located in the promoter, two distinct regions in the CpG island of the promoter were identified with high importance for gene silencing and outcome prediction. A logistic regression model (MGMT-STP27) comprising probes cg1243587 and cg12981137 provided good classification properties and prognostic value (kappa = 0.85; log-rank p < 0.001) using a training-set of 63 glioblastomas from homogenously treated patients, for whom MGMT methylation was previously shown to be predictive for outcome based on classification by methylation-specific PCR. MGMT-STP27 was successfully validated in an independent cohort of chemo-radiotherapy-treated glioblastoma patients (n = 50; kappa = 0.88; outcome, log-rank p < 0.001). Lower prevalence of MGMT methylation among CpG island methylator phenotype (CIMP) positive tumors was found in glioblastomas from The Cancer Genome Atlas than in low grade and anaplastic glioma cohorts, while in CIMP-negative gliomas MGMT was classified as methylated in approximately 50 % regardless of tumor grade. The proposed MGMT-STP27 prediction model allows mining of datasets derived on the HM-450K or HM-27K BeadChip to explore effects of distinct epigenetic context of MGMT methylation suspected to modulate treatment resistance in different tumor types.

Robust accelerated failure time regression

Robust estimators for accelerated failure time models with asymmetric (or symmetric) error distribution and censored observations are proposed. It is assumed that the error model belongs to a log-location-scale family of distributions and that the mean response is the parameter of interest. Since scale is a main component of mean, scale is not treated as a nuisance parameter. A three steps procedure is proposed. In the first step, an initial high breakdown point S estimate is computed. In the second step, observations that are unlikely under the estimated model are rejected or down weighted. Finally, a weighted maximum likelihood estimate is computed. To define the estimates, functions of censored residuals are replaced by their estimated conditional expectation given that the response is larger than the observed censored value. The rejection rule in the second step is based on an adaptive cut-off that, asymptotically, does not reject any observation when the data are generat ed according to the model. Therefore, the final estimate attains full efficiency at the model, with respect to the maximum likelihood estimate, while maintaining the breakdown point of the initial estimator. Asymptotic results are provided. The new procedure is evaluated with the help of Monte Carlo simulations. Two examples with real data are discussed.

DNA extraction using the QIAsymphony: Evaluation of PCR inhibitor removal

The goal of this study was to compare the quantity and purity of DNA extracted from biological tracesusing the QIAsymphony robot with that of the manual QIAamp DNA mini kit currently in use in ourlaboratory. We found that the DNA yield of robot was 1.6-3.5 times lower than that of the manualprotocol. This resulted in a loss of 8% and 29% of the alleles correctly scored when analyzing 1/400 and 1/800 diluted saliva samples, respectively. Specific tests showed that the QIAsymphony was at least 2-16times more efficient at removing PCR inhibitors. The higher purity of the DNA may therefore partlycompensate for the lower DNA yield obtained. No case of cross-contamination was observed amongsamples. After purification with the robot, DNA extracts can be automatically transferred in 96-wellsplates, which is an ideal format for subsequent RT-qPCR quantification and DNA amplification. Lesshands-on time and reduced risk of operational errors represent additional advantages of the robotic platform.

Concordance among digital gene expression, microarrays, and qPCR when measuring differential expression of microRNAs.

Profiling microRNA (miRNA) expression is of widespread interest given the critical role of miRNAs in many cellular functions. Profiling can be achieved via hybridization-based (microarrays), sequencing-based, or amplification-based (quantitative reverse transcription-PCR, qPCR) technologies. Among these, microarrays face the significant challenge of accurately distinguishing between mature and immature miRNA forms, and different vendors have developed different methods to meet this challenge. Here we measure differential miRNA expression using the Affymetrix, Agilent, and Illumina microarray platforms, as well as qPCR (Applied Biosystems) and ultra high-throughput sequencing (Illumina). We show that the differential expression measurements are more divergent when the three types of microarrays are compared than when the Agilent microarray, qPCR, and sequencing technology measurements are compared, which exhibit a good overall concordance.

Development of a real-time PCR for the specific detection of Waddlia chondrophila in clinical samples.

Waddlia chondrophila is considered as an emerging human pathogen likely involved in miscarriage and lower respiratory tract infections. Given the low sensitivity of cell culture to recover such an obligate intracellular bacteria, molecular-based diagnostic approaches are warranted. We thus developed a real-time PCR that amplifies Waddlia chondrophila DNA. Specific primers and probe were selected to target the 16S rRNA gene. The PCR specifically amplified W. chondrophila but did not amplify other related-bacteria such as Parachlamydia acanthamoebae, Simkania negevensis and Chlamydia pneumoniae. The PCR exhibited a good intra-run and inter-run reproducibility and a sensitivity of less than ten copies of the positive control. This real-time PCR was then applied to 32 nasopharyngeal aspirates taken from children with bronchiolitis not due to respiratory syncytial virus (RSV). Three samples revealed to be Waddlia positive, suggesting a possible role of this Chlamydia-related bacteria in this setting.

Blocking hypoxia-induced autophagy in tumors restores cytotoxic T-cell activity and promotes regression.

The relationship between hypoxic stress, autophagy, and specific cell-mediated cytotoxicity remains unknown. This study shows that hypoxia-induced resistance of lung tumor to cytolytic T lymphocyte (CTL)-mediated lysis is associated with autophagy induction in target cells. In turn, this correlates with STAT3 phosphorylation on tyrosine 705 residue (pSTAT3) and HIF-1α accumulation. Inhibition of autophagy by siRNA targeting of either beclin1 or Atg5 resulted in impairment of pSTAT3 and restoration of hypoxic tumor cell susceptibility to CTL-mediated lysis. Furthermore, inhibition of pSTAT3 in hypoxic Atg5 or beclin1-targeted tumor cells was found to be associated with the inhibition Src kinase (pSrc). Autophagy-induced pSTAT3 and pSrc regulation seemed to involve the ubiquitin proteasome system and p62/SQSTM1. In vivo experiments using B16-F10 melanoma tumor cells indicated that depletion of beclin1 resulted in an inhibition of B16-F10 tumor growth and increased tumor apoptosis. Moreover, in vivo inhibition of autophagy by hydroxychloroquine in B16-F10 tumor-bearing mice and mice vaccinated with tyrosinase-related protein-2 peptide dramatically increased tumor growth inhibition. Collectively, this study establishes a novel functional link between hypoxia-induced autophagy and the regulation of antigen-specific T-cell lysis and points to a major role of autophagy in the control of in vivo tumor growth.