Search

Global ocean biogeochemistry models currently employed in climate change projections use highly simplified representations of pelagic food webs. These food webs do not necessarily include critical pathways by which ecosystems interact with ocean biogeochemistry and climate. Here we present a global biogeochemical model which incorporates ecosystem dynamics based on the representation of ten plankton functional types (PFTs): six types of phytoplankton, three types of zooplankton, and heterotrophic procaryotes. We improved the representation of zooplankton dynamics in our model through (a) the explicit inclusion of large, slow-growing macrozooplankton (e.g. krill), and (b) the introduction of trophic cascades among the three zooplankton types. We use the model to quantitatively assess the relative roles of iron vs. grazing in determining phytoplankton biomass in the Southern Ocean high-nutrient low-chlorophyll (HNLC) region during summer. When model simulations do not include macrozooplankton grazing explicitly, they systematically overestimate Southern Ocean chlorophyll biomass during the summer, even when there is no iron deposition from dust. When model simulations include a slow-growing macrozooplankton and trophic cascades among three zooplankton types, the high-chlorophyll summer bias in the Southern Ocean HNLC region largely disappears. Our model results suggest that the observed low phytoplankton biomass in the Southern Ocean during summer is primarily explained by the dynamics of the Southern Ocean zooplankton community, despite iron limitation of phytoplankton community growth rates. This result has implications for the representation of global biogeochemical cycles in models as zooplankton faecal pellets sink rapidly and partly control the carbon export to the intermediate and deep ocean.

Ocean biogeochemistry (OBGC) models span a wide variety of complexities, including highly simplified nutrient-restoring schemes, nutrient–phytoplankton–zooplankton–detritus (NPZD) models that crudely represent the marine biota, models that represent a broader trophic structure by grouping organisms as plankton functional types (PFTs) based on their biogeochemical role (dynamic green ocean models) and ecosystem models that group organisms by ecological function and trait. OBGC models are now integral components of Earth system models (ESMs), but they compete for computing resources with higher resolution dynamical setups and with other components such as atmospheric chemistry and terrestrial vegetation schemes. As such, the choice of OBGC in ESMs needs to balance model complexity and realism alongside relative computing cost. Here we present an intercomparison of six OBGC models that were candidates for implementation within the next UK Earth system model (UKESM1). The models cover a large range of biological complexity (from 7 to 57 tracers) but all include representations of at least the nitrogen, carbon, alkalinity and oxygen cycles. Each OBGC model was coupled to the ocean general circulation model Nucleus for European Modelling of the Ocean (NEMO) and results from physically identical hindcast simulations were compared. Model skill was evaluated for biogeochemical metrics of global-scale bulk properties using conventional statistical techniques. The computing cost of each model was also measured in standardised tests run at two resource levels. No model is shown to consistently outperform all other models across all metrics. Nonetheless, the simpler models are broadly closer to observations across a number of fields and thus offer a high-efficiency option for ESMs that prioritise high-resolution climate dynamics. However, simpler models provide limited insight into more complex marine biogeochemical processes and ecosystem pathways, and a parallel approach of low-resolution climate dynamics and high-complexity biogeochemistry is desirable in order to provide additional insights into biogeochemistry–climate interactions.

