A new expectation-maximization statistical test for case-control association studies considering rare variants...

Area of Science:

• Genetics
• Bioinformatics
• Statistical genetics

Background:

• Genome-wide association studies (GWAS) identify common variants but leave heritability unexplained.
• Rare genetic variants may contribute significantly to disease risk but are often missed by current GWAS.
• Decreasing sequencing costs enable direct analysis of rare variants in case-control studies.

Purpose of the Study:

• To develop a novel statistical method for association testing of rare variants using sequencing data.
• To account for potential errors in sequencing reads during genotype calling.
• To evaluate the performance of the new method in identifying disease-associated rare variants.

Main Methods:

• Developed a test statistic for association testing using raw sequencing reads.
• Employed the expectation-maximization algorithm to determine genotype probabilities from base pair reads.
• Applied the SumStat procedure to aggregate association signals across multiple rare variant loci.
• Validated the method through simulations, assessing type I error rates and statistical power.

Main Results:

• The developed method maintains correct type I error rates, even with sequencing read errors.
• The SumStat procedure demonstrates good statistical power across various simulated scenarios.
• The method's power is comparable to or better than single-locus approaches.

Conclusions:

• The new statistical method effectively detects associations involving rare variants from sequencing data.
• This approach addresses limitations of traditional GWAS in explaining genetic heritability.
• The SumStat procedure offers a robust and powerful tool for rare variant association studies.