Search

Objectives: To evaluate the efficacy of a single direct epidural injection of tumor necrosis factor (TNF)-α inhibitor to reduce the pathological changes in nerve fiber injuries in a rat model of chronic compression of the dorsal root ganglion (CCD).Materials and Methods: A surgical procedure for CCD was performed in 17 adult male F-344 rats. The effects of the epidural TNF-α inhibitors on CCD-induced pathological changes were investigated. Three groups of rats (n = 17) were used: (1) CCD + saline (n = 4), (2) CCD + triamcinolone (n = 5), and (3) CCD + TNF-α inhibitors (n = 8). Their dorsal root ganglia and nerve roots were removed on postoperative day 14. The intraneural edema, demyelination, and Wallerian degeneration of all 17 rats were scored pathologically.Results: The pathology scores of the rats in the TNF-α inhibitor treatment group (1.38 ± 0.74) indicated a mild degree of intraneural edema compared to the saline treatment group (2.25 ± 0.50, p = 0.041). In addition, rats in the TNF-α inhibitor treatment group (2.13 ± 0.35) had a mild degree of demyelination compared to the saline treatment group (2.75 ± 0.50, p = 0.038) and the triamcinolone treatment group (2.80 ± 0.45, p = 0.019). The differences in the pathology scores for Wallerian degeneration were not statistically significant in all three study groups (p = 0.658).Conclusion: The epidural injection of a TNF-α inhibitor was more effective than a placebo and comparable to triamcinolone in reducing pathological nerve injury progression. [ABSTRACT FROM AUTHOR]

Control of hospital-associated Enterococcus faecium infection is a strenuous task due to the difficulty of identifying transmission routes and the persistence of this nosocomial pathogen despite the implementation of infection control measures that have been successful with other important nosocomial pathogens. This study provides a comprehensive analysis of over 100 E. faecium isolates collected from 66 cancer patients at the University of Arkansas for Medical Sciences (UAMS) between June 2018 and May 2019. In the top-down approach used in this study, we employed, in addition to the 106 E. faecium UAMS isolates, a filtered set of 2,167 E. faecium strains from the GenBank database to assess the current population structure of E. faecium species and, consequently, to identify the lineages associated with our clinical isolates. We then evaluated the antibiotic resistance and virulence profiles of hospital-associated strains from the species pool, focusing on antibiotics of last resort, to establish an updated classification of high-risk and multidrug-resistant nosocomial clones. Further investigation of the clinical isolates collected from UAMS patients using whole-genome sequencing analytical methodologies (core genome multilocus sequence typing [cgMLST], core single nucleotide polymorphism [coreSNP] analysis, and phylogenomics), with the addition of patient epidemiological data, revealed a polyclonal outbreak of three sequence types occurring simultaneously in different patient wards. The integration of genomic and epidemiological data collected from the patients increased our understanding of the relationships and transmission dynamics of the E. faecium isolates. Our study provides new insights into genomic surveillance of E. faecium to assist in monitoring and further limiting the spread of multidrug-resistant E. faecium. IMPORTANCE Enterococcus faecium is a member of the gastrointestinal microbiota. Although its virulence is low in healthy, immunocompetent individuals, E. faecium has become the third leading cause of health care-associated infections in the United States. This study provides a comprehensive analysis of over 100 E. faecium isolates collected from cancer patients at the University of Arkansas for Medical Sciences (UAMS). We employed a top-down analytical approach (from population genomics to molecular biology) to classify our clinical isolates into their genetic lineages and thoroughly evaluate their antibiotic resistance and virulence profiles. The addition of patient epidemiological data to the whole-genome sequencing analytical methodologies performed in the study allowed us to increase our understanding of the relationships and transmission dynamics of the E. faecium isolates. This study provides new insights into genomic surveillance of E. faecium to help monitor and further limit the spread of multidrug-resistant E. faecium., Competing Interests: The authors declare no conflict of interest.

Background: The relationships of patient factors and caregiver contribution to patients' self-care to different types of self-care have been rarely examined in Korean patients with heart failure., Objective: The aim of this study was to examine patient (ie, age, depressive symptoms, and self-care confidence) and caregiver (ie, caregiver contribution to self-care maintenance and self-care management, and caregiver confidence in contributing to self-care) factors related to different types of self-care (ie, self-care maintenance, symptom perception, and self-care management) and self-care confidence in Korean patients with heart failure., Methods: In this cross-sectional, correlational study, data from 41 pairs of patients (mean age, 68.0 years) and caregivers (mean age, 54.1 years) were collected and analyzed using multiple regression., Results: Higher levels of self-care confidence in patients were related to higher levels of self-care maintenance in patients. Higher levels of self-care confidence in patients were related to higher levels of symptom perception. Higher levels of self-care confidence in patients and caregiver contribution to self-care maintenance were related to higher levels of self-care management. Less severe depressive symptoms in patients and higher levels of caregiver confidence in contributing to self-care were related to higher levels of self-care confidence., Conclusion: Different patient and caregiver factors were related to different types of self-care and self-care confidence in Korean patients, but patients' self-care confidence was related to all types of self-care. Clinicians and researchers need to develop and deliver effective interventions to both patients and their caregivers to improve patients' self-care confidence and, in turn, self-care, considering different factors associated with each type of self-care., Competing Interests: The authors have no conflicts of interest to disclose., (Copyright © 2022 Wolters Kluwer Health, Inc. All rights reserved.)

Here, we report the complete genome sequence of Providencia rettgeri isolate PROV_UAMS_01, which was recovered in 2021 from a urine sample from a hospitalized patient in Arkansas, USA. The genome sequence of P. rettgeri isolate PROV_UAMS_01 comprises a single chromosomal replicon with a G+C content of 40.51% and a total of 3,887 genes.

Insertion sequences (ISs) and other transposable elements are associated with the mobilization of antibiotic resistance determinants and the modulation of pathogenic characteristics. In this work, we aimed to investigate the association between ISs and antibiotic resistance genes, and their role in the dissemination and modification of the antibiotic-resistant phenotype. To that end, we leveraged fully resolved Enterococcus faecium and Enterococcus faecalis genomes of isolates collected over 5 days from an inpatient with prolonged bacteraemia. Isolates from both species harboured similar IS family content but showed significant species-dependent differences in copy number and arrangements of ISs throughout their replicons. Here, we describe two inter-specific IS-mediated recombination events and IS-mediated excision events in plasmids of E. faecium isolates. We also characterize a novel arrangement of the ISs in a Tn1546-like transposon in E. faecalis isolates likely implicated in a vancomycin genotype-phenotype discrepancy. Furthermore, an extended analysis revealed a novel association between daptomycin resistance mutations in liaSR genes and a putative composite transposon in E. faecium , offering a new paradigm for the study of daptomycin resistance and novel insights into its dissemination. In conclusion, our study highlights the role ISs and other transposable elements play in the rapid adaptation and response to clinically relevant stresses such as aggressive antibiotic treatment in enterococci.

