Enhanced T cell transmigration across the murine liver sinusoidal endothelium is mediated by transcytosis and surface presentation of chemokines

Transmigration through the liver endothelium is a prerequisite for the homeostatic balance of intrahepatic T cells and a key regulator of inflammatory processes within the liver. Extravasation into the liver parenchyma is regulated by the distinct expression patterns of adhesion molecules and chemokines and their receptors on the lymphocyte and endothelial cell surface. In the present study, we investigated whether liver sinusoidal endothelial cells (LSEC) inhibit or support the chemokine-driven transmigration and differentially influence the transmigration of pro-inflammatory or anti-inflammatory CD4(+) T cells, indicating a mechanism of hepatic immunoregulation. Finally, the results shed light on the molecular mechanisms by which LSEC modulate chemokine-dependent transmigration. LSEC significantly enhanced the chemotactic effect of CXC-motif chemokine ligand 12 (CXCL12) and CXCL9, but not of CXCL16 or CCL20, on naive and memory CD4(+) T cells of a T helper 1, T helper 2, or interleukin-10-producing phenotype. In contrast, brain and lymphatic endothelioma cells and ex vivo isolated lung endothelia inhibited chemokine-driven transmigration. As for the molecular mechanisms, chemokine-induced activation of LSEC was excluded by blockage of G(i)-protein-coupled signaling and the use of knockout mice. After preincubation of CXCL12 to the basal side, LSEC took up CXCL12 and enhanced transmigration as efficiently as in the presence of the soluble chemokine. Blockage of transcytosis in LSEC significantly inhibited this effect, and this suggested that chemokines taken up from the basolateral side and presented on the luminal side of endothelial cells trigger T cell transmigration. CONCLUSION: Our findings demonstrate a unique capacity of LSEC to present chemokines to circulating lymphocytes and highlight the importance of endothelial cells for the in vivo effects of chemokines. Chemokine presentation by LSEC could provide a future therapeutic target for inhibiting lymphocyte immigration and suppressing hepatic inflammation.

Effect of reduced oxygen tension and long-term mechanical stimulation on chondrocyte-polymer constructs

We have investigated the influence of long-term confined dynamic compression and surface motion under low oxygen tension on tissue-engineered cell-scaffold constructs. Porous polyurethane scaffolds (8 mm x 4 mm) were seeded with bovine articular chondrocytes and cultured under normoxic (21% O(2)) or hypoxic (5% O(2)) conditions for up to 4 weeks. By means of our joint-simulating bioreactor, cyclic axial compression (10-20%; 0.5 Hz) was applied for 1 h daily with a ceramic ball, which simultaneously oscillated over the construct surface (+/-25 degrees; 0.5 Hz). Culture under reduced oxygen tension resulted in an increase in mRNA levels of type II collagen and aggrecan, whereas the expression of type I collagen was down-regulated at early time points. A higher glycosaminoglycan content was found in hypoxic than in normoxic constructs. Immunohistochemical analysis showed more intense type II and weaker type I collagen staining in hypoxic than in normoxic cultures. Type II collagen gene expression was slightly elevated after short-term loading, whereas aggrecan mRNA levels were not influenced by the applied mechanical stimuli. Of importance, the combination of loading and low oxygen tension resulted in a further down-regulation of collagen type I mRNA expression, contributing to the stabilization of the chondrocytic phenotype. Histological results confirmed the beneficial effect of mechanical loading on chondrocyte matrix synthesis. Thus, mechanical stimulation combined with low oxygen tension is an effective tool for modulating the chondrocytic phenotype and should be considered when chondrocytes or mesenchymal stem cells are cultured and differentiated with the aim of generating cartilage-like tissue in vitro.

Reanalysis-driven climate simulation over CORDEX North America domain using the Canadian Regional Climate Model, version 5: model performance evaluation

The performance of reanalysis-driven Canadian Regional Climate Model, version 5 (CRCM5) in reproducing the present climate over the North American COordinated Regional climate Downscaling EXperiment domain for the 1989–2008 period has been assessed in comparison with several observation-based datasets. The model reproduces satisfactorily the near-surface temperature and precipitation characteristics over most part of North America. Coastal and mountainous zones remain problematic: a cold bias (2–6 °C) prevails over Rocky Mountains in summertime and all year-round over Mexico; winter precipitation in mountainous coastal regions is overestimated. The precipitation patterns related to the North American Monsoon are well reproduced, except on its northern limit. The spatial and temporal structure of the Great Plains Low-Level Jet is well reproduced by the model; however, the night-time precipitation maximum in the jet area is underestimated. The performance of CRCM5 was assessed against earlier CRCM versions and other RCMs. CRCM5 is shown to have been substantially improved compared to CRCM3 and CRCM4 in terms of seasonal mean statistics, and to be comparable to other modern RCMs.

