Abstract:In comparisons between mutant and wild-type genotypes, transcriptome analysis can reveal the direct impacts of a mutation, together with the homeostatic responses of the biological system. Recent studies have highlighted that, when the effects of homozygosity for recessive mutations are studied in non-isogenic backgrounds, genes located proximal to the mutation on the same chromosome often appear over-represented among those genes identified as differentially expressed (DE). One hypothesis suggests that DE genes chromosomally linked to a mutation may not reflect functional responses to the mutation but, instead, result from an unequal distribution of expression quantitative trait loci (eQTLs) between sample groups of mutant or wild-type genotypes. This is problematic because eQTL expression differences are difficult to distinguish from genes that are DE due to functional responses to a mutation. Here we show that chromosomally co-located differentially expressed genes (CC-DEGs) are also observed in analyses of dominant mutations in heterozygotes. We define a method and a metric to quantify, in RNA-sequencing data, localised differential allelic representation (DAR) between those sample groups subjected to differential expression analysis. We show how the DAR metric can predict regions prone to eQTL-driven differential expression, and how it can improve functional enrichment analyses through gene exclusion or weighting-based approaches. Advantageously, this improved ability to identify probable eQTLs also reveals examples of CC-DEGs that are likely to be functionally related to a mutant phenotype. This supports a long-standing prediction that selection for advantageous linkage disequilibrium influences chromosome evolution. By comparing the genomes of zebrafish ( Danio rerio ) and medaka ( Oryzias latipes ), a teleost with a conserved ancestral karyotype, we find possible examples of chromosomal aggregation of CC-DEGs during evolution of the zebrafish lineage. Our method for DAR analysis requires only RNA-sequencing data, facilitating its application across new and existing datasets. Many human-relevant diseases result from genetic mutations that disrupt cellular functions. We can model these mutations in other organisms (e.g. mouse, zebrafish) and employ gene expression analysis (transcriptomics) to determine how mutations directly affect cells and how cells adjust expression of their genes to compensate for these mutations. In our transcriptome analyses of dominant disease-causative mutations in zebrafish, we identified an interesting phenomenon where a disproportionate number of differentially expressed genes reside on the same chromosome as a mutated gene. Here, we provide strong evidence supporting that the differential expression of some of these chromosomally co-located genes is not due to the mutation but is due to differential segregation of gene alleles with innately different expression levels (i.e. expression quantitative trait loci, eQTLs). We have developed a procedure to measure the likelihood of differential gene expression being due to an eQTL. This allows us to compensate for the presence of such eQTLs in bioinformatic analyses. Our procedure, Differential Allelic Representation (DAR) analysis, revealed evidence for aggregation of genes with related functions on the same chromosome over evolutionary timescales. DAR analysis allows disentanglement of eQTLs from mutation-dependent gene expression responses, thereby permitting more comprehensive investigation of transcriptome data.

Differential allelic representation (DAR) identifies candidate eQTLs and improves transcriptome analysis

Robust identification of regulatory variants (eQTLs) using a differential expression framework developed for RNA-sequencing

Uncovering genetic regulatory network divergence between duplicate genes using yeast eQTL landscape.

Haplotype-aware modeling of cis-regulatory effects highlights the gaps remaining in eQTL data

Probing the limits of cis-acting gene regulation using a model of allelic imbalance quantitative trait loci

Single-cell Eqtlgen Consortium: a Personalized Understanding of Disease

A New Approach for Discovering Functional Links Connecting Non-Coding Regulatory Variants to Gene Targets

Aberrant Gene Expression in Humans

SNP-Based and Kmer-Based eQTL Analysis Using Transcriptome Data

Structural Basis of Differential Gene Expression at eQTLs Loci from High-Resolution Ensemble Models of 3D Single-Cell Chromatin Conformations

Shared and distinct molecular effects of regulatory genetic variants provide insight into mechanisms of distal enhancer-promoter communication

The Role of Cis-Regulatory Motifs and Genetical Control of Expression in the Divergence of Yeast Duplicate Genes

Harnessing mobile genetic elements to explore gene regulation

Efficient and Effective Control of Confounding in Eqtl Mapping Studies Through Joint Differential Expression and Mendelian Randomization Analyses.

Alternative Strategies for Development of a Reference Transcriptome for Quantification of Allele Specific Expression in Organisms Having Sparse Genomic Resources.

Transcriptional Regulatory Pattern in Yeast Revealed Through Expression Quantitative Trait Locus Mapping

Gene Set Enrichment in Eqtl Data Identifies Novel Annotations and Pathway Regulators.

Allelic Expression of Deleterious Protein-Coding Variants Across Human Tissues

Distributed Eqtl Analysis with Auxiliary Information

Survey of Human Chromosome 21 Gene Expression Effects on Early Development in Danio rerio

A Powerful Approach Reveals Numerous Expression Quantitative Trait Haplotypes In Multiple Tissues