Is My Network Module Preserved and Reproducible?

In many applications, one is interested in determining which of the properties of a network module change across conditions. For example, to validate the existence of a module, it is desirable to show that it is reproducible (or preserved) in an independent test network. Here we study several types of network preservation statistics that do not require a module assignment in the test network. We distinguish network preservation statistics by the type of the underlying network. Some preservation statistics are defined for a general network (defined by an adjacency matrix) while others are only defined for a correlation network (constructed on the basis of pairwise correlations between numeric variables). Our applications show that the correlation structure facilitates the definition of particularly powerful module preservation statistics. We illustrate that evaluating module preservation is in general different from evaluating cluster preservation. We find that it is advantageous to aggregate multiple preservation statistics into summary preservation statistics. We illustrate the use of these methods in six gene co-expression network applications including 1) preservation of cholesterol biosynthesis pathway in mouse tissues, 2) comparison of human and chimpanzee brain networks, 3) preservation of selected KEGG pathways between human and chimpanzee brain networks, 4) sex differences in human cortical networks, 5) sex differences in mouse liver networks. While we find no evidence for sex specific modules in human cortical networks, we find that several human cortical modules are less preserved in chimpanzees. In particular, apoptosis genes are differentially co-expressed between humans and chimpanzees. Our simulation studies and applications show that module preservation statistics are useful for studying differences between the modular structure of networks. Data, R software and accompanying tutorials can be downloaded from the following webpage: http://www.genetics.ucla.edu/labs/horvath/CoexpressionNetwork/ModulePreservation

P2Y6 receptor potentiates pro-inflammatory responses in macrophages and exhibits differential roles in atherosclerotic lesion development.

BACKGROUND: P2Y(6), a purinergic receptor for UDP, is enriched in atherosclerotic lesions and is implicated in pro-inflammatory responses of key vascular cell types and macrophages. Evidence for its involvement in atherogenesis, however, has been lacking. Here we use cell-based studies and three murine models of atherogenesis to evaluate the impact of P2Y(6) deficiency on atherosclerosis. METHODOLOGY/PRINCIPAL FINDINGS: Cell-based studies in 1321N1 astrocytoma cells, which lack functional P2Y(6) receptors, showed that exogenous expression of P2Y(6) induces a robust, receptor- and agonist-dependent secretion of inflammatory mediators IL-8, IL-6, MCP-1 and GRO1. P2Y(6)-mediated inflammatory responses were also observed, albeit to a lesser extent, in macrophages endogenously expressing P2Y(6) and in acute peritonitis models of inflammation. To evaluate the role of P2Y(6) in atherosclerotic lesion development, we used P2Y(6)-deficient mice in three mouse models of atherosclerosis. A 43% reduction in aortic arch plaque was observed in high fat-fed LDLR knockout mice lacking P2Y(6) receptors in bone marrow-derived cells. In contrast, no effect on lesion development was observed in fat-fed whole body P2Y(6)xLDLR double knockout mice. Interestingly, in a model of enhanced vascular inflammation using angiotensin II, P2Y(6) deficiency enhanced formation of aneurysms and exhibited a trend towards increased atherosclerosis in the aorta of LDLR knockout mice. CONCLUSIONS: P2Y(6) receptor augments pro-inflammatory responses in macrophages and exhibits a pro-atherogenic role in hematopoietic cells. However, the overall impact of whole body P2Y(6) deficiency on atherosclerosis appears to be modest and could reflect additional roles of P2Y(6) in vascular disease pathophysiologies, such as aneurysm formation

The genetic architecture of NAFLD among inbred strains of mice

To identify genetic and environmental factors contributing to the pathogenesis of non-alcoholic fatty liver disease, we examined liver steatosis and related clinical and molecular traits in more than 100 unique inbred mouse strains, which were fed a diet rich in fat and carbohydrates. A >30-fold variation in hepatic TG accumulation was observed among the strains. Genome-wide association studies revealed three loci associated with hepatic TG accumulation. Utilizing transcriptomic data from the liver and adipose tissue, we identified several high-confidence candidate genes for hepatic steatosis, including Gde1, a glycerophosphodiester phosphodiesterase not previously implicated in triglyceride metabolism. We confirmed the role of Gde1 by in vivo hepatic over-expression and shRNA knockdown studies. We hypothesize that Gde1 expression increases TG production by contributing to the production of glycerol-3-phosphate. Our multi-level data, including transcript levels, metabolite levels, and gut microbiota composition, provide a framework for understanding genetic and environmental interactions underlying hepatic steatosis. DOI: http://dx.doi.org/10.7554/eLife.05607.00

Systems Genetics Analysis of Gene-by-Environment Interactions in Human Cells

Gene by environment (GxE) interactions are clearly important in many human diseases, but they have proven to be difficult to study on a molecular level. We report genetic analysis of thousands of transcript abundance traits in human primary endothelial cell (EC) lines in response to proinflammatory oxidized phospholipids implicated in cardiovascular disease. Of the 59 most regulated transcripts, approximately one-third showed evidence of GxE interactions. The interactions resulted primarily from effects of distal-, trans-acting loci, but a striking example of a local-GxE interaction was also observed for FGD6. Some of the distal interactions were validated by siRNA knockdown experiments, including a locus involved in the regulation of multiple transcripts involved in the ER stress pathway. Our findings add to the understanding of the overall architecture of complex human traits and are consistent with the possibility that GxE interactions are responsible, in part, for the failure of association studies to more fully explain common disease variation

Genetic Regulation of Adipose Gene Expression and Cardio-Metabolic Traits

Subcutaneous adipose tissue stores excess lipids and maintains energy balance. We performed expression quantitative trait locus (eQTL) analyses by using abdominal subcutaneous adipose tissue of 770 extensively phenotyped participants of the METSIM study. We identified cis-eQTLs for 12,400 genes at a 1% false-discovery rate. Among an approximately 680 known genome-wide association study (GWAS) loci for cardio-metabolic traits, we identified 140 coincident cis-eQTLs at 109 GWAS loci, including 93 eQTLs not previously described. At 49 of these 140 eQTLs, gene expression was nominally associated (p < 0.05) with levels of the GWAS trait. The size of our dataset enabled identification of five loci associated (p < 5 × 10(−8)) with at least five genes located >5 Mb away. These trans-eQTL signals confirmed and extended the previously reported KLF14-mediated network to 55 target genes, validated the CIITA regulation of class II MHC genes, and identified ZNF800 as a candidate master regulator. Finally, we observed similar expression-clinical trait correlations of genes associated with GWAS loci in both humans and a panel of genetically diverse mice. These results provide candidate genes for further investigation of their potential roles in adipose biology and in regulating cardio-metabolic traits