Search

Integrative structure modeling is increasingly used for determining the architectures of biological assemblies, especially those that are structurally heterogeneous. Recently, we reported on how to convert in vivo genetic interaction measurements into spatial restraints for structural modeling: first, phenotypic profiles are generated for each point mutation and thousands of gene deletions or environmental perturbations. Following, the phenotypic profile similarities are converted into distance restraints on the pairs of mutated residues. We illustrate the approach by determining the structure of the histone H3-H4 complex. The method is implemented in our open-source IMP program, expanding the structural biology toolbox by allowing structural characterization based on in vivo data without the need to purify the target system. We compare genetic interaction measurements to other sources of structural information, such as residue coevolution and deep-learning structure prediction of complex subunits. We also suggest that determining genetic interactions could benefit from new technologies, such as CRISPR-Cas9 approaches to gene editing, especially for mammalian cells. Finally, we highlight the opportunity for using genetic interactions to determine recalcitrant biomolecular structures, such as those of disordered proteins, transient protein assemblies, and host-pathogen protein complexes.

Seasonal shedding of winter hair at the start of summer is well studied in wild and domesticated populations. However, the genetic influences on this trait and their interactions are poorly understood. We use data from 13,364 cattle with 36,899 repeated phenotypes to investigate the relationship between hair shedding and environmental variables, single nucleotide polymorphisms, and their interactions to understand quantitative differences in seasonal shedding. Using deregressed estimated breeding values from a repeated records model in a genome-wide association analysis (GWAA) and meta-analysis of year-specific GWAA gave remarkably similar results. These GWAA identified hundreds of variants associated with seasonal hair shedding. There were especially strong associations between chromosomes 5 and 23. Genotype-by-environment interaction GWAA identified 1,040 day length-by-genotype interaction associations and 17 apparent temperature-by-genotype interaction associations with hair shedding, highlighting the importance of day length on hair shedding. Accurate genomic predictions of hair shedding were created for the entire dataset, Angus, Hereford, Brangus, and multibreed datasets. Loci related to metabolism and light-sensing have a large influence on seasonal hair shedding. This is one of the largest genetic analyses of a phenological trait and provides insight into both agriculture production and basic science. [ABSTRACT FROM AUTHOR]

Mediator of ERBB2-driven cell motility 1 (MEMO1) is an evolutionary conserved protein implicated in many biological processes; however, its primary molecular function remains unknown. Importantly, MEMO1 is overexpressed in many types of cancer and was shown to modulate breast cancer metastasis through altered cell motility. To better understand the function of MEMO1 in cancer cells, we analyzed genetic interactions of MEMO1 using gene essentiality data from 1028 cancer cell lines and found multiple iron-related genes exhibiting genetic relationships with MEMO1. We experimentally confirmed several interactions between MEMO1 and iron-related proteins in living cells, most notably, transferrin receptor 2 (TFR2), mitoferrin-2 (SLC25A28), and the global iron response regulator IRP1 (ACO1). These interactions indicate that cells with high-MEMO1 expression levels are hypersensitive to the disruptions in iron distribution. Our data also indicate that MEMO1 is involved in ferroptosis and is linked to iron supply to mitochondria. We have found that purified MEMO1 binds iron with high affinity under redox conditions mimicking intracellular environment and solved MEMO1 structures in complex with iron and copper. Our work reveals that the iron coordination mode in MEMO1 is very similar to that of iron-containing extradiol dioxygenases, which also display a similar structural fold. We conclude that MEMO1 is an iron-binding protein that modulates iron homeostasis in cancer cells.

Genome-wide association studies and experimental mouse models implicate the MIB1 and GATA6 genes in congenital heart disease (CHD). Their close physical proximity and conserved synteny suggest that these two genes might be involved in analogous cardiac developmental processes. Heterozygous Gata6 loss-of-function mutations alone or humanized Mib1 mutations in a NOTCH1-sensitized genetic background cause bicuspid aortic valve (BAV) and a membranous ventricular septal defect (VSD), consistent with MIB1 and NOTCH1 functioning in the same pathway. To determine if MIB1-NOTCH and GATA6 interact in valvular and septal development, we generated compound heterozygote mice carrying different Mib1 missense (Mib1K735R and Mib1V943F) or nonsense (Mib1R530X) mutations with the Gata6STOP/+ heterozygous null mutation. Combining Mib1R530X/+ or Mib1K735R/+ with Gata6STOP/+ does not affect Gata6STOP/+ single mutant phenotypes. In contrast, combining Mib1V943F/+ with Gata6STOP/+ decreases the incidence of BAV and VSD by 50%, suggesting a suppressive effect of Mib1V943F/+ on Gata6STOP/+. Transcriptomic and functional analyses revealed that while the EMT pathway term is depleted in the Gata6STOP/+ mutant, introducing the Mib1V943F variant robustly enriches this term, consistent with the Mib1V943F/+ phenotypic suppression of Gata6STOP/+. Interestingly, combined Notch1 and Gata6 insufficiency led to a nearly fully penetrant VSD but did not affect the BAV phenotype, underscoring the complex functional relationship between MIB1, NOTCH, and GATA6 in valvular and septal development.

Abstract Background Understanding the impact of gene interactions on disease phenotypes is increasingly recognised as a crucial aspect of genetic disease research. This trend is reflected by the growing amount of clinical research on oligogenic diseases, where disease manifestations are influenced by combinations of variants on a few specific genes. Although statistical machine-learning methods have been developed to identify relevant genetic variant or gene combinations associated with oligogenic diseases, they rely on abstract features and black-box models, posing challenges to interpretability for medical experts and impeding their ability to comprehend and validate predictions. In this work, we present a novel, interpretable predictive approach based on a knowledge graph that not only provides accurate predictions of disease-causing gene interactions but also offers explanations for these results. Results We introduce BOCK, a knowledge graph constructed to explore disease-causing genetic interactions, integrating curated information on oligogenic diseases from clinical cases with relevant biomedical networks and ontologies. Using this graph, we developed a novel predictive framework based on heterogenous paths connecting gene pairs. This method trains an interpretable decision set model that not only accurately predicts pathogenic gene interactions, but also unveils the patterns associated with these diseases. A unique aspect of our approach is its ability to offer, along with each positive prediction, explanations in the form of subgraphs, revealing the specific entities and relationships that led to each pathogenic prediction. Conclusion Our method, built with interpretability in mind, leverages heterogenous path information in knowledge graphs to predict pathogenic gene interactions and generate meaningful explanations. This not only broadens our understanding of the molecular mechanisms underlying oligogenic diseases, but also presents a novel application of knowledge graphs in creating more transparent and insightful predictors for genetic research.

Although it is well known that epistasis plays an important role in many evolutionary processes (e.g., speciation, evolution of sex), our knowledge on the frequency and prevalent sign of epistatic interactions is mainly limited to unicellular organisms or cell cultures of multicellular organisms. This is even more pronounced in regard to how the environment can influence genetic interactions. To broaden our knowledge in that respect we studied gene–gene interactions in a whole multicellular organism, Caenorhabditis elegans. We screened over one thousand gene interactions, each one in standard laboratory conditions, and under three different stressors: heat shock, oxidative stress, and genotoxic stress. Depending on the condition, between 7% and 22% of gene pairs showed significant genetic interactions and an overall sign of epistasis changed depending on the condition. Sign epistasis was quite common, but reciprocal sign epistasis was extremally rare. One interaction was common to all conditions, whereas 78% of interactions were specific to only one environment. Although epistatic interactions are quite common, their impact on evolutionary processes will strongly depend on environmental factors. [ABSTRACT FROM AUTHOR]

