# A TALEN-Exon Skipping Design for a Bethlem Myopathy Model in Zebrafish.

Zlatko Radev, Jean-Michel Hermel, Yannick Elipot, Sandrine Bretaud, Sylvain Arnould, Philippe Duchateau, Florence Ruggiero, Jean-Stéphane Joly, Frédéric Sohm

► To cite this version:

Zlatko Radev, Jean-Michel Hermel, Yannick Elipot, Sandrine Bretaud, Sylvain Arnould, et al.. A TALEN-Exon Skipping Design for a Bethlem Myopathy Model in Zebrafish.. PLoS ONE, 2015, 10 (7), pp.e0133986. 10.1371/journal.pone.0133986 . hal-01240441

HAL Id: hal-01240441

https://hal.science/hal-01240441

Submitted on 28 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
RESEARCH ARTICLE

# A TALEN-Exon Skipping Design for a Bethlem Myopathy Model in Zebrafish

Zlatko Radev¹,²✉, Jean-Michel Hermel³✉*, Yannick Elipot⁴, Sandrine Bretaud⁵, Sylvain Arnould⁶, Philippe Duchateau⁶, Florence Ruggiero⁵, Jean-Stéphane Joly⁴, Frédéric Sohm¹,²*

1 UMS 1374, AMAGEN, INRA, Jouy en Josas, Domaine de Vilvert, France, 2 UMS 3504, AMAGEN, CNRS, Gif-sur-Yvette, France, 3 UMR 9197, INRA-CASBAH team, NEURO-Psi, CNRS, Gif sur Yvette, France, 4 UMR 9197, DECA team, NEURO-Psi, CNRS, Gif sur Yvette, France, 5 UMR 5242, Institut de Génomique Fonctionnelle de Lyon, ENS de Lyon, CNRS, Université Lyon 1, Lyon, France, 6 CELLECTIS SA, Paris, France

✉ These authors contributed equally to this work.

✉ Current Address: Department of Biology & Centre for Translational Medicine and Cell Therapy, Varna Medical University, Varna, Bulgaria

* sohm@inaf.cnrs.gif.fr (FS); jmhermel@inaf.cnrs-gif.fr (JMH)

<figure><img src="image_2.png" /></figure>

CrossMark
click for updates

<figure><img src="image_3.png" /></figure>

## OPEN ACCESS

**Citation:** Radev Z, Hermel J-M, Elipot Y, Bretaud S, Arnould S, Duchateau P, et al. (2015) A TALEN-Exon Skipping Design for a Bethlem Myopathy Model in Zebrafish. PLoS ONE 10(7): e0133986. doi:10.1371/journal.pone.0133986

**Editor:** Bruce B Riley, Texas A&M University, UNITED STATES

**Received:** June 1, 2015

**Accepted:** July 3, 2015

**Published:** July 29, 2015

**Copyright:** © 2015 Radev et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Data Availability Statement:** All relevant data are within the paper and its Supporting Information files. Additional information is available from the Neuro-PSI, Institut des Neurosciences Paris-Saclay UMR 9197, CNRS- Université Paris-Sud, Bât. 32, CNRS, 91190 Gif-sur-Yvette, France.

**Funding:** Cellectis SA supported the project in the form of a participation to ZR salary and to the running cost of the project and provided research material (i.e. the TALEN). SA and PD participated to the design of the TALEN and the technical choices related to gene targeting and to mutation detection. SA and PD did not participate to the choice of the

## Abstract

Presently, human collagen VI-related diseases such as Ullrich congenital muscular dystrophy (UCMD) and Bethlem myopathy (BM) remain incurable, emphasizing the need to unravel their etiology and improve their treatments. In UCMD, symptom onset occurs early, and both diseases aggravate with ageing. In zebrafish fry, morpholinos reproduced early UCMD and BM symptoms but did not allow to study the late phenotype. Here, we produced the first zebrafish line with the human mutation frequently found in collagen VI-related disorders such as UCMD and BM. We used a transcription activator-like effector nuclease (TALEN) to design the *col6a1*<sup>ama605003</sup>-line with a mutation within an essential splice donor site, in intron 14 of the *col6a1* gene, which provoke an in-frame skipping of exon 14 in the processed mRNA. This mutation at a splice donor site is the first example of a template-independent modification of splicing induced in zebrafish using a targetable nuclease. This technique is readily expandable to other organisms and can be instrumental in other disease studies. Histological and ultrastructural analyzes of homozygous and heterozygous mutant fry and 3 months post-fertilization (mpf) fish revealed co-dominantly inherited abnormal myofibers with disorganized myofibrils, enlarged sarcoplasmic reticulum, altered mitochondria and misaligned sarcomeres. Locomotion analyzes showed hypoxia-response behavior in 9 mpf *col6a1* mutant unseen in 3 mpf fish. These symptoms worsened with ageing as described in patients with collagen VI deficiency. Thus, the *col6a1*<sup>ama605003</sup>-line is the first adult zebrafish model of collagen VI-related diseases; it will be instrumental both for basic research and drug discovery assays focusing on this type of disorders.
gene to target nor did they participate to data
collection or the analysis made to characterize the
isolated mutants. ZR's PhD thesis was financed by
Physiologie Animale et Système d'Élevage (PHASE)
department of Institut National de la Recherche
Agronomique (INRA; contrat number NA; http://www.
phase.inra.fr), CELLECTIS SA (#13000408; http://
www.cellectis.com/) and Animaux Modèles
Aquatiques et GENétique (AMAGEN, #439549; http://
www.inaf.cnrs.gif.fr/neuro-psi/amagen.html); YE was
supported by the Fondation de la Recherche
Médicale (FRM; #ING201401129350; http://www.frm.
org). The electron microscopy facilities of the Imagif
Cell Biology Unit of the Gif campus (www.imagif.cnrs.
fr) is supported by the Conseil Général de l'Essonne
(http://www.essonne.fr). The authors acknowledge
the France-Biolmaging infrastructure (http://france-
bioimaging.org) supported by the French National
Research Agency (ANR-10-INSB-04-01; http://www.
agence-nationale-recherche.fr). This work has
benefited from the support of TEFOR—
Investissement d'avenir (#ANR-II-INSB-0014; http://
tefor.net/main/pages/main/). The funders had no role
in study design, data collection and analysis, decision
to publish, or preparation of the manuscript. The
funder provided support in the form of salaries for ZR,
but did not have any additional role in the study
design, data collection and analysis, decision to
publish, or preparation of the manuscript. The specific
roles of these authors are articulated in the 'author
contributions' section.

**Competing Interests:** The authors have read the journal's policy and the authors of this manuscript have the following competing interests: SA is at the head of the Manufacturing Department of CELLECTIS SA, PD is CSO and board member of CELLECTIS SA. Cellectis SA supported the project in the form of a participation to ZR salary and to the running cost of the project and provided research material (i.e. the TALEN). SA and PD participated to the design of the TALEN and the technical choices related to gene targeting and to mutation detection. SA and PD did not participate to the choice of the gene to target nor did they participate to data collection or the analysis made to characterize the isolated mutants. SA and PD did not participate to the writing of the manuscript. This does not alter the authors' adherence to PLOS ONE policies on sharing data and materials.

**Introduction**

Collagen type VI proteins are found in the extracellular matrix; they are associated with inter-
stitial type I–III collagen fibers, in connective tissues of skin [1] and skeletal muscle [2–4]. Col-
lagen VI is involved in cell-cell attachment and interact with collagen IV of the basement
membrane [5], thus playing an important role in muscle maintenance. The predominant form
of collagen VI is a heterotrimer made of three polypeptide chains: alpha 1, alpha 2 and alpha
3 that are encoded by the COL6A1, COL6A2 and COL6A3 genes respectively. Each alpha chain
contains a short triple helical collagenous domain, flanked by N- and C-terminal globular
domains of variable size [6]. In the functional collagen VI, the three alpha chains combine
intracellularly into triple helical monomer that assembles into antiparallel dimer, and then two
dimers form a tetramer. The tetramers are secreted into the extracellular matrix where they
assemble into beaded microfibrils, the functional unit of collagen VI [7–10]. In muscle, intersti-
tial fibroblasts are the main source of collagen VI which is secreted in the endomysium sur-
rounding the muscle fibers [3].

