Genetics: Published Articles Ahead of Print, published on June 22, 2009 as 10.1534/genetics.109.104695

TITLE

FORK STALLING AND TEMPLATE SWITCHING (FOSTES)

AS A MECHANISM FOR POLY-ALANINE TRACT

EXPANSION AFFECTING THE DYC MUTANT OF HOXD13, A

NEW MURINE MODEL OF

Olivier Cocquempot1,2,*, Véronique Brault1,2, Charles Babinet3 and Yann Herault1,2,4,#

Affiliations

1 Université d’Orléans, UMR6218, Molecular Immunology and Embryology, Orléans, France

2 CNRS, UMR6218, MIE, 3B rue de la Férollerie, 45071 Orleans cedex 2, France

3 Unité de Biologie du Développement, URA 1960, CNRS, Institut Pasteur, 25 rue du Docteur Roux 75015 Paris, France

4 CNRS, UPS44, TAAM, Institut de Transgenose, Orléans, France;

* Present Address: Université de Limoges, Limoges, France

#Corresponding author: Yann Hérault

UMR6218, IEM, Uni Orléans, CNRS, UPS44, TAAM, Institut de Transgénose

3B rue de la Férollerie, 45071 Orléans cedex 2 France

Tel: +33 238 25 7976 Fax: +33 238 25 5450

email: [email protected]

1 ABSTRACT

Polyalanine expansion diseases are proposed to result from unequal cross-over of sister chromatids that increases the number of repeats. In this report we suggest an alternative mechanism we put forward while we investigated a new spontaneous mutant that we named “Dyc“ for “ in Y and Carpe“ phenotype. Phenotypic analysis revealed an abnormal limb patterning similar to that of the human inherited congenital disease synpolydactyly (SPD) and to the mouse mutant model Spdh. Both human SPD and mouse Spdh mutations affect the Hoxd13 gene within a 15-residue polyalanine encoding repeat in the first exon of the gene, leading to a dominant negative HOXD13. Genetic analysis of the Dyc mutant revealed a trinucleotide expansion in the polyalanine-encoding region of the Hoxd13 gene resulting in 7 alanine expansion. However, unlike the Spdh mutation, this expansion cannot result from a simple duplication of a short segment. Instead, we propose the Fork Stalling and Template Switching (FosTeS) described for generation of nonrecurrent genomic rearrangements as a possible mechanism for the Dyc polyalanine extension, as well as for other polyalanine expansions described in the literature and that could not be explained by unequal crossing-over.

2 INTRODUCTION

Highly conserved Hox genes encode transcription factors containing a homeobox and forming multimeric complexes with other specific DNA-binding partners to regulate the transcription of specific genes (MOENS and SELLERI 2006).

They play key roles in the control of positional information during body axis specification and limb patterning. In mammals, Hox genes are organized in 4 clusters

A to D in which they are expressed with a precise spatial and temporal pattern that reflects their genomic organisation, with 5’genes expressed late and in more posterior and distal regions (KMITA and DUBOULE 2003; DESCHAMPS and VAN

NES. 2005). In the appendices, only the more 5’ located genes (9-13) are expressed, with a predominant role of A- and D-cluster genes. Duplications and deletions analyses (HERAULT et al. 1998; KMITA et al. 2005) have revealed a complex network of interactions among Hox genes from the different clusters as well as a global regulation of genes within a cluster.

Despite their central role in vertebrate body patterning, human congenital abnormalities attributable to mutations in HOX genes are rare and have only mild effects (for review see (GOODMAN 2002)). Mutations in HOX genes have been mostly associated with some human congenital malformation syndromes of the developing limbs. Whereas mutations in HOXA13 cause Hand-Foot-Genitalia syndrome (HFGS,

OMIM#140000), characterized by short and halluces, and clinobrachydactyly of the second and fifth fingers (GOODMAN et al. 2000; WAROT et al. 1997), mutations in HOXD13 give rise to a wide spectrum of limb malformations with variable penetrance and expressivity. Missense mutations within the homeodomain of HOXD13 result in a loss of function that leads to type D (BDD, MIM

3 113200) or to brachydactyly type E (BDE, MIM 113300). Expansion (MURAGAKI et al. 1996; AKARSU et al. 1996; GOODMAN et al. 1997) and, in few cases, frameshift deletions (GOODMAN et al. 1998; CALABRESE et al. 2000) of the polyalanine tract in exon 1 of HOXD13 trigger synpolydactyly (SPD, MIM 186000), a semi-dominant limb malformation syndrome characterized by digit duplication and fusion in heterozygotes and a combination of , and shortening of the hands of feet in homozygotes (YUCEL et al. 2005; CARONIA et al. 2003).

A spontaneous mouse mutant of Hoxd13 was isolated by Johnson and collaborators (1998) that, in its homozygous state, phenotypically and molecularly modelled human SPD. The Spdh (synpolydactyly homolog) mutation is a 21 bp in- frame duplication within the polyalanine stretch of exon 1, expanding the number of alanines from 15 to 22 and corresponds to the major type of expansion found in human SPD. Spdh heterozygous mice have only a slight shortening of digits 2 and 5, whereas homozygotes present a severe shortening of all fingers from both forelimbs and hindlimbs, associated with a combination of syndactyly and polydactyly due to a substantial delay in the ossification of the limb bony elements. In addition, homozygous mice lack preputial glands and males are not fertile (BRUNEAU et al.

2001; ALBRECHT et al. 2002).

We identified a new spontaneous mutant originated in the BALB/c line that presented an alteration of the distal part of the limbs and was hence named Dyc for

“Digit in Y shape and Carpe”. This mutation appeared spontaneously in the last backcross of a transgenic line on the BALB/c genetic background, during the final intercross to generate mice homozygous for the transgene. The phenotype was similar to human synpolydactyly and to phenotypes observed in mutations affecting

