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

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.

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.

Behavioral and Molecular Genetics of Reading-Related AM and FM Detection Thresholds

Auditory detection thresholds for certain frequencies of both amplitude modulated (AM) and frequency modulated (FM) dynamic auditory stimuli are associated with reading in typically developing and dyslexic readers. We present the first behavioral and molecular genetic characterization of these two auditory traits. Two extant extended family datasets were given reading tasks and psychoacoustic tasks to determine FM 2 Hz and AM 20 Hz sensitivity thresholds. Univariate heritabilities were significant for both AM (h ²  = 0.20) and FM (h ²  = 0.29). Bayesian posterior probability of linkage (PPL) analysis found loci for AM (12q, PPL = 81 %) and FM (10p, PPL = 32 %; 20q, PPL = 65 %). Bivariate heritability analyses revealed that FM is genetically correlated with reading, while AM was not. Bivariate PPL analysis indicates that FM loci (10p, 20q) are not also associated with reading.

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.

Joint analysis of multiple phenotypes: summary of results and discussions from the Genetic Analysis Workshop 19

For Genetic Analysis Workshop 19, 2 extensive data sets were provided, including whole genome and whole exome sequence data, gene expression data, and longitudinal blood pressure outcomes, together with nongenetic covariates. These data sets gave researchers the chance to investigate different aspects of more complex relationships within the data, and the contributions in our working group focused on statistical methods for the joint analysis of multiple phenotypes, which is part of the research field of data integration. The analysis of data from different sources poses challenges to researchers but provides the opportunity to model the real-life situation more realistically.Our 4 contributions all used the provided real data to identify genetic predictors for blood pressure. In the contributions, novel multivariate rare variant tests, copula models, structural equation models and a sparse matrix representation variable selection approach were applied. Each of these statistical models can be used to investigate specific hypothesized relationships, which are described together with their biological assumptions.The results showed that all methods are ready for application on a genome-wide scale and can be used or extended to include multiple omics data sets. The results provide potentially interesting genetic targets for future investigation and replication. Furthermore, all contributions demonstrated that the analysis of complex data sets could benefit from modeling correlated phenotypes jointly as well as by adding further bioinformatics information.