Recorded quadrant incidence of female breast cancer in Great Britain suggests a disproportionate increase in the upper outer quadrant of the breast

Background: The upper outer quadrant (UOQ) of the breast is the most frequent site for incidence of breast cancel; but the reported disproportionate incidence in this quadrant appears to rise with year of publication. Materials and Methods: In order to determine whether this increasing incidence in the UOQ is an artifact of different study populations or is chronological, data have been analysed for annual quadrant incidence of female breast cancer recorded nationally in England and Wales between 1979 and 2000 and in Scotland between 1980 and 2001. Results: In England and Wales, the recorded incidence of female breast cancer in the UOQ rose front 47.9% in 1979 to 53.3% in 2000, and has done so linearly over tune with a con-elation coefficient R of +/- 0.71 +/- SD 0.01 (p < 0.001). Analysis of independent data front Scotland showed a similar trend in that recorded female breast cancer had also increased in the UOQ from 38.3% in 1980 to 54.7% in 2001, with a con-elation coefficient R for the linear annual increase of +0.80 +/- SD 0.03 (p < 0.001). Conclusion: These results are inconsistent with current views that the high level of UOQ breast cancer is due solely to a greater amount of target epithelial tissue in that region. Identification of the reasons for such a disproportionate site-specific increase could provide clues as to causative factors in breast cancer.

Comparative study of oestrogenic properties of eight phytoestrogens in MCF7 human breast cancer cells

Previous studies have compared the oestrogenic properties of phytoestrogens in a wide variety of disparate assays. Since not all phytoestrogens have been tested in each assay, this makes inter-study comparisons and ranking oestrogenic potency difficult. In this report, we have compared the oestrogen agonist and antagonist activity of eight phytoestrogens (genistein, daidzein, equol, miroestrol, deoxymiroestrol, 8-prenylnaringenin, coumestrol and resveratrol) in a range of assays all based within the same receptor and cellular context of the MCF7 human breast cancer cell line. The relative binding of each phytoestrogen to oestrogen receptor (ER) of MCF7 cytosol was calculated from the molar excess needed for 50 % inhibition of [H-3]oestradiol binding (IC50), and was in the order coumestrol (35x)/8-prenylnaringenin (45 x)/deoxymiroestrol (50 x) > miroestrol (260x) > genistein (1000x) > equol (4000x) > daidzein (not achieved: 40 % inhibition at 10(4)-fold molar excess) > resveratrol (not achieved: 10 % inhibition at 10(5)-fold molar excess). For cell-based assays, the rank order of potency (estimated in terms of the concentration needed to achieve a response equivalent to 50 % of that found with 17 beta-oestradiol (IC50)) remained very similar for all the assays whether measuring ligand ability to induce a stably transfected oestrogen-responsive ERE-CAT reporter gene, cell growth in terms of proliferation rate after 7 days or cell growth in terms of saturation density after 14 days. The IC50 values for these three assays in order were for 17 beta-oestradiol (1 x 10-(11) M, 1 x 10-(11) M, 2 x 10(-11) M), and in rank order of potency for the phytoestrogens, deoxymiroestrol (1 x 10(-10) M, 3 x 10(-11) M, 2 x 10(-11) M) > miroestrol (3 x 10(-10) M, 2 x 10(-11) M, 8 x 10(-11) M) > 8-prenylnaringenin (1 x 10(-9) M, 3 x 10(-10) M, 3 x 10(-10) M) > cournestrol (3 x 10(-8) M, 2 x 10(-8) M, 3 x 10(-8) M) > genistein (4 x 10(-8) M, 2 x 10(-8) M, 1 x 10(-8) M)/equol (1 x 10(-7) M, 3 x 10(-8) M, 2 x 10(-8) M) > daidzein (3 x 10(-7) M, 2 x 10(-7) M, 4 x 10(-8) M) > resveratrol (4 x 10(-6) M, not achieved, not achieved). Despite using the same receptor context of the MCF7 cells, this rank order differed from that determined from receptor binding. The most marked difference was for cournestrol and 8-prenylnaringenin which both displayed a relatively potent ability to displace [3H]oestradiol from cytosolic ER compared with their much lower activity in the cell-based assays. Albeit at varying concentrations, seven of the eight phytoestrogens (all except resveratrol) gave similar maximal responses to that given by 17 beta-oestradiol in cell-based assays which makes them full oestrogen agonists. We found no evidence for any oestrogen antagonist action of any of these phytoestrogens at concentrations of up to 10(-6) M on either reporter gene induction or on stimulation of cell growth. (c) 2005 Elsevier Ltd. All rights reserved.