The selective autophagic degradation of mitochondria via mitophagy is essential for preserving mitochondrial homeostasis and, thereby, disease maintenance and progression in acute myeloid leukemia (AML). Mitophagy is orchestrated by a variety of mitophagy receptors whose interplay is not well understood. Here, we established a pairwise multiplexed CRISPR screen targeting mitophagy receptors to elucidate redundancies and gain a deeper understanding of the functional interactome governing mitophagy in AML. We identified OPTN (optineurin) as sole non-redundant mitophagy receptor and characterized its unique role in AML. Knockdown and overexpression experiments demonstrated that OPTN expression is rate-limiting for AML cell proliferation. In a MN1-driven murine transplantation model, loss of OPTN prolonged overall median survival by 7 days (+21%). Mechanistically, we found broadly impaired mitochondrial respiration and function with increased mitochondrial ROS, that most likely caused the proliferation defect. Our results decipher the intertwined network of mitophagy receptors in AML for both ubiquitin-dependent and receptor-mediated mitophagy, identify OPTN as a non-redundant tool to study mitophagy in the context of leukemia and suggest OPTN inhibition as an attractive therapeutic strategy. Abbreviations: AML: acute myeloid leukemia; CRISPR: Clustered Regularly Interspaced Short Palindromic Repeats; CTRL: control; DFP: deferiprone; GI: genetic interaction; KD: knockdown; KO: knockout; ldMBM, lineage-depleted murine bone marrow; LFC: log2 fold change; LIR: LC3-interacting region; LSC: leukemic stem cell; MAGeCK: Model-based Analysis of Genome-wide CRISPR-Cas9 Knockout; MDIVI-1: mitochondrial division inhibitor 1; MOI: multiplicity of infection; MOM: mitochondrial outer membrane; NAC: N-acetyl-L-cysteine; OA: oligomycin-antimycin A; OCR: oxygen consumption rate; OE: overexpression; OPTN: optineurin; PINK1: PTEN induced putative kinase 1; ROS: reactive oxygen species; SEM: standard error of the mean; TCGA: The Cancer Genome Atlas; TEM: transmission electron microscopy; UBD: ubiquitin-binding domain; WT: wild type [ABSTRACT FROM AUTHOR]

Dysregulation of physiological processes may contribute to Alzheimer’s disease (AD) development. We previously found that an increase in the level of physiological dysregulation (PD) in the aging body is associated with declining resilience and robustness to major diseases. Also, our genome-wide association study found that genes associated with the age-related increase in PD frequently represented pathways implicated in axon guidance and synaptic function, which in turn were linked to AD and related traits (e.g., amyloid, tau, neurodegeneration) in the literature. Here, we tested the hypothesis that genes involved in PD and axon guidance/synapse function may jointly influence onset of AD. We assessed the impact of interactions between SNPs in such genes on AD onset in the Long Life Family Study and sought to replicate the findings in the Health and Retirement Study. We found significant interactions between SNPs in the UNC5C and CNTN6, and PLXNA4 and EPHB2 genes that influenced AD onset in both datasets. Associations with individual SNPs were not statistically significant. Our findings, thus, support a major role of genetic interactions in the heterogeneity of AD and suggest the joint contribution of genes involved in PD and axon guidance/synapse function (essential for the maintenance of complex neural networks) to AD development. [ABSTRACT FROM AUTHOR]

Hypoxia is a critical pathological element in many human diseases, including ischemic stroke, myocardial infarction, and solid tumors. Of particular significance and interest of ours are the cellular and molecular mechanisms that underlie susceptibility or tolerance to low O2. Previous studies have demonstrated that Notch signaling pathway regulates hypoxia tolerance in both Drosophila melanogaster and humans. However, the mechanisms mediating Notch-conferred hypoxia tolerance are largely unknown. In this study, we delineate the evolutionarily conserved mechanisms underlying this hypoxia tolerant phenotype. We determined the role of a group of conserved genes that were obtained from a comparative genomic analysis of hypoxia-tolerant D.melanogaster populations and human highlanders living at the high-altitude regions of the world (Tibetans, Ethiopians, and Andeans). We developed a novel dual-UAS/Gal4 system that allows us to activate Notch signaling in the Eaat1-positive glial cells, which remarkably enhances hypoxia tolerance in D.melanogaster, and, simultaneously, knock down a candidate gene in the same set of glial cells. Using this system, we discovered that the interactions between Notch signaling and bnl (fibroblast growth factor), croc (forkhead transcription factor C), or Mkk4 (mitogen-activated protein kinase kinase 4) are important for hypoxia tolerance, at least in part, through regulating neuronal development and survival under hypoxic conditions. Becausethese genetic mechanisms are evolutionarily conserved, this group of genes may serve as novel targets for developing therapeutic strategies and have a strong potential to be translated to humans to treat/prevent hypoxia-related diseases.

Background: Understanding the impact of gene interactions on disease phenotypes is increasingly recognised as a crucial aspect of genetic disease research. This trend is reflected by the growing amount of clinical research on oligogenic diseases, where disease manifestations are influenced by combinations of variants on a few specific genes. Although statistical machine-learning methods have been developed to identify relevant genetic variant or gene combinations associated with oligogenic diseases, they rely on abstract features and black-box models, posing challenges to interpretability for medical experts and impeding their ability to comprehend and validate predictions. In this work, we present a novel, interpretable predictive approach based on a knowledge graph that not only provides accurate predictions of disease-causing gene interactions but also offers explanations for these results. Results: We introduce BOCK, a knowledge graph constructed to explore disease-causing genetic interactions, integrating curated information on oligogenic diseases from clinical cases with relevant biomedical networks and ontologies. Using this graph, we developed a novel predictive framework based on heterogenous paths connecting gene pairs. This method trains an interpretable decision set model that not only accurately predicts pathogenic gene interactions, but also unveils the patterns associated with these diseases. A unique aspect of our approach is its ability to offer, along with each positive prediction, explanations in the form of subgraphs, revealing the specific entities and relationships that led to each pathogenic prediction. Conclusion: Our method, built with interpretability in mind, leverages heterogenous path information in knowledge graphs to predict pathogenic gene interactions and generate meaningful explanations. This not only broadens our understanding of the molecular mechanisms underlying oligogenic diseases, but also presents a novel application of knowledge graphs in creating more transparent and insightful predictors for genetic research. [ABSTRACT FROM AUTHOR]

Predicting mutational effects is essential for the control of antibiotic resistance (ABR). Predictions are difficult when there are strong genotype-by-environment (G × E), gene-by-gene (G × G or epistatic) or gene-by-gene-by-environment (G × G × E) interactions. We quantified G × G × E effects in Escherichia coli across environmental gradients. We created intergenic fitness landscapes using gene knock-outs and single-nucleotide ABR mutations previously identified to vary in the extent of G × E effects in our environments of interest. Then, we measured competitive fitness across a complete combinatorial set of temperature and antibiotic dosage gradients. In this way, we assessed the predictability of 15 fitness landscapes across 12 different but related environments. We found G × G interactions and rugged fitness landscapes in the absence of antibiotic, but as antibiotic concentration increased, the fitness effects of ABR genotypes quickly overshadowed those of gene knock-outs, and the landscapes became smoother. Our work reiterates that some single mutants, like those conferring resistance or susceptibility to antibiotics, have consistent effects across genetic backgrounds in stressful environments. Thus, although epistasis may reduce the predictability of evolution in benign environments, evolution may be more predictable in adverse environments. This article is part of the theme issue 'Interdisciplinary approaches to predicting evolutionary biology'. [ABSTRACT FROM AUTHOR]

The objective of the present work was to study the development of the aerial part and the root system and the nutritional situation in young rubber tree plants, selecting clones with greater vigor and desired characteristics in a rootstock, relating them to the reduction of production time of seedlings, increasing the degree of compatibility and consequently uniformity in order to know their potential for use as rootstock for commercial clones and guarantee the producer an adequate and uniform productivity during the bleeding period of his rubber plantation. The experiment was carried out at APTA Regional at Pindorama, São Paulo State, Brazil. The experiment was conducted in the nursery in plastic bags using a randomized block design with seven treatments and four replications. Each plot was composed of 120 plants, totaling 3360 useful plants. The treatments consisted of plants from the six clones PB 235, RRIM 600, GT 1, PR 261, IAN 873, Tjir 1 x Tjir 16 and unselected seeds (SNS). At 10 months all clones are ready for grafting. Clone IAN 873 presents grafting conditions at 8 months, with greater precocity compared to the others. Regarding the diameter of the stem, it can be said that the best results were respectively from the clones IAN 873, GT 1 and PB 235. That in the endosperm there was a higher concentration of zinc and copper. As for the concentration of calcium in the seed, a greater presence was observed in the shell, while the other macronutrients and total protein occurred in greater amounts in the endosperm. Within foliar nutrients, there is no significant difference between micronutrients. [ABSTRACT FROM AUTHOR]

Simple Summary: Epigenetic factors control how genes are expressed in different cell types. Members of the Polycomb Repressive Complexes (PRCs) are critical for epigenetic control of gene transcription in blood cells, and mutations and deletions in these factors are common in the blood cancer acute myeloid leukaemia (AML). This review article provides an overview of how the structure and function of PRCs are affected by genetic alterations in AML, with a primary focus on PRC2 core factors. We document how mutations and deletions in PRC2 factors are linked to other AML-associated genetic alterations and discuss how these observations might inform potential treatment avenues in future. Epigenetic dysregulation is a hallmark of many haematological malignancies and is very frequent in acute myeloid leukaemia (AML). A cardinal example is the altered activity of the Polycomb Repressive Complex 2 (PRC2) due to somatic mutations and deletions in genes encoding PRC2 core factors that are necessary for correct complex assembly. These genetic alterations typically lead to reduced histone methyltransferase activity that, in turn, has been strongly linked to poor prognosis and chemoresistance. In this review, we provide an overview of genetic alterations of PRC components in AML, with particular reference to structural and functional features of PRC2 factors. We further review genetic interactions between these alterations and other AML-associated mutations in both adult and paediatric leukaemias. Finally, we discuss reported prognostic links between PRC2 mutations and deletions and disease outcomes and potential implications for therapy. [ABSTRACT FROM AUTHOR]

Alzheimer's disease (AD) is one of the most challenging neurodegenerative diseases because of its complicated and progressive mechanisms, and multiple risk factors. Increasing research evidence demonstrates that genetics may be a key factor responsible for the occurrence of the disease. Although previous reports identified quite a few AD-associated genes, they were mostly limited owing to patient sample size and selection bias. There is a lack of comprehensive research aimed to identify AD-associated risk mutations systematically. To address this challenge, we hereby construct a large-scale AD mutation and co-mutation framework ('AD-Syn-Net'), and propose deep learning models named Deep-SMCI and Deep-CMCI configured with fully connected layers that are capable of predicting cognitive impairment of subjects effectively based on genetic mutation and co-mutation profiles. Next, we apply the customized frameworks to data sets to evaluate the importance scores of the mutations and identified mutation effectors and co-mutation combination vulnerabilities contributing to cognitive impairment. Furthermore, we evaluate the influence of mutation pairs on the network architecture to dissect the genetic organization of AD and identify novel co-mutations that could be responsible for dementia, laying a solid foundation for proposing future targeted therapy for AD precision medicine. Our deep learning model codes are available open access here: https://github.com/Pan-Bio/AD-mutation-effectors. [ABSTRACT FROM AUTHOR]

Powerful new technologies for perturbing genetic elements have recently expanded the study of genetic interactions in model systems ranging from yeast to human cell lines. However, technical artifacts can confound signal across genetic screens and limit the immense potential of parallel screening approaches. To address this problem, we devised a novel PCA-based method for correcting genome-wide screening data, bolstering the sensitivity and specificity of detection for genetic interactions. Applying this strategy to a set of 436 whole genome CRISPR screens, we report more than 1.5 million pairs of correlated "co-functional" genes that provide finer-scale information about cell compartments, biological pathways, and protein complexes than traditional gene sets. Lastly, we employed a gene community detection approach to implicate core genes for cancer growth and compress signal from functionally related genes in the same community into a single score. This work establishes new algorithms for probing cancer cell networks and motivates the acquisition of further CRISPR screen data across diverse genotypes and cell types to further resolve complex cellular processes.

Phenotypic complexity is caused by the contributions of environmental factors and multiple genetic loci, interacting or acting independently. Studies of yeast and Arabidopsis often find that the majority of natural variation across phenotypes is attributable to independent additive quantitative trait loci (QTL). Detected loci in these organisms explain most of the estimated heritable variation. By contrast, many heritable components underlying phenotypic variation in metazoan models remain undetected. Before the relative impacts of additive and interactive variance components on metazoan phenotypic variation can be dissected, high replication and precise phenotypic measurements are required to obtain sufficient statistical power to detect loci contributing to this missing heritability. Here, we used a panel of 296 recombinant inbred advanced intercross lines of Caenorhabditis elegans and a high-throughput fitness assay to detect loci underlying responses to 16 different toxins, including heavy metals, chemotherapeutic drugs, pesticides, and neuropharmaceuticals. Using linkage mapping, we identified 82 QTL that underlie variation in responses to these toxins, and predicted the relative contributions of additive loci and genetic interactions across various growth parameters. Additionally, we identified three genomic regions that impact responses to multiple classes of toxins. These QTL hotspots could represent common factors impacting toxin responses. We went further to generate near-isogenic lines and chromosome substitution strains, and then experimentally validated these QTL hotspots, implicating additive and interactive loci that underlie toxin-response variation.

The aim of this study is to investigate whether the genetic, epigenetic and environmental factors that give rise to supernumeraries in the maxillary incisor region and larger dimensions of the adjacent maxillary incisors are also associated with variations in the morphology of the mandibular incisors. If so, this would contribute to understanding the distribution and interactions of factors during dental development and how these can be modelled. The sample consisted of 34 patients with supernumerary teeth in the maxillary anterior region, matched for gender, age and White Caucasian ethnicity with 34 control subjects. The average ages of the supernumerary and control groups were 12.8 and 12.2 years, respectively. Study models of all subjects were constructed and imaged using a previously validated system. Using custom software, each of the mandibular incisor teeth were measured to obtain 17 parameters from the labial view and 17 from the occlusal view. Principal component analysis (PCA) was used to summarize the measurements into a smaller set representing distinct features of the clinical crowns, followed by a comparison between the supernumerary and control groups using 2-way ANOVA. Seven factors of tooth size of the mandibular central incisors and six factors of the mandibular lateral incisors were identified as major features of the clinical crowns. All parameters of both mandibular incisors were greater in the supernumerary group than in the control, with three of these, located in the incisal and cervical regions of the mandibular lateral incisors, being statistically significantly larger. The findings of this study indicate that the aetiological factors associated with supernumerary teeth in the maxillary anterior region also affect tooth crown dimensions of mandibular incisors. This new evidence enhances several models of the interactions of genetic, epigenetic and environmental components of dental development and supports a multi-model approach to increase understanding of this process and its variations. [ABSTRACT FROM AUTHOR]

