ALBERT

Description: We explore thermal convection of a fluid with a temperature-dependent viscosity in a basally heated 3-D spherical shell using linear stability analyses and numerical experiments, while considering the application of our results to terrestrial planets. The inner to outer radius ratio of the shell f assumed in the linear stability analyses is in the range of 0.11–0.88. The critical Rayleigh number R c for the onset of thermal convection decreases by two orders of magnitude as f increases from 0.11 to 0.88, when the viscosity depends sensitively on the temperature, as is the case for real mantle materials. Numerical simulations carried out in the range of f  = 0.11–0.55 show that a thermal boundary layer (TBL) develops both along the surface and bottom boundaries to induce cold and hot plumes, respectively, when f is 0.33 or larger. However, for smaller f values, a TBL develops only on the bottom boundary. Convection occurs in the stagnant-lid regime where the root mean square velocity on the surface boundary is less than 1 per cent of its maximum at depth, when the ratio of the viscosity at the surface boundary to that at the bottom boundary exceeds a threshold that depends on f . The threshold decreases from 10 6.5 at f  = 0.11 to 10 4 at f  = 0.55. If the viscosity at the base of the convecting mantle is 10 20 –10 21  Pa s, the Rayleigh number exceeds R c for Mars, Venus and the Earth, but does not for the Moon and Mercury; convection is unlikely to occur in the latter planets unless the mantle viscosity is much lower than 10 20  Pa s and/or the mantle contains a strong internal heat source.

Description: Subducting oceanic lithosphere is an example of a thin sheet-like object whose characteristic lateral dimension greatly exceeds its thickness. Here we exploit this property to derive a new hybrid boundary-integral/thin sheet (BITS) representation of subduction that combines in a single equation all the forces acting on the sheet: gravity, internal resistance to bending and stretching, and the tractions exerted by the ambient mantle. For simplicity, we limit ourselves to 2-D. We solve the BITS equations using a discrete Lagrangian approach in which the sheet is represented by a set of vertices connected by edges. Instantaneous solutions for the sinking speed of a slab attached to a trailing flat sheet obey a scaling law of the form V / V Stokes  = fct(St), where V Stokes is a characteristic Stokes sinking speed and St is the sheet's flexural stiffness. Time-dependent solutions for the evolution of the sheet's shape and thickness show that these are controlled by the viscosity ratio between the sheet and its surroundings. An important advantage of the BITS approach is the possibility of generalizing the sheet's rheology, either to a viscosity that varies along the sheet or to a non-Newtonian shear-thinning rheology.

Description: The aftershock distribution of the 1933 Sanriku-oki outer trench earthquake is estimated by using modern relocation methods and a newly developed velocity structure to examine the spatial extent of the source-fault and the possibility of a triggered interplate seismicity. In this study, we first examined the regional data quality of the 1933 earthquake based on smoked-paper records and then relocated the earthquakes by using the 3-D velocity structure and double-difference method. The improvements of hypocentre locations using these methods were confirmed by the examination of recent earthquakes that are accurately located based on ocean bottom seismometer data. The results show that the 1933 aftershocks occurred under both the outer- and inner-trench-slope regions. In the outer-trench-slope region, aftershocks are distributed in a ~280-km-long area and their depths are shallower than 50 km. Although we could not constrain the fault geometry from the hypocentre distribution, the depth distribution suggests the whole lithosphere is probably not under deviatoric tension at the time of the 1933 earthquake. The occurrence of aftershocks under the inner trench slope was also confirmed by an investigation of waveform frequency difference between outer and inner trench earthquakes as recorded at Mizusawa. The earthquakes under the inner trench slope were shallow (depth less double equals30 km) and the waveforms show a low-frequency character similar to the waveforms of recent, precisely located earthquakes in the same area. They are also located where recent activity of interplate thrust earthquakes is high. These suggest that the 1933 outer-trench-slope main shock triggered interplate earthquakes, which is an unusual case in the order of occurrence in contrast with the more common pairing of a large initial interplate shock with subsequent outer-slope earthquakes. The off-trench earthquakes are distributed about 80 km width in the trench perpendicular direction. This wide width cannot be explained from a single high-angle fault confined at a shallow depth (depth less double equals50 km). The upward motion of the 1933 tsunami waveform records observed at Sanriku coast also cannot be explained from a single high-angle west-dipping normal fault. If we consider additional fault, involvement of high-angle, east-dipping normal faults can better explain the tsunami first motion and triggering of the aftershock in a wide area under the outer trench slope. Therefore multiple off-trench normal faults may have activated during the 1933 earthquake. We also relocated recent (2001–2012) seismicity by the same method. The results show that the present seismicity in the outer-trench-slope region can be divided into several groups along the trench. Comparison of the 1933 rupture dimensions based on our aftershock relocations with the morphologies of fault scarps in the outer trench slope suggest that the rupture was limited to the region where fault scarps are largely trench parallel and cross cut the seafloor spreading fabric. These findings imply that bending geometry and structural segmentation of the incoming plate largely controls the spatial extent of the 1933 seismogenic faulting. In this shallow rupture model for this largest outer trench earthquake, triggered seismicity in the forearc and structural control of faulting represent an important deformation styles for off-trench and shallow megathrust zones.

Description: The 16S rRNA gene (16S rDNA) codes for RNA that plays a fundamental role during translation in the ribosome and is used extensively as a marker gene to establish relationships among bacteria. However, the complementary non-coding 16S rDNA (nc16S rDNA) has been ignored. An idea emerged in the course of analyzing bacterial 16S rDNA sequences in search for nucleotide composition and substitution patterns: Does the nc16S rDNA code? If so, what does it code for? More importantly: Does 16S rDNA evolution reflect its own evolution or the evolution of its counterpart nc16S rDNA? The objective of this minireview is to discuss these thoughts. nc strands often encode small RNAs (sRNAs), ancient components of gene regulation. nc16S rDNA sequences from different bacterial groups were used to search for possible matches in the Bacterial Small Regulatory RNA Database. Intriguingly, the sequence of one published sRNA obtained from Legionella pneumophila (GenBank: AE017354.1) showed high non-random similarity with nc16S rDNA corresponding in part to the V5 region especially from Legionella and relatives. While the target(s) of this sRNA is unclear at the moment, its mere existence might open up a new chapter in the use of the 16S rDNA to study relationships among bacteria.

Description: Thermotolerance of the fungus Fomes sp. EUM1 was evaluated in solid state fermentation (SSF). This thermotolerant strain improved both hyphal invasiveness (38%) and length (17%) in adverse thermal conditions exceeding 30°C and to a maximum of 40°C. In contrast, hyphal branching decreased by 46% at 45°C. The production of cellulases over corn stover increased 1.6-fold in 30°C culture conditions, xylanases increased 2.8-fold at 40°C, while laccase production improved 2.7-fold at 35°C. Maximum production of lignocellulolytic enzymes was obtained at elevated temperatures in shorter fermentation times (8–6 days), although the proteases appeared as a thermal stress response associated with a drop in lignocellulolytic activities. Novel and multiple isoenzymes of xylanase (four bands) and cellulase (six bands) were secreted in the range of 20–150 kDa during growth in adverse temperature conditions. However, only a single laccase isoenzyme (46 kDa) was detected. This is the first report describing the advantages of a thermotolerant white-rot fungus in SSF. These results have important implications for large-scale SSF, where effects of metabolic heat are detrimental to growth and enzyme production, which are severely affected by the formation of high temperature gradients.

Description: Geodetic, geologic and palaeomagnetic data reveal that Oregon (western USA) rotates clockwise at 0.3 to 1.0° Ma –1 (relative to North America) about an axis near the Idaho–Oregon–Washington border, while northeast Washington is relatively fixed. This rotation has been going on for at least 15 Ma. The Yakima fold and thrust belt (YFTB) forms the boundary between northern Oregon and central Washington where convergence of the clockwise-rotating Oregon block is apparently accommodated. North–south shortening across the YFTB has been thought to occur in a fan-like manner, increasing in rate to the west. We obtained high-accuracy, high-density geodetic GPS measurements in 2012–2014 that are used with earlier GPS measurements from the 1990s to characterize YFTB kinematics. The new results show that the deformation associated with the YFTB starts at the Blue Mountains Anticline in northern Oregon and extends north beyond the Frenchman Hills in Washington, past the epicentre of the 1872 M w 7.0 Entiat earthquake to 49°N. The north–south strain rate across the region is 2 to 3 x 10 –9 yr –1 between the volcanic arc and the eastern edge of the YFTB (241.0°E); east of there it drops to about 10 –9 yr –1 . At the eastern boundary of the YFTB, faults and earthquake activity are truncated by a north-trending, narrow zone of deformation that runs along the Pasco Basin and Moses Lake regions near 240.9°E. This zone, abutting the Department of Energy Hanford Nuclear Reservation, accommodates about 0.5 mm yr –1 of east to northeast shortening. A similar zone of N-trending transpression is seen along 239.9°E where there is a change in the strike of the Yakima folds. The modern deformation of the YFTB is about 600 km wide from south to north and internally may be controlled by pre-existing crustal structure.

Description: The Central Anatolian orogenic plateau is represented by young volcanism, rapid plateau uplift and distinctive (past and active) tectonic deformation. In this study, we consider observational data in terms of regional present-day geodynamics in the region. The residual topography of Central Anatolia was derived to define the regional isostatic conditions according to Airy isostasy and infer the potential role of ‘dynamic topography’. 2-D thermomechanical forward models for coupled mantle-lithosphere flow/deformation were conducted along an N–S directional profile through the region (e.g. northern/Pontides, interior and southern/Taurides). These models were based on seismic tomography data that provide estimates about the present-day mantle thermal structure beneath the Anatolian plate. We compare the modelling results with calculated residual topography and independent data sets of geological deformation, gravity and high surface heat flow/widespread geothermal activity. Model results suggest that there is ~1 km of mantle flow induced dynamic topography associated with the sublithospheric flow driven by the seismically inferred mantle structure. The uprising mantle may have also driven the asthenospheric source of volcanism in the north (e.g. Galatia volcanic province) and the Cappadocia volcanic province in the south while elevating the surface in the last 10 Myr. Our dynamic topography calculations emphasize the role of vertical forcing under other orogenic plateaux underlain by relatively thin crust and low-density asthenospheric mantle.

Description: In many geological systems, inversion of density stratification sets in Rayleigh–Taylor (RT) instabilities, leading to an ascent of relatively low-density materials through the high-density overburden in the form of diapirs. These diapirs often originate from dipping low-density layers. This study aims to show how the initial tilt of such source layers can control the ascent behaviour of diapirs initiated by RT instabilities. Using two-layer viscous models we produced RT instabilities in physical experiments, and investigated the effects of source-layer tilts ( β ). Our experiments suggest that these diapirs ascend with contrasting lateral spreading rates in the up and down slope directions, resulting in their axi-asymmetric geometry. However, their heads retain a circular outline on the horizontal top surface, where the upwelling axis is located away from their geometric centre in the upslope direction. In this paper, we present a series of experimental models to demonstrate the spectrum of axi-symmetric to -asymmetric geometrical transitions with increasing β . Our experiments also reveal that when β is large (〉4°) the diapirs become unstable, resulting in a continuous migration of their stems in the upslope direction. Using the volume of fluid method we ran computational fluid dynamic (CFD) simulations to study the underlying hydrodynamics of axi-asymmetric diapiric growth. The CFD simulations show that β 〉 0° conditions develop stronger flow vortices on the downslope side of an ascending diapir, leading to a pressure difference between the up- and downslope flanks. Such a differential pressure causes the diapir head to spread at a faster rate in the tilt direction. An estimate of the asymmetric spreading rates is given as a function of β . Our present study provides a fundamental understanding of the hydrodynamic flow structure responsible for the asymmetric growth of RT instabilities on tilted source layers, as applicable to a wide range of large-scale geological settings, such as sedimentary basins and subduction zones.

Description: The geomechanical analysis of a highly compartmentalized reservoir is performed to simulate the seafloor subsidence due to gas production. The available observations over the hydrocarbon reservoir consist of bathymetric surveys carried out before and at the end of a 10-yr production life. The main goal is the calibration of the reservoir compressibility c M , that is, the main geomechanical parameter controlling the surface response. Two conceptual models are considered: in one (i) c M varies only with the depth and the vertical effective stress (heterogeneity due to lithostratigraphic variability); in another (ii) c M varies also in the horizontal plane, that is, it is spatially distributed within the reservoir stratigraphic units. The latter hypothesis accounts for a possible partitioning of the reservoir due to the presence of sealing faults and thrusts that suggests the idea of a block heterogeneous system with the number of reservoir blocks equal to the number of uncertain parameters. The method applied here relies on an ensemble-based data assimilation (DA) algorithm (i.e. the ensemble smoother, ES), which incorporates the information from the bathymetric measurements into the geomechanical model response to infer and reduce the uncertainty of the parameter c M . The outcome from conceptual model (i) indicates that DA is effective in reducing the c M uncertainty. However, the maximum settlement still remains underestimated, while the areal extent of the subsidence bowl is overestimated. We demonstrate that the selection of the heterogeneous conceptual model (ii) allows to reproduce much better the observations thus removing a clear bias of the model structure. DA allows significantly reducing the c M uncertainty in the five blocks (out of the seven) characterized by large volume and large pressure decline. Conversely, the assimilation of land displacements only partially constrains the prior c M uncertainty in the reservoir blocks marginally contributing to the cumulative seafloor subsidence, that is, blocks with low pressure.

Description: This work aims to explore the ongoing tectonic activity of structures in the outermost sector of the Northern Apennines, which represents the active leading edge of the thrust belt and is dominated by compressive deformation. We have applied the Persistent Scatterer Interferometry (PSI) technique to obtain new insights into the present-day deformation pattern of the frontal area of the Northern Apennine. PSI has proved to be effective in detecting surface deformation of wide regions involved in low tectonic movements. We used 34 Envisat images in descending geometry over the period of time between 2004 and 2010, performing about 300 interferometric pairs. The analysis of the velocity maps and of the PSI time-series has allowed to observe ground deformation over the sector of the Po Plain between Piacenza and Reggio Emilia. The time-series of permanent GPS stations located in the study area, validated the results of the PSI technique, showing a good correlation with the PS time-series. The PS analysis reveals the occurrence of a well-known subsidence area on the rear of the Ferrara arc, mostly connected to the exploitation of water resources. In some instances, the PS velocity pattern reveals ground uplift (with mean velocities ranging from 1 to 2.8 mm yr –1 ) above active thrust-related anticlines of the Emilia and Ferrara folds, and part of the Pede-Apennine margin. We hypothesize a correlation between the observed uplift deformation pattern and the growth of the thrust-related anticlines. As the uplift pattern corresponds to known geological features, it can be used to constrain the seismo-tectonic setting, and a working hypothesis may involve that the active Emilia and Ferrara thrust folds would be characterized by interseismic periods possibly dominated by aseismic creep.

Description: The rhizobacterium Serratia plymuthica 4Rx13 emits the novel and unique volatile sodorifen (C 16 H 26 ), which has a polymethylated bicyclic structure. Transcriptome analysis revealed that gene SOD_c20750 (annotated as terpene cyclase) is involved in the biosynthesis of sodorifen. Here we show that this gene is located in a small cluster of four genes ( SOD_c20750 – SOD_c20780 ), and the analysis of the knockout mutants demonstrated that SOD_c20760 (annotated as methyltransferase) and SOD_c20780 (annotated as isopentenyl pyrophosphate (IPP) isomerase) are needed for the biosynthesis of sodorifen, while a sodorifen-negative phenotype was not achieved with the SOD_c20770 (annotated as deoxy-xylulose-5-phosphate (DOXP) synthase) mutant. Altogether, the function of this new gene cluster was assigned to the biosynthesis of this structurally unusual volatile compound sodorifen.

