Genetic architecture of ambulatory blood pressure in the general population: insights from cardiovascular gene-centric array.

Genetic determinants of blood pressure are poorly defined. We undertook a large-scale, gene-centric analysis to identify loci and pathways associated with ambulatory systolic and diastolic blood pressure. We measured 24-hour ambulatory blood pressure in 2020 individuals from 520 white European nuclear families (the Genetic Regulation of Arterial Pressure of Humans in the Community Study) and genotyped their DNA using the Illumina HumanCVD BeadChip array, which contains ≈50 000 single nucleotide polymorphisms in >2000 cardiovascular candidate loci. We found a strong association between rs13306560 polymorphism in the promoter region of MTHFR and CLCN6 and mean 24-hour diastolic blood pressure; each minor allele copy of rs13306560 was associated with 2.6 mm Hg lower mean 24-hour diastolic blood pressure (P=1.2×10(-8)). rs13306560 was also associated with clinic diastolic blood pressure in a combined analysis of 8129 subjects from the Genetic Regulation of Arterial Pressure of Humans in the Community Study, the CoLaus Study, and the Silesian Cardiovascular Study (P=5.4×10(-6)). Additional analysis of associations between variants in gene ontology-defined pathways and mean 24-hour blood pressure in the Genetic Regulation of Arterial Pressure of Humans in the Community Study showed that cell survival control signaling cascades could play a role in blood pressure regulation. There was also a significant overrepresentation of rare variants (minor allele frequency: <0.05) among polymorphisms showing at least nominal association with mean 24-hour blood pressure indicating that a considerable proportion of its heritability may be explained by uncommon alleles. Through a large-scale gene-centric analysis of ambulatory blood pressure, we identified an association of a novel variant at the MTHFR/CLNC6 locus with diastolic blood pressure and provided new insights into the genetic architecture of blood pressure.

Lineage relationship of prostate cancer cell types based on gene expression

Background: Prostate tumor heterogeneity is a major factor in disease management. Heterogeneity could be due to multiple cancer cell types with distinct gene expression. Of clinical importance is the so-called cancer stem cell type. Cell type-specific transcriptomes are used to examine lineage relationship among cancer cell types and their expression similarity to normal cell types including stem/progenitor cells. Methods: Transcriptomes were determined by Affymetrix DNA array analysis for the following cell types. Putative prostate progenitor cell populations were characterized and isolated by expression of the membrane transporter ABCG2. Stem cells were represented by embryonic stem and embryonal carcinoma cells. The cancer cell types were Gleason pattern 3 (glandular histomorphology) and pattern 4 (aglandular) sorted from primary tumors, cultured prostate cancer cell lines originally established from metastatic lesions, xenografts LuCaP 35 (adenocarcinoma phenotype) and LuCaP 49 (neuroendocrine/small cell carcinoma) grown in mice. No detectable gene expression differences were detected among serial passages of the LuCaP xenografts. Results: Based on transcriptomes, the different cancer cell types could be clustered into a luminal-like grouping and a non-luminal-like (also not basal-like) grouping. The non-luminal-like types showed expression more similar to that of stem/progenitor cells than the luminal-like types. However, none showed expression of stem cell genes known to maintain stemness. Conclusions: Non-luminal-like types are all representatives of aggressive disease, and this could be attributed to the similarity in overall gene expression to stem and progenitor cell types.

Gene expression relationship between prostate cancer cells of Gleason 3, 4 and normal epithelial cells as revealed by cell type-specific transcriptomes

