A Newly Described Bovine Type 2 Scurs Syndrome Segregates with a Frame-Shift Mutation in TWIST1

Abstract

The developmental pathways involved in horn development are complex and still poorly understood. Here we report the description of a new dominant inherited syndrome in the bovine Charolais breed that we have named type 2 scurs. Clinical examination revealed that, despite a strong phenotypic variability, all affected individuals show both horn abnormalities similar to classical scurs phenotype and skull interfrontal suture synostosis. Based on a genome-wide linkage analysis using Illumina BovineSNP50 BeadChip genotyping data from 57 half-sib and full-sib progeny, this locus was mapped to a 1.7 Mb interval on bovine chromosome 4. Within this region, the TWIST1 gene encoding a transcription factor was considered as a strong candidate gene since its haploinsufficiency is responsible for the human Saethre-Chotzen syndrome, characterized by skull coronal suture synostosis. Sequencing of the TWIST1 gene identified a c.148_157dup (p.A56RfsX87) frame-shift mutation predicted to completely inactivate this gene. Genotyping 17 scurred and 20 horned founders of our pedigree as well as 48 unrelated horned controls revealed a perfect association between this mutation and the type 2 scurs phenotype. Subsequent genotyping of 32 individuals born from heterozygous parents showed that homozygous mutated progeny are completely absent, which is consistent with the embryonic lethality reported in Drosophila and mouse suffering from TWIST1 complete insufficiency. Finally, data from previous studies on model species and a fine description of type 2 scurs symptoms allowed us to propose different mechanisms to explain the features of this syndrome. In conclusion, this first report on the identification of a potential causal mutation affecting horn development in cattle offers a unique opportunity to better understand horn ontogenesis.

About the Authors:

Aurélien Capitan

* E-mail: [email protected]

Affiliations Unité Mixte de Recherche 1313 Génétique Animale et Biologie Intégrative, Institut National de la Recherche Agronomique, Jouy-en-Josas, France, Service Génétique, Département Fédéral, Union Nationale des Coopératives d'Elevage et d'Insémination Animale, Paris, France

Cécile Grohs

Affiliation: Unité Mixte de Recherche 1313 Génétique Animale et Biologie Intégrative, Institut National de la Recherche Agronomique, Jouy-en-Josas, France

Bernard Weiss

Affiliation: Unité Mixte de Recherche 1313 Génétique Animale et Biologie Intégrative, Institut National de la Recherche Agronomique, Jouy-en-Josas, France

Marie-Noëlle Rossignol

Affiliation: Laboratoire d'Analyses Génétiques pour les Espèces Animales, Jouy-en-Josas, France

Patrick Reversé

Affiliation: Section Charolaise, Gènes Diffusion, Douai, France

André Eggen

Affiliation: Unité Mixte de Recherche 1313 Génétique Animale et Biologie Intégrative, Institut National de la Recherche Agronomique, Jouy-en-Josas, France

Introduction

Horns in bovine as in all members of the Cavicorn superfamily, are permanent and not ramified. They consist of a bony core covered by a corium producing the keratin sheath. Contrary to antlers in deer, the developmental pathways involved in horn formation have not been extensively studied and are still poorly understood.

Studies by Dove [1] contributed greatly to the comprehension of this complex process. Using tissue transplantation, Dove showed that: (i) the bony core is not an outgrowth of the skull but originates from a separated center of ossification located in the dermis and hypodermis of the calves' horn bud; (ii) the keratinization of the horn bud epidermis does not induce ossification of the underlying dermis and hypodermis and conversely, thus both phenomena are probably programmed during embryogenesis; (iii) the ossifying hypodermal tissue induces the frontal bone to grow upward and to form the base of the horn spike, then it fuses with the skull by dissolving it locally. (Figure S1). Thus, horn development is the result of the differentiation and remodeling of various tissues originating from two distinct germ layers: ectoderm and mesoderm.

Genetic abnormalities affecting horn development represent unique models to identify genes and pathways involved in this process. Two main approaches are generally used to achieve this goal: comparison between wild-type and affected horn buds gene expression (as recently used by Mariasegaram et al. [2]) or genetic mapping followed by candidate gene sequencing to identify the causal mutation. In this study, the latter approach was used to determine the genetic basis of the polled and scurs phenotypes in the French Charolais breed.

The polled phenotype is characterized by the complete absence of horns as well as of any type of corneous growth. On the contrary, scurs share similar shapes and locations with horns but they are generally smaller and characterized by an absence of fusion between the bony core and the skull [1], [3], [4]. Even if several exceptions have been reported (for a review see [5]), it is generally believed that the genetic determinism of these horn abnormalities involves the interaction of two autosomal biallelic loci: the polled and scurs loci. Indeed, the P allele of the polled locus is dominant and specifies the absence of wild type horns whereas the presence of scurs or the complete absence of appendage is determined by the Sc and sc alleles of the scurs locus, respectively [6]–[8]. Numerous studies have mapped the polled locus to the centromeric region of BTA01 in various breeds, but to date the causal mutation has not been identified and/or published [9]–[14]. However, only one study mapped the scurs locus on BTA19 in a crossbred pedigree [15] and we were not able to confirm this result in the French Charolais breed as reported in a previous study based on BTA19 microsatellites genotyping data [5].

