Microarray and comparative genomics-based identification of genes and gene regulatory regions of the mouse immune system

Bioinformatics analysis and validation of CSRNP1 as a key prognostic gene in non-small cell lung cancer

Cysteine and serine-rich nuclear protein 1 (CSRNP1) has shown prognostic significance in various cancers, but its role in non-small cell lung cancer (NSCLC) remains elusive. We investigated CSRNP1 expression in NSCLC cases using bioinformatics tools from the GEO public repository and validated our findings through RT-qPCR in tumor and adjacent normal tissues. KEGG and GO enrichment analyses were employed to unveil the significant deregulation in signaling pathways. Additionally, clinical significance of CSRNP1 in NSCLC was determined through receiver operating curve (ROC) analysis, and its impact on survival was assessed using Kaplan-Meier analysis. To explore the functional impact of CSRNP1, we silenced its expression in NSCLC cells and assessed the effects on cell viability, migration, and invasion using MTT, Transwell, and wound-healing assays, respectively. Additionally, we investigated the influence of CSRNP1 silencing on the phosphorylation patterns of critical signaling proteins such as p53, p-Akt, and p-MDM2. Our results demonstrated significantly lower CSRNP1 expression in NSCLC tumor tissues (P < 0.01). ROC analysis indicated that NSCLC patients with high CSRNP1 expression exhibited extended overall survival and disease-free survival. Furthermore, CSRNP1 silencing promoted NSCLC cells viability, migration, and invasion (P < 0.05). Mechanistically, CSRNP1 silencing led to increased phosphorylation of AKT and MDM2, along with a concurrent reduction in p53 protein expression, suggesting its impact on NSCLC through deregulated cell cycle processes. In conclusion, our study underscores the significance of CSRNP1 in NSCLC pathogenesis, offering insights for targeted therapeutic interventions of NSCLC.

The CSRNP Gene Family Serves as a Prognostic Biomarker in Clear Cell Renal Cell Carcinoma

The cysteine-serine-rich nuclear protein (CSRNP) family has prognostic value for various cancers. However, the association between this proteins and prognosis of clear cell renal cell carcinoma (ccRCC) remains unclear. This study aimed to determine the prognostic value of the CSRNP family for patients with ccRCC. Therefore, the gene expression profiling interactive analysis database was used to analyze the mRNA expression of CSRNP family members (CSRNPs) in relation with survival. Combined and independent prognostic values of CSRNPs were evaluated using SurvExpress and multivariate Cox regression analyses, respectively. Potential signaling pathways impacted by CSRNPs were evaluated using Metascape. Associations between the CSRNP family and immunocyte infiltration were determined from single-sample gene set enrichment analysis. Both cBioPortal and MethSurv were used to explore whether genomic and epidemic alterations might influence prognosis. We found that when both CSRNP1 and CSRNP3 had a low expression, patients with ccRCC had a worse overall survival (OS). Therefore, a prognostic signature was constructed as follows: risk score = -0.224 × exp_{mRNA of CSRNP1} + 0.820 × exp_{mRNA of CSRNP2} - 1.428 × exp_{mRNA of CSRNP3} . We found that OS was worse in patients from the high- than from the low-risk groups (AUC = 0.69). Moreover, this signature was an independent predictor after adjusting for clinical features. Functional enrichment analysis positively associated CSRNPs with the acute inflammatory response and humoral immune response pathways. This was validated by correlating each CSRNP with 28 types of immunocytes in tumor and normal tissues. A higher expression of CSRNP1 and CSRNP3 was associated with a better prognosis in both the high- and low-mutant burden groups. Cg19538674, cg07772537, and cg07811002 of CSRNP1, CSRNP2, and CSRNP3, respectively, were the predominant DNA methylation sites affecting OS. The CSRNP gene family signature may serve as a prognostic biomarker for predicting OS in patients with ccRCC. The association between CSRNPs and immune infiltration might offer future clinical treatment options.

Restoration of GATA4 expression impedes breast cancer progression by transcriptional repression of ReLA and inhibition of NF-κB signaling

There are increasing reports of aberrant expression of GATA4, correlated with oncogenesis and malignant progression in some solid tumors, but whether GATA4 functions as an oncogenic driver or a tumor suppressor in carcinogenesis remains controversial. Because the role and mechanism of GATA4 in breast cancer (BrCa) remain poorly understood, we focused on the expression of GATA4 in BrCa cell lines and tissues and its mechanism in breast oncogenesis. Semiquantitative real-time polymerase chain reaction (RT-PCR), quantitative RT-PCR, Western blot analysis, and immunohistochemistry were used to detect expression of GATA4 in BrCa cell lines and adjacent breast tissues. Methylation statuses of the GATA4 promoter were studied using methylation-specific PCR in BrCa cell lines.The effects of GATA4 on proliferation, invasion, and cell cycle were also analyzed. Compared with adjacent breast tissue, GATA4 expression in BrCa tissue and cell lines was obviously lower and low expression levels of GATA4 predicted poor outcome. Methylation of GATA4 occurred in almost all of BrCa cell lines . GATA4 overexpression decreased viability, invasion, migration, and epithelial-to-mesenchymal transition of MB-231 and BT549 cells, and markedly induced cell cycle arrest and apoptosis. Exogenous expression GATA4 accompanied a significant alteration of MMP2, MMP3, E-cadherin, and N-cadherin expression and induction of the caspase-8 pathway. Moreover, GATA4 could directly repress RelA (p65) transcription, reduce the nuclear phosphorylation-p65 and upregulate inhibitor kappa B expression. Altogether, GATA4 plays a tumor-suppressive role via repression of NF-κB signaling in BrCa cells. Our findings suggest that GATA4 is a potential prognostic biomarker and gene therapeutic target for human BrCa.

Screening for genes involved in antibody response to sheep red blood cells in the chicken, Gallus gallus

Antibody response, an important trait in both agriculture and biomedicine, plays a part in protecting animals from infection. Dissecting molecular basis of antibody response may improve artificial selection for natural disease resistance in livestock and poultry. A number of genetic markers associated with antibody response have been identified in the chicken and mouse by linkage-based association studies, which only define genomic regions by genetic markers but do not pinpoint genes for antibody response. In contrast, global expression profiling has been applied to define the molecular bases of a variety of biological traits through identification of differentially expressed genes (DEGs). Here, we employed Affimetrix GeneChip Chicken Genome Arrays to identify differentially expressed genes for antibody response to sheep red blood cells (SRBC) using chickens challenged with and without SRBC or chickens with high and low anti-SRBC titers. The DEGs include those with known (i.e., MHC class I and IgH genes) or unknown function in antibody response. Classification test of these genes suggested that the response of the chicken to intravenous injection of SRBC involved multiple biological processes, including response to stress or other different stimuli, sugar, carbohydrate or protein binding, and cell or soluble fraction, in addition to antibody response. This preliminary study thus provides an insight into molecular basis of antibody response to SRBC in the chicken.

