@article {525, title = {MendelProb: probability and sample size calculations for Mendelian studies of exome and whole genome sequence data.}, journal = {Bioinformatics}, volume = {35}, year = {2019}, month = {2019 02 01}, pages = {529-531}, abstract = {

Motivation: For the design of genetic studies, it is necessary to perform power calculations. Although for Mendelian traits the power of detecting linkage for pedigree(s) can be determined, it is also of great interest to determine the probability of identifying multiple pedigrees or unrelated cases with variants in the same gene. For many diseases, due to extreme locus heterogeneity this probability can be small. If only one family is observed segregating a variant classified as likely pathogenic or of unknown significance, the gene cannot be implicated in disease etiology. The probability of identifying several disease families or cases is dependent on the gene-specific disease prevalence and the sample size. The observation of multiple disease families or cases with variants in the same gene as well as evidence of pathogenicity from other sources, e.g. in silico prediction, expression and functional studies, can aid in implicating a gene in disease etiology. MendelProb can determine the probability of detecting a minimum number of families or cases with variants in the same gene. It can also calculate the probability of detecting genes with variants in different data types, e.g. identifying a variant in at least one family that can establish linkage and more the two additional families regardless of their size. Additionally, for a specified probability MendelProb can determine the number of probands which need to be screened to detect a minimum number of individuals with variants within the same gene.

Results: A single Mendelian disease family is not sufficient to implicate a gene in disease etiology. It is necessary to observe multiple families or cases with potentially pathogenic variants in the same gene. MendelProb, an R library, was developed to determine the probability of observing multiple families and cases with variants within a gene and to also establish the numbers of probands to screen to detect multiple observations of variants within a gene.

Availability and implementation: https://github.com/statgenetics/mendelprob.

}, keywords = {Exome, Genetic Linkage, Genomics, Humans, Pedigree, Probability, Sample Size, Software}, issn = {1367-4811}, doi = {10.1093/bioinformatics/bty542}, author = {He, Zongxiao and Wang, Lu and DeWan, Andrew T and Leal, Suzanne M} } @article {178, title = {Collapsed haplotype pattern method for linkage analysis of next-generation sequence data.}, journal = {Eur J Hum Genet}, volume = {23}, year = {2015}, month = {2015 Dec}, pages = {1739-43}, abstract = {

Recent advances in next-generation sequencing (NGS) make it possible to directly sequence genomes and exomes of individuals with Mendelian diseases and screen sequence data for causal variants. With the reduction in cost of NGS, DNA samples from entire families can be sequenced and linkage analysis can be performed directly using NGS data. Inspired by {\textquoteright}burden{\textquoteright} tests, which are used for complex trait rare variant association studies, we developed the collapsed haplotype pattern (CHP) method for linkage analysis. Using data from several deafness genes we demonstrate that the CHP method is substantially more powerful than analyzing individual variants. Unlike applying NGS data filtering approaches, the CHP method provides statistical evidence of a gene{\textquoteright}s involvement in disease etiology and is also less likely to exclude causal variants in the presence of phenocopies and/or reduced penetrance. The CHP method was implemented in the SEQLinkage software package, which can perform linkage analysis on NGS data or can generate data compatible with many linkage analysis programs, reviving them for use in NGS era.

}, keywords = {Chromosome Mapping, Genetic Linkage, Haplotypes, Hearing Loss, High-Throughput Nucleotide Sequencing, Humans, Sequence Analysis, DNA, Software}, issn = {1476-5438}, doi = {10.1038/ejhg.2015.64}, author = {Wang, Gao T and Zhang, Di and Li, Biao and Dai, Hang and Leal, Suzanne M} } @article {145, title = {Generation of sequence-based data for pedigree-segregating Mendelian or Complex traits.}, journal = {Bioinformatics}, volume = {31}, year = {2015}, month = {2015 Nov 15}, pages = {3706-8}, abstract = {

MOTIVATION: There is great interest in analyzing next generation sequence data that has been generated for pedigrees. However, unlike for population-based data there are only a limited number of rare variant methods to analyze pedigree data. One limitation is the ability to evaluate type I and II errors for family-based methods, due to lack of software that can simulate realistic sequence data for pedigrees.

SUMMARY: We developed RarePedSim (Rare-variant Pedigree-based Simulator), a program to simulate region/gene-level genotype and phenotype data for complex and Mendelian traits for any given pedigree structure. Using a genetic model, sequence variant data can be generated either conditionally or unconditionally on pedigree members{\textquoteright} qualitative or quantitative phenotypes. Additionally, qualitative or quantitative traits can be generated conditional on variant data. Sequence data can either be simulated using realistic population demographic models or obtained from sequence-based studies. Variant sites can be annotated with positions, allele frequencies and functionality. For rare variants, RarePedSim is the only program that can efficiently generate both genotypes and phenotypes, regardless of pedigree structure. Data generated by RarePedSim are in standard Linkage file (.ped) and Variant Call (.vcf) formats, ready to be used for a variety of purposes, including evaluation of type I error and power, for association methods including mixed models and linkage analysis methods.

AVAILABILITY AND IMPLEMENTATION: bioinformatics.org/simped/rare

CONTACT: sleal@bcm.edu.

}, keywords = {Base Sequence, Chromosome Segregation, Humans, Pedigree, Quantitative Trait, Heritable, Sequence Analysis, Software}, issn = {1367-4811}, doi = {10.1093/bioinformatics/btv412}, author = {Li, Biao and Wang, Gao T and Leal, Suzanne M} } @article {182, title = {Genetic linkage analysis in the age of whole-genome sequencing.}, journal = {Nat Rev Genet}, volume = {16}, year = {2015}, month = {2015 May}, pages = {275-84}, abstract = {

For many years, linkage analysis was the primary tool used for the genetic mapping of Mendelian and complex traits with familial aggregation. Linkage analysis was largely supplanted by the wide adoption of genome-wide association studies (GWASs). However, with the recent increased use of whole-genome sequencing (WGS), linkage analysis is again emerging as an important and powerful analysis method for the identification of genes involved in disease aetiology, often in conjunction with WGS filtering approaches. Here, we review the principles of linkage analysis and provide practical guidelines for carrying out linkage studies using WGS data.

}, keywords = {Algorithms, Female, Genetic Linkage, Genome, Human, Genome-Wide Association Study, Genotyping Techniques, Humans, Lod Score, Male, Pedigree, Penetrance, Phenotype, Polymorphism, Single Nucleotide, Sequence Analysis, DNA, Software}, issn = {1471-0064}, doi = {10.1038/nrg3908}, author = {Ott, Jurg and Wang, Jing and Leal, Suzanne M} } @article {250, title = {Power analysis and sample size estimation for sequence-based association studies.}, journal = {Bioinformatics}, volume = {30}, year = {2014}, month = {2014 Aug 15}, pages = {2377-8}, abstract = {

MOTIVATION: Statistical methods have been developed to test for complex trait rare variant (RV) associations, in which variants are aggregated across a region, which is typically a gene. Power analysis and sample size estimation for sequence-based RV association studies are challenging because of the necessity to realistically model the underlying allelic architecture of complex diseases within a suitable analytical framework to assess the performance of a variety of RV association methods in an unbiased manner.

SUMMARY: We developed SEQPower, a software package to perform statistical power analysis for sequence-based association data under a variety of genetic variant and disease phenotype models. It aids epidemiologists in determining the best study design, sample size and statistical tests for sequence-based association studies. It also provides biostatisticians with a platform to fairly compare RV association methods and to validate and assess novel association tests.

AVAILABILITY AND IMPLEMENTATION: The SEQPower program, source code, multi-platform executables, documentation, list of association tests, examples and tutorials are available at http://bioinformatics.org/spower.

}, keywords = {Alleles, Data Interpretation, Statistical, Disease, Genetic Association Studies, Genetic Variation, Humans, Phenotype, Sample Size, Sequence Analysis, DNA, Software}, issn = {1367-4811}, doi = {10.1093/bioinformatics/btu296}, author = {Wang, Gao T and Li, Biao and Santos-Cortez, Regie P Lyn and Peng, Bo and Leal, Suzanne M} } @article {249, title = {Variant association tools for quality control and analysis of large-scale sequence and genotyping array data.}, journal = {Am J Hum Genet}, volume = {94}, year = {2014}, month = {2014 May 01}, pages = {770-83}, abstract = {

Currently there is great interest in detecting associations between complex traits and rare variants. In this report, we describe Variant Association Tools (VAT) and the VAT pipeline, which implements best practices for rare-variant association studies. Highlights of VAT include variant-site and call-level quality control (QC), summary statistics, phenotype- and genotype-based sample selection, variant annotation, selection of variants for association analysis, and a collection of rare-variant association methods for analyzing qualitative and quantitative traits. The association testing framework for VAT is regression based, which readily allows for flexible construction of association models with multiple covariates and weighting themes based on allele frequencies or predicted functionality. Additionally, pathway analyses, conditional analyses, and analyses of gene-gene and gene-environment interactions can be performed. VAT is capable of rapidly scanning through data by using multi-process computation, adaptive permutation, and simultaneously conducting association analysis via multiple methods. Results are available in text or graphic file formats and additionally can be output to relational databases for further annotation and filtering. An interface to R language also facilitates user implementation of novel association methods. The VAT{\textquoteright}s data QC and association-analysis pipeline can be applied to sequence, imputed, and genotyping array, e.g., "exome chip," data, providing a reliable and reproducible computational environment in which to analyze small- to large-scale studies with data from the latest genotyping and sequencing technologies. Application of the VAT pipeline is demonstrated through analysis of data from the 1000 Genomes project.

}, keywords = {Genetic Association Studies, Genetic Variation, Genotyping Techniques, High-Throughput Nucleotide Sequencing, Humans, Multifactorial Inheritance, Quality Control, Software}, issn = {1537-6605}, doi = {10.1016/j.ajhg.2014.04.004}, author = {Wang, Gao T and Peng, Bo and Leal, Suzanne M} } @article {326, title = {SEQCHIP: a powerful method to integrate sequence and genotype data for the detection of rare variant associations.}, journal = {Bioinformatics}, volume = {28}, year = {2012}, month = {2012 Jul 01}, pages = {1745-51}, abstract = {

MOTIVATION: Next-generation sequencing greatly increases the capacity to detect rare-variant complex-trait associations. However, it is still expensive to sequence a large number of samples and therefore often small datasets are used. Given cost constraints, a potentially more powerful two-step strategy is to sequence a subset of the sample to discover variants, and genotype the identified variants in the remaining sample. If only cases are sequenced, directly combining sequence and genotype data will lead to inflated type-I errors in rare-variant association analysis. Although several methods have been developed to correct for the bias, they are either underpowered or theoretically invalid. We proposed a new method SEQCHIP to integrate genotype and sequence data, which can be used with most existing rare-variant tests.

RESULTS: It is demonstrated using both simulated and real datasets that the SEQCHIP method has controlled type-I errors, and is substantially more powerful than all other currently available methods.

AVAILABILITY: SEQCHIP is implemented in an R-Package and is available at http://linkage.rockefeller.edu/suzanne/seqchip/Seqchip.html.

}, keywords = {Adenoma, Case-Control Studies, Colorectal Neoplasms, Genetic Association Studies, Genetic Variation, Genotype, Humans, Phenotype, Sequence Analysis, DNA, Software}, issn = {1367-4811}, doi = {10.1093/bioinformatics/bts263}, author = {Liu, Dajiang J and Leal, Suzanne M} } @article {321, title = {SimRare: a program to generate and analyze sequence-based data for association studies of quantitative and qualitative traits.}, journal = {Bioinformatics}, volume = {28}, year = {2012}, month = {2012 Oct 15}, pages = {2703-4}, abstract = {

MOTIVATION: Currently, there is great interest in detecting complex trait rare variant associations using next-generation sequence data. On a monthly basis, new rare variant association methods are published. It is difficult to evaluate these methods because there is no standard to generate data and often comparisons are biased. In order to fairly compare rare variant association methods, it is necessary to generate data using realistic population demographic and phenotypic models.

RESULT: SimRare is an interactive program that integrates generation of rare variant genotype/phenotype data and evaluation of association methods using a unified platform. Variant data are generated for gene regions using forward-time simulation that incorporates realistic population demographic and evolutionary scenarios. Phenotype data can be obtained for both case-control and quantitative traits. SimRare has a user-friendly interface that allows for easy entry of genetic and phenotypic parameters. Novel rare variant association methods implemented in R can also be imported into SimRare, to evaluate their performance and compare results, e.g. power and Type I error, with other currently available methods both numerically and graphically.

}, keywords = {Genetic Association Studies, Genetic Variation, Humans, Phenotype, Quantitative Trait, Heritable, Software}, issn = {1367-4811}, doi = {10.1093/bioinformatics/bts499}, author = {Li, Biao and Wang, Gao and Leal, Suzanne M} } @article {316, title = {A unified method for detecting secondary trait associations with rare variants: application to sequence data.}, journal = {PLoS Genet}, volume = {8}, year = {2012}, month = {2012}, pages = {e1003075}, abstract = {

Next-generation sequencing has made possible the detection of rare variant (RV) associations with quantitative traits (QT). Due to high sequencing cost, many studies can only sequence a modest number of selected samples with extreme QT. Therefore association testing in individual studies can be underpowered. Besides the primary trait, many clinically important secondary traits are often measured. It is highly beneficial if multiple studies can be jointly analyzed for detecting associations with commonly measured traits. However, analyzing secondary traits in selected samples can be biased if sample ascertainment is not properly modeled. Some methods exist for analyzing secondary traits in selected samples, where some burden tests can be implemented. However p-values can only be evaluated analytically via asymptotic approximations, which may not be accurate. Additionally, potentially more powerful sequence kernel association tests, variable selection-based methods, and burden tests that require permutations cannot be incorporated. To overcome these limitations, we developed a unified method for analyzing secondary trait associations with RVs (STAR) in selected samples, incorporating all RV tests. Statistical significance can be evaluated either through permutations or analytically. STAR makes it possible to apply more powerful RV tests to analyze secondary trait associations. It also enables jointly analyzing multiple cohorts ascertained under different study designs, which greatly boosts power. The performance of STAR and commonly used RV association tests were comprehensively evaluated using simulation studies. STAR was also implemented to analyze a dataset from the SardiNIA project where samples with extreme low-density lipoprotein levels were sequenced. A significant association between LDLR and systolic blood pressure was identified, which is supported by pharmacogenetic studies. In summary, for sequencing studies, STAR is an important tool for detecting secondary-trait RV associations.

}, keywords = {Blood Pressure, Computer Simulation, Genetic Predisposition to Disease, Genetic Variation, Genome-Wide Association Study, High-Throughput Nucleotide Sequencing, Humans, Italy, Lipoproteins, LDL, Models, Genetic, Phenotype, Quantitative Trait Loci, Software}, issn = {1553-7404}, doi = {10.1371/journal.pgen.1003075}, author = {Liu, Dajiang J and Leal, Suzanne M} }