2019
DOI: 10.1016/j.fsigen.2019.01.008
|View full text |Cite
|
Sign up to set email alerts
|

Mitigating the effects of reference sequence bias in single-multiplex massively parallel sequencing of the mitochondrial DNA control region

Abstract: Highlights mtDNA control region of 101 diverse samples amplified in a single reaction as 10 overlapping amplicons and sequenced via MPS. Primers create reference bias, compromising ability to call variants or heteroplasmy in primer-binding regions. Bioinformatic selection of overarching reads bypasses effects of proprietary primers and mitigates bias. Data processing permits accurate calling of variants, and heteroplasmies down to 5% lev… Show more

Help me understand this report

Search citation statements

Order By: Relevance

Paper Sections

Select...
2
1
1
1

Citation Types

0
4
0

Year Published

2019
2019
2024
2024

Publication Types

Select...
7
1

Relationship

0
8

Authors

Journals

citations
Cited by 9 publications
(6 citation statements)
references
References 42 publications
0
4
0
Order By: Relevance
“…The amplicon settings, which define the data ranges for each of the assay amplicons, were initially set at 16013-16126, 16116-16225, 16224-16408, 16387-16486, 16474-30, 16555-152, 136-257, 246-364, 342-436, and 429-592, based on consultation with the vendor. These ranges sit just inside the primer regions for each of the 10 amplicons of the PowerSeq mtDNA CR Nested System, as analyzing data within the PCR primer binding regions can produce false point heteroplasmies [30]. Upon examination of the mixture data, the amplicon 3 settings were updated to 16223-16408 (for the reasons described in Section 3.3, below).…”
Section: Discussionmentioning
confidence: 99%
“…The amplicon settings, which define the data ranges for each of the assay amplicons, were initially set at 16013-16126, 16116-16225, 16224-16408, 16387-16486, 16474-30, 16555-152, 136-257, 246-364, 342-436, and 429-592, based on consultation with the vendor. These ranges sit just inside the primer regions for each of the 10 amplicons of the PowerSeq mtDNA CR Nested System, as analyzing data within the PCR primer binding regions can produce false point heteroplasmies [30]. Upon examination of the mixture data, the amplicon 3 settings were updated to 16223-16408 (for the reasons described in Section 3.3, below).…”
Section: Discussionmentioning
confidence: 99%
“…The remaining positions were usually of low quality and were counted as false negatives, because they were missed in at least one sequence sample (see below). We reasoned that a reference allele bias effect could impact the quality of variant calling, especially in the highly polymorphic and repetitive control region [35]. Therefore, we performed the same variant calling procedure using two alternative references: a Thoroughbred breed more closely related to H25 and H27, and the Twilight horse (a Thoroughbred) used for the equine reference genome.…”
Section: Mitochondrial Dna Variants Callingmentioning
confidence: 99%
“…Some approaches to MSA primer design, like the JCVI primer designer, often make use of reference sequences ( Li et al 2008 , 2012 ). One problem with this approach is reference sequences introduce a primer matching bias within interspecies comparisons favoring sequences with the most matches to the given reference ( Huszar et al 2019 ). Within the human microbiome, where intraspecies genetic diversity and mutation rates may be higher, the reference sequence may not remain a valid basis of comparison for any particular experiment ( Ramiro et al 2020 ).…”
Section: Resultsmentioning
confidence: 99%