Molecular bases of catalysis and ADP-ribose preference of human Mn2+-dependent ADP-ribose/CDP-alcohol diphosphatase and conversion by mutagenesis to a preferential cyclic ADP-ribose phosphohydrolase

Among metallo-dependent phosphatases, ADP-ribose/CDP-alcohol diphosphatases form a protein family (ADPRibase-Mn-like) mainly restricted, in eukaryotes, to vertebrates and plants, with preferential expression, at least in rodents, in immune cells. Rat and zebrafish ADPRibase-Mn, the only biochemically studied, are phosphohydrolases of ADP-ribose and, somewhat less efficiently, of CDP-alcohols and 2´,3´-cAMP. Furthermore, the rat but not the zebrafish enzyme displays a unique phosphohydrolytic activity on cyclic ADP-ribose. The molecular basis of such specificity is unknown. Human ADPRibase-Mn showed similar activities, including cyclic ADP-ribose phosphohydrolase, which seems thus common to mammalian ADPRibase-Mn. Substrate docking on a homology model of human ADPRibase-Mn suggested possible interactions of ADP-ribose with seven residues located, with one exception (Cys253), either within the metallo-dependent phosphatases signature (Gln27, Asn110, His111), or in unique structural regions of the ADPRibase-Mn family: s2s3 (Phe37 and Arg43) and h7h8 (Phe210), around the active site entrance. Mutants were constructed, and kinetic parameters for ADP-ribose, CDP-choline, 2´,3´-cAMP and cyclic ADP-ribose were determined. Phe37 was needed for ADP-ribose preference without catalytic effect, as indicated by the increased ADP-ribose Km and unchanged kcat of F37A-ADPRibase-Mn, while the Km values for the other substrates were little affected. Arg43 was essential for catalysis as indicated by the drastic efficiency loss shown by R43A-ADPRibase-Mn. Unexpectedly, Cys253 was hindering for cADPR phosphohydrolase, as indicated by the specific tenfold gain of efficiency of C253A-ADPRibase-Mn with cyclic ADP-ribose. This allowed the design of a triple mutant (F37A+L196F+C253A) for which cyclic ADP-ribose was the best substrate, with a catalytic efficiency of 3.5´104 M-1s-1 versus 4´103 M-1s-1 of the wild type.

No more non-model species: the promise of next generation sequencing for comparative immunology

Next generation sequencing (NGS) allows for the rapid, comprehensive and cost effective analysis of entire genomes and transcriptomes. NGS provides approaches for immune response gene discovery, profiling gene expression over the course of parasitosis, studying mechanisms of diversification of immune receptors and investigating the role of epigenetic mechanisms in regulating immune gene expression and/or diversification. NGS will allow meaningful comparisons to be made between organisms from different taxa in an effort to understand the selection of diverse strategies for host defence under different environmental pathogen pressures. At the same time, it will reveal the shared and unique components of the immunological toolkit and basic functional aspects that are essential for immune defence throughout the living world. In this review, we argue that NGS will revolutionize our understanding of immune responses throughout the animal kingdom because the depth of information it provides will circumvent the need to concentrate on a few "model" species.

Gene responses in the central nervous system of zebrafish embryos exposed to the neurotoxicant methyl mercury

Methyl mercury (MeHg) is a neurotoxicant with adverse effects on the development of the nervous system from fish to man. Despite a detailed understanding of the molecular mechanisms by which MeHg affects cellular homeostasis, it is still not clear how MeHg causes developmental neurotoxicity. We performed here a genome-wide transcriptional analysis of MeHg-exposed zebrafish embryos and combined this with a whole-mount in situ expression analysis of 88 MeHg-affected genes. The majority of the analyzed genes showed tissue- and region-restricted responses in various organs and tissues. The genes were linked to gene ontology terms like oxidative stress, transport and cell protection. Areas even within the central nervous system (CNS) are affected differently resulting in distinct cellular stress responses. Our study revealed an unexpected heterogeneity in gene responses to MeHg exposure in different tissues and neuronal subregions, even though the known molecular action of MeHg would predict a similar burden of exposed cells. The overall structure of the developing brain of MeHg-exposed embryos appeared normal, suggesting that the mechanism leading to differentiation of the CNS is not overtly affected by exposure to MeHg. We propose that MeHg disturbs the function of the CNS by disturbing the cellular homeostasis. As these cellular stress responses comprise genes that are also involved in normal neuronal activity and learning, MeHg may affect the developing CNS in a subtle manner that manifests itself in behavioral deficits.

Analysis of the membrane proteome of ciprofloxacin-resistant macrophages by stable isotope labeling with amino acids in cell culture (SILAC)

Overexpression of multidrug transporters is a well-established mechanism of resistance to chemotherapy, but other changes may be co-selected upon exposure to drugs that contribute to resistance. Using a model of J774 macrophages made resistant to the fluoroquinolone antibiotic ciprofloxacin and comparing it with the wild-type parent cell line, we performed a quantitative proteomic analysis using the stable isotope labeling with amino acids in cell culture technology coupled with liquid chromatography electrospray ionization Fourier transform tandem mass spectrometry (LC-ESI-FT-MS/MS) on 2 samples enriched in membrane proteins (fractions F1 and F2 collected from discontinuous sucrose gradient). Nine hundred proteins were identified with at least 3 unique peptides in these 2 pooled fractions among which 61 (F1) and 69 (F2) showed a significantly modified abundance among the 2 cell lines. The multidrug resistance associated protein Abcc4, known as the ciprofloxacin efflux transporter in these cells, was the most upregulated, together with Dnajc3, a protein encoded by a gene located downstream of Abcc4. The other modulated proteins are involved in transport functions, cell adhesion and cytoskeleton organization, immune response, signal transduction, and metabolism. This indicates that the antibiotic ciprofloxacin is able to trigger a pleiotropic adaptative response in macrophages that includes the overexpression of its efflux transporter.

An integrative computational systems biology approach identifies differentially regulated dynamic transcriptome signatures which drive the initiation of human T helper cell differentiation

Background

A proper balance between different T helper (Th) cell subsets is necessary for normal functioning of the adaptive immune system. Revealing key genes and pathways driving the differentiation to distinct Th cell lineages provides important insight into underlying molecular mechanisms and new opportunities for modulating the immune response. Previous computational methods to quantify and visualize kinetic differential expression data of three or more lineages to identify reciprocally regulated genes have relied on clustering approaches and regression methods which have time as a factor, but have lacked methods which explicitly model temporal behavior.

Results

We studied transcriptional dynamics of human umbilical cord blood T helper cells cultured in absence and presence of cytokines promoting Th1 or Th2 differentiation. To identify genes that exhibit distinct lineage commitment dynamics and are specific for initiating differentiation to different Th cell subsets, we developed a novel computational methodology (LIGAP) allowing integrative analysis and visualization of multiple lineages over whole time-course profiles. Applying LIGAP to time-course data from multiple Th cell lineages, we identified and experimentally validated several differentially regulated Th cell subset specific genes as well as reciprocally regulated genes. Combining differentially regulated transcriptional profiles with transcription factor binding site and pathway information, we identified previously known and new putative transcriptional mechanisms involved in Th cell subset differentiation. All differentially regulated genes among the lineages together with an implementation of LIGAP are provided as an open-source resource.

