An Annotated Draft Genome of the Mountain Hare (Lepus timidus).

João P Marques, Fernando A Seixas, Liliana Farelo, Colin M Callahan, Jeffrey M Good, W Ian Montgomery, Neil Reid, Paulo C Alves, Pierre Boursot, José Melo-Ferreira
Author Information
  1. João P Marques: CIBIO, Centro de Investigação em Biodiversidade e Recursos Genéticos, InBIO Laboratório Associado, Universidade do Porto, Vairão, Portugal.
  2. Fernando A Seixas: CIBIO, Centro de Investigação em Biodiversidade e Recursos Genéticos, InBIO Laboratório Associado, Universidade do Porto, Vairão, Portugal.
  3. Liliana Farelo: CIBIO, Centro de Investigação em Biodiversidade e Recursos Genéticos, InBIO Laboratório Associado, Universidade do Porto, Vairão, Portugal.
  4. Colin M Callahan: Division of Biological Sciences, University of Montana.
  5. Jeffrey M Good: Division of Biological Sciences, University of Montana.
  6. W Ian Montgomery: Institute of Global Food Security, School of Biological Sciences, Queen's University Belfast, United Kingdom.
  7. Neil Reid: Institute of Global Food Security, School of Biological Sciences, Queen's University Belfast, United Kingdom.
  8. Paulo C Alves: CIBIO, Centro de Investigação em Biodiversidade e Recursos Genéticos, InBIO Laboratório Associado, Universidade do Porto, Vairão, Portugal.
  9. Pierre Boursot: Institut des Sciences de l'Evolution Montpellier (ISEM), Univ Montpellier, CNRS, IRD, EPHE, France.
  10. José Melo-Ferreira: CIBIO, Centro de Investigação em Biodiversidade e Recursos Genéticos, InBIO Laboratório Associado, Universidade do Porto, Vairão, Portugal.

Abstract

Hares (genus Lepus) provide clear examples of repeated and often massive introgressive hybridization and striking local adaptations. Genomic studies on this group have so far relied on comparisons to the European rabbit (Oryctolagus cuniculus) reference genome. Here, we report the first de novo draft reference genome for a hare species, the mountain hare (Lepus timidus), and evaluate the efficacy of whole-genome re-sequencing analyses using the new reference versus using the rabbit reference genome. The genome was assembled using the ALLPATHS-LG protocol with a combination of overlapping pair and mate-pair Illumina sequencing (77x coverage). The assembly contained 32,294 scaffolds with a total length of 2.7 Gb and a scaffold N50 of 3.4 Mb. Re-scaffolding based on the rabbit reference reduced the total number of scaffolds to 4,205 with a scaffold N50 of 194 Mb. A correspondence was found between 22 of these hare scaffolds and the rabbit chromosomes, based on gene content and direct alignment. We annotated 24,578 protein coding genes by combining ab-initio predictions, homology search, and transcriptome data, of which 683 were solely derived from hare-specific transcriptome data. The hare reference genome is therefore a new resource to discover and investigate hare-specific variation. Similar estimates of heterozygosity and inferred demographic history profiles were obtained when mapping hare whole-genome re-sequencing data to the new hare draft genome or to alternative references based on the rabbit genome. Our results validate previous reference-based strategies and suggest that the chromosome-scale hare draft genome should enable chromosome-wide analyses and genome scans on hares.

Keywords

Associated Data

Dryad | 10.5061/dryad.x95x69pd8

