Inferring epidemic contact structure from phylogenetic trees.

Contact structure is believed to have a large impact on epidemic spreading and consequently using networks to model such contact structure continues to gain interest in epidemiology. However, detailed knowledge of the exact contact structure underlying real epidemics is limited. Here we address the question whether the structure of the contact network leaves a detectable genetic fingerprint in the pathogen population. To this end we compare phylogenies generated by disease outbreaks in simulated populations with different types of contact networks. We find that the shape of these phylogenies strongly depends on contact structure. In particular, measures of tree imbalance allow us to quantify to what extent the contact structure underlying an epidemic deviates from a null model contact network and illustrate this in the case of random mixing. Using a phylogeny from the Swiss HIV epidemic, we show that this epidemic has a significantly more unbalanced tree than would be expected from random mixing.

Quantitative comparison of reconstruction methods for intra-voxel fiber recovery from diffusion MRI.

Validation is arguably the bottleneck in the diffusion magnetic resonance imaging (MRI) community. This paper evaluates and compares 20 algorithms for recovering the local intra-voxel fiber structure from diffusion MRI data and is based on the results of the "HARDI reconstruction challenge" organized in the context of the "ISBI 2012" conference. Evaluated methods encompass a mixture of classical techniques well known in the literature such as diffusion tensor, Q-Ball and diffusion spectrum imaging, algorithms inspired by the recent theory of compressed sensing and also brand new approaches proposed for the first time at this contest. To quantitatively compare the methods under controlled conditions, two datasets with known ground-truth were synthetically generated and two main criteria were used to evaluate the quality of the reconstructions in every voxel: correct assessment of the number of fiber populations and angular accuracy in their orientation. This comparative study investigates the behavior of every algorithm with varying experimental conditions and highlights strengths and weaknesses of each approach. This information can be useful not only for enhancing current algorithms and develop the next generation of reconstruction methods, but also to assist physicians in the choice of the most adequate technique for their studies.

Enhancing actions of peptides derived from the γ-chain of fetal human hemoglobin on the immunostimulant activities of monophosphoryl lipid A.

Hemoglobin and its structures have been described since the 1990s to enhance a variety of biological activities of endotoxins (LPS) in a dose-dependent manner. To investigate the interaction processes in more detail, the system was extended by studying the interactions of newly designed peptides from the γ-chain of human hemoglobin with the adjuvant monophosphoryl lipid A (MPLA), a partial structure of lipid A lacking its 1-phosphate. It was found that some selected Hbg peptides, in particular two synthetic substructures designated Hbg32 and Hbg35, considerably increased the bioactivity of MPLA, which alone was only a weak activator of immune cells. These findings hold true for human mononuclar cells, monocytes and T lymphocytes. To understand the mechanisms of action in more detail, biophysical techniques were applied. These showed a peptide-induced change of the MPLA aggregate structure from multilamellar into a non-lamellar, probably inverted, cubic structure. Concomitantly, the peptides incorporated into the tightly packed MPLA aggregates into smaller units down to monomers. The fragmentation of the aggregates was an endothermic process, differing from a complex formation but rather typical for a catalytic reaction.

Reproductive conflicts and egg discrimination in a socially polymorphic ant

The ability to discriminate against competitors shapes cooperation and conflicts in all forms of social life. In insect societies, workers may detect and destroy eggs laid by other workers or by foreign queens, which can contribute to regulate reproductive conflicts among workers and queens. Variation in colony kin structure affects the magnitude of these conflicts and the diversity of cues used for discrimination, but the impact of the number of queens per colony on the ability of workers to discriminate between eggs of diverse origin has so far not been investigated. Here, we examined whether workers from the socially polymorphic ant Formica selysi distinguished eggs laid by nestmate workers from eggs laid by nestmate queens, as well as eggs laid by foreign queens from eggs laid by nestmate queens. Workers from single- and multiple-queen colonies discriminated worker-laid from queen-laid eggs, and eliminated the former. This suggests that workers collectively police each other in order to limit the colony-level costs of worker reproduction and not because of relatedness differences towards queens' and workers' sons. Workers from single-queen colonies discriminated eggs laid by foreign queens of the same social structure from eggs laid by nestmate queens. In contrast, workers from multiple-queen colonies did not make this distinction, possibly because cues on workers or eggs are more diverse. Overall, these data indicate that the ability of F. selysi workers to discriminate eggs is sufficient to restrain worker reproduction but does not permit discrimination between matrilines in multiple-queen colonies.