Conclusions

The LIGAP method is widely applicable to quantify differential time-course dynamics of many types of datasets and generalizes to any number of conditions. It summarizes all the time-course measurements together with the associated uncertainty for visualization and manual assessment purposes. Here we identified novel human Th subset specific transcripts as well as regulatory mechanisms important for the initiation of the Th cell subset differentiation.

Characterization of bovine FUT7 furthers understanding of FUT7 evolution in mammals

Background

The Sialyl-Lewis X (Sle^x) is a well-known glycan structure involved in leukocyte homing and recruitment to inflammatory sites. SLe^x is well conserved among species and is mainly synthesized by FucT-VII in vertebrates. The enzyme responsible for its biosynthesis in cattle was not known.

Results

We cloned a cDNA sequence encoding bovine α3-fucosyltransferase VII that shares 83% identity with its human counterpart. Located at the BTA 11 telomeric region, the 1029 bp open reading frame is spread over two different exons, E1 which also contains the unique 5’-untranslated region and E2 which includes the entire 3’-untranslated region. The bfut7 expression pattern is restricted to thymus and spleen. A single transcript leading to the synthesis of a 342 aa protein was identified. The encoded fucosyltransferase, produced as a recombinant enzyme in COS-1 cells, was shown to be specifically responsible for SLe^x synthesis in cattle. In addition, we showed that the gene promoter evolved from fish to mammals towards a complex system related to the immune system. But beyond the fact that the gene regulation seems to be conserved among mammals, we also identified 7 SNPs including 3 missense mutations in the coding region in a small panel of animals.

Conclusions

The FUT7 sequence was highly conserved as well as the specific activity of the encoded protein FucT-VII. In addition, our in silico promoter analysis and the high rate of polymorphism suggested that its function is evolving toward a complex system related to the immune system. Furthermore, comparing bovine to human and mouse sequences, it appeared that a decrease in gene regulation was correlated with an increase in mutation rate and wider tissue expression.

Characterization of Danio rerio Mn2+-dependent ADP-ribose/CDP-alcohol diphosphatase, the structural prototype of the ADPRibase-Mn-like protein family

The ADPRibase-Mn-like protein family, that belongs to the metallo-dependent phosphatase superfamily, has different functional and structural prototypes. The functional one is the Mn²⁺-dependent ADP-ribose/CDP-alcohol diphosphatase from Rattus norvegicus, which is essentially inactive with Mg²⁺ and active with low micromolar Mn²⁺ in the hydrolysis of the phosphoanhydride linkages of ADP-ribose, CDP-alcohols and cyclic ADP-ribose (cADPR) in order of decreasing efficiency. The structural prototype of the family is a Danio rerio protein with a known crystallographic structure but functionally uncharacterized. To estimate the structure-function correlation with the same protein, the activities of zebrafish ADPRibase-Mn were studied. Differences between zebrafish and rat enzymes are highlighted. The former showed a complex activity dependence on Mn²⁺, significant (≈25%) Mg²⁺-dependent activity, but was almost inactive on cADPR (150-fold less efficient than the rat counterpart). The low cADPR hydrolase activity agreed with the zebrafish genome lacking genes coding for proteins with significant homology with cADPR-forming enzymes. Substrate-docking to zebrafish wild-type protein, and characterization of the ADPRibase-Mn H97A mutant pointed to a role of His-97 in catalysis by orientation, and to a bidentate water bridging the dinuclear metal center as the potential nucleophile. Finally, three structural elements that delimit the active site entrance in the zebrafish protein were identified as unique to the ADPRibase-Mn-like family within the metallo-dependent phosphatase superfamily.

Reorganization of the host epigenome by a viral oncogene

Adenovirus small e1a oncoprotein causes ∼70% reduction in cellular levels of histone H3 lysine 18 acetylation (H3K18ac). It is unclear, however, where this dramatic reduction occurs genome-wide. ChIP-sequencing revealed that by 24 h after expression, e1a erases 95% of H3K18ac peaks in normal, contact-inhibited fibroblasts and replaces them with one-third as many at new genomic locations. The H3K18ac peaks at promoters and intergenic regions of genes with fibroblast-related functions are eliminated after infection, and new H3K18ac peaks are established at promoters of highly induced genes that regulate cell cycling and at new putative enhancers. Strikingly, the regions bound by the retinoblastoma family of proteins in contact-inhibited fibroblasts gain new peaks of H3K18ac in the e1a-expressing cells, including 55% of RB1-bound loci. In contrast, over half of H3K9ac peaks are similarly distributed before and after infection, independently of RB1. The strategic redistribution of H3K18ac by e1a highlights the importance of this modification for transcriptional activation and cellular transformation as well as functional differences between the RB-family member proteins.

Isolation and preparation of RNA from rat blood and lymphatic microvessels for use in microarray analysis

DNA microarray methodologies have proven to be an indispensable tool for genome-wide transcriptional profiling of organs, tissues, and cells. Here, we present a protocol for the optimized isolation and preparation of RNA from rat microvessels (including arteries, veins, and lymphatics) for subsequent use in two-color microarray analysis. The investigation of wide-ranging vessel sizes from all three vessel lineages necessitates an RNA isolation strategy that can effectively isolate high-quality RNA from varying and often very small quantities (<1 mg) of fibrous vessel tissue. Additionally, the lack of sample biomass necessitates the use of amplification strategies to generate enough RNA for use in microarray analysis. While the methods presented here were developed for use with two-color microarray analysis, the procedures and general concepts are applicable to most fluorescence-based microarray platforms.

Computational identification of transcriptional regulators in human endotoxemia

One of the great challenges in the post-genomic era is to decipher the underlying principles governing the dynamics of biological responses. As modulating gene expression levels is among the key regulatory responses of an organism to changes in its environment, identifying biologically relevant transcriptional regulators and their putative regulatory interactions with target genes is an essential step towards studying the complex dynamics of transcriptional regulation. We present an analysis that integrates various computational and biological aspects to explore the transcriptional regulation of systemic inflammatory responses through a human endotoxemia model. Given a high-dimensional transcriptional profiling dataset from human blood leukocytes, an elementary set of temporal dynamic responses which capture the essence of a pro-inflammatory phase, a counter-regulatory response and a dysregulation in leukocyte bioenergetics has been extracted. Upon identification of these expression patterns, fourteen inflammation-specific gene batteries that represent groups of hypothetically ‘coregulated’ genes are proposed. Subsequently, statistically significant cis-regulatory modules (CRMs) are identified and decomposed into a list of critical transcription factors (34) that are validated largely on primary literature. Finally, our analysis further allows for the construction of a dynamic representation of the temporal transcriptional regulatory program across the host, deciphering possible combinatorial interactions among factors under which they might be active. Although much remains to be explored, this study has computationally identified key transcription factors and proposed a putative time-dependent transcriptional regulatory program associated with critical transcriptional inflammatory responses. These results provide a solid foundation for future investigations to elucidate the underlying transcriptional regulatory mechanisms under the host inflammatory response. Also, the assumption that coexpressed genes that are functionally relevant are more likely to share some common transcriptional regulatory mechanism seems to be promising, making the proposed framework become essential in unravelling context-specific transcriptional regulatory interactions underlying diverse mammalian biological processes.

