Inter-individual gene expression variability implies stable regulation of brain-biased genes across organs in three ray-finned fishes

Phenotypic variation among individuals provides the raw material for evolution, and gene expression is a key mediator between genetic and phenotypic variation. As for phenotypes, the range of gene expression is limited and biased by evolutionary and developmental constraints. Observed expression variability due to biomolecular stochasticity and cell-to-cell heterogeneity has been well-studied in isogenic populations of unicellular organisms. However, for multicellular organisms with a diversity of cells and tissues sharing the same genetic background, the interplay between expression variability, gene and organ function, and gene regulation remains an open question. Here, we use highly multiplexed 3-end bulk RNA sequencing to generate transcriptome profiles spanning at least nine organs in outbred individuals of three ray-finned fish species: zebrafish, Northern pike, and spotted gar. Per organ, we quantify individual-to-individual gene expression variability independent of mean expression level. Lowly variable genes are enriched in cellular housekeeping functions whereas highly variable genes are enriched in stimulus-response functions. Furthermore, highly variable genes evolve under weaker purifying selection at the protein-coding sequence, indicating that intra-species expression variability predicts inter-species protein sequence divergence. Genes that are broadly expressed across organs are both highly expressed and lowly variable, whereas organ-biased genes are typically highly variable within their top organ. Among organ-biased genes, patterns of expression variance are dependent on the top organ. Specifically, brain- and gonads-biased genes have lowly variable expression across different organs, suggesting stabilizing selection. These patterns suggest that gene regulatory mechanisms evolve under organ-specific selective pressures.