Vegetables may serve as a reservoir for antibiotic resistant bacteria and resistance genes. AmpC β-lactamases and extended spectrum beta-lactamases (ESBL) inactivate commonly used β-lactam antibiotics, including penicillins and cephalosporins. In this study, we determined the prevalence of AmpC and ESBL-producing Enterobacterales in retail vegetables in the United States. A total of 88 vegetable samples were collected for the screening of AmpC and ESBL-producing Enterobacterales using CHROMagar ESBL agar. These vegetables included washed ready-to-eat salad (23), microgreens/sprouts (13), lettuce (11), herbs (11), spinach (5), mushrooms (5), brussels sprouts (4), kale (3), and other vegetable samples (13). AmpC and ESBL activity in these isolates were determined using double disk combination tests. Two vegetable samples (2.27%), organic basil and brussels sprouts, were positive for AmpC-producing Enterobacterales and eight samples (9.09%), including bean sprouts, organic parsley, organic baby spinach, and several mixed salads, were positive for ESBL-producing Enterobacterales. Whole genome sequencing was used to identify the bacterial species and resistance genes in these isolates. Genes encoding AmpC β-lactamases were found in Enterobacter hormaechei strains S43-1 and 74-2, which were consistent with AmpC production phenotypes. Multidrug-resistant E. hormaechei strains S11-1, S17-1, and S45-4 possess an ESBL gene, bla SHV66 , whereas five Serratia fonticola isolates contain genes encoding a minor ESBL, FONA-5. In addition, we used shotgun metagenomic sequencing approach to examine the microbiome and resistome profiles of three spinach samples. We found that Pseudomonas was the most prevalent bacteria genus in the spinach samples. Within the Enterobacteriaceae family, Enterobacter was the most abundant genus in the spinach samples. Moreover, antibiotic resistance genes encoding 12 major classes of antibiotics, including β-lactam antibiotics, aminoglycoside, macrolide, fluoroquinolone, and others, were found in these spinach samples. Therefore, vegetables can serve as an important vehicle for transmitting antibiotic resistance. The study highlights the need for antibiotic resistance surveillance in vegetable products., Competing Interests: Declarations of interest: none

In this chapter, we leverage a novel approach to assess the seamless population structure of Pseudomonas aeruginosa, using the full repertoire of genomes sequenced to date (GenBank, April 6, 2020). In order to assess the set of core functions that represents the species as well as the differences in these core functions among the phylogroups observed in the population structure analysis, we performed pangenome analyses at the species level and at the phylogroup level. The existence of the phylogroups described in the population structure analyses was supported by their different profiles of antibiotic-resistant determinants. Finally, we utilized a presence/absence matrix of protein families from the entire species to evaluate if P. aeruginosa phylogroups can be differentiated according to their accessory genomic content. Our analysis shows that the core genome of P. aeruginosa is approximately 62% of the average gene content for the species, and it is highly enriched with pathways related to the metabolism of carbohydrates and amino acids as well as cellular processes and cell maintenance. The analysis of the accessory genome of P. aeruginosa performed in this chapter confirmed not only the existence of the three phylogroups previously described in the population structure analysis, but also of 29 genetic substructures (subgroups) within the main phylogroups. Our work illustrates the utility of populations genomics pipelines to better understand highly complex bacterial species such as P. aeruginosa., (© 2022. The Author(s), under exclusive license to Springer Nature Switzerland AG.)

The acute radiation syndrome is defined in large part by radiation injury in the hematopoietic and gastrointestinal (GI) systems. To identify new pathways involved in radiation-induced GI injury, this study assessed dose- and time-dependent changes in plasma metabolites in a nonhuman primate model of whole abdominal irradiation. Male and female adult Rhesus monkeys were exposed to 6 MV photons to the abdomen at doses ranging between 8 and 14 Gy. At time points from 1 to 60 days after irradiation, plasma samples were collected and subjected to untargeted metabolomics. With the limited sample size of females, different discovery times after irradiation between males and females were observed in metabolomics pattern. Detailed analyses are restricted to only males for the discovery power. Radiation caused an increase in fatty acid oxidation and circulating levels of corticosteroids which may be an indication of physiological stress, and amino acids, indicative of a cellular repair response. The largest changes were observed at days 9 and 10 post-irradiation, with most returning to baseline at day 30. In addition, dysregulated metabolites involved in amino acid pathways, which might indicate changes in the microbiome, were detected. In conclusion, abdominal irradiation in a nonhuman primate model caused a plasma metabolome profile indicative of GI injury. These results point to pathways that may be targeted for intervention or used as early indicators of GI radiation injury. Moreover, our results suggest that effects are sex-specific and that interventions may need to be tailored accordingly.

Age-associated mitochondrial dysfunction and oxidative damage are primary causes for multiple health problems including sarcopenia and cardiovascular disease (CVD). Though the role of Nrf2, a transcription factor that regulates cytoprotective gene expression, in myopathy remains poorly defined, it has shown beneficial properties in both sarcopenia and CVD. Sulforaphane (SFN), a natural compound Nrf2-related activator of cytoprotective genes, provides protection in several disease states including CVD and is in various stages of clinical trials, from cancer prevention to reducing insulin resistance. This study aimed to determine whether SFN may prevent age-related loss of function in the heart and skeletal muscle. Cohorts of 2-month-old and 21- to 22-month-old mice were administered regular rodent diet or diet supplemented with SFN for 12 weeks. At the completion of the study, skeletal muscle and heart function, mitochondrial function, and Nrf2 activity were measured. Our studies revealed a significant drop in Nrf2 activity and mitochondrial functions, together with a loss of skeletal muscle and cardiac function in the old control mice compared to the younger age group. In the old mice, SFN restored Nrf2 activity, mitochondrial function, cardiac function, exercise capacity, glucose tolerance, and activation/differentiation of skeletal muscle satellite cells. Our results suggest that the age-associated decline in Nrf2 signaling activity and the associated mitochondrial dysfunction might be implicated in the development of age-related disease processes. Therefore, the restoration of Nrf2 activity and endogenous cytoprotective mechanisms by SFN may be a safe and effective strategy to protect against muscle and heart dysfunction due to aging., (© 2020 The Authors. Aging Cell published by the Anatomical Society and John Wiley & Sons Ltd.)

