Expression Levels Of Multiple Genes Research Articles

Long non-coding RNA regulation of gene expression during differentiation.

Transcriptome analysis of mammalian genomes has revealed widespread transcription, much of which does not encode protein. Long non-coding RNAs (lncRNAs) are a subset of the non-coding transcriptome that are emerging as critical regulators of various cellular processes. Differentiation of stem and progenitor cells requires a careful execution of specific genetic programs, and recent studies have revealed that lncRNA expression contributes to specification of cell identity. LncRNAs participate in regulating differentiation at multiple levels of gene expression through various mechanisms of action. In this review, functional roles of lncRNAs in regulating cellular differentiation of blood, muscle, skin, cardiomyocytes, adipocytes, and neurons are discussed.

Tumor-derived exosomes regulate expression of immune function-related genes in human T cell subsets.

Tumor cell-derived exosomes (TEX) suppress functions of immune cells. Here, changes in the gene profiles of primary human T lymphocytes exposed in vitro to exosomes were evaluated. CD4+ Tconv, CD8+ T or CD4+ CD39+ Treg were isolated from normal donors’ peripheral blood and co-incubated with TEX or exosomes isolated from supernatants of cultured dendritic cells (DEX). Expression levels of 24–27 immune response-related genes in these T cells were quantified by qRT-PCR. In activated T cells, TEX and DEX up-regulated mRNA expression levels of multiple genes. Multifactorial data analysis of ΔCt values identified T cell activation and the immune cell type, but not exosome source, as factors regulating gene expression by exosomes. Treg were more sensitive to TEX-mediated effects than other T cell subsets. In Treg, TEX-mediated down-regulation of genes regulating the adenosine pathway translated into high expression of CD39 and increased adenosine production. TEX also induced up-regulation of inhibitory genes in CD4+ Tconv, which translated into a loss of CD69 on their surface and a functional decline. Exosomes are not internalized by T cells, but signals they carry and deliver to cell surface receptors modulate gene expression and functions of human T lymphocytes.

Detection of colorectal cancer genes using a dye-free method combining barcode-base multiplex ligation-dependent probe amplification and pyrosequencing

Early detection of abnormal expression levels of cancer-related genes is crucial for reducing cancer mortality. Here, we describe a dye-free multiplex bioassay for simultaneous and quantitative analysis of the expression levels of multiple genes from one sample in a single assay, based on Sequence-tagged Multiplex ligationdependent probe Amplification (MLPA) coupled with pyrosequencing (termed as “SMAP”). Each pair of MLPA probes, containing a designed barcode, represents a gene of interest; thus, the use of various dyes to label different genes was avoided. The unique three-base barcode design on the probes, which can be decoded by pyrosequencing, enables individual quantification of the expression levels of six genes. Moreover, a new carryover contamination prevention system based on the use of restriction endonucleases was developed for PCR-based diagnostic screening assays. SMAP analysis revealed significant differences between the expression levels of CRC-related genes in the tumor tissues and normal tissues from a CRC patient. For PCR-based diagnostic screening assays, 0.5 U of the FokI restriction endonuclease was sufficient for the removal 0.01 pmol of PCR contamination. The ability to analyze the expression levels of a greater number of cancer-related genes would improve diagnostic sensitivity and efficacy. SMAP is amenable to the detection of an increased number of genes by lengthening the artificially designed barcodes; thus our method provides a promising means for cancer diagnostics and improving the treatment options available to cancer patients.

Assessment of a molecular diagnostic platform for integrated isolation and quantification of mRNA in whole blood.

Implementation of point-of-care tests may facilitate the health management of infectious diseases by reducing the timeframe on pathogen identification and host response measurements, allowing for immediate diagnosis and guided clinical intervention. In this feasibility study, a novel totally integrated and fully automated real-time polymerase chain reaction (PCR) platform (Idylla™, Biocartis) was assessed to determine the mRNA expression levels of multiple genes from 1 mL of whole blood. To this purpose, a sample-in result-out assay, including mRNA extraction and RT-qPCR-based detection, was ported to the platform. The genes used (matrix metallopeptidase 9, olfactomedin 4, NB1 glycoprotein and lipocalin 2) were previously identified as predictive for severity of disease caused by infection with respiratory syncytial virus (RSV). The reproducibility and robustness of the prototype assay was determined using the blood samples of 21 healthy donors. The data showed that the Idylla™ platform allows for a fast and user-friendly determination of the relative expression levels of the four selected mRNA markers.