Bayesian estimation of ancestral character states on phylogenies

Biologists frequently attempt to infer the character states at ancestral nodes of a phylogeny from the distribution of traits observed in contemporary organisms. Because phylogenies are normally inferences from data, it is desirable to account for the uncertainty in estimates of the tree and its branch lengths when making inferences about ancestral states or other comparative parameters. Here we present a general Bayesian approach for testing comparative hypotheses across statistically justified samples of phylogenies, focusing on the specific issue of reconstructing ancestral states. The method uses Markov chain Monte Carlo techniques for sampling phylogenetic trees and for investigating the parameters of a statistical model of trait evolution. We describe how to combine information about the uncertainty of the phylogeny with uncertainty in the estimate of the ancestral state. Our approach does not constrain the sample of trees only to those that contain the ancestral node or nodes of interest, and we show how to reconstruct ancestral states of uncertain nodes using a most-recent-common-ancestor approach. We illustrate the methods with data on ribonuclease evolution in the Artiodactyla. Software implementing the methods ( BayesMultiState) is available from the authors.

A novel spectral representation of electromyographic signals

Time/frequency and temporal analyses have been widely used in biomedical signal processing. These methods represent important characteristics of a signal in both time and frequency domain. In this way, essential features of the signal can be viewed and analysed in order to understand or model the physiological system. Historically, Fourier spectral analyses have provided a general method for examining the global energy/frequency distributions. However, an assumption inherent to these methods is the stationarity of the signal. As a result, Fourier methods are not generally an appropriate approach in the investigation of signals with transient components. This work presents the application of a new signal processing technique, empirical mode decomposition and the Hilbert spectrum, in the analysis of electromyographic signals. The results show that this method may provide not only an increase in the spectral resolution but also an insight into the underlying process of the muscle contraction.

Alignment of a model amyloid peptide fragment in bulk and at a solid surface

The alignment of model amyloid peptide YYKLVFFC is investigated in bulk and at a solid surface using a range of spectroscopic methods employing polarized radiation. The peptide is based on a core sequence of the amyloid beta (A beta) peptide, KLVFF. The attached tyrosine and cysteine units are exploited to yield information on alignment and possible formation of disulfide or dityrosine links. Polarized Raman spectroscopy on aligned stalks provides information on tyrosine orientation, which complements data from linear dichroism (LD) on aqueous solutions subjected to shear in a Couette cell. LD provides a detailed picture of alignment of peptide strands and aromatic residues and was also used to probe the kinetics of self-assembly. This suggests initial association of phenylalanine residues, followed by subsequent registry of strands and orientation of tyrosine residues. X-ray diffraction (XRD) data from aligned stalks is used to extract orientational order parameters from the 0.48 nm reflection in the cross-beta pattern, from which an orientational distribution function is obtained. X-ray diffraction on solutions subject to capillary flow confirmed orientation in situ at the level of the cross-beta pattern. The information on fibril and tyrosine orientation from polarized Raman spectroscopy is compared with results from NEXAFS experiments on samples prepared as films on silicon. This indicates fibrils are aligned parallel to the surface, with phenyl ring normals perpendicular to the surface. Possible disulfide bridging leading to peptide dimer formation was excluded by Raman spectroscopy, whereas dityrosine formation was probed by fluorescence experiments and was found not to occur except under alkaline conditions. Congo red binding was found not to influence the cross-beta XRD pattern.

European Union: Shadow WTO agricultural domestic support notifications