In humans, the deficiency of COL6 genes expression is associated with numerous diseases
[4,11–15]. Mutations in one of the three genes have been demonstrated to cause the onset of
muscular dystrophy [4,16,17,12,13]. Collagen VI-related disorders account for 7 to 19% of con-
genital muscular dystrophies as reported in several studies on cohorts of patients of different
ethnicities [18–20]. Bethlem myopathy (BM, MIM #158810) was first described in 1976 as a
benign myopathy with a childhood onset and a slow progression into adulthood [21]. BM rep-
resents the mildest form of the phenotypic spectrum of collagen VI-related disorders and Ull-
rich congenital muscular dystrophy (UCMD, MIM #254090) is at the most severe end of the
spectrum [4,13,14]. UCMD is characterized by clinico-pathological hallmarks such as kyphos-
coliosis, torticollis, follicular hyperkeratosis, excessive scar formation following skin trauma,
distal joint laxity, proximal joint contractures, protruding calcanei, scoliosis and early respira-
tory insufficiency [1,22,23]. BM and UCMD are rare diseases with reported incidences in
Northern England below 1:100,000 and 1:700,000 respectively [24]. Muscle biopsies of patients
with BM present characteristics of a myopathic rather than dystrophic histology i.e. muscle
fibers exhibit abnormally variable sizes, with an increase in the number of internal nuclei while
typically no fiber necrosis is found [2,25]. In UCMD, respiratory problems have been reported
at early disease onset, but this outcome remains infrequent in BM [4,26,27]. UCMD progresses
slowly and worsens with ageing. Ultimately more than half of the patients end up with severe
impairment of their walking ability [13,28]. Currently, BM and UCMD remain without cure.

BM and UCMD syndromes have been associated with mutations in one of the COL6A1, COL6A2 or COL6A3 genes [4,13]. Most of these mutations are substitutions that modify a splice site or create a premature stop codon. These mutations frequently disrupt the collagen protein structure by perturbing normal helix folding. In the case of BM, mutations frequently affect the correct assembly of collagen VI fibers or prevent the intracellular assembly and secretion of tetramers in the extracellular matrix, eventually diminishing the amount of functional collagen VI protein [4,13,29].

Mouse and zebrafish models reproducing mutations found in patients have contributed to a better understanding of the collagen VI-related diseases. In mouse lines, each of the three genes encoding collagen VI has been individually inactivated and the corresponding alterations in muscle were extensively studied in these mutants [30–35]. Similarly, zebrafish models with different forms of altered collagen VI chains were obtained by the injection of morpholinos [36,37]. However, the morphants only allowed the study of transient phenotypes in embryos up to 2 days post-fertilization (dpf). These studies have led to the identification of histopathological abnormalities and severe ultrastructural changes in skeletal muscles. In particular,
abnormal mitochondria and sarcoplasmic reticulum (SR) were observed as well as defective autophagy [36,38]. In addition, increased levels of apoptosis were detected in skeletal muscles of most animal models except in the *Col6a3* deficient mice [34]. Although potential treatment strategies have been tested on collagen VI-deficient mice as well as on zebrafish with a more severe UCMD-like phenotype, no pharmacological studies have been performed on zebrafish morphants modeling the milder BM-like conditions [36,39] and there were limited to very young fish.

In the present study, we used a transcription activator-like effector nuclease (TALEN [40]) to disrupt an essential splice site resulting in an in-frame skipping of exon 14 in the N-terminal portion of the collagenous domain of the Col6a1 protein chain. We reproduced in zebrafish a co-dominant mutation frequently found in BM patients that prevents the assembly of collagen VI dimers and tetramers and inhibits secretion of tetramers into the extracellular matrix [2, 4, 11, 29]. This mutation enabled us to generate and characterize the first mutant zebrafish line with altered collagen type VI alpha-1 chain.

The resulting *col6a1* mutants (*col6a1*<sup>ama605003</sup> hereafter and in ZFIN and ZIRC databases) represent the first model to study collagen VI-related disorders based on a stable zebrafish line and will give the opportunity to study the disease's mechanisms and progression from embryo to ageing fish. Moreover, the *col6a1*<sup>ama605003</sup> line will facilitate the search for a treatment for these disorders since zebrafish are particularly amenable to pharmacological screens.

# Materials and Methods

## Ethics Statement

Zebrafish of the strain TU were raised and bred according to standard procedures. The creation and maintenance of zebrafish with mutant collagen VI was approved by the local ethics committee. All anesthesia and lethal procedures were performed in tricaine (MS222) solutions as recommended for zebrafish, and all efforts were made to reduce suffering. We confirm that all ethical procedures presented in this study have been approved by the french Comité d'Étique en Expérimentation Animale n°59, Ministère de l'Enseignement Supérieur et de la Recherche (agreement number pending).

All chemicals were purchased from Sigma-Aldrich, unless specified otherwise.

## Generation of TALE nuclease coding sequences and molecular cloning

The steps performed to establish a zebrafish line with skipping of exon 14 in collagen VI are outlined in the flowchart (Fig 1). TALE nuclease (TALEN) coding sequences were obtained from CELLECTIS SA bioresearch. The DNA binding sites for the TALEN pair targeting *col6a1* were (5' to 3'): left site TAAAGGGTCACCAGGGC, right site TTCATTAGGTAACAGTA (Fig 2A). Coding sequences of TALE nucleases were subcloned in a derivative of the pSpe3-RfA vector for mRNA *in vitro* synthesis [41] between the 5' and 3' globin UTRs included in the vector. The integrity of the final constructs was verified by sequencing.

### *In vitro* mRNA synthesis

The final TALEN plasmids were linearized by digestion with *Kpn*I, the DNAs purified and used as template for *in vitro* mRNA synthesis by the mMESSAGE mMACHINE T3 Transcription Kit (Life Technologies). mRNAs were *in vitro* polyadenylated with the Poly(A) Tailing Kit (Life Technologies) and then purified using the NucleoSpin RNA Clean-up Kit (Macherey-Nagel).
Injection of zebrafish embryos

Zebrafish embryos at the one-cell stage were injected with three different doses of TALEN mRNAs (200-400-600 pg/pair), using a Picospritzer injector (Science Products). The toxicity of

<figure><img src="image_4.png" /><figcaption>**Fig 1.** Flowchart for the generation of TALEN-mediated targeted genomic modifications within the *col6a1* gene from the injection at the one-cell stage to the generation of the stable *col6a1*<sup>ama605003</sup> mutant zebrafish line. The steps performed to establish a zebrafish line with skipping of exon 14 in collagen VI are outlined.</figcaption></figure>

doi:10.1371/journal.pone.0133986.g001
the injected solutions was determined at 24 hpf by calculating the proportion of dead and malformed embryos for each dose.

## Lysis of embryos and fin clips

Normal zebrafish embryos of 2 dpf were euthanized with a lethal dose of tricaine (MS222) buffered with sodium bicarbonate, placed individually in 50 µl of lysis buffer with composition 10 mM Tris-HCl pH 8.0, 50 mM KCl, 0.3% Tween 20, 0.3% Igepal CA-630, 4 mM EDTA, supplemented with proteinase-K to 250 ng/µl just before use. Lysis was performed for 16 h at 55°C and proteinase-K was subsequently inactivated by incubation for 10 min at 95°C. Lysates were stored at -20°C. For fin clip lysates, adult fish (3 mpf) were anesthetized with 0.5 mg/mL tricaine buffered with sodium bicarbonate and biopsies from the caudal fin removed with a sharp blade. Fin clip biopsies were lysed in the same manner as embryos.

## PCR on embryo or fin clip lysates

PCRs encompassing the targeted region of *col6a1* were performed using the GoTaq Flexi DNA polymerase (Promega) in standard 50 µl reactions with 1 µl lysate per reaction as template. The primers used were as follows: forward primer TGCCACCATGAAGAAGAGTG and reverse primer TCAGATGTGAGTTGCTCAGAC (gp1, gp2, Fig 3A). The size of the amplified sequence was 246 bp (see Fig 3). Specificity and concentration of the PCR product was verified on a 2.5% agarose gel in 0.5 x TBE containing ethidium bromide at a final concentration of 0.5 µg/mL.

## T7 endonuclease I assay

For the T7 endonuclease I assay, 10 µl of non-purified PCR product (ca. 150 ng) were denatured and re-annealed in a thermocycler by incubation for 10 min at 95°C, 95°C-85°C at a rate

<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

**Fig 2.** Generation of somatic mutations in the *col6a1* gene in zebrafish. (A) Organization of exon-14 and intron-14 junction and positioning of the TALEN monomer binding sites (underlined). The splice donor site is boxed. (B) Toxicity of the TALEN injection in zebrafish. The proportion of dead, malformed and normal embryos at 2 dpf is shown for each injected concentration. The number of injected embryos is indicated in brackets. (C) T7 EI assay on a representative sample of the injected embryos at 2 dpf, for each condition. The proportion of positive and negative embryos is shown for each injected concentration. The number of injected embryos is indicated in brackets. The relative efficiency of the TALEN doses was not significantly different (p-value 0.1534 exact two-tailed Fisher's test).

doi:10.1371/journal.pone.0133986.g002
<figure><img src="image_10.png" /></figure>

<table><thead><tr><th rowspan="2">F0 Founder</th><th rowspan="2">F1 fin clips<br>genotyped by<br>direct<br>sequencing</th><th colspan="6">F1 fish heterozygous for the<br>corresponding mutation</th><th rowspan="2">Wild-type<br>F1 fish</th></tr><tr><th>M1</th><th>M2</th><th>M3</th><th>M4</th><th>M5</th><th>M6</th></tr></thead><tbody><tr><td>3</td><td>139</td><td>5</td><td>3</td><td>-</td><td>-</td><td>-</td><td>-</td><td>131</td></tr><tr><td>4</td><td>48</td><td>3</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>45</td></tr><tr><td>10</td><td>93</td><td>10</td><td>-</td><td>4</td><td>7</td><td>2</td><td>2</td><td>68</td></tr></tbody></table>

<table><thead><tr><th></th><th>EXON 14 col6a1</th><th>INTRON 14 col6a1</th></tr></thead><tbody><tr><td>WT:</td><td>GCTGTAAAGGGTCACCAGGGCCAGATGTGAGTATTACTGTTACCTAATGAAGTCA</td><td>col6a1 locus</td></tr><tr><td>ORM M1:</td><td>GCTGTAAAGGGTCACCAGGGCCAGATG--AGTATTACTGTTACCTAATGAAGTCA</td><td>(-2 bp)</td></tr><tr><td>M2:</td><td>GCTGTAAAGGGTCACCAGGGCCAG---GAGTATTACTGTTACCTAATGAAGTCA</td><td>(-4 bp)</td></tr><tr><td>*</td><td>GCTGTAAAGGGTCACCAGGGCCAGATGT--GTATTACTGTTACCTAATGAAGTCA</td><td>(-2 bp)</td></tr><tr><td>M4:</td><td>GCTGTAAAGGGTCACCAGGGCCAGA-GTGAGTATTACTGTTACCTAATGAAGTCA</td><td>(-1 bp)</td></tr><tr><td>M5:</td><td>GCTGTAAAGGGTCACCAGGGCCAGAT---AGTATTACTGTTACCTAATGAAGTCA</td><td>(-3 bp)</td></tr><tr><td>M6:</td><td>GCTGTAAAGGGTCACCAGGGCCA----------------GTTACCTAATGAAGTCA</td><td>(-15 bp)</td></tr></tbody></table>

<table><thead><tr><th></th><th>M1</th><th>M2</th><th>M3</th><th>M4</th><th>M5</th><th>M6</th><th>WT</th></tr></thead><tbody><tr><td>298 bp</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>244 bp</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

**Fig 3. Organization of the targeted locus and validation of the presence of alternatively spliced *col6a1* mRNA in the mutants.** (A) Organization of the genomic locus surrounding the *col6a1* exon 14. Scale in kb is shown on the left-hand side. Exons 10 to 19 are represented by vertical dashes. Introns are symbolized by a continuous line. Genotyping primers gp1 and gp2 are represented on top. (B) The mRNA structures of wild type (WT, 298 bp) and exon 14-skipped forms (244 bp) are represented. The scheme represents the mRNA with exon skipping of the mutant M2 (deleted bases are symbolized by stars). Exons are boxed. (C) Table of transmission of mutated alleles for the F0 founders 3, 4 and 10. The numbers of genotyped adult F1 fish, the number of times each different type of mutation (M1-M6) occurred and the counts for WT F1 fish are indicated. (D) Nucleotide sequence of the junction between exon and intron 14 and alignment of the mutations M1 to M6 targeted by the TALEN. Each deleted nucleotides is represented by a dash. The number of deleted bases is reported on the right-hand side of the alignment. M1 was an over-represented mutation (ORM) but was absent out of 5 males and 5 females wild type TU zebrafish (see S3 Fig). The *col6a1*<sup>ama605003</sup> line was generated from the M3 mutation (star). (E) Results of RT-PCRs on fin clips of heterozygous F1 fish containing one of the corresponding mutations (as indicated). The 298 bp band corresponds to the wild-type allele of the mRNA, the 244 bp band to the mutated allele. RT-PCR from a wild type (WT) fish is shown on the right hand side.

doi:10.1371/journal.pone.0133986.g003

of -2°C/s, 85°C-25°C at a rate of -0.3°C/s and then held at 4°C. The re-annealed DNAs were then treated with 5 U (0.5 µl) of T7 endonuclease I (New England Biolabs) upon addition of
1.5 µl buffer 2 (New England Biolabs) and 3.5 µl H₂O to a final volume of 15 µl. Reaction mixes were then incubated for 30 minutes at 37°C and resolved on agarose gels analogously to untreated PCR products. Mismatches at the targeted region of *col6a1* were expected to result in cleavage of heteroduplexes thus generating two additional fragments of ca. 136 bp and 110 bp.

## Confirmation of mutations by sequencing

To confirm mutations in the *col6a1* gene in mosaic F0 embryos or adult fish positive by T7 endonuclease I assay, PCRs on lysates encompassing the targeted site were cloned in the TA-pGEM-T Easy Vector (Promega) according to the manufacturer's protocol and transformed in One Shot TOP10 Chemically Competent bacteria (Life Technologies). Upon overnight blue-white antibiotic selection, PCRs were performed on white colonies using the GoTaq Flexi DNA polymerase (Promega) in standard 10 µl reactions. PCR products were purified using a Nucleo-Fast 96 PCR Plate and sent for sequencing by the reverse primer. Sequences were aligned to the reference sequence using the CLC Main Workbench program (Qiagen). Only insertions/deletions ≥ 2 bp were considered as mutations due to TALEN activity in order to eliminate possible errors generated in the amplification process.

For non-mosaic embryos/adult fish of generation F1 or higher, PCRs on lysates were directly purified and sequenced. If needed, TA cloning was performed as described above for sequences containing mixed traces, and sequences of individual alleles obtained.

## RT-PCR and detection of exon 14 skipping of *col6a1*

Total RNA extractions from embryos and fin clips were performed using the TRIzol reagent (Life Technologies) according to the manufacturer's protocol. Reverse transcription was performed on total RNA using MMLV transcriptase (Promega) and an oligo (dT) 15 primer (Promega) for 1 h at 42°C followed by inactivation of the enzyme for 5 min at 94°C. The cDNA was used as template in standard PCR reactions with GoTaq Flexi DNA polymerase (Promega). The primers used were: forward primer GGTCCAGTCGGTTACCAAGG (in exon 10, sp1, Fig 3B) and reverse primer CAGATGGTCCGTAGTTTCCAGG (in exon 16, sp2, Fig 3B) and the expected size of the PCR product was 298 bp for a WT allele and 244 bp for an allele with exon skipping (see Fig 3). PCR products were resolved on 2.5% agarose gels in 0.5x TBE containing ethidium bromide at a final concentration of 0.5 µg/mL and skipping of exon 14 was confirmed by sequencing.

## Histology and transmission electron microscopy

For hematoxylin-eosin-safran and Masson's trichrome (Novotec) stainings, muscle samples of 5 mpf fish were fixed in Bouin, embedded in paraffin, and 5 µm-thick sections were performed (Novotec, Lyon, France). Deparaffined sections were then stained with either one of the two coloration methods. For semi-thin sections and TEM procedures, fish at 2-day and 3-week post-fertilization were euthanized in lethal concentration of MS222 according to the animal care regulation, fixed in F1 (2% paraformaldehyde, 2.5% glutaraldehyde), embedded in agar 3%, vibratome-sectioned at 200 µm. Sections were post-fixed in 1% osmium plus 1% potassium ferricyanide, stained with uranyl-acetate en-bloc, dehydrated and embedded in epon (Embed-812, EMS). Three fish of each WT, HT and HM group of 3 mpf were immersed into lethal concentration of MS222 then perfused through the heart with F1 fixative. Two mm surgery blade transversal hand-sections of the whole body, posterior to the intestinal cavity, were embedded in agar 3% and processed as described above. Epon blocs were sectioned with an ultramicrotome (Leica EM UC6) equipped with a Jumbo-histo diamond knife (Diatome). Semi-thin sections (1 µm) were stained with Richardson's stain for histopathological analysis and imaged
with a Leica DMR microscope equipped with a Color View Soft Imaging camera. Ultrathin sections (80 nm) were prepared with the same diamond knife and collected on slim-bar (200 mesh) grid and 1 mm x 2 mm slot grid with poloform membrane (Agar Scientific) without further staining. Observations were performed using a transmission electron microscope (JEOL 1400). Images were acquired using a post-column high-resolution (11 megapixels) high-speed camera (SC1000 Orius; Gatan) and processed with Digital Micrograph (Gatan). In each fish, at least 3 myomeres were analyzed.

## TUNEL assay for detection of apoptosis

Fish at 3 week post fertilization (wpf) or 3 month post fertilization (mpf) were euthanized in lethal concentration of tricaine, fixed in 4% PFA and sagittal cryosections performed. TUNEL assay was performed with the DeadEnd Fluorometric TUNEL kit (Promega).

## Locomotion analysis

For the locomotion analysis, we used 3 mpf (n = 12, 16 and 11 for WT, HT and HM respectively) and 9 mpf fish (n = 15, 19 and 7 for WT, HT and HM respectively). Fish were individualized in a dedicated behavior room in tanks (8.8 x 11.8 x 4.6 cm for 3 mpf or 18 x 24 x 9.5 cm for 9 mpf). Tanks were placed on an infrared light box (Viewpoint, Lyon, France) and covered with a Plexiglas lid. Fish were placed, one per tank, in 200 mL (3 mpf) or 1 L (9 mpf) of fish facility water at 26°C and left to recover from stress overnight in the experimental setup. No physical, chemical or visual contact was possible between single fish. The light cycle in the behavior room was identical to the one of the fish facility (light switched on between 9 a.m. and 11 p.m.). On the next day, the spontaneous swimming activity of fish was recorded for one hour at noon with a DragonFly2 infrared digital camera (Point Grey Research, Richmond Canada) at 25 frames per second. Fish were tracked by the VirtualDub software (http://www.virtualdub.org) using a tracking plug-in developed by Yves Lhuillier. The total distance swum and the geographic coordinates of each fish in each frame were automatically determined. The distance swum between two consecutive frames was calculated from the geographical coordinates by triangulation. Instantaneous velocities were calculated for each second as the sum of the distances obtained every 25 frames and the maximum value of instantaneous speed was determined for each fish.

## Analysis of speed distribution

The total cumulated time (s) each fish swam at a given speed was calculated for each 1 cm/s increment from the geographical coordinates determined from the video analysis (see above and S7 Fig). The three speed activity profiles (SAP) were defined according to fish body length (bl.), with 3 mpf (about 2 cm) being half the length of 9 mpf fish (about 4 cm), as follows: rest SAP as the cumulated time the fish swam at a speed inferior to 1 cm/s for 3 mpf or 2 cm/s for 9 mpf, (speed < 0.5 bl./s); the moderate SAP as the cumulated time the fish swam at a speed comprise between 1 cm/s to 6 cm/s for 3 mpf or 2 to 12 cm/s for 9 mpf fish (0.5 bl./s < speed ≤ 3 bl./s); the fast SAP as the cumulated time the fish swam at a speed superior to 6 cm/s for 3 mpf or 12 cm/s for 9 mpf (speed ≥ 3 bl./s). The relative proportions of each of the 3 SAP for each genotype were calculated as a percentage for 3 mpf and 9 mpf fish.

## Occupancy of fish in the water column

Fish were conditioned one per tank overnight in tanks of 11.5 x 22.5 x 12.5 cm with 2 L of facility water and then video-recorded at 25 frames per second, in the vertical plane for 30 min.
The cumulated time each fish swam in the upper quarter of the tank was calculated from the videos (watched in accelerated mode at speed 4x). Vertical swimming trajectories were calculated from the geographical coordinates of fish tracked analogously to the locomotion analysis. For the diagrams, the position of the fish every 6 sec (1 frame every 150 frames) was plotted through the entire length of the films (30 min).

## Centrophobia analysis

To determine if fish at 9 mpf presented signs of centrophobia, we calculated the time swum by fish in the centre of the tank. For this, the centre was defined as the central horizontal zone, homothetic to the aquarium and covering 50% of the total surface of the tank. To realize this test, the films from the locomotion analysis (see above) were used.

## Respiration rate

Fish were conditioned individually overnight in tanks of 11.5 x 16.5 x 9 cm with 1.5 L of facility water. Then the respiration rate was determined by manually counting the oral/opercular movements of each fish for one minute in the fish facility to minimize stress conditions.

## Statistical analyses

For the behavioral studies on adult fish, data distributions were not tested to determine normal distribution, therefore they were analyzed through nonparametric tests. Chi-square tests were used to compare speed distributions and were performed with R. Mann-Whitney tests were used for all other pairwise comparisons and were performed with Statview. In the behavior analysis Figs (9 and 10), symbols *, ** and **** correspond to $p < 0.05$, $p < 0.01$ and $p < 0.0001$ respectively.

## Results

### Establishment of a stable zebrafish line using TALEN-mediated targeted mutagenesis of the exon 14 / intron 14 junction in the *col6a1* gene

In patients, skipping of exon 14 of the *col6a1* gene encoding the collagen type VI alpha-1 chain has been demonstrated to lead to the onset of BM, although a case of mild UCMD caused by this mutation has also been reported [42]. In zebrafish, exon 14 of the *col6a1* gene (genome assembly Zv9, NCBI Danio rerio annotation release 103), is homologous to exon 14 in humans. In zebrafish embryos, injection of morpholinos blocking the splice donor consensus sequence of intron 14 resulted in an in-frame skipping of exon 14 [36]. In the study of Telfer et al. [36], exon 14 is referred to as exon 13 since their report was based on Danio rerio genome assembly Zv8. In order to reproduce the skipping of exon 14 in a stable zebrafish line, we targeted the splice donor site of intron 14 with a TALEN directed against the exon/intron junction (Figs 1 and 2A). As natural polymorphism is a hindrance when trying to isolate NHEJ-induced mutations, we checked the targeted locus for polymorphism. The surrounding of the locus was found to be homozygous and to be of the same sequence as reported in the Zv9 genome in all fish tested (fin clip, PCR of 417 bp followed by direct sequencing of 5 males and 5 females non-injected TU, (data not shown). We injected three different doses of the TALEN mRNA pair (200, 400 and 600 pg/pair) at the zygote stage and counted normal, malformed and dead embryos (Fig 2B). The mutation rate as a function of the dose the TALEN mRNA was assessed by T7 endonuclease I (T7 EI) assay using the primers gp1 and gp2 (Figs 2C and 3A). The mutation rate of the three doses was not significantly different (exact two-tailed Fisher's test) ranging from 5% for the two lower doses to 25% for the highest dose. F0 fish injected with the
highest dose (600 pg) were raised to sexual maturity. Only normal embryos were retained for further studies and husbandry. We obtained 10 F0 fertile adult fish and analyzed fin biopsies of each of them for the presence of somatic mutations. The analysis was performed by sequencing of PCR products flanking the target site of the TALEN i.e. amplicons were individualized by TA cloning and 16 to 20 of them were sequenced (S2 Fig). We found that three fish carried somatic mutations at a rate of 5% or higher (fish # 1, 3 and 8; S2A Fig). Such high mosaicism reflects the cumulative effect of the fast development of zebrafish and the relatively high stability of the TALEN mRNA [43]. In these conditions, the TALEN induces stochastically new mutations in different cells at different times (S1 Fig). Indeed, mosaicism is not limited to somatic cells but is found also in the germline.

In the next step, we tested the progeny of each F0 fish for the presence of mutations. Briefly, each injected F0 fish was crossed with a wild type (WT) zebrafish. Of the resulting F1 eggs, 24 were individually lysed and served to amplify a DNA fragment encompassing the TALEN target site. The PCR fragments were then submitted to T7 EI assay. Positive samples were further analyzed by direct sequencing. The analysis of the F1 eggs showed that 7 F0 fish transmitted at least one mutation to their progeny (S2A Fig). One of the fish positive for the fin clip test failed to transmit any mutation. The number of positive embryos per clutch ranged from 1 to 11 out of 24. One over-represented mutation was present in the progeny of all 7 fish but absent from non-injected animals (S2 Fig). Three F0 fish transmitted more than one mutation to their progeny (fish 3, 4 and 10; S2 Fig). We chose to raise to adulthood only the progeny of those three fish.

In a second set of experiments, we obtained 280 F1 adult fish and analyzed a fin clip of each of them by direct sequencing of PCRs flanking the target site (gp1, gp2, Fig 3A) to identify heterozygous mutants (HT). Positive samples were thereafter TA-cloned and sequenced to determine the precise sequence of the mutation. Out of the 280 F1 fish, we isolated 36 adult heterozygous fish carrying one amongst six different mutations (small deletions M1-M6) at the targeted locus (Fig 3B, 3C and 3D). Five of the isolated mutations affected the splice donor consensus sequence at the beginning of intron 14: M1, M2, M3, M5, and M6. Two of those mutations also affected the end of exon 14: M2 and M6. Finally, one mutation, M4, did not affect the splice donor and consisted in a deletion of one nucleotide at the end of exon 14 introducing a frame-shift (Fig 3D).

## Correlation between isolated mutations in the *col6a1* gene and in-frame skipping of exon 14 in the spliced transcript

In order to confirm that the isolated mutations lead to exon skipping, we performed RT-PCR on total RNA from fin clip of F1 fish (Fig 3). We used primers sp1 and sp2 located within exon 10 and 16 of the *col6a1* gene (see Materials and methods and Fig 3B). The expected PCR product size was 298 bp for a wild type allele and 244 bp if exon 14 (54 bp long) was skipped (Fig 3B). In-frame exon skipping was confirmed for four mutations (M2, M3, M5 and M6) by gel electrophoresis of RT-PCR products (Fig 3E) and cDNA sequencing (data not shown). As expected, no exon skipping was observed for mutation M4 (see Fig 3E), however, cDNA sequencing confirmed that mutation M4 does indeed induce a frame-shift (data not shown). More surprisingly, we could not detect any exon skipping allele corresponding to mutation M1 (Fig 3E); similarly, sequencing did reveal no frame-shift either (data not shown). The mutation M1 was absent in the sequence of 10 wild type adults, 5 males and 5 females (S3 Fig). The most probable explanation for this result is nonsense-mediated mRNA decay [44] which is a well-established outcome of some splice site mutations and results in the selective degradation of the mutated mRNA.
<figure><img src="image_11.png" /></figure>

**Fig 4. Validation of the presence of *col6a1* mRNA by RT-PCR at different developmental stages in different organs of zebrafish.** (A) Confirmation of the presence of *col6a1* mRNA in different organs of adult WT zebrafish: *col6a1* mRNA is detected in the electrophoregram of organs from adult fish as indicated at the top of the panel. β-actin is used as internal control. (B) Confirmation of the presence of two forms of *col6a1* mRNA in wild-type (WT) and mutant (presenting exon 14 skipping) zebrafish embryos at 2 dpf from *col6a1*<sup>ama605003</sup> line. WT fish express only the wild type allele (298 bp); heterozygous mutant (HT) carry both wild type and mutated allele (244 bp) and homozygous fish (HM) express only the exon-skipped form.

doi:10.1371/journal.pone.0133986.g004

As we were in search for an exon-skipping mutation, we decided to establish the M3 muta-
tion as our working fish line (*col6a1*<sup>ama605003</sup>) as we had both male and female F1 heterozygous
fish available. This enabled us to generate directly F2 homozygous mutant (HM) fish. All
remaining experiments described in the paper were carried out with the *col6a1*<sup>ama605003</sup> fish
line.

col6a1 is expressed from embryo to adult in WT and zebrafish mutants

We performed RT-PCR on samples of different organs of adult WT fish that confirmed the
ubiquitous *col6a1* expression in intestine, caudal fin, jaw, spleen and in skeletal muscle (Fig
4A). The presence of the *col6a1* transcript in WT as well as mutant fish was confirmed by
RT-PCR analysis on whole zebrafish embryos at 2 dpf (Fig 4B).

*col6a1*<sup>ama605003</sup> embryos and fry were not morphologically different from
WT and fry displayed similar touch-evoked response

Mutant *col6a1*<sup>ama605003</sup> and WT sibling embryos and fry displayed neither gross phenotype dif-
ferences nor obvious behavioral changes (data not shown). Birefringence muscle imagery,
which is commonly used to identify dystrophic mutants [45,46], was not different in 3 and 4
dpf zebrafish mutants when compared to WT siblings (S4 Fig). Finally, actin staining with fluo-
rescent-conjugated phalloidin did not reveal significant structural alterations of muscle cells in
the trunk of 3 dpf mutants and WT siblings (S4 Fig). In order to assess locomotor activity of
the fry, the touch-evoked response was measured in 2 dpf mutant fry as previously described [46,47] and did not reveal significant differences between mutants and WT siblings (S5 Fig).

## Age-dependent progressive disorganization of myofibers of *col6a1*<sup>ama605003</sup> mutant fish without evidence of apoptosis

Histology of the trunk skeletal white muscle was evaluated by Richardson's staining from transversal and sagittal semi-thin sections. We studied HT and HM *col6a1*<sup>ama605003</sup> mutants and WT siblings at 2 days, 3 weeks and 4 months post-fertilization (2 dpf, 3 wpf, 4 mpf). At 2 dpf, we observed a mild disorganization of the skeletal muscle tissue in transverse sections of HT (Fig 5B1) and HM (Fig 5C1) *col6a1*<sup>ama605003</sup> fry as compared to WT (Fig 5A1): the vast majority of the myofibers were normal but a few of them already presented abnormal intracellular vacuoles. At 3 wpf, the phenotype worsened with noticeable abnormal intracellular vacuoles in the myofibers of mutants (red arrowheads, Fig 5B2-3 and 5C2-3). While the vacuoles were indeed present in myofibers of HT fish (Fig 5B2-3), they were clearly more numerous in the myofibers of HM fish (red arrowheads, Fig 5C2-5C3). At 4 mpf, myofibers with abnormal vacuoles persisted (Fig 5B4 and 5C4). We observed a relatively low number of altered myofibers in the representative sections we analyzed, which reflects the patchy distribution of abnormal myofibers amongst normal ones in mutant fish. In contrast, we were unable to find any abnormal vacuoles in myofibers in any muscle section from the WT fish we analyzed (Fig 5A1-4). Additionally, in *col6a1*<sup>ama605003</sup> mutants, cell-to-cell contact appeared to be weaker: we observed large gaps between pathological fibers and their neighbours in sagittal section (red arrowheads; Fig 5B2-5C2), although healthy myofibers were in close contact. This was probably due to a greater fragility of mutant myofibers towards the contraction-induced fixation artifact [48]. Generally, the muscular tissue of *col6a1*<sup>ama605003</sup> mutants appeared more fragile and the myofibers were not tethered to each other as well as in WT. The Richardson's staining appeared also weaker in the altered cells than in the healthy ones (Fig 5B and 5C). The presence of apoptosis was assessed by TUNEL on cryostat sections of trunk skeletal muscles at 3 wpf and 4 mpf on WT as well as on HT and HM *col6a1*<sup>ama605003</sup>. No obvious difference in TUNEL was found between WT and either *col6a1*<sup>ama605003</sup> mutant (data not shown). Finally, at 5 mpf, hematoxylin-eosin-safran or Masson's trichrome stainings of transversal paraffin sections of white muscles were performed to visualize cell nuclei and collagen content respectively. Indeed, *col6a1*<sup>ama605003</sup> mutant muscles showed a increase in the number of nuclei, most probably due to the presence of numerous fibroblasts, and associated abundant extracellular matrix indicating the development of fibrosis in damaged muscles (violet; Fig 6B, 6C, 6E and 6F). In HM, the increase in number of nuclei was amplified (arrows, Fig 6C) and was associated with the development of areas with accumulation of extracellular material (arrowheads, Fig 6C). Finally, we also observed in HM collagen-rich areas with numerous nuclei (star, Fig 6F) most probably due to fibrosis that was unseen in WT or in HT (Fig 6).

## Ultrastructural defects of muscle fibers in *col6a1*<sup>ama605003</sup> mutants were already present in embryos and worsened with age

In order to better characterize myofibre defects, we studied transverse and sagittal, ultra-thin sections of muscles from *col6a1*<sup>ama605003</sup> mutant and WT zebrafish, at 2 dpf, 3 wpf and 4 mpf (Figs 7-9) by transmission electron microscopy (TEM). These results confirm and complete those obtained by light microscopy (Figs 5 and 6) such as the patchy distribution of abnormal myofibers and the presence of abnormal intracellular vacuoles, and finally emphasize structural defects in mitochondria, myofibrils and the sarcoplasmic reticulum (SR).
<figure><img src="image_12.png" /><figcaption>**Fig 5.** Age-dependent progressive disorganization of muscle fibers of *col6a1*<sup>ama605003</sup> mutants. Light photomicrographs of Richardson's stained semi-thin 1-µm sections from wild type (WT, A1-4), heterozygous (HT, B1-4) and homozygous (HM, C1-4) *col6a1*<sup>ama605003</sup> mutants at 2 days, 3 weeks and 4 months post-fertilization (2 dpf, 3wpf, 4 mpf, respectively). Right from 2 dpf in HT (B1) as in HM (C1) mutants, we observed abnormal vacuoles (red arrowheads) in the centre of some of the muscle fibers. Muscle fibers adjacent to abnormal ones remained similar to the ones WT (A1). At 3 wpf and at 4 mpf, in HT (B2-4) as in HM (C2-4), the number and the size of abnormal vacuoles in myofibers varied from one area to another. Abnormal myofibers with vacuoles were scattered among unaffected myofibers identical to those observed in WT (A2-4). In sagittal sections of muscle from 3 wpf HT we observed breaks in fiber tethering (B2, red arrowhead) or cell in advanced degradation (C2, red arrowheads). Sections were cut according to transversal plane, except A2, B2, and C2 which were sagittal sections. Scale bars, 10 µm.</figcaption></figure>

doi:10.1371/journal.pone.0133986.g005
<figure><img src="image_13.png" /><figcaption>**Fig 6.** At 5 mpf, in *col6a1*<sup>ama605003</sup> mutants muscles, the number of nuclei was increased additionally HM muscle showed fibrosis. Light photomicrograph of 5 µm-thick paraffin section of trunk muscle of 5 mpf wild type (WT, A, D), heterozygous (HT, B, E) and homozygous (HM, C, F) *col6a1*<sup>ama605003</sup> mutants. The section were cut according to transversal plane and stained with hematoxylin-eosin-safran (A-C) or with Masson's trichrome (D-F) colorations. We observed a slight increase in nuclei in the fibrous septa of *col6a1*<sup>ama605003</sup> HT (B) and numerous nuclei in HM (C, arrows). Additionally, we observed unidentified amorphous material only in HM (C, arrowheads), and a figure of putative fibrosis (C, star). The Masson's trichrome confirmed the increasing number of nuclei (violet) from HT (E) to HM (F) as compared to WT (D). This staining also confirmed the presence of large collagen-rich areas most probably marking fibrosis (star, F) that was absent in WT (D) and HT (E). Scale bars 25µm.</figcaption></figure>

doi:10.1371/journal.pone.0133986.g006

In transverse ultrathin sections (Fig 7), an obvious disorganization of myofibrils was observed in the outer border of some myofibers, at all ages both in HT and HM *col6a1*<sup>ama605003</sup> (d, Fig 7B1-3 and 7C1-3 respectively). Large areas of disorganized myofibers were observed (d, Figs 7B and 7C and 8B and 8C). Some myofibers were completely devoid of myofibrils (asterisk, Fig 7C2) even though they were located beside normally organized ones. In these disorganized areas, the mitochondria appeared damaged (see below and arrows Fig 9) and most often bounded the limit between the Cohnheim's field of well-organized, crystal-like myofibrils and the disorganized ones (Figs 7B and 7C, 8B and 8C). In the same damaged areas, the SR was dilated and disorganized with globular enlargements. Swelling of the SR (arrowheads, Fig 7B1-2 and 7C1-3) was particularly obvious at the level of the terminal cisternae indicative of a possible defect in the excitation-contraction coupling. The section in myofibers of WT sibling presented consistently organized Cohnheim's field, mitochondria with finely defined cristae and regular SR (Fig 7A1-3).

In sagittal ultrathin sections (Fig 8), the *col6a1*<sup>ama605003</sup> mutants showed another striking defect. At all ages, the sarcomeres were misaligned both in HT (broken lines, Fig 8B1-3) and HM (broken lines, Fig 8C1-3). Within some myofibers, the Z-discs and M-lines of two adjacent sarcomeres were not aligned in *col6a1*<sup>ama605003</sup> mutants, contrary to those of WT zebrafish (Fig 8A1-3). Moreover, the cisternal dilations of the SR appeared more conspicuous in the sagittal
<figure><img src="image_14.png" /><figcaption>**Fig 7.** Disorganized Cohnheim's fields, swollen reticulum and altered mitochondria in *col6a1*<sup>ama605003</sup> mutant fish muscle. TEM pictures of transversal sections of muscle from wild type (WT, A1-3) and *col6a1*<sup>ama605003</sup> heterozygous (HT, B1-3) and homozygous (HM, C1-3) mutants at 2 dpf, 3 wpf and 4 mpf. At 2dpf (A1, B1, C1), 3 wpf (A2, B2, C2) and 4 mpf (A3, B3, C3) in HT and HM mutants, we observed in a few muscle fibers, abnormal mitochondria with areas of the matrix devoid of electron dense material (B1-3, C1-3, arrows). In these myofibers, the abnormal mitochondria were often located in close contact with enlarged sarcoplasmic reticulum (B1-2, C1-3, arrowheads). The altered mitochondria often marked the limit between crystal-like organized (o) Cohnheim's field and pathologic disorganized (d) myofibrils (B1-3, C1-2). In HM (C1), some fibers appeared more affected than in HT, with drastically disorganized (d) areas separated from organized myofibrils (o) with mitochondria having swollen cristae (C1-3, arrows). The sarcoplasm area of the more affected fibers appeared devoid of material (B3, C1-2, asterisks). (B3) In a few muscle fibers of HM, the reticulum appeared drastically dilated (asterisk) with numerous blebs or autophagic vacuoles localized nearby or within myofibrillar disarray foci (C2) that were also in close contact with mitochondria with enlarged cristae (C2-3, arrows). None abnormal mitochondrion, myofibril disorganization or enlarged reticulum was seen in WT (A1-3).</figcaption></figure>

doi:10.1371/journal.pone.0133986.g007

plane than in the transverse sections. In fact, the whole triads were completely disorganized in mutants with abnormal swelling of SR terminal cristae and T-tubules (arrowheads, Fig 8B and 8C). Mitochondria with swollen cristae were also found within the myofibrillar disarray foci (arrows, Fig 8B1 and 8C1). All the healthy mitochondria as seen in WT muscle, displayed properly delineated outer and inner membrane with well-disposed and defined cristae (Figs 7A1-2, 8A1 and 9A1-4) indicating well-preserved ultrastructure and excluding fixation artefact. These characteristics of mitochondria in good condition were absent in the myofibrillar disarray foci in HT as in HM *col6a1*<sup>ama605003</sup> mutants from all ages (arrows, Figs 7-9). Damage onset appeared to begin with a swelling of the mitochondria and the cristae (arrows, Fig 8B1
<figure><img src="image_15.png" /><figcaption>**Fig 8.** Misaligned sarcomeres in *col6a1*<sup>ama605003</sup> mutant fish muscle. TEM pictures of sagittal section of muscle from wild type (WT, A1-3) and *col6a1*<sup>ama605003</sup> heterozygous (HT, B1-3) and homozygous (HM, C1-3) mutants at 2 dpf (A1, B1, C1), 3 wpf (A2, B2, C2) and 4 mpf (A3, B3, C3). As shown in Fig 7 (transversal sections), we observed a swelling of the sarcoplasmic reticulum (B, C, arrowheads), the presence of abnormal mitochondria (B, C, arrows) and vacuoles with membrane blebbing of possibly autophagic nature (B and C, asterisks) in HT and HM mutants at all three ages. In the same myofibre, we also observed disorganized (d) regions amongst still well-organized (o) myofibrils. In sagittal section of HT (B) and HM (C) mutants muscles, the position of the Z-discs and M-bands revealed a conspicuous misalignment of adjacent sarcomeres (white lines). Finally, in HM at 3wpf and 4 mpf (C2-3), the actin filaments in some sarcomeres appeared detached from the Z-line where the vacuole/autophagic vesicles were present (asterisks). No ultrastructural abnormalities as the ones described above were seen in WT (A1-3).</figcaption></figure>

doi:10.1371/journal.pone.0133986.g008

and 8C1) then to evolve to a complete absence of matrix, with the detachment of the outer mitochondrial membrane from the inner membrane and ending with the hyper-condensation of the internal membrane (arrows, Fig 9B1-2). In what seemed to be an intermediate state, the mitochondria were slightly swollen and the cristae appeared blurred (arrows, Fig 9C1). In these areas, vesicles that resembled autophagic vacuoles or drastically enlarged T-tubules were frequently observed in close contact with altered mitochondria (arrowheads, Figs 7B3 and 7C1-2 and 8B1, 8B3 and 8C2-3). The altered mitochondria were most often located between organized and disorganized myofibrils and, for some of them, in direct contact with enlarged SR (arrows, Figs 7B1, 7B3 and 7C3, 8B2-3 and 8C2-3, 9B1, 9B3.2, 9C3 and 9C4). Finally, at 3 wpf,
<figure><img src="image_16.png" /><figcaption>**Fig 9.** Altered ultrastructure of mitochondria in *col6a1*<sup>ama605003</sup> mutant fish myofibers. TEM pictures of sagittal sections of muscles from wild type (WT, A1-4), heterozygous (HT, B1-4) and homozygous (HM, C1-4) *col6a1*<sup>ama605003</sup> mutants, at 2 dpf (A1, B1, C1), 3 wpf (A2, B2, C2), and 4 mpf (A3-4, B3-4, C3-4). In WT, the mitochondria constantly exhibited a dense and well-delineated lamination of the inner membrane and a well-defined outer membrane (A), attesting of the good quality of the fixation-inclusion procedures. In HT and HM muscles, the morphology of mitochondria was diversiform. Similarly to observations in Figs 7 and 8, some mitochondria presented a matrix partly cleared, or even devoid, of electron dense material (B, C arrow) and with severely dilated external membrane detached from the inner condensed membrane (B2-3, C3-4, arrows). These abnormalities were often associated with figures of vacuole/autophagic vesicles (C2, C4, asterisk). In the same area of HT or HM sections, we often observed a normal mitochondrion close to another one exhibiting either swollen cristae (B4, C4, arrows), or even cristae reduced to a few vesicles (B1-2, C2). In a few mitochondria the cristae were even ultra-condensed to a dense core separated from the outer membrane by a vesicle (arrows, B3.1-2, C1). The abnormal mitochondria most often bounded the limits between the crystal-like organized (o) and disorganized (d) myofibrils.</figcaption></figure>

doi:10.1371/journal.pone.0133986.g009
3 mpf

A
Total distance swum (m)

<figure><img src="image_17.png" /></figure>

<table><thead><tr><th>MW p-value</th><th>Total distance</th><th>Max. speed</th></tr></thead><tbody><tr><td>WT vs. HT</td><td>0.8527</td><td>0.425</td></tr><tr><td>WT vs. HM</td><td>0.9020</td><td>0.742</td></tr><tr><td>HT vs. HM</td><td>0.9214</td><td>0.132</td></tr></tbody></table>

B
Maximum speed (cm/s)

<figure><img src="image_18.png" /></figure>

9 mpf

C
Total distance swum (m)

<figure><img src="image_19.png" /></figure>

D
Maximum speed (cm/s)

<figure><img src="image_20.png" /></figure>

<table><thead><tr><th>MW p-value</th><th>Total distance</th><th>Max. speed</th></tr></thead><tbody><tr><td>WT vs. HT</td><td>0.1601</td><td>0.986</td></tr><tr><td>WT vs. HM</td><td><strong>0.0015</strong> **</td><td>0.827</td></tr><tr><td>HT vs. HM</td><td><strong>0.0077</strong> **</td><td>0.855</td></tr></tbody></table>

<figure><img src="image_21.png" /></figure>

<figure><img src="image_22.png" /></figure>

<figure><img src="image_23.png" /></figure>

<figure><img src="image_24.png" /></figure>

<figure><img src="image_25.png" /></figure>

<figure><img src="image_26.png" /></figure>

<figure><img src="image_27.png" /></figure>

<figure><img src="image_28.png" /></figure>

<figure><img src="image_29.png" /></figure>

<figure><img src="image_30.png" /></figure>

<figure><img src="image_31.png" /></figure>

<figure><img src="image_32.png" /></figure>

<table><thead><tr><th>Chi² p-value</th><th>p-value</th><th>p-value</th></tr></thead><tbody><tr><td>WT vs. HT 0.6315</td><td>WT vs. HM 0.4987</td><td>HT vs. HM 0.9014</td></tr></tbody></table>

<table><thead><tr><th>Chi² p-value</th><th>p-value</th><th>p-value</th></tr></thead><tbody><tr><td>WT vs. HT 0.3866</td><td>WT vs. HM <strong>9.4 x 10<sup>-11</sup></strong></td><td>HT vs. HM <strong>9.8 x 10<sup>-14</sup></strong></td></tr></tbody></table>

<figure><img src="image_33.png" /></figure>

<figure><img src="image_34.png" /></figure>

<figure><img src="image_35.png" /></figure>

<table><thead><tr><th>MW p-value</th><th>Rest</th><th>Middle</th><th>Fast</th></tr></thead><tbody><tr><td>WT vs. HT</td><td>0.9630</td><td>0.7103</td><td>0.7806</td></tr><tr><td>WT vs. HM</td><td>0.5796</td><td>0.3559</td><td>0.4237</td></tr><tr><td>HT vs. HM</td><td>0.5537</td><td>0.8823</td><td>0.8823</td></tr></tbody></table>

<figure><img src="image_36.png" /></figure>

<figure><img src="image_37.png" /></figure>

<figure><img src="image_38.png" /></figure>

<table><thead><tr><th>MW p-value</th><th>Rest</th><th>Middle</th><th>Fast</th></tr></thead><tbody><tr><td>WT vs. HT</td><td>0.1601</td><td>0.1934</td><td>0.6900</td></tr><tr><td>WT vs. HM</td><td><strong>0.0035</strong> **</td><td><strong>0.0035</strong> **</td><td>0.1405</td></tr><tr><td>HT vs. HM</td><td><strong>0.0182</strong> *</td><td><strong>0.0082</strong> **</td><td>0.4592</td></tr></tbody></table>

**Fig 10. The 9 mpf col6a1ama605003 HM fish swam around twice the distance of WT and HT fish.** We video-recorded free-swimming WT, HT and HM col6a1ama605003 fish at 3 mpf and 9 mpf for one hour in the horizontal plane, then we calculated by triangulation the total distance swum (A, C) and the maximum instantaneous speed (B, D). (A) At 3 mpf, there was no significant difference in the total distance swum between WT and HT, between WT and HM or between HT and HM. (C) At 9 mpf, there was a significant difference in the total swum distance between WT and HM (C, p-value 0.0015) and between HT and HM (C, p-value 0.0077), but not between WT and HT, indicating the late development of impairment in HM. There was no difference in the maximum speed between genotypes at the two ages studied (B, D). (E-J) We represented in pie-charts the speed activity profile (SAP) distributions in percentage of 3 classes and according to the fish bl. i.e. rest SAP (horizontal lines, speed slower than 1 cm/s for 3 mpf, E, and 2 cm/s for 9 mpf, H), middle range SAP
(squares, speed comprise between 1 to 6 cm/s for 3 mpf, F and 2 to 12 cm/s for 9 mpf, J) and fast SAP (open dots, speed superior to 6 cm/s for 3 mpf, G and to 12 cm/s for 9 mpf, J). This representation allowed us to show that at 9 mpf, there is a highly significant difference between WT versus HM (H vs. J, p-value 9.421x10⁻¹¹) and between HT and HM (I vs. J, 9.83x10⁻¹⁴). We further analyzed in a pairwise manner the histograms of the time (s) fish swum in the 3 SAP classes (speed distribution) described above for 3 mpf (K, L and M for rest, middle and fast SAP respectively) and 9 mpf (N, O and P for rest, middle and fast SAP respectively). This analysis showed that there was no difference in SAP within any of the 3 mpf groups. But the analysis showed that 9 mpf HM swam significantly more time in middle range SAP at the expense of their resting time i.e. in the rest SAP the difference is significant between WT and HM (N, p-value 0.0035); between HT and HM (N, p-value 0.0182). In the middle range SAP, the difference was significant between WT and HM (O, p-value 0.0035) and between HT and HM (O, p-value 0.0082). For 3 mpf fish, n = 12, 16 and 11 for WT, HT and HM respectively; for 9 mpf fish, n = 15, 19 and 7 for WT, HT and HM respectively. For each histogram and pie-chart, we performed either a Mann-Whitney (MW) or Chi-square test respectively, with * , ** and **** indicating p-values of <0.05, 0.01 and 0.0001 respectively.

doi:10.1371/journal.pone.0133986.g010

some myofibers presented sarcoplasmic blebs (S6A Fig) and a few ones showed the typical characteristics of dying cells (S6B Fig).

Altogether, our observations showed that muscle fibers of both HT and HM *col6a1<sup>ama605003</sup>* displayed disorganized myofibrils, enlarged SR, altered mitochondria morphology and misaligned sarcomeres. The defects started at embryonic stages and worsened as the fish aged. None of these alterations was observed in WT muscles. HM presented a more severe phenotype than HT as expected for a co-dominant mutation. Moreover, we found dying cells only in HM fish and only from 3 wpf onward. Most often, these alterations had a patchy pattern as they were observed in a few myofibers located within undamaged muscle regions. The presence of dying cells only in 3 wpf and 4 mpf HM fish emphasizes the progressive muscle alteration observed in the *col6a1<sup>ama605003</sup>* homozygous context.

## The *col6a1*<sup>ama605003</sup> mutants present symptoms of hypoxia at 9 months

Histopathological and TEM imagery showed that muscle tissue, myofibre ultrastructure and mitochondria were clearly altered in *col6a1*<sup>ama605003</sup> mutants, yet we did not observe any obvious locomotor impairment up to the age 3 mpf. Therefore we studied locomotion in older mutants at 9 mpf. HT, HM and WT siblings were placed individually into tanks and video-recorded in horizontal plane in a free swimming context. Briefly, after transfer into the experimental set-up, the fish were left overnight to recover from stress in a dedicated behavior laboratory (see *Materials and methods*). On the next day, the fish were video-recorded and the impact of their genotype on swimming behavior was analyzed.

At 3 mpf, no significant changes were observed between WT and *col6a1*<sup>ama605003</sup> mutant fish when the total distance swum or the maximum swimming speed (Fig 10A and 10B) was measured. All statistical analyses are indicated in the corresponding tables in Figs 10 and 11. At 9 mpf, maximal instantaneous speed remained similar between genotypes (Fig 10D). Strikingly, however, HM fish swam significantly longer distances, roughly double those swam by WT and HT (Fig 10C; p-value WT vs. HM, 0.0015; HT vs. HM, 0.0077, Mann-Whitney, MW). This result was rather unexpected since HM fish had worse myofibre defects than HT. Thus, we further analyzed the distribution of instantaneous speed from 1 to 21 cm/s of 3 and 9 mpf mutants and WT sibling fish (see S7 Fig). Since the average fish body length (bl.) of 3 mpf (about 2 cm) is the half of the bl. of the fish at 9 mpf (about 4 cm), we defined three classes of speed activity profile (SAP) according to the fish bl.: rest, middle range and fast SAP. Rest SAP was defined as speed below 0.5 bl./s, middle range SAP as speed ranging from 0.5 to 3 bl./s and fast SAP as speed greater than 3 bl./s (see *Materials and methods* and Fig 10E-10J).

Irrespective of their genotype, all 3 mpf fish groups displayed the same distribution between the three classes of SAP: fish swam roughly 70% of the time at speeds within the middle range and about 15% of the time at speeds within the rest or fast range respectively (Fig 10E-10G). Note that the distribution of instantaneous speed was different between 3 and 9 mpf fish (S7 Fig), thus SAP should only be used to compare free swimming behavior between fish of the
same age. At 9 mpf, HT and WT displayed similar SAP: they swam most of the time at speeds within the rest range (~70% of the time), most of the remaining time at speeds within the middle range (~25% of the time), and a little time in fast speed (~5% of the time; Fig 10H and 10I). However, 9 mpf HM fish displayed a strikingly different swimming behavior than their HT and WT siblings (Fig 10J). Thus, HM fish at 9 mpf rested much less (~5% of the time) and spent most of their time swimming at middle range speed (~75% of the time). This difference in SAP was highly significant between WT and HM (p-value 9.421 x 10⁻¹¹, Chi-square) and between HM and HT fish (p-value 9.83 x 10⁻¹⁴, Chi-square, Fig 10H–10J). Thus to better analyze this difference in SAP, we further compared in a pairwise manner, the corresponding histograms of swimming speed distributions for the three genotypes in each age group. At 3 mpf, there was no difference between mutant and control fish (Fig 10K–10M). But in the 9 mpf group, HM fish swam significantly more time at middle range activity than WT and HT fish (HM vs. WT p-value 0.082 and HM vs. HT p-value 0.0035, MW, Fig 10O) at the expense of time spent at rest speeds (HM vs. WT p-value 0.0182 and HM vs. HT p-value 0.0035, MW, Fig 10N). Since this type of increased activity or erratic behavior may indicate potential problems in oxygen intake [49] or suffering [50,51], we performed further studies on the respiration behavior of 9 mpf *col6a1*<sup>ama605003</sup> mutant fish.

To this end, we first tracked the trajectory in the vertical plane of single fish (see Materials and methods). We showed that *col6a1*<sup>ama605003</sup> mutant fish mostly swam in the upper part of the water column where oxygen is more abundant and defined by Kramer [52] as aquatic surface respiration (ASR; Fig 11A–11C). Hence, HT and HM mutant fish spent a significantly longer time in the 25% upper part of the tank than their WT siblings (HT vs. WT, p-value 0.0009, MW; HM vs. WT, p-value 0.0283, MW; Fig 11D). But unexpectedly, the cumulated times swum in the 25% upper part of the tank for HT fish was significantly longer (p-value 0.0472, MW) than the one of the HM (Fig 11D). Indeed, the main types of behavioral response of fish to decreased oxygen availability in water are activity profile changes, increase in respiration rate and surface respiration, and changes in vertical or horizontal habitat [49,53]. As a matter of fact, 9 mpf *col6a1*<sup>ama605003</sup> HM displayed a dramatic increase in respiration rate (Fig 11E) that is considered a direct response to hypoxia [53–55] i.e. the oral/opercular movement rate of HM fish was significantly higher than that of HT (p-value 0.0090) and WT (p-value 0.0090). There was no significant difference between the respiratory rate of WT and HT *col6a1*<sup>ama605003</sup> mutant fish. Since centrophobia is a stress index [56,57], we checked if the difference in respiratory behavior was due to stress by measuring the occupancy of the centre of the tank in the horizontal plane. We showed that there was no significant difference in occupancy of the centre of the tank between WT, HT and HM fish groups (Fig 11F). Therefore, we concluded that increased respiratory rates in HM were most probably a physiological response to hypoxia rather than to an increase in stress or anxiety.

Taken together, our data showed that as early as 2 dpf, *col6a1*<sup>ama605003</sup> mutant fish present alteration of muscle tissue with ultrastructural alteration of SR and of mitochondria. These histological defects progressed slowly in ageing 3 mpf fish, showing clear signs of muscle fibrosis at 5 mpf without any obvious behavioral changes at this age. Then, only at 9 mpf, *col6a1*<sup>ama605003</sup> fish displayed swimming behaviors typical of hypoxia. Interestingly, the observation of these simple physiological parameters in the behavior of free moving 9 mpf fish clearly discriminates between WT and *col6a1*<sup>ama605003</sup> mutants. Finally, the severity of the observed phenotype increased from HT to HM exemplifying the co-dominant character of the *col6a1*<sup>ama605003</sup> mutation.
Discussion

A TALEN designed for an exon skipping disease model

We generated several *col6a1* mutant zebrafish lines using a TALE nuclease. To the best of our knowledge, this is the first example in zebrafish which shows that template-independent NHEJ mediated repair can be used to produce an exon skipping mutation. In the zebrafish egg, the

<figure><img src="image_39.png" /></figure>

<figure><img src="image_40.png" /></figure>

<figure><img src="image_41.png" /></figure>

D
Cumulated time (s)
spent in the upper
quarter of the tank

<figure><img src="image_42.png" /></figure>

F Cumulated time (s)
in the half center area

<figure><img src="image_43.png" /></figure>

E
Respiratory rate (min⁻¹)

<figure><img src="image_44.png" /></figure>

**Fig 11.** The 9 mpf *col6a1*<sup>ama605003</sup> mutant fish presented a hypoxia behavior. We video-recorded for 30 min (25 frames/sec) and tracked the swimming trajectories in the vertical plane of WT, HT and HM *col6a1*<sup>ama605003</sup> fish at 9 mpf in a tank. The position of the fish every 6 sec (1 frame every 150 frames) was plotted. A representative trajectory of each genotype is presented, WT (A), HT (B) and HM (C). From these videos (n = 5 of each genotype) the cumulated time fish swam in the upper quarter of the tank (D) was calculated and was significantly different between WT and HT (p-value 0.009), WT and HM (p-value 0.0283) and HT and HM (p-value 0.0472). HT and HM mutants spent visibly more time in the oxygen-richer part of the tank. The respiratory rate (oral/ opercular movements, min⁻¹) was determined for mutants and controls, n = 5 (E) and this rate was significantly different between WT and HM (p-value 0.009) and HT and HM (p-value 0.009). As a test for centrophobia, we measured the time fish spent in the centre of the tank in the horizontal plane. The centre corresponds to the half centre area of the tank. There was no difference of occupancy between the three genotypes. n = 15, 19 and 7 for WT, HT and HM respectively. For each histogram, a Mann-Whitney test was performed; *, ** indicate p-values of < 0.05 and < 0.01 respectively.

doi:10.1371/journal.pone.0133986.g011
majority of generated double-strand breaks (DSB) are repaired by the error-prone NHEJ path-
way rather than by homology directed repair [58]. The use of a template-independent method is
thus expected to lead to high efficiency of splice site modification similar to NHEJ mediated gene
knock out. Short insertions or deletions generated within essential splice sites, are expected to
almost invariably modify the splicing. Indeed, the analysis of the various mutations isolated in
this study showed that four mutations out of six were located within the splice donor consensus
sequence of intron 14 of the gene *col6a1* and induced exon skipping. However, one mutation sit-
uated into the splice consensus sequence did not generate mRNA with skipping of exon 14 (M1).
There is no easy explanation why mutation M1 (Fig 3) failed to induce exon skipping while M3
(and M2, M5 and M6 for that matter) did so. M1 is clearly over-represented amongst the isolated
mutations. This mutation was absent of 10 non-injected wild type TU fish tested (S3 Fig). Simi-
larly, we failed to detect the mutated allele in fin clip of most of the founders (S2A Fig). Hence,
parental polymorphism seems unlikely. Sequencing artefact is unlikely as the fish transmitting
the M1 mutation were detected using T7 endonuclease assay prior sequencing. Therefore, it
seems the M1 mutation is actively favored in our fish. In fact, the total absence of mutated
mRNA in the M1 fish line may provide the best explanation for M1 over-representation. Indeed,
the absence of mutated mRNA is consistent with a mechanism of nonsense-mediated decay of
the mutated mRNA form [44]. The resulting phenotype, a reduced capability to produce Col6a1
protein, would be somehow milder than the structural defects introduced by the other mutations.
Thus, M1 mutation might avoid most of the counter-selection which hinders the deleterious co-
dominant mutations. Consequently, its apparent prevalence would increased.

Other outcomes of splice site mutations are possible such as the activation of alternative splice
sites upon deletion of the normally used ones. Thus, it is important to check that the isolated
mutation indeed induces the expected defect. Conversely, different mutations can produce the
same exon skipping as our results showed (see Fig 3D). Indeed, any mutation targeting one of the
two first nucleotides of the intron can induce the skipping of the preceding exon. Hence, target-
ing a splice site is as easy as realizing a knock-out with a TALE/CRISPR-Cas9 nuclease. Since
many human genetic diseases are caused by splicing defects, this simple and efficient method will
help to generate animal models for these diseases [44,59–61]. Application of this strategy may
also be useful to induce skipping of defective exons as therapy for other diseases [62,63].

Genotyping of injected animals and their progeny in order to find a heterozygous mutation
carrier may be a long and taxing procedure. Most frequently, only a few F0 and F1 animals are
analyzed. Mainly due to the mosaicism in F0 and limited sample size, mutations found in a F0
fish may be different from those found in its progeny (S1 Fig). In our study, we detected somatic
mutations in only three fish out of ten and from those three only two transmitted one or several
mutations to their progeny. Conversely, five fish which were negative for somatic mutations,
transmitted germline mutations to their progeny. In zebrafish, soma and germ-line cells are iso-
lated from each other early in development. Consequently, somatic and germ line nuclease-gen-
erated mutations derive most frequently than not from independent events. Hence, finding a
mutation in the fin clip of a F0 fish is hardly indicative of the presence of mutations in its germ
line. Direct screening of the F1 eggs should be preferred to screen for mutations.

Sarcoplasmic disarray and mitochondrial damage in *col6a1*<sup>ama60500</sup> muscles

The *col6a1*<sup>ama605003</sup> stable line we generated, allowed us to study at the functional and histolog-
ical level collagen VI-deficient fish up to 9 mpf for the first time. The *col6a1*<sup>ama605003</sup> zebrafish
mutants exhibited histological and ultrastructural muscle alterations identical to the ones
observed in muscle cell cultures [64–67] and biopsies [67,68] from patients with collagen VI
deficiency, in mice [30–35], as well as in the transient zebrafish model generated by morpholino-based knock-down [36,37].

The accumulation of extracellular matrix we observed in transversal section of white muscles (Fig 6) had the typical figure of fibrosis. The presence of fibrosis, essentially composed of collagen, a characteristic of all muscular dystrophies [69], is present in UCMD patient biopsies [23,67] and was also described in the *Col6a3* deficient mouse model [34]. At the ultrastructure level (Figs 7–9), the observed alterations included disorganized myofibrils, conspicuous abnormal dilation of the SR, terminal cisternae and T-tubules in the triads, abnormal mitochondria and misaligned adjacent sarcomeres. These ultrastructural muscle alterations match the clinical manifestations of BM [67,68,70] or UCMD [67,71] and were also observed in the muscles of the *mdx* mouse model of Duchenne muscular dystrophy [72] as well as in damaged muscle of exhausted athletes [73]. Moreover, the ultrastructural alterations we have observed from fry to young adult fish mutants are characteristic for the early ages of collagen VI-related myopathies in humans [67,68,70]. These findings suggest that our fish line is a good model to study the implication of mitochondria and SR in these diseases, which have not been studied extensively.

In our *col6a1*<sup>ama605003</sup> mutant line, we did not detect increased levels of apoptosis by TUNEL. Our data are in agreement with findings in other animal models with collagen VI deficiency [34], although in *col6a1*<sup>-/-</sup> mice, apoptosis was described as a hallmark of the model [31] and it was also present in the transient zebrafish morphants with exon skipping [36]. Nevertheless, we observed dying myofibers most frequently in mutants at 3 wpf and 4 mpf (S6 Fig). Our observation of the presence of autophagic vesicles associated with the aggregation and ultracondensation of mitochondria together with the absence of TUNEL-positive cells in mutant zebrafish suggest that *col6a1*<sup>ama605003</sup> deficient muscle cells do not generally die of apoptosis. However, further studies will be necessary to determine whether autophagy is altered in some way in the *col6a1*<sup>ama605003</sup> zebrafish mutants and could be stimulated as has been described in *Col6a1*<sup>-/-</sup> mouse models [31,38,74], since collagen VI can modulate autophagic signalling pathways [75]. Studies on animal models [36,38,76] and clinical trials on patients [65] strongly suggest that activation of autophagy and inhibition of the PTP of mitochondria should be considered as potential treatments for diseases caused by defects in collagen VI.

Moreover, we observed a general increase of muscle tissue alterations in HM compared to HT e.g. dying cells were observed at 3 dpf and 4 mpf in HM but not in HT and increased respiratory rate was observed only in 9 mpf HM. This allele-dependent gravity of the phenotype is consistent with a co-dominant inheritance and correlates with a protein alteration mechanism. This loss of function originates from perturbed assembly and secretion of collagen VI, which reduces the amount of functional protein in the extracellular matrix [4,33] and is probably responsible for the greater fragility of the muscle tissue to contraction-induced stress. The co-dominant mode of inheritance, together with late locomotor/hypoxia phenotype development is reminiscent of observations in patients with BM [4,17].

We also demonstrated an increase of phenotype severity with ageing at the histological ultrastructural and functional levels. In this regard, the *col6a1*<sup>ama605003</sup> zebrafish line allows the study of symptom progression and muscle alterations in ageing fish. We should emphasize the significance of a fish model which allows studying alterations in ageing muscle in analogy to health deterioration observed in some ageing patients with collagen VI deficiency [21,28,26,71].

## Hypoxia response behavior in aging *col6a1*<sup>ama605003</sup> fish, a possible link between collagen deficiency, mitochondria damage and respiration

The functional studies on the *col6a1*<sup>ama605003</sup> mutant line showed no alterations in normal touch-evoked response in 2 dpf fish fry (S5 Fig) and in swimming behavior of young 3 mpf
adults. With age progression, at 9 mpf, we observed three obvious changes in the behaviors of
mutants. First, compared to WT and HT, HM fish swam significantly more time in the range
of moderate SAP, spending much less time in the rest range (Fig 10 and S7 Fig). Secondly, com-
pared to WT and HT, HM fish had a significant increase in respiratory rate (Fig 11). Finally, 9
mpf HT and to a lesser extent HM fish swam preferentially in the upper part of the tanks
where oxygen, diffusing from air, is more concentrated (Fig 11). Moreover, beside their pre-
ferred aquatic surface respiration, the behaviors of HT and HM fish were different: the HT dra-
matically increased erratic swimming and the HM modified their swimming behavior slightly
but doubled their respiratory rate. These observations suggested a suffering behavior [50,51] or
more specifically hypoxia as they are consistent with reports of increased spontaneous locomo-
tor activity in fish exposed to oxygen deprivation [49,52,53,77]. This conclusion is further sup-
ported by the increase in respiratory rate observed in HM as this has been described for fish
lacking oxygen [49,55,78,79]. Finally, these behavior changes observed in *col6a1*<sup>ama605003</sup>
mutants were not induced by anxiety (Fig 11) since the mutant fish did not display thigmotaxis,
a centrophobic stress behavior [57]. Altogether these elements indicate the presence of a respi-
ratory adaptive response in *col6a1*<sup>ama605003</sup> mutant fish with a more pronounced phenotype in
HM. Indeed, while the change in swimming behavior seems to be sufficient for HT to compen-
sate for hypoxia, the more severely affected HM fish needed additionally to increase their respi-
ration rate. Since collagen VI is expressed in gills [80], the expression of the modified Col6a1
protein might provoke defects in the structure or the function of this organ resulting in an
altered intake of oxygen. Alternatively, the reason for mutants to suffer from hypoxia might be
attributed to the presence of degenerating mitochondria, which would interfere with the proper
cell respiration.

Moreover, interactions between hypoxia and collagen expression have been reported; for
review see Salminen et al. [81]. In normal conditions, collagen is stabilized by propyl and lysyl
hydroxylation hydrolases (PHD1-3) of the 2-oxoglutarate-dependent dioxygenases (2-OGDO)
family, the functions of which are oxygen dependent. PHD1-3-related, non-epigenetic path-
ways are linked to the hypoxia inducible factors (HIF) which are in turn targeted by PHD1-3.
Thus, oxygen controls these pathways through the PHDs that are oxygen sensors [81]. In nor-
mal normoxic conditions, PHDs are active and HIF are degraded, whereas in hypoxia PHD are
inactive and HIF accumulate inducing the expression of 300 genes involved in the response to
oxygen deprivation.

In turn, mitochondria play a key role in many diseases e.g. neurodegeneration, cancer, heart
and muscular diseases [82] to cite a few. Here, we showed that mitochondria are probably one of
the first organelles altered by the *col6a1* mutation in myofibers. Actually, on one hand the alter-
ation of mitochondria observed in many dystrophies has been linked to the Ca²⁺-overload due to
the SR calcium storage function alteration [83]. On the other hand, there is much evidence in the
literature that mitochondria control collagen gene expression through two interdependent epige-
netic and non-epigenetic mechanisms [81]. Both DNA and histone demethylases are 2-OGDO
enzymes. 2-oxoglutarate itself is a key metabolite in the Krebs cycle [81]. As an example, in
human fibroblasts, the *Col6A2* gene expression is reduced by CpG DNA methylation leading to
an alteration of the extracellular matrix [84]. Conversely, damaged extracellular matrix induces
the over-expression of DNMT3, one of the major proteins inducing DNA hypermethylation and
finely tunes the regulation of gene expression in rat smooth muscle cell [85]. Hence, it is rather
tempting to speculate that mitochondrial defects in *col6a1*<sup>ama605003</sup> line will create the condition
of a cross-talk between *col6a* genes expression and mitochondrial health. Then, in muscular dys-
trophy, collagen mutation might be expected to induce the alteration of SR and mitochondria,
which in turn will alter collagen expression thus amplifying in a vicious circle the initial mutation
effect. Such a process could account for the striking, patchy phenotype observed in histological
and TEM sections. Similarly, such process will amplify the mitochondrial alteration due to the ageing process [81] accounting at least partially for the late aggravation of the symptoms observed in many dystrophy diseases. Such pathway cross-talk could also be implicated in the age-related quality decrease of collagen, as it would trigger alterations of the matrix by activation of metalloproteinases [86].

## Conclusion

Our data show that the *col6a1*<sup>ama605003</sup> zebrafish line develops a mild and progressive form of collagen VI-related disease thus modeling perhaps best human BM. Adult zebrafish has become a widely used model organism for behavioral genetics [87] and for pharmacological studies due to the relatively low cost of maintenance as compared to mammals. Also, experiments on zebrafish embryos and fry allow the application of powerful statistical tests because of the simple production of high embryo numbers, which are easy to obtain and treat with a substance of interest by bathing. Data on long-term track record of BM patients remain scarce although the disease displays a great phenotypic variability [88]. Thus this first collagen VI mutant zebrafish line, modeling the condition in ageing animals might be used to study the pathological changes from as early as 2 dpf, to late 9 mpf, or even older fish. More generally, the *col6a1*<sup>ama605003</sup> line represents a suitable model for extensive drug tests to search for collagen VI-related disorder cures. We would gladly freely provide the *col6a1*<sup>ama605003</sup> line to any interested laboratory, please send your request to amagen@inaf.cnrs-gif.fr.

## S1 Materials and Methods

### Touch-evoked escape response

Dechorionated fish fry at 2 dpf were placed in the centre of a circle with a diameter of 10 mm and their behavior upon a gentle touch on the tail was recorded using a Zeiss axio zoom stereo-scope at one frame per 70 ms. The necessary time for each fry to exit the circle was measured.

### Birefringence

Fish at 3 dpf or 4 dpf raised in embryo medium with or without methylcellulose were anesthetized with 0.5 mg/mL tricaine (MS222) buffered with sodium bicarbonate. They were immobilized in 1% low-melting agarose on microscope slides and placed between two polarizers (3D lens, Taiwan), then photographed with a Zeiss axio zoom stereoscope or a Nikon AZ 100 zoom microscope.

### Fluorescent phalloidin staining of F-actin

Zebrafish fries at 3 dpf were euthanized in a lethal concentration of tricaine, fixed in 4% PFA overnight at 4°C, permeabilized in 5% Triton X-100/PBS followed by F-actin staining using phalloidin conjugated to rhodamine (Sigma) over night. Stained fries were mounted in 50% glycerol and observed under a Zeiss LSM780 confocal microscope.

## Supporting Information

S1 Fig. Origin of putative difference between mutations found in somatic and germline cells of zebrafish injected with TALENs. Upon injection of the TALEN-containing solution in fish embryos at the one-cell stage, different mutations are induced in some of the cells of the embryos (e. g. M1-M6) while no mutations arise in other cells (empty circles). The proportions of independently mutated and unmodified cells are variable and lead to the generation of mosaic fish. Lysis and genotyping of a sample of the injected embryos reveals the presence of
some of the mutations in the F0 embryos. Similarly, lysis and genotyping of a fin clip of adult F0 fish detects some of the mutations present in somatic cells of these fish. These somatic mutations may be the same or different as compared to the mutations present in the germline cells of the fish. Thus, lysis and genotyping of a sample of the progeny, F1 embryos, may reveal the presence of identical or different mutations as compared to the mutations found in somatic cells of F0 parents. Also, mutations in heterozygous non-mosaic F1 fish raised to adulthood for the establishment of fish lines to analyse (*) may differ or not from the mutations found in F1 embryos. On the whole, genotyping at each stage of the process often detects only a subset of the mutations present in the pool of available ones and can reveal the sequence of some of the induced mutations. Other mutations may remain undetected and may appear at later stages of the process of establishment of a fish line.

(TIF)

**S2 Fig. Validation of somatic and heritable mutations in the col6a1 gene in zebrafish. (A)** Number of mutations in the somatic cells of adult F0 fish and number of positive F1 fish found in a representative sample for each corresponding F0 fish. One of the mutations found in F1 embryos, M1 is overrepresented (ORM). Indicated are the numbers of F1 embryos containing the overrepresented mutation and the numbers of F1 embryos containing a different mutation. (B) Sequences of the five mutations found in the fin clips of F0 fish. The sequences of the over-represented mutation (F1 3-2) and the other three targeted mutations found in F1 embryos are also shown. The corresponding numbers of F0 and F1 fish are indicated at the left and the number of deleted base pairs is indicated at the right.

(TIF)

**S3 Fig. The M1 mutation is absent out of the sequence of 10 wild type adult.** The sequences of wild type adults, 5 females (top panel) and 5 males (bottom panel) are displayed together with their sequencing trace. The M1 over represented mutation, between vertical red lines (-2 bp, red box) is absent in these sequences. PCR primers used: (fwd) TCACTC CGTCTTCATT CAAAGATC, (rev) GTGATGGCAGCTTAAAGACG; sequencing primers (fwd) TCACTCCGTCTTCATTCAAAGATC.

(TIF)

**S4 Fig. No broad muscle disorganisation at 2 dpf in col6a1<sup>ama60500</sup> mutants.** Both wild type, HT and HM col6a1<sup>ama605003</sup> fish showed bright birefringence indicative of high muscle organi-sation. This was confirmed with phalloidin-rhodamine staining revealing actin myofilaments.

(TIF)

**S5 Fig. Similar touch-evoked response between WT and HT and HM col6a1<sup>ama605003</sup> mutants at 2 dpf.** To assess locomotor activity of the fry, the touch-evoked response was mea-sured in 2 dpf col6a1<sup>ama605003</sup> mutant fry video-recorded at 70 ms/frame. The time between the touch and escape of the centre circle of 1 cm diameter was measured. There was no significant difference in touch-evoked response between the 3 genotypes.

(TIF)

**S6 Fig. Dying cells in col6a1<sup>ama605003</sup> HM mutant at 5 mpf.** TEM pictures of col6a1<sup>ama605003</sup> HM mutant muscle showing (A) a dying myofibre with onset of cell membrane blebbing and disintegration. (B) Necrotic myofibre with ruptured plasma membrane and disorganized and degenerating organelles.

(TIF)

**S7 Fig. At 9 mpf, HM swam faster than WT and HT.** Distribution histogram of cumulated time (s) WT and HT and HM col6a1<sup>ama605003</sup> mutants swam at the indicated speed (S) in
increment of 1 cm/s. This distribution was calculated from the analysis of video-recorded for 1 h free-moving conditioned fish. The 9 mpf HM swam at higher speed for longer time than WT and HT. For 3 mpf, N = 12, 16 and 11 and for 9 mpf N = 15, 19 and 7 for WT, HT and HM respectively.

(TIF)

**S1 Supporting Information. The supporting information comprises a supporting Materials and Methods (below) and 6 supporting Figures (see S1-S7 Figs Captions below).**

(DOCX)

**Acknowledgments**

We would like to thank Emmanuel Dijoux, Céline Haye, Laurent Legendre and Benjamin Maurice for excellent animal care at the AMAGEN fish facility, and Sosthène Barbachou, Myriam Cousseau, Noémie de Crozé, Joanne Edouard, Emilie Guillon, Naïma El Khallouki, Emilie Menant and Yvette Monfils, for assistance and discussions. We are thankful to Valérie Allamand, Matthew Guille, Françoise Hermel and Sylvie Retaux for critical reading of the manuscript. We thank Yves Lhuillier for assistance with Virtualdub software (http://www.virtualdub.org). This work has benefited from the EM facilities and expertise of the Imagif Cell Biology Unit of the Gif campus (www.imagif.cnrs.fr).

**Author Contributions**

Conceived and designed the experiments: ZR JMH YE SB SA PD FR JSJ FS. Performed the experiments: ZR JMH YE SB. Analyzed the data: ZR JMH YE SB FR JSJ FS. Contributed reagents/materials/analysis tools: JMH YE SB SA PD FR JSJ FS. Wrote the paper: ZR JMH YE SB FR JSJ FS.

**References**

1. Kirschner J, Hausser I, Zou Y, Schreiber G, Christen H-J, Brown SC, et al. Ullrich congenital muscular dystrophy: Connective tissue abnormalities in the skin support overlap with Ehlers-Danlos syndromes. Am J Med Genet A. 2005; 132A: 296-301. doi: 10.1002/ajmg.a.30443 PMID: 15690374

2. Bertini E, Pepe G. Collagen type VI and related disorders: Bethlem myopathy and Ullrich scleroatonic muscular dystrophy. Eur J Paediatr Neurol EJPN Off J Eur Paediatr Neurol Soc. 2002; 6: 193-198.

3. Zou Y, Zhang R-Z, Sabatelli P, Chu M-L, Bönnemann CG. Muscle interstitial fibroblasts are the main source of collagen VI synthesis in skeletal muscle: implications for congenital muscular dystrophy types Ullrich and Bethlem. J Neuropathol Exp Neurol. 2008; 67: 144-154. doi: 10.1097/nen.0b013e3181634ef7 PMID: 18219255

4. Lampe AK, Bushby KMD. Collagen VI related muscle disorders. J Med Genet. 2005; 42: 673-685. doi: 10.1136/jmg.2002.002311 PMID: 16141002

5. Kuo H-J, Maslen CL, Keene DR, Glanville RW. Type VI Collagen Anchors Endothelial Basement Membranes by Interacting with Type IV Collagen. J Biol Chem. 1997; 272: 26522-26529. doi: 10.1074/jbc.272.42.26522 PMID: 9334230

6. Chu ML, Zhang RZ, Pan TC, Stokes D, Conway D, Kuo HJ, et al. Mosaic structure of globular domains in the human type VI collagen alpha 3 chain: similarity to von Willebrand factor, fibronectin, actin, salicylate proteins and aprotinin type protease inhibitors. EMBO J. 1990; 9: 385-393. PMID: 1689238

7. Furthmayr H, Wiedemann H, Timpl R, Odermatt E, Engel J. Electron-microscopical approach to a structural model of intima collagen. Biochem J. 1983; 211: 303-311. PMID: 6307276

8. Engvall E, Hessle H, Klier G. Molecular assembly, secretion, and matrix deposition of type VI collagen. J Cell Biol. 1986; 102: 703-710. PMID: 3456350

9. Baldock C, Sherratt MJ, Shuttleworth CA, Kielty CM. The supramolecular organization of collagen VI microfibrils. J Mol Biol. 2003; 330: 297-307. PMID: 12823969

10. Kuo HJ, Keene DR, Glanville RW. The macromolecular structure of type-VI collagen. Formation and stability of filaments. Eur J Biochem FEBS. 1995; 232: 364-372.
11. Kanagawa M, Toda T. The genetic and molecular basis of muscular dystrophy: roles of cell-matrix linkage in the pathogenesis. J Hum Genet. 2006; 51: 915-926. doi: 10.1007/s10038-006-0056-7 PMID: 16969582

12. Lisi M, Cohn R. Congenital muscular dystrophies: New aspects of an expanding group of disorders. Biochim Biophys Acta BBA-Mol Basis Dis. 2007; 1772: 159-172. doi: 10.1016/j.bbadis.2006.09.006

13. Bönnemann CG. The collagen VI-related myopathies: muscle meets its matrix. Nat Rev Neurol. 2011; 7: 379-390. doi: 10.1038/nrneurol.2011.81 PMID: 21691338

14. Bushby KMD, Collins J, Hicks D. Collagen type VI myopathies. Adv Exp Med Biol. 2014; 802: 185-199. doi: 10.1007/978-94-007-7893-1_12 PMID: 24443028

15. Hermanns-Lê T, Piérard GE, Piérard-Franchimont C, Delvenne P. Dermal Ultrastructure in Collagen VI Myopathy. Ultrastruct Pathol. 2014; 38: 164-166. doi: 10.3109/01913123.2013.829151 PMID: 24134684

16. Pan TC, Zhang RZ, Pericak-Vance MA, Tandan R, Fries T, Stajich JM, et al. Missense mutation in a von Willebrand factor type A domain of the alpha 3(VI) collagen gene (COL6A3) in a family with Bethlem myopathy. Hum Mol Genet. 1998; 7: 807-812. PMID: 9536084

17. Jöbsis GJ, Keizers H, Vreijling JP, de Visser M, Speer MC, Wolterman RA, et al. Type VI collagen mutations in Bethlem myopathy, an autosomal dominant myopathy with contractures. Nat Genet. 1996; 14: 113-115. doi: 10.1038/ng0996-113 PMID: 8782832

18. Okada M, Kawahara G, Noguchi S, Sugie K, Murayama K, Nonaka I, et al. Primary collagen VI deficiency is the second most common congenital muscular dystrophy in Japan. Neurology. 2007; 69: 1035-1042. doi: 10.1212/01.wnl.0000271387.10404.4e PMID: 17785673

19. Peat RA, Smith JM, Compton AG, Baker NL, Pace RA, Burkin DJ, et al. Diagnosis and etiology of congenital muscular dystrophy. Neurology. 2008; 71: 312-321. doi: 10.1212/01.wnl.0000284605.27654.5a PMID: 18160674

20. Clement EM, Feng L, Mein R, Sewry CA, Robb SA, Manzur AY, et al. Relative frequency of congenital muscular dystrophy subtypes: analysis of the UK diagnostic service 2001-2008. Neuromuscul Disord NMD. 2012; 22: 522-527. doi: 10.1016/j.nmd.2012.01.010 PMID: 22480491

21. Bethlem J, Wijngaarden GK. Benign myopathy, with autosomal dominant inheritance. A report on three pedigrees. Brain J Neurol. 1976; 99: 91-100.

22. Quijano-Roy S, Khirani S, Colella M, Ramirez A, Aloui S, Wehbi S, et al. Diaphragmatic dysfunction in Collagen VI myopathies. Neuromuscul Disord NMD. 2014; 24: 125-133. doi: 10.1016/j.nmd.2013.11.002 PMID: 24314752

23. Yonekawa T, Nishino I. Ullrich congenital muscular dystrophy: clinicopathological features, natural history and pathomechanism(s). J Neurol Neurosurg Psychiatry. 2014; doi: 10.1136/jnnp-2013-307052

24. Norwood FLM, Harling C, Chinnery PF, Eagle M, Bushby K, Straub V. Prevalence of genetic muscle disease in Northern England: in-depth analysis of a muscle clinic population. Brain J Neurol. 2009; 132: 3175-3186. doi: 10.1093/brain/awp236

25. Pepe G, de Visser M, Bertini E, Bushby K, Vanegas OC, Chu ML, et al. Bethlem myopathy (BETHLEM) 86th ENMC international workshop, 10-11 November 2000, Naarden, The Netherlands. Neuromuscul Disord NMD. 2002; 12: 296-305. PMID: 11801404

26. Haq RU, Speer MC, Chu ML, Tandan R. Respiratory muscle involvement in Bethlem myopathy. Neurology. 1999; 52: 174-176. PMID: 9921869

27. Van der Kooi AJ, de Voogt WG, Bertini E, Merlini L, Talim FB, Ben Yaou R, et al. Cardiac and pulmonary investigations in Bethlem myopathy. Arch Neurol. 2006; 63: 1617-1621. doi: 10.1001/archneur.63.11.1617 PMID: 17101832

28. Jöbsis GJ, Boers JM, Barth PG, de Visser M. Bethlem myopathy: a slowly progressive congenital muscular dystrophy with contractures. Brain J Neurol. 1999; 122 (Pt 4): 649-655.

29. Lampe AK, Zou Y, Sudano D, O'Brien KK, Hicks D, Laval SH, et al. Exon skipping mutations in collagen VI are common and are predictive for severity and inheritance. Hum Mutat. 2008; 29: 809-822. doi: 10.1002/humu.20704 PMID: 18366090

30. Bonaldo P, Braghetta P, Zanetti M, Piccolo S, Volpin D, Bressan GM. Collagen VI deficiency induces early onset myopathy in the mouse: an animal model for Bethlem myopathy. Hum Mol Genet. 1998; 7: 2135-2140. PMID: 9817932

31. Irwin WA, Bergamin N, Sabatelli P, Reggiani C, Megighian A, Merlini L, et al. Mitochondrial dysfunction and apoptosis in myopathic mice with collagen VI deficiency. Nat Genet. 2003; 35: 367-371. doi: 10.1038/ng1270 PMID: 14625552

32. Solares Perez A, Gartioux C, Beuvin M, Thao MV, Laine J, Medja F, et al. G.P.18 Muscle pathology and dysfunction in a novel mouse model of COLVI-myopathy. Neuromuscul Disord. 2012; 22: 827-828. doi: 10.1016/j.nmd.2012.06.088
33. Bernardi P, Bonaldo P. Mitochondrial Dysfunction and Defective Autophagy in the Pathogenesis of Collagen VI Muscular Dystrophies. Cold Spring Harb Perspect Biol. 2013; 5: a011387-a011387. doi: 10.1101/cshperspect.a011387 PMID: 23580791

34. Pan T-C, Zhang R-Z, Markova D, Arita M, Zhang Y, Bogdanovich S, et al. COL6A3 protein deficiency in mice leads to muscle and tendon defects similar to human collagen VI congenital muscular dystrophy. J Biol Chem. 2013; 288: 14320-14331. doi: 10.1074/jbc.M112.433078 PMID: 23564457

35. Pan T-C, Zhang R-Z, Arita M, Bogdanovich S, Adams SM, Gara SK, et al. A mouse model for dominant collagen VI disorders: heterozygous deletion of Col6a3 Exon 16. J Biol Chem. 2014; 289: 10293-10307. doi: 10.1074/jbc.M114.549311 PMID: 24563484

36. Telfer WR, Busta AS, Bonnemann CG, Feldman EL, Dowling JJ. Zebrafish models of collagen VI-related myopathies. Hum Mol Genet. 2010; 19: 2433-2444. doi: 10.1093/hmg/ddq126 PMID: 20338942

37. Ramanoudjame L, Rocancourt C, Lainé J, Lyphout L, Gartioux C, Schwartz ME, et al. G.P.19 Collagen VI genes in zebrafish skeletal muscle: Implications for collagen VI-myopathies. Neuromuscul Disord. 2012; 22: 828. doi: 10.1016/j.nmd.2012.06.089

38. Grumati P, Coletto L, Sabatelli P, Cescon M, Angelin A, Bertaggia E, et al. Autophagy is defective in collagen VI muscular dystrophies, and its reactivation rescues myofiber degeneration. Nat Med. 2010; 16: 1313-1320. doi: 10.1038/nm.2247 PMID: 21037586

39. Zulian A, Rizzo E, Schiavone M, Palma E, Tagliavini F, Blaauw B, et al. NIM811, a cyclophilin inhibitor without immunosuppressive activity, is beneficial in collagen VI congenital muscular dystrophy models. Hum Mol Genet. 2014; 23: 5353-5363. doi: 10.1093/hmg/ddu254 PMID: 24852368

40. Valton J, Cabaniols J-P, Galetto R, Delacote F, Duhamel M, Paris S, et al. Efficient strategies for TALEN-mediated genome editing in mammalian cell lines. Methods San Diego Calif. 2014; 69: 151-170. doi: 10.1016/j.ymeth.2014.06.013

41. Roure A, Rothbächer U, Robin F, Kalmar E, Ferone G, Lamy C, et al. A Multicassette Gateway Vector Set for High Throughput and Comparative Analyses in Ciona and Vertebrate Embryos. PLoS ONE. 2007; 2: e916. doi: 10.1371/journal.pone.0000916 PMID: 17878951

42. Pan T-C, Zhang R-Z, Sudano DG, Marie SK, Bönnemann CG, Chu M-L. New Molecular Mechanism for Ullrich Congenital Muscular Dystrophy: A Heterozygous In-Frame Deletion in the COL6A1 Gene Causes a Severe Phenotype. Am J Hum Genet. 2003; 73: 355-369. PMID: 12840783

43. Kimmel CB, Ballard WW, Kimmel SR, Ullmann B, Schilling TF. Stages of embryonic development of the zebrafish. Dev Dyn Off Publ Am Assoc Anat. 1995; 203: 253-310. doi: 10.1002/aja.1002030302

44. Carter MS, Li S, Wilkinson MF. A splicing-dependent regulatory mechanism that detects translation signals. EMBO J. 1996; 15: 5965-5975. PMID: 8918474

45. Gupta V, Kawahara G, Gundry SR, Chen AT, Lencer WI, Zhou Y, et al. The zebrafish dag1 mutant: a novel genetic model for dystroglycanopathies. Hum Mol Genet. 2011; 20: 1712-1725. doi: 10.1093/hmg/ddr047 PMID: 21296866

46. Charvet B, Ruggiero F, Le Guellec D. The development of the myotendinous junction. A review. Muscles Ligaments Tendons J. 2012; 2: 53. PMID: 23738275

47. Goody MF, Kelly MW, Reynolds CJ, Khalil A, Crawford BD, Henry CA. NAD+ biosynthesis ameliorates a zebrafish model of muscular dystrophy. PLoS Biol. 2012; 10: e1001409. doi: 10.1371/journal.pbio.1001409 PMID: 23109907

48. Carlsen F, Knappeis GG, Buchthal F. Ultrastructure of the resting and contracted striated muscle fiber at different degrees of stretch. J Biophys Biochem Cytol. 1961; 11: 95-117. PMID: 13876626

49. Kramer DL. Dissolved oxygen and fish behavior. Environ Biol Fishes. 1987; 18: 81-92.

50. Reilly SC, Quinn JP, Cossins AR, Sneddon LU. Behavioural analysis of a nociceptive event in fish: Comparisons between three species demonstrate specific responses. Appl Anim Behav Sci. 2008; 114: 248-259. doi: 10.1016/j.applanim.2008.01.016

51. Kalueff AV, Gebhardt M, Stewart AM, Cachat JM, Brimmer M, Chawla JS, et al. Towards a Comprehensive Catalog of Zebrafish Behavior 1.0 and Beyond. Zebrafish. 2013; 10: 70-86. doi: 10.1089/zeb.2012.0861 PMID: 23590400

52. Kramer DL. The evolutionary ecology of respiratory mode in fishes: an analysis based on the costs of breathing. Environ Biol Fishes. 1983; 9: 145-158.

53. Abdallah SJ, Thomas BS, Jonz MG. Aquatic surface respiration and swimming behaviour in adult and developing zebrafish exposed to hypoxia. J Exp Biol. 2015; doi: 10.1242/jeb.116343

54. Jonz MG. Development of oxygen sensing in the gills of zebrafish. J Exp Biol. 2005; 208: 1537-1549. doi: 10.1242/jeb.01564 PMID: 15802677
55. Shakarchi K, Zachar PC, Jonz MG. Serotonergic and cholinergic elements of the hypoxic ventilatory response in developing zebrafish. J Exp Biol. 2013; 216: 869–880. doi: 10.1242/jeb.079657 PMID: 23155078

56. Gerlai R, Lahav M, Guo S, Rosenthal A. Drinks like a fish: zebra fish (Danio rerio) as a behavior genetic model to study alcohol effects. Pharmacol Biochem Behav. 2000; 67: 773–782. doi: 10.1016/S0091-3057(00)00422-6 PMID: 11166068

57. Maximino C, de Brito TM, da Silva Batista AW, Herculano AM, Morato S, Gouveia A. Measuring anxiety in zebrafish: A critical review. Behav Brain Res. 2010; 214: 157–171. doi: 10.1016/j.bbr.2010.05.031 PMID: 20510300

58. Auer TO, Duroure K, De Cian A, Concordet J-P, Del Bene F. Highly efficient CRISPR/Cas9-mediated knock-in in zebrafish by homology-independent DNA repair. Genome Res. 2014; 24: 142–153. doi: 10.1101/gr.161638.113 PMID: 24179142

59. Ahlborn LB, Steffensen AY, Jønson L, Djursby M, Nielsen FC, Gerdes A-M, et al. Identification of a breast cancer family double heterozygote for RAD51C and BRCA2 gene mutations. Fam Cancer. 2014; doi: 10.1007/s10689-014-9747-y

60. Gardner JC, Liew G, Quan Y-H, Ermetal B, Ueyama H, Davidson AE, et al. Three different cone opsin gene array mutational mechanisms with genotype-phenotype correlation and functional investigation of cone opsin variants. Hum Mutat. 2014; 35: 1354–1362. doi: 10.1002/humu.22679 PMID: 25168334

61. Nozu K, Vorechovsky I, Kaito H, Fu XJ, Nakanishi K, Hashimura Y, et al. X-Linked Alport Syndrome Caused by Splicing Mutations in COL4A5. Clin J Am Soc Nephrol CJASN. 2014; doi: 10.2215/CJN.04140414

62. Lu Q-L, Yokota T, Takeda S, Garcia L, Muntoni F, Partridge T. The status of exon skipping as a therapeutic approach to duchenne muscular dystrophy. Mol Ther J Am Soc Gene Ther. 2011; 19: 9–15. doi: 10.1038/mt.2010.219

63. Pichavant C, Aartsma-Rus A, Clemens PR, Davies KE, Dickson G, Takeda S, et al. Current status of pharmaceutical and genetic therapeutic approaches to treat DMD. Mol Ther J Am Soc Gene Ther. 2011; 19: 830–840. doi: 10.1038/mt.2011.59

64. Merlini L, Bernardi P. Therapy of collagen VI-related myopathies (Bethlem and Ullrich). Neurother J Am Soc Exp Neurother. 2008; 5: 613–618. doi: 10.1016/j.nurt.2008.08.004

65. Maraldi NM, Sabatelli P, Columbaro M, Zamparelli A, Manzoli FA, Bernardi P, et al. Collagen VI myopathies: from the animal model to the clinical trial. Adv Enzyme Regul. 2009; 49: 197–211. doi: 10.1016/j.advenzreg.2008.12.009 PMID: 19162063

66. Angelin A, Tiepolo T, Sabatelli P, Grumati P, Bergamin N, Golfieri C, et al. Mitochondrial dysfunction in the pathogenesis of Ullrich congenital muscular dystrophy and prospective therapy with cyclosporins. Proc Natl Acad Sci. 2007; 104: 991–996. PMID: 17215366

67. Tagliavini F, Sardone F, Squarzoni S, Maraldi NM, Merlini L, Faldini C, et al. Ultrastructural changes in muscle cells of patients with collagen VI-related myopathies. Muscles Ligaments Tendons J. 2013; 3: 281–286. PMID: 24596691

68. Zamurs LK, Idoate MA, Hanssen E, Gomez-Ibañez A, Pastor P, Lamandé SR. Aberrant Mitochondria in a Bethlem Myopathy Patient with a Homozygous Amino Acid Substitution That Destabilizes the Colagen VI α2(VI) Chain. J Biol Chem. 2015; 290: 4272–4281. doi: 10.1074/jbc.M114.632208 PMID: 25533456

69. Mann CJ, Perdiguero E, Kharraz Y, Aguilar S, Pessina P, Serrano AL, et al. Aberrant repair and fibrosis development in skeletal muscle. Skelet Muscle. 2011; 1: 21–21. doi: 10.1186/2044-5040-1-21 PMID: 21798099

70. Bernardi P, Bonaldo P. Dysfunction of mitochondria and sarcoplasmic reticulum in the pathogenesis of collagen VI muscular dystrophies. Ann N Y Acad Sci. 2008; 1147: 303–311. doi: 10.1196/annals.1427.009 PMID: 19076452

71. Petrini S, D'Amico A, Sale P, Lucarini L, Sabatelli P, Tessa A, et al. Ullrich myopathy phenotype with secondary ColVI defect identified by confocal imaging and electron microscopy analysis. Neuromuscul Disord. 2007; 17: 587–596. doi: 10.1016/j.nmd.2007.04.010 PMID: 17588753

72. Rouède D, Coumailleau P, Schaub E, Bellanger J-J, Blanchard-Desce M, Tiaho F. Myofibrillar misalignment correlated to triad disappearance of mdx mouse gastrocnemius muscle probed by SHG microscopy. Biomed Opt Express. 2014; 5: 858. doi: 10.1364/BOE.5.000858 PMID: 24688819

73. Grobler LA. Skeletal muscle pathology in endurance athletes with acquired training intolerance. Br J Sports Med. 2004; 38: 697–703. doi: 10.1136/bjsm.2003.006502 PMID: 15562162

74. Gattazzo F, Molon S, Morbidoni V, Braghetta P, Blaauw B, Urciuolo A, et al. Cyclosporin A Promotes in vivo Myogenic Response in Collagen VI-Deficient Myopathic Mice. Front Aging Neurosci. 2014; 6: 244. doi: 10.3389/fnagi.2014.00244 PMID: 25309428
75. Neill T, Schaefer L, Iozzo RV. Instructive roles of extracellular matrix on autophagy. Am J Pathol. 2014; 184: 2146-2153. doi: 10.1016/j.ajpath.2014.05.010 PMID: 24976620

76. Tiepolo T, Angelin A, Palma E, Sabatelli P, Merlini L, Nicolosi L, et al. The cyclophilin inhibitor Debio 025 normalizes mitochondrial function, muscle apoptosis and ultrastructural defects in Col6a1-/- myopathic mice. Br J Pharmacol. 2009; 157: 1045-1052. doi: 10.1111/j.1476-5381.2009.00316.x PMID: 19519726

77. Kramer DL. Aquatic surface respiration in the fishes of Panama: distribution in relation to risk of hypoxia. Environ Biol Fishes. 1983; 8: 49-54.

78. Jonz MG, Nurse CA. Ontogenesis of oxygen chemoreception in aquatic vertebrates. Respir Physiol Neurobiol. 2006; 154: 139-152. doi: 10.1016/j.resp.2006.01.004 PMID: 16488670

79. Mcneil DG, Closs GP. Behavioural responses of a south-east Australian floodplain fish community to gradual hypoxia. Freshw Biol. 2007; 52: 412-420. doi: 10.1111/j.1365-2427.2006.01705.x

80. Nolte H, Konzer A, Ruhs A, Jungblut B, Braun T, Krüger M. Global protein expression profiling of zebrafish organs based on in vivo incorporation of stable isotopes. J Proteome Res. 2014; 13: 2162-2174. doi: 10.1021/pr5000335 PMID: 24611545

81. Salminen A, Kauppinen A, Hiltunen M, Kaarniranta K. Krebs cycle intermediates regulate DNA and histone methylation: Epigenetic impact on the aging process. Ageing Res Rev. 2014; 16: 45-65. doi: 10.1016/j.arr.2014.05.004 PMID: 24910305

82. Russell AP, Foletta VC, Snow RJ, Wadley GD. Skeletal muscle mitochondria: A major player in exercise, health and disease. Biochim Biophys Acta BBA-Gen Subj. 2014; 1840: 1276-1284. doi: 10.1016/j.bbagen.2013.11.016

83. Blake DJ, Weir A, Newey SE, Davies KE. Function and genetics of dystrophin and dystrophin-related proteins in muscle. Physiol Rev. 2002; 82: 291-329. PMID: 11917091

84. Kopp MU, Winterhalter KH, Trueb B. DNA methylation accounts for the inhibition of collagen VI expression in transformed fibroblasts. Eur J Biochem. 1997; 249: 489-496. PMID: 9370358

85. Jiang J-X, Aitken KJ, Sotiropoulos C, Kirwan T, Panchal T, Zhang N, et al. Phenotypic Switching Induced by Damaged Matrix Is Associated with DNA Methyltransferase 3A (DNMT3A) Activity and Nuclear Localization in Smooth Muscle Cells (SMC). Pasterkamp G, editor. PLoS ONE. 2013; 8: e69089. doi: 10.1371/journal.pone.0069089 PMID: 24282625

86. Aitken KJ, Block G, Lorenzo A, Herz D, Sabha N, Dessouki O, et al. Mechanotransduction of Extracellular Signal-Regulated Kinases 1 and 2 Mitogen-Activated Protein Kinase Activity in Smooth Muscle Is Dependent on the Extracellular Matrix and Regulated by Matrix Metalloproteinases. Am J Pathol. 2006; 169: 459-470. doi: 10.2353/ajpath.2006.050969 PMID: 16877348

87. Norton W, Bally-Cuif L. Adult zebrafish as a model organism for behavioural genetics. BMC Neurosci. 2010; 11: 90. doi: 10.1186/1471-2202-11-90 PMID: 20678210

88. Deconinck N, Richard P, Allamand V, Behin A, Lafôret P, Ferreiro A, et al. Bethlem myopathy: long-term follow-up identifies COL6 mutations predicting severe clinical evolution. J Neurol Neurosurg Psychiatry. 2014; doi: 10.1136/jnnp-2013-307245