Description: Among staphylococci Staphylococcus saprophyticus is the only species that is typically uropathogenic and an important cause of urinary tract infections in young women. The amino acid D-serine occurs in relatively high concentrations in human urine and has a bacteriostatic or toxic effect on many bacteria. In uropathogenic Escherichia coli and S. saprophyticus , the amino acid regulates the expression of virulence factors and can be used as a nutrient. The ability of uropathogens to respond to or to metabolize D-serine has been suggested as a factor that enables colonization of the urinary tract. Until now nothing is known about D-serine transport in S.   saprophyticus . We generated mutants of putative transporter genes in S.   saprophyticus 7108 that show homology to the D-serine transporter cyc A of E. coli and tested them in a D-serine depletion assay to analyze the D-serine uptake rate of the cells. The mutant of SPP1070 showed a strong decrease in D-serine uptake. Therefore, SSP1070 was identified as a major D-serine transporter in S. saprophyticus 7108 and was named D-serine transporter A (DstA). D-serine caused a prolonged lag phase of S. saprophyticus in a chemically defined medium. This negative effect was dependent on the presence of DstA.

Description: In Europe, common input data types for seismic hazard evaluation include earthquake catalogues, seismic zonation models and ground motion models, all with well-constrained epistemic uncertainties. In contrast, neotectonic deformation models and their related uncertainties are rarely considered in earthquake forecasting and seismic hazard studies. In this study, for the first time in Europe, we developed a seismic hazard model based exclusively on active fault and geodynamic deformation models. We applied it to the External Dinarides, a slow-deforming fold-and-thrust belt in the Central Mediterranean. The two deformation models furnish consistent long-term earthquake rates above the M w 4.7 threshold on a latitude/longitude grid with 0.2° spacing. Results suggest that the use of deformation models is a valid alternative to empirical-statistical approaches in earthquake forecasting in slow-deforming regions of Europe. Furthermore, we show that the variability of different deformation models has a comparable effect on the peak ground motion acceleration uncertainty as do the ground motion prediction equations.

Description: An important real world application of doublet flow occurs in well design of both geothermal and hydrocarbon reservoirs. A guiding principle for fluid management of injection and extraction wells is that mass balance is commonly assumed between the injected and produced fluid. Because the doublets are considered closed loops, the injection fluid is assumed to eventually reach the producer well and all the produced fluid ideally comes from stream tubes connected to the injector of the well pair making up the doublet. We show that when an aquifer background flow occurs, doublets will rarely retain closed loops of fluid recirculation. When the far-field flow rate increases relative to the doublet's strength, the area occupied by the doublet will diminish and eventually vanishes. Alternatively, rather than using a single injector (source) and single producer (sink), a linear array of multiple injectors separated by some distance from a parallel array of producers can be used in geothermal energy projects as well as in waterflooding of hydrocarbon reservoirs. Fluid flow in such an arrangement of parallel source-sink arrays is shown to be macroscopically equivalent to that of a line doublet. Again, any far-field flow that is strong enough will breach through the line doublet, which then splits into two vortices. Apart from fundamental insight into elementary flow dynamics, our new results provide practical clues that may contribute to improve the planning and design of doublets and direct line drives commonly used for flow management of groundwater, geothermal and hydrocarbon reservoirs.

Description: Genes encoding fluoride transporters have been identified in bacterial and archaeal species. The genome sequence of the cariogenic  Streptococcus mutans  bacteria suggests the presence of a putative fluoride transporter, which is referred to as a chloride channel permease. Two homologues of this gene (GenBank locus tags SMU_1290c and SMU_1289c) reside in tandem in the genome of  S. mutans . The aim of this study was to determine whether the chloride channel permeases contribute to fluoride resistance. We constructed SMU_1290c- and SMU_1289c-knockout  S. mutans  UA159 strains. We also constructed a double-knockout strain lacking both genes. SMU_1290c or SMU_1289c was transformed into a fluoride transporter- disrupted  Escherichia coli strain. All bacterial strains were cultured under appropriate conditions with or without sodium fluoride, and fluoride resistance was evaluated. All three gene-knockout  S. mutans  strains showed lower resistance to sodium fluoride than did the wild-type strain. No significant changes in resistance to other sodium halides were recognized between the wild-type and double-knockout strains. Both SMU_1290c and SMU_1289c transformation rescued fluoride transporter-disrupted  E. coli  cell from fluoride toxicity. We conclude that the chloride channel permeases contribute to fluoride resistance in  S. mutans .

Description: Clostridium thermocellum and Thermoanaerobacterium saccharolyticum are bacteria under investigation for production of biofuels from plant biomass. Thermoanaerobacterium saccharolyticum has been engineered to produce ethanol at high yield (〉90% of theoretical) and titer (〉70 g/l). Efforts to engineer C. thermocellum have not, to date, been as successful, and efforts are underway to transfer the ethanol production pathway from T. saccharolyticum to C. thermocellum . One potential challenge in transferring metabolic pathways is the possibility of incompatible levels of nicotinamide cofactors. These cofactors (NAD + , NADH, NADP + and NADPH) and their oxidation state are important in the context of microbial redox metabolism. In this study we directly measured the concentrations and reduced oxidized ratios of these cofactors in a number of strains of C. thermocellum and T. saccharolyticum by using acid/base extraction and enzymatic assays. We found that cofactor ratios are maintained in a fairly narrow range, regardless of the metabolic network modifications considered. We have found that the ratios are similar in both organisms, which is a relevant observation in the context of transferring the T. saccharolyticum ethanol production pathway to C. thermocellum .

Description: In November 2010, intense seismic activity including 29 events with a magnitude above 5.0, started in the western part of the Gulf of Aden, where the structure of the oceanic spreading ridge is characterized by a series of N115°-trending slow-spreading segments set within an EW-trending rift. Using signals recorded by permanent and temporary networks in Djibouti and Yemen, we located 1122 earthquakes, with a magnitude ranging from 2.1 to 5.6 from 2010 November 1 to 2011 March 31. By looking in detail at the space–time distribution of the overall seismicity, and both the frequency and the moment tensor of large earthquakes, we re-examine the chronology of this episode. In addition, we also interpret the origin of the activity using high-resolution bathymetric data, as well as from observations of seafloor cable damage caused by high temperatures and lava flows. The analysis allows us to identify distinct active areas. First, we interpret that this episode is mainly related to a diking event along a specific ridge segment, located at E044°. In light of previous diking episodes in nearby subaerial rift segments, for which field constraints and both seismic and geodetic data exist, we interpret the space–time evolution of the seismicity of the first few days. Migration of earthquakes suggests initial magma ascent below the segment centre. This is followed by a southeastward dike propagation below the rift immediately followed by a northwestward dike propagation below the rift ending below the northern ridge wall. The cumulative seismic moment associated with this sequence reaches 9.1 x 10 17 Nm, and taking into account a very low seismic versus geodetic moment, we estimate a horizontal opening of ~0.58–2.9 m. The seismic activity that followed occurred through several bursts of earthquakes aligned along the segment axis, which are interpreted as short dike intrusions implying fast replenishment of the crustal magma reservoir feeding the dikes. Over the whole period, the opening is estimated to be ~1.76–8.8 m across the segment. A striking feature of this episode is that the seismicity remained confined within one individual segment, whereas the adjacent en-echelon segments were totally quiescent, suggesting that the magma supply system of one segment is disconnected from those of the neighbouring segments. Second, we identify activity induced by the first intrusion with epicentres aligned along an N035°E-trending, ~30 km long at the northwestern end of the active opening segment. This group encompasses more than seven earthquakes with magnitude larger than 5.0, and with strike-slip focal mechanisms consistent with the faults identified in the bathymetry and the structural pattern of the area. We propose that a transform fault is currently in formation which indicates an early stage of the ridge segmentation, at the locus of the trend change of the spreading ridge, which also corresponds to the boundary between a clear oceanic lithosphere and the zone of transform between continental and oceanic crust.

Description: We challenge the perspective that seismicity could contribute to polar motion by arguing quantitatively that, in first approximation and on the average, interseismic deformations can compensate for it. This point is important because what we must simulate and observe in Earth Orientation Parameter time-series over intermediate timescales of decades or centuries is the residual polar motion resulting from the two opposing processes of coseismic and interseismic deformations. In this framework, we first simulate the polar motion caused by only coseismic deformations during the longest period available of instrumental seismicity, from 1900 to present, using both the CMT and ISC-GEM catalogues. The instrumental seismicity covering a little longer than one century does not represent yet the average seismicity that we should expect on the long term. Indeed, although the simulation shows a tendency to move the Earth rotation pole towards 133°E at the average rate of 16.5 mm yr –1 , this trend is still sensitive to individual megathrust earthquakes, particularly to the 1960 Chile and 1964 Alaska earthquakes. In order to further investigate this issue, we develop a global seismicity model (GSM) that is independent from any earthquake catalogue and that describes the average seismicity along plate boundaries on the long term by combining information about present-day plate kinematics with the Anderson theory of faulting, the seismic moment conservation principle and a few other assumptions. Within this framework, we obtain a secular polar motion of 8 mm yr –1 towards 112.5°E that is comparable with that estimated from 1900 to present using the earthquake catalogues, although smaller by a factor of 2 in amplitude and different by 20° in direction. Afterwards, in order to reconcile the idea of a secular polar motion caused by earthquakes with our simplest understanding of the seismic cycle, we adapt the GSM in order to account for interseismic deformations and we use it to quantify, for the first time ever, their contribution to polar motion. Taken together, coseismic and interseismic deformations make the rotation pole wander around the north pole with maximum polar excursions of about 1 m. In particular, the rotation pole moves towards about Newfoundland when the interseismic contribution dominates over the coseismic ones (i.e. during phases of low seismicity or, equivalently, when most of the fault system associated with plate boundaries is locked). When megathrust earthquakes occur, instead, the rotation pole is suddenly shifted in an almost opposite direction, towards about 133°E.

Description: Cadmium is a widespread environmental pollutant and poses some potential risks to human health. However, the signaling events controlling cadmium toxicity are not fully understood. In this study, we examined the effect of cadmium chloride on cell viability and the intracellular nitric oxide (NO) level in yeast cells. The results showed that exposure of yeast cells to cadmium (0–100 μM) could induce cell killing with significantly increased intracellular NO levels. Morphological analysis of the nuclei with 4 ' ,6-diamidino-2-phenylindole staining and DNA strand breaks analysis showed that cadmium at 50 μM can induce cell apoptosis in yeast cells. Treatment of yeast cells with cadmium (50 μM) and the nitric oxide scavenger c-PTIO [2-(4-carboxyphenyl)-4,4,5,5-teramethylimidazoline-1-oxyl-3-oxide; 0.2 mM] showed that c-PTIO attenuated the cadmium-induced cell killing. Our findings indicated that cadmium-induced yeast cell killing is mediated by a directly increased intracellular NO level.

Description: While elasticity is a defining characteristic of the Earth's lithosphere, it is often ignored in numerical models of long-term tectonic processes in favour of a simpler viscoplastic description. Here we assess the consequences of this assumption on a well-studied geodynamic problem: the growth of normal faults at an extensional plate boundary. We conduct 2-D numerical simulations of extension in elastoplastic and viscoplastic layers using a finite difference, particle-in-cell numerical approach. Our models simulate a range of faulted layer thicknesses and extension rates, allowing us to quantify the role of elasticity on three key observables: fault-induced topography, fault rotation, and fault life span. In agreement with earlier studies, simulations carried out in elastoplastic layers produce rate-independent lithospheric flexure accompanied by rapid fault rotation and an inverse relationship between fault life span and faulted layer thickness. By contrast, models carried out with a viscoplastic lithosphere produce results that may qualitatively resemble the elastoplastic case, but depend strongly on the product of extension rate and layer viscosity U x L . When this product is high, fault growth initially generates little deformation of the footwall and hanging wall blocks, resulting in unrealistic, rigid block-offset in topography across the fault. This configuration progressively transitions into a regime where topographic decay associated with flexure is fully accommodated within the numerical domain. In addition, high U x L favours the sequential growth of multiple short-offset faults as opposed to a large-offset detachment. We interpret these results by comparing them to an analytical model for the fault-induced flexure of a thin viscous plate. The key to understanding the viscoplastic model results lies in the rate-dependence of the flexural wavelength of a viscous plate, and the strain rate dependence of the force increase associated with footwall and hanging wall bending. This behaviour produces unrealistic deformation patterns that can hinder the geological relevance of long-term rifting models that assume a viscoplastic rheology.

Description: Effective uptake of fermentable substrates is a fundamentally important aspect of any fermentation process. The solventogenic bacterium Clostridium beijerinckii is noted for its ability to ferment a wide range of carbohydrates, yet few of its sugar transport systems have been characterized. In common with other anaerobes, C. beijerinckii shows a marked dependence on the PEP-dependent phosphotransferase system (PTS) for sugar accumulation. In this study, the gene cbe0751 encoding the sugar-specific domains of a phosphotransferase belonging to the glucose family was cloned into an Escherichia coli strain lacking the ability to take up and phosphorylate glucose. Transformants gained ability to ferment glucose, and also mannose, and further analysis of a selected transformant demonstrated that it could take up and phosphorylate glucose, confirming that cbe0751 encodes a glucose PTS which also recognizes mannose as a substrate. RT-PCR analysis showed that cbe0751 was expressed in cultures grown on both substrates, but also to varying extents during growth on some other carbon sources. Although analogue inhibition studies suggested that Cbe0751 is not the only glucose PTS in C. beijerinckii , this system should nevertheless be regarded as a potential target for metabolic engineering to generate a strain showing improved sugar fermentation properties.

Description: Heavy metals are dense chemicals with dual biological role as micronutrients and intoxicants. A few hypersaline environmental systems are naturally enriched with heavy metals, while most metal-contaminated sites are a consequence of human activities. Numerous halotolerant and moderately halophilic Bacteria possess metal tolerance, whereas a few archaeal counterparts share similar features. The main mechanisms underlying heavy metal resistance in halophilic Bacteria and Archaea include extracellular metal sequestration by biopolymers, metal efflux mediated by specific transporters and enzymatic detoxification. Biotransformation of metals by halophiles has implications both for trace metal turnover in natural saline ecosystems and for development of novel bioremediation strategies.

Description: We screened for a gene that inhibits streptomycin production in Streptomyces griseus when it is introduced on a high-copy-number plasmid pIJ702, and obtained a plasmid pKM545. The introduction of pKM545 abolished streptomycin production on all media tested including YMP-sugar and Nutrient broth. S1 protection analysis demonstrated that the introduction of this plasmid downregulated the transcriptional activity of the promoter preceding strR , the pathway-specific transcriptional regulator for streptomycin biosynthesis. The 2.8-kb Bam HI fragment cloned onto pKM545 contained two coding sequences SGR_5442 and 5443. These coding sequences and the two downstream ones (SGR_5444 and 5445) constituted a possible operon structure designated to be rspABCD (regulation of streptomycin production). RspB and RspC exhibited a marked similarity with an ATP-binding domain and a membrane-associating domain of an ABC-2 type transporter, respectively, suggesting that the Rsp proteins comprise a membrane exporter. The gene cluster consisting of the rsp operon and the upstream divergent small coding sequence (SGR_5441) was widely distributed to Streptomyces genome. An rspB mutant of S. griseus produced 3-fold streptomycin of the parental strain in YMP liquid medium. The evidence implies that the Rsp translocator is involved in the export of a substance that specifies the expression level of streptomycin biosynthesis genes in S. griseus .

Description: We previously demonstrated that in Streptomyces coelicolor two-component system AfsQ1/Q2 activates the production of the yellow-colored coelimycin P2 (also named as yCPK) on glutamate-supplemented minimal medium, and the response regulator AfsQ1 could specifically bind to the intergenic region between two structural genes, cpkA and cpkD . Here, a more in-depth investigation was performed to elucidate the mechanism underlying the role of AfsQ1/Q2 in regulating coelimycin P2 biosynthesis. Deletion of afsQ1/Q2 resulted in markedly decreased expression of the whole coelimycin P2 biosynthetic gene cluster. Electrophoretic mobility shift assays revealed that AfsQ1 bound only to the target site identified previously, but not to any other promoters in the gene cluster. Mutations of AfsQ1-binding motif only resulted in drastically reduced transcription of the cpkA/B/C operon (encoding three type I polyketide synthases) and intriguingly, led to enhanced expression of some coelimcyin P2 genes, particularly accA1 and scF . These results suggested the direct role of AfsQ1/Q2 in regulating coelimycin production, which is directly mediated by the structural genes, but not the cluster-situated regulatory genes, and also implied that other unknown mechanisms may be involved in AfsQ1/Q2-mediated regulation of coelimycin P2 biosynthesis.

