Cis-regulatory Changes Research Articles

Analysis of cis-regulatory changes underlying phenotype divergence shaped by domestication in pigs

BackgroundCis-regulatory elements (CREs) are regions of DNA that regulate the expression of nearby genes. Changes in these elements can lead to phenotypic variations and adaptations in different populations. However, the regulatory dynamics underlying the local adaptation of traits remain poorly understood in Chinese and Western pigs. By comparing the chromatin accessibility profiles of skeletal muscle, liver, and fat between these two pig populations, we aimed to identify key regulatory elements that could explain phenotypic differences observed between the two groups.ResultsOur results revealed that the genome-wide chromatin accessibility profiles were largely similar at a qualitative level within tissues. However, we also identified local regions that exhibited quantitative differences, most of which occurred in liver tissue. Interestingly, we found that most of the increased chromatin accessibility in the livers of Chinese pigs was associated with tissue-specific openness. Furthermore, we observed a positive correlation between the ATAC-seq signal at the transcript start site (TSS) and the expression levels of nearby genes. Motif enrichment analysis revealed NR2F1 as a key regulator in Chinese pigs. Differentially expressed genes (DEGs) in Chinese pigs showed enrichment for NR2F1 response targets. One of the genes regulated by NR2F1 in Chinese pigs, NPC1, harbored a high alternative allelic frequency in the intron region.ConclusionOverall, our study provides valuable insights into the regulatory dynamics underlying phenotypic variation in pigs. By elucidating the role of CREs in driving phenotypic variation, we can better understand the genetic basis of complex traits and potentially identify targets for genetic improvement in livestock breeding programs.

Evolution of the regulation of developmental gene expression in blind Mexican cavefish.

Developmental evolution and diversification of morphology can arise through changes in the regulation of gene expression or protein-coding sequence. To unravel mechanisms underlying early developmental evolution in cavefish of the species Astyanax mexicanus, we compared transcriptomes of surface-dwelling and blind cave-adapted morphs at the end of gastrulation. Twenty percent of the transcriptome was differentially expressed. Allelic expression ratios in cave X surface hybrids showed that cis-regulatory changes are the quasi-exclusive contributors to inter-morph variations in gene expression. Among a list of 108 genes with change at the cis-regulatory level, we explored the control of expression of rx3, which is a master eye gene. We discovered that cellular rx3 levels are cis-regulated in a cell-autonomous manner, whereas rx3 domain size depends on non-autonomous Wnt and Bmp signalling. These results highlight how uncoupled mechanisms and regulatory modules control developmental gene expression and shape morphological changes. Finally, a transcriptome-wide search for fixed coding mutations and differential exon use suggested that variations in coding sequence have a minor contribution. Thus, during early embryogenesis, changes in gene expression regulation are the main drivers of cavefish developmental evolution.

Genomic signatures associated with recurrent scale loss in cyprinid fish.

Scale morphology represents a fundamental feature of fish and a key evolutionary trait underlying fish diversification. Despite frequent and recurrent scale loss throughout fish diversification, comprehensive genome-wide analyses of the genomic signatures associated with scale loss in divergent fish lineages remain scarce. In the current study, we investigated genome-wide signatures, specifically convergent protein-coding gene loss, amino acid substitutions, and cis-regulatory sequence changes, associated with recurrent scale loss in two divergent Cypriniformes lineages based on large-scale genomic, transcriptomic, and epigenetic data. Results demonstrated convergent changes in many genes related to scale formation in divergent scaleless fish lineages, including loss of P/Q-rich scpp genes (e.g. scpp6 and scpp7), accelerated evolution of non-coding elements adjacent to the fgf and fgfr genes, and convergent amino acid changes in genes (e.g. snap29) under relaxed selection. Collectively, these findings highlight the existence of a shared genetic architecture underlying recurrent scale loss in divergent fish lineages, suggesting that evolutionary outcomes may be genetically repeatable and predictable in the convergence of scale loss in fish.

Abstract P15: Aberrant Non-canonical NF-κB Signalling Reprograms the Epigenome Landscape to Drive Oncogenic Transcriptomes in Multiple Myeloma

