Web-based Gene Pathogenicity Analysis (WGPA): a web platform to interpret gene pathogenicity from personal genome data

UNLABELLED: As the volume of patient-specific genome sequences increases the focus of biomedical research is switching from the detection of disease-mutations to their interpretation. To this end a number of techniques have been developed that use mutation data collected within a population to predict whether individual genes are likely to be disease-causing or not. As both sequence data and associated analysis tools proliferate, it becomes increasingly difficult for the community to make sense of these data and their implications. Moreover, no single analysis tool is likely to capture all relevant genomic features that contribute to the gene's pathogenicity. Here, we introduce Web-based Gene Pathogenicity Analysis (WGPA), a web-based tool to analyze genes impacted by mutations and rank them through the integration of existing prioritization tools, which assess different aspects of gene pathogenicity using population-level sequence data. Additionally, to explore the polygenic contribution of mutations to disease, WGPA implements gene set enrichment analysis to prioritize disease-causing genes and gene interaction networks, therefore providing a comprehensive annotation of personal genomes data in disease. AVAILABILITY AND IMPLEMENTATION: wgpa.systems-genetics.net

Genome-wide co-expression analysis in multiple tissues

Expression quantitative trait loci (eQTLs) represent genetic control points of gene expression, and can be categorized as cis- and trans-acting, reflecting local and distant regulation of gene expression respectively. Although there is evidence of co-regulation within clusters of trans-eQTLs, the extent of co-expression patterns and their relationship with the genotypes at eQTLs are not fully understood. We have mapped thousands of cis- and trans-eQTLs in four tissues (fat, kidney, adrenal and left ventricle) in a large panel of rat recombinant inbred (RI) strains. Here we investigate the genome-wide correlation structure in expression levels of eQTL transcripts and underlying genotypes to elucidate the nature of co-regulation within cis- and trans-eQTL datasets. Across the four tissues, we consistently found statistically significant correlations of cis-regulated gene expression to be rare (<0.9% of all pairs tested). Most (>80%) of the observed significant correlations of cis-regulated gene expression are explained by correlation of the underlying genotypes. In comparison, co-expression of trans-regulated gene expression is more common, with significant correlation ranging from 2.9%-14.9% of all pairs of trans-eQTL transcripts. We observed a total of 81 trans-eQTL clusters (hot-spots), defined as consisting of > or =10 eQTLs linked to a common region, with very high levels of correlation between trans-regulated transcripts (77.2-90.2%). Moreover, functional analysis of large trans-eQTL clusters (> or =30 eQTLs) revealed significant functional enrichment among genes comprising 80% of the large clusters. The results of this genome-wide co-expression study show the effects of the eQTL genotypes on the observed patterns of correlation, and suggest that functional relatedness between genes underlying trans-eQTLs is reflected in the degree of co-expression observed in trans-eQTL clusters. Our results demonstrate the power of an integrative, systematic approach to the analysis of a large gene expression dataset to uncover underlying structure, and inform future eQTL studies

Unique Regulatory Properties of Mesangial Cells Are Genetically Determined in the Rat

Mesangial cells are glomerular cells of stromal origin. During immune complex mediated crescentic glomerulonephritis (Crgn), infiltrating and proliferating pro-inflammatory macrophages lead to crescent formation. Here we have hypothesised that mesangial cells, given their mesenchymal stromal origin, show similar immunomodulatory properties as mesenchymal stem cells (MSCs), by regulating macrophage function associated with glomerular crescent formation. We show that rat mesangial cells suppress conA-stimulated splenocyte proliferation in vitro, as previously shown for MSCs. We then investigated mesangial cell-macrophage interaction by using mesangial cells isolated from nephrotoxic nephritis (NTN)-susceptible Wistar Kyoto (WKY) and NTN-resistant Lewis (LEW) rats. We first determined the mesangial cell transcriptome in WKY and LEW rats and showed that this is under marked genetic control. Supernatant transfer results show that WKY mesangial cells shift bone marrow derived macrophage (BMDM) phenotype to M1 or M2 according to the genetic background (WKY or LEW) of the BMDMs. Interestingly, these effects were different when compared to those of MSCs suggesting that mesangial cells can have unique immunomodulatory effects in the kidney. These results demonstrate the importance of the genetic background in the immunosuppressive effects of cells of stromal origin and specifically of mesangial cell-macrophage interactions in the pathophysiology of crescentic glomerulonephritis

Influence of shear stress magnitude and direction on atherosclerotic plaque composition

British Heart Foundation (BHF) grants (no. RG/11/13/29055 and PG/15/49/31595), awarded to R.K. and E.P. who are employed by Bioengineering and the MRC Clinical Sciences Centre at Imperial College London. The grants also supported R.M.P. and S.M.B. V.V.M. was supported by a BHF. PhD studentship

Inducing persistent flow disturbances accelerates atherogenesis and promotes thin cap fibroatheroma development in D374Y-PCSK9 hypercholesterolemic minipigs

