Why Mapping Windows Need an Explicit Definition
A molecular QTL analysis tests whether genetic variants are associated with a molecular phenotype, such as expression of a gene. Before testing, the analyst must decide which variant–phenotype pairs are in scope. A search near a gene is commonly called cis; a search farther away, potentially across the genome, is commonly called trans. These labels describe a mapping rule, not a proven biological mechanism.
There is no single universal distance at which a variant changes from cis to trans. A window is an operational choice: it specifies the genomic locations to search for a given molecular phenotype. That choice affects which associations can be discovered, what “cis” means in the reported results, how much computation is required, and how many statistical tests contribute to the search.
How to Specify a Cis Window
A common rule centers the window on a gene’s transcription start site (TSS). Let \(t_g\) be the TSS position for gene \(g\), and let \(W\) be a chosen distance in base pairs. A symmetric TSS-centered window is the interval from \(t_g-W\) to \(t_g+W\). Under an inclusive endpoint convention, a variant at either boundary is included.
Another rule starts from the gene body, whose annotated start and end positions are \(a_g\) and \(b_g\). A gene-body window can extend \(L\) base pairs upstream and \(R\) base pairs downstream, producing an interval from \(a_g-L\) to \(b_g+R\). The upstream and downstream directions refer to transcription, so the genomic coordinates used for those flanks depend on whether the gene is on the forward or reverse strand.
These definitions are not interchangeable. A TSS-centered window is anchored to one point; a gene-body window includes the span of the annotated gene and can be asymmetric. A protocol should state which is used, the chosen distance or flank lengths, and how endpoints are handled. The TSS also depends on the transcript annotation: genes can have multiple transcripts with different start sites. An analysis should identify the annotation and transcript-selection rule rather than silently treating every gene as having one unambiguous TSS.
Coordinates must be tied to a reference genome build. A position such as 10,000,000 is meaningful only in relation to a particular genome assembly and coordinate convention. For example, one may specify that variant positions and gene coordinates use the same build and 1-based reference coordinates. If variants and annotations use different builds, positions must be harmonized before applying the window rule. Otherwise, a variant can be assigned to the wrong interval even when the stated distance rule is clear.
Worked Example: Applying a TSS-Centered Cis Window
Suppose a hypothetical gene has a TSS at position 10,000,000 on a specified reference build. The analysis defines cis as within 500,000 base pairs of the TSS, inclusive. Classify variants at positions 9,500,000, 10,420,000, and 10,500,001.
The interval is calculated as \(10{,}000{,}000 \pm 500{,}000\), giving \([9{,}500{,}000,\ 10{,}500{,}000]\). The variant at 9,500,000 is at the lower boundary and is included. The variant at 10,420,000 is 420,000 base pairs from the TSS and is included. The variant at 10,500,001 is 500,001 base pairs away, so it falls outside this cis window.
This classification follows the stated positional rule. It does not establish that either included variant directly regulates the gene, nor that the excluded variant has no biological effect on it.
Defining the Trans Search
Once the cis rule is defined, a straightforward pair-specific convention calls variants outside that gene’s cis interval trans for that gene. This includes distant variants on the same chromosome as well as variants on other chromosomes. Some studies instead restrict “trans” to different-chromosome pairs or impose a separate minimum distance. Those conventions can be useful, but they produce different search sets and must be reported explicitly.
Classification is relative to the target gene. A variant can be inside the cis window for one gene and outside the cis window for another. Overlapping windows make this especially clear: the same variant may be tested as cis for one phenotype and trans for a different phenotype. Thus “cis variant” is often shorthand; the more precise statement is “cis for gene \(g\) under the specified window rule.”
Worked Example: Classifying Same-Chromosome and Different-Chromosome Pairs
Consider a hypothetical gene on chromosome 4 with a TSS at position 20,000,000. Its cis window is defined as \(\pm 250{,}000\) base pairs. Three variants are considered: one at position 20,180,000 on chromosome 4, one at position 21,000,000 on chromosome 4, and one at position 8,000,000 on chromosome 12.
The cis interval is \([19{,}750{,}000,\ 20{,}250{,}000]\). The chromosome 4 variant at 20,180,000 is 180,000 base pairs from the TSS and is cis for this gene. The variant at 21,000,000 is 1,000,000 base pairs from the TSS, so it is outside the cis interval and is trans under the pair-specific complement rule. The chromosome 12 variant is also outside the interval and is trans under that rule.
If the study defines trans to mean “on a different chromosome,” only the chromosome 12 pair qualifies as trans. The distant chromosome 4 pair would instead be outside the cis search but not included in that particular trans category. Reporting the rule prevents these distinct classifications from being confused.
Search-Space Size, Computation, and Power
For each molecular phenotype, the number of variant–gene tests depends on the chosen window and on which variants are available. A narrow cis interval usually contains fewer variants than a broad interval, reducing the number of tests and the computational work required. A genome-wide trans search can involve far more pairs: in a simple design, each variant may be tested against many molecular phenotypes. The number of pairs can become large even when the number of donors is unchanged.
Search-space size also affects statistical power when an analysis controls the chance of any false positive across the tests. To see the basic relationship, suppose \(M\) null hypotheses are tested and each test has a valid p-value. If each test uses a threshold \(\alpha/M\), the probability of at least one false positive is at most \(\alpha\). This follows because the probability of the union of false-positive events is no greater than the sum of their probabilities:
This bound does not require tests to be independent, though it can be conservative when tests are correlated. It shows why, under this particular family-wise error rule, a larger search requires a smaller per-test threshold. A smaller threshold makes it harder for a true association to be declared significant, all else equal. The precise power still depends on sample size, effect size, phenotype variability, variant frequency, linkage disequilibrium, and the chosen analysis procedure. A narrow window does not automatically have higher power if it excludes the causal variant or a useful proxy.
Worked Example: Comparing the Burden of Two Searches
Suppose a hypothetical analysis tests 20,000 variants for one gene in a cis search and 1,000,000 variant–gene pairs in a trans search. For illustration, use the Bonferroni threshold \(\alpha/M\) with a family-wise level \(\alpha=0.05\). Assume the p-values are valid under their null hypotheses.
For the cis search, the threshold is \(0.05/20{,}000=0.0000025=2.5\times10^{-6}\). For the trans search, it is \(0.05/1{,}000{,}000=0.00000005=5.0\times10^{-8}\). A hypothetical p-value of \(1.0\times10^{-7}\) is below the cis threshold, so it meets this illustrative cis criterion. It is above the trans threshold, so it does not meet the illustrative trans criterion.
This comparison isolates the effect of the number of tests under the chosen correction; it is not a claim that the two searches have identical test dependence or biological relevance. Correlation among variants can affect the effective information in a search, while the Bonferroni bound remains valid when the individual p-values are calibrated. The broader search has a greater testing burden under this rule, but it may be needed to detect long-range or cross-chromosome associations.
Biological Interpretation Depends on the Boundary
A broader cis window can include regulatory elements that are far from a TSS, including elements between genes or in neighboring genes. This can increase the chance of including a variant that is biologically relevant to the target phenotype. It can also make a “cis association” less local in the everyday sense and complicate interpretation when nearby genes share regulatory regions or are in linkage disequilibrium. An associated variant in a broad window does not necessarily act through the nearest gene.
A narrower window offers a more restrictive definition of local association, often with fewer tests and a simpler positional description. But a fixed boundary can exclude a functional variant beyond it. Such a variant may then be labelled trans by the operational rule even if it acts through a nearby regulatory element. Conversely, being within a cis window does not prove a local mechanism. Genomic distance is a useful search convention, not a causal test.
Window choices also affect comparisons across studies. Two studies can use the same word, “cis,” while testing different intervals, reference builds, or transcript annotations. Their reported counts of cis QTLs are not directly comparable unless those choices are aligned or their differences are made explicit. A transparent report gives the build, annotation, TSS or gene-body rule, window size or flank lengths, endpoint convention, and trans definition.
Common Mistakes
- Treating one distance as a biological law. A variant just outside a chosen boundary is not known to be biologically unrelated to the gene. The boundary is an operational convention.
- Calling every trans pair different-chromosome. Some definitions include distant variants on the same chromosome; others do not. State whether trans is the complement of the cis window or uses a separate rule.
- Ignoring strand and transcript annotation. A TSS is transcript-specific, and upstream/downstream flanks depend on strand. Specify the annotation and how the TSS is selected.
- Comparing windows across incompatible coordinates. Variant and gene positions must use the same genome build and coordinate convention before distances are calculated.
- Equating a cis association with direct regulation. Position alone does not establish mechanism, identify the causal variant, or prove that the target gene mediates an effect.
Check Your Understanding
Answer each question using a clearly stated mapping rule.
- A gene has a TSS at position 5,000,000 and uses an inclusive \(\pm 100{,}000\)-base-pair window. What is the interval, and how is a variant at 4,900,000 classified?
- Why can the same variant be cis for one gene and trans for another?
- Give two distinct conventions for defining a trans search, and explain why they may yield different results.
- How can widening a cis window affect biological interpretation, computation, and power under a fixed family-wise error procedure?
- What coordinate and annotation details should be reported so that another analyst can reproduce a window assignment?