Hox genes. Phenotypic and genetic characterization of Dyc mice demonstrated that

4 Dyc is a new allele of Hoxd13 affecting the polyalanine cluster. Even though the resulting mutant allele contains, like the Spdh mutant, a 7-alanine expansion, sequencing analysis revealed that the nucleotide expansion in the Dyc mutant is different from that of the Spdh. Spdh polyalanine expansion is a straightforward reduplication of a short segment of the imperfect trinucleotide repeat, whereas the

Dyc expansion seems to results from two smaller duplications. Unlike classical long unstable trinucleotide repeat expansions that have been described to result from polymerase slippage during DNA replication (see KOVTUN et al. 2004 for review), polyalanine expansions are short, imperfect and stable when transmitted from generation to generation and hence have been proposed to result from an unequal crossing-over between two mispaired wild-type alleles (WARREN 1997). However, the Dyc expansion, like a few other alanine tract expansions reported in the literature do not all fit with the model of unequal-crossing over (GOODMAN et al. 1997; DE

BAERE et al. 2003; TROCHET et al. 2005; ROBINSON et al. 2005). Until now, no mechanism other than replication slippage, which does not fit with the “nature” of the polyalanine repeats, was proposed (TROCHET et al. 2007). In the manuscript, we describe the phenotype and mutation of the Dyc allele of Hoxd13, apply the mechanism of Fork Stalling and Template Switching (FosTeS) (LEE et al. 2007), based on the switching forward or backward of the replication fork along the DNA template using microhomology sequences, to generate the Dyc expansion and propose FosTeS as a new mechanism for polyalanine extensions.

5 MATERIAL AND METHODS

Mouse Lines and complementation test

The Dyc spontaneous mutation was isolated during the intercross of a mutant line backcrossed at N10 on the BALB/c genetic background. C57BL/6J HoxD

(ZAKANY et al. 1997) mice were obtained from D. Duboule through the European

Mouse Mutant Archive (EMMA; www.emmanet.org). The lines were maintained in

Specific Pathogen Free conditions and all the experiments were carried in accordance with the European and French regulation. Complementation test was done in B6C hybrid background.

Skeletal preparations

Preparations of skeletons were realized for adults, new-borns and embryos, according to JOHNSON et al. (1998). Briefly, the adults and new-born are eviscerated, skinned and fixed in ethanol. Staining of the bones is obtained with alizarin red and of the cartilage (new borns) with alcian blue. The flesh are clarified in glycerol. The embryos are fixed in Bouin solution, stained with alcian blue and the flesh is clarified in methyl salicylate.

In situ hybridization

Whole-mount in situ hybridization was performed using digoxigenin-labeled probes for Hoxd13, Hoxd12, Hoxd11 and Hoxd10 genes as described previously (HERAULT et al. 1998) and embryos of stages 12.5 days post coïtum (dpc) obtained by time matings from Dyc/+ x Dyc/+. Embryos were collected, fixed in 4% paraformaldehyde at 4°C, dehydrated and stored in methanol at -20°C until use. After rehydratation,

6 bleaching with H2O2 and a pre-treatment with proteinase K (10µg/ml), the embryos were incubated over night at 70°C with the probes (1µg/ml). They were then washed several times, incubated with anti-DIG alkaline phosphatase antibody (Roche

Biochemicals) at the dilution of 1:5000 overnight at 4°C. After washing of the antibody, staining was performed with BM purple (Roche).

Exon sequencing of Hoxd13 to Hoxd10 genes

To screen for mutation in the coding sequence of Hoxd10-13 genes, specific primers were designed with the DsGene program for amplification of the exons of those genes by PCR (JOHNSON et al. 1998). Exons were amplified by classical PCR.

Then, both strands of the PCR fragments were sequenced using the Big Dye

Terminator Reaction Mix (Applera France) on an ABI310 sequencer. The sequence was analysed with the DSGene Software (Acceleris, UK).

7 RESULTS

Dyc mutants show digits malformations and agenesis of preputial glands

Mice with abnormal digits in fore- and hind-limbs feet were discovered in the BALB/c line during a backcross of a transgene to generate homozygous animals. In the progeny of this intercross, 71% of the animals (n=109 produced progeny) had a digit alteration, with 48% having a slight phenotype affecting digits 2 and 5, and 23% having a stronger phenotype with shortening of all the fingers. This was indicative of the appearance of a semi-dominant mutation with the group of animals bearing the mild phenotype being heterozygote and the group with the strong phenotype being homozygote and suggested that the parents were heterozygote. Looking back at the parents, they were indeed found to have the mild digit phenotype.

In heterozygous adult animals, skeletal analysis (Fig. 1A) revealed a shortening of digits 2 and 5. Homozygous animals exhibited a complete alteration of the distal part of the five digits and malformation of all the bones from the ankle (metacarpal, carpal, metatarsal and tarsal bones). These alterations were associated with a syndactyly for digits 4 and 5 and the appearance of an extra sixth postaxial digit (polydactyly). The hindlimbs had an additional strong deformity of the metatarsal bone of digit 1. This new mutant was called “Dyc” for “Doigt en Y et Carpe” (Y-shaped finger and carpe).

No other alteration of the was found (n=10 wt and 10 Dyc/Dyc ; 12 week old). Necropsy analyses (n= 3 wt and 4 Dyc/Dyc) did not reveal any other drastic organ malformation and a haematological analysis (n= 3 wt and 4 Dyc/Dyc;

12 week old) indicated a normal globular numeration. Homozygous Dyc males are not fertile. Analysis of the uro-genital apparatus of these males revealed that prepucial glands were absent. Hence, the Dyc mutation produces two major