Recruitment of EB1, a master regulator of microtubule dynamics, to the surface of the Theileria annulata schizont

The apicomplexan parasite Theileria annulata transforms infected host cells, inducing uncontrolled proliferation and clonal expansion of the parasitized cell population. Shortly after sporozoite entry into the target cell, the surrounding host cell membrane is dissolved and an array of host cell microtubules (MTs) surrounds the parasite, which develops into the transforming schizont. The latter does not egress to invade and transform other cells. Instead, it remains tethered to host cell MTs and, during mitosis and cytokinesis, engages the cell's astral and central spindle MTs to secure its distribution between the two daughter cells. The molecular mechanism by which the schizont recruits and stabilizes host cell MTs is not known. MT minus ends are mostly anchored in the MT organizing center, while the plus ends explore the cellular space, switching constantly between phases of growth and shrinkage (called dynamic instability). Assuming the plus ends of growing MTs provide the first point of contact with the parasite, we focused on the complex protein machinery associated with these structures. We now report how the schizont recruits end-binding protein 1 (EB1), a central component of the MT plus end protein interaction network and key regulator of host cell MT dynamics. Using a range of in vitro experiments, we demonstrate that T. annulata p104, a polymorphic antigen expressed on the schizont surface, functions as a genuine EB1-binding protein and can recruit EB1 in the absence of any other parasite proteins. Binding strictly depends on a consensus SxIP motif located in a highly disordered C-terminal region of p104. We further show that parasite interaction with host cell EB1 is cell cycle regulated. This is the first description of a pathogen-encoded protein to interact with EB1 via a bona-fide SxIP motif. Our findings provide important new insight into the mode of interaction between Theileria and the host cell cytoskeleton.

The influence of bedrock orientation on the landscape evolution, surface morphology and denudation (10Be) at the Niesen, Switzerland

Landscape evolution and surface morphology in mountainous settings are a function of the relative importance between sediment transport processes acting on hillslopes and in channels, modulated by climate variables. The Niesen nappe in the Swiss Penninic Prealps presents a unique setting in which opposite facing flanks host basins underlain by identical lithologies, but contrasting litho-tectonic architectures where lithologies either dip parallel to the topographic slope or in the opposite direction (i.e. dip slope and non-dip slope). The north-western facing Diemtigen flank represents such a dip slope situation and is characterized by a gentle topography, low hillslope gradients, poorly dissected channels, and it hosts large landslides. In contrast, the south-eastern facing Frutigen side can be described as non-dip slope flank with deeply incised bedrock channels, high mean hillslope gradients and high relief topography. Results from morphometric analysis reveal that noticeable differences in morphometric parameters can be related to the contrasts in the relative importance of the internal hillslope-channel system between both valley flanks. While the contrasting dip-orientations of the underlying flysch bedrock has promoted hillslope and channelized processes to contrasting extents and particularly the occurrence of large landslides on the dip slope flank, the flank averaged beryllium-10 (10Be)-derived denudation rates are very similar and range between 0.20 and 0.26 mm yr−1. In addition, our denudation rates offer no direct relationship to basin's slope, area, steepness or concavity index, but reveal a positive correlation to mean basin elevation that we interpret as having been controlled by climatically driven factors such as frost-induced processes and orographic precipitation. Our findings illustrate that while the landscape properties in this part of the northern Alpine border can mainly be related to the tectonic architecture of the underlying bedrock, the denudation rates have a strong orographic control through elevation dependent mean annual temperature and precipitation.

MODELING SPORADIC TUMOR FORMATION DRIVEN BY TELOMERE DYSFUNCTION IN THE GASTROINTESTINAL TRACT