Transcriptomes of the B and T lineages compared by multiplatform microarray profiling

T and B lymphocytes are developmentally and functionally related cells of the immune system, representing the two major branches of adaptive immunity. Although originating from a common precursor, they play very different roles: T cells contribute to and drive cell-mediated immunity, whereas B cells secrete Abs. Because of their functional importance and well-characterized differentiation pathways, T and B lymphocytes are ideal cell types with which to understand how functional differences are encoded at the transcriptional level. Although there has been a great deal of interest in defining regulatory factors that distinguish T and B cells, a truly genomewide view of the transcriptional differences between these two cells types has not yet been taken. To obtain a more global perspective of the transcriptional differences underlying T and B cells, we exploited the statistical power of combinatorial profiling on different microarray platforms, and the breadth of the Immunological Genome Project gene expression database, to generate robust differential signatures. We find that differential expression in T and B cells is pervasive, with the majority of transcripts showing statistically significant differences. These distinguishing characteristics are acquired gradually, through all stages of B and T differentiation. In contrast, very few T versus B signature genes are uniquely expressed in these lineages, but are shared throughout immune cells.

Comparative analysis of acute and chronic corticosteroid pharmacogenomic effects in rat liver: transcriptional dynamics and regulatory structures

Background

Comprehensively understanding corticosteroid pharmacogenomic effects is an essential step towards an insight into the underlying molecular mechanisms for both beneficial and detrimental clinical effects. Nevertheless, even in a single tissue different methods of corticosteroid administration can induce different patterns of expression and regulatory control structures. Therefore, rich in vivo datasets of pharmacological time-series with two dosing regimens sampled from rat liver are examined for temporal patterns of changes in gene expression and their regulatory commonalities.

Results

The study addresses two issues, including (1) identifying significant transcriptional modules coupled with dynamic expression patterns and (2) predicting relevant common transcriptional controls to better understand the underlying mechanisms of corticosteroid adverse effects. Following the orientation of meta-analysis, an extended computational approach that explores the concept of agreement matrix from consensus clustering has been proposed with the aims of identifying gene clusters that share common expression patterns across multiple dosing regimens as well as handling challenges in the analysis of microarray data from heterogeneous sources, e.g. different platforms and time-grids in this study. Six significant transcriptional modules coupled with typical patterns of expression have been identified. Functional analysis reveals that virtually all enriched functions (gene ontologies, pathways) in these modules are shown to be related to metabolic processes, implying the importance of these modules in adverse effects under the administration of corticosteroids. Relevant putative transcriptional regulators (e.g. RXRF, FKHD, SP1F) are also predicted to provide another source of information towards better understanding the complexities of expression patterns and the underlying regulatory mechanisms of those modules.

Conclusions

We have proposed a framework to identify significant coexpressed clusters of genes across multiple conditions experimented from different microarray platforms, time-grids, and also tissues if applicable. Analysis on rich in vivo datasets of corticosteroid time-series yielded significant insights into the pharmacogenomic effects of corticosteroids, especially the relevance to metabolic side-effects. This has been illustrated through enriched metabolic functions in those transcriptional modules and the presence of GRE binding motifs in those enriched pathways, providing significant modules for further analysis on pharmacogenomic corticosteroid effects.

Gene expression profiles identify inflammatory signatures in dendritic cells

Dendritic cells (DCs) constitute a heterogeneous group of antigen-presenting leukocytes important in activation of both innate and adaptive immunity. We studied the gene expression patterns of DCs incubated with reagents inducing their activation or inhibition. Total RNA was isolated from DCs and gene expression profiling was performed with oligonucleotide microarrays. Using a supervised learning algorithm based on Random Forest, we generated a molecular signature of inflammation from a training set of 77 samples. We then validated this molecular signature in a testing set of 38 samples. Supervised analysis identified a set of 44 genes that distinguished very accurately between inflammatory and non inflammatory samples. The diagnostic performance of the signature genes was assessed against an independent set of samples, by qRT-PCR. Our findings suggest that the gene expression signature of DCs can provide a molecular classification for use in the selection of anti-inflammatory or adjuvant molecules with specific effects on DC activity.

Genome-wide analysis of immune activation in human T and B cells reveals distinct classes of alternatively spliced genes

Alternative splicing of pre-mRNA is a mechanism that increases the protein diversity of a single gene by differential exon inclusion/exclusion during post-transcriptional processing. While alternative splicing is established to occur during lymphocyte activation, little is known about the role it plays during the immune response. Our study is among the first reports of a systematic genome-wide analysis of activated human T and B lymphocytes using whole exon DNA microarrays integrating alternative splicing and differential gene expression. Purified human CD2⁺ T or CD19⁺ B cells were activated using protocols to model the early events in post-transplant allograft immunity and sampled as a function of time during the process of immune activation. Here we show that 3 distinct classes of alternatively spliced and/or differentially expressed genes change in an ordered manner as a function of immune activation. We mapped our results to function-based canonical pathways and demonstrated that some are populated by only one class of genes, like integrin signaling, while other pathways, such as purine metabolism and T cell receptor signaling, are populated by all three classes of genes. Our studies augment the current view of T and B cell activation in immunity that has been based exclusively upon differential gene expression by providing evidence for a large number of molecular networks populated as a function of time and activation by alternatively spliced genes, many of which are constitutively expressed.

Lyase to live by: sphingosine phosphate lyase as a therapeutic target

BACKGROUND

Sphingosine 1-phosphate (S1P) is a bioactive lipid that regulates cell proliferation, survival and migration and plays an essential role in angiogenesis and lymphocyte trafficking. S1P levels in the circulation and tissues are tightly regulated for proper cell functioning, and dysregulation of this system may contribute to the pathophysiology of certain human diseases. Sphingosine phosphate lyase (SPL) irreversibly degrades S1P and thereby acts as a gatekeeper that regulates S1P signaling by modulating intracellular S1P levels and the chemical S1P gradient that exists between lymphoid organs and circulating blood and lymph. However, SPL also generates biochemical products that may be relevant in human disease. SPL has been directly implicated in various physiological and pathological processes, including cell stress responses, cancer, immunity, hematopoietic function, muscle homeostasis, inflammation and development.

OBJECTIVE/METHODS

This review summarizes the current know-ledge of SPL structure, function and regulation, its involvement in various disease states and currently available small molecules known to modulate SPL activity.

RESULTS/CONCLUSION

This review provides evidence that SPL is a potential target for pharmacological manipulation for the treatment of malignant, autoimmune, inflammatory and other diseases.

Hydrolysis of the phosphoanhydride linkage of cyclic ADP-ribose by the Mn(2+)-dependent ADP-ribose/CDP-alcohol pyrophosphatase

