Transcriptome profiling of tendon fibroblasts at the onset of embryonic muscle contraction reveals novel force-responsive genes

  1. Pavan K Nayak
  2. Arul Subramanian
  3. Thomas F Schilling  Is a corresponding author
  1. Department of Developmental and Cell Biology, University of California, United States
8 figures, 1 table and 8 additional files

Figures

Figure 1 with 1 supplement
Onset of embryonic muscle contraction regulates transcription in tenocytes.

(A–D) Diagrams depicting changes in tenocyte distribution and morphology during onset of trunk and cranial muscle contractions, (A) 36 hpf when twitching movements are sporadic and (B) 48 hpf when embryos become free swimming. Lateral views of 36 (A) and 48 hpf embryos (B). Insets show lateral and transverse views of migrating tenocyte progenitors (A) and differentiated tenocytes at somite boundaries with polarized, branched projections (B). Ventral views of the embryonic head in 60 hpf (C) and 72 hpf (D) embryos just prior to and during the onset of jaw movements. Cartilage (green), tenocytes (magenta), and muscles (cyan) showing tenocyte elongation, particularly in the sternohyoid tendon (sht) and condensation, as well as the mandibulohyoid junction (mhj). (E) Heatmaps from bulk RNA-sequencing (RNA-seq) showing the top 1000 differentially expressed genes (DEGs) between 36 and 48 hpf. p < 0.05. (F) Elevated expression of tenocyte marker genes mkxa, tnmd, and egr1 and extracellular matrix (ECM) genes col2a1a, col9a1a in RNA-seq experiments at 48 hpf. Datapoints represent normalized read counts of single biological replicates at each color-coded timepoint (n = 7 for 36 hpf, n = 4 for 48 hpf). (G) Elevated expression of cartilage marker genes col2a1a and col9a1a in 48 hpf samples. (H) PCA of individual replicates showing separation of experimental conditions by timepoint. (I) GO analysis using Biological Process (BP) terms of top 2788 DEGs by adjusted p-value.

Figure 1—figure supplement 1
FACS gating thresholds for mCherry+ cells.

(Left to right for each) Forward Scatter A (FSC-A) versus Side Scatter A (SSC-A) shows P1 threshold. FSC-A versus Forward Scatter H (FSC-H) shows P2 threshold to select for single cells. FSC-A versus mCherry-A shows P3 fluorescence gating for mCherry+ cells. (A) P3 selection gating allowed selection of cells with strong mCherry signal based on Negative control AB (WT) sample versus mCherry expressing tenocytes (48 hpf) FACS gating. Established P3 gating selected for mCherry-positive cells in all 36 and 48 hpf Tg(scxa:mCherry) samples. (B) Thresholds used in 36 and 48 hpf Tg(scxa:mCherry) samples for FACS prior to bulk RNA-sequencing (RNA-seq).

Figure 2 with 2 supplements
Expression of matn1, klf2a, and mxra5b with scxa in cranial and trunk tenocytes.

Ventral cranial (A–F, K–V) and lateral trunk (G–J) views of 72 hpf (A–F, K–V) and 48 hpf (G–J) embryos showing isHCR of matn1 (A, C–F), klf2a (G, I–K, M–P), and mxra5b (Q, S–V) in combination with scxa (B–F, H–J, L–P, R–V). (D–F, J, N–P, T–V) Higher magnification views of tenocyte nuclei in marked ROI. (C, D, M, N, S, T) ROI and panels outlined in magenta show magnified views of 3D volumes of tenocytes associated with imt. (I, J) ROI and panels outlined in cyan show magnified views of 3D volume of VMS tenocytes. (C, E, M, O, S, V) ROI and panels outlined in royal blue show magnified views of 3D volume of tenocytes associated with sht enthesis. (C, F, M, P) ROI and panels outlined in yellow show magnified views of 3D volumes of tenocytes associated with ch-pqt. (S, U) ROI and panels outlined in green show magnified views of 3D volumes of tenocytes associated with mhj. Each magnified view of ROI displays a translucent outline of the nuclear 3D volume with white puncta representing voxel colocalizations of isHCR as depicted by the colocalization function in Imaris (see Methods). mc – Meckel’s cartilage, pq – palatoquadrate, ch – ceratohyal, bh – basihyal cartilage, ima – intermadibularis anterior tendon, mhj – mandibulohyoid junction, sht – sternohyoideus tendon, ch-pqt – ceratohyal-palatoquadrate tendon, sb – somite boundary. Scale bars = 20 µm.