In order to fine-map both loci, we performed Illumina BovineSNP50 genotyping on a French Charolais pedigree consisting of 323 individuals (73 horned, 153 scurred and 97 polled) representing 40 paternal and 35 maternal half-sib families (unpublished data). After haplotype reconstruction for the BTA01 centromeric region, two different haplotypes were identified among the polled individuals but absent among the horned individuals. To avoid potential bias due to different interactions between the scurs locus and two different polled mutations, we classified the polled and scurred individuals into two groups, according to their polled haplotype at BTA01, before performing the mapping of the scurs locus within each group. Interestingly, several scurred individuals could not be classified into these two groups. In other words, those animals were scurred without exhibiting one of the two identified polled haplotypes on BTA01. A pedigree analysis revealed that these animals are related to the same sire over a maximum of six generations and that the scurs phenotype is transmitted in a pattern consistent with autosomal dominant inheritance. However this transmission occured independently from the BTA01 haplotype pointing to a different etiology than the common scurred phenotype, the expression of which is fully dependent on the presence of the P allele from the polled locus [7], [8]. Based on these evidences, this new genetic disorder affecting horn development was called type 2 scurs.

In the study reported here, our objective was (i) to describe more precisely the type 2 scurs phenotype and (ii) to fine-map this locus and identify the causal mutation, in order to better understand the developmental pathways involved in bovine horn formation.

Results

Clinical findings

Visual examinations showed a strong phenotypic variability between individuals and genders among affected animals:

1. The size of scurs varies from small scabs to 15 cm-long appendages in adult females (Figure 1) whereas in adult males scurs are systematically massive (more than 10 cm-long) and often less mobile, but still not completely fused to the skull. Furthermore, scurs grow earlier in life in males: they are usually detected at first examination (between 4 and 6 months) unlike in females where they become visible at second examination (between 9 and 18 months) or later.

2. The structure of the keratin sheath also presents different levels of alteration depending on the scurs size: tiny scurs develop as scabs made of scaly patches, whereas small scurs and the terminal portion of long scurs are covered by irregular keratin sheets (Figure 2).

3. In addition, affected individuals show mild to pronounced acrocephaly and a ridge-shaped extra bone deposition along the interfrontal suture, which both appear to be negatively correlated with the size of scurs (Figure 1). These pathologies are also attested by the particular shape of the poll and the denser calcification of the interfrontal suture of the type 2 scurred vs horned skull radiographs (Figure 3).

[Figure omitted. See PDF.]

Figure 1. Phenotypic description of type 2 scurs syndrome.

Ten-(A), two-(B) and four-(C) year old females affected by type 2 scurs syndrome. Note the marked phenotypic variability and the negative correlation between the size of scurs and the importance of acrocephaly and ridge-shaped extra bone deposition along the interfrontal suture. (D). Four-year old horned female.

https://doi.org/10.1371/journal.pone.0022242.g001

[Figure omitted. See PDF.]

Figure 2. Details of the type 2 scurs keratin sheath.

(A). Scaly patches. (B). Small scurs (∼2 cm) with irregular keratin sheath. (C and D). Long scurs (∼15 cm) with an irregular keratin sheath at their end (see arrows). (E). End of a normal horn (∼25 cm) with a regular keratin sheath.

https://doi.org/10.1371/journal.pone.0022242.g002

[Figure omitted. See PDF.]

Figure 3. Frontal radiographs of type 2 scurred and horned skulls.

(A). Control skull of a five-year old horned female. (B). Skull of a four-year old female affected by the type 2 scurs and carrying long scurs. Note: (1) the particular shape of the poll, (2) the denser calcification of the interfrontal suture, (3) the absence of fusion between the frontal bone and the horn bony core, (4) the absence of frontal bone drawing up and (5) the non-pneumatization of the bony core in affected vs. horned skulls.

https://doi.org/10.1371/journal.pone.0022242.g003

Radiographs also revealed that (i) the frontal bone of affected individuals is not drawn up to form the basilar portion of the horn spike (ii) the scurs bony core is not pneumatized and (iii) the space between the skull and the bony core is filled with soft tissues.

In summary, despite a strong phenotypic variability, all affected individuals show both horn abnormalities similar to classical scurs phenotype and skull interfrontal suture synostosis, constituting a single pathological entity (i.e. a syndrome). These animals are readily distinguishable from their wild-type relatives and from the wild-type controls: (i) horns of wild-type male and female born from affected dams are always visible at first examination and often already fused to the skull at this time; (ii) in adults they are firmly attached to the skull, their keratin sheath is regular and their size is longer than 15 cm; (iii) finally none of the wild-type animals show evidence of acrocephaly and extra bone deposition along the skull interfrontal suture.

Mapping of the type 2 scurs gene

