Fabricated devices for performing bacterial-fungal interaction experiments across scales.

Diverse and complex microbiomes are found in virtually every environment on Earth. Bacteria and fungi often co-dominate environmental microbiomes, and there is growing recognition that bacterial-fungal interactions (BFI) have significant impacts on the functioning of their associated microbiomes, environments, and hosts. Investigating BFI in vitro remains a challenge, particularly when attempting to examine interactions at multiple scales of system complexity. Fabricated devices can provide control over both biotic composition and abiotic factors within an experiment to enable the characterization of diverse BFI phenotypes such as modulation of growth rate, production of biomolecules, and alterations to physical movements. Engineered devices ranging from microfluidic chips to simulated rhizosphere systems have been and will continue to be invaluable to BFI research, and it is anticipated that such devices will continue to be developed for diverse applications in the field. This will allow researchers to address specific questions regarding the nature of BFI and how they impact larger microbiome and environmental processes such as biogeochemical cycles, plant productivity, and overall ecosystem resilience. Devices that are currently used for experimental investigations of bacteria, fungi, and BFI are discussed herein along with some of the associated challenges and several recommendations for future device design and applications.

Atp1a2 and Kcnj9 are candidate genes underlying sensitivity to oxycodone-induced locomotor activation and withdrawal-induced anxiety-like behaviors in C57BL/6 substrains.

Opioid use disorder is heritable, yet its genetic etiology is largely unknown. C57BL/6J and C57BL/6NJ mouse substrains exhibit phenotypic diversity in the context of limited genetic diversity which together can facilitate genetic discovery. Here, we found C57BL/6NJ mice were less sensitive to oxycodone (OXY)-induced locomotor activation versus C57BL/6J mice in a conditioned place preference paradigm. Narrow-sense heritability was estimated at 0.22-0.31, implicating suitability for genetic analysis. Quantitative trait locus (QTL) mapping in an F2 cross identified a chromosome 1 QTL explaining 7-12% of the variance in OXY locomotion and anxiety-like withdrawal in the elevated plus maze. A second QTL for EPM withdrawal behavior on chromosome 5 near Gabra2 (alpha-2 subunit of GABA-A receptor) explained 9% of the variance. To narrow the chromosome 1 locus, we generated recombinant lines spanning 163-181 Mb, captured the QTL for OXY locomotor traits and withdrawal, and fine-mapped a 2.45-Mb region (170.16-172.61 Mb). Transcriptome analysis identified five, localized striatal cis-eQTL transcripts and two were confirmed at the protein level (KCNJ9, ATP1A2). Kcnj9 codes for a potassium channel (GIRK3) that is a major effector of mu opioid receptor signaling. Atp1a2 codes for a subunit of a Na+/K+ ATPase enzyme that regulates neuronal excitability and shows functional adaptations following chronic opioid administration. To summarize, we identified two candidate genes underlying the physiological and behavioral properties of opioids, with direct preclinical relevance to investigators employing these widely used substrains and clinical relevance to human genetic studies of opioid use disorder.

A Practical Approach to Using the Genomic Standards Consortium MIxS Reporting Standard for Comparative Genomics and Metagenomics.

Comparative analysis of (meta)genomes necessitates aggregation, integration, and synthesis of well-annotated data using standards. The Genomic Standards Consortium (GSC) collaborates with the research community to develop and maintain the Minimum Information about any (x) Sequence (MIxS) reporting standard for genomic data. To facilitate the use of the GSC's MIxS reporting standard, we provide a description of the structure and terminology, how to navigate ontologies for required terms in MIxS, and demonstrate practical usage through a soil metagenome example.

Associated bacterial communities, confrontation studies, and comparative genomics reveal important interactions between Morchella with Pseudomonas spp.

