Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences

Profiling phylogenetic marker genes, such as the 16S rRNA gene, is a key tool for studies of microbial communities but does not provide direct evidence of a community’s functional capabilities. Here we describe PICRUSt (Phylogenetic Investigation of Communities by Reconstruction of Unobserved States), a computational approach to predict the functional composition of a metagenome using marker gene data and a database of reference genomes. PICRUSt uses an extended ancestral-state reconstruction algorithm to predict which gene families are present and then combines gene families to estimate the composite metagenome. Using 16S information, PICRUSt recaptures key findings from the Human Microbiome Project and accurately predicts the abundance of gene families in host-associated and environmental communities, with quantifiable uncertainty. Our results demonstrate that phylogeny and function are sufficiently linked that this ‘predictive metagenomic’ approach should provide useful insights into the thousands of uncultivated microbial communities for which only marker gene surveys are currently available

Estimating Impact Forces of Tail Club Strikes by Ankylosaurid Dinosaurs

BACKGROUND: It has been assumed that the unusual tail club of ankylosaurid dinosaurs was used actively as a weapon, but the biological feasibility of this behaviour has not been examined in detail. Ankylosaurid tail clubs are composed of interlocking vertebrae, which form the handle, and large terminal osteoderms, which form the knob. METHODOLOGY/PRINCIPAL FINDINGS: Computed tomographic (CT) scans of several ankylosaurid tail clubs referred to Dyoplosaurus and Euoplocephalus, combined with measurements of free caudal vertebrae, provide information used to estimate the impact force of tail clubs of various sizes. Ankylosaurid tails are modeled as a series of segments for which mass, muscle cross-sectional area, torque, and angular acceleration are calculated. Free caudal vertebrae segments had limited vertical flexibility, but the tail could have swung through approximately 100 degrees laterally. Muscle scars on the pelvis record the presence of a large M. longissimus caudae, and ossified tendons alongside the handle represent M. spinalis. CT scans showed that knob osteoderms were predominantly cancellous, which would have lowered the rotational inertia of the tail club and made it easier to wield as a weapon. CONCLUSIONS/SIGNIFICANCE: Large knobs could generate sufficient force to break bone during impacts, but average and small knobs could not. Tail swinging behaviour is feasible in ankylosaurids, but it remains unknown whether the tail was used for interspecific defense, intraspecific combat, or both

Substrate Type Determines Metagenomic Profiles from Diverse Chemical Habitats

Environmental parameters drive phenotypic and genotypic frequency variations in microbial communities and thus control the extent and structure of microbial diversity. We tested the extent to which microbial community composition changes are controlled by shifting physiochemical properties within a hypersaline lagoon. We sequenced four sediment metagenomes from the Coorong, South Australia from samples which varied in salinity by 99 Practical Salinity Units (PSU), an order of magnitude in ammonia concentration and two orders of magnitude in microbial abundance. Despite the marked divergence in environmental parameters observed between samples, hierarchical clustering of taxonomic and metabolic profiles of these metagenomes showed striking similarity between the samples (>89%). Comparison of these profiles to those derived from a wide variety of publically available datasets demonstrated that the Coorong sediment metagenomes were similar to other sediment, soil, biofilm and microbial mat samples regardless of salinity (>85% similarity). Overall, clustering of solid substrate and water metagenomes into discrete similarity groups based on functional potential indicated that the dichotomy between water and solid matrices is a fundamental determinant of community microbial metabolism that is not masked by salinity, nutrient concentration or microbial abundance

The Comparative Osteology of the Petrotympanic Complex (Ear Region) of Extant Baleen Whales (Cetacea: Mysticeti)

Anatomical comparisons of the ear region of baleen whales (Mysticeti) are provided through detailed osteological descriptions and high-resolution photographs of the petrotympanic complex (tympanic bulla and petrosal bone) of all extant species of mysticete cetaceans. Salient morphological features are illustrated and identified, including overall shape of the bulla, size of the conical process of the bulla, morphology of the promontorium, and the size and shape of the anterior process of the petrosal. We place our comparative osteological observations into a phylogenetic context in order to initiate an exploration into petrotympanic evolution within Mysticeti.The morphology of the petrotympanic complex is diagnostic for individual species of baleen whale (e.g., sigmoid and conical processes positioned at midline of bulla in Balaenoptera musculus; confluence of fenestra cochleae and perilymphatic foramen in Eschrichtius robustus), and several mysticete clades are united by derived characteristics. Balaenids and neobalaenids share derived features of the bulla, such as a rhomboid shape and a reduced anterior lobe (swelling) in ventral aspect, and eschrichtiids share derived morphologies of the petrosal with balaenopterids, including loss of a medial promontory groove and dorsomedial elongation of the promontorium. Monophyly of Balaenoidea (Balaenidae and Neobalaenidae) and Balaenopteroidea (Balaenopteridae and Eschrichtiidae) was recovered in phylogenetic analyses utilizing data exclusively from the petrotympanic complex.This study fills a major gap in our knowledge of the complex structures of the mysticete petrotympanic complex, which is an important anatomical region for the interpretation of the evolutionary history of mammals. In addition, we introduce a novel body of phylogenetically informative characters from the ear region of mysticetes. Our detailed anatomical descriptions, illustrations, and comparisons provide valuable data for current and future studies on the phylogenetic relationships, evolution, and auditory physiology of mysticetes and other cetaceans throughout Earth's history

Precision measurement of the top quark mass from dilepton events at CDF II

