Evaluation of the reconstruction limits of a frequency-independent crosshole georadar waveform inversion scheme in the presence of dispersion

Waveform tomographic imaging of crosshole georadar data is a powerful method to investigate the shallow subsurface because of its ability to provide images of pertinent petrophysical parameters with extremely high spatial resolution. All current crosshole georadar waveform inversion strategies are based on the assumption of frequency-independent electromagnetic constitutive parameters. However, in reality, these parameters are known to be frequency-dependent and complex and thus recorded georadar data may show significant dispersive behavior. In this paper, we evaluate synthetically the reconstruction limits of a recently published crosshole georadar waveform inversion scheme in the presence of varying degrees of dielectric dispersion. Our results indicate that, when combined with a source wavelet estimation procedure that provides a means of partially accounting for the frequency-dependent effects through an "effective" wavelet, the inversion algorithm performs remarkably well in weakly to moderately dispersive environments and has the ability to provide adequate tomographic reconstructions.

Decision analysis for the genotype designation in low-template-DNA profiles

What genotype should the scientist specify for conducting a database search to try to find the source of a low-template-DNA (lt-DNA) trace? When the scientist answers this question, he or she makes a decision. Here, we approach this decision problem from a normative point of view by defining a decision-theoretic framework for answering this question for one locus. This framework combines the probability distribution describing the uncertainty over the trace's donor's possible genotypes with a loss function describing the scientist's preferences concerning false exclusions and false inclusions that may result from the database search. According to this approach, the scientist should choose the genotype designation that minimizes the expected loss. To illustrate the results produced by this approach, we apply it to two hypothetical cases: (1) the case of observing one peak for allele xi on a single electropherogram, and (2) the case of observing one peak for allele xi on one replicate, and a pair of peaks for alleles xi and xj, i ≠ j, on a second replicate. Given that the probabilities of allele drop-out are defined as functions of the observed peak heights, the threshold values marking the turning points when the scientist should switch from one designation to another are derived in terms of the observed peak heights. For each case, sensitivity analyses show the impact of the model's parameters on these threshold values. The results support the conclusion that the procedure should not focus on a single threshold value for making this decision for all alleles, all loci and in all laboratories.

Reverse-transcriptase-associated RNaseH activity mediates template switching during reverse transcription in vitro.

During the first steps of reverse transcription of the retroviral genome, sequences present at the extremities of the RNA are used to reconstitute a host cell PolII promoter. The assembly of the promoter occurs by template switching, which takes advantage of a direct repeat at the ends of the RNA molecule. These steps are catalysed by the viral reverse transcriptase, which carries an intrinsic RNaseH activity that is probably also involved therein. To study the role of the RNaseH activity in this first template-switching event, an in vitro system has been developed based on primer extensions of synthetic RNAs. When an RNA was reverse transcribed with wild-type reverse transcriptase in the presence of a second RNA the 3' part of which was repeated at the 5' end of the first one, extension products could be observed corresponding to a chimeric cDNA comprising both RNA species. This template switching could not be detected when a mutant reverse transcriptase lacking the RNaseH activity was used. The results show that the RNaseH activity is needed to remove the 5' RNA sequences from the cDNA:RNA hybrid thereby enabling its translocation to another RNA containing an appropriate complementary target sequence.

Rapid Site-Directed Mutagenesis Using Two-PCR-Generated DNA Fragments Reproducing the Plasmid Template.

We describe a new rapid and efficient polymerase chain reaction (PCR)-based site-directed mutagenesis method. This procedure is effective with any plasmid and it employs four oligonucleotide primers. One primer contains the desired mutation, the second is oriented in the opposite direction (one of these two primers should be phosphorylated), and the third and fourth should be coding in complementary fashion for a unique restriction site to be introduced in a nonessential region. The method consists of two simultaneous PCR reactions; the PCR products are digested with the enzyme that recognizes the newly introduced unique restriction site and then ligased and used to transform competent bacteria. Additionally, the use of Dpn I facilitates the elimination of template DNA. The newly introduced restriction site is essential for ligation in the correct orientation of the two-PCR products and is further used for mutant screening. Resulting plasmids carry both the new restriction site and the desired mutation. Using this method, more than 20 mutants have already been generated (using two different kinds of templates); all these mutants were sequenced for the desired mutation and transfected into AtT-20 cells and the expressed mutant proteins encoded by the vector were assayed.