Abstract Activation of the NF-κB pathway is recognized to be necessary for the survival and progression of haematological cancers including Multiple Myeloma (MM). In particular, recurrent genetic lesions resulting in the chronic induction of non-canonical NF-κB (ncNF-κB) signalling have been exclusively linked to MM. Such mutations among patients are associated with significant proliferative and metastatic advantages and have been causally linked to treatment resistance. How the constitutive activity of ncNF-κB transcription factors can support such adaptations in this common haematological malignancy has been a fundamental question in cancer. In this study, we captured the epigenomic landscape of patient-derived MM cell lines harbouring recurrent NF-κB activating mutations through perturbation of nfkb2/p52 and integrated these cis-regulatory changes with MM patient epigenome data. Here, we find that sustained NF-κB/p52 levels resulting from such mutations favours the recruitment of enhancers beyond the normal B-cell repertoire. Furthermore, through targeted disruption of p52, we characterise how such enhancers are complicit in the formation of super-enhancers and the establishment of cis-regulatory interactions with myeloma dependencies during constitutive activation of p52. Finally, we functionally validate the pathological impact of these cis-regulatory modules on cell and tumour phenotypes using in vitro and in vivo models, confirming RGS1 as a novel p52-dependent myeloma driver. We conclude that the divergent epigenomic reprogramming enforced by aberrant non-canonical NF-κB signalling potentiates transcriptional programs beneficial for multiple myeloma progression. Citation Format: Daniel A. Ang, Jean-Michel Carter, Kamalakshi Deka, Joel H.L. Tan, Jianbiao Zhou, Qingfeng Chen, Wee Joo Chng, Nathan Harmston, Yinghui Li. Aberrant Non-canonical NF-κB Signalling Reprograms the Epigenome Landscape to Drive Oncogenic Transcriptomes in Multiple Myeloma [abstract]. In: Proceedings of Frontiers in Cancer Science; 2023 Nov 6-8; Singapore. Philadelphia (PA): AACR; Cancer Res 2024;84(8_Suppl):Abstract nr P15.

Cell-type-specific cis-regulatory divergence in gene expression and chromatin accessibility revealed by human-chimpanzee hybrid cells.

Although gene expression divergence has long been postulated to be the primary driver of human evolution, identifying the genes and genetic variants underlying uniquely human traits has proven to be quite challenging. Theory suggests that cell-type-specific cis-regulatory variants may fuel evolutionary adaptation due to the specificity of their effects. These variants can precisely tune the expression of a single gene in a single cell-type, avoiding the potentially deleterious consequences of trans-acting changes and non-cell type-specific changes that can impact many genes and cell types, respectively. It has recently become possible to quantify human-specific cis-acting regulatory divergence by measuring allele-specific expression in human-chimpanzee hybrid cells-the product of fusing induced pluripotent stem (iPS) cells of each species in vitro. However, these cis-regulatory changes have only been explored in a limited number of cell types. Here, we quantify human-chimpanzee cis-regulatory divergence in gene expression and chromatin accessibility across six cell types, enabling the identification of highly cell-type-specific cis-regulatory changes. We find that cell-type-specific genes and regulatory elements evolve faster than those shared across cell types, suggesting an important role for genes with cell-type-specific expression in human evolution. Furthermore, we identify several instances of lineage-specific natural selection that may have played key roles in specific cell types, such as coordinated changes in the cis-regulation of dozens of genes involved in neuronal firing in motor neurons. Finally, using novel metrics and a machine learning model, we identify genetic variants that likely alter chromatin accessibility and transcription factor binding, leading to neuron-specific changes in the expression of the neurodevelopmentally important genes FABP7 and GAD1. Overall, our results demonstrate that integrative analysis of cis-regulatory divergence in chromatin accessibility and gene expression across cell types is a promising approach to identify the specific genes and genetic variants that make us human.

Cell-type-specific cis-regulatory divergence in gene expression and chromatin accessibility revealed by human-chimpanzee hybrid cells

