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.

VARIABILITY OF PHENOTYPIC, PROTEOMIC AND TRANSCRIPTOMIC EXPRESSION OF STAPHYLOCOCCUS AUREUS SURFACE ADHESINS

Staphylococcus aureus is a highly successful pathogen responsible of a wide variety of diseases, from minor skin infection to life-threatening sepsis or infective endocarditis, as well as food poisoning and toxic shock syndrome. This heterogeneity of infections and the ability of S. aureus to develop antibiotic-resistance to virtually any available drugs reflect its extraordinary capacity to adapt and survive in a great variety of environments. The pathogenesis of S. aureus infection involves a wide range of cell wall-associated adhesins and extracellular toxins that promote host colonization and invasion. In addition, S. aureus is extremely well equipped with regulatory systems that sense environmental conditions and respond by fine tuning the expression of metabolic and virulence determinants. Surface adhesins referred to MSCRAMMs - for Microbial Surface Component Recognizing Adherence Matrix Molecules - mediate binding to the host extracellular matrix or serum components, including fibrinogen, fibronectin, collagen and elastin, and promote tissue colonization and invasion. Major MSCRAMMs include a family of surface-attached proteins covalently bound to the cell wall peptidoglycan via a conserved LPXTG motif. Genomic analyses indicate that S. aureus contain up to 22 LPXTG surface proteins, which could potentially act individually or in synergy to promote infection. In the first part of this study we determined the range of adherence phenotypes to fibrinogen and fibronectin among 30 carriage isolates of S. aureus and compared it to the adherence phenotypes of 30 infective endocarditis and 30 blood culture isolates. Overall there were great variations in in vitro adherence, but no differences were observed between carriage and infection strains. We further determined the relation between in vitro adherence and in vivo infectivity in a rat model of experimental endocarditis, using 4 isolates that displayed either extremely low or high adherence phenotypes. Unexpectedly, no differences were observed between the in vivo infectivity of isolates that were poorly and highly adherent in vitro. We concluded that the natural variability of in vitro adherence to fibrinogen and fibronectin did not correlate with in vivo infectivity, and thus that pathogenic differences between various strains might only be expressed in in vivo conditions, but not in vitro. Therefore, considering the importance of adhesins expression for infection, direct measurement of those adhesins present on the bacterial surface were made by proteomic approach. 5 In the second series of experiments we assessed the physical presence of the LPXTG species at the staphylococcal surface, as measured at various time points during growth in different culture media. S. aureus Newman was grown in either tryptic soy broth (TSB) or in Roswell Park Memorial Institute (RPMI) culture medium, and samples were removed from early exponential growth phase to late stationary phase. Experiments were performed with mutants in the global accessory-gene regulator (agr), surface protein A (Spa) and clumping factor A (ClfA). Peptides of surface proteins were recovered by "trypsin-shaving" of live bacteria, and semi-quantitative proteomic analysis was performed by tandem liquid-chromatography and mass-spectrometry (LC-MS). We also determined in parallel the mRNA expression by microarrays analysis, as well as the phenotypic adherence of the bacteria to fibrinogen in vitro. The surface proteome was highly complex and contained numerous proteins theoretically not belonging to the bacterial envelope, including ribosomal proteins and metabolic enzymes. Sixteen of the 21 known LPXTG species were detected, but were differentially expressed. As expected, 9 known agr-regulated proteins (e.g. including Spa, FnBPA, ClfA, IsdA, IsdB, SasH, SasD, SasG and FmtB) increased up to the late exponential growth phase, and were abrogated in agr-negative mutants. However, only Spa and SasH modified their proteomic and mRNA profiles in parallel in the parent and its agr negative mutant, while all other LPXTG proteins modified their proteomic profiles independently of their mRNA. Moreover, ClfA became highly transcribed and active in in vitro fibrinogen adherence tests during late growth (24h), whereas it remained poorly detected by proteomics. Differential expression was also detected in iron-rich TSB versus iron-poor RPMI. Proteins from the iron-regulated surface determinant (isd) system, including IsdA, IsdB and IsdH were barely expressed in iron-rich TSB, whereas they increased their expression by >10 time in iron-poor RPMI. We conclude that semi-quantitative proteomic analysis of specific protein species is feasible in S. aureus and that proteomic, transcriptomic and adherence phenotypes demonstrated differential profiles in S. aureus. Furthermore, peptide signatures released by trypsin shaving suggested differential protein domain exposures in various environments, which might be relevant for antiadhesins vaccines. A comprehensive understanding of the S. aureus physiology should integrate all these approaches.