Members of the fungal genus Morchella are widely known for their important ecological roles and significant economic value. In this study, we used amplicon and genome sequencing to characterize bacterial communities associated with sexual fruiting bodies from wild specimens, as well as vegetative mycelium and sclerotia obtained from Morchella isolates grown in vitro. These investigations included diverse representatives from both Elata and Esculenta Morchella clades. Unique bacterial community compositions were observed across the various structures examined, both within and across individual Morchella isolates or specimens. However, specific bacterial taxa were frequently detected in association with certain structures, providing support for an associated core bacterial community. Bacteria from the genus Pseudomonas and Ralstonia constituted the core bacterial associates of Morchella mycelia and sclerotia, while other genera (e.g., Pedobacter spp., Deviosa spp., and Bradyrhizobium spp.) constituted the core bacterial community of fruiting bodies. Furthermore, the importance of Pseudomonas as a key member of the bacteriome was supported by the isolation of several Pseudomonas strains from mycelia during in vitro cultivation. Four of the six mycelial-derived Pseudomonas isolates shared 16S rDNA sequence identity with amplicon sequences recovered directly from the examined fungal structures. Distinct interaction phenotypes (antagonistic or neutral) were observed in confrontation assays between these bacteria and various Morchella isolates. Genome sequences obtained from these Pseudomonas isolates revealed intriguing differences in gene content and annotated functions, specifically with respect to toxin-antitoxin systems, cell adhesion, chitinases, and insecticidal toxins. These genetic differences correlated with the interaction phenotypes. This study provides evidence that Pseudomonas spp. are frequently associated with Morchella and these associations may greatly impact fungal physiology.

The endohyphal microbiome: current progress and challenges for scaling down integrative multi-omic microbiome research.

As microbiome research has progressed, it has become clear that most, if not all, eukaryotic organisms are hosts to microbiomes composed of prokaryotes, other eukaryotes, and viruses. Fungi have only recently been considered holobionts with their own microbiomes, as filamentous fungi have been found to harbor bacteria (including cyanobacteria), mycoviruses, other fungi, and whole algal cells within their hyphae. Constituents of this complex endohyphal microbiome have been interrogated using multi-omic approaches. However, a lack of tools, techniques, and standardization for integrative multi-omics for small-scale microbiomes (e.g., intracellular microbiomes) has limited progress towards investigating and understanding the total diversity of the endohyphal microbiome and its functional impacts on fungal hosts. Understanding microbiome impacts on fungal hosts will advance explorations of how "microbiomes within microbiomes" affect broader microbial community dynamics and ecological functions. Progress to date as well as ongoing challenges of performing integrative multi-omics on the endohyphal microbiome is discussed herein. Addressing the challenges associated with the sample extraction, sample preparation, multi-omic data generation, and multi-omic data analysis and integration will help advance current knowledge of the endohyphal microbiome and provide a road map for shrinking microbiome investigations to smaller scales. Video Abstract.

A centralized resource for bacterial-fungal interactions research.

Research on bacterial-fungal interactions (BFIs) has revealed that fungi and bacteria frequently interact with one another within diverse ecosystems and microbiomes. Assessing the current state of knowledge within the field of BFI research, particularly with respect to what interactions between bacteria and fungi have been previously described, is very challenging and time consuming. This is largely due to a lack of any centralized resource, with reports of BFIs being spread across publications in numerous journals using non-standardized text to describe the relationships. To address this issue, we have developed the BFI Research Portal, a publicly accessible database of previously reported interactions between bacterial and fungal taxa to serve as a centralized resource for the field. Users can query bacterial or fungal taxa to see what members from the other kingdom have been observed as interaction partners. Search results are accompanied by interactive and intuitive visual outputs, and the database is a dynamic resource that will be updated as new BFIs are reported.

Advances and Challenges in Fluorescence in situ Hybridization for Visualizing Fungal Endobacteria.

Several bacteria have long been known to interact intimately with fungi, but molecular approaches have only recently uncovered how cosmopolitan these interactions are in nature. Currently, bacterial-fungal interactions (BFI) are inferred based on patterns of co-occurrence in amplicon sequencing investigations. However, determining the nature of these interactions, whether the bacteria are internally or externally associated, remains a grand challenge in BFI research. Fluorescence in situ hybridization (FISH) is a robust method that targets unique sequences of interest which can be employed for visualizing intra-hyphal targets, such as mitochondrial organelles or, as in this study, bacteria. We evaluate the challenges and employable strategies to resolve intra-hyphal BFI to address pertinent criteria in BFI research, such as culturing media, spatial distribution of bacteria, and abundance of bacterial 16S rRNA copies for fluorescent labeling. While these experimental factors influence labeling and detection of endobacteria, we demonstrate how to overcome these challenges thorough permeabilization, appropriate media choice, and targeted amplification using hybridization chain reaction FISH. Such microscopy imaging approaches can now be utilized by the broader research community to complement sequence-based investigations and provide more conclusive evidence on the nature of specific bacterial-fungal relationships.

