New statistical method identifies cytokines that distinguish stool microbiomes.

Regressing an outcome or dependent variable onto a set of input or independent variables allows the analyst to measure associations between the two so that changes in the outcome can be described by and predicted by changes in the inputs. While there are many ways of doing this in classical statistics, where the dependent variable has certain properties (e.g., a scalar, survival time, count), little progress on regression where the dependent variable are microbiome taxa counts has been made that do not impose extremely strict conditions on the data. In this paper, we propose and apply a new regression model combining the Dirichlet-multinomial distribution with recursive partitioning providing a fully non-parametric regression model. This model, called DM-RPart, is applied to cytokine data and microbiome taxa count data and is applicable to any microbiome taxa count/metadata, is automatically fit, and intuitively interpretable. This is a model which can be applied to any microbiome or other compositional data and software (R package HMP) available through the R CRAN website

Evaluation of 16S rRNA gene sequencing for species and strain-level microbiome analysis.

The 16S rRNA gene has been a mainstay of sequence-based bacterial analysis for decades. However, high-throughput sequencing of the full gene has only recently become a realistic prospect. Here, we use in silico and sequence-based experiments to critically re-evaluate the potential of the 16S gene to provide taxonomic resolution at species and strain level. We demonstrate that targeting of 16S variable regions with short-read sequencing platforms cannot achieve the taxonomic resolution afforded by sequencing the entire (~1500 bp) gene. We further demonstrate that full-length sequencing platforms are sufficiently accurate to resolve subtle nucleotide substitutions (but not insertions/deletions) that exist between intragenomic copies of the 16S gene. In consequence, we argue that modern analysis approaches must necessarily account for intragenomic variation between 16S gene copies. In particular, we demonstrate that appropriate treatment of full-length 16S intragenomic copy variants has the potential to provide taxonomic resolution of bacterial communities at species and strain level

Extensive Gene Amplification as a Mechanism for Piperacillin-Tazobactam Resistance in Escherichia coli.

Although the TEM-1 β-lactamase (BlaTEM-1) hydrolyzes penicillins and narrow-spectrum cephalosporins, organisms expressing this enzyme are typically susceptible to β-lactam/β-lactamase inhibitor combinations such as piperacillin-tazobactam (TZP). However, our previous work led to the discovery of 28 clinical isolates of Escherichia coli resistant to TZP that contained only blaTEM-1 One of these isolates, E. coli 907355, was investigated further in this study. E. coli 907355 exhibited significantly higher β-lactamase activity and BlaTEM-1 protein levels when grown in the presence of subinhibitory concentrations of TZP. A corresponding TZP-dependent increase in blaTEM-1 copy number was also observed, with as many as 113 copies of the gene detected per cell. These results suggest that TZP treatment promotes an increase in blaTEM-1 gene dosage, allowing BlaTEM-1 to reach high enough levels to overcome inactivation by the available tazobactam in the culture. To better understand the nature of the blaTEM-1 copy number proliferation, whole-genome sequence (WGS) analysis was performed on E. coli 907355 in the absence and presence of TZP. The WGS data revealed that the blaTEM-1 gene is located in a 10-kb genomic resistance module (GRM) that contains multiple resistance genes and mobile genetic elements. The GRM was found to be tandemly repeated at least 5 times within a p1ESCUM/p1ECUMN-like plasmid when bacteria were grown in the presence of TZP.IMPORTANCE Understanding how bacteria acquire resistance to antibiotics is essential for treating infected patients effectively, as well as preventing the spread of resistant organisms. In this study, a clinical isolate of E. coli was identified that dedicated more than 15% of its genome toward tandem amplification of a ~10-kb resistance module, allowing it to escape antibiotic-mediated killing. Our research is significant in that it provides one possible explanation for clinical isolates that exhibit discordant behavior when tested for antibiotic resistance by different phenotypic methods. Our research also shows that GRM amplification is difficult to detect by short-read WGS technologies. Analysis of raw long-read sequence data was required to confirm GRM amplification as a mechanism of antibiotic resistance. MBio 2018 Apr 24; 9(2):e00583-18

Global epidemiology and clinical outcomes of carbapenem-resistant Pseudomonas aeruginosa and associated carbapenemases (POP): a prospective cohort study

