|
|
|
|
Published online before print
August 9, 2006, 10.1101/gr.5268806 Genome Res. 16:1084-1090, 2006 ©2006 by Cold Spring Harbor Laboratory Press; ISSN 1088-9051/06 $5.00
Letter An ~140-kb deletion associated with feline spinal muscular atrophy implies an essential LIX1 function for motor neuron survival1Laboratory of Comparative Medical Genetics, Department of Microbiology & Molecular Genetics, College of Veterinary Medicine, Michigan State University, East Lansing, Michigan 48824, USA; 2Laboratory of Genomic Diversity, National Cancer InstituteFrederick, Frederick, Maryland 21702, USA; 3Institute of Human Genetics, Institute of Genetics, and Center for Molecular Medicine Cologne, University of Cologne, 50931 Cologne, Germany; 4National Center for Biotechnology Information, National Institutes of Health, Department of Health and Human Services, Bethesda, Maryland 20894, USA; 5Department of Veterinary Integrative Biosciences, College of Veterinary Medicine and Biomedical Sciences, Texas A&M University, College Station, Texas 77843, USA; 6Department of Biochemistry & Molecular Biology, Michigan State University, East Lansing, Michigan 48824, USA
The leading genetic cause of infant mortality is spinal muscular atrophy (SMA), a clinically and genetically heterogeneous group of disorders. Previously we described a domestic cat model of autosomal recessive, juvenile-onset SMA similar to human SMA type III. Here we report results of a whole-genome scan for linkage in the feline SMA pedigree using recently developed species-specific and comparative mapping resources. We identified a novel SMA gene candidate, LIX1, in an ~140-kb deletion on feline chromosome A1q in a region of conserved synteny to human chromosome 5q15. Though LIX1 function is unknown, the predicted secondary structure is compatible with a role in RNA metabolism. LIX1 expression is largely restricted to the central nervous system, primarily in spinal motor neurons, thus offering explanation of the tissue restriction of pathology in feline SMA. An exon sequence screen of 25 human SMA cases, not otherwise explicable by mutations at the SMN1 locus, failed to identify comparable LIX1 mutations. Nonetheless, a LIX1-associated etiology in feline SMA implicates a previously undetected mechanism of motor neuron maintenance and mandates consideration of LIX1 as a candidate gene in human SMA when SMN1 mutations are not found.
The spinal muscular atrophies (SMAs) are a genetically heterogeneous group of disorders that vary in clinical severity, from lethal in infancy to onset of mild weakness in adulthood, but all are characterized by neurogenic muscle atrophy due to degeneration of lower motor neurons of the spinal cord (Talbot and Davies 2001
The high incidence of human SMA is mainly due to a human-specific (Rochette et al. 2001
We recently described a domestic cat model of SMA exhibiting autosomal recessive, juvenile-onset skeletal muscle atrophy and weakness with onset at ~12 wk of age due to lower motor neuron loss (He et al. 2005
Whole-genome scan and linkage analysis We genotyped 23 privately owned and 45 purpose-bred members of a domestic cat SMA pedigree (Fig. 1) with previously mapped, highly polymorphic short tandem repeat (STR) markers (Menotti-Raymond et al. 2003a
Comparison of the published order of feline markers (Murphy et al. 2000
Deletion abrogating LIX1 and LNPEP expression
Because the LIX1 microsatellite was deleted in affected cats, we avoided the usual challenge of choosing and sequencing positional candidate genes for mutations, but were instead faced with the challenge of determining the endpoints of the mutation and supplying further proof that the deletion disrupts LIX1 expression. Fifty-two additional sequence-tagged sites (STS) in LIX1 exons and at regular genomic intervals 3' of LIX1 were analyzed in affected cat DNA, including the downstream genes LNPEP, LRAP, ARTS-1, and CAST. A deletion was confirmed, and the breakpoints were narrowed to intron 3 of LIX1 and intron 1 of LNPEP (Fig. 2A). LIX1 and LNPEP are transcribed from opposite strands, their most 3' exons adjacent, so the data indicated loss of exons 46 of LIX1 and all but exon 1 of LNPEP. Based on comparison to the homologous region in the dog genome by BLAT analysis (http://genome.ucsc.edu/cgi-bin/hgBlat?), we estimated the deletion at the cat SMA disease locus to be ~140 kbp. We identified the precise deletion breakpoints by sequencing the PCR product amplified from affected cat genomic DNA with primers flanking the deletion. Two copies of the short sequence, AGTTTA, flanked the deletion site, a finding compatible with a deletion mechanism of nonhomologous recombination (Woods-Samuels et al. 1991
We confirmed that all 31 affected cats were homozygous for the deleted allele and all 28 obligate carriers were heterozygous by multiplex PCR assay amplifying products from each of the deletion breakpoints on the normal and deleted alleles of FCA A1q (Fig. 2B). One hundred fifty-two clinically normal purebred Maine coon cats, the breed from which the feline SMA family was derived, were heterozygous or homozygous for the normal allele. RTPCR analysis demonstrated LIX1 and LNPEP mRNA expression in cervical spinal ventral horn gray matter of genetically normal and heterozygous carrier cats, but no expression of either gene was detected in the affected cats (Fig. 2C). Aberrant RNA transcripts were not observed in affected cat RNA by Northern blotting (data not shown). These data thus also confirmed that feline SMA in this kindred is an autosomal recessive trait caused by disruption of the LIX1/LNPEP locus.
Characterization of LIX1
The nucleotide sequence of LIX1 cDNA amplified from cat spinal cord demonstrated an open reading frame (ORF) of 846 bp (GenBank accession no. DQ250154 [GenBank] ). The deduced amino acid sequence was 99%, 97%, 97%, 92%, 85%, 59%, and 51% identical to those of dog, cow, human, mouse, chicken, Danio rerio, and Drosophila melanogaster, respectively, and amino acid substitutions among these species are highly conservative. Protein sequence conservation across this broad range of species suggested that LIX1 function is essential. Of particular note was that computerized secondary structure analysis predicted amino acid residues 2999 of LIX1 fold into a double-stranded RNA-binding domain, a domain common in proteins of RNA metabolism or transport and previously predicted in CG13139, the Drosophila LIX1 ortholog (Saunders and Barber 2003
We observed additional conservation outside of the protein-coding sequence in LIX1 mRNA compatible with conservation of essential post-transcriptional cis regulatory elements. Sequence identity was highly concentrated in the ~400 bp before the start translation codon (5' UTR), where the cat sequence was 97%, 95%, 95%, 88%, 69%, 52%, and 38% identical to the dog, cow, human, mouse, frog, chicken, and fruitfly sequences, respectively. In each of the mammals, the degree of nucleotide sequence conservation in the 5' UTR was identical to or greater than that found in the ORF of LIX1. Included in this conserved segment are sequences predicted (Zuker 2003
The LIX1 3' UTR sequence is also highly conserved, further suggesting that post-transcriptional regulation may be an important element of LIX1 expression. The cat 3' UTR sequence was 83%, 75%, 75%, 37%, and 39% identical to the dog, cow, human, mouse, and chicken, respectively, after excluding an Alu-J element from the human and two ruminant-specific SINEs from the cow sequences. The aligned sequences were 2.652.93 kb, except for the mouse sequence, which was 1.98 kb. In all but the mouse and chicken, there were four similarly placed poly(A) addition signals. We demonstrated the alternative use of two poly(A) addition signals in LIX1 mRNA isolated from cat lumbar ventral horn by RTPCR and 3'-rapid amplification of cDNA ends (RACE) with polyadenylation beginning 1186 and 2642 bp 3' of the stop translation codon, respectively. Among the five AAUAAA sequences in the cat LIX1 3' UTR, the two we observed being used were those that were followed by a GU-rich region immediately after the polyadenylation site, conforming thereby to the bipartite consensus for efficiently used polyadenylation signals (Proudfoot 1991
Survey for human LIX1 mutations
To assess the possibility that mutations in LIX1 might also occur in cases of human SMA, we sequenced the LIX1-coding regions of 25 unrelated SMA patients (23 exhibited the type III phenotype) previously determined to be free of SMN1 deletions or subtle mutations (Wirth et al. 1999
We report here the molecular basis of type III SMA without SMN mutation in a cat model, the determination of which is the first discovery of a novel disease gene in the domestic cat by whole-genome scan and positional reasoning. Identification of LIX1 as the putative feline SMA gene and having function that is apparently required for motor neuron maintenance are findings of potential relevance to human SMA. This success demonstrates the current maturity of gene-mapping resources available for this species (Murphy et al. 2000
To date, there are reports of ~250 spontaneously occurring pathologies of likely genetic origin in cats (http://omia.angis.org.au/), and almost all have human counterparts. With the current sophistication of medical monitoring for pet cats, of which there are ~90 million in the United States alone (http://www.appma.org/press_industrytrends.asp), and the inbred nature of many cat breeds, new disease entities and new cases of previously described pathologies continue to be ascertained. With the development of genomic resources in the cat and the application of complementary comparative tools developed in other species, the domestic cat is emerging as a promising resource of pheno-typically defined genetic variation of biomedical significance (O'Brien et al. 2002
Strictly speaking, the question still remains as to which of the deleted genes is the SMA disease gene or whether loss of both, or perhaps of an unidentified nonprotein coding sequence is required to produce the disease. Based on reported functionality and/or expression profiles, we suggest that LIX1 alone is the disease gene. LNPEP is an aminopeptidase most highly expressed in placenta, heart, kidney, and small intestine that is variously secreted or a type II membrane protein, both forms having activity that degrade various peptides, including hormones such as oxytocin and vasopressin (Laustsen et al. 1997
While LNPEP does not present a compelling case for significant expression or essential function in spinal cord, LIX1 expression appears to be highest in the spinal cord of humans (http://genome.ucsc.edu/cgi-bin/hgGene?hgsid=70251085&db=hg18&hgg_gene=NM_153234). LIX1 was first reported as transiently expressed in the chicken developing limb bud and facial primordia, the gene identifier standing for limb expression 1 (Swindell et al. 2001
The high evolutionary conservation of LIX1 suggests that function of the gene product is essential, and the cat SMA phenotype caused by LIX1 disruption points to an important role in motor neuron development and/or maintenance. While LIX1 function remains to be determined, this putative disease gene identified by investigation of the cat SMA model exhibits an expression profile more restricted to the site of human and cat pathology than is the expression of SMN1, the most commonly altered gene in human SMA phenotypes. If the predicted double-stranded RNA-binding domain is a functional domain in this protein, it suggests that LIX1 may participate in RNA metabolism or transport, as does SMN (Meister et al. 2002
Animals and patients Privately owned, purebred Maine coon cats were investigated as previously described (He et al. 2005
Genotyping
Linkage analysis
Radiation hybrid mapping
Mutation and gene expression analysis
Gene expression was analyzed by RTPCR as previously described (He et al. 2005
The LIX1 mRNA sequence was determined by standard automated dideoxy chain termination cycle sequencing methods of RTPCR products amplified from cat spinal cord RNA. We assembled overlapping sequences and carried out sequence analysis with the Lasergene suite of programs (DNAstar). Northern blots were prepared by standard protocols (Sambrook and Russell 2001
Primary and secondary protein structure prediction
Human LIX1 sequencing
RTPCR analysis of RNA isolated from EBV-transformed lymphoblastoid cell lines of SMA patients was performed using the primers 5'-GATCCGGCTCTAGTCTTCAAA-3 and 5'-CCT CAGGGCCATTCGTAGC-3' in exons 1 and 6, respectively. Positive controls included RTPCR amplification of GAPDH and SFRS10 (Htra2-
We thank The Broad Institute and Agencourt Bioscience Corporation for public access to the 2x sequence assembly of the domestic cat genome. This work was supported in part by grants from the Association Française contre les Myopathies, National Institutes of Health (HD39888), Michigan State University Companion Animal Fund, Deutsche Forschungsgemeinschaft (Wi-945-12), Center for Molecular Medicine Cologne (TV98), and the Intramural Programs of the NIH, NCI, and NLM. The content of this publication does not necessarily reflect the views or policies of the Department of Health and Human Services, nor does mention of trade names imply endorsement by the U.S. Government.
7 These two authors contributed equally to this work.
E-mail fyfe{at}cvm.msu.edu; fax (517) 353-8957. Supplemental material is available online at www.genome.org. The sequence data from this study has been submitted to GenBank under accession no. DQ250154. Article published online before print. Article and publication date are at http://www.genome.org/cgi/doi/10.1101/gr.5268806.
Agarwala R., Applegate D.L., Maglott D., Schuler G.D., Schäffer A.A. 2000. A fast and scalable radiation hybrid map construction and integration strategy. Genome Res. 10: 350364. Altschul S.F., Madden T.L., Schäffer A.A., Zhang J., Zhang Z., Miller W., Lipman D.J. 1997. Gapped BLAST and PSI-BLAST: A new generation of protein database search programs. Nucleic Acids Res. 25: 33893402. Applegate D., Bixby R., Chvátal V., Cook W. 1998. On the solution of traveling salesman problems. In Documenta Math. III pp. 645656. International Congress of Mathematics III. Bassell G.J. and Kelic S. 2004. Binding proteins for mRNA localization and local translation, and their dysfunction in genetic neurological disease. Curr. Opin. Neurobiol. 14: 574581.[CrossRef][Medline] Ben-Dor A. and Chor B. 1997. On constructing radiation hybrid maps. J. Comput. Biol. 4: 517533.[Medline] Coovert D.D., Le T.T., McAndrew P.E., Strasswimmer J., Crawford T.O., Mendell J.R., Coulson S.E., Androphy E.J., Prior T.W., Burghes A.H.M. 1997. The survival motor neuron protein in spinal muscular atrophy. Hum. Mol. Genet. 6: 12051214. Fishelson M. and Geiger D. 2004. Optimizing exact genetic linkage computations. J. Comput. Biol. 11: 263275.[CrossRef][Medline] He Q., Lowrie C., Shelton G.D., Castellani R.J., Menotti-Raymond M., Murphy W., O'Brien S.J., Swanson W.F., Fyfe J.C. 2005. Inherited motor neuron disease in domestic cats: A model of spinal muscular atrophy. Pediatr. Res. 57: 324330.[CrossRef][Medline] Keller S.R., Davis A.C., Clairmont K.B. 2002. Mice deficient in the insulin-regulated membrane aminopeptidase show substantial decreases in glucose transporter GLUT4 levels but maintain normal glucose homeostasis. J. Biol. Chem. 277: 1767717686. Koh I.Y.Y., Eyrich V.A., Marti-Renom M.A., Przybylski D., Madhusudhan M.S., Eswar N., Graña O., Pazos F., Valencia A., Sali A. et al. 2003. EVA: Evaluation of protein structure prediction servers. Nucleic Acids Res. 31: 33113315. Laustsen P.G., Rasmussen T.E., Petersen K., Pedraza-Diaz S., Moestrup S.K., Gliemann J., Sottrup-Jensen L., Kristensen T. 1997. The complete amino acid sequence of human placental oxytocinase. Biochim. Biophys. Acta 1352: 17.[Medline] Lefebvre S., Burglen L., Reboullet S., Clermont O., Burlet P., Viollet L., Benichou B., Cruaud C., Millasseau P., Zeviani M. et al. 1995. Identification and characterization of a spinal muscular atrophy-determining gene. Cell 80: 155165.[CrossRef][Medline] Lefebvre S., Burlet P., Liu Q., Bertrandy S., Clermont O., Munnich A., Dreyfuss G., Melki J. 1997. Correlation between severity and SMN protein level in spinal muscular atrophy. Nat. Genet. 16: 265269.[CrossRef][Medline] Lin L., Faraco J., Li R., Kadotani H., Rogers W., Lin X., Qiu X., de Jong P.J., Nishino S., Mignot E. 1999. The sleep disorder canine narcolepsy is caused by a mutation in the hypocretin (orexin) receptor 2 gene. Cell 98: 365376.[CrossRef][Medline] Lorson C.L., Hahnen E., Androphy E.J., Wirth B. 1999. A single nucleotide in the SMN gene regulates splicing and is responsible for spinal muscular atrophy. Proc. Natl. Acad. Sci. 96: 63076311. Meister G., Eggert C., Fischer U. 2002. SMN-mediated assembly of RNPs: A complex story. Trends Cell Biol. 12: 472478.[CrossRef][Medline] Mendell J.T. and Dietz H.C. 2001. When the message goes awry: Disease-producing mutations that influence mRNA content and performance. Cell 107: 411414.[CrossRef][Medline] Menotti-Raymond M., David V.A., Agarwala R., Schäffer A.A., Stephens R., O'Brien S.J., Murphy W.J. 2003a. Radiation hybrid mapping of 304 novel microsatellites in the domestic cat genome. Cytogenet. Genome Res. 102: 272276.[CrossRef][Medline] Menotti-Raymond M., David V.A., Roelke M.E., Chen Z.Q., Menotti K.A., Sun S., Schäffer A.A., Tomlin J.F., Agarwala R., O'Brien S.J. et al. 2003b. Second-generation integrated genetic linkage/radiation hybrid maps of the domestic cat (Felis catus). J. Hered. 94: 95106. Moeller C., Yaylaoglu M.B., Alvarez-Bolado G., Thaller C., Eichele G. 2002. Murine Lix1, a novel marker for substantia nigra, cortical layer 5, and hindbrain structures. Brain Res. Gene Expr. Patterns 1: 199203.[CrossRef][Medline] Monani U.R., Lorson C.L., Parsons D.W., Prior T.W., Androphy E.J., Burghes A.H.M., McPherson J.D. 1999. A single nucleotide difference that alters splicing patterns distinguishes the SMA gene SMN1 from the copy gene. SMN2. Hum. Mol. Genet. 8: 11771183. Monani U.R., Coovert D.D., Burghes A.H.M. 2000. Animal models of spinal muscular atrophy. Hum. Mol. Genet. 9: 24512457. Munsat T.L. and Davies K.E. 1992. The International SMA consortium meeting. Neuromuscul. Disord. 2: 423428.[Medline] Murphy W.J., Mennotti-Raymond M., Lyons L.A., Thompson M.A., O'Brien S.J. 1999. Development of a feline whole-genome radiation hybrid panel and comparative mapping of human chromosome 12 and 22 loci. Genomics 57: 18.[CrossRef][Medline] Murphy W.J., Sun S., Chen Z., Yuhki N., Hirschmann D., Menotti-Raymond M., O'Brien S.J. 2000. A radiation hybrid map of the cat genome: Implications for comparative mapping. Genome Res. 10: 691702. O'Brien S.J., Menotti-Raymond M., Murphy W.J., Yuhki N. 2002. The Feline Genome Project. Annu. Rev. Genet. 36: 657686.[CrossRef][Medline] Ogino S. and Wilson R.B. 2002. Genetic testing and risk assessment for spinal muscular atrophy (SMA). Hum. Genet. 111: 477500.[CrossRef][Medline] Pearn J. 1978. Incidence, prevalence, and gene frequency studies of chronic childhood spinal muscular atrophy. J. Med. Genet. 15: 409428.[Abstract] Proudfoot N. 1991. Poly(A) signals. Cell 64: 671674.[CrossRef][Medline] Rochette C.F., Gilbert N., Simard L.R. 2001. SMN gene duplication and the emergence of the SMN2 gene occurred in distinct hominids: SMN2 is unique to Homo sapiens . Hum. Genet. 108: 255266.[CrossRef][Medline] Sambrook J. and Russell D.W. 2001. Molecular cloning: A laboratory manual. 3rd ed. Cold Spring Harbor Laboratory Press, Cold Spring Harbor, NY. Saunders L.R. and Barber G.N. 2003. The dsRNA binding protein family: Critical roles, diverse cellular functions. FASEB J. 17: 961983. Stevens G., Yawitch T., Rodda J., Verhaart S., Krause A. 1999. Different molecular basis for spinal muscular atrophy in South African black patients. Am. J. Med. Genet. 86: 420426.[CrossRef][Medline] Swanson W.F. 2006. Application of assisted reproduction for population management in felids: The potential and reality for conservation of small cats. Theriogenology 66: 4958.[CrossRef][Medline] Swindell E.C., Moeller C., Thaller C., Eichele G. 2001. Cloning and expression analysis of chicken Lix1, a founding member of a novel gene family. Mech. Dev. 109: 405408.[CrossRef][Medline] Talbot K. and Davies K.E. 2001. Spinal muscular atrophy. Semin. Neurol. 21: 189197. Wirth B. 2000. An update of the mutation spectrum of the survival motor neuron gene (SMN1) in autosomal recessive spinal muscular atrophy (SMA). Hum. Mutat. 15: 228237.[CrossRef][Medline] Wirth B., Schmidt T., Hahnen E., Rudnik-Schoneborn S., Krawczak M., Muller-Myhsok B., Schonling J., Zerres K. 1997. De novo rearrangements found in 2% of index patients with spinal muscular atrophy: Mutational mechanisms, parental origin, mutation rate, and implications for genetic counseling. Am. J. Hum. Genet. 61: 11021111.[CrossRef][Medline] Wirth B., Herz M., Wetter A., Moskau S., Hahnen E., Rudnik-Schoneborn S., Wienker T., Zerres K. 1999. Quantitative analysis of survival motor neuron copies: Identification of subtle SMN1 mutations in patients with spinal muscular atrophy, genotype-phenotype correlation, and implications for genetic counseling. Am. J. Hum. Genet. 64: 13401356.[CrossRef][Medline] Wirth B., Brichta L., Schrank B., Lochmüller H., Blick S., Baasner A., Heller R. 2006. Mildly affected patients with spinal muscular atrophy are partially protected by an increased SMN2 copy number. Hum. Genet. 119: 422428.[CrossRef][Medline] Woods-Samuels P., Kazazian Jr. H.H., Antonarakis S.E. 1991. Nonhomologous recombination on the human genome: Deletions in the human factor VIII gene. Genomics 10: 94101.[CrossRef][Medline] Zhang Z., Schwartz S., Wagner L., Miller W. 2000. A greedy algorithm for aligning DNA sequences. J. Comput. Biol. 7: 203214.[CrossRef][Medline] Zuker M. 2003. Mfold web server for nucleic acid folding and hybridization prediction. Nucleic Acids Res. 31: 34063415.
Received February 27, 2006; accepted in revised format June 6, 2006. This article has been cited by other articles:
|
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||