Colorectal cancer is a complex disease that is thought to arise when cells accumulate mutations that allow for uncontrolled growth. There are several recognized mechanisms for generating such mutations in sporadic colon cancer; one of which is chromosomal instability (CIN). One hypothesized driver of CIN in cancer is the improper repair of dysfunctional telomeres. Telomeres comprise the linear ends of chromosomes and play a dual role in cancer. Its length is maintained by the ribonucleoprotein, telomerase, which is not a normally expressed in somatic cells and as cells divide, telomeres continuously shorten. Critically shortened telomeres are considered dysfunctional as they are recognized as sites of DNA damage and cells respond by entering into replicative senescence or apoptosis, a process that is p53-dependent and the mechanism for telomere-induced tumor suppression. Loss of this checkpoint and improper repair of dysfunctional telomeres can initiate a cycle of fusion, bridge and breakage that can lead to chromosomal changes and genomic instability, a process that can lead to transformation of normal cells to cancer cells. Mouse models of telomere dysfunction are currently based on knocking out the telomerase protein or RNA component; however, the naturally long telomeres of mice require multiple generational crosses of telomerase null mice to achieve critically short telomeres. Shelterin is a complex of six core proteins that bind to telomeres specifically. Pot1a is a highly conserved member of this complex that specifically binds to the telomeric single-stranded 3’ G-rich overhang. Previous work in our lab has shown that Pot1a is essential for chromosomal end protection as deletion of Pot1a in murine embryonic fibroblasts (MEFs) leads to open telomere ends that initiate a DNA damage response mediated by ATR, resulting in p53-dependent cellular senescence. Loss of Pot1a in the background of p53 deficiency results in increased aberrant homologous recombination at telomeres and elevated genomic instability, which allows Pot1a-/-, p53-/- MEFs to form tumors when injected into SCID mice. These phenotypes are similar to those seen in cells with critically shortened telomeres. In this work, we created a mouse model of telomere ysfunction in the gastrointestinal tract through the conditional deletion of Pot1a that recapitulates the microscopic features seen in severe telomere attrition. Combined intestinal loss of Pot1a and p53 lead to formation of invasive adenocarcinomas in the small and large intestines. The tumors formed with long latency, low multiplicity and had complex genomes due to chromosomal instability, features similar to those seen in sporadic human colorectal cancers. Taken together, we have developed a novel mouse model of intestinal tumorigenesis based on genomic instability driven by telomere dysfunction.

Surface charge of polymer coated SPIONs influences the serum protein adsorption, colloidal stability and subsequent cell interaction in vitro

It is known that the nanoparticle-cell interaction strongly depends on the physicochemical properties of the investigated particles. In addition, medium density and viscosity influence the colloidal behaviour of nanoparticles. Here, we show how nanoparticle-protein interactions are related to the particular physicochemical characteristics of the particles, such as their colloidal stability, and how this significantly influences the subsequent nanoparticle-cell interaction in vitro. Therefore, different surface charged superparamagnetic iron oxide nanoparticles were synthesized and characterized. Similar adsorbed protein profiles were identified following incubation in supplemented cell culture media, although cellular uptake varied significantly between the different particles. However, positively charged nanoparticles displayed a significantly lower colloidal stability than neutral and negatively charged particles while showing higher non-sedimentation driven cell-internalization in vitro without any significant cytotoxic effects. The results of this study strongly indicate therefore that an understanding of the aggregation state of NPs in biological fluids is crucial in regards to their biological interaction(s).

Magma-Driven Multiple Dike Propagation and Fracture Toughness of Crustal Rocks

Dike swarms consisting of tens to thousands of subparallel dikes are commonly observed at Earth's surface, raising the possibility of simultaneous propagation of two or more dikes at various stages of a swarm's development. The behavior of multiple propagating dikes differs from that of a single dike owing to the interacting stress fields associated with each dike. We analyze an array of parallel, periodically spaced dikes that grow simultaneously from an overpressured source into a semi-infinite, linear elastic host rock. To simplify the analysis, we assume steady state (constant velocity) magma flow and dike propagation. We use a perturbation method to analyze the coupled, nonlinear problem of multiple dike propagation and magma transport. The stress intensity factor at the dike tips and the opening displacements of the dike surfaces are calculated. The numerical results show that dike spacing has a profound effect on the behavior of dike propagation. The stress intensity factors at the tips of parallel dikes decrease with a decrease in dike spacing and are significantly smaller than that for a single dike with the same length. The reduced stress intensity factor indicates that, compared to a single dike, propagation of parallel dikes is more likely to be arrested under otherwise the same conditions. It also implies that fracture toughness of the host rock in a high confining pressure environment may not be as high as inferred from the propagation of a single dike. Our numerical results suggest fracture toughness values on the order of 100 MPa root m. The opening displacements for parallel dikes are smaller than that for a single dike, which results in higher magma pressure gradients in parallel dikes and lower flux of magma transport.