BACKGROUND: Carbapenem-resistant Pseudomonas aeruginosa (CRPA) is a global threat, but the distribution and clinical significance of carbapenemases are unclear. The aim of this study was to define characteristics and outcomes of CRPA infections and the global frequency and clinical impact of carbapenemases harboured by CRPA. METHODS: We conducted an observational, prospective cohort study of CRPA isolated from bloodstream, respiratory, urine, or wound cultures of patients at 44 hospitals (10 countries) between Dec 1, 2018, and Nov 30, 2019. Clinical data were abstracted from health records and CRPA isolates were whole-genome sequenced. The primary outcome was 30-day mortality from the day the index culture was collected. We compared outcomes of patients with CRPA infections by infection type and across geographic regions and performed an inverse probability weighted analysis to assess the association between carbapenemase production and 30-day mortality. FINDINGS: We enrolled 972 patients (USA n=527, China n=171, south and central America n=127, Middle East n=91, Australia and Singapore n=56), of whom 581 (60%) had CRPA infections. 30-day mortality differed by infection type (bloodstream 21 [30%] of 69, respiratory 69 [19%] of 358, wound nine [14%] of 66, urine six [7%] of 88; p=0·0012) and geographical region (Middle East 15 [29%] of 52, south and central America 20 [27%] of 73, USA 60 [19%] of 308, Australia and Singapore three [11%] of 28, China seven [6%] of 120; p=0·0002). Prevalence of carbapenemase genes among CRPA isolates also varied by region (south and central America 88 [69%] of 127, Australia and Singapore 32 [57%] of 56, China 54 [32%] of 171, Middle East 27 [30%] of 91, USA ten [2%] of 527; p\u3c0·0001). KPC-2 (n=103 [49%]) and VIM-2 (n=75 [36%]) were the most common carbapenemases in 211 carbapenemase-producing isolates. After excluding USA patients, because few US isolates had carbapenemases, patients with carbapenemase-producing CRPA infections had higher 30-day mortality than those with non-carbapenemase-producing CRPA infections in both unadjusted (26 [22%] of 120 vs 19 [12%] of 153; difference 9%, 95% CI 3-16) and adjusted (difference 7%, 95% CI 1-14) analyses. INTERPRETATION: The emergence of different carbapenemases among CRPA isolates in different geographical regions and the increased mortality associated with carbapenemase-producing CRPA infections highlight the therapeutic challenges posed by these organisms. FUNDING: National Institutes of Health

Epidemiología mundial y resultados clínicos de Pseudomonas aeruginosa resistente a carbapenemes y carbapenemasas asociadas (POP): un estudio prospectivo de cohortes

