Genetic Analysis Workshop Research Articles

Homozygosity disequilibrium and its gene regulation.

Homozygosity disequilibrium (HD) describes a nonrandom pattern of sizable runs of homozygosity (ROH) that deviated from a random distribution of homozygotes and heterozygotes in the genome. In this study, we developed a double-weight local polynomial model for estimating homozygosity intensity. This new estimation method enables considering the local property and genetic information of homozygosity in the human genome when detecting regions of HD. By using this new method, we estimated whole-genome homozygosity intensities by analyzing real whole-genome sequencing data of 959 related individuals from 20 large pedigrees provided by Genetic Analysis Workshop 19 (GAW19). Through the analysis, we derived the distribution of HD in the human genome and provided evidence for the genetic component of natural variation in HD. Generalized estimating equation analysis for 855 related individuals was performed to identify regions of HD associated with diastolic blood pressure (DBP), systolic blood pressure, and hypertension (HTN), with concomitant adjustment for age and sex. We identified one DBP-associated and 2 HTN-associated regions of HD. We also studied the gene regulation of HD by analyzing the real whole-genome transcription data of 647 individuals. A set of gene expressions regulated by the DBP- and HTN-associated regions of HD was identified. Finally, we conducted simulation studies to evaluate the performance of our homozygosity association test. The results showed that the association test had a high power and that type 1 error was controlled. The methods have been integrated into our developed Loss-of-Heterozygosity Analysis Suite software, which can be downloaded at http://www.stat.sinica.edu.tw/hsinchou/genetics/loh/LOHAS.htm.

Genome-wide association of trajectories of systolic blood pressure change.

BackgroundThere is great interindividual variation in systolic blood pressure (SBP) as a result of the influences of several factors, including sex, ancestry, smoking status, medication use, and, especially, age. The majority of genetic studies have examined SBP measured cross-sectionally; however, SBP changes over time, and not necessarily in a linear fashion. Therefore, this study conducted a genome-wide association (GWA) study of SBP change trajectories using data available through the Genetic Analysis Workshop 19 (GAW19) of 959 individuals from 20 extended Mexican American families from the San Antonio Family Studies with up to 4 measures of SBP. We performed structural equation modeling (SEM) while taking into account potential genetic effects to identify how, if at all, to include covariates in estimating the SBP change trajectories using a mixture model based latent class growth modeling (LCGM) approach for use in the GWA analyses.ResultsThe semiparametric LCGM approach identified 5 trajectory classes that captured SBP changes across age. Each LCGM identified trajectory group was ranked based on the average number of cumulative years as hypertensive. Using a pairwise comparison of these classes the heritability estimates range from 12 to 94 % (SE = 17 to 40 %).ConclusionThese identified trajectories are significantly heritable, and we identified a total of 8 promising loci that influence one’s trajectory in SBP change across age. Our results demonstrate the potential utility of capitalizing on extant genetic data and longitudinal SBP assessments available through GAW19 to explore novel analytical methods with promising results.

Constrained multivariate association with longitudinal phenotypes.

BackgroundThe incorporation of longitudinal data into genetic epidemiological studies has the potential to provide valuable information regarding the effect of time on complex disease etiology. Yet, the majority of research focuses on variables collected from a single time point. This aim of this study was to test for main effects on a quantitative trait across time points using a constrained maximum-likelihood measured genotype approach. This method simultaneously accounts for all repeat measurements of a phenotype in families. We applied this method to systolic blood pressure (SBP) measurements from three time points using the Genetic Analysis Workshop 19 (GAW19) whole-genome sequence family simulated data set and 200 simulated replicates. Data consisted of 849 individuals from 20 extended Mexican American pedigrees. Comparisons were made among 3 statistical approaches: (a) constrained, where the effect of a variant or gene region on the mean trait value was constrained to be equal across all measurements; (b) unconstrained, where the variant or gene region effect was estimated separately for each time point; and (c) the average SBP measurement from three time points. These approaches were run for nine genetic variants with known effect sizes (>0.001) for SBP variability and a known gene-centric kernel (MAP4)-based test under the GAW19 simulation model across 200 replicates.ResultsWhen compared to results using two time points, the constrained method utilizing all 3 time points increased power to detect association. Averaging SBP was equally effective when the variant has a large effect on the phenotype, but less powerful for variants with lower effect sizes. However, averaging SBP was far more effective than either the constrained or unconstrained approaches when using a gene-centric kernel-based test.ConclusionWe determined that this constrained multivariate approach improves genetic signal over the bivariate method. However, this method is still only effective in those variants that explain a moderate to large proportion of the phenotypic variance but is not as effective for gene-centric tests.

