Standardized metadata for human pathogen/vector genomic sequences

High throughput sequencing has accelerated the determination of genome sequences for thousands of human infectious disease pathogens and dozens of their vectors. The scale and scope of these data are enabling genotype-phenotype association studies to identify genetic determinants of pathogen virulence and drug/insecticide resistance, and phylogenetic studies to track the origin and spread of disease outbreaks. To maximize the utility of genomic sequences for these purposes, it is essential that metadata about the pathogen/vector isolate characteristics be collected and made available in organized, clear, and consistent formats. Here we report the development of the GSCID/BRC Project and Sample Application Standard, developed by representatives of the Genome Sequencing Centers for Infectious Diseases (GSCIDs), the Bioinformatics Resource Centers (BRCs) for Infectious Diseases, and the U.S. National Institute of Allergy and Infectious Diseases (NIAID), part of the National Institutes of Health (NIH), informed by interactions with numerous collaborating scientists. It includes mapping to terms from other data standards initiatives, including the Genomic Standards Consortium's minimal information (MIxS) and NCBI's BioSample/BioProjects checklists and the Ontology for Biomedical Investigations (OBI). The standard includes data fields about characteristics of the organism or environmental source of the specimen, spatial-temporal information about the specimen isolation event, phenotypic characteristics of the pathogen/vector isolated, and project leadership and support. By modeling metadata fields into an ontology-based semantic framework and reusing existing ontologies and minimum information checklists, the application standard can be extended to support additional project-specific data fields and integrated with other data represented with comparable standards. The use of this metadata standard by all ongoing and future GSCID sequencing projects will provide a consistent representation of these data in the BRC resources and other repositories that leverage these data, allowing investigators to identify relevant genomic sequences and perform comparative genomics analyses that are both statistically meaningful and biologically relevant

Genetic Variation In Vitro and In Vivo of an Attenuated Lassa Vaccine Candidate

The attenuated Lassa vaccine candidate ML29 is a laboratory-produced reassortant between Lassa and Mopeia viruses, two Old World arenaviruses that differ by 40% in nucleic acid sequence. In our previous studies, ML29 elicited sterilizing immunity against Lassa virus challenge in guinea pigs and marmosets and virus-specific cell-mediated immunity in both simian immunodeficiency virus (SIV)-infected and uninfected rhesus macaques. Here, we show that ML29 is stable after 12 passages in vitro without losing its plaque morphology or its attenuated phenotype in suckling mice. Additionally, we used deep sequencing to characterize the viral population comprising the original stock of ML29, the stock of ML29 after 12 passages in Vero cells, and the ML29 isolates obtained from vaccinated animals. Twenty-seven isolates bore approximately 77 mutations that exceeded 20% of the single-nucleotide polymorphism (SNP) changes at any single locus. Of these 77 mutations, 5 appeared to be host specific, for example, appearing in mice but not in primates. None of these mutations were reversions of ML29 to the sequences of the parental Lassa and Mopeia viruses. The host-specific mutations indicate viral adaptations to virus-host interactions, and such interactions make reasonable targets for antiviral approaches. Variants capable of chronic infection did not emerge from any of the primate infections, even in immune-deficient animals, indicating that the ML29 reassortant is reasonably stable in vivo. In conclusion, the preclinical studies of ML29 as a Lassa virus vaccine candidate have been advanced, showing high levels of protection in nonhuman primates and acceptable stability both in vitro and in vivo. © 2014, American Society for Microbiology

Evolution of Mycobacterium abscessus in the human lung: Cumulative mutations and genomic rearrangement of porin genes in patient isolates

