CAIM: Coverage-based Analysis for Identification of Microbiome.

Accurate taxonomic profiling of microbial taxa in a metagenomic sample is vital to gain insights into microbial ecology. Recent advancements in sequencing technologies have contributed tremendously toward understanding these microbes at species resolution through a whole shotgun metagenomic (WMS) approach. In this study, we developed a new bioinformatics tool, CAIM, for accurate taxonomic classification and quantification within both long- and short-read metagenomic samples using an alignment-based method. CAIM depends on two different containment techniques to identify species in metagenomic samples using their genome coverage information to filter out false positives rather than the traditional approach of relative abundance. In addition, we propose a nucleotide-count based abundance estimation, which yield lesser root mean square error than the traditional read-count approach. We evaluated the performance of CAIM on 28 metagenomic mock communities and 2 synthetic datasets by comparing it with other top-performing tools. CAIM maintained a consitently good performance across datasets in identifying microbial taxa and in estimating relative abundances than other tools. CAIM was then applied to a real dataset sequenced on both Nanopore (with and without amplification) and Illumina sequencing platforms and found high similality of taxonomic profiles between the sequencing platforms. Lastly, CAIM was applied to fecal shotgun metagenomic datasets of 232 colorectal cancer patients and 229 controls obtained from 4 different countries and primary 44 liver cancer patients and 76 controls. The predictive performance of models using the genome-coverage cutoff was better than those using the relative-abundance cutoffs in discriminating colorectal cancer and primary liver cancer patients from healthy controls with a highly confident species markers.

Short-Term Ingestion of Essential Amino Acid Based Nutritional Supplements or Whey Protein Improves the Physical Function of Older Adults Independently of Gut Microbiome.

SCOPE: Dietary proteins and essential amino acids (EAAs) are the major nutritional supplements that support the growth and activity of gut microbes contributing to the wellbeing of their host. This study hypothesizes that daily supplementation of the diet with either EAAs or whey protein for 12 weeks would improve the gut microbiome of older adults. METHODS AND RESULTS: The stool samples are processed and subjected to Illumina-based 16S ribosomal ribonucleic acid (rRNA) gene amplicon sequencing. In both groups, the most abundant families are found in order of relative abundance included: Bacteroidaceae, Lachnospiraceae, Ruminococcaceae, Prevotellaceae, Rikenellaceae, Enterobacteriaceae, Oscillospiraceae, Tannerellaceae, and Akkermansiaceae, which indicate that these subjects are able to maintain a same healthy microbial diversity in their guts. A significant finding is a reduction of proinflammatory cytokine, interleukin-18 (IL-18) in the EAAs group. It also uses the standard 6-min walking test (6MWT) as a measure of cardiopulmonary fitness. At the end of the study, the subjects in the EAAs group perform significantly better in the 6MWT as compared to the whey group. CONCLUSION: It seems plausible that the improved physical performance and reduced proinflammatory cytokine, IL-18 seen in the EAAs group, are independent of changes in gut microbiota.

Refining the identity of mesenchymal cell types associated with murine periosteal and endosteal bone.

Single-cell RNA-seq has led to novel designations for mesenchymal cells associated with bone as well as multiple designations for what appear to be the same cell type. The main goals of this study were to increase the amount of single-cell RNA sequence data for osteoblasts and osteocytes, to compare cells from the periosteum to those inside bone, and to clarify the major categories of cell types associated with murine bone. We created an atlas of murine bone-associated cells by harmonizing published datasets with in-house data from cells targeted by Osx1-Cre and Dmp1-Cre driver strains. Cells from periosteal bone were analyzed separately from those isolated from the endosteum and trabecular bone. Over 100,000 mesenchymal cells were mapped to reveal 11 major clusters designated fibro-1, fibro-2, chondrocytes, articular chondrocytes, tenocytes, adipo-Cxcl12 abundant reticular (CAR), osteo-CAR, preosteoblasts, osteoblasts, osteocytes, and osteo-X, the latter defined in part by periostin expression. Osteo-X, osteo-CAR, and preosteoblasts were closely associated with osteoblasts at the trabecular bone surface. Wnt16 was expressed in multiple cell types from the periosteum but not in cells from endocortical or cancellous bone. Fibro-2 cells, which express markers of stem cells, localized to the periosteum but not trabecular bone in adult mice. Suppressing bone remodeling eliminated osteoblasts and altered gene expression in preosteoblasts but did not change the abundance or location of osteo-X or osteo-CAR cells. These results provide a framework for identifying bone cell types in murine single-cell RNA-seq datasets and suggest that osteoblast progenitors reside near the surface of remodeling bone.

