Search

A modified Chromium 10x droplet-based protocol that subsamples cells for both short-read and long-read (nanopore) sequencing together with a new computational pipeline (FLAMES) is developed to enable isoform discovery, splicing analysis, and mutation detection in single cells. We identify thousands of unannotated isoforms and find conserved functional modules that are enriched for alternative transcript usage in different cell types and species, including ribosome biogenesis and mRNA splicing. Analysis at the transcript level allows data integration with scATAC-seq on individual promoters, improved correlation with protein expression data, and linked mutations known to confer drug resistance to transcriptome heterogeneity. Supplementary Information The online version contains supplementary material available at 10.1186/s13059-021-02525-6.

Impaired oxidative capacity and mitochondrial function contribute to the dystrophic pathology in muscles of patients with Duchenne muscular dystrophy (DMD) and in relevant mouse models of the disease. Emerging evidence suggests an association between disrupted core clock expression and mitochondrial quality control, but this has not been established in muscles lacking dystrophin. We examined the diurnal regulation of muscle core clock and mitochondrial quality control expression in dystrophin-deficient C57BL/10ScSn- Dmdmdx ( mdx) mice, an established model of DMD. Male C57BL/10 (BL/10; n = 18) and mdx mice ( n = 18) were examined every 4 h beginning at the dark cycle. Throughout the entire light-dark cycle, extensor digitorum longus (EDL) muscles from mdx mice had decreased core clock mRNA expression ( Arntl, Cry1, Cry2, Nr1d2; P < 0.05) and disrupted mitochondrial quality control mRNA expression related to biogenesis (decreased; Ppargc1a, Esrra; P < 0.05), fission (increased; Dnm1l; P < 0.01), fusion (decreased; Opa1, Mfn1; P < 0.05), and autophagy/mitophagy (decreased: Bnip3; P < 0.05; increased: Becn1; P < 0.05). Cosinor analysis revealed a decrease in the rhythmicity parameters mesor and amplitude for Arntl, Cry1, Cry2, Per2, and Nr1d1 ( P < 0.001) in mdx mice. Diurnal oscillations in Esrra, Sirt1, Map1lc3b, and Sqstm1 were absent in mdx mice, along with decreased mesor and amplitude of Ppargc1a mRNA expression ( P < 0.01). The expression of proteins involved in mitochondrial biogenesis (decreased: PPARGC1A, P < 0.05) and autophagy/mitophagy (increased: MAP1LC3BII, SQSTM1, BNIP3; P < 0.05) were also dysregulated in tibialis anterior muscles of mdx mice. These findings suggest that dystrophin deficiency in mdx mice impairs the regulation of the core clock and mitochondrial quality control, with relevance to DMD and related disorders.

The development of the musculoskeletal system and its maintenance depends on the reciprocal relationship between muscle and bone. The size of skeletal muscles and the forces generated during muscle contraction are potent sources of mechanical stress on the developing skeleton, and they shape bone structure during growth. This is particularly evident in hypermuscular global myostatin (Mstn)‐null mice, where larger muscles during development increase bone mass and alter bone shape. However, whether muscle hypertrophy can similarly influence the shape of bones after the embryonic and prepubertal period is unknown. To address this issue, bone structure was assessed after inducing muscle hypertrophy in the lower hindlimbs of young‐adult C57BL/6J male mice by administering intramuscular injections of recombinant adeno‐associated viral vectors expressing follistatin (FST), a potent antagonist of Mstn. Two FST isoforms were used: the full‐length 315 amino acid isoform (FST‐315) and a truncated 288 amino acid isoform (FST‐288). In both FST‐treated cohorts, muscle hypertrophy was observed, and the anterior crest of the tibia, adjacent to the tibialis anterior muscle, was lengthened. Hypertrophy of the muscles surrounding the tibia caused the adjacent cortical shell to recede inward toward the central axis: an event driven by bone resorption adjacent to the hypertrophic muscle. The findings reveal that inducing muscle hypertrophy in mice can confer changes in bone shape in early adulthood. © 2021 The Authors. JBMR Plus published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research.

