Local and global error models to improve uncertainty quantification

In groundwater applications, Monte Carlo methods are employed to model the uncertainty on geological parameters. However, their brute-force application becomes computationally prohibitive for highly detailed geological descriptions, complex physical processes, and a large number of realizations. The Distance Kernel Method (DKM) overcomes this issue by clustering the realizations in a multidimensional space based on the flow responses obtained by means of an approximate (computationally cheaper) model; then, the uncertainty is estimated from the exact responses that are computed only for one representative realization per cluster (the medoid). Usually, DKM is employed to decrease the size of the sample of realizations that are considered to estimate the uncertainty. We propose to use the information from the approximate responses for uncertainty quantification. The subset of exact solutions provided by DKM is then employed to construct an error model and correct the potential bias of the approximate model. Two error models are devised that both employ the difference between approximate and exact medoid solutions, but differ in the way medoid errors are interpolated to correct the whole set of realizations. The Local Error Model rests upon the clustering defined by DKM and can be seen as a natural way to account for intra-cluster variability; the Global Error Model employs a linear interpolation of all medoid errors regardless of the cluster to which the single realization belongs. These error models are evaluated for an idealized pollution problem in which the uncertainty of the breakthrough curve needs to be estimated. For this numerical test case, we demonstrate that the error models improve the uncertainty quantification provided by the DKM algorithm and are effective in correcting the bias of the estimate computed solely from the MsFV results. The framework presented here is not specific to the methods considered and can be applied to other combinations of approximate models and techniques to select a subset of realizations

Unedited in vivo detection and quantification of γ-aminobutyric acid in the occipital cortex using short-TE MRS at 3 T.

Short-TE MRS has been proposed recently as a method for the in vivo detection and quantification of γ-aminobutyric acid (GABA) in the human brain at 3 T. In this study, we investigated the accuracy and reproducibility of short-TE MRS measurements of GABA at 3 T using both simulations and experiments. LCModel analysis was performed on a large number of simulated spectra with known metabolite input concentrations. Simulated spectra were generated using a range of spectral linewidths and signal-to-noise ratios to investigate the effect of varying experimental conditions, and analyses were performed using two different baseline models to investigate the effect of an inaccurate baseline model on GABA quantification. The results of these analyses indicated that, under experimental conditions corresponding to those typically observed in the occipital cortex, GABA concentration estimates are reproducible (mean reproducibility error, <20%), even when an incorrect baseline model is used. However, simulations indicate that the accuracy of GABA concentration estimates depends strongly on the experimental conditions (linewidth and signal-to-noise ratio). In addition to simulations, in vivo GABA measurements were performed using both spectral editing and short-TE MRS in the occipital cortex of 14 healthy volunteers. Short-TE MRS measurements of GABA exhibited a significant positive correlation with edited GABA measurements (R = 0.58, p < 0.05), suggesting that short-TE measurements of GABA correspond well with measurements made using spectral editing techniques. Finally, within-session reproducibility was assessed in the same 14 subjects using four consecutive short-TE GABA measurements in the occipital cortex. Across all subjects, the average coefficient of variation of these four GABA measurements was 8.7 ± 4.9%. This study demonstrates that, under some experimental conditions, short-TE MRS can be employed for the reproducible detection of GABA at 3 T, but that the technique should be used with caution, as the results are dependent on the experimental conditions. Copyright © 2013 John Wiley & Sons, Ltd.

The SOS screen in Arabidopsis: a search for functions involved in DNA metabolism.

