Two novel quantitative trait linkage analysis statistics based on the posterior probability of linkage: application to the COGA families

BACKGROUND: In this paper we apply two novel quantitative trait linkage statistics based on the posterior probability of linkage (PPL) to chromosome 4 from the GAW 14 COGA dataset. Our approaches are advantageous since they use the full likelihood, use full phenotypic information, do not assume normality at the population level or require population/sample parameter estimates; and like other forms of the PPL, they are specifically tailored to accumulate linkage evidence, either for or against linkage, across multiple sets of heterogeneous data. RESULTS: The first statistic uses all quantitative trait (QT) information from the pedigree (QT-posterior probability of linkage, PPL); we applied the QT-PPL to the trait ecb21 (resting electroencephalogram). The second statistic allows simultaneous incorporation of dichotomous trait data into the QT analysis via a threshold model (QTT-PPL); we applied the QTT-PPL to combined data on ecb21 and ALDX1. We obtained a QT-PPL of 96% at GABRB1 and a QT-PPL of 18% at FABP2 while the QTT-PPL was 4% and 2% at the same two loci, respectively. By comparison, the variance-components (VC) method, as implemented in SOLAR, yielded multipoint VC LOD scores of 2.05 and 2.21 at GABRB1 and FABP2, respectively; no other VC LODs were greater than 2. CONCLUSION: The QTT-PPL was only 4% at GABARB1, which might suggest that the underlying ecb21 gene does not also cause ALDX1, although features of the data complicate interpretation of this result

Evidence for the multiple hits genetic theory for inherited language impairment: a case study

Communication disorders have complex genetic origins, with constellations of relevant gene markers that vary across individuals. Some genetic variants are present in healthy individuals as well as those affected by developmental disorders. Growing evidence suggests that some variants may increase susceptibility to these disorders in the presence of other pathogenic gene mutations. In the current study, we describe eight children with specific language impairment and four of these children had a copy number variant in one of these potential susceptibility regions on chromosome 15. Three of these four children also had variants in other genes previously associated with language impairment. Our data support the theory that 15q11.2 is a susceptibility region for developmental disorders, specifically language impairment.University of Nebraska. Health Research ConsortiumBarkley Memorial Trus

De novo protein design:How do we expand into the universe of possible protein structures?

Protein scientists are paving the way to a new phase in protein design and engineering. Approaches and methods are being developed that could allow the design of proteins beyond the confines of natural protein structures. This possibility of designing entirely new proteins opens new questions: What do we build? How do we build into protein-structure space where there are few, if any, natural structures to guide us? To what uses can the resulting proteins be put? And, what, if anything, does this pursuit tell us about how natural proteins fold, function and evolve? We describe the origins of this emerging area of fully de novo protein design, how it could be developed, where it might lead, and what challenges lie ahead

Addressing the unmet need for visualizing Conditional Random Fields in Biological Data

Background: The biological world is replete with phenomena that appear to be ideally modeled and analyzed by one archetypal statistical framework - the Graphical Probabilistic Model (GPM). The structure of GPMs is a uniquely good match for biological problems that range from aligning sequences to modeling the genome-to-phenome relationship. The fundamental questions that GPMs address involve making decisions based on a complex web of interacting factors. Unfortunately, while GPMs ideally fit many questions in biology, they are not an easy solution to apply. Building a GPM is not a simple task for an end user. Moreover, applying GPMs is also impeded by the insidious fact that the complex web of interacting factors inherent to a problem might be easy to define and also intractable to compute upon. Discussion: We propose that the visualization sciences can contribute to many domains of the bio-sciences, by developing tools to address archetypal representation and user interaction issues in GPMs, and in particular a variety of GPM called a Conditional Random Field(CRF). CRFs bring additional power, and additional complexity, because the CRF dependency network can be conditioned on the query data. Conclusions: In this manuscript we examine the shared features of several biological problems that are amenable to modeling with CRFs, highlight the challenges that existing visualization and visual analytics paradigms induce for these data, and document an experimental solution called StickWRLD which, while leaving room for improvement, has been successfully applied in several biological research projects.Comment: BioVis 2014 conferenc

Validation of a Cost-Efficient Multi-Purpose SNP Panel for Disease Based Research

BACKGROUND: Here we present convergent methodologies using theoretical calculations, empirical assessment on in-house and publicly available datasets as well as in silico simulations, that validate a panel of SNPs for a variety of necessary tasks in human genetics disease research before resources are committed to larger-scale genotyping studies on those samples. While large-scale well-funded human genetic studies routinely have up to a million SNP genotypes, samples in a human genetics laboratory that are not yet part of such studies may be productively utilized in pilot projects or as part of targeted follow-up work though such smaller scale applications require at least some genome-wide genotype data for quality control purposes such as DNA "barcoding" to detect swaps or contamination issues, determining familial relationships between samples and correcting biases due to population effects such as population stratification in pilot studies. PRINCIPAL FINDINGS: Empirical performance in classification of relative types for any two given DNA samples (e.g., full siblings, parental, etc) indicated that for outbred populations the panel performs sufficiently to classify relationship in extended families and therefore also for smaller structures such as trios and for twin zygosity testing. Additionally, familial relationships do not significantly diminish the (mean match) probability of sharing SNP genotypes in pedigrees, further indicating the uniqueness of the "barcode." Simulation using these SNPs for an African American case-control disease association study demonstrated that population stratification, even in complex admixed samples, can be adequately corrected under a range of disease models using the SNP panel. CONCLUSION: The panel has been validated for use in a variety of human disease genetics research tasks including sample barcoding, relationship verification, population substructure detection and statistical correction. Given the ease of genotyping our specific assay contained herein, this panel represents a useful and economical panel for human geneticists