Although gene expression divergence has long been postulated to be the primary driver of human evolution, identifying the genes and genetic variants underlying uniquely human traits has proven to be quite challenging. Theory suggests that cell-type-specific cis-regulatory variants may fuel evolutionary adaptation due to the specificity of their effects. These variants can precisely tune the expression of a single gene in a single cell-type, avoiding the potentially deleterious consequences of trans-acting changes and non-cell type-specific changes that can impact many genes and cell types, respectively. It has recently become possible to quantify human-specific cis-acting regulatory divergence by measuring allele-specific expression in human-chimpanzee hybrid cells—the product of fusing induced pluripotent stem (iPS) cells of each species in vitro. However, these cis-regulatory changes have only been explored in a limited number of cell types. Here, we quantify human-chimpanzee cis-regulatory divergence in gene expression and chromatin accessibility across six cell types, enabling the identification of highly cell-type-specific cis-regulatory changes. We find that cell-type-specific genes and regulatory elements evolve faster than those shared across cell types, suggesting an important role for genes with cell-type-specific expression in human evolution. Furthermore, we identify several instances of lineage-specific natural selection that may have played key roles in specific cell types, such as coordinated changes in the cis-regulation of dozens of genes involved in neuronal firing in motor neurons. Finally, using novel metrics and a machine learning model, we identify genetic variants that likely alter chromatin accessibility and transcription factor binding, leading to neuron-specific changes in the expression of the neurodevelopmentally important genes FABP7 and GAD1. Overall, our results demonstrate that integrative analysis of cis-regulatory divergence in chromatin accessibility and gene expression across cell types is a promising approach to identify the specific genes and genetic variants that make us human.

Tissue-specific and cis-regulatory changes underlie parallel, adaptive gene expression evolution in house mice.

Changes in gene regulation have long been appreciated as a driving force of adaptive evolution, however the relative contributions of cis- and trans-acting changes to gene regulation over short evolutionary timescales remain unclear. Instances of recent, parallel phenotypic evolution provide an opportunity to assess whether parallel patterns are seen at the level of gene expression, and to assess the relative contribution of cis- and trans- changes to gene regulation in the early stages of divergence. Here, we studied gene expression in liver and brown adipose tissue in two wild-derived strains of house mice that independently adapted to cold, northern environments, and we compared them to a strain of house mice from a warm, tropical environment. To investigate gene regulatory evolution, we studied expression in parents and allele-specific expression in F1 hybrids of crosses between warm-adapted and cold-adapted strains. First, we found that the different cold-adapted mice showed both unique and shared changes in expression, but that the proportion of shared changes (i.e. parallelism) was greater than expected by chance. Second, we discovered that expression evolution occurred largely at tissue-specific and cis-regulated genes, and that these genes were over-represented in parallel cases of evolution. Finally, we integrated the expression data with scans for selection in natural populations and found substantial parallelism in the two northern populations for genes under selection. Furthermore, selection outliers were associated with cis-regulated genes more than expected by chance; cis-regulated genes under selection influenced phenotypes such as body size, immune functioning, and activity level. These results demonstrate that parallel patterns of gene expression in mice that have independently adapted to cold environments are driven largely by tissue-specific and cis-regulatory changes, providing insight into the mechanisms of adaptive gene regulatory evolution at the earliest stages of divergence.

Gene loss and cis-regulatory novelty shaped core histone gene evolution in the apiculate yeast Hanseniaspora uvarum.

Core histone genes display a remarkable diversity of cis-regulatory mechanisms despite their protein sequence conservation. However, the dynamics and significance of this regulatory turnover are not well understood. Here, we describe the evolutionary history of core histone gene regulation across 400 million years in budding yeasts. We find that canonical mode of core histone regulation-mediated by the trans-regulator Spt10-is ancient, likely emerging between 320 and 380 million years ago and is fixed in the majority of extant species. Unexpectedly, we uncovered the emergence of a novel core histone regulatory mode in the Hanseniaspora genus, from its fast-evolving lineage, which coincided with the loss of 1 copy of its paralogous core histone genes. We show that the ancestral Spt10 histone regulatory mode was replaced, via cis-regulatory changes in the histone control regions, by a derived Mcm1 histone regulatory mode and that this rewiring event occurred with no changes to the trans-regulator, Mcm1, itself. Finally, we studied the growth dynamics of the cell cycle and histone synthesis in genetically modified Hanseniaspora uvarum. We find that H. uvarum divides rapidly, with most cells completing a cell cycle within 60 minutes. Interestingly, we observed that the regulatory coupling between histone and DNA synthesis was lost in H. uvarum. Our results demonstrate that core histone gene regulation was fixed anciently in budding yeasts, however it has greatly diverged in the Hanseniaspora fast-evolving lineage.

