SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Amplification-based gene characterization allows for sensitive and specific quantification of functional genes. There is often a large diversity of genes represented for functional gene groups, and multiple primers may be necessary to target associated genes.
Genes belonging to the same functional group may include numerous and variable gene sequences, making characterizing and quantifying difficult. Therefore, high-throughput design tools are needed to simultaneously create primers for improved quantification of target genes. We developed MetaFunPrimer, a bioinformatic pipeline, to design primers for numerous genes of interest. This tool also enables gene target prioritization based on ranking the presence of genes in user-defined references, such as environment-specific metagenomes. Given inputs of protein and nucleotide sequences for gene targets of interest and an accompanying set of reference metagenomes or genomes, MetaFunPrimer generates primers for ranked genes of interest. To demonstrate the usage and benefits of MetaFunPrimer, a total of 78 primer pairs were designed to target observed ammonia monooxygenase subunit A (amoA) genes of ammonia-oxidizing bacteria (AOB) in 1,550 publicly available soil metagenomes. We demonstrate computationally that these amoA-AOB primers can cover 94% of the amoA-AOB genes observed in the 1,550 soil metagenomes compared with a 49% estimated coverage by previously published primers. Finally, we verified the utility of these primer sets in incubation experiments that used long-term nitrogen fertilized or unfertilized soils. High-throughput quantitative PCR (qPCR) results and statistical analyses showed significant differences in relative quantification patterns between the two soils, and subsequent absolute quantifications also confirmed that target genes enumerated by six selected primer pairs were significantly more abundant in the nitrogen-fertilized soils. This new tool gives microbial ecologists a new approach to assess functional gene abundance and related microbial community dynamics quickly and affordably.
Fixed or frozen testicular tissues from five rats per group were analyzed by: subjective and quantitative evaluations of spermatogenesis; Northern-blot analysis for expression of selected genes; quantification of testosterone and receptors for LH; and morphometric analysis of Leydig cells. Based on observations of fixed tissue, it was evident that some rats in the flight and vivarium groups had testicular abnormalities unassociated with treatment, and probably existing when they were assigned randomly to the four treatment groups; the simulated-launch group contained no abnormal rat. Lesions induced in testes of caudal-elevation rats precluded discernment of any pre-existing abnormality. Considering rats without pre-existing abnormalities, diameter of seminiferous tubules and numbers of germ cells per tubule cross section were lower (E less than 0.05) in flight rats than in simulated-launch or vivarium rats. However, ratios of germ cells to each other, or to Sertoli cells, and number of homogenization-resistant spermatids did not differ from values for simulated-launch or vivarium controls. There was no effect of flight on normal expression of testis-specific hsp gene products, or evidence for production of stress-inducible transcripts of the hsp70 or hsp90 genes. Concentration of receptors for rLH in testicular tissue, and surface densities of smooth endoplasmic reticulum and peroxisomes in Leydig cells, were similar in flight and simulated-launch rats. However, concentrations of testosterone in testicular tissue or peripheral blood plasma were reduced (P less than 0.05) in flight rats to less than 20 percent of values for simulated-launch or vivarium controls. Thus, spermatogenesis was essentially normal in flight rats, but production of testosterone was severely depressed. Sequela of reduced androgen production on turnover of muscle and bone should be considered when interpreting data from mammals exposed to microgravity.
Xanthohumol (XN), a polyphenol from hops (Humulus lupulus), exhibits antioxidant, anti-inflammatory, antihyperlipidemic, and chemo-preventive activity. Preclinical evidence suggests gut microbiota are critical to mediating some of these bioactivities. Nevertheless, its precise impact on human gut microbiota, particularly at supplemental doses, remains poorly characterized. We evaluated 200 mg/day XN for 3 weeks on human gut microbiota in a eubiotic and dysbiotic model using the Simulator of the Human Intestinal Microbial Ecosystem (SHIME®). Functional assessments of microbiota included quantification of XN metabolites, short-chain fatty acids (SCFAs), and untargeted metabolomics of the digestive metabolome. Bacterial composition was assessed by 16S rRNA gene sequencing. XN reduced alpha-diversity and short-chain fatty acid production in both models, as well as altered taxa abundance variably between models. XN disrupted bile acid metabolism through inhibition of microbial bile salt hydrolase (BSH). The modulation of bile acid metabolism has important implications for host-level bioactivity of XN.
Advances in RNA sequencing technologies have achieved an unprecedented accuracy in the quantification of mRNA isoforms, but our knowledge of isoform-specific functions has lagged behind. There is a need to understand the functional consequences of differential splicing, which could be supported by the generation of accurate and comprehensive isoform-specific gene ontology annotations. We present isoform interpretation, a method that uses expectation–maximization to infer isoform-specific functions based on the relationship between sequence and functional isoform similarity. We predicted isoform-specific functional annotations for 85 617 isoforms of 17 900 protein-coding human genes spanning a range of 17 430 distinct gene ontology terms. Comparison with a gold-standard corpus of manually annotated human isoform functions showed that isoform interpretation significantly outperforms state-of-the-art competing methods. We provide experimental evidence that functionally related isoforms predicted by isoform interpretation show a higher degree of domain sharing and expression correlation than functionally related genes. We also show that isoform sequence similarity correlates better with inferred isoform function than with gene-level function.
Abstract Protein activity, abundance, and stability can be regulated by post-translational modification including ubiquitination. Ubiquitination is conserved among eukaryotes and plays a central role in modulating cellular function; yet, we lack comprehensive catalogs of proteins that are modified by ubiquitin in plants. In this study, we describe an antibody-based approach to enrich ubiquitinated peptides coupled with isobaric labeling to enable quantification of up to 18-multiplexed samples. This approach identified 17,940 ubiquitinated lysine sites arising from 6,453 proteins from Arabidopsis (Arabidopsis thaliana) primary roots, seedlings, and rosette leaves. Gene ontology analysis indicated that ubiquitinated proteins are associated with numerous biological processes including hormone signaling, plant defense, protein homeostasis, and metabolism. We determined ubiquitinated lysine residues that directly regulate the stability of three transcription factors, CRYPTOCHROME-INTERACTING BASIC-HELIX-LOOP-HELIX 1 (CIB1), CIB1 LIKE PROTEIN 2 (CIL2), and SENSITIVE TO PROTON RHIZOTOXICITY1 (STOP1) using in vivo degradation assays. Furthermore, codon mutation of CIB1 to create a K166R conversion to prevent ubiquitination, via CRISPR/Cas9-derived adenosine base editing, led to an early flowering phenotype and increased expression of FLOWERING LOCUS T (FT). These comprehensive site-level ubiquitinome profiles provide a wealth of data for future functional studies related to modulation of biological processes mediated by this post-translational modification in plants.
The biomolecular networks underpinning cell function exhibit canalization, or the buffering of fluctuations required to function in a noisy environment. We present a new major release of $\tt{CANA}$, v1.0.0, an open-source Python package for understanding canalization in automata network models, discrete dynamical systems in which activation of biomolecular entities (e.g. transcription of genes) is modeled as the activity of coupled automata. One understudied putative mechanism for canalization is the functional equivalence of biomolecular regulators (e.g. among the transcription factors for a gene). We study this mechanism using the theory of symmetry in discrete functions. We present a new exact method, $\tt{schematodes}$, for finding maximal symmetry groups among the inputs to discrete functions, and integrate it into $\tt{CANA}$. The $\tt{schematodes}$ method substantially outperforms the inexact method of previous $\tt{CANA}$ versions both in speed and accuracy. We apply $\tt{CANA}$ v1.0.0 to study symmetry in 74 experimentally supported automata network models from the Cell Collective (CC) repository. The symmetry distribution is significantly different in the CC than in random automata with the same in-degree (connectivity) and bias (average output) (Kolmogorov–Smirnov test, P ≪ .001). Its spread is much wider than in a null model (IQR 0.31 versus IQR 0.20 with equal medians), demonstrating that the CC is enriched in functions with extreme symmetry or asymmetry.
Established research has illustrated that moderate exposure to stress in the womb influences both adult phonotype and genotype for several physiological pathways, especially in males. Proposed explanations include adaptions made by the fetus resulting from a limited supply of nutrients, referred to as the thrifty phenotype. In this study, we examine this fetal programming effect on the appetite control and energy expenditure pathways in prenatally stressed adult male offspring. Subjects were male rats born from time-mated female rats exposed to unpredictable, variable prenatal stress (UVPS) throughout gestation. An analysis of the adult male rat offspring genetic expression of epididymal fat pads and the plasma concentrations of hormones involved in appetite control and energy expenditure pathways showed a significantly diminished expression of leptin and adiponectin compared to unstressed controls. Leptin and adiponectin are both major hormones involved in the appetite control and energy expenditure pathways, with leptin regulating energy balance due to its function as an inhibitor of hunger, and adiponectin modulating glucose levels and fatty acid breakdown. We observed higher leptin concentrations within the prenatally stressed male plasma, and lower expression of leptin (OB) and adiponectin (ADIPOQ) genes from the epididymal fat pads. We suggest that elevated leptin in the plasma elicited a negative feedback effect on OB expression levels, decreasing their quantification compared to control animals. Further analysis will include plasma quantification of insulin and glucose, as well as expression of ghrelin, a peptide which acts on the central nervous system and the bodys perception of hunger.
Persistent drought affects global crop production and is becoming more severe in many parts of the world in recent decades. Deciphering how plants respond to drought will facilitate the development of flexible mitigation strategies. Sorghum bicolor L. Moench (sorghum), a major cereal crop and an emerging bioenergy crop, exhibits remarkable resilience to drought. To better understand the molecular traits that underlie sorghum's remarkable drought tolerance, we undertook a large-scale sorghum gene expression profiling effort, totaling nearly 1500 transcriptome profiles, across a 3-year field study with replicated plots in California's Central Valley. This study included time-resolved gene expression data from roots and leaves of two sorghum genotypes, BTx642 and RTx430, with different pre-flowering and post-flowering drought-tolerance adaptations under control and drought conditions. Quantification of genotype-specific drought tolerance effects was enabled by de novo sequencing, assembly, and annotation of both BTx642 and RTx430 genomes. These reference-quality genomes were used to construct a pangene set for characterizing conserved and genotype-specific expression. By integrating time-resolved transcriptomic responses to drought in the field across three consecutive years, we identified a set of 726 drought-responsive genes that responded similarly in all 3 years of our field study. Functional enrichment analysis identified abiotic stress, secondary cell wall-related processes and metabolism as particularly affected under both types of drought stress. We also found that some glyoxylate cycle pathway genes, including malate synthase and isocitrate lyase, are differentially regulated particularly during post-flowering drought stress, implicating this pathway as potentially important for drought responsiveness. This expansive dataset represents a unique resource for sorghum and drought research communities and provides a methodological framework for the integration of multi-faceted time-resolved transcriptomic datasets.
High-resolution imaging with secondary ion mass spectrometry (nanoSIMS) has become a standard method in systems biology and environmental biogeochemistry and is broadly used to decipher ecophysiological traits of environmental microorganisms, metabolic processes in plant and animal tissues, and cross-kingdom symbioses. When combined with stable isotope-labeling—an approach we refer to as nanoSIP—nanoSIMS imaging offers a distinctive means to quantify net assimilation rates and stoichiometry of individual cell-sized particles in both low- and high-complexity environments. While the majority of nanoSIP studies in environmental and microbial biology have focused on nitrogen and carbon metabolism (using 15 N and 13 C tracers), multiple advances have pushed the capabilities of this approach in the past decade. The development of a high-brightness oxygen ion source has enabled high-resolution metal analyses that are easier to perform, allowing quantification of metal distribution in cells and environmental particles. New preparation methods, tools for automated data extraction from large data sets, and analytical approaches that push the limits of sensitivity and spatial resolution have allowed for more robust characterization of populations ranging from marine archaea to fungi and viruses. Further, NanoSIMS studies continue to be enhanced by correlation with orthogonal imaging and ‘omics approaches; when linked to molecular visualization methods, such as in situ hybridization and antibody labeling, these techniques enable in situ function to be linked to microbial identity and gene expression. Here we present an updated description of the primary materials, methods, and calculations used for nanoSIP, with an emphasis on recent advances in nanoSIMS applications, key methodological steps, and potential pitfalls.
The evaluation of epigenetic features such as DNA methylation is becoming increasingly important in many biochemical processes like gene expression and transcription as well as in several diseases like schizophrenia or diabetes. Here, in this work, we report that self-assembled nanomechanical resonators entirely composed of DNA molecules can be used to explore gross changes in DNA methylation levels (0–25–50%), while careful control of tensile stress is needed to reduce the variability of resonance frequency for rigorous quantification. The effect of the tensile stress retained by the suspended DNA nanoresonators on the application of the technique is extensively explored using a combination of laser Doppler vibrometry and atomic force spectroscopy. DNA nanoresonators are real-time, label-free sensors and could avoid chemical functionalization and sample amplification. Therefore, they may represent in the future a key complementary routine tool for global DNA methylation analysis needed to evaluate the consequences of environmental stresses on the human genome.
OBJECTIVES AND METHODS This project addresses the challenge of understanding and predicting individual radiation sensitivity by integrating genetics, demographics and biomarker characteristics across species (mice and humans). We hypothesize that ex vivo DNA repair response to GCR components is a central determinant of cancer risk from space radiation and can serve as a biomarker of radiation risk in combination with genetics. Automated image quantification of 53BP1+ radiation-induced foci (RIF) during the first 4-48 h post-irradiation was performed as a function of dose and LET in non-immortalized primary skin fibroblasts derived from 76 mice across 15 strains (5 inbred reference strains and 10 collaborative-cross strains) exposed to X rays (0.1, 1 and 4 Gy), 350 MeV/n 40Ar and 600 MeV/n 56Fe (1.1 and 3 particles/100μm2), as well as in peripheral blood mononuclear cells (PBMCs) from 768 healthy donors (matched ethnicity, 50/50 male/female, 18-70 years old) exposed to gamma rays (0.1 and 1 Gy), 350 MeV/n 28Si, 350 MeV/n 40Ar and 600 MeV/n 56Fe (1.1 and 3 particles/100μm2). QUANTIFICATION OF 53BP1+ FOCI IN VITRO AND ASSOCIATIONS TO IN VIVO RADIATION SUSCEPTIBILITY IN 15 MOUSE STRAINS We reported in vitro repair kinetic and repairable fractions of RIF for the 15 mouse strains and introduced a mathematical model for RIF as a function of time, dose and LET. We noted that the metabolic activity of cells modulates the RIF response, and we introduced the open access tool terRIFic (Tool for Enhanced Results of RIF In Cells, https://radbiolab.shinyapps.io/terrific/) to correct for such bias using confluence level. Notably, at 4h post-irradiation, RIF/Gy decreased with dose or LET: as the dose or LET increases, so does the proximity of DNA double-strand-breaks (DSB) and our data suggest that proximal DSBs are brought together inside isolated RIF for repair. The RIF/Gy trend was inverted at 24h, suggesting RIF with high DSB content are more difficult to repair. We showed that in vitro metrics correlate with in vivo measurements in the same 15 mouse strains, such as survival levels of immune cells or spontaneous cancer incidence, suggesting a relationship between the efficiency of DSB repair and cancer risk or radiation toxicity. In addition to the efficiency of repair and persistent RIF, the amount of spontaneous foci before irradiation was also found to be strain dependent. Finally, we performed genome-wide association study in the same 15 mouse strains using all RIF phenotypes measured in vitro, identifying genes of interest and validating RIF as an ideal biomarker for individual radiation sensitivity. BASELINE 53BP1+ FOCI PREDICTS INDIVIDUAL HUMAN RESPONSE TO GCR COMPONENTS Based on the analysis of radiation responses of 576 donor PBMCs (using quantification of 53BP1+ foci, oxidative stress and cell death), we observed a wide variability of subject- and LET-dependent radiation responses, with radiation-induced DNA repair foci increasing with LET, though oxidative stress being notably reduced by high-LET irradiation, potentially due to a switch between hydrogen peroxide and oxygen radical-based mechanisms. We identified a relationship between few spontaneous DNA foci at baseline and increased DNA repair after irradiation, accompanied by an alteration in immunoregulatory cytokine secretion, which might be adapted as biomarkers to predict ionizing radiation sensitivity. Among demographic variables, only latent cytomegalovirus infection and age were predictive of high baseline foci formation. Finally, we have performed low-throughput whole genome sequencing of all samples and are currently in the process of identifying the genes and pathways associated with low and high-LET ionizing radiation sensitivity in humans.
Surface phytoplankton communities were linked with the carbon they export into the deep ocean by comparing 18 S rRNA gene sequence communities from surface seawater and individually isolated sinking particles. Particles were collected in sediment traps deployed at locations in the North Pacific subtropical gyre and the California Current. DNA was isolated from individual particles, bulk-collected trap particles, and the surface seawater. The relative sequence abundance of exported phytoplankton taxa in the surface water varied across functional groups and ecosystems. Of the sequences detected in sinking particles, about half were present in large (>300 μm), individually isolated particles and primarily belonged to taxa with small cell sizes (<50 μm). Exported phytoplankton taxa detected only in bulk trap samples, and thus presumably packaged in the smaller sinking size fraction, contained taxa that typically have large cell sizes (>500 m). The effect of particle degradation on the detectable 18 S rRNA gene community differed across taxa, and differences in community composition among individual particles from the same location largely reflected differences in relative degradation state. Using these data and particle imaging, we present an approach that incorporates genetic diversity into mechanistic models of the ocean's biological carbon pump, which will lead to better quantification of the ocean’s carbon cycle.
One-carbon (C 1 ) metabolism, centered on the universal methyl donor S-adenosyl methionine (SAM), plays critical roles in biosynthesis, redox regulation, and stress responses across plants and microbes. A recently proposed photosynthetic C 1 pathway links SAM methyl groups directly to RuBisCO-mediated CO 2 assimilation and integrates with nitrogen and sulfur metabolism. Light-dependent SAM synthesis may regulate the methylation of biopolymers and specialized metabolites and help mitigate photorespiratory stress under elevated temperature and drought. Phylogenetic analysis of two core enzymes suggests evolutionary continuity from methylotrophic microbes to land plants, supporting microbial origins via endosymbiotic gene transfer. Beyond intracellular roles, C 1 metabolism drives biosphere–atmosphere exchange via gases such as methane, methanol, formic acid, and formaldehyde, and numerous specialized volatiles synthesized through SAM methylation. S-methylmethionine, a mobile C 1 metabolite, may mediate phloem transport of reduced sulfur, nitrogen, and methyl groups, linking above- and belowground C 1 cycling in plants. Advances in real-time gas sensing now allow the high-frequency quantification of C 1 fluxes from leaves, stems, and soils, highlighting C 1 metabolism as a critical and underrecognized component of terrestrial carbon and nutrient cycling. Given its microbial ancestry and the production of diverse volatile biosignatures, C 1 metabolism may also offer unique insights into life's origins and biosignature detection on exoplanets.
This project addresses the challenge of understanding and predicting individual radiation sensitivity by integrating genetics, demographics and biomarker characteristics across species (mice and humans). We hypothesize that ex vivo DNA repair response to GCR components is a central determinant of cancer risk from space radiation and can serve as a biomarker of radiation risk in combination with genetics. Automated image quantification of 53BP1+ radiation-induced foci (RIF) during the first 4-48 h post-irradiation was performed as a function of dose and LET in non-immortalized primary skin fibroblasts derived from 76 mice across 15 strains (5 inbred reference strains and 10 collaborative-cross strains) exposed to X rays (0.1, 1 and 4 Gy), 350 MeV/n 40Ar and 600 MeV/n 56Fe (1.1 and 3 particles/100sq. μm), as well as in peripheral blood mononuclear cells (PBMCs) from 768 healthy donors (matched ethnicity, 50/50 male/female, 18-70 years old) exposed to gamma rays (0.1 and 1 Gy), 350 MeV/n 28Si, 350 MeV/n 40Ar and 600 MeV/n 56Fe (1.1 and 3 particles/100sq. μm). A genome-wide association study (GWAS) was performed on the mouse strains between DNA damage responses to space radiation and single nucleotide polymorphisms (SNPs). We found SNPs, which were significantly associated to the RIF phenotype, mapped to genes and pathways that are functionally linked to health hazards for deep space exploration (e.g. carcinogenesis, nervous system damage and immune dysfunction). Some of these SNPs were located within protein coding regions, potentially interfering with protein functions and providing promising genetic targets for countermeasures. We also found correlations between both spontaneous and radiation-induced DNA damage and SNPs mapped to pathways associated with cellular metabolism. GWAS is undergoing for the human data. All data have been made available via the NASA Space Biology Open-Science database (genelab.nasa.gov) and we will discuss how various genomic and transcriptomic datasets can be accessed for modeling and integrated using machine learning methods for discovering new radiation biology.
Phages are one of the key ecological drivers of microbial community dynamics, function, and evolution. Despite their importance in bacterial ecology and evolutionary processes, phage genes are poorly characterized, hampering their usage in a variety of biotechnological applications. Methods to characterize such genes, even those critical to the phage life cycle, are labor intensive and are generally phage specific. Here, we develop a systematic gene essentiality mapping method scalable to new phage–host combinations that facilitate the identification of nonessential genes. As a proof of concept, we use an arrayed genome-wide CRISPR interference (CRISPRi) assay to map gene essentiality landscape in the canonical coliphages λ and P1. Results from a single panel of CRISPRi probes largely recapitulate the essential gene roster determined from decades of genetic analysis for lambda and provide new insights into essential and nonessential loci in P1. We present evidence of how CRISPRi polarity can lead to false positive gene essentiality assignments and recommend caution towards interpreting CRISPRi data on gene essentiality when applied to less studied phages. Finally, we show that we can engineer phages by inserting DNA barcodes into newly identified inessential regions, which will empower processes of identification, quantification, and tracking of phages in diverse applications.
Aims/hypothesis Progression to type 1 diabetes is associated with genetic factors, the presence of autoantibodies and a decline in beta cell insulin secretion in response to glucose. Very little is known regarding the molecular changes that occur in human insulin-secreting beta cells prior to the onset of type 1 diabetes. Herein, we applied an unbiased proteomics approach to identify changes in proteins and potential mechanisms of islet dysfunction in islet-autoantibody-positive organ donors with pre-symptomatic stage 1 type 1 diabetes (HbA1c ≤42 mmol/mol [6.0%]). We aimed to identify pathways in islets that are indicative of beta cell dysfunction. Methods Multiple islet sections were collected through laser microdissection of frozen pancreatic tissues from organ donors positive for single or multiple islet autoantibodies (AAb + , n=5), and age (±2 years)- and sex-matched non-diabetic (ND) control donors (n=5) obtained from the Network for Pancreatic Organ donors with Diabetes (nPOD). Islet sections were subjected to MS-based proteomics and analysed with label-free quantification followed by pathway and functional annotations. Results Analyses resulted in ~4500 proteins identified with low false discovery rate (<1%), with 2165 proteins reliably quantified in every islet sample. We observed large inter-donor variations that presented a challenge for statistical analysis of proteome changes between donor groups. We therefore focused on only the donors with stage 1 type 1 diabetes who were positive for multiple autoantibodies (mAAb + , n=3) and genetic risk compared with their matched ND controls (n=3) for the final statistical analysis. Approximately 10% of the proteins (n=202) were significantly different (unadjusted p<0.025, q<0.15) for mAAb + vs ND donor islets. The significant alterations clustered around major functions for upregulation in the immune response and glycolysis, and downregulation in endoplasmic reticulum (ER) stress response as well as protein translation and synthesis. The observed proteome changes were further supported by several independent published datasets, including a proteomics dataset from in vitro proinflammatory cytokine-treated human islets and single-cell RNA-seq datasets from AAb + individuals. Conclusions/interpretation In situ human islet proteome alterations in stage 1 type 1 diabetes centred around several major functional categories, including an expected increase in immune response genes (elevated antigen presentation/HLA), with decreases in protein synthesis and ER stress response, as well as compensatory metabolic response. The dataset serves as a proteomics resource for future studies on beta cell changes during type 1 diabetes progression and pathogenesis. Data availability The LC-MS raw datasets that support the findings of this study have been deposited in the online repository: MassIVE (https://massive.ucsd.edu/ProteoSAFe/static/massive.jsp) with accession no. MSV000090212.
Abstract Background The endocannabinoid system (ECS) is highly integrated with seemingly all physiological and pathophysiological processes in the body. There is increasing interest in utilizing bioactive plant compounds, for promoting health and improving production in livestock. Given the established interaction between phytochemicals and the ECS, there are many opportunities for identification and development of therapies to address a range of diseases and disorders. However, the ECS has not been thoroughly characterized in cattle, especially in the gastrointestinal tract. The objective of this study was to characterize the distribution and transcriptional abundance of genes associated with the endocannabinoid system in bovine tissues. Methods Tissues including brain, spleen, thyroid, lung, liver, kidney, mesenteric vein, tongue, sublingual mucosa, rumen, omasum, duodenum, jejunum, ileum and colon were collected from 10-mo old Holstein steers (n = 6). Total RNA was extracted and gene expression was measured using absolute quantification real time qPCR. Gene expression of endocannabinoid receptorsCNR1andCNR2, synthesis enzymesDAGLA,DAGLBandNAPEPLD, degradation enzymesMGLLandFAAH, and transient receptor potential vanilloidsTRPV3andTRPV6was measured. Data were analyzed in R using a Kruskal-Wallis followed by a Wilcoxon rank-sum test. Results are reported as the median copy number/20 ng of equivalent cDNA (CN) with interquartile range (IQR). Results The greatest expression ofCNR1andCNR2was in the brain and spleen, respectively. Expression of either receptor was not detected in any gastrointestinal tissues, however there was a tendency (P = 0.095) forCNR2to be expressed above background in rumen. Expression of endocannabinoid synthesis and degradation enzymes varied greatly across tissues. Brain tissue had the greatestDAGLAexpression at 641 CN (IQR 52;P ≤ 0.05).DAGLBwas detected in all tissues, with brain and spleen having the greatest expression (P ≤ 0.05). Expression ofNAPEPLDin the gastrointestinal tract was lowest in tongue and sublingual mucosal. There was no difference in expression ofNAPEPLDbetween hindgut tissues, however these tissues collectively had 592% greater expression than rumen and omasum (P ≤ 0.05). WhileMGLLwas found to be expressed in all tissues, expression ofFAAHwas only above the limit of detection in brain, liver, kidney, jejunum and ileum.TRPV3was expressed above background in tongue, rumen, omasum and colon. Although not different from each other, thyroid and duodenum had the greatest expression ofTRPV6, with 285 (IQR 164) and 563 (IQR 467) CN compared to all other tissues (P < 0.05). Conclusions These data demonstrate the complex distribution and variation of the ECS in bovine tissues. Expression patterns suggest that regulatory functions of this system are tissue dependent, providing initial insight into potential target tissues for manipulation of the ECS.