The notification of the level of domestic support to the World Trade Organization (WTO) is intended to reflect compliance with obligations entered into at the time of the Uruguay Round. WTO members have often been slow to provide notification of domestic support levels. This makes the process of notification less useful as an indicator of the degree to which changes in policy have or have not benefited the trade system as a whole and exporting countries in particular. The notification of domestic support in the E.U. illustrates the value of a measure that reflects current policies and can therefore act as a basis for negotiation of further disciplines where these are necessary. The E.U. has made major changes in its Common Agricultural Policy (CAP) over the period since 1992 when the MacSharry reforms were implemented. Payments originally notified in the blue box (related to supply control) have over time been changed until in their present form they are unrelated to current production or price levels, and hence can satisfy the criteria for the green box. The E.U. has therefore much more latitude in trade talks to agree to reductions in the allowable trade-distorting support. This paper reproduced the E.U. notifications relating to 2003/04 and extends these with official statistics to the year 2006/07. It then projects forward the components of domestic support until the year 2013/14, based on forecasts of future production and estimates of policy parameters. The impact of a successful Doha Round is simulated, showing that the constraints envisaged in the WTO draft modalities document of May 19, 2008, would be binding by the year 2013, at about the time the next budget cycle in the E.U. starts. Without the Doha Round constraints, further reform might still happen for domestic reasons, but the framework provided by the WTO for domestic policy spending would be less relevant. In that case, much could hinge on the legitimacy of the Single Farm Payment system under the current rules governing the green box.

Perspex machine X: software development

The Perspex Machine arose from the unification of computation with geometry. We now report significant redevelopment of both a partial C compiler that generates perspex programs and of a Graphical User Interface (GUI). The compiler is constructed with standard compiler-generator tools and produces both an explicit parse tree for C and an Abstract Syntax Tree (AST) that is better suited to code generation. The GUI uses a hash table and a simpler software architecture to achieve an order of magnitude speed up in processing and, consequently, an order of magnitude increase in the number of perspexes that can be manipulated in real time (now 6,000). Two perspex-machine simulators are provided, one using trans-floating-point arithmetic and the other using transrational arithmetic. All of the software described here is available on the world wide web. The compiler generates code in the neural model of the perspex. At each branch point it uses a jumper to return control to the main fibre. This has the effect of pruning out an exponentially increasing number of branching fibres, thereby greatly increasing the efficiency of perspex programs as measured by the number of neurons required to implement an algorithm. The jumpers are placed at unit distance from the main fibre and form a geometrical structure analogous to a myelin sheath in a biological neuron. Both the perspex jumper-sheath and the biological myelin-sheath share the computational function of preventing cross-over of signals to neurons that lie close to an axon. This is an example of convergence driven by similar geometrical and computational constraints in perspex and biological neurons.

The distribution of solar wind speeds during solar minimum: calibration for numerical solar wind modeling constraints on the source of the slow solar wind

It took the solar polar passage of Ulysses in the early 1990s to establish the global structure of the solar wind speed during solar minimum. However, it remains unclear if the solar wind is composed of two distinct populations of solar wind from different sources (e.g., closed loops which open up to produce the slow solar wind) or if the fast and slow solar wind rely on the superradial expansion of the magnetic field to account for the observed solar wind speed variation. We investigate the solar wind in the inner corona using the Wang-Sheeley-Arge (WSA) coronal model incorporating a new empirical magnetic topology–velocity relationship calibrated for use at 0.1 AU. In this study the empirical solar wind speed relationship was determined by using Helios perihelion observations, along with results from Riley et al. (2003) and Schwadron et al. (2005) as constraints. The new relationship was tested by using it to drive the ENLIL 3-D MHD solar wind model and obtain solar wind parameters at Earth (1.0 AU) and Ulysses (1.4 AU). The improvements in speed, its variability, and the occurrence of high-speed enhancements provide confidence that the new velocity relationship better determines the solar wind speed in the outer corona (0.1 AU). An analysis of this improved velocity field within the WSA model suggests the existence of two distinct mechanisms of the solar wind generation, one for fast and one for slow solar wind, implying that a combination of present theories may be necessary to explain solar wind observations.

Elastic-net prefiltering for two-class classification

