Fig 2: Sucrose profiles and global translation in patients 1, 2, and 3. Comparison of ribosomal fractions from fibroblasts of healthy control parent (A) and affected patient 1 (P1) heterozygous for the p.A178D RPL13 mutation (B). (C) Overlay of sucrose gradient profiles derived from healthy control parent (red line) and affected individual P1 (black line). Ribosomal subunit ratio (60S:80S) is indicated in the top right corner. (D) Box plot showing the results of the 60S:80S from fibroblasts of healthy controls (parents of patient 1) and mutation-positive (MP) subjects (patients 1 to 3 and mother of patient 2), each of whom is heterozygous for an RPL13 mutation. An increase median value for the 60S:80S ratio is observed in the MP group. The p value from one-sample t test comparing mutation-positive group against the median of the control samples. (E) Relative translation ratio from OP-puro translation assay. Global protein rate is reduced in the MP subjects (patients 1 to 3 and mother of patient 2) compared with controls (n = 2). Data are presented as mean and interquartile range. Experiments were performed twice, in duplicate or triplicate/donor/experiment, and analyzed using Mann–Whitney U test.

Fig 3: Evaluation of the rpl13 L191Afs/L191Afs zebrafish. (A) Sanger sequencing confirmed the genotype of the homozygous rpl13 L191Afs/L191Afs (rpl13 mutant) fish. (B) qPCR performed on RNA extracted from three WT and three rpl13 L191Afs / L191Afs reveals a partial activation of nonsense-mediated mRNA decay in mutant zebrafish compared with WT. The p value from Student's t test. (C) WB performed on protein extracts from three WT and three rpl13 L191Afs / L191Afs fish reveals reduced expression of eL13 in mutant fish compared with WT. The p value from Student's t test. (D) Alcian blue staining of control (WT) and rpl13 L191Af / L191Afs zebrafish at 5 dpf, 7 dpf, and 1 mpf. Compared with WT, cartilage deformities in the head of the rpl13 L191Afs/L191Afs are evident. Red lines indicate the ceratohyal (CH) angle, yellow line indicates the distance between the tip of CH and the tip of Meckel's (MK) cartilages, and green line indicates the MK width. Scale bar at 5 and 7 dpf = 200 µm; scale bar at 1 mpf = 500 µm. (E) At 5 dpf, rpl13 L191Af/L191Afs show a significant increase of the CH angle compared with WT and (F) at 7 dpf a significant reduction in the distance between the CH and MK cartilages is evident in mutants compared with WT. The p values from Student's t test. (G) At 1 mpf, the frequency of fish with short distance (<0.53 mm) between the CH and MK cartilages and (H) the frequency of fish with narrower MK cartilage were significantly higher in mutants compared with WT. The p values from chi-square test. (I) Standard length (SL) distribution is significantly different in WT compared with rpl13 mutants at 1 mpf, with a higher number of long fish (>6.5 mm) in the WT population compared with mutants. The p value from chi-square test. (J) A higher frequency of fish with high SL/HAA (height at anterior margin of anal fin) ratio is found in the mutant group compared with WT group. Morphometric parameters (panels G-J) were categorized by global quartiles, presented as number of fish (%) in each category, and frequencies analyzed by chi-square test. (K) Swim bladder inflation in control (WT) and rpl13 L191Af/L191Afs zebrafish at 7 dpf. Arrows indicate a partially inflated lobe in mutants compared with WT. Scale bar = 1 mm. Quantitative analysis shows a significant delay in the swim bladder inflation of the first lobe in rpl13 L191Af/L191Afs compared with WT. The p values from crosstabs with chi-square test for frequencies.

Fig 4: Expression of the eL13 protein in dermal fibroblasts from mutation-positive and control subjects. (A) Ribbon representation of eL13 bound to the human 60S subunit (based on pdb code 4V6X). Residues in eL13 targeted by the two missense mutations are highlighted in black boxes. The starting position of the affected region caused by the intronic mutation at N159 (p.N159_V160ins18) (black arrow) and the potentially disrupted interaction between Y161 and nucleotide A509 are highlighted. The start of the frameshift deletion harbored by our zebrafish model, p.L191Afs*32, which occurs at the point where the C-terminal a-helix of eL13 binds to the rRNA, is also marked (red arrow). Ribosomal proteins that interact with eL13 (red) include eL33 (blue), eL36 (purple), or eL27 (gray) are shown. Figure was generated using the PyMOL molecular graphics system, version 2.0.7 (Schrödinger, LLC). (B) Representative immunofluorescence images of dermal fibroblasts from patients and healthy controls showing staining for eL13 (column 1), structural dyes Hoechst and Phalloidin (column 2), and merged (column 3). Scale bar = 30 µm. (C) Western immunoblotting (WB) demonstrated expression of the eL13 protein in dermal fibroblasts from mutation-positive subjects at similar levels as in fibroblasts from control subjects. Acetylated tubulin was used as a loading control. (D) Densitometry analysis of the immunoblot confirmed the lack of significant differences in eL13 expression (mean ± SD, Student's t test). (E) Dermal fibroblasts from patients and controls showed similar growth kinetics over time in culture. Control 1 = male, age 29 years; control 2 = female age 30 years; control 3 = male age 17 years; control 4 = female age 31 years.

Fig 5: Genetic findings in the RPL13 gene. Genetic variants and pedigrees of the four families. RPL13 mutations were identified by WGS in family 1 (A), by Sanger sequencing in families 2 and 3 (B, C), and by WES followed by Sanger sequencing in family 4 (D). Schematic illustration of the RPL13 gene showing all the variants identified in this report (marked above the gene) and previously reported (below the gene) (E).