@swbarnes2 cool. never realized there was a difference.
@oyvindbusk thanks, I also tried this and ended up with similar thresholds (male < 40% and female > 50%). I didn't try to filter out pseudoautosomal regions since their coordinates differ across species and assembly versions (based on PAR coordinates at:
).
Looking at 322 CCLE samples, 233 were called Male, 73 Female, and 10 Unknown (which is >= 40% and <= 50%). Out of the 233 Male, only 5 would have been called differently with your thresholds. I will see if I can check the thresholds against a different approach. Also, a lot of the CCLE cells have copy number amplifications / deletions, so these results might be skewed by that.
Here is the distribution of nHet / nHomo for chrX in CCLE samples (I used this instead of nHet/(nHet+nHomo)). The 2 vertical blue lines are equivalent to 40% and 50% thresholds, and the 30% threshold is the red line.
Seqanswers Leaderboard Ad
Collapse
Announcement
Collapse
No announcement yet.
X
-
Gender = biological sex + culture. You don't care about people's gender, you care about their sex. (And even the biology is not black and white 100% of the time)
Leave a comment:
-
How about using % heterozygosity on X (without the pseudoautosomal regions (X:60000-2699520 and X:154931043-155260560). In our lab, male = < 30 % and female = > 50 %.
Leave a comment:
-
Testing sex determination with CCLE samples
I've tried using [num reads mapped to chrX] / [num reads mapped to chrY]
to determine sex in some CCLE exome-seq samples. The ratios turned out to be:
9.4 -- s1
304.6 -- s2
272.9 -- s3
168.3 -- s4
220.6 -- s5
297.8 -- s6
226.1 -- s7
257.1 -- s8
241.9 -- s9
287.0 -- s10
278.6 -- s11
260.3 -- s12
9.7 -- s13
8.7 -- s14
261.2 -- s15
279.3 -- s16
9.0 -- s17
8.5 -- s18
260.7 -- s19
297.4 -- s20
8.7 -- s21
261.8 -- s22
189.0 -- s23
147.4 -- s24
291.2 -- s25
[num reads mapped to chrX] / [num reads mapped to chrY] is < 10 for all male samples and > 100 for all female samples.
Still, I'm not sure whether these thresholds are stable across exome-seq kits, gene panels, etc. I wonder if there's a more robust way to determine sex.Last edited by bw.; 02-13-2014, 11:07 AM.
Leave a comment:
-
in plink there is an option to check sex using X chr, I hope you are looking for this..
plink --bfile data --impute-sex --make-bed --out newfile
Leave a comment:
-
I would second that -- map reads against a panel of Y-chromosome genes/exons.
Leave a comment:
-
gender check using sequencing data
Hi All,
I have exome sequencing data from siblings and would like to confirm their gender and relationship using genetics to make sure they are in fact siblings. I also want to confirm that the gender information is correct. Are there tools out there that can do this using exome seq data?
Thanks,Tags: None
Latest Articles
Collapse
-
by seqadmin
Metagenomics has improved the way researchers study microorganisms across diverse environments. Historically, studying microorganisms relied on culturing them in the lab, a method that limits the investigation of many species since most are unculturable1. Metagenomics overcomes these issues by allowing the study of microorganisms regardless of their ability to be cultured or the environments they inhabit. Over time, the field has evolved, especially with the advent...-
Channel: Articles
09-23-2024, 06:35 AM -
-
by seqadmin
During the COVID-19 pandemic, scientists observed that while some individuals experienced severe illness when infected with SARS-CoV-2, others were barely affected. These disparities left researchers and clinicians wondering what causes the wide variations in response to viral infections and what role genetics plays.
Jean-Laurent Casanova, M.D., Ph.D., Professor at Rockefeller University, is a leading expert in this crossover between genetics and infectious...-
Channel: Articles
09-09-2024, 10:59 AM -
ad_right_rmr
Collapse
News
Collapse
Topics | Statistics | Last Post | ||
---|---|---|---|---|
Started by seqadmin, 10-02-2024, 04:51 AM
|
0 responses
13 views
0 likes
|
Last Post
by seqadmin
10-02-2024, 04:51 AM
|
||
Started by seqadmin, 10-01-2024, 07:10 AM
|
0 responses
21 views
0 likes
|
Last Post
by seqadmin
10-01-2024, 07:10 AM
|
||
Started by seqadmin, 09-30-2024, 08:33 AM
|
0 responses
25 views
0 likes
|
Last Post
by seqadmin
09-30-2024, 08:33 AM
|
||
Started by seqadmin, 09-26-2024, 12:57 PM
|
0 responses
18 views
0 likes
|
Last Post
by seqadmin
09-26-2024, 12:57 PM
|
Leave a comment: