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.

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.

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.

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.

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.

Primary open-angle glaucoma (POAG), a leading cause of irreversible blindness globally, shows disparity in prevalence and manifestations across ancestries. We perform meta-analysis across 15 biobanks (of the Global Biobank Meta-analysis Initiative) (n = 1,487,441: cases = 26,848) and merge with previous multi-ancestry studies, with the combined dataset representing the largest and most diverse POAG study to date (n = 1,478,037: cases = 46,325) and identify 17 novel significant loci, 5 of which were ancestry specific. Gene-enrichment and transcriptome-wide association analyses implicate vascular and cancer genes, a fifth of which are primary ciliary related. We perform an extensive statistical analysis of SIX6 and CDKN2B-AS1 loci in human GTEx data and across large electronic health records showing interaction between SIX6 gene and causal variants in the chr9p21.3 locus, with expression effect on CDKN2A/B. Our results suggest that some POAG risk variants may be ancestry specific, sex specific, or both, and support the contribution of genes involved in programmed cell death in POAG pathogenesis.

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.

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.

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.

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.
CONCLUSIONS: 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.