A novel cis-regulatory element drives early expression of Nkx3.2 in the gnathostome primary jaw joint

  1. Jake Leyhr
  2. Laura Waldmann
  3. Beata Filipek-Górniok
  4. Hanqing Zhang
  5. Amin Allalou
  6. Tatjana Haitina  Is a corresponding author
  1. Subdepartment of Evolution and Development, Department of Organismal Biology, Uppsala University, Sweden
  2. Science for Life Laboratory Genome Engineering Zebrafish Facility, Department of Organismal Biology, Uppsala University, Sweden
  3. Division of Visual Information and Interaction, Department of Information Technology, Uppsala University, Sweden
  4. Science for Life Laboratory BioImage Informatics Facility, Sweden
8 figures, 1 video and 1 additional file

Figures

Gene synteny around Nkx3.2 is well conserved in vertebrate genomes.

Phylogenetic tree of vertebrates based on commonly accepted topology. Pointed boxes represent gene orientation with gene names indicated on top. Red hexagons mark the position of conserved non-coding element (CNE) (jaw joint regulatory sequence 1, JRS1) downstream of Nkx3.2, where filled hexagons mark CNEs selected for in vivo functional characterization in this study. The corresponding chromosome/contig number of each locus is indicated below the gene order schematic of each species. Daggers indicate extinct species.

Figure 2 with 2 supplements
A conserved non-coding element, jaw joint regulatory sequence 1 (JRS1), identified using mVISTA and MEME.

(A) mVISTA alignment of vertebrate Nkx3.2 loci, using the spotted gar locus as reference. Peaks indicate conserved sequences >50% identity, coloured peaks indicate >70% identity. Dark blue peaks indicate conserved exon sequences, pink indicates conserved non-coding sequences. (B) Shared sequence motifs (1–4) in the core ~245 bp sequence of JRS1 in different species identified with MEME analysis. p values indicated per species (horizontal) and per motif (vertical). (C) Relevant transcription factor-binding sites predicted by Tomtom at motifs 1–4 with associated p values for each match.

Figure 2—figure supplement 1
Multiple sequence alignment, genome coordinates, and percent identity matrix of the jaw joint regulatory sequence 1 (JRS1) conserved core.

Asterisks mark sites in the alignment conserved in 6/6 species, dots mark sites conserved in ≥4/6 species.

Figure 2—figure supplement 2
Jaw joint regulatory sequence 1 (JRS1) is differentially bound by transcription factors in the first and second embryonic mouse branchial arches.

Schematic view of the transcription factor binding at JRS1 in E11.5 mouse BA1 and BA2, determined by ChIP-seq assays performed by Amin et al., 2015 and Donaldson et al., 2012. The height of each bar represents the maximum enrichment at the summit of the enrichment peak, and the width of the bar represents the width of the summit region (200 bp). The genomic location of each enrichment summit region in relation to JRS1 is precise and to scale, while the Nkx3.2 and Rab28 gene icons are purely schematic to indicate the strand direction of JRS1.

Figure 2—figure supplement 2—source data 1

Genomic coordinates of transcription factor-binding enrichment at jaw joint regulatory sequence 1 (JRS1) in E11.5 mouse branchial arch (BA) 1 and 2 determined by ChIP-seq assays.

These data were obtained from the publicly available supplementary information of Amin et al., 2015 and Donaldson et al., 2012 . Genomic coordinates correspond to chromosome 5 in the mm9 mouse assembly.

https://cdn.elifesciences.org/articles/75749/elife-75749-fig2-figsupp2-data1-v1.xlsx
Jaw joint regulatory sequence 1 (JRS1) is absent in most acanthopterygian fish.

Phylogenetic tree based on Hughes et al., 2018 and Ghezelayagh et al., 2022. Pointed boxes represent gene orientation with gene names indicated on top. The slimehead nkx3.2 box is displayed as slightly grayed out as only the 3’ UTR was present in the contig. Red hexagons mark the presence and position of JRS1, where filled hexagons mark enhancers selected for in vivo functional characterization in this study. Small red rectangles mark the position of JRS1 fragments. Below the gene order schematic of each species is marked the chromosome or contig containing this region. The plot of BLASTN E-Values represents the E-Values of the top hits to bichir (circles) and spotted gar (triangles) JRS1 query sequences in the genomes of each species. Shapes filled with white or black indicate that the top BLASTN hit using the bichir or spotted gar JRS1 query sequence, respectively, was found between the nkx3.2 and rab28 genes, while grey-filled shapes indicate the top hit was found in a different locus. The x-axis is truncated on the right side to a minimum of 1e−30 for ease of comparison, and the red dotted line indicates an E-Value of 1e−04. Taxonomic orders and higher clade classifications are shown on the far right.