Analysis of gene expression to predict dynamics of future hypertension incidence in type 2 diabetic patients.

BackgroundThe main focus of the Genetic Analysis Workshop 19 (GAW19) is identification of genes related to the occurrence of hypertension in the cohort of patients with type 2 diabetes mellitus (T2DM). The aim of our study was to predict dynamics of the future hypertension incidence, based on gene expression profiles, systolic and diastolic blood pressure changes in time, sex, baseline age, and cigarette smoking status.We analyzed data made available to GAW19 participants, which included gene expression profiles of peripheral blood mononuclear cells (PBMCs) from the diabetic members of 20 Mexican American families.MethodsOn the basis of mid blood pressure measurements at several time points, the coefficient of regression (slope) was calculated for each individual. We corrected the slope value in patients treated with antihypertensive medications. Feature preprocessing methods were used to remove highly correlated probes and linear dependencies between them. Subsequently, multiple linear regression model was used to associate gene expression with the regression coefficient calculated for each T2DM patient. Tenfold cross-validation was used to validate the model. We used linear mixed effects model and kinship coefficients to account for the family structure. All calculations were performed in R.ResultsThis analysis allowed us to identify 6 well-annotated genes: RTP4, FXYD6, GDF11, IFNAR1, NOX3, and HLA-DQ2, associated with dynamics of future hypertension incidence. Two of them, IFNAR1 and NOX3 were previously implicated in pathogenesis of hypertension.ConclusionsThere is no obvious mechanism that links all detected genes with dynamics of hypertension incidence. Identification of possible connection with hypertension needs further investigation.

Identification of interactions using model-based multifactor dimensionality reduction.

BackgroundCommon complex traits may involve multiple genetic and environmental factors and their interactions. Many methods have been proposed to identify these interaction effects, among them several machine learning and data mining methods. These are attractive for identifying interactions because they do not rely on specific genetic model assumptions. To handle the computational burden arising from an exhaustive search, including all possible combinations of factors, filter methods try to select promising factors in advance.MethodsModel-based multifactor dimensionality reduction (MB-MDR), a semiparametric machine learning method allowing adjustment for confounding variables and lower level effects, is applied to Genetic Analysis Workshop 19 (GAW19) data to identify interaction effects on different traits. Several filtering methods based on the nearest neighbor algorithm are assessed in terms of compatibility with MB-MDR.ResultsSingle nucleotide polymorphism (SNP) rs859400 shows a significant interaction effect (corrected p value <0.05) with age on systolic blood pressure (SBP). We identified 23 SNP–SNP interaction effects on hypertension status (HS), 42 interaction effects on SBP, and 26 interaction effects on diastolic blood pressure (DBP). Several of these SNPs are in strong linkage disequilibrium (LD). Three of the interaction effects on HS are identified in filtered subsets.ConclusionsThe considered filtering methods seem not to be appropriate to use with MB-MDR. LD pruning is further quality control to be incorporated, which can reduce the combinatorial burden by removing redundant SNPs.

Comparing family-based rare variant association tests for dichotomous phenotypes.

BackgroundIt has been repeatedly stressed that family-based samples suffer less from genetic heterogeneity and that association analyses with family-based samples are expected to be powerful for detecting susceptibility loci for rare disease. Various approaches for rare-variant analysis with family-based samples have been proposed.MethodsIn this report, performances of the existing methods were compared with the simulated data set provided as part of Genetic Analysis Workshop 19 (GAW19). We considered the rare variant transmission disequilibrium test (RV-TDT), generalized estimating equations-based kernel association (GEE-KM) test, an extended combined multivariate and collapsing test for pedigree data (known as Pedigree Combined Multivariate and Collapsing [PedCMC]), gene-level kernel and burden association tests with disease status for pedigree data (PedGene), and the family-based rare variant association test (FARVAT).ResultsThe results show that PedGene and FARVAT are usually the most efficient, and the optimal test statistic provided by FARVAT is robust under different disease models. Furthermore, FARVAT was implemented with C++, which is more computationally faster than other methods.ConclusionsConsidering both statistical and computational efficiency, we conclude that FARVAT is a good choice for rare-variant analysis with extended families.

Network-guided interaction mining for the blood pressure phenotype of unrelated individuals in genetic analysis workshop 19.