SNORD116 and SNORD115 change expression of multiple genes and modify each other's activity

The loss of two gene clusters encoding small nucleolar RNAs, SNORD115 and SNORD116 contribute to Prader–Willi syndrome (PWS), the most common syndromic form of obesity in humans. SNORD115 and SNORD116 are considered to be orphan C/D box snoRNAs (SNORDs) as they do not target rRNAs or snRNAs. SNORD115 exhibits sequence complementarity towards the serotonin receptor 2C, but SNORD116 shows no extended complementarities to known RNAs. To identify molecular targets, we performed genome-wide array analysis after overexpressing SNORD115 and SNORD116 in HEK 293T cells, either alone or together. We found that SNORD116 changes the expression of over 200 genes. SNORD116 mainly changed mRNA expression levels. Surprisingly, we found that SNORD115 changes SNORD116's influence on gene expression. In similar experiments, we compared gene expression in post-mortem hypothalamus between individuals with PWS and aged-matched controls. The synopsis of these experiments resulted in 23 genes whose expression levels were influenced by SNORD116. Together our results indicate that SNORD115 and SNORD116 influence expression levels of multiple genes and modify each other activity.

Digital detection of multiple minority mutants and expression levels of multiple colorectal cancer-related genes using digital-PCR coupled with bead-array.

To simultaneously analyze mutations and expression levels of multiple genes on one detection platform, we proposed a method termed “multiplex ligation-dependent probe amplification–digital amplification coupled with hydrogel bead-array” (MLPA–DABA) and applied it to diagnose colorectal cancer (CRC). CRC cells and tissues were sampled to extract nucleic acid, perform MLPA with sequence-tagged probes, perform digital emulsion polymerase chain reaction (PCR), and produce a hydrogel bead-array to immobilize beads and form a single bead layer on the array. After hybridization with fluorescent probes, the number of colored beads, which reflects the abundance of expressed genes and the mutation rate, was counted for diagnosis. Only red or green beads occurred on the chips in the mixed samples, indicating the success of single-molecule PCR. When a one-source sample was analyzed using mixed MLPA probes, beads of only one color occurred, suggesting the high specificity of the method in analyzing CRC mutation and gene expression. In gene expression analysis of a CRC tissue from one CRC patient, the mutant percentage was 3.1%, and the expression levels of CRC-related genes were much higher than those of normal tissue. The highly sensitive MLPA–DABA succeeds in the relative quantification of mutations and gene expressions of exfoliated cells in stool samples of CRC patients on the same chip platform. MLPA–DABA coupled with hydrogel bead-array is a promising method in the non-invasive diagnosis of CRC.

MAP kinase phosphatase-1 expression is regulated by 15-deoxy-Δ12,14-prostaglandin J2 via a HuR-dependent post-transcriptional mechanism

In the present study, we demonstrate a mechanism through which 15-deoxy-Δ12,14-prostaglandin J2 (15d-PGJ2) induces MKP-1 expression in rat primary astrocytes, leading to the regulation of inflammatory responses. We show that 15d-PGJ2 enhances the efficiency of MKP-1 pre-mRNA processing (constitutive splicing and 3′-end processing) and increases the stability of the mature mRNA. We further report that this occurs via the RNA-binding protein, Hu antigen R (HuR). Our experiments show that HuR knockdown abrogates the 15d-PGJ2-induced increases in the pre-mRNA processing and mature mRNA stability of MKP-1, whereas HuR overexpression further enhances the 15d-PGJ2-induced increases in these parameters. Using cysteine (Cys)-mutated HuR proteins, we show that the Cys-245 residue of HuR (but not Cys-13 or Cys-284) is critical for the direct binding of HuR with 15d-PGJ2 and the effects downstream of this interaction. Collectively, our data show that HuR is a novel target of 15d-PGJ2 and reveal HuR-mediated pre-mRNA processing and mature mRNA stabilization as important regulatory steps in the 15d-PGJ2-induced expression of MKP-1. The potential to use a small molecule such as 15d-PGJ2 to regulate the induction of MKP-1 at multiple levels of gene expression could be exploited as a novel therapeutic strategy aimed at combating a diverse range of MKP-1-associated pathologies.

Mapping the genetic architecture of gene regulation in whole blood.