Dietary factors modulate interactions between the microbiome, metabolome, and immune system. Sulforaphane (SFN) exerts effects on aging, cancer prevention and reducing insulin resistance. This study investigated effects of SFN on the gut microbiome and metabolome in old mouse model compared with young mice. Young (6-8 weeks) and old (21-22 months) male C57BL/6J mice were provided regular rodent chow ± SFN for 2 months. We collected fecal samples before and after SFN administration and profiled the microbiome and metabolome. Multi-omics datasets were analyzed individually and integrated to investigate the relationship between SFN diet, the gut microbiome, and metabolome. The SFN diet restored the gut microbiome in old mice to mimic that in young mice, enriching bacteria known to be associated with an improved intestinal barrier function and the production of anti-inflammatory compounds. The tricarboxylic acid cycle decreased and amino acid metabolism-related pathways increased. Integration of multi-omic datasets revealed SFN diet-induced metabolite biomarkers in old mice associated principally with the genera, Oscillospira , Ruminococcus , and Allobaculum . Collectively, our results support a hypothesis that SFN diet exerts anti-aging effects in part by influencing the gut microbiome and metabolome. Modulating the gut microbiome by SFN may have the potential to promote healthier aging.

Genomic DNA is the best "unique identifier" for organisms. Alignment-free phylogenomic analysis, simple, fast, and efficient method to compare genome sequences, relies on looking at the distribution of small DNA sequence of a particular length, referred to as k-mer. The k-mer approach has been explored as a basis for sequence analysis applications, including assembly, phylogenetic tree inference, and classification. Although this approach is not novel, selecting the appropriate k-mer length to obtain the optimal resolution is rather arbitrary. However, it is a very important parameter for achieving the appropriate resolution for genome/sequence distances to infer biologically meaningful phylogenetic relationships. Thus, there is a need for a systematic approach to identify the appropriate k-mer from whole-genome sequences. We present K-mer-length Iterative Selection for UNbiased Ecophylogenomics (KITSUNE), a tool for assessing the empirically optimal k-mer length of any given set of genomes of interest for phylogenomic analysis via a three-step approach based on (1) cumulative relative entropy (CRE), (2) average number of common features (ACF), and (3) observed common features (OCF). Using KITSUNE, we demonstrated the feasibility and reliability of these measurements to obtain empirically optimal k-mer lengths of 11, 17, and ∼34 from large genome datasets of viruses, bacteria, and fungi, respectively. Moreover, we demonstrated a feature of KITSUNE for accurate species identification for the two de novo assembled bacterial genomes derived from error-prone long-reads sequences, and for a published yeast genome. In addition, KITSUNE was used to identify the shortest species-specific k-mer accurately identifying viruses. KITSUNE is freely available at https://github.com/natapol/kitsune., (Copyright © 2020 Pornputtapong, Acheampong, Patumcharoenpol, Jenjaroenpun, Wongsurawat, Jun, Yongkiettrakul, Chokesajjawatee and Nookaew.)

Due to animal welfare issues, European Union has banned the use of conventional cages (CC) and non-EU countries including the US are also under constant public pressure to restrict their use in egg production. Very limited information is available on the composition of the microbial community of hens raised in different housing environments. This study was conducted to determine the effects of CC and enriched colony cages (EC) on cecal microbiota of two commercial laying hen strains, Hy-Line W36 (W36) and Hy-Line Brown (HB) during the late production stage (53, 58, 67, and 72 weeks of age). Cecal microbiota was studied by analyzing 16S rRNA gene sequences with Quantitative Insights Into Microbial Ecology (QIIME) 2 ver. 2018.8. Differentially abundant taxa were identified by Linear discriminant analysis Effect Size (LEfSe) analysis ( P < 0.05, LDA score > 2.0). At phylum level, Actinobacteria was significantly enriched in W36 at all time points while Synergistetes (53 weeks), Spirochaetes (58 weeks), and Synergistetes and Spirochaetes (67 weeks) were significantly higher in HB. At genus level, Bifidobacterium (at all time points) and butyric acid producing genera such as Butyricicoccus and Subdoligranulum (58 and 72 weeks) were significantly higher in W36 as compared to HB. Moreover, Proteobacteria (72 weeks) and its associated genus Campylobacter (67 and 72 weeks) were significantly enriched in EC as compared to CC. Alpha diversity was significantly higher in HB (at all time points) and in EC (67 weeks) as compared to W36 and CC, respectively. Similarly, there was a significant difference in community structure (beta diversity) between W36 and HB (all time points) as well as between EC and CC (67 weeks). The effect of housing and strains was not only seen at the bacterial composition and structure but also reflected at their functional level. Notably, KEGG metabolic pathways predicted to be involved in carbohydrates degradation and amino acids biosynthesis by PICRUSt analysis were significantly different between W36 and HB housed at CC and EC. In sum, cecal microbiota composition, diversities, and their functional pathways were affected by housing type which further varied between two commercial laying hen strains, HB and W36. This suggests that both housing and genetic strains of laying hens should be considered for selection of the alternative housing systems such as enriched colony cage., (Copyright © 2020 Adhikari, Jun, Kwon, Kiess and Adhikari.)

Here, we present complete genome sequences of four Enterococcus faecium isolates, obtained from two patients with apparent vancomycin-resistant Enterococcus faecium bacteremia; these isolates also carried two mutations known to be associated with daptomycin resistance. Sequences were obtained using de novo and hybrid assembly of Oxford Nanopore and Illumina sequence data., (Copyright © 2020 Jenjaroenpun et al.)