The SOS screen, as originally described by Perkins et al. (1999) [7], was setup with the aim of identifying Arabidopsis functions that might potentially be involved in the DNA metabolism. Such functions, when expressed in bacteria, are prone to disturb replication and thus trigger the SOS response. Consistently, expression of AtRAD51 and AtDMC1 induced the SOS response in bacteria, even affecting E. coli viability. 100 SOS-inducing cDNAs were isolated from a cDNA library constructed from an Arabidopsis cell suspension that was found to highly express meiotic genes. A large proportion of these SOS(+) candidates are clearly related to the DNA metabolism, others could be involved in the RNA metabolism, while the remaining cDNAs encode either totally unknown proteins or proteins that were considered as irrelevant. Seven SOS(+) candidate genes are induced following gamma irradiation. The in planta function of several of the SOS-inducing clones was investigated using T-DNA insertional mutants or RNA interference. Only one SOS(+) candidate, among those examined, exhibited a defined phenotype: silenced plants for DUT1 were sensitive to 5-fluoro-uracil (5FU), as is the case of the leaky dut-1 mutant in E. coli that are affected in dUTPase activity. dUTPase is essential to prevent uracil incorporation in the course of DNA replication.

The functions of DNA methylation by CcrM in Caulobacter crescentus: a global approach.

DNA methylation is involved in a diversity of processes in bacteria, including maintenance of genome integrity and regulation of gene expression. Here, using Caulobacter crescentus as a model, we exploit genome-wide experimental methods to uncover the functions of CcrM, a DNA methyltransferase conserved in most Alphaproteobacteria. Using single molecule sequencing, we provide evidence that most CcrM target motifs (GANTC) switch from a fully methylated to a hemi-methylated state when they are replicated, and back to a fully methylated state at the onset of cell division. We show that DNA methylation by CcrM is not required for the control of the initiation of chromosome replication or for DNA mismatch repair. By contrast, our transcriptome analysis shows that >10% of the genes are misexpressed in cells lacking or constitutively over-expressing CcrM. Strikingly, GANTC methylation is needed for the efficient transcription of dozens of genes that are essential for cell cycle progression, in particular for DNA metabolism and cell division. Many of them are controlled by promoters methylated by CcrM and co-regulated by other global cell cycle regulators, demonstrating an extensive cross talk between DNA methylation and the complex regulatory network that controls the cell cycle of C. crescentus and, presumably, of many other Alphaproteobacteria.

Quantification of the local heartwall motion by magnetic resonance myocardial tagging.

Sophisticated magnetic resonance tagging techniques provide powerful tools for the non-invasive assessment of the local heartwall motion towards a deeper fundamental understanding of local heart function. For the extraction of motion data from the time series of magnetic resonance tagged images and for the visualization of the local heartwall motion a new image analysis procedure has been developed. New parameters have been derived which allows quantification of the motion patterns and are highly sensitive to any changes in these patterns. The new procedure has been applied for heart motion analysis in healthy volunteers and in patient collectives with different heart diseases. The achieved results are summarized and discussed.

Simulations of electrophoretic collisions of DNA knots with gel obstacles

Gel electrophoresis can be used to separate nicked circular DNA molecules of equal length but forming different knot types. At low electric fields, complex knots drift faster than simpler knots. However, at high electric field the opposite is the case and simpler knots migrate faster than more complex knots. Using Monte Carlo simulations we investigate the reasons of this reversal of relative order of electrophoretic mobility of DNA molecules forming different knot types. We observe that at high electric fields the simulated knotted molecules tend to hang over the gel fibres and require passing over a substantial energy barrier to slip over the impeding gel fibre. At low electric field the interactions of drifting molecules with the gel fibres are weak and there are no significant energy barriers that oppose the detachment of knotted molecules from transverse gel fibres.

Automated DNA extraction using the QIAsymphony platform: Estimation of DNA recovery from simulated forensic stains