Gene regulation and speciation in a migratory divide between songbirds

Behavioral variation abounds in nature. This variation is important for adaptation and speciation, but its molecular basis remains elusive. Here, we use a hybrid zone between two subspecies of songbirds that differ in migration – an ecologically important and taxonomically widespread behavior---to gain insight into this topic. We measure gene expression in five brain regions. Differential expression between migratory states was dominated by circadian genes in all brain regions. The remaining patterns were largely brain-region specific. For example, expression differences between the subspecies that interact with migratory state likely help maintain reproductive isolation in this system and were documented in only three brain regions. Contrary to existing work on regulatory mechanisms underlying species-specific traits, two lines of evidence suggest that trans- (vs. cis) regulatory changes underlie these patterns – no evidence for allele-specific expression in hybrids and minimal associations between genomic differentiation and expression differences. Additional work with hybrids shows expression levels were often distinct (transgressive) from parental forms. Behavioral contrasts and functional enrichment analyses allowed us to connect these patterns to mitonuclear incompatibilities and compensatory responses to stress that could exacerbate selection on hybrids and contribute to speciation.

Ascertainment Bias in the Genomic Test of Positive Selection on Regulatory Sequences.

Evolution of gene expression mediated by cis-regulatory changes is thought to be an important contributor to organismal adaptation, but identifying adaptive cis-regulatory changes is challenging due to the difficulty in knowing the expectation under no positive selection. A new approach for detecting positive selection on transcription factor binding sites (TFBSs) was recently developed, thanks to the application of machine learning in predicting transcription factor (TF) binding affinities of DNA sequences. Given a TFBS sequence from a focal species and the corresponding inferred ancestral sequence that differs from the former at n sites, one can predict the TF-binding affinities of many n-step mutational neighbors of the ancestral sequence and obtain a null distribution of the derived binding affinity, which allows testing whether the binding affinity of the real derived sequence deviates significantly from the null distribution. Applying this test genomically to all experimentally identified binding sites of 3 TFs in humans, a recent study reported positive selection for elevated binding affinities of TFBSs. Here, we show that this genomic test suffers from an ascertainment bias because, even in the absence of positive selection for strengthened binding, the binding affinities of known human TFBSs are more likely to have increased than decreased in evolution. We demonstrate by computer simulation that this bias inflates the false positive rate of the selection test. We propose several methods to mitigate the ascertainment bias and show that almost all previously reported positive selection signals disappear when these methods are applied.

EPCO-08. CONVERGENCE OF NON-CODING SOMATIC MUTATIONS ON CANCER DRIVER PLEXUSES IN PEDIATRIC HIGH-GRADE GLIOMA

Abstract Central nervous system tumors are the most common pediatric malignancies after leukemia, and the most common cause of cancer death among persons ages 0 – 14 [1,2]. Despite its low incidence rate, pediatric high-grade gliomas (pHGGs) account for over 40% of all childhood brain tumor deaths [3, 4]. As such, it remains the focus of many research teams across the globe. The discovery of recurrent mutations in genes encoding histone H3.3 [5] was a breakthrough in the genetic basis of pHGGs and highlighted the relevance of epigenetic reconfiguration for oncogene activation. However, much work is still needed to uncover driver alterations responsible for tumor initiation or progression. pHGGs have frequent epigenetic dysregulation [6,7] and abundant somatic mutations in the noncoding regions of the genome [8], but the functional relevance of the noncoding somatic mutations has not been evaluated systematically. Here, we assess correlations between noncoding somatic mutations in pediatric brain-specific cis-regulatory elements (CREs) and changes in gene expression. We demonstrate the sparsity of recurrent mutations in CREs. Then illustrate the power of graph database in consolidating multi-omics datasets for efficient queries. Our network approach for charting the genomic, epigenetic, and transcriptomic landscape of pHGG allow us to 1) identify clusters of CREs associated with dysregulated gene expression, 2) nominate transcription factors whose binding affinities are recurrently impacted by somatic mutations found in pHGGs, and 3) establish regulatory networks potentially guiding pHGG tumor growth and progression.

