Progress in methods for rare variant association

Emerging roles of rare and low-frequency genetic variants in type 1 diabetes mellitus

Type 1 diabetes mellitus (T1DM) is defined as an autoimmune disorder and has enormous complexity and heterogeneity. Although its precise pathogenic mechanisms are obscure, this disease is widely acknowledged to be precipitated by environmental factors in individuals with genetic susceptibility. To date, the known susceptibility loci, which have mostly been identified by genome-wide association studies, can explain 80%–85% of the heritability of T1DM. Researchers believe that at least a part of its missing genetic component is caused by undetected rare and low-frequency variants. Most common variants have only small to modest effect sizes, which increases the difficulty of dissecting their functions and restricts their potential clinical application. Intriguingly, many studies have indicated that rare and low-frequency variants have larger effect sizes and play more significant roles in susceptibility to common diseases, including T1DM, than common variants do. Therefore, better recognition of rare and low-frequency variants is beneficial for revealing the genetic architecture of T1DM and for providing new and potent therapeutic targets for this disease. Here, we will discuss existing challenges as well as the great significance of this field and review current knowledge of the contributions of rare and low-frequency variants to T1DM.

A quantitative trait rare variant nonparametric linkage method with application to age-at-onset of Alzheimer's disease

To analyze pedigrees with quantitative trait (QT) and sequence data, we developed a rare variant (RV) quantitative nonparametric linkage (QNPL) method, which evaluates sharing of minor alleles. RV-QNPL has greater power than the traditional QNPL that tests for excess sharing of minor and major alleles. RV-QNPL is robust to population substructure and admixture, locus heterogeneity, and inclusion of nonpathogenic variants and can be readily applied outside of coding regions. When QNPL was used to analyze common variants, it often led to loci mapping to large intervals, e.g., >40 Mb. In contrast, when RVs are analyzed, regions are well defined, e.g., a gene. Using simulation studies, we demonstrate that RV-QNPL is substantially more powerful than applying traditional QNPL methods to analyze RVs. RV-QNPL was also applied to analyze age-at-onset (AAO) data for 107 late-onset Alzheimer's disease (LOAD) pedigrees of Caribbean Hispanic and European ancestry with whole-genome sequence data. When AAO of AD was analyzed regardless of APOE ε4 status, suggestive linkage (LOD = 2.4) was observed with RVs in KNDC1 and nominally significant linkage (p < 0.05) was observed with RVs in LOAD genes ABCA7 and IQCK. When AAO of AD was analyzed for APOE ε4 positive family members, nominally significant linkage was observed with RVs in APOE, while when AAO of AD was analyzed for APOE ε4 negative family members, nominal significance was observed for IQCK and ADAMTS1. RV-QNPL provides a powerful resource to analyze QTs in families to elucidate their genetic etiology.

Gene-based analyses of the maternal genome implicate maternal effect genes as risk factors for conotruncal heart defects

Congenital heart defects (CHDs) affect approximately 1% of newborns. Epidemiological studies have identified several genetically-mediated maternal phenotypes (e.g., pregestational diabetes, chronic hypertension) that are associated with the risk of CHDs in offspring. However, the role of the maternal genome in determining CHD risk has not been defined. We present findings from gene-level, genome-wide studies that link CHDs to maternal effect genes as well as to maternal genes related to hypertension and proteostasis. Maternal effect genes, which provide the mRNAs and proteins in the oocyte that guide early embryonic development before zygotic gene activation, have not previously been implicated in CHD risk. Our findings support a role for and suggest new pathways by which the maternal genome may contribute to the development of CHDs in offspring.

Gene-based genome-wide association studies and meta-analyses of conotruncal heart defects

