Convergent Balancing Selection on the Mu-Opioid Receptor in Primates

Variation in the mu-opioid receptor gene (OPRM1) moderates the influence of maternal sensitivity on child attachment

The endogenous opioid system is thought to play an important role in mother-infant attachment. In infant rhesus macaques, variation in the μ-opioid receptor gene (OPRM1) is related to differences in attachment behavior that emerges following repeated separation from the mother; specifically, infants carrying at least one copy of the minor G allele of the OPRM1 C77G polymorphism show heightened and more persistent separation distress, as well as a pattern of increased contact-seeking behavior directed towards the mother during reunions (at the expense of affiliation with other group members). Research in adult humans has also linked the minor G allele of the analogous OPRM1 A118G polymorphism with greater interpersonal sensitivity. Adopting an interactionist approach, we examined whether OPRM1 A118G genotype and maternal (in)sensitivity are associated with child attachment style, predicting that children carrying the G allele may be more likely to develop an ambivalent attachment pattern in response to less sensitive maternal care. The sample consisted of 191 mothers participating with their children (n = 223) in the Maternal Adversity, Vulnerability and Neurodevelopment (MAVAN) project, a community-based, birth cohort study of Canadian mothers and their children assessed longitudinally across the child's development. Maternal sensitivity was coded from at-home mother-child interactions videotaped when the child was 18 months of age. Child attachment was assessed at 36 months using the Strange Situation paradigm. As predicted, G allele carriers, but not AA homozygotes, showed increasing odds of being classified as ambivalently attached with decreasing levels of maternal sensitivity. Paralleling earlier non-human animal research, this work provides support for the theory that endogenous opioids contribute to the expression of attachment behaviors in humans.

OPRM1 and BDNF polymorphisms associated with a compensatory neurophysiologic signature in knee osteoarthritis patients

OBJECTIVE

The present study investigated the relationship between three genetic polymorphisms of OPRM1 (rs1799971 - A118G and rs1799972 - C17T) and BDNF (rs6265 - C196T) and EEG-measured brain oscillations in Knee Osteoarthritis (KOA) patients.

MATERIALS AND METHODS

We performed a cross-sectional analysis of a cohort study (DEFINE cohort), KOA arm, with 66 patients, considering demographic (age, sex, and education), clinical (pain intensity and duration), OPRM1 (rs1799971 - A118G and rs1799972 - C17T) and BDNF (rs6265 - C196T) genotypes, and electrophysiological measures. Brain oscillations relative power from Delta, Theta, Alpha, Low Alpha, High Alpha, Beta, Low Beta and High Beta oscillations were measured during resting state EEG. Multivariate regression models were used to explore the main brain oscillation predictors of the three genetic polymorphisms.

RESULTS

Our findings demonstrate that Theta and Low Beta oscillations are associated with the variant allele of OPRM1-rs1799971 (A118G) on left frontal and left central regions, respectively, while Alpha brain oscillation is associated with variant genotypes (CT/TT) of BDNF-rs6265 on frontal (decrease of oscillation power) and left central (increase of oscillation power) regions. No significant model was found for OPRM1-rs1799972 (C17T) in addition to the inclusion of pain intensity as a significant predictor of this last model.

CONCLUSION

One potential interpretation for these findings is that polymorphisms of OPRM1 - that is involved with endogenous pain control - lead to increased compensatory oscillatory mechanisms, characterized by increased theta oscillations. Along the same line, polymorphisms of the BDNF lead to decreased alpha oscillations in the frontal area, likely also reflecting the disruption of resting states to also compensate for the increased injury associated with knee OA. It is possible that these polymorphisms require additional brain adaption to the knee OA related injury.

Tapentadol shows lower intrinsic efficacy at µ receptor than morphine and oxycodone

