Appearance
Multi-Hit Method (MULTI-HIT)
MULTI-HIT is a statistical method that allows for the examination of the fit of models accommodating multiple instantaneous nucleotide substitutions in sequence data. The method expands on traditional models by incorporating the complexities involved with multiple hits that can occur at a single nucleotide position in evolutionary processes.
Citation
If you utilize the MULTI-HIT method in your research or analysis, please cite:
Input Parameters
MULTI-HIT requires various parameters for its analysis:
Required Inputs
- Alignment File: A file containing the multiple sequence alignment (supported formats:
.fasta
,.phy
, etc.). - Genetic Code: The genetic code to use for the analysis (default: "Universal").
Optional Inputs
- Triple Islands: Toggle to specify if synonymous triple-hit substitutions should be accounted for (options: "Yes" or "No").
- Rate Classes: The number of omega rate classes to be included in the model (default: 3, range: 1-10).
Outputs
Summary
The MULTI-HIT analysis generates a JSON file that includes:
- Job Metadata: Overview of the parameters used and a timestamp of analysis completion.
- Site-Level Results: Detailed information on each site in the alignment, including statistical evidence for multiple hits and their significance.
Detailed Site Output
Each site analyzed is accompanied by:
- Evidence Ratios: Estimates comparing the occurrence of multiple hits versus other types of substitutions at that site.
- P-values: Significance levels indicating whether the observed multiple substitutions significantly deviate from the null hypothesis.
Visualization
Results from the MULTI-HIT analysis can be visualized using the integrated visualization tools available on the platform.
- Circos Plots: Interactive visual representations that map codons and their relationships based on evidence ratios.
- Data Tables: Clear, sortable tables that summarize the results for each site, allowing users to explore the analysis efficiently.
Example Workflow
Here’s a step-by-step guide to using the MULTI-HIT tool via the website:
Upload Data:
- Navigate to the MULTI-HIT page and upload your alignment file with phylogenetic tree.
- Select the appropriate genetic code from the dropdown menu.
Set Parameters:
- Adjust optional parameters such as rate classes and triple islands if needed.
Run Analysis:
- Click the "Run Analysis" button to initiate the MULTI-HIT process.
- Monitor the progress on the interface. Once completed, you will receive a notification if opted in.
Review Results:
- Access your results through a generated link or through your notifications.
- Review both the summary statistics and detailed site-level data.
- Use the visualization tools to explore your findings interactively.
Download Outputs:
- Optionally, download the output files for storage or further analysis.
Example CLI Usage
To run the FitMultiModel analysis using HyPhy, you can use the following command:
bash
/path/to/hyphy/hyphy $LIB_PATH/res/TemplateBatchFiles/SelectionMethods/FitMultiModel.bf
Parameters
- --code: Which genetic code should be used (default: "Universal").
- --alignment: An in-frame codon alignment in a supported format (required).
- --tree: A phylogenetic tree file (optional).
- --rates: The number of omega rate classes to include in the model (default: 3).
- --triple-islands: Use a separate rate parameter for synonymous triple-hit substitutions ("Yes" or "No").
- --output: File path to save the resulting JSON (default: same path as alignment file + 'FITTER.json').
- --save-fit: File path to save model fit files, with extensions (.MODEL_NAME.bf); default is NOT to save.
Full Example Usage
Here’s a complete command with parameters set:
bash
/path/to/hyphy/hyphy $LIB_PATH/res/TemplateBatchFiles/SelectionMethods/FitMultiModel.bf --code Universal --alignment alignment.fasta --tree tree.nwk --rates 3 --triple-islands Yes --output results.json --save-fit dev/null
Minimal Example Command
For a minimal example without default parameters:
bash
/path/to/hyphy/hyphy $LIB_PATH/res/TemplateBatchFiles/SelectionMethods/FitMultiModel.bf --alignment alignment.fasta