The aim of this study was to evaluate the forensic protocol recently developed by Qiagen for the QIAsymphony automated DNA extraction platform. Samples containing low amounts of DNA were specifically considered, since they represent the majority of samples processed in our laboratory. The analysis of simulated blood and saliva traces showed that the highest DNA yields were obtained with the maximal elution volume available for the forensic protocol, that is 200 ml. Resulting DNA extracts were too diluted for successful DNA profiling and required a concentration. This additional step is time consuming and potentially increases inversion and contamination risks. The 200 ml DNA extracts were concentrated to 25 ml, and the DNA recovery estimated with real-time PCR as well as with the percentage of SGM Plus alleles detected. Results using our manual protocol, based on the QIAamp DNA mini kit, and the automated protocol were comparable. Further tests will be conducted to determine more precisely DNA recovery, contamination risk and PCR inhibitors removal, once a definitive procedure, allowing the concentration of DNA extracts from low yield samples, will be available for the QIAsymphony.

Numerical simulation of gel electrophoresis of DNA knots in weak and strong electric fields.

Gel electrophoresis allows one to separate knotted DNA (nicked circular) of equal length according to the knot type. At low electric fields, complex knots, being more compact, drift faster than simpler knots. Recent experiments have shown that the drift velocity dependence on the knot type is inverted when changing from low to high electric fields. We present a computer simulation on a lattice of a closed, knotted, charged DNA chain drifting in an external electric field in a topologically restricted medium. Using a Monte Carlo algorithm, the dependence of the electrophoretic migration of the DNA molecules on the knot type and on the electric field intensity is investigated. The results are in qualitative and quantitative agreement with electrophoretic experiments done under conditions of low and high electric fields.

One billion years of bZIP transcription factor evolution: conservation and change in dimerization and DNA-binding site specificity.

The genomic era has revealed that the large repertoire of observed animal phenotypes is dependent on changes in the expression patterns of a finite number of genes, which are mediated by a plethora of transcription factors (TFs) with distinct specificities. The dimerization of TFs can also increase the complexity of a genetic regulatory network manifold, by combining a small number of monomers into dimers with distinct functions. Therefore, studying the evolution of these dimerizing TFs is vital for understanding how complexity increased during animal evolution. We focus on the second largest family of dimerizing TFs, the basic-region leucine zipper (bZIP), and infer when it expanded and how bZIP DNA-binding and dimerization functions evolved during the major phases of animal evolution. Specifically, we classify the metazoan bZIPs into 19 families and confirm the ancient nature of at least 13 of these families, predating the split of the cnidaria. We observe fixation of a core dimerization network in the last common ancestor of protostomes-deuterostomes. This was followed by an expansion of the number of proteins in the network, but no major dimerization changes in interaction partners, during the emergence of vertebrates. In conclusion, the bZIPs are an excellent model with which to understand how DNA binding and protein interactions of TFs evolved during animal evolution.

Simultaneous quantification of selective serotonin reuptake inhibitors and metabolites in human plasma by liquid chromatography-electrospray mass spectrometry for therapeutic drug monitoring.

A simple and sensitive liquid chromatography-electrospray ionization mass spectrometry method was developed for the simultaneous quantification in human plasma of all selective serotonin reuptake inhibitors (citalopram, fluoxetine, fluvoxamine, paroxetine and sertraline) and their main active metabolites (desmethyl-citalopram and norfluoxetine). A stable isotope-labeled internal standard was used for each analyte to compensate for the global method variability, including extraction and ionization variations. After sample (250μl) pre-treatment with acetonitrile (500μl) to precipitate proteins, a fast solid-phase extraction procedure was performed using mixed mode Oasis MCX 96-well plate. Chromatographic separation was achieved in less than 9.0min on a XBridge C18 column (2.1×100mm; 3.5μm) using a gradient of ammonium acetate (pH 8.1; 50mM) and acetonitrile as mobile phase at a flow rate of 0.3ml/min. The method was fully validated according to Société Française des Sciences et Techniques Pharmaceutiques protocols and the latest Food and Drug Administration guidelines. Six point calibration curves were used to cover a large concentration range of 1-500ng/ml for citalopram, desmethyl-citalopram, paroxetine and sertraline, 1-1000ng/ml for fluoxetine and fluvoxamine, and 2-1000ng/ml for norfluoxetine. Good quantitative performances were achieved in terms of trueness (84.2-109.6%), repeatability (0.9-14.6%) and intermediate precision (1.8-18.0%) in the entire assay range including the lower limit of quantification. Internal standard-normalized matrix effects were lower than 13%. The accuracy profiles (total error) were mainly included in the acceptance limits of ±30% for biological samples. The method was successfully applied for routine therapeutic drug monitoring of more than 1600 patient plasma samples over 9 months. The β-expectation tolerance intervals determined during the validation phase were coherent with the results of quality control samples analyzed during routine use. This method is therefore precise and suitable both for therapeutic drug monitoring and pharmacokinetic studies in most clinical laboratories.