Figure 2—figure supplement 1
matn1, klf2a, and mxra5b are expressed in musculoskeletal tissues of developing embryos.

Lateral (A, B, D–I) and ventral (C) views of embryos showing expression of matn1 (A–C), klf2a (D–F), and mxra5b (G, H). (A–C) 48 hpf embryos show matn1 expression in cartilage progenitors at 48 hpf and in differentiated cartilages (and associated tenocytes) at 60 hpf (B, C). Lateral views of 36 hpf (D), 48 hpf (E), and 60 hpf (F) embryos show klf2a expression in pharyngeal mesenchyme (D), skeletal progenitors and in tenocytes along VMS (E, F). Lateral views of 36 hpf (G), 48 hpf (H), and 60 hpf (I) embryos show mxra5b expression in tenocytes along horizontal myosepta (HMS) along the notochord and VMS. Scale bars = 100 μm. Abbreviations: abc = anterior basicranial commissure, ch = ceratohyal cartilage, ep = ethmoid plate, hs = hyosymplectic cartilage, mc = meckel’s cartilage, nc = notochord, pf = pectoral fin, pq = palatoquadrate cartilage, sb = somite boundaries, t = trabeculae cartilage.

Figure 2—figure supplement 2
matn1 is expressed in differentiating cranial tenocytes.

Ventral view of the developing mandibular arch in a 51 hpf embryo showing in situ Hybridization Chain Reaction (isHCR) of matn1 (A, C, D) and scxa (B–D). (D) magnified view of yellow ROI (C) shows outline of tenocyte nuclear 3D volume with white puncta representing voxel colocalizations of matn1 and scxa as depicted by colocalization using Imaris (see methods). ep = ethmoid plate cartilage, ch = ceratohyal cartilage, mc = meckel’s cartilage, pq = palatoquadrate cartilage, ch-pqt = ceratohyal-palatoquadrate tendon, ima = intermandibularis anterior tendon, mhj = mandibulohyoid junction tendon, sht = sternohyoideus tendon. Scale bars = 20 µm.

Figure 3 with 1 supplement
Paralysis regulates tenocyte gene expression in developing musculoskeletal system.

(A) Heatmap of differentially expressed genes (DEGs) from bulk RNA-sequencing (RNA-seq) between WT and aBTX-injected (aBTX-inj) paralyzed 48 hpf embryos (force perturbed). (B) PCA of individual replicates WT versus aBTX-inj embryos’ RNA-seq separate by experimental condition. (C) Venn diagram shows overlap of genes between developmental time-point and force perturbed RNA-seq experiments. (D) Comparison of normalized read counts between replicates of matn1, klf2a, and mxra5b in 36 versus 48 hpf and WT versus aBTX RNA-seq experiments. (E) KEGG pathway analysis plot shows enrichment of overlapping genes from (C). ns = not significant, *p < 0.05, ***p < 0.001.

Figure 3—figure supplement 1
Paralysis regulates gene expression of matn1, klf2a, and mxra5b in developing embryos.

Bar plots showing global changes in relative expression from RT-qPCR of matn1, klf2a, and mxra5b genes between WT and cacnb1−/− mutant embryos at 48 hpf (A) and 72 hpf (B). Bar plots show global changes in relative expression of matn1, klf2a, and mxra5b between 48 hpf uninjected WT controls (blue) and aBTX-injected paralyzed (green) embryos (C), aBTX-injected paralyzed (green) and aBTX-injected ‘Twitching’ (partially recovered, magenta) embryos (D), and between aBTX-injected paralyzed (green) and aBTX-injected, ‘Recovered’ (cyan) embryos (E) at 48 hpf (right). ns = not significant, *p < 0.05, **p < 0.01, ***p < 0.001.

