OMeta: an ontology-based, data-driven metadata tracking system.

BACKGROUND: The development of high-throughput sequencing and analysis has accelerated multi-omics studies of thousands of microbial species, metagenomes, and infectious disease pathogens. Omics studies are enabling genotype-phenotype association studies which identify genetic determinants of pathogen virulence and drug resistance, as well as phylogenetic studies designed to track the origin and spread of disease outbreaks. These omics studies are complex and often employ multiple assay technologies including genomics, metagenomics, transcriptomics, proteomics, and metabolomics. To maximize the impact of omics studies, it is essential that data be accompanied by detailed contextual metadata (e.g., specimen, spatial-temporal, phenotypic characteristics) in clear, organized, and consistent formats. Over the years, many metadata standards developed by various metadata standards initiatives have arisen; the Genomic Standards Consortium's minimal information standards (MIxS), the GSCID/BRC Project and Sample Application Standard. Some tools exist for tracking metadata, but they do not provide event based capabilities to configure, collect, validate, and distribute metadata. To address this gap in the scientific community, an event based data-driven application, OMeta, was created that allows users to quickly configure, collect, validate, distribute, and integrate metadata. RESULTS: A data-driven web application, OMeta, has been developed for use by researchers consisting of a browser-based interface, a command-line interface (CLI), and server-side components that provide an intuitive platform for configuring, capturing, viewing, and sharing metadata. Project and sample metadata can be set based on existing standards or based on projects goals. Recorded information includes details on the biological samples, procedures, protocols, and experimental technologies, etc. This information can be organized based on events, including sample collection, sample quantification, sequencing assay, and analysis results. OMeta enables configuration in various presentation types: checkbox, file, drop-box, ontology, and fields can be configured to use the National Center for Biomedical Ontology (NCBO), a biomedical ontology server. Furthermore, OMeta maintains a complete audit trail of all changes made by users and allows metadata export in comma separated value (CSV) format for convenient deposition of data into public databases. CONCLUSIONS: We present, OMeta, a web-based software application that is built on data-driven principles for configuring and customizing data standards, capturing, curating, and sharing metadata.

Gastrointestinal microbial populations can distinguish pediatric and adolescent Acute Lymphoblastic Leukemia (ALL) at the time of disease diagnosis.

BACKGROUND: An estimated 15,000 children and adolescents under the age of 19 years are diagnosed with leukemia, lymphoma and other tumors in the USA every year. All children and adolescent acute leukemia patients will undergo chemotherapy as part of their treatment regimen. Fortunately, survival rates for most pediatric cancers have improved at a remarkable pace over the past three decades, and the overall survival rate is greater than 90 % today. However, significant differences in survival rate have been found in different age groups (94 % in 1-9.99 years, 82 % in ≥10 years and 76 % in ≥15 years). ALL accounts for about three out of four cases of childhood leukemia. Intensive chemotherapy treatment coupled with prophylactic or therapeutic antibiotic use could potentially have a long-term effect on the resident gastrointestinal (GI) microbiome. The composition of GI microbiome and its changes upon chemotherapy in pediatric and adolescent leukemia patients is poorly understood. In this study, using 16S rRNA marker gene sequences we profile the GI microbial communities of pediatric and adolescent acute leukemia patients before and after chemotherapy treatment and compare with the microbiota of their healthy siblings. RESULTS: Our study cohort consisted of 51 participants, made up of matched pediatric and adolescent patients with ALL and a healthy sibling. We elucidated and compared the GI microbiota profiles of patients and their healthy sibling controls via analysis of 16S rRNA gene sequencing data. We assessed the GI microbiota composition in pediatric and adolescent patients with ALL during the course of chemotherapy by comparing stool samples taken before chemotherapy with stool samples collected at varying time points during the chemotherapeutic treatment. The microbiota profiles of both patients and control sibling groups are dominated by members of Bacteroides, Prevotella, and Faecalibacterium. At the genus level, both groups share many taxa in common, but the microbiota diversity of the patient group is significantly lower than that of the control group. It was possible to distinguish between the patient and control groups based on their microbiota profiles. The top taxa include Anaerostipes, Coprococcus, Roseburia, and Ruminococcus2 with relatively higher abundance in the control group. The observed microbiota changes are likely the result of several factors including a direct influence of therapeutic compounds on the gut flora and an indirect effect of chemotherapy on the immune system, which, in turn, affects the microbiota. CONCLUSIONS: This study provides significant information on GI microbiota populations in immunocompromised children and opens up the potential for developing novel diagnostics based on stool tests and therapies to improve the dysbiotic condition of the microbiota at the time of diagnosis and in the earliest stages of chemotherapy.