8 alterations: a severe malformation of the autopod and the lack of preputial glands in males. These two phenotypes were similar to those of the Spdh mutant mouse, a model for human synpolydactyly characterized by an expansion of unstable repeats within the 5’-end of the Hoxd13 gene (JOHNSON et al. 1998).

Effect of the Dyc mutation on the ossification

We went further to check whether the origin of the phenotype of the Dyc mutant during bone and cartilage development was similar to the Spdh mutant. Alcian blue/alizarin red staining of the skeleton was performed at several time points during embryonic and postnatal development to compare cartilage and bone development between Dyc/Dyc mutants and wild-type littermates (Suppl. Fig.1). Alcian blue staining of autopods at E14.5 dpc reveals a developmental delay in individualisation of the fingers and cartilage formation in the Dyc/Dyc autopod. Supernumerary fingers are present. At 15.5 dpc, deformation at the level of the metatarsus and metacarpus are visible. At birth (0.5 days post partum or dpp), ossification centres are detected in all wild-type digits but not in Dyc/Dyc mutant forelimbs. Polydactily is accompanied by syndactily. At 4.5 dpp, while ossification of the phalanx is completed in wild-type autopods. in Dyc/Dyc autopods, the ossification centres location is not clearly defined and there are defects in segmentation of the finger elements, making it hard to distinguish metacarpus and phalanx.

Dyc mutation is an allele of Hoxd genes

Limb phenotypes associated to the Dyc mutation are more severe than those found the Hoxd13 mutation obtained by genetic inactivation (KMITA et al. 2002), but similar to the alterations found in mutants of the posterior Hoxd genes, controlling limb

9 patterning (ZAKANY et al. 1997, KMITA et al. 2002). Heterozygote HoxDDel3/+ mutant mice (ZAKANY et al. 1997) bearing a deletion that inactivates the expression of

Hoxd11, Hoxd12 and Hoxd13 genes simultaneously have a shortening of digits 2 and

5 like Dyc heterozygous animals (Fig. 1b vs Fig. 1a). HoxDDel3/Del3 homozygous for this deficiency showed reduction in the length of all the digits () with a rougher and stiffer aspect of the digits and polydactyly/syndactyly in the hindlimbs

(Fig. 1b). The Dyc mutation in its homozygous form seems hence to be less severe than the Del3 one. We carried out a complementation test by crossing Dyc animals with HoxDDel3 mutant mice. HoxDDel3/Dyc trans-heterozygous mice displayed digit alteration similar to the ones present in HoxDDel3/Del3 homozygotes (Fig. 1b). Both have a shortening of all fingers with a malformation of toe 1 and polydactyly in the hindlimbs. Interestingly synpolydactily is rarely observed in transheterozygote animals. These results indicate that the Dyc mutation is an allele of one of the genes from the HoxD complex.

Posterior genes from the HoxD do not have their expression perturbed by the

Dyc mutation

We checked whether the Dyc mutation affected the transcriptional control of posterior genes from the HoxD complex in the limb bud. Only the more 5’ located genes (9-13) are expressed in the limbs. They are expressed in overlapping domains in more or less anterior to posterior gradients, are involved in the proper patterning of proximo- distal elements and impact the development of the antero-posterior axis of the limbs.

Comparing the expression patterns of Hoxd13, Hoxd12, Hoxd11 and Hoxd10 genes in the limb buds of wild-type controls and Dyc/Dyc mutant embryos at 12.5 dpc using whole-mount in situ hybridization revealed identical pattern of expression between

10 mutant and control animals in the distal part of the autopod (Suppl. fig. 2). Additional hybridization of Hoxd13 at 11.5 dpc and 13.5 dpc and of Hoxa13 at 13.5 dpc did not reveal any abnormal expression of these markers whereas physical alteration is already visible (data not shown).

The Dyc mutation affects the Hoxd13 protein sequence

After having checked that expression patterns of the posterior genes from the HoxD complex were not perturbed, we looked for the presence of mutations within the coding sequences of those genes. Sequencing analysis of Hoxd10, Hoxd11 and

Hoxd12 did not show any mutation within those genes in Dyc mutants. However,

PCR amplification of the 5’ part of exon 1 of Hoxd13 indicated an increase in the size of the amplified segment compared to BALB/c and C57BL/6J controls (data not shown). Sequencing of this fragment revealed a duplicated region responsible for the increased size (Fig. 2a). This amplification corresponds to 21 base pairs containing mainly a trinucleotide repeats, coding for a series of 15 alanine residues in the wild- type control and for 22 alanine residues in the Dyc mutant. Hence, the Dyc mutation corresponds to polyalanine expansion similar to the Hoxd13Spdh mutation that creates a gain-of-function of HOXD13 protein. However, alignment between Spdh and Dyc polyalanine stretches (Fig. 2b) indicates that the nature of the expansion is different between the two allelic mutants. Whereas the Spdh expansion appears to be a straightforward reduplication of the nucleotide repeats corresponding to alanines 1 to

7, the Dyc expansion is more complex, with a duplication of nucleotide repeats coding for alanines 1 to 5 and the two first nucleotides of alanine 7, whereas a small duplication corresponding of the sequence AGCA shifted in front of the first duplication instead of being after as in the wild-type sequence.

11 DISCUSSION

The Dyc mutation is a semi-dominant mutation causing severe limb malformations that are strickingly similar to that of the Spdh mutant, a 21 bp in-frame duplication within the polyalanine-encoding region of the Hoxd13 gene (JOHNSON et al. 1998). Like Spdh/Spdh mutants, homozygous Dyc mice have delayed distal limb development and altered chondrogenesis resulting in reduced size of all the fingers associated with polydactyly and syndactyly. Heterozygous animals are less affected, with only a slight shortening of fingers 2 and 5, indicating that the Hoxd13 wild-type allele is able to partly compensate the effect of the Dyc allele. Like in the Spdh mutant, alterations found in Dyc homozygous mice are more severe than that of null mutants of Hoxd13 (DOLLE et al. 1993; DAVIS and CAPECCHI, 1996) and but are similar to those found in the loss-of-function of the 3 most posterior genes of the

