Plasma glycomics predict cardiovascular disease in patients with ART-controlled HIV infections.

Despite effective control of HIV infection with antiretroviral drugs, individuals with HIV have high incidences of secondary diseases. These sequelae, such as cardiovascular disease (CVD), are poorly understood and represent a major health burden. To date, predictive biomarkers of HIV-associated secondary disease have been elusive, making preventative clinical management essentially impossible. Here, we applied a newly developed and easy to deploy, multitarget, and high-throughput glycomic analysis to banked HIV+ human plasma samples to determine whether the glycome may include biomarkers that predict future HIV-associated cardiovascular events or CVD diagnoses. Using 324 patient samples, we identified a glycomic fingerprint that was predictive of future CVD events but independent of CD4 counts, diabetes, age, and birth sex, suggesting that the plasma glycome may serve as a biomarker for specific HIV-associated sequelae. Our findings constitute the discovery of novel glycan biomarkers that could classify patients with HIV with elevated risk for CVD and reveal the untapped prognostic potential of the plasma glycome in human disease.-Oswald, D. M., Sim, E. S., Baker, C., Farhan, O., Debanne, S. M., Morris, N. J., Rodriguez, B. G., Jones, M. B., Cobb, B. A. Plasma glycomics predict cardiovascular disease in patients with ART-controlled HIV infections.

Heritability and longitudinal outcomes of spelling skills in individuals with histories of early speech and language disorders.

This study examined the spelling skills in middle childhood and adolescence in individuals with histories of early childhood speech sound disorders (SSD) with and without language impairment (LI). Youth without such histories were also included (No SSD/LI group). The heritability of spelling skills at each age level was estimated. Children with SSD were classified as SSD-only, SSD with LI but without childhood apraxia of speech (SSD + LI/ No CAS), and CAS and LI (CAS + LI). The SSD-only group did not differ in spelling from the No SSD/LI group, suggesting that SSD-only did not increase risk for poor spelling. The SSD + LI/No CAS and CAS + LI groups had poorer spelling skills than the SSD-only and No SSD/LI groups. Spelling was associated with phonological awareness in the middle childhood and adolescent samples and with rapid automatized naming in the adolescent sample. Heritability of spelling skills was stronger in adolescence than in middle childhood. Differences in the correlates of spelling and in heritability at the two ages suggest developmental changes in the factors contributing to spelling.

strum: an R package for structural modeling of latent variables for general pedigrees.

BACKGROUND: Structural equation modeling (SEM) is an extremely general and powerful approach to account for measurement error and causal pathways when analyzing data, and it has been used in wide range of applied sciences. There are many commercial and freely available software packages for SEM. However, it is difficult to use any of the packages to analyze general pedigree data, and SEM packages for genetics are limited in their application. RESULTS: We present the new R package strum to serve the need of a suitable SEM software tool for genetic analysis. It implements a general framework for SEM within the context of general pedigree data. This context requires specialized considerations such as familial correlations and ascertainment. Our package is an extraordinarily flexible tool capable of modeling genetic association, linkage analysis, polygenic effects, shared environment, and ascertainment combined with confirmatory factor analysis and general SEM. It also provides a convenient tool for model visualization, and integrates tools for simulating pedigree data. The various features of this package are tested through a simulation study to evaluate performance, and our results show that strum is very reliable and robust in terms of the accuracy and coverage of parameter estimates. CONCLUSIONS: strum is a valuable new tool for genetic analysis. It can be easily used with general pedigree data, incorporating both measurement and structural models, giving it some significant advantages over other software packages. It also includes a built-in approach for handling ascertainment, a helpful integrated tool for genetic data simulation, and built-in tools for model visualization, providing a significant addition to biomedical research.

Power of single- vs. multi-marker tests of association.