Stool microbiota composition is associated with the prospective risk of Plasmodium falciparum infection.

BACKGROUND: In humans it is unknown if the composition of the gut microbiota alters the risk of Plasmodium falciparum infection or the risk of developing febrile malaria once P. falciparum infection is established. Here we collected stool samples from a cohort composed of 195 Malian children and adults just prior to an intense P. falciparum transmission season. We assayed these samples using massively parallel sequencing of the 16S ribosomal RNA gene to identify the composition of the gut bacterial communities in these individuals. During the ensuing 6-month P. falciparum transmission season we examined the relationship between the stool microbiota composition of individuals in this cohort and their prospective risk of both P. falciparum infection and febrile malaria. RESULTS: Consistent with prior studies, stool microbial diversity in the present cohort increased with age, although the overall microbiota profile was distinct from cohorts in other regions of Africa, Asia and North America. Age-adjusted Cox regression analysis revealed a significant association between microbiota composition and the prospective risk of P. falciparum infection; however, no relationship was observed between microbiota composition and the risk of developing febrile malaria once P. falciparum infection was established. CONCLUSIONS: These findings underscore the diversity of gut microbiota across geographic regions, and suggest that strategic modulation of gut microbiota composition could decrease the risk of P. falciparum infection in malaria-endemic areas, potentially as an adjunct to partially effective malaria vaccines.

TIGRFAMs and Genome Properties in 2013.

TIGRFAMs, available online at http://www.jcvi.org/tigrfams is a database of protein family definitions. Each entry features a seed alignment of trusted representative sequences, a hidden Markov model (HMM) built from that alignment, cutoff scores that let automated annotation pipelines decide which proteins are members, and annotations for transfer onto member proteins. Most TIGRFAMs models are designated equivalog, meaning they assign a specific name to proteins conserved in function from a common ancestral sequence. Models describing more functionally heterogeneous families are designated subfamily or domain, and assign less specific but more widely applicable annotations. The Genome Properties database, available at http://www.jcvi.org/genome-properties, specifies how computed evidence, including TIGRFAMs HMM results, should be used to judge whether an enzymatic pathway, a protein complex or another type of molecular subsystem is encoded in a genome. TIGRFAMs and Genome Properties content are developed in concert because subsystems reconstruction for large numbers of genomes guides selection of seed alignment sequences and cutoff values during protein family construction. Both databases specialize heavily in bacterial and archaeal subsystems. At present, 4284 models appear in TIGRFAMs, while 628 systems are described by Genome Properties. Content derives both from subsystem discovery work and from biocuration of the scientific literature.

CharProtDB: a database of experimentally characterized protein annotations.