Asymmetric Evolution of Protein Domains in the Leucine-Rich Repeat Receptor-Like Kinase Family of Plant Signaling Proteins.

The coding sequences of developmental genes are expected to be deeply conserved, with cis-regulatory change driving the modulation of gene function. In contrast, proteins with roles in defense are expected to evolve rapidly, in molecular arms races with pathogens. However, some gene families include both developmental and defense genes. In these families, does the tempo and mode of evolution differ between genes with divergent functions, despite shared ancestry and structure? The leucine-rich repeat receptor-like kinase (LRR-RLKs) protein family includes members with roles in plant development and defense, thus providing an ideal system for answering this question. LRR-RLKs are receptors that traverse plasma membranes. LRR domains bind extracellular ligands; RLK domains initiate intracellular signaling cascades in response to ligand binding. In LRR-RLKs with roles in defense, LRR domains evolve faster than RLK domains. To determine whether this asymmetry extends to LRR-RLKs that function primarily in development, we assessed evolutionary rates and tested for selection acting on 11 subfamilies of LRR-RLKs, using deeply sampled protein trees. To assess functional evolution, we performed heterologous complementation assays in Arabidopsis thaliana (Arabidopsis). We found that the LRR domains of all tested LRR-RLK proteins evolved faster than their cognate RLK domains. All tested subfamilies of LRR-RLKs had strikingly similar patterns of molecular evolution, despite divergent functions. Heterologous transformation experiments revealed that multiple mechanisms likely contribute to the evolution of LRR-RLK function, including escape from adaptive conflict. Our results indicate specific and distinct evolutionary pressures acting on LRR versus RLK domains, despite diverse organismal roles for LRR-RLK proteins.

Environmentally robust cis-regulatory changes underlie rapid climatic adaptation

Changes in gene expression are thought to play a major role in adaptive evolution. While it is known that gene expression is highly sensitive to the environment, very few studies have determined the influence of genetic and environmental effects on adaptive gene expression differences in natural populations. Here, we utilize allele-specific expression to characterize cis and trans gene regulatory divergence in temperate and tropical house mice in two metabolic tissues under two thermal conditions. First, we show that gene expression divergence is pervasive between populations and across thermal conditions, with roughly 5 to 10% of genes exhibiting genotype-by-environment interactions. Second, we found that most expression divergence was due to cis-regulatory changes that were stable across temperatures. In contrast, patterns of expression plasticity were largely attributable to trans-effects, which showed greater sensitivity to temperature. Nonetheless, we found a small subset of temperature-dependent cis-regulatory changes, thereby identifying loci underlying expression plasticity. Finally, we performed scans for selection in wild house mice to identify genomic signatures of rapid adaptation. Genomic outliers were enriched in genes with evidence for cis-regulatory divergence. Notably, these genes were associated with phenotypes that affected body weight and metabolism, suggesting that cis-regulatory changes are a possible mechanism for adaptive body size evolution between populations. Our results show that gene expression plasticity, largely controlled in trans, may facilitate the colonization of new environments, but that evolved changes in gene expression are largely controlled in cis, illustrating the genetic and nongenetic mechanisms underlying the establishment of populations in new environments.

Analysis of the genetic loci of pigment pattern evolution in vertebrates.