Figure 4 with 1 supplement
Mechanical force differentially regulates expression of matn1 in ima enthesis tenocytes.

Ventral views of Meckel’s cartilage and associated tenocytes showing in situ Hybridization Chain Reaction (isHCR) of matn1 (green) and anti-mCherry immunofluorescence (magenta) marking the tenocytes in Tg(scxa:mCherry) embryos at 72 hpf in WT uninjected (WT) (A–E), aBTX-inj (Paralyzed) (F–J), partially recovered aBTX-inj (Twitching) (K–O), and completely recovered aBTX-inj (Full Recovery) (P–T) conditions at ima enthesis. (D, I, N, S) Grayscale images showing nuclei stained with DAPI with ROIs showing isolated 3D volumes of chondrocytes (green) and enthesis tenocytes (magenta) based on DAPI signal. (E, J, O, T) Insets showing magnified views of the 3D volumes of tenocytes associated with ima enthesis depicting expression of matn1 and stained for mCherry. (U) Violin plot showing changes in mean fluorescence intensity of matn1 in ima enthesis tenocyte nuclei between WT (n = 8), Paralyzed (n = 8), Twitching (n = 6), and Full Recovery (n = 7) with ~8 nuclei measured per embryo. p-value calculated with linear mixed effects model with Tukey post hoc test. *p < 0.05, **p < 0.01, ***p < 0.001. Scale bars = 20 µm.

Figure 4—source data 1

Measurements of matn1 isHCR signal intensity in ima enthesis tenocytes.

https://cdn.elifesciences.org/articles/105802/elife-105802-fig4-data1-v2.xlsx
Figure 4—figure supplement 1
Mechanical force regulates expression of matn1 and klf2a in sht enthesis tenocytes.

Ventral views of ceratohyal (ch) cartilage and associated tenocytes showing in situ Hybridization Chain Reaction (isHCR) of matn1 (green) (A–L) and klf2a (green) (N–Y) and anti-mCherry immunofluorescence (magenta) marking the tenocytes in Tg(scxa:mCherry) embryos at 72 hpf in WT uninjected (WT) (A–C, N–P), aBTX-inj (Paralyzed) (D–F, Q–S), partially recovered aBTX-inj (Twitching) (G–I, T–V), and completely recovered aBTX-inj (Full Recovery) (J–L, W–Y) conditions at sht enthesis. (B, E, H, K, O, R, U, X) Grayscale images showing nuclei stained with DAPI with ROIs showing isolated 3D volumes of chondrocytes (green) and sternohyoideus-ceratohyal (sht) enthesis tenocytes (magenta) based on DAPI signal. (C, F, I, L, P, S, V, Y) Insets showing magnified views of the 3D volumes of tenocytes associated with sht enthesis depicting expression of matn1 and stained for mCherry. (M) Violin plot showing changes in mean fluorescence intensity of matn1 in sht enthesis tenocyte nuclei between WT (n = 8), Paralyzed (n = 8), Twitching (n = 6), and Full Recovery (n = 7) with ~8 nuclei measured per embryo. (Z) Violin plot showing changes in mean fluorescence intensity of klf2a in sht enthesis tenocyte nuclei between WT (n = 15), Paralyzed (n = 16), Twitching (n = 14), and Full Recovery (n = 11) with ~8 nuclei measured per embryo. p-values calculated with linear mixed effects model with Tukey post hoc test. ***p < 0.001. Scale bars = 20 µm.

Figure 4—figure supplement 1—source data 1

Measurements of matn1 and klf2a isHCR signal intensity in sht enthesis tenocytes.

https://cdn.elifesciences.org/articles/105802/elife-105802-fig4-figsupp1-data1-v2.xlsx
Figure 5 with 3 supplements
Mechanical force differentially regulates expression of mxra5b in sht enthesis tenocytes.