ABSTRACTBackground Mycobacterium abscessus subspecies massiliense (M. massiliense) is increasingly recognized as an emerging bacterial pathogen, particularly in cystic fibrosis (CF) patients and CF centres’ respiratory outbreaks. We characterized genomic and phenotypic changes in 15 serial isolates from two CF patients (1S and 2B) with chronic pulmonary M. massiliense infection leading to death, as well as four isolates from a CF centre outbreak in which patient 2B was the index case.Results Comparative genomic analysis revealed the mutations affecting growth rate, metabolism, transport, lipids (loss of glycopeptidolipids), antibiotic susceptibility (macrolides and aminoglycosides resistance), and virulence factors. Mutations in 23S rRNA, mmpL4, porin locus and tetR genes occurred in isolates from both CF patients. Interestingly, we identified two different spontaneous mutation events at the mycobacterial porin locus: a fusion of two tandem porin paralogs in patient 1S and a partial deletion of the first porin paralog in patient 2B. These genomic changes correlated with reduced porin protein expression, diminished 14C-glucose uptake, slower bacterial growth rates, and enhanced TNF-α induction in mycobacteria-infected THP-1 human cells. Porin gene complementation of porin mutants partly restored 14C-glucose uptake, growth rate and TNF-α levels to those of intact porin strains.Conclusions We hypothesize that specific mutations accumulated and maintained over time in M. massiliense, including mutations shared among transmissible strains, collectively lead to more virulent, host adapted lineages in CF patients and other susceptible hosts

Targeted enrichment outperforms other enrichment techniques and enables more multi-species RNA-Seq analyses.

Enrichment methodologies enable the analysis of minor members in multi-species transcriptomic data. We compared the standard enrichment of bacterial and eukaryotic mRNA to a targeted enrichment using an Agilent SureSelect (AgSS) capture for Brugia malayi, Aspergillus fumigatus, and the Wolbachia endosymbiont of B. malayi (wBm). Without introducing significant systematic bias, the AgSS quantitatively enriched samples, resulting in more reads mapping to the target organism. The AgSS-enriched libraries consistently had a positive linear correlation with their unenriched counterparts (r2 = 0.559-0.867). Up to a 2,242-fold enrichment of RNA from the target organism was obtained following a power law (r2 = 0.90), with the greatest fold enrichment achieved in samples with the largest ratio difference between the major and minor members. While using a single total library for prokaryote and eukaryote enrichment from a single RNA sample could be beneficial for samples where RNA is limiting, we observed a decrease in reads mapping to protein coding genes and an increase in multi-mapping reads to rRNAs in AgSS enrichments from eukaryotic total RNA libraries compared to eukaryotic poly(A)-enriched libraries. Our results support a recommendation of using AgSS targeted enrichment on poly(A)-enriched libraries for eukaryotic captures, and total RNA libraries for prokaryotic captures, to increase the robustness of multi-species transcriptomic studies

Targeted enrichment outperforms other enrichment techniques and enables more multi-species RNA-Seq analyses.

Enrichment methodologies enable the analysis of minor members in multi-species transcriptomic data. We compared the standard enrichment of bacterial and eukaryotic mRNA to a targeted enrichment using an Agilent SureSelect (AgSS) capture for Brugia malayi, Aspergillus fumigatus, and the Wolbachia endosymbiont of B. malayi (wBm). Without introducing significant systematic bias, the AgSS quantitatively enriched samples, resulting in more reads mapping to the target organism. The AgSS-enriched libraries consistently had a positive linear correlation with their unenriched counterparts (r2 = 0.559-0.867). Up to a 2,242-fold enrichment of RNA from the target organism was obtained following a power law (r2 = 0.90), with the greatest fold enrichment achieved in samples with the largest ratio difference between the major and minor members. While using a single total library for prokaryote and eukaryote enrichment from a single RNA sample could be beneficial for samples where RNA is limiting, we observed a decrease in reads mapping to protein coding genes and an increase in multi-mapping reads to rRNAs in AgSS enrichments from eukaryotic total RNA libraries compared to eukaryotic poly(A)-enriched libraries. Our results support a recommendation of using AgSS targeted enrichment on poly(A)-enriched libraries for eukaryotic captures, and total RNA libraries for prokaryotic captures, to increase the robustness of multi-species transcriptomic studies