Population subdivision in marine environments: the contributions of biogeography, geographic distance, and discontinuous habitat to genetic differentiation in a blennioid

The relative importance of factors that may promote genetic differentiation in marine organisms is largely unknown. Here, contributions to population structure from biogeography, habitat distribution, and isolation by distance were investigated in Axoclinus nigricaudus, a small subtidal rock reef fish, throughout its range in the Gulf of California. A 408 basepair fragment of the mitochondrial control region was sequenced from 105 individuals. Variation was significantly partitioned between many pairs of populations. Phylogenetic analyses, hierarchical analyses of variance, and general linear models substantiated a major break between two putative biogeographic regions. This genetic discontinuity coincides with an abrupt change in ecological characteristics (including temperature and salinity) but does not coincide with known oceanographic circulation patterns. Geographic distance and the nature of habitat separating populations (continuous habitat along a shoreline, discontinuous habitat along a shoreline, and open water) also contributed to population structure in general linear model analyses. To verify that local populations are genetically stable over time, one population was resampled on four occasions over eighteen months; it showed no evidence of a temporal component to diversity. These results indicate that having a planktonic life stage does not preclude geographically partitioned genetic variation over relatively small geographic distances in marine environments. Moreover, levels of genetic differentiation among populations of Axoclinus nigricaudus cannot be explained by a single factor, but are due to the combined influences of a biogeographic boundary, habitat, and geographic distance.

Riddling and invariance for discontinuous maps preserving Lebesgue measure

In this paper we use the mixture of topological and measure-theoretic dynamical approaches to consider riddling of invariant sets for some discontinuous maps of compact regions of the plane that preserve two-dimensional Lebesgue measure. We consider maps that are piecewise continuous and with invertible except on a closed zero measure set. We show that riddling is an invariant property that can be used to characterize invariant sets, and prove results that give a non-trivial decomposion of what we call partially riddled invariant sets into smaller invariant sets. For a particular example, a piecewise isometry that arises in signal processing (the overflow oscillation map), we present evidence that the closure of the set of trajectories that accumulate on the discontinuity is fully riddled. This supports a conjecture that there are typically an infinite number of periodic orbits for this system.

Magma Flow Instabilities in a Volcanic Conduit: Implications for Long-Period Seismicity

Silicic volcanic eruptions are typically accompanied by repetitive Long-Period (LP) seismicity that originates from a small region of the upper conduit. These signals have the capability to advance eruption prediction, since they commonly precede a change in the eruption vigour. Shear bands forming along the conduit wall, where the shear stresses are highest, have been linked to providing the seismic trigger. However, existing computational models are unable to generate shear bands at the depths where the LP signals originate using simple magma strength models. Presented here is a model in which the magma strength is determined from a constitutive relationship dependent upon crystallinity and pressure. This results in a depth-dependent magma strength, analogous to planetary lithospheres. Hence, in shallow highly-crystalline regions a macroscopically discontinuous brittle type of deformation will prevail, whilst in deeper crystal-poor regions there will be a macroscopically continuous plastic deformation mechanism. This will result in a depth where the brittle-ductile transition occurs, and here shear bands disconnected from the free-surface may develop. We utilize the Finite Element Method and use axi-symmetric coordinates to model magma flow as a viscoplastic material, simulating quasi-static shear bands along the walls of a volcanic conduit. Model results constrained to the Soufrière Hills Volcano, Montserrat, show the generation of two types of shear bands: upper-conduit shear bands that form between the free-surface to a few 100 metres below it and discrete shear bands that form at the depths where LP seismicity is measured to occur corresponding to the brittle-ductile transition and the plastic shear region. It is beyond the limitation of the model to simulate a seismic event, although the modelled viscosity within the discrete shear bands suggests a failure and healing cycle time that supports the observed LP seismicity repeat times. However, due to the paucity of data and large parameter space available these results can only be considered to be qualitative rather than quantitative at this stage.

Fast evolution strategies