Background: Prostate cancer cells in primary tumors have been typed CD10(-)/CD13(-)/CD24(hi)/CD26(+)/CD38(lo)/CD44(-)/CD104(-). This CD phenotype suggests a lineage relationship between cancer cells and luminal cells. The Gleason grade of tumors is a descriptive of tumor glandular differentiation. Higher Gleason scores are associated with treatment failure. Methods: CD26(+) cancer cells were isolated from Gleason 3+3 (G3) and Gleason 4+4 (G4) tumors by cell sorting, and their gene expression or transcriptome was determined by Affymetrix DNA array analysis. Dataset analysis was used to determine gene expression similarities and differences between G3 and G4 as well as to prostate cancer cell lines and histologically normal prostate luminal cells. Results: The G3 and G4 transcriptomes were compared to those of prostatic cell types of non-cancer, which included luminal, basal, stromal fibromuscular, and endothelial. A principal components analysis of the various transcriptome datasets indicated a closer relationship between luminal and G3 than luminal and G4. Dataset comparison also showed that the cancer transcriptomes differed substantially from those of prostate cancer cell lines. Conclusions: Genes differentially expressed in cancer are potential biomarkers for cancer detection, and those differentially expressed between G3 and G4 are potential biomarkers for disease stratification given that G4 cancer is associated with poor outcomes. Differentially expressed genes likely contribute to the prostate cancer phenotype and constitute the signatures of these particular cancer cell types.

Alterations in gene expression and activity during squamous cell carcinoma development

This study focuses on characterizing the genetic and biological alterations associated with squamous cell carcinoma development. Normal human epidermal keratinocytes (HEKs), cells isolated from a preneoplastic lesion (IEC-1), and two neoplastic cell lines, SCC-25 and COLD-16, were grown as raft cultures, and their gene expression profiles were screened using cDNA arrays. Our data indicated that the expression levels of at least 37 genes were significantly (P less than or equal to 0.05; 1.9% of genes screened) altered in neoplastic cells compared with normal cells. Of these genes, 10 genes were up-regulated and 27 genes were down-regulated in the neoplastic cells. In addition, 51% of the genes altered in the neoplastic cells were already altered in the preneoplastic IEC-1 cells. Immunohistochemical staining of patient tumors was used to verify the cDNA array analysis. Our analysis indicated that alterations in genes associated with extracellular matrix production and apoptosis are disrupted in preneoplastic cells, whereas later stages of neoplasia are associated with alterations in gene expression for genes involved in DNA repair or epidermal growth factor (EGF) receptor/mitogen-activated protein kinase kinase (MAPKK)/MAPK/activator protein-1 (AP-1) signaling. Subsequent functional analysis of the alterations in expression of the EGF receptor/MAPKK/MAPK/AP-1 genes suggested they did not contribute to the neoplastic phenotype.

The Interleukin-1 receptor antagonist is a direct target gene of PPARalpha in liver.

BACKGROUND/AIMS: The Peroxisome Proliferator-Activated Receptor (PPAR) alpha belongs to the superfamily of Nuclear Receptors and plays an important role in numerous cellular processes, including lipid metabolism. It is known that PPARalpha also has an anti-inflammatory effect, which is mainly achieved by down-regulating pro-inflammatory genes. The objective of this study was to further characterize the role of PPARalpha in inflammatory gene regulation in liver. RESULTS: According to Affymetrix micro-array analysis, the expression of various inflammatory genes in liver was decreased by treatment of mice with the synthetic PPARalpha agonist Wy14643 in a PPARalpha-dependent manner. In contrast, expression of Interleukin-1 receptor antagonist (IL-1ra), which was acutely stimulated by LPS treatment, was induced by PPARalpha. Up-regulation of IL-1ra by LPS was lower in PPARalpha -/- mice compared to Wt mice. Transactivation and chromatin immunoprecipitation studies identified IL-1ra as a direct positive target gene of PPARalpha with a functional PPRE present in the promoter. Up-regulation of IL-1ra by PPARalpha was conserved in human HepG2 hepatoma cells and the human monocyte/macrophage THP-1 cell line. CONCLUSIONS: In addition to down-regulating expression of pro-inflammatory genes, PPARalpha suppresses the inflammatory response by direct up-regulation of genes with anti-inflammatory properties.

Classification of human astrocytic gliomas on the basis of gene expression: a correlated group of genes with angiogenic activity emerges as a strong predictor of subtypes.

