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.

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

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

Divergence Evolution of Duplicated Yeast Genes in Protein Interaction Network

In situ dissecting the evolution of gene duplication with different histone modification patterns based on high-throughput data analysis in <i>Arabidopsis thaliana</i>

Comparative Gene Retention Analysis in Barley, Wild Emmer, and Bread Wheat Pangenome Lines Reveals Factors Affecting Gene Retention Following Gene Duplication

Histone Modification Pattern Evolution after Yeast Gene Duplication.

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

Cis-regulatory Change and Expression Divergence Between Duplicate Genes Formed by Genome Duplication of Arabidopsis Thaliana

Coordinated Evolution Of Transcriptional And Post-Transcriptional Regulation For Mitochondrial Functions In Yeast Strains

Genomic Background Predicts the Fate of Duplicated Genes: Evidence from the Yeast Genome

Single-cell eQTL mapping in yeast reveals a tradeoff between growth and reproduction

Statistical Analysis Reveals Co-Expression Patterns of Many Pairs of Genes in Yeast Are Jointly Regulated by Interacting Loci

The 3-Dimensional Genome Drives the Evolution of Asymmetric Gene Duplicates via Enhancer Capture-Divergence

Asymmetric evolution of human transcription factor regulatory networks.

Duplicate genes increase gene expression diversity within and between species

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

Ohnologs and SSD Paralogs Differ in Genomic and Expression Features Related to Dosage Constraints

The Role of the 3-Dimensional Genome in New Gene Evolution

Distant regulatory effects of genetic variation in multiple human tissues

Aberrant Gene Expression in Humans

Epistasis And Destabilizing Mutations Shape Gene Expression Variability In Humans Via Distinct Modes Of Action