Dynamics and genetic structure of Argentine ant supercolonies in their native range.

Some introduced ant populations have an extraordinary social organization, called unicoloniality, whereby individuals mix freely within large supercolonies. We investigated whether this mode of social organization also exists in native populations of the Argentine ant Linepithema humile. Behavioral analyses revealed the presence of 11 supercolonies (width 1 to 515 m) over a 3-km transect. As in the introduced range, there was always strong aggression between but never within supercolonies. The genetic data were in perfect agreement with the behavioral tests, all nests being assigned to identical supercolonies with the different methods. There was strong genetic differentiation between supercolonies but no genetic differentiation among nests within supercolonies. We never found more than a single mitochondrial haplotype per supercolony, further supporting the view that supercolonies are closed breeding units. Genetic and chemical distances between supercolonies were positively correlated, but there were no other significant associations between geographic, genetic, chemical, and behavioral distances. A comparison of supercolonies sampled in 1999 and 2005 revealed a very high turnover, with about one-third of the supercolonies being replaced yearly. This dynamic is likely to involve strong competition between supercolonies and thus act as a potent selective force maintaining unicoloniality over evolutionary time.

Using composite estimators to improve both domain and total area estimation

In this article we propose using small area estimators to improve the estimatesof both the small and large area parameters. When the objective is to estimateparameters at both levels accurately, optimality is achieved by a mixed sampledesign of fixed and proportional allocations. In the mixed sample design, oncea sample size has been determined, one fraction of it is distributedproportionally among the different small areas while the rest is evenlydistributed among them. We use Monte Carlo simulations to assess theperformance of the direct estimator and two composite covariant-freesmall area estimators, for different sample sizes and different sampledistributions. Performance is measured in terms of Mean Squared Errors(MSE) of both small and large area parameters. It is found that the adoptionof small area composite estimators open the possibility of 1) reducingsample size when precision is given, or 2) improving precision for a givensample size.

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.

Assessing an intuitive condition for stability under a range of traffic conditions via a generalised Lu-Kumar network

We argue the importance both of developing simple sufficientconditions for the stability of general multiclass queueing networks and also of assessing such conditions under a range of assumptions on the weight of the traffic flowing between service stations. To achieve the former, we review a peak-rate stability condition and extend its range of application and for the latter, we introduce a generalisation of the Lu-Kumar network on which the stability condition may be tested for a range of traffic configurations. The peak-rate condition is close to exact when the between-station traffic is light, but degrades as this traffic increases.

The three phases of the ensemble forecasting of niche models: geographic range and shifts in climatically suitable areas of Utetheisa ornatrix (Lepidoptera, Arctiidae)

Species' geographic ranges are usually considered as basic units in macroecology and biogeography, yet it is still difficult to measure them accurately for many reasons. About 20 years ago, researchers started using local data on species' occurrences to estimate broad scale ranges, thereby establishing the niche modeling approach. However, there are still many problems in model evaluation and application, and one of the solutions is to find a consensus solution among models derived from different mathematical and statistical models for niche modeling, climatic projections and variable combination, all of which are sources of uncertainty during niche modeling. In this paper, we discuss this approach of ensemble forecasting and propose that it can be divided into three phases with increasing levels of complexity. Phase I is the simple combination of maps to achieve a consensual and hopefully conservative solution. In Phase II, differences among the maps used are described by multivariate analyses, and Phase III consists of the quantitative evaluation of the relative magnitude of uncertainties from different sources and their mapping. To illustrate these developments, we analyzed the occurrence data of the tiger moth, Utetheisa ornatrix (Lepidoptera, Arctiidae), a Neotropical moth species, and modeled its geographic range in current and future climates.

Prospective study evaluating the predictability of need for retreatment with intravitreal ranibizumab for age-related macular degeneration.