HoxD complex (Hoxd13, Hoxd12 and Hoxd11) (ZAKANY and DUBOULE 1996).

However, normal pattern expression of those genes during limb development of

Dyc/Dyc embryos indicates that inactivation does not occur at the level of the transcription as it could have been expected from the transcriptional control loops between Hox genes that have been described in the literature (GAVALAS et al.

2003).

Mapping of the Dyc mutation revealed an expansion of the trinucleotide repeat in the upstream exon of Hoxd13 similar to that found in the Spdh spontaneous mutant that was discovered in a mouse colony with the B6C3Fe genetic background.

Both Dyc and Spdh alleles give an additional 7 alanines at the 15 alanines stretch, conferring a dominant-negative property to the HOXD13 protein, probably by interfering with the activity of the wild-type HOXD13 and other HOX proteins. This

12 suggests that the polyalanine tract is important either for interaction with other partners or for formation of a functionally important structure of the protein and that whatever the genetic background, there is a threshold length of polyalanine extension above which the extension becomes pathologic and digital anomalies appear. Polyalanine tract expansion is also the common mutation found in the typical form of human SPD (MURAGAKI et al. 1996; AKARASU et al. 1996; HORSNELL et al. 2006), with additional alanine residues varying from 7 to 14 (GOODMAN et al.

1997). Like in the mouse, a minimum of 7 additional alanine residues is needed in human HOXD13 for a phenotype to appear and the genetic background does not seem to influence the severity of the phenotype that is influenced by the length of the alanine expansion (GOODMAN et al. 1997). But human and mouse differ in their heterozygous and homozygous states. Whereas digits fusions and duplications are present in human heterozygotes, this phenotype is only visible in mouse homozygotes indicating a difference in molecular interplay between human and mouse.

Polyalanine tract expansions are found in at least nine disorders (AMIEL et al.

2004; BROWN and BROWN 2004), mostly in genes that are coding for transcription factors involved in developmental processes (KARLIN et al. 2002; LAVOIE et al.

2003). Unlike classical long unstable trinucleotide repeat expansions found in genes implicated in hereditary neurodegenerative diseases, the polyalanine expansions are imperfect, short and stable through meiosis and mitosis. This suggests a different mechanism of mutation from strand slippage occurring during DNA replication, repair or recombination in the other diseases of unstable repeat expansion (SINDEN et al.

2002; PEARSON et al. 2005; KAPLAN et al. 2007; KOVTUN and MCMURRAY

2008). Indeed, the expanded polyalanine tract in SPD is quite stable (no increase in

13 mutation size) when transmitted from one generation to another. In addition, polyalanine tracts are short, with the longest non-pathogenic tract being of 20 residues (LAVOIE et al. 2003), whereas tracts of approximately 35 perfect trinucleotide repeats are required for instability and expansion (EICHLER et al. 1994;

KUNST and WARREN 1994). WARREN (1997) suggested that polyalanine expansion found in the mouse Spdh allele and human SPD may have resulted from unequal crossing over in the HOXD13 gene during meiosis due to misalignment between triplet repeats. This mechanism could also explain expansions found in other polyalanine disorders (ARAI et al. 2007; ROBINSON et al. 2005). However, unlike the Spdh mutation, the complex duplication found in the Dyc mutant is not a straightforward reduplication and cannot be explained by unequal crossing-over.

Looking at the 7-alanine expansion tract (Fig. 2b), the additional six last alanines can correspond to a reduplication of alanines 1-6, but there is an additional GCA codon in front of those 6 alanines whose origin cannot be explained. We checked for codon polymorphism in the B6C3Fe and BALB/c lines that could explain this atypical insertion, but could not find any. In addition, even if this additional GCA codon is ignored, the 6-alanine expansion cannot result from simple unequal allelic homologous recombination. Other mutations have been reported in the literature

(GOODMAN et al. 1997; TROCHET et al. 2007; DE BAERE et al. 2003; ROBINSON et al. 2005), which do not fit the unequal allelic recombination model for mutational mechanism described by WARREN (1997). TROCHET and collaborators (2007) already refuted the unequal crossing-over mechanism as the explanation of polyalanine expansions. However, they only suggested that the replication slippage mechanism proposed for diseases with unstable triplet expansions was also applying to polyalanine extensions.

14 We therefore tried to find another mechanism that could lead to the Dyc expansion. Among the mechanisms that have kept our attention is the FoSTeS described by Lee et al. (2007). FoSTeS was proposed to explain some nonrecurrent chromosomal rearrangements (duplications, deletions…) causing genomic disorders and that could not be explained by the usual nonallelic homologous recombination

(NAHR) and nonhomologous end joining (NHEJ) (SHAW and LUPSKI 2004; see GU et al. 2008) for review of the three mechanisms). Performing a breakpoint sequence analysis of nonrecurrent duplications, they uncovered complex arrays of normal, duplicated and triplicated sequences at the junctions of the duplications and identified few base-pair microhomologies that could act as “bridges” for the DNA replication fork to skip along the chromosome to create these complex arrays. They proposed that the DNA fork could skip forward or backward along the chromosome using these short sequence homologies when encountering a complex genomic architecture or a