The development of targeted treatment strategies adapted to individual patients requires identification of the different tumor classes according to their biology and prognosis. We focus here on the molecular aspects underlying these differences, in terms of sets of genes that control pathogenesis of the different subtypes of astrocytic glioma. By performing cDNA-array analysis of 53 patient biopsies, comprising low-grade astrocytoma, secondary glioblastoma (respective recurrent high-grade tumors), and newly diagnosed primary glioblastoma, we demonstrate that human gliomas can be differentiated according to their gene expression. We found that low-grade astrocytoma have the most specific and similar expression profiles, whereas primary glioblastoma exhibit much larger variation between tumors. Secondary glioblastoma display features of both other groups. We identified several sets of genes with relatively highly correlated expression within groups that: (a). can be associated with specific biological functions; and (b). effectively differentiate tumor class. One prominent gene cluster discriminating primary versus nonprimary glioblastoma comprises mostly genes involved in angiogenesis, including VEGF fms-related tyrosine kinase 1 but also IGFBP2, that has not yet been directly linked to angiogenesis. In situ hybridization demonstrating coexpression of IGFBP2 and VEGF in pseudopalisading cells surrounding tumor necrosis provided further evidence for a possible involvement of IGFBP2 in angiogenesis. The separating groups of genes were found by the unsupervised coupled two-way clustering method, and their classification power was validated by a supervised construction of a nearly perfect glioma classifier.

Infrequent promoter methylation of the MGMT gene in liver metastases from uveal melanoma.

Uveal melanoma is associated with a high mortality rate once metastases occur, with over >90% of metastatic patients dying within less than 1 year from metastases to the liver. The intraarterial hepatic (iah) administration of the alkylating agent fotemustine holds some promise with response rates of 36% and median survival of 15 months. Here, we investigated whether the DNA-repair-protein MGMT may be involved in the variability of response to fotemustine and temozolomide in uveal melanoma. Epigenetic inactivation of MGMT has been demonstrated to be a predictive marker for benefit from alkylating agent therapy in glioblastoma. We found a methylated MGMT promoter in 6% of liver metastases from 34 uveal melanoma patients. The mean MGMT activity measured in liver metastases with negligible liver tissue content was significantly lower than in liver tissue (146 versus 523 fmol/mg protein, p = 0.002). Expression of the MGMT protein was detectable in 50% of 88 metastases by immunohistochemistry on a tissue microarray. Expression was heterogeneous, and in accordance with MGMT activity data, usually lower than in the surrounding liver. Differential MGMT activity/expression between metastasis and liver tissue and more efficient depletion of MGMT with higher doses of alkylating agent therapy using iah delivery may provide the pharmacologic window for the higher response rate. However, these results do not support MGMT methylation status or protein expression as predictive markers for treatment outcome to iah chemotherapy with alkylating agents.

Genome-wide analysis of cancer/testis gene expression.

Cancer/Testis (CT) genes, normally expressed in germ line cells but also activated in a wide range of cancer types, often encode antigens that are immunogenic in cancer patients, and present potential for use as biomarkers and targets for immunotherapy. Using multiple in silico gene expression analysis technologies, including twice the number of expressed sequence tags used in previous studies, we have performed a comprehensive genome-wide survey of expression for a set of 153 previously described CT genes in normal and cancer expression libraries. We find that although they are generally highly expressed in testis, these genes exhibit heterogeneous gene expression profiles, allowing their classification into testis-restricted (39), testis/brain-restricted (14), and a testis-selective (85) group of genes that show additional expression in somatic tissues. The chromosomal distribution of these genes confirmed the previously observed dominance of X chromosome location, with CT-X genes being significantly more testis-restricted than non-X CT. Applying this core classification in a genome-wide survey we identified >30 CT candidate genes; 3 of them, PEPP-2, OTOA, and AKAP4, were confirmed as testis-restricted or testis-selective using RT-PCR, with variable expression frequencies observed in a panel of cancer cell lines. Our classification provides an objective ranking for potential CT genes, which is useful in guiding further identification and characterization of these potentially important diagnostic and therapeutic targets.