Description: Sulfate-reducing bacteria (SRBs) gain their energy by coupling the oxidation of organic substrate to the reduction of sulfate to sulfide. Several SRBs are able to use alternative terminal electron acceptors to sulfate such as nitrate. Nitrate-reducing SRBs have been isolated from a diverse range of environments. In order to be able to understand the significance of nitrate reduction in SRBs, we need to examine the ecology and physiology of the nitrate-reducing SRB isolates.

Description: The ADP-ribosylating enzymes are encoded in many pathogenic bacteria in order to affect essential functions of the host. In this study, we show that Neisseria gonorrhoeae possess a locus that corresponds to the ADP-ribosyltransferase NarE, a previously characterized enzyme in N. meningitidis . The 291 bp coding sequence of gonococcal narE shares 100% identity with part of the coding sequence of the meningococcal narE gene due to a frameshift previously described, thus leading to a 49-amino-acid deletion at the N-terminus of gonococcal NarE protein. However, we found a promoter region and a GTG start codon, which allowed expression of the protein as demonstrated by RT-PCR and western blot analyses. Using a gonococcal NarE–6xHis fusion protein, we demonstrated that the gonococcal enzyme underwent auto-ADP-ribosylation but to a lower extent than meningococcal NarE. We also observed that gonoccocal NarE exhibited ADP-ribosyltransferase activity using agmatine and cell-free host proteins as ADP-ribose acceptors, but its activity was inhibited by human β-defensins. Taken together, our results showed that NarE of Neisseria gonorrhoeae is a functional enzyme that possesses key features of bacterial ADP-ribosylating enzymes.

Description: The effect of fructose 1,6-bisphosphate (Fru 1,6-P 2 ) on the regulatory enzymes of pentose phosphate pathway of Escherichia coli was examined. Fru 1,6-P 2 inhibited E. coli transaldolase (EC 2.2.1.2) competitively against fructose 6-phosphate and uncompetitively against erythrose 4-phosphate, whereas Fru 1,6-P 2 did not affect glucose 6-phosphate dehydrogenase (EC 1.1.1.49) and 6-phosphogluconate dehydrogenase (EC 1.1.1.44). Kinetic results can be explained by assuming that transaldolase has two kinds of binding sites for Fru 1,6-P 2 : a competitive binding site for fructose 6-phosphate and a second binding site on the enzyme-erythrose 4-phosphate complex. Fru 1,6-P 2 increased resulting from the stimulation of glycolysis, can inhibit transaldolase and further participates in the elevation of the concentration of ribose 5-phosphate that can be preferentially utilized for anabolic reaction in exponential phase of E. coli .

Description: In the context of the 2014 realization of the International Terrestrial Reference Frame, the International DORIS (Doppler Orbitography Radiopositioning Integrated by Satellite) Service (IDS) has delivered to the IERS a set of 1140 weekly SINEX files including station coordinates and Earth orientation parameters, covering the time period from 1993.0 to 2015.0. From this set of weekly SINEX files, the IDS combination centre estimated a cumulative DORIS position and velocity solution to obtain mean horizontal and vertical motion of 160 stations at 71 DORIS sites. The main objective of this study is to validate the velocities of the DORIS sites by comparison with external models or time-series. Horizontal velocities are compared with two recent global plate models (GEODVEL 2010 and NNR-MORVEL56). Prior to the comparisons, DORIS horizontal velocities were corrected for Global Isostatic Adjustment from the ICE-6G (VM5a) model. For more than half of the sites, the DORIS horizontal velocities differ from the global plate models by less than 2–3 mm yr –1 . For five of the sites (Arequipa, Dionysos/Gavdos, Manila and Santiago) with horizontal velocity differences with respect to these models larger than 10 mm yr –1 , comparisons with GNSS estimates show the veracity of the DORIS motions. Vertical motions from the DORIS cumulative solution are compared with the vertical velocities derived from the latest GPS cumulative solution over the time span 1995.0–2014.0 from the University of La Rochelle solution at 31 co-located DORIS-GPS sites. These two sets of vertical velocities show a correlation coefficient of 0.83. Vertical differences are larger than 2 mm yr –1 at 23 percent of the sites. At Thule, the disagreement is explained by fine-tuned DORIS discontinuities in line with the mass variations of outlet glaciers. Furthermore, the time evolution of the vertical time-series from the DORIS station in Thule show similar trends to the GRACE equivalent water height.

Description: We have used ensemble averages of satellite-derived free-air gravity anomaly data, together with inverse modelling techniques, to determine the effective elastic thickness, T e , of circum-Pacific subducting oceanic lithosphere and its relationship to plate age. Synthetic modelling tests show that T e can be recovered best using gravity anomaly, rather than bathymetry, data and profiles that are at least 750 km long. Inverse modelling based on a uniform T e elastic plate suggests that T e increases with age of the subducting oceanic lithosphere and is given approximately by the depth to the 390 ± 10 °C oceanic isotherm based on a cooling plate model. Misfits between the observed and calculated gravity anomalies are significantly improved if a mechanically weak zone is included between the trench axis and the outer rise. This weak zone is coincident with observations of bend-faulting and seismicity. Inverse modelling shows that T e landward of the outer rise is generally 40–65 per cent less than the T e seaward of the outer rise. Both landward and seaward T e increases with age of the lithosphere and are given by the depth to the 342–349 °C and 671–714 °C oceanic isotherm, respectively. A dependence of T e on age is consistent with models for the cooling of oceanic lithosphere as it moves away from a mid-ocean ridge and the temperature-dependent ductile creep of oceanic lithospheric minerals such as olivine. By comparing the observed T e to the predicted T e based on laboratory-derived yield strength envelopes and an assumption of elastic-perfectly plastic deformation, we have attempted to constrain the rheology of oceanic lithosphere. Regardless of the assumed friction coefficient, the dry-olivine low-temperature plasticity flow laws of Goetze, Evans & Goetze, Raterron et al . and Mei et al . all provide quite a good fit to the observed T e at circum-Pacific subduction zones. This result contrasts with the Hawaiian Islands, where these flow laws are generally too strong to fit the observations. The discrepancy in rheology within Pacific plate may be caused by differences in the timescale of loading and therefore the amount of viscoelastic stress relaxation that has occurred. Other possibilities include thermal rejuvenation and magma-assisted flexure at the Hawaiian Islands.

Description: Large gaps and inconsistencies remain in published estimates of Nubia–Somalia plate motion based on reconstructions of seafloor spreading data around Africa. Herein, we use newly available reconstructions of the Southwest Indian Ridge at ~1-Myr intervals since 20 Ma to estimate Nubia–Somalia plate motion farther back in time than previously achieved and with an unprecedented degree of temporal resolution. At the northern end of the East African rift, our new estimates of Nubia–Somalia motion for six times from 0.78 Ma to 5.2 Ma differ by only 2 per cent from the rift-normal component of motion that is extrapolated from a recently estimated GPS angular velocity. The rate of rift-normal extension thus appears to have remained steady since at least 5.2 Ma. Our new rotations indicate that the two plates have moved relative to each other since at least 16 Ma and possibly longer. Motion has either been steady since at least 16 Ma or accelerated modestly between 6 and 5.2 Ma. Our Nubia–Somalia rotations predict 42.5 ± 3.8 km of rift-normal extension since 10.6 Ma across the well-studied, northern segment of the Main Ethiopian Rift, consistent with 40–50 km estimates for extension since 10.6 Myr based on seismological surveys of this narrow part of the plate boundary. Nubia–Somalia rotations are also derived by combining newly estimated Somalia–Arabia rotations that reconstruct the post-20-Ma opening of the Gulf of Aden with Nubia–Arabia rotations estimated via a probabilistic analysis of plausible opening scenarios for the Red Sea. These rotations predict Nubia–Somalia motion since 5.2 Myr that is consistent with that determined from Southwest Indian Ridge data and also predict 40 ± 3 km of rift-normal extension since 10.6 Ma across the Main Ethiopian Rift, consistent with our 42.5 ± 3.8 km Southwest Indian Ridge estimate. Our new rotations exclude at high confidence level previous estimates of 12 ± 13 and 123 ± 14 km for rift-normal extensions across the Main Ethiopian Rift since 10.6 Ma based on reconstructions of Chron 5n.2 along the Southwest Indian Ridge. Sparse coverage of magnetic reversals older than 16 Ma along the western third of the Southwest Indian Ridge precludes reliable determinations of Nubia–Somalia plate motion before 16 Ma, leaving unanswered the key question of when the motion between the two plates began.

Description: The TAMMAR segment of the Mid-Atlantic Ridge forms a classic propagating system centred about two degrees south of the Kane Fracture Zone. The segment is propagating to the south at a rate of 14 mm yr –1 , 15 per cent faster than the half-spreading rate. Here, we use seismic refraction data across the propagating rift, sheared zone and failed rift to investigate the crustal structure of the system. Inversion of the seismic data agrees remarkably well with crustal thicknesses determined from gravity modelling. We show that the crust is thickened beneath the highly magmatic propagating rift, reaching a maximum thickness of almost 8 km along the seismic line and an inferred (from gravity) thickness of about 9 km at its centre. In contrast, the crust in the sheared zone is mostly 4.5–6.5 km thick, averaging over 1 km thinner than normal oceanic crust, and reaching a minimum thickness of only 3.5 km in its NW corner. Along the seismic line, it reaches a minimum thickness of under 5 km. The PmP reflection beneath the sheared zone and failed rift is very weak or absent, suggesting serpentinisation beneath the Moho, and thus effective transport of water through the sheared zone crust. We ascribe this increased porosity in the sheared zone to extensive fracturing and faulting during deformation. We show that a bookshelf-faulting kinematic model predicts significantly more crustal thinning than is observed, suggesting that an additional mechanism of deformation is required. We therefore propose that deformation is partitioned between bookshelf faulting and simple shear, with no more than 60 per cent taken up by bookshelf faulting.

Description: The Hdr (heterodisulfide reductase)-like enzyme is predicted, from gene transcript profiling experiments previously published, to be essential in oxidative sulfur metabolism in a number of bacteria and archaea. Nevertheless, no biochemical and physicochemical data are available so far about this enzyme. Genes coding for it were identified in Aquifex aeolicus , a Gram-negative, hyperthermophilic, chemolithoautotrophic and microaerophilic bacterium that uses inorganic sulfur compounds as electron donor to grow. We provide biochemical evidence that this Hdr-like enzyme is present in this sulfur-oxidizing prokaryote (cultivated with thiosulfate or elemental sulfur). We demonstrate, by immunolocalization and cell fractionation, that Hdr-like enzyme is associated, presumably monotopically, with the membrane fraction. We show by co-immunoprecipitation assay or partial purification, that the Hdr proteins form a stable complex composed of at least five subunits, HdrA, HdrB1, HdrB2, HdrC1 and HdrC2, present in two forms of high molecular mass on native gel (~240 and 450 kDa). These studies allow us to propose a revised model for dissimilatory sulfur oxidation pathways in A. aeolicus , with Hdr predicted to generate sulfite.

Description: Geobacter sulfurreducens is an anaerobic soil bacterium that is involved in biogeochemical cycles of elements such as Fe and Mn. Although significant progress has been made in the understanding of the electron transfer processes in G. sulfurreducens , little is known about the regulatory mechanisms involved in their control. To expand the study of gene regulation in G. sulfurreducens , we carried out a genome-wide identification of transcription start sites (TSS) by 5'RACE and by deep RNA sequencing of primary mRNAs in two growth conditions. TSSs were identified along G. sulfurreducens genome and over 50% of them were located in the upstream region of the associated gene, and in some cases we detected genes with more than one TSS. Our global mapping of TSSs contributes with valuable information, which is needed for the study of transcript structure and transcription regulation signals and can ultimately contribute to the understanding of transcription initiation phenomena in G. sulfurreducens .

Description: DnpA, a putative de- N -acetylase of the PIG-L superfamily, is required for antibiotic tolerance in Pseudomonas aeruginosa . Exactly how dnpA (gene locus PA5002) directs the formation of antibiotic-tolerant persister cells is currently unknown. Previous research provided evidence for a role in surface-associated process(es), possibly in lipopolysaccharide biosynthesis. In silico sequence analysis of DnpA predicts a single transmembrane domain and N in /C out orientation of DnpA. In contrast, we here show that DnpA is an integral inner membrane protein containing two transmembrane domains, with the major C-terminal part located at the cytoplasmic face. Correct insertion into the inner membrane is necessary for DnpA to promote fluoroquinolone tolerance. The membrane localization of DnpA further supports its role in cell envelope-associated process(es). In addition to shedding light on the biological role of DnpA, this study highlights the risks of overreliance on the predictive value of bioinformatics tools and the importance of rigorous experimental validation of in silico predictions.

Description: Cell wall impermeability and active efflux of drugs are among the primary reasons for drug resistance in mycobacteria. Efflux pumps are tripartite membrane localized transport proteins that expel drug molecules outside the cells. Several of such efflux pumps are annotated in mycobacteria, but few have been characterized, like MSMEG_2991, a putative efflux pump permease of Mycobacterium smegmatis . To substantiate this, we overexpressed MSMEG_2991 protein in Escherichia coli 2443. Expression of MSMEG_2991 elevated the resistance towards structurally unrelated groups of antibiotics. An active antibiotic efflux pump nature of MSMEG_2991 was revealed by assessing the acquisition of ciprofloxacin in the absence and presence of the efflux pump inhibitor, carbonyl cyanide m-chlorophenyl hydrazone, indicating the involvement of proton-motive force (pmf) during the efflux activity. MSMEG_2991 expression elevated biofilm formation in E. coli by 4-fold, keeping parity to some of the earlier reported efflux pumps. In silico analysis suggested the presence of 12 transmembrane helices in MSMEG_2991 resembling EmrD efflux pump of E. coli . Based on in vivo and in silico analyses, MSMEG_2991 may be designated as a pmf-mediated multidrug efflux pump protein that expels diverse groups of antibiotics and might as well be involved in the biofilm enhancement.

Description: Bacterial small RNAs (sRNAs) play essential roles in the post-transcriptional control of gene expression. To improve their detection by conventional microarrays, we designed a custom microarray containing a group of probes targeting known and some putative Escherichia coli sRNAs. To assess its potential in detection of sRNAs, RNA profiling experiments were performed with total RNA extracted from E. coli MG1655 cells exponentially grown in rich (Luria–Bertani) and minimal (M9/glucose) media. We found that many sRNAs could yield reasonably strong and statistically significant signals corresponding to nearly all sRNAs annotated in the EcoCyc database. Besides differential expression of two sRNAs (GcvB and RydB), expression of other sRNAs was less affected by the composition of the growth media. Other examples of the differentially expressed sRNAs were revealed by comparing gene expression of the wild-type strain and its isogenic mutant lacking functional poly(A) polymerase I ( pcnB ). Further, northern blot analysis was employed to validate these data and to assess the existence of new putative sRNAs. Our results suggest that the use of custom microarrays with improved capacities for detection of sRNAs can offer an attractive opportunity for efficient gene expression profiling of sRNAs and their target mRNAs at the whole transcriptome level.

Description: Escherichia coli DedA/Tvp38 family proteins YghB and YqjA are putative membrane transporters with 62% amino acid identity and overlapping functions. An E. coli strain (BC202) with nonpolar yghB and yqjA mutations displays cell-division defects and temperature sensitivity and is sensitive to antibiotics and alkaline pH. In this study, we performed site-directed mutagenesis on conserved, charged amino acids of YqjA and YghB. We discovered two conserved predicted membrane-embedded arginines (R130 and R136) that are critical for function in both proteins as defined by their ability to complement BC202 phenotypes, when expressed from a plasmid. Lysine can substitute for arginine at position R130 indicating a charge dependence at this position, but could not substitute at R136. In light of the established role that arginine plays in the translocation mechanism of numerous membrane transporters, we hypothesize that these amino acids play a role in the transport mechanism of these DedA/Tvp38 family proteins.

