unable to run Edit_dbSNP_Files function

[laurentierney]

Hi there,

I am unable to run this Edit_dbSNP_Files function.

I get this error : Error in file(file, "rt") : cannot open the connection

I downloaded a gft file per chromosome and they are stored in a directory called eSNP. The file path to an individual gtf files is like the following, with just the number on the end ranging from 1-24:

/Users/laurentierney/Desktop/eSNP/dbSNP_database_snp151_UCSC_01.gtf

this is the code I'm trying to run: edited <- Edit_dbSNP_Files("/Users/laurentierney/Desktop/eSNP", File_Name = "dbSNP_database_snp151_UCSC", Organism = "Human")

I have noticed when the error message comes back it is adding the number "1" to the end of this file path:

cannot open file '/Users/laurentierney/Desktop/eSNPdbSNP_database_snp151_UCSC1': No such file or directoryError in file(file, "rt") : cannot open the connection

Sorry if this is very basic but I can not understand what I am doing wrong

thanks

Lauren

[BenvenLab]

Hi Lauren, I'm Attaching a much more detailed manual for the eSNP-Kryotyping. Try following the steps there and let me know if the problem is still occurring.

Thanks, Uri

‫בתאריך יום ד׳, 17 במרץ 2021 ב-8:09 מאת ‪Lauren‬‏ @.*** ‬‏>:‬

Hi there,

I am unable to run this Edit_dbSNP_Files function.

I get this error : Error in file(file, "rt") : cannot open the connection

I downloaded a gft file per chromosome and they are stored in a directory called eSNP. The file path to an individual gtf files is like the following, with just the number on the end ranging from 1-24:

/Users/laurentierney/Desktop/eSNP/dbSNP_database_snp151_UCSC_01.gtf

this is the code I'm trying to run: edited <- Edit_dbSNP_Files("/Users/laurentierney/Desktop/eSNP", File_Name = "dbSNP_database_snp151_UCSC", Organism = "Human")

I have noticed when the error message comes back it is adding the number "1" to the end of this file path:

cannot open file '/Users/laurentierney/Desktop/eSNPdbSNP_database_snp151_UCSC1': No such file or directoryError in file(file, "rt") : cannot open the connection

Sorry if this is very basic but I can not understand what I am doing wrong

thanks

Lauren

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/BenvenLab/eSNPKaryotyping/issues/4, or unsubscribe https://github.com/notifications/unsubscribe-auth/AD7XRUTI4DTP7WL65VLXARLTECLZPANCNFSM4ZKO2B2A .

[laurentierney]

Hi Uri,

thank you so much for getting back to me, but a colleague had figured out the issues, it was because my db files I had downloaded file extensions.

thank you,

Lauren

---