DNA lesion. Since the proposition of FoSTeS as a mechanism to explain the non- recurrent rearrangements in PMD patients, another genomic disorder have been analysed for submicroscopic rearrangements, revealing the potential role of FoSTeS in an increasing number of complex pathologic rearrangements (CARVALHO et al.

2009). We checked if this long-range template switching mechanism could be applied to the Dyc mutation and more generally to polyalanine extensions via short-range template skipping. We made multiple alignment of the wild-type Hoxd13 polyalanine tract sequence to the Dyc sequence at microhomologies to try to reconstitute the Dyc expansion (Fig. 2c). We could identify two microhomology regions of five base-pairs each that could have made the replication fork skipped two times forward along the wild-type polyalanine sequence to create the Dyc extension. We then listed the polyalanine expansions referenced in the literature that could not be explained by

15 unequal crossing-over (Table 1) and we tried the same microhomology alignment method between the expanded alleles and their normal control alleles in order to determine the process of nucleotide amplification. We encounter the same problem as GOODMAN et al. (1997) to interpret the 9-alanine expansion of pedigree P. This one can either be interpreted as a simple reduplication of alanines 1-9 via the mechanism of unequal crossing-over or of FoSTeS (Fig. 3) only if there is a polymorphism or a mutation in the triplet coding for alanine 8 (GCT to GCG). Both the

15-bp (+ 5 alanines) and 18-bp (+ 6 alanines) insertions found in PHOX2B can be interpreted as a skipping “backward” of the replication fork followed by a skipping

“forward”. The 21-bp insertion (+ 7 alanines) could be interpreted as a single skipping

“backward”. However, the microhomology found is only of a single nucleotide. The 15 alanines extension in FOXL2 can be recovered by three skipping “backward” of the replication fork along the normal polyalanine tract. Finally, the +7 alanine extension in the PABPN1 polyalanine can be the result of a first skipping “backward” at a 2-bp homology followed directly by a second skipping “backward” of three nucleotides.

Hence, FoSTeS was tested on seven atypical polyalanine extensions and allowed to reconstitute the pathologic expansion for at least five of them. We therefore propose FoSTeS as an alternative mutational mechanism to unequal crossing-over in the generation of polyalanine expansions. Our hypothesis is further supported by the fact that FoSTeS is supposed to be engendered, or at least facilitated, by a surrounding high-order genomic architecture that contains cruciform or other non-B structures (Lee et al., 2007). Such alternative DNA structures are facilitated by specific sequence motifs such as low copy repeats (LCR), symmetrical features and short, direct or inverted, repeats, and the polyalanine trinucleotide repeat sequence possesses such sequence motifs and symmetry elements. Looking

16 for specific features flanking the breakpoints of complex rearrangements occurring in the MECP2 region and proposed to result from FoSTeS, researchers found increased frequency of the sequences 5’-CTG-3’/5’-CAG-3’ (CARVALHO et al.

2009).They suggest that these motifs might represent a cis acting sequence that may be a recognition site for proteins involved in priming DNA replication in eukaryotes.

Interestingly, 5’-CTG-3’/5’-CAG-3’ are motifs that are also found around the breakpoints in the polyalanine expansion tracts (see figures 2c and 3).

Hence, after having been described as a mechanism for large genomic rearrangements, our findings that FoSTeS could also be implicated in monogenic traits via very short stretches of sequences duplications reveals that a variety of genomic rearrangements ranging from gross DNA changes of several kilobases or megabases to rearrangements involving only few nucleotides expansions within single genes could originate from the same molecular mechanism.

Literature Cited

AKARSU, A. N., I. STOILOV, E. YILMAZ, B. S. SAYLI and M. SARFARAZI, 1996 Genomic structure of HOXD13 gene: A nine polyalanine duplication causes synpolydactyly in two unrelated families. Human Molecular Genetics 5: 945-952.

ALBRECHT, A. N., G. C. SCHWABE, S. STRICKER, A. BODDRICH, E. E. WANKER et al., 2002 The synpolydactyly homolog (spdh) mutation in the mouse - a defect in patterning and growth of limb cartilage elements. Mechanisms of Development 112: 53-67.

AMIEL, J., B. LAUDIER, T. ATTIE-BITACH, H. TRANG, L. DE PONTUAL et al., 2003 Polyalanine expansion and frameshift mutations of the paired-like homeobox gene PHOX2B in congenital central hypoventilation syndrome. Nature Genetics 33: 459- 461.

AMIEL, J., D. TROCHET, M. CLEMENT-ZIZA, A. MUNNICH and S. LYONNET, 2004 Polyalanine expansions in human. Human Molecular Genetics 13: R235-R243.

ARAI, H., T. OTAGIRI, A. SASAKI, T. HASHIMOTO, K. UMETSU et al., 2007 De novo polyalanine expansion of PHOX2B in congenital central hypoventilation syndrome:

17 unequal sister chromatid exchange during paternal gametogenesis. Journal of Human Genetics 52: 921-925.

BROWN, L. Y., and S. A. BROWN, 2004 Alanine tracts: the expanding story of human illness and trinucleotide repeats. Trends in Genetics 20: 51-58.

BROWN, L. Y., S. E. HODGE, W. G. JOHNSON, S. G. GUY, J. S. NYE et al., 2001 Possible association of NTDs with a poly-histidine tract polymorphism in the ZIC2 gene. American Journal of Human Genetics 69: 546-546.

BRUNEAU, S., K. R. JOHNSON, M. YAMAMOTO, A. KUROIWA and D. DUBOULE, 2001 The mouse Hoxd13(spdh) mutation, a polyalanine expansion similar to human type II synpolydactyly (SPD), disrupts the function but not the expression of other Hoxd genes. Developmental Biology 237: 345-353.

CALABRESE, O., S. BIGONI, F. GUALANDI, C. TRABANELLI, G. CAMERA et al., 2000 A new mutation in HOXD13 associated with foot pre-postaxial polydactyly. European Journal of Human Genetics 8: 140.

CARONIA, G., F. R. GOODMAN, C. M. E. MCKEOWN, P. J. SCAMBLER and V. ZAPPAVIGNA, 2003 An 147L substitution in the HOXD13 homeodomain causes a novel human limb malformation by producing a selective loss of function. Development 130: 1701-1712.

