Modern simulation utilities for genetic analysis.

BACKGROUND: Statistical geneticists employ simulation to estimate the power of proposed studies, test new analysis tools, and evaluate properties of causal models. Although there are existing trait simulators, there is ample room for modernization. For example, most phenotype simulators are limited to Gaussian traits or traits transformable to normality, while ignoring qualitative traits and realistic, non-normal trait distributions. Also, modern computer languages, such as Julia, that accommodate parallelization and cloud-based computing are now mainstream but rarely used in older applications. To meet the challenges of contemporary big studies, it is important for geneticists to adopt new computational tools. RESULTS: We present TraitSimulation, an open-source Julia package that makes it trivial to quickly simulate phenotypes under a variety of genetic architectures. This package is integrated into our OpenMendel suite for easy downstream analyses. Julia was purpose-built for scientific programming and provides tremendous speed and memory efficiency, easy access to multi-CPU and GPU hardware, and to distributed and cloud-based parallelization. TraitSimulation is designed to encourage flexible trait simulation, including via the standard devices of applied statistics, generalized linear models (GLMs) and generalized linear mixed models (GLMMs). TraitSimulation also accommodates many study designs: unrelateds, sibships, pedigrees, or a mixture of all three. (Of course, for data with pedigrees or cryptic relationships, the simulation process must include the genetic dependencies among the individuals.) We consider an assortment of trait models and study designs to illustrate integrated simulation and analysis pipelines. Step-by-step instructions for these analyses are available in our electronic Jupyter notebooks on Github. These interactive notebooks are ideal for reproducible research. CONCLUSION: The TraitSimulation package has three main advantages. (1) It leverages the computational efficiency and ease of use of Julia to provide extremely fast, straightforward simulation of even the most complex genetic models, including GLMs and GLMMs. (2) It can be operated entirely within, but is not limited to, the integrated analysis pipeline of OpenMendel. And finally (3), by allowing a wider range of more realistic phenotype models, TraitSimulation brings power calculations and diagnostic tools closer to what investigators might see in real-world analyses.

Correction to: A candidate gene study of intermediate histopathological phenotypes in HIV-associated neurocognitive disorders.

Due to a production error data in Table 1 were not presented correctly.

A candidate gene study of intermediate histopathological phenotypes in HIV-associated neurocognitive disorders.

HIV-associated neurocognitive disorders (HAND) describe a spectrum of neuropsychological impairment caused by HIV-1 infection. While the sequence of cellular and physiological events that lead to HAND remains obscure, it likely involves chronic neuroinflammation. Host genetic markers that increase the risk for HAND have been reported, but replication of such studies is lacking, possibly due to inconsistent application of a behavioral phenotype across studies. In the current study, we used histopathologic phenotypes in order to validate putative risk alleles for HAND. The National NeuroAIDS Tissue Consortium, a longitudinal study of the neurologic manifestations of HIV. Data and specimens were obtained from 175 HIV-infected adults. After determining several potential covariates of neurocognitive functioning, we quantified levels of six histopathological markers in the frontal lobe in association with neurocognitive functioning: SYP, MAP 2, HLA-DR, Iba1, GFAP, and ß-amyloid. We then determined alleles of 15 candidate genes for their associations with neurocognitive functioning and histopathological markers. Finally, we identified the most plausible causal pathway based on our data using a multi-stage linear regression-based mediation analysis approach. None of the genetic markers were associated with neurocognitive functioning. Of the histopathological markers, only MAP 2 and SYP were associated with neurocognitive functioning; however, MAP 2 and SYP did not vary as a function of genotype. Mediation analysis suggests a causal pathway in which presynaptic degeneration (SYP) leads to somatodendritic degeneration (MAP 2) and ultimately neurocognitive impairment. This study did not support the role of host genotype in the histopathology underlying HAND. The findings lend further support for synaptodendritic degeneration as the proximal underlying neuropathological substrate of HAND.

OPENMENDEL: a cooperative programming project for statistical genetics.

Statistical methods for genome-wide association studies (GWAS) continue to improve. However, the increasing volume and variety of genetic and genomic data make computational speed and ease of data manipulation mandatory in future software. In our view, a collaborative effort of statistical geneticists is required to develop open source software targeted to genetic epidemiology. Our attempt to meet this need is called the OPENMENDEL project (https://openmendel.github.io). It aims to (1) enable interactive and reproducible analyses with informative intermediate results, (2) scale to big data analytics, (3) embrace parallel and distributed computing, (4) adapt to rapid hardware evolution, (5) allow cloud computing, (6) allow integration of varied genetic data types, and (7) foster easy communication between clinicians, geneticists, statisticians, and computer scientists. This article reviews and makes recommendations to the genetic epidemiology community in the context of the OPENMENDEL project.

Activity of Genes with Functions in Human Williams-Beuren Syndrome Is Impacted by Mobile Element Insertions in the Gray Wolf Genome.

In canines, transposon dynamics have been associated with a hyper-social behavioral syndrome, although the functional mechanism has yet to be described. We investigate the epigenetic and transcriptional consequences of these behavior-associated mobile element insertions (MEIs) in dogs and Yellowstone gray wolves. We posit that the transposons themselves may not be the causative feature; rather, their transcriptional regulation may exert the functional impact. We survey four outlier transposons associated with hyper-sociability, with the expectation that they are targeted for epigenetic silencing. We predict hyper-methylation of MEIs, suggestive that the epigenetic silencing of and not the MEIs themselves may be driving dysregulation of nearby genes. We found that transposon-derived sequences are significantly hyper-methylated, regardless of their copy number or species. Further, we have assessed transcriptome sequence data and found evidence that MEIs impact the expression levels of six genes (WBSCR17, LIMK1, GTF2I, WBSCR27, BAZ1B, and BCL7B), all of which have known roles in human Williams-Beuren syndrome due to changes in copy number, typically hemizygosity. Although further evidence is needed, our results suggest that a few insertions alter local expression at multiple genes, likely through a cis-regulatory mechanism that excludes proximal methylation.