Epistasis, 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.
We 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.
This 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.
Our 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.
•Smoking-related phenotypes represent the most common comorbidities and risk indicators of OUD.•Genes associated with OUD phenotypes are also heavily associated with smoking-related phenotypes.•The ...rates of association between OUD and psychiatric disorders are similar across genes and observationally in the literature.•This review highlights the importance of multi-disorder comorbidity in both identifying and treating those prone to OUD.
Opioid use disorder (OUD) creates significant public health and economic burdens worldwide. Therefore, understanding the risk factors that lead to the development of OUD is fundamental to reducing both its prevalence and its impact. Significant sources of OUD risk include co-occurring lifetime and current diagnoses of both psychiatric disorders, primarily mood disorders, and other substance use disorders, and unique and shared genetic factors. Although there appears to be pleiotropy between OUD and both mood and substance use disorders, this aspect of OUD risk is poorly understood. In this review, we describe the prevalence and clinical significance of addictive and affective comorbidities as risk factors for OUD development as a basis for rational opioid prescribing and OUD treatment and to improve efforts to prevent the disorder. We also review the genetic variants that have been associated with OUD and other addictive and affective disorders to highlight targets for future study and risk assessment protocols.
Neutrophil migration is critical to the initiation and resolution of inflammation. Macrophage-1 antigen (Mac-1; CD11b/CD18, αMβ2) is a leukocyte integrin essential for firm adhesion to endothelial ...ICAM-1 (intercellular adhesion molecule 1) and migration of neutrophils in the shear forces of the circulation. PDI (protein disulfide isomerase) has been reported to influence neutrophil adhesion and migration. We aimed to elucidate the molecular mechanism of PDI control of Mac-1 affinity for ICAM-1 during neutrophil migration under fluid shear.
Neutrophils isolated from whole blood were perfused over microfluidic chips coated with ICAM-1. Colocalization of Mac-1 and PDI on neutrophils was visualized by fluorescently labeled antibodies and confocal microscopy. The redox state of Mac-1 disulfide bonds was mapped by differential cysteine alkylation and mass spectrometry. Wild-type or disulfide mutant Mac-1 was expressed recombinantly in Baby Hamster Kidney cells to measure ligand affinity. Mac-1 conformations were measured by conformation-specific antibodies and molecular dynamics simulations. Neutrophils crawling on immobilized ICAM-1 were measured in presence of oxidized or reduced PDI, and the effect of PDI inhibition using isoquercetin on neutrophil crawling on inflamed endothelial cells was examined. Migration indices in the X- and Y-direction were determined and the crawling speed was calculated.
PDI colocalized with high-affinity Mac-1 at the trailing edge of stimulated neutrophils when crawling on ICAM-1 under fluid shear. PDI cleaved 2 allosteric disulfide bonds, C169-C176 and C224-C264, in the βI domain of the β2 subunit, and cleavage of the C224-C264 disulfide bond selectively controls Mac-1 disengagement from ICAM-1 under fluid shear. Molecular dynamics simulations and conformation-specific antibodies reveal that cleavage of the C224-C264 bond induces conformational change and mechanical stress in the βI domain. This allosterically alters the exposure of an αI domain epitope associated with a shift of Mac-1 to a lower-affinity state. These molecular events promote neutrophil motility in the direction of flow at high shear stress. Inhibition of PDI by isoquercetin reduces neutrophil migration in the direction of flow on endothelial cells during inflammation.
Shear-dependent PDI cleavage of the neutrophil Mac-1 C224-C264 disulfide bond triggers Mac-1 de-adherence from ICAM-1 at the trailing edge of the cell and enables directional movement of neutrophils during inflammation.
Environments often vary across a life cycle, imposing fluctuating natural selection across development. Such fluctuating selection can favor different phenotypes in different life stages, but ...stage-specific evolutionary responses will depend on genetic variance, covariance, and their interaction across development and across environments. Thus, quantifying how genetic architecture varies with plastic responses to the environment and across development is vital to predict whether stage-specific adaptation will occur in nature. Additionally, the interaction of genetic variation and environmental plasticity (GxE) may be stage-specific, leading to a three-way interaction between genotype, environment, and development or GxDxE. To test for these patterns, we exposed larvae and adults of Drosophila melanogaster isogenic lines derived from a natural population to extreme heat and cold stress after developmental acclimation to cool (18 °C) and warm (25 °C) conditions and measured genetic variance for thermal hardiness. We detected significant GxE that was specific to larvae and adults for cold and heat hardiness (GxDxE), but no significant genetic correlation across development for either trait at either acclimation temperature. However, cross-development phenotypic correlations for acclimation responses suggest that plasticity itself may be developmentally constrained, though rigorously testing this hypothesis requires more experimentation. These results illustrate the potential for stage-specific adaptation within a complex life cycle and demonstrate the importance of measuring traits at appropriate developmental stages and environmental conditions when predicting evolutionary responses to changing climates.
In follicular lymphoma (FL), surface immunoglobulin (sIg) carries mandatory N-glycosylation sites in the variable regions, inserted during somatic hypermutation. These glycosylation sites are ...tumor-specific, indicating a critical function in FL. Added glycan unexpectedly terminates at high mannose (Mann) and confers capability for sIg-mediated interaction with local macrophage-expressed DC-SIGN lectin resulting in low-level activation of upstream B-cell receptor signaling responses. Here we show that despite being of low-level, DC-SIGN induces a similar downstream transcriptional response to anti-IgM in primary FL cells, characterized by activation of pathways associated with B-cell survival, proliferation and cell-cell communication. Lectin binding was also able to engage post-transcriptional receptor cross-talk pathways since, like anti-IgM, DC-SIGN down-modulated cell surface expression of CXCR4. Importantly, pre-exposure of a FL-derived cell line expressing sIgM-Mann or primary FL cells to DC-SIGN, which does not block anti-IgM binding, reversibly paralyzed the subsequent Ca
response to anti-IgM. These novel findings indicate that modulation of sIg function occurs in FL via lectin binding to acquired mannoses. The B-cell receptor alternative engagement described here provides two advantages to lymphoma cells: (i) activation of signaling, which, albeit of low-level, is sufficient to trigger canonical lymphoma-promoting responses, and (ii) protection from exogenous antigen by paralyzing anti-IgM-induced signaling. Blockade of this alternative engagement could offer a new therapeutic strategy.
How proteins harness mechanical force to control function is a significant biological question. Here we describe a human cell surface receptor that couples ligand binding and force to trigger a ...chemical event which controls the adhesive properties of the receptor. Our studies of the secreted platelet oxidoreductase, ERp5, have revealed that it mediates release of fibrinogen from activated platelet αIIbβ3 integrin. Protein chemical studies show that ligand binding to extended αIIbβ3 integrin renders the βI-domain Cys177-Cys184 disulfide bond cleavable by ERp5. Fluid shear and force spectroscopy assays indicate that disulfide cleavage is enhanced by mechanical force. Cell adhesion assays and molecular dynamics simulations demonstrate that cleavage of the disulfide induces long-range allosteric effects within the βI-domain, mainly affecting the metal-binding sites, that results in release of fibrinogen. This coupling of ligand binding, force and redox events to control cell adhesion may be employed to regulate other protein-protein interactions.
The opioid epidemic continues to contribute to loss of life through overdose and significant social and economic burdens. Many individuals who develop problematic opioid use (POU) do so after being ...exposed to prescribed opioid analgesics. Therefore, it is important to accurately identify and classify risk factors for POU. In this review, we discuss the etiology of POU and highlight novel approaches to identifying its risk factors. These approaches include the application of polygenic risk scores (PRS) and diverse machine learning (ML) algorithms used in tandem with data from electronic health records (EHR), clinical notes, patient demographics, and digital footprints. The implementation and synergy of these types of data and approaches can greatly assist in reducing the incidence of POU and opioid-related mortality by increasing the knowledge base of patient-related risk factors, which can help to improve prescribing practices for opioid analgesics.
As organisms age the environment fluctuates, exerting differential selection across ontogeny. In particular, highly seasonal environments expose life stages to often drastically different thermal ...environments. This developmental variation is particularly striking in organisms with complex life cycles, wherein life history stages also exhibit distinct morphologies, physiologies, and behaviors. Genes acting pleiotropically on thermal responses may produce genetic correlations across ontogeny, constraining the independent evolution of each life stage to their respective thermal environments. To investigate whether developmental genetic correlations constrain the evolution thermal hardiness of the fly Drosophila melanogaster, we applied quantitative genetic analyses to cold hardiness measured in both larvae and adults from isogenic lines of the Drosophila Genetic Reference Panel (DGRP), using survival at stressful low temperatures as the phenotypic metric. Using full genome resequencing data for the DGRP, we also implemented genome-wide association (GWA) analysis using Bayesian Sparse Linear Mixed Models (BSLMMs) to estimate associations between naturally segregating variation and cold hardiness for both larvae and adults. Quantitative genetic analyses revealed no significant genetic correlation for cold hardiness between life stages, suggesting complete genetic decoupling of thermal hardiness across the metamorphic boundary. Both quantitative genetic and GWA analyses suggested that polygenic variation underlies cold hardiness in both stages, and that associated loci largely affected one stage or the other, but not both. However, reciprocal enrichment tests and correlations between BSLMM parameters for each life stage support some shared physiological mechanisms that may reflect common cellular thermal response pathways. Overall, these results suggest no developmental genetic constraints on cold hardiness across metamorphosis in D. melanogaster, an important consideration in evolutionary models of responses to changing climates. Genetic correlations for environmental sensitivity across ontogeny remains largely unexplored in other organisms, thus assessing the generality of genetic decoupling will require further quantitative or population genetic analysis in additional species.
The αIIbβ3 integrin receptor coordinates platelet adhesion, activation, and mechanosensing in thrombosis and hemostasis. Using differential cysteine alkylation and mass spectrometry, we have ...identified a disulfide bond in the αIIb subunit linking cysteines 490 and 545 that is missing in ∼1 in 3 integrin molecules on the resting and activated human platelet surface. This alternate covalent form of αIIbβ3 is predetermined as it is also produced by human megakaryoblasts and baby hamster kidney fibroblasts transfected with recombinant integrin. From coimmunoprecipitation experiments, the alternate form selectively partitions into focal adhesions on the activated platelet surface. Its function was evaluated in baby hamster kidney fibroblast cells expressing a mutant integrin with an ablated C490-C545 disulfide bond. The disulfide mutant integrin has functional outside-in signaling but extended residency time in focal adhesions due to a reduced rate of clathrin-mediated integrin internalization and recycling, which is associated with enhanced affinity of the αIIb subunit for clathrin adaptor protein 2. Molecular dynamics simulations indicate that the alternate covalent form of αIIb requires higher forces to transition from bent to open conformational states that is in accordance with reduced affinity for fibrinogen and activation by manganese ions. These findings indicate that the αIIbβ3 integrin receptor is produced in various covalent forms that have different cell surface distribution and function. The C490, C545 cysteine pair is conserved across all 18 integrin α subunits, and the disulfide bond in the αV and α2 subunits in cultured cells is similarly missing, suggesting that the alternate integrin form and function are also conserved.
•The platelet αIIbβ3 integrin receptor is produced in various covalent forms that have different functions.•One-third of αIIbβ3 molecules are missing an αIIb disulfide bond that changes distribution and function of the receptor.
Display omitted
β2-Glycoprotein I (β2GPI) is an evolutionary conserved, abundant circulating protein. Although its function remains uncertain, accumulated evidence points toward interactions with endothelial cells ...and components of the coagulation system, suggesting a regulatory role in vascular biology. Our group has shown that thioredoxin 1 (TRX-1) generates free thiols in β2GPI, a process that may have a regulatory role in platelet adhesion. This report extends these studies and shows for the first time evidence of β2GPI with free thiols in vivo in both multiple human and murine serum samples. To explore how the vascular surface may modulate the redox status of β2GPI, unstimulated human endothelial cells and EAhy926 cells are shown to be capable of amplifying the effect of free thiol generation within β2GPI. Multiple oxidoreductase enzymes, such as endoplasmic reticulum protein 46 (ERp 46) and TRX-1 reductase, in addition to protein disulfide isomerase are secreted on the surface of endothelial cells. Furthermore, one or more of these generated free thiols within β2GPI are also shown to be nitrosylated. Finally, the functional significance of these findings is explored, by showing that free thiol–containing β2GPI has a powerful effect in protecting endothelial cells and EAhy926 cells from oxidative stress–induced cell death.