BEAGLE: an application programming interface and high-performance computing library for statistical phylogenetics

Phylogenetic inference is fundamental to our understanding of most aspects of the origin and evolution of life, and in recent years, there has been a concentration of interest in statistical approaches such as Bayesian inference and maximum likelihood estimation. Yet, for large data sets and realistic or interesting models of evolution, these approaches remain computationally demanding. High-throughput sequencing can yield data for thousands of taxa, but scaling to such problems using serial computing often necessitates the use of nonstatistical or approximate approaches. The recent emergence of graphics processing units (GPUs) provides an opportunity to leverage their excellent floating-point computational performance to accelerate statistical phylogenetic inference. A specialized library for phylogenetic calculation would allow existing software packages to make more effective use of available computer hardware, including GPUs. Adoption of a common library would also make it easier for other emerging computing architectures, such as field programmable gate arrays, to be used in the future. We present BEAGLE, an application programming interface (API) and library for high-performance statistical phylogenetic inference. The API provides a uniform interface for performing phylogenetic likelihood calculations on a variety of compute hardware platforms. The library includes a set of efficient implementations and can currently exploit hardware including GPUs using NVIDIA CUDA, central processing units (CPUs) with Streaming SIMD Extensions and related processor supplementary instruction sets, and multicore CPUs via OpenMP. To demonstrate the advantages of a common API, we have incorporated the library into several popular phylogenetic software packages. The BEAGLE library is free open source software licensed under the Lesser GPL and available from http://beagle-lib.googlecode.com. An example client program is available as public domain software.     

Excitation energy transfer and trapping dynamics in the core complex of the filamentous photosynthetic bacterium Roseiflexus castenholzii

The light-harvesting core complex of the thermophilic filamentous anoxygenic phototrophic bacterium Roseiflexus castenholzii is intrinsic to the cytoplasmic membrane and intimately bound to the reaction center (RC). Using ultrafast transient absorption and time-resolved fluorescence spectroscopy with selective excitation, energy transfer, and trapping dynamics in the core complex have been investigated at room temperature in both open and closed RCs. Results presented in this report revealed that the excited energy transfer from the BChl 800 to the BChl 880 band of the antenna takes about 2?ps independent of the trapping by the RC. The time constants for excitation quenching in the core antenna BChl 880 by open and closed RCs were found to be 60 and 210?ps, respectively. Assuming that the light harvesting complex is generally similar to LH1 of purple bacteria, the possible structural and functional aspects of this unique antenna complex are discussed. The results show that the core complex of Roseiflexus castenholzii contains characteristics of both purple bacteria and Chloroflexus aurantiacus.     

Stability of helix-rich proteins at high concentrations

A number of techniques, including circular dichroism, FTIR, front face fluorescence, and UV absorption spectrophotometries, dynamic light scattering, and DSC, were used to directly measure the colloidal and conformational stability of proteins in highly concentrated solutions. Using bovine serum albumin (BSA), chicken egg white lysozyme, human hemoglobin A0, and bovine fibrinogen as model proteins, the thermal transition temperatures of proteins in dilute and concentrated solutions were compared. At 10 degrees C, no significant differences in both secondary and tertiary structures were detected for proteins at different concentrations. When temperature was introduced as a variable, however, hemoglobin and fibrinogen demonstrated higher transition midpoints (T(m)s) in concentrated rather than in dilute solutions (deltaT(m) approximately 2-10 degrees C). In contrast, lysozyme and BSA in concentrated solutions exhibit a lower T(m) than in dilute solutions (deltaT(m) approximately 2-20 degrees C). From these studies, it appears that a variety of factors determine the effect of high concentrations on the colloidal and conformational stability of a particular protein. While the prediction of excluded volume theory is that high concentrations should conformationally stabilize proteins, other factors such as pH, kinetics, protein dynamics, and intermolecular charge-charge effects may affect the overall stability of proteins at high concentrations under certain conditions.     

Dispersal traits linked to range size through range location, not dispersal ability, in Western Australian angiosperms

Aim  We examine the relative importance of seed dispersal mode in determining the range size and range placement in 524 species from six focal plant families (Agavaceae, Euphorbiaceae, Malvacaeae, Sapindaceae, Proteaceae and Fabaceae ( Acacia )).
Location  Western Australia.
Methods  Taxa were categorized by dispersal mode and life-form and their distributions modelled using maxent . Geographical range size was compared amongst dispersal mode, life-form and biome using phylogenetically independent contrasts. Geographical range placement was considered in a similar manner.
Results  Range size did not vary with dispersal mode (ant versus wind and vertebrate dispersal) or life-form, and instead varied primarily as a function of the biogeographical region in which a species was found. Range placement, however, did vary among dispersal modes, with the consequence that diversity of wind- and ant-dispersed plants increased with latitude while the diversity of vertebrate-dispersed plants was more evenly distributed.
Main conclusions  For the taxa studied, range sizes were a function of the biogeographical region in which species were found. Although differences in range size may exist among species differing in dispersal modes, they are likely to be far smaller than differences among species from different biogeographical regions. The trait most likely to affect species geographical range size, and hence rarity and risks associated with other threats, may simply be the geographical region in which that species has evolved.     