The intraspecies genomic diversity of the single-strand RNA (+) virus species hepatitis A virus (Hepatovirus), hepatitis C virus (Hepacivirus), and hepatitis E virus (Orthohepevirus) was compared. These viral species all can cause liver inflammation (hepatitis), but share no gene similarity. The codon usage of human hepatitis A virus (HAV) is suboptimal for replication in its host, a characteristic it shares with taxonomically related rodent, simian, and bat hepatitis A virus species. We found this codon usage to be strikingly similar to that of Triatoma virus that infects blood-sucking kissing bugs. The codon usage of that virus is well adapted to its insect host. The codon usage of HAV is also similar to other invertebrate viruses of various taxonomic families. An evolutionary ancestor of HAV and related virus species is hypothesized to be an insect virus that underwent a host jump to infect mammals. The similarity between HAV and invertebrate viruses goes beyond codon usage, as they also share amino acid composition characteristics, while not sharing direct sequence homology. In contrast, hepatitis C virus and hepatitis E virus are highly similar in codon usage preference, nucleotide composition, and amino acid composition, and share these characteristics with Human pegivirus A, West Nile virus, and Zika virus. We present evidence that these observations are only partly explained by differences in nucleotide composition of the complete viral codon regions. We consider the combination of nucleotide composition, amino acid composition, and codon usage preference suitable to provide information on possible evolutionary similarities between distant virus species that cannot be investigated by phylogeny., (© 2019 The Authors. MicrobiologyOpen published by John Wiley & Sons Ltd.)

In 2016, a year-long large-scale mumps outbreak occurred in Arkansas among a highly-vaccinated population. A total of 2954 mumps cases were identified during this outbreak. The majority of cases (1676 (57%)) were school-aged children (5-17 years), 1536 (92%) of these children had completed the mumps vaccination schedule. To weigh the possibility that the mumps virus evaded vaccine-induced immunity in the affected Arkansas population, we established a pipeline for genomic characterization of the outbreak strains. Our pipeline produces whole-genome sequences along with phylogenetic analysis of the outbreak mumps virus strains. We collected buccal swab samples of patients who tested positive for the mumps virus during the 2016 Arkansas outbreak, and used the portable Oxford Nanopore Technology to sequence the extracted strains. Our pipeline identified the genotype of the Arkansas mumps strains as genotype G and presented a genome-based phylogenetic tree with superior resolution to a standard small hydrophobic (SH) gene-based tree. We phylogenetically compared the Arkansas whole-genome sequences to all publicly available mumps strains. While these analyses show that the Arkansas mumps strains are evolutionarily distinct from the vaccine strains, we observed no correlation between vaccination history and phylogenetic grouping. Furthermore, we predicted potential B-cell epitopes encoded by the Arkansas mumps strains using a random forest prediction model trained on antibody-antigen protein structures. Over half of the predicted epitopes of the Jeryl-Lynn vaccine strains in the Hemagglutinin-Neuraminidase (HN) surface glycoprotein (a major target of neutralizing antibodies) region are missing in the Arkansas mumps strains. In-silico analyses of potential epitopes may indicate that the Arkansas mumps strains display antigens with reduced immunogenicity, which may contribute to reduced vaccine effectiveness. However, our in-silico findings should be assessed by robust experiments such as cross neutralization assays. Metadata analysis showed that vaccination history had no effect on the evolution of the Arkansas mumps strains during this outbreak. We conclude that the driving force behind the spread of the mumps virus in the 2016 Arkansas outbreak remains undetermined., (Copyright © 2019 Elsevier B.V. All rights reserved.)

Vancomycin-resistant Enterococcus faecium (VREfm) is a major cause of nosocomial infections of the bloodstream and urinary tract. Here, we report the draft genome sequences of 48 vancomycin-resistant E. faecium isolates recovered from inpatients exhibiting clinical signs of bacteremia at the University of Arkansas for Medical Sciences (UAMS)., (Copyright © 2019 Udaondo et al.)

Background: Here we investigated the influence of different stabilization and storage strategies on the quality and composition of the fecal microbial community. Namely, same-day isolated murine DNA was compared to samples stored for 1 month in air at ambient temperature, with or without preservative buffers (i.e. EDTA and lysis buffer), different temperatures (i.e. 4 °C, - 20 °C, and - 80 °C), and hypoxic conditions., Results: Only storage in lysis buffer significantly reduced DNA content, yet without integrity loss. Storage in EDTA affected alpha diversity the most, which was also reflected in cluster separation. Distinct changes were also seen in the phyla and bacterial species abundance per storage strategy. Metabolic function analysis showed 22 pathways not significantly affected by storage conditions, whereas the tyrosine metabolism pathway was significantly changed in all strategies except by EDTA., Conclusion: Each long-term storage strategy introduced a unique post-collection bias, which is important to take into account when interpreting data.

Infections due to Clostridioides difficile (previously known as Clostridium difficile) are a major problem in hospitals, where cases can be caused by community-acquired strains as well as by nosocomial spread. Whole genome sequences from clinical samples contain a lot of information but that needs to be analyzed and compared in such a way that the outcome is useful for clinicians or epidemiologists. Here, we compare 663 public available complete genome sequences of C. difficile using average amino acid identity (AAI) scores. This analysis revealed that most of these genomes (640, 96.5%) clearly belong to the same species, while the remaining 23 genomes produce four distinct clusters within the Clostridioides genus. The main C. difficile cluster can be further divided into sub-clusters, depending on the chosen cutoff. We demonstrate that MLST, either based on partial or full gene-length, results in biased estimates of genetic differences and does not capture the true degree of similarity or differences of complete genomes. Presence of genes coding for C. difficile toxins A and B (ToxA/B), as well as the binary C. difficile toxin (CDT), was deduced from their unique PfamA domain architectures. Out of the 663 C. difficile genomes, 535 (80.7%) contained at least one copy of ToxA or ToxB, while these genes were missing from 128 genomes. Although some clusters were enriched for toxin presence, these genes are variably present in a given genetic background. The CDT genes were found in 191 genomes, which were restricted to a few clusters only, and only one cluster lacked the toxin A/B genes consistently. A total of 310 genomes contained ToxA/B without CDT (47%). Further, published metagenomic data from stools were used to assess the presence of C. difficile sequences in blinded cases of C. difficile infection (CDI) and controls, to test if metagenomic analysis is sensitive enough to detect the pathogen, and to establish strain relationships between cases from the same hospital. We conclude that metagenomics can contribute to the identification of CDI and can assist in characterization of the most probable causative strain in CDI patients.