Cyclic ADP-ribose (cADPR) metabolism in mammals is catalyzed by NAD glycohydrolases (NADases) that, besides forming ADP-ribose, form and hydrolyze the N(1)-glycosidic linkage of cADPR. Thus far, no cADPR phosphohydrolase was known. We tested rat ADP-ribose/CDP-alcohol pyrophosphatase (ADPRibase-Mn) and found that cADPR is an ADPRibase-Mn ligand and substrate. ADPRibase-Mn activity on cADPR was 65-fold less efficient than on ADP-ribose, the best substrate. This is similar to the ADP-ribose/cADPR formation ratio by NADases. The product of cADPR phosphohydrolysis by ADPRibase-Mn was N(1)-(5-phosphoribosyl)-AMP, suggesting a novel route for cADPR turnover.

Mn2+-dependent ADP-ribose/CDP-alcohol pyrophosphatase: a novel metallophosphoesterase family preferentially expressed in rodent immune cells

ADPRibase-Mn (Mn2+-dependent ADP-ribose/CDP-alcohol pyrophosphatase) was earlier isolated from rat liver supernatants after separation from ADPRibase-I and ADPRibase-II (Mg2+-activated ADP-ribose pyrophosphatases devoid of CDP-alcohol pyrophosphatase activity). The last mentioned are putative Nudix hydrolases, whereas the molecular identity of ADPRibase-Mn is unknown. MALDI (matrix-assisted laser-desorption ionization) MS data from rat ADPRibase-Mn pointed to a hypothetical protein that was cloned and expressed and showed the expected specificity. It is encoded by the RGD1309906 rat gene, which so far has been annotated simply as 'hydrolase'. ADPRibase-Mn is not a Nudix hydrolase, but it shows the sequence and structural features typical of the metallophosphoesterase superfamily. It may constitute a protein family of its own, the members of which appear to be specific to vertebrates, plants and algae. ADP-ribose was successfully docked to a model of rat ADPRibase-Mn, revealing its putative active centre. Microarray data from the GEO (Gene Expression Omnibus) database indicated that the mouse gene 2310004I24Rik, an orthologue of RGD1309906, is preferentially expressed in immune cells. This was confirmed by Northern-blot and activity assay of ADPRibase-Mn in rat tissues. A possible role of ADPRibase-Mn in immune cell signalling is suggested by the second-messenger role of ADP-ribose, which activates TRPM2 (transient receptor potential melastatin channel-2) ion channels as a mediator of oxidative/nitrosative stress, and by the signalling function assigned to many of the microarray profile neighbours of 2310004I24Rik. Furthermore, the influence of ADPRibase-Mn on the CDP-choline or CDP-ethanolamine pathways of phospholipid biosynthesis cannot be discounted.

Sphingosine-1-phosphate lyase in development and disease: sphingolipid metabolism takes flight

Sphingosine-1-phosphate lyase (SPL) is a highly conserved enzyme that catalyses the final step of sphingolipid degradation, namely the irreversible cleavage of the carbon chain at positions 2-3 of a long-chain base phosphate (LCBP), thereby yielding a long-chain aldehyde and phosphoethanolamine. LCBPs are potent signaling molecules involved in cell proliferation, survival, migration, cell-cell interactions and cell stress responses. Therefore, tight regulation of LCBP signaling is required for proper cell function, and perturbations of this system can lead to alterations in biological processes including development, reproduction and physiology. SPL is a key enzyme in regulating the intracellular and circulating levels of LCBPs and is, therefore, gaining attention as a putative target for pharmacological intervention. This review provides an overview of our current understanding of SPL structure and function, mechanisms involved in SPL regulation and the role of SPL in development and disease.

Novel insights into the relationships between dendritic cell subsets in human and mouse revealed by genome-wide expression profiling

Genome-wide expression profiling of mouse and human leukocytes reveal conserved transcriptional programs of plasmacytoid or conventional dendritic cell subsets.

Background

Dendritic cells (DCs) are a complex group of cells that play a critical role in vertebrate immunity. Lymph-node resident DCs (LN-DCs) are subdivided into conventional DC (cDC) subsets (CD11b and CD8α in mouse; BDCA1 and BDCA3 in human) and plasmacytoid DCs (pDCs). It is currently unclear if these various DC populations belong to a unique hematopoietic lineage and if the subsets identified in the mouse and human systems are evolutionary homologs. To gain novel insights into these questions, we sought conserved genetic signatures for LN-DCs and in vitro derived granulocyte-macrophage colony stimulating factor (GM-CSF) DCs through the analysis of a compendium of genome-wide expression profiles of mouse or human leukocytes.

Results

We show through clustering analysis that all LN-DC subsets form a distinct branch within the leukocyte family tree, and reveal a transcriptomal signature evolutionarily conserved in all LN-DC subsets. Moreover, we identify a large gene expression program shared between mouse and human pDCs, and smaller conserved profiles shared between mouse and human LN-cDC subsets. Importantly, most of these genes have not been previously associated with DC function and many have unknown functions. Finally, we use compendium analysis to re-evaluate the classification of interferon-producing killer DCs, lin^-CD16⁺HLA-DR⁺ cells and in vitro derived GM-CSF DCs, and show that these cells are more closely linked to natural killer and myeloid cells, respectively.

Conclusion

Our study provides a unique database resource for future investigation of the evolutionarily conserved molecular pathways governing the ontogeny and functions of leukocyte subsets, especially DCs.

Two polymorphisms of presenilin-2 gene (PSEN2) 5' regulatory region are not associated with Alzheimer's disease (AD) in the Polish population

Presenilin 2 gene (PSEN2) is one of the causative genes for familial Alzheimer's disease. A delA polymorphism located in PSEN2 promoter was proposed to be a risk factor for early-onset AD. We examined association between AD and PSEN2 polymorphisms located in two 5'UTR regions in group of 217 late-onset AD patients, 109 mild cognitive impairment patients, and 225 non-demented control subjects. No significant differences for genotype and allele distributions of a delA and a novel insAC polymorphisms in the studied groups as compared to controls were observed. Univariate and multivariate risk estimation shows that neither delA, insAC alleles nor the genotypes are risk factors for AD. No significant interaction between the APOE4 and PSEN2 polymorphisms was found. A bioinformatic analysis showed that delA polymorphism influences binding sites of transcription factors involved in the cellular processes related to AD. The rare variants identified in exon 3 of the PSEN2 could have a potential influence on PSEN2 transcript splicing.

Increased Transcriptional Activity of Milk-Related Genes following the Active Phase of Experimental Autoimmune Encephalomyelitis and Multiple Sclerosis