Inversion of crosshole seismic data in heterogenous environments : Comparison of waveform and ray-based approaches

High-resolution tomographic imaging of the shallow subsurface is becoming increasingly important for a wide range of environmental, hydrological and engineering applications. Because of their superior resolution power, their sensitivity to pertinent petrophysical parameters, and their far reaching complementarities, both seismic and georadar crosshole imaging are of particular importance. To date, corresponding approaches have largely relied on asymptotic, ray-based approaches, which only account for a very small part of the observed wavefields, inherently suffer from a limited resolution, and in complex environments may prove to be inadequate. These problems can potentially be alleviated through waveform inversion. We have developed an acoustic waveform inversion approach for crosshole seismic data whose kernel is based on a finite-difference time-domain (FDTD) solution of the 2-D acoustic wave equations. This algorithm is tested on and applied to synthetic data from seismic velocity models of increasing complexity and realism and the results are compared to those obtained using state-of-the-art ray-based traveltime tomography. Regardless of the heterogeneity of the underlying models, the waveform inversion approach has the potential of reliably resolving both the geometry and the acoustic properties of features of the size of less than half a dominant wavelength. Our results do, however, also indicate that, within their inherent resolution limits, ray-based approaches provide an effective and efficient means to obtain satisfactory tomographic reconstructions of the seismic velocity structure in the presence of mild to moderate heterogeneity and in absence of strong scattering. Conversely, the excess effort of waveform inversion provides the greatest benefits for the most heterogeneous, and arguably most realistic, environments where multiple scattering effects tend to be prevalent and ray-based methods lose most of their effectiveness.

The apical localization of transcribing RNA polymerases on supercoiled DNA prevents their rotation around the template.

The interaction of Escherichia coli RNA polymerase with supercoiled DNA was visualized by cryo-electron microscopy of vitrified samples and by classical electron microscopy methods. We observed that when E. coli RNA polymerase binds to a promoter on supercoiled DNA, this promoter becomes located at an apical loop of the interwound DNA molecule. During transcription RNA polymerase shifts the apical loop along the DNA, always remaining at the top of the moving loop. This relationship between RNA polymerase and the supercoiled template precludes circling of the RNA polymerase around the DNA and prevents the growing RNA transcript from becoming entangled with the template DNA.

Source wavelet estimation for waveform inversion of crosshole georadar data

AbstractFor a wide range of environmental, hydrological, and engineering applications there is a fast growing need for high-resolution imaging. In this context, waveform tomographic imaging of crosshole georadar data is a powerful method able to provide images of pertinent electrical properties in near-surface environments with unprecedented spatial resolution. In contrast, conventional ray-based tomographic methods, which consider only a very limited part of the recorded signal (first-arrival traveltimes and maximum first-cycle amplitudes), suffer from inherent limitations in resolution and may prove to be inadequate in complex environments. For a typical crosshole georadar survey the potential improvement in resolution when using waveform-based approaches instead of ray-based approaches is in the range of one order-of- magnitude. Moreover, the spatial resolution of waveform-based inversions is comparable to that of common logging methods. While in exploration seismology waveform tomographic imaging has become well established over the past two decades, it is comparably still underdeveloped in the georadar domain despite corresponding needs. Recently, different groups have presented finite-difference time-domain waveform inversion schemes for crosshole georadar data, which are adaptations and extensions of Tarantola's seminal nonlinear generalized least-squares approach developed for the seismic case. First applications of these new crosshole georadar waveform inversion schemes on synthetic and field data have shown promising results. However, there is little known about the limits and performance of such schemes in complex environments. To this end, the general motivation of my thesis is the evaluation of the robustness and limitations of waveform inversion algorithms for crosshole georadar data in order to apply such schemes to a wide range of real world problems.One crucial issue to making applicable and effective any waveform scheme to real-world crosshole georadar problems is the accurate estimation of the source wavelet, which is unknown in reality. Waveform inversion schemes for crosshole georadar data require forward simulations of the wavefield in order to iteratively solve the inverse problem. Therefore, accurate knowledge of the source wavelet is critically important for successful application of such schemes. Relatively small differences in the estimated source wavelet shape can lead to large differences in the resulting tomograms. In the first part of my thesis, I explore the viability and robustness of a relatively simple iterative deconvolution technique that incorporates the estimation of the source wavelet into the waveform inversion procedure rather than adding additional model parameters into the inversion problem. Extensive tests indicate that this source wavelet estimation technique is simple yet effective, and is able to provide remarkably accurate and robust estimates of the source wavelet in the presence of strong heterogeneity in both the dielectric permittivity and electrical conductivity as well as significant ambient noise in the recorded data. Furthermore, our tests also indicate that the approach is insensitive to the phase characteristics of the starting wavelet, which is not the case when directly incorporating the wavelet estimation into the inverse problem.Another critical issue with crosshole georadar waveform inversion schemes which clearly needs to be investigated is the consequence of the common assumption of frequency- independent electromagnetic constitutive parameters. This is crucial since in reality, these parameters are known to be frequency-dependent and complex and thus recorded georadar data may show significant dispersive behaviour. In particular, in the presence of water, there is a wide body of evidence showing that the dielectric permittivity can be significantly frequency dependent over the GPR frequency range, due to a variety of relaxation processes. The second part of my thesis is therefore dedicated to the evaluation of the reconstruction limits of a non-dispersive crosshole georadar waveform inversion scheme in the presence of varying degrees of dielectric dispersion. I show that the inversion algorithm, combined with the iterative deconvolution-based source wavelet estimation procedure that is partially able to account for the frequency-dependent effects through an "effective" wavelet, performs remarkably well in weakly to moderately dispersive environments and has the ability to provide adequate tomographic reconstructions.