Widespread bacterial diversity within the bacteriome of fungi.

Knowledge of associations between fungal hosts and their bacterial associates has steadily grown in recent years as the number and diversity of examinations have increased, but current knowledge is predominantly limited to a small number of fungal taxa and bacterial partners. Here, we screened for potential bacterial associates in over 700 phylogenetically diverse fungal isolates, representing 366 genera, or a tenfold increase compared with previously examined fungal genera, including isolates from several previously unexplored phyla. Both a 16 S rDNA-based exploration of fungal isolates from four distinct culture collections spanning North America, South America and Europe, and a bioinformatic screen for bacterial-specific sequences within fungal genome sequencing projects, revealed that a surprisingly diverse array of bacterial associates are frequently found in otherwise axenic fungal cultures. We demonstrate that bacterial associations with diverse fungal hosts appear to be the rule, rather than the exception, and deserve increased consideration in microbiome studies and in examinations of microbial interactions.

A quantitative trait variant in Gabra2 underlies increased methamphetamine stimulant sensitivity.

Psychostimulant (methamphetamine, cocaine) use disorders have a genetic component that remains mostly unknown. We conducted genome-wide quantitative trait locus (QTL) analysis of methamphetamine stimulant sensitivity. To facilitate gene identification, we employed a Reduced Complexity Cross between closely related C57BL/6 mouse substrains and examined maximum speed and distance traveled over 30 min following methamphetamine (2 mg/kg, i.p.). For maximum methamphetamine-induced speed following the second and third administration, we identified a single genome-wide significant QTL on chromosome 11 that peaked near the Cyfip2 locus (LOD = 3.5, 4.2; peak = 21 cM [36 Mb]). For methamphetamine-induced distance traveled following the first and second administration, we identified a genome-wide significant QTL on chromosome 5 that peaked near a functional intronic indel in Gabra2 coding for the alpha-2 subunit of the GABA-A receptor (LOD = 3.6-5.2; peak = 34-35 cM [66-67 Mb]). Striatal cis-expression QTL mapping corroborated Gabra2 as a functional candidate gene underlying methamphetamine-induced distance traveled. CRISPR/Cas9-mediated correction of the mutant intronic deletion on the C57BL/6J background to the wild-type C57BL/6NJ allele was sufficient to reduce methamphetamine-induced locomotor activity toward the wild-type C57BL/6NJ-like level, thus validating the quantitative trait variant (QTV). These studies show the power and efficiency of Reduced Complexity Crosses in identifying causal variants underlying complex traits. Functionally restoring Gabra2 expression decreased methamphetamine stimulant sensitivity and supports preclinical and human genetic studies implicating the GABA-A receptor in psychostimulant addiction-relevant traits. Importantly, our findings have major implications for studying psychostimulants in the C57BL/6J strain-the gold standard strain in biomedical research.

Sex Differences in Behavioral and Brainstem Transcriptomic Neuroadaptations following Neonatal Opioid Exposure in Outbred Mice.

The opioid epidemic led to an increase in the number of neonatal opioid withdrawal syndrome (NOWS) cases in infants born to opioid-dependent mothers. Hallmark features of NOWS include weight loss, severe irritability, respiratory problems, and sleep fragmentation. Mouse models provide an opportunity to identify brain mechanisms that contribute to NOWS. Neonatal outbred Swiss Webster Cartworth Farms White (CFW) mice were administered morphine (15 mg/kg, s.c.) twice daily from postnatal day 1 (P1) to P14, an approximation of the third trimester of human gestation. Female and male mice underwent behavioral testing on P7 and P14 to determine the impact of opioid exposure on anxiety and pain sensitivity. Ultrasonic vocalizations (USVs) and daily body weights were also recorded. Brainstems containing pons and medulla were collected during morphine withdrawal on P14 for RNA sequencing. Morphine induced weight loss from P2 to P14, which persisted during adolescence (P21) and adulthood (P50). USVs markedly increased at P7 in females, emerging earlier than males. On P7 and P14, both morphine-exposed female and male mice displayed hyperalgesia on the hot plate and tail-flick assays, with females showing greater hyperalgesia than males. Morphine-exposed mice exhibited increased anxiety-like behavior in the open-field arena on P21. Transcriptome analysis of the brainstem, an area implicated in opioid withdrawal and NOWS, identified pathways enriched for noradrenergic signaling in females and males. We also found sex-specific pathways related to mitochondrial function and neurodevelopment in females and circadian entrainment in males. Sex-specific transcriptomic neuroadaptations implicate unique neurobiological mechanisms underlying NOWS-like behaviors.