Vertebrate pigmentation patterns are amongst the best characterised model systems for studying the genetic basis of adaptive evolution. The wealth of available data on the genetic basis for pigmentation evolution allows for analysis of trends and quantitative testing of evolutionary hypotheses. We employed Gephebase, a database of genetic variants associated with natural and domesticated trait variation, to examine trends in how cis-regulatory and coding mutations contribute to vertebrate pigmentation phenotypes, as well as factors that favour one mutation type over the other. We found that studies with lower ascertainment bias identified higher proportions of cis-regulatory mutations, and that cis-regulatory mutations were more common amongst animals harbouring a higher number of pigment cell classes. We classified pigmentation traits firstly according to their physiological basis and secondly according to whether they affect colour or pattern, and identified that carotenoid-based pigmentation and variation in pattern boundaries are preferentially associated with cis-regulatory change. We also classified genes according to their developmental, cellular, and molecular functions. We found a greater proportion of cis-regulatory mutations in genes implicated in upstream developmental processes compared to those involved in downstream cellular functions, and that ligands were associated with a higher proportion of cis-regulatory mutations than their respective receptors. Based on these trends, we discuss future directions for research in vertebrate pigmentation evolution.

Evolution of Spatial and Temporal cis-Regulatory Divergence in Sticklebacks.

Cis-regulatory changes are thought to play a major role in adaptation. Threespine sticklebacks have repeatedly colonized freshwater habitats in the Northern Hemisphere, where they have evolved a suite of phenotypes that distinguish them from marine populations, including changes in physiology, behavior, and morphology. To understand the role of gene regulatory evolution in adaptive divergence, here we investigate cis-regulatory changes in gene expression between marine and freshwater ecotypes through allele-specific expression (ASE) in F1 hybrids. Surveying seven ecologically relevant tissues, including three sampled across two developmental stages, we identified cis-regulatory divergence affecting a third of genes, nearly half of which were tissue-specific. Next, we compared allele-specific expression in dental tissues at two timepoints to characterize cis-regulatory changes during development between marine and freshwater fish. Applying a genome-wide test for selection on cis-regulatory changes, we find evidence for lineage-specific selection on several processes between ecotypes, including the Wnt signaling pathway in dental tissues. Finally, we show that genes with ASE, particularly those that are tissue-specific, are strongly enriched in genomic regions of repeated marine-freshwater divergence, supporting an important role for these cis-regulatory differences in parallel adaptive evolution of sticklebacks to freshwater habitats. Altogether, our results provide insight into the cis-regulatory landscape of divergence between stickleback ecotypes across tissues and during development, and support a fundamental role for tissue-specific cis-regulatory changes in rapid adaptation to new environments.

Natural overexpression of CAROTENOID CLEAVAGE DIOXYGENASE 4 in tomato alters carotenoid flux.

Carotenoids and apocarotenoids function as pigments and flavor volatiles in plants that enhance consumer appeal and offer health benefits. Tomato (Solanum lycopersicum.) fruit, especially those of wild species, exhibit a high degree of natural variation in carotenoid and apocarotenoid contents. Using positional cloning and an introgression line (IL) of Solanum habrochaites "LA1777', IL8A, we identified carotenoid cleavage dioxygenase 4 (CCD4) as the factor responsible for controlling the dark orange fruit color. CCD4b expression in ripe fruit of IL8A plants was ∼8,000 times greater than that in the wild type, presumably due to 5' cis-regulatory changes. The ShCCD4b-GFP fusion protein localized in the plastid. Phytoene, ζ-carotene, and neurosporene levels increased in ShCCD4b-overexpressing ripe fruit, whereas trans-lycopene, β-carotene, and lutein levels were reduced, suggestive of feedback regulation in the carotenoid pathway by an unknown apocarotenoid. Solid-phase microextraction-gas chromatography-mass spectrometry analysis showed increased levels of geranylacetone and β-ionone in ShCCD4b-overexpressing ripe fruit coupled with a β-cyclocitral deficiency. In carotenoid-accumulating Escherichia coli strains, ShCCD4b cleaved both ζ-carotene and β-carotene at the C9-C10 (C9'-C10') positions to produce geranylacetone and β-ionone, respectively. Exogenous β-cyclocitral decreased carotenoid synthesis in the ripening fruit of tomato and pepper (Capsicum annuum), suggesting feedback inhibition in the pathway. Our findings will be helpful for enhancing the aesthetic and nutritional value of tomato and for understanding the complex regulatory mechanisms of carotenoid and apocarotenoid biogenesis.

Accounting for cis-regulatory constraint prioritizes genes likely to affect species-specific traits