Waveform inversion of crosshole georadar data: influence of source wavelet variability and the suitability of a single wavelet assumption

Waveform-based tomographic imaging of crosshole georadar data is a powerful method to investigate the shallow subsurface because of its ability to provide images of electrical properties in near-surface environments with unprecedented spatial resolution. A critical issue with waveform inversion is the a priori unknown source signal. Indeed, the estimation of the source pulse is notoriously difficult but essential for the effective application of this method. Here, we explore the viability and robustness of a recently proposed deconvolution-based procedure to estimate the source pulse during waveform inversion of crosshole georadar data, where changes in wavelet shape with location as a result of varying near-field conditions and differences in antenna coupling may be significant. Specifically, we examine whether a single, average estimated source current function can adequately represent the pulses radiated at all transmitter locations during a crosshole georadar survey, or whether a separate source wavelet estimation should be performed for each transmitter gather. Tests with synthetic and field data indicate that remarkably good tomographic reconstructions can be obtained using a single estimated source pulse when moderate to strong variability exists in the true source signal with antenna location. Only in the case of very strong variability in the true source pulse are tomographic reconstructions clearly improved by estimating a different source wavelet for each transmitter location.

A template-assembled synthetic protein surface mimetic of the von Willebrand factor A1 domain inhibits botrocetin-induced platelet aggregation.

Platelet adhesion, the initial step of platelet activation, is mediated by the interaction of von Willebrand factor (VWF) with its platelet receptor, the GPIb-IX complex. The binding of VWF to GPIb-IX is induced either by increased shear stress or by exogenous modulators, such as botrocetin. At a molecular level, this interaction takes place between the A1 domain of VWF and the GPIb alpha chain of the GPIb-IX complex. We report here the design and functional characteristics of a VWF template-assembled synthetic protein (TASP), a chimeric four-helix-bundle TASP scaffold mimicking the surface of the A1 domain. Twelve residues located on helices alpha 3 and alpha 4 in the native A1 domain were grafted onto a surface formed by two neighboring helices of the TASP. VWF TASP was found to inhibit specifically botrocetin-induced platelet aggregation and to bind both botrocetin and GPIb alpha. However, in contrast to the native A1 domain, VWF TASP did not bind simultaneously to both ligands. Modeling studies revealed that the relative orientation of the alpha helices in VWF TASP led to a clash of bound botrocetin and GPIb alpha. These results demonstrate that a chimeric four-helix-bundle TASP as a scaffold offers a suitable surface for presenting crucial residues of the VWF A1 domain; the potential of the TASP approach for de novo protein design and mimicry is thereby illustrated.

Intracranial hypertension: what additional information can be derived from ICP waveform after head injury?

OBJECTIVE: Although intracranial hypertension is one of the important prognostic factors after head injury, increased intracranial pressure (ICP) may also be observed in patients with favourable outcome. We have studied whether the value of ICP monitoring can be augmented by indices describing cerebrovascular pressure-reactivity and pressure-volume compensatory reserve derived from ICP and arterial blood pressure (ABP) waveforms. METHOD: 96 patients with intracranial hypertension were studied retrospectively: 57 with fatal outcome and 39 with favourable outcome. ABP and ICP waveforms were recorded. Indices of cerebrovascular reactivity (PRx) and cerebrospinal compensatory reserve (RAP) were calculated as moving correlation coefficients between slow waves of ABP and ICP, and between slow waves of ICP pulse amplitude and mean ICP, respectively. The magnitude of 'slow waves' was derived using ICP low-pass spectral filtration. RESULTS: The most significant difference was found in the magnitude of slow waves that was persistently higher in patients with a favourable outcome (p<0.00004). In patients who died ICP was significantly higher (p<0.0001) and cerebrovascular pressure-reactivity (described by PRx) was compromised (p<0.024). In the same patients, pressure-volume compensatory reserve showed a gradual deterioration over time with a sudden drop of RAP when ICP started to rise, suggesting an overlapping disruption of the vasomotor response. CONCLUSION: Indices derived from ICP waveform analysis can be helpful for the interpretation of progressive intracranial hypertension in patients after brain trauma.

Validity of an arterial pressure waveform analysis device: does the puncture site play a role in the agreement with intermittent pulmonary artery catheter thermodilution measurements?

OBJECTIVE: The measurement of cardiac output is a key element in the assessment of cardiac function. Recently, a pulse contour analysis-based device without need for calibration became available (FloTrac/Vigileo, Edwards Lifescience, Irvine, CA). This study was conducted to determine if there is an impact of the arterial catheter site and to investigate the accuracy of this system when compared with the pulmonary artery catheter using the bolus thermodilution technique (PAC). DESIGN: Prospective study. SETTING: The operating room of 1 university hospital. PARTICIPANTS: Twenty patients undergoing cardiac surgery. INTERVENTIONS: CO was determined in parallel by the use of the Flotrac/Vigileo systems in the radial and femoral position (CO_rad and CO_fem) and by PAC as the reference method. Data triplets were recorded at defined time points. The primary endpoint was the comparison of CO_rad and CO_fem, and the secondary endpoint was the comparison with the PAC. MEASUREMENTS AND MAIN RESULTS: Seventy-eight simultaneous data recordings were obtained. The Bland-Altman analysis for CO_fem and CO_rad showed a bias of 0.46 L/min, precision was 0.85 L/min, and the percentage error was 34%. The Bland-Altman analysis for CO_rad and PAC showed a bias of -0.35 L/min, the precision was 1.88 L/min, and the percentage error was 76%. The Bland-Altman analysis for CO_fem and PAC showed a bias of 0.11 L/min, the precision was 1.8 L/min, and the percentage error was 69%. CONCLUSION: The FloTrac/Vigileo system was shown to not produce exactly the same CO data when used in radial and femoral arteries, even though the percentage error was close to the clinically acceptable range. Thus, the impact of the introduction site of the arterial catheter is not negligible. The agreement with thermodilution was low.

The immune space: a concept and template for rationalizing vaccine development.

Abstract Empirical testing of candidate vaccines has led to the successful development of a number of lifesaving vaccines. The advent of new tools to manipulate antigens and new methods and vectors for vaccine delivery has led to a veritable explosion of potential vaccine designs. As a result, selection of candidate vaccines suitable for large-scale efficacy testing has become more challenging. This is especially true for diseases such as dengue, HIV, and tuberculosis where there is no validated animal model or correlate of immune protection. Establishing guidelines for the selection of vaccine candidates for advanced testing has become a necessity. A number of factors could be considered in making these decisions, including, for example, safety in animal and human studies, immune profile, protection in animal studies, production processes with product quality and stability, availability of resources, and estimated cost of goods. The "immune space template" proposed here provides a standardized approach by which the quality, level, and durability of immune responses elicited in early human trials by a candidate vaccine can be described. The immune response profile will demonstrate if and how the candidate is unique relative to other candidates, especially those that have preceded it into efficacy testing and, thus, what new information concerning potential immune correlates could be learned from an efficacy trial. A thorough characterization of immune responses should also provide insight into a developer's rationale for the vaccine's proposed mechanism of action. HIV vaccine researchers plan to include this general approach in up-selecting candidates for the next large efficacy trial. This "immune space" approach may also be applicable to other vaccine development endeavors where correlates of vaccine-induced immune protection remain unknown.