Description: Candida albicans is a dimorphic yeast commonly found on human mucosal membranes that switches from yeast to hyphal morphology in response to environmental factors. The change to hyphal growth requires histone H3 modifications by the yeast-specific histone acetyltransferase Rtt109. In addition to its role in morphogenesis, Rtt109-dependent acetylation of histone H3 lysine residues 9 and 56 has regulatory functions during DNA replication and repair. Boric acid (BA) is a broad-spectrum agent that specifically inhibits C. albicans hyphal growth, locking the fungus in its harmless commensal yeast state. The present study characterizes the effect of BA on C. albicans histone acetylation in respect to specificity, time-course and significance. We demonstrate that sublethal concentrations of BA reduce H3K9/H3K56 acetylation, both on a basal level and in response to genotoxic stress. Acetylation at other selected histone sites were not affected by BA. qRT-PCR expression analysis of the DNA repair gene Rad51 indicated no elevated level of genotoxic stress during BA exposure. A forward-mutation analysis demonstrated the BA does not increase spontaneous or induced mutations . The findings suggest that DNA repair remains effective even when histone H3 acetylation decreases and dispels the notion that BA treatment impairs genome integrity in yeast.

Description: The impact of gene mutations on the growth of the cells can be studied using pure cultures. However, the importance of certain proteins and pathways can be also examined via co-culturing wild type and its mutant derivative. Here, the relative fitness of a mutant strain that lacks the global nitrogen regulator, CodY, was examined in Bacillus cereus , a food poisoning Gram-positive bacterium. Fitness measurements revealed that the codY strain was outcompeted when cocultured with the wild-type ATCC 14579 under various rich laboratory medium, and also when inoculated in certain beverages. In nutrient-poor minimal medium, the codY mutant had comparable fitness to the wild-type strain. Interestingly, the relative fitness of the codY strain was antagonistic when it was cultivated in apple or orange juices due to unknown properties of these beverages, highlighting the importance of chemical composition of the test medium during the bacterial fitness measurements.

Description: Out of 8000 candidates from a genetic screening for Pseudomonas putida KT2442 mutants showing defects in biofilm formation, 40 independent mutants with diminished levels of biofilm were analyzed. Most of these mutants carried insertions in genes of the lap cluster, whose products are responsible for synthesis, export and degradation of the adhesin LapA. All mutants in this class were strongly defective in biofilm formation. Mutants in the flagellar regulatory genes fleQ and flhF showed similar defects to that of the lap mutants . On the contrary, transposon insertions in the flagellar structural genes fliP and flgG , that also impair flagellar motility, had a modest defect in biofilm formation. A mutation in gacS , encoding the sensor element of the GacS/GacA two-component system, also had a moderate effect on biofilm formation. Additional insertions targeted genes involved in cell envelope function: PP3222, encoding the permease element of an ABC-type transporter and tolB , encoding the periplasmic component of the Tol-OprL system required for outer membrane stability. Our results underscore the central role of LapA, suggest cross-regulation between motility and adhesion functions and provide insights on the role of cell envelope trafficking and maintenance for biofilm development in P. putida .

Description: Horizontal gene transfer plays a significant role in bacterial evolution and has major clinical importance. Thus, it is vital to understand the mechanisms and kinetics of genetic transformations. Natural transformation is the driving mechanism for horizontal gene transfer in diverse genera of bacteria. Our study introduces a simple and rapid method for the investigation of natural transformation. This highly sensitive system allows the detection of a transformation event directly from a bacterial population without any separation step or selection of cells. The system is based on the bacterial luciferase operon from Photorhabdus luminescens . The studied molecular tools consist of the functional modules luxCDE and luxAB , which involve a replicative plasmid and an integrative gene cassette. A well-established host for bacterial genetic investigations, Acinetobacter baylyi ADP1, is used as the model bacterium. We show that natural transformation followed by homologous recombination or plasmid recircularization can be readily detected in both actively growing and static biofilm-like cultures, including very rare transformation events. The system allows the detection of natural transformation within 1 h of introducing sample DNA into the culture. The introduced method provides a convenient means to study the kinetics of natural transformation under variable conditions and perturbations.

Description: The legume–rhizobial interaction results in the formation of symbiotic nodules in which rhizobia fix nitrogen. During the process of symbiosis, reactive oxygen species (ROS) are generated. Thus, the response of rhizobia to ROS is important for successful nodulation and nitrogen fixation. In this study, we investigated how Azorhizobium caulinodans , a rhizobium that forms both root and stem nodules on its host plant, regulates ROS resistance. We found that in-frame deletions of a gene encoding the putative catalase-peroxidase katG or a gene encoding a LysR-family regulatory protein, oxyR , exhibited increased sensitivity to H 2 O 2 . We then showed that OxyR positively regulated katG expression in an H 2 O 2 -independent fashion. Furthermore, we found that deletion of katG or oxyR led to significant reduction in the number of stem nodules and decrease of nitrogen fixation capacities in symbiosis. Our results revealed that KatG and OxyR are not only critical for antioxidant defense in vitro , but also important for nodule formation and nitrogen fixation during interaction with plant hosts.

Description: A Neogene magmatic reactivation of the Manihiki Plateau, a large igneous province (LIP) in the central Pacific, is studied using seismic reflection data. Igneous diapirs have been identified exclusively within a narrow WNW–ESE striking corridor in the southern High Plateau (HP), which is parallel to the Neogene Pacific Plate motion and overlaps with an extrapolation of the Society Islands Hotspot (SIH) path. The igneous diapirs are characterized by a narrow width (〉5 km), penetration of the Neogene sediments, and they become progressively younger towards the East (23–10 Ma). The magmatic source appears to be of small lateral extent, which leads to the conclusion that the diapirs represent Neogene hotspot volcanism within a LIP, and thus may be an older, previously unknown extension of the SIH track (〉4.5 Ma). Comparing hotspot volcanism within oceanic and continental lithosphere, we further conclude that hotspot volcanism within LIP crust has similarities to tectonically faulted continental crust.

Description: Peptidoglycan (also known as murein) is an important envelope component of bacteria, and its turnover usually takes place at considerable levels during normal growth. Amino sugars and murein tripeptide resulting from murein degradation are used for resynthesis of peptidoglycan or as self-generated nutrients or energy sources for cell growth. PgrR (regulator of peptide glycan recycling; formerly YcjZ) was recently identified as a repressor of several genes participating in uptake and degradation of murein tripeptide. In this study, we identified the ycjG gene involved in murein tripeptide degradation as a new direct target of PgrR. The expression of PgrR-regulated genes including ycjY , mppA , mpaA and ycjG was repressed in the presence of a good nitrogen source, but their expression increased under poor nitrogen conditions. Under nitrogen starvation, the pgrR mutant cells exhibited faster growth than wild-type cells, implying that derepression of genes under the control of PgrR may help cells overcome nitrogen limitation. Therefore, these results suggest that nitrogen starvation induces derepression of PgrR-controlled genes involved in uptake and degradation of murein tripeptide, and this may stimulate the utilization of murein tripeptide as a nitrogen source.

Description: Many mechanisms control the state of stress within Earth plates. First-order well-known mechanisms include stresses induced by lateral variations of lithospheric density structure, sublithospheric tractions, ridge push and subduction pull. In this study, we attempt to quantify the influence of these mechanisms to understand the origin of stresses in the lithosphere, choosing the African plate (TAP) as an example. A finite-element based suite, Proshell, was developed to combine several data sets, to estimate the gravitational potential energy (GPE) of the lithosphere and to calculate stresses acting on the real (non-planar) geometry of TAP. We introduce several quantitative parameters to measure the degree of fit between the model and observations. Our modelling strategy involves nine series of numerical experiments. We start with the simplest possible model and then, step by step, build it up to be a more physically realistic model, all the while discussing the influence of each additional component. The starting (oversimplified) model series (1) is based on the CRUST2 data set for the crust and a half-space-cooling approximation of the lithospheric mantle. We then describe models (series 2–5) that account for lithospheric mantle density heterogeneities to build a more reliable GPE model. The consecutive series involve basal traction from the convective mantle (series A, C) and the rheological heterogeneity of the TAP via variations in its effective elastic thickness (series B, C). The model quality reflects the increase in complexity between series with an improving match to observed stress regimes and directions. The most complex model (series D) also accounts for the bending stresses in the elastic lithosphere and achieves a remarkably good fit to observations. All of our experiments were based on the iteration of controlling parameters in order to achieve the best fit between modelled and observed stresses, always considering physically feasible values. This gives us confidence that our methodology appropriately models the stress pattern of Africa, and that it may be further applied to other plates on Earth. Our modelling approach allows us to quantify several important features controlling the lithospheric stress pattern. Even though the initial oversimplified model does not fit the observations satisfactorily, it shows how ridge push may create significant compressive stresses in the lithosphere. More complex models show the importance of the density structure of the lithosphere, specifically in the subcrustal lithosphere. The stress regime within the TAP mainly results from a global balance of masses and mass moments between continental and oceanic parts of the plate. The orientation of stresses, in turn, is influenced more by local features expressed by topographic and crustal density variations, whereas existent subcrustal density variations appear to be smoothed by the crust above. The models show that accounting separately for either basal tractions or rheological heterogeneities brings moderate improvement, but the combination of these two mechanisms results in a substantially better match between model and observations. The bending stresses caused by isostatical re-adjustment improve the model match, but they have to be analysed with caution because of their depth-dependent nature.

Description: The Red Sea and Gulf of Aden represent two young basins that formed between Africa and Arabia since the early Oligocene, floored by oceanic crust or by transitional and thinned continental crust. While in the easternmost Gulf of Aden, the rift–drift transition can be dated chron C6 (~20.1 Ma), here we show that in the Red Sea the first pulse of seafloor spreading occurred during chron C3n.2n (~4.6 Ma) around ~17.1°N (present-day coordinates) and propagated southwards from this location, separating the Danakil microplate from Arabia. It is also shown that seafloor spreading between Arabia and Nubia started later, around chron 2A (~2.58 Ma), and propagated northwards. At present, there is no magnetic evidence for the existence of a linear spreading centre in the northern Red Sea at latitudes higher than ~24°N and in the southern Red Sea below ~14.8°N. The present-day plate kinematics of this region can be described with high accuracy by a network of five interacting plates (Nubia, Arabia, Somalia, Sinai and Danakil) and six triple junctions. For times older than anomaly 2A (~2.58 Ma) and up to anomaly 3, the absence of marine magnetic anomalies between Arabia and Nubia prevents a rigorous kinematic description of the five-plates system. However, there is strong evidence that the unique changes in plate motions during the last 5 Myr were a dramatic slowdown at chron C2 (~1.77 Ma) in the spreading or extension rates along the ridge and rift axes, thereby a good representation of the real plate motions can be obtained anyway by backward extension of the oldest Arabia–Nubia and Arabia–Danakil stage rotations determined on the basis of marine magnetic anomalies, respectively, C2–C2A and C2A–C3. The proposed kinematic reconstructions are accompanied by a geodynamic explanation for the genesis of large continent–continent fracture zones at the rift–drift transition and by an analysis of the strain associated with plate motions in Afar, northeastern Egypt and Sinai.

Description: On 2016 February 6 the South Taiwan earthquake ( M w 6.4) occurred in the Meinong District of Kaohsiung, southern Taiwan, at a depth of 17 km. It caused 117 fatalities and widespread damage to infrastructures, especially in the Tainan city. To clarify the generating mechanism of this damaging earthquake, we determined high-resolution 3-D images of P - and S -wave velocity ( V P , V S ) and Poisson's ratio ( ) in the epicentral area. We used 105 712 P - and 61 250 S -wave arrival times of 8279 local earthquakes (1.5 ≤ M ≤ 6.4) recorded at 41 seismic stations in South Taiwan during 2000–2011. In the upper crust (depth ≤ 10 km), the most remarkable feature is low- V P , low- V S and high- anomalies in areas with known active faults in the southwestern and easternmost parts of Taiwan. In contrast, high- V P , high- V S and low- anomalies become dominant in the lower crust. The hypocentre of the 2016 South Taiwan earthquake is located in a boundary zone where seismic velocity and Poisson's ratio change drastically in both the horizontal and vertical directions. Furthermore, the hypocentre is underlain by a vertically elongated high- anomaly at depths of 23–40 km, which may reflect ascending fluids from the upper (or uppermost) mantle. The low- V and high- anomalies in the upper crust coincide with areas of low heat flow, negative Bouguer gravity anomaly, and low magnetotelluric resistivity, which may reflect crustal fluids contained in the young fold-and-thrust belt. These results suggest that the 2016 South Taiwan earthquake was triggered by ascending fluids from dehydration of the subducting Eurasian slab, invading into active faults with a high loading rate.

Description: Crystallization experiments in the dendritic regime have been carried out in hypergravity conditions (from 1 to 1300 g) from an ammonium chloride solution (NH 4 Cl and H 2 O). A commercial centrifuge was equipped with a slip ring so that electric power (needed for a Peltier device and a heating element), temperature and ultrasonic signals could be transmitted between the experimental setup and the laboratory. Ultrasound measurements (2–6 MHz) were used to detect the position of the front of the mushy zone and to determine attenuation in the mush. Temperature measurements were used to control a Peltier element extracting heat from the bottom of the setup and to monitor the evolution of crystallization in the mush and in the liquid. A significant increase of solid fraction and attenuation in the mush is observed as gravity is increased. Kinetic undercooling is significant in our experiments and has been included in a macroscopic mush model. The other ingredients of the model are conservation of energy and chemical species, along with heat/species transfer between the mush and the liquid phase: boundary-layer exchanges at the top of the mush and bulk convection within the mush (formation of chimneys). The outputs of the model compare well with our experiments. We have then run the model in a range of parameters suitable for the Earth's inner core. This has shown the role of bulk mush convection for the inner core and the reason why a solid fraction very close to unity should be expected. We have also run melting experiments: after crystallization of a mush, the liquid has been heated from above until the mush started to melt, while the bottom cold temperature was maintained. These melting experiments were motivated by the possible local melting at the inner core boundary that has been invoked to explain the formation of the anomalously slow F-layer at the bottom of the outer core or inner core hemispherical asymmetry. Oddly, the consequences of melting are an increase in solid fraction and a decrease in attenuation. It is hence possible that surface seismic velocity and attenuation of the inner core are strongly affected by melting.

Description: The coupling that exists between surface processes and deformation within both the shallow crust and the deeper mantle-lithosphere has stimulated the development of computational geodynamic models that incorporate a free surface boundary condition. We introduce a treatment of this boundary condition that is suitable for staggered grid, finite difference schemes employing a structured Eulerian mesh. Our interface capturing treatment discretizes the free surface boundary condition via an interface that conforms with the edges of control volumes (e.g. a ‘staircase’ representation) and requires only local stencil modifications to be performed. Comparisons with analytic solutions verify that the method is first-order accurate. Additional intermodel comparisons are performed between known reference models to further validate our free surface approximation. Lastly, we demonstrate the applicability of a multigrid solver to our free surface methodology and demonstrate that the local stencil modifications do not strongly influence the convergence of the iterative solver.

Description: Ellesmere Island in Arctic Canada displays a complex geological evolution. The region was affected by two distinct orogenies, the Palaeozoic Ellesmerian orogeny (the Caledonian equivalent in Arctic Canada and Northern Greenland) and the Palaeogene Eurekan orogeny, related to the opening of Baffin Bay and the consequent convergence of the Greenland plate. The details of this complex evolution and the present-day deep structure are poorly constrained in this remote area and deep geophysical data are sparse. Receiver function analysis of seven temporary broad-band seismometers of the Ellesmere Island Lithosphere Experiment complemented by two permanent stations provides important data on the crustal velocity structure of Ellesmere Island. The crustal expression of the northernmost tectonic block of Ellesmere Island (~82°–83°N), Pearya, which was accreted during the Ellesmerian orogeny, is similar to that at the southernmost part, which is part of the Precambrian Laurentian (North America-Greenland) craton. Both segments have thick crystalline crust (~35–36 km) and comparable velocity–depth profiles. In contrast, crustal thickness in central Ellesmere Island decreases from ~24–30 km in the Eurekan fold and thrust belt (~79.7°–80.6°N) to ~16–20 km in the Hazen Stable Block (HSB; ~80.6°–81.4°N) and is covered by a thick succession of metasediments. A deep crustal root (~48 km) at ~79.6°N is interpreted as cratonic crust flexed beneath the Eurekan fold and thrust belt. The Carboniferous to Palaeogene sedimentary succession of the Sverdrup Basin is inferred to be up to 1–4 km thick, comparable to geologically-based estimates, near the western margin of the HSB.