Antecedentes: La Pseudomonas aeruginosa resistente a los carbapenemes (CRPA) es una amenaza mundial, pero la distribución y la importancia clínica de las carbapenemasas no están claras. El objetivo de este estudio fue definir las características y los resultados de las infecciones por CRPA, así como la frecuencia global y el impacto clínico de las carbapenemasas albergadas por CRPA. Métodos: Llevamos a cabo un estudio de cohortes observacional y prospectivo de CRPA aislados de cultivos de torrente sanguíneo, respiratorio, orina o heridas de pacientes en 44 hospitales (10 países) entre el 1 de diciembre de 2018 y el 30 de noviembre de 2019. Los datos clínicos se extrajeron de los registros de salud y los aislados de CRPA se secuenciaron en todo el genoma. El resultado primario fue la mortalidad a 30 días a partir del día en que se recolectó el cultivo índice. Se compararon los resultados de los pacientes con infecciones por CRPA por tipo de infección y entre regiones geográficas y se realizó un análisis ponderado de probabilidad inversa para evaluar la asociación entre la producción de carbapenemasas y la mortalidad a 30 días. Resultados: Se incluyeron 972 pacientes (EE.UU. n=527, China n=171, América del Sur y Central n=127, Oriente Medio n=91, Australia y Singapur n=56), de los cuales 581 (60%) tenían infecciones por CRPA. La mortalidad a los 30 días difería según el tipo de infección (torrente sanguíneo 21 [30%] de 69, respiratoria 69 [19%] de 358, herida nueve [14%] de 66, orina seis [7%] de 88; p=0-0012) y la región geográfica (Oriente Medio 15 [29%] de 52, América del Sur y Central 20 [27%] de 73, EE.UU. 60 [19%] de 308, Australia y Singapur tres [11%] de 28, China siete [6%] de 120; p=0-0002). La prevalencia de genes carbapenemasa entre los aislados CRPA también varió según la región (América del Sur y Central 88 [69%] de 127, Australia y Singapur 32 [57%] de 56, China 54 [32%] de 171, Oriente Medio 27 [30%] de 91, EE.UU. diez [2%] de 527; p<0-0001). KPC-2 (n=103 [49%]) y VIM-2 (n=75 [36%]) fueron las carbapenemasas más comunes en 211 aislados productores de carbapenemasas. Después de excluir a los pacientes de EE.UU., porque pocos aislados de EE.UU. tenían carbapenemasas, los pacientes con infecciones por CRPA productoras de carbapenemasas tuvieron una mayor mortalidad a los 30 días que aquellos con infecciones por CRPA no productoras de carbapenemasas, tanto en los análisis no ajustados (26 [22%] de 120 frente a 19 [12%] de 153; diferencia 9%, IC 95% 3-16) como ajustados (diferencia 7%, IC 95% 1-14). Interpretación: La aparición de diferentes carbapenemasas entre los aislados de CRPA en diferentes regiones geográficas y el aumento de la mortalidad asociada a las infecciones por CRPA productores de carbapenemasas ponen de manifiesto los retos terapéuticos que plantean estos organismos. Financiación: Institutos Nacionales de Salud.Background: Carbapenem-resistant Pseudomonas aeruginosa (CRPA) is a global threat, but the distribution and clinical significance of carbapenemases are unclear. The aim of this study was to define characteristics and outcomes of CRPA infections and the global frequency and clinical impact of carbapenemases harboured by CRPA. Methods: We conducted an observational, prospective cohort study of CRPA isolated from bloodstream, respiratory, urine, or wound cultures of patients at 44 hospitals (10 countries) between Dec 1, 2018, and Nov 30, 2019. Clinical data were abstracted from health records and CRPA isolates were whole-genome sequenced. The primary outcome was 30-day mortality from the day the index culture was collected. We compared outcomes of patients with CRPA infections by infection type and across geographic regions and performed an inverse probability weighted analysis to assess the association between carbapenemase production and 30-day mortality. Findings: We enrolled 972 patients (USA n=527, China n=171, south and central America n=127, Middle East n=91, Australia and Singapore n=56), of whom 581 (60%) had CRPA infections. 30-day mortality differed by infection type (bloodstream 21 [30%] of 69, respiratory 69 [19%] of 358, wound nine [14%] of 66, urine six [7%] of 88; p=0·0012) and geographical region (Middle East 15 [29%] of 52, south and central America 20 [27%] of 73, USA 60 [19%] of 308, Australia and Singapore three [11%] of 28, China seven [6%] of 120; p=0·0002). Prevalence of carbapenemase genes among CRPA isolates also varied by region (south and central America 88 [69%] of 127, Australia and Singapore 32 [57%] of 56, China 54 [32%] of 171, Middle East 27 [30%] of 91, USA ten [2%] of 527; p<0·0001). KPC-2 (n=103 [49%]) and VIM-2 (n=75 [36%]) were the most common carbapenemases in 211 carbapenemase-producing isolates. After excluding USA patients, because few US isolates had carbapenemases, patients with carbapenemase-producing CRPA infections had higher 30-day mortality than those with non-carbapenemase-producing CRPA infections in both unadjusted (26 [22%] of 120 vs 19 [12%] of 153; difference 9%, 95% CI 3–16) and adjusted (difference 7%, 95% CI 1–14) analyses. Interpretation: The emergence of different carbapenemases among CRPA isolates in different geographical regions and the increased mortality associated with carbapenemase-producing CRPA infections highlight the therapeutic challenges posed by these organisms. Funding: National Institutes of Health

Ceruloplasmin Deficiency Reduces Levels of Iron and BDNF in the Cortex and Striatum of Young Mice and Increases Their Vulnerability to Stroke