Systems genetic analysis of binge-like eating in a C57BL/6J x DBA/2J-F2 cross.

Binge eating is a heritable trait associated with eating disorders and refers to the rapid consumption of a large quantity of energy-dense food that is, associated with loss of control and negative affect. Binge eating disorder is the most common eating disorder in the United States; however, the genetic basis is unknown. We previously identified robust mouse inbred strain differences between C57BL/6J and DBA/2J in binge-like eating of sweetened palatable food in an intermittent access, conditioned place preference paradigm. To map the genetic basis of changes in body weight and binge-like eating (BLE) and to identify candidate genes, we conducted quantitative trait locus (QTL) analysis in 128 C57BL/6J x DBA/2J-F2 mice combined with PheQTL and trait covariance analysis in GeneNetwork2 using legacy BXD-RI trait datasets. We identified a QTL on Chromosome 18 influencing changes in body weight across days in females (log of the odds [LOD] = 6.3; 1.5-LOD: 3-12 cM) that contains the candidate gene Zeb1. We also identified a sex-combined QTL influencing initial palatable food intake on Chromosome 5 (LOD = 5.8; 1.5-LOD: 21-28 cM) that contains the candidate gene Lcorl and a second QTL influencing escalated palatable food intake on Chromosome 6 in males (LOD = 5.4; 1.5-LOD: 50-59 cM) that contains the candidate genes Adipor2 and Plxnd1. Finally, we identified a suggestive QTL in females for slope of BLE on distal Chromosome 18 (LOD = 4.1; p = 0.055; 1.5-LOD: 23-35 cM). Future studies will use BXD-RI strains to fine map loci and support candidate gene nomination for gene editing.

Democratization of fungal highway columns as a tool to investigate bacteria associated with soil fungi.

Bacteria-fungi interactions (BFIs) are essential in ecosystem functioning. These interactions are modulated not only by local nutritional conditions but also by the physicochemical constraints and 3D structure of the environmental niche. In soils, the unsaturated and complex nature of the substrate restricts the dispersal and activity of bacteria. Under unsaturated conditions, some bacteria engage with filamentous fungi in an interaction (fungal highways) in which they use fungal hyphae to disperse. Based on a previous experimental device to enrich pairs of organisms engaging in this interaction in soils, we present here the design and validation of a modified version of this sampling system constructed using additive printing. The 3D printed devices were tested using a novel application in which a target fungus, the common coprophilous fungus Coprinopsis cinerea, was used as bait to recruit and identify bacterial partners using its mycelium for dispersal. Bacteria of the genera Pseudomonas, Sphingobacterium and Stenotrophomonas were highly enriched in association with C. cinerea. Developing and producing these new easy-to-use tools to investigate how bacteria overcome dispersal limitations in cooperation with fungi is important to unravel the mechanisms by which BFIs affect processes at an ecosystem scale in soils and other unsaturated environments.

Challenges in Bioinformatics Workflows for Processing Microbiome Omics Data at Scale.

The nascent field of microbiome science is transitioning from a descriptive approach of cataloging taxa and functions present in an environment to applying multi-omics methods to investigate microbiome dynamics and function. A large number of new tools and algorithms have been designed and used for very specific purposes on samples collected by individual investigators or groups. While these developments have been quite instructive, the ability to compare microbiome data generated by many groups of researchers is impeded by the lack of standardized application of bioinformatics methods. Additionally, there are few examples of broad bioinformatics workflows that can process metagenome, metatranscriptome, metaproteome and metabolomic data at scale, and no central hub that allows processing, or provides varied omics data that are findable, accessible, interoperable and reusable (FAIR). Here, we review some of the challenges that exist in analyzing omics data within the microbiome research sphere, and provide context on how the National Microbiome Data Collaborative has adopted a standardized and open access approach to address such challenges.

Multiple Cases of Bacterial Sequence Erroneously Incorporated Into Publicly Available Chloroplast Genomes.