Description: Joint analysis of the seismic velocities and geoid, gravity and gravity gradients are used to constrain the viscosity profile within the mantle as well as the lateral density variations. Recent ESA's Gravity field and steady-state Ocean Circulation Explorer measurements of the second-order derivatives of the Earth's gravity potential give new possibilities to determine these mantle properties. Using a simple mantle model and seismic tomography results, we investigate how the gravitational potential, the three components of the gravity vector and the gravity gradients can bring information on the radial viscosity profile and on the mantle mass anomalies. We start with lateral density variations in the Earth's mantle based either on slab history or deduced from seismic tomography. The main uncertainties are: for the latter case, the relationship between seismic velocity and density—the so-called density/velocity scaling factor—and for the former case, the variation with depth of the density contrast between the cold slabs and the surrounding mantle. We perform a Monte Carlo search for the viscosity and the density/velocity scaling factor profiles within the mantle, which allows to fit the observed geoid, gravity and gradients of gravity. We compute the posterior probability distribution of the unknown parameters, and find that the gravity gradients improve the estimate of the scaling factor within the upper mantle, because of their sensitivity to the masses within the upper mantle, whereas the geoid and the gravity better constrain the scaling factor in the lower mantle. In the upper mantle, it is less than 0.02 in the upper part and about 0.08–0.14 in the lower part, and it is significantly larger for depths greater than 1200 km (about 0.32–0.34). In any case, the density/velocity scaling factor between 670 and 1150 km depth is not well constrained. We show that the viscosity of the upper part of the mantle is strongly correlated with the viscosity of the lower part of the mantle and that the viscosity profile is characterized by a decrease in the lower part of the upper mantle (about 10 20 –2  x  10 20 Pa s) and by an increase (about 10 23 –2  x  10 23 Pa s) at the top of the lower mantle (between 670 and 1150 km). The viscosity of the mantle below 1150 km depth is well estimated in our Monte Carlo search and is about 10 22 –4  x  10 22 Pa s.

Description: Large-scale topography may be due to several causes, including (1) variations in crustal thickness and density structure, (2) oceanic lithosphere age differences, (3) subcrustal density variations in the continental lithosphere and (4) convective flow in the mantle beneath the lithosphere. The last contribution in particular may change with time and be responsible for continental inundations; distinguishing between these contributions is therefore important for linking Earth's history to its observed geological record. As a step towards this goal, this paper aims at such distinction for the present-day topography: the approach taken is deriving a ‘model’ topography due to contributions (3) and (4), along with a model geoid, using a geodynamic mantle flow model. Both lithosphere thickness and density anomalies beneath the lithosphere are inferred from seismic tomography. Density anomalies within the continental lithosphere are uncertain, because they are probably due to variations in composition and temperature, making a simple scaling from seismic to density anomalies inappropriate. Therefore, we test a number of different assumptions regarding these. As a reality check, model topography is compared, in terms of both correlation and amplitude ratio, to ‘residual’ topography, which follows from observed topography after subtracting contributions (1) and (2). The model geoid is compared to observations as well. Comparatively good agreement is found if there is either an excess density of 0.2 per cent in the lithosphere above 150 km depth, with anomalies below as inferred from tomography, or if the excess density is 0.4 per cent in the entire lithosphere. Further, a good fit is found for viscosity 10 20 Pa s in the asthenosphere, increasing to 10 23 Pa s in the lower mantle above D'. Results are quite dependent on which tomography models they are based on; for some recent ones, topography correlation is 0.6, many smaller scale features are matched, topography amplitude is less than 30 per cent too large, while geoid variance reduction exceeds 70 per cent—overall a considerable improvement compared to previous models. Correlation becomes less if smaller scale features (corresponding to spherical harmonic degrees 15 and higher), which are probably largely due to anomalies in the lithosphere, are also considered. Comparison of results with different viscosity structures, and a regional comparison of amplitude ratios, indicates that lateral viscosity variations can be quite strong, but only leading to moderate variations in model topography of a factor probably less than two.

Description: The central metabolite acetyl phosphate (acP) has long been proposed to influence transcription regulation by directly transferring its phosphoryl group to a number of response regulators in many bacterial species. Here, we provide in vitro evidence for this proposition and demonstrate, using an in vitro transcription system, that acP-dependent phosphorylation of aspartate 51 of CpxR induces transcription of one of its regulon members in E. coli , cpxP . We also used this in vitro transcription system to extend our previously reported in vivo data that hypothesized that acetylation of RNA polymerase (RNAP) influences acP-dependent cpxP transcription, using glutamine as a genetic mimic for acetylated arginine 291 of the carboxy-terminal domain of RNAP α subunit. The data we present here lend strong support to the hypothesis that acP has a direct effect on transcription regulation in E. coli via phosphorylation of CpxR, and that RNAP acetylation can modulate this response.

Description: Chaperones are a diverse class of ubiquitous proteins that assist other cellular proteins in folding correctly and maintaining their native structure. Many different chaperones cooperate to constitute the ‘proteostasis’ machinery in the cells. It has been proposed earlier that archaeal organisms could be ideal model systems for deciphering the basic functioning of the ‘protein folding machinery’ in higher eukaryotes. Several chaperone families have been characterized in archaea over the years but mostly one protein at a time, making it difficult to decipher the composition and mechanistics of the protein folding system as a whole. In order to deal with these lacunae, we have developed a database of all archaeal chaperone proteins, CrAgDb ( C haperone r epertoire in A rchaeal g enomes). The data have been presented in a systematic way with intuitive browse and search facilities for easy retrieval of information. Access to these curated datasets should expedite large-scale analysis of archaeal chaperone networks and significantly advance our understanding of operation and regulation of the protein folding machinery in archaea. Researchers could then translate this knowledge to comprehend the more complex protein folding pathways in eukaryotic systems. The database is freely available at http://14.139.227.92/mkumar/cragdb/ .

Description: Oblique convergence across Tibet leads to slip partitioning with the coexistence of strike-slip, normal and thrust motion on major fault systems. A key point is to understand and model how faults interact and accumulate strain at depth. Here, we extract ground deformation across the Haiyuan Fault restraining bend, at the northeastern boundary of the Tibetan plateau, from Envisat radar data spanning the 2001–2011 period. We show that the complexity of the surface displacement field can be explained by the partitioning of a uniform deep-seated convergence. Mountains and sand dunes in the study area make the radar data processing challenging and require the latest developments in processing procedures for Synthetic Aperture Radar interferometry. The processing strategy is based on a small baseline approach. Before unwrapping, we correct for atmospheric phase delays from global atmospheric models and digital elevation model errors. A series of filtering steps is applied to improve the signal-to-noise ratio across high ranges of the Tibetan plateau and the phase unwrapping capability across the fault, required for reliable estimate of fault movement. We then jointly invert our InSAR time-series together with published GPS displacements to test a proposed long-term slip-partitioning model between the Haiyuan and Gulang left-lateral Faults and the Qilian Shan thrusts. We explore the geometry of the fault system at depth and associated slip rates using a Bayesian approach and test the consistency of present-day geodetic surface displacements with a long-term tectonic model. We determine a uniform convergence rate of 10 [8.6–11.5] mm yr –1 with an N89 [81–97]°E across the whole fault system, with a variable partitioning west and east of a major extensional fault-jog (the Tianzhu pull-apart basin). Our 2-D model of two profiles perpendicular to the fault system gives a quantitative understanding of how crustal deformation is accommodated by the various branches of this thrust/strike-slip fault system and demonstrates how the geometry of the Haiyuan fault system controls the partitioning of the deep secular motion.

Description: In this paper, we study stability of the rate, state and temperature friction (RSTF) model. The Segall and Rice approach is used to model heat transfer at the sliding interface with its surroundings. The effect of pore pressure is not considered in the model to avoid the complex expression for critical stiffness. Linear stability analysis of the spring-mass sliding system is carried out with the ageing law under the quasistatic conditions in order to determine the critical stiffness above which sliding behaviour changes from unstable to stable or vice versa. Our numerical simulations establish that critical stiffness of the heated surface may increase or decrease from corresponding to the critical stiffness of the unheated surface depending on the relative values of two contradictory parameters related with velocity effect and temperature effect. Parametric studies are also carried out to understand shear velocity and temperature of the sliding surface dependence of steady friction. The RSTF model is also used to study the gravity induced failure of a creeping rock slope and the results are justified.

Description: In bacteria, cysteine can be synthesized from serine by two steps involving an L-serine O-acetyltransferase (SAT) and a cysteine synthase (CysK). While CysK is found in the publicly available annotated genome from Lactobacillus casei ATCC 334, a gene encoding SAT ( cysE ) is missing. In this study, we found that various strains of L. casei grew in a chemically defined medium containing sulfide as the sole sulfur source, indicating the presence of a serine O-acetyltransferase. The gene lying upstream of cysK is predicted to encode a homoserine trans-succinylase ( metA ). To study the function of this gene, it was cloned from L. casei FAM18110. The purified, recombinant protein did not acylate L-homoserine in vitro . Instead, it catalyzed the formation of O-acetyl serine from L-serine and acetyl-CoA. Furthermore, the plasmid expressing the L. casei gene complemented an Escherichia coli cysE mutant strain but not an E. coli metA mutant. This clearly demonstrated that the gene annotated as metA in fact encodes the SAT function and should be annotated as cysE .

Description: The Siple Coast region of Antarctica contains a number of fast-flowing ice streams, which control the dynamics and mass balance of the region. These ice streams are known to undergo stagnation and reactivation cycles, which lead to ice thickness changes that may be sufficient to excite a viscous solid Earth response (glacial isostatic adjustment; GIA). This study aims to quantify Siple Coast ice thickness changes during the last 2000 yr in order to determine the degree to which they might contribute to GIA and associated present-day bedrock uplift rates. This is important because accurate modelling of GIA is necessary to determine the rate of present-day ice-mass change from satellite gravimetry. Recently-published reconstructions of ice-stream variability were used to create a suite of kinematic models for the stagnation-related thickening of Kamb Ice Stream since ~1850 AD, and a GIA model was used to predict present-day deformation rates in response to this thickening. A number of longer-term loading scenarios, which include the stagnation and reactivation of ice streams across the Siple Coast over the past 2000 yr, were also constructed, and used to investigate the longer term GIA signal in the region. Uplift rates for each of the ice loading histories, based on a range of earth models, were compared with regional GPS-observed uplift rates and an empirical GIA estimate. We estimate Kamb Ice Stream to have thickened by 70–130 m since stagnation ~165 years ago. Modelled present-day vertical motion in response to this load increase peaks at –17 mm yr –1 (i.e. 17 mm yr –1 subsidence) for the weakest earth models tested here. Comparison of the solid Earth response to ice load changes throughout the last glacial cycle, including ice stream stagnation and reactivation across the Siple Coast during the last 2000 yr, with an empirical GIA estimate suggests that the upper mantle viscosity of the region is greater than 1 x 10 20 Pa s. When upper mantle viscosity values of 1 x 10 20  Pa s or smaller are combined with our suite of ice-load scenarios we predict uplift rates across Siple Coast that are at least 4 mm yr –1 smaller than those predicted by the empirical GIA estimate. GPS data are unable to further constrain model parameters due to the distance of the GPS sites from the study area. Our results demonstrate that Late Holocene ice load changes related to the stagnation and reactivation of ice streams on the Siple Coast may play a dominant role in defining the present-day uplift signal. However, both the detailed Earth structure and deglacial history of the region need to be better constrained in order to reduce uncertainties associated with the GIA signal of this region.

Description: The thickness and equivalent global sea level contribution of an improved model of the central and northern Laurentide Ice Sheet is constrained by 24 relative sea level histories and 18 present-day GPS-measured vertical land motion rates. The final model, termed Laur16, is derived from the ICE-5G model by holding the timing history constant and iteratively adjusting the thickness history, in four regions of northern Canada. In the final model, the last glacial maximum (LGM) thickness of the Laurentide Ice Sheet west of Hudson Bay was ~3.4–3.6 km. Conversely, east of Hudson Bay, peak ice thicknesses reached ~4 km. The ice model thicknesses inferred for these two regions represent, respectively, a ~30 per cent decrease and an average ~20–25 per cent increase to the load thickness relative to the ICE-5G reconstruction, which is generally consistent with other recent studies that have focussed on Laurentide Ice Sheet history. The final model also features peak ice thicknesses of 1.2–1.3 km in the Baffin Island region, a modest reduction relative to ICE-5G and unchanged thicknesses for a region in the central Canadian Arctic Archipelago west of Baffin Island. Vertical land motion predictions of the final model fit observed crustal uplift rates well, after an adjustment is made for the elastic crustal response to present-day ice mass changes of regional ice cover. The new Laur16 model provides more than a factor of two improvement of the fit to the RSL data ( 2 measure of misfit) and a factor of nine improvement to the fit of the GPS data (mean squared error measure of fit), compared to the ICE-5G starting model. Laur16 also fits the regional RSL data better by a factor of two and gives a slightly better fit to GPS uplift rates than the recent ICE-6G model. The volume history of the Laur16 reconstruction corresponds to an up to 8 m reduction in global sea level equivalent compared to ICE-5G at LGM.

Description: The 2011 October 23 M W 7.1 Van earthquake in eastern Turkey caused ~600 deaths and caused widespread damage and economic loss. The seismogenic rupture was restricted to 10–25 km in depth, but aseismic surface creep, coincident with outcrop fault exposures, was observed in the hours to months after the earthquake. We combine observations from radar interferometry, seismology, geomorphology and Quaternary dating to investigate the geological slip rate and seismotectonic context of the Van earthquake, and assess the implications for continuing seismic hazard in the region. Transient post-seismic slip on the upper Van fault started immediately following the earthquake, and decayed over a period of weeks; it may not fully account for our long-term surface slip-rate estimate of ≥0.5 mm yr –1 . Post-seismic slip on the Bostanici splay fault initiated several days to weeks after the main shock, and we infer that it may have followed the M W  5.9 aftershock on the 9th November. The Van earthquake shows that updip segmentation can be important in arresting seismic ruptures on dip-slip faults. Two large, shallow aftershocks show that the upper 10 km of crust can sustain significant earthquakes, and significant slip is observed to have reached the surface in the late Quaternary, so there may be a continuing seismic hazard from the upper Van fault and the associated splay. The wavelength of folding in the hanging wall of the Van fault is dominated by the structure in the upper 10 km of the crust, masking the effect of deeper seismogenic structures. Thus, models of subsurface faulting based solely on surface folding and faulting in regions of reverse faulting may underestimate the full depth extent of seismogenic structures in the region. In measuring the cumulative post-seismic offsets to anthropogenic structures, we show that Structure-from-Motion can be rapidly deployed to create snapshots of post-seismic displacement. We also demonstrate the utility of declassified Corona mission imagery (1960s–1970s) for geomorphic mapping in areas where recent urbanization has concealed the geomorphic markers.

Description: We purified a chitinase from pilei extractions of Coprinopsis cinerea fruiting bodies by ammonium sulfate precipitation and CM Sepharose cation exchange chromatography. MALDI-TOF/TOF MS analysis characterized this purified chitinase as a putative class V chitinase, ChiB1. ChiB1 hydrolyzed colloidal chitin and chitosan, whereas it did not hydrolyze chitin powder. ChiB1 cleaved only pNP-(GlcNAc) 2 , rather than pNP-GlcNAc or pNP-(Glc-NAc) 3 , to release nitrophenol. ChiB1 preferably and progressively released (GlcNAc)2 from (GlcNAc)6 and digested (GlcNAc)6 to two molecules of (GlcNAc)3 in a small proportion, but did not split (GlcNAc) 2 , so it is an exochitinase. ChiB1 has an optimum temperature range of 35°C to 40°C and an optimum pH of 5.0. ChiB1 exhibited Km and Vmax values of 2.63 mg ml –1 and 2.31 μmol min –1  mg protein –1 for colloidal chitin, respectively. The ChiB1 gene, along with another putative endochitinase (class III chitinase gene), was expressed dominantly among eight predicted chitinase genes in the genome, and its expression level increased with the maturation of fruiting bodies. ChiB1 incubation released a large amount of soluble β-glucan fractions from alkali-insoluble cell wall fractions of C. cinerea fruiting bodies, thereby it may promote the degradation of cell walls in synergy with the β-1,3-glucanases during pileus autolysis.