BackgroundWe aimed to assess whether whole blood expression quantitative trait loci (eQTLs) with effects in cis and trans are robust and can be used to identify regulatory pathways affecting disease susceptibility.Materials and MethodsWe performed whole-genome eQTL analyses in 890 participants of the KORA F4 study and in two independent replication samples (SHIP-TREND, N = 976 and EGCUT, N = 842) using linear regression models and Bonferroni correction.ResultsIn the KORA F4 study, 4,116 cis-eQTLs (defined as SNP-probe pairs where the SNP is located within a 500 kb window around the transcription unit) and 94 trans-eQTLs reached genome-wide significance and overall 91% (92% of cis-, 84% of trans-eQTLs) were confirmed in at least one of the two replication studies. Different study designs including distinct laboratory reagents (PAXgene™ vs. Tempus™ tubes) did not affect reproducibility (separate overall replication overlap: 78% and 82%). Immune response pathways were enriched in cis- and trans-eQTLs and significant cis-eQTLs were partly coexistent in other tissues (cross-tissue similarity 40–70%). Furthermore, four chromosomal regions displayed simultaneous impact on multiple gene expression levels in trans, and 746 eQTL-SNPs have been previously reported to have clinical relevance. We demonstrated cross-associations between eQTL-SNPs, gene expression levels in trans, and clinical phenotypes as well as a link between eQTLs and human metabolic traits via modification of gene regulation in cis.ConclusionsOur data suggest that whole blood is a robust tissue for eQTL analysis and may be used both for biomarker studies and to enhance our understanding of molecular mechanisms underlying gene-disease associations.

Use of toxicogenomics to predict the potential toxic effect of Benzo(a)pyrene on zebrafish embryos: Ocular developmental toxicity

Benzo(a)pyrene (BaP) is a representative polycyclic aromatic hydrocarbon (PAH), which is ubiquitous in the environment. The toxic effects of BaP on fish embryos have been described in detail, but some potentially toxic effects of BaP might have been neglected owing to the limitations of traditional techniques. In the present research, global transcriptional patterns were used to study the potentially toxic effects of BaP, as well as its underlying toxicological mechanisms. The expression levels of multiple genes were significantly changed by BaP exposure. The results of ontology assignments and cluster analysis showed that BaP could affect the processes of photoreceptor maintenance and phototransduction. We also conducted an experiment on phototactic response and found that larvae exposed to BaP displayed a decreasing response to light. In addition, BaP exposure decreased the cellular density of the ganglion cell layer (GCL) significantly. These results suggested that BaP exposure induced visual system developmental defects and dysfunction by perturbation of photoreceptor development related genes. Our results were helpful for an understanding of the toxicity of BaP. This study also indicated that microarray analysis was effective for predicting the potential toxicity of chemicals with high sensitivity and accuracy.

An Internal Ribosome Entry Site (IRES) Mutant Library for Tuning Expression Level of Multiple Genes in Mammalian Cells

A set of mutated Encephalomyocarditis virus (EMCV) internal ribosome entry site (IRES) elements with varying strengths is generated by mutating the translation initiation codons of 10th, 11th, and 12th AUG to non-AUG triplets. They are able to control the relative expression of multiple genes over a wide range in mammalian cells in both transient and stable transfections. The relative strength of each IRES mutant remains similar in different mammalian cell lines and is not gene specific. The expressed proteins have correct molecular weights. Optimization of light chain over heavy chain expression by these IRES mutants enhances monoclonal antibody expression level and quality in stable transfections. Uses of this set of IRES mutants can be extended to other applications such as synthetic biology, investigating interactions between proteins and its complexes, cell engineering, multi-subunit protein production, gene therapy, and reprogramming of somatic cells into stem cells.

Cross-Platform Prediction of Gene Expression Signatures

Gene expression signatures can predict the activation of oncogenic pathways and other phenotypes of interest via quantitative models that combine the expression levels of multiple genes. However, as the number of platforms to measure genome-wide gene expression proliferates, there is an increasing need to develop models that can be ported across diverse platforms. Because of the range of technologies that measure gene expression, the resulting signal values can vary greatly. To understand how this variation can affect the prediction of gene expression signatures, we have investigated the ability of gene expression signatures to predict pathway activation across Affymetrix and Illumina microarrays. We hybridized the same RNA samples to both platforms and compared the resultant gene expression readings, as well as the signature predictions. Using a new approach to map probes across platforms, we found that the genes in the signatures from the two platforms were highly similar, and that the predictions they generated were also strongly correlated. This demonstrates that our method can map probes from Affymetrix and Illumina microarrays, and that this mapping can be used to predict gene expression signatures across platforms.