References

  1. Heredity (Edinb). 2006 Nov;97(5):355-65 [PMID: 16896342]
  2. Gigascience. 2017 Dec 1;6(12):1-5 [PMID: 29099922]
  3. Bioinformatics. 2011 Feb 15;27(4):578-9 [PMID: 21149342]
  4. J Hered. 2016 Jul;107(4):295-308 [PMID: 26921276]
  5. Bioinformatics. 2014 May 1;30(9):1236-40 [PMID: 24451626]
  6. Genome Biol. 2018 Jul 30;19(1):91 [PMID: 30056805]
  7. Proc Natl Acad Sci U S A. 2019 Nov 26;116(48):24150-24156 [PMID: 31712446]
  8. Proc Natl Acad Sci U S A. 2011 Jan 25;108(4):1513-8 [PMID: 21187386]
  9. Genome Biol. 2013 May 27;14(5):R47 [PMID: 23710727]
  10. Sci Data. 2019 May 22;6(1):68 [PMID: 31118413]
  11. Nucleic Acids Res. 2017 Jan 4;45(D1):D183-D189 [PMID: 27899595]
  12. Bioinformatics. 2014 Aug 1;30(15):2114-20 [PMID: 24695404]
  13. BMC Bioinformatics. 2014 Jun 12;15:182 [PMID: 24925680]
  14. Bioinformatics. 2015 Oct 1;31(19):3210-2 [PMID: 26059717]
  15. G3 (Bethesda). 2019 Jun 5;9(6):1785-1793 [PMID: 31000506]
  16. Genome Res. 2008 Jan;18(1):188-96 [PMID: 18025269]
  17. Mol Ecol. 2007 Feb;16(3):605-18 [PMID: 17257116]
  18. Cell Rep. 2015 Jan 6;10(1):112-22 [PMID: 25565328]
  19. Mol Ecol. 2009 Jun;18(12):2643-58 [PMID: 19457181]
  20. Science. 2018 Jun 22;360(6395):1355-1358 [PMID: 29930138]
  21. Cytogenet Genome Res. 2002;96(1-4):223-7 [PMID: 12438803]
  22. Bioinformatics. 2011 Nov 1;27(21):2987-93 [PMID: 21903627]
  23. Nucleic Acids Res. 2016 Jan 4;44(D1):D457-62 [PMID: 26476454]
  24. PeerJ. 2018 Jun 4;6:e4958 [PMID: 29888139]
  25. Nat Methods. 2012 Mar 04;9(4):357-9 [PMID: 22388286]
  26. Bioinformatics. 2015 Aug 1;31(15):2443-51 [PMID: 25810435]
  27. Science. 2014 Aug 29;345(6200):1074-1079 [PMID: 25170157]
  28. Nucleic Acids Res. 2016 Jan 4;44(D1):D279-85 [PMID: 26673716]
  29. Syst Biol. 2004 Jun;53(3):433-47 [PMID: 15503672]
  30. Nucleic Acids Res. 2018 Jan 4;46(D1):D649-D655 [PMID: 29145629]
  31. Bioinformatics. 2011 Mar 15;27(6):764-70 [PMID: 21217122]
  32. Curr Protoc Bioinformatics. 2013;43:11.10.1-11.10.33 [PMID: 25431634]
  33. Genome Res. 2012 Mar;22(3):557-67 [PMID: 22147368]
  34. Bioinformatics. 2011 Mar 15;27(6):863-4 [PMID: 21278185]
  35. Bioinformatics. 2009 Aug 15;25(16):2078-9 [PMID: 19505943]
  36. Mol Ecol. 2005 Sep;14(10):3005-16 [PMID: 16101770]
  37. Mol Ecol. 1997 May;6(5):463-7 [PMID: 9161014]
  38. Bioinformatics. 2016 Apr 1;32(7):1088-90 [PMID: 26614127]
  39. Philos Trans R Soc Lond B Biol Sci. 2008 Sep 12;363(1505):2831-9 [PMID: 18508749]
  40. Nature. 2011 Jul 13;475(7357):493-6 [PMID: 21753753]
  41. Mol Phylogenet Evol. 2003 Apr;27(1):70-80 [PMID: 12679072]
  42. BMC Bioinformatics. 2009 Dec 15;10:421 [PMID: 20003500]

MeSH Term

Animals
Female
Genome
Genomics
Hares
Molecular Sequence Annotation
Transcriptome

Word Cloud

Created with Highcharts 10.0.0genomeharereferencerabbitLepusdraftwhole-genomeusingnewscaffoldsbaseddatadenovotimidusre-sequencinganalysessequencingassemblytotalscaffoldN50transcriptomehare-specificharesHaresgenusprovideclearexamplesrepeatedoftenmassiveintrogressivehybridizationstrikinglocaladaptationsGenomicstudiesgroupfarreliedcomparisonsEuropeanOryctolaguscuniculusreportfirstspeciesmountainevaluateefficacyversusassembledALLPATHS-LGprotocolcombinationoverlappingpairmate-pairIllumina77xcoveragecontained32294length27 Gb34 MbRe-scaffoldingreducednumber4205194 Mbcorrespondencefound22chromosomesgenecontentdirectalignmentannotated24578proteincodinggenescombiningab-initiopredictionshomologysearch683solelyderivedthereforeresourcediscoverinvestigatevariationSimilarestimatesheterozygosityinferreddemographichistoryprofilesobtainedmappingalternativereferencesresultsvalidatepreviousreference-basedstrategiessuggestchromosome-scaleenablechromosome-widescansAnnotatedDraftGenomeMountainHareLagomorphaLeporidsannotation

Similar Articles

Cited By