Conotruncal heart defects (CTDs) are among the most common and severe groups of congenital heart defects. Despite evidence of an inherited genetic contribution to CTDs, little is known about the specific genes that contribute to the development of CTDs. We performed gene-based genome-wide analyses using microarray-genotyped and imputed common and rare variants data from two large studies of CTDs in the United States. We performed two case-parent trio analyses (N = 640 and 317 trios), using an extension of the family-based multi-marker association test, and two case-control analyses (N = 482 and 406 patients and comparable numbers of controls), using a sequence kernel association test. We also undertook two meta-analyses to combine the results from the analyses that used the same approach (i.e. family-based or case-control). To our knowledge, these analyses are the first reported gene-based, genome-wide association studies of CTDs. Based on our findings, we propose eight CTD candidate genes (ARF5, EIF4E, KPNA1, MAP4K3, MBNL1, NCAPG, NDFUS1 and PSMG3). Four of these genes (ARF5, KPNA1, NDUFS1 and PSMG3) have not been previously associated with normal or abnormal heart development. In addition, our analyses provide additional evidence that genes involved in chromatin-modification and in ribonucleic acid splicing are associated with congenital heart defects.

MULTILAYER KNOCKOFF FILTER: CONTROLLED VARIABLE SELECTION AT MULTIPLE RESOLUTIONS

We tackle the problem of selecting from among a large number of variables those that are "important" for an outcome. We consider situations where groups of variables are also of interest. For example, each variable might be a genetic polymorphism, and we might want to study how a trait depends on variability in genes, segments of DNA that typically contain multiple such polymorphisms. In this context, to discover that a variable is relevant for the outcome implies discovering that the larger entity it represents is also important. To guarantee meaningful results with high chance of replicability, we suggest controlling the rate of false discoveries for findings at the level of individual variables and at the level of groups. Building on the knockoff construction of Barber and Candès [Ann. Statist. 43 (2015) 2055-2085] and the multilayer testing framework of Barber and Ramdas [J. Roy. Statist. Soc. Ser. B 79 (2017) 1247-1268], we introduce the multilayer knockoff filter (MKF). We prove that MKF simultaneously controls the FDR at each resolution and use simulations to show that it incurs little power loss compared to methods that provide guarantees only for the discoveries of individual variables. We apply MKF to analyze a genetic dataset and find that it successfully reduces the number of false gene discoveries without a significant reduction in power.

A Zoom-Focus algorithm (ZFA) to locate the optimal testing region for rare variant association tests

Motivation

Increasing amounts of whole exome or genome sequencing data present the challenge of analysing rare variants with extremely small minor allele frequencies. Various statistical tests have been proposed, which are specifically configured to increase power for rare variants by conducting the test within a certain bin, such as a gene or a pathway. However, a gene may contain from several to thousands of markers, and not all of them are related to the phenotype. Combining functional and non-functional variants in an arbitrary genomic region could impair the testing power.

Results

We propose a Zoom-Focus algorithm (ZFA) to locate the optimal testing region within a given genomic region. It can be applied as a wrapper function in existing rare variant association tests to increase testing power. The algorithm consists of two steps. In the first step, Zooming, a given genomic region is partitioned by an order of two, and the best partition is located. In the second step, Focusing, the boundaries of the zoomed region are refined. Simulation studies showed that ZFA substantially increased the statistical power of rare variants' tests, including the SKAT, SKAT-O, burden test and the W-test. The algorithm was applied on real exome sequencing data of hypertensive disorder, and identified biologically relevant genetic markers to metabolic disorders that were undetectable by a gene-based method. The proposed algorithm is an efficient and powerful tool to enhance the power of association study for whole exome or genome sequencing data.

Availability and Implementation

The ZFA software is available at: http://www2.ccrb.cuhk.edu.hk/statgene/software.html.

Contact

maggiew@cuhk.edu.hk or bzee@cuhk.edu.hk.

Supplementary information

Supplementary data are available at Bioinformatics online.

Four Genome-Wide Association Studies Identify New Extreme Longevity Variants

