Epistatic Relationship Research Articles

Synergistic effects of tRNA modification defects in Escherichia coli K12.

tRNAs are the central adaptor molecule in translation and require a wide variety of post-transcriptional modifications to fulfill their functions. The model gram-negative Escherichia coli K-12 is one of the handfuls of bacteria where all the tRNA modifications and corresponding genes have been characterized. This work dissects epistatic relationships between tRNA modification genes in E. coli by conducting a synthetic lethal screen revealing 5 pairs of modifications that cannot be deleted in combination when cells are grown in rich media, and 15 pairs of modifications that lead to growth defects when deleted in combination. The truA gene involved in the insertion of Psi residues at positions 38 to 40 in multiple tRNAs gave the highest number of synthetic lethality phenotypes that could be complemented by the expression of the truA gene in trans and in some cases suppressed by the overexpression of target tRNAs. A pilot phenotype screen of strains lacking two tRNA modifications genes viable in rich media lead to the identification of growth conditions that exhibited poor growth. This work lays the foundation to dissect the role of modifications in tRNA quality control.

Combinatorial effector targeting (COMET) for transcriptional modulation and locus-specific biochemistry.

Understanding how human gene expression is coordinately regulated by functional units of proteins across the genome remains a major biological goal. Here, we present COMET, a high-throughput screening platform for combinatorial effector targeting for the identification of transcriptional modulators. We generate libraries of combinatorial dCas9-based fusion proteins, containing two to six effector domains, allowing us to systematically investigate more than 110,000 combinations of effector proteins at endogenous human loci for their influence on transcription. Importantly, we keep full proteins or domains intact, maintaining catalytic cores and surfaces for protein-protein interactions. We observe more than 5800 significant hits that modulate transcription, we demonstrate cell type specific transcriptional modulation, and we further investigate epistatic relationships between our effector combinations. We validate unexpected combinations as synergistic or buffering, emphasizing COMET as both a method for transcriptional effector discovery, and as a functional genomics tool for identifying novel domain interactions and directing locus-specific biochemistry.

Chromatin remodeller Chd7 is developmentally regulated in the neural crest by tissue-specific transcription factors.

Neurocristopathies such as CHARGE syndrome result from aberrant neural crest development. A large proportion of CHARGE cases are attributed to pathogenic variants in the gene encoding CHD7, chromodomain helicase DNA binding protein 7, which remodels chromatin. While the role for CHD7 in neural crest development is well documented, how this factor is specifically up-regulated in neural crest cells is not understood. Here, we use epigenomic profiling of chick and human neural crest to identify a cohort of enhancers regulating Chd7 expression in neural crest cells and other tissues. We functionally validate upstream transcription factor binding at candidate enhancers, revealing novel epistatic relationships between neural crest master regulators and Chd7, showing tissue-specific regulation of a globally acting chromatin remodeller. Furthermore, we find conserved enhancer features in human embryonic epigenomic data and validate the activity of the human equivalent CHD7 enhancers in the chick embryo. Our findings embed Chd7 in the neural crest gene regulatory network and offer potentially clinically relevant elements for interpreting CHARGE syndrome cases without causative allocation.

The epistatic relationship of Drosophila melanogaster CtIP and Rif1 in homology-directed repair of DNA double-strand breaks.