From: BenvenLab @.> Sent: 22 March 2021 13:38 To: BenvenLab/eSNPKaryotyping @.> Cc: Lauren Tierney @.>; Manual @.> Subject: Re: [BenvenLab/eSNPKaryotyping] unable to run Edit_dbSNP_Files function (#4)

Hi Lauren, I'm Attaching a much more detailed manual for the eSNP-Kryotyping. Try following the steps there and let me know if the problem is still occurring.

Thanks, Uri

‫בתאריך יום ד׳, 17 במרץ 2021 ב-8:09 מאת ‪Lauren‬‏ @.*** ‬‏>:‬

Hi there,

I am unable to run this Edit_dbSNP_Files function.

I get this error : Error in file(file, "rt") : cannot open the connection

I downloaded a gft file per chromosome and they are stored in a directory called eSNP. The file path to an individual gtf files is like the following, with just the number on the end ranging from 1-24:

/Users/laurentierney/Desktop/eSNP/dbSNP_database_snp151_UCSC_01.gtf

this is the code I'm trying to run: edited <- Edit_dbSNP_Files("/Users/laurentierney/Desktop/eSNP", File_Name = "dbSNP_database_snp151_UCSC", Organism = "Human")

I have noticed when the error message comes back it is adding the number "1" to the end of this file path:

cannot open file '/Users/laurentierney/Desktop/eSNPdbSNP_database_snp151_UCSC1': No such file or directoryError in file(file, "rt") : cannot open the connection

Sorry if this is very basic but I can not understand what I am doing wrong

thanks

Lauren

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/BenvenLab/eSNPKaryotyping/issues/4, or unsubscribe https://github.com/notifications/unsubscribe-auth/AD7XRUTI4DTP7WL65VLXARLTECLZPANCNFSM4ZKO2B2A .

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHubhttps://github.com/BenvenLab/eSNPKaryotyping/issues/4#issuecomment-804067882, or unsubscribehttps://github.com/notifications/unsubscribe-auth/AO3JUULHSWCLE7ZX63IL7SDTE5B5NANCNFSM4ZKO2B2A.

[BenvenLab]

Great! I'm glad to hear that! I hope this tool will be helpful to you! Uri

‫בתאריך יום ב׳, 22 במרץ 2021 ב-10:51 מאת ‪Lauren‬‏ <‪ @.***‬‏>:‬

Hi Uri,

thank you so much for getting back to me, but a colleague has just this minute figured out the issues, it was because my files had file extensions.

I will make a comment on your github page to explain my error.

thank you,

Lauren

---

From: BenvenLab @.> Sent: 22 March 2021 13:38 To: BenvenLab/eSNPKaryotyping @.> Cc: Lauren Tierney @.>; Manual @.> Subject: Re: [BenvenLab/eSNPKaryotyping] unable to run Edit_dbSNP_Files function (#4)

Hi Lauren, I'm Attaching a much more detailed manual for the eSNP-Kryotyping. Try following the steps there and let me know if the problem is still occurring.

Thanks, Uri

‫בתאריך יום ד׳, 17 במרץ 2021 ב-8:09 מאת ‪Lauren‬‏ @.*** ‬‏>:‬

Hi there,

I am unable to run this Edit_dbSNP_Files function.

I get this error : Error in file(file, "rt") : cannot open the connection

I downloaded a gft file per chromosome and they are stored in a directory called eSNP. The file path to an individual gtf files is like the following, with just the number on the end ranging from 1-24:

/Users/laurentierney/Desktop/eSNP/dbSNP_database_snp151_UCSC_01.gtf

this is the code I'm trying to run: edited <- Edit_dbSNP_Files("/Users/laurentierney/Desktop/eSNP", File_Name = "dbSNP_database_snp151_UCSC", Organism = "Human")

I have noticed when the error message comes back it is adding the number "1" to the end of this file path:

cannot open file '/Users/laurentierney/Desktop/eSNPdbSNP_database_snp151_UCSC1': No such file or directoryError in file(file, "rt") : cannot open the connection

Sorry if this is very basic but I can not understand what I am doing wrong

thanks

Lauren

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/BenvenLab/eSNPKaryotyping/issues/4, or unsubscribe < https://github.com/notifications/unsubscribe-auth/AD7XRUTI4DTP7WL65VLXARLTECLZPANCNFSM4ZKO2B2A>

.

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub< https://github.com/BenvenLab/eSNPKaryotyping/issues/4#issuecomment-804067882>, or unsubscribe< https://github.com/notifications/unsubscribe-auth/AO3JUULHSWCLE7ZX63IL7SDTE5B5NANCNFSM4ZKO2B2A>.

— You are receiving this because you commented. Reply to this email directly, view it on GitHub https://github.com/BenvenLab/eSNPKaryotyping/issues/4#issuecomment-804123419, or unsubscribe https://github.com/notifications/unsubscribe-auth/AD7XRUTB6STNOUBUEOLJLATTE5KPNANCNFSM4ZKO2B2A .

[laurentierney]

Hi Uri,

I am really keen to apply your eSNP karyotyping protocol to my research however, I am struggling with section 3.2.2

After running cufflinks, I'm struggling with how to produce the merged table you show in the manuscript, with Gene_id, Gene_name, Chromosome number, gene start position and then FPKM for each of the samples. How did you go from running cufflinks to producing this table?

Any help would be greatly appreciated,

thank you,

Lauren

---

From: BenvenLab @.> Sent: 22 March 2021 15:17 To: BenvenLab/eSNPKaryotyping @.> Cc: Lauren Tierney @.>; Manual @.> Subject: Re: [BenvenLab/eSNPKaryotyping] unable to run Edit_dbSNP_Files function (#4)

Great! I'm glad to hear that! I hope this tool will be helpful to you! Uri

‫בתאריך יום ב׳, 22 במרץ 2021 ב-10:51 מאת ‪Lauren‬‏ <‪ @.***‬‏>:‬

Hi Uri,

thank you so much for getting back to me, but a colleague has just this minute figured out the issues, it was because my files had file extensions.

I will make a comment on your github page to explain my error.

thank you,

Lauren

---

From: BenvenLab @.> Sent: 22 March 2021 13:38 To: BenvenLab/eSNPKaryotyping @.> Cc: Lauren Tierney @.>; Manual @.> Subject: Re: [BenvenLab/eSNPKaryotyping] unable to run Edit_dbSNP_Files function (#4)

Hi Lauren, I'm Attaching a much more detailed manual for the eSNP-Kryotyping. Try following the steps there and let me know if the problem is still occurring.

Thanks, Uri

‫בתאריך יום ד׳, 17 במרץ 2021 ב-8:09 מאת ‪Lauren‬‏ @.*** ‬‏>:‬

Hi there,

I am unable to run this Edit_dbSNP_Files function.

I get this error : Error in file(file, "rt") : cannot open the connection

I downloaded a gft file per chromosome and they are stored in a directory called eSNP. The file path to an individual gtf files is like the following, with just the number on the end ranging from 1-24:

/Users/laurentierney/Desktop/eSNP/dbSNP_database_snp151_UCSC_01.gtf

this is the code I'm trying to run: edited <- Edit_dbSNP_Files("/Users/laurentierney/Desktop/eSNP", File_Name = "dbSNP_database_snp151_UCSC", Organism = "Human")

I have noticed when the error message comes back it is adding the number "1" to the end of this file path:

cannot open file '/Users/laurentierney/Desktop/eSNPdbSNP_database_snp151_UCSC1': No such file or directoryError in file(file, "rt") : cannot open the connection

Sorry if this is very basic but I can not understand what I am doing wrong

thanks

Lauren

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/BenvenLab/eSNPKaryotyping/issues/4, or unsubscribe < https://github.com/notifications/unsubscribe-auth/AD7XRUTI4DTP7WL65VLXARLTECLZPANCNFSM4ZKO2B2A>

.

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub< https://github.com/BenvenLab/eSNPKaryotyping/issues/4#issuecomment-804067882>, or unsubscribe< https://github.com/notifications/unsubscribe-auth/AO3JUULHSWCLE7ZX63IL7SDTE5B5NANCNFSM4ZKO2B2A>.

— You are receiving this because you commented. Reply to this email directly, view it on GitHub https://github.com/BenvenLab/eSNPKaryotyping/issues/4#issuecomment-804123419, or unsubscribe https://github.com/notifications/unsubscribe-auth/AD7XRUTB6STNOUBUEOLJLATTE5KPNANCNFSM4ZKO2B2A .

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHubhttps://github.com/BenvenLab/eSNPKaryotyping/issues/4#issuecomment-804143429, or unsubscribehttps://github.com/notifications/unsubscribe-auth/AO3JUULJDBANXW5AQHGUN7LTE5NO7ANCNFSM4ZKO2B2A.

Please consider the environment before printing this email. This email and its contents may be legally privileged and /or confidential. If it has come to you in error you should not copy or show it to anyone; nor should you take any action based on it, other than to reply to the sender to notify them of the error immediately and to then delete the email from your inbox and deleted items. Keele University staff and students are required to abide by the University's IT Policies when sending and receiving email. Keele University email is hosted by a cloud provider and may be stored outside of the UK.

[BenvenLab]

Hi Lauren,

For eSNP-Karyotyping you don't need to obtain the gene-expression levels because it uses the BAM files to measure the SNP ratios. For e-Karyotyping you do need to use cufflinks since the detection of chromosomal aberration is based on the changes in gene expression. In case you want to use e-Karyotyping, you have to use cufflinks or any other program to get the expression values of each gene in each sample (then I used cufflinks, now I use a program called stringtie...). Then you have to merge the table you have for each sample, to one big sample.You can do this step in Excel you can do it using VLOOKUP function or in R using the 'merge' function... I don't really know where you got stuck... where you were able to get the table of gene expression per sample and the problem is in the merging?

Uri

‫בתאריך יום ג׳, 4 במאי 2021 ב-5:55 מאת ‪Lauren‬‏ @.*** ‬‏>:‬

Hi Uri,

I am really keen to apply your eSNP karyotyping protocol to my research however, I am struggling with section 3.2.2

After running cufflinks, I'm struggling with how to produce the merged table you show in the manuscript, with Gene_id, Gene_name, Chromosome number, gene start position and then FPKM for each of the samples. How did you go from running cufflinks to producing this table?

Any help would be greatly appreciated,

thank you,

Lauren

---

From: BenvenLab @.> Sent: 22 March 2021 15:17 To: BenvenLab/eSNPKaryotyping @.> Cc: Lauren Tierney @.>; Manual @.> Subject: Re: [BenvenLab/eSNPKaryotyping] unable to run Edit_dbSNP_Files function (#4)

Great! I'm glad to hear that! I hope this tool will be helpful to you! Uri

‫בתאריך יום ב׳, 22 במרץ 2021 ב-10:51 מאת ‪Lauren‬‏ <‪ @.***‬‏>:‬

Hi Uri,

thank you so much for getting back to me, but a colleague has just this minute figured out the issues, it was because my files had file extensions.

I will make a comment on your github page to explain my error.

thank you,

Lauren

---

From: BenvenLab @.> Sent: 22 March 2021 13:38 To: BenvenLab/eSNPKaryotyping @.> Cc: Lauren Tierney @.>; Manual @.> Subject: Re: [BenvenLab/eSNPKaryotyping] unable to run Edit_dbSNP_Files function (#4)

Hi Lauren, I'm Attaching a much more detailed manual for the eSNP-Kryotyping. Try following the steps there and let me know if the problem is still occurring.

Thanks, Uri

‫בתאריך יום ד׳, 17 במרץ 2021 ב-8:09 מאת ‪Lauren‬‏ @.*** ‬‏>:‬

Hi there,

I am unable to run this Edit_dbSNP_Files function.

I get this error : Error in file(file, "rt") : cannot open the connection

I downloaded a gft file per chromosome and they are stored in a directory called eSNP. The file path to an individual gtf files is like the following, with just the number on the end ranging from 1-24:

/Users/laurentierney/Desktop/eSNP/dbSNP_database_snp151_UCSC_01.gtf

this is the code I'm trying to run: edited <- Edit_dbSNP_Files("/Users/laurentierney/Desktop/eSNP", File_Name = "dbSNP_database_snp151_UCSC", Organism = "Human")

I have noticed when the error message comes back it is adding the number "1" to the end of this file path:

cannot open file '/Users/laurentierney/Desktop/eSNPdbSNP_database_snp151_UCSC1': No such file or directoryError in file(file, "rt") : cannot open the connection

Sorry if this is very basic but I can not understand what I am doing wrong

thanks

Lauren

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/BenvenLab/eSNPKaryotyping/issues/4, or unsubscribe <

https://github.com/notifications/unsubscribe-auth/AD7XRUTI4DTP7WL65VLXARLTECLZPANCNFSM4ZKO2B2A>

.

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub<

https://github.com/BenvenLab/eSNPKaryotyping/issues/4#issuecomment-804067882>,

or unsubscribe<

https://github.com/notifications/unsubscribe-auth/AO3JUULHSWCLE7ZX63IL7SDTE5B5NANCNFSM4ZKO2B2A>.

— You are receiving this because you commented. Reply to this email directly, view it on GitHub < https://github.com/BenvenLab/eSNPKaryotyping/issues/4#issuecomment-804123419>,

or unsubscribe < https://github.com/notifications/unsubscribe-auth/AD7XRUTB6STNOUBUEOLJLATTE5KPNANCNFSM4ZKO2B2A>

.

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub< https://github.com/BenvenLab/eSNPKaryotyping/issues/4#issuecomment-804143429>, or unsubscribe< https://github.com/notifications/unsubscribe-auth/AO3JUULJDBANXW5AQHGUN7LTE5NO7ANCNFSM4ZKO2B2A>.

Please consider the environment before printing this email. This email and its contents may be legally privileged and /or confidential. If it has come to you in error you should not copy or show it to anyone; nor should you take any action based on it, other than to reply to the sender to notify them of the error immediately and to then delete the email from your inbox and deleted items. Keele University staff and students are required to abide by the University's IT Policies when sending and receiving email. Keele University email is hosted by a cloud provider and may be stored outside of the UK.

— You are receiving this because you commented. Reply to this email directly, view it on GitHub https://github.com/BenvenLab/eSNPKaryotyping/issues/4#issuecomment-831821328, or unsubscribe https://github.com/notifications/unsubscribe-auth/AD7XRUXIMO4OTUZ3B2AALKDTL7AAZANCNFSM4ZKO2B2A .

[laurentierney]

Hi Uri,

Apologies for the long email, but I am in need of your expertise!

I am using your eSNP karyotyping method, all appears to works fine when using your sample files from the paper but I am now running the same protocol on my own files and there have been various errors

firstly the accepted_hits.bam file was truncated, but this was resolved by changing the order of the functions, so reorderSam was ran first and then add or remove read groups.

However, I then ran into another problem where the sex chromosomes appear to be missed out of the VCF file (well there is only one row from the X chromsome and its is appended with characters like so ... chrX_K12709...

I'm not sure why the sex chromsomes are missing from the VCF file, do you have any thoughts?

I should also mention, then when I ran the protocol using your sample files from the paper (SRR3501465 & SRR3501466), I was not able to use your functions such as CreateVCF, instead I copied the underlying code into R and ran each command separately and that appeared to work, producing all the outputs, but I'm not sure if my outputs are correct. Is it possible to gain access to your intermediate files such as the accepted_hits.bam and VCF files, deletion.txt and variant table to try and compare?

thanks

Lauren

---

From: BenvenLab @.> Sent: 05 May 2021 17:35 To: BenvenLab/eSNPKaryotyping @.> Cc: Lauren Tierney @.>; Manual @.> Subject: Re: [BenvenLab/eSNPKaryotyping] unable to run Edit_dbSNP_Files function (#4)

Hi Lauren,

For eSNP-Karyotyping you don't need to obtain the gene-expression levels because it uses the BAM files to measure the SNP ratios. For e-Karyotyping you do need to use cufflinks since the detection of chromosomal aberration is based on the changes in gene expression. In case you want to use e-Karyotyping, you have to use cufflinks or any other program to get the expression values of each gene in each sample (then I used cufflinks, now I use a program called stringtie...). Then you have to merge the table you have for each sample, to one big sample.You can do this step in Excel you can do it using VLOOKUP function or in R using the 'merge' function... I don't really know where you got stuck... where you were able to get the table of gene expression per sample and the problem is in the merging?

Uri

‫בתאריך יום ג׳, 4 במאי 2021 ב-5:55 מאת ‪Lauren‬‏ @.*** ‬‏>:‬

Hi Uri,

I am really keen to apply your eSNP karyotyping protocol to my research however, I am struggling with section 3.2.2

After running cufflinks, I'm struggling with how to produce the merged table you show in the manuscript, with Gene_id, Gene_name, Chromosome number, gene start position and then FPKM for each of the samples. How did you go from running cufflinks to producing this table?

Any help would be greatly appreciated,

thank you,

Lauren

---

From: BenvenLab @.> Sent: 22 March 2021 15:17 To: BenvenLab/eSNPKaryotyping @.> Cc: Lauren Tierney @.>; Manual @.> Subject: Re: [BenvenLab/eSNPKaryotyping] unable to run Edit_dbSNP_Files function (#4)

Great! I'm glad to hear that! I hope this tool will be helpful to you! Uri

‫בתאריך יום ב׳, 22 במרץ 2021 ב-10:51 מאת ‪Lauren‬‏ <‪ @.***‬‏>:‬

Hi Uri,

thank you so much for getting back to me, but a colleague has just this minute figured out the issues, it was because my files had file extensions.

I will make a comment on your github page to explain my error.

thank you,

Lauren

---

From: BenvenLab @.> Sent: 22 March 2021 13:38 To: BenvenLab/eSNPKaryotyping @.> Cc: Lauren Tierney @.>; Manual @.> Subject: Re: [BenvenLab/eSNPKaryotyping] unable to run Edit_dbSNP_Files function (#4)

Hi Lauren, I'm Attaching a much more detailed manual for the eSNP-Kryotyping. Try following the steps there and let me know if the problem is still occurring.

Thanks, Uri

‫בתאריך יום ד׳, 17 במרץ 2021 ב-8:09 מאת ‪Lauren‬‏ @.*** ‬‏>:‬

Hi there,

I am unable to run this Edit_dbSNP_Files function.

I get this error : Error in file(file, "rt") : cannot open the connection

I downloaded a gft file per chromosome and they are stored in a directory called eSNP. The file path to an individual gtf files is like the following, with just the number on the end ranging from 1-24:

/Users/laurentierney/Desktop/eSNP/dbSNP_database_snp151_UCSC_01.gtf

this is the code I'm trying to run: edited <- Edit_dbSNP_Files("/Users/laurentierney/Desktop/eSNP", File_Name = "dbSNP_database_snp151_UCSC", Organism = "Human")

I have noticed when the error message comes back it is adding the number "1" to the end of this file path:

cannot open file '/Users/laurentierney/Desktop/eSNPdbSNP_database_snp151_UCSC1': No such file or directoryError in file(file, "rt") : cannot open the connection

Sorry if this is very basic but I can not understand what I am doing wrong

thanks

Lauren

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/BenvenLab/eSNPKaryotyping/issues/4, or unsubscribe <

https://github.com/notifications/unsubscribe-auth/AD7XRUTI4DTP7WL65VLXARLTECLZPANCNFSM4ZKO2B2A>

.

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub<

https://github.com/BenvenLab/eSNPKaryotyping/issues/4#issuecomment-804067882>,

or unsubscribe<

https://github.com/notifications/unsubscribe-auth/AO3JUULHSWCLE7ZX63IL7SDTE5B5NANCNFSM4ZKO2B2A>.

— You are receiving this because you commented. Reply to this email directly, view it on GitHub < https://github.com/BenvenLab/eSNPKaryotyping/issues/4#issuecomment-804123419>,

or unsubscribe < https://github.com/notifications/unsubscribe-auth/AD7XRUTB6STNOUBUEOLJLATTE5KPNANCNFSM4ZKO2B2A>

.

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub< https://github.com/BenvenLab/eSNPKaryotyping/issues/4#issuecomment-804143429>, or unsubscribe< https://github.com/notifications/unsubscribe-auth/AO3JUULJDBANXW5AQHGUN7LTE5NO7ANCNFSM4ZKO2B2A>.

Please consider the environment before printing this email. This email and its contents may be legally privileged and /or confidential. If it has come to you in error you should not copy or show it to anyone; nor should you take any action based on it, other than to reply to the sender to notify them of the error immediately and to then delete the email from your inbox and deleted items. Keele University staff and students are required to abide by the University's IT Policies when sending and receiving email. Keele University email is hosted by a cloud provider and may be stored outside of the UK.

— You are receiving this because you commented. Reply to this email directly, view it on GitHub https://github.com/BenvenLab/eSNPKaryotyping/issues/4#issuecomment-831821328, or unsubscribe https://github.com/notifications/unsubscribe-auth/AD7XRUXIMO4OTUZ3B2AALKDTL7AAZANCNFSM4ZKO2B2A .

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHubhttps://github.com/BenvenLab/eSNPKaryotyping/issues/4#issuecomment-832837605, or unsubscribehttps://github.com/notifications/unsubscribe-auth/AO3JUULRTFKFENAW7WXYR63TMFXVFANCNFSM4ZKO2B2A.

---

Please consider the environment before printing this email.This email and its contents may be legally privileged and /or confidential.If it has come to you in error you should not copy or show it to anyone; nor should you take any action based on it, other than to reply to the sender to notify them of the error immediately and to then delete the email from your inbox and deleted items. Keele University staff and students are required to abide by the University's IT Policies when sending and receiving email. Keele University email is hosted by a cloud provider and may be stored outside of the UK.

[laurentierney]

Hi Uri,

I'm very sorry to keep bothering you, but I'm still struggling to run your protocol successfully with my data files, more details on the query below:

I am using your eSNP karyotyping method, all appears to works fine when using your sample files from the paper but I am now running the same protocol on my own files and there have been various issues

Firstly, the accepted_hits.bam produced file was truncated, but this was resolved by changing the order of the commands (within the createVCF function), so reorderSam was ran first and then add or remove read groups, this allowed me to progress with subsequent steps but I noticed there were no allelic ratios produced for the X and Y chromsomes.

Then the code failed when trying to create deletion tables using the function from the R package, it happened when indexing the accepted hits bam file, an error appeared after reaching chromosome 23.

for one of the files, it run through to completion but again the plots did not show information for X or Y chromosomes.

Do you know why this might be?

additional info I am using the UCSC hg38 reference genome to map the files from your paper and the sequencing files from my project and I noticed the VCF from one of my files, there is only one row from the X chromsome and its is appended with characters like so ... chrX_K12709, which I think means its an ALT chromosome)

thanks

Lauren

---

From: Lauren Tierney @.> Sent: 01 September 2021 13:32 To: BenvenLab/eSNPKaryotyping @.> Subject: Re: [BenvenLab/eSNPKaryotyping] unable to run Edit_dbSNP_Files function (#4)

Hi Uri,

Apologies for the long email, but I am in need of your expertise!

I am using your eSNP karyotyping method, all appears to works fine when using your sample files from the paper but I am now running the same protocol on my own files and there have been various errors

firstly the accepted_hits.bam file was truncated, but this was resolved by changing the order of the functions, so reorderSam was ran first and then add or remove read groups.

However, I then ran into another problem where the sex chromosomes appear to be missed out of the VCF file (well there is only one row from the X chromsome and its is appended with characters like so ... chrX_K12709...

I'm not sure why the sex chromsomes are missing from the VCF file, do you have any thoughts?

I should also mention, then when I ran the protocol using your sample files from the paper (SRR3501465 & SRR3501466), I was not able to use your functions such as CreateVCF, instead I copied the underlying code into R and ran each command separately and that appeared to work, producing all the outputs, but I'm not sure if my outputs are correct. Is it possible to gain access to your intermediate files such as the accepted_hits.bam and VCF files, deletion.txt and variant table to try and compare?

thanks

Lauren

---

From: BenvenLab @.> Sent: 05 May 2021 17:35 To: BenvenLab/eSNPKaryotyping @.> Cc: Lauren Tierney @.>; Manual @.> Subject: Re: [BenvenLab/eSNPKaryotyping] unable to run Edit_dbSNP_Files function (#4)

Hi Lauren,

For eSNP-Karyotyping you don't need to obtain the gene-expression levels because it uses the BAM files to measure the SNP ratios. For e-Karyotyping you do need to use cufflinks since the detection of chromosomal aberration is based on the changes in gene expression. In case you want to use e-Karyotyping, you have to use cufflinks or any other program to get the expression values of each gene in each sample (then I used cufflinks, now I use a program called stringtie...). Then you have to merge the table you have for each sample, to one big sample.You can do this step in Excel you can do it using VLOOKUP function or in R using the 'merge' function... I don't really know where you got stuck... where you were able to get the table of gene expression per sample and the problem is in the merging?

Uri

‫בתאריך יום ג׳, 4 במאי 2021 ב-5:55 מאת ‪Lauren‬‏ @.*** ‬‏>:‬

Hi Uri,

I am really keen to apply your eSNP karyotyping protocol to my research however, I am struggling with section 3.2.2

After running cufflinks, I'm struggling with how to produce the merged table you show in the manuscript, with Gene_id, Gene_name, Chromosome number, gene start position and then FPKM for each of the samples. How did you go from running cufflinks to producing this table?

Any help would be greatly appreciated,

thank you,

Lauren

---

From: BenvenLab @.> Sent: 22 March 2021 15:17 To: BenvenLab/eSNPKaryotyping @.> Cc: Lauren Tierney @.>; Manual @.> Subject: Re: [BenvenLab/eSNPKaryotyping] unable to run Edit_dbSNP_Files function (#4)

Great! I'm glad to hear that! I hope this tool will be helpful to you! Uri

‫בתאריך יום ב׳, 22 במרץ 2021 ב-10:51 מאת ‪Lauren‬‏ <‪ @.***‬‏>:‬

Hi Uri,

thank you so much for getting back to me, but a colleague has just this minute figured out the issues, it was because my files had file extensions.

I will make a comment on your github page to explain my error.

thank you,

Lauren

---

From: BenvenLab @.> Sent: 22 March 2021 13:38 To: BenvenLab/eSNPKaryotyping @.> Cc: Lauren Tierney @.>; Manual @.> Subject: Re: [BenvenLab/eSNPKaryotyping] unable to run Edit_dbSNP_Files function (#4)

Hi Lauren, I'm Attaching a much more detailed manual for the eSNP-Kryotyping. Try following the steps there and let me know if the problem is still occurring.

Thanks, Uri

‫בתאריך יום ד׳, 17 במרץ 2021 ב-8:09 מאת ‪Lauren‬‏ @.*** ‬‏>:‬

Hi there,

I am unable to run this Edit_dbSNP_Files function.

I get this error : Error in file(file, "rt") : cannot open the connection

I downloaded a gft file per chromosome and they are stored in a directory called eSNP. The file path to an individual gtf files is like the following, with just the number on the end ranging from 1-24:

/Users/laurentierney/Desktop/eSNP/dbSNP_database_snp151_UCSC_01.gtf

this is the code I'm trying to run: edited <- Edit_dbSNP_Files("/Users/laurentierney/Desktop/eSNP", File_Name = "dbSNP_database_snp151_UCSC", Organism = "Human")

I have noticed when the error message comes back it is adding the number "1" to the end of this file path:

cannot open file '/Users/laurentierney/Desktop/eSNPdbSNP_database_snp151_UCSC1': No such file or directoryError in file(file, "rt") : cannot open the connection

Sorry if this is very basic but I can not understand what I am doing wrong

thanks

Lauren

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/BenvenLab/eSNPKaryotyping/issues/4, or unsubscribe <

https://github.com/notifications/unsubscribe-auth/AD7XRUTI4DTP7WL65VLXARLTECLZPANCNFSM4ZKO2B2A>

.

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub<

https://github.com/BenvenLab/eSNPKaryotyping/issues/4#issuecomment-804067882>,

or unsubscribe<

https://github.com/notifications/unsubscribe-auth/AO3JUULHSWCLE7ZX63IL7SDTE5B5NANCNFSM4ZKO2B2A>.

— You are receiving this because you commented. Reply to this email directly, view it on GitHub < https://github.com/BenvenLab/eSNPKaryotyping/issues/4#issuecomment-804123419>,

or unsubscribe < https://github.com/notifications/unsubscribe-auth/AD7XRUTB6STNOUBUEOLJLATTE5KPNANCNFSM4ZKO2B2A>

.

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub< https://github.com/BenvenLab/eSNPKaryotyping/issues/4#issuecomment-804143429>, or unsubscribe< https://github.com/notifications/unsubscribe-auth/AO3JUULJDBANXW5AQHGUN7LTE5NO7ANCNFSM4ZKO2B2A>.

Please consider the environment before printing this email. This email and its contents may be legally privileged and /or confidential. If it has come to you in error you should not copy or show it to anyone; nor should you take any action based on it, other than to reply to the sender to notify them of the error immediately and to then delete the email from your inbox and deleted items. Keele University staff and students are required to abide by the University's IT Policies when sending and receiving email. Keele University email is hosted by a cloud provider and may be stored outside of the UK.

— You are receiving this because you commented. Reply to this email directly, view it on GitHub https://github.com/BenvenLab/eSNPKaryotyping/issues/4#issuecomment-831821328, or unsubscribe https://github.com/notifications/unsubscribe-auth/AD7XRUXIMO4OTUZ3B2AALKDTL7AAZANCNFSM4ZKO2B2A .

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHubhttps://github.com/BenvenLab/eSNPKaryotyping/issues/4#issuecomment-832837605, or unsubscribehttps://github.com/notifications/unsubscribe-auth/AO3JUULRTFKFENAW7WXYR63TMFXVFANCNFSM4ZKO2B2A.

---

Please consider the environment before printing this email.This email and its contents may be legally privileged and /or confidential.If it has come to you in error you should not copy or show it to anyone; nor should you take any action based on it, other than to reply to the sender to notify them of the error immediately and to then delete the email from your inbox and deleted items. Keele University staff and students are required to abide by the University's IT Policies when sending and receiving email. Keele University email is hosted by a cloud provider and may be stored outside of the UK.