Search

Understanding the molecular parameters that regulate cross-species transmission and host adaptation of potential pathogens is crucial to control emerging infectious disease. Although microbial pathotype diversity is conventionally associated with gene gain or loss, the role of pathoadaptive nonsynonymous single-nucleotide polymorphisms (nsSNPs) has not been systematically evaluated. Here, our genome-wide analysis of core genes within Salmonella enterica serovar Typhimurium genomes reveals a high degree of allelic variation in surface-exposed molecules, including adhesins that promote host colonization. Subsequent multinomial logistic regression, MultiPhen and Random Forest analyses of known/suspected adhesins from 580 independent Typhimurium isolates identifies distinct host-specific nsSNP signatures. Moreover, population and functional analyses of host-associated nsSNPs for FimH, the type 1 fimbrial adhesin, highlights the role of key allelic residues in host-specific adherence in vitro. Together, our data provide the first concrete evidence that functional differences between allelic variants of bacterial proteins likely contribute to pathoadaption to diverse hosts., One of the key aspects for controlling infectious diseases is understanding how pathogens cross host species. Here the authors conduct a genome-wide analysis of Salmonella and show a high degree of variation, enabling host-adapted colonization among Salmonella intestinal and systemic serovars.

Background: Samples collected from CF patient airways often contain large amounts of host-derived nucleic acids that interfere with recovery and purification of microbial and viral nucleic acids. This study describes metagenomic and metatranscriptomic methods that address these issues. Methods: Microbial and viral metagenomes, and microbial metatranscriptomes, were successfully prepared from sputum samples from five adult CF patients. Results: Contaminating host DNA was dramatically reduced in the metagenomes. Each CF patient presented a unique microbiome; in some Pseudomonas aeruginosa was replaced by other opportunistic bacteria. Even though the taxonomic composition of the microbiomes is very different, the metabolic potentials encoded by the community are very similar. The viral communities were dominated by phages that infect major CF pathogens. The metatranscriptomes reveal differential expression of encoded metabolic potential with changing health status. Conclusions: Microbial and viral metagenomics combined with microbial transcriptomics characterize the dynamic polymicrobial communities found in CF airways, revealing both the taxa present and their current metabolic activities. These approaches can facilitate the development of individualized treatment plans and novel therapeutic approaches.

The advent of high-throughput sequencing, coupled with advances in computational methods, has enabled genome-wide dissection of genetics, evolution, and disease, with nucleotide resolution. The discoveries derived from genomics promise benefits to basic research, biotechnology, and medicine; however, the speed and affordability of sequencing has resulted in a flood of “big data” in the life sciences. In addition, the current heterogeneity of sequencing platforms and diversity of applications complicate the development of tools for analysis, and this has slowed widespread adoption of the technology. Making sense of the data and delivering actionable insight requires improved computational infrastructure, new methods for interpreting the data, and unique collaborative approaches. Here we review the role of big data in genomics, its impact on the development of tools for collaborative analysis of genomes, and successes and ongoing challenges in coping with big data.

Microbial communities in the lungs of patients with cystic fibrosis (CF) and chronic obstructive pulmonary disease (COPD) have been shown to be spatially heterogeneous. Viral communities may also vary spatially, leading to localized viral populations and infections. Here, we characterized viral communities from multiple areas of the lungs of two patients with late-stage CF using metagenomics, that is, the explanted lungs from a transplant patient and lungs acquired postmortem. All regions harbored eukaryotic viruses that may infect the human host, notably herpesviruses, anelloviruses, and papillomaviruses. In the highly diseased apical lobes of explant lungs, viral diversity was extremely low, and only eukaryotic viruses were present. The absence of phage suggests that CF-associated microbial biofilms may escape top-down controls by phage predation. The phages present in other lobes of explant lungs and in all lobes of postmortem lungs comprised distinct communities, and encoded genes for clinically important microbial phenotypes, including small colony variants and antibiotic resistance. Based on the these observations, we postulate that viral communities in CF lungs are spatially distinct and contribute to CF pathology by augmenting the metabolic potential of resident microbes, as well as by directly damaging lung tissue via carcinomas and herpesviral outbreaks.