Surface "Waves" on Byrd Glacier, Antarctica

Byrd Glacier has one of the largest ice catchment areas in Antarctica, delivers more ice to the Ross Ice Shelf than any other ice stream, and is the fastest of these ice streams. A force balance, combined with a mass balance, demonstrates that stream flow in Byrd Glacier is transitional from sheet flow in East Antarctica to shelf flow in the Ross Ice Shelf. The longitudinal pulling stress, calculated along an ice flowband from the force balance, is linked to variations of ice thickness, to the ratio of the basal water pressure to the ice overburden pressure where Byrd Glacier is grounded, and is reduced by an ice-shelf buttressing stress where Byrd Glacier is floating. Longitudinal tension peaks at pressure-ratio maxima in grounded ice and close to minima in the ratio of the pulling stress to the buttressing stress in floating ice. The longitudinal spacing of these tension peaks is rather uniform and, for grounded ice, the peaks occur at maxima in surface slope that have no clear relation to the bed slope. This implies that the maxima in surface slope constitute a "wave train" that is related to regular variations in ice-bed coupling, not primarily to bed topography. It is unclear whether these surface "waves" are "standing waves" or are migrating either upslope or downslope, possibly causing the grounding line to either retreat or advance. Deciding which is the case will require obtaining bed topography in the map plane, a new map of surface topography, and more sophisticated modeling that includes ice flow linked to subglacial hydrology in the map plane.

Experimental, theoretical and numerical investigation of the nonlinear micromechanical properties of bone

Aging societies suffer from an increasing incidence of bone fractures. Bone strength depends on the amount of mineral measured by clinical densitometry, but also on the micromechanical properties of the bone hierarchical organization. A good understanding has been reached for elastic properties on several length scales, but up to now there is a lack of reliable postyield data on the lower length scales. In order to be able to describe the behavior of bone at the microscale, an anisotropic elastic-viscoplastic damage model was developed using an eccentric generalized Hill criterion and nonlinear isotropic hardening. The model was implemented as a user subroutine in Abaqus and verified using single element tests. A FE simulation of microindentation in lamellar bone was finally performed show-ing that the new constitutive model can capture the main characteristics of the indentation response of bone. As the generalized Hill criterion is limited to elliptical and cylindrical yield surfaces and the correct shape for bone is not known, a new yield surface was developed that takes any convex quadratic shape. The main advantage is that in the case of material identification the shape of the yield surface does not have to be anticipated but a minimization results in the optimal shape among all convex quadrics. The generality of the formulation was demonstrated by showing its degeneration to classical yield surfaces. Also, existing yield criteria for bone at multiple length scales were converted to the quadric formulation. Then, a computational study to determine the influence of yield surface shape and damage on the in-dentation response of bone using spherical and conical tips was performed. The constitutive model was adapted to the quadric criterion and yield surface shape and critical damage were varied. They were shown to have a major impact on the indentation curves. Their influence on indentation modulus, hardness, their ratio as well as the elastic to total work ratio were found to be very well described by multilinear regressions for both tip shapes. For conical tips, indentation depth was not a significant fac-tor, while for spherical tips damage was insignificant. All inverse methods based on microindentation suffer from a lack of uniqueness of the found material properties in the case of nonlinear material behavior. Therefore, monotonic and cyclic micropillar com-pression tests in a scanning electron microscope allowing a straightforward interpretation comple-mented by microindentation and macroscopic uniaxial compression tests were performed on dry ovine bone to identify modulus, yield stress, plastic deformation, damage accumulation and failure mecha-nisms. While the elastic properties were highly consistent, the postyield deformation and failure mech-anisms differed between the two length scales. A majority of the micropillars showed a ductile behavior with strain hardening until failure by localization in a slip plane, while the macroscopic samples failed in a quasi-brittle fashion with microcracks coalescing into macroscopic failure surfaces. In agreement with a proposed rheological model, these experiments illustrate a transition from a ductile mechanical behavior of bone at the microscale to a quasi-brittle response driven by the growth of preexisting cracks along interfaces or in the vicinity of pores at the macroscale. Subsequently, a study was undertaken to quantify the topological variability of indentations in bone and examine its relationship with mechanical properties. Indentations were performed in dry human and ovine bone in axial and transverse directions and their topography measured by AFM. Statistical shape modeling of the residual imprint allowed to define a mean shape and describe the variability with 21 principal components related to imprint depth, surface curvature and roughness. The indentation profile of bone was highly consistent and free of any pile up. A few of the topological parameters, in particular depth, showed significant correlations to variations in mechanical properties, but the cor-relations were not very strong or consistent. We could thus verify that bone is rather homogeneous in its micromechanical properties and that indentation results are not strongly influenced by small de-viations from the ideal case. As the uniaxial properties measured by micropillar compression are in conflict with the current literature on bone indentation, another dissipative mechanism has to be present. The elastic-viscoplastic damage model was therefore extended to viscoelasticity. The viscoelastic properties were identified from macroscopic experiments, while the quasistatic postelastic properties were extracted from micropillar data. It was found that viscoelasticity governed by macroscale properties has very little influence on the indentation curve and results in a clear underestimation of the creep deformation. Adding viscoplasticity leads to increased creep, but hardness is still highly overestimated. It was possible to obtain a reasonable fit with experimental indentation curves for both Berkovich and spherical indenta-tion when abandoning the assumption of shear strength being governed by an isotropy condition. These results remain to be verified by independent tests probing the micromechanical strength prop-erties in tension and shear. In conclusion, in this thesis several tools were developed to describe the complex behavior of bone on the microscale and experiments were performed to identify its material properties. Micropillar com-pression highlighted a size effect in bone due to the presence of preexisting cracks and pores or inter-faces like cement lines. It was possible to get a reasonable fit between experimental indentation curves using different tips and simulations using the constitutive model and uniaxial properties measured by micropillar compression. Additional experimental work is necessary to identify the exact nature of the size effect and the mechanical role of interfaces in bone. Deciphering the micromechanical behavior of lamellar bone and its evolution with age, disease and treatment and its failure mechanisms on several length scales will help preventing fractures in the elderly in the future.