Adverse growth conditions can lead to decreased plant growth, productivity, and survival, resulting in poor yields or failure of crops and biofeedstocks. In some cases, the microbial community associated with plants has been shown to alleviate plant stress and increase plant growth under suboptimal growing conditions. A systematic understanding of how the microbial community changes under these conditions is required to understand the contribution of the microbiome to water utilization, nutrient uptake, and ultimately yield. Using a microbiome inoculation strategy, we studied how the belowground microbiome of Populus deltoides changes in response to diverse environmental conditions, including water limitation, light limitation (shading), and metal toxicity. While plant responses to treatments in terms of growth, photosynthesis, gene expression and metabolite profiles were varied, we identified a core set of bacterial genera that change in abundance in response to host stress. The results of this study indicate substantial structure in the plant microbiome community and identify potential drivers of the phytobiome response to stress. IMPORTANCE The identification of a common "stress microbiome" indicates tightly controlled relationships between the plant host and bacterial associates and a conserved structure in bacterial communities associated with poplar trees under different growth conditions. The ability of the microbiome to buffer the plant from extreme environmental conditions coupled with the conserved stress microbiome observed in this study suggests an opportunity for future efforts aimed at predictably modulating the microbiome to optimize plant growth.

Background: Zika virus (ZIKV) is an emerging human pathogen. Since its arrival in the Western hemisphere, from Africa via Asia, it has become a serious threat to pregnant women, causing microcephaly and other neuropathies in developing fetuses. The mechanisms behind these teratogenic effects are unknown, although epidemiological evidence suggests that microcephaly is not associated with the original, African lineage of ZIKV. The sequences of 196 published ZIKV genomes were used to assess whether recently proposed mechanistic explanations for microcephaly are supported by molecular level changes that may have increased its virulence since the virus left Africa. For this we performed phylogenetic, recombination, adaptive evolution and tetramer frequency analyses, and compared protein sequences for the presence of protease cleavage sites, Pfam domains, glycosylation sites, signal peptides, trans-membrane protein domains, and phosphorylation sites., Results: Recombination events within or between Asian and Brazilian lineages were not observed, and likewise there were no differences in protease cleavage, glycosylation sites, signal peptides or trans-membrane domains between African and Brazilian strains. The frequency of Retinoic Acid Response Element (RARE) sequences was increased in Brazilian strains. Genetic adaptation was also apparent by tetramer signatures that had undergone major changes in the past but has stabilized in the Brazilian lineage despite subsequent geographic spread, suggesting the viral population presently propagates in the same host species in various regions. Evidence for selection pressure was recognized for several amino acid sites in the Brazilian lineage compared to the African lineage, mainly in nonstructural proteins, especially protein NS4B. A number of these positively selected mutations resulted in an increased potential to be phosphorylated in the Brazilian lineage compared to the African linage, which may have increased their potential to interfere with neural fetal development., Conclusions: ZIKV seems to have adapted to a limited number of hosts, including humans, during which its virulence increased. Its protein NS4B, together with NS4A, has recently been shown to inhibit Akt-mTOR signaling in human fetal neural stem cells, a key pathway for brain development. We hypothesize that positive selection of novel phosphorylation sites in the protein NS4B of the Brazilian lineage could interfere with phosphorylation of Akt and mTOR, impairing Akt-mTOR signaling and this may result in an increased risk for developmental neuropathies.

We present here the complete genome sequences of three mumps virus (MuV) strains isolated from patients who tested positive for the mumps virus during a mumps outbreak in Springdale, AR (USA), in 2016. The virus genomes, sequenced with Oxford Nanopore technology, belong to genotype G and have an average length of 15,342 nucleotides (nt)., (Copyright © 2017 Alkam et al.)

Background: Functional annotation of bacterial genomes is an obligatory and crucially important step of information processing from the genome sequences into cellular mechanisms. However, there is a lack of computational methods to evaluate the quality of functional assignments., Results: We developed a genome-scale model that assigns Bayesian probability to each gene utilizing a known property of functional similarity between neighboring genes in bacteria., Conclusions: Our model clearly distinguished true annotation from random annotation with Bayesian annotation probability >0.95. Our model will provide a useful guide to quantitatively evaluate functional annotation methods and to detect gene sets with reliable annotations.

The development of rapid, economical genome sequencing has shed new light on the classification of viruses. As of October 2016, the National Center for Biotechnology Information (NCBI) database contained >2 million viral genome sequences and a reference set of ~4000 viral genome sequences that cover a wide range of known viral families. Whole-genome sequences can be used to improve viral classification and provide insight into the viral "tree of life". However, due to the lack of evolutionary conservation amongst diverse viruses, it is not feasible to build a viral tree of life using traditional phylogenetic methods based on conserved proteins. In this study, we used an alignment-free method that uses k-mers as genomic features for a large-scale comparison of complete viral genomes available in RefSeq. To determine the optimal feature length, k (an essential step in constructing a meaningful dendrogram), we designed a comprehensive strategy that combines three approaches: (1) cumulative relative entropy, (2) average number of common features among genomes, and (3) the Shannon diversity index. This strategy was used to determine k for all 3,905 complete viral genomes in RefSeq. The resulting dendrogram shows consistency with the viral taxonomy of the ICTV and the Baltimore classification of viruses.

It has been 30 years since the initial emergence and subsequent rapid global spread of multidrug-resistant Salmonella entericaserovar Typhimurium DT104 (MDR DT104). Nonetheless, its origin and transmission route have never been revealed. We used whole-genome sequencing (WGS) and temporally structured sequence analysis within a Bayesian framework to reconstruct temporal and spatial phylogenetic trees and estimate the rates of mutation and divergence times of 315S Typhimurium DT104 isolates sampled from 1969 to 2012 from 21 countries on six continents. DT104 was estimated to have emerged initially as antimicrobial susceptible in ∼1948 (95% credible interval [CI], 1934 to 1962) and later became MDR DT104 in ∼1972 (95% CI, 1972 to 1988) through horizontal transfer of the 13-kb Salmonella genomic island 1 (SGI1) MDR region into susceptible strains already containing SGI1. This was followed by multiple transmission events, initially from central Europe and later between several European countries. An independent transmission to the United States and another to Japan occurred, and from there MDR DT104 was probably transmitted to Taiwan and Canada. An independent acquisition of resistance genes took place in Thailand in ∼1975 (95% CI, 1975 to 1990). In Denmark, WGS analysis provided evidence for transmission of the organism between herds of animals. Interestingly, the demographic history of Danish MDR DT104 provided evidence for the success of the program to eradicate Salmonellafrom pig herds in Denmark from 1996 to 2000. The results from this study refute several hypotheses on the evolution of DT104 and suggest that WGS may be useful in monitoring emerging clones and devising strategies for prevention of Salmonella infections., (Copyright © 2016 Leekitcharoenphon et al.)