Prophylactic and therapeutic use of antibodies for protection against respiratory infection with Francisella tularensis

The role of Abs in protection against respiratory infection with the intracellular bacterium Francisella tularensis is not clear. To investigate the ability of Abs to clear bacteria from the lungs and prevent systemic spread, immune serum was passively administered i.p. to naive mice before intranasal F. tularensis live vaccine strain infection. It was found that immune serum treatment provided 100% protection against lethal challenge while normal serum or Ig-depleted immune serum provided no protection. Protective efficacy was correlated with increased clearance of bacteria from the lung and required expression of FcgammaR on phagocytes, including macrophages and neutrophils. However, complement was not required for protection. In vitro experiments demonstrated that macrophages were more readily infected by Ab-opsonized bacteria but became highly efficient in killing upon activation by IFN-gamma. Consistent with this finding, in vivo Ab-mediated protection was found to be dependent upon IFN-gamma. SCID mice were not protected by passive Ab transfer, suggesting that T cells but not NK cells serve as the primary source for IFN-gamma. These data suggest that a critical interaction of humoral and cellular immune responses is necessary to provide sterilizing immunity against F. tularensis. Of considerable interest was the finding that serum Abs were capable of conferring protection against lethal respiratory tularemia when given 24-48 h postexposure. Thus, this study provides the first evidence for the therapeutic use of Abs in Francisella-infected individuals.     

Engineering of Saccharomyces cerevisiae for efficient anaerobic alcoholic fermentation of L-arabinose

For cost-effective and efficient ethanol production from lignocellulosic fractions of plant biomass, the conversion of not only major constituents, such as glucose and xylose, but also less predominant sugars, such as l-arabinose, is required. Wild-type strains of Saccharomyces cerevisiae, the organism used in industrial ethanol production, cannot ferment xylose and arabinose. Although metabolic and evolutionary engineering has enabled the efficient alcoholic fermentation of xylose under anaerobic conditions, the conversion of l-arabinose into ethanol by engineered S. cerevisiae strains has previously been demonstrated only under oxygen-limited conditions. This study reports the first case of fast and efficient anaerobic alcoholic fermentation of l-arabinose by an engineered S. cerevisiae strain. This fermentation was achieved by combining the expression of the structural genes for the l-arabinose utilization pathway of Lactobacillus plantarum, the overexpression of the S. cerevisiae genes encoding the enzymes of the nonoxidative pentose phosphate pathway, and extensive evolutionary engineering. The resulting S. cerevisiae strain exhibited high rates of arabinose consumption (0.70 g h(-1) g [dry weight](-1)) and ethanol production (0.29 g h(-1) g [dry weight](-1)) and a high ethanol yield (0.43 g g(-1)) during anaerobic growth on l-arabinose as the sole carbon source. In addition, efficient ethanol production from sugar mixtures containing glucose and arabinose, which is crucial for application in industrial ethanol production, was achieved.     

Habitat fragmentation,genetic diversity,and inbreeding depression in a threatened grassland legume: is genetic rescue necessary?

Kincaid’s lupine (Lupinus oreganus), a threatened perennial legume of western Oregon grasslands, is composed of small, fragmented populations that have consistently low natural seed set, suggesting they may have accumulated high enough levels of genetic load to be candidates for genetic rescue. We used simple sequence repeat (SSR) loci, both nuclear DNA and chloroplast DNA, to screen populations throughout the species’ range for evidence of severe inbreeding and recent genetic bottlenecks due to habitat fragmentation. After genotyping about 40% of the known populations, only one of 24 populations had strong statistical evidence for a recent genetic bottleneck (H _e > H _eq). Both mean nSSR fixation coefficients and genetic diversity did not statistically differ between very small, small, medium, and large lupine population size classes. Within population chloroplast DNA haplotype number was high for an animal pollinated species, ≈4.2 haplotypes/population, and within population haplotype diversity was also relatively evenly distributed. Within population patterns of nSSR and cpSSR genetic diversity suggest that genetic diversity has not been lost over the last century of habitat fragmentation. With genet lifespan thought to exceed 100 years, overlap of several to many generations, and substantial reductions in seed set from inbreeding depression that shifts cohort composition towards those generated by outcrossing events, Kincaid’s lupine is likely maintain the currently high levels of within population genetic diversity. The case of Kincaid’s lupine provides an example of how the assumptions of severe inbreeding depression with small population size and habitat fragmentation can be inaccurate.     

Tricyclic aminopyrimidine histamine H4 receptor antagonists

This report discloses the development of a series of tricyclic histamine H(4) receptor antagonists. Starting with a low nanomolar benzofuranopyrimidine HTS hit devoid of pharmaceutically acceptable properties, we navigated issues with metabolism and solubility to furnish a potent, stable and water soluble tricyclic histamine H(4) receptor antagonist with desirable physiochemical parameters which demonstrated efficacy a mouse ova model.