We analyzed global transcriptional changes in the lymph nodes of mice with experimental autoimmune encephalomyelitis in a longitudinal fashion. Most of the transcriptional activity was observed between 3 and 5 days postimmunization. After that period, gene expression changes decayed sharply back to baseline levels. A comparison of transcriptional profiles between immunized and control mice at the time of peak disease activity revealed 266 transcripts, mostly involved in cell-cell interaction and protein synthesis. When the same comparison was performed at the time of recovery from an attack, increased expression of genes coding for milk components were identified. Specifically, casein alpha (Csn1s1), beta (Csn2), gamma (Csn1s2a), and kappa (Csn3), in addition to lactoalbumin alpha and extracellular proteinase were elevated >3-fold in immunized animals compared with CFA-injected controls. We confirmed these findings by quantitative RT-PCR and immunostaining of Csn3. Interestingly, the expression of Csn3 was also found elevated in the blood of multiple sclerosis (MS) patients after a relapse. Altogether, our data suggest that increased production of milk-related transcripts in the lymph nodes and blood succeeds an inflammatory event in experimental autoimmune encephalomyelitis and MS. The potential role of lactogenic hormones in MS is discussed.

Dynamic SPR monitoring of yeast nuclear protein binding to a cis-regulatory element

Gene expression is controlled by protein complexes binding to short specific sequences of DNA, called cis-regulatory elements. Expression of most eukaryotic genes is controlled by dozens of these elements. Comprehensive identification and monitoring of these elements is a major goal of genomics. In pursuit of this goal, we are developing a surface plasmon resonance (SPR) based assay to identify and monitor cis-regulatory elements. To test whether we could reliably monitor protein binding to a regulatory element, we immobilized a 16bp region of Saccharomyces cerevisiae chromosome 5 onto a gold surface. This 16bp region of DNA is known to bind several proteins and thought to control expression of the gene RNR1, which varies through the cell cycle. We synchronized yeast cell cultures, and then sampled these cultures at a regular interval. These samples were processed to purify nuclear lysate, which was then exposed to the sensor. We found that nuclear protein binds this particular element of DNA at a significantly higher rate (as compared to unsynchronized cells) during G1 phase. Other time points show levels of DNA-nuclear protein binding similar to the unsynchronized control. We also measured the apparent association complex of the binding to be 0.014s(-1). We conclude that (1) SPR-based assays can monitor DNA-nuclear protein binding and that (2) for this particular cis-regulatory element, maximum DNA-nuclear protein binding occurs during G1 phase.

Characterization of a family of novel cysteine- serine-rich nuclear proteins (CSRNP)

Gene array analysis has been widely used to identify genes induced during T cell activation. Our studies identified an immediate early gene that is strongly induced in response to IL-2 in mouse T cells which we named cysteine- serine-rich nuclear protein-1 (CSRNP-1). The human ortholog was previously identified as an AXIN1 induced gene (AXUD1). The protein does not contain sequence defined domains or motifs annotated in public databases, however the gene is a member of a family of three mammalian genes that share conserved regions, including cysteine- and serine-rich regions and a basic domain, they encode nuclear proteins, possess transcriptional activation domain and bind the sequence AGAGTG. Consequently we propose the nomenclature of CSRNP-1, -2 and -3 for the family. To elucidate the physiological functions of CSRNP-1, -2 and -3, we generated mice deficient for each of these genes by homologous recombination in embryonic stem cells. Although the CSRNP proteins have the hallmark of transcription factors and CSRNP-1 expression is highly induced by IL-2, deletion of the individual genes had no obvious consequences on normal mouse development, hematopoiesis or T cell functions. However, combined deficiencies cause partial neonatal lethality suggesting that the genes have redundant functions.

Post hoc pattern matching: assigning significance to statistically defined expression patterns in single channel microarray data

Background

Researchers using RNA expression microarrays in experimental designs with more than two treatment groups often identify statistically significant genes with ANOVA approaches. However, the ANOVA test does not discriminate which of the multiple treatment groups differ from one another. Thus, post hoc tests, such as linear contrasts, template correlations, and pairwise comparisons are used. Linear contrasts and template correlations work extremely well, especially when the researcher has a priori information pointing to a particular pattern/template among the different treatment groups. Further, all pairwise comparisons can be used to identify particular, treatment group-dependent patterns of gene expression. However, these approaches are biased by the researcher's assumptions, and some treatment-based patterns may fail to be detected using these approaches. Finally, different patterns may have different probabilities of occurring by chance, importantly influencing researchers' conclusions about a pattern and its constituent genes.

Results

We developed a four step, post hoc pattern matching (PPM) algorithm to automate single channel gene expression pattern identification/significance. First, 1-Way Analysis of Variance (ANOVA), coupled with post hoc 'all pairwise' comparisons are calculated for all genes. Second, for each ANOVA-significant gene, all pairwise contrast results are encoded to create unique pattern ID numbers. The # genes found in each pattern in the data is identified as that pattern's 'actual' frequency. Third, using Monte Carlo simulations, those patterns' frequencies are estimated in random data ('random' gene pattern frequency). Fourth, a Z-score for overrepresentation of the pattern is calculated ('actual' against 'random' gene pattern frequencies). We wrote a Visual Basic program (StatiGen) that automates PPM procedure, constructs an Excel workbook with standardized graphs of overrepresented patterns, and lists of the genes comprising each pattern. The visual basic code, installation files for StatiGen, and sample data are available as supplementary material.

Conclusion

The PPM procedure is designed to augment current microarray analysis procedures by allowing researchers to incorporate all of the information from post hoc tests to establish unique, overarching gene expression patterns in which there is no overlap in gene membership. In our hands, PPM works well for studies using from three to six treatment groups in which the researcher is interested in treatment-related patterns of gene expression. Hardware/software limitations and extreme number of theoretical expression patterns limit utility for larger numbers of treatment groups. Applied to a published microarray experiment, the StatiGen program successfully flagged patterns that had been manually assigned in prior work, and further identified other gene expression patterns that may be of interest. Thus, over a moderate range of treatment groups, PPM appears to work well. It allows researchers to assign statistical probabilities to patterns of gene expression that fit a priori expectations/hypotheses, it preserves the data's ability to show the researcher interesting, yet unanticipated gene expression patterns, and assigns the majority of ANOVA-significant genes to non-overlapping patterns.

Sphingosine-1-phosphate lyase in immunity and cancer: silencing the siren

Sphingosine-1-phosphate (S1P) is a bioactive lipid that promotes cell survival, proliferation and migration, platelet aggregation, mediates ischemic preconditioning, and is essential for angiogenesis and lymphocyte trafficking. Sphingosine-1-phosphate lyase (SPL) is the enzyme responsible for the irreversible degradation of S1P and is, thus, in a strategic position to regulate these same processes by removing available S1P signaling pools, that is, silencing the siren. In fact, recent studies have implicated SPL in the regulation of immunity, cancer surveillance and other physiological processes. Here, we summarize the current understanding of SPL function and regulation, and discuss how SPL might facilitate cancer chemoprevention and serve as a target for modulation of immune responses in transplantation settings and in the treatment of autoimmune disease.

IFN Regulatory Factor-2 Regulates Macrophage Apoptosis through a STAT1/3- and Caspase-1-Dependent Mechanism