Description: Acidithiobacillus ferrivorans is an acidophile implicated in low-temperature biomining for the recovery of metals from sulfide minerals. Acidithiobacillus ferrivorans obtains its energy from the oxidation of inorganic sulfur compounds, and genes encoding several alternative pathways have been identified. Next-generation sequencing of At. ferrivorans RNA transcripts identified the genes coding for metabolic and electron transport proteins for energy conservation from tetrathionate as electron donor. RNA transcripts suggested that tetrathionate was hydrolyzed by the tetH1 gene product to form thiosulfate, elemental sulfur and sulfate. Despite two of the genes being truncated, RNA transcripts for the SoxXYZAB complex had higher levels than for thiosulfate quinone oxidoreductase ( doxDA genes). However, a lack of heme-binding sites in soxX suggested that DoxDA was responsible for thiosulfate metabolism. Higher RNA transcript counts also suggested that elemental sulfur was metabolized by heterodisulfide reductase ( hdr genes) rather than sulfur oxygenase reductase ( sor ). The sulfite produced as a product of heterodisulfide reductase was suggested to be oxidized by a pathway involving the sat gene product or abiotically react with elemental sulfur to form thiosulfate. Finally, several electron transport complexes were involved in energy conservation. This study has elucidated the previously unknown At. ferrivorans tetrathionate metabolic pathway that is important in biomining.

Description: TRP1 is a frequently used auxotrophic marker for genetic modifications and selections in trp – budding yeast strains, including the commonly used wild-type strain W303a. However, we found that introduction of the TRP1 gene into a trp – strain significantly affected vegetative growth at low and high temperatures. Therefore, caution should be needed when working in a trp – background strain and using the TRP1 marker to study stress response phenotypes, particularly when analyzing temperature sensitivities.

Description: The factors driving the dominance of the Calvin–Benson–Bassham cycle (CBB) or reductive citric acid cycle (rCAC) in autotrophic microorganisms in different habitats are debated. Based on costs for synthesizing a few metabolic intermediates, it has been suggested that the CBB poses a disadvantage due to higher metabolic cost. The purpose of this study was to extend this estimate of cost from metabolite synthesis to biomass synthesis. For 12 gammaproteobacteria (CBB) and five epsilonproteobacteria (rCAC), the amount of ATP to synthesize a gram of biomass from CO 2 was calculated from genome sequences via metabolic maps. The eleven central carbon metabolites needed to synthesize biomass were all less expensive to synthesize via the rCAC (66%–89% of the ATP needed to synthesize them via CBB). Differences in cell compositions did result in differing demands for metabolites among the organisms, but the differences in cost to synthesize biomass were small among organisms that used a particular pathway (e.g. rCAC), compared to the difference between pathways (rCAC versus CBB). The rCAC autotrophs averaged 0.195 moles ATP per g biomass, while their CBB counterparts averaged 0.238. This is the first in silico estimate of the relative expense of both pathways to generate biomass.

Description: Within a slab experiencing present-day lateral break-off, a particular type of earthquakes is expected to cluster at the detachment horizon tip: namely, events generated by reverse faulting, with the approximately horizontal compression involved acting along the strike of the slab. Such a cluster of moderate magnitude earthquakes (4.7 ≤  m b ≤ 5.0) was identified in this study at the 160–175 km depth range of the Vrancea seismogenic body, in the Southeast Carpathians mountains collision environment. The corresponding cluster epicentres were systematically positioned at the boundary between a region being subject ( cf. published GPS records), to present-day upward movements, and another one that underwent present-day subsidence. Such an overall setting seems to suggest that a lateral break-off is currently developing at the indicated depth within the Vrancea slab, leading to topographic uplift above the already detached slab section, and to enhanced subsidence above the section to which the gravitational slab pull was being transferred. In addition, by taking into account some systematic time correspondence which we documented between moderate magnitude events of the 160–175 km depth cluster and the subsequent strong Vrancea shocks ( M w ≥ 6.9), it appears that the latter, although occurring at much shallower depths (roughly, in the 80–140 km range), were also controlled by the break-off progress.

Description: Bajada del Diablo is located in the Northern Patagonian Massif, Chubut, Argentina. The study area includes several circular structures found in Miocene olivine basalts of the Quiñelaf Eruptive Complex and in the Late Pliocene/Early Pleistocene Pampa Sastre conglomerates. An impact origin has been proposed for these circular structures. With the aim of further investigate the proposed impact origin, topographic, gravimetric, magnetic, resistivity, palaeomagnetic and electromagnetic surveys in two circular structures (‘8’ and ‘G’) located in Pampa Sastre conglomerates and in basalts of the Quiñelaf Eruptive Complex were carried out. The new geophysical results support the hypothesis of an impact origin. However, the confirmation of such an origin through the findings of shock metamorphism evidences and/or the recovery of meteorites remains elusive.

Description: Bacteria carry a number of genes that cause cell growth arrest or cell lysis upon expression. Notably, defective prophages retain many lysis proteins. Here, we identified a novel lytic gene, ydfD , on the Qin prophage segment of the Escherichia coli genome. YdfD lyses 99.9% of cells within 2 h of its induction. The co-expression of the upstream gene, dicB , encoding a cell division inhibitor, as well as sulA , encoding another cell division inhibitor, abolished YdfD-induced cell lysis. These results imply that YdfD-induced lysis is a cell division-dependent event. We further found that by deleting the hydrophobic 22-residue N-terminal domain, the resulting 42-residue C-terminal domain was still toxic to cause cell lysis. We propose that YdfD, associated with the cytoplasmic membrane, inhibits an essential cellular process(s).

Description: Estimating in situ stress based on hydraulic fracturing data typically depends on interpretation of the breakdown, secondary breakdown (‘reopening’) and shut-in pressure. While it has been recognized that the near-wellbore stress field should be taken into account and that the compressibility of the injection system and the viscous flow of the fluid can diminish the accuracy of stress estimates, these issues have not been well quantified. A coupled numerical model that includes the compressibility of the injection system and the flow of a viscous fluid in a plane-strain hydraulic fracture extending from a wellbore, in an impermeable rock, and in the presence of a non-isotropic in situ stress field provides a basic tool for estimating the order of the error associated with hydraulic fracturing stress measurements under non-ideal conditions. The main findings of this work are model-based guidelines on the values of relevant dimensionless parameter groups to ensure sufficient accuracy of stress estimates that use idealized models. When these guidelines cannot be met under field conditions, the model can be further applied to obtain first-order corrections that account for compressibility, viscosity and near-wellbore effects.

Description: We consider the problem of seismic velocity change estimation using ambient noise recordings. Motivated by Zhan et al. , we study how the velocity change estimation is affected by seasonal fluctuations in the noise sources. More precisely, we consider a numerical model and introduce spatio-temporal seasonal fluctuations in the noise sources. We show that indeed, as pointed out by Zhan et al. , the stretching method is affected by these fluctuations and produces misleading apparent velocity variations which reduce dramatically the signal to noise ratio of the method. We also show that these apparent velocity variations can be eliminated by an adequate normalization of the cross-correlation functions. Theoretically we expect our approach to work as long as the seasonal fluctuations in the noise sources are uniform, an assumption which holds for closely located seismic stations. We illustrate with numerical simulations in homogeneous and scattering media that the proposed normalization significantly improves the accuracy of the velocity change estimation. Similar behaviour is also observed with real data recorded in the Aegean volcanic arc. We study in particular the volcano of Santorini during the seismic unrest of 2011–2012 and observe a decrease in the velocity of seismic waves which is correlated with GPS measured elevation.

Description: Many bacteria are capable of accumulating intracellular granules of polyhydroxyalkanoates (PHA). In this work, we developed confocal microscopy analysis of bacterial cells to study changes in the diameters of cells as well as PHA granules during growth and PHA accumulation in the bacterium Cupriavidus necator H16 (formerly Ralstonia eutropha ). The cell envelope was stained by DiD ® fluorescent probe and PHA granules by Nile Red. Signals from both probes were separated based on their spectral and fluorescence life-time properties. During growth and PHA accumulation, bacterial cells increased their length but the width of the cells remained constant. The volume fraction of PHA granules in cells increased during PHA accumulation, nevertheless, its value did not exceed 40 vol. % regardless of the PHA weight content. It seems that bacterial cultures lengthen the cells in order to control the PHA volume portion. However, since similar changes in cell length were also observed in a PHA non-accumulating mutant, it seems that there is no direct control mechanism, which regulates the prolongation of the cells with respect to PHA granules volume. It is more likely that PHA biosynthesis and the length of cells are influenced by the same external stimuli such as nutrient limitation.

Description: Three putative regulatory genes, namely ttmRI, ttmRII and ttmRIII , which are present in the tetramycin ( ttm ) biosynthetic gene cluster, were found in Streptomyces ahygroscopicus . Disruption of ttmRI, ttmRII or ttmRIII reduced tetramycin production, and their complementation restored production to varying degrees. Gene expression analysis of the wild-type (WT) and mutant strains through reverse transcriptase–polymerase chain reaction (RT-PCR) of the ttm gene cluster showed that the expression levels of most of the biosynthetic genes were reduced in ttmRI , ttmRII and ttmRIII . Electrophoretic mobility shift assays demonstrated that TtmRI, TtmRII and TtmRIII bound the promoters of several genes in the ttm gene cluster. This study found that these three proteins are a group of positive regulators that activate the transcription of the ttm gene cluster in S. ahygroscopicus . In addition, ttmRII had a reduced degree of grey pigmentation. Thus, TtmRII has a pleiotropic regulatory function in the tetramycin biosynthetic pathway and in development.

Description: Sulfide:quinone oxidoreductase (SQR) is the primary sulfide-oxidizing enzyme found in all three domains of life. Of the six phylogenetically distinct types of SQR, four have representatives that have been biochemically characterized. The genome of Chlorobaculum tepidum encodes three SQR homologs. One of these, encoded by CT1087, is a type VI SQR that has been previously shown to be required for growth at high sulfide concentrations and to be expressed in sulfide-dependent manner. Therefore, CT1087 was hypothesized to be a high sulfide adapted SQR. CT1087 was expressed in Escherichia coli with an N-terminal His-tag (CT1087NHis 6 ) and purified by Ni-NTA chromatography. CT1087NHis 6 was active and contained FAD as a strongly bound cofactor. The measured kinetic parameters for CT1087NHis 6 indicate a low affinity for sulfide and a high enzymatic turnover rate consistent with the hypothesis for its function inferred from genetic and expression data. These are the first kinetic data for a type VI SQR and have implications for structure-function analyses of all SQR's.

Description: Chloroflexus aggregans is an unbranched multicellular filamentous bacterium having the ability of gliding motility. The filament moves straightforward at a constant rate, ~3 μm sec –1 on solid surface and occasionally reverses the moving direction. In this study, we successfully detected movements of glass beads on the cell-surface along long axis of the filament indicating that the cell-surface movement was the direct force for gliding. Microscopic analyses found that the cell-surface movements were confined to a cell of the filament, and each cell independently moved and reversed the direction. To understand how the cellular movements determine the moving direction of the filament, we proposed a discrete-time stochastic model; sum of the directions of the cellular movements determines the moving direction of the filament only when the filament pauses, and after moving, the filament keeps the same directional movement until all the cells pause and/or move in the opposite direction. Monte Carlo simulation of this model showed that reversal frequency of longer filaments was relatively fixed to be low, but the frequency of shorter filaments varied widely. This simulation result appropriately explained the experimental observations. This study proposed the relevant mechanism adequately describing the motility of the multicellular filament in C. aggregans .

Description: Members of the Gram-positive lactic acid bacteria (LAB) are well-known for their beneficial properties as starter cultures and probiotics. Many LAB species produce ribosomally synthesized proteinaceous antibiotics (bacteriocins). Weissella confusa MBF8-1 is a strain isolated from a fermented soybean product that not only produces useful exopolysaccharides but also exhibits bacteriocin activity, which we call weissellicin MBF. Here, we show that bacteriocin production by W. confusa MBF8-1 is specified by a large plasmid, pWcMBF8-1. Plasmid pWcMBF8-1 (GenBank accession number KR350502), which was identified from the W. confusa MBF8-1 draft genome sequence, is 17 643 bp in length with a G + C content of 34.8% and contains 25 open reading frames (ORFs). Six ORFs constitute the weissellicin MBF locus, encoding three putative double-glycine-motif peptides (Bac1, Bac2, Bac3), an ABC transporter complex (BacTE) and a putative immunity protein (BacI). Two ORFs encode plasmid partitioning and mobilization proteins, suggesting that pWcMBF8-1 is transferable to other hosts. To the best of our knowledge, plasmid pWcMBF8-1 not only represents the first large Weissella plasmid to be sequenced but also the first to be associated with bacteriocin production in W. confusa .

Description: Passive seismic experiment was carried out at the SW contact of the Dinarides and Pannonian basin to determine the crustal structure and velocity discontinuities. The aim of the experiment was to define the relationship between the Adriatic microplate and the Pannonian segment as a part of the European plate. Most of the temporary seismic stations were deployed in Croatia along the Alp07 profile—a part of the active-source ALP 2002 project. About 300-km-long profile stretches from Istra peninsula to the Drava river, in a WSW–ESE direction. Teleseismic events recorded on 13 temporary seismic stations along the profile were analysed by P -receiver function method. Two types of characteristic receiver functions (RF) have been identified, belonging to Dinaridic and Pannonian crusts as defined on the Alp07 profile, while in transitional zone there are both types. Three major crustal discontinuities can be identified for the Dinaridic type: sedimentary basement, intracrustal discontinuity and Mohorovičić discontinuity, whereas the Pannonian type revealed only two discontinuities. The intracrustal discontinuity was not observed in the Pannonian type, thus pointing to a single-layered crust in the Pannonian basin. Two interpretation methods were applied: forward modelling of the receiver functions and H – stacking method, and the results were compared with the active-source seismic data at deep refraction profile Alp07. The receiver function modelling has given reliable results of the Moho depths that are in accordance with the seismic refraction results at the end of the Alp07 profile, that is in the area of Pannonian crust characterized by simple crustal structure and low seismic velocities ( Vp between 5.9 and 6.2 km s –1 ). In the Dinarides and its peripheral parts, receiver function modelling regularly gives greater Moho depths, up to +15 per cent, due to more complex crustal structure. The depths of the Moho calculated by the H – stacking method vary within wide limits (±13 km), due to band limited data of short-period stations. The results at five stations have to be rejected because of huge deviations in comparison with all previous results, while at the other seven stations the Moho depths vary within ±15 per cent around the Moho discontinuity of the Alp07 profile.

Description: Megathrust earthquakes of magnitude close to 9 are followed by large-scale (thousands of km) and long-lasting (decades), significant crustal and mantle deformation. This deformation can be observed at the surface and quantified with GPS measurements. Here we report on deformation observed during the 5 yr time span after the 2010 M w 8.8 Maule Megathrust Earthquake (2010 February 27) over the whole South American continent. With the first 2 yr of those data, we use finite element modelling (FEM) to relate this deformation to slip on the plate interface and relaxation in the mantle, using a realistic layered Earth model and Burgers rheologies. Slip alone on the interface, even up to large depths, is unable to provide a satisfactory fit simultaneously to horizontal and vertical displacements. The horizontal deformation pattern requires relaxation both in the asthenosphere and in a low-viscosity channel along the deepest part of the plate interface and no additional low-viscosity wedge is required by the data. The vertical velocity pattern (intense and quick uplift over the Cordillera) is well fitted only when the channel extends deeper than 100 km. Additionally, viscoelastic relaxation alone cannot explain the characteristics and amplitude of displacements over the first 200 km from the trench and aseismic slip on the fault plane is needed. This aseismic slip on the interface generates stresses, which induce additional relaxation in the mantle. In the final model, all three components (relaxation due to the coseismic slip, aseismic slip on the fault plane and relaxation due to aseismic slip) are taken into account. Our best-fit model uses slip at shallow depths on the subduction interface decreasing as function of time and includes (i) an asthenosphere extending down to 200 km, with a steady-state Maxwell viscosity of 4.75 x 10 18 Pa s; and (ii) a low-viscosity channel along the plate interface extending from depths of 55–135 km with viscosities below 10 18 Pa s.

