Our work
Publications
Lab publications
Click a title to see the authors, abstract and a link to the article. Siddiq Lab members highlighted in bold.
Preprints
Haque T, Siddiq MA, Duveau F, Wittkopp PJ
bioRxiv (preprint)
Abstract
Genetically identical cells grown in the same environment show variation in gene expression known as expression noise. Expression noise can be heritable and impact fitness, making it subject to natural selection. Increasing expression noise for the Saccharomyces cerevisiae TDH3 gene was shown to be beneficial in glucose-based media when mean TDH3 expression was far from the fitness optimum but deleterious when it was close to this optimum. Here, we show that growth on different carbon sources alters the effects of new mutations on TDH3 expression noise and examine the fitness effects of changing expression noise. In galactose-based media, we observed the same relationship between expression noise and fitness seen in glucose-based media, but in glycerol- and ethanol-based media, we observed the opposite relationship or no significant relationship, respectively. Using simulations of single-cell organisms, we found that these differences were most likely explained by environment-specific relationships between gene expression and fitness. We also found that, far from the optimum, the fitness effects of noise were greatest when expression was highly heritable between mother and daughter cells. The empirical observations and simulations reported in this study show how environments influence both the production of expression noise and its impacts on fitness.
Ertl HA, Bayala EX, Siddiq MA, Wittkopp PJ
bioRxiv (preprint)
Abstract
Pioneer factors are critical for gene regulation and development because they bind chromatin and make DNA more accessible for binding by other transcription factors. The pioneer factor Grainy head (Grh) is present across metazoans and has been shown to retain a role in epithelium development in fruit flies, nematodes, and mice despite extensive divergence in both amino acid sequence and length. Here, we investigate the evolution of Grh function by comparing the effects of the fly ( Drosophila melanogaster) and worm ( Caenorhabditis elegans) Grh orthologs on chromatin accessibility, gene expression, embryonic development, and viability in transgenic D. melanogaster. We found that the Caenorhabditis elegans ortholog rescued cuticle development but not full embryonic viability in Drosophila melanogaster grh null mutants. At the molecular level, the C. elegans ortholog only partially rescued chromatin accessibility and gene expression. Divergence in the disordered N-terminus of the Grh protein contributes to these differences in embryonic viability and molecular phenotypes. These data show how pioneer factors can diverge in sequence and function at the molecular level while retaining conserved developmental functions at the organismal level. Despite divergence in a disordered region that affects function at both molecular and organismal levels, the Caenorhabditis elegans Grainy head (Grh) protein rescued cuticle morphology in D. melanogaster embryos.
Siddiq Lab (2026–present)
Siddiq MA, Kania HP, Brown NJ, Wittkopp PJ
Molecular Biology and Evolution 43(8): msag203 (2026)
Abstract
Changes in regulatory sequences controlling the timing and activity of gene products underlie much of natural phenotypic variation. Yet, identifying which divergent sites matter and how they impact gene expression remains challenging. Here, we investigate how transcriptional activity and homeostatic responsiveness of orthologous promoters of the metabolic gene TDH3 evolved among Saccharomyces yeast. We found that promoter expression level increased specifically in the Saccharomyces cerevisiae lineage and that a substantial part of this increase was caused by genetic variants located between the well-characterized, conserved binding sites for two direct transcriptional regulators. These nucleotide changes altered the promoters’ expression levels while leaving the expression dynamics conserved. Further, the effects of these nucleotide changes were only seen in the presence of a third transcription factor, TYE7p, which is known to be recruited by the other transcription factors through protein–protein interactions. These results suggest that the cis-regulatory changes act through their influence on the collective assembly/activation of a transcription factor complex and that changes acting through such a mechanism can allow distinct parts of gene expression, such as expression level and dynamics, to separately diverge.
Rapala JR, Siddiq M, Wittkopp PJ, O’Meara MJ, O’Meara TR
Nature Communications 17: 4593 (2026)
Abstract
A central tenet of biology is that protein structure mediates the sequence-function relationship. Recently, there has been excitement about the promise of advances in protein structure modeling to generate hypotheses about sequence-structure-function relationships. Here, we leverage structural similarity to identify rapidly evolving proteasome assembly chaperones and characterize their function in Candidozyma (Candida) auris. Despite extensive sequence divergence, we demonstrate conservation of function, corroborating that specific folds, and not sequences, are required for function. This theoretical premise suggests that protein structures with certain properties should be functionally interchangeable, even if they were not products of a common evolutionary history. To reduce this theory to practice, we performed structure-informed protein design, exploring sequence space that is not accessible via stepwise evolution, and mutated more than 40 residues in the Poc4 proteasome assembly chaperone to demonstrate that artificial proteins can rescue complex biological processes in the context of the whole cell. This sequence-structure-function relationship expands our ability to use structure to identify deep evolutionary relationships between proteins and generate hypotheses about gene function in non-model organisms. Overall, this helps to define and understand functional constraints on protein evolution, with important implications for both future protein design and retrospective function prediction.
Earlier work (before 2026)
Siddiq MA, Duveau F, Wittkopp PJ
Nature Ecology & Evolution 8(12): 2184–2194 (2024)
Abstract
The environment influences how an organism’s genotype determines its phenotype and how this phenotype affects its fitness. Here, to better understand this dual role of environment in the production and selection of phenotypic variation, we determined genotype-phenotype-fitness relationships for mutant strains of Saccharomyces cerevisiae in four environments. Specifically, we measured how promoter mutations of the metabolic gene TDH3 modified expression level and affected growth for four different carbon sources. In each environment, we observed a clear relationship between TDH3 expression level and fitness, but this relationship differed among environments. Mutations with similar effects on expression in different environments often had different effects on fitness and vice versa. Such environment-specific relationships between phenotype and fitness can shape the evolution of phenotypic plasticity. We also found that mutations disrupting binding sites for transcription factors had more variable effects on expression among environments than those disrupting the TATA box, which is part of the core promoter. However, mutations with the most environmentally variable effects on fitness were located in the TATA box, because of both the lack of plasticity of TATA box mutations and environment-specific fitness functions. This observation suggests that mutations affecting different molecular mechanisms contribute unequally to regulatory sequence evolution in changing environments.
Vande Zande P, Siddiq MA, Hodgins-Davis A, Kim L, Wittkopp PJ
PLoS Genetics 19(12): e1011078 (2023)
Abstract
Genetic networks are surprisingly robust to perturbations caused by new mutations. This robustness is conferred in part by compensation for loss of a gene’s activity by genes with overlapping functions, such as paralogs. Compensation occurs passively when the normal activity of one paralog can compensate for the loss of the other, or actively when a change in one paralog’s expression, localization, or activity is required to compensate for loss of the other. The mechanisms of active compensation remain poorly understood in most cases. Here we investigate active compensation for the loss or reduction in expression of the Saccharomyces cerevisiae gene TDH3 by its paralog TDH2. TDH2 is upregulated in a dose-dependent manner in response to reductions in TDH3 by a mechanism requiring the shared transcriptional regulators Gcr1p and Rap1p. TDH1, a second and more distantly related paralog of TDH3, has diverged in its regulation and is upregulated by another mechanism. Other glycolytic genes regulated by Rap1p and Gcr1p show changes in expression similar to TDH2, suggesting that the active compensation by TDH3 paralogs is part of a broader homeostatic response mediated by shared transcriptional regulators.
Dulchavsky M, Mitra R, Wu K, Li J, Boer K, Liu X, Zhang Z, Vasquez C, Clark CT, Funckes K, Shankar K, Bonnet-Zahedi S, Siddiq M, Sepulveda Y, Suhandynata RT, Momper JD, Calabrese AN, George O, Stull F, Bardwell JCA
Nature Chemical Biology 19(11): 1406–1414 (2023)
Abstract
The flavoenzyme nicotine oxidoreductase (NicA2) is a promising injectable treatment to aid in the cessation of smoking, a behavior responsible for one in ten deaths worldwide. NicA2 acts by degrading nicotine in the bloodstream before it reaches the brain. Clinical use of NicA2 is limited by its poor catalytic activity in the absence of its natural electron acceptor CycN. Without CycN, NicA2 is instead oxidized slowly by dioxygen (O2), necessitating unfeasibly large doses in a therapeutic setting. Here, we report a genetic selection strategy that directly links CycN-independent activity of NicA2 to growth of Pseudomonas putida S16. This selection enabled us to evolve NicA2 variants with substantial improvement in their rate of oxidation by O2. The encoded mutations cluster around a putative O2 tunnel, increasing flexibility and accessibility to O2 in this region. These mutations further confer desirable clinical properties. A variant form of NicA2 is tenfold more effective than the wild type at degrading nicotine in the bloodstream of rats.
Raicu AM, Kadiyala D, Niblock M, Jain A, Yang Y, Bird KM, Bertholf K, Seenivasan A, Siddiq M, Arnosti DN
Molecular Biology and Evolution 40(2): msad003 (2023)
Abstract
Evolution of sequence-specific transcription factors clearly drives lineage-specific innovations, but less is known about how changes in the central transcriptional machinery may contribute to evolutionary transformations. In particular, transcriptional regulators are rich in intrinsically disordered regions that appear to be magnets for evolutionary innovation. The C-terminal Binding Protein (CtBP) is a transcriptional corepressor derived from an ancestral lineage of alpha hydroxyacid dehydrogenases; it is found in mammals and invertebrates, and features a core NAD-binding domain as well as an unstructured C-terminus (CTD) of unknown function. CtBP can act on promoters and enhancers to repress transcription through chromatin-linked mechanisms. Our comparative phylogenetic study shows that CtBP is a bilaterian innovation whose CTD of about 100 residues is present in almost all orthologs. CtBP CTDs contain conserved blocks of residues and retain a predicted disordered property, despite having variations in the primary sequence. Interestingly, the structure of the C-terminus has undergone radical transformation independently in certain lineages including flatworms and nematodes. Also contributing to CTD diversity is the production of myriad alternative RNA splicing products, including the production of “short” tailless forms of CtBP in Drosophila. Additional diversity stems from multiple gene duplications in vertebrates, where up to five CtBP orthologs have been observed. Vertebrate lineages show fewer major modifications in the unstructured CTD, possibly because gene regulatory constraints of the vertebrate body plan place specific constraints on this domain. Our study highlights the rich regulatory potential of this previously unstudied domain of a central transcriptional regulator.
Siddiq MA, Wittkopp PJ
Current Opinion in Genetics & Development 77: 101998 (2022)
Abstract
Studies of regulatory variation in yeast — at the level of new mutations, polymorphisms within a species, and divergence between species — have provided great insight into the molecular and evolutionary processes responsible for the evolution of gene expression in eukaryotes. The increasing ease with which yeast genomes can be manipulated and expression quantified in a high-throughput manner has recently accelerated mechanistic studies of cis– and trans-regulatory variation at multiple evolutionary timescales. These studies have, for example, identified differences in the properties of cis– and trans-acting mutations that affect their evolutionary fate, experimentally characterized the molecular mechanisms through which cis– and trans-regulatory variants act, and illustrated how regulatory networks can diverge between species with or without changes in gene expression.
Duveau F, Vande Zande P, Metzger BPH, Diaz CJ, Walker EA, Tryban S, Siddiq MA, Yang B, Wittkopp PJ
eLife 10: e67806 (2021)
Abstract
Heritable variation in a gene’s expression arises from mutations impacting cis– and trans-acting components of its regulatory network. Here, we investigate how trans-regulatory mutations are distributed within the genome and within a gene regulatory network by identifying and characterizing 69 mutations with trans-regulatory effects on expression of the same focal gene in Saccharomyces cerevisiae. Relative to 1766 mutations without effects on expression of this focal gene, we found that these trans-regulatory mutations were enriched in coding sequences of transcription factors previously predicted to regulate expression of the focal gene. However, over 90% of the trans-regulatory mutations identified mapped to other types of genes involved in diverse biological processes including chromatin state, metabolism, and signal transduction. These data show how genetic changes in diverse types of genes can impact a gene’s expression in trans, revealing properties of trans-regulatory mutations that provide the raw material for trans-regulatory variation segregating within natural populations.
Siddiq MA, Thornton JW
Proceedings of the National Academy of Sciences USA 116(43): 21634–21640 (2019)
Abstract
Polymorphism in the alcohol dehydrogenase (ADH) protein of Drosophila melanogaster, like genetic variation in many other enzymes, has long been hypothesized to be maintained by a selective trade-off between thermostability and enzyme activity. Two major Adh variants, named Fast and Slow, are distributed along latitudinal clines on several continents. The balancing selection trade-off hypothesis posits that Fast is favored at high latitudes because it metabolizes alcohol faster, whereas Slow is favored at low latitudes because it is more stable at high temperatures. Here we use biochemical and physiological assays of precisely engineered genetic variants to directly test this hypothesis. As predicted, the Fast protein has higher catalytic activity than Slow, and both the Fast protein and regulatory variants linked to it confer greater ethanol tolerance on transgenic animals. But we found no evidence of a temperature-mediated trade-off: The Fast protein is not less stable or active at high temperatures, and Fast alleles increase ethanol tolerance and survivorship at all temperatures tested. Further, analysis of a population genomic dataset reveals no signature of balancing selection in the Adh gene. These results provide strong evidence against balancing selection driven by a stability/activity trade-off in Adh and they justify caution about this hypothesis for other enzymes except those for which it has been directly tested. Our findings tentatively suggest that environment-specific selection for the Fast allele, coupled with demographic history, may have produced the observed pattern of Adh variation.
Siddiq MA, Hochberg GKA, Thornton JW
Current Opinion in Structural Biology 47: 113–122 (2017)
Abstract
Specific interactions between proteins and their molecular partners drive most biological processes, so understanding how these interactions evolve is an important question for biochemists and evolutionary biologists alike. It is often thought that ancestral proteins were systematically more promiscuous than modern proteins and that specificity usually evolves after gene duplication by partitioning and refining the activities of multifunctional ancestors. However, recent studies using ancestral protein reconstruction (APR) have found that ligand-specific functions in some modern protein families evolved de novo from ancestors that did not already have those functions. Further, the new specific interactions evolved by simple mechanisms, with just a few mutations changing classically recognized biochemical determinants of specificity, such as steric and electrostatic complementarity. Acquiring new specific interactions during evolution therefore appears to be neither difficult nor rare. Rather, it is likely that proteins continually gain and lose new activities over evolutionary time as mutations cause subtle but consequential changes in the shape and electrostatics of interaction interfaces. Only a few of these activities, however, are incorporated into the biological processes that contribute to fitness before they are lost to the ravages of further mutation.
Siddiq MA, Loehlin DW, Montooth KL, Thornton JW
Nature Ecology & Evolution 1(2): 25 (2017)
Abstract
Identifying the genetic basis for adaptive differences between species requires explicit tests of historical hypotheses concerning the effects of past changes in gene sequence on molecular function, organismal phenotype and fitness. We address this challenge by combining ancestral protein reconstruction with biochemical experiments and physiological analysis of transgenic animals that carry ancestral genes. We tested a widely held hypothesis of molecular adaptation — that changes in the alcohol dehydrogenase protein (ADH) along the lineage leading to Drosophila melanogaster increased the catalytic activity of the enzyme and thereby contributed to the ethanol tolerance and adaptation of the species to its ethanol-rich ecological niche. Our experiments strongly refute the predictions of the adaptive ADH hypothesis and caution against accepting intuitively appealing accounts of historical molecular adaptation that are based on correlative evidence. The experimental strategy we employed can be used to decisively test other adaptive hypotheses and the claims they entail about past biological causality.
Hoekstra LA, Siddiq MA, Montooth KL
Genetics 195(3): 1129–1139 (2013)
Abstract
Interactions between mitochondrial and nuclear gene products that underlie eukaryotic energy metabolism can cause the fitness effects of mutations in one genome to be conditional on variation in the other genome. In ectotherms, the effects of these interactions are likely to depend upon the thermal environment, because increasing temperature accelerates molecular rates. We find that temperature strongly modifies the pleiotropic phenotypic effects of an incompatible interaction between a Drosophila melanogaster polymorphism in the nuclear-encoded, mitochondrial tyrosyl-transfer (t)RNA synthetase and a D. simulans polymorphism in the mitochondrially encoded tRNATyr. The incompatible mitochondrial–nuclear genotype extends development time, decreases larval survivorship, and reduces pupation height, indicative of decreased energetic performance. These deleterious effects are ameliorated when larvae develop at 16° and exacerbated at warmer temperatures, leading to complete sterility in both sexes at 28°. The incompatible genotype has a normal metabolic rate at 16° but a significantly elevated rate at 25°, consistent with the hypothesis that inefficient energy metabolism extends development in this genotype at warmer temperatures. Furthermore, the incompatibility decreases metabolic plasticity of larvae developed at 16°, indicating that cooler development temperatures do not completely mitigate the deleterious effects of this genetic interaction. Our results suggest that the epistatic fitness effects of metabolic mutations may generally be conditional on the thermal environment. The expression of epistatic interactions in some environments, but not others, weakens the efficacy of selection in removing deleterious epistatic variants from populations and may promote the accumulation of incompatibilities whose fitness effects will depend upon the environment in which hybrids occur.
Meiklejohn CD, Holmbeck MA, Siddiq MA, Abt DN, Rand DM, Montooth KL
PLoS Genetics 9(1): e1003238 (2013)
Abstract
Mitochondrial transcription, translation, and respiration require interactions between genes encoded in two distinct genomes, generating the potential for mutations in nuclear and mitochondrial genomes to interact epistatically and cause incompatibilities that decrease fitness. Mitochondrial-nuclear epistasis for fitness has been documented within and between populations and species of diverse taxa, but rarely has the genetic or mechanistic basis of these mitochondrial-nuclear interactions been elucidated, limiting our understanding of which genes harbor variants causing mitochondrial-nuclear disruption and of the pathways and processes that are impacted by mitochondrial-nuclear coevolution. Here we identify an amino acid polymorphism in the Drosophila melanogaster nuclear-encoded mitochondrial tyrosyl-tRNA synthetase that interacts epistatically with a polymorphism in the D. simulans mitochondrial-encoded tRNATyr to significantly delay development, compromise bristle formation, and decrease fecundity. The incompatible genotype specifically decreases the activities of oxidative phosphorylation complexes I, III, and IV that contain mitochondrial-encoded subunits. Combined with the identity of the interacting alleles, this pattern indicates that mitochondrial protein translation is affected by this interaction. Our findings suggest that interactions between mitochondrial tRNAs and their nuclear-encoded tRNA synthetases may be targets of compensatory molecular evolution. Human mitochondrial diseases are often genetically complex and variable in penetrance, and the mitochondrial-nuclear interaction we document provides a plausible mechanism to explain this complexity.