Helical repeat of DNA in the region of homologous pairing.

The process of DNA strand exchange during general genetic recombination is initiated within protein-stabilized synaptic filaments containing homologous regions of interacting DNA molecules. The RecA protein in bacteria and its analogs in eukaryotic organisms start this process by forming helical filamentous complexes on single-stranded or partially single-stranded DNA molecules. These complexes then progressively bind homologous double-stranded DNA molecules so that homologous regions of single- and double-stranded DNA molecules become aligned in register while presumably winding around common axis. The topological assay presented herein allows us to conclude that in synaptic complexes containing homologous single- and double-stranded DNA molecules, all three DNA strands have a helicity of approximately 19 nt per turn.

DNA topology: feeling the pulse of a topoisomerase.

The action of individual type II DNA topoisomerases has been followed in real time by observing the elastic response of single DNA molecules to sequential strand passage events. Micromanipulation methods provide a complementary approach to biochemical studies for investigating the mechanism of DNA topoisomerases.

Associations of serum EBV DNA and gammopathy with post-transplant lymphoproliferative disease.

BACKGROUND: Post-transplant lymphoproliferative disease (PTLD) is a life-threatening complication of immunosuppression following transplantation. Epstein-Barr virus (EBV) and gammopathy in serum are associated with PTLD, but these two parameters have not been evaluated in parallel for their association with PTLD. METHODS: We evaluated the incidence of EBV load positivity, gammopathy, and protein expression in sera from all PTLD patients diagnosed at our hospital during the past seven yr. Results were compared with those of a control group including matched transplanted patients who did not develop PTLD. RESULTS: Seven of 10 PTLD patients presented EBV(+) PTLD, for which five patients had detectable serum EBV DNA levels compared with none of 38 controls (RR between two groups =121, p &lt; 0.0001). Five out of 10 patients had gammopathy at PTLD diagnosis compared with 5/38 controls (RR between two groups = 6.6, p = 0.022). Additionally, protein serum analysis by high-resolution two-dimensional gel electrophoresis and image examination failed to evidence specific abnormality in patients with PTLD compared with controls. CONCLUSIONS: Our results confirm an association between EBV in sera and gammopathy with PTLD, and highlight the high specificity of the former analysis. Whether a combination of both analyses will improve the clinical detection of PTLD remains to be evaluated in a larger prospective cohort study.

Scattering of repetitive DNA sequences in the albumin and vitellogenin gene loci of Xenopus laevis.

We have analyzed middle repetitive DNA in the albumin and vitellogenin gene families of Xenopus laevis. Mapping specific repetitive DNA sequences derived from introns of the A1 vitellogenin gene reveals that these sequences are scattered within and around the four vitellogenin genes (A1, A2, B1 and B2) and the two albumin genes (74 kd and 68 kd). Three repetitive DNA elements present in the A1 vitellogenin transcriptional unit are also located in introns of the 74 kd albumin gene. This apparently random distribution of middle repetitive DNA in the two gene families suggests that the analyzed sequences are not involved in gene regulation, but rather that they might represent unstable genetic elements. This hypothesis is further supported by the finding that size polymorphism in the A1 vitellogenin gene and in the 74 kd albumin gene is correlated with the presence or absence of repetitive DNA.