Double-strand breaks (DSBs) are genotoxic DNA lesions that pose significant threats to genomic stability, necessitating precise and efficient repair mechanisms to prevent cell death or mutations. DSBs are repaired through nonhomologous end-joining (NHEJ) or homology-directed repair (HDR), which includes homologous recombination (HR) and single-strand annealing (SSA). CtIP and Rif1 are conserved proteins implicated in DSB repair pathway choice, possibly through redundant roles in promoting DNA end-resection required for HDR. Although the roles of these proteins have been well-established in other organisms, the role of Rif1 and its potential redundancies with CtIP in Drosophila melanogaster remain elusive. To examine the roles of DmCtIP and DmRif1 in DSB repair, this study employed the direct repeat of white (DR-white) assay, tracking across indels by decomposition (TIDE) analysis, and P{wIw_2 kb 3'} assay to track repair outcomes in HR, NHEJ, and SSA, respectively. These experiments were performed in DmCtIPΔ/Δ single mutants, DmRif1Δ/Δ single mutants, and DmRif1Δ/Δ; DmCtIPΔ/Δ double mutants. This work demonstrates significant defects in both HR and SSA repair in DmCtIPΔ/Δ and DmRif1Δ/Δ single mutants. However, experiments in DmRif1Δ/Δ; DmCtIPΔ/Δ double mutants reveal that DmCtIP is epistatic to DmRif1 in promoting HDR. Overall, this study concludes that DmRif1 and DmCtIP do not perform their activities in a redundant pathway, but rather DmCtIP is the main driver in promoting HR and SSA, most likely through its role in end resection.

Addressing epistasis in the design of protein function

Mutations in protein active sites can dramatically improve function. The active site, however, is densely packed and extremely sensitive to mutations. Therefore, some mutations may only be tolerated in combination with others in a phenomenon known as epistasis. Epistasis reduces the likelihood of obtaining improved functional variants and dramatically slows natural and lab evolutionary processes. Research has shed light on the molecular origins of epistasis and its role in shaping evolutionary trajectories and outcomes. In addition, sequence- and AI-based strategies that infer epistatic relationships from mutational patterns in natural or experimental evolution data have been used to design functional protein variants. In recent years, combinations of such approaches and atomistic design calculations have successfully predicted highly functional combinatorial mutations in active sites. These were used to design thousands of functional active-site variants, demonstrating that, while our understanding of epistasis remains incomplete, some of the determinants that are critical for accurate design are now sufficiently understood. We conclude that the space of active-site variants that has been explored by evolution may be expanded dramatically to enhance natural activities or discover new ones. Furthermore, design opens the way to systematically exploring sequence and structure space and mutational impacts on function, deepening our understanding and control over protein activity.

Molecular definition of the endogenous Toll-like receptor signalling pathways.

Innate immune pattern recognition receptors, such as the Toll-like receptors (TLRs), are key mediators of the immune response to infection and central to our understanding of health and disease1. After microbial detection, these receptors activate inflammatory signal transduction pathways that involve IκB kinases, mitogen-activated protein kinases, ubiquitin ligases and other adaptor proteins. The mechanisms that connect the proteins in the TLR pathways are poorly defined. To delineate TLR pathway activities, we engineered macrophages to enable microscopy and proteomic analysis of the endogenous myddosome constituent MyD88. We found that myddosomes form transient contacts with activated TLRs and that TLR-free myddosomes are dynamic in size, number and composition over the course of 24 h. Analysis using super-resolution microscopy revealed that, within most myddosomes, MyD88 forms barrel-like structures that function as scaffolds for effector protein recruitment. Proteomic analysis demonstrated that myddosomes contain proteins that act at all stages and regulate all effector responses of the TLR pathways, and genetic analysis defined the epistatic relationship between these effector modules. Myddosome assembly was evident in cells infected with Listeria monocytogenes, but these bacteria evaded myddosome assembly and TLR signalling during cell-to-cell spread. On the basis of these findings, we propose that the entire TLR signalling pathway is executed from within the myddosome.

Differential requirement for RecFOR pathway components in Thermus thermophilus.

Recombinational repair is an important mechanism that allows DNA replication to overcome damaged templates, so the DNA is duplicated timely and correctly. The RecFOR pathway is one of the common ways to load RecA, while the RuvABC complex operates in the resolution of DNA intermediates. We have generated deletions of recO, recR and ruvB genes in Thermus thermophilus, while a recF null mutant could not be obtained. The recO deletion was in all cases accompanied by spontaneous loss of function mutations in addA or addB genes, which encode a helicase-exonuclease also key for recombination. The mutants were moderately affected in viability and chromosome segregation. When we generated these mutations in a Δppol/addAB strain, we observed that the transformation efficiency was maintained at the typical level of Δppol/addAB, which is 100-fold higher than that of the wild type. Most mutants showed increased filamentation phenotypes, especially ruvB, which also had DNA repair defects. These results suggest that in T. thermophilus (i) the components of the RecFOR pathway have differential roles, (ii) there is an epistatic relationship of the AddAB complex over the RecFOR pathway and (iii) that neither of the two pathways or their combination is strictly required for viability although they are necessary for normal DNA repair and chromosome segregation.

Genome-scale analysis of interactions between genetic perturbations and natural variation

Interactions between genetic perturbations and segregating loci can cause perturbations to show different phenotypic effects across genetically distinct individuals. To study these interactions on a genome scale in many individuals, we used combinatorial DNA barcode sequencing to measure the fitness effects of 8046 CRISPRi perturbations targeting 1721 distinct genes in 169 yeast cross progeny (or segregants). We identified 460 genes whose perturbation has different effects across segregants. Several factors caused perturbations to show variable effects, including baseline segregant fitness, the mean effect of a perturbation across segregants, and interacting loci. We mapped 234 interacting loci and found four hub loci that interact with many different perturbations. Perturbations that interact with a given hub exhibit similar epistatic relationships with the hub and show enrichment for cellular processes that may mediate these interactions. These results suggest that an individual’s response to perturbations is shaped by a network of perturbation-locus interactions that cannot be measured by approaches that examine perturbations or natural variation alone.

ETV2 induces endothelial, but not hematopoietic, lineage specification in birds.

Cardiovascular system develops from the lateral plate mesoderm. Its three primary cell lineages (hematopoietic, endothelial, and muscular) are specified by the sequential actions of conserved transcriptional factors. ETV2, a master regulator of mammalian hemangioblast development, however, is absent in the chicken genome and acts downstream of NPAS4L in zebrafish. Here, we investigated the epistatic relationship between NPAS4L and ETV2 in avian hemangioblast development. We showed that ETV2 is deleted in all 363 avian genomes analyzed. Mouse ETV2 induced LMO2, but not NPAS4L or SCL, expression in chicken mesoderm. Squamate (lizards, geckos, and snakes) genomes contain both NPAS4L and ETV2 In Madagascar ground gecko, both genes were expressed in developing hemangioblasts. Gecko ETV2 induced only LMO2 in chicken mesoderm. We propose that both NPAS4L and ETV2 were present in ancestral amniote, with ETV2 acting downstream of NPAS4L in endothelial lineage specification. ETV2 may have acted as a pioneer factor by promoting chromatin accessibility of endothelial-specific genes and, in parallel with NPAS4L loss in ancestral mammals, has gained similar function in regulating blood-specific genes.

Genomic Prediction from Multi-Environment Trials of Wheat Breeding.

Genomic prediction relates a set of markers to variability in observed phenotypes of cultivars and allows for the prediction of phenotypes or breeding values of genotypes on unobserved individuals. Most genomic prediction approaches predict breeding values based solely on additive effects. However, the economic value of wheat lines is not only influenced by their additive component but also encompasses a non-additive part (e.g., additive × additive epistasis interaction). In this study, genomic prediction models were implemented in three target populations of environments (TPE) in South Asia. Four models that incorporate genotype × environment interaction (G × E) and genotype × genotype (GG) were tested: Factor Analytic (FA), FA with genomic relationship matrix (FA + G), FA with epistatic relationship matrix (FA + GG), and FA with both genomic and epistatic relationship matrices (FA + G + GG). Results show that the FA + G and FA + G + GG models displayed the best and a similar performance across all tests, leading us to infer that the FA + G model effectively captures certain epistatic effects. The wheat lines tested in sites in different TPE were predicted with different precisions depending on the cross-validation employed. In general, the best prediction accuracy was obtained when some lines were observed in some sites of particular TPEs and the worse genomic prediction was observed when wheat lines were never observed in any site of one TPE.

Combinatorial transcriptomic and genetic dissection of insulin/IGF-1 signaling-regulated longevity in Caenorhabditis elegans.

Classical genetic analysis is invaluable for understanding the genetic interactions underlying specific phenotypes, but requires laborious and subjective experiments to characterize polygenic and quantitative traits. Contrarily, transcriptomic analysis enables the simultaneous and objective identification of multiple genes whose expression changes are associated with specific phenotypes. Here, we conducted transcriptomic analysis of genes crucial for longevity using datasets with daf-2/insulin/IGF-1 receptor mutant Caenorhabditis elegans. Our analysis unraveled multiple epistatic relationships at the transcriptomic level, in addition to verifying genetically established interactions. Our combinatorial analysis also revealed transcriptomic changes associated with longevity conferred by daf-2 mutations. In particular, we demonstrated that the extent of lifespan changes caused by various mutant alleles of the longevity transcription factor daf-16/FOXO matched their effects on transcriptomic changes in daf-2 mutants. We identified specific aging-regulating signaling pathways and subsets of structural and functional RNA elements altered by different genes in daf-2 mutants. Lastly, we elucidated the functional cooperation between several longevity regulators, based on the combination of transcriptomic and molecular genetic analysis. These data suggest that different biological processes coordinately exert their effects on longevity in biological networks. Together our work demonstrates the utility of transcriptomic dissection analysis for identifying important genetic interactions for physiological processes, including aging and longevity.

Genome-Enabled Analysis of Population Dynamics and Virulence-Associated Loci in the Oat Crown Rust Fungus Puccinia coronata f. sp. avenae.

Puccinia coronata f. sp. avenae (Pca) is an important fungal pathogen causing crown rust that impacts oat production worldwide. Genetic resistance for crop protection against Pca is often overcome by the rapid virulence evolution of the pathogen. This study investigated the factors shaping adaptive evolution of Pca using pathogen populations from distinct geographic regions within the United States and South Africa. Phenotypic and genome-wide sequencing data of these diverse Pca collections, including 217 isolates, uncovered phylogenetic relationships and established distinct genetic composition between populations from northern and southern regions from the United States and South Africa. The population dynamics of Pca involve a bidirectional movement of inoculum between northern and southern regions of the United States and contributions from clonality and sexuality. The population from South Africa is solely clonal. A genome-wide association study (GWAS) employing a haplotype-resolved Pca reference genome was used to define 11 virulence-associated loci corresponding to 25 oat differential lines. These regions were screened to determine candidate Avr effector genes. Overall, the GWAS results allowed us to identify the underlying genetic factors controlling pathogen recognition in an oat differential set used in the United States to assign pathogen races (pathotypes). Key GWAS findings support complex genetic interactions in several oat lines, suggesting allelism among resistance genes or redundancy of genes included in the differential set, multiple resistance genes recognizing genetically linked Avr effector genes, or potentially epistatic relationships. A careful evaluation of the composition of the oat differential set accompanied by the development or implementation of molecular markers is recommended. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Interaction models matter: an efficient, flexible computational framework for model-specific investigation of epistasis

PurposeEpistasis, the interaction between two or more genes, is integral to the study of genetics and is present throughout nature. Yet, it is seldom fully explored as most approaches primarily focus on single-locus effects, partly because analyzing all pairwise and higher-order interactions requires significant computational resources. Furthermore, existing methods for epistasis detection only consider a Cartesian (multiplicative) model for interaction terms. This is likely limiting as epistatic interactions can evolve to produce varied relationships between genetic loci, some complex and not linearly separable.MethodsWe present new algorithms for the interaction coefficients for standard regression models for epistasis that permit many varied models for the interaction terms for loci and efficient memory usage. The algorithms are given for two-way and three-way epistasis and may be generalized to higher order epistasis. Statistical tests for the interaction coefficients are also provided. We also present an efficient matrix based algorithm for permutation testing for two-way epistasis. We offer a proof and experimental evidence that methods that look for epistasis only at loci that have main effects may not be justified. Given the computational efficiency of the algorithm, we applied the method to a rat data set and mouse data set, with at least 10,000 loci and 1,000 samples each, using the standard Cartesian model and the XOR model to explore body mass index.ResultsThis study reveals that although many of the loci found to exhibit significant statistical epistasis overlap between models in rats, the pairs are mostly distinct. Further, the XOR model found greater evidence for statistical epistasis in many more pairs of loci in both data sets with almost all significant epistasis in mice identified using XOR. In the rat data set, loci involved in epistasis under the XOR model are enriched for biologically relevant pathways.ConclusionOur results in both species show that many biologically relevant epistatic relationships would have been undetected if only one interaction model was applied, providing evidence that varied interaction models should be implemented to explore epistatic interactions that occur in living systems.

Determining Photoreceptor Cell Identity: Rod Versus Cone Fate Governed by tbx2b Opposing nrl.

NRL is an influential transcription factor and central to animal modeling in ophthalmology. Disrupting NRL abrogates rod development and produces an excess of S-cones (also known as "UV cones" or "short-wavelength-sensitive1 [SWS1] cones"). Strikingly, mutations in zebrafish tbx2b produce the exact opposite phenotypes (excess rods and loss of SWS1 cones). We sought to define what genetic relationship exists, if any, between these transcription factors. We also infer whether these two phenotypes (altered rod abundance and altered SWS1 cone abundance) are independent versus inter-related. Zebrafish mutants were bred to disrupt nrl and tbx2b in concert. Rods and SWS1 cones were quantified and characterized at ultrastructural and transcriptional levels. Considering single mutant zebrafish, we confirmed previously established phenotypes and noted that the number of rods lost in nrl-/- mutants is reflected by a concomitant increase in SWS1 cone abundance. The tbx2b-/- mutants present the opposite phenotype(s) but exhibit a similar trade-off in cell abundances, with lots of rods and a concomitant decrease in SWS1 cones. Double mutant nrl-/-;tbx2b-/- zebrafish recapitulate the nrl-/- mutant phenotype(s). The tbx2b is thought to be required for producing SWS1 cones in zebrafish, but this can be over-ridden when nrl is absent. Regarding the altered cell abundances observed in either tbx2b-/- or nrl-/- mutants, the alterations in rod and SWS1 cones appear to not be two separate phenotypes but are instead a single intertwined outcome. The tbx2b and nrl are in an epistatic relationship, with nrl phenotypes dominating, implying that tbx2b is upstream of nrl in photoreceptor cell fate determination.

Abstract B017: VCP extracts the chromatin remodeler SNF2H from nascent DNA to stabilize stressed replication forks

Abstract Valosin-containing protein (VCP) is an evolutionarily conserved AAA+ ATPase which plays pleiotropic roles in global proteostasis by extracting polyubiquitinated proteins from various cellular organelles and structures and facilitating their turnover. Our recent work demonstrated that nuclear VCP activity during DNA damage response specifically depends on ATM/ATR/DNA-PK-mediated Ser784 phosphorylation. Ser784phosphorylation increases VCP ability to facilitate chromatin-associated protein degradation and is required for DNA repair, checkpoint signaling, and cell survival. Clinically, high nuclear pSer784-VCP levels are significantly associated with poor survival outcome among chemotherapy-treated breast and pancreatic cancer patients. Thus, pSer784-VCP is a broadly important and therapeutically relevant genome stabilizer. Here, we focused specifically on the role of pSer784-VCP in DNA replication fork stability during stress, a known source of DNA damage if not properly regulated. By performing single-molecule DNA fiber analysis using human cell lines with genetic gain and loss-of-function of VCP, we observed that VCP protects nascent DNA from hydroxyurea (HU)-induced over-resection and this activity requires Ser784 phosphorylation. Using proximity ligation assay (PLA), we detected HU-induced pSer784-VCP on nascent DNA, suggesting that it may protect stressed forks by extracting functionally important regulators. Consistent with this theory, we detected an interaction of pSer784-VCP with SNF2H, a nucleosome-sliding enzyme which we recently found essential for nascent DNA over-resection at unprotected forks. In support of the hypothesis that SNF2H is a novel substrate of pSer784-VCP at stressed forks, we observed SNF2H retention on nascent DNA upon VCP knockdown in HU-treated cells, and the functional rescue of SNF2H retention by VCP re-expression requires Ser784phosphorylation. In support of their epistatic relationship, SNF2H depletion fully rescued fork over-resection caused by VCP loss. Similarly, blocking fork reversal by depleting SMARCAL1 or FBH1 also rescued fork over-resection caused by VCP loss. Importantly, mutating the ATPase activity of SNF2H abolished its ability to promote fork over-resection, suggesting that nucleosome remodeling by SNF2H underlies its effect on fork stability and may be a key molecular event enabling fork reversal upstream of fork degradation. Collectively, our data suggest that pSer784-VCP is a previously unrecognized stabilizer of stressed DNA replication forks, and it does so at least in part via its physical extraction of SNF2H from nascent DNA and the consequent attenuation of nucleosome remodeling needed for fork reversal and over-resection. More experiments are currently underway to further dissect the mechanistic aspects of this working model as well as to evaluate the potential of developing pSer784-VCP into a clinically relevant chemo-predictive biomarker and chemo-sensitizing target for different types of cancer. Citation Format: Jieya Shao, Mari Iwase, Rong Xu, Shuyang Lin. VCP extracts the chromatin remodeler SNF2H from nascent DNA to stabilize stressed replication forks [abstract]. In: Proceedings of the AACR Special Conference in Cancer Research: DNA Damage Repair: From Basic Science to Future Clinical Application; 2024 Jan 9-11; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2024;84(1 Suppl):Abstract nr B017.

Complex epistatic interactions between ELF3, PRR9, and PRR7 regulate the circadian clock and plant physiology.

Circadian clocks are endogenous timekeeping mechanisms that coordinate internal physiological responses with the external environment. EARLY FLOWERING3 (ELF3), PSEUDO RESPONSE REGULATOR (PRR9), and PRR7 are essential components of the plant circadian clock and facilitate entrainment of the clock to internal and external stimuli. Previous studies have highlighted a critical role for ELF3 in repressing the expression of PRR9 and PRR7. However, the functional significance of activity in regulating circadian clock dynamics and plant development is unknown. To explore this regulatory dynamic further, we first employed mathematical modeling to simulate the effect of the prr9/prr7 mutation on the elf3 circadian phenotype. These simulations suggested that simultaneous mutations in prr9/prr7 could rescue the elf3 circadian arrhythmia. Following these simulations, we generated all Arabidopsis elf3/prr9/prr7 mutant combinations and investigated their circadian and developmental phenotypes. Although these assays could not replicate the results from the mathematical modeling, our results have revealed a complex epistatic relationship between ELF3 and PRR9/7 in regulating different aspects of plant development. ELF3 was essential for hypocotyl development under ambient and warm temperatures, while PRR9 was critical for root thermomorphogenesis. Finally, mutations in prr9 and prr7 rescued the photoperiod-insensitive flowering phenotype of the elf3 mutant. Together, our results highlight the importance of investigating the genetic relationship among plant circadian genes.

A dual sgRNA library design to probe genetic modifiers using genome-wide CRISPRi screens

Mapping genetic interactions is essential for determining gene function and defining novel biological pathways. We report a simple to use CRISPR interference (CRISPRi) based platform, compatible with Fluorescence Activated Cell Sorting (FACS)-based reporter screens, to query epistatic relationships at scale. This is enabled by a flexible dual-sgRNA library design that allows for the simultaneous delivery and selection of a fixed sgRNA and a second randomized guide, comprised of a genome-wide library, with a single transduction. We use this approach to identify epistatic relationships for a defined biological pathway, showing both increased sensitivity and specificity than traditional growth screening approaches.

The Saccharomyces cerevisiae acetyltransferase Gcn5 exerts antagonistic pleiotropic effects on chronological ageing.

Compared to replicative lifespan, epigenetic regulation of chronological lifespan (CLS) is less well understood in yeast. Here, by screening all the viable mutants of histone acetyltransferase (HAT) and histone deacetylase (HDAC), we demonstrate that Gcn5, functioning in the HAT module of the SAGA/SLIK complex, exhibits an epistatic relationship with the HDAC Hda1 to control the expression of starvation-induced stress response and respiratory cell growth. Surprisingly, the gcn5Δ mutants lose their colony-forming potential early in the stationary phase but display a longer maximum CLS than their WT counterparts, suggesting the contradictory roles of Gcn5 in lifespan regulation. Integrative analyses of the transcriptome, metabolome and ChIP assays reveal that Gcn5 is necessary for the activation of two regulons upon glucose starvation: the Msn2/4-/Gis1-dependent stress response and the Cat8-/Adr1-mediated metabolic reprogramming, to enable pro-longevity characteristics, including redox homeostasis, stress resistance and maximal storage of carbohydrates. The activation of Cat8-/Adr1-dependent regulon also promotes the pyruvate dehydrogenase (PDH) bypass, leading to acetyl-CoA synthesis, global and targeted H3K9 acetylation. Global H3K9 acetylation levels mediated by Gcn5 and Hda1 during the transition into stationary phase are positively correlated with senescent cell populations accumulated in the aged cell cultures. These data suggest that Gcn5 lies in the centre of a feed-forward loop between histone acetylation and starvation-induced gene expression, enabling stress resistance and homeostasis but also promoting chronological ageing concomitantly.

Idiosyncratic and dose-dependent epistasis drives variation in tomato fruit size.

Epistasis between genes is traditionally studied with mutations that eliminate protein activity, but most natural genetic variation is in cis-regulatory DNA and influences gene expression and function quantitatively. In this study, we used natural and engineered cis-regulatory alleles in a plant stem-cell circuit to systematically evaluate epistatic relationships controlling tomato fruit size. Combining a promoter allelic series with two other loci, we collected over 30,000 phenotypic data points from 46 genotypes to quantify how allele strength transforms epistasis. We revealed a saturating dose-dependent relationship but also allele-specific idiosyncratic interactions, including between alleles driving a step change in fruit size during domestication. Our approach and findings expose an underexplored dimension of epistasis, in which cis-regulatory allelic diversity within gene regulatory networks elicits nonlinear, unpredictable interactions that shape phenotypes.

Impact of piperaquine resistance in Plasmodium falciparum on malaria treatment effectiveness in The Guianas: a descriptive epidemiological study

Plasmodium falciparum is an apicomplexan parasite responsible for lethal cases of malaria. According to WHO recommendations, P falciparum cases are treated with artemisinin-based combination therapy including dihydroartemisinin-piperaquine. However, the emergence of resistant parasites against dihydroartemisinin-piperaquine was reported in southeast Asia in 2008 and, a few years later, suspected in South America. To characterise resistance emergence, a treatment efficacy study was performed on the reported patients infected with P falciparum and treated with dihydroartemisinin-piperaquine in French Guiana (n=6, 2016-18). Contemporary isolates collected in French Guiana were genotyped for P falciparum chloroquine resistance transporter (pfCRT; n=845) and pfpm2 and pfpm3 copy number (n=231), phenotyped using the in vitro piperaquine survival assay (n=86), and analysed through genomic studies (n=50). Additional samples from five Amazonian countries and one outside the region were genotyped (n=1440). In field isolates, 40 (47%) of 86 (95% CI 35·9-57·1) were resistant to piperaquine in vitro; these phenotypes were more associated with pfCRTC350R (ie, Cys350Arg) and pfpm2 and pfpm3 amplifications (Dunn test, p<0·001). Those markers were also associated with dihydroartemisinin-piperaquine treatment failure (n=3 [50%] of 6). A high prevalence of piperaquine resistance markers was observed in Suriname in 19 (83%) of 35 isolates and in Guyana in 579 (73%) of 791 isolates. The pfCRTC350R mutation emerged before pfpm2 and pfpm3 amplification in a temporal sequence different from southeast Asia, and in the absence of artemisinin partial resistance, suggesting a geographically distinctive epistatic relationship between these genetic markers. The high prevalence of piperaquine resistance markers in parasite populations of the Guianas, and the risk of associated therapeutic failures calls for caution on dihydroartemisinin-piperaquine use in the region. Furthermore, greater attention should be given to potential differences in genotype to phenotype mapping across genetically distinct parasite populations from different continents. Pan American Health Organization and WHO, French Ministry for Research, European Commission, Santé publique France, Agence Nationale de la Recherche, Fundação de Amparo à Pesquisa do Estado do Amazonas, Ministry of Health of Brazil, Oswaldo Cruz Foundation, and National Institutes of Health. For the French and Portuguese translations of the abstract see Supplementary Materials section.