The human oropharynx is a reservoir for many potential pathogens, including streptococcal species that cause endocarditis. Although oropharyngeal microbes have been well described, viral communities are essentially uncharacterized. We conducted a metagenomic study to determine the composition of oropharyngeal DNA viral communities (both phage and eukaryotic viruses) in healthy individuals and to evaluate oropharyngeal swabs as a rapid method for viral detection. Viral DNA was extracted from 19 pooled oropharyngeal swabs and sequenced. Viral communities consisted almost exclusively of phage, and complete genomes of several phage were recovered, including Escherichia coli phage T3, Propionibacterium acnes phage PA6, and Streptococcus mitis phage SM1. Phage relative abundances changed dramatically depending on whether samples were chloroform treated or filtered to remove microbial contamination. pblA and pblB genes of phage SM1 were detected in the metagenomes. pblA and pblB mediate the attachment of S. mitis to platelets and play a significant role in S. mitis virulence in the endocardium, but have never previously been detected in the oral cavity. These genes were also identified in salivary metagenomes from three individuals at three time points and in individual saliva samples by PCR. Additionally, we demonstrate that phage SM1 can be induced by commonly ingested substances. Our results indicate that the oral cavity is a reservoir for pblA and pblB genes and for phage SM1 itself. Further studies will determine the association between pblA and pblB genes in the oral cavity and the risk of endocarditis.

We identified a novel rhabdovirus, American bat vesiculovirus, from postmortem tissue samples from 120 rabies-negative big brown bats with a history of human contact. Five percent of the tested bats were infected with this virus. The extent of zoonotic exposure and possible health effects in humans from this virus are unknown.

Understanding the molecular parameters that regulate cross-species transmission and host adaptation of potential pathogens is crucial to control emerging infectious disease. Although microbial pathotype diversity is conventionally associated with gene gain or loss, the role of pathoadaptive nonsynonymous single-nucleotide polymorphisms (nsSNPs) has not been systematically evaluated. Here, our genome-wide analysis of core genes within Salmonella enterica serovar Typhimurium genomes reveals a high degree of allelic variation in surface-exposed molecules, including adhesins that promote host colonization. Subsequent multinomial logistic regression, MultiPhen and Random Forest analyses of known/suspected adhesins from 580 independent Typhimurium isolates identifies distinct host-specific nsSNP signatures. Moreover, population and functional analyses of host-associated nsSNPs for FimH, the type 1 fimbrial adhesin, highlights the role of key allelic residues in host-specific adherence in vitro. Together, our data provide the first concrete evidence that functional differences between allelic variants of bacterial proteins likely contribute to pathoadaption to diverse hosts.

A novel targeted massive parallel sequencing approach identified genetic variation in eight known or predicted fimbrial adhesins for 46 Salmonella strains. The results highlight associations between specific adhesin alleles, host species, and antimicrobial resistance. The differentiation of allelic variants has potential applications for diagnostic microbiology and epidemiological investigations.

Nature Communications 6 Article: 8754 (2015); Published: 30 Octomber 2015; Updated: 8 August 2017 In Fig. 3 of this Article, the numbers of isolates studied for the fimH41_Newp, fimH44_Newp and fimH45_Newp alleles were inadvertently swapped. The correct version of Fig. 3 appears below as Fig. 1.

Summary: Here, we present riboPicker, a robust framework for the rapid, automated identification and removal of ribosomal RNA sequences from metatranscriptomic datasets. The results can be exported for subsequent analysis, and the databases used for the web-based version are updated on a regular basis. riboPicker categorizes rRNA-like sequences and provides graphical visualizations and tabular outputs of ribosomal coverage, alignment results and taxonomic classifications. Availability and implementation: This open-source application was implemented in Perl and can be used as stand-alone version or accessed online through a user-friendly web interface. The source code, user help and additional information is available at http://ribopicker.sourceforge.net/. Contact: rschmied@sciences.sdsu.edu; rschmied@sciences.sdsu.edu Supplementary information: Supplementary data are available at Bioinformatics online.