Objective: Effective adjuvant therapeutic strategies are urgently needed to overcome MAPK inhibitor (MAPKi) resistance, which is one of the most common forms of resistance that has emerged in many types of cancers. Here, we aimed to systematically identify the genetic interactions underlying MAPKi resistance, and to further investigate the mechanisms that produce the genetic interactions that generate synergistic MAPKi resistance. Methods: We conducted a comprehensive pair-wise sgRNA-based high-throughput screening assay to identify synergistic interactions that sensitized cancer cells to MAPKi, and validated 3 genetic combinations through competitive growth, cell viability, and spheroid formation assays. We next conducted Kaplan-Meier survival analysis based on The Cancer Genome Atlas database and conducted immunohistochemistry to determine the clinical relevance of these synergistic combinations. We also investigated the MAPKi resistance mechanisms of these validated synergistic combinations by using co-immunoprecipitation, Western blot, qRT-PCR, and immunofluorescence assays. Results: We constructed a systematic interaction network of MAPKi resistance and identified 3 novel synergistic combinations that effectively targeted MAPKi resistance (ITGB3 + IGF1R, ITGB3 + JNK, and HDGF + LGR5). We next analyzed their clinical relevance and the mechanisms by which they sensitized cancer cells to MAPKi exposure. Specifically, we discovered a novel protein complex, HDGF-LGR5, that adaptively responded to MAPKi to enhance cancer cell stemness, which was up- or downregulated by the inhibitors of ITGB3 + JNK or ITGB3 + IGF1R. Conclusions: Pair-wise sgRNA library screening provided systematic insights into elucidating MAPKi resistance in cancer cells. ITGB3-+ IGF1R-targeting drugs (cilengitide + linsitinib) could be used as an effective therapy for suppressing the adaptive formation of the HDGF-LGR5 protein complex, which enhanced cancer stemness during MAPKi stress.

Seminal yeast studies have established the value of comprehensively mapping genetic interactions (GIs) for inferring gene function. Efforts in human cells using focused gene sets underscore the utility of this approach, but the feasibility of generating large-scale, diverse human GI maps remains unresolved. We developed a CRISPR interference platform for large-scale quantitative mapping of human GIs. We systematically perturbed 222,784 gene pairs in two cancer cell lines. The resultant maps cluster functionally related genes, assigning function to poorly characterized genes, including TMEM261, a new electron transport chain component. Individual GIs pinpoint unexpected relationships between pathways, exemplified by a specific cholesterol biosynthesis intermediate whose accumulation induces deoxynucleotide depletion, causing replicative DNA damage and a synthetic-lethal interaction with the ATR/9-1-1 DNA repair pathway. Our map provides a broad resource, establishes GI maps as a high-resolution tool for dissecting gene function, and serves as a blueprint for mapping the genetic landscape of human cells.

Chemotherapy is used to treat most cancer patients, yet our understanding of factors that dictate response and resistance to such drugs remains limited. We report the generation of a quantitative chemical-genetic interaction map in human mammary epithelial cells charting the impact of the knockdown of 625 genes related to cancer and DNA repair on sensitivity to 29 drugs, covering all classes of chemotherapy. This quantitative map is predictive of interactions maintained in other cell lines, identifies DNA-repair factors, predicts cancer cell line responses to therapy, and prioritizes synergistic drug combinations. We identify that ARID1A loss confers resistance to PARP inhibitors in cells and ovarian cancer patients and that loss of GPBP1 causes resistance to cisplatin and PARP inhibitors through the regulation of genes involved in homologous recombination. This map helps navigate patient genomic data and optimize chemotherapeutic regimens by delineating factors involved in the response to specific types of DNA damage.

The metabolic pathways fueling tumor growth have been well characterized, but the specific impact of transforming events on network topology and enzyme essentiality remains poorly understood. To this end, we performed combinatorial CRISPR-Cas9 screens on a set of 51 carbohydrate metabolism genes that represent glycolysis and the pentose phosphate pathway (PPP). This high-throughput methodology enabled systems-level interrogation of metabolic gene dispensability, interactions, and compensation across multiple cell types. The metabolic impact of specific combinatorial knockouts was validated using 13C and 2H isotope tracing, and these assays together revealed key nodes controlling redox homeostasis along the KEAP-NRF2 signaling axis. Specifically, targeting KEAP1 in combination with oxidative PPP genes mitigated the deleterious effects of these knockouts on growth rates. These results demonstrate how our integrated framework, combining genetic, transcriptomic, and flux measurements, can improve elucidation of metabolic network alterations and guide precision targeting of metabolic vulnerabilities based on tumor genetics.

The natural ends of the linear eukaryotic chromosomes are protected by telomeres, which also play an important role in aging and cancer development. Telomere length varies between species, but it is strictly controlled in all organisms. The process of Telomere Length Maintenance (TLM) involves many pathways, protein complexes and interactions that were first discovered in budding and fission yeast model organisms (Saccharomyces cerevisiae, Schizosaccharomyces pombe). In particular, large-scale systematic genetic screens in budding yeast uncovered a network of -500 genes that, when mutated, cause telomeres to lengthen or to shorten. In contrast, the TLM network in fission yeast remains largely unknown and systematic data is still lacking. In this work we try to close this gap and develop a unified interpretable machine learning framework for TLM gene discovery and phenotype prediction in both species. We demonstrate the utility of our framework in pinpointing the pathways by which TLM homeostasis is maintained and predicting novel TLM genes in fission yeast. The results of this study could be used for better understanding of telomere biology and serve as a step towards the adaptation of computational methods based on telomeric data for human prognosis. [ABSTRACT FROM AUTHOR]

Background: Maternal genetic variations, including variations in mitochondrial biogenesis (MB) and oxidative phosphorylation (OP), are associated with placental abruption (PA). However, the role of maternal-fetal genetic interactions (MFGI) and parent-of-origin (imprinting) effects in PA remain unknown. Objective: To investigate MFGI in MB-OP, and imprinting effects in relation to risk of PA. Methods: Among Peruvian mother-infant pairs (503 PA cases and 1052 controls), independent single nucleotide polymorphisms (SNPs), with linkage-disequilibrium coefficient <0.80, were selected to characterize genetic variations in MB-OP (78 SNPs in 24 genes) and imprinted genes (2713 SNPs in 73 genes). For each MB-OP SNP, four multinomial models corresponding to fetal allele effect, maternal allele effect, maternal and fetal allele additive effect, and maternal-fetal allele interaction effect were fit under Hardy-Weinberg equilibrium, random mating, and rare disease assumptions. The Bayesian information criterion (BIC) was used for model selection. For each SNP in imprinted genes, imprinting effect was tested using a likelihood ratio test. Bonferroni corrections were used to determine statistical significance (p-value < 6.4e-4 for MFGI and p-value < 1.8e-5 for imprinting). Abruption cases were more likely to experience preeclampsia, have shorter gestational age, and deliver infants with lower birthweight compared with controls. Models with MFGI effects provided improved fit than models with only maternal and fetal genotype main effects for SNP rs12530904 (p-value = 1.2e-04) in calcium/calmodulin-dependent protein kinase [CaM kinase] II beta (CAMK2B), and, SNP rs73136795 (p-value = 1.9e-04) in peroxisome proliferator-activated receptor-gamma (PPARG), both MB genes. We identified 320 SNPs in 45 maternally-imprinted genes (including potassium voltage-gated channel subfamily Q member 1 [KCNQ1], neurotrimin [NTM], and, ATPase phospholipid transporting 10 A [ATP10A]) associated with abruption. Top hits included rs2012323 (p-value = 1.6E-16) and rs12221520 (p-value1.3e-13) in KCNQ1, rs8036892 (p-value = 9.3E-17) and rs188497582 in ATP10A, rs12589854 (p-value = 2.9E-11) and rs80203467 (p-value = 4.6e-11) in maternally expressed 8, small nucleolar RNA host (MEG8), and rs138281088 in solute carrier family 22 member 2 (SLC22A2) (p-value = 6.8e-9). We identified novel PA-related maternal-fetal MB gene interactions and imprinting effects that highlight the role of the fetus in PA risk development. Findings can inform mechanistic investigations to understand the pathogenesis of PA. [ABSTRACT FROM AUTHOR]

