Metagenomic Analysis of Respiratory Tract DNA Viral Communities in Cystic Fibrosis and Non-Cystic Fibrosis Individuals

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

Assessing Coral Reefs on a Pacific-Wide Scale Using the Microbialization Score

The majority of the world’s coral reefs are in various stages of decline. While a suite of disturbances (overfishing, eutrophication, and global climate change) have been identified, the mechanism(s) of reef system decline remain elusive. Increased microbial and viral loading with higher percentages of opportunistic and specific microbial pathogens have been identified as potentially unifying features of coral reefs in decline. Due to their relative size and high per cell activity, a small change in microbial biomass may signal a large reallocation of available energy in an ecosystem; that is the microbialization of the coral reef. Our hypothesis was that human activities alter the energy budget of the reef system, specifically by altering the allocation of metabolic energy between microbes and macrobes. To determine if this is occurring on a regional scale, we calculated the basal metabolic rates for the fish and microbial communities at 99 sites on twenty-nine coral islands throughout the Pacific Ocean using previously established scaling relationships. From these metabolic rate predictions, we derived a new metric for assessing and comparing reef health called the microbialization score. The microbialization score represents the percentage of the combined fish and microbial predicted metabolic rate that is microbial. Our result

Confirmatory sequencing of full-genome amplicons generated by PCR.

<p>A. For each amplicons, the sequencing was carried out with primers other than those used for full-genome amplification. The linear sequencing template (black line), the PCR primers are shown at both ends of the amplicons. Bleu lines (amplicons 671F-892R), green lines (amplicons 6838F-6972R), red lines (amplicons 6838F-7076R) represent the sequenced regions attached to the corresponding sequencing primers. Numbers included in primer IDs represent the position relative to HPV_SD2 genome. B. Alignment of sequences generated by Sanger method along the HPV_SD2 generated by 454 sequencing. Bars in blue, green and red correspond to sequences shown in <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0058404#pone-0058404-g006" target="_blank">Figure 6A</a>. The newly generated sequences covered 4,622 nucleotides representing 63.3% of the total HPV_SD2 genome. C. Comparison of the sequences generated by Sanger method and HPV_SD2 contig generated by 454 sequencing. Size of each Sanger sequence and percentage identity are shown.</p

Composition of sample pools and methods for DNA and RNA extraction.

<p>Composition of sample pools and methods for DNA and RNA extraction.</p

Summary of known viruses detected in nasopharyngeal samples.

<p>Summary of known viruses detected in nasopharyngeal samples.</p

Genomic organization of the HPV_SD2 virus.

<p>Open reading frames (L2, L1, E6, E7, E1, E2, E4) and a 530 bp non-coding long control region (LCR) are shown. B. Details of the LCR region showing the TATA box (TATAAA, positions 3735–3740), a polyadenylated site (AATAAA, positions 3336–3341), 3 palindrome sites (ACCG-N₄-CGGT; positions 3483–3494, 3650–3661, 3720–3731) and 1 degenerate palindrome (ACC-N₆-GGT, positions 3524–3535). C. Metal-binding domains in deduced E6 and E7 proteins.</p

Confirmatory PCRs to verify HPV_SD2 is a circular double-stranded genome.

<p>A. Graphical representation of the binding sites of primers 671F-892R on the putative circular structure of HPV_SD2, and the predicted PCR product sizes (I: 222 bp, II: 7,591 bp, III: 14,890 bp). The predicted short band (I) indicates the amplification of the proximal region between primers 671F and 892R. The large band (II) indicates that <i>Taq</i> DNA polymerase would amplify the region between the primers 671F and 892R by making the full circle of the HPV genome. The large band (III) indicates the <i>Taq</i> DNA polymerase would make 2 full circle around the HPV genome. PCRs were also performed using primer sets 6838F-6972R, 6838F-7076R. B. Agarose gel (0.5%) showing the amplified HPV_SD2. Primer sets used are shown: primer sets 6838F-6972R, 6838F-7076R. For each PCR, the same sample pool was tested at different concentrations [1∶1 (lanes 1, 4, 7), 1∶10 (lanes 2, 5, 8) and 1∶50 (lanes (3, 6, 9)] using 1 µl per reaction. L: DNA ladder. Amplicons can be seen at the expected sizes.</p