Estimating and evaluating the statistics of gapped local-alignment scores

We present a novel maximum-likelihood-based algorithm for estimating the distribution of alignment scores from the scores of unrelated sequences in a database search. Using a new method for measuring the accuracy of p-values, we show that our maximum-likelihood-based algorithm is more accurate than existing regression-based and lookup table methods. We explore a more sophisticated way of modeling and estimating the score distributions (using a two-component mixture model and expectation maximization), but conclude that this does not improve significantly over simply ignoring scores with small E-values during estimation. Finally, we measure the classification accuracy of p-values estimated in different ways and observe that inaccurate p-values can, somewhat paradoxically, lead to higher classification accuracy. We explain this paradox and argue that statistical accuracy, not classification accuracy, should be the primary criterion in comparisons of similarity search methods that return p-values that adjust for target sequence length.

Chemical and structural disorder in eumelanins: A possible explanation for broadband absorbance

We report the results of an experimental and theoretical study of the electronic and structural properties of a key eumelanin precursor-5,6,-dihydroxyindole-2-carboxylic acid ( DHICA) - and its dimeric forms. We have used optical spectroscopy to follow the oxidative polymerization of DHICA to eumelanin and observe red shifting and broadening of the absorption spectrum as the reaction proceeds. First principles density functional theory calculations indicate that DHICA oligomers ( possible reaction products of oxidative polymerization) have the gap between the highest occupied molecular orbital and the lowest unoccupied molecular orbital red-shifted gaps with respect to the monomer. Furthermore, different bonding configurations ( leading to oligomers with different structures) produce a range of gaps. These experimental and theoretical results lend support to the chemical disorder model where the broadband monotonic absorption characteristic of all melanins is a consequence of the superposition of a large number of nonhomogeneously broadened Gaussian transitions associated with each of the components of a melanin ensemble. These results suggest that the traditional model of eumelanin as an amorphous organic semiconductor is not required to explain its optical properties and should be thoroughly reexamined. These results have significant implications for our understanding of the physics, chemistry, and biological function of these important biological macromolecules. Indeed, one may speculate that the robust functionality of melanins in vitro is a direct consequence of its heterogeneity, i.e., chemical disorder is a "low cost" natural resource in these systems

Simulations of thermal Bose fields in the classical limit

We demonstrate that the time-dependent projected Gross-Pitaevskii equation (GPE) derived earlier [M. J. Davis, R. J. Ballagh, and K. Burnett, J. Phys. B 34, 4487 (2001)] can represent the highly occupied modes of a homogeneous, partially-condensed Bose gas. Contrary to the often held belief that the GPE is valid only at zero temperature, we find that this equation will evolve randomized initial wave functions to a state describing thermal equilibrium. In the case of small interaction strengths or low temperatures, our numerical results can be compared to the predictions of Bogoliubov theory and its perturbative extensions. This demonstrates the validity of the GPE in these limits and allows us to assign a temperature to the simulations unambiguously. However, the GPE method is nonperturbative, and we believe it can be used to describe the thermal properties of a Bose gas even when Bogoliubov theory fails. We suggest a different technique to measure the temperature of our simulations in these circumstances. Using this approach we determine the dependence of the condensate fraction and specific heat on temperature for several interaction strengths, and observe the appearance of vortex networks. Interesting behavior near the critical point is observed and discussed.

Characterization of a distinct plasma membrane macrodomain in differentiated adipocytes

Caveolae are small invaginations of the cell surface that are abundant in mature adipocytes. A recent study (Kanzaki, M., and Pessin, J. E. (2002) J. Biol Chem 277, 25867-25869) described novel caveolin- and actin-containing structures associated with the adipocyte cell surface that contain specific signaling proteins. We have characterized these structures, here termed caves, using light and electron microscopy and observe that they represent surface-connected wide invaginations of the basal plasma membrane that are sometimes many micrometers in diameter. Rather than simply a caveolar domain, these structures contain all elements of the plasma membrane including clathrin-coated pits, lipid raft markers, and non-raft markers. GLUT4 is recruited to caves in response to insulin stimulation. Caves can occupy a significant proportion of the plasma membrane area and are surrounded by cortical actin. Caveolae density in caves is similar to that on the bulk plasma membrane, but because these structures protrude much deeper into the plane of focus of the light microscope molecules such as caveolin and other plasma membrane proteins appear more concentrated in caves. We conclude that the adipocyte surface membrane contains numerous wide invaginations that do not represent novel caveolar structures but rather large surface caves.

The structure of quantum Lie algebras for the classical series, B-l, C-l and D-l