Public sequencing databases are invaluable resources to biological researchers, but assessing data veracity as well as the curation and maintenance of such large collections of data can be challenging. Genomes of eukaryotic organelles, such as chloroplasts and other plastids, are particularly susceptible to assembly errors and misrepresentations in these databases due to their close evolutionary relationships with bacteria, which may co-occur within the same environment, as can be the case when sequencing plants. Here, based on sequence similarities with bacterial genomes, we identified several suspicious chloroplast assemblies present in the National Institutes of Health (NIH) Reference Sequence (RefSeq) collection. Investigations into these chloroplast assemblies reveal examples of erroneous integration of bacterial sequences into chloroplast ribosomal RNA (rRNA) loci, often within the rRNA genes, presumably due to the high similarity between plastid and bacterial rRNAs. The bacterial lineages identified within the examined chloroplasts as the most likely source of contamination are either known associates of plants, or co-occur in the same environmental niches as the examined plants. Modifications to the methods used to process untargeted 'raw' shotgun sequencing data from whole genome sequencing efforts, such as the identification and removal of bacterial reads prior to plastome assembly, could eliminate similar errors in the future.

Development of Anti-Yersinia pestis Human Antibodies with Features Required for Diagnostic and Therapeutic Applications.

BACKGROUND: Yersinia pestis is a category A infective agent that causes bubonic, septicemic, and pneumonic plague. Notably, the acquisition of antimicrobial or multidrug resistance through natural or purposed means qualifies Y. pestis as a potential biothreat agent. Therefore, high-quality antibodies designed for accurate and sensitive Y. pestis diagnostics, and therapeutics potentiating or replacing traditional antibiotics are of utmost need for national security and public health preparedness. METHODS: Here, we describe a set of human monoclonal immunoglobulins (IgG1s) targeting Y. pestis fraction 1 (F1) antigen, previously derived from in vitro evolution of a phage-display library of single-chain antibodies (scFv). We extensively characterized these antibodies and their effect on bacterial and mammalian cells via: ELISA, flow cytometry, mass spectrometry, spectroscopy, and various metabolic assays. RESULTS: Two of our anti-F1 IgG (αF1Ig 2 and αF1Ig 8) stood out for high production yield, specificity, and stability. These two antibodies were additionally attractive in that they displayed picomolar affinity, did not compete when binding Y. pestis, and retained immunoreactivity upon chemical derivatization. Most importantly, these antibodies detected <1,000 Y. pestis cells in sandwich ELISA, did not harm respiratory epithelial cells, induced Y. pestis agglutination at low concentration (350 nM), and caused apparent reduction in cell growth when radiolabeled at a nonagglutinating concentration (34 nM). CONCLUSION: These antibodies are amenable to the development of accurate and sensitive diagnostics and immuno/radioimmunotherapeutics.

The effect of the demyelinating agent cuprizone on binge-like eating of sweetened palatable food in female and male C57BL/6 substrains.

Binge eating is a heritable symptom of eating disorders with an unknown genetic etiology. Rodent models for binge-like eating (BLE) of palatable food permit the study of genetic and biological mechanisms. We previously genetically mapped a coding mutation in Cyfip2 associated with increased BLE of sweetened palatable food in the C57BL/6NJ versus C57BL/6J substrain. The increase in BLE in C57BL/6NJ mice was associated with a decrease in transcription of genes enriched for myelination in the striatum. Here, we tested the hypothesis that decreasing myelin levels with the demyelinating agent cuprizone would enhance BLE. Mice were treated with a 0.3% cuprizone home cage diet for two weeks. Cuprizone induced similar weight loss in both substrains and sexes that recovered within 48 h after removal of cuprizone. Following a three-week recovery period, mice were trained for BLE in an intermittent, limited access procedure. Surprisingly, cuprizone significantly reduced BLE in male but not female C57BL/6NJ mice while having no effect in C57BL/6J mice. Cuprizone also reduced myelin basic protein (MBP) at seven weeks post-cuprizone removal while having no effect on myelin-associated glycoprotein at this time point. C57BL/6NJ mice also showed less MBP than C57BL/6J mice. There were no statistical interactions of Treatment with Sex on MBP levels, indicating that differences in MBP reduction are unlikely to account for sex differences in BLE. To summarize, cuprizone induced an unexpected, significant reduction in BLE in C57BL/6NJ males, which could indicate genotype-dependent sex differences in the biological mechanisms of BLE.

Complete Genome Sequence of a Parabacteroides distasonis Strain (CavFT hAR46) Isolated from a Gut Wall-Cavitating Microlesion in a Patient with Severe Crohn's Disease.