As DNA sequencing becomes faster and cheaper, genomics-based approaches are being explored for their use in personalized diagnoses and treatments. Here, we provide a proof of principle for disease monitoring using personal metagenomic sequencing and traditional clinical microbiology by focusing on three adults with cystic fibrosis (CF). The CF lung is a dynamic environment that hosts a complex ecosystem composed of bacteria, viruses, and fungi that can vary in space and time. Not surprisingly, the microbiome data from the induced sputum samples we collected revealed a significant amount of species diversity not seen in routine clinical laboratory cultures. The relative abundances of several species changed as clinical treatment was altered, enabling the identification of the climax and attack communities that were proposed in an earlier work. All patient microbiomes encoded a diversity of mechanisms to resist antibiotics, consistent with the characteristics of multidrug-resistant microbial communities that are commonly observed in CF patients. The metabolic potentials of these communities differed by the health status and recovery route of each patient. Thus, this pilot study provides an example of how metagenomic data might be used with clinical assessments for the development of treatments tailored to individual patients.

Recent evidence showing host specificity of colonizing bacteria supports the view that multicellular organisms are holobionts comprised of the macroscopic host in synergistic interdependence with a heterogeneous and host-specific microbial community. Whereas host-bacteria interactions have been extensively investigated, comparatively little is known about host-virus interactions and viral contribution to the holobiont. We sought to determine the viral communities associating with different Hydra species, whether these viral communities were altered with environmental stress, and whether these viruses affect the Hydra-associated holobiont. Here we show that each species of Hydra harbors a diverse host-associated virome. Primary viral families associated with Hydra are Myoviridae, Siphoviridae, Inoviridae, and Herpesviridae. Most Hydra-associated viruses are bacteriophages, a reflection of their involvement in the holobiont. Changes in environmental conditions alter the associated virome, increase viral diversity, and affect the metabolism of the holobiont. The specificity and dynamics of the virome point to potential viral involvement in regulating microbial associations in the Hydra holobiont. While viruses are generally regarded as pathogenic agents, our study suggests an evolutionary conserved ability of viruses to function as holobiont regulators and, therefore, constitutes an emerging paradigm shift in host-microbe interactions.

The impaired mucociliary clearance in individuals with Cystic Fibrosis (CF) enables opportunistic pathogens to colonize CF lungs. Here we show that Rothia mucilaginosa is a common CF opportunist that was present in 83% of our patient cohort, almost as prevalent as Pseudomonas aeruginosa (89%). Sequencing of lung microbial metagenomes identified unique R. mucilaginosa strains in each patient, presumably due to evolution within the lung. The de novo assembly of a near-complete R. mucilaginosa (CF1E) genome illuminated a number of potential physiological adaptations to the CF lung, including antibiotic resistance, utilization of extracellular lactate, and modification of the type I restriction-modification system. Metabolic characteristics predicted from the metagenomes suggested R. mucilaginosa have adapted to live within the microaerophilic surface of the mucus layer in CF lungs. The results also highlight the remarkable evolutionary and ecological similarities of many CF pathogens; further examination of these similarities has the potential to guide patient care and treatment.