Ventral views of ceratohyal (ch) cartilage and associated tenocytes showing in situ Hybridization Chain Reaction (isHCR) of mxra5b (green) and anti-mCherry immunofluorescence (magenta) marking the tenocytes in Tg(scxa:mCherry) embryos at 72 hpf in WT uninjected (WT) (A–E), aBTX-inj paralyzed (F–J), partially recovered aBTX-inj (Twitching) (K–O), and completely recovered aBTX-inj (Full Recovery) (P–T) conditions at sht enthesis. (D, I, N, S) Grayscale images showing nuclei stained with DAPI with ROIs showing isolated 3D volumes of chondrocytes (green) and sht enthesis tenocytes (magenta) based on DAPI signal. (E, J, O, T) Insets showing magnified views of the 3D volumes of tenocytes associated with sht enthesis depicting expression of mxra5b and stained for mCherry. (U) Violin plot showing changes in mean fluorescence intensity of mxra5b in sht enthesis tenocyte nuclei between WT (n = 7), Paralyzed (n = 8), Twitching (n = 8), and Full Recovery (n = 4) with ~8 nuclei measured per embryo. p-value calculated with linear mixed effects model with Tukey post hoc test. *p < 0.05, **p < 0.01. Scale bars = 20 µm.

Figure 5—source data 1

Measurements of mxra5b isHCR signal intensity in sht enthesis tenocytes.

https://cdn.elifesciences.org/articles/105802/elife-105802-fig5-data1-v2.xlsx
Figure 5—figure supplement 1
Mechanical force differentially regulates expression of mxra5b and klf2a in ima enthesis tenocytes.

Ventral views of Meckels cartilage and associated tenocytes showing in situ Hybridization Chain Reaction (isHCR) of klf2a (green) (A–L) and mxra5b (green) (N–Y) and anti-mCherry immunofluorescence (magenta) marking the tenocytes in Tg(scxa:mCherry) embryos at 72 hpf in WT uninjected (WT) (A–C, N–P), aBTX-inj (Paralyzed) (D–F, Q–S), partially recovered aBTX-inj (Twitching) (G–I, T–V), and completely recovered aBTX-inj (Full Recovery) (J–L, W–Y) conditions at ima enthesis. (B, E, H, K, O, R, U, X) Grayscale images showing nuclei stained with DAPI with ROIs showing isolated 3D volumes of chondrocytes (green) and ima enthesis tenocytes (magenta) based on DAPI signal. (C, F, I, L, P, S, V, Y) Insets showing magnified views of the 3D volumes of tenocytes associated with ima enthesis depicting expression of mxra5b and klf2a and stained for mCherry. (M) Violin plot showing changes in mean fluorescence intensity of klf2a in ima enthesis tenocyte nuclei between WT (n = 15), Paralyzed (n = 16), Twitching (n = 14), and Full Recovery (n = 11) with ~8 nuclei measured per embryo. (Z) Violin plot showing changes in mean fluorescence intensity of mxra5b in ima enthesis tenocyte nuclei between WT (n = 7), Paralyzed (n = 8), Twitching (n = 8), and Full Recovery (n = 4) with ~8 nuclei measured per embryo. p-values calculated with linear mixed effects model with Tukey post hoc test. ***p < 0.001. Scale bars = 20 µm.

Figure 5—figure supplement 1—source data 1

Measurements of klf2a and mxra5b isHCR signal intensity in ima enthesis tenocytes.

https://cdn.elifesciences.org/articles/105802/elife-105802-fig5-figsupp1-data1-v2.xlsx
Figure 5—figure supplement 2
Mechanical force regulates expression of mxra5b in mhj myotendinous junction tenocytes.