PURPOSE: To investigate the rhythm and predictability of the need for retreatment with intravitreal injections of ranibizumab for neovascular age-related macular degeneration (nAMD). METHODS: This prospective study enrolled 39 patients with treatment-naïve nAMD. After three loading doses of intravitreal ranibizumab, patients underwent an intensified follow-up for 12 months (initially weekly, then with stepwise increases to every 2 weeks and to monthly after each injection). Patients were retreated on an as-needed basis if any fluid or increased central retinal thickness (CRT) (>50μm) was found on spectral domain optical coherence tomography (OCT). Statistical analysis included patients who received at least two retreatments (five injections). RESULTS: A mean of 7.5 injections (range 0-12) were given between months 3 and 15. The mean visual acuity increased by 13.1 and 12.6 ETDRS letters at months 12 and 15 respectively. Two or more injection-retreatment intervals were found in 31 patients. The variability of their intra-individual intervals up to 14 weeks was small (SD 0-2.13 weeks), revealing a high regularity of the retreatment rhythm. The SD was correlated with the mean interval duration (r = 0.89, p < 0.001). The first interval was a good predictor of the following intervals (regression coefficient =0.81). One retreatment criterion was stable in 97 % of patients (cysts or subretinal fluid). CONCLUSION: The results of this study demonstrate a high intra-individual predictability of retreatment need with ranibizumab injections for nAMD. These findings may be helpful for developing individualized treatment plans for maintained suppression of disease activity with a minimum of injections and visits.

GPS analysis of human locomotion: further evidence for long-range correlations in stride-to-stride fluctuations of gait parameters.

During free walking, gait is automatically adjusted to provide optimal mechanical output and minimal energy expenditure; gait parameters, such as cadence, fluctuate from one stride to the next around average values. It was described that this fluctuation exhibited long-range correlations and fractal-like patterns. In addition, it was suggested that these long-range correlations disappeared if the participant followed the beep of metronome to regulate his or her pace. Until now, these fractal fluctuations were only observed for stride interval, because no technique existed to adequately analyze an extended time of free walking. The aim of the present study was to measure walking speed (WS), step frequency (SF) and step length (SL) with high accuracy (<1 cm) satellite positioning method (global positioning system or GPS) in order to detect long-range correlations in the stride-to-stride fluctuations. Eight participants walked 30 min under free and constrained (metronome) conditions. Under free walking conditions, DFA (detrended fluctuation analysis) and surrogate data tests showed that the fluctuation of WS, SL and SF exhibited a fractal pattern (i.e., scaling exponent alpha: 0.5 < alpha < 1) in a large majority of participants (7/8). Under constrained conditions (metronome), SF fluctuations became significantly anti-correlated (alpha < 0.5) in all participants. However, the scaling exponent of SL and WS was not modified. We conclude that, when the walking pace is controlled by an auditory signal, the feedback loop between the planned movement (at supraspinal level) and the sensory inputs induces a continual shifting of SF around the mean (persistent anti-correlation), but with no effect on the fluctuation dynamics of the other parameters (SL, WS).

The domain and interpretation of utility functions: An exploration

This paper proposes an exploration of the methodology of utilityfunctions that distinguishes interpretation from representation. Whilerepresentation univocally assigns numbers to the entities of the domainof utility functions, interpretation relates these entities withempirically observable objects of choice. This allows us to makeexplicit the standard interpretation of utility functions which assumesthat two objects have the same utility if and only if the individual isindifferent among them. We explore the underlying assumptions of suchan hypothesis and propose a non-standard interpretation according towhich objects of choice have a well-defined utility although individualsmay vary in the way they treat these objects in a specific context.We provide examples of such a methodological approach that may explainsome reversal of preferences and suggest possible mathematicalformulations for further research.

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.

Molecular characterization of HLA-C incompatibilities in HLA-ABDR-matched unrelated bone marrow donor-recipient pairs. Sequence of two new Cw alleles (Cw*02023 and Cw*0707) and recognition by cytotoxic T lymphocytes.