As shown in Figure 4, the genome-wide scan revealed a significant linkage (maximum LOD score of 7.2) between the type 2 scurs phenotype and several clusters of markers located on chromosome 4. The 95% confidence interval spanned 1.7 Mb (from marker ARS-BFGL-NGS-57582 to BTB-01114634) encompassing six different genes: SNX13, PRPS1L1, HDAC9, UBE2D4, TWIST1 and FERD3L. Among them TWIST1 was the most compelling candidate gene since it encodes a basic helix-loop-helix (bHLH) transcription factor regulating many processes including cranial suture patterning and fusion [16], [17]. Numerous mutations in this gene have been reported to be responsible for the human Saethre-Chotzen syndrome (SCS; OMIM#101400) also known as Acrocephalosyndactyly type III. The hallmark of this autosomal-dominant syndrome is acrocephaly (due to premature fusion –synostosis– of the skull coronal suture) associated with variable additional features such as mild craniofacial and limb deformities (for a review see [18], [19]).

[Figure omitted. See PDF.]

Figure 4. Genome-wide multipoint linkage analysis of type 2 scurs syndrome.

The genome-wide scan reveals a significant linkage (maximum LOD score = 7.2) with clusters of markers located on chromosome 4. The 95% confidence interval spans 1.7 Mb encompassing six different genes.

https://doi.org/10.1371/journal.pone.0022242.g004

Mutation analysis

Sequencing the entire TWIST1 gene in a trio consisting of two affected females and one unaffected male allowed us to identify a 10-bp duplication (c.148_157dup) in a GC-rich fragment of exon 1. Subsequent genotyping of this mutation on a broader panel of animals revealed a perfect association with the type 2 scurs phenotype: all 17 affected founders of our pedigree were heterozygous whereas all 20 non-affected founders and 48 unrelated controls were homozygous for the wild-type allele.

Assumed consequence of the TWIST1 c.148_157dup mutation

TWIST1 has two highly conserved domains: the basic helix–loop–helix (bHLH) domain and the tryptophan–arginine (WR) domain [20] (Figure 5). The first domain is a bipartite domain for DNA binding (basic motif) and protein-protein interactions (helix-loop-helix motif) shared by numerous transcription factors [23]–[26]. The second one, also known as the TWIST box since it is specific to this subfamily, has been shown to inhibit the function of the Runx2 DNA binding domain [27] and to be essential for the transactivating function of TWIST1 in mice [28].

[Figure omitted. See PDF.]

Figure 5. Characterization of the TWIST1 c.148_157dup mutation.

(A) TWIST1 gene's organization scheme; (B) DNA sequencing chromatograms showing the c.148_157dup mutation of type 2 scurs affected animals and the wild type allele; (C) putative p.A56RfsX87 mutated protein and multispecies alignment of the TWIST1 protein sequence using CLUSTALW [21]. Cattle (Bta), human (Hsa), mouse (Mmu), chicken (Gga), anolis (Aca), xenopus (Xtr) and zebrafish (Dre) sequences accession numbers are respectively DAA30767, NP_000465, AAH33434, NP_990070, DAA06059, AAH74558, ABC73066 in Genbank. The frame-shift change in the bovine p.A56RfsX87 mutated protein (Mut.) is underlined. The basic, helix-loop-helix and tryptophan–arginine domains are respectively highlighted with dashed, solid and mixed dashed boxes [22].

https://doi.org/10.1371/journal.pone.0022242.g005

As shown in Figure 5, the c.148_157dup mutation is predicted to cause a frame-shift change of TWIST1 with Alanine-56 as the first affected amino acid (p.A56RfsX87). This frame-shift is assumed to completely inactivate this gene because the 142-amino acid mutated protein lacks both above mentioned functional domains.

Since targeted twist-null mutations are embryonic lethal in Drosophila [29], [30] and mouse [31], we genotyped 32 individuals born from heterozygous parents to look for distorted Mendelian transmission. As presented in Table 1, no homozygous individual for the c.148_157dup was observed. Moreover the observed genotype distribution fitted its expectation under the hypothesis of lethality of this mutation in the homozygous state. Thus, we conclude that the TWIST1 c.148_157dup mutation associated with type 2 scurs may be homozygous lethal.

[Figure omitted. See PDF.]

Table 1. Results of TWIST1 c.148_157dup mutation genotyping in 32 progeny born from heterozygous parents.

https://doi.org/10.1371/journal.pone.0022242.t001

Discussion

In this study, we describe the identification, the characterization and the fine-mapping of a new genetic defect affecting both horn development and skull interfrontal suture fusion in cattle.

This new syndrome was named type 2 scurs for its striking similarity with the already known scurs phenotype. In both syndromes, horn appendages are indistinguishable at the macroscopic scale: (i) the frontal bone is not drawn up to form the basilar portion of the horn spike; (ii) the space between the skull and the bony core is filled with soft tissues; (iii) the bony core is densely ossified and covered by an irregular keratin sheath, and finally (iv) there is a marked phenotypic variability between individual and gender [1], [3]. Such similarity suggests a close etiology between these syndromes: genes involved in the same metabolic pathway might be responsible for these abnormalities.

Our study shows that the similarity between these two disorders could have interfered with the identification of the genetic determinism of the “classical” scurs and the type 2 scurs, explaining the numerous exceptions [5] reported to the scurs genetic determinism proposed by Long and Gregory [7] and Brem et al. [8].

In addition, we identified a frameshift mutation (p.A56RfsX87) predicted to inactivate TWIST1 and demonstrated that a perfect association exists between this mutation and the type 2 scurs syndrome.

Although we cannot provide functional proof of the causality of the p.A56RfsX87 mutation at this time, the large amount of functional data available for the TWIST1 gene strongly supports this hypothesis.

Indeed, more than 80 mutations in the TWIST1 gene have been reported in humans as the cause of skull coronal suture synostosis, a symptom of the Saethre-Chotzen syndrome [18], [19], [32]–[36]. Functional studies have identified TWIST1 haploinsufficiency as the disease-causing mechanism of this syndrome and unravelled the major role played by this gene in the regulation of cranial suture patterning and fusion [16], [17], [27], [37]–[43]. Namely TWIST1 deficiency inhibits osteogenic stem cells proliferation and leads to premature osteoblast differentiation altering the balance between these two phenomena which is essential for normal sutural growth [44]. Therefore, we assume that the p.A56RfsX87 mutation causes TWIST1 haploinsufficiency which in turn is responsible for the skull interfrontal suture synostosis observed in cattle affected by the type 2 scurs syndrome.

Moreover the absence of homozygous individuals for the mutation in our second pedigree is consistent with the embryonic lethality reported in Drosophila and mouse suffering from TWIST1 complete insufficiency [29]–[31]. Contrary to the Drosophila mutant, twist1 -/- mouse embryos undergo normal gastrulation. However, later they display severe defects in cranial neural tube closure, head mesenchyme, somites and limb buds and finally die at E10.5–11 [31]. These experiments have underlined the critical role played by TWIST1 in diverse developmental pathways during embryogenesis such as specification of the mesodermal somites' derivatives and neural crest cell migration and differentiation (for a review see [45]–[48]). Thus, we believe that homozygosity for the bovine p.A56RfsX87 mutation causes TWIST1 complete insufficiency and is embryonic lethal.

Since all available model animals are not horned, it is not possible to infer from previous studies the possible mode of action of the TWIST1 p.A56RfsX87 mutation on the development of scurs. However, fine examination of the affected animals reveals interesting clues on the underlying mechanisms.

As reported in the human Saethre-Chotzen syndrome, there is a marked phenotypic variability (including horn development) among individuals carrying the same TWIST1 mutation. Moreover, there is a positive correlation between the level of horn abnormality and the importance of craniosynostosis symptoms. These facts suggest that the mechanism causing the scurs phenotype is the same than for craniosynostosis, i.e. TWIST1 haploinsufficiency. Nevertheless contrary to craniosynostosis, the antiosteogenic function of TWIST1 cannot be the only cause of horn abnormalities as attested hereafter.

Despite the marked phenotypic variability, all p.A56RfsX87/+ animals present both bony core and keratin sheath abnormalities. Since the keratinization of the normal horn bud epidermis and the ossification of the underlying dermis and hypodermis are both programmed during embryogenesis and not induced by each other after birth [1], this observation rules out a possible action of an abnormal bony core on the organization of the overlying corium. Rather, it advocates for three possible etiologies: (i) an early role of TWIST1 in horn bud cells programming during embryogenesis or fetal life; (ii) postnatal requirement of TWIST1 for the modification of normally programmed horn bud epidermis and underlying tissues in the corium and bony core respectively; or (iii) both.

The first hypothesis is consistent with the above-mentioned critical role played by TWIST1 in cell lineage specification and differentiation during embryogenesis. The second is supported by the unique role played by TWIST1 in promoting cell dedifferentiation, migration and proliferation in processes like cancer or Epithelial-Mesenchymal Transition (EMT) [49], [50]. Interestingly, Mariasegaram et al. [2] have observed a marked enrichment of gene networks relating to EMT by studying the differentially expressed genes between horn buds from 1 to 2 week-old polled, scurred and horned calves using a bovine gene expression microarray. However, they did not report a significant difference in TWIST1 expression between the three categories. This result suggests an earlier involvement of TWIST1 in processes leading to epidermis keratinisation and bony core ossification (hypothesis 1).

In conclusion, we describe a new autosomal dominant inherited syndrome characterised by horn development anomalies, craniosynostosis and an absence of homozygous affected calves. In addition, we have identified the p.A56RfsX87 mutation in the TWIST1 gene as the candidate causative mutation and propose different mechanisms involving TWIST1 haploinsufficiency in diverse developmental pathways to explain the three main features of this syndrome. To our knowledge, this is the first report on the identification of a potential causal mutation affecting horn development in cattle. Better, among Bovinae, type 2 scurs would be the only genetic disorder affecting horn development explained by a simple mechanism possibly involving loss of gene function. This makes type 2 scurs an ideal model to study horn ontogenesis. The detailed involvement of TWIST1 in horn development remains to be investigated by functional studies.

Materials and Methods

Ethics statement

Experiments reported in this work comply with the French National Institute for Agricultural Research (INRA) ethical guidelines. Animals were extremely well cared. Approval by the INRA Ethical Committee was not necessary for blood sampling, sperm sampling, radiographs and routine husbandry procedures. Blood was collected by the following agricultural technicians licensed by the French Etablissements Départementaux de l'Elevage (EDE): Rémi Bierbaum, Christophe Caron, Vincent Colas, Michel Dewaele, Bruno Elmanowsky, Denis Faradji, Sébastien Landemaine, Jean-Marie Moinel, Arnaud Poilvert, Bernard Raimbault, Stéphane Thibaux and Arnaud Tranier. Sperm was obtained from semen straws generously provided by Gènes Diffusion and UCATRC breeding companies. Radiographs were performed by Dominique Rémy and Guillaume Belbis (licensed veterinarians). All the samples were obtained with the permission of the French Polled Charolais Program.

Animals

Seventeen scurred dams, 20 horned Artificial Insemination (AI) sires and their 57 progeny (40 scurred and 17 horned) were genotyped for linkage analysis. These individuals belong to the French Polled Charolais Program which aims at producing high genetic value polled sires by mating the best horned AI sires to polled and scurred cows for several generations. All the scurred dams are related to the same sire over a maximum of five generations.

Sequencing of the candidate gene was performed on a trio consisting of a scurred dam, a horned bull and the scurred heifer. In addition, 48 unrelated horned Charolais were recruited as controls. Finally, 32 additional individuals born from heterozygous parents (five sires and 15 dams) were used to study the Mendelian transmission of the candidate mutation.

Phenotypes

The progeny were phenotyped twice as described in Capitan et al. [5] whereas founders and control individuals were phenotyped once at adulthood. All types of corneous growths that were loosely attached to the skull were considered as scurs [1], [3]–[5]. Following the identification of type 2 scurs, most of the affected individuals were re-examined by visual inspection to refine the phenotype and to search for other associated abnormalities. To complete this study, the skulls of a four-year old female carrying long scurs and a five-year old horned control were radiographed using a Gierth HF 80 Plus (Vtrade international, Fernelmont, Belgium) with the following parameters: 50 volts, 10 mA and 80 cm. Finally, a survey was carried out among breeders to collect all past observations noted among such animals and their ancestors.

Samples

DNA was extracted from blood using the Wizard® Genomic DNA purification Kit (Promega, Charbonnières-les-bains, France) or from sperm using a standard phenol-chloroform method.

Linkage analysis

DNA samples were genotyped with the Illumina BovineSNP50 BeadChip [51]. Marker order and map distances were based on the bovine sequence assembly Btau_4.0, assuming 1 Mb for 1 cM. A genome-wide multipoint linkage analysis was performed using MERLIN software (version 1.1.2) [52] and assuming a dominant model of disease inheritance with an allele frequency of 0.0001. Penetrance values were set at 0.01, 0.99 and 0.99 for the homozygous wild-type, heterozygous and homozygous affected individuals, respectively. The marker-marker linkage disequilibrium was modeled using the – rsq 0.1 option. Finally, genotypes for BTAX markers of the bull-calves sires were set as missing to avoid Mendelian errors.

Mutation analysis

PCR primers covering the whole TWIST1 gene were designed from the bovine genome sequence assembly Btau_4.0 with Primer3 sofware [53] (Table S1). PCR reactions were performed using the Go-Taq Flexi (Promega, Charbonnières-les-bains, France) or the GC-RICH PCR System (Roche, Meylan, France) according to the manufacturer's instructions on a PTC-100 thermocycler (BioRad, Marnes-la-Coquette, France). The resulting amplicons were purified on MultiScreen PCR96 Filter Plates (Millipore, Molsheim, France) and bidirectionnally sequenced by Qiagen (Hilden, Germany) using conventional Sanger sequencing. Polymorphism was detected with the NovoSNP software [54]. The candidate mutation was subsequently genotyped by PCR-sequencing of exon 1 under the same conditions.

Supporting Information

[Figure omitted. See PDF.]

Figure S1.

Horn development stages adapted from Dove's (1935) experimental report.

https://doi.org/10.1371/journal.pone.0022242.s001

(PDF)

Table S1.

TWIST1 primers sequences.

https://doi.org/10.1371/journal.pone.0022242.s002

(PDF)

Acknowledgments

The authors would like to thank the Gènes Diffusion and UCATRC breeding companies as well as the Charolais breeders for generously providing samples and phenotypic information. We are particularly grateful to Marie-Thérèse Blanchon (Gènes Diffusion), Bernard Boiret and GAEC Humeau for the excellent collaboration and the warm welcome. We acknowledge François Ménissier (INRA, retired) for exciting discussions and for sharing his valuable archives on the founder stocks of the present French polled Charolais population. The contributions of Carine Bernard-Capel (Institut de l'Elevage) and Didier Boichard (INRA) in writing the Apis Gene-Hornout funding convention and of Arnaud Delpeuch (Institut de l'Elevage) in building the Prosanco phenotype database are highly appreciated. We also wish to thank Christophe Richard and the staff of the UCEA INRA experimental farm in Bressonvilliers (France) for maintaining three type 2 scurs affected cows, as well as Dominique Rémy and Guillaume « Gaston » Belbis (Ecole Nationale Vétérinaire d'Alfort) for having generously radiographed these animals. We finally acknowledge Dominique Rocha and Didier Boichard (INRA) for revising the paper and Hélène Hayes (INRA) for English correction of the manuscript.

Author Contributions

Conceived and designed the experiments: AC AE. Performed the experiments: AC CG M-NR. Analyzed the data: AC. Contributed reagents/materials/analysis tools: BW PR. Wrote the paper: AC.

Citation: Capitan A, Grohs C, Weiss B, Rossignol M-N, Reversé P, Eggen A (2011) A Newly Described Bovine Type 2 Scurs Syndrome Segregates with a Frame-Shift Mutation in TWIST1. PLoS ONE6(7): e22242. https://doi.org/10.1371/journal.pone.0022242

References

1. Dove WF (1935) The physiology of horn growth: a study of the morphogenesis, the interaction of tissues and the evolutionary processes of a Mendelian recessive character by means of transplantation of tissues. J exp Zool 69: 347–405.WF Dove1935The physiology of horn growth: a study of the morphogenesis, the interaction of tissues and the evolutionary processes of a Mendelian recessive character by means of transplantation of tissues.J exp Zool69347405

2. Mariasegaram M, Reverter A, Barris W, Lehnert SA, Dalrymple B, et al. (2010) Transcription profiling provides insights into gene pathways involved in horn and scurs development in cattle. BMC Genomics 11: 370.M. MariasegaramA. ReverterW. BarrisSA LehnertB. Dalrymple2010Transcription profiling provides insights into gene pathways involved in horn and scurs development in cattle.BMC Genomics11370

3. Cole LJ (1924) The Wisconsin experiment in cross-breeding cattle. Proc World's Dairy Congress 2: 1383–1388.LJ Cole1924The Wisconsin experiment in cross-breeding cattle.Proc World's Dairy Congress213831388

4. Williams HD, Williams T (1952) The inheritance of horns and their modifications in polled Hereford cattle. J Hered 43: 267–272.HD WilliamsT. Williams1952The inheritance of horns and their modifications in polled Hereford cattle.J Hered43267272

5. Capitan A, Grohs C, Gautier M, Eggen A (2009) The scurs inheritance: new insights from the French Charolais breed. BMC Genet 10: 33.A. CapitanC. GrohsM. GautierA. Eggen2009The scurs inheritance: new insights from the French Charolais breed.BMC Genet1033

6. White WT, Ibsen HL (1936) Horn inheritance in Galloway-Holstein Cattle crosses. J Genet 32: 33–49.WT WhiteHL Ibsen1936Horn inheritance in Galloway-Holstein Cattle crosses.J Genet323349

7. Long CR, Gregory KE (1978) Inheritance of the horned, scurred and polled condition in cattle. J Hered 69: 395–400.CR LongKE Gregory1978Inheritance of the horned, scurred and polled condition in cattle.J Hered69395400

8. Brem G, Karnbaum B, Rosenberg E (1982) Zur Vererbung der Hornlosigkeit beim Fleckvieh. Bayer. Landwirtsch. Jahrb 59: 688–695.G. BremB. KarnbaumE. Rosenberg1982Zur Vererbung der Hornlosigkeit beim Fleckvieh. Bayer. Landwirtsch.Jahrb59688695

9. Georges M, Drinkwater R, King T, Mishra A, Moore SS, et al. (1993) Microsatellite mapping of a gene affecting horn development in Bos Taurus. Nat Genet 4: 206–210.M. GeorgesR. DrinkwaterT. KingA. MishraSS Moore1993Microsatellite mapping of a gene affecting horn development in Bos Taurus.Nat Genet4206210

10. Schmutz SM, Marquess FL, Berryere TG, Moker JS (1995) DNA marker-assisted selection of the polled condition in Charolais cattle. Mamm Genome 6: 710–713.SM SchmutzFL MarquessTG BerryereJS Moker1995DNA marker-assisted selection of the polled condition in Charolais cattle.Mamm Genome6710713

11. Brenneman RA, Davis SK, Sanders JO, Burns BM, Wheeler TC, et al. (1996) The polled locus maps to BTA1 in a Bos indicus X Bos taurus cross. J Hered 87: 156–161.RA BrennemanSK DavisJO SandersBM BurnsTC Wheeler1996The polled locus maps to BTA1 in a Bos indicus X Bos taurus cross.J Hered87156161

12. Harlizius B, Tammen I, Eichler K, Eggen A, Hetzel DJ (1997) New markers on bovine chromosome 1 are closely linked to the polled gene in Simmental and Pinzgauer cattle. Mamm Genome 8: 255–257.B. HarliziusI. TammenK. EichlerA. EggenDJ Hetzel1997New markers on bovine chromosome 1 are closely linked to the polled gene in Simmental and Pinzgauer cattle.Mamm Genome8255257

13. Drögemüller C, Wöhlke A, Leeb T, Distl O (2005a) A 4 Mb high resolution BAC contig on bovine chromosome 1q12 and comparative analysis with human chromosome 21q22. Comp Funct Genomics 6: 194–203.C. DrögemüllerA. WöhlkeT. LeebO. Distl2005aA 4 Mb high resolution BAC contig on bovine chromosome 1q12 and comparative analysis with human chromosome 21q22.Comp Funct Genomics6194203

14. Drögemüller C, Wöhlke A, Mömke S, Distl O (2005b) Fine mapping of the polled locus to a 1-Mb region on bovine chromosome 1q12. Mamm Genome 16: 613–620.C. DrögemüllerA. WöhlkeS. MömkeO. Distl2005bFine mapping of the polled locus to a 1-Mb region on bovine chromosome 1q12.Mamm Genome16613620

15. Asai M, Berryere TG, Schmutz SM (2004) The scurs locus in cattle maps to bovine chromosome 19. Anim Genet 35: 34–39.M. AsaiTG BerryereSM Schmutz2004The scurs locus in cattle maps to bovine chromosome 19.Anim Genet353439

16. Connerney J, Andreeva V, Leshem Y, Muentener C, Mercado MA, et al. (2006) Twist1 dimer selection regulates cranial suture patterning and fusion. Dev Dyn 235: 1345–1357.J. ConnerneyV. AndreevaY. LeshemC. MuentenerMA Mercado2006Twist1 dimer selection regulates cranial suture patterning and fusion.Dev Dyn23513451357

17. Connerney J, Andreeva V, Leshem Y, Mercado MA, Dowell K, et al. (2008) Twist1 homodimers enhance FGF responsiveness of the cranial sutures and promote suture closure. Dev Biol 318: 323–334.J. ConnerneyV. AndreevaY. LeshemMA MercadoK. Dowell2008Twist1 homodimers enhance FGF responsiveness of the cranial sutures and promote suture closure.Dev Biol318323334

18. Rose CS, Malcolm S (1997) A TWIST in development. Trends Genet 13: 384–387.CS RoseS. Malcolm1997A TWIST in development.Trends Genet13384387

19. Gripp KW, Zackai EH, Stolle CA (2000) Mutations in the human TWIST gene. Hum Mutat 15: 150–155.KW GrippEH ZackaiCA Stolle2000Mutations in the human TWIST gene.Hum Mutat15150155

20. Spring J, Yanze N, Middel AM, Stierwald M, Groger H, et al. (2000) The mesoderm specification factor twist in the life cycle of jellyfish. Dev Biol 228: 363–375.J. SpringN. YanzeAM MiddelM. StierwaldH. Groger2000The mesoderm specification factor twist in the life cycle of jellyfish.Dev Biol228363375

21. Thompson JD, Higgins DG, Gibson TJ (1994) CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position specific gap penalties and weight matrix choice. Nucleic Acids Res 22: 4673–4680.JD ThompsonDG HigginsTJ Gibson1994CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position specific gap penalties and weight matrix choice.Nucleic Acids Res2246734680

22. UniProt Consortium (2010) The Universal Protein Resource (UniProt) in 2010. Nucleic Acids Res 38: D142–148.UniProt Consortium2010The Universal Protein Resource (UniProt) in 2010.Nucleic Acids Res38D142148Available: http://www.uniprot.org. Available: http://www.uniprot.org.

23. Murre C, McCaw PS, Baltimore D (1989) A new DNA binding and dimerization motif in immunoglobulin enhancer binding, daughterless, MyoD, and myc proteins. Cell 56: 777–783.C. MurrePS McCawD. Baltimore1989A new DNA binding and dimerization motif in immunoglobulin enhancer binding, daughterless, MyoD, and myc proteins.Cell56777783

24. Murre C, McCaw PS, Vaessin H, Caudy M, Jan LY, et al. (1989) Interactions between heterologous helix-loop-helix proteins generate complexes that bind specifically to a common DNA sequence. Cell 58: 537–544.C. MurrePS McCawH. VaessinM. CaudyLY Jan1989Interactions between heterologous helix-loop-helix proteins generate complexes that bind specifically to a common DNA sequence.Cell58537544

25. Murre C, Bain G, van Dijk MA, Engel I, Furnari BA, et al. (1994) Structure and function of helix-loop-helix proteins. Biochim Biophys Acta 1218: 129–135.C. MurreG. BainMA van DijkI. EngelBA Furnari1994Structure and function of helix-loop-helix proteins.Biochim Biophys Acta1218129135

26. Atchley WR, Fitch WM (1997) A natural classification of the basic helix-loop-helix class of transcription factors. Proc Natl Acad Sci U S A 94: 5172–5176.WR AtchleyWM Fitch1997A natural classification of the basic helix-loop-helix class of transcription factors.Proc Natl Acad Sci U S A9451725176

27. Bialek P, Kern B, Yang X, Schrock M, Sosic D, et al. (2004) A twist code determines the onset of osteoblast differentiation. Dev Cell 6: 423–435.P. BialekB. KernX. YangM. SchrockD. Sosic2004A twist code determines the onset of osteoblast differentiation.Dev Cell6423435

28. Laursen KB, Mielke E, Iannaccone P, Füchtbauer EM (2007) Mechanism of transcriptional activation by the proto-oncogene Twist1. J Biol Chem 282: 34623–34633.KB LaursenE. MielkeP. IannacconeEM Füchtbauer2007Mechanism of transcriptional activation by the proto-oncogene Twist1.J Biol Chem2823462334633

29. Simpson P (1983) Maternal-Zygotic Gene Interactions during Formation of the Dorsoventral Pattern in Drosophila Embryos. Genetics 105: 615–632.P. Simpson1983Maternal-Zygotic Gene Interactions during Formation of the Dorsoventral Pattern in Drosophila Embryos.Genetics105615632

30. Nüsslein-Volhard C, Wieschaus E, Kluding H (1984) Mutations affecting the pattern of the larval cuticle in Drosophila melanogaster. I. Zygotic loci on the second chromosome. Roux's Arch Dev Biol 193: 267–282.C. Nüsslein-VolhardE. WieschausH. Kluding1984Mutations affecting the pattern of the larval cuticle in Drosophila melanogaster. I. Zygotic loci on the second chromosome.Roux's Arch Dev Biol193267282

31. Chen ZF, Behringer RR (1995) Twist is required in head mesenchyme for cranial neural tube morphogenesis. Genes Dev 9: 686–699.ZF ChenRR Behringer1995Twist is required in head mesenchyme for cranial neural tube morphogenesis.Genes Dev9686699

32. Howard TD, Paznekas WA, Green ED, Chiang LC, Ma N, et al. (1997) Mutations in TWIST, a basic helix-loop-helix transcription factor, in Saethre-Chotzen syndrome. Nat Genet 15: 36–41.TD HowardWA PaznekasED GreenLC ChiangN. Ma1997Mutations in TWIST, a basic helix-loop-helix transcription factor, in Saethre-Chotzen syndrome.Nat Genet153641

33. El Ghouzzi V, Le Merrer M, Perrin-Schmitt F, Lajeunie E, Benit P, et al. (1997) Mutations of the TWIST gene in the Saethre-Chotzen syndrome. Nat Genet 15: 42–46.V. El GhouzziM. Le MerrerF. Perrin-SchmittE. LajeunieP. Benit1997Mutations of the TWIST gene in the Saethre-Chotzen syndrome.Nat Genet154246

34. Rose CS, Patel P, Reardon W, Malcolm S, Winter RM (1997) The TWIST gene, although not disrupted in Saethre-Chotzen patients with apparently balanced translocations of 7p21, is mutated in familial and sporadic cases. Hum Mol Genet 6: 1369–1373.CS RoseP. PatelW. ReardonS. MalcolmRM Winter1997The TWIST gene, although not disrupted in Saethre-Chotzen patients with apparently balanced translocations of 7p21, is mutated in familial and sporadic cases.Hum Mol Genet613691373

35. Johnson D, Horsley SW, Moloney DM, Oldridge M, Twigg SR, et al. (1998) A comprehensive screen for TWIST mutations in patients with craniosynostosis identifies a new microdeletion syndrome of chromosome band 7p21.1. Am J Hum Genet 63: 1282–1293.D. JohnsonSW HorsleyDM MoloneyM. OldridgeSR Twigg1998A comprehensive screen for TWIST mutations in patients with craniosynostosis identifies a new microdeletion syndrome of chromosome band 7p21.1.Am J Hum Genet6312821293

36. Paznekas WA, Cunningham ML, Howard TD, Korf BR, Lipson MH, et al. (1998) Genetic heterogeneity of Saethre-Chotzen syndrome, due to TWIST and FGFR mutations. Am J Hum Genet 62: 1370–1380.WA PaznekasML CunninghamTD HowardBR KorfMH Lipson1998Genetic heterogeneity of Saethre-Chotzen syndrome, due to TWIST and FGFR mutations.Am J Hum Genet6213701380

37. El Ghouzzi V, Legeai-Mallet L, Aresta S, Benoist C, Munnich A, et al. (2000) Saethre–Chotzen mutations cause TWIST protein degradation or impaired nuclear location. Hum Mol Genet 9: 813–819.V. El GhouzziL. Legeai-MalletS. ArestaC. BenoistA. Munnich2000Saethre–Chotzen mutations cause TWIST protein degradation or impaired nuclear location.Hum Mol Genet9813819

38. El Ghouzzi V, Legeai-Mallet L, Benoist-Lasselin C, Lajeunie E, Renier D, et al. (2001) Mutations in the basic domain and the loop–helix II junction of TWIST abolish DNA binding in Saethre–Chotzen syndrome. FEBS Lett 492: 112–118.V. El GhouzziL. Legeai-MalletC. Benoist-LasselinE. LajeunieD. Renier2001Mutations in the basic domain and the loop–helix II junction of TWIST abolish DNA binding in Saethre–Chotzen syndrome.FEBS Lett492112118

39. El Ghouzzi V, Bonaventure J, Munnich A (2001) TWIST: un nouvel acteur de l'ossification des os plats. Médecine/sciences 17: 1281–1288.V. El GhouzziJ. BonaventureA. Munnich2001TWIST: un nouvel acteur de l'ossification des os plats.Médecine/sciences1712811288

40. Yousfi M, Lasmoles F, El Ghouzzi V, Marie PJ (2002) Twist haploinsufficiency in Saethre-Chotzen syndrome induces calvarial osteoblast apoptosis due to increased TNFalpha expression and caspase-2 activation. Hum Mol Genet 11: 359–369.M. YousfiF. LasmolesV. El GhouzziPJ Marie2002Twist haploinsufficiency in Saethre-Chotzen syndrome induces calvarial osteoblast apoptosis due to increased TNFalpha expression and caspase-2 activation.Hum Mol Genet11359369

41. Firulli BA, Krawchuk D, Centonze VE, Vargesson N, Virshup DM, et al. (2005) Altered Twist1 and Hand2 dimerization is associated with Saethre-Chotzen syndrome and limb abnormalities. Nat Genet 37: 373–381.BA FirulliD. KrawchukVE CentonzeN. VargessonDM Virshup2005Altered Twist1 and Hand2 dimerization is associated with Saethre-Chotzen syndrome and limb abnormalities.Nat Genet37373381

42. Funato N, Twigg SR, Higashihori N, Ohyama K, Wall SA, et al. (2005) Functional analysis of natural mutations in two TWIST protein motifs. Hum Mutat 25: 550–556.N. FunatoSR TwiggN. HigashihoriK. OhyamaSA Wall2005Functional analysis of natural mutations in two TWIST protein motifs.Hum Mutat25550556

43. Guenou H, Kaabeche K, Mée SL, Marie PJ (2005) A role for fibroblast growth factor receptor-2 in the altered osteoblast phenotype induced by Twist haploinsufficiency in the Saethre-Chotzen syndrome. Hum Mol Genet 14: 1429–1439.H. GuenouK. KaabecheSL MéePJ Marie2005A role for fibroblast growth factor receptor-2 in the altered osteoblast phenotype induced by Twist haploinsufficiency in the Saethre-Chotzen syndrome.Hum Mol Genet1414291439

44. Johnson D, Iseki S, Wilkie AO, Morriss-Kay GM (2000) Expression patterns of Twist and Fgfr1, -2 and -3 in the developing mouse coronal suture suggest a key role for twist in suture initiation and biogenesis. Mech Dev 91: 341–345.D. JohnsonS. IsekiAO WilkieGM Morriss-Kay2000Expression patterns of Twist and Fgfr1, -2 and -3 in the developing mouse coronal suture suggest a key role for twist in suture initiation and biogenesis.Mech Dev91341345

45. O'Rourke MP, Tam PP (2002) Twist functions in mouse development. Int J Dev Biol 46: 401–413.MP O'RourkePP Tam2002Twist functions in mouse development.Int J Dev Biol46401413

46. Castanon I, Baylies MK (2002) A Twist in fate: evolutionary comparison of Twist structure and function. Gene 287: 11–22.I. CastanonMK Baylies2002A Twist in fate: evolutionary comparison of Twist structure and function.Gene2871122

47. Barnes RM, Firulli AB (2009) A twist of insight - the role of Twist-family bHLH factors in development. Int J Dev Biol 53: 909–924.RM BarnesAB Firulli2009A twist of insight - the role of Twist-family bHLH factors in development.Int J Dev Biol53909924

48. Cakouros D, Raices RM, Gronthos S, Glackin CA (2010) Twist-ing cell fate: mechanistic insights into the role of twist in lineage specification/differentiation and tumorigenesis. J Cell Biochem 110: 1288–1298.D. CakourosRM RaicesS. GronthosCA Glackin2010Twist-ing cell fate: mechanistic insights into the role of twist in lineage specification/differentiation and tumorigenesis.J Cell Biochem11012881298

49. Yang J, Mani SA, Donaher JL, Ramaswamy S, Itzykson RA, et al. (2004) Twist, a master regulator of morphogenesis, plays an essential role in tumor metastasis. Cell 117: 927–939.J. YangSA ManiJL DonaherS. RamaswamyRA Itzykson2004Twist, a master regulator of morphogenesis, plays an essential role in tumor metastasis.Cell117927939

50. Karreth F, Tuveson DA (2004) Twist induces an epithelial-mesenchymal transition to facilitate tumor metastasis. Cancer Biol Ther 3: 1058–1059.F. KarrethDA Tuveson2004Twist induces an epithelial-mesenchymal transition to facilitate tumor metastasis.Cancer Biol Ther310581059

51. Matukumalli LK, Lawley CT, Schnabel RD, Taylor JF, Allan MF, et al. (2009) Development and characterization of a high density SNP genotyping assay for cattle. PLoS One 4: e5350.LK MatukumalliCT LawleyRD SchnabelJF TaylorMF Allan2009Development and characterization of a high density SNP genotyping assay for cattle.PLoS One4e5350

52. Abecasis GR, Cherny SS, Cookson WO, Cardon LR (2002) Merlin-rapid analysis of dense genetic maps using sparse gene flow trees. Nat Genet 30: 97–101.GR AbecasisSS ChernyWO CooksonLR Cardon2002Merlin-rapid analysis of dense genetic maps using sparse gene flow trees.Nat Genet3097101

53. Rozen S, Skaletsky HJ (2000) Primer3 on the WWW for general users and for biologist programmers. Methods Molecular Biol 132: 365–386.S. RozenHJ Skaletsky2000Primer3 on the WWW for general users and for biologist programmers.Methods Molecular Biol132365386

54. Weckx S, Del-Favero J, Rademakers R, Claes L, Cruts M, et al. (2005) NovoSNP, a novel computational tool for sequence variation discovery. Genome Res 15: 436–442.S. WeckxJ. Del-FaveroR. RademakersL. ClaesM. Cruts2005NovoSNP, a novel computational tool for sequence variation discovery.Genome Res15436442