An eQTL biological data visualization challenge and approaches from the visualization community

In 2011, the IEEE VisWeek conferences inaugurated a symposium on Biological Data Visualization. Like other domain-oriented Vis symposia, this symposium's purpose was to explore the unique characteristics and requirements of visualization within the domain, and to enhance both the Visualization and Bio/Life-Sciences communities by pushing Biological data sets and domain understanding into the Visualization community, and well-informed Visualization solutions back to the Biological community. Amongst several other activities, the BioVis symposium created a data analysis and visualization contest. Unlike many contests in other venues, where the purpose is primarily to allow entrants to demonstrate tour-de-force programming skills on sample problems with known solutions, the BioVis contest was intended to whet the participants' appetites for a tremendously challenging biological domain, and simultaneously produce viable tools for a biological grand challenge domain with no extant solutions. For this purpose expression Quantitative Trait Locus (eQTL) data analysis was selected. In the BioVis 2011 contest, we provided contestants with a synthetic eQTL data set containing real biological variation, as well as a spiked-in gene expression interaction network influenced by single nucleotide polymorphism (SNP) DNA variation and a hypothetical disease model. Contestants were asked to elucidate the pattern of SNPs and interactions that predicted an individual's disease state. 9 teams competed in the contest using a mixture of methods, some analytical and others through visual exploratory methods. Independent panels of visualization and biological experts judged entries. Awards were given for each panel's favorite entry, and an overall best entry agreed upon by both panels. Three special mention awards were given for particularly innovative and useful aspects of those entries. And further recognition was given to entries that correctly answered a bonus question about how a proposed "gene therapy" change to a SNP might change an individual's disease status, which served as a calibration for each approaches' applicability to a typical domain question. In the future, BioVis will continue the data analysis and visualization contest, maintaining the philosophy of providing new challenging questions in open-ended and dramatically underserved Bio/Life Sciences domains

Robust Online Hamiltonian Learning

In this work we combine two distinct machine learning methodologies, sequential Monte Carlo and Bayesian experimental design, and apply them to the problem of inferring the dynamical parameters of a quantum system. We design the algorithm with practicality in mind by including parameters that control trade-offs between the requirements on computational and experimental resources. The algorithm can be implemented online (during experimental data collection), avoiding the need for storage and post-processing. Most importantly, our algorithm is capable of learning Hamiltonian parameters even when the parameters change from experiment-to-experiment, and also when additional noise processes are present and unknown. The algorithm also numerically estimates the Cramer-Rao lower bound, certifying its own performance.Comment: 24 pages, 12 figures; to appear in New Journal of Physic

A cGAS-dependent response links DNA damage and senescence in alveolar epithelial cells:A potential drug target in IPF

Alveolar epithelial cell (AEC) senescence is implicated in the pathogenesis of idiopathic pulmonary fibrosis (IPF). Mitochondrial dysfunction including release of mitochondrial DNA (mtDNA) is a feature of senescence, which led us to investigate the role of the DNA-sensing GMP-AMP synthase (cGAS) in IPF, with a focus on AEC senescence. cGAS expression in fibrotic tissue from lungs of IPF patients was detected within cells immunoreactive for epithelial cell adhesion molecule (EpCAM) and p21, epithelial and senescence markers respectively. Submerged primary cultures of AECs isolated from lung tissue of IPF patients (IPF-AECs, n=5) exhibited higher baseline senescence than AECs from control donors (Ctrl-AECs, n=5-7), as assessed by increased nuclear histone 2AXγ phosphorylation, p21 mRNA and expression of senescence-associated secretory phenotype (SASP) cytokines. Pharmacological cGAS inhibition using RU.521 diminished IPF-AEC senescence in culture and attenuated induction of Ctrl-AEC senescence following etoposide-induced DNA damage. Short interfering RNA (siRNA) knockdown of cGAS also attenuated etoposide-induced senescence of the AEC line, A549. Higher levels of mtDNA were detected in the cytosol and culture supernatants of primary IPF- and etoposide-treated Ctrl-AECs when compared to Ctrl-AECs at baseline. Furthermore, ectopic mtDNA augmented cGAS-dependent senescence of Ctrl-AECs, whereas DNAse I treatment diminished IPF-AEC senescence. This study provides evidence that a self-DNA driven, cGAS-dependent response augments AEC senescence, identifying cGAS as a potential therapeutic target for IPF