Simultaneous rapid sequencing of multiple RNA virus genomes

AbstractComparing sequences of archived viruses collected over many years to the present allows the study of viral evolution and contributes to the design of new vaccines. However, the difficulty, time and expense of generating full-length sequences individually from each archived sample have hampered these studies. Next generation sequencing technologies have been utilized for analysis of clinical and environmental samples to identify viral pathogens that may be present. This has led to the discovery of many new, uncharacterized viruses from a number of viral families. Use of these sequencing technologies would be advantageous in examining viral evolution. In this study, a sequencing procedure was used to sequence simultaneously and rapidly multiple archived samples using a single standard protocol. This procedure utilized primers composed of 20 bases of known sequence with 8 random bases at the 3′-end that also served as an identifying barcode that allowed the differentiation each viral library following pooling and sequencing. This conferred sequence independence by random priming both first and second strand cDNA synthesis. Viral stocks were treated with a nuclease cocktail to reduce the presence of host nucleic acids. Viral RNA was extracted, followed by single tube random-primed double-stranded cDNA synthesis. The resultant cDNAs were amplified by primer-specific PCR, pooled, size fractionated and sequenced on the Ion Torrent PGM platform. The individual virus genomes were readily assembled by both de novo and template-assisted assembly methods. This procedure consistently resulted in near full length, if not full-length, genomic sequences and was used to sequence multiple bovine pestivirus and coronavirus isolates simultaneously

Draft Genome Sequence of a Mycobacterium avium Complex Isolate from a Broadbill Bird

Mycobacterium avium complex (MAC) organisms cause opportunistic infections in humans, yet their epidemiology remains poorly understood. They are slowly growing environmental and animal-associated mycobacteria that have little notoriety except for the strains that cause disseminated infections in HIV- infected humans (1). Most MAC organisms are classified taxonomically as a single species, M. avium, which is divided into at least four subspecies, M. avium subsp. avium, M. avium subsp. hominissuis, M. avium subsp. paratuberculosis, and M. avium subsp. silvaticum (2). The only other species in this group is M. intracellulare. Genotyping of this diverse bacterial group has been achieved using intergenic spacers (3) and rpoB sequence analysis (4, 5)

Estimation of viral richness from shotgun metagenomes using a frequency count approach