Purpose: To evaluate the feasibility of using an internally cooled directional electrode (ICDE) to create sufficient directional ablation in the desired area and to investigate whether use of the ICDE reduces thermal injury to the adjacent organs., Materials and Methods: The animal care and use committee approved this animal study. In ex vivo studies of bovine livers and colons, temperature was measured in both ablation and nonablation directions with the ICDE, temperature at the liver-colon interface and the presence of burns on the colon were evaluated with each ICDE and internally cooled conventional electrode (ICCE), and the ablation area with the use of three ICDEs in the switching multichannel mode was evaluated. In an in vivo study of 10 pigs, 11 ablations were performed by using three ICDEs in the switching multichannel mode. In addition, the difference in thermal injury to the gallbladder and stomach was compared between use of the ICDE and the ICCE. Mann-Whitney test and the Fisher exact test were used for statistical analysis., Results: In the ex vivo study, the hepatic temperature in the nonablation direction was lower than 50°C; temperature higher than 47°C was maintained for less than 3 minutes, even 0.5 cm away from the ICDE; and the temperature at the liver-colon interface did not increase above 50°C. In the switching multichannel mode, a fused coagulation zone was created. In the in vivo study, confluent ablation was created in the inner circle of three ICDEs with a mean interelectrode distance of less than 2.6 cm. Use of the ICDE did not create a moderate-to-severe gallbladder wall injury in any of the pigs, but use of the ICCE caused injury in three of the five pigs (P > .05). The ICDE caused moderate-to-severe gastric wall injury at one of eight ablations, and the ICCE caused injury at eight of nine ablations (P < .05)., Conclusion: Use of an ICDE can create a sufficient ablation in the desired direction and can decrease thermal injury to the adjacent organs., (© RSNA, 2015.)

The Pseudomonas genus contains a metabolically versatile group of organisms that are known to occupy numerous ecological niches, including the rhizosphere and endosphere of many plants. Their diversity influences the phylogenetic diversity and heterogeneity of these communities. On the basis of average amino acid identity, comparative genome analysis of >1,000 Pseudomonas genomes, including 21 Pseudomonas strains isolated from the roots of native Populus deltoides (eastern cottonwood) trees resulted in consistent and robust genomic clusters with phylogenetic homogeneity. All Pseudomonas aeruginosa genomes clustered together, and these were clearly distinct from other Pseudomonas species groups on the basis of pangenome and core genome analyses. In contrast, the genomes of Pseudomonas fluorescens were organized into 20 distinct genomic clusters, representing enormous diversity and heterogeneity. Most of our 21 Populus-associated isolates formed three distinct subgroups within the major P. fluorescens group, supported by pathway profile analysis, while two isolates were more closely related to Pseudomonas chlororaphis and Pseudomonas putida. Genes specific to Populus-associated subgroups were identified. Genes specific to subgroup 1 include several sensory systems that act in two-component signal transduction, a TonB-dependent receptor, and a phosphorelay sensor. Genes specific to subgroup 2 contain hypothetical genes, and genes specific to subgroup 3 were annotated with hydrolase activity. This study justifies the need to sequence multiple isolates, especially from P. fluorescens, which displays the most genetic variation, in order to study functional capabilities from a pangenomic perspective. This information will prove useful when choosing Pseudomonas strains for use to promote growth and increase disease resistance in plants., (Copyright © 2015 Jun et al.)

The bacterial microbiota of plants is diverse, with 1000s of operational taxonomic units (OTUs) associated with any individual plant. In this work, we used phenotypic analysis, comparative genomics, and metabolic models to investigate the differences between 19 sequenced Pseudomonas fluorescens strains. These isolates represent a single OTU and were collected from the rhizosphere and endosphere of Populus deltoides. While no traits were exclusive to either endosphere or rhizosphere P. fluorescens isolates, multiple pathways relevant for plant-bacterial interactions are enriched in endosphere isolate genomes. Further, growth phenotypes such as phosphate solubilization, protease activity, denitrification and root growth promotion are biased toward endosphere isolates. Endosphere isolates have significantly more metabolic pathways for plant signaling compounds and an increased metabolic range that includes utilization of energy rich nucleotides and sugars, consistent with endosphere colonization. Rhizosphere P. fluorescens have fewer pathways representative of plant-bacterial interactions but show metabolic bias toward chemical substrates often found in root exudates. This work reveals the diverse functions that may contribute to colonization of the endosphere by bacteria and are enriched among closely related isolates.