Interactions between genes are an important part of the genetic architecture of complex diseases. In this paper, we use literature-guided individual genes known to be associated with type 2 diabetes (referred to as “seed genes”) to create a larger list of genes that share implied or direct networks with these seed genes. This larger list of genes are known to interact with each other, but whether they interact in ways to influence hypertension in individuals presents an interesting question. Using Genetic Analysis Workshop data on individuals with diabetes, for which only case-control labels of hypertension are known, we offer a foray into identification of diabetes-related gene interactions that are associated with hypertension. We use the approach of Lo et al. (Proc Natl Acad Sci U S A 105: 12387-12392, 2008), which creates a score to identify pairwise significant gene associations. We find that the genes GCK and PAX4, formerly known to be found within similar coexpression and pathway networks but without specific direct interactions, do, in fact, show significant joint interaction effects for hypertension.

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

BackgroundThe 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.MethodsGAW19 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.

Type I error rates of rare single nucleotide variants are inflated in tests of association with non-normally distributed traits using simple linear regression methods.

In this study, the effects of (a) the minor allele frequency of the single nucleotide variant (SNV), (b) the degree of departure from normality of the trait, and (c) the position of the SNVs on type I error rates were investigated in the Genetic Analysis Workshop (GAW) 19 whole exome sequence data. To test the distribution of the type I error rate, 5 simulated traits were considered: standard normal and gamma distributed traits; 2 transformed versions of the gamma trait (log10 and rank-based inverse normal transformations); and trait Q1 provided by GAW 19. Each trait was tested with 313,340 SNVs. Tests of association were performed with simple linear regression and average type I error rates were determined for minor allele frequency classes. Rare SNVs (minor allele frequency < 0.05) showed inflated type I error rates for non–normally distributed traits that increased as the minor allele frequency decreased. The inflation of average type I error rates increased as the significance threshold decreased. Normally distributed traits did not show inflated type I error rates with respect to the minor allele frequency for rare SNVs. There was no consistent effect of transformation on the uniformity of the distribution of the location of SNVs with a type I error.

Genetic association analysis based on a joint model of gene expression and blood pressure.

Recent work on genetic association studies suggests that much of the heritable variation in complex traits is unexplained, which indicates a need for using more biologically meaningful modeling approaches and appropriate statistical methods. In this study, we propose a biological framework and a corresponding statistical model incorporating multilevel biological measures, and illustrate it in the analysis of the real data provided by the Genetic Analysis Workshop (GAW) 19, which contains whole genome sequence (WGS), gene expression (GE), and blood pressure (BP) data. We investigate the direct effect of single-nucleotide variants (SNVs) on BP and GE, while considering the non-directional dependence between BP and GE, by using copula functions to jointly model BP and GE conditional on SNVs. We implement the method for analysis on a genome-wide scale, and illustrate it within an association analysis of 68,727 SNVs on chromosome 19 that lie in or around genes with available GE measures. Although there is no indication for inflated type I errors under the proposed method, our results show that the association tests have smaller p values than tests under univariate models for common and rare variants using single-variant tests and gene-based multimarker tests. Hence, considering multilevel biological measures and modeling the dependence structure between these measures by using a plausible graphical approach may lead to more informative findings than standard univariate tests of common variants and well-recognized gene-based rare variant tests.

A LASSO penalized regression approach for genome-wide association analyses using related individuals: application to the Genetic Analysis Workshop 19 simulated data.

We propose a novel LASSO (least absolute shrinkage and selection operator) penalized regression method used to analyze samples consisting of (potentially) related individuals. Developed in the context of linear mixed models, our method models the relatedness of individuals in the sample through a random effect whose covariance structure is a linear function of known matrices with elements combinations of the condensed coefficients of identity between the individuals in the sample. We implement our method to analyze the simulated family data provided by the 19th Genetic Analysis Workshop in an effort to identify loci regulating the simulated trait of systolic blood pressure. The analyses were performed with full knowledge of the simulation model. Our findings demonstrate that we can significantly reduce the rate of false positive signals by incorporating the relatedness of the study participants.

Genetic complexity at expression quantitative trait loci.

BackgroundIdentifying variants that regulate gene expression and delineating their genetic architecture is a critical next step in our endeavors to better understand the genetic etiology of complex diseases. The appropriate genomic tools are in place, and preliminary analytic strategies have been developed.MethodsHere we used Genetic Analysis Workshop (GAW) 19 data to investigate the genetic complexity of expression quantitative trait loci (eQTL), chromosomal regions likely to harbor regulatory elements responsible for gene expression. For this investigation, we analyzed the lymphocyte expression profiles of 653 individuals in 20 pedigrees who were also genotyped by single nucleotide polymorphism (SNP) arrays, followed by sequencing and imputation. We used these data to examine the degree of allelic heterogeneity, a contributor to genetic complexity at eQTL, by sequentially conditioning on the most significantly associated SNPs.ResultSOLAR (Sequential Oligogenic Linkage Analysis Routines)-MGA (measured genotype approach) and FaST-LMM (Factored Spectrally Transformed Linear Mixed Model) software allowed us to analyze pedigree data. The power and Type 1 error rates for single SNP association testing and multiple SNP sequential association testing were consistent for these programs. Sequential conditioning of the real expression data revealed substantial levels of allelic heterogeneity at the 2 eQTL examined, illustrating this feature of genetic complexity.ConclusionseQTL exhibit substantial genetic complexity among and within pedigrees.

