Comparison of the Equine Reference Sequence with Its Sanger Source Data and New Illumina Reads.
The reference assembly for the domestic horse, EquCab2, published in 2009, was built using approximately 30 million Sanger reads from a Thoroughbred mare named Twilight. Contiguity in the assembly was facilitated using nearly 315 thousand BAC end sequences from Twilight's half brother Bravo. Si...
Main Authors: | , , , , , , |
---|---|
Format: | Article |
Language: | English |
Published: |
Public Library of Science (PLoS)
2015-01-01
|
Series: | PLoS ONE |
Online Access: | http://europepmc.org/articles/PMC4479572?pdf=render |
id |
doaj-90e2856a8abb4f91887f16007f30752b |
---|---|
record_format |
Article |
spelling |
doaj-90e2856a8abb4f91887f16007f30752b2020-11-24T20:51:47ZengPublic Library of Science (PLoS)PLoS ONE1932-62032015-01-01106e012685210.1371/journal.pone.0126852Comparison of the Equine Reference Sequence with Its Sanger Source Data and New Illumina Reads.Jovan Rebolledo-MendezMatthew S HestandStephen J ColemanZheng ZengLudovic OrlandoJames N MacLeodTed KalbfleischThe reference assembly for the domestic horse, EquCab2, published in 2009, was built using approximately 30 million Sanger reads from a Thoroughbred mare named Twilight. Contiguity in the assembly was facilitated using nearly 315 thousand BAC end sequences from Twilight's half brother Bravo. Since then, it has served as the foundation for many genome-wide analyses that include not only the modern horse, but ancient horses and other equid species as well. As data mapped to this reference has accumulated, consistent variation between mapped datasets and the reference, in terms of regions with no read coverage, single nucleotide variants, and small insertions/deletions have become apparent. In many cases, it is not clear whether these differences are the result of true sequence variation between the research subjects' and Twilight's genome or due to errors in the reference. EquCab2 is regarded as "The Twilight Assembly." The objective of this study was to identify inconsistencies between the EquCab2 assembly and the source Twilight Sanger data used to build it. To that end, the original Sanger and BAC end reads have been mapped back to this equine reference and assessed with the addition of approximately 40X coverage of new Illumina Paired-End sequence data. The resulting mapped datasets identify those regions with low Sanger read coverage, as well as variation in genomic content that is not consistent with either the original Twilight Sanger data or the new genomic sequence data generated from Twilight on the Illumina platform. As the haploid EquCab2 reference assembly was created using Sanger reads derived largely from a single individual, the vast majority of variation detected in a mapped dataset comprised of those same Sanger reads should be heterozygous. In contrast, homozygous variations would represent either errors in the reference or contributions from Bravo's BAC end sequences. Our analysis identifies 720,843 homozygous discrepancies between new, high throughput genomic sequence data generated for Twilight and the EquCab2 reference assembly. Most of these represent errors in the assembly, while approximately 10,000 are demonstrated to be contributions from another horse. Other results are presented that include the binary alignment map file of the mapped Sanger reads, a list of variants identified as discrepancies between the source data and resulting reference, and a BED annotation file that lists the regions of the genome whose consensus was likely derived from low coverage alignments.http://europepmc.org/articles/PMC4479572?pdf=render |
collection |
DOAJ |
language |
English |
format |
Article |
sources |
DOAJ |
author |
Jovan Rebolledo-Mendez Matthew S Hestand Stephen J Coleman Zheng Zeng Ludovic Orlando James N MacLeod Ted Kalbfleisch |
spellingShingle |
Jovan Rebolledo-Mendez Matthew S Hestand Stephen J Coleman Zheng Zeng Ludovic Orlando James N MacLeod Ted Kalbfleisch Comparison of the Equine Reference Sequence with Its Sanger Source Data and New Illumina Reads. PLoS ONE |
author_facet |
Jovan Rebolledo-Mendez Matthew S Hestand Stephen J Coleman Zheng Zeng Ludovic Orlando James N MacLeod Ted Kalbfleisch |
author_sort |
Jovan Rebolledo-Mendez |
title |
Comparison of the Equine Reference Sequence with Its Sanger Source Data and New Illumina Reads. |
title_short |
Comparison of the Equine Reference Sequence with Its Sanger Source Data and New Illumina Reads. |
title_full |
Comparison of the Equine Reference Sequence with Its Sanger Source Data and New Illumina Reads. |
title_fullStr |
Comparison of the Equine Reference Sequence with Its Sanger Source Data and New Illumina Reads. |
title_full_unstemmed |
Comparison of the Equine Reference Sequence with Its Sanger Source Data and New Illumina Reads. |
title_sort |
comparison of the equine reference sequence with its sanger source data and new illumina reads. |
publisher |
Public Library of Science (PLoS) |
series |
PLoS ONE |
issn |
1932-6203 |
publishDate |
2015-01-01 |
description |
The reference assembly for the domestic horse, EquCab2, published in 2009, was built using approximately 30 million Sanger reads from a Thoroughbred mare named Twilight. Contiguity in the assembly was facilitated using nearly 315 thousand BAC end sequences from Twilight's half brother Bravo. Since then, it has served as the foundation for many genome-wide analyses that include not only the modern horse, but ancient horses and other equid species as well. As data mapped to this reference has accumulated, consistent variation between mapped datasets and the reference, in terms of regions with no read coverage, single nucleotide variants, and small insertions/deletions have become apparent. In many cases, it is not clear whether these differences are the result of true sequence variation between the research subjects' and Twilight's genome or due to errors in the reference. EquCab2 is regarded as "The Twilight Assembly." The objective of this study was to identify inconsistencies between the EquCab2 assembly and the source Twilight Sanger data used to build it. To that end, the original Sanger and BAC end reads have been mapped back to this equine reference and assessed with the addition of approximately 40X coverage of new Illumina Paired-End sequence data. The resulting mapped datasets identify those regions with low Sanger read coverage, as well as variation in genomic content that is not consistent with either the original Twilight Sanger data or the new genomic sequence data generated from Twilight on the Illumina platform. As the haploid EquCab2 reference assembly was created using Sanger reads derived largely from a single individual, the vast majority of variation detected in a mapped dataset comprised of those same Sanger reads should be heterozygous. In contrast, homozygous variations would represent either errors in the reference or contributions from Bravo's BAC end sequences. Our analysis identifies 720,843 homozygous discrepancies between new, high throughput genomic sequence data generated for Twilight and the EquCab2 reference assembly. Most of these represent errors in the assembly, while approximately 10,000 are demonstrated to be contributions from another horse. Other results are presented that include the binary alignment map file of the mapped Sanger reads, a list of variants identified as discrepancies between the source data and resulting reference, and a BED annotation file that lists the regions of the genome whose consensus was likely derived from low coverage alignments. |
url |
http://europepmc.org/articles/PMC4479572?pdf=render |
work_keys_str_mv |
AT jovanrebolledomendez comparisonoftheequinereferencesequencewithitssangersourcedataandnewilluminareads AT matthewshestand comparisonoftheequinereferencesequencewithitssangersourcedataandnewilluminareads AT stephenjcoleman comparisonoftheequinereferencesequencewithitssangersourcedataandnewilluminareads AT zhengzeng comparisonoftheequinereferencesequencewithitssangersourcedataandnewilluminareads AT ludovicorlando comparisonoftheequinereferencesequencewithitssangersourcedataandnewilluminareads AT jamesnmacleod comparisonoftheequinereferencesequencewithitssangersourcedataandnewilluminareads AT tedkalbfleisch comparisonoftheequinereferencesequencewithitssangersourcedataandnewilluminareads |
_version_ |
1716801208259182592 |