Tutorials › High-Dimensional Genomic Data Analysis › Y-Chromosome Calls for Sex Checks

Variant-Level and Sample-Level Quality Control · Tutorial 33 of 795

Y-Chromosome Calls for Sex Checks

Learn to summarize Y-chromosome calls as supporting evidence while separating informative genotypes from missing, low-quality, or otherwise unreliable data.

Intermediate 11 min read

What You'll Learn

  • Define an eligible Y-chromosome marker set for a sample-level check.
  • Calculate a call rate using only markers that pass the stated site-level criteria.
  • Distinguish reference and alternate genotypes from missing and low-quality calls.
  • Interpret patterns across multiple non-pseudoautosomal Y markers.
  • Recognize suspicious heterozygous calls and other signs that a pattern needs investigation.
  • Explain why Y-call patterns are supporting evidence, not a stand-alone determination.

Why Y-Chromosome Calls Can Support a Sex Check

Genetic sex checks compare observed chromosome-specific genotype patterns with an expectation tied to the presence or absence of particular sex chromosomes. The earlier material on X-chromosome dosage introduced the complementary X-based signal. Here the focus is the Y chromosome: a pattern of reliable calls at appropriate Y-chromosome loci can support evidence that Y-chromosome DNA is present in a sample.

The useful signal is not simply the number of alternate alleles. At a well-chosen Y-chromosome variant, a high-quality reference call and a high-quality alternate call are both evidence that the locus was observed. The allele may vary among individuals; what matters for this basic check is whether the sample has credible genotype calls at multiple suitable Y loci. By contrast, a missing genotype gives no allele observation at that locus.

Definition: For a Y-chromosome sex check, an informative Y call is a non-missing genotype call that passes the specified quality criteria at an eligible, reliably assayed Y-chromosome locus. A missing call has no genotype observation. A low-quality call is recorded but fails the stated criteria, so it is not counted as informative for the check.

Choose Loci and Define Callability Before Summarizing

A useful summary starts with a carefully defined set of loci. For a simple presence check, use Y-specific loci outside the pseudoautosomal regions, where X and Y chromosomes share homologous sequence. The selected loci should also be suitable for the genotyping platform and have reliable chromosome annotation and mapping. Poorly mapped or otherwise unreliable loci can produce apparent calls that do not represent the intended Y-chromosome signal.

Set quality rules before comparing samples. For example, an analysis might require a genotype quality score of at least 20 and a read depth of at least 10. Those values are illustrative, not universal recommendations: appropriate thresholds depend on the assay and its validated quality measures. A site that fails sample-level quality criteria is not informative for that sample, even if a genotype value appears in the file. Similarly, a site that was excluded globally should not be counted as an eligible opportunity for any sample.

Definition: Suppose \(K\) eligible Y loci remain after applying the stated site-level rules. Let \(C\) be the number of those loci at which the sample has a non-missing genotype that passes the sample-level quality rules. The sample’s Y call rate is \(C/K\). Report the numerator and denominator, as well as the rate, so that a low rate is not mistaken for strong evidence of chromosome absence.

This ratio describes how much of the chosen panel was observed well enough to interpret; it is not a probability that the sample carries a Y chromosome. The denominator must be appropriate to the question. If a locus is not included in the assay, or was removed for failing a site-level rule, it does not belong in \(K\). If it is eligible but missing or low quality in one sample, it remains in \(K\) but does not contribute to \(C\).

Worked Example: Counting Calls at Eligible Y Loci

Consider a hypothetical sample and a panel of eight eligible, non-pseudoautosomal Y loci. The sample has six passing genotype calls: four are reference calls and two are alternate calls. One additional locus has a genotype that fails the pre-specified quality threshold, and one is missing.

Only the six passing calls count as informative. The two alternate calls are not the only evidence: the four high-quality reference calls also show that those loci were successfully observed. The low-quality genotype and the missing genotype do not count.

The call rate is

$$ \frac{C}{K}=\frac{6}{8}=0.750 $$

Thus, 75.0% of the eligible loci were observed at sufficient quality in this hypothetical sample. This supports a description of the sample as having a moderately complete set of Y-locus calls, but the call rate alone does not establish whether a Y chromosome is present. The pattern, quality, and suitability of the loci must also be considered.

Interpreting the Pattern Across Loci

For a typical single-copy region outside the pseudoautosomal regions, an individual with a Y chromosome has one allele at each locus. Genotype files may represent that allele as a haploid call or, depending on software and conventions, as a homozygous diploid call. Check the data representation before interpreting apparent genotype classes. A reference allele at one locus and an alternate allele at another can both be compatible with a Y-chromosome pattern; they are not, by themselves, evidence of a problem.

A pattern of passing calls across multiple suitable Y loci is more useful than a single isolated call. Multiple calls reduce dependence on any one marker, which might be affected by local assay problems or a mapping error. But the calls are not necessarily independent: nearby loci may share technical issues, and sample-wide coverage or assay problems can affect many loci at once. Therefore, a large count of calls should be read alongside the locus set and quality measures, not treated as a formal probability calculation.

Worked Example: Reference Calls Are Still Evidence

In a hypothetical panel of 20 eligible Y loci, sample A has 18 passing calls: 13 reference and 5 alternate. Sample B has 2 passing calls, both alternate, with 18 loci missing or failing quality rules.