Analysis of an iterative deconvolution approach for estimating the source wavelet during waveform inversion of crosshole ground-penetrating radar data

A major issue in the application of waveform inversion methods to crosshole georadar data is the accurate estimation of the source wavelet. Here, we explore the viability and robustness of incorporating this step into a time-domain waveform inversion procedure through an iterative deconvolution approach. Our results indicate that, at least in non-dispersive electrical environments, such an approach provides remarkably accurate and robust estimates of the source wavelet even in the presence of strong heterogeneity in both the dielectric permittivity and electrical conductivity. Our results also indicate that the proposed source wavelet estimation approach is relatively insensitive to ambient noise and to the phase characteristics of the starting wavelet. Finally, there appears to be little-to-no trade-off between the wavelet estimation and the tomographic imaging procedures.

Source-wavelet estimation during full-waveform inversion of ground-penetrating radar data

A major issue in the application of waveform inversion methods to crosshole ground-penetrating radar (GPR) data is the accurate estimation of the source wavelet. Here, we explore the viability and robustness of incorporating this step into a recently published time-domain inversion procedure through an iterative deconvolution approach. Our results indicate that, at least in non-dispersive electrical environments, such an approach provides remarkably accurate and robust estimates of the source wavelet even in the presence of strong heterogeneity of both the dielectric permittivity and electrical conductivity. Our results also indicate that the proposed source wavelet estimation approach is relatively insensitive to ambient noise and to the phase characteristics of the starting wavelet. Finally, there appears to be little to no trade-off between the wavelet estimation and the tomographic imaging procedures.

Full-waveform inversion of cross-hole ground-penetrating radar data to characterize a gravel aquifer close to the Thur River, Switzerland

Cross-hole radar tomography is a useful tool for mapping shallow subsurface electrical properties viz. dielectric permittivity and electrical conductivity. Common practice is to invert cross-hole radar data with ray-based tomographic algorithms using first arrival traveltimes and first cycle amplitudes. However, the resolution of conventional standard ray-based inversion schemes for cross-hole ground-penetrating radar (GPR) is limited because only a fraction of the information contained in the radar data is used. The resolution can be improved significantly by using a full-waveform inversion that considers the entire waveform, or significant parts thereof. A recently developed 2D time-domain vectorial full-waveform crosshole radar inversion code has been modified in the present study by allowing optimized acquisition setups that reduce the acquisition time and computational costs significantly. This is achieved by minimizing the number of transmitter points and maximizing the number of receiver positions. The improved algorithm was employed to invert cross-hole GPR data acquired within a gravel aquifer (4-10 m depth) in the Thur valley, Switzerland. The simulated traces of the final model obtained by the full-waveform inversion fit the observed traces very well in the lower part of the section and reasonably well in the upper part of the section. Compared to the ray-based inversion, the results from the full-waveform inversion show significantly higher resolution images. At either side, 2.5 m distance away from the cross-hole plane, borehole logs were acquired. There is a good correspondence between the conductivity tomograms and the natural gamma logs at the boundary of the gravel layer and the underlying lacustrine clay deposits. Using existing petrophysical models, the inversion results and neutron-neutron logs are converted to porosity. Without any additional calibration, the values obtained for the converted neutron-neutron logs and permittivity results are very close and similar vertical variations can be observed. The full-waveform inversion provides in both cases additional information about the subsurface. Due to the presence of the water table and associated refracted/reflected waves, the upper traces are not well fitted and the upper 2 m in the permittivity and conductivity tomograms are not reliably reconstructed because the unsaturated zone is not incorporated into the inversion domain.