A NOTCH3-CXCL12-driven myeloma-tumor niche signaling axis promotes chemoresistance in multiple myeloma.

Multiple myeloma (MM) remains incurable due to disease relapse and drug resistance. Notch signals from the tumor microenvironment (TME) confer chemoresistance, but the cellular and molecular mechanisms are not entirely understood. Using clinical and transcriptomic datasets, we found that NOTCH3 is upregulated in CD138+ cells from newly diagnosed MM (NDMM) patients compared to healthy individuals and increased in progression/relapsed MM (PRMM) patients. Further, NDMM patients with high NOTCH3 expression exhibited worse responses to Bortezomib (BOR)-based therapies. Cells of the TME, including osteocytes, upregulated NOTCH3 in MM cells and protected them from apoptosis induced by BOR. NOTCH3 activation (NOTCH3OE) in MM cells decreased BOR anti-MM efficacy and its ability to improve survival in in vivo myeloma models. Molecular analyses revealed that NDMM and PRMM patients with high NOTCH3 exhibit CXCL12 upregulation. TME cells upregulated CXCL12 and activated the CXCR4 pathway in MM cells in a NOTCH3-dependent manner. Moreover, genetic or pharmacologic inhibition of CXCL12 in NOTCH3OE MM cells restored sensitivity to BOR regimes in vitro and in human bones bearing NOTCH3OE MM tumors cultured ex vivo. Our clinical and preclinical data unravel a novel NOTCH3-CXCL12 pro-survival signaling axis in the TME and suggest that osteocytes transmit chemoresistance signals to MM cells.

Exploiting nanopore sequencing for characterization and grading of IDH-mutant gliomas.

The 2021 WHO Classification of Central Nervous System Tumors recommended evaluation of cyclin-dependent kinase inhibitor 2A/B (CDKN2A/B) deletion in addition to codeletion of 1p/19q to characterize IDH-mutant gliomas. Here, we demonstrated the use of a nanopore-based copy-number variation sequencing (nCNV-seq) approach to simultaneously identify deletions of CDKN2A/B and 1p/19q. The nCNV-seq approach was initially evaluated on three distinct glioma cell lines and then applied to 19 IDH-mutant gliomas (8 astrocytomas and 11 oligodendrogliomas) from patients. The whole-arm 1p/19q codeletion was detected in all oligodendrogliomas with high concordance among nCNV-seq, FISH, DNA methylation profiling, and whole-genome sequencing. For the CDKN2A/B deletion, nCNV-seq detected the loss in both astrocytoma and oligodendroglioma, with strong correlation with the CNV profiles derived from whole-genome sequencing (Pearson correlation (r) = 0.95, P < 2.2 × 10-16 to r = 0.99, P < 2.2 × 10-16 ) and methylome profiling. Furthermore, nCNV-seq can differentiate between homozygous and hemizygous deletions of CDKN2A/B. Taken together, nCNV-seq holds promise as a new, alternative approach for a rapid and simultaneous detection of the molecular signatures of IDH-mutant gliomas without capital expenditure for a sequencer.

A framework for defining mesenchymal cell types associated with murine periosteal and endosteal bone.

Single-cell RNA sequencing has led to numerous novel designations for mesenchymal cell types associated with bone. Consequently, there are now multiple designations for what appear to be the same cell type. In addition, existing datasets contain relatively small numbers of mature osteoblasts and osteocytes and there has been no comparison of periosteal bone cells to those at the endosteum and trabecular bone. The main goals of this study were to increase the amount of single cell RNA sequence data for osteoblasts and osteocytes, to compare cells from the periosteum to those inside bone, and to clarify the major categories of cell types associated with murine bone. To do this, we created an atlas of murine bone-associated cells by harmonizing published datasets with in-house data from cells targeted by Osx1-Cre and Dmp1-Cre driver strains. Cells from periosteal bone were analyzed separately from those isolated from the endosteum and trabecular bone. Over 100,000 mesenchymal cells were mapped to reveal 11 major clusters designated fibro-1, fibro-2, chondrocytes, articular chondrocytes, tenocytes, adipo-CAR, osteo-CAR, pre-osteoblasts, osteoblasts, osteocytes, and osteo-X, the latter defined in part by Postn expression. Osteo-X, osteo-CAR, and pre-osteoblasts were closely associated with osteoblasts at the trabecular bone surface. Wnt16 was expressed in multiple cell types from the periosteum but not in any cells from endocortical or cancellous bone. Fibro-2 cells, which express markers of skeletal stem cells, localized to the periosteum but not trabecular bone in adult mice. Suppressing bone remodeling eliminated osteoblasts and altered gene expression in pre-osteoblasts but did not change the abundance or location of osteo-X or osteo-CAR cells. These results provide a framework for identifying bone cell types in murine single cell RNA sequencing datasets and suggest that osteoblast progenitors reside near the surface of remodeling bone.