CARVALHO, C.M.B., F. ZHANG, P. LIU, A. PATEL, T. SAHOO et al., 2009 Complex rearrangements in patients with duplications of MECP2 can occur by Fork Stalling and Template Switching. Human Molecular Genetics [Epub ahead of print].

CRISPONI, L., M. DEIANA, A. LOI, F. CHIAPPE, M. UDA et al., 2001 The putative forkhead transcription factor FOXL2 is mutated in blepharophimosis/ptosis/epicanthus inversus syndrome. Nature Genetics 27: 159-166.

DAVIS, A. P., and M. R. CAPECCHI, 1996 A mutational analysis of the 5' HoxD genes: Dissection of genetic interactions during limb development in the mouse. Development 122: 1175-1185.

DE BAERE, E., D. BEYSEN, C. OLEY, B. LORENZ, J. COCQUET et al., 2003 FOXL2 and BPES: Mutational hotspots, phenotypic variability, and revision of the genotype- phenotype correlation. American Journal of Human Genetics 72: 478-487.

DESCHAMPS, J., and J. VAN NES, 2005 Developmental regulation of the Hox genes during axial morphogenesis in the mouse. Development 132: 2931-2942.

DOLLE, P., A. DIERICH, M. LEMEUR, T. SCHIMMANG, B. SCHUHBAUR et al., 1993 DISRUPTION OF THE HOXD-13 GENE INDUCES LOCALIZED HETEROCHRONY LEADING TO MICE WITH NEOTENIC LIMBS. Cell 75: 431- 441.

18 EICHLER, E. E., J. J. A. HOLDEN, B. W. POPOVICH, A. L. REISS, K. SNOW et al., 1994 LENGTH OF UNINTERRUPTED CGG REPEATS DETERMINES INSTABILITY IN THE FMR1 GENE. Nature Genetics 8: 88-94.

GAVALAS, A., C. RUHRBERG, J. LIVET, C. E. HENDERSON and R. KRUMLAUF, 2003 Neuronal defects in the hindbrain of Hoxa1, Hoxb1 and Hoxb2 mutants reflect regulatory interactions among these Hox genes. Development 130: 5663-5679.

GOODMAN, F., M. L. GIOVANNUCCI-UZIELLI, C. HALL, W. REARDON, R. WINTER et al., 1998 Deletions in HOXD13 segregate with an identical, novel foot malformation in two unrelated families. American Journal of Human Genetics 63: 992-1000.

GOODMAN, F. R., 2003 Congenital abnormalities of body patterning: embryology revisited. Lancet 362: 651-662.

GOODMAN, F.R., 2002 Limb malformations and the human HOX genes. Am J Med Genet 112: 256-265.

GOODMAN, F. R., C. BACCHELLI, A. F. BRADY, L. A. BRUETON, J. P. FRYNS et al., 2000 Novel HOXA13 mutations and the phenotypic spectrum of hand-foot-genital syndrome. American Journal of Human Genetics 67: 197-202.

GOODMAN, F. R., S. MUNDLOS, Y. MURAGAKI, D. DONNAI, M. L. GIOVANNUCCIUZIELLI et al., 1997 Synpolydactyly phenotypes correlate with size of expansions in HOXD13 polyalanine tract. Proceedings of the National Academy of Sciences of the United States of America 94: 7458-7463.

GU, W., F. ZHANG and J. R. LUPSKI, 2008 Mechanisms for human genomic rearrangements. Pathogenetics 1: 4.

HERAULT, Y., M. RASSOULZADEGAN, F. CUZIN and D. DUBOULE, 1998 Engineering chromosomes in mice through targeted meiotic recombination (TAMERE). Nature Genetics 20: 381-384.

HORSNELL, K., M. ALI, S. MALIK, L. WILSON, C. HALL et al., 2006 Clinical phenotype associated with homozygosity for a HOXD13 7-residue polyalanine tract expansion. European Journal of 49: 396-401.

JOHNSON, K. R., H. O. SWEET, L. R. DONAHUE, P. WARD-BAILEY, R. T. BRONSON et al., 1998 A new spontaneous mouse mutation of Hoxd13 with a polyalanine expansion and phenotype similar to human synpolydactyly. Human Molecular Genetics 7: 1033-1038.

KAPLAN, S., S. ITZKOVITZ and E. SHAPIRO, 2007 A universal mechanism ties genotype to phenotype in trinucleotide diseases. Plos Computational Biology 3: 2291-2298.

KARLIN, S., L. BROCCHIERI, A. BERGMAN, J. MRAZEK and A. J. GENTLES, 2002 Amino acid runs in eukaryotic proteomes and disease associations. Proceedings of the National Academy of Sciences of the United States of America 99: 333-338.

19 KMITA, M., and D. DUBOULE, 2003 Organizing axes in time and space; 25 years of colinear tinkering. Science 301: 331-333.

KMITA, M., N. FRAUDEAU, Y. HERAULT and D. DUBOULE, 2002a Serial deletions and duplications suggest a mechanism for the collinearity of Hoxd genes in limbs. Nature 420: 145-150.

KMITA, M., B. TARCHINI, D. DUBOULE and Y. HERAULT, 2002b Evolutionary conserved sequences are required for the insulation of the vertebrate Hoxd complex in neural cells. Development 129: 5521-5528.

KMITA, M., B. TARCHINI, J. ZAKANY, M. LOGAN, C. J. TABIN et al., 2005 Early developmental arrest of mammalian limbs lacking HoxA/HoxD gene function. Nature 435: 1113-1116.

KOVTUN, I. V., and C. T. MCMURRAY, 2008 Features of trinucleotide repeat instability in vivo. Cell Research 18: 198-213.