Ventral views of mandibulohyoid (mhj) myotendinous junction (MTJ) associated tenocytes showing in situ Hybridization Chain Reaction (isHCR) of mxra5b (green) and anti-mCherry immunofluorescence (magenta) marking the tenocytes in Tg(scxa:mCherry) embryos at 72 hpf in WT uninjected (WT) (A–C), aBTX-inj (Paralyzed) (D–F), partially recovered aBTX-inj (Twitching) (G–I), and completely recovered aBTX-inj (Full Recovery) (J–L) conditions. (B, E, H, K) Grayscale images showing nuclei stained with DAPI with ROIs showing isolated 3D volumes of mhj tenocytes (magenta) based on DAPI signal. (C, F, I, L) Insets showing magnified views of the 3D volumes of tenocytes associated with mhj MTJ depicting expression of mxra5b and stained for mCherry. (M) Violin plot showing changes in mean fluorescence intensity of mxra5b in mhj MTJ tenocyte nuclei between WT (n = 7), Paralyzed (n = 8), Twitching (n = 8), and Full Recovery (n = 4) with ~10 nuclei measured per embryo. p-value calculated with linear mixed effects model with Tukey post hoc test. ***p < 0.001. Scale bars = 20 µm.

Figure 5—figure supplement 2—source data 1

Measurements of mxra5b isHCR signal intensity in mhj MTJ tenocytes.

https://cdn.elifesciences.org/articles/105802/elife-105802-fig5-figsupp2-data1-v2.xlsx
Figure 5—figure supplement 3
Mechanical force differentially regulates expression of mxra5b and klf2a in sht myotendinous junction tenocytes.

Ventral views of sht and associated myotendinous junction (MTJ) Tenocytes showing in situ Hybridization Chain Reaction (isHCR) of mxra5b (green) (A–L) and klf2a (green) (N–Y) and anti-mCherry immunofluorescence (magenta) marking the tenocytes in Tg(scxa:mCherry) embryos at 72 hpf in WT uninjected (WT) (A–C, N–P), aBTX-inj (Paralyzed) (D–F, Q–S), partially recovered aBTX-inj (Twitching) (G–I, T–V), and completely recovered aBTX-inj (Full Recovery) (J–L, W–Y) conditions at sht MTJ. (B, E, H, K, O, R, U, X) Grayscale images showing nuclei stained with DAPI with ROIs showing isolated 3D volumes of ima enthesis tenocytes (magenta) based on DAPI signal. (C, F, I, L, P, S, V, Y) Insets showing magnified views of the 3D volumes of tenocytes associated with sht enthesis depicting expression of mxra5b and klf2a and stained for mCherry. (M) Violin plot showing changes in mean fluorescence intensity of mxra5b in sht MTJ tenocyte nuclei between WT (n = 7), Paralyzed (n = 8), Twitching (n = 8), and Full Recovery (n = 4) with ~4 nuclei measured per embryo. (Z) Violin plot showing changes in mean fluorescence intensity of klf2a in ima enthesis tenocyte nuclei between WT (n = 15), Paralyzed (n = 16), Twitching (n = 14), and Full Recovery (n = 11) with ~4 nuclei measured per embryo. p-values calculated with linear mixed effects model with Tukey post hoc test. *p < 0.05, **p < 0.01. Scale bars = 20 µm.

Figure 5—figure supplement 3—source data 1

Measurements of mxra5b and klf2a isHCR signal intensity in sht MTJ tenocytes.

https://cdn.elifesciences.org/articles/105802/elife-105802-fig5-figsupp3-data1-v2.xlsx
Mechanical force regulates expression of klf2a in mhj myotendinous junction tenocytes.

Ventral views of mandibulohyoid junction (mhj), myotendinous junction (MTJ) associated tenocytes showing in situ Hybridization Chain Reaction (isHCR) of klf2a (green) and anti-mCherry immunofluorescence (magenta) marking the tenocytes in Tg(scxa:mCherry) embryos at 72 hpf in WT uninjected (WT) (A–E), aBTX-inj (Paralyzed) (F–J), partially recovered aBTX-inj (Twitching) (K–O), and completely recovered aBTX-inj (Full Recovery) (P–T) conditions. (D, I, N, S) Grayscale images showing nuclei stained with DAPI with ROIs showing isolated 3D volumes of mhj tenocytes (magenta) based on DAPI signal. (E, J, O, T) Insets showing magnified views of the 3D volumes of tenocytes associated with mhj MTJ depicting expression of klf2a and stained for mCherry. (U) Violin plot showing changes in mean fluorescence intensity of klf2a in mhj MTJ tenocyte nuclei between WT (n = 17), Paralyzed (n = 15), Twitching (n = 14), and Full Recovery (n = 11) with ~10 nuclei measured per embryo. p-value calculated with linear mixed effects model with Tukey post hoc test. *p < 0.05, ***p < 0.001. Scale bars = 20 µm.