High NEK2 expression in myeloid progenitors suppresses T cell immunity in multiple myeloma.

Multiple myeloma (MM) growth is supported by an immune-tolerant bone marrow microenvironment. Here, we find that loss of Never in mitosis gene A (NIMA)-related kinase 2 (NEK2) in tumor microenvironmental cells is associated with MM growth suppression. The absence of NEK2 leads to both fewer tumor-associated macrophages (TAMs) and inhibitory T cells. NEK2 expression in myeloid progenitor cells promotes the generation of functional TAMs when stimulated with MM conditional medium. Clinically, high NEK2 expression in MM cells is associated with increased CD8+ T effector memory cells, while low NEK2 is associated with an IFN-γ gene signature and activated T cell response. Inhibition of NEK2 upregulates PD-L1 expression in MM cells and myeloid cells. In a mouse model, the combination of NEK2 inhibitor INH154 with PD-L1 blockade effectively eliminates MM cells and prolongs survival. Our results provide strong evidence that NEK2 inhibition may overcome tumor immune escape and support its further clinical development.

Long-term benefit of DAAs on gut dysbiosis and microbial translocation in HCV-infected patients with and without HIV coinfection.

Long-term effect of Direct-acting antivirals (DAAs) on gut microbiota, short-chain fatty acids (SCFAs) and microbial translocation in patients with hepatitis C virus (HCV) infection who achieve sustained virological response (SVR) were limited. A longitudinal study of 50 patients with HCV monoinfection and 19 patients with HCV/HIV coinfection received DAAs were conducted. Fecal specimens collected at baseline and at week 72 after treatment completion (FUw72) were analyzed for 16S rRNA sequencing and the butyryl-CoA:acetateCoA transferase (BCoAT) gene expression using real-time PCR. Plasma lipopolysaccharide binding protein (LBP) and intestinal fatty acid binding protein (I-FABP) were quantified by ELISA assays. SVR rates in mono- and coinfected patients were comparable (94% vs. 100%). The improvement of gut dysbiosis and microbial translocation was found in responders but was not in non-responders. Among responders, significant restoration of alpha-diversity, BCoAT and LBP were observed in HCV patients with low-grade fibrosis (F0-F1), while HCV/HIV patients exhibited partial improvement at FUw72. I-FABP did not decline significantly in responders. Treatment induced microbiota changes with increasing abundance of SCFAs-producing bacteria, including Blautia, Fusicatenibacter, Subdoligranulum and Bifidobacterium. In conclusion, long-term effect of DAAs impacted the restoration of gut dysbiosis and microbial translocation. However, early initiation of DAAs required for an alteration of gut microbiota, enhanced SCFAs-producing bacteria, and could reduce HCV-related complications.

Reduced osteoprotegerin expression by osteocytes may contribute to rebound resorption after denosumab discontinuation.

Denosumab is an anti-RANKL Ab that potently suppresses bone resorption, increases bone mass, and reduces fracture risk. Discontinuation of denosumab causes rapid rebound bone resorption and bone loss, but the molecular mechanisms are unclear. We generated humanized RANKL mice and treated them with denosumab to examine the cellular and molecular conditions associated with rebound resorption. Denosumab potently suppressed both osteoclast and osteoblast numbers in cancellous bone in humanized RANKL mice. The decrease in osteoclast number was not associated with changes in osteoclast progenitors in bone marrow. Long-term, but not short-term, denosumab administration reduced osteoprotegerin (OPG) mRNA in bone. Localization of OPG expression revealed that OPG mRNA is produced by a subpopulation of osteocytes. Long-term denosumab administration reduced osteocyte OPG mRNA, suggesting that OPG expression declines as osteocytes age. Consistent with this, osteocyte expression of OPG was more prevalent near the surface of cortical bone in humans and mice. These results suggest that new osteocytes are an important source of OPG in remodeling bone and that suppression of remodeling reduces OPG abundance by reducing new osteocyte formation. The lack of new osteocytes and the OPG they produce may contribute to rebound resorption after denosumab discontinuation.

