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.

Novel mutation in the tyrosine kinase domain of FGFR2 in a patient with Pfeiffer syndrome.

Mutations in the fibroblast growth factor receptor 2 (FGFR2) cause a variety of craniosynostosis syndromes. The mutational spectrum tends to be narrow with the majority of mutations occurring in either exon IIIa or IIIc or in the intronic sequence preceding exon IIIc. Mutations outside of this hotspot are uncommon and the few identified mutations have demonstrated wide clinical variability, making it difficult to establish a clear-cut genotype-phenotype correlation. To better delineate the clinical picture associated with these unusual mutations, we describe a severely affected patient with Pfeiffer syndrome and a missense mutation in the tyrosine kinase (TK) domain of FGFR2.

The death domain-containing protein Unc5CL is a novel MyD88-independent activator of the pro-inflammatory IRAK signaling cascade.

The family of death domain (DD)-containing proteins are involved in many cellular processes, including apoptosis, inflammation and development. One of these molecules, the adapter protein MyD88, is a key factor in innate and adaptive immunity that integrates signals from the Toll-like receptor/interleukin (IL)-1 receptor (TLR/IL-1R) superfamily by providing an activation platform for IL-1R-associated kinases (IRAKs). Here we show that the DD-containing protein Unc5CL (also known as ZUD) is involved in a novel MyD88-independent mode of IRAK signaling that culminates in the activation of the transcription factor nuclear factor kappa B (NF-κB) and c-Jun N-terminal kinase. Unc5CL required IRAK1, IRAK4 and TNF receptor-associated factor 6 but not MyD88 for its ability to activate these pathways. Interestingly, the protein is constitutively autoproteolytically processed, and is anchored by its N-terminus specifically to the apical face of mucosal epithelial cells. Transcriptional profiling identified mainly chemokines, including IL-8, CXCL1 and CCL20 as Unc5CL target genes. Its prominent expression in mucosal tissues, as well as its ability to induce a pro-inflammatory program in cells, suggests that Unc5CL is a factor in epithelial inflammation and immunity as well as a candidate gene involved in mucosal diseases such as inflammatory bowel disease.

Pathogenic Vibrio activate NLRP3 inflammasome via cytotoxins and TLR/nucleotide-binding oligomerization domain-mediated NF-kappa B signaling.

Vibrio vulnificus and Vibrio cholerae are Gram-negative pathogens that cause serious infectious disease in humans. The beta form of pro-IL-1 is thought to be involved in inflammatory responses and disease development during infection with these pathogens, but the mechanism of beta form of pro-IL-1 production remains poorly defined. In this study, we demonstrate that infection of mouse macrophages with two pathogenic Vibrio triggers the activation of caspase-1 via the NLRP3 inflammasome. Activation of the NLRP3 inflammasome was mediated by hemolysins and multifunctional repeat-in-toxins produced by the pathogenic bacteria. NLRP3 activation in response to V. vulnificus infection required NF-kappaB activation, which was mediated via TLR signaling. V. cholerae-induced NLRP3 activation also required NF-kappaB activation but was independent of TLR stimulation. Studies with purified V. cholerae hemolysin revealed that toxin-stimulated NLRP3 activation was induced by TLR and nucleotide-binding oligomerization domain 1/2 ligand-mediated NF-kappaB activation. Our results identify the NLRP3 inflammasome as a sensor of Vibrio infections through the action of bacterial cytotoxins and differential activation of innate signaling pathways acting upstream of NF-kappaB.

Recurrent dominant mutations affecting two adjacent residues in the motor domain of the monomeric kinesin KIF22 result in skeletal dysplasia and joint laxity.