A two-stage linear-in-the-parameter model construction algorithm is proposed aimed at noisy two-class classification problems. The purpose of the first stage is to produce a prefiltered signal that is used as the desired output for the second stage which constructs a sparse linear-in-the-parameter classifier. The prefiltering stage is a two-level process aimed at maximizing a model's generalization capability, in which a new elastic-net model identification algorithm using singular value decomposition is employed at the lower level, and then, two regularization parameters are optimized using a particle-swarm-optimization algorithm at the upper level by minimizing the leave-one-out (LOO) misclassification rate. It is shown that the LOO misclassification rate based on the resultant prefiltered signal can be analytically computed without splitting the data set, and the associated computational cost is minimal due to orthogonality. The second stage of sparse classifier construction is based on orthogonal forward regression with the D-optimality algorithm. Extensive simulations of this approach for noisy data sets illustrate the competitiveness of this approach to classification of noisy data problems.

Analysis of full-waveform LiDAR data for classification of an orange orchard scene

Full-waveform laser scanning data acquired with a Riegl LMS-Q560 instrument were used to classify an orange orchard into orange trees, grass and ground using waveform parameters alone. Gaussian decomposition was performed on this data capture from the National Airborne Field Experiment in November 2006 using a custom peak-detection procedure and a trust-region-reflective algorithm for fitting Gauss functions. Calibration was carried out using waveforms returned from a road surface, and the backscattering coefficient c was derived for every waveform peak. The processed data were then analysed according to the number of returns detected within each waveform and classified into three classes based on pulse width and c. For single-peak waveforms the scatterplot of c versus pulse width was used to distinguish between ground, grass and orange trees. In the case of multiple returns, the relationship between first (or first plus middle) and last return c values was used to separate ground from other targets. Refinement of this classification, and further sub-classification into grass and orange trees was performed using the c versus pulse width scatterplots of last returns. In all cases the separation was carried out using a decision tree with empirical relationships between the waveform parameters. Ground points were successfully separated from orange tree points. The most difficult class to separate and verify was grass, but those points in general corresponded well with the grass areas identified in the aerial photography. The overall accuracy reached 91%, using photography and relative elevation as ground truth. The overall accuracy for two classes, orange tree and combined class of grass and ground, yielded 95%. Finally, the backscattering coefficient c of single-peak waveforms was also used to derive reflectance values of the three classes. The reflectance of the orange tree class (0.31) and ground class (0.60) are consistent with published values at the wavelength of the Riegl scanner (1550 nm). The grass class reflectance (0.46) falls in between the other two classes as might be expected, as this class has a mixture of the contributions of both vegetation and ground reflectance properties.

Some results from an intercomparison of the climates simulated by 14 atmospheric general circulation models

Some climatological information from 14 atmospheric general circulation models is presented and compared in order to assess the ability of a broad group of models to simulate current climate. The quantities considered are cross sections of temperature, zonal wind, and meridional stream function together with latitudinal distributions of mean sea level pressure and precipitation rate. The nature of the deficiencies in the simulated climates that are common to all models and those which differ among models is investigated; the general improvement in the ability of models to simulate certain aspects of the climate is shown; consideration is given to the effect of increasing resolution on simulated climate; and approaches to understanding and reducing model deficiencies are discussed. The information presented here is a subset of a more voluminous compilation which is available in report form (Boer et al., 1991). This report contains essentially the same text, but results from all 14 models are presented together with additional results in the form of geographical distributions of surface variables and certain difference statistics.

An intercomparison of the climates simulated by 14 atmospheric general circulation models

Climatological information from fourteen atmospheric general circulation models is presented and compared in order to assess the ability of a broad group of models to simulate current climate. The quantities considered are cross sections of temperature, zonal wind and meridional stream function together with latitudinal distributions of mean sea-level pressure and precipitation rate. The nature of the deficiencies in the simulated climates that are common to all models and those which differ among models is investigated, general improvement in the ability of models to simulate certain aspects of the climate is shown, consideration is given to the effect of increasing resolution on simulated climate and approaches to the understanding and reduction of model deficiencies are discussed.

Climate impacts on river flow: projections for the Medway catchment UK with UKCP09 and CATCHMOD