While the influence of HLA-AB and -DRB1 matching on the outcome of bone marrow transplantation (BMT) with unrelated donors is clear, the evaluation of HLA-C has been hampered by its poor serological definition. Because the low resolution of standard HLA-C typing could explain the significant number of positive cytotoxic T lymphocyte precursor frequency (CTLpf) tests found among HLA-AB-subtype, DRB1/B3/B5-subtype matched patient/donor pairs, we have identified by sequencing the incompatibilities recognized by CD8+ CTL clones obtained from such positive CTLpf tests. In most cases the target molecules were HLA-C antigens that had escaped detection by serology (e.g. Cw*1601, 1502 or 0702). Direct recognition of HLA-C by a CTL clone was demonstrated by lysis of the HLA class I-negative 721.221 cell line transfected with Cw*1601 cDNA. Because of the functional importance of Cw polymorphism, a PCR-SSO oligotyping procedure was set up allowing the resolution of 29 Cw alleles. Oligotyping of a panel of 382 individuals (including 101 patients and their 272 potential unrelated donors, 5 related donors and 4 platelet donors) allowed to determine HLA-C and HLA A-B-Cw-DRB1 allelic frequencies, as well as a number of A-Cw, B-Cw, and DRB1-Cw associations. Two new HLA-Cw alleles (Cw*02023 and Cw*0707) were identified by DNA sequencing of PCR-amplified exon 2-intron 2-exon 3 amplicons. Furthermore, we determined the degree of HLA-C compatibility in 287 matched pairs that could be formed from 73 patients and their 184 potential unrelated donors compatible for HLA-AB by serology and for HLA-DRB1/ B3/B5 by oligotyping. Cw mismatches were identified in 42.1% of these pairs, and AB-subtype oligotyping showed that 30% of these Cw-incompatible pairs were also mismatched for A or B-locus subtype. The degree of HLA-C incompatibility was strongly influenced by the linkage with B alleles and by the ABDR haplotypes. Cw alleles linked with B*4403, B*5101, B18, and B62 haplotypes were frequently mismatched. Apparently high resolution DNA typing for HLA-AB does not result in full matching at locus C. Since HLA-C polymorphism is recognized by alloreactive CTLs, such incompatibilities might be as relevant as AB-subtype mismatches in clinical transplantation.

What is what?: A simple time-domain test of long-memory vs. structural breaks

This paper proposes a new time-domain test of a process being I(d), 0 < d = 1, under the null, against the alternative of being I(0) with deterministic components subject to structural breaks at known or unknown dates, with the goal of disentangling the existing identification issue between long-memory and structural breaks. Denoting by AB(t) the different types of structural breaks in the deterministic components of a time series considered by Perron (1989), the test statistic proposed here is based on the t-ratio (or the infimum of a sequence of t-ratios) of the estimated coefficient on yt-1 in an OLS regression of ?dyt on a simple transformation of the above-mentioned deterministic components and yt-1, possibly augmented by a suitable number of lags of ?dyt to account for serial correlation in the error terms. The case where d = 1 coincides with the Perron (1989) or the Zivot and Andrews (1992) approaches if the break date is known or unknown, respectively. The statistic is labelled as the SB-FDF (Structural Break-Fractional Dickey- Fuller) test, since it is based on the same principles as the well-known Dickey-Fuller unit root test. Both its asymptotic behavior and finite sample properties are analyzed, and two empirical applications are provided.

Influence of learning on range expansion and adaptation to novel habitats.

Learning has been postulated to 'drive' evolution, but its influence on adaptive evolution in heterogeneous environments has not been formally examined. We used a spatially explicit individual-based model to study the effect of learning on the expansion and adaptation of a species to a novel habitat. Fitness was mediated by a behavioural trait (resource preference), which in turn was determined by both the genotype and learning. Our findings indicate that learning substantially increases the range of parameters under which the species expands and adapts to the novel habitat, particularly if the two habitats are separated by a sharp ecotone (rather than a gradient). However, for a broad range of parameters, learning reduces the degree of genetically-based local adaptation following the expansion and facilitates maintenance of genetic variation within local populations. Thus, in heterogeneous environments learning may facilitate evolutionary range expansions and maintenance of the potential of local populations to respond to subsequent environmental changes.