Testing Rare-Variant Association without Calling Genotypes Allows for Systematic Differences in Sequencing between Cases and Controls.
Next-generation sequencing of DNA provides an unprecedented opportunity to discover rare genetic variants associated with complex diseases and traits. However, the common practice of first calling underlying genotypes and then treating the called values as known is prone to false positive findings,...
| Main Authors: | , , , , |
|---|---|
| Format: | Article |
| Language: | English |
| Published: |
Public Library of Science (PLoS)
2016-05-01
|
| Series: | PLoS Genetics |
| Online Access: | http://europepmc.org/articles/PMC4859496?pdf=render |
| _version_ | 1818158771178831872 |
|---|---|
| author | Yi-Juan Hu Peizhou Liao H Richard Johnston Andrew S Allen Glen A Satten |
| author_facet | Yi-Juan Hu Peizhou Liao H Richard Johnston Andrew S Allen Glen A Satten |
| author_sort | Yi-Juan Hu |
| collection | DOAJ |
| description | Next-generation sequencing of DNA provides an unprecedented opportunity to discover rare genetic variants associated with complex diseases and traits. However, the common practice of first calling underlying genotypes and then treating the called values as known is prone to false positive findings, especially when genotyping errors are systematically different between cases and controls. This happens whenever cases and controls are sequenced at different depths, on different platforms, or in different batches. In this article, we provide a likelihood-based approach to testing rare variant associations that directly models sequencing reads without calling genotypes. We consider the (weighted) burden test statistic, which is the (weighted) sum of the score statistic for assessing effects of individual variants on the trait of interest. Because variant locations are unknown, we develop a simple, computationally efficient screening algorithm to estimate the loci that are variants. Because our burden statistic may not have mean zero after screening, we develop a novel bootstrap procedure for assessing the significance of the burden statistic. We demonstrate through extensive simulation studies that the proposed tests are robust to a wide range of differential sequencing qualities between cases and controls, and are at least as powerful as the standard genotype calling approach when the latter controls type I error. An application to the UK10K data reveals novel rare variants in gene BTBD18 associated with childhood onset obesity. The relevant software is freely available. |
| first_indexed | 2024-12-11T15:35:23Z |
| format | Article |
| id | doaj.art-4fa59e839d544eb2b396fef81bc0eb19 |
| institution | Directory Open Access Journal |
| issn | 1553-7390 1553-7404 |
| language | English |
| last_indexed | 2024-12-11T15:35:23Z |
| publishDate | 2016-05-01 |
| publisher | Public Library of Science (PLoS) |
| record_format | Article |
| series | PLoS Genetics |
| spelling | doaj.art-4fa59e839d544eb2b396fef81bc0eb192022-12-22T00:59:57ZengPublic Library of Science (PLoS)PLoS Genetics1553-73901553-74042016-05-01125e100604010.1371/journal.pgen.1006040Testing Rare-Variant Association without Calling Genotypes Allows for Systematic Differences in Sequencing between Cases and Controls.Yi-Juan HuPeizhou LiaoH Richard JohnstonAndrew S AllenGlen A SattenNext-generation sequencing of DNA provides an unprecedented opportunity to discover rare genetic variants associated with complex diseases and traits. However, the common practice of first calling underlying genotypes and then treating the called values as known is prone to false positive findings, especially when genotyping errors are systematically different between cases and controls. This happens whenever cases and controls are sequenced at different depths, on different platforms, or in different batches. In this article, we provide a likelihood-based approach to testing rare variant associations that directly models sequencing reads without calling genotypes. We consider the (weighted) burden test statistic, which is the (weighted) sum of the score statistic for assessing effects of individual variants on the trait of interest. Because variant locations are unknown, we develop a simple, computationally efficient screening algorithm to estimate the loci that are variants. Because our burden statistic may not have mean zero after screening, we develop a novel bootstrap procedure for assessing the significance of the burden statistic. We demonstrate through extensive simulation studies that the proposed tests are robust to a wide range of differential sequencing qualities between cases and controls, and are at least as powerful as the standard genotype calling approach when the latter controls type I error. An application to the UK10K data reveals novel rare variants in gene BTBD18 associated with childhood onset obesity. The relevant software is freely available.http://europepmc.org/articles/PMC4859496?pdf=render |
| spellingShingle | Yi-Juan Hu Peizhou Liao H Richard Johnston Andrew S Allen Glen A Satten Testing Rare-Variant Association without Calling Genotypes Allows for Systematic Differences in Sequencing between Cases and Controls. PLoS Genetics |
| title | Testing Rare-Variant Association without Calling Genotypes Allows for Systematic Differences in Sequencing between Cases and Controls. |
| title_full | Testing Rare-Variant Association without Calling Genotypes Allows for Systematic Differences in Sequencing between Cases and Controls. |
| title_fullStr | Testing Rare-Variant Association without Calling Genotypes Allows for Systematic Differences in Sequencing between Cases and Controls. |
| title_full_unstemmed | Testing Rare-Variant Association without Calling Genotypes Allows for Systematic Differences in Sequencing between Cases and Controls. |
| title_short | Testing Rare-Variant Association without Calling Genotypes Allows for Systematic Differences in Sequencing between Cases and Controls. |
| title_sort | testing rare variant association without calling genotypes allows for systematic differences in sequencing between cases and controls |
| url | http://europepmc.org/articles/PMC4859496?pdf=render |
| work_keys_str_mv | AT yijuanhu testingrarevariantassociationwithoutcallinggenotypesallowsforsystematicdifferencesinsequencingbetweencasesandcontrols AT peizhouliao testingrarevariantassociationwithoutcallinggenotypesallowsforsystematicdifferencesinsequencingbetweencasesandcontrols AT hrichardjohnston testingrarevariantassociationwithoutcallinggenotypesallowsforsystematicdifferencesinsequencingbetweencasesandcontrols AT andrewsallen testingrarevariantassociationwithoutcallinggenotypesallowsforsystematicdifferencesinsequencingbetweencasesandcontrols AT glenasatten testingrarevariantassociationwithoutcallinggenotypesallowsforsystematicdifferencesinsequencingbetweencasesandcontrols |