Spondyloepimetaphyseal dysplasia with joint laxity, leptodactylic type (lepto-SEMDJL, aka SEMDJL, Hall type), is an autosomal dominant skeletal disorder that, in spite of being relatively common among skeletal dysplasias, has eluded molecular elucidation so far. We used whole-exome sequencing of five unrelated individuals with lepto-SEMDJL to identify mutations in KIF22 as the cause of this skeletal condition. Missense mutations affecting one of two adjacent amino acids in the motor domain of KIF22 were present in 20 familial cases from eight families and in 12 other sporadic cases. The skeletal and connective tissue phenotype produced by these specific mutations point to functions of KIF22 beyond those previously ascribed functions involving chromosome segregation. Although we have found Kif22 to be strongly upregulated at the growth plate, the precise pathogenetic mechanisms remain to be elucidated.

ZNRD1 (zinc ribbon domain-containing 1) is a host cellular factor that influences HIV-1 replication and disease progression.

BACKGROUND: Human immunodeficiency virus (HIV) takes advantage of multiple host proteins to support its own replication. The gene ZNRD1 (zinc ribbon domain-containing 1) has been identified as encoding a potential host factor that influenced disease progression in HIV-positive individuals in a genomewide association study and also significantly affected HIV replication in a large-scale in vitro short interfering RNA (siRNA) screen. Genes and polymorphisms identified by large-scale analysis need to be followed up by means of functional assays and resequencing efforts to more precisely map causal genes. METHODS: Genotyping and ZNRD1 gene resequencing for 208 HIV-positive subjects (119 who experienced long-term nonprogression [LTNP] and 89 who experienced normal disease progression) was done by either TaqMan genotyping assays or direct sequencing. Genetic association analysis was performed with the SNPassoc package and Haploview software. siRNA and short hairpin RNA (shRNA) specifically targeting ZNRD1 were used to transiently or stably down-regulate ZNRD1 expression in both lymphoid and nonlymphoid cells. Cells were infected with X4 and R5 HIV strains, and efficiency of infection was assessed by reporter gene assay or p24 assay. RESULTS: Genetic association analysis found a strong statistically significant correlation with the LTNP phenotype (single-nucleotide polymorphism rs1048412; [Formula: see text]), independently of HLA-A10 influence. siRNA-based functional analysis showed that ZNRD1 down-regulation by siRNA or shRNA impaired HIV-1 replication at the transcription level in both lymphoid and nonlymphoid cells. CONCLUSION: Genetic association analysis unequivocally identified ZNRD1 as an independent marker of LTNP to AIDS. Moreover, in vitro experiments pointed to viral transcription as the inhibited step. Thus, our data strongly suggest that ZNRD1 is a host cellular factor that influences HIV-1 replication and disease progression in HIV-positive individuals.

Variscan lamprophyres in the Lower Penninic domain (Central Alps): age and tectonic significance

Lamprophyre dykes have been recently discovered in blocks of gneiss embedded in a calcschist formation of wildflysch type that forms the top of the Mesozoic-Tertiary metasedimentary cover of the Antigorio nappe (the Teggiolo zone) in the Val Bavona (Lower Penninic, NW Ticino, Switzerland). The presence of the lamprophyres gives a clue to the possible source of these blocks. Similar dykes occur in the N part of the Maggia nappe where they are intruded into the Matorello granite and the surrounding gneisses. We studied these lamprophyres at two localities in the Teggiolo zone (Tamierpass and Lago del Zott) and at one locality in the Maggia nappe (Laghetti). Detailed mineralogical and geochemical investigations confirm their great similarity, particularly between the Tamier and Laghetti dykes. They all recrystallized during Alpine metamorphism under amphibolite facies conditions and lost their primary mineral assemblages and textures. The chemistry reveals a calc-alkaline affinity, a limited differentiation range, features of mineral accumulation and intense remobilization in some cases. The lamprophyres are characterized by a high mg# and relatively low contents in REE and other incompatible elements. In situ SHRIMP and LA-ICPMS U-Pb zircon dating yielded ages of 284.8 +/- 1.7 Ma (Tamier), 290.0 +/- 1.3 Ma (Zott) and 290.5 +/- 3.7 Ma (Laghetti). These ages are compatible with the general late- to post-Variscan magmatic evolution of the Helvetic and Lower Penninic domains. The lamprophyres are considered as melts derived from the lithospheric Variscan mantle, variously hybridized and differentiated at the contact with crustal material during late- to post-orogenic extension. These lamprophyres are chemically distinct from earlier lamprophyres of Visean age, emplaced together with their associated granites in transcurrent fault zones during the Variscan orogenic compression. The similarity of these different dykes suggests that the front of the Maggia nappe is a likely source of the gneissic blocks embedded in the calcschists at the top of the Teggiolo zone. They would have been provided by the advancing Maggia nappe during its thrusting over the Antigorio nappe and simultaneous closure of the Teggiolo sedimentary basin.