Current genome-wide association studies still heavily rely on a single-marker strategy, in which each single nucleotide polymorphism (SNP) is tested individually for association with a phenotype. Although methods and software packages that consider multimarker models have become available, they have been slow to become widely adopted and their efficacy in real data analysis is often questioned. Based on conducting extensive simulations, here we endeavor to provide more insights into the performance of simple multimarker association tests as compared to single-marker tests. The results reveal the power advantage as well as disadvantage of the two- vs. the single-marker test. Power differentials depend on the correlation structure among tag SNPs, as well as that between tag SNPs and causal variants. A two-marker test has relatively better performance than single-marker tests when the correlation of the two adjacent markers is high. However, using HapMap data, two-marker tests tended to have a greater chance of being less powerful than single-marker tests, due to constraints on the number of actual possible haplotypes in the HapMap data. Yet, the average power difference was small whenever the one-marker test is more powerful, while there were many situations where the two-marker test can be much more powerful. These findings can be useful to guide analyses of future studies.

How meaningful are heritability estimates of liability?

It is commonly acknowledged that estimates of heritability from classical twin studies have many potential shortcomings. Despite this, in the post-GWAS era, these heritability estimates have come to be a continual source of interest and controversy. While the heritability estimates of a quantitative trait are subject to a number of biases, in this article we will argue that the standard statistical approach to estimating the heritability of a binary trait relies on some additional untestable assumptions which, if violated, can lead to badly biased estimates. The ACE liability threshold model assumes at its heart that each individual has an underlying liability or propensity to acquire the binary trait (e.g., disease), and that this unobservable liability is multivariate normally distributed. We investigated a number of different scenarios violating this assumption such as the existence of a single causal diallelic gene and the existence of a dichotomous exposure. For each scenario, we found that substantial asymptotic biases can occur, which no increase in sample size can remove. Asymptotic biases as much as four times larger than the true value were observed, and numerous cases also showed large negative biases. Additionally, regions of low bias occurred for specific parameter combinations. Using simulations, we also investigated the situation where all of the assumptions of the ACE liability model are met. We found that commonly used sample sizes can lead to biased heritability estimates. Thus, even if we are willing to accept the meaningfulness of the liability construct, heritability estimates under the ACE liability threshold model may not accurately reflect the heritability of this construct. The points made in this paper should be kept in mind when considering the meaningfulness of a reported heritability estimate for any specific disease.

A variance component based multi-marker association test using family and unrelated data.

BACKGROUND: Incorporating family data in genetic association studies has become increasingly appreciated, especially for its potential value in testing rare variants. We introduce here a variance-component based association test that can test multiple common or rare variants jointly using both family and unrelated samples. RESULTS: The proposed approach implemented in our R package aggregates or collapses the information across a region based on genetic similarity instead of genotype scores, which avoids the power loss when the effects are in different directions or have different association strengths. The method is also able to effectively leverage the LD information in a region and it can produce a test statistic with an adaptively estimated number of degrees of freedom. Our method can readily allow for the adjustment of non-genetic contributions to the familial similarity, as well as multiple covariates. CONCLUSIONS: We demonstrate through simulations that the proposed method achieves good performance in terms of Type I error control and statistical power. The method is implemented in the R package "fassoc", which provides a useful tool for data analysis and exploration.

Single-marker and two-marker association tests for unphased case-control genotype data, with a power comparison.

In case-control single nucleotide polymorphism (SNP) data, the allele frequency, Hardy Weinberg Disequilibrium, and linkage disequilibrium (LD) contrast tests are three distinct sources of information about genetic association. While all three tests are typically developed in a retrospective context, we show that prospective logistic regression models may be developed that correspond conceptually to the retrospective tests. This approach provides a flexible framework for conducting a systematic series of association analyses using unphased genotype data and any number of covariates. For a single stage study, two single-marker tests and four two-marker tests are discussed. The true association models are derived and they allow us to understand why a model with only a linear term will generally fit well for a SNP in weak LD with a causal SNP, whatever the disease model, but not for a SNP in high LD with a non-additive disease SNP. We investigate the power of the association tests using real LD parameters from chromosome 11 in the HapMap CEU population data. Among the single-marker tests, the allelic test has on average the most power in the case of an additive disease, but for dominant, recessive, and heterozygote disadvantage diseases, the genotypic test has the most power. Among the four two-marker tests, the Allelic-LD contrast test, which incorporates linear terms for two markers and their interaction term, provides the most reliable power overall for the cases studied. Therefore, our result supports incorporating an interaction term as well as linear terms in multi-marker tests.