Tapentadol is a centrally acting analgesic with a dual mechanism of action. It acts as an agonist at the µ receptor and inhibitor of noradrenaline reuptake. Clinical trials suggest similar analgesic efficacy of tapentadol, oxycodone, and morphine in acute and chronic pain. Given the limited information about the molecular actions of tapentadol at the µ receptor, we investigated the intrinsic efficacy of tapentadol and compared it with other opioids. β-chlornaltrexamine (β-CNA, 100 nM, 20 min) was used to deplete spare receptors in AtT20 cells stably transfected with human µ receptor wild-type (WT). Opioid-mediated changes in membrane potential were measured in real-time using a membrane potential-sensitive fluorescent dye. Using Black and Leff's operational model, intrinsic efficacy relative to DAMGO was calculated for each opioid. Tapentadol (0.05 ± 0.01) activated the GIRK channel with lesser intrinsic efficacy than morphine (0.17 ± 0.02) and oxycodone (0.16 ± 0.02). We further assessed the signaling of tapentadol in the common µ receptor variants (N40D and A6V) which are associated with altered receptor signaling. We found no difference in the response of tapentadol between these receptor variants.

Population-specific genetic background for the OPRM1 variant rs1799971 (118A>G): implications for genomic medicine and functional analysis

The mu-opioid receptor (MOR, OPRM1) has important roles in diverse functions including reward, addiction, and response to pain treatment. SNP rs1799971 (118A > G, N40D) which occur at a high frequency (40-60%) in Asia and moderate frequency (15%) in samples of European ancestry, is the only common coding variant in the canonical transcript, in non-African populations. Despite extensive studies, the molecular consequences of this variation remained unresolved. The aim of this study was to determine the genetic background of the OPRM1 region of 118G in four representative populations and to assess its potential modulatory effect. Seven common haplotypes with distinct population distribution were identified based on seven SNPs. Three haplotypes carry the 118G and additional highly linked regulatory SNPs (e.g., rs9383689) that could modulate the effect of 118G. Extended analysis in the 1000 Genomes database (n = 2504) revealed a common East Asian-specific haplotype with a different genetic background in which there are no variant alleles for an upstream LD block tagged by the eQTL rs9397171. The major European haplotype specifically includes the eQTL intronic SNP rs62436463 that must have arisen after the split between European and Asian populations. Differentiating between the effect of 118G and these SNPs requires specific experimental approaches. The analysis also revealed a significant increase in two 118A haplotypes with eQTL SNPs associated with drug addiction (rs510769) and obesity (rs9478496) in populations with native Mexican ancestry. Future studies are required to assess the clinical implication of these findings.

From Pharmacology to Physiology: Endocrine Functions of μ-Opioid Receptor Networks

The steady rise in opioid users and abusers has uncovered multiple detrimental health consequences of perturbed opioid receptor signaling, thereby creating the need to better understand the biology of these systems. Among endogenous opioid networks, μ-receptors have received special attention due to their unprecedented biological complexity and broad implications in homeostatic functions. Here, we review the origin, molecular biology, and physiology of endogenous opioids with a special focus on μ-opioid receptor networks within the endocrine system. Moreover, we summarize the current evidence supporting an involvement of the latter in regulating distinct endocrine functions. Finally, we combine these insights to present an integrated perspective on μ-opioid receptor biology and provide an outlook on future studies and unresolved questions in this field.

Flexible Mixture Model Approaches That Accommodate Footprint Size Variability for Robust Detection of Balancing Selection