Bone strength is controlled by both bone mass, and the organization and quality of the bone material. The current standard method for measuring bone mass in mouse and rat studies is micro-computed tomography. This method typically uses a single threshold to identify bone material in the cortical and trabecular regions. However, this single threshold method obscures information about the mineral content of the bone material and depends on normal morphology to separately analyze cortical and trabecular structures. To extend this method to identify bone mass at multiple density levels, we have established a protocol for unbiased selection and application of multiple thresholds using a standard laboratory-based micro-computed tomography instrument. This non-invasive method can be applied to longitudinal studies and archived samples and provides additional information about bone structure and strength.

Alternative splicing shapes the phenotype of cells in development and disease. Long-read RNA-sequencing recovers full-length transcripts but has limited throughput at the single-cell level. Here we developed single-cell full-length transcript sequencing by sampling (FLT-seq), together with the computational pipeline FLAMES to overcome these issues and perform isoform discovery and quantification, splicing analysis and mutation detection in single cells. With FLT-seq and FLAMES, we performed the first comprehensive characterization of the full-length isoform landscape in single cells of different types and species and identified thousands of unannotated isoforms. We found conserved functional modules that were enriched for alternative transcript usage in different cell populations, including ribosome biogenesis and mRNA splicing. Analysis at the transcript-level allowed data integration with scATAC-seq on individual promoters, improved correlation with protein expression data and linked mutations known to confer drug resistance to transcriptome heterogeneity. Our methods reveal previously unseen isoform complexity and provide a better framework for multi-omics data integration.

Recycling of internalized receptors from endosomal compartments is essential for the receptors' cell-surface homeostasis. Sorting nexin 27 (SNX27) cooperates with the retromer complex in the recycling of proteins containing type I PSD95-Dlg-ZO1 (PDZ)-binding motifs. Here we define specific acidic amino acid sequences upstream of the PDZ-binding motif required for high-affinity engagement of the human SNX27 PDZ domain. However, a subset of SNX27 ligands, such as the β2 adrenergic receptor and N-methyl-D-aspartate (NMDA) receptor, lack these sequence determinants. Instead, we identified conserved sites of phosphorylation that substitute for acidic residues and dramatically enhance SNX27 interactions. This newly identified mechanism suggests a likely regulatory switch for PDZ interaction and protein transport by the SNX27-retromer complex. Defining this SNX27 binding code allowed us to classify more than 400 potential SNX27 ligands with broad functional implications in signal transduction, neuronal plasticity and metabolite transport.

Nephrogenic syndrome of inappropriate antidiuresis (NSIAD) is a genetic disease first described in 2 unrelated male infants with severe symptomatic hyponatremia. Despite undetectable arginine vasopressin levels, patients have inappropriately concentrated urine resulting in hyponatremia, hypoosmolality, and natriuresis. Here, we describe and functionally characterize a novel vasopressin type 2 receptor (V2R) gain-of-function mutation. An L312S substitution in the seventh transmembrane domain was identified in a boy presenting with water-induced hyponatremic seizures at the age of 5.8 years. We show that, compared with wild-type V2R, the L312S mutation results in the constitutive production of cAMP, indicative of the gain-of-function NSIAD profile. Interestingly, like the previously described F229V and I130N NSIAD-causing mutants, this appears to both occur in the absence of notable constitutive β-arrestin2 recruitment and can be reduced by the inverse agonist Tolvaptan. In addition, to understand the effect of various V2R substitutions on the full receptor "life-cycle," we have used and further developed a bioluminescence resonance energy transfer intracellular localization assay using multiple localization markers validated with confocal microscopy. This allowed us to characterize differences in the constitutive and ligand-induced localization and trafficking profiles of the novel L312S mutation as well as for previously described V2R gain-of-function mutants (NSIAD; R137C and R137L), loss-of-function mutants (nephrogenic diabetes insipidus; R137H, R181C, and M311V), and a putative silent V266A V2R polymorphism. In doing so, we describe differences in trafficking between unique V2R substitutions, even at the same amino acid position, therefore highlighting the value of full and thorough characterization of receptor function beyond simple signaling pathway analysis.