A framework for structural equation models in general pedigrees.

BACKGROUND/AIMS: Structural Equation Modeling (SEM) is an analysis approach that accounts for both the causal relationships between variables and the errors associated with the measurement of these variables. In this paper, a framework for implementing structural equation models (SEMs) in family data is proposed. METHODS: This framework includes both a latent measurement model and a structural model with covariates. It allows for a wide variety of models, including latent growth curve models. Environmental, polygenic and other genetic variance components can be included in the SEM. Kronecker notation makes it easy to separate the SEM process from a familial correlation model. A limited information method of model fitting is discussed. We show how missing data and ascertainment may be handled. We give several examples of how the framework may be used. RESULTS: A simulation study shows that our method is computationally feasible, and has good statistical properties. CONCLUSION: Our framework may be used to build and compare causal models using family data without any genetic marker data. It also allows for a nearly endless array of genetic association and/or linkage tests. A preliminary Matlab program is available, and we are currently implementing a more complete and user-friendly R package.

Calculating asymptotic significance levels of the constrained likelihood ratio test with application to multivariate genetic linkage analysis.

The asymptotic distribution of the multivariate variance component linkage analysis likelihood ratio test has provoked some contradictory accounts in the literature. In this paper we confirm that some previous results are not correct by deriving the asymptotic distribution in one special case. It is shown that this special case is a good approximation to the distribution in many situations. We also introduce a new approach to simulating from the asymptotic distribution of the likelihood ratio test statistic in constrained testing problems. It is shown that this method is very efficient for small p-values, and is applicable even when the constraints are not convex. The method is related to a multivariate integration problem. We illustrate how the approach can be applied to multivariate linkage analysis in a simulation study. Some more philosophical issues relating to one-sided tests in variance components linkage analysis are discussed.

Model-Free Linkage Analysis of a Quantitative Trait.

Model-free methods of linkage analysis for quantitative traits are a class of easily implemented, computationally efficient and statistically robust approaches to searching for linkage to a quantitative trait. By "model-free" we refer to methods of linkage analysis that do not fully specify a genetic model (i.e., the causal allele frequency, and penetrance functions). In this chapter we briefly survey the methods that are available, and then we discuss the necessary steps to implement an analysis using the programs GENIBD, SIBPAL and RELPAL in the S.A.G.E. (Statistical Analysis for Genetic Epidemiology) software suite.

Structural Equation Modeling.

Structural equation modeling (SEM) is a multivariate statistical framework that is used to model complex relationships between directly observed and indirectly observed (latent) variables. SEM is a general framework that involves simultaneously solving systems of linear equations and encompasses other techniques such as regression, factor analysis, path analysis, and latent growth curve modeling. Recently, SEM has gained popularity in the analysis of complex genetic traits because it can be used to better analyze the relationships between correlated variables (traits), to model genes as latent variables as a function of multiple observed genetic variants, and to assess the association between multiple genetic variants and multiple correlated phenotypes of interest. Though the general SEM framework only allows for the analysis of independent observations, recent work has extended SEM for the analysis of data on general pedigrees. Here, we review the theory of SEM for both unrelated and family data, describe the available software for SEM, and provide examples of SEM analysis.

Genome-wide association study identifies three novel loci in Fuchs endothelial corneal dystrophy.

The structure of the cornea is vital to its transparency, and dystrophies that disrupt corneal organization are highly heritable. To understand the genetic aetiology of Fuchs endothelial corneal dystrophy (FECD), the most prevalent corneal disorder requiring transplantation, we conducted a genome-wide association study (GWAS) on 1,404 FECD cases and 2,564 controls of European ancestry, followed by replication and meta-analysis, for a total of 2,075 cases and 3,342 controls. We identify three novel loci meeting genome-wide significance (P<5 × 10-8): KANK4 rs79742895, LAMC1 rs3768617 and LINC00970/ATP1B1 rs1200114. We also observe an overwhelming effect of the established TCF4 locus. Interestingly, we detect differential sex-specific association at LAMC1, with greater risk in women, and TCF4, with greater risk in men. Combining GWAS results with biological evidence we expand the knowledge of common FECD loci from one to four, and provide a deeper understanding of the underlying pathogenic basis of FECD.