MicroRNA‐214 Protects Cardiac Myocytes Against H2O2‐Induced Injury

Reactive oxygen species (ROS)-induced cardiac myocyte injury resulting from changes in the expression levels of multiple genes plays a critical role in the pathogenesis of numerous heart diseases. The purpose of this study was to determine the potential roles of microRNA-214 (miR-214) in hydrogen peroxide (H2O2)-mediated gene regulation in cardiac myocytes. In this study, we used quantitative real-time RT-PCR (qRT-PCR) to demonstrate that miR-214 was upregulated in cardiac myocytes after treatment with H2O2. We transfected cells with pre-miR-214 to upregulate miR-214 expression and transfected cells with a miR-214 inhibitor (anti-miR-214) to downregulate miR-214 expression. H2O2-induced cardiac cell apoptosis was detected by flow cytometry. The level of apoptosis was increased by the miR-214 inhibitor and decreased by pre-miR-214. Therefore, we believe that miR-214 plays a positive role in H2O2-induced cardiac cell apoptosis. Phosphatase and tensin homolog deleted on chromosome 10 (PTEN) is constitutively active and is considered to be the primary downregulator of the pro-oncogenic PI3K/Akt pathway. Western blot analysis revealed that the expression of the PTEN protein in cardiac myocytes decreased after H2O2 induction. Anti-miR-214 increased PTEN protein expression level, in contrast, pre-miR-214 decreased the PTEN protein expression level in cultured cardiac myocytes. These results indicate that PTEN is regulated by miR-214 and serves as an important target of miR-214 in cardiac myocytes. In conclusion, miR-214 is sensitive to H2O2 stimulation, and miR-214 protects cardiac myocytes against H2O2-induced injury via one of its targets, PTEN.

Positional information, in bits.

Cells in a developing embryo have no direct way of "measuring" their physical position. Through a variety of processes, however, the expression levels of multiple genes come to be correlated with position, and these expression levels thus form a code for "positional information." We show how to measure this information, in bits, using the gap genes in the Drosophila embryo as an example. Individual genes carry nearly two bits of information, twice as much as would be expected if the expression patterns consisted only of on/off domains separated by sharp boundaries. Taken together, four gap genes carry enough information to define a cell's location with an error bar of ~1 along the anterior/posterior axis of the embryo. This precision is nearly enough for each cell to have a unique identity, which is the maximum information the system can use, and is nearly constant along the length of the embryo. We argue that this constancy is a signature of optimality in the transmission of information from primary morphogen inputs to the output of the gap gene network.

Coupling of gene expression in medial prefrontal cortex and nucleus accumbens after neonatal ventral hippocampal lesions accompanies deficits in sensorimotor gating and auditory processing in rats