We report a measurement of the top quark mass, M_t, in the dilepton decay channel of $t\bar{t}\to b\ell'^{+}\nu_{\ell'}\bar{b}\ell^{-}\bar{\nu}_{\ell}$ using an integrated luminosity of 1.0 fb^{-1} of p\bar{p} collisions collected with the CDF II detector. We apply a method that convolutes a leading-order matrix element with detector resolution functions to form event-by-event likelihoods; we have enhanced the leading-order description to describe the effects of initial-state radiation. The joint likelihood is the product of the likelihoods from 78 candidate events in this sample, which yields a measurement of M_{t} = 164.5 \pm 3.9(\textrm{stat.}) \pm 3.9(\textrm{syst.}) \mathrm{GeV}/c^2, the most precise measurement of M_t in the dilepton channel.Comment: 7 pages, 2 figures, version includes changes made prior to publication by journa

Measurement of the Ratios of Branching Fractions B(Bs -> Ds pi pi pi) / B(Bd -> Dd pi pi pi) and B(Bs -> Ds pi) / B(Bd -> Dd pi)

Using 355 pb^-1 of data collected by the CDF II detector in \ppbar collisions at sqrt{s} = 1.96 TeV at the Fermilab Tevatron, we study the fully reconstructed hadronic decays B -> D pi and B -> D pi pi pi. We present the first measurement of the ratio of branching fractions B(Bs -> Ds pi pi pi) / B(Bd -> Dd pi pi pi) = 1.05 pm 0.10 (stat) pm 0.22 (syst). We also update our measurement of B(Bs -> Ds pi) / B(Bd -> Dd pi) to 1.13 pm 0.08 (stat) pm 0.23 (syst) improving the statistical uncertainty by more than a factor of two. We find B(Bs -> Ds pi) = [3.8 pm 0.3 (stat) pm 1.3 (syst)] \times 10^{-3} and B(Bs -> Ds pi pi pi) = [8.4 pm 0.8 (stat) pm 3.2 (syst)] \times 10^{-3}.Comment: 7 pages, 2 figure

Cross Section Measurements of High-pTp_T Dilepton Final-State Processes Using a Global Fitting Method

We present a new method for studying high-$p_T$ dilepton events ($e^{\pm}e^{\mp}$, $\mu^{\pm}\mu^{\mp}$, $e^{\pm}\mu^{\mp}$) and simultaneously extracting the production cross sections of $p\bar{p} \to t\bar{t}$, $p\bar{p} \to W^+W^-$, and p\bar{p} \to \ztt at a center-of-mass energy of $\sqrt{s} = 1.96$ TeV. We perform a likelihood fit to the dilepton data in a parameter space defined by the missing transverse energy and the number of jets in the event. Our results, which use $360 {\rm pb^{-1}}$ of data recorded with the CDF II detector at the Fermilab Tevatron Collider, are $\sigma(t\bar{t}) = 8.5_{-2.2}^{+2.7}$ pb, $\sigma(W^+W^-) = 16.3^{+5.2}_{-4.4}$ pb, and \sigma(\ztt) =291^{+50}_{-46} pb.Comment: 20 pages, 2 figures, to be submitted to PRD-R

Top Quark Mass Measurement from Dilepton Events at CDF II with the Matrix-Element Method

We describe a measurement of the top quark mass using events with two charged leptons collected by the CDF II detector from $p\bar{p}$ collisions with $\sqrt s = 1.96$ TeV at the Fermilab Tevatron. The likelihood in top mass is calculated for each event by convoluting the leading order matrix element describing $q\bar{q} \to t\bar{t} \to b\ell\nu_{\ell}\bar{b}\ell'\nu_{\ell'}$ with detector resolution functions. The presence of background events in the data sample is modeled using similar calculations involving the matrix elements for major background processes. In a data sample with integrated luminosity of 340 pb$^{-1}$, we observe 33 candidate events and measure $M_{top} = 165.2 \pm 6.1(\textrm{stat.}) \pm 3.4(\textrm{syst.}) \mathrm{~GeV}/c^2.$ This measurement represents the first application of this method to events with two charged leptons and is the most precise single measurement of the top quark mass in this channel.Comment: 21 pages, 14 figure

Search for New Physics in Lepton + Photon + X Events with L=305 pb-1 of ppbar Collisions at roots=1.96 TeV

We present results of a search for anomalous production of events containing a charged lepton (either electron or muon) and a photon, both with high transverse momentum, accompanied by additional signatures, X, including missing transverse energy (MET) and additional leptons and photons. We use the same kinematic selection criteria as in a previous CDF search, but with a substantially larger data set, 305 pb-1, a ppbar collision energy of 1.96 TeV, and the upgraded CDF II detector. We find 42 Lepton+Photon+MET events versus a standard model expectation of 37.3 +- 5.4 events. The level of excess observed in Run I, 16 events with an expectation of 7.6 +- 0.7 events (corresponding to a 2.7 sigma effect), is not supported by the new data. In the signature of Multi-Lepton+Photon+X we observe 31 events versus an expectation of 23.0 +- 2.7 events. In this sample we find no events with an extra photon or MET and so find no events like the one ee+gg+MET event observed in Run I.Comment: 7 pages, 3 figures, 1 table. Accepted to PR

Measurement of the Lambda_b Lifetime in Lambda_b --> J/psi Lambda0 in p-pbar Collisions at sqrt(s)=1.96 TeV

We report a measurement of the Lambda_b lifetime in the exclusive decay Lambda_b --> J/psi Lambda0 in p-pbar collisions at sqrt(s) = 1.96 TeV using an integrated luminosity of 1.0 fb^{-1} of data collected by the CDF II detector at the Fermilab Tevatron. Using fully reconstructed decays, we measure tau(Lambda_b) = 1.593 ^{+0.083}_{-0.078} (stat.) +- 0.033 (syst.) ps. This is the single most precise measurement of tau(Lambda_b) and is 3.2 sigma higher than the current world average.Comment: 7 Pages, 2 Figures, 1 Table. Submitted to Phys. Rev. Let