Long-term balancing selection typically leaves narrow footprints of increased genetic diversity, and therefore most detection approaches only achieve optimal performances when sufficiently small genomic regions (i.e., windows) are examined. Such methods are sensitive to window sizes and suffer substantial losses in power when windows are large. Here, we employ mixture models to construct a set of five composite likelihood ratio test statistics, which we collectively term B statistics. These statistics are agnostic to window sizes and can operate on diverse forms of input data. Through simulations, we show that they exhibit comparable power to the best-performing current methods, and retain substantially high power regardless of window sizes. They also display considerable robustness to high mutation rates and uneven recombination landscapes, as well as an array of other common confounding scenarios. Moreover, we applied a specific version of the B statistics, termed B2, to a human population-genomic data set and recovered many top candidates from prior studies, including the then-uncharacterized STPG2 and CCDC169-SOHLH2, both of which are related to gamete functions. We further applied B2 on a bonobo population-genomic data set. In addition to the MHC-DQ genes, we uncovered several novel candidate genes, such as KLRD1, involved in viral defense, and SCN9A, associated with pain perception. Finally, we show that our methods can be extended to account for multiallelic balancing selection and integrated the set of statistics into open-source software named BalLeRMix for future applications by the scientific community.

The μ-opioid receptor gene A118G polymorphism is associated with insecure attachment in children with disruptive mood regulation disorder and their mothers

BACKGROUND

The A118G single nucleotide polymorphism (SNP) of the μ-opioid receptor gene, with high expression of the A allele and low expression of the G allele, has been associated with emotional/behavioral dysregulation and depressive disorders and is recognized as a mediator of affiliative behavior. No study has thus far investigated this SNP in school-age children with disruptive mood regulation disorder (DMDD). This study compared a sample of healthy children and their mothers with a sample of children with DMDD and their mothers, evaluating whether insecure attachment and psychopathological symptoms are associated with A allele- or G allele-carrying mothers and children and whether caregiving capacities are associated with A allele- or G allele-carrying mothers.

METHODS

For evaluation of their psychopathological symptoms and attachment styles, mothers filled out the CBCL/6-18, the SCL-90-R, and the ECR. To evaluate the types of relationship children were experiencing with their mothers, children filled out the ECR-revised child version and the PBI. Genotypic analyses were conducted on DNA samples obtained by buccal swabbing from children and mothers.

RESULTS

An insecure attachment style was more frequent in mothers and children carrying the G allele (G/G + A/G genotypes). In the clinical sample, G allele-carrying children scored higher than homozygous A/A ones on the subscales of Withdrawal and Conduct Problems. G-carrying mothers showed higher interpersonal sensitivity, depression, hostility, and paranoid ideation and provided less care than A/A mothers.

CONCLUSIONS

This study offers new insights into the associations between the A118G SNP of the μ-opioid receptor gene and emotional/behavioral functioning, attachment style in children, and psychopathology and caregiving ability in mothers.

In Vitro Effects of Ligand Bias on Primate Mu Opioid Receptor Downstream Signaling

Interest has emerged in biased agonists at the mu opioid receptor (MOR) as a possible means for maintaining potent analgesis with reduced side effect profiles. While approaches measuring in vitro biased agonism are used in the development of these compounds, their therapeutic utility will ultimately be determined by in vivo functional effects. Nonhuman primates (NHPs) are the most translational model for evaluating the behavioral effects of candidate medications, but biased signaling of these drugs at NHP MOR receptors has been unstudied. The goal of the current work was to characterize MOR ligand bias in rhesus macaques, focusing on agonists that have previously been reported to show different patterns of biased agonism in rodents and humans. Downstream signaling pathways that responded to MOR activation were identified using a luciferase reporter array. Concentration-response curves for specific pathways (cAMP, NF-ĸB, MAPK/JNK) were generated using six agonists previously reported to differ in terms of signaling bias at rodent and human MORs. Using DAMGO as a reference ligand, relative cAMP, NF-ĸB and MAPK/JNK signaling by morphine, endomorphin-1, and TRV130 were found to be comparable between species. Further, the bias patterns of across ligands for NF-ĸB and MAPK/JNK were largely similar between species. There was a high degree of concordance between rhesus macaque and human MOR receptor signaling bias for all agonists tested, further demonstrating their utility for future translational behavioral studies.

Endogenous opiates and behavior: 2017

