- Research article
- Open Access
Arachnomelia syndrome in Simmental cattle is caused by a homozygous 2-bp deletion in the molybdenum cofactor synthesis step 1 gene (MOCS1)
BMC Genetics volume 12, Article number: 11 (2011)
Arachnomelia syndrome is an autosomal recessive inherited disease in cattle. Affected calves die around birth and show malformations of the skeleton mainly affecting the legs, the spinal column and the skull. A number of arachnomelia syndrome affected Simmental calves were recently detected by a surveillance system of anomalies with a peak of more than 120 recorded cases in the year 2006. The causative mutation was previously mapped to a 9 cM-region on bovine chromosome 23. We herein report the fine-mapping and identification of the gene causing arachnomelia syndrome in Simmental cattle.
By using a dense set of markers, the arachnomelia syndrome linked region could be refined to 1.5 cM harbouring three protein coding genes. Comparative sequencing of these genes revealed a two-bp-deletion in the bovine MOCS1 gene resulting in a frame-shift and a premature termination codon. We genotyped affected calves and their ancestors and found that all affected were homozygous for the deletion whereas all carriers were heterozygous. Furthermore, cattle from the same population, but not directly related to known carriers mostly showed the wild type genotype.
MOCS1 encodes two proteins that are involved in the first synthesis step of molybdenum cofactor. A non functional sulfite-oxydase, one of the enzymes requiring molybdenum cofactor, leads to a similar pathology in Brown Swiss cattle. In combination the perfect association of the mutation with the phenotype and the obvious disruption of protein translation provide strong evidence for the causality of the MOCS1 mutation. Our results are the first example for an oligogenic lethal inherited disease in cattle. Furthermore, they show the potential involvement of sulfite metabolism in aberrant bone development.
The congenital arachnomelia syndrome (AS, OMIA Phene ID 139, Group 000059) in Holstein Friesian, Red Holstein and Simmental cattle was initially described by Rieck and Schade . The cardinal pathological changes are skeletal malformations of the skull, legs and the spinal column . The facial deformation leads to brachygnathia inferior and concave rounding of the maxilla forming a dent ('pointer-head'). The legs show abnormally thin diaphyses of the long bones and stiffened, hyperextended fetlocks ('spider-legs', dolichostenomelia), leading to frequent fractures of the metacarpus and metatarsus in the course of forced birth assistance. The vertebral formation failure results in kyphosis and scoliosis.
In the 1980s, the syndrome was dispersed in European Brown Swiss cattle by the use of American Brown Swiss sires  and in 2005 AS reappeared in Simmental. There were more than 150 confirmed cases in this breed with a peak in 2006 . It could be shown by pedigree analysis that AS is most probably a monogenic recessive inherited disease with an estimated allele frequency of 3% in the current cow population . Subsequently, AS was mapped to a 9 cM-region on bovine chromosome 23 by linkage analyses of families including 89 affected Simmental calves . An indirect gene test was developed and applied to several hundred sires that had a pedigree-based risk to be a carrier of AS (results are publicly available ). The availability of the indirect genetic test and the awareness of the disease by the artificial insemination centres and breeders led to a sharp decline in cases after 2007 . AS is the first oligogenic Mendelian inherited disease found in cattle, since it was mapped to a different location (on bovine chromosome 5) in Brown Swiss .
In this study, we present the fine mapping and identification of the putative disease causing mutation of AS in Simmental cattle thereby providing insights into the genetic basis of AS and the biology of bone development.
Fine mapping of the AS region
In a previous work we mapped the AS-condition in Simmental cattle to a ~9 cM-region. The corresponding region was covered with additional markers to further narrow down the AS locus. Published  as well as newly developed microsatellites were tested within the AS families. Finally, 11 informative additional microsatellites were chosen and the AS families were genotyped (Table 1, Additional file 1). The 86 cases available for fine mapping came from 4 half sib families. Therefore, the paternal and, subsequently, the maternal haplotype could unambiguously be deduced in most cases. 17-marker haplotypes were constructed for all genotyped cases. We analysed the minimal region without recombination and the region with complete homozygosity within the AS-affected animals. Both analyses led to the identical minimal candidate region of AS between the markers RM033 and LFL014 (Figure 1, Table 2). This region is about 0.7 cM or 1.5 Mb in length. These results excluded some prominent candidate genes that are located close to the region in question: RUNX2 and CUL7 that had been associated with diseases of the skeleton [8, 9].
Investigating candidate genes within the region carrying the AS mutation
The 1.5 Mb-core candidate region on chromosome 23 harboured one pseudogene, a TAF9 RNA polymerase II, one UAU-encoding tRNA gene, and three protein-coding genes, DAAM2, MOCS1, and LRFN2, that were annotated in BTAU4.0 (Figure 1). Based on their functional annotation none of these genes was a straightforward candidate with a known function in bone or cartilage development. The pseudogene is similar to TAF2G encoding an ubiquitously expressed transcription factor that is essential for cell viability . DAAM2 encodes the dishevelled associated activator of morphogenesis 2 that is involved in the wnt signalling pathway  implying a role in early cell development. Nevertheless, the gene function is not very well understood and there are no data from disease associations or animal models. MOCS1 encodes the molybdenum cofactor synthesis 1 protein that is involved in the synthesis of the molybdenum cofactor (Moco). Moco is required by at least three enzymes, sulfite-oxidase, xanthine-dehydrogenase, and aldehyde-oxidase. Known mutations within MOCS1 lead to an autosomal recessive disease, the molybdenum cofactor deficiency, with neurological symptoms in humans . The last gene, LRFN2 encoding the leucine-rich repeat and fibronectin type III domain-containing protein 2, is known to interact with different receptors in brain tissue and plays a role in t-cell development and hematopoiesis . Nevertheless, the protein seems to interact with the NMDA receptors  that, among other functions, are involved in bone development [15, 16]. Since there was no convincing candidate, we systematically sequenced the exonic and part of the intergenic, 5' and 3' regions of all three genes. The genomic organization of the bovine DAAM2 and MOCS1 was derived from the BTAU4.0 annotation, whereas the exons of LRFN2 were deduced from human and mouse EST and RNA data, since the first exon was not annotated. In total 22.4 kb were sequenced and twelve mutations, including one deletion, were detected within the three genes by comparative sequencing of one unrelated, one AS-carrier, and two AS-affected animals (Table 3).
A panel of additional animals was tested for these alleles either by direct sequencing, allele-specific PCR or PCR-RFLP assays. Excluding (either an AS-affected calf being not homozygote or a known carrier being not heterozygous for the mutation and the wt allele) genotypes were found for all mutations, except for the deletion within MOCS1.
Mutation of the bovine MOCS1 gene
The two bp-CA-deletion (Figure 2) is located at RNA (GI:261490660) position 1224-1225 (numbering starting at the A of the translation start codon) that is located in exon 11 of MOCS1 (c.1224-1225delCA). A genotyping system based on allele specific PCR was designed for the mutation and all cases as well as the family members available were genotyped (Table 4). As expected, all affected calves were homozygote for the mutation (del/del), whereas all parents were heterozygote (wt/del). As a control, cattle sampled randomly from the Bavarian Simmental population excluding first-degree relatives of known carriers were genotyped. None of them was homozygote for the c.1224-1225delCA mutation, but about 2.8% were heterozygote. The mutation was not observed in 120 cattle from other breeds (German Gelbvieh, Holstein Friesian, Belgian Blue and Braunvieh). We analysed the pedigrees of the heterozygous Simmental cattle from the random sample and found that, as previously observed for the affected calves, all could be traced back to SEMPER or his sire SENAT, confirming the hypothesis, that SENAT may be the founder and that most of the carriers inherited the mutation via SEMPER .
Finally, 67 sires that were known to be putative carriers of the mutation, either by pedigree analysis or by results from genetic testing , were genotyped. These sires include ROMEL and REXON, that produced a large number of daughters. All 67 sires were heterozygous for the deletion. DNA from the suspected founder, SEMPER, was not available, but DNA from one of his sons (SENDER) and brothers (SET) were genotyped. These were homozygous for the wild type allele. This result is in agreement with the finding, that all pedigrees available can be traced back to SEMPER via one of his daughters, but not to one of his sons. Two affected calves, that were considered as phenocopies since there was no connection of their pedigrees with SEMPER or SENAT [2, Table 4] were homozygote for the wild type allele.
MOCS1 is one of the rare genes known to express a polycistronic RNA. It produces two different enzymes (MOCS1A and MOCS1B) from non-overlapping open-reading-frames from a bicistronic transcript . The RNA is known to occur in different splice types in vertebrates including mouse, human, cattle and invertebrates . At least two splice types seem to exist in cattle, encoding either MOCS1A or a functional MOCS1B that is translated as a fusion with an inactive MOCS1A protein (Figure 3).
The reference sequence (GI:261490661) of the bovine MOCS1 is 633 amino acids in length and the C-terminus encoding the Moa-C domain is highly conserved within mammals (Figure 4). The 2 bp deletion leads to a predicted frame shift beginning at amino acid position 24 of MOCS1B. The altered protein is 73 amino acids in length and the complete MoaC domain is skipped. Therefore, the altered protein is unlikely to be functional.
The emergence of the bovine AS caused serious concern among European Simmental breeders in 2006, since the disease was transmitted by some top sires that had been extensively used. Therefore, we mapped the AS locus and an indirect gene test was made available in summer 2007 that was an effective tool to control the disease by genotyping sires with conspicuous pedigrees.
To explore the disease mechanism and to avoid the drawbacks of an indirect gene test, the causative mutation should be identified based on the initial mapping data. Additional regional markers were applied and the AS region was finally narrowed down to a 1.5 Mb-region containing five genes. Since none of them could convincingly be associated with the phenotype of AS, three genes were comparatively sequenced. As a result, exactly one obviously functional mutation in MOCS1 could be identified. Genotyping of a large panel of affected, obligate carrier and unrelated cattle provided additional evidence, that we identified the genuine AS mutation for Simmental cattle.
MOCS1 encodes a bicistronic RNA that is translated into MOCS1A and MOCS1B that are involved in the first step of the synthesis of Moco. They are analogous to the MoaC and MoaA enzymes of Escherichia coli and promote the synthesis of precursor Z (or cyclic pyranopterin monophosphate, cPMP) from a 5'-GTP derivate. Mutations of MOCS1 that lead to an inability of synthesising functional precursor Z cause Moco-deficiency type A (MIM 252150) in human, an autosomal recessive inherited disease that is mainly characterised by neurological symptoms due to brain dysmorphisms (loss of white matter in the CNS) . Nevertheless, mutations of an other gene leading to a disturbed sulphur metabolism, CBX, cause homocystinuria that is associated with different signs including abnormally developed phalanges (arachnodactyly) .
Homozygous MOCS1 knock-out mice show curly whiskers, die between day 1 and 11 of age and are smaller than wt mice, but do not show any other morphological abnormalities or significant changes in CNS . MOCS1-/- mice show the biochemical characteristics of sulfite and xanthine intoxication, e.g. elevated sulfite and xanthine levels and undetectable levels of uric acid in the urine. The main pathological mechanism seems to be mediated by the sulfite oxidase deficiency (SOD), since, probably due to redundancy, known isolated deficiencies of other molybdenum dependent enzymes (aldehyde oxidase and xandthine oxidoreductase) do not result in clinical symptoms in humans [21, 22]. The precise mechanisms of this toxicity are not completely understood. In mice, cell death appears to be triggered by elevated sulfite levels .
All splice forms of MOCS1 are found in a constant ratio in a wide range of tissues, probably reflecting its universal use . In cattle the main predilection site of Moco-deficiency is the developing bone. Even though bone malformations are not common signs of Moco-deficiency in human and mice, rarely microcephaly and a prominent forehead are observed in patients with Moco-  and SOD-deficiency .
Mice and children with Moco-deficiency are usually born alive. It was postulated that sulfite is ameliorated in affected embryos by maternal clearance and therefore damage seems to be modulated till birth . Most probably this is not the case in cattle due to the non-invasive implantation. Therefore, (bone-)tissue is not protected from damage from elevated sulfite levels.
Finally, the causative role of MOCS1 in bovine AS is strongly supported by the recent discovery of a SUOX mutation leading to a similar phenotype in Brown Swiss cattle . Therefore, AS is the first known example for an oligogenic inherited disease in cattle caused by mutations within genes of an enzyme and its cofactor. These results provide impetus for investigating the possible role of elevated sulfite levels in aberrant bone development.
Moco can be treated by purified precursor Z in mouse . In cattle it will not be possible to treat Moco in near future, since the critical insult occurs during fetal development and affected calves usually are not viable postnatally. Instead, AS is controlled by excluding carriers of the mutation from breeding. Replacing the indirect by a direct gene test will accelerate the eradication of AS from the Simmental population. Since the direct gene test is not only more efficient but also easier and cheaper than the indirect gene test, it will become feasible to test larger numbers of breeding cows. This is important due to the comparatively high frequency of carriers in the current cow population.
Our results strongly support the hypothesis that a missing molybdenum cofactor, leading to non functional Moco-dependent enzymes, namely sulfite oxidase, lead to aberrant bone development in bovine AS.
There are 4 lines of evidence for the mutation of MOCS1 being causative for AS. 1) AS is an autosomal recessive inherited lethal disease. Therefore, it could be expected that a loss of function mutation is the cause of AS. 2) The genotyping results of additional animals show complete concordance of the deletion genotype with the AS phenotype. 3) Mutations of MOCS1 or the dependant SUOX can also lead to bone malformations in human. 4) The corresponding mutation of SUOX leads to a similar phenotype in Brown Swiss cattle.
Collection and DNA extraction
The collection of AS-affected Simmental calves had been described previously . In short, affected calves were identified in the course of the Bavarian surveillance programme for inherited congenital malformations. Only calves with AS status confirmed by pathological investigation were included in the study. DNA from tissue was prepared using the DNA tissue kit (Qiagen). DNA from blood was prepared using a blood kit (PEQLAB Biotechnology GmbH). Pedigrees were extracted from the joint German and Austrian Simmental pedigree data used for routine genetic evaluation.
For determination of allele frequencies Simmental cattle were collected from the performance testing station reflecting the breeding population of the years 2002 - 2008 when the AS mutation was spread in the population by some prominent sires . Cattle with parents being known carriers were excluded from the analyses. Holstein Friesian and Braunvieh cattle were male calves bought from local breeders in Bavaria.
DNA from pedigree members were genotyped with 17 microsatellites covering the region 13-40 cM on chromosome 23 (Table 1, Additional File 1). These were selected from the map published by Ihara  or developed based on the bovine RefSeq assembly BTAU4.0 (Table 1). One primer of each pair was fluorescently labelled and microsatellites were amplified in 10 μl-reactions and separated on an ABI Prism 310 Genetic Analyzer (Applied Biosystems) as described . Fragment sizes were measured and allele tables were generated using the GENEMAPPER (V. 3.0, Applied Biosystems) software.
Exonic as well as part of the intergenic, 5'- and 3'-regions of three genes were amplified in 10 μl reaction volume using 2 μl of genomic DNA, 0.05 μM of each primer (Supplement 1), 50 μM of each dNTP, 0.5 units of HotStar-taq polymerase (Qiagen, Hilden, Germany), and reaction buffer containing 1.5 mM MgCl2. Cycling-conditions on a Biometra T gradient 96-well thermocycler were: 15 min at 95°C, [0.5 min at 95°C, 1 min at 60°C, 0.5 min at 72°C]35 × and 10 min, 60°C final extension. Before direct sequencing PCR-reactions were purified using the QIAquick PCR Purification Kit. Sequencing was performed using BigDye V 3.0 terminator cycle sequencing kit (Applied Biosystems, Foster City, CA). Sequencing analysis was run on an ABI PRISM® 310 genetic analyser. Single nucleotide variation was analysed using PolyPhred [28–30].
Allele specific PCR genotyping system
Primers (5'-cctgacatgaacaggggaac-3', 5'- ccaggtgggaactgaagtgt-3', 5'-tggagaaactcactctgtgtctc-3', 5'- gtgttcaggttctgagacagagt-3') were designed in order to amplify the two alleles of the 2-bp deletion in an tetra-ARMS reaction. 1.5 pmol of each primer were used to amplify 30 ng genomic DNA in a 10 μl PCR reaction with 0.5 units of HotStar-taq polymerase (Qiagen, Hilden, Germany) under standard conditions. Cycling-conditions on a Biometra T gradient 96-well thermocycler were: 15 min at 95°C, [0.5 min at 95°C, 40 sec at 60°C, 1 min at 72°C]29×. The fragments were separated on a 2% agarose gel in TAE-buffer at 150 V for 40 min. The PCR reaction resulted in fragments of 208 bp and 142 bp in length for the wt and AS allele, respectively.
Rieck GW, Schade W: Arachnomelia (spider limbs), a new hereditary fatal malformation syndrome of cattle. [in German]. Dtsch Tierärztl Wochenschr. 1975, 82: 342-347.
Buitkamp J, Luntz B, Emmerling R, Reichenbach HD, Weppert M, Schade B, Meier N, Götz KU: Syndrome of arachnomelia in Simmental cattle. BMC Vet Res. 2008, 4: 39-10.1186/1746-6148-4-39.
König H, Galliard C, Chavaz J, Hunziker F, Tontis A: Prüfung von Schweizer Braunvieh-Bullen auf das vererbte Syndrom der Arachnomelie und Arthrogrypose (SAA) durch Untersuchung der Nachkommen im Fetalstadium. Tierärztl Umsch. 1987, 42: 692-697.
Buitkamp J, Kühn C, Semmer J, Götz KU: Assignment of the locus for arachnomelia syndrome to bovine chromosome 23 in Simmental cattle. Anim Genet. 2009, 40: 894-899. 10.1111/j.1365-2052.2009.01933.x.
Bayerische Zuchtwert-Informationen, BaZI-Rind. [http://www.lfl.bayern.de/bazi-rind/]
Drögemüller C, Rossi M, Gentile A, Testoni S, Jorg H, Stranzinger G, Drögemüller M, Glowatzki-Mullis ML, Leeb T: Arachnomelia in Brown Swiss cattle maps to chromosome 5. Mamm Genome. 2009, 20: 53-59.
Ihara N, Takasuga A, Mizoshita K, Takeda H, Sugimoto M, Mizoguchi Y, Hirano T, Itoh T, Watanabe T, Reed KM: A comprehensive genetic map of the cattle genome based on 3802 microsatellites. Genome Res. 2004, 14: 1987-1998. 10.1101/gr.2741704.
Huber C, Dias-Santagata D, Glaser A, O'Sullivan J, Brauner R, Wu K, Xu X, Pearce K, Wang R, Uzielli ML: Identification of mutations in CUL7 in 3-M syndrome. Nat Genet. 2005, 37: 1119-1124. 10.1038/ng1628.
Mundlos S, Otto F, Mundlos C, Mulliken JB, Aylsworth AS, Albright S, Lindhout D, Cole WG, Henn W, Knoll JH: Mutations involving the transcription factor CBFA1 cause cleidocranial dysplasia. Cell. 1997, 89: 773-779. 10.1016/S0092-8674(00)80260-3.
Frontini M, Soutoglou E, Argentini M, Bole-Feysot C, Jost B, Scheer E, Tora L: TAF9b (formerly TAF9L) is a bona fide TAF that has unique and overlapping roles with TAF9. Mol Cell Biol. 2005, 25: 4638-4649. 10.1128/MCB.25.11.4638-4649.2005.
Katoh M: Networking of WNT, FGF, Notch, BMP, and Hedgehog signaling pathways during carcinogenesis. Stem Cell Rev. 2007, 3: 30-38. 10.1007/s12015-007-0006-6.
Reiss J, Cohen N, Dorche C, Mandel H, Mendel RR, Stallmeyer B, Zabot MT, Dierks T: Mutations in a polycistronic nuclear gene associated with molybdenum cofactor deficiency. Nat Genet. 1998, 20: 51-53. 10.1038/1706.
Castellanos A, Lang G, Frampton J, Weston K: Regulation of erythropoiesis by the neuronal transmembrane protein Lrfn2. Exp Hematol. 2007, 35: 724-734. 10.1016/j.exphem.2007.02.004.
Wang CY, Chang K, Petralia RS, Wang YX, Seabold GK, Wenthold RJ: A novel family of adhesion-like molecules that interacts with the NMDA receptor. J Neurosci. 2006, 26: 2174-2183. 10.1523/JNEUROSCI.3799-05.2006.
Merle B, Itzstein C, Delmas PD, Chenu C: NMDA glutamate receptors are expressed by osteoclast precursors and involved in the regulation of osteoclastogenesis. J Cell Biochem. 2003, 90: 424-436. 10.1002/jcb.10625.
Itzstein C, Espinosa L, Delmas PD, Chenu C: Specific antagonists of NMDA receptors prevent osteoclast sealing zone formation required for bone resorption. Biochem Biophys Res Commun. 2000, 268: 201-209. 10.1006/bbrc.2000.2097.
Gray TA, Nicholls RD: Diverse splicing mechanisms fuse the evolutionarily conserved bicistronic MOCS1A and MOCS1B open reading frames. Rna. 2000, 6: 928-936. 10.1017/S1355838200000182.
Mize C, Johnson JL, Rajagopalan KV: Defective molybdopterin biosynthesis: clinical heterogeneity associated with molybdenum cofactor deficiency. J Inherit Metab Dis. 1995, 18: 283-290. 10.1007/BF00710416.
Sebastio G, Sperandeo MP, Panico M, de Franchis R, Kraus JP, Andria G: The molecular basis of homocystinuria due to cystathionine beta-synthase deficiency in Italian families, and report of four novel mutations. Am J Hum Genet. 1995, 56: 1324-1333.
Lee HJ, Adham IM, Schwarz G, Kneussel M, Sass JO, Engel W, Reiss J: Molybdenum cofactor-deficient mice resemble the phenotype of human patients. Hum Mol Genet. 2002, 11: 3309-3317. 10.1093/hmg/11.26.3309.
Johnson JL, Duran M: Molybdenum cofactor deficiency and isolated sulfite oxidase deficiency. The metabolic and molecular bases of inherited disease. Edited by: Scriver C, Beaudet A, Sly W, Valle D. 2001, New York: McGraw-Hill, 3163-3177.
Simmonds HA, Reiter S, Nishino T: Hereditary xanthinuria. The metabolic and molecular bases of inherited disease. Edited by: Scriver C, Beaudet A, Sly W, Valle D. 1995, New York: McGraw-Hill, 1781-1797.
Reiss J, Bonin M, Schwegler H, Sass JO, Garattini E, Wagner S, Lee HJ, Engel W, Riess O, Schwarz G: The pathogenesis of molybdenum cofactor deficiency, its delay by maternal clearance, and its expression pattern in microarray analysis. Mol Genet Metab. 2005, 85: 12-20. 10.1016/j.ymgme.2005.01.008.
Gross-Hardt S, Reiss J: The bicistronic MOCS1 gene has alternative start codons on two mutually exclusive exons. Mol Genet Metab. 2002, 76: 340-343. 10.1016/S1096-7192(02)00100-2.
Endres W, Shin YS, Gunther R, Ibel H, Duran M, Wadman SK: Report on a new patient with combined deficiencies of sulphite oxidase and xanthine dehydrogenase due to molybdenum cofactor deficiency. Eur J Pediatr. 1988, 148: 246-249. 10.1007/BF00441412.
Hobson EE, Thomas S, Crofton PM, Murray AD, Dean JC, Lloyd D: Isolated sulphite oxidase deficiency mimics the features of hypoxic ischaemic encephalopathy. Eur J Pediatr. 2005, 164: 655-659. 10.1007/s00431-005-1729-5.
Drögemüller C, Tetens J, Sigurdsson S, Gentile A, Testoni S, Lindblad-Toh K, Leeb T: Identification of the bovine arachnomelia mutation by massively parallel sequencing implicates sulfite oxidase (SUOX) in bone development. PLoS Genet. 2010, 6: e1001079-
Ewing B, Hillier L, Wendl MC, Green P: Base-calling of automated sequencer traces using phred. I. Accuracy assessment. Genome Res. 1998, 8: 175-185.
Gordon D, Abajian C, Green P: Consed: a graphical tool for sequence finishing. Genome Res. 1998, 8: 195-202.
Nickerson DA, Tobe VO, Taylor SL: PolyPhred: automating the detection and genotyping of single nucleotide substitutions using fluorescence-based resequencing. Nucleic Acids Res. 1997, 25: 2745-2751. 10.1093/nar/25.14.2745.
Rizzi M, Schindelin H: Structural biology of enzymes involved in NAD and molybdenum cofactor biosynthesis. Curr Opin Struct Biol. 2002, 12: 709-720. 10.1016/S0959-440X(02)00385-8.
We thank the farmers for sending their calves and blood samples and the AI-stations for providing semen samples for genetic analyses. Tissue samples and phenotypic data were provided by Drs. Schade and Meier from the TGD, Grub.
JB drafted the manuscript, designed the mapping and sequencing strategy and analyzed the data. Genotyping and sequencing was done by JS for wet-lab portions.
K-UG participated in study design and coordination and critically revised the manuscript. All authors read and approved the final manuscript.