Figure 3—source data 1

The E-values of the top BLASTN hits using the bichir and spotted gar jaw joint regulatory sequence 1 (JRS1) query sequences in the genomes of each species.

https://cdn.elifesciences.org/articles/75749/elife-75749-fig3-data1-v1.xlsx
Zebrafish jaw joint regulatory sequence 1 (JRS1) enhancer drives mCherry reporter gene expression in jaw joint-forming chondroprogenitor cells and partly overlaps with sox10:egfp expressing cells.

(A, A’) At 2 dpf jaw joint progenitor cells express both GFP and mCherry. (B, B’) By 3 dpf cells start to differentiate and nkx3.2(JRS1):mCherry activity is restricted to jaw joint-forming interzone, overlapping with sox10:egfp. Single-labelled mCherry-expressing cells are surrounding the joint-forming region. (C, C’) At 5 dpf mCherry-expressing cells are restricted to the articulation forming area between Meckel’s cartilage (m) and the palatoquadrate (pq). Double mCherry/GFP-expressing cells are restricted to posterior Meckel’s cartilage and anterior palatoquadrate. (D, D’) At 14 dpf a clear joint cavity is visible. nkx3.2(JRS1):mCherry activity is restricted to the joint cavity and to both lateral and medial palatoquadrate. Dashed box in (A–D) is magnified in (A’–D’). (A–D) Represents maximum projection of confocal Z-stack, and (A’–D’) represents a single confocal image. In nkx3.2−/− mutants nkx3.2(JRS1):mCherry marks the cells outside of the fused jaw joint at 3 dpf (E, E’) and 5 dpf (F, F’). Dashed boxes in (E, F) are magnified in (E’, F’). (E, F) and (E’, F’) represent maximum projection of confocal Z-stack. Asterisks indicate the approximate location of the fusion site. Scale bars: 100 μm (A–F) and 25 μm (A’–F’).

Functional conservation of the jaw joint regulatory sequence 1 (JRS1) enhancer within tested gnathostome species.

(A–C, G–I) Maximum projection images of 3 dpf transgenic zebrafish embryos (ventral view) driving mCherry expression in jaw joint and mandibular arch elements under the control of the JRS1 sequence of (A) human Homo sapiens, (B) mouse Mus musculus, (C) frog Xenopus tropicalis, (G) zebrafish Danio rerio, (H) bichir Polypterus senegalus and (I) elephant shark Callorhinchus milii. Asterisks mark a jaw joint. (D–F and J–L) Maximum projection images of 3 dpf nkx3.2(JRS1):mCherry transgenic zebrafish driving mCherry expression under the control of species-specific enhancer sequences with fli1a:egfp background reveals mCherry expression in GFP-labelled perichondrium cells. Dashed boxes are magnified in (D’–F’) and (J’–L’) as single confocal images. White arrowheads mark mCherry expression in the perichondrium. m: Meckel’s cartilage; pq: palatoquadrate. Scale bars: 75 μm (A–L), 25 μm (D’–F’ and J’–L’).

Figure 6 with 1 supplement
Homozygous jaw joint regulatory sequence 1 (JRS1) enhancer deletion results in early jaw joint dysmorphology.

(A) Schematic of the JRS1 deletion allele generation. (B–J) Alcian blue and alizarin red-stained jaw joints of representative wild-type, nkx3.2+/ΔJRS1, and nkx3.2ΔJRS1/ΔJRS1 zebrafish at 5, 14, and 30 dpf. Phenotypes are quantified in the top right of each panel. Standard lengths (mm) are given for 14 and 30 dpf juveniles. White arrowhead marks the partial fusion between Meckel’s cartilage (m) and the palatoquadrate (pq). Grey arrowheads indicate the rounded, non-bulbous anteroposterior process of Meckel’s cartilage. Black arrowheads mark the reduced retroarticular process (ra). Asterisks mark the ossifying retroarticular bone. Scale bars: 25 µm. (K–N) Representative embryos following in situ hybridization staining for nkx3.2 in the ventral craniofacial region at 48 and 56 hpf. Arrows indicate the reduced expression in the jaw joint domains in nkx3.2ΔJRS1/ΔJRS1 embryos. (O) nkx3.2 expression levels at 48 and 56 hpf in wild-type and nkx3.2+/ΔJRS1 versus nkx3.2ΔJRS1/ΔJRS1 zebrafish embryos, measured as pixel intensity of individual jaw joint expression domains in in situ stained embryos (N = 16, 18, 18, 17, respectively). (P) Relative whole-body nkx3.2 expression levels at 6 dpf in wild-type, nkx3.2+/ΔJRS1, and nkx3.2ΔJRS1/ΔJRS1 zebrafish, determined by qPCR. NS. denotes p > 0.05, ** p < 0.01, *** p < 0.001. Error bars represent mean ± standard error of the mean (SEM).