CharProtDB (http://www.jcvi.org/charprotdb/) is a curated database of biochemically characterized proteins. It provides a source of direct rather than transitive assignments of function, designed to support automated annotation pipelines. The initial data set in CharProtDB was collected through manual literature curation over the years by analysts at the J. Craig Venter Institute (JCVI) [formerly The Institute of Genomic Research (TIGR)] as part of their prokaryotic genome sequencing projects. The CharProtDB has been expanded by import of selected records from publicly available protein collections whose biocuration indicated direct rather than homology-based assignment of function. Annotations in CharProtDB include gene name, symbol and various controlled vocabulary terms, including Gene Ontology terms, Enzyme Commission number and TransportDB accession. Each annotation is referenced with the source; ideally a journal reference, or, if imported and lacking one, the original database source.

Draft genome sequence determination for cystic fibrosis and chronic granulomatous disease Burkholderia multivorans isolates.

Burkholderia multivorans is a Gram-negative bacterium and a member of the Burkholderia cepacia complex, which is frequently associated with respiratory infections in people with cystic fibrosis (CF) and chronic granulomatous disease (CGD). We are reporting the genome sequences of 4 B. multivorans strains, 2 from CF patients and 2 from CGD patients.

The Protein Naming Utility: a rules database for protein nomenclature.

Generation of syntactically correct and unambiguous names for proteins is a challenging, yet vital task for functional annotation processes. Proteins are often named based on homology to known proteins, many of which have problematic names. To address the need to generate high-quality protein names, and capture our significant experience correcting protein names manually, we have developed the Protein Naming Utility (PNU, http://www.jcvi.org/pn-utility). The PNU is a web-based database for storing and applying naming rules to identify and correct syntactically incorrect protein names, or to replace synonyms with their preferred name. The PNU allows users to generate and manage collections of naming rules, optionally building upon the growing body of rules generated at the J. Craig Venter Institute (JCVI). Since communities often enforce disparate conventions for naming proteins, the PNU supports grouping rules into user-managed collections. Users can check their protein names against a selected PNU rule collection, generating both statistics and corrected names. The PNU can also be used to correct GenBank table files prior to submission to GenBank. Currently, the database features 3080 manual rules that have been entered by JCVI Bioinformatics Analysts as well as 7458 automatically imported names.

Pathema: a clade-specific bioinformatics resource center for pathogen research.

Pathema (http://pathema.jcvi.org) is one of the eight Bioinformatics Resource Centers (BRCs) funded by the National Institute of Allergy and Infectious Disease (NIAID) designed to serve as a core resource for the bio-defense and infectious disease research community. Pathema strives to support basic research and accelerate scientific progress for understanding, detecting, diagnosing and treating an established set of six target NIAID Category A-C pathogens: Category A priority pathogens; Bacillus anthracis and Clostridium botulinum, and Category B priority pathogens; Burkholderia mallei, Burkholderia pseudomallei, Clostridium perfringens and Entamoeba histolytica. Each target pathogen is represented in one of four distinct clade-specific Pathema web resources and underlying databases developed to target the specific data and analysis needs of each scientific community. All publicly available complete genome projects of phylogenetically related organisms are also represented, providing a comprehensive collection of organisms for comparative analyses. Pathema facilitates the scientific exploration of genomic and related data through its integration with web-based analysis tools, customized to obtain, display, and compute results relevant to ongoing pathogen research. Pathema serves the bio-defense and infectious disease research community by disseminating data resulting from pathogen genome sequencing projects and providing access to the results of inter-genomic comparisons for these organisms.

Draft genome sequence of Bacteroides vulgatus PC510, a strain isolated from human feces.

Although Bacteroides vulgatus is one of the most prevalent microorganisms in the human gastrointestinal tract, little is known about the genetic potential of this species. Here, we describe the annotated draft genome sequence of B. vulgatus PC510 isolated from human feces.

Draft genome sequence of Turicibacter sanguinis PC909, isolated from human feces.

While the microbiota resident in the human gut is now known to provide a range of functions relevant to host health, many of the microbial members of the community have not yet been cultured or are represented by a limited number of isolates. We describe here the draft genome sequence of Turicibacter sanguinis PC909, isolated from a pooled healthy human fecal sample as part of the Australian Human Gut Microbiome Project.

Advances in Microbiome Research for Animal Health.

Host-associated microbiomes contribute in many ways to the homeostasis of the metaorganism. The microbiome's contributions range from helping to provide nutrition and aiding growth, development, and behavior to protecting against pathogens and toxic compounds. Here we summarize the current knowledge of the diversity and importance of the microbiome to animals, using representative examples of wild and domesticated species. We demonstrate how the beneficial ecological roles of animal-associated microbiomes can be generally grouped into well-defined main categories and how microbe-based alternative treatments can be applied to mitigate problems for both economic and conservation purposes and to provide crucial knowledge about host-microbiota symbiotic interactions. We suggest a Customized Combination of Microbial-Based Therapies to promote animal health and contribute to the practice of sustainable husbandry. We also discuss the ecological connections and threats associated with animal biodiversity loss, microorganism extinction, and emerging diseases, such as the COVID-19 pandemic.

Longitudinal Study of Oral Microbiome Variation in Twins.

Humans are host to a multitude of microorganisms that rapidly populate the body at birth, subject to a complex interplay that is dependent on host genetics, lifestyle, and environment. The host-associated microbiome, including the oral microbiome, presents itself in a complex ecosystem important to health and disease. As the most common chronic disease globally, dental caries is induced by host-microbial dysbiosis in children and adults. Multiple biological and environmental factors are likely to impact disease predisposition, onset, progression, and severity, yet longitudinal studies able to capture these influences are missing. To investigate how host genetics and environment influenced the oral microbial communities over time, we profiled supragingival plaque microbiomes of dizygotic and monozygotic twins during 3 visits over 12-months. Dental plaque DNA samples were amplified by targeting the 16S rRNA gene V4 region, and microbial findings were correlated with clinical, diet and genetic metadata. We observed that the oral microbiome variances were shaped primarily by the environment when compared to host genetics. Among the environmental factors shaping microbial changes of our subjects, significant metadata included age of the subject, and the age by which subjects initiated brushing habits, and the types of actions post-brushing. Relevant heritability of the microbiome included Actinomyces and Capnocytophaga in monozygotic twins and Kingella in dizygotic twins. Corynebacterium and Veillonella abundances were associated with age, whereas Aggregatibacter was associated with younger subjects. Streptococcus abundance showed an inverse association over time, and Selenomonas abundances increased with brushing frequency per day. Unraveling the exact biological mechanisms in caries has the potential to reveal novel host-microbial biomarkers, pathways, and targets important to effective preventive measures, and early disease control in children.

Xylitol and sorbitol effects on the microbiome of saliva and plaque.

Chewing gum containing xylitol may help prevent caries by reducing levels of mutans streptococci (MS) and lactobacilli in saliva and plaque. Very little is known about other species which are possibly beneficial to oral health. In this study, we employed high-throughput sequencing of the 16S rRNA gene to profile microbial communities of saliva and plaque following short-term consumption of xylitol and sorbitol containing chewing gum. Participants (n = 30) underwent a washout period and were randomly assigned to one of two groups. Each group chewed either xylitol or sorbitol gum for three weeks, before undergoing a second four-week washout period after which they switched to the alternate gum for three weeks. Analysis of samples collected before and after each intervention identified distinct plaque and saliva microbial communities that altered dependent on the order in which gum treatments were given. Neither the xylitol nor sorbitol treatments significantly affected the bacterial composition of plaque. Lactobacilli were undetected and the number of Streptococcus mutans sequence reads was very low and unaffected by either xylitol or sorbitol. However, sorbitol affected several other streptococcal species in saliva including increasing the abundance of S. cristatus, an oral commensal shown to inhibit bacteria associated with chronic periodontitis.

Insights into the environmental resistance gene pool from the genome sequence of the multidrug-resistant environmental isolate Escherichia coli SMS-3-5.

The increasing occurrence of multidrug-resistant pathogens of clinical and agricultural importance is a global public health concern. While antimicrobial use in human and veterinary medicine is known to contribute to the dissemination of antimicrobial resistance, the impact of microbial communities and mobile resistance genes from the environment in this process is not well understood. Isolated from an industrially polluted aquatic environment, Escherichia coli SMS-3-5 is resistant to a record number of antimicrobial compounds from all major classes, including two front-line fluoroquinolones (ciprofloxacin and moxifloxacin), and in many cases at record-high concentrations. To gain insights into antimicrobial resistance in environmental bacterial populations, the genome of E. coli SMS-3-5 was sequenced and compared to the genome sequences of other E. coli strains. In addition, selected genetic loci from E. coli SMS-3-5 predicted to be involved in antimicrobial resistance were phenotypically characterized. Using recombinant vector clones from shotgun sequencing libraries, resistance to tetracycline, streptomycin, and sulfonamide/trimethoprim was assigned to a single mosaic region on a 130-kb plasmid (pSMS35_130). The remaining plasmid backbone showed similarity to virulence plasmids from avian-pathogenic E. coli (APEC) strains. Individual resistance gene cassettes from pSMS35_130 are conserved among resistant bacterial isolates from multiple phylogenetic and geographic sources. Resistance to quinolones was assigned to several chromosomal loci, mostly encoding transport systems that are also present in susceptible E. coli isolates. Antimicrobial resistance in E. coli SMS-3-5 is therefore dependent both on determinants acquired from a mobile gene pool that is likely available to clinical and agricultural pathogens, as well, and on specifically adapted multidrug efflux systems. The association of antimicrobial resistance with APEC virulence genes on pSMS35_130 highlights the risk of promoting the spread of virulence through the extensive use of antibiotics.

A host subtraction database for virus discovery in human cell line sequencing data.

The human cell lines HepG2, HuH-7, and Jurkat are commonly used for amplification of the RNA viruses present in environmental samples. To assist with assays by RNAseq, we sequenced these cell lines and developed a subtraction database that contains sequences expected in sequence data from uninfected cells. RNAseq data from cell lines infected with Sendai virus were analyzed to test host subtraction. The process of mapping RNAseq reads to our subtraction database vastly reduced the number non-viral reads in the dataset to allow for efficient secondary analyses.

A draft genome sequence for the Ixodes scapularis cell line, ISE6.

Background: The tick cell line ISE6, derived from Ixodes scapularis, is commonly used for amplification and detection of arboviruses in environmental or clinical samples. Methods: To assist with sequence-based assays, we sequenced the ISE6 genome with single-molecule, long-read technology. Results: The draft assembly appears near complete based on gene content analysis, though it appears to lack some instances of repeats in this highly repetitive genome. The assembly appears to have separated the haplotypes at many loci. DNA short read pairs, used for validation only, mapped to the cell line assembly at a higher rate than they mapped to the Ixodes scapularis reference genome sequence. Conclusions: The assembly could be useful for filtering host genome sequence from sequence data obtained from cells infected with pathogens.

Supragingival Plaque Microbiome Ecology and Functional Potential in the Context of Health and Disease.

To address the question of how microbial diversity and function in the oral cavities of children relates to caries diagnosis, we surveyed the supragingival plaque biofilm microbiome in 44 juvenile twin pairs. Using shotgun sequencing, we constructed a genome encyclopedia describing the core supragingival plaque microbiome. Caries phenotypes contained statistically significant enrichments in specific genome abundances and distinct community composition profiles, including strain-level changes. Metabolic pathways that are statistically associated with caries include several sugar-associated phosphotransferase systems, antimicrobial resistance, and metal transport. Numerous closely related previously uncharacterized microbes had substantial variation in central metabolism, including the loss of biosynthetic pathways resulting in auxotrophy, changing the ecological role. We also describe the first complete Gracilibacteria genomes from the human microbiome. Caries is a microbial community metabolic disorder that cannot be described by a single etiology, and our results provide the information needed for next-generation diagnostic tools and therapeutics for caries.IMPORTANCE Oral health has substantial economic importance, with over $100 billion spent on dental care in the United States annually. The microbiome plays a critical role in oral health, yet remains poorly classified. To address the question of how microbial diversity and function in the oral cavities of children relate to caries diagnosis, we surveyed the supragingival plaque biofilm microbiome in 44 juvenile twin pairs. Using shotgun sequencing, we constructed a genome encyclopedia describing the core supragingival plaque microbiome. This unveiled several new previously uncharacterized but ubiquitous microbial lineages in the oral microbiome. Caries is a microbial community metabolic disorder that cannot be described by a single etiology, and our results provide the information needed for next-generation diagnostic tools and therapeutics for caries.

Analysis of the Aedes albopictus C6/36 genome provides insight into cell line utility for viral propagation.

Background: The 50-year-old Aedes albopictus C6/36 cell line is a resource for the detection, amplification, and analysis of mosquito-borne viruses including Zika, dengue, and chikungunya. The cell line is derived from an unknown number of larvae from an unspecified strain of Aedes albopictus mosquitoes. Toward improved utility of the cell line for research in virus transmission, we present an annotated assembly of the C6/36 genome. Results: The C6/36 genome assembly has the largest contig N50 (3.3 Mbp) of any mosquito assembly, presents the sequences of both haplotypes for most of the diploid genome, reveals independent null mutations in both alleles of the Dicer locus, and indicates a male-specific genome. Gene annotation was computed with publicly available mosquito transcript sequences. Gene expression data from cell line RNA sequence identified enrichment of growth-related pathways and conspicuous deficiency in aquaporins and inward rectifier K+ channels. As a test of utility, RNA sequence data from Zika-infected cells were mapped to the C6/36 genome and transcriptome assemblies. Host subtraction reduced the data set by 89%, enabling faster characterization of nonhost reads. Conclusions: The C6/36 genome sequence and annotation should enable additional uses of the cell line to study arbovirus vector interactions and interventions aimed at restricting the spread of human disease.

Initial genome sequencing of the sugarcane CP 96-1252 complex hybrid.

The CP 96-1252 cultivar of sugarcane is a complex hybrid of commercial importance. DNA was extracted from lab-grown leaf tissue and sequenced. The raw Illumina DNA sequencing results provide 101 Gbp of genome sequence reads. The dataset is available from https://www.ncbi.nlm.nih.gov/bioproject/PRJNA345486/.

Type 1 Diabetes: Urinary Proteomics and Protein Network Analysis Support Perturbation of Lysosomal Function.

While insulin replacement therapy restores the health and prevents the onset of diabetic complications (DC) for many decades, some T1D patients have elevated hemoglobin A1c values suggesting poor glycemic control, a risk factor of DC. We surveyed the stool microbiome and urinary proteome of a cohort of 220 adolescents and children, half of which had lived with T1D for an average of 7 years and half of which were healthy siblings. Phylogenetic analysis of the 16S rRNA gene did not reveal significant differences in gut microbial alpha-diversity comparing the two cohorts. The urinary proteome of T1D patients revealed increased abundances of several lysosomal proteins that correlated with elevated HbA1c values. In silico protein network analysis linked such proteins to extracellular matrix components and the glycoprotein LRG1. LRG1 is a prominent inflammation and neovascularization biomarker. We hypothesize that these changes implicate aberrant glycation of macromolecules that alter lysosomal function and metabolism in renal tubular epithelial cells, cells that line part of the upper urinary tract.