BACKGROUND: Viruses are important drivers of ecosystem functions, yet little is known about the vast majority of viruses. Viral shotgun metagenomics enables the investigation of broad ecological questions in phage communities. One ecological characteristic is species richness, which is the number of different species in a community. Viruses do not have a phylogenetic marker analogous to the bacterial 16S rRNA gene with which to estimate richness, and so contig spectra are employed to measure the number of virus taxa in a given community. A contig spectrum is generated from a viral shotgun metagenome by assembling the random sequence reads into groups of sequences that overlap (contigs) and counting the number of sequences that group within each contig. Current tools available to analyze contig spectra to estimate phage richness are limited by relying on rank-abundance data. RESULTS: We present statistical estimates of virus richness from contig spectra. The program CatchAll (http://www.northeastern.edu/catchall/) was used to analyze contig spectra in terms of frequency count data rather than rank-abundance, thus enabling formal statistical analyses. Also, the influence of potentially spurious low-frequency counts on richness estimates was minimized by two methods, empirical and statistical. The results show greater estimates of viral richness than previous calculations in nearly all environments analyzed, including swine feces and reclaimed fresh water. CONCLUSIONS: CatchAll yielded consistent estimates of richness across viral metagenomes from the same or similar environments. Additionally, analysis of pooled viral metagenomes from different environments via mixed contig spectra resulted in greater richness estimates than those of the component metagenomes. Using CatchAll to analyze contig spectra will improve estimations of richness from viral shotgun metagenomes, particularly from large datasets, by providing statistical measures of richness

Early antibody response against Mycobacterium avium subspecies paratuberculosis antigens in subclinical cattle

<p>Abstract</p> <p>Background</p> <p>Our laboratories have previously reported on the experimental infection of cattle with <it>Mycobacterium avium </it>subsp <it>paratuberculosis </it>(<it>M. paratuberculosis</it>) using an intratonsillar infection model. In addition, we have recently developed a partial protein array representing 92 <it>M. paratuberculosis </it>coding sequences. These combined tools have enabled a unique look at the temporal analysis of <it>M. paratuberculosis </it>antigens within the native host. The primary objective of this study was to identify <it>M. paratuberculosis </it>antigens detected by cattle early during infection. A secondary objective was to evaluate the humoral immune response in cattle during the initial year of infection.</p> <p>Results</p> <p>Sera from two experimentally infected cattle, taken pre-inoculation and at day 70, 194 and 321 post infection, identified dynamic antibody reactivity among antigens with some showing an increased response over time and others showing declining levels of reactivity over the same time period. A <it>M. paratuberculosis </it>specific protein, encoded by MAP0862, was strongly detected initially, but the antibody response became weaker with time. The most reactive protein was a putative surface antigen encoded by MAP1087. A second protein, MAP1204, implicated in virulence, was also strongly detected by day 70 in both cattle. Subsequent experiments showed that these two proteins were detected with sera from 5 of 9 naturally infected cattle in the subclinical stage of Johne's disease.</p> <p>Conclusion</p> <p>Collectively these results demonstrate that <it>M. paratuberculosis </it>proteins are detected by sera from experimentally infected cattle as early as 70 days after exposure. These data further suggest at least two antigens may be useful in the early diagnosis of <it>M. paratuberculosis </it>infections. Finally, the construction and use of a protein array in this pilot study has led to a novel approach for discovery of <it>M. paratuberculosis </it>antigens.</p

The Application and Performance of Single Nucleotide Polymorphism Markers for Population Genetic Analyses of Lepidoptera

Microsatellite markers are difficult to apply within lepidopteran studies due to the lack of locus-specific PCR amplification and the high proportion of “null” alleles, such that erroneous estimations of population genetic parameters often result. Herein single nucleotide polymorphism (SNP) markers are developed from Ostrinia nubilalis (Lepidoptera: Crambidae) using next generation expressed sequence tag (EST) data. A total of 2742 SNPs were predicted within a reference assembly of 7414 EST contigs, and a subset of 763 were incorporated into 24 multiplex PCR reactions. To validate this pipeline, 5 European and North American sample sites were genotyped at 178 SNP loci, which indicated 84 (47.2%) were in Hardy–Weinberg equilibrium. Locus-by-locus FST, analysis of molecular variance, and STRUCTURE analyses indicate significant genetic differentiation may exist between European and North American O. nubilalis. The observed genetic diversity was significantly lower among European sites, which may result from genetic drift, natural selection, a genetic bottleneck, or ascertainment bias due to North American origin of EST sequence data. SNPs are an abundant source of mutation data for molecular genetic marker development in non-model species, with shared ancestral SNPs showing application within closely related species. These markers offer advantages over microsatellite markers for genetic and genomic analyses of Lepidoptera, but the source of mutation data may affect the estimation of population parameters and likely need to be considered in the interpretation of empirical data

The In-Feed Antibiotic Carbadox Induces Phage Gene Transcription in the Swine Gut Microbiome

Carbadox is a quinoxaline-di-N-oxide antibiotic fed to over 40% of young pigs in the United States that has been shown to induce phage DNA transduction in vitro; however, the effects of carbadox on swine microbiome functions are poorly understood. We investigated the in vivo longitudinal effects of carbadox on swine gut microbial gene expression (fecal metatranscriptome) and phage population dynamics (fecal dsDNA viromes). Microbial metagenome, transcriptome, and virome sequences were annotated for taxonomic inference and gene function by using FIGfam (isofunctional homolog sequences) and SEED subsystems databases. When the beta diversities of microbial FIGfam annotations were compared, the control and carbadox communities were distinct 2 days after carbadox introduction. This effect was driven by carbadox-associated lower expression of FIGfams (n = 66) related to microbial respiration, carbohydrate utilization, and RNA metabolism (q \u3c 0.1), suggesting bacteriostatic or bactericidal effects within certain populations. Interestingly, carbadox treatment caused greater expression of FIGfams related to all stages of the phage lytic cycle 2 days following the introduction of carbadox (q ≤0.07), suggesting the carbadox-mediated induction of prophages and phage DNA recombination. These effects were diminished by 7 days of continuous carbadox in the feed, suggesting an acute impact. Additionally, the viromes included a few genes that encoded resistance to tetracycline, aminoglycoside, and beta-lactam antibiotics but these did not change in frequency over time or with treatment. The results show decreased bacterial growth and metabolism, prophage induction, and potential transduction of bacterial fitness genes in swine gut bacterial communities as a result of carbadox administration

Inactivation of the Pta-AckA Pathway Causes Cell Death in \u3ci\u3eStaphylococcus aureus\u3c/i\u3e

During growth under conditions of glucose and oxygen excess, Staphylococcus aureus predominantly accumulates acetate in the culture medium, suggesting that the phosphotransacetylase-acetate kinase (Pta-AckA) pathway plays a crucial role in bacterial fitness. Previous studies demonstrated that these conditions also induce the S. aureus CidR regulon involved in the control of cell death. Interestingly, the CidR regulon is comprised of only two operons, both encoding pyruvate catabolic enzymes, suggesting an intimate relationship between pyruvate metabolism and cell death. To examine this relationship, we introduced ackA and pta mutations in S. aureus and tested their effects on bacterial growth, carbon and energy metabolism, cid expression, and cell death. Inactivation of the Pta-AckA pathway showed a drastic inhibitory effect on growth and caused accumulation of dead cells in both pta and ackA mutants. Surprisingly, inactivation of the Pta-AckA pathway did not lead to a decrease in the energy status of bacteria, as the intracellular concentrations of ATP, NAD+, and NADH were higher in the mutants. However, inactivation of this pathway increased the rate of glucose consumption, led to a metabolic block at the pyruvate node, and enhanced carbon flux through both glycolysis and the tricarboxylic acid (TCA) cycle. Intriguingly, disruption of the Pta-AckA pathway also induced the CidR regulon, suggesting that activation of alternative pyruvate catabolic pathways could be an important survival strategy for the mutants. Collectively, the results of this study demonstrate the indispensable role of the Pta-AckA pathway in S. aureus for maintaining energy and metabolic homeostasis during overflow metabolism

Genome sequencing of ovine isolates of Mycobacterium avium subspecies paratuberculosis offers insights into host association

<p>Abstract</p> <p>Background</p> <p>The genome of <it>Mycobacterium avium </it>subspecies <it>paratuberculosis </it>(<it>MAP</it>) is remarkably homogeneous among the genomes of bovine, human and wildlife isolates. However, previous work in our laboratories with the bovine K-10 strain has revealed substantial differences compared to sheep isolates. To systematically characterize all genomic differences that may be associated with the specific hosts, we sequenced the genomes of three U.S. sheep isolates and also obtained an optical map.</p> <p>Results</p> <p>Our analysis of one of the isolates, <it>MAP </it>S397, revealed a genome 4.8 Mb in size with 4,700 open reading frames (ORFs). Comparative analysis of the <it>MAP </it>S397 isolate showed it acquired approximately 10 large sequence regions that are shared with the human <it>M. avium </it>subsp. <it>hominissuis </it>strain 104 and lost 2 large regions that are present in the bovine strain. In addition, optical mapping defined the presence of 7 large inversions between the bovine and ovine genomes (~ 2.36 Mb). Whole-genome sequencing of 2 additional sheep strains of <it>MAP </it>(JTC1074 and JTC7565) further confirmed genomic homogeneity of the sheep isolates despite the presence of polymorphisms on the nucleotide level.</p> <p>Conclusions</p> <p>Comparative sequence analysis employed here provided a better understanding of the host association, evolution of members of the <it>M. avium </it>complex and could help in deciphering the phenotypic differences observed among sheep and cattle strains of <it>MAP</it>. A similar approach based on whole-genome sequencing combined with optical mapping could be employed to examine closely related pathogens. We propose an evolutionary scenario for <it>M. avium </it>complex strains based on these genome sequences.</p