Accurate assessments of anthropogenic carbon dioxide (CO2) emissions and their redistribution among the atmosphere, ocean, and terrestrial biosphere is important to better understand the global carbon cycle, support the climate policy process, and project future climate change. Present-day analysis requires the combination of a range of data, algorithms, statistics and model estimates and their interpretation by a broad scientific community. Here we describe datasets and a methodology developed by the global carbon cycle science community to quantify all major components of the global carbon budget, including their uncertainties. We discuss changes compared to previous estimates, consistency within and among components, and methodology and data limitations. CO2 emissions from fossil fuel combustion and cement production (EFF) are based on energy statistics, while emissions from Land-Use Change (ELUC), including deforestation, are based on combined evidence from land cover change data, fire activity in regions undergoing deforestation, and models. The global atmospheric CO2 concentration is measured directly and its rate of growth (GATM) is computed from the concentration. The mean ocean CO2 sink (SOCEAN) is based on observations from the 1990s, while the annual anomalies and trends are estimated with ocean models. Finally, the global residual terrestrial CO2 sink (SLAND) is estimated by the difference of the other terms. For the last decade available (2002–2011), EFF was 8.3 ± 0.4 PgC yr−1, ELUC 1.0 ± 0.5 PgC yr−1, GATM 4.3 ± 0.1PgC yr−1, SOCEAN 2.5 ± 0.5 PgC yr−1, and SLAND 2.6 ± 0.8 PgC yr−1. For year 2011 alone, EFF was 9.5 ± 0.5 PgC yr−1, 3.0 percent above 2010, reflecting a continued trend in these emissions; ELUC was 0.9 ± 0.5 PgC yr−1, approximately constant throughout the decade; GATM was 3.6 ± 0.2 PgC yr−1, SOCEAN was 2.7 ± 0.5 PgC yr−1, and SLAND was 4.1 ± 0.9 PgC yr−1. GATM was low in 2011 compared to the 2002–2011 average because of a high uptake by the land probably in response to natural climate variability associated to La Niña conditions in the Pacific Ocean. The global atmospheric CO2 concentration reached 391.31 ± 0.13 ppm at the end of year 2011. We estimate that EFF will have increased by 2.6% (1.9–3.5%) in 2012 based on projections of gross world product and recent changes in the carbon intensity of the economy. All uncertainties are reported as ±1 sigma (68% confidence assuming Gaussian error distributions that the real value lies within the given interval), reflecting the current capacity to characterise the annual estimates of each component of the global carbon budget. This paper is intended to provide a baseline to keep track of annual carbon budgets in the future. All data presented here can be downloaded from the Carbon Dioxide Information Analysis Center (doi:10.3334/CDIAC/GCP_V2013).

Phage therapy recently passed a key milestone with success of the first regulated clinical trial using systemic administration. In this single-arm non-comparative safety study, phages were administered intravenously to patients with invasive Staphylococcus aureus infections with no adverse reactions reported. Here, we examined features of 78 lytic S. aureus phages, most of which were propagated using a S. carnosus host modified to be broadly susceptible to staphylococcal phage infection. Use of this host eliminates the threat of contamination with staphylococcal prophage — the main vector of S. aureus horizontal gene transfer. We determined the host range of these phages against an international collection of 185 S. aureus isolates with 56 different multilocus sequence types that included multiple representatives of all epidemic MRSA and MSSA clonal complexes. Forty of our 78 phages were able to infect > 90% of study isolates, 15 were able to infect > 95%, and two could infect all 184 clinical isolates, but not a phage-resistant mutant generated in a previous study. We selected the 10 phages with the widest host range for in vitro characterization by planktonic culture time-kill analysis against four isolates:- modified S. carnosus strain TM300H, methicillin-sensitive isolates D329 and 15981, and MRSA isolate 252. Six of these 10 phages were able to rapidly kill, reducing cell numbers of at least three isolates. The four best-performing phages, in this assay, were further shown to be highly effective in reducing 48 h biofilms on polystyrene formed by eight ST22 and eight ST36 MRSA isolates. Genomes of 22 of the widest host-range phages showed they belonged to the Twortvirinae subfamily of the order Caudovirales in three main groups corresponding to Silviavirus, and two distinct groups of Kayvirus. These genomes assembled as single-linear dsDNAs with an average length of 140 kb and a GC content of c. 30%. Phages that could infect > 96% of S. aureus isolates were found in all three groups, and these have great potential as therapeutic candidates if, in future studies, they can be formulated to maximize their efficacy and eliminate emergence of phage resistance by using appropriate combinations.