Figure 6—source data 1

Relative pixel intensity values obtained from in situ images of nkx3.2 expression in zebrafish at 48 and 56 hpf, and statistical analyses.

https://cdn.elifesciences.org/articles/75749/elife-75749-fig6-data1-v1.xlsx
Figure 6—source data 2

qPCR raw values.

https://cdn.elifesciences.org/articles/75749/elife-75749-fig6-data2-v1.xlsx
Figure 6—figure supplement 1
Alignment of wild-type and uu3731 jaw joint regulatory sequence 1 (JRS1) deletion alleles with annotated MEME motifs, primers, and CRISPR gRNA targets.

The conserved core of the zebrafish JRS1 enhancer sequence identified by MEME is labelled in grey. Motifs 1–3 are coloured according to Figure 2B. Scissors mark the beginning and end of the deletion location, and an asterisk marks the chance insertion of a random 7-bp sequence to the deletion location during DNA repair.

Figure 7 with 1 supplement
Geometric morphometric analysis of jaw joint regulatory sequence 1 (JRS1) deletion phenotypes at 9 dpf.

(A) Principal components analysis of geometric morphometric comparison of posterior Meckel’s cartilage shape in wild-type, nkx3.2+/ΔJRS1, and nkx3.2ΔJRS1/ΔJRS1 zebrafish at 9 dpf. Thin-plate splines display the extremes of shape along PC1 and PC2. Inset is a histogram showing the percentage of variance explained by PCs 1–5. (B–G) Lateral and dorsal 3D renderings of the left jaw joint from averaged optical projection tomography (OPT) models of wild-type (N = 12), nkx3.2+/ΔJRS1 (N = 10), and nkx3.2ΔJRS1/ΔJRS1 (N = 8) zebrafish at 9 dpf. White arrowheads mark the partial fusion between Meckel’s cartilage (m) and the palatoquadrate (pq). Black arrowhead marks the reduced retroarticular process (ra). Scale bars: 25 µm.

Figure 7—figure supplement 1
Optical projection tomography of 9 dpf jaw joint regulatory sequence 1 (JRS1) deletion zebrafish.

(A–C) Lateral and (D–F) ventral maximum projection views of the average head skeleton of 9 dpf wild-type (N = 12), nkx3.2+/ΔJRS1 (N = 10), and nkx3.2ΔJRS1/ΔJRS1 (N = 8) zebrafish. White arrowhead indicates the subtle jaw joint dysmorphology seen in nkx3.2ΔJRS1/ΔJRS1 mutants. Scale bars: 150 µM.

Graphical summary of this study.

The conserved jaw joint regulatory sequence 1 (JRS1) enhancer was first identified as a conserved non-coding element (CNE) in silico, then confirmed with transgenesis experiments of a range of species-specific JRS1 sequences driving fluorescent reporter expression in the jaw joint. Finally, the JRS1 enhancer was deleted from the zebrafish genome to reveal a transient jaw joint dysmorphology and partial fusion. m: Meckel’s cartilage; pq: palatoquadrate; ra: retroarticular process.

Videos

Video 1
Jaw joint regulatory sequence 1 (JRS1) enhancer drives reporter expression in the early notochord and jaw joints.

Lateral view of nkx3.2(JRS1):mCherry/sox10:egfp zebrafish embryo developing from 44 to 71 hpf. Inset in situ hybridization image of nkx3.2 expression in a wild-type 56 hpf zebrafish embryo. White arrowhead indicates jaw joint expression domain while black arrowheads indicate expression in the notochord.

Additional files

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. Jake Leyhr
  2. Laura Waldmann
  3. Beata Filipek-Górniok
  4. Hanqing Zhang
  5. Amin Allalou
  6. Tatjana Haitina
(2022)
A novel cis-regulatory element drives early expression of Nkx3.2 in the gnathostome primary jaw joint
eLife 11:e75749.
https://doi.org/10.7554/eLife.75749