KOVTUN, I. V., C. SPIRO and C. T. MCMURRAY, 2004 Triplet repeats and DNA repair - Germ cell and somatic cell instability in transgenic mice. Methods in Molecular Biology: 309-319.

KUNST, C. B., and S. T. WARREN, 1994 CRYPTIC AND POLAR VARIATION OF THE FRAGILE-X REPEAT COULD RESULT IN PREDISPOSING NORMAL ALLELES. Cell 77: 853-861.

LAUMONNIER, F., N. RONCE, B. C. J. HAMEL, P. THOMAS, J. LESPINASSE et al., 2002 Transcription factor SOX3 is involved in X-linked mental retardation with growth hormone deficiency. American Journal of Human Genetics 71: 1450-1455.

LAVOIE, H., F. DEBEANE, Q. D. TRINH, J. F. TURCOTTE, L. P. CORBEIL-GIRARD et al., 2003 Polymorphism, shared functions and convergent evolution of genes with sequences coding for polyalanine domains. Human Molecular Genetics 12: 2967- 2979.

LEE, J. A., C. M. B. CARVALHO and J. R. LUPSKI, 2007 A DNA replication mechanism for generating nonrecurrent rearrangements associated with genomic disorders. Cell 131: 1235-1247.

PROTEIN. Molecular and Cellular Biology 14: 4057-4066.

MOENS, C. B., and L. SELLERI, 2006 Hox cofactors in vertebrate development. Developmental Biology 291: 193-206.

MUNDLOS, S., F. OTTO, C. MUNDLOS, J. B. MULLIKEN, A. S. AYLSWORTH et al., 1997 Mutations involving the transcription factor CBFA1 cause cleidocranial dysplasia. Cell 89: 773-779.

20 MURAGAKI, Y., S. MUNDLOS, J. UPTON and B. R. OLSEN, 1996 Altered growth and branching patterns in synpolydactyly caused by mutations in HOXD13. Science 272: 548-551.

PEARSON, C. E., K. N. EDAMURA and J. D. CLEARY, 2005 Repeat instability: Mechanisms of dynamic mutations. Nature Reviews Genetics 6: 729-742.

ROBINSON, D. O., S. R. HAMMANS, S. P. READ and J. SILLIBOURNE, 2005 Oculopharyngeal muscular dystrophy (OPMD): analysis of the PABPN1 gene expansion sequence in 86 patients reveals 13 different expansion types and further evidence for unequal recombination as the mutational mechanism. Human Genetics 116: 267-271.

SHAW, C. J., and J. R. LUPSKI, 2004 Implications of human genome architecture for rearrangement-based disorders: the genomic basis of disease. Human Molecular Genetics 13: R57-R64.

SINDEN, R. R., V. N. POTAMAN, E. A. OUSSATCHEVA, C. E. PEARSON, Y. L. LYUBCHENKO et al., 2002 Triplet repeat DNA structures and human genetic disease: dynamic mutations from dynamic DNA. Journal of Biosciences 27: 53-65.

STROMME, P., M. E. MANGELSDORF, M. A. SHAW, K. M. LOWER, S. M. E. LEWIS et al., 2002 Mutations in the human ortholog of Aristaless cause X-linked mental retardation and epilepsy. Nature Genetics 30: 441-445.

TROCHET, D., S. J. HONG, J. K. LIM, J. F. BRUNET, A. MUNNICH et al., 2005 Molecular consequences of PHOX2B missense, frameshift and alanine expansion mutations leading to autonomic dysfunction. Human Molecular Genetics 14: 3697- 3708.

TROCHET, D., L. DE PONTUAL, B. KEREN, A. MUNNICH, M. VEKEMANS et al., 2007 Polyalanline expansions might not result from unequal crossing-over. Human Mutation 28: 1043-1044.

UTSCH, B., K. BECKER, D. BROCK, M. J. LENTZE, F. BIDLINGMAIER et al., 2002 A novel stable polyalanine [poly(A)] expansion in the HOXA13 gene associated with hand-foot-genital syndrome: proper function of poly(A)-harbouring transcription factors depends on a critical repeat length? Human Genetics 110: 488-494.

WAROT, X., C. FROMENTAL-RAMAIN, V. FRAULOB, P. CHAMBON and P. DOLLE, 1997 Gene dosage-dependent effects of the Hoxa-13 and Hoxd-13 mutations on morphogenesis of the terminal parts of the digestive and urogenital tracts. Development 124: 4781-4791.

WARREN, S. T., 1997 Polyalanine expansion in synpolydactyly might result from unequal crossing-over of HOXD13. Science 275: 408-409.

YUCEL, A., I. KURU, M. E. BOZAN, M. ACAR and M. SOLAK, 2005 Radiographic evaluation and unusual bone formations in different genetic patterns in synpolydactyly. Skeletal Radiology 34: 468-476.

21 ZAKANY, J., and D. DUBOULE, 1996 Synpolydactyly in mice with a targeted deficiency in the HoxD complex. Nature 384: 69-71.

ZAKANY, J., C. FROMENTALRAMAIN, X. WAROT and D. DUBOULE, 1997 Regulation of number and size of digits by posterior Hox genes: A dose-dependent mechanism with potential evolutionary implications. Proceedings of the National Academy of Sciences of the United States of America 94: 13695-13700.

Tables