Capsular polysaccharides enable clinically important clones of Klebsiella pneumoniae to cause severe systemic infections in susceptible hosts. Phage-encoded capsule depolymerases have the potential to provide an alternative treatment paradigm in patients when multiple drug resistance has eroded the efficacy of conventional antibiotic chemotherapy. An investigation of 164 K. pneumoniae from intensive care patients in Thailand revealed a large number of distinct K types in low abundance but four (K2, K51, K1, K10) with a frequency of at least 5%. To identify depolymerases with the capacity to degrade capsules associated with these common K-types, 62 lytic phage were isolated from Thai hospital sewage water using K1, K2 and K51 isolates as hosts; phage plaques, without exception, displayed halos indicative of the presence of capsule-degrading enzymes. Phage genomes ranged in size from 41–348 kb with between 50 and 535 predicted coding sequences (CDSs). Using a custom phage protein database we were successful in applying annotation to 30 - 70% (mean = 58%) of these CDSs. The largest genomes, of so-called jumbo phage, carried multiple tRNAs as well as CRISPR repeat and spacer sequences. One of the smaller phage genomes was found to contain a putative Cas type 1E gene, indicating a history of host DNA acquisition in these obligate lytic phage. Whole-genome sequencing (WGS) indicated that some phage displayed an extended host range due to the presence of multiple depolymerase genes; in total, 42 candidate depolymerase genes were identified with up to eight in a single genome. Seven distinct virions were selected for further investigation on the basis of host range, phage morphology and WGS. Candidate genes for K1, K2 and K51 depolymerases were expressed and purified as his6-tagged soluble protein and enzymatic activity demonstrated against K. pneumoniae capsular polysaccharides by gel electrophoresis and Anton-Paar rolling ball viscometry. Depolymerases completely removed the capsule in K-type-specific fashion from K. pneumoniae cells. We conclude that broad-host range phage carry multiple enzymes, each with the capacity to degrade a single K-type, and any future use of these enzymes as therapeutic agents will require enzyme cocktails for utility against a range of K. pneumoniae infections.

Nineteen bacteriophages against five main capsular types of multidrug-resistant Acinetobacter baumannii were isolated from tertiary care hospital sewage. Eight representative phages from each capsular type were characterized and tested for their biological properties. The biological features revealed that phages T1245, T444, and T515 had a large burst size of more than 420 pfu/mL, together with a short latent period lasting less than 6 min, and were readily adsorbed to a bacterial host within 10 min. Moreover, these phages demonstrated host specificity and stability over a broad range of temperatures (−20 to 60 °C) and pH (5.0–9.0). A whole-genome analysis of six lytic and two temperate phages revealed high genomic similarity with double-stranded DNA between 40 and 50 kb and G + C content of 38–39%. The protein compositions disclosed the absence of toxin-coding genes. The phylogenic results, together with morphological micrographs, confirmed that three selected phages (T1245, T444, and T515) belong to the Podoviridae family within the order Caudovirales. The biological data and bioinformatics analysis indicated that these novel A. baumannii phages possess important enzymes, including depolymerase and endolysin, which could be further developed as promising alternative antibacterial agents to control A. baumannii infections.