Molecular role of Cx37 in advanced atherosclerosis: a micro-array study.

Recently, we showed that connexin37 (Cx37) protects against early atherosclerotic lesion development by regulating monocyte adhesion. The expression of this gap junction protein is altered in mouse and human atherosclerotic lesions; it is increased in macrophages newly recruited to the lesions and disappears from the endothelium of advanced plaques. To obtain more insight into the molecular role of Cx37 in advanced atherosclerosis, we used micro-array analysis for gene expression profiling in aortas of ApoE(-/-) and Cx37(-/-)ApoE(-/-) mice before and after 18 weeks of cholesterol-rich diet. Out of >15,000 genes, 106 genes were significantly differentially expressed in young mice before diet (P-value of <0.05, fold change of >0.7 or <-0.7, and intensity value >2.2 times background). Ingenuity pathway analysis (IPA) revealed differences in genes involved in cell-to-cell signaling and interaction, cellular compromise and nutritional disease. In addition, we identified 100 genes that were significantly perturbed after the cholesterol-rich diet. Similar to the analysis on 10-week-old mice, IPA revealed differences in genes involved in cell-to-cell signaling and interaction as well as to immuno-inflammatory disease. Furthermore, we found important changes in genes involved in vascular calcification and matrix degradation, some of which were confirmed at protein level by (immuno-)histochemistry. In conclusion, we suggest that Cx37 deficiency alters the global differential gene expression profiles in young mice towards a pro-inflammatory phenotype, which are then further influenced in advanced atherosclerosis. The results provide new insights into the significance of Cx37 in plaque calcification.

Test of four colon cancer risk-scores in formalin fixed paraffin embedded microarray gene expression data.

BACKGROUND: Prognosis prediction for resected primary colon cancer is based on the T-stage Node Metastasis (TNM) staging system. We investigated if four well-documented gene expression risk scores can improve patient stratification. METHODS: Microarray-based versions of risk-scores were applied to a large independent cohort of 688 stage II/III tumors from the PETACC-3 trial. Prognostic value for relapse-free survival (RFS), survival after relapse (SAR), and overall survival (OS) was assessed by regression analysis. To assess improvement over a reference, prognostic model was assessed with the area under curve (AUC) of receiver operating characteristic (ROC) curves. All statistical tests were two-sided, except the AUC increase. RESULTS: All four risk scores (RSs) showed a statistically significant association (single-test, P < .0167) with OS or RFS in univariate models, but with HRs below 1.38 per interquartile range. Three scores were predictors of shorter RFS, one of shorter SAR. Each RS could only marginally improve an RFS or OS model with the known factors T-stage, N-stage, and microsatellite instability (MSI) status (AUC gains < 0.025 units). The pairwise interscore discordance was never high (maximal Spearman correlation = 0.563) A combined score showed a trend to higher prognostic value and higher AUC increase for OS (HR = 1.74, 95% confidence interval [CI] = 1.44 to 2.10, P < .001, AUC from 0.6918 to 0.7321) and RFS (HR = 1.56, 95% CI = 1.33 to 1.84, P < .001, AUC from 0.6723 to 0.6945) than any single score. CONCLUSIONS: The four tested gene expression-based risk scores provide prognostic information but contribute only marginally to improving models based on established risk factors. A combination of the risk scores might provide more robust information. Predictors of RFS and SAR might need to be different.

Comparison of differential gene expression to water stress among bacteria with relevant pollutant-degradation properties.