Measuring allele-specific expression in interspecies hybrids is a powerful way to detect cis-regulatory changes underlying adaptation. However, it remains difficult to identify genes most likely to explain species-specific traits. Here, we outline a simple strategy that leverages population-scale allele-specific RNA-seq data to identify genes that show constrained cis-regulation within species yet show divergence between species. Applying this strategy to data from human-chimpanzee hybrid cortical organoids, we identify signatures of lineage-specific selection on genes related to saccharide metabolism, neurodegeneration, and primary cilia. We also highlight cis-regulatory divergence in CUX1 and EDNRB that may shape the trajectory of human brain development.

Trans-regulatory changes underpin the evolution of the Drosophila immune response.

When an animal is infected, the expression of a large suite of genes is changed, resulting in an immune response that can defend the host. Despite much evidence that the sequence of proteins in the immune system can evolve rapidly, the evolution of gene expression is comparatively poorly understood. We therefore investigated the transcriptional response to parasitoid wasp infection in Drosophila simulans and D. sechellia. Although these species are closely related, there has been a large scale divergence in the expression of immune-responsive genes in their two main immune tissues, the fat body and hemocytes. Many genes, including those encoding molecules that directly kill pathogens, have cis regulatory changes, frequently resulting in large differences in their expression in the two species. However, these changes in cis regulation overwhelmingly affected gene expression in immune-challenged and uninfected animals alike. Divergence in the response to infection was controlled in trans. We argue that altering trans-regulatory factors, such as signalling pathways or immune modulators, may allow natural selection to alter the expression of large numbers of immune-responsive genes in a coordinated fashion.

Evolution of stickleback spines through independent cis-regulatory changes at HOXDB

Understanding the mechanisms leading to new traits or additional features in organisms is a fundamental goal of evolutionary biology. We show that HOXDB regulatory changes have been used repeatedly in different fish genera to alter the length and number of the prominent dorsal spines used to classify stickleback species. In Gasterosteus aculeatus (typically ‘three-spine sticklebacks’), a variant HOXDB allele is genetically linked to shortening an existing spine and adding an additional spine. In Apeltes quadracus (typically ‘four-spine sticklebacks’), a variant HOXDB allele is associated with lengthening a spine and adding an additional spine in natural populations. The variant alleles alter the same non-coding enhancer region in the HOXDB locus but do so by diverse mechanisms, including single-nucleotide polymorphisms, deletions and transposable element insertions. The independent regulatory changes are linked to anterior expansion or contraction of HOXDB expression. We propose that associated changes in spine lengths and numbers are partial identity transformations in a repeating skeletal series that forms major defensive structures in fish. Our findings support the long-standing hypothesis that natural Hox gene variation underlies key patterning changes in wild populations and illustrate how different mutational mechanisms affecting the same region may produce opposite gene expression changes with similar phenotypic outcomes.

Convergent copy number increase of genes associated with freshwater colonization in fishes.

Copy number variation (CNV) can cause phenotypic changes. However, in contrast to amino acid substitutions and cis-regulatory changes, little is known about the functional categories of genes in which CNV is important for adaptation to novel environments. It is also unclear whether the same genes repeatedly change the copy numbers for adapting to similar environments. Here, we investigate CNV associated with freshwater colonization in fishes, which was observed multiple times across different lineages. Using 48 ray-finned fishes across diverse orders, we identified 23 genes whose copy number increases were associated with freshwater colonization. These genes showed enrichment for peptide receptor activity, hexosyltransferase activity and unsaturated fatty acid metabolism. We further revealed that three of the genes showed copy number increases in freshwater populations compared to marine ancestral populations of the stickleback genus Gasterosteus. These results indicate that copy number increases of genes involved in fatty acid metabolism (FADS2), immune function (PSMB8a) and thyroid hormone metabolism (UGT2) may be important for freshwater colonization at both the inter-order macroevolutionary scale and at the intra-genus microevolutionary scale. Further analysis across diverse taxa will help to understand the role of CNV in the adaptation to novel environments. This article is part of the theme issue 'Genetic basis of adaptation and speciation: from loci to causative mutations'.