The structure constants of quantum Lie algebras depend on a quantum deformation parameter q and they reduce to the classical structure constants of a Lie algebra at q = 1. We explain the relationship between the structure constants of quantum Lie algebras and quantum Clebsch-Gordan coefficients for adjoint x adjoint --> adjoint We present a practical method for the determination of these quantum Clebsch-Gordan coefficients and are thus able to give explicit expressions for the structure constants of the quantum Lie algebras associated to the classical Lie algebras B-l, C-l and D-l. In the quantum case the structure constants of the Cartan subalgebra are non-zero and we observe that they are determined in terms of the simple quantum roots. We introduce an invariant Killing form on the quantum Lie algebras and find that it takes values which are simple q-deformations of the classical ones.

Millennial and orbital variations of El Nino/Southern Oscillation and high-latitude climate in the last glacial period

The El Nino/Southern Oscillation (ENSO) phenomenon is believed to have operated continuously over the last glacial interglacial cycle(1). ENSO variability has been suggested to be linked to millennial-scale oscillations in North Atlantic climate during that time(2,3), but the proposals disagree on whether increased frequency of El Nino events, the warm phase of ENSO, was linked to North Atlantic warm or cold periods. Here we present a high-resolution record of surface moisture, based on the degree of peat humification and the ratio of sedges to grass, from northern Queensland, Australia, covering the past 45,000 yr. We observe millennial-scale dry periods, indicating periods of frequent El Nino events ( summer precipitation declines in El Nino years in northeastern Australia). We find that these dry periods are correlated to the Dansgaard - Oeschger events - millennial-scale warm events in the North Atlantic climate record - although no direct atmospheric connection from the North Atlantic to our site can be invoked. Additionally, we find climatic cycles at a semiprecessional timescale (, 11,900 yr). We suggest that climate variations in the tropical Pacific Ocean on millennial as well as orbital timescales, which determined precipitation in northeastern Australia, also exerted an influence on North Atlantic climate through atmospheric and oceanic teleconnections.

Nutrient-induced perturbations to delta C-13 and delta N-15 in symbiotic dinoilagellates and their coral hosts

Inorganic nutrients play a critical role in determining benthic community structure in tropical seas. This study examined the impact of adding inorganic nutrients (ammonium and phosphate) on the isotopic composition of 2 reef-building corals, Pocillopora damicornis and Heliofungia actiniformis, on the southern Great Barrier Reef. The addition of elevated nutrients to patch reefs that pond at low tide did not perturb the C:N ratio of either species or their symbiotic dinoflagellates. The C:N ratios were significantly higher in material extracted from the skeleton (14.8 +/- 1.50 and 10.8 +/- 1.42) than either host (7.6 +/- 0.87 and 6.0 +/- 0.71) or symbiotic dinoflagellates (5.7 +/- 0.48 and 6.9 +/- 0.66) (P. damicornis and H. actiniformis respectively; 95 confidence intervals). The ratio of acquired N to background N suggests that the added dissolved inorganic nitrogen (DIN) accounted for 50 to 100% of total nitrogen within the tissues of P. damicornis and H. actiniformis at the end of the experiment. The addition of the isotopically depleted nutrients (delta(15) N = 0parts per thousand) to patch reefs significantly decreased delta(15)N from control values of between 3 and 4 to values to below 1 in the case of all compartments, while delta(13)C values were relatively unresponsive to nutrient treatments. These findings suggest that coral delta(15)N has the potential to provide a historical record of the delta(15)N of dissolved nitrogen surrounding reef-building corals and their symbiotic dinoflagellates.

The selection of intended actions and the observation of others' actions: A time-resolved fMRI study

Whenever we plan, imagine, or observe an action, the motor systems that would be involved in preparing and executing that action are similarly engaged. The way in which such common motor activation is formed, however, is likely to differ depending on whether it arises from our own intentional selection of action or from the observation of another's action. In this study, we use time-resolved event-related functional MRI to tease apart neural processes specifically related to the processing of observed actions, the selection of our own intended actions, the preparation for movement, and motor response execution. Participants observed a finger gesture movement or a cue indicating they should select their own finger gesture to perform, followed by a 5-s delay period; participants then performed the observed or self-selected action. During the preparation and readiness for action, prior to initiation, we found activation in a common network of higher motor areas, including dorsal and ventral premotor areas and the pre-supplementary motor area (pre-SMA); the more caudal SMA showed greater activation during movement execution. Importantly, the route to this common motor activation differed depending on whether participants freely selected the actions to perform or whether they observed the actions performed by another person. Observation of action specifically involved activation of inferior and superior parietal regions, reflecting involvement of the dorsal visual pathway in visuomotor processing required for planning the action. In contrast, the selection of action specifically involved the dorsal lateral prefrontal and anterior cingulate cortex, reflecting the role of these prefrontal areas in attentional selection and guiding the selection of responses. (c) 2005 Elsevier Inc. All rights reserved.

Marine macroalgal biodiversity hotspots: why is there high species richness and endemism in southern Australian marine benthic flora?

