- Research article
- Open Access
The wooly mutation (wly) on mouse chromosome 11 is associated with a genetic defect in Fam83g
BMC Research Notesvolume 6, Article number: 189 (2013)
Mice homozygous for the spontaneous wooly mutation (abbreviated wly) are recognized as early as 3–4 weeks of age by the rough or matted appearance of their coats. Previous genetic analysis has placed wly in a 5.9 Mb interval on Chromosome 11 that contains over 200 known genes. Assignment of wly to one of these genes is needed in order to provide probes that would ultimately facilitate a complete molecular analysis of that gene’s role in the normal and disrupted development of the mammalian integument.
Here, a large intraspecific backcross family was used to genetically map wly to a smaller (0.8 Mb) span on mouse Chromosome 11 that includes fewer than 20 genes. DNA sequencing of the coding regions in two of these candidates known to be expressed in skin has revealed a 955 bp, wly-specific deletion. This deletion, which lies within the coordinates of both Slc5a10 [for solute carrier family 5 (sodium/glucose cotransporter), member 10] and Fam83g (for family with sequence similarity 83, member G), alters the splicing of mutant Fam83g transcripts only, and is predicted to result in a severely truncated (probably non-functional) protein product.
We suggest that this mutation in Fam83g is the likely basis of the mouse wooly phenotype.
The spontaneous wooly (or woolly) mutation (abbreviated wly) was initially identified at The Jackson Laboratory (Bar Harbor, ME, USA) among a litter of inbred NOD/ShiLtJ mice . Mutants are recognized as early as 3–4 weeks of age by the rough or matted appearance of their coats (see Figure 1), but—in spite of this presentation—all hair types examined (auchene, guard, zigzag and vibrissae) appear microscopically normal , and histological examination of skin has revealed no marked anomaly compared to normal-coated (heterozygous) littermates . When wly was mapped to mouse Chromosome (Chr) 11, it was immediately tested for genetic complementation in crosses with the waved 2 (wa2) mutant , but since no affected progeny were produced, these recessive variants must be due to defects in distinct, albeit syntenic, genes . Indeed, homozygosity mapping based on 54 affected F2 animals has placed wly between D11Mit313 and D11Mit261, a 5.9 Mb interval that contains over 200 known genes , but does not include wa2 or any other obvious gene candidates.
To associate this mutation with a causative molecular defect, wly was fine-mapped to a genetic region where fewer than 20 genes are located—only a few of which are known to be expressed in skin. Direct sequence analysis of the coding regions in two of these candidates has identified a mutant-specific defect in Fam83g (for family with sequence similarity 83, member G) that we propose to be the likely genetic basis of the wly mutation.
Standard inbred strains A/J, C57BL/6J, and NOD/ShiLtJ; and NOD/ShiLtJ-wly/J mutant mice were obtained from The Jackson Laboratory (Bar Harbor, ME, USA). Mice homozygous for the mutant wly allele were reliably identified by the matted appearance of their coats which is first evident by 3–4 weeks of age and persists throughout life. Both male and female wly/wly homozygotes appear to be fully fertile, and we have maintained the NOD/ShiLtJ-wly/J line since 2009 by crossing homozygotes. The treatment and use of all mice in this study were compliant with protocols approved by the Institutional Animal Care and Use Committee (IACUC) at Central Connecticut State University (New Britain, CT, USA).
Genomic DNA was isolated from 2–4 mm tail-tip biopsies taken from two-week-old mice, using Nucleospin® kits from BD Biosciences (Palo Alto, CA, USA), as directed. DNA samples from standard inbred and mutant strains that we do not routinely maintain in our colony were purchased from The Jackson Laboratory’s Mouse DNA Resource.
The polymerase chain reaction (PCR) was performed in 13 ul reactions using the Titanium® PCR kit from BD Biosciences, as directed. Oligonucleotide primers for PCR were designed and synthesized by Invitrogen (Carlsbad, CA, USA), based on sequence information available online [3, 4]. In addition to standard microsatellite markers  on Chr 11, six DNA markers based on single-nucleotide-polymorphisms previously reported to differ between strains A/J and NOD/ShiLtJ [3, 4] were also scored. These markers (herein designated SNP1-6) are described in detail in Additional file 1 & Additional file 2. To visualize PCR product sizes, reactions plus 2 ul loading buffer (bromophenol blue in 20% Tris-buffered sucrose) were electrophoresed through 3.5% NuSieve® agarose (Lonza, Rockland, ME, USA) gels. Gels were stained with ethidium bromide (0.5ug/mL) and photographed under ultraviolet light. For sequence analysis, about 1.5 ug of individual PCR amplimers were purified and concentrated into a 30 ul volume using QIAquick® PCR Purification kits (Qiagen, Valencia, CA, USA). Amplimers were shipped to SeqWright, Inc. (Houston, TX, USA) for primer-extension sequencing.
Total RNA was isolated from tail skin samples taken from 3-month-old mice using the Nucleospin® RNA L kit by Macherey-Nagel (Easton, PA, USA). cDNA was generated using the SMARTer™ RACE cDNA amplification kit (Clontech Laboratories). To amplify Slc5a10-specific cDNA, primers 1F (5’ TGTTCCGGGACCCTTCCACAGGAGACCT 3’), taken from Exon 9, and 1R (5’ ATGACCAGCCGTCCCACCAGCAGCAACT 3’), taken from Exons12 and 13, were used in a “step-down” PCR reaction. The products of this initial reaction were diluted 1:10 in Tricine-KOH buffer (10 mM, pH 8.5) plus 1 mM EDTA, and were amplified again in a standard PCR reaction using a nested primer pair: 2F (5’ AGCGGTCCCTGTCTGCCCGGAACTTGAA 3’), taken from Exon 10, and 2R (5’ TGGGCATCAGCTCCATGACCAGCTTCGGGT 3’), taken from Exon 11. To amplify Fam83g-specific cDNA, primers 3F (5’ ACGGGCAGCCGCACATCAAGGAAGTGGT 3’), taken from Exon 1, and 3R (5’ AGCACAATGGGCTCTGGCTCTGGCTCCT 3’), taken from Exon 4, were used in a standard PCR reaction. The products of this initial reaction were diluted (as above), and were amplified again with a nested primer pair: 4F (5’ TGCGCAAGATGGTCAGCCAGGCGCAGAA 3’), taken from Exon 1, and 4R (5’ ATGGGCTCTGGCTCTGGCTCCTTCTCCA 3’), taken from Exon 4. Final (second-round) amplimers were purified (as described above) and shipped to SeqWright, Inc., for primer-extension sequencing.
Meiotic fine-mapping of wly
To more precisely locate wly on proximal Chr 11, (A/J x NOD/ShiLtJ-wly/J)F1, wly/+ females were bred back to NOD/ShiLtJ-wly/wly males, producing a large family of 1,679 backcross progeny that segregated for wly (and numerous other molecular markers). These progeny were typed for wly and six, PCR-scorable, microsatellite markers on Chr 11, as summarized in Figure 2. These data are in agreement with the 1 wild type : 1 mutant ratio expected for a testcross (χ2 < 0.05; P > 0.8), suggesting that wly/wly mice are fully viable, at least compared to heterozygous wly/+ littermates. Furthermore, this analysis indicated that wly is located between D11Mit208 and D11Mit242, a region of about 4.8 Mb and fully consistent with wly’s previously-defined mapping interval between D11Mit313 and D11Mit216.
To further restrict the physical position of wly on Chr 11, mice with a meiotic crossover between D11Mit208 and D11Mit242 were typed for six, single-nucleotide polymorphisms (herein designated SNP1-6) known to lie in this “critical region”. This analysis identified two recombinants that carried a crossover between SNP2 and wly, and one with a crossover between wly and SNP6. (No crossovers separated wly from D11Mit260, or from SNP3, 4, or 5.) Thus, wly must be located between SNP2 and SNP6, an interval that measures less than 0.8 Mb and includes fewer than 20 genes or predicted genes (see Figure 3a).
Sequence analysis of gene candidates from the wly-critical interval
For some of the genes located in the SNP2 to SNP6 interval (including Slc47a1, Epn2, Grap, Ulk2, Mapk7, and Akap10) loss-of-function alleles have been engineered (by others), and null-allele homozygotes have been reported to display no apparent changes to hair development or coat texture, making them unlikely candidates for the gene responsible for the wly mutation. For several others (including four predicted genes not shown in Figure 3a), available expression data [3, 4] failed to suggest any obvious functional role in skin. By contrast, two other genes in the critical interval, the overlapping Slc5a10 [for solute carrier family 5 (sodium/glucose cotransporter), member 10] and Fam83g (for family with sequence similarity 83, member G) genes, were isolated from cDNA libraries derived from mouse skin . To investigate these “primary” candidates as the possible basis of the mutant wooly phenotype, the coding regions of both genes were sequenced in genomic DNA from A/J, C57BL/6J, NOD/ShiLtJ, and mutant NOD/ShiLtJ-wly/J mice.
While this analysis revealed several sites that are polymorphic among these four strains—some new and some previously reported (as summarized in Additional file 3 and shown in detail in Additional file 4 & Additional file 5)—the only sequence distinction found between the coisogenic NOD/ShiLtJ and NOD/ShiLtJ-wly/J strains was a 955 bp deletion in the mutant strain (see Additional file 6a) that lies in Intron 10–11 of Slc5a10 and removes part of Intron 2–3 and Exon 3 of Fam83g (see Figure 3b). Because none of 17 other inbred strains tested (including nine strains of Swiss origin, and therefore related to or derived from the NOD strain) show this sequence alteration, while wly mutants tested from three distinct colonies do (Additional file 6b), we suggest that this 955 bp deletion is specifically associated with the wly mutation.
A Fam83g deletion is likely to be the molecular basis of the mutant wooly phenotype
Because this deletion lies in Intron 10–11 of the Slc5a10 gene, we amplified sequences between Exons 9 and 13 from cDNA templates based on total RNA isolated from wly/wly or wild type skin to determine if this mutation could affect splicing of Slc5a10 transcripts. Amplimers copied from wly/wly and wild type NOD/ShiLtJ cDNA templates were identical in length, and sequencing verified that Exons 10 and 11 of Slc5a10 are spliced normally in wly/wly mutants (data not shown).
Because this deletion removes the splice acceptor site at the 5’ end of Exon 3 in Fam83g, it was predicted that the mutant transcript would likely be spliced such that Exon 2 is joined with Exon 4, rather than Exon 3 (see Figure 3b). This prediction was tested by PCR amplification of sequences between Exon 1 and Exon 4 from cDNA templates based on total RNA isolated from mutant or wild type skin. Sequencing of the 359 bp, wly-specific product and the 484 bp, wild type product (see Additional file 6c) demonstrated that the predicted, aberrant Exon 2–4 splice does occur in wly/wly mutants; while the expected Exon 1-2-3-4 joining occurs in wild type mice. Because the skipped Exon 3 contains 125 nucleotides, the inappropriate junction of Exon 2 with 4 causes a frameshift after Codon 224 that is predicted to introduce 24 novel amino acids in the mutant Fam83g gene product before terminating translation quite early in Exon 4 (see Figure 4). Such a severe truncation in the Fam83g protein sequence (from an expected 812 a.a. for the normal protein to only 524 a.a. in the mutant) is very likely to negatively impact protein function, and suggests further that this mutation could be the basis of the mutant wooly phenotype.
The eight members of the Fam83 family are mostly uncharacterized proteins, in both mouse and man. Fam83a, also known as tumor antigen BJ-TSA-9 , is a novel, tumor-specific protein highly expressed in human lung adenocarcinoma cells. Fam83d, also referred to as spindle protein CHICA , is a cell-cycle-regulated spindle component which localizes to the mitotic spindle and is both up regulated and phosphorylated during mitosis. Defects in the gene encoding Fam83h cause autosomal dominant hypocalicified ameliogenesis imperfecta (ADHCAI) [15–19]. Fam83b, c, f, and g are uncharacterized proteins present across vertebrates, while Fam83e is an uncharacterized protein found only in mammals. Some proteins with structural similarity to Fam83g’s N-terminal domain of unknown function (Pfam domain DUF1669, [20, 21]; superfamily domain SSF5624, [22, 23]; conserved domain cd 09119, [24, 25]; see Figure 4) are known to be phospholipases, but this domain in Fam83g shows only trace similarity to the phospholipase D catalytic domain and lacks the functionally-important histidine residue [26, 27], so while Fam83g may share a similar 3-dimensional fold with some phospholipase D-like enzymes, it is unlikely to display phospholipase D-like activity. An additional proline-rich domain has been identified from amino acid 470 to 511 (Prosite profiles, PS5099, [28, 29]; see Figure 4), but again no function has yet been assigned. Whatever their functional significance, both of these conserved protein domains would be partially or entirely removed by the Fam83g mutation in wooly mice (see Figure 4). We predict that this mutant Fam83g protein is, therefore, likely to be non-functional, consistent with our suggestion that the lack of normal Fam83g product in wly/wly homozygotes may be responsible for the mutant wooly phenotype.
More direct evidence for a causal link between this deletion and the wooly phenotype would require, for example, transgenic rescue of mutant homozygotes, or recreation of the wooly phenotype in engineered Fam83g “knock-out” homozygotes. While such single-addition and single-subtraction experiments are beyond our laboratory’s ability to perform, we anticipate testing for complementation between NOD/ShiLtJ- wly/J and a recessive, Fam83g knock-out variant (as soon as one becomes available). Non-complementation (i.e., production of phenotypically mutant wly/Fam83gk.o. offspring) would provide definitive proof that homozygousity for Fam83g defects, alone, is the molecular basis of the mutant wooly phenotype.
No other similar mutations affecting skin or hair have been described in this region on mouse Chr 11 , nor are we aware of any human conditions involving skin or hair in the orthologous region of human Chr 17p11-12 . In man, the term “woolly hair” (WH) is used to describe a group of inherited hair shaft disorders characterized by fine and tightly curled hair , but—based on their distinct phenotypes, chromosomal locations, and, in some cases, known molecular bases [see [32, 33]—none of these described disorders appears related to mouse wooly. We therefore suggest that the NOD/ShiLtJ- wly/J mouse strain may provide a unique animal resource, the study of which will be crucial to any future investigation of Fam83g and its functional role in the normal or disrupted development of the mammalian integument. For example, it will be interesting to learn where Fam83g is expressed in normal and mutant hair follicles (e.g., in the dermal papilla, the epithelial lineage, or other hair-follicle-associated structures), and whether its expression varies during the anagen, catagen and telogen phases of the hair follicle growth cycle. While one previous histological evaluation failed to uncover a microscopic manifestation of the wooly phenotype , we hope—especially with molecular probes made possible by the likely genetic assignment of wly to Fam83g—that in situ-based expression or immunohistological analyses can now be approached to finally reveal the cellular basis of the wooly phenotype.
The 955 bp deletion we describe appears to be specifically associated with the wly mutation; it alters the splicing of mutant Fam83g transcripts; and is predicted to generate a severely truncated, mutant Fam83g protein. We therefore suggest that this defect is likely to be the molecular basis of the mutant wooly phenotype.
TRK is a professor in the Biomolecular Sciences Department at Central Connecticut State University (New Britain, CT). LAR, KMC, EBA, DVS, and AMF were undergraduates majoring in Biomolecular Sciences or Biochemistry when they conducted this research.
Harris B, Ward-Bailey PF, Bronson RT: Wooly: a new hair mutation on mouse chromosome 11. Mouse mutant resources Web site. 2003, Bar Harbor, Maine: The Jackson Laboratory, MGI Direct Data Submission http://mousemutant.jax.org/articles/MMRmutantwooly.html
Luetteke NC, Phillips HK, Qiu TH, Copeland NG, Earp HS, Jenkins NA, Lee DC: The mouse waved-2 phenotype results from a point mutation in the EGF receptor tyrosine kinase. Genes Dev. 1994, 8: 399-413. 10.1101/gad.8.4.399.
Ensembl M: Mouse genome sequencing consortium: the European Bioinformatics Institute (EBI) and the Welcome Trust Sanger Institute (WTSI). 2012, Release 65, Available at http://www.ensembl.org/Mus_musculus/
Mouse Genome Database (MGD): Mouse genome database group: the mouse genome informatics website. 2012, Bar Harbor, ME: The Jackson Laboratory, Available at http://www.informatics.jax.org
Dietrich WF, Miller J, Steen R, Merchant MA, Damron-Boles D, Husain Z, Dredge R, Daly MJ, Ingalls KA, O'Connor TJ: A comprehensive genetic map of the mouse genome. Nature. 1996, 380: 149-152. 10.1038/380149a0.
Nakamura T, Yonezawa A, Hashimoto S, Katsura T, Inui K: Disruption of multidrug and toxin extrusion MATE1 potentiates cisplatin-induced nephrotoxicity. Biochem Pharmacology. 2010, 80: 1762-1767. 10.1016/j.bcp.2010.08.019.
Chen H, Ko G, Zatti A, Di Giacomo G, Lui L, Raiteri E, Perucco E, Collesi C, Min W, Zeiss C, De Camilli P, Cremona O: Embryonic arrest at midgestation and disruption of notch signaling produced by absence of both epsin 1 and epsin 2 in mice. Proc Natl Acad Sci USA. 2009, 106: 13838-13843. 10.1073/pnas.0907008106.
Shen R, Ouyang YB, Qu CK, Alonso A, Sperzel L, Mustelin T, Kaplan MH, Feng GS: Grap negatively regulates T-cell receptor-elicited lymphocyte proliferation and interleukin-2 induction. Mol Cell Biol. 2002, 22: 3230-3236. 10.1128/MCB.22.10.3230-3236.2002.
Cheong H, Lindsten T, Wu J, Lu C, Thompson CB: Ammonia-induced autophagy is independent of ULK1/ULK2 kinases. Proc Natl Acad Sci USA. 2011, 108: 11121-11126. 10.1073/pnas.1107969108.
Hayashi M, Kim SW, Imanaka-Yosjida K, Yoshida T, Abel ED, Eliceiri B, Yang Y, Ulevitch RJ, Lee JD: Targeted deletion of BMK1/ERK5 in adult mice perturbs vascular integrity and leads to endothelial failure. J Clin Invest. 2004, 113: 1138-1148.
Tingley WG, Pawlikowska L, Zaroff JG, Kim T, Nguyen T, Young SG, Vranizan K, Kwok P-Y Whooley MA, Conklin BR: Gene-trapped mouse embryonic stem cell-derived cardiac myocytes and human genetics implicate AKAP10 in heart rhythm regulation. Proc Natl Acad Sci USA. 2007, 104: 8461-8466. 10.1073/pnas.0610393104.
Kawai J, Shinagawa A, Shibata K, Yoshino M, Itoh M, Ishii Y, Arakawa T, Hara A, Fukunishi Y, Konno H, Adachi J, Fukuda S, Aizawa K, Izawa M, Nishi K, Kiyosawa H, Kondo S, Yamanaka I, Saito T, Okazaki Y, Gojobori T, Bono H, Kasukawa T, Saito R, Kadota K, Matsuda H, Ashburner M, Batalov S, Casavant T, Fleischmann W: Functional annotation of a full-length mouse cDNA collection. Nature. 2001, 409: 685-690. 10.1038/35055500.
Li Y, Dong X, Yin Y, Su Y, Xu Q, Zhang Y, Pang X, Zhang Y, Chen W: BJ-TSA-9, a novel human tumor-specific gene, has potential as a biomarker of lung cancer. Neoplasia. 2005, 7: 1073-1080. 10.1593/neo.05406.
Santamaria A, Nagel S, Sillje HH, Nigg EA: The spindle protein CHICA mediates localization of the chromokinesin Kid to the mitotic spindle. Curr Biol. 2008, 18: 723-729. 10.1016/j.cub.2008.04.041.
Hyun HK, Lee SK, Lee KE, Kamg HY, Kim EJ, Choung PH, Kim JW: Identification of a novel FAM83H mutation and microhardness of an affected molar in autosomal dominant hypocalcified amelogenesis imperfecta. Int Endod J. 2009, 42: 1039-1043. 10.1111/j.1365-2591.2009.01617.x.
Ding Y, Estrella MR, Hu YY, Chan HL, Zhang HD, Kin JW, Simmer JP, Hu JC: Fam83h Is associated with intracellular vesicles and ADHCAI. J Dent Res. 2009, 88: 991-996. 10.1177/0022034509349454.
El-Sayed W, Shore RC, Parry DA, Inglehearn CF, Mighell AJ: Ultrastructural analyses of deciduous teeth affected by hypocalcified amelogenesis imperfecta from a family with a novel Y458X FAM83H nonsense mutation. Cell Tissues Organs. 2010, 191: 235-239. 10.1159/000252801.
Wright JT, Frazier-Bowers S, Simmons D, Alexander K, Crawford P, Han ST, Hart PS, Hart TC: Phenotypic variation in FAM83H-associated amelogenesis imperfecta. J Dent Res. 2009, 88: 356-360. 10.1177/0022034509333822.
Hart PS, Becertik S, Cogulu D, Emingil G, Ozdemir-Ozenen D, Han ST, Sulima PP, Firatli E, Hart TC: Novel FAM83H mutations in Turkish families with autosomal dominant hypocalcified amelogenesis imperfecta. Clin Genet. 2009, 75: 401-404.
Pfam: Wellcome trust Sanger institute. Release 26.0. 2012, Available at http://pfam.sanger.ac.uk/
Bateman A, Coggill P, Finn RD: DUFs: families in search of a function. Acta Crystallogr Sect F Struct Biol Cryst Commun. 2010, 66: 1148-1152. 10.1107/S1744309110001685.
Structural Classification of Proteins Database (SCOP): MRC laboratory of molecular biology and centre for protein engineering. Release 1.75. 2012, Available at http://scop.mrc-lmb.cam.ac.uk/scop/
Andreva A, Howorth D, Chandonia JM, Brenner SE, Hubbard TJP, Chothia C, Murzin AG: Data growth and its impact on the SCOP database: New developments. Nucleic Acids Res. 2008, 36: D419-D425.
Conserved Domain Database (CDD): National center for biotechnology information. 2012, Available at http://www.ncbi.nlm.nih.gov/cdd/
Marchler-Bauer A, Lu S, Anderson JB, Chitsaz F, Derbyshire MK, DeWeese-Scott C, Fong JH, Geer LY, Geer RC, Gonzales NR, Gwadz M, Hurwitz DI, Jackson JD, Ke Z, Lanczycki CJ, Lu F, Marchler GH, Mullokandov M, Omelchenko MV, Robertson CL, Song JS, Thanki N, Yamashita RA, Zhang D, Zhang N, Zheng C, Bryant SH: CDD: a conserved domain database for the functional annotation of proteins. Nucleic Acids Res. 2011, 39: 225-229. 10.1093/nar/gkq769.
Raikwar NS, Bowen RF, Deeg MA: Mutating His29, His125, His133 or His158 abolishes glycosylphosphatidylinositol-specific Phospholipase D catalytic activity. Biochem J. 2005, 391: 285-289. 10.1042/BJ20050656.
Stuckey JA, Dixon JE: Crystal structure of a Phospholipase D family member. Nat Struct Biol. 1999, 6: 278-284. 10.1038/6716.
PROSITE: Swiss institute bioinformtics. Release 26.79. 2012, Available at http://www.expasy.org/prosite/PROSITE
Hulo N, Bairoch A, Bulliard V, Cerutti L, Cuche B, De Castro E, Lachaize C, Langendijk-Genevaux PS, Sigrist CJA: The 20 years of PROSITE. Nucleic Acids Res. 2007, 36: D245-D249. 10.1093/nar/gkm977.
Online Medelian Inheritance in Man (OMIM): Johns Hopkins university (Baltimore, MD) and national center for biotechnology information. 2012, Bethesda, MD: National Library of Medicine, Available at http://www.ncbi.nlm.nih.gov/omim/
Chien AJ, Valentine MC, Sybert VP: Hereditary woolly hair and keratosis pilaris. J Am Acad Dermatol. 2006, 54: S35-S39. 10.1016/j.jaad.2005.01.092.
Shimomura Y, Wajid M, Petukhova L, Shapiro L, Christiano AM: Mutations in the lipase H gene underlie autosomal recessive woolly hair/hypotrichosis. J Invest Dermatol. 2009, 129: 622-628. 10.1038/jid.2008.290.
Shimomura Y, Wajid M, Ishii Y, Shapiro L, Petukhova L, Gordon D, Christiano AM: Disruption of P2RY5, and orphan G protein-coupled receptor, underlies autosomal recessive woolly hair. Nat Genet. 2008, 40: 335-339. 10.1038/ng.100.
The authors thank CCSU undergraduates Nisrine Dagamseh, Michelle Kayser, Andrea Nixon, Andy Morrison, Lori LeBlanc, Randy Taylor, Amarilis Perez, Jose Horak; and all the high-school interns (Edie Tinker, Mariam Hasan, Somaly Chhean, Rick Deschenes, Alicia Davis, Joshua Wrice, Marcelino Thillet, Ananda Thillet, Danerick Peralta, Diego Peralta, Ashley Feliciano, Krystal Garcia, Kayla Garcia, Jodalis Montalvo, Milagros Molina, Brianna Cirinna, and Rebecca Fuentes) for help with DNA isolation and marker typing. We also thank Mary Mantzaris for excellent animal care. This work was supported by NIH AREA Grant 1R15AR059572.
The authors declare that they have no competing interests.
LAR and KMC led all aspects of this study, including experimental design, data acquisition and interpretation. DVS and EBA made substantial contributions to the genetic analysis. EBA and AMF participated in the cDNA analysis. TRK conceived of the study, carried out all procedures involving mice, and drafted the manuscript. All authors read, edited, and approved the final manuscript.
Legairre A Radden, Kevin M Child contributed equally to this work.