Description: We provide an updated present-day stress map for the Italian territory. Following the World Stress Map (WSM) Project guidelines, we list the different stress indicators, explaining the criteria used to select data. We discuss the data, which will also be included in the 2016 release of the WSM, highlighting the areas for which we have added stress information. Our map displays the minimum horizontal stress orientations inferred from crustal stress indicators down to 40 km depth using data of A–C quality, updated for earthquakes until December 2015. We have completely reviewed all data, and the data set now contains 855 entries, in contrast to the previous 715. The number of data with A–C quality of 630 corresponds to an increase of 26 per cent relative to the previous data set. In particular, the new data set contains the results of the analysis of borehole breakouts, critically reviewed data from earthquake focal mechanisms, data concerning active faults, formal inversions of focal mechanisms of seismic sequences or of restricted areas and one stress determination from overcoring. The new data set defines the stress field in areas not well covered by the previous data: the region north to the Po Plain and the central Adriatic sea, both characterized by a thrust- and strike-faulting regime, the northern Sicilian belt with a prevailing normal-faulting regime, and the Ionian sea with a strike-slip regime.

Description: Small non-coding regulatory RNAs (sRNAs) are key players in post-transcriptional regulation of gene expression. Hundreds of sRNAs have been identified in Sinorhizobium meliloti , but their biological function remains unknown for most of them. In this study, we characterized the expression pattern of the gene encoding the 77-nt sRNA MmgR in S. meliloti strain 2011. A chromosomal transcriptional reporter fusion (P mmgR - gfp ) showed that the mmgR promoter is active along different stages of the interaction with alfalfa roots. In pure cultures, P mmgR - gfp activity paralleled the sRNA abundance indicating that mmgR expression is primarily controlled at the level of transcriptional initiation. P mmgR - gfp activity was higher during growth in rhizobial defined medium (RDM) than in TY medium. Furthermore, P mmgR - gfp was induced at 60 min after shifting growing cells from TY to RDM medium, i.e. shorter than the cell doubling time. In defined RDM medium containing NO 3 – , both P mmgR - gfp and MmgR level were repressed by the addition of tryptone or single amino acids, suggesting that mmgR expression depends on the cellular nitrogen (N) status. In silico analysis failed to detect conserved motifs upstream the promoter RNA polymerase binding site, but revealed a strongly conserved motif centered at –28 that may be linked to the observed regulatory pattern by the N source.

Description: A β-glycoside hydrolase was isolated from cell walls material in Coprinopsis cinerea elongating stipes. By analysis of SDS-PAGE, MALDI-TOF/TOF MS and substrate specificity, this enzyme was characterized as an extracellular β-glucosidase which is a trimer consisting of three homosubunits. β-Glucosidase did not degrade β-glucans with modified ends, whereas it hydrolyzed various β-glucans with free ends and related oligosaccharides with β-1,3-, β-1,4- or β-1,6-linkages. Although this β-glucosidase possesses glycosyltransferase activity on laminarioligosaccharides, it did not transfer glucose residues from laminaritriose to β-glucan in stipe cell walls to produce larger β-glucan molecules; instead, it caused a decrease in the molecular size of stipe wall β-glucan by removing glucose. Relatively, the molecular size of wall β-glucans in the elongating apical stipe was less than that found in the non-elongating basal stipes, and this β-glucosidase was more highly expressed in the elongating apical stipe than in non-elongating basal regions. Therefore, we propose that β-glucosidase functions by trimming or cutting the β-glucan side chains on the β-1,3-glucan backbone to prevent them from forming longer branches, keeping the wall plastic to promote diffuse wall growth.

Description: Measurements of the velocity field associated with plumes rising through a viscous fluid are performed using stereoscopic Particle-Image Velocimetry in the Rayleigh number range 4.4 x 10 5 –6.4 x 10 5 . The experimental model is analogous to a mantle plume rising from the core–mantle boundary to the base of the lithosphere. The behaviour of the plume is studied throughout its life cycle, which is broken up into four stages; the Formation Stage, when the plume forms; the Rising Stage, when the plume rises through the fluid; the Spreading Stage, when the plume reaches the surface and spreads; and finally the Declining Stage, when the heat source has been removed and the plume weakens. The latter three stages are examined in terms of the Finite-Time Lyapunov Exponent fields and the advection of passive tracers throughout the flow. The temperature at the heater and near the fluid surface are measured using thermocouples to infer how the presence of a mantle plume would produce excess temperature near the lithosphere throughout the various stages of its life cycle. In all experiments, a time lag is observed between the removal of the heat source and the decline in the excess temperature near the surface, which is proportional to the rise time. A simple analytical model is presented, which suggests that under mantle conditions (i.e. negligible thermal diffusion), the relationship between the time lag and the rise time is robust and independent of the Rayleigh number; however, the constant of proportionality is closer to unity in the absence of diffusion. Once the heat source is removed, the excess temperature near the surface declines exponentially at a rate that is inversely proportional to the rise time. The implications of this result are discussed in terms of the decline in volcanism in the Louisville hotspot chain over the past 20 Ma. The rise velocity of material in the plume is examined; the rise velocity is found to vary significantly with the plume height in a manner that is inconsistent with many of the common semi-analytical models of thermal plumes in the literature. It is also argued that this height dependency will cause estimates of the rise velocity based on the decay series of uranium isotopes to significantly underestimate the true value.

Description: Streptococcus mutans harbours an intracellular, human DPP IV-analogous enzyme Xaa-Pro dipeptidyl-peptidase (EC 3.4.14.11). According to previous reports, an extracellular isozyme in S. gordonii and S. suis has been associated with virulence. Speculating that even an intracellular form may aid in virulence of S. mutans , we have tried to purify, characterize and evaluate enzyme inhibition by specific inhibitors. The native enzyme was partially purified by ion-exchange and gel filtration chromatography. Owing to low yield, the enzyme was overexpressed in Lactococcus lactis and purified by affinity chromatography. The recombinant enzyme (rSm-XPDAP) had a specific activity of 1070 U mg –1 , while the V max and K m were 7 μM min –1 and 89 ± 7 μM ( n = 3), respectively. The serine protease inhibitor phenylmethylsulphonyl fluoride and a DPP IV-specific inhibitor diprotin A proved to be active against rSm-XPDAP. As a novel approach, the evaluation of the effect of anti-human DPP IV (AHD) drugs on rSm-XPDAP activity found saxagliptin to be effective to some extent ( K i = 129 ± 16 μM), which may lead to the synthesis and development of a new class of antimicrobial agents.

Description: The Makran subduction zone is one of the last convergent margins to be investigated using space-based geodesy. While there is a lack of historical and modern instrumentation in the region, a sparse sampling of continuous and campaign measurements over the past decade has allowed us to make the first estimates of convergence rates. We combine GPS measurements from 20 stations located in Iran, Pakistan and Oman along with hypocentral locations from the International Seismological Centre to create a preliminary 3-D estimate of the geometry of the megathrust, along with a preliminary fault-coupling model for the Makran subduction zone. Using a convergence rate which is strongly constrained by measurements from the incoming Arabia plate along with the backslip method of Savage, we find the Makran subduction zone appears to be locked to a depth of at least 38 km and accumulating strain.We also find evidence for a segmentation of plate coupling, with a 300 km long section of reduced plate coupling. The range of acceptable locking depths from our modelling and the 900 km along-strike length for the megathrust, makes the Makran subduction zone capable of earthquakes up to M w  = 8.8. In addition, we find evidence for slow-slip-like transient deformation events on two GPS stations. These observations are suggestive of transient deformation events observed in Cascadia, Japan and elsewhere.

Description: We present numerical models of mantle dynamics forced by plate velocities history in the last 450 Ma. The lower-mantle rheology and the thickness of a dense basal layer are systematically varied and several initial procedures are considered for each case. For some cases, the dependence on the mantle convection vigour is also examined. The resulting evolution of the CMB heat flux is analysed in terms of criteria to promote or inhibit reversals inferred from numerical dynamos. Most models present a rather dynamic lower mantle with the emergence of two thermochemical piles towards present-day. Only a small minority of models present two stationary piles over the last 450 Myr. At present-day, the composition field obtained in our models is found to correlate better with tomography than the temperature field. In addition, the temperature field immediately at the CMB (and thus the heat flux pattern) slightly differs from the average temperature field over the 100-km thick mantle layer above it. The evolution of the mean CMB heat flux or of the amplitude of heterogeneity seldom presents the expected correlation with the evolution of the palaeomagnetic reversal frequency suggesting these effects cannot explain the observations. In contrast, our analysis favours ‘inertial control’ on the geodynamo associated with polar cooling and in some cases break of Taylor columns in the outer core as sources of increased reversal frequency. Overall, the most likely candidates among our mantle dynamics models involve a viscosity increase in the mantle equal or smaller than 30: models with a discontinuous viscosity increase at the transition zone tend to agree better at present-day with observations of seismic tomography, but models with a gradual viscosity increase agree better with some of the criteria proposed to affect reversal frequency.

Description: Lateral viscosity variations (LVVs) in the mantle influence geodynamic processes and their surface expressions. With the observed long-wavelength geoid, free-air anomaly, gravity gradient in three directions and discrete, high-accuracy residual topography, we invert for depth- and temperature-dependent and tectonically regionalized mantle viscosity with a mantle flow model. The inversions suggest that long-wavelength gravitational and topographic signals are mainly controlled by the radial viscosity profile; the pre-Cambrian lithosphere viscosity is slightly (~ one order of magnitude) higher than that of oceanic and Phanerozoic lithosphere; plate margins are substantially weaker than plate interiors; and viscosity has only a weak apparent, dependence on temperature, suggesting either a balancing between factors or a smoothing of actual higher amplitude, but short wavelength, LVVs. The predicted large-scale lithospheric stress regime (compression or extension) is consistent with the world stress map (thrust or normal faulting). Both recent compiled high-accuracy residual topography and the predicted dynamic topography yield ~1 km amplitude long-wavelength dynamic topography, inconsistent with recent studies suggesting amplitudes of ~100 to ~500 m. Such studies use a constant, positive admittance (transfer function between topography and gravity), in contrast to the evidence which shows that the earth has a spatially and wavelength-dependent admittance, with large, negative admittances between ~4000 and ~10 4 km wavelengths.

Description: The hydrocarbonoclastic bacterium Acinetobacter venetianus RAG-1 has attracted substantial attention due to its powerful oil-degrading capabilities and its potential to play an important ecological role in the cleanup of alkanes. In this study, we compare the transcriptome of the strain RAG-1 grown in dodecane, the corresponding alkanol (dodecanol), and sodium acetate for the characterization of genes involved in dodecane uptake and utilization. Comparison of the transcriptional responses of RAG-1 grown on dodecane led to the identification of 1074 genes that were differentially expressed relative to sodium acetate. Of these, 622 genes were upregulated when grown in dodecane. The highly upregulated genes were involved in alkane catabolism, along with stress response. Our data suggest AlkMb to be primarily involved in dodecane oxidation. Transcriptional response of RAG-1 grown on dodecane relative to dodecanol also led to the identification of permease, outer membrane protein and thin fimbriae coding genes potentially involved in dodecane uptake. This study provides the first model for key genes involved in alkane uptake and metabolism in A. venetianus RAG-1.

Description: Dissimilatory sulfate reduction is the central microbial metabolism in global sulfur cycling. Understanding the importance of sulfate reduction to Earth's biogeochemical S cycle requires aggregating single-cell processes with geochemical signals. For sulfate reduction, these signals include the ratio of stable sulfur isotopes preserved in minerals, as well as the hydrogen isotope ratios and structures of microbial membrane lipids preserved in organic matter. In this study, we cultivated the model sulfate reducer, Desulfovibrio vulgaris DSM 644 T , to investigate how these parameters were perturbed by changes in expression of the protein DsrC. DsrC is critical to the final metabolic step in sulfate reduction to sulfide. S and H isotopic fractionation imposed by the wild type was compared to three mutants. Discrimination against 34 S in sulfate, as calculated from the residual reactant, did not discernibly differ among all strains. However, a closed-system sulfur isotope distillation model, based on accumulated sulfide, produced inconsistent results in one mutant strain IPFG09. Lipids produced by IPFG09 were also slightly enriched in 2 H. These results suggest that DsrC alone does not have a major impact on sulfate-S, though may influence sulfide-S and lipid-H isotopic compositions. While intriguing, a mechanistic explanation requires further study under continuous culture conditions.

Description: We present a new, computationally efficient numerical method to simulate global seismic wave propagation in realistic 3-D Earth models. We characterize the azimuthal dependence of 3-D wavefields in terms of Fourier series, such that the 3-D equations of motion reduce to an algebraic system of coupled 2-D meridian equations, which is then solved by a 2-D spectral element method (SEM). Computational efficiency of such a hybrid method stems from lateral smoothness of 3-D Earth models and axial singularity of seismic point sources, which jointly confine the Fourier modes of wavefields to a few lower orders. We show novel benchmarks for global wave solutions in 3-D structures between our method and an independent, fully discretized 3-D SEM with remarkable agreement. Performance comparisons are carried out on three state-of-the-art tomography models, with seismic period ranging from 34 s down to 11 s. It turns out that our method has run up to two orders of magnitude faster than the 3-D SEM, featured by a computational advantage expanding with seismic frequency.

Description: We report on a study to explore the deep electrical conductivity structure of the Dead Sea Basin (DSB) using magnetotelluric (MT) data collected along a transect across the DSB where the left lateral strike-slip Dead Sea transform (DST) fault splits into two fault strands forming one of the largest pull-apart basins of the world. A very pronounced feature of our 2-D inversion model is a deep, subvertical conductive zone beneath the DSB. The conductor extends through the entire crust and is sandwiched between highly resistive structures associated with Precambrian rocks of the basin flanks. The high electrical conductivity could be attributed to fluids released by dehydration of the uppermost mantle beneath the DSB, possibly in combination with fluids released by mid- to low-grade metamorphism in the lower crust and generation of hydrous minerals in the middle crust through retrograde metamorphism. Similar high conductivity zones associated with fluids have been reported from other large fault systems. The presence of fluids and hydrous minerals in the middle and lower crust could explain the required low friction coefficient of the DST along the eastern boundary of the DSB and the high subsidence rate of basin sediments. 3-D inversion models confirm the existence of a subvertical high conductivity structure underneath the DSB but its expression is far less pronounced. Instead, the 3-D inversion model suggests a deepening of the conductive DSB sediments off-profile towards the south, reaching a maximum depth of approximately 12 km, which is consistent with other geophysical observations. At shallower levels, the 3-D inversion model reveals salt diapirism as an upwelling of highly resistive structures, localized underneath the Al-Lisan Peninsula. The 3-D model furthermore contains an E–W elongated conductive structure to the northeast of the DSB. More MT data with better spatial coverage are required, however, to fully constrain the robustness of the above-mentioned off-profile features.

Description: Estimating how topography is maintained provides insights into the different factors responsible for surface deformations and their relative roles. Here, we develop a new and simple approach to assess the degree of isostatic compensation of continental topography at regional scale from GOCE gravity gradients. We calculate the ratio between the radial gradient observed by GOCE and that calculated from topography only. From analytical and statistical formulations, simple relationships between this ratio and the degree of compensation are obtained under the Airy–Heiskanen isostasy hypothesis. Then, a value of degree of compensation at each point of study area can be easily deduced. We apply our method to the Alaska-Canada Cordillera and validate our results by comparison with a standard isostatic gravity anomaly model and additional geophysical information for this area. Both our GOCE-based results and the isostatic anomaly show that Airy–Heiskanen isostasy prevails for the Yukon Plateau whereas additional mechanisms are required to support topography below the Northwest Territories Craton and the Yakutat collision zone.