IFN regulatory factor (IRF)-2(-/-) mice are significantly more resistant to LPS challenge than wild-type littermates, and this was correlated with increased numbers of apoptotic Kupffer cells. To assess the generality of this observation, and to understand the role of IRF-2 in apoptosis, responses of peritoneal macrophages from IRF-2(+/+) and IRF-2(-/-) mice to apoptotic stimuli, including the fungal metabolite, gliotoxin, were compared. IRF-2(-/-) macrophages exhibited a consistently higher incidence of apoptosis that failed to correlate with caspase-3/7 activity. Using microarray gene expression profiling of liver RNA samples derived from IRF-2(+/+) and IRF-2(-/-) mice treated with saline or LPS, we identified >40 genes that were significantly down-regulated in IRF-2(-/-) mice, including Stat3, which has been reported to regulate apoptosis. Compared with IRF-2(+/+) macrophages, STAT3alpha mRNA was up-regulated constitutively or after gliotoxin treatment of IRF-2(-/-) macrophages, whereas STAT3beta mRNA was down-regulated. Phospho-Y705-STAT3, phospho-S727-STAT1, and phospho-p38 protein levels were also significantly higher in IRF-2(-/-) than control macrophages. Activation of the STAT signaling pathway has been shown to elicit expression of CASP1 and apoptosis. IRF-2(-/-) macrophages exhibited increased basal and gliotoxin-induced caspase-1 mRNA expression and enhanced caspase-1 activity. Pharmacologic inhibition of STAT3 and caspase-1 abolished gliotoxin-induced apoptosis in IRF-2(-/-) macrophages. A novel IFN-stimulated response element, identified within the murine promoter of Casp1, was determined to be functional by EMSA and supershift analysis. Collectively, these data support the hypothesis that IRF-2 acts as a transcriptional repressor of Casp1, and that the absence of IRF-2 renders macrophages more sensitive to apoptotic stimuli in a caspase-1-dependent process.

Toxicogenomics in the assessment of immunotoxicity

Microarray analysis is used for simultaneous measurement of expression of thousands of genes in a given sample and as such extends and deepens our understanding of biological processes. Application of the technique in toxicology is referred to as toxicogenomics. The examples of assessment of immunotoxicity by gene expression profiling presented and discussed here, show that microarray analysis is able to detect known and novel effects of a wide range of immunomodulating agents. Besides the elucidation of mechanisms of action, toxicogenomics is also applied to predict consequences of exposing biological systems to toxic agents. Successful attempts to classify compounds using signature gene expression profiles have been reported. These did, however, not specifically focus on immunotoxicity. Databases containing expression profiles can facilitate the applications of toxicogenomics. Platforms and methodologies for gene expression profiling may vary, however, hampering data compiling across different laboratories. Therefore, attention is paid to standardization of the generation, reporting, and management of microarray data. Obtained gene expression profiles should be anchored to pathological and functional endpoints for correct interpretation of results. These issues are also important when using toxicogenomics in risk assessment. The application of toxicogenomics in evaluation of immunotoxicity is thus not yet without challenges. It already contributes to the understanding of immunotoxic processes and the development of in vitro screening assays, though, and is therefore expected to be of value for mechanistic insight into immunotoxicity and hazard identification of existing and novel compounds.

Selective gene silencing in activated leukocytes by targeting siRNAs to the integrin lymphocyte function-associated antigen-1

Silencing gene expression by RNAi is a powerful method for exploring gene function and validating drug targets and potentially for therapy. Lymphocytes and other primary blood cells are resistant to lipid-based transfection in vitro and are difficult to target in vivo. We show here that antibody-protamine fusion proteins targeting the human integrin lymphocyte function-associated antigen-1 (LFA-1) efficiently deliver siRNAs and specifically induce silencing in primary lymphocytes, monocytes, and dendritic cells. Moreover, a fusion protein constructed from an antibody that preferentially recognizes activation-dependent conformational changes in LFA-1 selectively targets activated leukocytes and can be used to suppress gene expression and cell proliferation only in activated lymphocytes. The siRNA-fusion protein complexes do not cause lymphocyte activation or induce IFN responses. K562 cells expressing latent WT or constitutively activated LFA-1 engrafted in the lungs of SCID mice are selectively targeted by intravenously injected fusion protein-siRNA complexes, demonstrating the potential in vivo applicability of LFA-1-directed siRNA delivery.

Molecular characterization of the immune system: emergence of proteins, processes, and domains

Many genes and proteins are required to carry out the processes of innate and adaptive immunity. For many studies, including systems biology, it is necessary to have a clear and comprehensive definition of the immune system, including the genes and proteins that take part in immunological processes. We have identified and cataloged a large portion of the human immunology-related genes, which we call the essential immunome. The 847 identified genes and proteins were annotated, and their chromosomal localizations were compared to the mouse genome. Relation to disease was also taken into account. We identified numerous pseudogenes, many of which are expressed, and found two putative new genes. We also carried out an evolutionary analysis of immune processes based on gene orthologs to gain an overview of the evolutionary past and molecular present of the human immune system. A list of genes and proteins were compiled. A comprehensive characterization of the member genes and proteins, including the corresponding pseudogenes is presented. Immunome genes were found to have three types of emergence in independent studies of their ontologies, domains, and functions.

Gene Expression Alterations in Immune System Pathways following Exposure to Immunosuppressive Chemicals

Exposure to environmental agents can affect a number of adverse immunological outcomes, including changes in the incidence of infectious disease. Diethylstilbestrol (DES), dexamethasone (DEX), cyclophosphamide, and 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) are immunosuppressive chemicals that can induce similar pathophysiological end points in the thymus; however, the mechanism of toxicity is different for each compound. We examined differential gene expression in the spleen and thymus following chemical exposure and correlated these changes with alterations in functional immune end points and our knowledge of the known mechanisms of action. RNA from the spleen and thymus has been analyzed using Illumina Sentrix arrays and BeadStudio software. Preliminary data suggest that DES induced the greatest number of gene changes in the spleen, while DEX induced the most changes in the thymus. In both spleen and thymus, genomic analysis revealed gene expression changes that were common to multiple chemicals and that may be associated with xenobiotic-induced immune system perturbations, including alterations in genes associated with apoptosis, antigen processing and presentation, and response to biotic stimulus. This was particularly evident in the thymus, where there were many similarities in the expression profiles, as well as gene alterations unique to a single compound. In contrast, expression profiles in spleen were more distinct. The category of genes most profoundly affected by all four chemicals was response to biotic stimulus: there were both clusters of genes modulated by multiple chemicals and genes altered by a single chemical. The distinct gene profiles may specifically relate to cellular targets and mechanism of action.

Multiple regulatory regions and tissue-specific transcription initiation mediate the expression of NEMO/IKKgamma gene