Ceruloplasmin (Cp) is an essential ferroxidase that plays important roles in cellular iron trafficking. Previous findings suggest that the proper regulation and subcellular localization of iron are very important in brain cell function and viability. Brain iron dyshomeostasis is observed during normal aging, as well as in several neurodegenerative disorders such as Alzheimer's, Parkinson's and Huntington's diseases, coincident with areas more susceptible to insults. Because of their high metabolic demand and electrical excitability, neurons are particularly vulnerable to ischemic injury and death. We therefore set out to look for abnormalities in the brain of young adult mice that lack Cp. We found that iron levels in the striatum and cerebral cortex of these young animals are significantly lower than wild-type (WT) controls. Also mRNA levels of the neurotrophin brain derived neurotrophic factor (BDNF), known for its role in maintenance of cell viability, were decreased in these brain areas. Chelator-mediated depletion of iron in cultured neural cells resulted in reduced BDNF expression by a posttranscriptional mechanism, suggesting a causal link between low brain iron levels and reduced BDNF expression. When the mice were subjected to middle cerebral artery occlusion, a model of focal ischemic stroke, we found increased brain damage in Cp-deficient mice compared to WT controls. Our data indicate that lack of Cp increases neuronal susceptibility to ischemic injury by a mechanism that may involve reduced levels of iron and BDNF

Multi-model seascape genomics identifies distinct environmental drivers of selection among sympatric marine species

Background As global change and anthropogenic pressures continue to increase, conservation and management increasingly needs to consider species’ potential to adapt to novel environmental conditions. Therefore, it is imperative to characterise the main selective forces acting on ecosystems, and how these may influence the evolutionary potential of populations and species. Using a multi-model seascape genomics approach, we compare putative environmental drivers of selection in three sympatric southern African marine invertebrates with contrasting ecology and life histories: Cape urchin (Parechinus angulosus), Common shore crab (Cyclograpsus punctatus), and Granular limpet (Scutellastra granularis). Results Using pooled (Pool-seq), restriction-site associated DNA sequencing (RAD-seq), and seven outlier detection methods, we characterise genomic variation between populations along a strong biogeographical gradient. Of the three species, only S. granularis showed significant isolation-by-distance, and isolation-by-environment driven by sea surface temperatures (SST). In contrast, sea surface salinity (SSS) and range in air temperature correlated more strongly with genomic variation in C. punctatus and P. angulosus. Differences were also found in genomic structuring between the three species, with outlier loci contributing to two clusters in the East and West Coasts for S. granularis and P. angulosus, but not for C. punctatus. Conclusion The findings illustrate distinct evolutionary potential across species, suggesting that species-specific habitat requirements and responses to environmental stresses may be better predictors of evolutionary patterns than the strong environmental gradients within the region. We also found large discrepancies between outlier detection methodologies, and thus offer a novel multi-model approach to identifying the principal environmental selection forces acting on species. Overall, this work highlights how adding a comparative approach to seascape genomics (both with multiple models and species) can elucidate the intricate evolutionary responses of ecosystems to global change

The genomic landscape of balanced cytogenetic abnormalities associated with human congenital anomalies

Despite the clinical significance of balanced chromosomal abnormalities (BCAs), their characterization has largely been restricted to cytogenetic resolution. We explored the landscape of BCAs at nucleotide resolution in 273 subjects with a spectrum of congenital anomalies. Whole-genome sequencing revised 93% of karyotypes and demonstrated complexity that was cryptic to karyotyping in 21% of BCAs, highlighting the limitations of conventional cytogenetic approaches. At least 33.9% of BCAs resulted in gene disruption that likely contributed to the developmental phenotype, 5.2% were associated with pathogenic genomic imbalances, and 7.3% disrupted topologically associated domains (TADs) encompassing known syndromic loci. Remarkably, BCA breakpoints in eight subjects altered a single TAD encompassing MEF2C, a known driver of 5q14.3 microdeletion syndrome, resulting in decreased MEF2C expression. We propose that sequence-level resolution dramatically improves prediction of clinical outcomes for balanced rearrangements and provides insight into new pathogenic mechanisms, such as altered regulation due to changes in chromosome topology

Global, regional, and national incidence, prevalence, and years lived with disability for 328 diseases and injuries for 195 countries, 1990–2016: a systematic analysis for the Global Burden of Disease Study 2016

As mortality rates decline, life expectancy increases, and populations age, non-fatal outcomes of diseases and injuries are becoming a larger component of the global burden of disease. The Global Burden of Diseases, Injuries, and Risk Factors Study 2016 (GBD 2016) provides a comprehensive assessment of prevalence, incidence, and years lived with disability (YLDs) for 328 causes in 195 countries and territories from 1990 to 2016