The potential impact of climate change on areas of strategic importance for water resources remains a concern. Here, river flow projections for the River Medway, above Teston in southeast England are presented, which is just such an area of strategic importance. The river flow projections use climate inputs from the Hadley Centre Regional Climate Model (HadRM3) for the time period 1960–2080 (a subset of the early release UKCP09 projections). River flow predictions are calculated using CATCHMOD, the main river flow prediction tool of the Environment Agency (EA) of England and Wales. In order to use this tool in the best way for climate change predictions, model setup and performance are analysed using sensitivity and uncertainty analysis. The model's representation of hydrological processes is discussed and the direct percolation and first linear storage constant parameters are found to strongly affect model results in a complex way, with the former more important for low flows and the latter for high flows. The uncertainty in predictions resulting from the hydrological model parameters is demonstrated and the projections of river flow under future climate are analysed. A clear climate change impact signal is evident in the results with a persistent lowering of mean daily river flows for all months and for all projection time slices. Results indicate that a projection of lower flows under future climate is valid even taking into account the uncertainties considered in this modelling chain exercise. The model parameter uncertainty becomes more significant under future climate as the river flows become lower. This has significant implications for those making policy decisions based on such modelling results. Copyright © 2010 John Wiley & Sons, Ltd.

The effect of model configuration on modelled hillslope–riparian interactions

The transfer of hillslope water to and through the riparian zone forms a research area of importance in hydrological investigations. Numerical modelling schemes offer a way to visualise and quantify first-order controls on catchment runoff response and mixing. We use a two-dimensional Finite Element model to assess the link between model setup decisions (e.g. zero-flux boundary definitions, soil algorithm choice) and the consequential hydrological process behaviour. A detailed understanding of the consequences of model configuration is required in order to produce reliable estimates of state variables. We demonstrate that model configuration decisions can determine effectively the presence or absence of particular hillslope flow processes and, the magnitude and direction of flux at the hillslope–riparian interface. If these consequences are not fully explored for any given scheme and application, the resulting process inference may well be misleading.

Open solar flux estimates from near-Earth measurements of the interplanetary magnetic field: comparison of the first two perihelion passes of the Ulysses spacecraft

Results from all phases of the orbits of the Ulysses spacecraft have shown that the magnitude of the radial component of the heliospheric field is approximately independent of heliographic latitude. This result allows the use of near- Earth observations to compute the total open flux of the Sun. For example, using satellite observations of the interplanetary magnetic field, the average open solar flux was shown to have risen by 29% between 1963 and 1987 and using the aa geomagnetic index it was found to have doubled during the 20th century. It is therefore important to assess fully the accuracy of the result and to check that it applies to all phases of the solar cycle. The first perihelion pass of the Ulysses spacecraft was close to sunspot minimum, and recent data from the second perihelion pass show that the result also holds at solar maximum. The high level of correlation between the open flux derived from the various methods strongly supports the Ulysses discovery that the radial field component is independent of latitude. We show here that the errors introduced into open solar flux estimates by assuming that the heliospheric field’s radial component is independent of latitude are similar for the two passes and are of order 25% for daily values, falling to 5% for averaging timescales of 27 days or greater. We compare here the results of four methods for estimating the open solar flux with results from the first and second perehelion passes by Ulysses. We find that the errors are lowest (1–5% for averages over the entire perehelion passes lasting near 320 days), for near-Earth methods, based on either interplanetary magnetic field observations or the aa geomagnetic activity index. The corresponding errors for the Solanki et al. (2000) model are of the order of 9–15% and for the PFSS method, based on solar magnetograms, are of the order of 13–47%. The model of Solanki et al. is based on the continuity equation of open flux, and uses the sunspot number to quantify the rate of open flux emergence. It predicts that the average open solar flux has been decreasing since 1987, as Correspondence to: M. Lockwood (m.lockwood@rl.ac.uk) is observed in the variation of all the estimates of the open flux. This decline combines with the solar cycle variation to produce an open flux during the second (sunspot maximum) perihelion pass of Ulysses which is only slightly larger than that during the first (sunspot minimum) perihelion pass.