BACKGROUND: -Although disturbed flow is thought to play a central role in the development of advanced coronary atherosclerotic plaques, no causal relationship has been established. We evaluated whether inducing disturbed flow would cause the development of advanced coronary plaques, including thin cap fibroatheroma (TCFA). METHODS AND RESULTS: -D374Y-PCSK9 hypercholesterolemic minipigs (N=5) were instrumented with an intracoronary shear-modifying stent (SMS). Frequency-domain optical coherence tomography was obtained at baseline, immediately post-stent, 19, and 34 weeks and used to compute shear stress metrics of disturbed flow. At 34 weeks, plaque type was assessed within serially-collected histological sections and co-registered to the distribution of each shear metric. The SMS caused a flow-limiting stenosis and blood flow exiting the SMS caused regions of increased shear stress on the outer curvature and large regions of low and multidirectional shear stress on the inner curvature of the vessel. As a result, plaque burden was ~3-fold higher downstream of the SMS compared to both upstream of the SMS and in the control artery (p<0.001). Advanced plaques were also primarily observed downstream of the SMS, in locations initially exposed to both low (p<0.002) and multidirectional (p<0.002) shear stress. TCFA regions demonstrated significantly lower shear stress that persisted over the duration of the study compared to other plaque types (p<0.005). CONCLUSIONS: -These data support a causal role for lowered and multidirectional shear stress in the initiation of advanced coronary atherosclerotic plaques. Persistently lowered shear stress appears to be the principal flow disturbance needed for the formation of TCFA

Intra- and inter-individual genetic differences in gene expression

Genetic variation is known to influence the amount of mRNA produced by a gene. Given that the molecular machines control mRNA levels of multiple genes, we expect genetic variation in the components of these machines would influence multiple genes in a similar fashion. In this study we show that this assumption is correct by using correlation of mRNA levels measured independently in the brain, kidney or liver of multiple, genetically typed, mice strains to detect shared genetic influences. These correlating groups of genes (CGG) have collective properties that account for 40-90% of the variability of their constituent genes and in some cases, but not all, contain genes encoding functionally related proteins. Critically, we show that the genetic influences are essentially tissue specific and consequently the same genetic variations in the one animal may up-regulate a CGG in one tissue but down-regulate the same CGG in a second tissue. We further show similarly paradoxical behaviour of CGGs within the same tissues of different individuals. The implication of this study is that this class of genetic variation can result in complex inter- and intra-individual and tissue differences and that this will create substantial challenges to the investigation of phenotypic outcomes, particularly in humans where multiple tissues are not readily available.&#xd;&#xa;&#xd;&#xa

Investigating the footprint of post-domestication dispersal on the diversity of modern European, African and Asian goats

Background: Goats were domesticated in the Fertile Crescent about 10,000&nbsp;years before present (YBP) and subsequently spread across Eurasia and Africa. This dispersal is expected to generate a gradient of declining genetic diversity with increasing distance from the areas of early livestock management. Previous studies have reported the existence of such genetic cline in European goat populations, but they were based on a limited number of microsatellite markers. Here, we have analyzed data generated by the AdaptMap project and other studies. More specifically, we have used the geographic coordinates and estimates of the observed (Ho) and expected (He) heterozygosities of 1077 European, 1187 African and 617 Asian goats belonging to 38, 43 and 22 different breeds, respectively, to find out whether genetic diversity and distance to Ganj Dareh, a Neolithic settlement in western Iran for which evidence of an early management of domestic goats has been obtained, are significantly correlated. Results: Principal component and ADMIXTURE analyses revealed an incomplete regional differentiation of European breeds, but two genetic clusters representing Northern Europe and the British-Irish Isles were remarkably differentiated from the remaining European populations. In African breeds, we observed five main clusters: (1) North Africa, (2) West Africa, (3) East Africa, (4) South Africa, and (5) Madagascar. Regarding Asian breeds, three well differentiated West Asian, South Asian and East Asian groups were observed. For European and Asian goats, no strong evidence of significant correlations between Ho and He and distance to Ganj Dareh was found. In contrast, in African breeds we detected a significant gradient of diversity, which decreased with distance to Ganj Dareh. Conclusions: The detection of a genetic cline associated with distance to the Ganj Dareh in African but not in European or Asian goat breeds might reflect differences in the post-domestication dispersal process and subsequent migratory movements associated with the management of caprine populations from these three continents

Expression quantitative trait loci are highly sensitive to cellular differentiation state

Blood cell development from multipotent hematopoietic stem cells to specialized blood cells is accompanied by drastic changes in gene expression for which the triggers remain mostly unknown. Genetical genomics is an approach linking natural genetic variation to gene expression variation, thereby allowing the identification of genomic loci containing gene expression modulators (eQTLs). In this paper, we used a genetical genomics approach to analyze gene expression across four developmentally close blood cell types collected from a large number of genetically different but related mouse strains. We found that, while a significant number of eQTLs (365) had a consistent “static” regulatory effect on gene expression, an even larger number were found to be very sensitive to cell stage. As many as 1,283 eQTLs exhibited a “dynamic” behavior across cell types. By looking more closely at these dynamic eQTLs, we show that the sensitivity of eQTLs to cell stage is largely associated with gene expression changes in target genes. These results stress the importance of studying gene expression variation in well-defined cell populations. Only such studies will be able to reveal the important differences in gene regulation between different ce