The southern Australian marine macroalgal flora has the highest levels of species richness and endemism of any regional macroalgal flora in the world. Analyses of species composition and distributions for the southern Australian flora have identified four different floristic elements, namely the southern Australian endemic element, the widely distributed temperate element, the tropical element and a cold water element. Within the southern Australian endemic element, four species distribution patterns are apparent, thought to largely result from the Jurassic to Oligocene fragmentation of East Gondwana, the subsequent migration of Tethyan ancestors from the west Australian coast and the later invasion of high latitude Pacific species. Climatic deterioration from the late Eocene to the present is thought responsible for the replacement of the previous tropical south coast flora by an endemic temperate flora which has subsequently diversified in response to fluctuating environmental conditions, abundant rocky substrata and substantial habitat heterogeneity. High levels of endemism are attributed to Australia's long isolation and maintained, as is the high species richness, by the lack of recent mass extinction events. The warm water Leeuwin Current has had profound influence in the region since the Eocene, flowing to disperse macroalgal species onto the south coast as well as ameliorating the local environment. It is now evident that the high species richness and endemism we now observe in the southern Australian marine macroalgal flora can be attributed to a complex interaction of biogeographical, ecological and phylogenetic processes over the last 160 million years.

Increase of fat oxidation and weight loss in obese mice caused by chronic treatment with human growth hormone or a modified C-terminal fragment

OBJECTIVE: To observe the chronic effects of human growth hormone (hGH) and AOD9604 (a C-terminal fragment of hGH) on body weight, energy balance, and substrate oxidation rates in obese (ob/ob) and lean C57BL/6Jmice. In vitro assays were used to confirm whether the effects of AOD9604 are mediated through the hGH receptor, and if this peptide is capable of cell proliferation via the hGH receptor. METHOD: Obese and lean mice were treated with hGH, AOD or saline for 14 days using mini-osmotic pumps. Body weight, caloric intake, resting energy expenditure, fat oxidation, glucose oxidation, and plasma glucose, insulin and glycerol were measured before and after treatment. BaF-BO3 cells transfected with the hGH receptor were used to measure in Vitro I-125-hGH receptor binding and cell proliferation. RESULTS: Both hGH and AOD significantly reduced body weight gain in obese mice. This was associated with increased in vivo fat oxidation and increased plasma glycerol levels (an index of lipolysis). Unlike hGH, however, AOD9604 did not induce hyperglycaemia or reduce insulin secretion. AOD9604 does not compete for the hGH receptor and nor does it induce cell proliferation, unlike hGH. CONCLUSIONS: Both hGH and its C-terminal fragment reduce body weight gain, increase fat oxidation, and stimulate lipolysis in obese mice, yet AOD9604 does not interact with the hGH receptor. Thus, the concept of hGH behaving as a pro-hormone is further confirmed. This data shows that fragments of hGH can act in a manner novel to traditional hGH-stimulated pathways.

Being breastfed in infancy and breast cancer incidence in adult life: Results from the two Nurses' Health studies

Events during perinatal and early life may influence the incidence of breast cancer in adult life, and some case-control studies suggest that having been breastfed may reduce breast cancer risk. The authors studied this association among premenopausal and postmenopausal women by using data from the two Nurses' Health Studies, the Nurses' Health Study (using data from 1992 to 1996) and the Nurses' Health Study II (using data from 1991 to 1997). A history of being breastfed was self-reported by the study participants. During a total of 695,655 person-years, 1,073 cases of invasive breast cancer were diagnosed. The authors did not observe any important overall association between having been breastfed and the development of breast cancer later in life among premenopausal women (covariate-adjusted relative risk = 0.97, 95% confidence interval (CI): 0.78, 1.20) or postmenopausal women (covariate-adjusted relative risk = 1.12, 95% CI: 0.92, 1.37). No significant trend was observed with increasing duration of breastfeeding. The authors also used data on breastfeeding retrospectively collected from 2,103 mothers of participants of the two Nurses' Health Studies. With the mothers' reports, the covariate-adjusted odds ratio of breast cancer was 1.11 (95% CI: 0.88, 1.39) for women who were breastfed compared with those who were not. Data from these two large cohorts do not support the hypothesis that being breastfed confers protection against subsequent breast cancer.

Causal and localizable quantum operations

We examine constraints on quantum operations imposed by relativistic causality. A bipartite superoperator is said to be localizable if it can be implemented by two parties (Alice and Bob) who share entanglement but do not communicate, it is causal if the superoperator does not convey information from Alice to Bob or from Bob to Alice. We characterize the general structure of causal complete-measurement superoperators, and exhibit examples that are causal but not localizable. We construct another class of causal bipartite superoperators that are not localizable by invoking bounds on the strength of correlations among the parts of a quantum system. A bipartite superoperator is said to be semilocalizable if it can be implemented with one-way quantum communication from Alice to Bob, and it is semicausal if it conveys no information from Bob to Alice. We show that all semicausal complete-measurement superoperators are semi localizable, and we establish a general criterion for semicausality. In the multipartite case, we observe that a measurement superoperator that projects onto the eigenspaces of a stabilizer code is localizable.