Resistance to semi-dry environments has been considered a crucial trait for superior growth and survival of strains used for bioaugmentation in contaminated soils. In order to compare water stress programmes, we analyse differential gene expression among three phylogenetically different strains capable of aromatic compound degradation: Arthrobacter chlorophenolicus A6, Sphingomonas wittichii RW1 and Pseudomonas veronii 1YdBTEX2. Standardized laboratory-induced water stress was imposed by shock exposure of liquid cultures to water potential decrease, induced either by addition of solutes (NaCl, solute stress) or by addition of polyethylene glycol (matric stress), both at absolute similar stress magnitudes and at those causing approximately similar decrease of growth rates. Genome-wide differential gene expression was recorded by micro-array hybridizations. Growth of P. veronii 1YdBTEX2 was the most sensitive to water potential decrease, followed by S. wittichii RW1 and A. chlorophenolicus A6. The number of genes differentially expressed under decreasing water potential was lowest for A. chlorophenolicus A6, increasing with increasing magnitude of the stress, followed by S. wittichii RW1 and P. veronii 1YdBTEX2. Gene inspection and gene ontology analysis under stress conditions causing similar growth rate reduction indicated that common reactions among the three strains included diminished expression of flagellar motility and increased expression of compatible solutes (which were strain-specific). Furthermore, a set of common genes with ill-defined function was found between all strains, including ABC transporters and aldehyde dehydrogenases, which may constitute a core conserved response to water stress. The data further suggest that stronger reduction of growth rate of P. veronii 1YdBTEX2 under water stress may be an indirect result of the response demanding heavy NADPH investment, rather than the presence or absence of a suitable stress defence mechanism per se.

Genetic analysis of the porcine group B rotavirus NSP2 gene from wild-type Brazilian strains

Group B rotaviruses (RV-B) were first identified in piglet feces, being later associated with diarrhea in humans, cattle, lambs, and rats. In human beings, the virus was only described in China, India, and Bangladesh, especially infecting adults. Only a few studies concerning molecular analysis of the RV-B NSP2 gene have been conducted, and porcine RV-B has not been characterized. In the present study, three porcine wild-type RV-B strains from piglet stool samples collected from Brazilian pig herds were used for analysis. PAGE results were inconclusive for those samples, but specific amplicons of the RV-B NSP2 gene (segment 8) were obtained in a semi-nested PCR assay. The three porcine RV-B strains showed the highest nucleotide identity with the human WH1 strain and the alignments with other published sequences resulted in three groups of strains divided according to host species. The group of human strains showed 92.4 to 99.7% nucleotide identity while the porcine strains of the Brazilian RV-B group showed 90.4 to 91.8% identity to each other. The identity of the Brazilian porcine RV-B strains with outer sequences consisting of group A and C rotaviruses was only 35.3 to 38.8%. A dendrogram was also constructed to group the strains into clusters according to host species: human, rat, and a distinct third cluster consisting exclusively of the Brazilian porcine RV-B strains. This is the first study of the porcine RV-B NSP2 gene that contributes to the partial characterization of this virus and demonstrates the relationship among RV-B strains from different host species.

Down-regulation of the CSLF6 gene results in decreased (1,3;1,4)-beta-D-glucan in endosperm of wheat

(1,3;1,4)-beta-d-Glucan (beta-glucan) accounts for 20% of the total cell walls in the starchy endosperm of wheat (Triticum aestivum) and is an important source of dietary fiber for human nutrition with potential health benefits. Bioinformatic and array analyses of gene expression profiles in developing caryopses identified the CELLULOSE SYNTHASE-LIKE F6 (CSLF6) gene as encoding a putative beta-glucan synthase. RNA interference constructs were therefore designed to down-regulate CSLF6 gene expression and expressed in transgenic wheat under the control of a starchy endosperm-specific HMW subunit gene promoter. Analysis of wholemeal flours using an enzyme-based kit and by high-performance anion-exchange chromatography after digestion with lichenase showed decreases in total beta-glucan of between 30% and 52% and between 36% and 53%, respectively, in five transgenic lines compared to three control lines. The content of water-extractable beta-glucan was also reduced by about 50% in the transgenic lines, and the M(r) distribution of the fraction was decreased from an average of 79 to 85 x 10(4) g/mol in the controls and 36 to 57 x 10(4) g/mol in the transgenics. Immunolocalization of beta-glucan in semithin sections of mature and developing grains confirmed that the impact of the transgene was confined to the starchy endosperm with little or no effect on the aleurone or outer layers of the grain. The results confirm that the CSLF6 gene of wheat encodes a beta-glucan synthase and indicate that transgenic manipulation can be used to enhance the health benefits of wheat products.