Asthma with severe exacerbation is one of the most common causes of hospitalization among young children. Exacerbations are typically triggered by respiratory infections, but the host factors causing recurrent infections and exacerbations in some children are poorly understood. As a result, current treatment options and preventive measures are inadequate. We sought to identify genetic interaction associated with the development of childhood asthma. We performed an exhaustive search for pairwise interaction between genetic single nucleotide polymorphisms using 1204 cases of a specific phenotype of early childhood asthma with severe exacerbations in patients aged 2 to 6 years combined with 5328 nonasthmatic controls. Replication was attempted in 3 independent populations, and potential underlying immune mechanisms were investigated in the COPSAC2010 and COPSAC2000 birth cohorts. We found evidence of interaction, including replication in independent populations, between the known childhood asthma loci CDHR3 and GSDMB. The effect of CDHR3 was dependent on the GSDMB genotype, and this interaction was more pronounced for severe and early onset of disease. Blood immune analyses suggested a mechanism related to increased IL-17A production after viral stimulation. We found evidence of interaction between CDHR3 and GSDMB in development of early childhood asthma, possibly related to increased IL-17A response to viral infections. This study demonstrates the importance of focusing on specific disease subtypes for understanding the genetic mechanisms of asthma. [Display omitted] [ABSTRACT FROM AUTHOR]

Although U2AF1 S34F is a recurrent splicing factor mutation in lung adenocarcinoma (ADC), U2AF1 S34F alone is insufficient for producing tumors in previous models. Because lung ADCs with U2AF1 S34F frequently have co-occurring KRAS mutations and smoking histories, we hypothesized that tumor-forming potential arises from U2AF1 S34F interacting with oncogenic KRAS and environmental stress. To elucidate the effect of U2AF1 S34F co-occurring with a second mutation, we generated human bronchial epithelial cells (HBEC3kt) with co-occurring U2AF1 S34F and KRAS G12V . Transcriptome analysis revealed that co-occurring U2AF1 S34F and KRAS G12V differentially impacts inflammatory, cell cycle, and KRAS pathways. Subsequent phenotyping found associated suppressed cytokine production, increased proliferation, anchorage-independent growth, and tumors in mouse xenografts. Interestingly, HBEC3kts harboring only U2AF1 S34F display increased splicing in stress granule protein genes and viability in cigarette smoke concentrate. Our results suggest that U2AF1 S34F may potentiate transformation by granting precancerous cells survival advantage in environmental stress, permitting accumulation of additional mutations like KRAS G12V , which synergize with U2AF1 S34F to transform the cell., Competing Interests: DECLARATION OF INTERESTS A.N.B. is a consultant for Remix Therapeutics, Inc.

The natural ends of the linear eukaryotic chromosomes are protected by telomeres, which also play an important role in aging and cancer development. Telomere length varies between species, but it is strictly controlled in all organisms. The process of Telomere Length Maintenance (TLM) involves many pathways, protein complexes and interactions that were first discovered in budding and fission yeast model organisms (Saccharomyces cerevisiae, Schizosaccharomyces pombe). In particular, large-scale systematic genetic screens in budding yeast uncovered a network of ≈500 genes that, when mutated, cause telomeres to lengthen or to shorten. In contrast, the TLM network in fission yeast remains largely unknown and systematic data is still lacking. In this work we try to close this gap and develop a unified interpretable machine learning framework for TLM gene discovery and phenotype prediction in both species. We demonstrate the utility of our framework in pinpointing the pathways by which TLM homeostasis is maintained and predicting novel TLM genes in fission yeast. The results of this study could be used for better understanding of telomere biology and serve as a step towards the adaptation of computational methods based on telomeric data for human prognosis.

Transfer RNA variants increase the frequency of mistranslation, the misincorporation of an amino acid not specified by the "standard" genetic code, to frequencies approaching 10% in yeast and bacteria. Cells cope with these variants by having multiple copies of each tRNA isodecoder and through pathways that deal with proteotoxic stress. In this study, we define the genetic interactions of the gene encoding tRNASer UGG, G26A, which mistranslates serine at proline codons. Using a collection of yeast temperature-sensitive alleles, we identify negative synthetic genetic interactions between the mistranslating tRNA and 109 alleles representing 91 genes, with nearly half of the genes having roles in RNA processing or protein folding and turnover. By regulating tRNA expression, we then compare the strength of the negative genetic interaction for a subset of identified alleles under differing amounts of mistranslation. The frequency of mistranslation correlated with the impact on cell growth for all strains analyzed; however, there were notable differences in the extent of the synthetic interaction at different frequencies of mistranslation depending on the genetic background. For many of the strains, the extent of the negative interaction with tRNASer UGG, G26A was proportional to the frequency of mistranslation or only observed at intermediate or high frequencies. For others, the synthetic interaction was approximately equivalent at all frequencies of mistranslation. As humans contain similar mistranslating tRNAs, these results are important when analyzing the impact of tRNA variants on disease, where both the individual's genetic background and the expression of the mistranslating tRNA variant need to be considered. [ABSTRACT FROM AUTHOR]

Abstract Background The impact of genetic interaction networks on evolution is a fundamental issue. Previous studies have demonstrated that the topology of the network is determined by the properties of the cellular machinery. Functionally related genes frequently interact with one another, and they establish modules, e.g., modules of protein complexes and biochemical pathways. In this study, we experimentally tested the hypothesis that compensatory evolutionary modifications, such as mutations and transcriptional changes, occur frequently in genes from perturbed modules of interacting genes. Results Using Saccharomyces cerevisiae haploid deletion mutants as a model, we investigated two modules lacking COG7 or NUP133, which are evolutionarily conserved genes with many interactions. We performed laboratory evolution experiments with these strains in two genetic backgrounds (with or without additional deletion of MSH2), subjecting them to continuous culture in a non-limiting minimal medium. Next, the evolved yeast populations were characterized through whole-genome sequencing and transcriptome analyses. No obvious compensatory changes resulting from inactivation of genes already included in modules were identified. The supposedly compensatory inactivation of genes in the evolved strains was only rarely observed to be in accordance with the established fitness effect of the genetic interaction network. In fact, a substantial majority of the gene inactivations were predicted to be neutral in the experimental conditions used to determine the interaction network. Similarly, transcriptome changes during continuous culture mostly signified adaptation to growth conditions rather than compensation of the absence of the COG7, NUP133 or MSH2 genes. However, we noticed that for genes whose inactivation was deleterious an upregulation of transcription was more common than downregulation. Conclusions Our findings demonstrate that the genetic interactions and the modular structure of the network described by others have very limited effects on the evolutionary trajectory following gene deletion of module elements in our experimental conditions and has no significant impact on short-term compensatory evolution. However, we observed likely compensatory evolution in functionally related (albeit non-interacting) genes.