BackgroundAfter neonatal ventral hippocampal lesions (NVHLs), adult rats exhibit evidence of neural processing deficits relevant to schizophrenia, including reduced prepulse inhibition (PPI) of acoustic startle and impaired sensory processing. In intact rats, the regulation of PPI by the ventral hippocampus (VH) is mediated via interactions with medial prefrontal cortex (mPFC) and nucleus accumbens (NAC). We assessed PPI, auditory-evoked responses and expression of 7 schizophrenia-related genes in mPFC and NAC, in adult rats after sham- or real NVHLs. MethodsMale inbred Buffalo (BUF) rat pups (d7; n = 36) received either vehicle or ibotenic acid infusion into the VH. PPI and auditory-evoked dentate gyrus local field potentials (LFPs) were measured on d56 and d66, respectively. Brains were processed for RT-PCR measures of mPFC and NAC Comt, Erbb4, Grid2, Ncam1, Slc1a2, Nrg1 and Reln. ResultsNVHL rats exhibited significant deficits in PPI (p = 0.005) and LFPs (p < 0.015) proportional to lesion size. Sham vs. NVHL rats did not differ in gene expression levels in mPFC or NAC. As we previously reported, multiple gene expression levels were highly correlated within- (mean r's ≈ 0.5), but not across-brain regions (mean r's ≈ 0). However, for three genes – Comt, Slc1a2 and Ncam1 – after NVHLs, expression levels became significantly correlated, or “coupled,” across the mPFC and NAC (p's < 0.03, 0.002 and 0.05, respectively), and the degree of “coupling” increased with VH lesion size. ConclusionsAfter NVHLs that disrupt PPI and auditory processing, specific gene expression levels suggest an abnormal functional coupling of the mPFC and NAC. This model of VH–mPFC–NAC network dysfunction after NVHLs may have implications for understanding the neural basis for PPI- and related sensory processing deficits in schizophrenia patients.

Insights into the Regular Patterns of Protein Abundance Regulation and Distribution

Protein is one of the fundamental building blocks of life, and its amount in one cell is closely related to biological function. Therefore, protein abundance is regulated strictly at multiple levels of gene expression. The determinants of protein abundance include its corresponding mRNA abundance, the rate of protein synthesis and degradation, so the final cellular protein abundance is determined by the balance of these three factors. Benefitting from the rapid advancement of quantitative proteomics and the outputting of large-scale protein abundance data, we can devote to uncovering the relationship between protein abundance and its intrinsic characteristics (such as evolutionary characteristics, structure characteristics and functional categories), which is of great significance for the deep understanding of the basic composition principle of life systems. Here we summarize the recent research advances of protein abundance regulation, and the regular patterns of the relationship between protein abundance and its intrinsic characteristics. In addition, we also discuss the biological significance underlying them.

Mamld1 Deficiency Significantly Reduces mRNA Expression Levels of Multiple Genes Expressed in Mouse Fetal Leydig Cells but Permits Normal Genital and Reproductive Development

Although mastermind-like domain containing 1 (MAMLD1) (CXORF6) on human chromosome Xq28 has been shown to be a causative gene for 46,XY disorders of sex development with hypospadias, the biological function of MAMLD1/Mamld1 remains to be elucidated. In this study, we first showed gradual and steady increase of testicular Mamld1 mRNA expression levels in wild-type male mice from 12.5 to 18.5 d postcoitum. We then generated Mamld1 knockout (KO) male mice and revealed mildly but significantly reduced testicular mRNA levels (65-80%) of genes exclusively expressed in Leydig cells (Star, Cyp11a1, Cyp17a1, Hsd3b1, and Insl3) as well as grossly normal testicular mRNA levels of genes expressed in other cell types or in Leydig and other cell types. However, no demonstrable abnormality was identified for cytochrome P450 17A1 and 3β-hydroxysteroid dehydrogenase (HSD3B) protein expression levels, appearance of external and internal genitalia, anogenital distance, testis weight, Leydig cell number, intratesticular testosterone and other steroid metabolite concentrations, histological findings, in situ hybridization findings for sonic hedgehog (the key molecule for genital tubercle development), and immunohistochemical findings for anti-Müllerian hormone (Sertoli cell marker), HSD3B (Leydig cell marker), and DEAD (Asp-Glu-Ala-Asp) box polypeptide 4 (germ cell marker) in the KO male mice. Fertility was also normal. These findings imply that Mamld1 deficiency significantly reduces mRNA expression levels of multiple genes expressed in mouse fetal Leydig cells but permits normal genital and reproductive development. The contrastive phenotypic findings between Mamld1 KO male mice and MAMLD1 mutation positive patients would primarily be ascribed to species difference in the fetal sex development.

Identification of genes induced by dibutyl phthalate in human thyroid follicular FTC-238 cells

Dibutyl phthalates (DBP) are widely used as plasticizers to soften and increase the flexibility of polyvinyl chloride plastics, although they can leach into the surrounding environment. Many studies have reported that changes in thyroid function are linked to a variety of chemicals in the environment. Although the mechanism is not known, DBP decreases T3 and T4, and downregulates the sodium/iodide symporter promoter in the thyroid. We undertook a genome-wide analysis of transcriptional responses in human thyroid follicular carcinoma FTC-238 cells to non-toxic doses of DBP (100 μM). Whole human genome microarrays detected a total of 4,480 genes that exhibited altered expression after exposure to non-toxic doses of DBP compared with the vehicle control. We found that 2,218 and 2,262 genes showed different expression levels following treatment with non-toxic doses of DBP (1.5-fold (p <0.01)). Significant changes in the expression levels of multiple genes potentially critical in the regulation of DBP-induced thyroid disease were found. A DBP exposure gene network was identified and established by computational analysis. The possible roles of responsive genes and their related pathways during DBP treatment are discussed. This study details how cells respond to DBP exposure via transcriptional regulation. Information from our network could provide clues for further mechanistic studies that explain how DBP exposure leads to thyroid toxicity and disease.

Requirements for ion and solute transport, and pH regulation during enamel maturation

Transcellular bicarbonate transport is suspected to be an important pathway used by ameloblasts to regulate extracellular pH and support crystal growth during enamel maturation. Proteins that play a role in amelogenesis include members of the ABC transporters (SLC gene family and CFTR). A number of carbonic anhydrases (CAs) have also been identified. The defined functions of these genes are likely interlinked during enamel mineralization. The purpose of this study is to quantify relative mRNA levels of individual SLC, Cftr, and CAs in enamel cells obtained from secretory and maturation stages on rat incisors. We also present novel data on the enamel phenotypes for two animal models, a mutant porcine (CFTR-ΔF508) and the NBCe1-null mouse. Our data show that two SLCs (AE2 and NBCe1), Cftr, and Car2, Car3, Car6, and Car12 are all significantly up-regulated at the onset of the maturation stage of amelogenesis when compared to the secretory stage. The remaining SLCs and CA gene transcripts showed negligible expression or no significant change in expression from secretory to maturation stages. The enamel of CFTR-ΔF508 adult pigs was hypomineralized and showed abnormal crystal growth. NBCe1-null mice enamel was structurally defective and had a marked decrease in mineral content relative to wild-type. These data demonstrate the importance of many non-matrix proteins to amelogenesis and that the expression levels of multiple genes regulating extracellular pH are modulated during enamel maturation in response to an increased need for pH buffering during hydroxyapatite crystal growth.

Benzo[a]pyrene exposure influences the cardiac development and the expression of cardiovascular relative genes in zebrafish (Daniorerio) embryos

It is reported that the most abundant polycyclic aromatic hydrocarbons (PAHs) in weathered crude oils are cardiotoxic. However, the action mechanism of PAHs on vertebrate cardiovascular development and disease is unclear. In the present study, the cardiac morphology and functioning of zebrafish embryos exposed to benzo[a]pyrene [B(a)P], as a high-ring PAHs, for 72h were observed and determined. The results showed that B(a)P exposure resulted in cardiac developmental defects in zebrafish embryos. Significant changes in expression level of multiple genes potentially critical for regulating the B(a)P-induced cardiovascular developmental defects were also found. A gene network regulating cardiac development perturbed by B(a)P exposure was identified and established by computational analysis and employment of some databases. The information from the network could provide a clue for further mechanistic studies explaining molecular events regulating B(a)P-mediated cardiovascular defects and consequences.

Genetical genomics of Th1 and Th2 immune response in a baboon model of atherosclerosis risk factors

Objective CD4 + T-cells mediate inflammation in atherosclerosis, but additive genetic effects on associated pathways of Th1 and Th2 immune response have not been described. We sought to characterize heritability, pleiotropy, and QTL effects on the expression of genes implicated in Th1 and Th2 immune response in a baboon model of risk factors for atherosclerosis. Methods We employed a maximum likelihood-based variance decomposition approach to estimate additive genetic effects on transcript levels generated from a gene expression profile of lymphocytes in 499 pedigreed baboons maintained on a basal diet. Transcript levels for 57 genes implicated in Th1 and Th2 immune response were selected for analysis based on significant heritability in this profile. Multipoint whole genome scans were conducted on heritable transcript levels to localize QTLs influencing these measures. To evaluate pleiotropic effects on transcript levels, we estimated genetic and phenotypic correlations among transcript measures, and assessed their correspondence using a Mantel test. Network analysis using GeneGo's MetaCore™ software was conducted to characterize known interaction among coded proteins. Results Heritabilities for candidate gene transcript levels ranged from 0.092–0.786 (median h 2 = 0.278, P = 4.72 × 10 −4). Linkage analyses yielded significant evidence (LOD ≥ 2.73) for 14 eQTLs (LOD score range 2.76–14.87, genome-wide P = 4.9 × 10 −2–1.03 × 10 −14). Estimates of genetic correlation supported shared additive genetic effects incorporating all 57 transcripts (null hypothesis of ρ G = 0 rejected at FDR ≤ 0.05 for 522 of 1596 estimates), and accounted for most of the observed phenotypic correlation among transcripts (Mantel test, r [ ρ P], [ ρ G] = 0.781, P < 0.0001). Network analysis revealed direct interactions among 54 of the 57 coded proteins. Conclusions We conclude that major genetic effects influence expression levels of multiple genes implicated in Th1 and Th2 immune response. Additionally, we find that expression levels of these candidate genes are characterized by extensive pleiotropy, consistent with known interaction among their coded proteins, many of which are independently associated with atherosclerosis.