Structural equation modeling with latent variables for longitudinal blood pressure traits using general pedigrees.

Structural equation modeling (SEM) has been used in a wide range of applied sciences including genetic analysis. The recently developed R package, strum, implements a framework for SEM for general pedigree data. We explored different SEM techniques using strum to analyze the multivariate longitudinal data and to ultimately test the association of genotypes on blood pressure traits. The quantitative blood pressure (BP) traits, systolic BP (SBP) and diastolic BP (DBP) were analyzed as the main traits of interest with age, sex, and smoking status as covariates. The single nucleotide polymorphism (SNP) genotype information from genome-wide association studies (GWAS) data was used for the test of association. The adjustment for hypertension treatment effect was done by the censored regression approach. Two different longitudinal data models, autoregressive model and latent growth curve model, were used to fit the longitudinal BP traits. The test of association for SNP was done using a novel score test within the SEM framework of strum. We found the 10 SNPs within the GWAS suggestive P value level, and among those 10, the most significant top 3 SNPs agreed in rank in both analysis models. The general SEM framework in strum is very useful to model and test for the association with massive genotype data and complex systems of multiple phenotypes with general pedigree data.

IV-treated pulmonary exacerbations in the prior year: An important independent risk factor for future pulmonary exacerbation in cystic fibrosis.

BACKGROUND: Single-center analyses have suggested that the number of CF pulmonary exacerbations (PEx) treated with intravenous antibiotics an individual has experienced in the prior year is significantly associated with their future PEx hazard. METHODS: We studied Prior-year PEx association with future PEx hazard by Cox proportional hazards regression among CF Foundation Patient Registry patients who experienced PEx after Jan 1, 2010. RESULTS: Among 13,579 patients, those with 1, 2, 3, or ≥4 Prior-year PEx treated with intravenous antibiotics were at 1.8, 2.9, 4.8, and 8.7 higher PEx hazard vs those without (P<.0001). Adjustment with significant demographic and clinical covariates (univariate P≤.0001) reduced Prior-year PEx hazard ratios to 1.6, 2.4, 3.6, and 6.0 (P<.0001). No other covariates had adjusted hazard ratios of >1.7. CONCLUSIONS: Prior-year PEx strongly associate with future PEx hazard and should be accounted for in prospective trials where treatment-associated change in PEx hazard is an efficacy outcome.

Model-free linkage analysis of a quantitative trait.

Model-free methods of linkage analysis for quantitative traits are a class of easily implemented, computationally efficient, and statistically robust approaches to searching for linkage to a quantitative trait. By "model-free" we refer to methods of linkage analysis that do not fully specify a genetic model (i.e., the causal allele frequency and penetrance functions). In this chapter, we briefly survey the methods that are available, and then we discuss the necessary steps to implement an analysis using the programs GENIBD, SIBPAL, and RELPAL in the Statistical Analysis for Genetic Epidemiology (S.A.G.E.) software suite.

A Parametric Survival Model When a Covariate is Subject to Left-Censoring.

PROBLEM STATEMENT: Modeling survival data with a set of covariates usually assumes that the values of the covariates are fully observed. However, in a variety of applications, some values of a covariate may be left-censored due to inadequate instrument sensitivity to quantify the biospecimen. When data are left-censored, the true values are missing but are known to be smaller than the detection limit. The most commonly used ad-hoc method to deal with nondetect values is to substitute the nondetect values by the detection limit. Such ad-hoc analysis of survival data with an explanatory variable subject to left-censoring may provide biased and inefficient estimators of hazard ratios and survivor functions. METHOD: We consider a parametric proportional hazards model to analyze time-to-event data. We propose a likelihood method for the estimation and inference of model parameters. In this likelihood approach, instead of replacing the nondetect values by the detection limit, we adopt a numerical integration technique to evaluate the observed data likelihood in the presence of a left-censored covariate. Monte Carlo simulations were used to demonstrate various properties of the proposed regression estimators including the consistency and efficiency. RESULTS: The simulation study shows that the proposed likelihood approach provides approximately unbiased estimators of the model parameters. The proposed method also provides estimators that are more efficient than those obtained under the ad-hoc method. Also, unlike the ad-hoc estimators, the coverage probabilities of the proposed estimators are at their nominal level. Analysis of a large cohort study, genetic and inflammatory marker of sepsis study, shows discernibly different results based on the proposed method. CONCLUSION: Naive use of detection limit in a parametric survival model may provide biased and inefficient estimators of hazard ratios and survivor functions. The proposed likelihood approach provides approximately unbiased and efficient estimators of hazard ratios and survivor functions.