Integrating multiple genomic data: sparse representation based biomarker selection for blood pressure.

BackgroundAlthough many genes have been implicated as hypertension candidates, to date, few studies have integrated different types of genomic data for the purpose of biomarker selection.MethodsApplying a newly proposed sparse representation based variable selection (SRVS) method to the Genetic Analysis Workshop19 data, we analyzed a combined data set consisting of 11522 gene expressions and 354893 single-nucleotide polymorphisms (SNPs) from 397 subjects (case/control: 151/246), with the aim to identify potential biomarkers for blood pressure using both gene expression measures and SNP data.ResultsAmong the top 1000 variables (SNPs/gene expressions = 575/425) selected, the bioinformatics analysis showed that 302 were plausibly associated with blood pressure. In addition, we identified 173 variables that were associated with body weight and 84 associated with left ventricular contractility. Together, 55.9 % of the top 1000 variables showed associations with blood pressure related phenotypes(SNP/gene expression =348/211).ConclusionsOur results support the feasibility of the SRVS algorithm in integrating multiple data sets of different structure for comprehensive analysis.

Comparison of 2 models for gene-environment interactions: an example of simulated gene-medication interactions on systolic blood pressure in family-based data.

BackgroundNearly half of adults in the United States who are diagnosed with hypertension use blood-pressure-lowering medications. Yet there is a large interindividual variability in the response to these medications. Two complementary gene–environment interaction methods have been published and incorporated into publicly available software packages to examine interaction effects, including whether genetic variants modify the association between medication use and blood pressure. The first approach uses a gene–environment interaction term to measure the change in outcome when both the genetic marker and medication are present (the “interaction model”). The second approach tests for effect-size differences between strata of an environmental exposure (the “med-diff” approach). However, no studies have quantitatively compared how these methods perform with respect to 1 or 2 degree of freedom (DF) tests or in family-based data sets. We evaluated these 2 approaches using simulated genotype–medication response interactions at 3 single nucleotide polymorphisms (SNPs) across a range of minor allele frequencies (MAFs 0.1–5.4 %) using the Genetic Analysis Workshop 19 family sample.ResultsThe estimated interaction effect sizes were on average larger in the interaction model approach compared to the med-diff approach. The true positive proportion was higher for the med-diff approach for SNPs less than 1 % MAF, but higher for the interaction model when common variants were evaluated (MAF >5 %). The interaction model produced lower false-positive proportions than expected (5 %) across a range of MAFs for both the 1DF and 2DF tests. In contrast, the med-diff approach produced higher but stable false-positive proportions around 5 % across MAFs for both tests.ConclusionsAlthough the 1DF tests both performed similarly for common variants, the interaction model estimated true interaction effects with less bias and higher true positive proportions than the med-diff approach. However, if rare variation (MAF <5 %) is of interest, our findings suggest that when convergence is achieved, the med-diff approach may estimate true interaction effects more conservatively and with less variability.

A general method for combining different family-based rare-variant tests of association to improve power and robustness of a wide range of genetic architectures.

Current rare-variant, gene-based tests of association often suffer from a lack of statistical power to detect genotype–phenotype associations as a result of a lack of prior knowledge of genetic disease models combined with limited observations of extremely rare causal variants in population-based samples. The use of pedigree data, in which rare variants are often more highly concentrated than in population-based data, has been proposed as 1 possible method for enhancing power. Methods for combining multiple gene-based tests of association into a single summary p value are a robust approach to different genetic architectures when little a priori knowledge is available about the underlying genetic disease model. To date, however, little consideration has been given to combining gene-based tests of association for the analysis of pedigree data. We propose a flexible framework for combining any number of family-based rare-variant tests of association into a single summary statistic and for assessing the significance of that statistic. We show that this approach maintains type I error and improves the robustness, to different genetic architectures, of the statistical power of family- and gene-based rare-variant tests through application to simulated phenotype data from Genetic Analysis Workshop 19.

Testing rare variants for hypertension using family-based tests with different weighting schemes.