The endocytic protein SNX27 functions to link the parathyroid hormone receptor (PTHR) to the retromer trafficking complex. Loss of SNX27 in mice leads to overactive PTHR signaling and reduced osteoblastic bone formation during postnatal bone growth. Thus SNX27 is a new modulator of PTHR signaling., The parathyroid hormone 1 receptor (PTHR) is central to the process of bone formation and remodeling. PTHR signaling requires receptor internalization into endosomes, which is then terminated by recycling or degradation. Here we show that sorting nexin 27 (SNX27) functions as an adaptor that couples PTHR to the retromer trafficking complex. SNX27 binds directly to the C-terminal PDZ-binding motif of PTHR, wiring it to retromer for endosomal sorting. The structure of SNX27 bound to the PTHR motif reveals a high-affinity interface involving conserved electrostatic interactions. Mechanistically, depletion of SNX27 or retromer augments intracellular PTHR signaling in endosomes. Osteoblasts genetically lacking SNX27 show similar disruptions in PTHR signaling and greatly reduced capacity for bone mineralization, contributing to profound skeletal deficits in SNX27-knockout mice. Taken together, our data support a critical role for SNX27-retromer mediated transport of PTHR in normal bone development.

Parathyroid hormone-related protein (PTHrP) expression in breast cancer is enriched in bone metastases compared to primary tumors. Human MCF7 breast cancer cells "home" to the bones of immune deficient mice following intracardiac inoculation, but do not grow well and stain negatively for Ki67, thus serving as a model of breast cancer dormancy in vivo. We have previously shown that PTHrP overexpression in MCF7 cells overcomes this dormant phenotype, causing them to grow as osteolytic deposits, and that PTHrP-overexpressing MCF7 cells showed significantly lower expression of genes associated with dormancy compared to vector controls. Since early work showed a lack of cyclic AMP (cAMP) response to parathyroid hormone (PTH) in MCF7 cells, and cAMP is activated by PTH/PTHrP receptor (PTHR1) signaling, we hypothesized that the effects of PTHrP on dormancy in MCF7 cells occur through non-canonical (i.e., PTHR1/cAMP-independent) signaling. The data presented here demonstrate the lack of cAMP response in MCF7 cells to full length PTHrP(1-141) and PTH(1-34) in a wide range of doses, while maintaining a response to three known activators of adenylyl cyclase: calcitonin, prostaglandin E2 (PGE2), and forskolin. PTHR1 mRNA was detectable in MCF7 cells and was found in eight other human breast and murine mammary carcinoma cell lines. Although PTHrP overexpression in MCF7 cells changed expression levels of many genes, RNAseq analysis revealed that PTHR1 was unaltered, and only 2/32 previous PTHR1/cAMP responsive genes were significantly upregulated. Instead, PTHrP overexpression in MCF7 cells resulted in significant enrichment of the calcium signaling pathway. We conclude that PTHR1 in MCF7 breast cancer cells is not functionally linked to activation of the cAMP pathway. Gene expression responses to PTHrP overexpression must, therefore, result from autocrine or intracrine actions of PTHrP independent of PTHR1, through signals emanating from other domains within the PTHrP molecule.

Parathyroid hormone (PTH)-related protein (PTHrP) (gene name Pthlh) was discovered as the factor responsible for the humoral hypercalcemia of malignancy. It shares such sequence similarity with PTH in the amino-terminal region that the two are equally able to act through a single G protein-coupled receptor, PTH1R. A number of biological activities are ascribed to domains of PTHrP beyond the amino-terminal domain. PTH functions as a circulating hormone, but PTHrP is generated locally in many tissues including bone, where it acts as a paracrine factor on osteoblasts and osteocytes. The present study compares how PTH and PTHrP influence cyclic AMP (cAMP) formation through adenylyl cyclase, the first event in cell activation through PTH1R. Brief exposure to full length PTHrP(1-141) in several osteoblastic cell culture systems was followed by sustained adenylyl cyclase activity for more than an hour after ligand washout. This effect was dose-dependent and was not found with shorter PTHrP or PTH peptides even though they were fully able to activate adenylyl cyclase with acute treatment. The persistent activation response to PTHrP(1-141) was seen also with later events in the cAMP/PKA pathway, including persistent activation of CRE-luciferase and sustained regulation of several CREB-responsive mRNAs, up to 24 h after the initial exposure. Pharmacologic blockade of endocytosis prevented the persistent activation of cAMP and gene responses. We conclude that full length PTHrP, the likely local physiological effector in bone, differs in intracellular action to PTH by undergoing endosomal translocation to induce a prolonged adenylyl cyclase activation in its target cells.