Scripts used for characterizing human gut bacteriophages in the following manuscript:

Camarillo-Guerrero LF, Almeida A, Rangel-Pineros G, Finn RD, Lawley TD (2020) [Massive expansion of human gut bacteriophage diversity]

Associated data can also be found in our FTP server

Neural network that distinguishes phages from integrative and conjugative elements (ICEs).

Requirements:

• Python (tested v3.6.7)
• TensorFlow (tested v1.10)
• Keras (tested v.2.2.4)

Usage:

classifier.py <input_features_file.txt>

Notes:

• input_features_file.txt: It contains a feature vector of 1026 dimensions: fraction of hypothetical proteins (1), gene density (1), 5-kmer signature (1024) that represents a phage or an ICE (1 feature vector per line)
  
• classifier/classifier_demo.py: It runs a demo of the classifier with 50 examples of phages and ICEs each
  

Input features generation files:

getGeneDensity.py: This function takes in a GFF3 file and returns the number of genes / kb.
getHypothetical.py: This function takes in a GFF3 file and returns the fraction of hypothetical proteins.
getKmer.py: This function takes in a DNA sequence and counts the proportion of each of the 1024 possible 5mers.

Usage:

getGeneDensity(<gff3_file_name>)
getHypothetical(<gff3_file_name>)
getSignature_hash(<DNA_sequence>)

figures/

• 'Figure 1.py': Distribution of MIUViG scores from CheckV analysis
• 'Figure 2.py': Viral diversity patterns across gut bacteria genera and broad host range VCs
• 'Figure 3.py': Gut phageome profiling across human populations and correlation with gut bacteria enterotypes
• 'Figure 4.py': Crass-like family global distribution and host-phage network of globally distributed VCs
• 'Figure 5.py': Phylogenetic structure of the pX phage and global distribution
• 'Figure S1.py': Quality control assessment of GPD
• 'Figure S2.py': Viral diversity patterns across gut bacteria phyla and host range analysis of gut phages
• 'Figure S3.py': Correlation between sequencing depth and number of phages detected in a sample
• 'Figure S4.py': Host range analysis of globally distributed phages