Structural equation modeling.

Structural equation modeling (SEM) is a multivariate statistical framework that is used to model complex relationships between directly and indirectly observed (latent) variables. SEM is a general framework that involves simultaneously solving systems of linear equations and encompasses other techniques such as regression, factor analysis, path analysis, and latent growth curve modeling. Recently, SEM has gained popularity in the analysis of complex genetic traits because it can be used to better analyze the relationships between correlated variables (traits), to model genes as latent variables as a function of multiple observed genetic variants, and assess the association between multiple genetic variants and multiple correlated phenotypes of interest. Though the general SEM framework only allows for the analysis of independent observations, recent work has extended SEM for the analysis of general pedigrees. Here, we review the theory of SEM for both unrelated and family data, the available software for SEM, and provide an example of SEM analysis.

Testing gene-environment interactions in gene-based association studies.

Gene-based and single-nucleotide polymorphism (SNP) set association studies provide an important complement to SNP analysis. Kernel-based nonparametric regression has recently emerged as a powerful and flexible tool for this purpose. Our goal is to explore whether this approach can be extended to incorporate and test for interaction effects, especially for genes containing rare variant SNPs. Here, we construct nonparametric regression models that can be used to include a gene-environment interaction effect under the framework of the least-squares kernel machine and examine the performance of the proposed method on the Genetic Analysis Workshop 17 unrelated individuals data set. Two hundred simulated replicates were used to explore the power for detecting interaction. We demonstrate through a genome scan of the quantitative phenotype Q1 that the simulated gene-environment interaction effect in the data can be detected with reasonable power by using the least-squares kernel machine method.

Mendelian randomization in family data.

The phrase "mendelian randomization" has become associated with the use of genetic polymorphisms to uncover causal relationships between phenotypic variables. The statistical methods useful in mendelian randomization are known as instrumental variable techniques. We present an approach to instrumental variable estimation that is useful in family data and is robust to the use of weak instruments. We illustrate our method to measure the causal influence of low-density lipoprotein on high-density lipoprotein, body mass index, triglycerides, and systolic blood pressure. We use the Framingham Heart Study data as distributed to participants in the Genetics Analysis Workshop 16.

Comparison of univariate and multivariate linkage analysis of traits related to hypertension.

Complex traits are often manifested by multiple correlated traits. One example of this is hypertension (HTN), which is measured on a continuous scale by systolic blood pressure (SBP). Predisposition to HTN is predicted by hyperlipidemia, characterized by elevated triglycerides (TG), low-density lipids (LDL), and high-density lipids (HDL). We hypothesized that the multivariate analysis of TG, LDL, and HDL would be more powerful for detecting HTN genes via linkage analysis compared with univariate analysis of SBP. We conducted linkage analysis of four chromosomal regions known to contain genes associated with HTN using SBP as a measure of HTN in univariate Haseman-Elston regression and using the correlated traits TG, LDL, and HDL in multivariate Haseman-Elston regression. All analyses were conducted using the Framingham Heart Study data. We found that multivariate linkage analysis was better able to detect chromosomal regions in which the angiotensinogen, angiotensin receptor, guanine nucleotide-binding protein 3, and prostaglandin I2 synthase genes reside. Univariate linkage analysis only detected the AGT gene. We conclude that multivariate analysis is appropriate for the analysis of multiple correlated phenotypes, and our findings suggest that it may yield new linkage signals undetected by univariate analysis.