The Effect of Mutation Subtypes on the Allele Frequency Spectrum and Population Genetics Inference

Analysis for paper: Effect of Mutation Subtypes on the Allele Frequency Spectrum and Population Genetics Inference

Step 1) Use tar -xvzf afs_analysis.tar to unzip file. Major files in directories include:

a) ./data

• ./data/mst_sfs/: Text file for each of 96 subtypes. File all_subtypes_afs.txt contains all 96 AFS in one file
• ./gw_sites_info.txt: Text file containing each variant info genome wide used in analysis
• ./all_MST_100kb_counts.txt: Text file containing data for 100Kb window analysis (in each window: D, mst counts, etc)
• ./data/DaDi/bootstrap_SFS: 100 bootstrapped AFS for each subtype
• ./data/simulated_neutral/all_*.txt: Simulated neutral AFS using theta estimates for A_C.AAA and C_T.ACG

b) ./output

• sigma_mat.txt: Large sigma matrix used for covariance calculations in D-2 statistic
• ./output/D2_stat/by_subtype: Output of D-2 statistic for each subtypes AFS
• ./output/DaDi: Output from running dadi 10 times on each subtypes AFS.

c) ./to_submit_scripts

• entire_analysis.R: R script containing entire pipeline sourcing each relevant analysis file
• afs_paper_functions.R: R script with functions used

Step 2) Install software needed to run all analysis completely from scratch.

• R libraries needed can be found in ./scripts/libraries_needed.R
• DaDi - python 3 version (https://dadi.readthedocs.io/en/latest/user-guide/installation/)
• SLURM like parallel computing envrionment

Step 3) Run analysis by going through entire_analysis.R in sections.

• Create AFS for each of 96 mutation subtypes. Outputs separate txt file for each and a single with all 96 subtypes.
• Create dataframe of genome wide AFS statistics for each of 96 subtypes. Includes D, D-2, % singles, doubles, triples, etc
• AFS heterogeneity across subtypes and signals of recurrent mutations and gene conversion shaping the AFS
  • Signals of mutation rate/recurrent mutations on AFS. Look at singleton/doubleton ratio by mutation rate
  • Signals of gBGC on AFS. Split into WS, SW, and indifferent mutation types. Compare D and D-2 across subtypes
• DaDi Analysis: Effect of genome wide AFS heterogeneity across subtypes on demographic inference
• Regional 100Kb Analysis: Effect of subtype composition and local genomic factors on the regional AFS