ACL reconstruction with physiological graft tension by intraoperative adjustment of the anteroposterior translation to the uninjured contralateral knee.

PURPOSE Fixation of anterior cruciate ligament (ACL) substitutes with non-physiological anteroposterior translation (APT) worsens outcome. The aim was to present a technique for physiological APT adjustment of the transplant in ACL reconstruction and its outcome at midterm. METHODS In a consecutive series of 28 patients (age 32 ± 11 years, 24 male), chronic ACL deficiency was treated by bone-patella-tendon-bone reconstruction. Transplant APT was adjusted to that of the contralateral uninjured ACL, measured 3, 6, and 12 months postoperatively using the Rolimeter. At a median follow-up of 5.3 years (3-8 years), 82% of the patients were re-evaluated with APT measurement and using IKDC-, Tegner-, Lysholm-Scores, conventional radiographs and MRI. RESULTS No differences in APT (mean ± SD) between uninjured and reconstructed knees were observed after adjustment (6 ± 1 versus 6 ± 1 mm, n.s.). Three months postoperatively, a statistically significant increase in APT (7 ± 1 mm) and a further increase at midterm (9 ± 2 mm) were observed. Patients scored "normal" or "nearly normal", respectively, in 79% (IKDC) and 4 (3-9) points (Tegner; median, range) or 89 ± 9 points (Lysholm; mean ± SD). Radiological evaluation showed no, minimal or moderate joint degeneration in 5, 20 and 75% of patients, respectively. MRI confirmed intact ACL transplants in all patients. CONCLUSION ACL reconstruction using the presented technique was considered successful, as patients did not suffer from subjective instability, radiographic analysis did not provide evidence for graft rupture at midterm. However, APT increase and occurrence of degenerative changes in reconstructed knees at the midterm might not be prevented even by restoration of a physiological APT in ACL reconstruction. The Rolimeter can be used for quick and easy intraoperative indirect control of the applied tension to the ACL transplant by measuring the APT to obtain physiological tensioning resulting in a satisfying outcome at midterm. LEVEL OF EVIDENCE IV.

