Merging microsatellite data: enhanced methodology and software to combine genotype data for linkage and association analysis

<p>Abstract</p> <p>Background</p> <p>Correctly merged data sets that have been independently genotyped can increase statistical power in linkage and association studies. However, alleles from microsatellite data sets genotyped with different experimental protocols or platforms cannot be accurately matched using base-pair size information alone. In a previous publication we introduced a statistical model for merging microsatellite data by matching allele frequencies between data sets. These methods are implemented in our software MicroMerge version 1 (v1). While MicroMerge v1 output can be analyzed by some genetic analysis programs, many programs can not analyze alignments that do not match alleles one-to-one between data sets. A consequence of such alignments is that codominant genotypes must often be analyzed as phenotypes. In this paper we describe several extensions that are implemented in MicroMerge version 2 (v2).</p> <p>Results</p> <p>Notably, MicroMerge v2 includes a new one-to-one alignment option that creates merged pedigree and locus files that can be handled by most genetic analysis software. Other features in MicroMerge v2 enhance the following aspects of control: 1) optimizing the algorithm for different merging scenarios, such as data sets with very different sample sizes or multiple data sets, 2) merging small data sets when a reliable set of allele frequencies are available, and 3) improving the quantity and 4) quality of merged data. We present results from simulated and real microsatellite genotype data sets, and conclude with an association analysis of three familial dyslipidemia (FD) study samples genotyped at different laboratories. Independent analysis of each FD data set did not yield consistent results, but analysis of the merged data sets identified strong association at locus D11S2002.</p> <p>Conclusion</p> <p>The MicroMerge v2 features will enable merging for a variety of genotype data sets, which in turn will facilitate meta-analyses for powering association analysis.</p

Integrating Genetic and Network Analysis to Characterize Genes Related to Mouse Weight

Systems biology approaches that are based on the genetics of gene expression have been fruitful in identifying genetic regulatory loci related to complex traits. We use microarray and genetic marker data from an F2 mouse intercross to examine the large-scale organization of the gene co-expression network in liver, and annotate several gene modules in terms of 22 physiological traits. We identify chromosomal loci (referred to as module quantitative trait loci, mQTL) that perturb the modules and describe a novel approach that integrates network properties with genetic marker information to model gene/trait relationships. Specifically, using the mQTL and the intramodular connectivity of a body weight–related module, we describe which factors determine the relationship between gene expression profiles and weight. Our approach results in the identification of genetic targets that influence gene modules (pathways) that are related to the clinical phenotypes of interest

Niche adaptation by expansion and reprogramming of general transcription factors

Experimental analysis of TFB family proteins in a halophilic archaeon reveals complex environment-dependent fitness contributions. Gene conversion events among these proteins can generate novel niche adaptation capabilities, a process that may have contributed to archaeal adaptation to extreme environments

Neural G0:a quiescent-like state found in neuroepithelial-derived cells and glioma

Single‐cell RNA sequencing has emerged as a powerful tool for resolving cellular states associated with normal and maligned developmental processes. Here, we used scRNA‐seq to examine the cell cycle states of expanding human neural stem cells (hNSCs). From these data, we constructed a cell cycle classifier that identifies traditional cell cycle phases and a putative quiescent‐like state in neuroepithelial‐derived cell types during mammalian neurogenesis and in gliomas. The Neural G0 markers are enriched with quiescent NSC genes and other neurodevelopmental markers found in non‐dividing neural progenitors. Putative glioblastoma stem‐like cells were significantly enriched in the Neural G0 cell population. Neural G0 cell populations and gene expression are significantly associated with less aggressive tumors and extended patient survival for gliomas. Genetic screens to identify modulators of Neural G0 revealed that knockout of genes associated with the Hippo/Yap and p53 pathways diminished Neural G0 in vitro, resulting in faster G1 transit, down‐regulation of quiescence‐associated markers, and loss of Neural G0 gene expression. Thus, Neural G0 represents a dynamic quiescent‐like state found in neuroepithelial‐derived cells and gliomas

A nonsynonymous SNP within PCDH15 is associated with lipid traits in familial combined hyperlipidemia

Familial combined hyperlipidemia (FCHL) is a common lipid disorder characterized by the presence of multiple lipoprotein phenotypes that increase the risk of premature coronary heart disease. In a previous study, we identified an intragenic microsatellite marker within the protocadherin 15 (PCDH15) gene to be associated with high triglycerides (TGs) in Finnish dyslipidemic families. In this study we analyzed all four known nonsynonymous SNPs within PCDH15 in 1,268 individuals from Finnish and Dutch multigenerational families with FCHL. Association analyses of quantitative traits for SNPs were performed using the QTDT test. The nonsynonymous SNP rs10825269 resulted in a P = 0.0006 for the quantitative TG trait. Additional evidence for association was observed with the same SNP for apolipoprotein B levels (apo-B) (P = 0.0001) and total cholesterol (TC) levels (P = 0.001). None of the other three SNPs tested showed a significant association with any lipid-related trait. We investigated the expression of PCDH15 in different human tissues and observed that PCDH15 is expressed in several tissues including liver and pancreas. In addition, we measured the plasma lipid levels in mice with loss-of-function mutations in Pcdh15 (Pcdh15av-Tg and Pcdh15av-3J) to investigate possible abnormalities in their lipid profile. We observed a significant difference in plasma TG and TC concentrations for the Pcdh15av-3J carriers when compared with the wild type (P = 0.013 and P = 0.044, respectively). Our study suggests that PCDH15 is associated with lipid abnormalities