Figure 6—source data 1

Measurements of klf2a isHCR signal intensity in mhj MTJ tenocytes.

https://cdn.elifesciences.org/articles/105802/elife-105802-fig6-data1-v2.xlsx
Loss of mxra5b function affects somite boundary structure.

(A) Lateral views of WT and mxra5b multiplex CRISPants at 48 and 72 hpf Tg(scx:mCherry) embryos stained with anti-mCherry to show tenocytes at the somite boundary (SB). (B) Quantification of somite boundary angle measurements of 48 hpf WT or mxra5b multiplex CRISPant embryos. p-value calculated with Watson’s U2 test. *p < 0.05.

Model depicting role of mechanical force in regulating expression of genes in tenocytes during onset of active muscle contraction.

(A) Cartoon showing role of force in regulating tenocyte morphogenesis and gene expression in tenocytes between 36 and 48 hpf stages correlating with onset of active swimming The variability in gene expression is related to increase in both magnitude and persistence of muscle contraction force. (B) Representative model summarizing the multifaceted role of muscle contractile force on expression dynamics of matn1, klf2a, and mxra5b genes in cranial tendon attachments. (C) Force-responsive gene expression is more nuanced than a binary on/off control.

Tables

Key resources table
Reagent type (species) or resourceDesignationSource or referenceIdentifiersAdditional information
Strain (Danio rerio)ABSchilling labRRID:NCBITaxon_7955
Genetic reagent (Danio rerio)Tg(scxa:mCherry)Galloway labfb301Tg; RRID:ZFIN_ZDB-GENO-180925-6scx BAC transgenic in AB background
Genetic reagent (Danio rerio)cacnb1−/−; Tg(scxa:mCherry)Schilling labIr1092;fb301; RRID:ZFIN_ZDB-ALT-191023-1cacnb1 mutant in Tg(scx:mCherry) background
Sequence-based reagentT7 sequence-tagged primersThis paperSupplementary file 72 mM final concentration
Commercial assay or kitProtoscript II first strand cDNA synthesis kitNew England BiolabsCat # E6560
Commercial assay or kitT7 RNA polymeraseMillipore Sigma (Roche)Cat # 10881767001
Commercial assay or kitMonarch Total RNA Miniprep kitNew England BiolabsCat # T2010S
Commercial assay or kitDIG RNA labeling mixMillipore Sigma (Roche)Cat # 11277073910
Commercial assay or kitMEGAshortscript T7 transcription kitThermo Fisher Scientific (Invitrogen)Cat # AM1354
Commercial assay or kitLuna Universal qPCR master mixNew England BiolabsCat # M3003S
Commercial assay or kitZirconium beadsBenchmark ScientificCat # D1032-10
Commercial assay or kitRNEasy Micro KitQIAGENCat # 74004
Commercial assay or kit40 µm filterPluriselect-USACat # 43-10040-50
Commercial assay or kitHCR Buffers (v3.0)Molecular InstrumentsHybridization buffer, Wash buffer, Amplifier buffer
AntibodyAnti-Digoxigenin-AP, Fab fragmentsMillipore Sigma (Roche)Cat # 11093274910 RRID:AB_5144971:2000
AntibodyRat monoclonal anti-mCherry antibodyInvitrogen (Thermo Fisher Scientific)Cat # M11217 RRID:AB_25366111:500
AntibodyChicken polyclonal anti-GFP antibodyabcamCat # ab13970 RRID:AB_3007981:1000
AntibodyMouse monoclonal anti-myosin heavy chain antibodyDevelopmental Studies Hybridoma Bank (DHSB)Cat # A4.1025 RRID:AB_5283561:200
AntibodyAlexa Fluor 594 AffiniPure F(ab’)2 Fragment Donkey polyclonal anti-Rat IgG (H+L)Jackson ImmunoResearch LaboratoriesCat # 712-586-153 RRID:AB_23406911:1000
AntibodyAlexa Fluor 488 AffiniPure F(ab’)2 Fragment Donkey polyclonal anti-Chicken IgY IgG (H+L)Jackson ImmunoResearch LaboratoriesCat # 703-546-155 RRID:AB_23403761:1000
AntibodyAlexa Fluor 647 AffiniPure F(ab’)2 Fragment Donkey polyclonal anti-Mouse IgG (H+L)Jackson ImmunoResearch LaboratoriesCat # 715-606-151 RRID:AB_23408661:1000
Chemical compound, drugNitro Blue Tetrazolium chloride solution (NBT)Millipore Sigma (Roche)Cat # 11383213001
PubChem CID: 9281
Chemical compound, drug5-Bromo-4-chloro-3-indolyl phosphate solutionMillipore Sigma (Roche)Cat # 11383221001
PubChem CID: 81059
Chemical compound, drugEthylenediaminetetraacetic acid disodium saltMillipore Sigma (Roche)Cat # E5134
PubChem CID: 8759
Chemical compound, drugCalcium chloride hexahydrateMillipore Sigma (Roche)Cat # 21108
PubChem CID: 6093252
Chemical compound, drugDulbecco’s phosphate-buffered saline (DPBS) 1×Thermo Fisher Scientific (Gibco)Cat # 14190144
Chemical compound, drugAgarose low gelling temperatureMillipore Sigma (Sigma-Aldrich)Cat # A9414
Chemical compound, drugSSC buffer 20×Millipore Sigma (Sigma-Aldrich)Cat # S6639-1L
Chemical compound, drugDAPIMillipore Sigma (Sigma-Aldrich)Cat # D9542
PubChem CID: 2954
Sequence-based reagentmatn1-B1Molecular instrumentsNM_001099740.220 probe set
Sequence-based reagentmxra5b-B1Molecular instrumentsXM_017357865.220 probe set
Sequence-based reagentmxra5b-B1Molecular instrumentsXM_017357865.220 probe set
Sequence-based reagentklf2a-B3Molecular instrumentsNM_131856.320 probe set
Sequence-based reagentscxa-B2Molecular instrumentsNM_00108306918 probe set
Sequence-based reagentB1-h1&h2- Alexa Fluor 488 amplifier hairpinsMolecular instrumentsHCR RNA-FISH (v3.0)
Sequence-based reagentB2-h1&h2- Alexa Fluor 546 amplifier hairpinsMolecular instrumentsHCR RNA-FISH (v3.0)
Sequence-based reagentB3-h1&h2- Alexa Fluor 647 amplifier hairpinsMolecular instrumentsHCR RNA-FISH (v3.0)
OtherBD FACSAria II Cell SorterBecton, Dickinson and CompanyRRID:SCR_018934
OtherBioanalyzer 2100 instrumentAgilentRRID:SCR_018043
OtherHiSeq 4000 sequencing systemIlluminaRRID:SCR_016386
OtherNextSeq 550 systemIlluminaRRID:SCR_016381
OtherLightCycler 480 Real Time PCR SystemRocheRRID:SCR_018626
OtherSP8 Lightning Confocal microscopeLeicaRRID:SCR_018169
OtherZeiss Axioplan 2 imaging systemZeissRRID:SCR_020918
OtherMicroPublisher color RTV-5.0 CCD cameraQImaging
OtherBeadBug 3 microtube homogenizerBenchmark ScientificCat # D1030
Peptide, recombinant proteinProtease (Subtilisin Carlsberg) from Bacillus licheniformisMillipore Sigma (Sigma-Aldrich)Cat # P5380
UniProtKB: P00780.SUBC_BACLI
Peptide, recombinant proteinCollagenase Type IV from Hathewaya histolytica (Clostridium histolyticum)Thermo Fisher Scientific (Gibco Life technologies)Cat # 17104019
Peptide, recombinant proteinDeoxyribonuclease I (DNase I) from bovine pancreasMillipore Sigma (Roche)Cat # 10104159001 UniProtKB: P00639.DNAS1_BOVIN
Peptide, recombinant proteinBovine serum albumin stock solution (10%)Miltenyi BiotecCat # 130-091-376
Recombinant DNA reagentpmtb-t7-alpha-bungarotoxinAddgene (Megason lab)Cat # 69542
RRID:Addgene_69542
Software, algorithmSpliced Transcripts Alignment to a Reference (STAR) v2.5.2aDobin labRRID:SCR_004463
Software, algorithmSmart-seq2 single sample pipelineBroad InstituteRRID:SCR_021228
Software, algorithmRSEM v1.2.31Dewey labRRID:SCR_000262
Software, algorithmDESeq2 v1.30.1Anders labRRID:SCR_015687
Software, algorithmClustVisVilo labRRID:SCR_017133
Software, algorithmClusterProfiler R packageQing-Yu labRRID:SCR_016884
Software, algorithmShinyGOGe labRRID:SCR_019213
Software, algorithmVennDiagram v1.7.3Boutros labRRID:SCR_002414
Software, algorithmGeneOverlap v1.26.0Shen labRRID:SCR_018419
Software, algorithmLightCycler SoftwareRocheRRID:SCR_012155
Software, algorithmZeiss Zen Microscopy softwareZeissRRID:SCR_013672
Software, algorithmLeica Application Suite XLeicaRRID:SCR_013673
Software, algorithmImarisBitplaneRRID:SCR_007370
OtherOptical Biology Core at UCIDepartment of Developmental Biology, UCIRRID:SCR_026614Core facility
OtherGenomics Research and Technology Hub Core at UCIDepartment of Biological Chemistry, UCIRRID:SCR_026615Core facility
OtherFlow Cytometry Core at UCIStem Cell Research Center, UCIRRID:SCR_026616Core facility