GEDI: a user-friendly toolbox for analysis of large-scale gene expression data

Abstract Background Several mathematical and statistical methods have been proposed in the last few years to analyze microarray data. Most of those methods involve complicated formulas, and software implementations that require advanced computer programming skills. Researchers from other areas may experience difficulties when they attempting to use those methods in their research. Here we present an user-friendly toolbox which allows large-scale gene expression analysis to be carried out by biomedical researchers with limited programming skills. Results Here, we introduce an user-friendly toolbox called GEDI (Gene Expression Data Interpreter), an extensible, open-source, and freely-available tool that we believe will be useful to a wide range of laboratories, and to researchers with no background in Mathematics and Computer Science, allowing them to analyze their own data by applying both classical and advanced approaches developed and recently published by Fujita et al. Conclusion GEDI is an integrated user-friendly viewer that combines the state of the art SVR, DVAR and SVAR algorithms, previously developed by us. It facilitates the application of SVR, DVAR and SVAR, further than the mathematical formulas present in the corresponding publications, and allows one to better understand the results by means of available visualizations. Both running the statistical methods and visualizing the results are carried out within the graphical user interface, rendering these algorithms accessible to the broad community of researchers in Molecular Biology.

Exome analysis of HIV patients submitted to dendritic cells therapeutic vaccine reveals an association of CNOT1 gene with response to the treatment

INTRODUCTION: With the aim of searching genetic factors associated with the response to an immune treatment based on autologous monocyte-derived dendritic cells pulsed with autologous inactivated HIV, we performed exome analysis by screening more than 240,000 putative functional exonic variants in 18 HIV-positive Brazilian patients that underwent the immune treatment. METHODS: Exome analysis has been performed using the ILLUMINA Infinium HumanExome BeadChip. zCall algorithm allowed us to recall rare variants. Quality control and SNP-centred analysis were done with GenABEL R package. An in-house implementation of the Wang method permitted gene-centred analysis. RESULTS: CCR4-NOT transcription complex, subunit 1 (CNOT1) gene (16q21), showed the strongest association with the modification of the response to the therapeutic vaccine (p=0.00075). CNOT1 SNP rs7188697 A/G was significantly associated with DC treatment response. The presence of a G allele indicated poor response to the therapeutic vaccine (p=0.0031; OR=33.00; CI=1.74-624.66), and the SNP behaved in a dominant model (A/A vs. A/G+G/G p=0.0009; OR=107.66; 95% CI=3.85-3013.31), being the A/G genotype present only in weak/transient responders, conferring susceptibility to poor response to the immune treatment. DISCUSSION: CNOT1 is known to be involved in the control of mRNA deadenylation and mRNA decay. Moreover, CNOT1 has been recently described as being involved in the regulation of inflammatory processes mediated by tristetraprolin (TTP). The TTP-CCR4-NOT complex (CNOT1 in the CCR4-NOT complex is the binding site for TTP) has been reported as interfering with HIV replication, through post-transcriptional control. Therefore, we can hypothesize that genetic variation occurring in the CNOT1 gene could impair the TTP-CCR4-NOT complex, thus interfering with HIV replication and/or host immune response. CONCLUSIONS: Being aware that our findings are exclusive to the 18 patients studied with a need for replication, and that the genetic variant of CNOT1 gene, localized at intron 3, has no known functional effect, we propose a novel potential candidate locus for the modulation of the response to the immune treatment, and open a discussion on the necessity to consider the host genome as another potential variant to be evaluated when designing an immune therapy study