Quadrant/octant sequencing and the role of coherent structures in bed load sediment entrainment

To permit the tracking of turbulent flow structures in an Eulerian frame from single-point measurements, we make use of a generalization of conventional two-dimensional quadrant analysis to three-dimensional octants. We characterize flow structures using the sequences of these octants and show how significance may be attached to particular sequences using statistical mull models. We analyze an example experiment and show how a particular dominant flow structure can be identified from the conditional probability of octant sequences. The frequency of this structure corresponds to the dominant peak in the velocity spectra and exerts a high proportion of the total shear stress. We link this structure explicitly to the propensity for sediment entrainment and show that greater insight into sediment entrainment can be obtained by disaggregating those octants that occur within the identified macroturbulence structure from those that do not. Hence, this work goes beyond critiques of Reynolds stress approaches to bed load entrainment that highlight the importance of outward interactions, to identifying and prioritizing the quadrants/octants that define particular flow structures. Key Points <list list-type=''bulleted'' id=''jgrf20196-list-0001''> <list-item id=''jgrf20196-li-0001''>A new method for analysing single point velocity data is presented <list-item id=''jgrf20196-li-0002''>Flow structures are identified by a sequence of flow states (termed octants) <list-item id=''jgrf20196-li-0003''>The identified structure exerts high stresses and causes bed-load entrainment

Pizgrischite, (Cu,Fe)Cu14PbBi17S35, a new sulfosalt from the Swiss Alps: Description, crystal structure and occurrence

Pizgrischite, (Cu,Fe)Cu14PbBi17S35, is a new mineral species named after the type locality, Piz Grisch Mountain, Val Ferrera, Graubunden, Switzerland. This sulfosalt occurs as thin, striated, metallic lead-grey blades measuring up to I cm in length, embedded in quartz and associated with tetrahedrite, chalcopyrite, pyrite, sphalerite, emplectite and derivatives of the aikinite-bismuthinite series. In plane-polarized light, the new species is brownish grey with no perceptible pleochroism; under crossed nicols in oil immersion, it presents a weak anisotropy with dark brown tints. Minimum and maximum reflectance values (in %) in air are: 40.7-42.15 (470 nm), 41.2-43.1 (546 nm), 41.2-43.35 (589 nm) and 40.7-43.3 (650 nm). Cleavage is perfect along 001 I and well developed on {010}. Abundant polysynthetic twinning is observed on (010). The mean micro-indentation hardness is 190 kg/mm(2) (Mohs hardness 3.3), and the calculated density is 6.58 g/cm(3). Electron-microprobe analyses yield (wt%; mean result of seven analyses): Cu 16.48, Pb 2.10, Fe 0.77, Bi 60.70, Sb 0.35, S 19.16, Se 0.04, total 99.60. The resulting empirical chemical formula is (Cu15.24Fe0.80Pb0.60)(Sigma 16.64)(Bi17.07Sb0.17)(Sigma 17.24)(S35.09Se0.03)(Sigma 35.12), in accordance with the formula derived from the single-crystal refinement of the structure, (Cu,Fe)Cu14PbBi17S35. Pizgrischite is monoclinic, space group C2/m, with the following unit-cell parameters: a 35.054(2), b3.91123(I), c43.192(2) angstrom, beta 96.713(4)degrees, V5881.24 angstrom(3), Z=4. The strongest seven X-ray powder-diffraction lines [d in angstrom (I)(hkl)] are: 5.364(40)((6) over bar 04), 4.080(50)((8) over bar 05), 3.120(40)(118), 3.104(68)((3) over bar 18), 2.759(53) ((9) over bar 11),2.752(44)(910) and 1.956(100)(020). The crystal structure is an expanded monoclinic derivative of kupcikite. Pizgrischite belongs to the cuprobismutite series of bismuth sulfosalts but, sensu stricto, it is not a homologue of cuprobismutite. At the type locality. pizarischite is the result of the Alpine metamorphism under greenschist-facies conditions of pre-Tertiary hydrothermal Cu-Bi mineralization.

