Carotid plaque age is a feature of plaque stability inversely related to levels of plasma insulin.

BACKGROUND: The stability of atherosclerotic plaques determines the risk for rupture, which may lead to thrombus formation and potentially severe clinical complications such as myocardial infarction and stroke. Although the rate of plaque formation may be important for plaque stability, this process is not well understood. We took advantage of the atmospheric (14)C-declination curve (a result of the atomic bomb tests in the 1950s and 1960s) to determine the average biological age of carotid plaques. METHODOLOGY/PRINCIPAL FINDING: The cores of carotid plaques were dissected from 29 well-characterized, symptomatic patients with carotid stenosis and analyzed for (14)C content by accelerator mass spectrometry. The average plaque age (i.e. formation time) was 9.6±3.3 years. All but two plaques had formed within 5-15 years before surgery. Plaque age was not associated with the chronological ages of the patients but was inversely related to plasma insulin levels (p = 0.0014). Most plaques were echo-lucent rather than echo-rich (2.24±0.97, range 1-5). However, plaques in the lowest tercile of plaque age (most recently formed) were characterized by further instability with a higher content of lipids and macrophages (67.8±12.4 vs. 50.4±6.2, p = 0.00005; 57.6±26.1 vs. 39.8±25.7, p<0.0005, respectively), less collagen (45.3±6.1 vs. 51.1±9.8, p<0.05), and fewer smooth muscle cells (130±31 vs. 141±21, p<0.05) than plaques in the highest tercile. Microarray analysis of plaques in the lowest tercile also showed increased activity of genes involved in immune responses and oxidative phosphorylation. CONCLUSIONS/SIGNIFICANCE: Our results show, for the first time, that plaque age, as judge by relative incorporation of (14)C, can improve our understanding of carotid plaque stability and therefore risk for clinical complications. Our results also suggest that levels of plasma insulin might be involved in determining carotid plaque age.

Multi-organ expression profiling uncovers a gene module in coronary artery disease involving transendothelial migration of leukocytes and LIM domain binding 2: the Stockholm Atherosclerosis Gene Expression (STAGE) study.

Environmental exposures filtered through the genetic make-up of each individual alter the transcriptional repertoire in organs central to metabolic homeostasis, thereby affecting arterial lipid accumulation, inflammation, and the development of coronary artery disease (CAD). The primary aim of the Stockholm Atherosclerosis Gene Expression (STAGE) study was to determine whether there are functionally associated genes (rather than individual genes) important for CAD development. To this end, two-way clustering was used on 278 transcriptional profiles of liver, skeletal muscle, and visceral fat (n = 66/tissue) and atherosclerotic and unaffected arterial wall (n = 40/tissue) isolated from CAD patients during coronary artery bypass surgery. The first step, across all mRNA signals (n = 15,042/12,621 RefSeqs/genes) in each tissue, resulted in a total of 60 tissue clusters (n = 3958 genes). In the second step (performed within tissue clusters), one atherosclerotic lesion (n = 49/48) and one visceral fat (n = 59) cluster segregated the patients into two groups that differed in the extent of coronary stenosis (P = 0.008 and P = 0.00015). The associations of these clusters with coronary atherosclerosis were validated by analyzing carotid atherosclerosis expression profiles. Remarkably, in one cluster (n = 55/54) relating to carotid stenosis (P = 0.04), 27 genes in the two clusters relating to coronary stenosis were confirmed (n = 16/17, P<10(-27 and-30)). Genes in the transendothelial migration of leukocytes (TEML) pathway were overrepresented in all three clusters, referred to as the atherosclerosis module (A-module). In a second validation step, using three independent cohorts, the A-module was found to be genetically enriched with CAD risk by 1.8-fold (P<0.004). The transcription co-factor LIM domain binding 2 (LDB2) was identified as a potential high-hierarchy regulator of the A-module, a notion supported by subnetwork analysis, by cellular and lesion expression of LDB2, and by the expression of 13 TEML genes in Ldb2-deficient arterial wall. Thus, the A-module appears to be important for atherosclerosis development and, together with LDB2, merits further attention in CAD research.

Reverse engineering of gene networks with LASSO and nonlinear basis functions.

The quest to determine cause from effect is often referred to as reverse engineering in the context of cellular networks. Here we propose and evaluate an algorithm for reverse engineering a gene regulatory network from time-series and steady-state data. Our algorithmic pipeline, which is rather standard in its parts but not in its integrative composition, combines ordinary differential equations, parameter estimations by least angle regression, and cross-validation procedures for determining the in-degrees and selection of nonlinear transfer functions. The result of the algorithm is a complete directed network, in which each edge has been assigned a score from a bootstrap procedure. To evaluate the performance, we submitted the outcome of the algorithm to the reverse engineering assessment competition DREAM2, where we used the data corresponding to the InSilico1 and InSilico2 networks as input. Our algorithm outperformed all other algorithms when inferring one of the directed gene-to-gene networks.

Transcriptional profiling uncovers a network of cholesterol-responsive atherosclerosis target genes.

Despite the well-documented effects of plasma lipid lowering regimes halting atherosclerosis lesion development and reducing morbidity and mortality of coronary artery disease and stroke, the transcriptional response in the atherosclerotic lesion mediating these beneficial effects has not yet been carefully investigated. We performed transcriptional profiling at 10-week intervals in atherosclerosis-prone mice with human-like hypercholesterolemia and a genetic switch to lower plasma lipoproteins (Ldlr(-/-)Apo(100/100)Mttp(flox/flox) Mx1-Cre). Atherosclerotic lesions progressed slowly at first, then expanded rapidly, and plateaued after advanced lesions formed. Analysis of lesion expression profiles indicated that accumulation of lipid-poor macrophages reached a point that led to the rapid expansion phase with accelerated foam-cell formation and inflammation, an interpretation supported by lesion histology. Genetic lowering of plasma cholesterol (e.g., lipoproteins) at this point all together prevented the formation of advanced plaques and parallel transcriptional profiling of the atherosclerotic arterial wall identified 37 cholesterol-responsive genes mediating this effect. Validation by siRNA-inhibition in macrophages incubated with acetylated-LDL revealed a network of eight cholesterol-responsive atherosclerosis genes regulating cholesterol-ester accumulation. Taken together, we have identified a network of atherosclerosis genes that in response to plasma cholesterol-lowering prevents the formation of advanced plaques. This network should be of interest for the development of novel atherosclerosis therapies.

Evidence of highly regulated genes (in-Hubs) in gene networks of Saccharomyces cerevisiae.

Uncovering interactions between genes, gene networks, is important to increase our understanding of intrinsic cellular processes and responses to external stimuli such as drugs. Gene networks can be computationally inferred from repeated measurements of gene expression, using algorithms, which assume that each gene is controlled by only a small number of other proteins. Here, by extending the transcription network with cofactors (defined from protein-protein binding data) as active regulators, we identified the effective gene network, providing evidence of in-hubs in the gene regulatory networks of yeast. Then, using the notion that in-hub genes will be differentially expressed over several experimental conditions, we designed an algorithm, the HubDetector, enabling identification of in-hubs directly from gene expression data. Applying the HubDetector to 488 genome-wide expression profiles from two independent datasets, we identified putative in-hubs overlapping significantly with in-hubs in the effective gene network.