Description: Bayesian sampling based inversions require many thousands or even millions of forward models, depending on how nonlinear or non-unique the inverse problem is, and how many unknowns are involved. The result of such a probabilistic inversion is not a single ‘best-fit’ model, but rather a probability distribution that is represented by the entire model ensemble. Often, a geophysical inverse problem is non-unique, and the corresponding posterior distribution is multimodal, meaning that the distribution consists of clusters with similar models that represent the observations equally well. In these cases, we would like to visualize the characteristic model properties within each of these clusters of models. However, even for a moderate number of inversion parameters, a manual appraisal for a large number of models is not feasible. This poses the question whether it is possible to extract end-member models that represent each of the best-fit regions including their uncertainties. Here, I show how a machine learning tool can be used to characterize end-member models, including their uncertainties, from a complete model ensemble that represents a posterior probability distribution. The model ensemble used here results from a nonlinear geodynamic inverse problem, where rheological properties of the lithosphere are constrained from multiple geophysical observations. It is demonstrated that by taking vertical cross-sections through the effective viscosity structure of each of the models, the entire model ensemble can be classified into four end-member model categories that have a similar effective viscosity structure. These classification results are helpful to explore the non-uniqueness of the inverse problem and can be used to compute representative data fits for each of the end-member models. Conversely, these insights also reveal how new observational constraints could reduce the non-uniqueness. The method is not limited to geodynamic applications and a generalized MATLAB code is provided to perform the appraisal analysis.

Description: Inference of the mantle viscosity from observations for glacial isostatic adjustment (GIA) process has usually been conducted through the analyses based on the simple three-layer viscosity model characterized by lithospheric thickness, upper- and lower-mantle viscosities. Here, we examine the viscosity structures for the simple three-layer viscosity model and also for the two-layer lower-mantle viscosity model defined by viscosities of 670, D (670- D km depth) and D ,2891 ( D -2891 km depth) with D -values of 1191, 1691 and 2191 km. The upper-mantle rheological parameters for the two-layer lower-mantle viscosity model are the same as those for the simple three-layer one. For the simple three-layer viscosity model, rate of change of degree-two zonal harmonics of geopotential due to GIA process (GIA-induced 2 ) of –(6.0–6.5)  x  10 –11 yr –1 provides two permissible viscosity solutions for the lower mantle, (7–20)  x  10 21 and (5–9)  x  10 22  Pa s, and the analyses with observational constraints of the 2 and Last Glacial Maximum (LGM) sea levels at Barbados and Bonaparte Gulf indicate (5–9)  x  10 22  Pa s for the lower mantle. However, the analyses for the 2 based on the two-layer lower-mantle viscosity model only require a viscosity layer higher than (5–10)  x  10 21  Pa s for a depth above the core–mantle boundary (CMB), in which the value of (5–10)  x  10 21  Pa s corresponds to the solution of (7–20)  x  10 21  Pa s for the simple three-layer one. Moreover, the analyses with the 2 and LGM sea level constraints for the two-layer lower-mantle viscosity model indicate two viscosity solutions: 670,1191  〉 3  x  10 21 and 1191,2891  ~ (5–10)  x  10 22  Pa s, and 670,1691  〉 10 22 and 1691,2891  ~ (5–10)  x  10 22  Pa s. The inferred upper-mantle viscosity for such solutions is (1–4)  x  10 20  Pa s similar to the estimate for the simple three-layer viscosity model. That is, these analyses require a high viscosity layer of (5–10)  x  10 22  Pa s at least in the deep mantle, and suggest that the GIA-based lower-mantle viscosity structure should be treated carefully in discussing the mantle dynamics related to the viscosity jump at ~670 km depth. We also preliminarily put additional constraints on these viscosity solutions by examining typical relative sea level (RSL) changes used to infer the lower-mantle viscosity. The viscosity solution inferred from the far-field RSL changes in the Australian region is consistent with those for the 2 and LGM sea levels, and the analyses for RSL changes at Southport and Bermuda in the intermediate region for the North American ice sheets suggest the solution of 670, D  〉 10 22 , D ,2891  ~ (5–10)  x  10 22  Pa s ( D  = 1191 or 1691 km) and upper-mantle viscosity higher than 6  x  10 20  Pa s.

Description: Lacticin Q is an unmodified leaderless bacteriocin produced by Lactococcus lactis QU 5. It has been revealed that the production and self-immunity of lacticin Q are facilitated by a gene cluster lnqQBCDEF . The gene for a putative TetR-family transcriptional regulator, termed lnqR , was found nearby the lnqQBCDEF cluster, but its involvement in lacticin Q biosynthesis remained unknown. In this study, we created an LnqR-overexpressing QU 5 recombinant by using lactococcal constitutive promoter P 32 . The recombinant QU 5 showed enhanced production of and self-immunity to lacticin Q. RT-PCR analysis has revealed that an overexpression of LnqR increases the amounts of lnqQBCDEF transcripts, and these six genes are transcribed as an operon in a single transcriptional unit. Interestingly, LnqR expression and thus lacticin Q production by L. lactis QU 5 was found temperature dependent, while LnzR, an LnqR-homologue, in L. lactis QU 14 was expressed in a similar but not identical manner to LnqR, resulting in dissimilar bacteriocin productivities by these strains. This report demonstrates LnqR as the first TetR-family transcriptional regulator involved in LAB bacteriocin biosynthesis and that, as an exceptional case of TetR-family regulators, LnqR positively regulates the transcription of these biosynthetic genes.

Description: Mid-crustal conductors are a common phenomenon in magnetotelluric studies. In the Andean Cordillera of southern Chile, they appear to concentrate along major fault zones. A high-resolution, broad-band magnetotelluric survey including 31 stations has been carried out along two profiles perpendicular to (1) the Liquiñe-Ofqui Fault Systems (LOFS) and (2) the Villarrica-Quetrupillán-Lanín volcanic lineament running parallel to the Mocha-Villarrica Fault Zone (MVFZ). The survey aimed at tracing one of the known conductors from mid-crustal depth to near-surface along these faults. Directionality and dimensionality were analysed using tensor decomposition. Phase tensors and induction arrows reveal two major geoelectric strike directions following the strike of LOFS and MVFZ. 2-D inversion shows low resistivity zones along both fault systems down to a depth of 〉10 km, where the brittle-ductile transition is expected. Along the LOFS, the two anomalies are linked to (1) Lake Caburgua, where the LOFS broadens to about 2 km of lateral extension and seems to represent a pull-apart structure, and (2) the intersection with the Villarrica-Quetrupillán-Lanín volcanic lineament, where seismic activity was observed during the latest eruption in March 2015. A connection of the mid-crustal conductor to the ESE-WNW-striking fault zones is indicated from the presented data.

Description: Cytosolic, globular droplets with an average diameter of 50 nm were observed in vitrified Thermoplasma acidophilum cells by means of cryo-electron tomography. These droplets were isolated by column chromatography and immunoprecipitation protein purification methods. Subsequent chemical and biochemical analyses identified lipid and protein components, respectively. Two major lipid components, comigrating menaquinones at the solvent front and the slower migrating Thermoplasma polar lipid U4, were detected by TLC experiments. The major protein component was identified as the 153 amino acid long Ta0547 vitellogenin-N domain protein. This domain has been found so far exclusively in large lipid transport proteins of vertebrates and non-vertebrates. Blast protein database homology searches with Ta0547 did not return any eukaryal hits; homologous sequences were found only in thermo-acidophilic archaeons. However, a profile-sequence domain search performed with the vitellogenin-N domain (PF01347) hmm-profile against the T. acidophilum proteome returned Ta0547 as hit. Electron microscopy appearance of isolated droplets resembled to lipoprotein particles. However, no (tetraether) lipid layer could be detected on the droplets surface, rather hydrophobic compounds of the electron dense lumen were surrounded by a denser discontinuous protein boundary. Based on described features, these particles qualify for a novel lipoprotein particle category, what we nominated Thermoplasma Quinone Droplet.

Description: We present new rotations that describe the relative positions and velocities of the Pacific and North America plates at 22 times during the past 19.7 Myr, offering 1-Myr temporal resolution for studies of the geotectonic evolution of western North America and other plate boundary locations. Derived from 18 000 magnetic reversal, fracture zone and transform fault identifications from the Pacific–Antarctic–Nubia–North America plate circuit and the velocities of 935 GPS sites on the Pacific and North America plates, the new rotations and GPS-derived angular velocity indicate that the rate of motion between the two plates increased by 70 per cent from 19.7 to 9±1 Ma, but changed by less than 2 per cent since 8 Ma and even less since 4.2 Ma. The rotations further suggest that the relative plate direction has rotated clockwise for most of the past 20 Myr, with a possible hiatus from 9 to 5 Ma. This conflicts with previously reported evidence for a significant clockwise change in the plate direction at 8–6 Ma. Our new rotations indicate that Pacific plate motion became obliquely convergent with respect to the San Andreas Fault of central California at 5.2–4.2 Ma, in agreement with geological evidence for a Pliocene onset of folding and faulting in central California. Our reconstruction of the northern Gulf of California at 6.3 Ma differs by only 15–30 km from structurally derived reconstructions after including 3–4 km Myr –1 of geodetically measured slip between the Baja California Peninsula and Pacific plate. This implies an approximate 15–30 km upper bound for plate non-rigidity integrated around the global circuit at 6.3 Ma. A much larger 200±54 km discrepancy between our reconstruction of the northern Gulf of California at 12 Ma and that estimated from structural and marine geophysical observations suggests that faults in northwestern Mexico or possibly west of the Baja California Peninsula accommodated large amounts of obliquely divergent dextral shear from 12–6.3 Ma. Pacific–North America plate motion since 16 Myr estimated with our new rotations agrees well with structurally summed deformation along two transects of western North America between the Colorado Plateau and western California, with a difference as small as 40 km out of 760 km of margin-parallel motion. A strong resemblance between a 20-Myr-to-present flow line reconstructed with our new rotations and the traces of the 700-km-long Queen Charlotte Fault and continental slope west of Canada suggests that the plate margin geometry was influenced by the passage of the Pacific plate and Yakutat block. The new rotations also suggest that (1) oblique convergence west of Canada initiated at 12–11 Ma, 5–8 Myr earlier than previously estimated, (2) no significant margin-normal shortening has occurred in areas of Canada located east of the Haida Gwaii archipelago since 20 Ma and (3) Pacific plate underthrusting of Haida Gwaii has accommodated the margin-normal component of plate motion since 12–11 Ma. Our rotations suggest an 70 per cent increase in the rate that the Pacific plate has been consumed by subduction beneath the Aleutian arc since 19.7 Ma, with still-unknown consequences for the rate of arc magmatism.

Description: Glyoxal (GO) belongs to the reactive electrophilic species generated in vivo in all organisms. In order to identify targets of GO and their response mechanisms, we attempted to screen for GO-sensitive mutants by random insertions of Tn phoA -132. The genes responsible for GO susceptibility were functionally classified as the following: (i) tRNA modification; trmE , gidA and truA , (ii) DNA repair; recA and recC , (iii) toxin–antitoxin; mqsA and (iv) redox metabolism; yqhD and caiC . In addition, an insertion in the crp gene, encoding the cAMP responsive transcription factor, exhibits a GO-resistant phenotype, which is consistent with the phenotype of adenylate cyclase ( cya ) mutant showing GO resistance. This suggests that global regulation involving cAMP is operated in a stress response to GO. To further characterize the CRP-regulated genes directly associated with GO resistance, we created double mutants deficient in both crp and one of the candidate genes including yqhD , gloA and sodB . The results indicate that these genes are negatively regulated by CRP as confirmed by real-time RT-PCR. We propose that tRNA as well as DNA are the targets of GO and that toxin/antitoxin, antioxidant and cAMP are involved in cellular response to GO.

Description: Described are results of laboratory experiments which revealed regularities of gradual transition from stick-slip mode to aseismic creep. The behaviour of model gouge-filled fault was investigated with experimental setup of the spring-bock model. It was experimentally proven that small variations of a percentage of materials with velocity strengthening and velocity weakening properties in the fault principal slip zone may result in significant variation of the portion of seismic energy radiated during a fault slip event. The tests simulated different modes of interblock sliding whose characteristic values of scaled kinetic energy varied by several orders of magnitude, while differences in contact strength and shear stress drop remained relatively small. The obtained results led to the conclusion that the earthquake radiation efficiency and the fault slip mode are governed by the ratio of two parameters—maximum fault slip-weakening rate and shear stiffness of the enclosing massif. The ratio can be essentially changed by small variations of the material composition of the fault principal slip zone.

Description: Among the outstanding tectonic questions regarding the convergence between the Tien Shan and Tarim basin in northwestern China are the manner in which deformation is accommodated within their lithospheres, and the extent that the Tarim lithosphere underthrusts the Tien Shan. In particular, the amount and type of deformation within the Tarim basin is poorly understood. It is also uncertain if the convergence between the Tarim and the Tien Shan takes place mainly along a discrete boundary, or if the Tarim lithosphere simply indents into the Kazach shield, forming the Tien Shan through crustal thickening accommodated by a distributed series of thrust faults. In this study we use hypocentres from published earthquake catalogues and waveforms recorded by regional seismic networks to determine earthquake source parameters through regional centroid moment tensor inversion. The entire dataset consists of 160 earthquakes that occurred between 1969 and 2009 and with moment magnitudes between 3.5 and 7 distributed throughout the central Tien Shan and northwestern Tarim Basin. The estimated focal depths of these earthquakes range from the near-surface to about 44 km. Focal mechanisms throughout much of the Tien Shan indicate active deformation accommodated by thrust faults from at least the upper crust to 30 km depth. South of the Tien Shan, the Jia-shi earthquake sequence within the Tarim basin suggests that both crustal shortening and localized flexure are part of a complicated process involving rotational convergence. Inside the Tarim basin, two earthquakes with thrust faulting mechanisms near the crust–mantle boundary beneath the Bachu uplift imply a brittle rheology of the lower crust. High-angle thrust events occur broadly across the Tien Shan, suggesting that the Tarim lithosphere as a whole is strong and indents into the Kazach shield to create the mountain range.

Description: The subduction of lithospheric plates is partitioned between subducting plate motion and lateral slab migration (i.e. trench retreat and advance). We use 3-D, dynamic models of subduction to address the role of a power-law mantle viscosity on subduction dynamics and, in particular, rates of trench retreat. For all numerical models tested, we find that a power-law rheology results in reduced rates of trench retreat, and elevated slab dip angles, relative to the equivalent isoviscous mantle model. We analyse the asthenospheric pressure distribution and the style of mantle flow, which exhibits only limited variability as a function of mantle rheology, in order to compute estimates of the mantle forces associated with subduction. The inclusion of a power-law rheology reduces the mantle shear force (which resists subducting plate motion) to a greater degree than it reduces the dynamic pressure gradient across the slab (which resists trench retreat). Therefore, the inclusion of a power-law mantle rheology favours a shift towards a subduction mode with a reduced trench retreat component, typically a relative reduction of order 25 per cent in our 3-D models. We suggest that this mechanism may be of importance for reducing the high trench retreat rates observed in many previous models to levels more in line with the average subduction partitioning observed on Earth at present (i.e. trench velocity ≤ plate velocity), for most absolute plate motion reference frames.

Description: Pulverized rocks (PR) are extremely incohesive and highly fractured rocks found within the damage zones of several large strike-slip faults around the world. They maintain their crystal structure, show little evidence of shearing or chemical alteration, and are believed to be produced by strong tensile forces. Several mechanisms for pulverization have been proposed based on simple qualitative analyses or laboratory experiments under simplified loading conditions. Numerical modelling, however, can offer new insights into what is needed to produce PR and likely conditions of formation. We perform dynamic rupture simulations of different earthquakes, varying the magnitude, the slip distribution, and the rupture speed (supershear and subshear), while measuring the stresses produced away from the fault. To contextualize our results, a basic threshold of 10 MPa is set as the tensile strength of the rock mass and recordings are made of where, when, and by how much this threshold is exceeded for each earthquake type. Guided by field observations, we discern that a large (〉 M w 7.1) subshear earthquake along a bimaterial fault produces a pulverized rock distribution most consistent with observations. The damage is asymmetric with the majority on the stiffer side of the fault extending out for several hundred metres. Within this zone there is a large and sudden volumetric expansion in all directions as the rupture passes. We propose that such an extreme tensile stress state, repeated for every earthquake, eventually produces the PR seen in the field.