CRISPR interference provides increased cell type-specificity compared to the Cre-loxP system.

Cre-mediated recombination is frequently used for cell type-specific loss of function (LOF) studies. A major limitation of this system is recombination in unwanted cell types. CRISPR interference (CRISPRi) has been used effectively for global LOF in mice. However, cell type-specific CRISPRi, independent of recombination-based systems, has not been reported. To test the feasibility of cell type-specific CRISPRi, we produced two novel knock-in mouse models that achieve gene suppression when used together: one expressing dCas9::KRAB under the control of a cell type-specific promoter and the other expressing a single guide RNA from a safe harbor locus. We then compared the phenotypes of mice in which the same gene was targeted by either CRISPRi or the Cre-loxP system, with cell specificity conferred by Dmp1 regulatory elements in both cases. We demonstrate that CRISPRi is effective for cell type-specific LOF and that it provides improved cell type-specificity compared to the Cre-loxP system.

Growth and dormancy control of myeloma cells by mesenchymal stem cells.

Bone marrow mesenchymal stem cells (MSCs) may have contrasting impacts on the progression of multiple myeloma (MM). Priming normal MSCs, by culturing them with MM cells, mimics the MSC-induced MM growth. We studied the contrasting effects of conditioned medium (CM) from unprimed or primed MSCs on growth of MM cells from newly diagnosed cases. We elucidated potential molecular pathways using global gene expression profiling and focused on the role of the mTOR2 component, RICTOR, as a novel mediator of dormancy in MM. Primed MSCs CM consistently increased proportions of proliferating cells and supported MM growth in 3-day (n = 20) and 10-day (n = 12) cultures, effects that were partially mediated through the IGF1 axis. In contrast, unprimed MSCs CM inhibited growth of MM cells in cases mainly from stages I/II MM. The genes most overexpressed in MM cells treated with primed MSCs CM were associated with cell cycle, DNA-damage repair, and proliferation; genes most overexpressed in MM cells treated with unprimed MSCs CM were associated with dormancy pathways including RICTOR (mTOR2 pathway), CXCR4, and BCL2. RICTOR protein level was induced by unprimed MSCs CM and was lower in KI67+ proliferating MM cells treated with primed MSCs CM. RICTOR was underexpressed in clinical relapse samples compared with baseline samples of the same patients. Inhibiting RICTOR expression in primary MM cells promoted their growth, and enforced expression of RICTOR in MM cell lines inhibited their growth. Our findings suggest that, after prolonged interactions with MM cells, bone marrow MSCs shift from MM-repressive to MM-permissive. AVAILABILITY OF DATA AND MATERIALS: Our institutional GEP data of MM cells from newly diagnosed patients used to show RICTOR expression have been deposited at Gene Expression Omnibus (GEO: GSE2658, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE2658).

The NAD salvage pathway in mesenchymal cells is indispensable for skeletal development in mice.

NAD is an essential co-factor for cellular energy metabolism and multiple other processes. Systemic NAD+ deficiency has been implicated in skeletal deformities during development in both humans and mice. NAD levels are maintained by multiple synthetic pathways but which ones are important in bone forming cells is unknown. Here, we generate mice with deletion of Nicotinamide Phosphoribosyltransferase (Nampt), a critical enzyme in the NAD salvage pathway, in all mesenchymal lineage cells of the limbs. At birth, NamptΔPrx1 exhibit dramatic limb shortening due to death of growth plate chondrocytes. Administration of the NAD precursor nicotinamide riboside during pregnancy prevents the majority of in utero defects. Depletion of NAD post-birth also promotes chondrocyte death, preventing further endochondral ossification and joint development. In contrast, osteoblast formation still occurs in knockout mice, in line with distinctly different microenvironments and reliance on redox reactions between chondrocytes and osteoblasts. These findings define a critical role for cell-autonomous NAD homeostasis during endochondral bone formation.

Reversal of the diabetic bone signature with anabolic therapies in mice.