Additional files

Supplementary file 1

Differentially expressed gene list of bulk RNA-seq of sorted mCherry+ tenocytes from 36 hpf vs 48 hpf embryos.

https://cdn.elifesciences.org/articles/105802/elife-105802-supp1-v2.csv
Supplementary file 2

ShinyGO analysis of 36 vs. 48 hpf bulk RNA-seq differentially expressed genes.

https://cdn.elifesciences.org/articles/105802/elife-105802-supp2-v2.csv
Supplementary file 3

DAVID analysis of 36 vs. 48 hpf bulk RNA-seq differentially expressed genes.

https://cdn.elifesciences.org/articles/105802/elife-105802-supp3-v2.xlsx
Supplementary file 4

Differentially expressed gene list from bulk RNA-seq of sorted mCherry+ tenocytes from 48 hpf WT vs aBTX-injected embryos.

https://cdn.elifesciences.org/articles/105802/elife-105802-supp4-v2.csv
Supplementary file 5

List of differentially expressed genes overlapping between 36 hpf vs 48 hpf bulk RNA-seq and 48 hpf WT vs a-BTX injected paralysis bulk RNA-seq.

https://cdn.elifesciences.org/articles/105802/elife-105802-supp5-v2.csv
Supplementary file 6

ShinyGO analysis of differentially expressed genes overlapping between 36 hpf vs 48 hpf bulk RNA-seq and 48 hpf WT vs a-BTX injected paralysis bulk RNA-seq.

https://cdn.elifesciences.org/articles/105802/elife-105802-supp6-v2.csv
Supplementary file 7

List of primers used for chromogenic in situ hybridizations and RT-qPCRs.

https://cdn.elifesciences.org/articles/105802/elife-105802-supp7-v2.xlsx
MDAR checklist
https://cdn.elifesciences.org/articles/105802/elife-105802-mdarchecklist1-v2.pdf

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Pavan K Nayak
  2. Arul Subramanian
  3. Thomas F Schilling
(2025)
Transcriptome profiling of tendon fibroblasts at the onset of embryonic muscle contraction reveals novel force-responsive genes
eLife 14:e105802.
https://doi.org/10.7554/eLife.105802