Crohn's disease (CD) is a chronic inflammatory bowel disease (IBD) of the digestive tract in humans. There is evidence that Parabacteroides distasonis could contribute to IBD. Here, we present the complete genome sequence of a strain designated CavFT-hAR46, which was isolated from a gut intramural cavernous fistulous tract (CavFT) microlesion in a CD patient.

Cyfip1 Haploinsufficiency Increases Compulsive-Like Behavior and Modulates Palatable Food Intake in Mice: Dependence on Cyfip2 Genetic Background, Parent-of Origin, and Sex.

Binge eating (BE) is a heritable trait associated with eating disorders and involves episodes of rapid, large amounts of food consumption. We previously identified cytoplasmic FMR1-interacting protein 2 (Cyfip2) as a genetic factor underlying compulsive-like BE in mice. CYFIP2 is a homolog of CYFIP1 which is one of four paternally-deleted genes in patients with Type I Prader-Willi Syndrome (PWS), a neurodevelopmental disorder whereby 70% of cases involve paternal 15q11-q13 deletion. PWS symptoms include hyperphagia, obesity (if untreated), cognitive deficits, and obsessive-compulsive behaviors. We tested whether Cyfip1 haploinsufficiency (+/-) would enhance compulsive-like behavior and palatable food (PF) intake in a parental origin- and sex-dependent manner on two Cyfip2 genetic backgrounds, including the BE-prone C57BL/6N (Cyfip2N/N) background and the BE-resistant C57BL/6J (Cyfip2J/J) background. Cyfip1+/- mice showed increased compulsive-like behavior on both backgrounds and increased PF intake on the Cyfip2N/N background. In contrast, maternal Cyfip1 haploinsufficiency on the BE-resistant Cyfip2J/J background induced a robust escalation in PF intake in wild-type Cyfip1J/J males while having no effect in Cyfip1J/- males. Notably, induction of behavioral phenotypes in wild-type males following maternal Fmr1+/- has previously been reported. In the hypothalamus, there was a paternally-enhanced reduction in CYFIP1 protein whereas in the nucleus accumbens, there was a maternally-enhanced reduction in CYFIP1 protein. Nochange in FMR1 protein (FMRP) was observed in Cyfip1+/- mice, regardless of parental origin. To summarize, Cyfip1 haploinsufficiency increased compulsive-like behavior and induced genetic background-dependent, sex-dependent, and parent-of-origin-dependent effects on PF consumption and CYFIP1 expression that could have relevance for neurodevelopmental and neuropsychiatric disorders.

C57BL/6 substrain differences in inflammatory and neuropathic nociception and genetic mapping of a major quantitative trait locus underlying acute thermal nociception.

Sensitivity to different pain modalities has a genetic basis that remains largely unknown. Employing closely related inbred mouse substrains can facilitate gene mapping of nociceptive behaviors in preclinical pain models. We previously reported enhanced sensitivity to acute thermal nociception in C57BL/6J (B6J) versus C57BL/6N (B6N) substrains. Here, we expanded on nociceptive phenotypes and observed an increase in formalin-induced inflammatory nociceptive behaviors and paw diameter in B6J versus B6N mice (Charles River Laboratories). No strain differences were observed in mechanical or thermal hypersensitivity or in edema following the Complete Freund's Adjuvant model of inflammatory pain, indicating specificity in the inflammatory nociceptive stimulus. In the chronic constrictive nerve injury, a model of neuropathic pain, no strain differences were observed in baseline mechanical threshold or in mechanical hypersensitivity up to one month post-chronic constrictive nerve injury. We replicated the enhanced thermal nociception in the 52.5°C hot plate test in B6J versus B6N mice from The Jackson Laboratory. Using a B6J × B6N-F2 cross (N = 164), we mapped a major quantitative trait locus underlying hot plate sensitivity to chromosome 7 that peaked at 26 Mb (log of the odds [LOD] = 3.81, p < 0.01; 8.74 Mb-36.50 Mb) that was more pronounced in males. Genes containing expression quantitative trait loci associated with the peak nociceptive marker that are implicated in pain and inflammation include Ryr1, Cyp2a5, Pou2f2, Clip3, Sirt2, Actn4, and Ltbp4 (false discovery rate < 0.05). Future studies involving positional cloning and gene editing will determine the quantitative trait gene(s) and potential pleiotropy of this locus across pain modalities.