Crystallization and preliminary crystallographic characterization of an SH3 domain from the IB1 scaffold protein.

IB1 is a mammalian scaffold protein that interacts with components of the c-Jun N-terminal kinase (JNK) signal-transduction pathway mainly via its protein-protein interaction domains. Crystallization of the key Src homology 3 (SH3) domain of IB1 has been achieved. Crystallization experiments with unmodified protein and deliberately oxidized protein have led to different crystal forms. X-ray data have been collected to 3.0 A resolution from a crystal form with rectangular prism morphology. These crystals are orthorhombic (P2(1)2(1)2(1)), with unit-cell parameters a = 45.9, b = 57.0, c = 145.5 A. These are the first crystallographic data on a scaffold molecule such as IB1 to be reported.

A pseudo-spectral method for the simulation of poro-elastic seismic wave propagation in 2D polar coordinates using domain decomposition

We present a novel numerical approach for the comprehensive, flexible, and accurate simulation of poro-elastic wave propagation in 2D polar coordinates. An important application of this method and its extensions will be the modeling of complex seismic wave phenomena in fluid-filled boreholes, which represents a major, and as of yet largely unresolved, computational problem in exploration geophysics. In view of this, we consider a numerical mesh, which can be arbitrarily heterogeneous, consisting of two or more concentric rings representing the fluid in the center and the surrounding porous medium. The spatial discretization is based on a Chebyshev expansion in the radial direction and a Fourier expansion in the azimuthal direction and a Runge-Kutta integration scheme for the time evolution. A domain decomposition method is used to match the fluid-solid boundary conditions based on the method of characteristics. This multi-domain approach allows for significant reductions of the number of grid points in the azimuthal direction for the inner grid domain and thus for corresponding increases of the time step and enhancements of computational efficiency. The viability and accuracy of the proposed method has been rigorously tested and verified through comparisons with analytical solutions as well as with the results obtained with a corresponding, previously published, and independently bench-marked solution for 2D Cartesian coordinates. Finally, the proposed numerical solution also satisfies the reciprocity theorem, which indicates that the inherent singularity associated with the origin of the polar coordinate system is adequately handled.

The SM protein Vps33 and the t-SNARE H(abc) domain promote fusion pore opening.

Intracellular membrane fusion proceeds via distinct stages of membrane docking, hemifusion and fusion pore opening and depends on interacting families of Rab, SNARE and SM proteins. Trans-SNARE complexes dock the membranes in close apposition. Efficient fusion requires further SNARE-associated proteins. They might increase the number of trans-SNARE complexes or the fusogenic potential of a single SNARE complex. We investigated the contributions of the SM protein Vps33 to hemifusion and pore opening between yeast vacuoles. Mutations in Vps33 that weaken its interactions with the SNARE complex allowed normal trans-SNARE pairing and lipid mixing but retarded content mixing. Deleting the H(abc) domain of the vacuolar t-SNARE Vam3, which interacts with Vps33, had the same effect. This suggests that SM proteins promote fusion pore opening by enhancing the fusogenic activity of a SNARE complex. They should thus be considered integral parts of the fusion machinery.

Homologous DNA pairing domain peptides of RecA protein: intrinsic propensity to form beta-structures and filaments.