NEMO/IKKgamma gene, which is responsible of two allelic diseases in human, EDA-ID and IP, encodes for a protein with a central regulatory role in the activation of the NF-kB pathway. We here provide insights into the molecular mechanism governing NEMO/IKKgamma expression. We mapped 4 distinctive NEMO/IKKgamma transcription start sites each corresponding to an alternative first exon, controlled by two conserved promoters. A distal promoter, named promoter A, located 10 kb upstream of the coding region and a proximal promoter, promoter B, with strong bi-directional activity driving also the transcription of G6PD gene in the opposite direction. The promoter B is housekeeping, it is embedded in a CpG island, required for proper expression and it is down-regulated by methylation. The promoter A is active in cells of hepatic origin and it directs transcription of the main NEMO/IKKgamma 5' UTR alternative transcript in liver, which starts at a tissue-specific site. Qualitative and quantitative expression analysis revealed that each NEMO/IKKgamma 5' UTR alternative transcript has different expression profiles indicating that the control of NEMO/IKKgamma expression is mediated through tissue-specific transcription initiation sites and multiple regulatory regions.

Immunotoxicogenomics: the potential of genomics technology in the immunotoxicity risk assessment process

Evaluation of xenobiotic-induced changes in gene expression as a method to identify and classify potential toxicants is being pursued by industry and regulatory agencies worldwide. A workshop was held at the Research Triangle Park campus of the Environmental Protection Agency to discuss the current state-of-the-science of "immunotoxicogenomics" and to explore the potential role of genomics techniques for immunotoxicity testing. The genesis of the workshop was the current lack of widely accepted triggering criteria for Tier 1 immunotoxicity testing in the context of routine toxicity testing data, the realization that traditional screening methods would require an inordinate number of animals and are inadequate to handle the number of chemicals that may need to be screened (e.g., high production volume compounds) and the absence of an organized effort to address the state-of-the-science of toxicogenomics in the identification of immunotoxic compounds. The major focus of the meeting was on the theoretical and practical utility of genomics techniques to (1) replace or supplement current immunotoxicity screening procedures, (2) provide insight into potential modes or mechanisms of action, and (3) provide data suitable for immunotoxicity hazard identification or risk assessment. The latter goal is of considerable interest to a variety of stakeholders as a means to reduce animal use and to decrease the cost of conducting and interpreting standard toxicity tests. A number of data gaps were identified that included a lack of dose response and kinetic data for known immunotoxic compounds and a general lack of data correlating genomic alterations to functional changes observed in vivo. Participants concluded that a genomics approach to screen chemicals for immunotoxic potential or to generate data useful to risk assessors holds promise but that routine use of these methods is years in the future. However, recent progress in molecular immunology has made mode and mechanism of action studies much more practical. Furthermore, a variety of published immunotoxicity studies suggest that microarray analysis is already a practical means to explore pathway-level changes that lead to altered immune function. To help move the science of immunotoxicogenomics forward, a partnership of industry, academia, and government was suggested to address data gaps, validation, quality assurance, and protocol development.

Gene expression microarrays: glimpses of the immunological genome

Successful microarray experimentation can generate enormous amounts of data, potentially very rich but also very unwieldy. Bold outlooks and new methods for data analysis and presentation should yield additional insight into the complexities of the immune system.

Construction of quantitative proteome reference maps of mouse spleen and lymph node based on two-dimensional gel electrophoresis

Quantitative features of the proteome are extremely useful for studying cellular processes at a molecular level. In this study, we attempted to construct quantitative reference proteome maps of the mouse spleen and lymph node based on 2-DE followed by protein identification using MS. We analyzed more than 1000 spots on the 2-DE images and consequently were able to determine that 919 spots were derived from 328 different genes. To obtain statistically reliable information of the protein levels from these 2-DE images, we measured the volumes of the respective spots on 2-DE images obtained by four to six independent experimental runs. These measurements were used to calculate the variability of the volumes of the respective spots on 2-DE following subcellular fractionation, which enabled us to discriminate differentially produced proteins from those within the range of intrinsic variability. More importantly, while the 2-DE data have been traditionally collected in a gel image-based manner, the resultant quantitative 2-DE data could be analyzed using the same procedure as that for mRNA expression profiles. This greatly assists in bridging the gap between the analyses of transcriptomes and proteomes and enables the integration of this data on the same informational platform.

Functional genomic responses to cystic fibrosis transmembrane conductance regulator (CFTR) and CFTR(delta508) in the lung

Cystic fibrosis (CF), a common lethal pulmonary disorder in Caucasians, is caused by mutations in the cystic fibrosis transmembrane conductance regulator gene (CFTR) that disturbs fluid homeostasis and host defense in target organs. The effects of CFTR and delta508-CFTR were assessed in transgenic mice that 1) lack CFTR expression (Cftr-/-); 2) express the human delta508 CFTR (CFTR(delta508)); 3) overexpress the normal human CFTR (CFTR(tg)) in respiratory epithelial cells. Genes were selected from Affymetrix Murine Gene-Chips analysis and subjected to functional classification, k-means clustering, promoter cis-elements/modules searching, literature mining, and pathway exploring. Genomic responses to Cftr-/- were not corrected by expression of CFTR(delta508). Genes regulating host defense, inflammation, fluid and electrolyte transport were similarly altered in Cftr-/- and CFTR(delta508) mice. CFTR(delta508) induced a primary disturbance in expression of genes regulating redox and antioxidant systems. Genomic responses to CFTR(tg) were modest and were not associated with lung pathology. CFTR(tg) and CFTR(delta508) induced genes encoding heat shock proteins and other chaperones but did not activate the endoplasmic reticulum-associated degradation pathway. RNAs encoding proteins that directly interact with CFTR were identified in each of the CFTR mouse models, supporting the hypothesis that CFTR functions within a multiprotein complex whose members interact at the level of protein-protein interactions and gene expression. Promoters of genes influenced by CFTR shared common regulatory elements, suggesting that their co-expression may be mediated by shared regulatory mechanisms. Genes and pathways involved in the response to CFTR may be of interest as modifiers of CF.

CisMols Analyzer: identification of compositionally similar cis-element clusters in ortholog conserved regions of coordinately expressed genes

Combinatorial interactions of sequence-specific trans-acting factors with localized genomic cis-element clusters are the principal mechanism for regulating tissue-specific and developmental gene expression. With the emergence of expanding numbers of genome-wide expression analyses, the identification of the cis-elements responsible for specific patterns of transcriptional regulation represents a critical area of investigation. Computational methods for the identification of functional cis-regulatory modules are difficult to devise, principally because of the short length and degenerate nature of individual cis-element binding sites and the inherent complexity that is generated by combinatorial interactions within cis-clusters. Filtering candidate cis-element clusters based on phylogenetic conservation is helpful for an individual ortholog gene pair, but combining data from cis-conservation and coordinate expression across multiple genes is a more difficult problem. To approach this, we have extended an ortholog gene-pair database with additional analytical architecture to allow for the analysis and identification of maximal numbers of compositionally similar and phylogenetically conserved cis-regulatory element clusters from a list of user-selected genes. The system has been successfully tested with a series of functionally related and microarray profile-based co-expressed ortholog pairs of promoters and genes using known regulatory regions as training sets and co-expressed genes in the olfactory and immunohematologic systems as test sets. CisMols Analyzer is accessible via a Web interface at .