Single cell analysis reveals satellite cell heterogeneity for proinflammatory chemokine expression

Background: The expression of proinflammatory signals at the site of muscle injury are essential for efficient tissue repair and their dysregulation can lead to inflammatory myopathies. Macrophages, neutrophils, and fibroadipogenic progenitor cells residing in the muscle are significant sources of proinflammatory cytokines and chemokines. However, the inducibility of the myogenic satellite cell population and their contribution to proinflammatory signaling is less understood.Methods: Mouse satellite cells were isolated and exposed to lipopolysaccharide (LPS) to mimic sterile skeletal muscle injury and changes in the expression of proinflammatory genes was examined by RT-qPCR and single cell RNA sequencing. Expression patterns were validated in skeletal muscle injured with cardiotoxin by RT-qPCR and immunofluorescence.Results: Satellite cells in culture were able to express Tnfa, Ccl2, and Il6, within 2 h of treatment with LPS. Single cell RNA-Seq revealed seven cell clusters representing the continuum from activation to differentiation. LPS treatment led to a heterogeneous pattern of induction of C-C and C-X-C chemokines (e.g., Ccl2, Ccl5, and Cxcl0) and cytokines (e.g., Tgfb1, Bmp2, Il18, and Il33) associated with innate immune cell recruitment and satellite cell proliferation. One cell cluster was enriched for expression of the antiviral interferon pathway genes under control conditions and LPS treatment. Activation of this pathway in satellite cells was also detectable at the site of cardiotoxin induced muscle injury.Conclusion: These data demonstrate that satellite cells respond to inflammatory signals and secrete chemokines and cytokines. Further, we identified a previously unrecognized subset of satellite cells that may act as sensors for muscle infection or injury using the antiviral interferon pathway

COL4A1 Mutations Cause Ocular Dysgenesis, Neuronal Localization Defects, and Myopathy in Mice and Walker-Warburg Syndrome in Humans

Muscle-eye-brain disease (MEB) and Walker Warburg Syndrome (WWS) belong to a spectrum of autosomal recessive diseases characterized by ocular dysgenesis, neuronal migration defects, and congenital muscular dystrophy. Until now, the pathophysiology of MEB/WWS has been attributed to alteration in dystroglycan post-translational modification. Here, we provide evidence that mutations in a gene coding for a major basement membrane protein, collagen IV alpha 1 (COL4A1), are a novel cause of MEB/WWS. Using a combination of histological, molecular, and biochemical approaches, we show that heterozygous Col4a1 mutant mice have ocular dysgenesis, neuronal localization defects, and myopathy characteristic of MEB/WWS. Importantly, we identified putative heterozygous mutations in COL4A1 in two MEB/WWS patients. Both mutations occur within conserved amino acids of the triple-helix-forming domain of the protein, and at least one mutation interferes with secretion of the mutant proteins, resulting instead in intracellular accumulation. Expression and posttranslational modification of dystroglycan is unaltered in Col4a1 mutant mice indicating that COL4A1 mutations represent a distinct pathogenic mechanism underlying MEB/WWS. These findings implicate a novel gene and a novel mechanism in the etiology of MEB/WWS and expand the clinical spectrum of COL4A1-associated disorders

A Systems Genetics Approach Implicates USF1, FADS3, and Other Causal Candidate Genes for Familial Combined Hyperlipidemia

We hypothesized that a common SNP in the 3' untranslated region of the upstream transcription factor 1 (USF1), rs3737787, may affect lipid traits by influencing gene expression levels, and we investigated this possibility utilizing the Mexican population, which has a high predisposition to dyslipidemia. We first associated rs3737787 genotypes in Mexican Familial Combined Hyperlipidemia (FCHL) case/control fat biopsies, with global expression patterns. To identify sets of co-expressed genes co-regulated by similar factors such as transcription factors, genetic variants, or environmental effects, we utilized weighted gene co-expression network analysis (WGCNA). Through WGCNA in the Mexican FCHL fat biopsies we identified two significant Triglyceride (TG)-associated co-expression modules. One of these modules was also associated with FCHL, the other FCHL component traits, and rs3737787 genotypes. This USF1-regulated FCHL-associated (URFA) module was enriched for genes involved in lipid metabolic processes. Using systems genetics procedures we identified 18 causal candidate genes in the URFA module. The FCHL causal candidate gene fatty acid desaturase 3 (FADS3) was associated with TGs in a recent Caucasian genome-wide significant association study and we replicated this association in Mexican FCHL families. Based on a USF1-regulated FCHL-associated co-expression module and SNP rs3737787, we identify a set of causal candidate genes for FCHL-related traits. We then provide evidence from two independent datasets supporting FADS3 as a causal gene for FCHL and elevated TGs in Mexicans