title: Mapping whole genome shotgun sequence and variant calling in mammalian species without their reference genomes.
abstract: Genomics research in mammals has produced reference genome sequences that are essential for identifying variation associated with disease.  High quality reference genome sequences are now available for humans, model species, and economically important agricultural animals.  Comparisons between these species have provided unique insights into mammalian gene function.  However, the number of species with reference genomes is small compared to those needed for studying molecular evolutionary relationships in the tree of life.  For example, among the even-toed ungulates there are approximately 300 species whose phylogenetic relationships have been calculated in the 10k trees project.  Only six of these have reference genomes:  cattle, swine, sheep, goat, water buffalo, and bison.  Although reference sequences will eventually be developed for additional hoof stock, the resources in terms of time, money, infrastructure and expertise required to develop a quality reference genome may be unattainable for most species for at least another decade.  In this work we mapped 35 Gb of next generation sequence data of a Katahdin sheep to its own species' reference genome ( Ovis aries Oar3.1) and to that of a species that diverged 15 to 30 million years ago ( Bos taurus UMD3.1).  In total, 56% of reads covered 76% of UMD3.1 to an average depth of 6.8 reads per site, 83 million variants were identified, of which 78 million were homozygous and likely represent interspecies nucleotide differences. Excluding repeat regions and sex chromosomes, nearly 3.7 million heterozygous sites were identified in this animal vs. bovine UMD3.1, representing polymorphisms occurring in sheep.  Of these, 41% could be readily mapped to orthologous positions in ovine Oar3.1 with 80% corroborated as heterozygous.  These variant sites, identified via interspecies mapping could be used for comparative genomics, disease association studies, and ultimately to understand mammalian gene function.
authors.affiliation: Department of Biochemistry and Molecular Biology, School of Medicine, University of Louisville, Louisville, KY, 40202, USA ; Intrepid Bioinformatics, Louisville, KY, 40202, USA.
authors.forename: Ted
authors.initials: T
authors.lastname: Kalbfleisch
authors.affiliation: USDA Meat Animal Research Center, Clay Center, Nebraska, 68933, USA.
authors.forename: Michael P
authors.initials: MP
authors.lastname: Heaton
references.citation: Genome Res. 2013 Jul;23(7):1063-8
references.pmid: 23817047
references.citation: Nature. 2007 Jun 14;447(7146):799-816
references.pmid: 17571346
references.citation: Bioinformatics. 2010 Mar 1;26(5):589-95
references.pmid: 20080505
references.citation: Genome Res. 2010 Sep;20(9):1297-303
references.pmid: 20644199
references.citation: BMC Vet Res. 2010 Apr 30;6:23
references.pmid: 20433741
references.citation: Anim Genet. 2010 Oct;41(5):449-53
references.pmid: 20809919
references.citation: Genome Biol. 2009;10(4):R42
references.pmid: 19393038
references.citation: Nature. 2012 Sep 6;489(7414):57-74
references.pmid: 22955616
references.citation: PLoS Biol. 2012 Feb;10(2):e1001258
references.pmid: 22346734
references.citation: Bioinformatics. 2011 Aug 1;27(15):2156-8
references.pmid: 21653522
references.citation: J Anim Sci. 1991 Mar;69(3):993-9
references.pmid: 2061269
references.citation: Genome Res. 2002 Apr;12(4):656-64
references.pmid: 11932250
references.citation: Nat Biotechnol. 2011 Jan;29(1):24-6
references.pmid: 21221095
references.citation: Am J Hum Genet. 2007 Sep;81(3):559-75
references.pmid: 17701901
references.citation: Bioinformatics. 2009 Aug 15;25(16):2078-9
references.pmid: 19505943
references.citation: Curr Protoc Bioinformatics. 2009 Dec;Chapter 1:Unit1.4
references.pmid: 19957273
references.citation: BMC Vet Res. 2008 Jul 14;4:25
references.pmid: 18625065
references.citation: Nature. 2012 Nov 1;491(7422):56-65
references.pmid: 23128226