For sample A, the call rate is \(18/20=0.900\). The broad set of passing calls, including reference calls, provides substantially more usable Y-locus information than a summary that counted alternate calls alone would suggest. For sample B, the call rate is \(2/20=0.100\); two alternate calls are observations, but most of the panel is not interpretable.

It would be wrong to say that sample B provides stronger Y evidence simply because both of its observed calls are alternate. The samples differ greatly in how much of the eligible panel was observed. A low call rate leaves the Y-based check inconclusive; it does not show that the chromosome is absent.

This interpretation rests on basic assumptions: the loci are correctly assigned to the Y chromosome, the chosen regions are appropriate for the check, and passing calls are sufficiently reliable to represent the sample’s genotype. If those assumptions are doubtful, the summary may reflect annotation, mapping, or assay behavior rather than chromosome content.

Missingness and Quality Failures Are Not Alleles

A no-call is not a reference genotype and not an alternate genotype. It records that the data do not provide a usable genotype at that locus. A low-quality genotype is also not a reliable allele observation for this purpose. Both may arise from inadequate coverage, assay failure, or other technical causes. The reason may vary, but neither category should be silently recoded as evidence for or against a Y chromosome.

The distinction between a globally excluded site and a sample-specific failure is important. A globally excluded site was not considered suitable for this check and should be removed from the denominator. A sample-specific missing or failed call at an otherwise eligible site stays in the denominator and lowers that sample’s call rate. Keeping these cases separate makes the summary reproducible and helps identify whether low callability is specific to a sample or shared across the panel.

Worked Example: Separating Excluded Sites From Failed Calls

A hypothetical panel begins with 10 Y loci. Two fail a site-level rule for all samples, so the eligible denominator is \(K=8\). In one sample, five loci have passing calls, one has a recorded genotype below the sample-level quality threshold, and two have no genotype.

The low-quality call and the two missing calls are all non-informative for this sample, but the two globally excluded loci are not part of its denominator. Therefore,

$$ \text{Y call rate}=\frac{5}{8}=0.625 $$

The sample has passing calls at 62.5% of the eligible loci. Dividing by 10 would mix in sites that were not eligible for any sample, while counting the low-quality genotype would overstate usable information. The arithmetic is meaningful only after both kinds of site selection have been stated.

Diagnostics and Common Pitfalls

A useful review considers both the call-rate summary and the pattern that produced it. Check whether the eligible loci are correctly annotated, whether calls pass the same stated quality rules, and whether missingness is widespread or concentrated in particular samples or loci. Previously introduced missingness diagnostics can help describe such patterns. A high call rate is not a guarantee that every call is correct; a low call rate is a warning that the sample-level Y pattern may be uninformative.

At a non-pseudoautosomal, single-copy Y locus, an apparent heterozygous genotype is unexpected under the simple single-allele representation. It should be treated as a diagnostic flag, not as a routine Y call to interpret in the same way as a clear haploid or software-encoded homozygous call. Review whether the locus is truly non-pseudoautosomal, how the file encodes ploidy, and whether the genotype passes quality checks. Possible technical explanations include mapping ambiguity or sample mixture; the observation alone does not identify its cause.

  • Mistake: treating every missing call as a reference call. A missing genotype has no observed allele. Keep it missing and exclude it from the informative-call count.
  • Mistake: treating low call rate as evidence that no Y chromosome is present. Low call rate means that the panel provides limited usable information. It may reflect technical failure, so it cannot be interpreted as chromosome absence on its own.
  • Mistake: counting only alternate calls as Y evidence. A passing reference call is also a genotype observation at a Y locus. Summarize passing calls, not just alternate alleles.
  • Mistake: interpreting an apparent heterozygous call without checking representation and locus type. First verify ploidy encoding, location, and quality; otherwise a technical or annotation issue can be mistaken for a biological pattern.

Y-chromosome calls are supporting evidence within quality control, not a stand-alone determination about an individual. A sound summary states which loci were eligible, how call quality was defined, how many calls passed, and how much of the panel was missing or failed. It avoids converting technical absence of data into a biological conclusion.

Check Your Understanding

Answer each question, then compare your reasoning with the solutions below.

  1. A sample has 12 eligible loci, with 9 passing calls, 2 missing calls, and 1 low-quality call. What is its call rate, and which observations count?
  2. Why is a passing reference genotype at a suitable Y locus relevant to a Y-based check?
  3. A panel has 15 eligible loci. A sample has one passing alternate call and 14 missing calls. Is this strong evidence that the sample lacks a Y chromosome? Explain.
  4. At a stated non-pseudoautosomal Y locus, a sample has an apparent heterozygous genotype. Give two checks to make before interpreting it.

Solutions. (1) The call rate is \(9/12=0.750\). Only the nine passing calls count; the two missing and one low-quality genotype do not. (2) A passing reference genotype is still an observed genotype at that locus, so it contributes to evidence that the locus was successfully assayed. The allele need not be alternate for the call to be informative. (3) No. The call rate is \(1/15\approx0.0667\), so most eligible loci provide no usable observation. The pattern is largely uninformative and does not establish chromosome absence. (4) Verify that the locus is correctly annotated outside the pseudoautosomal regions and check how the genotype file encodes ploidy; also review the call’s quality and mapping reliability before treating it as meaningful.

Clinical Biostats