The twist, writhe and overall shape of supercoiled DNA change during counterion-induced transition from a loosely to a tightly interwound superhelix. Possible implications for DNA structure in vivo.

A cryo-electron microscopy study of supercoiled DNA molecules freely suspended in cryo-vitrified buffer was combined with Monte Carlo simulations and gel electrophoretic analysis to investigate the role of intersegmental electrostatic repulsion in determining the shape of supercoiled DNA molecules. It is demonstrated here that a decrease of DNA-DNA repulsion by increasing concentrations of counterions causes a higher fraction of the linking number deficit to be partitioned into writhe. When counterions reach concentrations likely to be present under in vivo conditions, naturally supercoiled plasmids adopt a tightly interwound conformation. In these tightly supercoiled DNA molecules the opposing segments of interwound superhelix seem to directly contact each other. This form of supercoiling, where two DNA helices interact laterally, may represent an important functional state of DNA. In the particular case of supercoiled minicircles (178 bp) the delta Lk = -2 topoisomers undergo a sharp structural transition from almost planar circles in low salt buffers to strongly writhed "figure-eight" conformations in buffers containing neutralizing concentrations of counterions. Possible implications of this observed structural transition in DNA are discussed.

Estimation of the lateral correlation structure of subsurface water content from surface-based ground-penetrating radar reflection images

Over the past decade, significant interest has been expressed in relating the spatial statistics of surface-based reflection ground-penetrating radar (GPR) data to those of the imaged subsurface volume. A primary motivation for this work is that changes in the radar wave velocity, which largely control the character of the observed data, are expected to be related to corresponding changes in subsurface water content. Although previous work has indeed indicated that the spatial statistics of GPR images are linked to those of the water content distribution of the probed region, a viable method for quantitatively analyzing the GPR data and solving the corresponding inverse problem has not yet been presented. Here we address this issue by first deriving a relationship between the 2-D autocorrelation of a water content distribution and that of the corresponding GPR reflection image. We then show how a Bayesian inversion strategy based on Markov chain Monte Carlo sampling can be used to estimate the posterior distribution of subsurface correlation model parameters that are consistent with the GPR data. Our results indicate that if the underlying assumptions are valid and we possess adequate prior knowledge regarding the water content distribution, in particular its vertical variability, this methodology allows not only for the reliable recovery of lateral correlation model parameters but also for estimates of parameter uncertainties. In the case where prior knowledge regarding the vertical variability of water content is not available, the results show that the methodology still reliably recovers the aspect ratio of the heterogeneity.

Constraints on the permeability structure of alluvial aquifers from P-wave sonic logs

Modern sonic logging tools designed for shallow environmental and engineering applications allow for P-wave phase velocity measurements over a wide frequency band. Methodological considerations indicate that, for saturated unconsolidated sediments in the silt to sand range and source frequencies ranging from approximately 1 to 30 kHz, the observable poro-elastic P-wave velocity dispersion is sufficiently pronounced to allow for reliable first-order estimations of the underlying permeability structure. These predictions have been tested on and verified for a surficial alluvial aquifer. Our results indicate that, even without any further calibration, the thus obtained permeability estimates as well as their variabilities within the pertinent lithological units are remarkably close to those expected based on the corresponding granulometric characteristics.

Use of an in-house approach to study the three-dimensional structures of various outer membrane proteins: structure of the alcaligin outer membrane transporter FauA from Bordetella pertussis.