Mating aggregations and mating success in the flower thrips, Frankliniella schultzei (Thysanoptera: Thripidae), and a possible role for pheromones

Aggregations of Frankliniella schultzei males were observed on the corollas of Hibiscus rosasinensis and Gossypium hirsutum flowers in southeast Queensland. Aggregations were seen only on the upper surfaces of corollas but may have occurred on other flower parts, which were hidden from view. Conspecific females entered aggregations and a small proportion of them mated [18% (n = 163), H. rosasinensis; 30% (n = 181), G. hirsutum]. Most females (87 and 72%, respectively) that did not mate in aggregations walked to other flower parts. Behavior was difficult to observe on these parts, but mating was sometimes observed there. The number of females that landed within aggregations on the upper surfaces of both H. rosasinensis and G. hirsutum corollas was highly correlated with the number of males (r = 0.88, r = 0.93, respectively; P < 0.001). Significantly more mating pairs were observed in high-density aggregations (mean +/- SE, 1.10 +/- 0.22 and 4.44 +/- 0.48, respectively) than in low-density aggregations (0.37 +/- 0.11 and 1.67 +/- 0.29, respectively) (P < 0.05) on flowers of both species. More F. schultzei females were attracted to sticky traps baited with live conspecific males set among flowering Ipomoea indica (mean +/- SE, 8.83 +/- 0.32) and G. hirsutum (10.90 +/- 0.79) plants than to control traps (0.10 +/- 0.05 and 0.70 +/- 0.25, respectively)( P < 0.05), presumably in response to male-produced pheromones. Significantly more females were attracted to traps with high male densities than to traps with low densities. We found no statistical evidence that aggregation size influenced mating success (proportion males that mated). Mating success, however, should be evaluated with respect to mating on all flower parts and not just the upper surfaces of corollas. The results of this study constitute the first behavioral evidence for an attractant sex pheromone in thrips.

A curious experiment: the paradigm switch from observation and speculation to experimentation, in the understanding of neuromuscular function and disease

The four-link chain of the motor unit represents the contemporary end-point of some two millennia of evolving knowledge in neuroscience. The paradigm shift in neuromuscular epistemology occurred in the mid-17th century. In 1666, the newly graduated Dutch doctor, Jan Swammerdam (1637-1680) published his former investigations of dissected nerve-muscle preparations. These experiments comprised the quantum leap from observation and speculation, to that of experimentation in the field of neuroanatomy and neurophysiology. In what he termed 'A Curious Experiment' he also described the phenomenon of intrinsic muscle excitability - I cannot observe that the muscle in the living animal ever absolutely ceases from all motion. Eighty years later (1752), von Haller demonstrated experimentally that irritability (contractility) was an intrinsic property of all muscular tissue; and distinguished between the sensibility of nerve impulses and the irritability of muscular contraction. This experimental progression from Swammerdam to von Haller culminated in 1850, when Claude Bernard's studies in experimental pharmacology confirmed that muscle was a functional unit, independent of any electrical innervation via its supplying nerve. This account comprises an audit of Swammerdam's work in the perspective of neuromuscular knowledge. (C) 2002 Elsevier Science B.V. All rights reserved.

Numerical simulations of wind and temperature structure within an Alpine lake basin, Lake Tekapo, New Zealand

High-resolution numerical model simulations have been used to study the local and mesoscale thermal circulations in an Alpine lake basin. The lake (87 km(2)) is situated in the Southern Alps, New Zealand and is located in a glacially excavated rock basin surrounded by mountain ranges that reach 3000 m in height. The mesoscale model used (RAMS) is a three-dimensional non-hydrostatic model with a level 2.5 turbulence closure scheme. The model demonstrates that thermal forcing at local (within the basin) and regional (coast-to-basin inflow) scales drive the observed boundary-layer airflow in the lake basin during clear anticyclonic summertime conditions. The results show that the lake can modify (perturb) both the local and regional wind systems. Following sunrise, local thermal circulations dominate, including a lake breeze component that becomes embedded within the background valley wind system. This results in a more divergent flow in the basin extending across the lake shoreline. However, a closed lake breeze circulation is neither observed nor modelled. Modelling results indicate that in the latter part of the day when the mesoscale (coast-to-basin) inflow occurs, the relatively cold pool of lake air in the basin can cause the intrusion to decouple from the surface. Measured data provide qualitative and quantitative support for the model results.