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.     

A mouse-adapted SARS-coronavirus causes disease and mortality in BALB/c mice

No single animal model for severe acute respiratory syndrome (SARS) reproduces all aspects of the human disease. Young inbred mice support SARS-coronavirus (SARS-CoV) replication in the respiratory tract and are available in sufficient numbers for statistical evaluation. They are relatively inexpensive and easily accessible, but their use in SARS research is limited because they do not develop illness following infection. Older (12- to 14-mo-old) BALB/c mice develop clinical illness and pneumonitis, but they can be hard to procure, and immune senescence complicates pathogenesis studies. We adapted the SARS-CoV (Urbani strain) by serial passage in the respiratory tract of young BALB/c mice. Fifteen passages resulted in a virus (MA15) that is lethal for mice following intranasal inoculation. Lethality is preceded by rapid and high titer viral replication in lungs, viremia, and dissemination of virus to extrapulmonary sites accompanied by lymphopenia, neutrophilia, and pathological changes in the lungs. Abundant viral antigen is extensively distributed in bronchial epithelial cells and alveolar pneumocytes, and necrotic cellular debris is present in airways and alveoli, with only mild and focal pneumonitis. These observations suggest that mice infected with MA15 die from an overwhelming viral infection with extensive, virally mediated destruction of pneumocytes and ciliated epithelial cells. The MA15 virus has six coding mutations associated with adaptation and increased virulence; when introduced into a recombinant SARS-CoV, these mutations result in a highly virulent and lethal virus (rMA15), duplicating the phenotype of the biologically derived MA15 virus. Intranasal inoculation with MA15 reproduces many aspects of disease seen in severe human cases of SARS. The availability of the MA15 virus will enhance the use of the mouse model for SARS because infection with MA15 causes morbidity, mortality, and pulmonary pathology. This virus will be of value as a stringent challenge in evaluation of the efficacy of vaccines and antivirals.     

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.     

Homogeneous and heterogeneous tertiary structure ensembles of amyloid-β peptides

The interplay of modern molecular simulation and high-quality nuclear magnetic resonance (NMR) experiments has reached a fruitful stage for quantitative characterization of structural ensembles of disordered peptides. Amyloid-β 1-42 (Aβ42), the primary peptide associated with Alzheimer's disease, and fragments such as Aβ21-30 are both classified as intrinsically disordered peptides (IDPs). We use a variety of NMR observables to validate de novo molecular dynamics simulations in explicit water to characterize the tertiary structure ensemble of Aβ42 and Aβ21-30 from the perspective of their classification as IDPs. Unlike the Aβ21-30 fragment that conforms to expectations of an IDP that is primarily extended, we find that Aβ42 samples conformations reflecting all possible secondary structure categories and spans the range of IDP classifications from collapsed structured states to highly extended conformations, making it an IDP with a far more heterogeneous tertiary ensemble.