Evolution strategies are a class of general optimisation algorithms which are applicable to functions that are multimodal, nondifferentiable, or even discontinuous. Although recombination operators have been introduced into evolution strategies, the primary search operator is still mutation. Classical evolution strategies rely on Gaussian mutations. A new mutation operator based on the Cauchy distribution is proposed in this paper. It is shown empirically that the new evolution strategy based on Cauchy mutation outperforms the classical evolution strategy on most of the 23 benchmark problems tested in this paper. The paper also shows empirically that changing the order of mutating the objective variables and mutating the strategy parameters does not alter the previous conclusion significantly, and that Cauchy mutations with different scaling parameters still outperform the Gaussian mutation with self-adaptation. However, the advantage of Cauchy mutations disappears when recombination is used in evolution strategies. It is argued that the search step size plays an important role in determining evolution strategies' performance. The large step size of recombination plays a similar role as Cauchy mutation.

Observation of a cubic phase in mesoporous silicate films grown at the air/water interface

X-Ray diffraction is reported from mesoporous silicate films grown at the air/water interface. The films were studied both as powdered films, and oriented on silicon or mica sheets. At early stages of growth we observe Bragg diffraction from a highly ordered cubic phase, with both long and short d-spacing peaks. We have assigned this as a discontinuous micellar Pm3n phase in which the silica is partly ordered. Later films retain only the known hexagonal p6m peaks and have lost any order both at short d-spacings and the longer d-spacing Bragg peaks characteristic of the cubic structure. The silica framework is considerably expanded from that in bulk amorphous silica, average Si Si distances are some 30% greater. Incorporation of glycerol or polyethylene glycol preserves the earlier cubic structure. To be consistent with earlier, in situ, X-ray and neutron reflectivity data we infer that both structures are produced after a phase transition from a less-ordered him structure late in the induction phase. The structural relations between the film Pm3n and p6m phase(s) and the known bulk SBA-1 and MCM-41 phases are briefly discussed.

Identification of a new ligand binding domain in the al subunit of the inhibitory glycine receptor

Four discontinuous extracellular sequence domains have been proposed to form the ligand binding sites of the ligand-gated ion channel receptor superfamily. In this study, we investigated the role of 12 contiguous residues of the inhibitory glycine receptor that define the proposed loop A ligand binding domain; Using the techniques of site-directed mutagenesis and patch-clamp electrophysiology, four of the 12 residues were shown to have impaired ligand binding. Three mutants, I93A, A101H, and N102A, resulted in significant (17-44-fold) increases in the agonist EC50 values as compared with the wild-type glycine receptor, whereas Hill coefficients, I-max values, and antagonist affinity remained largely unaffected. Consideration of receptor efficacy values indicates that these residues are involved in ligand binding rather than channel activation. A fourth mutant, W94A, failed to give rise to any glycine-activated currents, although cell-surface expression was observed, suggesting that this residue may also be involved in agonist binding. These data provide the most extensive characterization of the loop A ligand binding domain available to date and define two new residue locations, Ile(93) and Asn(102), as contributing to the four-loop model of ligand binding.

Low intraspecific variation in the rRNA internal transcribed spacer 2 (ITS2) of the Australian paralysis tick, Ixodes holocyclus

Ixodes holocyclus has a narrow, discontinuous distribution along the east coast of Australia. We studied ticks from 17 localities throughout the geographic range of this tick. The ITS2 of I. holocyclus is 793 bp long. We found nucleotide variation at eight of the 588 nucleotide positions (1.4%) that were compared for all ticks. There were eight different nucleotide sequences. Most sequences were not restricted to a particular geographic region. However, sequences F, G and H, which had an adenine at position 197, were found only in the far north of Queensland - all other ticks had a guanine at this position. The low level of intraspecific variation in this tick (0.7%) contrasts with the sequence divergence between L holocyclus and its close relative, I. cornuatus (13.1 %). These data indicate that L holocyclus does not contain cryptic species despite possible geographic isolation of some populations. We conclude that variation in the ITS2 is likely to be informative about the phylogeny of the group.