Phages (viruses that infect bacteria) have gained significant attention because of their abundance, diversity and important ecological roles. However, the lack of a universal gene shared by all phages presents a challenge for phage identification and characterization, especially in environmental samples where it is difficult to culture phage-host systems. Homologous conserved genes (or "signature genes") present in groups of closely-related phages can be used to explore phage diversity and define evolutionary relationships amongst these phages. Bioinformatic approaches are needed to identify candidate signature genes and design PCR primers to amplify those genes from environmental samples; however, there is currently no existing computational tool that biologists can use for this purpose. Here we present PhiSiGns, a web-based and standalone application that performs a pairwise comparison of each gene present in user-selected phage genomes, identifies signature genes, generates alignments of these genes, and designs potential PCR primer pairs. PhiSiGns is available at ( http://www.phantome.org/phisigns/ ; http://phisigns.sourceforge.net/ ) with a link to the source code. Here we describe the specifications of PhiSiGns and demonstrate its application with a case study. PhiSiGns provides phage biologists with a user-friendly tool to identify signature genes and design PCR primers to amplify related genes from uncultured phages in environmental samples. This bioinformatics tool will facilitate the development of novel signature genes for use as molecular markers in studies of phage diversity, phylogeny, and evolution.

Motivation: Metagenomes are often characterized by high levels of unknown sequences. Reads derived from known microorganisms can easily be identified and analyzed using fast homology search algorithms and a suitable reference database, but the unknown sequences are often ignored in further analyses, biasing conclusions. Nevertheless, it is possible to use more data in a comparative metagenomic analysis by creating a cross-assembly of all reads, i.e. a single assembly of reads from different samples. Comparative metagenomics studies the interrelationships between metagenomes from different samples. Using an assembly algorithm is a fast and intuitive way to link (partially) homologous reads without requiring a database of reference sequences. Results: Here, we introduce crAss, a novel bioinformatic tool that enables fast simple analysis of cross-assembly files, yielding distances between all metagenomic sample pairs and an insightful image displaying the similarities. Availability and implementation: crAss is available as a web server at http://edwards.sdsu.edu/crass/, and the Perl source code can be downloaded to run as a stand-alone command line tool. Contact: dutilh@cmbi.ru.nl Supplementary information: Supplementary data are available at Bioinformatics online.

Cystic fibrosis (CF) is a common fatal genetic disorder with mortality most often resulting from microbial infections of the lungs. Culture-independent studies of CF-associated microbial communities have indicated that microbial diversity in the CF airways is much higher than suggested by culturing alone. However, these studies have relied on indirect methods to sample the CF lung such as expectorated sputum and bronchoalveolar lavage (BAL). Here, we characterize the diversity of microbial communities in tissue sections from anatomically distinct regions of the CF lung using barcoded 16S amplicon pyrosequencing. Microbial communities differed significantly between different areas of the lungs, and few taxa were common to microbial communities in all anatomical regions surveyed. Our results indicate that CF lung infections are not only polymicrobial, but also spatially heterogeneous suggesting that treatment regimes tailored to dominant populations in sputum or BAL samples may be ineffective against infections in some areas of the lung.

The health of the coral reefs of the Abrolhos Bank (Southwestern Atlantic) was characterized with a holistic approach using measurements of four ecosystem components: (i) inorganic and organic nutrient concentrations, [1] fish biomass, [1] macroalgal and coral cover and (iv) microbial community composition and abundance. The possible benefits of protection from fishing were particularly evaluated by comparing sites with varying levels of protection. Two reefs within the well-enforced no-take area of the National Marine Park of Abrolhos (Parcel dos Abrolhos and California) were compared with two unprotected coastal reefs (Sebastião Gomes and Pedra de Leste) and one legally protected but poorly enforced coastal reef (the "paper park" of Timbebas Reef). The fish biomass was lower and the fleshy macroalgal cover was higher in the unprotected reefs compared with the protected areas. The unprotected and protected reefs had similar seawater chemistry. Lower vibrio CFU counts were observed in the fully protected area of California Reef. Metagenome analysis showed that the unprotected reefs had a higher abundance of archaeal and viral sequences and more bacterial pathogens, while the protected reefs had a higher abundance of genes related to photosynthesis. Similar to other reef systems in the world, there was evidence that reductions in the biomass of herbivorous fishes and the consequent increase in macroalgal cover in the Abrolhos Bank may be affecting microbial diversity and abundance. Through the integration of different types of ecological data, the present study showed that protection from fishing may lead to greater reef health. The data presented herein suggest that protected coral reefs have higher microbial diversity, with the most degraded reef (Sebastião Gomes) showing a marked reduction in microbial species richness. It is concluded that ecological conditions in unprotected reefs may promote the growth and rapid evolution of opportunistic microbial pathogens.

International audience; The bacterioplankton diversity of coastal waters along a latitudinal gradient between Puerto Rico and Argentina was analyzed using a total of 134,197 high-quality sequences from the V6 hypervariable region of the small-subunit ribosomal RNA gene (16S rRNA) (mean length of 60 nt). Most of the OTUs were identified into Proteobacteria, Bacteriodetes, Cyanobacteria, and Actinobacteria, corresponding to approx. 80% of the total number of sequences. The number of OTUs corresponding to species varied between 937 and 1946 in the seven locations. Proteobacteria appeared at high frequency in the seven locations. An enrichment of Cyanobacteria was observed in Puerto Rico, whereas an enrichment of Bacteroidetes was detected in the Argentinian shelf and Uruguayan coastal lagoons. The highest number of sequences of Actinobacteria and Acidobacteria were obtained in the Amazon estuary mouth. The rarefaction curves and Good coverage estimator for species diversity suggested a significant coverage, with values ranging between 92 and 97% for Good coverage. Conserved taxa corresponded to aprox. 52% of all sequences. This study suggests that human-contaminated environments may influence bacterioplankton diversity

Background Sequencing metagenomes that were pre-amplified with primer-based methods requires the removal of the additional tag sequences from the datasets. The sequenced reads can contain deletions or insertions due to sequencing limitations, and the primer sequence may contain ambiguous bases. Furthermore, the tag sequence may be unavailable or incorrectly reported. Because of the potential for downstream inaccuracies introduced by unwanted sequence contaminations, it is important to use reliable tools for pre-processing sequence data. Results TagCleaner is a web application developed to automatically identify and remove known or unknown tag sequences allowing insertions and deletions in the dataset. TagCleaner is designed to filter the trimmed reads for duplicates, short reads, and reads with high rates of ambiguous sequences. An additional screening for and splitting of fragment-to-fragment concatenations that gave rise to artificial concatenated sequences can increase the quality of the dataset. Users may modify the different filter parameters according to their own preferences. Conclusions TagCleaner is a publicly available web application that is able to automatically detect and efficiently remove tag sequences from metagenomic datasets. It is easily configurable and provides a user-friendly interface. The interactive web interface facilitates export functionality for subsequent data processing, and is available at http://edwards.sdsu.edu/tagcleaner.

Metagenomic studies characterize both the composition and diversity of uncultured viral and microbial communities. BLAST-based comparisons have typically been used for such analyses; however, sampling biases, high percentages of unknown sequences, and the use of arbitrary thresholds to find significant similarities can decrease the accuracy and validity of estimates. Here, we present Genome relative Abundance and Average Size (GAAS), a complete software package that provides improved estimates of community composition and average genome length for metagenomes in both textual and graphical formats. GAAS implements a novel methodology to control for sampling bias via length normalization, to adjust for multiple BLAST similarities by similarity weighting, and to select significant similarities using relative alignment lengths. In benchmark tests, the GAAS method was robust to both high percentages of unknown sequences and to variations in metagenomic sequence read lengths. Re-analysis of the Sargasso Sea virome using GAAS indicated that standard methodologies for metagenomic analysis may dramatically underestimate the abundance and importance of organisms with small genomes in environmental systems. Using GAAS, we conducted a meta-analysis of microbial and viral average genome lengths in over 150 metagenomes from four biomes to determine whether genome lengths vary consistently between and within biomes, and between microbial and viral communities from the same environment. Significant differences between biomes and within aquatic sub-biomes (oceans, hypersaline systems, freshwater, and microbialites) suggested that average genome length is a fundamental property of environments driven by factors at the sub-biome level. The behavior of paired viral and microbial metagenomes from the same environment indicated that microbial and viral average genome sizes are independent of each other, but indicative of community responses to stressors and environmental conditions., Author Summary Metagenomics uses DNA or RNA sequences isolated directly from the environment to determine what viruses or microorganisms exist in natural communities and what metabolic activities they encode. Typically, metagenomic sequences are compared to annotated sequences in public databases using the BLAST search tool. Our methods, implemented in the Genome relative Abundance and Average Size (GAAS) software, improve the way BLAST searches are processed to estimate the taxonomic composition of communities and their average genome length. GAAS provides a more accurate picture of community composition by correcting for a systematic sampling bias towards larger genomes, and is useful in situations where organisms with small genomes are abundant, such as disease outbreaks caused by small RNA viruses. Microbial average genome length relates to environmental complexity and the distribution of genome lengths describes community diversity. A study of the average genome length of viruses and microorganisms in four different biomes using GAAS on 169 metagenomes showed significantly different average genome sizes between biomes, and large variability within biomes as well. This also revealed that microbial and viral average genome sizes in the same environment are independent of each other, which reflects the different ways that microorganisms and viruses respond to stress and environmental conditions.

The human respiratory tract is constantly exposed to a wide variety of viruses, microbes and inorganic particulates from environmental air, water and food. Physical characteristics of inhaled particles and airway mucosal immunity determine which viruses and microbes will persist in the airways. Here we present the first metagenomic study of DNA viral communities in the airways of diseased and non-diseased individuals. We obtained sequences from sputum DNA viral communities in 5 individuals with cystic fibrosis (CF) and 5 individuals without the disease. Overall, diversity of viruses in the airways was low, with an average richness of 175 distinct viral genotypes. The majority of viral diversity was uncharacterized. CF phage communities were highly similar to each other, whereas Non-CF individuals had more distinct phage communities, which may reflect organisms in inhaled air. CF eukaryotic viral communities were dominated by a few viruses, including human herpesviruses and retroviruses. Functional metagenomics showed that all Non-CF viromes were similar, and that CF viromes were enriched in aromatic amino acid metabolism. The CF metagenomes occupied two different metabolic states, probably reflecting different disease states. There was one outlying CF virome which was characterized by an over-representation of Guanosine-5'-triphosphate,3'-diphosphate pyrophosphatase, an enzyme involved in the bacterial stringent response. Unique environments like the CF airway can drive functional adaptations, leading to shifts in metabolic profiles. These results have important clinical implications for CF, indicating that therapeutic measures may be more effective if used to change the respiratory environment, as opposed to shifting the taxonomic composition of resident microbiota.

Summary: The affordability of high-throughput sequencing has created an unprecedented surge in the use of genomic data in basic, translational and clinical research. The rapid evolution of sequencing technology, coupled with its broad adoption across biology and medicine, necessitates fast, collaborative interdisciplinary discussion. SEQanswers provides a real-time knowledge-sharing resource to address this need, covering experimental and computational aspects of sequencing and sequence analysis. Developers of popular analysis tools are among the >4000 active members, and ~40 peer-reviewed publications have referenced SEQanswers. Availability: The SEQanswers community is freely accessible at http://SEQanswers.com/ Contact: david.mittelman@vt.edu; ecolivares@gmail.com Supplementary information: Supplementary data are available at Bioinformatics online.

Summary: Here, we present PRINSEQ for easy and rapid quality control and data preprocessing of genomic and metagenomic datasets. Summary statistics of FASTA (and QUAL) or FASTQ files are generated in tabular and graphical form and sequences can be filtered, reformatted and trimmed by a variety of options to improve downstream analysis. Availability and Implementation: This open-source application was implemented in Perl and can be used as a stand alone version or accessed online through a user-friendly web interface. The source code, user help and additional information are available at http://prinseq.sourceforge.net/. Contact: rschmied@sciences.sdsu.edu; redwards@cs.sdsu.edu

Mass culture of algae for the production of biofuels is a developing technology designed to offset the depletion of fossil fuel reserves. However, large scale culture of algae in open ponds can be challenging because of incidences of infestation with algal parasites. Without knowledge of the identity of the specific parasite and how to control these pests, algal-based biofuel production will be limited. We have characterized a eukaryotic parasite of Scenedesmus dimorphus growing in outdoor ponds used for biofuel production. We demonstrated that as the genomic DNA of parasite FD01 increases, the concentration of S. dimorphus cells decreases; consequently, this is a highly destructive pathogen. Techniques for culture of the parasite and host were developed, and the endoparasite was identified as the Aphelidea, Amoeboaphelidium protococcarum. Phylogenetic analysis of ribosomal sequences revealed that parasite FD01 placed within the recently described Cryptomycota, a poorly known phylum based on two species of Rozella and environmental samples. Transmission electron microscopy demonstrated that aplanospores of the parasite produced filose pseudopodia, which contained fine fibers the diameter of actin microfilaments. Multiple lipid globules clustered and were associated with microbodies, mitochondria and a membrane cisternae, an arrangement characteristic of the microbody-lipid globule complex of chytrid zoospores. After encystment and attachment to the host cells, the parasite injected its protoplast into the host between the host cell wall and plasma membrane. At maturity the unwalled parasite occupied the entire host cell. After cleavage of the protoplast into aplanospores, a vacuole and lipids remained in the host cell. Amoeboaphelidium protococcarum isolate FD01 is characteristic of the original description of this species and is different from strain X-5 recently characterized. Our results help put a face on the Cryptomycota, revealing that the phylum is more diverse than previously understood and include some of the Aphelidea as well as Rozella species and potentially Microsporidia.

Viruses are the most abundant and diverse genetic entities on Earth; however, broad surveys of viral diversity are hindered by the lack of a universal assay for viruses and the inability to sample a sufficient number of individual hosts. This study utilized vector-enabled metagenomics (VEM) to provide a snapshot of the diversity of DNA viruses present in three mosquito samples from San Diego, California. The majority of the sequences were novel, suggesting that the viral community in mosquitoes, as well as the animal and plant hosts they feed on, is highly diverse and largely uncharacterized. Each mosquito sample contained a distinct viral community. The mosquito viromes contained sequences related to a broad range of animal, plant, insect and bacterial viruses. Animal viruses identified included anelloviruses, circoviruses, herpesviruses, poxviruses, and papillomaviruses, which mosquitoes may have obtained from vertebrate hosts during blood feeding. Notably, sequences related to human papillomaviruses were identified in one of the mosquito samples. Sequences similar to plant viruses were identified in all mosquito viromes, which were potentially acquired through feeding on plant nectar. Numerous bacteriophages and insect viruses were also detected, including a novel densovirus likely infecting Culex erythrothorax. Through sampling insect vectors, VEM enables broad survey of viral diversity and has significantly increased our knowledge of the DNA viruses present in mosquitoes.

Genome sequencing has become routine, however genome assembly still remains a challenge despite the computational advances in the last decade. In particular, the abundance of repeat elements in genomes makes it difficult to assemble them into a single complete sequence. Identical repeats shorter than the average read length can generally be assembled without issue. However, longer repeats such as ribosomal RNA operons cannot be accurately assembled using existing tools. The application Scaffold_builder was designed to generate scaffolds – super contigs of sequences joined by N-bases – based on the similarity to a closely related reference sequence. This is independent of mate-pair information and can be used complementarily for genome assembly, e.g. when mate-pairs are not available or have already been exploited. Scaffold_builder was evaluated using simulated pyrosequencing reads of the bacterial genomes Escherichia coli 042, Lactobacillus salivarius UCC118 and Salmonella enterica subsp. enterica serovar Typhi str. P-stx-12. Moreover, we sequenced two genomes from Salmonella enterica serovar Typhimurium LT2 G455 and Salmonella enterica serovar Typhimurium SDT1291 and show that Scaffold_builder decreases the number of contig sequences by 53% while more than doubling their average length. Scaffold_builder is written in Python and is available at http://edwards.sdsu.edu/scaffold_builder. A web-based implementation is additionally provided to allow users to submit a reference genome and a set of contigs to be scaffolded.