Bordetella pertussis is the bacterial agent of whooping cough in humans. Under iron-limiting conditions, it produces the siderophore alcaligin. Released to the extracellular environment, alcaligin chelates iron, which is then taken up as a ferric alcaligin complex via the FauA outer membrane transporter. FauA belongs to a family of TonB-dependent outer membrane transporters that function using energy derived from the proton motive force. Using an in-house protocol for membrane-protein expression, purification and crystallization, FauA was crystallized in its apo form together with three other TonB-dependent transporters from different organisms. Here, the protocol used to study FauA is described and its three-dimensional structure determined at 2.3 A resolution is discussed.

Hepatitis C virus proteins: from structure to function.

Great progress has been made over the past years in elucidating the structure and function of the hepatitis C virus (HCV) proteins, most of which are now actively being pursued as antiviral targets. The structural proteins, which form the viral particle, include the core protein and the envelope glycoproteins E1 and E2. The nonstructural proteins include the p7 viroporin, the NS2 protease, the NS3-4A complex harboring protease and NTPase/RNA helicase activities, the NS4B and NS5A proteins, and the NS5B RNA-dependent RNA polymerase. NS4B is a master organizer of replication complex formation while NS5A is a zinc-containing phosphoprotein involved in the regulation of HCV RNA replication versus particle production. Core to NS2 make up the assembly module while NS3 to NS5B represent the replication module (replicase). However, HCV proteins exert multiple functions during the viral life cycle, and these may be governed by different structural conformations and/or interactions with viral and/or cellular partners. Remarkably, each viral protein is anchored to intracellular membranes via specific determinants that are essential to protein function in the cell. This review summarizes current knowledge of the structure and function of the HCV proteins and highlights recent advances in the field.

Examining the Agreement of Structural Loss Measured With Heidelberg Retinal Tomograhy and Functional Loss With Standard Automated Perimetry (SAP; Octopus), Pulsar Perimetry (PP), and Moorfields Motion Displacement Test (MDT) Perimetry

Purpose: To examine the relationship of functional measurements with structural measures. Methods: 146 eyes of 83 test subjects underwent Heidelberg Retinal Tomography (HRTIII) (disc area<2.43, mphsd<40), and perimetry testing with Octopus (SAP; Dynamic), Pulsar (PP; TOP) and Moorfields MDT (ESTA). Glaucoma was defined as progressive structural or functional loss (20 eyes). Perimetry test points were grouped into 6 sectors based on the estimated optic nerve head angle into which the associated nerve fiber bundle enters (Garway-Heath map). Perimetry summary measures (PSM) (MD SAP/ MD PP/ PTD MDT) were calculated from the average total deviation of each measured threshold from the normal for each sector. We calculated the 95% significance level of the sectorial PSM from the respective normative data. We calculated the percentage agreement with group1 (G1), healthy on HRT and within normal perimetric limits, and group 2 (G2), abnormal on HRT and outside normal perimetric limits. We also examined the relationship of PSM and rim area (RA) in those sectors classified as abnormal by MRA (Moorfields Regression Analysis) of HRT. Results: The mean age was 65 (range= [37, 89]). The global sensitivity versus specificity of each instrument in detecting glaucomatous eyes was: MDT 80% vs. 88%, SAP 80% vs. 80%, PP 70% vs. 89% and HRT 80% vs. 79%. Highest percentage agreement of HRT (respectively G1, G2, sector) with PSM were MDT (89%, 57%, nasal superior), SAP (83%, 74%, temporal superior), PP (74%, 63%, nasal superior). Globally percentage agreement (respectively G1, G2) was MDT (92%, 28%), SAP (87%, 40%) and PP (77%, 49%). Linear regression showed there was no significant trend globally associating RA and PSM. However, sectorally the supero-nasal sector had a statistically significant (p<0.001) trend with each instrument, the associated r2 coefficients are (MDT 0.38 SAP 0.56 and PP 0.39). Conclusions: There were no significant differences in global sensitivity or specificity between instruments. Structure-function relationships varied significantly between instruments and were consistently strongest supero-nasally. Further studies are required to investigate these relationships in detail.