aswarren / clustertracker_va

Code for a webpage displaying introductions and resulting clusters of coronavirus across the USA. Customized for CA Dept. of Public Health.

Geek Repo:Geek Repo

Github PK Tool:Github PK Tool

Calfornia Big Tree Cluster-Tracker

This repo contains code to generate a webpage displaying SARS-CoV-2 clusters and introductions inferred via matUtils introduce.

Two versions are included here. The primary version of the California Big Tree Cluster Tracker (in the "cdph" subfolder), displays SARS-CoV-2 introductions between California counties and U.S. states. Introductions are displayed geographically on an interactive map and in a separate table below the map, with links to view and explore the data in Taxonium and the California Big Tree Investigator. A second version (in the "example" subfolder) is a generalizeable version that is set up to display introductions between U.S. states, and uses sequences from the global public SARS-CoV-2 phylogenetic tree.

The California Big Tree Cluster Tracker is itself a customized version of the UC Santa Cruz Genomics Institute's Cluster Tracker (GitHub repo here) and was developed to use data produced from the Calfornia Big Tree effort for the California Department of Public Health.

Development and Testing: Notes for development and testing of the California Big Tree Cluster Tracker can be found in the README in the "cdph" directory. The notes in that file address our particular data processing pipeline and architecture (i.e., Terra and GCP).

Customization: Information on how to create a customized version of Cluster Tracker can be found in the README in the "examples" directory. A python script and a number of data files would be required to preprocess your data if you wish to use data not in the global public SARS-CoV-2 phylogenetic tree.

Contents

Screenshot and Features

CA Big Tree Cluster Tracker screen shot

Features

  • Displays on a map the total number of introductions into a region, and introductions from one region into another.
  • Numbers of introductions can be displayed with log scaling applied or as actual numbers of introductions.
  • Introductions on the map can be filtered to show just those from the past 3 months, past 6 months, past 12 months, or the whole pandemic.
  • Users can toggle between two levels of analysis: California counties or California state. The county-level introductions display introductions to and from Calfornia counties as well as other U.S. states. The state-level introductions display introductions between all of California (as a single region) and other states.
  • A table below the map displays information about clusters in the selected region and includes:
    • the number of samples in the cluster, the range of dates for the cluster, clade and lineage information, a list of potential origins for the cluster and an index value indicating the confidence associated with the potential origin, the cluster's growth score (an importance estimate based on cluster size and age)
    • a link to view the cluster in Taxonium
    • a link to view the samples in the California Big Tree Investigator tool, where users can join PHI to the phylogentic tree information (for authorized users only)
    • all sample names in the cluster
    • for CDPH samples, a list of all associated specimen IDs
  • All data in the table can be sorted in ascending or descending order. Mulitple column sorting can be achieved by pressing the Shift key while clicking on a second column.
  • A search tool allows the user to filter the table contents.
    • Search on one or more columns for text containing a search string. Multiple phrases can be searched for using logical "and" or "or" criteria.
    • The table can also be filtered by date range, cluster size, or growth score.
  • Table data can be copied by selecting a range of table cells and using Ctrl-C to copy the contents to the computer's clipboard.
  • The Taxonium JSONL file with phylogenetic tree information and a tab-separated file with information for all introductions can be downloaded for further analysis.

Aproach

This site uses Python to perform backend setup and vanilla JavaScript for website rendering. We use the protobuf file format to store mutation-annotated phylogentic tree information, and tab-separated text files (TSV) for metadata and text-based output files. For our primary version, Python scripts are used to process a mix of samples from California State and public repositories and create a phylogenetic tree containing global and California sequences (which we call the California Big Tree). From sequence metadata we identify the region where each sample was collected (i.e., which county in California or which U.S. state) and then we calculate introductions into each region. We use the matUtils suite of tools for manipulating the protobuf files and calculating the introductions into each region. TaxoniumTools is used to convert protobuf files to JSONL for use in the Taxonium viewer. We use the Terra platform as our primary data pipeline, so our python scripts are modularized to be be compatible with WDL, but can also be run on the user's desktop computer. A GeoJSON file with introductions for each region is created for input into a Leaflet-based map. JSON files supply a table below the map with cluster details.

Quickstart

For instructions on how to set up a simplified version of Cluster Tracker using public data, see the instructions in the example directory README.

Instructions for setting up the Calfornia Big Tree Cluster Tracker on UCSC's infrastructure can be found in the "cdph" README and notes in the Pathogen Genomics Google Drive.

Software Requirements: For either the CDPH or basic version of Cluster Tracker, you will need to have the UShER software suite (version 0.6.2 or later) and TaxoniumTools installed and available on your path. You should have python 3 installed on your computer. Verify that the python "dateutil" package is included as some versions of python may be missing the dateutil standard package; it can be installed if needed via conda.

Further Details

Overview of Repo Folder Structure

  • cdph: Contains source code for data processing and the web app for the California Big Tree Cluster Tracker.
    • data: Contains python source code for data processing as well as extra files needed during creation and transfer of output files to the Google Cloud storage bucket, and a copy of the Taxonium config file for customizing the Taxonium interface.
    • www: Contains files (except for the data files) to be deployed via Google App Engine for the CA Big Tree Cluster Tracker web site. (Note: data files are hosted via GCP.)
    • README.md: Contains data processing and deployment notes for the California Big Tree Cluster Tracker.
  • example: Contains source code for data processing and the web app for creating a basic implementation of Cluster Tracker. Note that this basic version simplifies the California Big Tree Cluster Tracker by eliminating the option for California-county-based introductions.
    • data: Contains python source code for data processing using the global public SARS-CoV-2 phylogenetic tree and creating US-state-based introductions.
    • www: Contains all files needed for the web app except for data files and a customized Taxonium backend.
    • README.md: Contains notes on how to customize Cluster Tracker for a different geographic region.
  • src/js: Contains source code for the JS libraries used in the California Big Tree Cluster Tracker app.

About

Code for a webpage displaying introductions and resulting clusters of coronavirus across the USA. Customized for CA Dept. of Public Health.


Languages

Language:JavaScript 30.9%Language:Python 26.0%Language:CSS 22.0%Language:HTML 21.2%