The mechanisms underlying the bone disease induced by diabetes are complex and not fully understood; and antiresorptive agents, the current standard of care, do not restore the weakened bone architecture. Herein, we reveal the diabetic bone signature in mice at the tissue, cell, and transcriptome levels and demonstrate that three FDA-approved bone-anabolic agents correct it. Diabetes decreased bone mineral density (BMD) and bone formation, damaged microarchitecture, increased porosity of cortical bone, and compromised bone strength. Teriparatide (PTH), abaloparatide (ABL), and romosozumab/anti-sclerostin antibody (Scl-Ab) all restored BMD and corrected the deteriorated bone architecture. Mechanistically, PTH and more potently ABL induced similar responses at the tissue and gene signature levels, increasing both formation and resorption with positive balance towards bone gain. In contrast, Scl-Ab increased formation but decreased resorption. All agents restored bone architecture, corrected cortical porosity, and improved mechanical properties of diabetic bone; and ABL and Scl-Ab increased toughness, a fracture resistance index. Remarkably, all agents increased bone strength over the healthy controls even in the presence of severe hyperglycemia. These findings demonstrate the therapeutic value of bone anabolic agents to treat diabetes-induced bone disease and suggest the need for revisiting the approaches for the treatment of bone fragility in diabetes.

Dissecting Phenotype from Genotype with Clinical Isolates of SARS-CoV-2 First Wave Variants.

The emergence and availability of closely related clinical isolates of SARS-CoV-2 offers a unique opportunity to identify novel nonsynonymous mutations that may impact phenotype. Global sequencing efforts show that SARS-CoV-2 variants have emerged and then been replaced since the beginning of the pandemic, yet we have limited information regarding the breadth of variant-specific host responses. Using primary cell cultures and the K18-hACE2 mouse, we investigated the replication, innate immune response, and pathology of closely related, clinical variants circulating during the first wave of the pandemic. Mathematical modeling of the lung viral replication of four clinical isolates showed a dichotomy between two B.1. isolates with significantly faster and slower infected cell clearance rates, respectively. While isolates induced several common immune host responses to infection, one B.1 isolate was unique in the promotion of eosinophil-associated proteins IL-5 and CCL11. Moreover, its mortality rate was significantly slower. Lung microscopic histopathology suggested further phenotypic divergence among the five isolates showing three distinct sets of phenotypes: (i) consolidation, alveolar hemorrhage, and inflammation, (ii) interstitial inflammation/septal thickening and peribronchiolar/perivascular lymphoid cells, and (iii) consolidation, alveolar involvement, and endothelial hypertrophy/margination. Together these findings show divergence in the phenotypic outcomes of these clinical isolates and reveal the potential importance of nonsynonymous mutations in nsp2 and ORF8.

Single cell cortical bone transcriptomics define novel osteolineage gene sets altered in chronic kidney disease.

Introduction: Due to a lack of spatial-temporal resolution at the single cell level, the etiologies of the bone dysfunction caused by diseases such as normal aging, osteoporosis, and the metabolic bone disease associated with chronic kidney disease (CKD) remain largely unknown. Methods: To this end, flow cytometry and scRNAseq were performed on long bone cells from Sost-cre/Ai9+ mice, and pure osteolineage transcriptomes were identified, including novel osteocyte-specific gene sets. Results: Clustering analysis isolated osteoblast precursors that expressed Tnc, Mmp13, and Spp1, and a mature osteoblast population defined by Smpd3, Col1a1, and Col11a1. Osteocytes were demarcated by Cd109, Ptprz1, Ramp1, Bambi, Adamts14, Spns2, Bmp2, WasI, and Phex. We validated our in vivo scRNAseq using integrative in vitro promoter occupancy via ATACseq coupled with transcriptomic analyses of a conditional, temporally differentiated MSC cell line. Further, trajectory analyses predicted osteoblast-to-osteocyte transitions via defined pathways associated with a distinct metabolic shift as determined by single-cell flux estimation analysis (scFEA). Using the adenine mouse model of CKD, at a time point prior to major skeletal alterations, we found that gene expression within all stages of the osteolineage was disturbed. Conclusion: In sum, distinct populations of osteoblasts/osteocytes were defined at the single cell level. Using this roadmap of gene assembly, we demonstrated unrealized molecular defects across multiple bone cell populations in a mouse model of CKD, and our collective results suggest a potentially earlier and more broad bone pathology in this disease than previously recognized.

Revisiting chloroplast genomic landscape and annotation towards comparative chloroplast genomes of Rhamnaceae.