ABSTRACT Pseudomonas aeruginosa is an antibiotic-refractory pathogen with a large genome and extensive genotypic diversity. Historically, P. aeruginosa has been a major model system for understanding the molecular mechanisms underlying type I clustered regularly interspaced short palindromic repeat (CRISPR) and CRISPR-associated protein (CRISPR-Cas)-based bacterial immune system function. However, little information on the phylogenetic distribution and potential role of these CRISPR-Cas systems in molding the P. aeruginosa accessory genome and antibiotic resistance elements is known. Computational approaches were used to identify and characterize CRISPR-Cas systems within 672 genomes, and in the process, we identified a previously unreported and putatively mobile type I-C P. aeruginosa CRISPR-Cas system. Furthermore, genomes harboring noninhibited type I-F and I-E CRISPR-Cas systems were on average ~300 kb smaller than those without a CRISPR-Cas system. In silico analysis demonstrated that the accessory genome (n = 22,036 genes) harbored the majority of identified CRISPR-Cas targets. We also assembled a global spacer library that aided the identification of difficult-to-characterize mobile genetic elements within next-generation sequencing (NGS) data and allowed CRISPR typing of a majority of P. aeruginosa strains. In summary, our analysis demonstrated that CRISPR-Cas systems play an important role in shaping the accessory genomes of globally distributed P. aeruginosa isolates. IMPORTANCE P. aeruginosa is both an antibiotic-refractory pathogen and an important model system for type I CRISPR-Cas bacterial immune systems. By combining the genome sequences of 672 newly and previously sequenced genomes, we were able to provide a global view of the phylogenetic distribution, conservation, and potential targets of these systems. This analysis identified a new and putatively mobile P. aeruginosa CRISPR-Cas subtype, characterized the diverse distribution of known CRISPR-inhibiting genes, and provided a potential new use for CRISPR spacer libraries in accessory genome analysis. Our data demonstrated the importance of CRISPR-Cas systems in modulating the accessory genomes of globally distributed strains while also providing substantial data for subsequent genomic and experimental studies in multiple fields. Understanding why certain genotypes of P. aeruginosa are clinically prevalent and adept at horizontally acquiring virulence and antibiotic resistance elements is of major clinical and economic importance.

Infections caused by community-acquired (CA)-methicillin resistant Staphylococcus aureus (MRSA) have been reported worldwide. We assessed whether any common genetic markers existed among 117 CA-MRSA isolates from the United States, France, Switzerland, Australia, New Zealand, and Western Samoa by performing polymerase chain reaction for 24 virulence factors and the methicillin-resistance determinant. The genetic background of the strain was analyzed by pulsed-field gel electrophoresis (PFGE) and multi-locus sequence typing (MLST). The CA-MRSA strains shared a type IV SCCmec cassette and the Panton-Valentine leukocidin locus, whereas the distribution of the other toxin genes was quite specific to the strains from each continent. PFGE and MLST analysis indicated distinct genetic backgrounds associated with each geographic origin, although predominantly restricted to the agr3 background. Within each continent, the genetic background of CA-MRSA strains did not correspond to that of the hospital-acquired MRSA.

There is a need for a regional assessment of the frequency and diversity of MRSA to determine major circulating clones and the extent to which community and healthcare MRSA reservoirs have mixed. We conducted a prospective cohort study of inpatients in Orange County, California, systematically collecting clinical MRSA isolates from 30 hospitals, to assess MRSA diversity and distribution. All isolates were characterized by spa typing, with selective PFGE and MLST to relate spa types with major MRSA clones. We collected 2,246 MRSA isolates from hospital inpatients. This translated to 91/10,000 inpatients with MRSA and an Orange County population estimate of MRSA inpatient clinical cultures of 86/100,000 people. spa type genetic diversity was heterogeneous between hospitals, and relatively high overall (72%). USA300 (t008/ST8), USA100 (t002/ST5) and a previously reported USA100 variant (t242/ST5) were the dominant clones across all Orange County hospitals, representing 83% of isolates. Fifteen hospitals isolated more t008 (USA300) isolates than t002/242 (USA100) isolates, and 12 hospitals isolated more t242 isolates than t002 isolates. The majority of isolates were imported into hospitals. Community-based infection control strategies may still be helpful in stemming the influx of traditionally community-associated strains, particularly USA300, into the healthcare setting.

Methicillin-resistant Staphylococcus aureus (MRSA) with reduced vancomycin susceptibility (VISA, vancomycin-intermediate S. aureus) has been reported from many countries. Whether resistance is evolving regularly in different genetic backgrounds or in a single clone with a genetic predisposition, as early results suggest, is unclear. We have studied 101 MRSA with reduced vancomycin susceptibility from nine countries by multilocus sequence typing (MLST) and characterization of SCCmec (staphylococcal chromosomal cassette mec) and agr (accessory gene regulator). We found nine genotypes by MLST, with isolates within all five major hospital MRSA lineages. Most isolates (88/101) belonged to two of the earliest MRSA clones that have global prevalence. Our results show that reduced susceptibility to vancomycin has emerged in many successful epidemic lineages with no clear clonal disposition. Increasing antimicrobial resistance in genetically distinct pandemic clones may lead to MRSA infections that will become increasingly difficult to treat.