The 20 amino acid residue peptides derived from RecA loop L2 have been shown to be the pairing domain of RecA. The peptides bind to ss- and dsDNA, unstack ssDNA, and pair the ssDNA to its homologous target in a duplex DNA. As shown by circular dichroism, upon binding to DNA the disordered peptides adopt a beta-structure conformation. Here we show that the conformational change of the peptide from random coil to beta-structure is important in binding ss- and dsDNA. The beta-structure in the DNA pairing peptides can be induced by many environmental conditions such as high pH, high concentration, and non-micellar sodium dodecyl sulfate (6 mM). This behavior indicates an intrinsic property of these peptides to form a beta-structure. A beta-structure model for the loop L2 of RecA protein when bound to DNA is thus proposed. The fact that aromatic residues at the central position 203 strongly modulate the peptide binding to DNA and subsequent biochemical activities can be accounted for by the direct effect of the aromatic amino acids on the peptide conformational change. The DNA-pairing domain of RecA visualized by electron microscopy self-assembles into a filamentous structure like RecA. The relevance of such a peptide filamentous structure to the structure of RecA when bound to DNA is discussed.

The alternatively spliced domain TnFnIII A1A2 of the extracellular matrix protein tenascin-C suppresses activation-induced T lymphocyte proliferation and cytokine production.

Several lines of evidences have suggested that T cell activation could be impaired in the tumor environment, a condition referred to as tumor-induced immunosuppression. We have previously shown that tenascin-C, an extracellular matrix protein highly expressed in the tumor stroma, inhibits T lymphocyte activation in vitro, raising the possibility that this molecule might contribute to tumor-induced immunosuppression in vivo. However, the region of the protein mediating this effect has remained elusive. Here we report the identification of the minimal region of tenascin-C that can inhibit T cell activation. Recombinant fragments corresponding to defined regions of the molecule were tested for their ability to inhibit in vitro activation of human peripheral blood T cells induced by anti-CD3 mAbs in combination with fibronectin or IL-2. A recombinant protein encompassing the alternatively spliced fibronectin type III domains of tenascin-C (TnFnIII A-D) vigorously inhibited both early and late lymphocyte activation events including activation-induced TCR/CD8 down-modulation, cytokine production, and DNA synthesis. In agreement with this, full length recombinant tenascin-C containing the alternatively spliced region suppressed T cell activation, whereas tenascin-C lacking this region did not. Using a series of smaller fragments and deletion mutants issued from this region, we have identified the TnFnIII A1A2 domain as the minimal region suppressing T cell activation. Single TnFnIII A1 or A2 domains were no longer inhibitory, while maximal inhibition required the presence of the TnFnIII A3 domain. Altogether, these data demonstrate that the TnFnIII A1A2 domain mediate the ability of tenascin-C to inhibit in vitro T cell activation and provide insights into the immunosuppressive activity of tenascin-C in vivo.

PIDD death-domain phosphorylation by ATM controls prodeath versus prosurvival PIDDosome signaling.

Biochemical evidence implicates the death-domain (DD) protein PIDD as a molecular switch capable of signaling cell survival or death in response to genotoxic stress. PIDD activity is determined by binding-partner selection at its DD: whereas recruitment of RIP1 triggers prosurvival NF-κB signaling, recruitment of RAIDD activates proapoptotic caspase-2 via PIDDosome formation. However, it remains unclear how interactor selection, and thus fate decision, is regulated at the PIDD platform. We show that the PIDDosome functions in the "Chk1-suppressed" apoptotic response to DNA damage, a conserved ATM/ATR-caspase-2 pathway antagonized by Chk1. In this pathway, ATM phosphorylates PIDD on Thr788 within the DD. This phosphorylation is necessary and sufficient for RAIDD binding and caspase-2 activation. Conversely, nonphosphorylatable PIDD fails to bind RAIDD or activate caspase-2, and engages prosurvival RIP1 instead. Thus, ATM phosphorylation of the PIDD DD enables a binary switch through which cells elect to survive or die upon DNA injury.