Gene Repeat Name Reference Sequence HOXD13 15 ala normal allele GCGGCGGCGGCGGCAGCGGCGGCTGCGGCGGCGGCGGCGGCAGCC GCGGCGGCGGCGGCAGCGGCGGCTGCGGCGGCGGCGGCGGCAGCGGCGGC HOXD13 .+ 9 ala pedigree P Goodman et al. 1997 GGCGGCGGCGGCGGCGGCAGCC GCAGCAGCAGCGGCGGCGGCCGCGGCAGCGGCGGCGGCGGCAGCGGCAGC PHOX2B 20 ala normal allele GGCGGCAGCT GCAGCAGCAGCGGCGGCGGCCGCGGCAGCGGCGGCGGCGGCAGCGGCAGC PHOX2B .+ 5 ala Trochet et al. 2005/2007 GGCGGCGGCGGCAGCGGCAGCGGCT GCAGCAGCAGCGGCGGCGGCCGCGGCAGCGGCGGCGGCGGCAGCGGCGGC PHOX2B .+ 6 ala Trochet et al. 2005/2007 GGCCGCGGCAGCGGCGGCGGCGGCAGCT GCAGCAGCAGCGGCGGCGGGCGCAGCAGCAGCGGCGGCGGCCGCGGCAGC PHOX2B .+ 7 ala Trochet et al. 2005/2007 GGCGGCGGCGGCAGCGGCAGCGGCGGCAGCT

FOXL2 14 ala normal allele GCGGCAGCCGCAGCGGCTGCAGCAGCTGCGGCTGCAGCCGCG GCGGCAGCCGCAGCGGCTGCAGCAGCTGCGGCTGCAGCAGCTGCGGC FOXL2 .+ 15 ala 921–935trip15 De Baere et al. 2003 TGCAGCAGCTGCGGCTGCAGCAGCTGCGGCTGCAGCCGCG PABPN1 10 ala normal allele GCGGCGGCGGCGGCGGCGGCAGCAGCAGCG GCGGCGGCGGCGGCGGCGGCGGCGGCGGCGGCGGCGGCGGC PABPN1 .+ 7 ala (GCG)13 Robinson et al. 2005 AGCAGCAGCG

Table 1: List of the alanine track expansions that cannot be explained by the unequal recombination mechanism. These expansions were all found in human pedigrees.

Figure legends

Figure 1: Alizarin red colored skeletal preparations of adult distal forelimbs and hindlimbs of (a) wild-type (+/+), heterozygous (Dyc/+) and homozygous

(Dyc/Dyc) mutants and (b) heterozygous Del3, homozygous Del3 and transheterozygous Dyc/Del3 mice. Fingers are numbered from 1 to 5, from the to the little finger and shortening of fingers is shown by arrows. Dyc/+ and

22 Del3/+ autopods have shortening of fingers 2 and 5. Dyc/Dyc have shortening and malformation of all fingers and appearance of a supernumerary finger, noted 6, in post-axial position. Syndactily between finger 3 and 4 with fusion at the metacarpus or metatarsus level is also present. In hindlimbs, the hallux is shortened with the metatarsus having a characteristic form. Del3/Del3 have shortening of all fingers and high deformation of toe 1, syndactily of toe 4 and 5 (arrow), and supernumerary toe 6 in the hindlimb. Transheterozygote Dyc/Del3 animals show shortening of all fingers with a higher deformation of toe 1 (arrow).

Figure 2: Analysis of the Dyc mutation. (A) Alignment between BALB/c (sequence of reference) and Dyc/Dyc showing the amplification in 5’ of Hoxd13 Exon1. (B)

Alignment between the BALB/c, C57Bl/6J, Spdh and Dyc sequences. The duplicated nucleotides in the two mutants are shown in red and purple and the corresponding fragments in the normal sequence are underlined in the same colour. The arrows show the straightforward tandem duplication in the Spdh mutant and the imperfect amplification in Dyc where the unequal allelic homologous recombination model is at fault. (C) Multiple alignments of the BALB/c sequence to the Dyc (yellow) at microhomologies (red letters) that are boxed. The extended sequence in the Dyc mutant is underlined. Aligned sequences are in blue, green and purple, non-aligned ones in grey.

Figure 3: Application of FoSTeS to atypical polyalanine extensions. Multiple alignments of the normal human sequences of HOXD13, PHOX2B, FOXL2 and

PABPN1 to the sequences with atypical polyalanine extensions (yellow) at microhomologies (red letters) that are boxed. The extended sequence in the mutants

23 is underlined. Aligned sequences are in blue, green and purple, non-aligned ones in grey.

Suppl. figure 1: Autopod development in wild-type and Dyc/Dyc mutant mice. In each panel, a wild-type paw (left) is compared to Dyc/Dyc one (right). (a ) At 14.5 dpc, wild-type autopods have well individualized fingers with visible cartilaginous elements at the level of the carpus metacarpus, tarsi and metatarsi, whereas

Dyc/Dyc autopod development is delayed: cartilaginous blastema have only formed in the first phalanx and are not visible in the more distal parts of the future digits that are shorter. (b)At 15.5 dpc, of the cartilaginous elements are visible in both wild-type and mutant autopods, but the digits in Dyc/Dyc embryos reveals shortened and supernumerary digits (6+7) with deformation at the level of the metatarsus and metacarpus. (c) Alizarin red staining of the skeleton at 0.5 dpp shows ossification centres in all wild-type digits but not in Dyc/Dyc mutant forelimbs. In the anterior limbs polydactily is accompanied by syndactily of digits 2, 3, 4 and 5, 6. In posterior limbs, there is also polydactily with syndactily of digits 2, 3 and 4, 5. At 4.5 dpp, ossification of the phalanx is completed in wild-type autopods. In Dyc/Dyc autopods, some ossification has occurred, but localisation of the ossifications centres is not clearly defined and there are defects in segmentation of the finger elements, making it hard to distinguish metacarpus and phalanx..

Suppl. figure 2: Expression of Hox genes in wild-type and homozygous (Dyc/Dyc) forelimbs at 12.5 dpc. Whole-mount in situ hybridization of Hoxd10, Hoxd11, Hoxd12,

Hoxd13 and Hoxa13 probes indicates that there is no difference in expression profiles of theses genes between wild type and homozygote.

24