Staphylococcus aureus bacteraemia (SAB) is commonly complicated by metastatic infection or relapse after treatment. Objectives. The study aim was to determine the role of bacterial, host, and management factors in development of complicated SAB. Methods. A prospectively-conducted observational study gathered data on predisposition, management and outcome of 100 consecutive SAB cases. Antibiotic susceptibilities and genetic lineage of bacterial isolates were determined. Further clinical and microbiological data were gathered on two retrospective series from 1999–2000 (n=57) and 2004 (n=116). Results. In the prospective cases, 27% met our definition of complicated disease. Expressed as RR and 95% CI, complicated disease was associated with diabetes (1.58, 1.00–2.48), injecting-drug use (5.48, 0.88–33.49), community-onset of symptoms (1.4, 1.02–1.92), and symptom duration ≥48 hours prior to starting effective antibiotic therapy (2.10, 1.22–3.61). Uncomplicated disease was associated with the presence of a central line (0.69, 0.55–0.88) and prompt removal of a primary focus (0.71, 0.57–0.90). Neither methicillin resistance nor genetic lineage was associated with complicated disease, but methicillin resistance was associated with higher mortality. Conclusions. This study demonstrates that clinical rather than microbial factors are the major determinants of SAB outcome and underscores the importance of early treatment.

The methicillin-resistant Staphylococcus aureus (MRSA) clonal group known as ST239-MRSA-III is notable for its hybrid origin and for causing sustained hospital epidemics worldwide since the late 1970s. We studied the population structure of this MRSA clonal group using a sample of 111 isolates that were collected over 34 years from 29 countries. Genetic variation was assessed using typing methods and novel ascertainment methods, resulting in approximately 15 kb of sequence from 32 loci for all isolates. A single most parsimonious tree, free of homoplasy, partitioned 28 haplotypes into geographically-associated clades, including prominent European, Asian, and South American clades. The rate of evolution was estimated to be approximately 100x faster than standard estimates for bacteria, and dated the most recent common ancestor of these isolates to the mid-20th century. Associations were discovered between the ST239 phylogeny and the ccrB and dru loci of the methicillin resistance genetic element, SCCmec type III, but not with the accessory components of the element that are targeted by multiplex PCR subtyping tools. In summary, the evolutionary history of ST239 can be characterized by rapid clonal diversification that has left strong evidence of geographic and temporal population structure. SCCmec type III has remained linked to the ST239 chromosome during clonal diversification, but it has undergone homoplasious losses of accessory components. These results provide a population genetics framework for the precise identification of emerging ST239 variants, and invite a re-evaluation of the markers used for subtyping SCCmec.

We determined the agr type, multilocus sequence type, protein A gene type (spa typing), toxin gene profile, and antimicrobial drug resistance profile of 469 isolates of Panton-Valentine leukocidin–positive community-acquired methicillin-resistant Staphylococcus aureus isolates (PVL-positive CA-MRSA). The isolates had been collected from around the world from 1999 through 2005 by the French National Reference Center for Staphylococci. We found that some continent-specific clones described in 2003, such as clone ST8, have now spread all over the world. Likewise, some PVL-positive CA-MRSA have spread to several countries on various continents. New clones have emerged (e.g., ST377) on new genetic backgrounds. PVL-positive CA-MRSA that were usually susceptible to most antistaphylococcal antimicrobial agents have acquired new resistance determinants (e.g., to gentamicin) in certain countries. The major trait shared by all these clones is a short staphylococcal chromosomal cassette mec element of type IV or V.