The ability to sense and respond to osmotic fluctuations is critical for the maintenance of cellular integrity. We used gene co-essentiality analysis to identify an unappreciated relationship between TSC22D2 , WNK1 , and NRBP1 in regulating cell volume homeostasis. All of these genes have paralogs and are functionally buffered for osmo-sensing and cell volume control. Within seconds of hyperosmotic stress, TSC22D, WNK, and NRBP family members physically associate into biomolecular condensates, a process that is dependent on intrinsically disordered regions (IDRs). A close examination of these protein families across metazoans revealed that TSC22D genes evolved alongside a domain in NRBPs that specifically binds to TSC22D proteins, which we have termed NbrT (NRBP binding region with TSC22D), and this co-evolution is accompanied by rapid IDR length expansion in WNK-family kinases. Our study reveals that TSC22D , WNK , and NRBP genes evolved in metazoans to co-regulate rapid cell volume changes in response to osmolarity. [Display omitted] • TSC22D2 , WNK1 , and NRBP1 are functionally dominant gene family members • TSC22D, WNK, and NRBP protein families rapidly form reversible condensates • TSC22Ds mediate the translocation of NRBP1 into condensates • Association of TSC22Ds and NRBPs co-evolved with the expansion of IDRs in WNK kinases Xiao et al. identify components of metazoan cell volume regulation machinery through gene co-essentiality analysis from large-scale CRISPR datasets. TSC22D , WNK , and NRBP genes sense rapid cell volume changes to form condensates, and disordered TSC22D proteins mediate the translocation of NRBP pseudokinases, which co-evolved in metazoans to regulate cell volume recovery. [ABSTRACT FROM AUTHOR]