Adsorption in slit-like pores of activated carbons: Improvement of the Horvath and Kawazoe method

Anew thermodynamic approach has been developed in this paper to analyze adsorption in slitlike pores. The equilibrium is described by two thermodynamic conditions: the Helmholtz free energy must be minimal, and the grand potential functional at that minimum must be negative. This approach has led to local isotherms that describe adsorption in the form of a single layer or two layers near the pore walls. In narrow pores local isotherms have one step that could be either very sharp but continuous or discontinuous benchlike for a definite range of pore width. The latter reflects a so-called 0 --> 1 monolayer transition. In relatively wide pores, local isotherms have two steps, of which the first step corresponds to the appearance of two layers near the pore walls, while the second step corresponds to the filling of the space between these layers. All features of local isotherms are in agreement with the results obtained from the density functional theory and Monte Carlo simulations. The approach is used for determining pore size distributions of carbon materials. We illustrate this with the benzene adsorption data on activated carbon at 20, 50, and 80 degreesC, argon adsorption on activated carbon Norit ROX at 87.3 K, and nitrogen adsorption on activated carbon Norit R1 at 77.3 K.

Stress correlation function evolution in lattice solid elasto-dynamic models of shear and fracture zones and earthquake prediction

It has been argued that power-law time-to-failure fits for cumulative Benioff strain and an evolution in size-frequency statistics in the lead-up to large earthquakes are evidence that the crust behaves as a Critical Point (CP) system. If so, intermediate-term earthquake prediction is possible. However, this hypothesis has not been proven. If the crust does behave as a CP system, stress correlation lengths should grow in the lead-up to large events through the action of small to moderate ruptures and drop sharply once a large event occurs. However this evolution in stress correlation lengths cannot be observed directly. Here we show, using the lattice solid model to describe discontinuous elasto-dynamic systems subjected to shear and compression, that it is for possible correlation lengths to exhibit CP-type evolution. In the case of a granular system subjected to shear, this evolution occurs in the lead-up to the largest event and is accompanied by an increasing rate of moderate-sized events and power-law acceleration of Benioff strain release. In the case of an intact sample system subjected to compression, the evolution occurs only after a mature fracture system has developed. The results support the existence of a physical mechanism for intermediate-term earthquake forecasting and suggest this mechanism is fault-system dependent. This offers an explanation of why accelerating Benioff strain release is not observed prior to all large earthquakes. The results prove the existence of an underlying evolution in discontinuous elasto-dynamic, systems which is capable of providing a basis for forecasting catastrophic failure and earthquakes.

Epitope map for a growth hormone receptor agonist monoclonal antibody, MAb 263

Monoclonal antibody (MAb) 263 is a widely used monoclonal antibody that recognizes the extracellular domain (ECD) of the GH receptor. It has been shown to act as a GH agonist both in vitro and in vivo, and we report here that it must be divalent to exert its effect on the full-length receptor. To understand the mechanism of its agonist action, we have determined the precise epitope for this antibody using a novel random PCR mutagenesis approach together with expression screening in yeast. A library of 5200 clones of rabbit GH receptor ECD mutants were screened both with MAb 263 and with an anticarboxy-tag antibody to verify complete ECD expression. Sequencing for clones that expressed complete ECD but were not MAb 263 positive identified 20 epitope residues distributed in a discontinuous manner throughout the ECD. The major part of the epitope, as revealed after mapping onto the crystal structure model of the ECD molecule, was located on the side and upper portion of domain 1, particularly within the D - E strand disulfide loop 79 - 96. Molecular dynamics docking of an antibody of the same isotype as MAb 263 was used to dock the bivalent antibody to the 1528-Angstrom(2) epitope and to visualize the likely consequences of MAb binding. The minimized model enables the antibody to grasp two receptors in a pincer-like movement from opposite sides, facilitating alignment of the receptor dimerization domains in a manner similar to, but not identical with, GH.