BACKGROUND: Massive parallel sequencing technologies have enabled the elucidation of plant phylogenetic relationships from chloroplast genomes at a high pace. These include members of the family Rhamnaceae. The current Rhamnaceae phylogenetic tree is from 13 out of 24 Rhamnaceae chloroplast genomes, and only one chloroplast genome of the genus Ventilago is available. Hence, the phylogenetic relationships in Rhamnaceae remain incomplete, and more representative species are needed. RESULTS: The complete chloroplast genome of Ventilago harmandiana Pierre was outlined using a hybrid assembly of long- and short-read technologies. The accuracy and validity of the final genome were confirmed with PCR amplifications and investigation of coverage depth. Sanger sequencing was used to correct for differences in lengths and nucleotide bases between inverted repeats because of the homopolymers. The phylogenetic trees reconstructed using prevalent methods for phylogenetic inference were topologically similar. The clustering based on codon usage was congruent with the molecular phylogenetic tree. The groups of genera in each tribe were in accordance with tribal classification based on molecular markers. We resolved the phylogenetic relationships among six Hovenia species, three Rhamnus species, and two Ventilago species. Our reconstructed tree provides the most complete and reliable low-level taxonomy to date for the family Rhamnaceae. Similar to other higher plants, the RNA editing mostly resulted in converting serine to leucine. Besides, most genes were subjected to purifying selection. Annotation anomalies, including indel calling errors, unaligned open reading frames of the same gene, inconsistent prediction of intergenic regions, and misannotated genes, were identified in the published chloroplast genomes used in this study. These could be a result of the usual imperfections in computational tools, and/or existing errors in reference genomes. Importantly, these are points of concern with regards to utilizing published chloroplast genomes for comparative genomic analysis. CONCLUSIONS: In summary, we successfully demonstrated the use of comprehensive genomic data, including DNA and amino acid sequences, to build a reliable and high-resolution phylogenetic tree for the family Rhamnaceae. Additionally, our study indicates that the revision of genome annotation before comparative genomic analyses is necessary to prevent the propagation of errors and complications in downstream analysis and interpretation.

CReSIL: accurate identification of extrachromosomal circular DNA from long-read sequences.

Extrachromosomal circular DNA (eccDNA) of chromosomal origin is found in many eukaryotic species and cell types, including cancer, where eccDNAs with oncogenes drive tumorigenesis. Most studies of eccDNA employ short-read sequencing for their identification. However, short-read sequencing cannot resolve the complexity of genomic repeats, which can lead to missing eccDNA products. Long-read sequencing technologies provide an alternative to constructing complete eccDNA maps. We present a software suite, Construction-based Rolling-circle-amplification for eccDNA Sequence Identification and Location (CReSIL), to identify and characterize eccDNA from long-read sequences. CReSIL's performance in identifying eccDNA, with a minimum F1 score of 0.98, is superior to the other bioinformatic tools based on simulated data. CReSIL provides many useful features for genomic annotation, which can be used to infer eccDNA function and Circos visualization for eccDNA architecture investigation. We demonstrated CReSIL's capability in several long-read sequencing datasets, including datasets enriched for eccDNA and whole genome datasets from cells containing large eccDNA products. In conclusion, the CReSIL suite software is a versatile tool for investigating complex and simple eccDNA in eukaryotic cells.

Piezo1 stimulates mitochondrial function via cAMP signaling.

Mechanical signals stimulate mitochondrial function but the molecular mechanisms are not clear. Here, we show that the mechanically sensitive ion channel Piezo1 plays a critical role in mitochondrial adaptation to mechanical stimulation. The activation of Piezo1 induced mitochondrial calcium uptake and oxidative phosphorylation (OXPHOS). In contrast, loss of Piezo1 reduced the mitochondrial oxygen consumption rate (OCR) and adenosine triphosphate (ATP) production in calvarial cells and these changes were associated with increased expression of the phosphodiesterases Pde4a and lower cyclic AMP (cAMP) levels. In addition, Piezo1 increased cAMP production and the activation of a cAMP-responsive transcriptional reporter. Consistent with this, cAMP was sufficient to increase mitochondrial OCR and the inhibition of phosphodiesterases augmented the increase in OCR induced by Piezo1. Moreover, the inhibition of cAMP production or activity of protein kinase A, a kinase activated by cAMP, prevented the increase in OCR induced by Piezo1. These results demonstrate that cAMP signaling contributes to the increase in mitochondrial OXPHOS induced by activation of Piezo1.