3D Intervertebral Disc Localization and Segmentation from MR Images by Data-Driven Regression and Classification

In this paper we propose a new fully-automatic method for localizing and segmenting 3D intervertebral discs from MR images, where the two problems are solved in a unified data-driven regression and classification framework. We estimate the output (image displacements for localization, or fg/bg labels for segmentation) of image points by exploiting both training data and geometric constraints simultaneously. The problem is formulated in a unified objective function which is then solved globally and efficiently. We validate our method on MR images of 25 patients. Taking manually labeled data as the ground truth, our method achieves a mean localization error of 1.3 mm, a mean Dice metric of 87%, and a mean surface distance of 1.3 mm. Our method can be applied to other localization and segmentation tasks.

Localization and Segmentation of 3D Intervertebral Discs in MR Images by Data Driven Estimation

This paper addresses the problem of fully-automatic localization and segmentation of 3D intervertebral discs (IVDs) from MR images. Our method contains two steps, where we first localize the center of each IVD, and then segment IVDs by classifying image pixels around each disc center as foreground (disc) or background. The disc localization is done by estimating the image displacements from a set of randomly sampled 3D image patches to the disc center. The image displacements are estimated by jointly optimizing the training and test displacement values in a data-driven way, where we take into consideration both the training data and the geometric constraint on the test image. After the disc centers are localized, we segment the discs by classifying image pixels around disc centers as background or foreground. The classification is done in a similar data-driven approach as we used for localization, but in this segmentation case we are aiming to estimate the foreground/background probability of each pixel instead of the image displacements. In addition, an extra neighborhood smooth constraint is introduced to enforce the local smoothness of the label field. Our method is validated on 3D T2-weighted turbo spin echo MR images of 35 patients from two different studies. Experiments show that compared to state of the art, our method achieves better or comparable results. Specifically, we achieve for localization a mean error of 1.6-2.0 mm, and for segmentation a mean Dice metric of 85%-88% and a mean surface distance of 1.3-1.4 mm.

A 49 year hindcast of surface winds over the Iberian Peninsula

The use of hindcast climatic data is quite extended for multiple applications. However, this approach needs the support of a validation process to allow its drawbacks and, therefore, confidence levels to be assessed. In this work, the strategy relies on an hourly wind database resulting from a dynamical downscaling experiment, with a spatial resolution of 10 km, covering the Iberian Peninsula (IP), driven by the ERA40 reanalysis (1959–2001) extended by European Centre for Medium-Range Weather Forecast (ECMWF) analysis (2002–2007) and comprising two main steps. Initially, the skill of the simulation is evaluated comparing the quality-tested observational database (Lorente-Plazas et al., 2014) at local and regional scales. The results show that the model is able to portray the main features of the wind over the IP: annual cycles, wind roses, spatial and temporal variability, as well as the response to different circulation types. In addition, there is a significant added value of the simulation with respect to driving conditions, especially in regions with a complex orography. However, some problems are evident, the major drawback being the systematic overestimation of the wind speed, which is mainly attributed to a missrepresentation of frictional forces. The model skill is also lower along the Mediterranean coast and for the Pyrenees. In a second phase, the high spatio-temporal resolution of the pseudo-real wind database is used to explore the limitations of the observational database. It is shown that missing values do not affect the characterisation of the wind climate over the IP, while the length of the observational period (6 years) is sufficient for most regions, with only a few exceptions. The spatial distribution of the observational sampling schemes should be enhanced to improve the correct assessment of all IP wind regimes, particularly in some mountainous areas.

Numerical study of the laser-tip coupling in surface plasmon assisted stacked-double-gate field emitter arrays

Recently, sub-wavelength-pitch stacked double-gate metal nanotip arrays have been proposed to realize high current, high brightness electron bunches for ultrabright cathodes for x-ray free-electron laser applications. With the proposed device structure, ultrafast field emission of photoexcited electrons is efficiently driven by vertical incident near infrared laser pulses, via near field coupling of the surface plasmon polariton resonance of the gate electrodes with the nanotip apex. In this work, in order to gain insight in the underlying physical processes, the authors report detailed numerical studies of the proposed device. The results indicate the importance of the interaction of the double-layer surface plasmon polariton, the position of the nanotip, as well as the incident angle of the near infrared laser pulses.