This paper is the fortieth consecutive installment of the annual anthological review of research concerning the endogenous opioid system, summarizing articles published during 2017 that studied the behavioral effects of molecular, pharmacological and genetic manipulation of opioid peptides and receptors as well as effects of opioid/opiate agonists and antagonists. The review is subdivided into the following specific topics: molecular-biochemical effects and neurochemical localization studies of endogenous opioids and their receptors (1), the roles of these opioid peptides and receptors in pain and analgesia in animals (2) and humans (3), opioid-sensitive and opioid-insensitive effects of nonopioid analgesics (4), opioid peptide and receptor involvement in tolerance and dependence (5), stress and social status (6), learning and memory (7), eating and drinking (8), drug abuse and alcohol (9), sexual activity and hormones, pregnancy, development and endocrinology (10), mental illness and mood (11), seizures and neurologic disorders (12), electrical-related activity and neurophysiology (13), general activity and locomotion (14), gastrointestinal, renal and hepatic functions (15), cardiovascular responses (16), respiration and thermoregulation (17), and immunological responses (18).

Fine-Scale Characterization of Genomic Structural Variation in the Human Genome Reveals Adaptive and Biomedically Relevant Hotspots

Genomic structural variants (SVs) are distributed nonrandomly across the human genome. The "hotspots" of SVs have been implicated in evolutionary innovations, as well as medical conditions. However, the evolutionary and biomedical features of these hotspots remain incompletely understood. Here, we analyzed data from 2,504 genomes to construct a refined map of 1,148 SV hotspots in human genomes. We confirmed that segmental duplication-related nonallelic homologous recombination is an important mechanistic driver of SV hotspot formation. However, to our surprise, we also found that a majority of SVs in hotspots do not form through such recombination-based mechanisms, suggesting diverse mechanistic and selective forces shaping hotspots. Indeed, our evolutionary analyses showed that the majority of SV hotspots are within gene-poor regions and evolve under relaxed negative selection or neutrality. However, we still found a small subset of SV hotspots harboring genes that are enriched for anthropologically crucial functions and evolve under geography-specific and balancing adaptive forces. These include two independent hotspots on different chromosomes affecting alpha and beta hemoglobin gene clusters. Biomedically, we found that the SV hotspots coincide with breakpoints of clinically relevant, large de novo SVs, significantly more often than genome-wide expectations. For example, we showed that the breakpoints of multiple large SVs, which lead to idiopathic short stature, coincide with SV hotspots. Therefore, the mutational instability in SV hotpots likely enables chromosomal breaks that lead to pathogenic structural variation formations. Overall, our study contributes to a better understanding of the mutational and adaptive landscape of the genome.

Common marmoset (Callithrix jacchus) personality, subjective well-being, hair cortisol level and AVPR1a, OPRM1, and DAT genotypes

We studied personality, subjective well-being, and hair cortisol level, in common marmosets Callithrix jacchus, a small, cooperatively breeding New World monkey, by examining their associations with one another and genotypes. Subjects were 68 males and 9 females that lived in the RIKEN Center for Life Science Technologies. Personality and subjective well-being were assessed by keeper ratings on two questionnaires, hair samples were obtained to assay cortisol level and buccal swabs were used to assess AVPR1a, OPRM1 and DAT genotypes. Three personality domains-Dominance, Sociability, and Neuroticism-were identified. Consistent with findings in other species, Sociability and Neuroticism were related to higher and lower subjective well-being, respectively. Sociability was also associated with higher hair cortisol levels. The personality domains and hair cortisol levels were heritable and associated with genotypes: the short form of AVPR1a was associated with lower Neuroticism and the AA genotype of the A111T SNP of OPRM1 was related to lower Dominance, lower Neuroticism, and higher hair cortisol level. Some genetic associations were not in directions that one would expect given findings in other species. These findings provide insights into the proximate and ultimate bases of personality in common marmosets, other primates and humans.