Genetic interactions (GIs) refer to two altered genes having a combined effect that is not seen individually. They play a crucial role in influencing drug efficacy. We utilized CGIdb 2.0 (http://www.medsysbio.org/CGIdb2/), an updated database of comprehensively published GIs information, encompassing synthetic lethality (SL), synthetic viability (SV), and chemical-genetic interactions. CGIdb 2.0 elucidates GIs relationships between or within protein complex models by integrating protein-protein physical interactions. Additionally, we introduced GENIUS (GENetic Interactions mediated drUg Signature) to leverage GIs for identifying the response signature of immune checkpoint inhibitors (ICIs). GENIUS identified high MAP4K4 expression as a resistant signature and high HERC4 expression as a sensitive signature for ICIs treatment. Melanoma patients with high expression of MAP4K4 were associated with decreased efficacy and poorer survival following ICIs treatment. Conversely, overexpression of HERC4 in melanoma patients correlated with a positive response to ICIs. Notably, HERC4 enhances sensitivity to immunotherapy by facilitating antigen presentation. Analyses of immune cell infiltration and single-cell data revealed that B cells expressing MAP4K4 may contribute to resistance to ICIs in melanoma. Overall, CGIdb 2.0, provides integrated GIs data, thus serving as a crucial tool for exploring drug effects. • CGIdb 2.0 comprehensively collects synthetic lethal, synthetic viable, and chemical-genetic interactions in cancer. • GENIUS predicts response signatures of immune checkpoint inhibitors by genetic interactions in cancer cells. • High expression of HERC4 mediate sensitivity and high expression of MAP4K4 induce resistance to immunotherapy in melanoma. [ABSTRACT FROM AUTHOR]

This paper presents the analytical and numerical comparison of two methods of estimation of additive × additive × additive (aaa) interaction of QTL effects. The first method takes into account only the plant phenotype, while in the second we also included genotypic information from molecular marker observation. Analysis was made on 150 doubled haploid (DH) lines of barley derived from cross Steptoe × Morex and 145 DH lines from Harrington × TR306 cross. In total, 153 sets of observation was analyzed. In most cases, aaa interactions were found with an exert effect on QTL. Results also show that with molecular marker observations, obtained estimators had smaller absolute values than phenotypic estimators. [ABSTRACT FROM AUTHOR]

We study local adaptation of a peripheral population by investigating the fate of new mutations using a haploid two-locus two-allele continent-island migration model. We explore how linkage, epistasis, and maladaptive gene flow affect the invasion probability of weakly beneficial de-novo mutations that arise on the island at an arbitrary physical distance to a locus that already maintains a stable migration-selection polymorphism. By assuming a slightly supercritical branching process, we deduce explicit conditions on the parameters that permit a positive invasion probability and we derive approximations for it. They show how the invasion probability depends on the additive and epistatic effects of the mutant, on its linkage to the polymorphism, and on the migration rate. We use these approximations together with empirically motivated distributions of epistatic effects to analyze the influence of epistasis on the expected invasion probability if mutants are drawn randomly from such a distribution and occur at a random physical distance to the existing polymorphism. We find that the invasion probability generally increases as the epistasis parameter increases or the migration rate decreases, but not necessarily as the recombination rate decreases. Finally, we shed light on the size of emerging genomic islands of divergence by exploring the size of the chromosomal neighborhood of the already established polymorphism in which 50% or 90% of the successfully invading mutations become established. These 'window sizes' always decrease in a reverse sigmoidal way with stronger migration and typically increase with increasing epistatic effect. [ABSTRACT FROM AUTHOR]

Mice expressing human amyloid precursor protein (APP) containing the dominant Swedish and Iberian mutations (AppNL–F) or also Arctic mutation (AppNL–G–F) show neuropathology and hippocampus-dependent cognitive impairments pertinent to Alzheimer's disease (AD) in mouse models at 18 and 6 months of age, respectively. Apolipoprotein E, involved in cholesterol metabolism, plays an important role in maintaining the brain. There are three human apolipoprotein E isoforms: E2, E3, and E4. Compared to E3, E4 increases while E2 protects against AD risk. At 6 months of age, prior to the onset of plaque pathology, E3, but not E4, protected against hAPP/Aβ-induced impairments in spatial memory retention in the Morris water maze. However, these earlier studies were limited as hapoE was not expressed outside the brain and E3 or E4 was not expressed under control of an apoE promotor, E2 was often not included, hAPP was transgenically overexpressed and both mouse and hAPP were present. Therefore, to determine whether apoE has isoform-dependent effects on hAPP/Aβ-induced behavioral alterations and cognitive impairments in adult female and male mice at 6 and 18 months of age, we crossed AppNL–G–F and AppNL–F mice with E2, E3, and E4 mice. To distinguish whether genotype differences seen at either time point were due to main effects of hAPP, hapoE, or hAPP × hapoE genetic interactions, we also behavioral and cognitively tested E2, E3, and E4 female and male mice at 6 and 18 months of age. We also compared behavioral and cognitive performance of 18-month-old AppNL–G–F and AppNL–F female and male mice on a murine apoE background along with that of age—and sex-matched C57BL/6J wild-type mice. For many behavioral measures at both time points there were APP × APOE interactions, supporting that apoE has isoform-dependent effects on hAPP/Aβ-induced behavioral and cognitive performance. NL-G-F/E3, but not NL-G-F/E2, mice had lower cortical insoluble Aβ42 levels than NL-G-F/E4 mice. NL-F/E3 and NL-F/E2 mice had lower cortical insoluble Aβ42 levels than NL-F/E4 mice. These results demonstrate that there are apoE isoform-dependent effects on hAPP/Aβ-induced behavioral alterations and cognitive impairments and cortical insoluble Aβ42 levels in mouse models containing only human APP and apoE. [ABSTRACT FROM AUTHOR]

Poly (ADPribose) polymerase inhibitors (PARPis) are clinically approved drugs designed according to the concept of synthetic lethality (SL) interaction. It is crucial to expand the scale of patients who can benefit from PARPis, and overcome drug resistance associated with it. Genetic interactions (GIs) include SL and synthetic viability (SV) that participate in drug response in cancer cells. Based on the hypothesis that mutated genes with SL or SV interactions with PARP1/2/3 are potential sensitive or resistant PARPis biomarkers, respectively, we developed a novel computational method to identify them. We analyzed fitness variation of cell lines to identify PARP1/2/3-related GIs according to CRISPR/Cas9 and RNA interference functional screens. Potential resistant/sensitive mutated genes were identified using pharmacogenomic datasets. We identified 41 candidate resistant and 130 candidate sensitive PARPi-response related genes, and observed that EGFR with gain-of-function mutation induced PARPi resistance, and predicted a combination therapy with PARP inhibitor (veliparib) and EGFR inhibitor (erlotinib) for lung cancer. We also revealed that a resistant gene set (TNN, PLEC, and TRIP12) in lower grade glioma and a sensitive gene set (BRCA2, TOP3A, and ASCC3) in ovarian cancer, which were associated with prognosis. Thus, cancer genome-derived GIs provide new insights for identifying PARPi biomarkers and a new avenue for precision therapeutics.

Mice expressing human amyloid precursor protein (APP) containing the dominant Swedish and Iberian mutations (AppNL–F) or also Arctic mutation (AppNL–G–F) show neuropathology and hippocampus-dependent cognitive impairments pertinent to Alzheimer’s disease (AD) in mouse models at 18 and 6 months of age, respectively. Apolipoprotein E, involved in cholesterol metabolism, plays an important role in maintaining the brain. There are three human apolipoprotein E isoforms: E2, E3, and E4. Compared to E3, E4 increases while E2 protects against AD risk. At 6 months of age, prior to the onset of plaque pathology, E3, but not E4, protected against hAPP/Aβ-induced impairments in spatial memory retention in the Morris water maze. However, these earlier studies were limited as hapoE was not expressed outside the brain and E3 or E4 was not expressed under control of an apoE promotor, E2 was often not included, hAPP was transgenically overexpressed and both mouse and hAPP were present. Therefore, to determine whether apoE has isoform-dependent effects on hAPP/Aβ-induced behavioral alterations and cognitive impairments in adult female and male mice at 6 and 18 months of age, we crossed AppNL–G–F and AppNL–F mice with E2, E3, and E4 mice. To distinguish whether genotype differences seen at either time point were due to main effects of hAPP, hapoE, or hAPP × hapoE genetic interactions, we also behavioral and cognitively tested E2, E3, and E4 female and male mice at 6 and 18 months of age. We also compared behavioral and cognitive performance of 18-month-old AppNL–G–F and AppNL–F female and male mice on a murine apoE background along with that of age—and sex-matched C57BL/6J wild-type mice. For many behavioral measures at both time points there were APP × APOE interactions, supporting that apoE has isoform-dependent effects on hAPP/Aβ-induced behavioral and cognitive performance. NL-G-F/E3, but not NL-G-F/E2, mice had lower cortical insoluble Aβ42 levels than NL-G-F/E4 mice. NL-F/E3 and NL-F/E2 mice had lower cortical insoluble Aβ42 levels than NL-F/E4 mice. These results demonstrate that there are apoE isoform-dependent effects on hAPP/Aβ-induced behavioral alterations and cognitive impairments and cortical insoluble Aβ42 levels in mouse models containing only human APP and apoE.

Genome-scale expression studies and comprehensive loss-of-function genetic screens have focused almost exclusively on the highest confidence candidate genes. Here, we describe a strategy for characterizing the lower confidence candidates identified by such approaches. We interrogated 177 genes that we classified as essential for the proliferation of cancer cells exhibiting constitutive β-catenin activity and integrated data for each of the candidates, derived from orthogonal short hairpin RNA (shRNA) knockdown and clustered regularly interspaced short palindromic repeats (CRISPR)-Cas9-mediated gene editing knockout screens, to yield 69 validated genes. We then characterized the relationships between sets of these genes using complementary assays: medium-throughput stable isotope labeling by amino acids in cell culture (SILAC)-based mass spectrometry, yielding 3,639 protein-protein interactions, and a CRISPR-mediated pairwise double knockout screen, yielding 375 combinations exhibiting greater- or lesser-than-additive phenotypic effects indicating genetic interactions. These studies identify previously unreported regulators of β-catenin, define functional networks required for the survival of β-catenin-active cancers, and provide an experimental strategy that may be applied to define other signaling networks.

Abstract Background Genome-wide association studies (GWAS) have proven successful in predicting genetic risk of disease using single-locus models; however, identifying single nucleotide polymorphism (SNP) interactions at the genome-wide scale is limited due to computational and statistical challenges. We addressed the computational burden encountered when detecting SNP interactions for survival analysis, such as age of disease-onset. To confront this problem, we developed a novel algorithm, called the Efficient Survival Multifactor Dimensionality Reduction (ES-MDR) method, which used Martingale Residuals as the outcome parameter to estimate survival outcomes, and implemented the Quantitative Multifactor Dimensionality Reduction method to identify significant interactions associated with age of disease-onset. Methods To demonstrate efficacy, we evaluated this method on two simulation data sets to estimate the type I error rate and power. Simulations showed that ES-MDR identified interactions using less computational workload and allowed for adjustment of covariates. We applied ES-MDR on the OncoArray-TRICL Consortium data with 14,935 cases and 12,787 controls for lung cancer (SNPs = 108,254) to search over all two-way interactions to identify genetic interactions associated with lung cancer age-of-onset. We tested the best model in an independent data set from the OncoArray-TRICL data. Results Our experiment on the OncoArray-TRICL data identified many one-way and two-way models with a single-base deletion in the noncoding region of BRCA1 (HR 1.24, P = 3.15 × 10–15), as the top marker to predict age of lung cancer onset. Conclusions From the results of our extensive simulations and analysis of a large GWAS study, we demonstrated that our method is an efficient algorithm that identified genetic interactions to include in our models to predict survival outcomes.

Abstract Essential genes tend to be highly conserved across eukaryotes, but, in some cases, their critical roles can be bypassed through genetic rewiring. From a systematic analysis of 728 different essential yeast genes, we discovered that 124 (17%) were dispensable essential genes. Through whole‐genome sequencing and detailed genetic analysis, we investigated the genetic interactions and genome alterations underlying bypass suppression. Dispensable essential genes often had paralogs, were enriched for genes encoding membrane‐associated proteins, and were depleted for members of protein complexes. Functionally related genes frequently drove the bypass suppression interactions. These gene properties were predictive of essential gene dispensability and of specific suppressors among hundreds of genes on aneuploid chromosomes. Our findings identify yeast's core essential gene set and reveal that the properties of dispensable essential genes are conserved from yeast to human cells, correlating with human genes that display cell line‐specific essentiality in the Cancer Dependency Map (DepMap) project.

Prediction of phenotypes from genotypes is an important objective to fulfill the promises of genomics, precision medicine and agriculture. Although it’s now possible to account for the majority of genetic variation through model fitting, prediction of phenotypes remains a challenge, especially across populations that have diverged in the past. In this study, we designed simulation experiments to specifically investigate the role of genetic interactions in failure of polygenic prediction. We found that non-additive genetic interactions can significantly reduce the accuracy of polygenic prediction. Our study demonstrated the importance of considering genetic interactions in genetic prediction.

Genome-wide association studies and experimental mouse models implicate the MIB1 and GATA6 genes in congenital heart disease (CHD). Their close physical proximity and conserved synteny suggest that these two genes might be involved in analogous cardiac developmental processes. Heterozygous Gata6 loss-of-function mutations alone or humanized Mib1 mutations in a NOTCH1-sensitized genetic background cause bicuspid aortic valve (BAV) and a membranous ventricular septal defect (VSD), consistent with MIB1 and NOTCH1 functioning in the same pathway. To determine if MIB1-NOTCH and GATA6 interact in valvular and septal development, we generated compound heterozygote mice carrying different Mib1 missense ( Mib1 K735R and Mib1 V943F ) or nonsense ( Mib1 R530X ) mutations with the Gata6 STOP/+ heterozygous null mutation. Combining Mib1 R530X/+ or Mib1 K735R/+ with Gata6 STOP/+ does not affect Gata6 STOP/+ single mutant phenotypes. In contrast, combining Mib1 V943F/+ with Gata6 STOP/+ decreases the incidence of BAV and VSD by 50%, suggesting a suppressive effect of Mib1 V943F/+ on Gata6 STOP/+ . Transcriptomic and functional analyses revealed that while the EMT pathway term is depleted in the Gata6 STOP/+ mutant, introducing the Mib1 V943F variant robustly enriches this term, consistent with the Mib1 V943F/+ phenotypic suppression of Gata6 STOP/+ . Interestingly, combined Notch1 and Gata6 insufficiency led to a nearly fully penetrant VSD but did not affect the BAV phenotype, underscoring the complex functional relationship between MIB1, NOTCH, and GATA6 in valvular and septal development.

Formation of the corpus callosum (CC), anterior commissure (AC), and postoptic commissure (POC), connecting the left and right cerebral hemispheres, is crucial for cerebral functioning. Collapsin response mediator protein 2 (CRMP2) has been suggested to be associated with the mechanisms governing this formation, based on knockout studies in mice and knockdown/knockout studies in zebrafish. Previously, we reported two cases of non-synonymous CRMP2 variants with S14R and R565C substitutions. Among the, the R565C substitution (p.R565C) was caused by the novel CRMP2 mutation c.1693C > T, and the patient presented with intellectual disability accompanied by CC hypoplasia. In this study, we demonstrate that crmp2 mRNA could rescue AC and POC formation in crmp2-knockdown zebrafish, whereas the mRNA with the R566C mutation could not. Zebrafish CRMP2 R566C corresponds to human CRMP2 R565C. Further experiments with transfected cultured cells indicated that CRMP2 with the R566C mutation could not bind to kinesin light chain 1 (KLC1). Knockdown of klc1a in zebrafish resulted in defective AC and POC formation, revealing a genetic interaction with crmp2. These findings suggest that the CRMP2 R566C mutant fails to bind to KLC1, preventing axonal elongation and leading to defective AC and POC formation in zebrafish and CC formation defects in humans. Our study highlights the importance of the interaction between CRMP2 and KLC1 in the formation of the forebrain commissures, revealing a novel mechanism associated with CRMP2 mutations underlying human neurodevelopmental abnormalities., (© 2024 Wiley Periodicals LLC.)

Dissecting the genetic mechanisms underlying agronomic traits is of great importance for crop breeding. Agronomic traits are usually controlled by multiple quantitative trait loci (QTLs) and genetic interactions, and mapping the underlying causal genes is still labor-intensive and time-consuming. Here, we present a genetic tool for directly targeting the specific causal genes by using a single-gene resolution linkage map that was constructed from 3756 F2 rice plants via targeted sequencing technology and Tukey-Kramer multiple comparisons test. Three large- and moderate-effect QTLs, qHD6-2, qGL3-1, and qGW5-2, were successfully mapped to their specific causal genes, Hd1, GS3, and GW5, respectively. A complex genetic interaction network containing 30 QTL–QTL interactions was constructed, revealing that the alternative allele of hub QTL, qHD6-2, can hide or release the genetic contributions of the alleles at interacting loci. Moreover, arranging genetic interactions in the models lead to more accurate phenotypic predictions. These results provide a community resource and new feasible strategy for deciphering the genetic mechanisms of complex agronomic traits and accelerating crop breeding. [ABSTRACT FROM AUTHOR]

A major goal of aging research is identifying genetic targets that could be used to slow or reverse aging – changes in the body and extend limits of human lifespan. However, majority of genes that showed the anti-aging and pro-survival effects in animal models were not replicated in humans, with few exceptions. Potential reasons for this lack of translation include a highly conditional character of genetic influence on lifespan, and its heterogeneity, meaning that better survival may be result of not only activity of individual genes, but also gene–environment and gene–gene interactions, among other factors. In this paper, we explored associations of genetic interactions with human lifespan. We selected candidate genes from well-known aging pathways (IGF1/FOXO growth signaling, P53/P16 apoptosis/senescence, and mTOR/SK6 autophagy and survival) that jointly decide on outcomes of cell responses to stress and damage, and so could be prone to interactions. We estimated associations of pairwise statistical epistasis between SNPs in these genes with survival to age 85+ in the Atherosclerosis Risk in Communities study, and found significant (FDR < 0.05) effects of interactions between SNPs in IGF1R, TGFBR2, and BCL2 on survival 85+. We validated these findings in the Cardiovascular Health Study sample, with P < 0.05, using survival to age 85+, and to the 90th percentile, as outcomes. Our results show that interactions between SNPs in genes from the aging pathways influence survival more significantly than individual SNPs in the same genes, which may contribute to heterogeneity of lifespan, and to lack of animal to human translation in aging research.

Understanding how genes interact is a central challenge in biology. Experimental evolution provides a useful, but underutilized, tool for identifying genetic interactions, particularly those that involve non-loss-of-function mutations or mutations in essential genes. We previously identified a strong positive genetic interaction between specific mutations in KEL1 (P344T) and HSL7 (A695fs) that arose in an experimentally evolved Saccharomyces cerevisiae population. Because this genetic interaction is not phenocopied by gene deletion, it was previously unknown. Using "evolutionary replay" experiments, we identified additional mutations that have positive genetic interactions with the kel1 -P344T mutation. We replayed the evolution of this population 672 times from six timepoints. We identified 30 populations where the kel1 -P344T mutation reached high frequency. We performed whole-genome sequencing on these populations to identify genes in which mutations arose specifically in the kel1 -P344T background. We reconstructed mutations in the ancestral and kel1 -P344T backgrounds to validate positive genetic interactions. We identify several genetic interactors with KEL1 , we validate these interactions by reconstruction experiments, and we show these interactions are not recapitulated by loss-of-function mutations. Our results demonstrate the power of experimental evolution to identify genetic interactions that are positive, allele specific, and not readily detected by other methods, shedding light on an underexplored region of the yeast genetic interaction network. [ABSTRACT FROM AUTHOR]

We have achieved a residue-level resolution of genetic interaction mapping - a technique that measures how the function of one gene is affected by the alteration of a second gene - by analyzing point mutations. Here, we describe how to interpret point mutant genetic interactions, and outline key applications for the approach, including interrogation of protein interaction interfaces and active sites, and examination of post-translational modifications. Genetic interaction analysis has proven effective for characterizing cellular processes; however, to date, systematic high-throughput genetic interaction screens have relied on gene deletions or knockdowns, which limits the resolution of gene function analysis and poses problems for multifunctional genes. Our point mutant approach addresses these issues, and further provides a tool for in vivo structure-function analysis that complements traditional biophysical methods. We also discuss the potential for genetic interaction mapping of point mutations in human cells and its application to personalized medicine.