Next-generation sequencing technology makes directly testing rare variants possible. However, existing statistical methods to detect common variants may not be optimal for testing rare variants because of allelic heterogeneity as well as the extreme rarity of individual variants. Recently, several statistical methods to detect associations of rare variants were developed, including population-based and family-based methods. Compared with population-based methods, family-based methods have more power and can prevent bias induced by population substructure. Both population-based and family-based methods for rare variant association studies are essentially testing the effect of a weighted combination of variants or its function. How to model the weights is critical for the testing power because the number of observations for any given rare variant is small and the multiple-test correction is more stringent for rare variants. We propose 4 weighting schemes for the family-based rare variants test (FBAT-v) to test for the effects of both rare and common variants across the genome. Applying FBAT-v with the proposed weighting schemes on the Genetic Analysis Workshop 19 family data indicates that the power of FBAT-v can be comparatively enhanced in most circumstances.

Prioritization of family member sequencing for the detection of rare variants.

BackgroundThe advent of affordable sequencing has enabled researchers to discover many variants contributing to disease, including rare variants. There are methods for determining the most informative individuals for sequencing, but the application of these methods is more complex when working with families. Sets of large families can be beneficial in finding rare variants, but it may be unfeasible to sequence all members of these family sets.MethodsUsing simulated data from the Genetic Analysis Workshop 19, we apply multiple regression to identify cases and controls. To find the best controls for each case, we used kinship coefficients to match within families. Selected cases and controls were analyzed for rare variants, collapsed by gene, associated with hypertension using the family-based rare variant association test (FARVAT).ResultsThe gene with the strongest simulated effect, MAP4, did not meet the Bonferroni corrected significance threshold. However, analysis of cases and controls using our selection method substantially improved the significance of MAP4, despite the reduction in sample size.ConclusionsTaking the additional steps to select the optimal cases and controls from large family data sets can help ensure that only informative individuals are included in analysis and may improve the ability to detect rare variants.

Identifying regions of disease-related variants in admixed populations with the summation partition approach.

We propose a new method for identifying disease-related regions of single nucleotide variants in recently admixed populations. We use principal component analysis to derive both global and local ancestry information. We then use the summation partition approach to search for disease-related regions based on both rare variants and the local ancestral information of each region. We demonstrate this method using individuals with high systolic blood pressure from a sample of unrelated Mexican American subjects provided in the 19th Genetic Analysis Workshop.

Estimating relationships between phenotypes and subjects drawn from admixed families.

BackgroundEstimating relationships among subjects in a sample, within family structures or caused by population substructure, is complicated in admixed populations. Inaccurate allele frequencies can bias both kinship estimates and tests for association between subjects and a phenotype. We analyzed the simulated and real family data from Genetic Analysis Workshop 19, and were aware of the simulation model.ResultsWe found that kinship estimation is more accurate when marker data include common variants whose frequencies are less variable across populations. Estimates of heritability and association vary with age for longitudinally measured traits. Accounting for local ancestry identified different true associations than those identified by a traditional approach. Principal components aid kinship estimation and tests for association, but their utility is influenced by the frequency of the markers used to generate them.ConclusionsAdmixed families can provide a powerful resource for detecting disease loci, as well as analytical challenges. Allele frequencies, although difficult to adequately estimate in admixed populations, have a strong impact on the estimation of kinship, ancestry, and association with phenotypes. Approaches that acknowledge population structure in admixed families outperform those which ignore it.

A multistep approach to single nucleotide polymorphism-set analysis: an evaluation of power and type I error of gene-based tests of association after pathway-based association tests.

The aggregation of functionally associated variants given a priori biological information can aid in the discovery of rare variants associated with complex diseases. Many methods exist that aggregate rare variants into a set and compute a single p value summarizing association between the set of rare variants and a phenotype of interest. These methods are often called gene-based, rare variant tests of association because the variants in the set are often all contained within the same gene. A reasonable extension of these approaches involves aggregating variants across an even larger set of variants (eg, all variants contained in genes within a pathway). Testing sets of variants such as pathways for association with a disease phenotype reduces multiple testing penalties, may increase power, and allows for straightforward biological interpretation. However, a significant variant-set association test does not indicate precisely which variants contained within that set are causal. Because pathways often contain many variants, it may be helpful to follow-up significant pathway tests by conducting gene-based tests on each gene in that pathway to narrow in on the region of causal variants. In this paper, we propose such a multistep approach for variant-set analysis that can also account for covariates and complex pedigree structure. We demonstrate this approach on simulated phenotypes from Genetic Analysis Workshop 19. We find generally better power for the multistep approach when compared to a more conventional, single-step approach that simply runs gene-based tests of association on each gene across the genome. Further work is necessary to evaluate the multistep approach on different data sets with different characteristics.