Package: glypniro
GlypNirOA automated script for processing and combining Byonic and PD standard output.RequirementsPythons 3.9+ with the following packages installed.pandassequaluniprotparserrequestsxlrdopenpyxlscipyInstallationThe script can be installed using the following command.pip install glypniroUsageThe program can be run from the location of the script using the commandglypniro. The following parameters can be used for operating the script within the commandline.ParemetersDescriptions-i,--input-fileFilepath to an xlsx file describing the experiment.-o,--outputFilepath to the output xlsx file for the analysis.-s,--score-cutoff(Optional) Default=200. Cutoff score for filtering of Byonic output-t,--trust-byonic(Optional) Instruct the script to trust glycan position assignment and used them for area under the curve calculation.-d,--debug(Optional) In conjunction to the final output, the script would also create debug files that contain the unique PSM selected for calculation of the data in the final output.-p,--parse-uniprot(Optional) Attempt to parse UniProt accession ID using regular expression and use them as master id.-g,--get-uniprot(Optional) Using theuniprotparsermodule to access and parse protein name from uniprot accession ids of the proteins of those within the dataset.Input file formatThe input file used in the-iparameter should have the following format.Ex:condition_idreplicate_idfilenamearea_filenameDepleted_Plasma_HCC1\data\NoMutliFuc_Nlink_10_20ppm_05Da_Depleted_Plasma_HCC_1.raw_Byonic.xlsx\data\Depleted_Plasma_HCC_1_MSnSpectrumInfo.txtDepleted_Plasma_HCC2\data\NoMutliFuc_Nlink_10_20ppm_05Da_Depleted_Plasma_HCC_2.raw_Byonic.xlsx\data\Depleted_Plasma_HCC_2_MSnSpectrumInfo.txtDepleted_Plasma_HCC3\data\NoMutliFuc_Nlink_10_20ppm_05Da_Depleted_Plasma_HCC_3.raw_Byonic.xlsx\data\Depleted_Plasma_HCC_3_MSnSpectrumInfo.txtDepleted_Plasma_Nor1\data\NoMutliFuc_Nlink_10_20ppm_05Da_Depleted_Plasma_Nor_1.raw_Byonic.xlsx\data\Depleted_Plasma_Nor_1_MSnSpectrumInfo.txtDepleted_Plasma_Nor2\data\NoMutliFuc_Nlink_10_20ppm_05Da_Depleted_Plasma_Nor_2.raw_Byonic.xlsx\data\Depleted_Plasma_Nor_2_MSnSpectrumInfo.txtDepleted_Plasma_Nor3\data\NoMutliFuc_Nlink_10_20ppm_05Da_Depleted_Plasma_Nor_3.raw_Byonic.xlsx\data\Depleted_Plasma_Nor_3_MSnSpectrumInfo.txtNon_Depleted_Plasma_HCC1\data\NoMutliFuc_Nlink_10_20ppm_05Da_Non_Depleted_Plasma_HCC_1.raw_Byonic.xlsx\data\Non_Depleted_Plasma_HCC_1_MSnSpectrumInfo.txtNon_Depleted_Plasma_HCC2\data\NoMutliFuc_Nlink_10_20ppm_05Da_Non_Depleted_Plasma_HCC_2.raw_Byonic.xlsx\data\Non_Depleted_Plasma_HCC_2_MSnSpectrumInfo.txtNon_Depleted_Plasma_HCC3\data\NoMutliFuc_Nlink_10_20ppm_05Da_Non_Depleted_Plasma_HCC_3.raw_Byonic.xlsx\data\Non_Depleted_Plasma_HCC_3_MSnSpectrumInfo.txtNon_Depleted_Plasma_Nor1\data\NoMutliFuc_Nlink_10_20ppm_05Da_Non_Depleted_Plasma_Nor_1.raw_Byonic.xlsx\data\Non_Depleted_Plasma_Nor_1_MSnSpectrumInfo.txtNon_Depleted_Plasma_Nor2\data\NoMutliFuc_Nlink_10_20ppm_05Da_Non_Depleted_Plasma_Nor_2.raw_Byonic.xlsx\data\Non_Depleted_Plasma_Nor_2_MSnSpectrumInfo.txtNon_Depleted_Plasma_Nor3\data\NoMutliFuc_Nlink_10_20ppm_05Da_Non_Depleted_Plasma_Nor_3.raw_Byonic.xlsx\data\Non_Depleted_Plasma_Nor_3_MSnSpectrumInfo.txtThe rows within the input file should not have to follow any particular order however the columns have to contain the necessary content:condition_idcondition label of the experimentreplicate_idreplicate label for the experimentfilenamefilepath to the Byonic .xlsx output file of the experiment.area_filenamefilepath to the PD tabulated output file of the experiment.Output for default execution modeThe output of the script is a single .xlsx file with 4 sheets.The first one is the output where we calculate the proportion by combining both the glycosylated and unglycosylated glycoform data.The second one is the output where the proportion of glycosylated was calculated without the unglycosylated data while the unglycosylated data was calculated similar to the first sheet.The third one is the only the filter of the glycosylated output from the second sheet.The forth one is the filter for only of the unglycosylated output from the first sheet.Exampleglypniro -i test_experiment.xlsx -o test_output.xlsx -t -gThe above command would instruct the script to use thetest_experiment.xlsxfile as input file and output astest_output.xlsx.
Inclusion of-twould mean that we trust Byonic assignment of glycan position and shall use them for calculation of that specific glycoform AUC within the proteins.
Inclusion of-gwould instruct the script to connect to the UniProt online database and attempt to parse protein name from the UniProt accession id contain in the protein name within the Byonic file.glypniro -i test_experiment.xlsx -o test_output.xlsxThe above command would instruct the script to use thetest_experiment.xlsxfile as input file and output astest_output.xlsx.
Without-toptional parameter, we only use the information of what glycans were found but not assigning them any positions. The AUC will only be combined for those PSMs with the same peptide sequence and glycan combination.Note for using with PeakView SWATH dataThe commandglypniro-reformatshould be used to generate appropriate input for GlypNirO from SWATH and Byonic data.glypniro-reformat -b byonic.xlsx -p peakview_peptide.xlsx -o description_peakview.xlsxThe command above will generate input files from experiment information and peptides information from GlypNirO and SWATH and create a description file that can be directly input into the main GlypNirO script.
For optimal result, the SWATH library used for the PeakView should be constructed from Byonic identification data.After that, theglypnirocommand can be used to process the data with the additional argument-s 0to instruct the script to ignore the Byonic score cutoff since the combining process does not parse Byonic score and only substitute with value 1.Note for using with Skyline dataThe commandglypniro-reformatsshould be used to generate appropriate input for GlypNirO from Skyline and Byonic data.glypniro-reformats -b byonic.xlsx -s skyline.csv -o description_skyline.xlsxThe command above will generate input files from experiment information and peptides information from GlypNirO and Skyline and create a description file that can be directly input into the main GlypNirO script.