The search for the genetic determinants of extreme human longevity has been challenged by the phenotype's rarity and its nonspecific definition by investigators. To address these issues, we established a consortium of four studies of extreme longevity that contributed 2,070 individuals who survived to the oldest one percentile of survival for the 1900 U.S. birth year cohort. We conducted various analyses to discover longevity-associated variants (LAV) and characterized those LAVs that differentiate survival to extreme age at death (eSAVs) from those LAVs that become more frequent in centenarians because of mortality selection (eg, survival to younger years). The analyses identified new rare variants in chromosomes 4 and 7 associated with extreme survival and with reduced risk for cardiovascular disease and Alzheimer's disease. The results confirm the importance of studying truly rare survival to discover those combinations of common and rare variants associated with extreme longevity and longer health span.

Omics-squared: human genomic, transcriptomic and phenotypic data for genetic analysis workshop 19

Background

The Genetic Analysis Workshops (GAW) are a forum for development, testing, and comparison of statistical genetic methods and software. Each contribution to the workshop includes an application to a specified data set. Here we describe the data distributed for GAW19, which focused on analysis of human genomic and transcriptomic data.

Methods

GAW19 data were donated by the T2D-GENES Consortium and the San Antonio Family Heart Study and included whole genome and exome sequences for odd-numbered autosomes, measures of gene expression, systolic and diastolic blood pressures, and related covariates in two Mexican American samples. These two samples were a collection of 20 large families with whole genome sequence and transcriptomic data and a set of 1943 unrelated individuals with exome sequence. For each sample, simulated phenotypes were constructed based on the real sequence data. ‘Functional’ genes and variants for the simulations were chosen based on observed correlations between gene expression and blood pressure. The simulations focused primarily on additive genetic models but also included a genotype-by-medication interaction. A total of 245 genes were designated as ‘functional’ in the simulations with a few genes of large effect and most genes explaining < 1 % of the trait variation. An additional phenotype, Q1, was simulated to be correlated among related individuals, based on theoretical or empirical kinship matrices, but was not associated with any sequence variants. Two hundred replicates of the phenotypes were simulated. The GAW19 data are an expansion of the data used at GAW18, which included the family-based whole genome sequence, blood pressure, and simulated phenotypes, but not the gene expression data or the set of 1943 unrelated individuals with exome sequence.

Associations Between Common and Rare Exonic Genetic Variants and Serum Levels of 20 Cardiovascular-Related Proteins: The Tromsø Study

Supplemental Digital Content is available in the text.

Background—

Genetic variation can be used to study causal relationships between biomarkers and diseases. Here, we identify new common and rare genetic variants associated with cardiovascular-related protein levels (protein quantitative trait loci [pQTLs]). We functionally annotate these pQTLs, predict and experimentally confirm a novel molecular interaction, and determine which pQTLs are associated with diseases and physiological phenotypes.

Methods and Results—

As part of a larger case–control study of venous thromboembolism, serum levels of 51 proteins implicated in cardiovascular diseases were measured in 330 individuals from the Tromsø Study. Exonic genetic variation near each protein’s respective gene (cis) was identified using sequencing and arrays. Using single site and gene-based tests, we identified 27 genetic associations between pQTLs and the serum levels of 20 proteins: 14 associated with common variation in cis, of which 6 are novel (ie, not previously reported); 7 associations with rare variants in cis, of which 4 are novel; and 6 associations in trans. Of the 20 proteins, 15 were associated with single sites and 7 with rare variants. cis-pQTLs for kallikrein and F12 also show trans associations for proteins (uPAR, kininogen) known to be cleaved by kallikrein and with NTproBNP. We experimentally demonstrate that kallikrein can cleave proBNP (NTproBNP precursor) in vitro. Nine of the pQTLs have previously identified associations with 17 disease and physiological phenotypes.

Conclusions—

We have identified cis and trans genetic variation associated with the serum levels of 20 proteins and utilized these pQTLs to study molecular mechanisms underlying disease and physiological phenotypes.