Background: For decades there has been increasing interest in understanding the relationships between microbial communities and ecosystem functions. Current DNA sequencing technologies allows for the exploration of microbial communities in two principle ways: targeted rRNA gene surveys and shotgun metagenomics. For large study designs, it is often still prohibitively expensive to sequence metagenomes at both the breadth and depth necessary to statistically capture the true functional diversity of a community. Although rRNA gene surveys provide no direct evidence of function, they do provide a reasonable estimation of microbial diversity, while being a very cost-effective way to screen samples of interest for later shotgun metagenomic analyses. However, there is a great deal of 16S rRNA gene survey data currently available from diverse environments, and thus a need for tools to infer functional composition of environmental samples based on 16S rRNA gene survey data., Results: We present a computational method called pangenome-based functional profiles (PanFP), which infers functional profiles of microbial communities from 16S rRNA gene survey data for Bacteria and Archaea. PanFP is based on pangenome reconstruction of a 16S rRNA gene operational taxonomic unit (OTU) from known genes and genomes pooled from the OTU's taxonomic lineage. From this lineage, we derive an OTU functional profile by weighting a pangenome's functional profile with the OTUs abundance observed in a given sample. We validated our method by comparing PanFP to the functional profiles obtained from the direct shotgun metagenomic measurement of 65 diverse communities via Spearman correlation coefficients. These correlations improved with increasing sequencing depth, within the range of 0.8-0.9 for the most deeply sequenced Human Microbiome Project mock community samples. PanFP is very similar in performance to another recently released tool, PICRUSt, for almost all of survey data analysed here. But, our method is unique in that any OTU building method can be used, as opposed to being limited to closed-reference OTU picking strategies against specific reference sequence databases., Conclusions: We developed an automated computational method, which derives an inferred functional profile based on the 16S rRNA gene surveys of microbial communities. The inferred functional profile provides a cost effective way to study complex ecosystems through predicted comparative functional metagenomes and metadata analysis. All PanFP source code and additional documentation are freely available online at GitHub ( https://github.com/srjun/PanFP ).

The 2014 Ebola outbreak in West Africa is the largest documented for this virus. To examine the dynamics of this genome, we compare more than 100 currently available ebolavirus genomes to each other and to other viral genomes. Based on oligomer frequency analysis, the family Filoviridae forms a distinct group from all other sequenced viral genomes. All filovirus genomes sequenced to date encode proteins with similar functions and gene order, although there is considerable divergence in sequences between the three genera Ebolavirus, Cuevavirus and Marburgvirus within the family Filoviridae. Whereas all ebolavirus genomes are quite similar (multiple sequences of the same strain are often identical), variation is most common in the intergenic regions and within specific areas of the genes encoding the glycoprotein (GP), nucleoprotein (NP) and polymerase (L). We predict regions that could contain epitope-binding sites, which might be good vaccine targets. This information, combined with glycosylation sites and experimentally determined epitopes, can identify the most promising regions for the development of therapeutic strategies.This manuscript has been authored by UT-Battelle, LLC under Contract No. DE-AC05-00OR22725 with the U.S. Department of Energy. The United States Government retains and the publisher, by accepting the article for publication, acknowledges that the United States Government retains a non-exclusive, paid-up, irrevocable, world-wide license to publish or reproduce the published form of this manuscript, or allow others to do so, for United States Government purposes. The Department of Energy will provide public access to these results of federally sponsored research in accordance with the DOE Public Access Plan (http://energy.gov/downloads/doe-public-access-plan)., (© FEMS 2015.)

Since the first two complete bacterial genome sequences were published in 1995, the science of bacteria has dramatically changed. Using third-generation DNA sequencing, it is possible to completely sequence a bacterial genome in a few hours and identify some types of methylation sites along the genome as well. Sequencing of bacterial genome sequences is now a standard procedure, and the information from tens of thousands of bacterial genomes has had a major impact on our views of the bacterial world. In this review, we explore a series of questions to highlight some insights that comparative genomics has produced. To date, there are genome sequences available from 50 different bacterial phyla and 11 different archaeal phyla. However, the distribution is quite skewed towards a few phyla that contain model organisms. But the breadth is continuing to improve, with projects dedicated to filling in less characterized taxonomic groups. The clustered regularly interspaced short palindromic repeats (CRISPR)-Cas system provides bacteria with immunity against viruses, which outnumber bacteria by tenfold. How fast can we go? Second-generation sequencing has produced a large number of draft genomes (close to 90 % of bacterial genomes in GenBank are currently not complete); third-generation sequencing can potentially produce a finished genome in a few hours, and at the same time provide methlylation sites along the entire chromosome. The diversity of bacterial communities is extensive as is evident from the genome sequences available from 50 different bacterial phyla and 11 different archaeal phyla. Genome sequencing can help in classifying an organism, and in the case where multiple genomes of the same species are available, it is possible to calculate the pan- and core genomes; comparison of more than 2000 Escherichia coli genomes finds an E. coli core genome of about 3100 gene families and a total of about 89,000 different gene families. Why do we care about bacterial genome sequencing? There are many practical applications, such as genome-scale metabolic modeling, biosurveillance, bioforensics, and infectious disease epidemiology. In the near future, high-throughput sequencing of patient metagenomic samples could revolutionize medicine in terms of speed and accuracy of finding pathogens and knowing how to treat them.

Background: More than 80% of the microbial genomes in GenBank are of 'draft' quality (12,553 draft vs. 2,679 finished, as of October, 2013). We have examined all the microbial DNA sequences available for complete, draft, and Sequence Read Archive genomes in GenBank as well as three other major public databases, and assigned quality scores for more than 30,000 prokaryotic genome sequences., Results: Scores were assigned using four categories: the completeness of the assembly, the presence of full-length rRNA genes, tRNA composition and the presence of a set of 102 conserved genes in prokaryotes. Most (~88%) of the genomes had quality scores of 0.8 or better and can be safely used for standard comparative genomics analysis. We compared genomes across factors that may influence the score. We found that although sequencing depth coverage of over 100x did not ensure a better score, sequencing read length was a better indicator of sequencing quality. With few exceptions, most of the 30,000 genomes have nearly all the 102 essential genes., Conclusions: The score can be used to set thresholds for screening data when analyzing "all published genomes" and reference data is either not available or not applicable. The scores highlighted organisms for which commonly used tools do not perform well. This information can be used to improve tools and to serve a broad group of users as more diverse organisms are sequenced. Unexpectedly, the comparison of predicted tRNAs across 15,000 high quality genomes showed that anticodons beginning with an 'A' (codons ending with a 'U') are almost non-existent, with the exception of one arginine codon (CGU); this has been noted previously in the literature for a few genomes, but not with the depth found here.

Purpose: To evaluate the feasibility and efficacy of antegrade renal cooling during radiofrequency (RF) ablation to protect the renal collecting system in pigs., Materials and Methods: With the approval of the institutional animal care and use committee, RF ablation was planned bilaterally after nephrostomy in 20 kidneys from 10 domestic pigs. A coaxial assembly consisting of a 7-F pigtail catheter and a 12-F arterial sheath was placed at the renal pelvis. Cooled 5% dextrose in water solution was perfused in only one of the kidneys in each pig. A 17-gauge internally cooled tip RF electrode with a 3-cm tip was placed near the renal pelvis with ultrasonographic guidance. On postprocedural days 1 and 7, computed tomography (CT) was performed to evaluate the RF ablation effect, the presence of a fluid collection, and the leakage of contrast material. Ex vivo retrograde pyelography was performed to assess the integrity of the urinary tract after kidney harvest. Histologic findings related to the urothelial damage were scored semiquantitatively., Results: The mean maximal diameter of the RF ablation area measured at CT and in specimens was not significantly different between cooled and noncooled kidneys. A fluid collection around the RF ablation area was seen in eight of the 10 noncooled kidneys on day 7 CT scans compared with two of the 10 cooled kidneys; the difference was statistically significant (P = .035). Ex vivo retrograde pyelography depicted leakage of contrast material in only the noncooled group. The mean histologic damage score was lower in the cooled group than in the noncooled group., Conclusion: RF ablation after placement of an antegrade cooling system in pig kidneys was feasible. Antegrade cold dextrose infusion of the urinary tract during RF ablation is effective in protecting the renal collecting system without compromising the RF ablation effect., ((c) RSNA, 2010.)

We present a whole-proteome phylogeny of prokaryotes constructed by comparing feature frequency profiles (FFPs) of whole proteomes. Features are l-mers of amino acids, and each organism is represented by a profile of frequencies of all features. The selection of feature length is critical in the FFP method, and we have developed a procedure for identifying the optimal feature lengths for inferring the phylogeny of prokaryotes, strictly speaking, a proteome phylogeny. Our FFP trees are constructed with whole proteomes of 884 prokaryotes, 16 unicellular eukaryotes, and 2 random sequences. To highlight the branching order of major groups, we present a simplified proteome FFP tree of monophyletic class or phylum with branch support. In our whole-proteome FFP trees (i) Archaea, Bacteria, Eukaryota, and a random sequence outgroup are clearly separated; (ii) Archaea and Bacteria form a sister group when rooted with random sequences; (iii) Planctomycetes, which possesses an intracellular membrane compartment, is placed at the basal position of the Bacteria domain; (iv) almost all groups are monophyletic in prokaryotes at most taxonomic levels, but many differences in the branching order of major groups are observed between our proteome FFP tree and trees built with other methods; and (v) previously "unclassified" genomes may be assigned to the most likely taxa. We describe notable similarities and differences between our FFP trees and those based on other methods in grouping and phylogeny of prokaryotes.

Ten complete mammalian genome sequences were compared by using the "feature frequency profile" (FFP) method of alignment-free comparison. This comparison technique reveals that the whole nongenic portion of mammalian genomes contains evolutionary information that is similar to their genic counterparts--the intron and exon regions. We partitioned the complete genomes of mammals (such as human, chimp, horse, and mouse) into their constituent nongenic, intronic, and exonic components. Phylogenic species trees were constructed for each individual component class of genome sequence data as well as the whole genomes by using standard tree-building algorithms with FFP distances. The phylogenies of the whole genomes and each of the component classes (exonic, intronic, and nongenic regions) have similar topologies, within the optimal feature length range, and all agree well with the evolutionary phylogeny based on a recent large dataset, multispecies, and multigene-based alignment. In the strictest sense, the FFP-based trees are genome phylogenies, not species phylogenies. However, the species phylogeny is highly related to the whole-genome phylogeny. Furthermore, our results reveal that the footprints of evolutionary history are spread throughout the entire length of the whole genome of an organism and are not limited to genes, introns, or short, highly conserved, nongenic sequences that can be adversely affected by factors (such as a choice of sequences, homoplasy, and different mutation rates) resulting in inconsistent species phylogenies.

The vast sequence divergence among different virus groups has presented a great challenge to alignment-based sequence comparison among different virus families. Using an alignment-free comparison method, we construct the whole-proteome phylogeny for a population of viruses from 11 viral families comprising 142 large dsDNA eukaryote viruses. The method is based on the feature frequency profiles (FFP), where the length of the feature (l-mer) is selected to be optimal for phylogenomic inference. We observe that (i) the FFP phylogeny segregates the population into clades, the membership of each has remarkable agreement with current classification by the International Committee on the Taxonomy of Viruses, with one exception that the mimivirus joins the phycodnavirus family; (ii) the FFP tree detects potential evolutionary relationships among some viral families; (iii) the relative position of the 3 herpesvirus subfamilies in the FFP tree differs from gene alignment-based analysis; (iv) the FFP tree suggests the taxonomic positions of certain "unclassified" viruses; and (v) the FFP method identifies candidates for horizontal gene transfer between virus families.

For comparison of whole-genome (genic + nongenic) sequences, multiple sequence alignment of a few selected genes is not appropriate. One approach is to use an alignment-free method in which feature (or l-mer) frequency profiles (FFP) of whole genomes are used for comparison-a variation of a text or book comparison method, using word frequency profiles. In this approach it is critical to identify the optimal resolution range of l-mers for the given set of genomes compared. The optimum FFP method is applicable for comparing whole genomes or large genomic regions even when there are no common genes with high homology. We outline the method in 3 stages: (i) We first show how the optimal resolution range can be determined with English books which have been transformed into long character strings by removing all punctuation and spaces. (ii) Next, we test the robustness of the optimized FFP method at the nucleotide level, using a mutation model with a wide range of base substitutions and rearrangements. (iii) Finally, to illustrate the utility of the method, phylogenies are reconstructed from concatenated mammalian intronic genomes; the FFP derived intronic genome topologies for each l within the optimal range are all very similar. The topology agrees with the established mammalian phylogeny revealing that intron regions contain a similar level of phylogenic signal as do coding regions.

The extrahepatic biliary system is one of the unusual locations of extrapulmonary neuroendocrine carcinoma. To our knowledge, there has been no report describing the imaging findings of this rare malignancy of the biliary tree. We report 3 large-cell neuroendocrine carcinomas and a small-cell carcinoma arising from the biliary system. Biliary tumors showing on computed tomography an aggressive invasion to the liver or extensive lymph node involvement should raise the suspicion of their being high-grade neuroendocrine tumors.

We have constructed a map of the "protein structure space" by using the pairwise structural similarity scores calculated for all nonredundant protein structures determined experimentally. As expected, proteins with similar structures clustered together in the map and the overall distribution of structural classes of this map followed closely that of the map of the "protein fold space" we have reported previously. Consequently, proteins sharing similar molecular functions also were found to colocalize in the protein structure space map, pointing toward a previously undescribed scheme for structure-based functional inference for remote homologues based on the proximity in the map of the protein structure space. We found that this scheme consistently outperformed other predictions made by using either the raw scores or normalized Z-scores of pairwise DALI structure alignment.

The ORF sequences of the gene encoding sepiapterin reductase were cloned from the genomic DNAs of Chlorobium tepidum and Chlorobium limicola, which are known to produce L-threo- and L-erythro-tetrahydrobiopterin (BH4)-N-acetylglucosamine, respectively. The deduced amino acid sequence of C. limicola consists of 241 residues, while C. tepidum SR has three residues more at the C-terminal. The overall protein sequence identity was 87.7%. Both recombinant proteins generated from Escherichia coli were identified to catalyze reduction of diketo compound 6-pyruvoyltetrahydropterin to L-threo-BH4. This result suggests that C. limicola needs an additional enzyme for L-erythro-BH4 synthesis to yield its glycoside. The catalytic activity of Chlorobium SRs also supports the previously proposed mechanism of two consecutive reductions of C1' carbonyl group of 6-pyruvoyltetrahydropterin via isomerization reaction.