A Mixture model with random-effects components for clustering correlated gene-expression profiles

Motivation: The clustering of gene profiles across some experimental conditions of interest contributes significantly to the elucidation of unknown gene function, the validation of gene discoveries and the interpretation of biological processes. However, this clustering problem is not straightforward as the profiles of the genes are not all independently distributed and the expression levels may have been obtained from an experimental design involving replicated arrays. Ignoring the dependence between the gene profiles and the structure of the replicated data can result in important sources of variability in the experiments being overlooked in the analysis, with the consequent possibility of misleading inferences being made. We propose a random-effects model that provides a unified approach to the clustering of genes with correlated expression levels measured in a wide variety of experimental situations. Our model is an extension of the normal mixture model to account for the correlations between the gene profiles and to enable covariate information to be incorporated into the clustering process. Hence the model is applicable to longitudinal studies with or without replication, for example, time-course experiments by using time as a covariate, and to cross-sectional experiments by using categorical covariates to represent the different experimental classes. Results: We show that our random-effects model can be fitted by maximum likelihood via the EM algorithm for which the E(expectation) and M(maximization) steps can be implemented in closed form. Hence our model can be fitted deterministically without the need for time-consuming Monte Carlo approximations. The effectiveness of our model-based procedure for the clustering of correlated gene profiles is demonstrated on three real datasets, representing typical microarray experimental designs, covering time-course, repeated-measurement and cross-sectional data. In these examples, relevant clusters of the genes are obtained, which are supported by existing gene-function annotation. A synthetic dataset is considered too.

The full Bayesian significance test for mixture models: results in gene expression clustering

Gene clustering is a useful exploratory technique to group together genes with similar expression levels under distinct cell cycle phases or distinct conditions. It helps the biologist to identify potentially meaningful relationships between genes. In this study, we propose a clustering method based on multivariate normal mixture models, where the number of clusters is predicted via sequential hypothesis tests: at each step, the method considers a mixture model of m components (m = 2 in the first step) and tests if in fact it should be m - 1. If the hypothesis is rejected, m is increased and a new test is carried out. The method continues (increasing m) until the hypothesis is accepted. The theoretical core of the method is the full Bayesian significance test, an intuitive Bayesian approach, which needs no model complexity penalization nor positive probabilities for sharp hypotheses. Numerical experiments were based on a cDNA microarray dataset consisting of expression levels of 205 genes belonging to four functional categories, for 10 distinct strains of Saccharomyces cerevisiae. To analyze the method's sensitivity to data dimension, we performed principal components analysis on the original dataset and predicted the number of classes using 2 to 10 principal components. Compared to Mclust (model-based clustering), our method shows more consistent results.

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.

The C242T polymorphism of the p22-phox gene (CYBA) is associated with higher left ventricular mass in Brazilian hypertensive patients

Background: Reactive oxygen species have been implicated in the physiopathogenesis of hypertensive end-organ damage. This study investigated the impact of the C242T polymorphism of the p22-phox gene (CYBA) on left ventricular structure in Brazilian hypertensive subjects. Methods: We cross-sectionally evaluated 561 patients from 2 independent centers [Campinas (n = 441) and Vitoria (n = 120)] by clinical history, physical examination, anthropometry, analysis of metabolic and echocardiography parameters as well as p22-phox C242T polymorphism genotyping. In addition, NADPH-oxidase activity was quantified in peripheral mononuclear cells from a subgroup of Campinas sample. Results: Genotype frequencies in both samples were consistent with the Hardy-Weinberg equilibrium. Subjects with the T allele presented higher left ventricular mass/height(2.7) than those carrying the CC genotype in Campinas (76.8 +/- 1.6 vs 70.9 +/- 1.4 g/m(2.7); p = 0.009), and in Vitoria (45.6 +/- 1.9 vs 39.9 +/- 1.4 g/m(2.7); p = 0.023) samples. These results were confirmed by stepwise regression analyses adjusted for age, gender, blood pressure, metabolic variables and use of anti-hypertensive medications. In addition, increased NADPH-oxidase activity was detected in peripheral mononuclear cells from T allele carriers compared with CC genotype carriers (p = 0.03). Conclusions: The T allele of the p22-phox C242T polymorphism is associated with higher left ventricular mass/height(2.7) and increased NADPH-oxidase activity in Brazilian hypertensive patients. These data suggest that genetic variation within NADPH-oxidase components may modulate left ventricular remodeling in subjects with systemic hypertension.

Evaluating gene by sex and age interactions on cardiovascular risk factors in Brazilian families

Background: In family studies, it is important to evaluate the impact of genes and environmental factors on traits of interest. In particular, the relative influences of both genes and the environment may vary in different strata of the population of interest, such as young and old individuals, or males and females. Methods: In this paper, extensions of the variance components model are used to evaluate heterogeneity in the genetic and environmental variance components due to the effects of sex and age (the cutoff between young and old was 43 yrs). The data analyzed were from 81 Brazilian families (1,675 individuals) of the Baependi Family Heart Study. Results: The models allowing for heterogeneity of variance components by sex suggest that genetic and environmental variances are not different in males and females for diastolic blood pressure, LDL-cholesterol, and HDL-cholesterol, independent of the covariates included in the models. However, for systolic blood pressure, fasting glucose and triglycerides, the evidence for heterogeneity was dependent on the covariates in the model. For instance, in the presence of sex and age covariates, heterogeneity in the genetic variance component was suggested for fasting glucose. But, for systolic blood pressure, there was no evidence of heterogeneity in any of the two variance components. Except for the LDL-cholesterol, models allowing for heterogeneity by age provide evidence of heterogeneity in genetic variance for triglycerides and systolic and diastolic blood pressure. There was evidence of heterogeneity in environmental variance in fasting glucose and HDL-cholesterol. Conclusions: Our results suggest that heterogeneity in trait variances should not be ignored in the design and analyses of gene-finding studies involving these traits, as it may generate additional information about gene effects, and allow the investigation of more sophisticated models such as the model including sex-specific oligogenic variance components.

Selective Decrease of Components of the Creatine Kinase System and ATP Synthase Complex in Chronic Chagas Disease Cardiomyopathy

Background: Chronic Chagas disease cardiomyopathy (CCC) is an inflammatory dilated cardiomyopathy with a worse prognosis than other cardiomyopathies. CCC occurs in 30 % of individuals infected with Trypanosoma cruzi, endemic in Latin America. Heart failure is associated with impaired energy metabolism, which may be correlated to contractile dysfunction. We thus analyzed the myocardial gene and protein expression, as well as activity, of key mitochondrial enzymes related to ATP production, in myocardial samples of end-stage CCC, idiopathic dilated (IDC) and ischemic (IC) cardiomyopathies. Methodology/Principal Findings: Myocardium homogenates from CCC (N = 5), IC (N = 5) and IDC (N = 5) patients, as well as from heart donors (N = 5) were analyzed for protein and mRNA expression of mitochondrial creatine kinase (CKMit) and muscular creatine kinase (CKM) and ATP synthase subunits aplha and beta by immunoblotting and by real-time RT-PCR. Total myocardial CK activity was also assessed. Protein levels of CKM and CK activity were reduced in all three cardiomyopathy groups. However, total CK activity, as well as ATP synthase alpha chain protein levels, were significantly lower in CCC samples than IC and IDC samples. CCC myocardium displayed selective reduction of protein levels and activity of enzymes crucial for maintaining cytoplasmic ATP levels. Conclusions/Significance: The selective impairment of the CK system may be associated to the loss of inotropic reserve observed in CCC. Reduction of ATP synthase alpha levels is consistent with a decrease in myocardial ATP generation through oxidative phosphorylation. Together, these results suggest that the energetic deficit is more intense in the myocardium of CCC patients than in the other tested dilated cardiomyopathies.

Intranasal vaccination with messenger RNA as a new approach in gene therapy: Use against tuberculosis

Background: mRNAs are highly versatile, non-toxic molecules that are easy to produce and store, which can allow transient protein expression in all cell types. The safety aspects of mRNA-based treatments in gene therapy make this molecule one of the most promising active components of therapeutic or prophylactic methods. The use of mRNA as strategy for the stimulation of the immune system has been used mainly in current strategies for the cancer treatment but until now no one tested this molecule as vaccine for infectious disease. Results: We produce messenger RNA of Hsp65 protein from Mycobacterium leprae and show that vaccination of mice with a single dose of 10 mu g of naked mRNA-Hsp65 through intranasal route was able to induce protection against subsequent challenge with virulent strain of Mycobacterium tuberculosis. Moreover it was shown that this immunization was associated with specific production of IL-10 and TNF-alpha in spleen. In order to determine if antigen presenting cells (APCs) present in the lung are capable of capture the mRNA, labeled mRNA-Hsp65 was administered by intranasal route and lung APCs were analyzed by flow cytometry. These experiments showed that after 30 minutes until 8 hours the populations of CD11c(+), CD11b(+) and CD19(+) cells were able to capture the mRNA. We also demonstrated in vitro that mRNA-Hsp65 leads nitric oxide (NO) production through Toll-like receptor 7 (TLR7). Conclusions: Taken together, our results showed a novel and efficient strategy to control experimental tuberculosis, besides opening novel perspectives for the use of mRNA in vaccines against infectious diseases and clarifying the mechanisms involved in the disease protection we noticed as well.

Probing the hormonal activity of fractionated molecular humic components in tomato auxin mutants

We progressively reduced the complexity of humic matter by a mild sequential removal of unbound or free components, weakly, and strongly bound molecules. The auxin-like response of residues from each step was tested using tomato (cv. Micro-Tom) seedlings expressing DR5 auxin synthetic promoter fused to the beta-glucuronidase (GUS) reporter gene and the low auxin-sensitivity diageotropica (dgt) mutant. Both exogenous auxin and humic matter promoted lateral root emergence in the control, but failed to induce lateral roots in the dgt mutant. When strongly bound components were removed from humic matter by breaking the ester and ether bonds, the humic residues lost their ability to induce the DR5

Transcription regulation of the PbGP43 gene by nitrogen in the human pathogen Paracoccidioides brasiliensis

We show indirect evidences for the possible involvement of NIT-2-like binding motifs in transcription modulation of the PbGP43 gene, which codes for an important antigen from the human fungal pathogen Paracoccidioides brasiliensis. This investigation was motivated by the finding of 23 NIT2-like sites within the proximal -2047 nucleotides of the PbGP43 5` intergenic region from the Pb339 isolate. They compose four clusters, two of them identical. We found four NIT2-containing probes that were positive in electrophoretic mobility shift assays and further analyzed them. PbGP43 could be modulated by nitrogen primary sources in Pb339, Pb3 and Pb18 isolates, as observed by reverse transcription (RT) real time-PCR. Gene reporter assays conducted in Aspergillus nidulans suggested that the minimal fragment responsible for nitrogen modulation lies within -480 bp of the PbGP43 gene. This is the first report on PbGP43 transcription modulation in response to nitrogen primary sources, which might help understand its regulation during infection. (C) 2008 Elsevier Inc. All rights reserved.

Identification of novel "pathologs" (human disease-related gene candidates) from the RIKEN full-length mouse cDNA data set

The majority of common diseases such as cancer, allergy, diabetes, or heart disease are characterized by complex genetic traits, in which genetic and environmental components contribute to disease susceptibility. Our knowledge of the genetic factors underlying most of such diseases is limited. A major goal in the post-genomic era is to identify and characterize disease susceptibility genes and to use this knowledge for disease treatment and prevention. More than 500 genes are conserved across the invertebrate and vertebrate genomes. Because of gene conservation, various organisms including yeast, fruitfly, zebrafish, rat, and mouse have been used as genetic models.

Genetic interactions of the Ashergillus nidulans atmA(ATM) homolog with different components of the DNA damage response pathway

Ataxia telangiectasia mutated (ATM) is a phosphatidyl-3-kinase-related protein kinase that functions as a central regulator of the DNA damage response in eukaryotic cells. In humans, mutations in ATM cause the devastating neurodegenerative disease ataxia telangiectasia. Previously, we characterized the homolog of ATM (AtmA) in the filamentous fungus Aspergillus nidulans. In addition to its expected role in the DNA damage response, we found that AtmA is also required for polarized hyphal growth. Here, we extended these studies by investigating which components of the DNA damage response pathway are interacting with AtmA. The AtmA(ATM) loss of function caused synthetic lethality when combined with mutation in UvsB(ATR). Our results suggest that AtmA and UvsB are interacting and they are probably partially redundant in terms of DNA damage sensing and/or repairing and polar growth. We identified and inactivated A. nidulans chkA(CHK1) and chkB(CHK2) genes. These genes are also redundantly involved in A. nidulans DNA damage response. We constructed several combinations of double mutants for Delta atmA, Delta uvsB, Delta chkA, and Delta chkB. We observed a complex genetic relationship with these mutations during the DNA replication checkpoint and DNA damage response. Finally, we observed epistatic and synergistic interactions between AtmA, and bimE(APCI), ankA(WEE1) and the cdc2-related kinase npkA, at S-phase checkpoint and in response to DNA-damaging agents.

Molecular profiling of isolated histological components of Wilms tumor implicates a common role for the Wnt signaling pathway in kidney and tumor development

Wilms tumor (WT), a tumor composed of three histological components - blastema (BL), epithelia and stroma - is considered an appropriate model system to study the biological relationship between differentiation and tumorigenesis. To investigate molecular associations between nephrogenesis and WT, the gene expression pattern of individual cellular components was analyzed, using a customized platform containing 4,608 genes. WT gene expression patterns were compared to genes regulated during kidney differentiation. BL had a closer gene expression pattern to the earliest stage of normal renal development. The BL gene expression pattern was compared to that of fetal kidney (FK) and also between FK and mature kidney, identifying 25 common de-regulated genes supposedly involved in the earliest events of WT onset. Quantitative RT-PCR was performed, confirming the difference in expression levels for 13 of 16 genes (81.2%) in the initial set and 8 of 13 (61.5%) in an independent set of samples. An overrepresentation of genes belonging to the Wnt signaling pathway was identified, namely PLCG2, ROCK2 and adenomatous polyposis coli (APC). Activation of the Wnt pathway was confirmed in WT, using APC at protein level and PLCG2 at mRNA and protein level. APC showed positive nuclear immunostaining for an independent set of WT samples, similarly to the FK in week 11. Lack of PLCG2 expression was confirmed in WT and in FK until week 18. Taken together, these results provided molecular evidence of the recapitulation of the embryonic kidney by WT as well as involvement of the Wnt pathway in the earliest events of WT onset. Copyright (C) 2008 S. Karger AG, Basel.

T allele of-344C/T polymorphism in aldosterone synthase gene is not associated with resistant hypertension

Resistant hypertension (RH) is the maintenance of elevated blood pressure concurrent with the use of three different antihypertensive drugs, one of which is a diuretic. The Renin-Angiotensin-Aldosterone System plays a major role in volume-dependent hypertension. Therefore, its components are interesting targets for genetic association studies. This work focused on the -344 C/T polymorphism in the CYP11b2 gene, which encodes aldosterone synthase. This work evaluates the association between T allele and resistance to anti-hypertensive treatment. Genotyping analysis included 88 subjects with RH, 142 who were responsive to anti-hypertensive treatment and 110 subjects as a control group. Plasmatic concentrations of aldosterone, renin and cortisol, carotid intima-media thickness and carotid-femoral pulse wave velocity were assessed in a smaller subset of hypertensive patients. An association was found between T allele and hypertension (P < 0.005), but there was no difference in allele frequencies between both hypertensive groups. There was no difference in plasmatic parameters either, in remodeling indicators between the genotypic groups.

Rhipicephalus (Boophilus) microplus: Clotting time in tick-infested skin varies according to local inflammation and gene expression patterns in tick salivary glands

Ticks deposit saliva at the site of their attachment to a host in order to inhibit haemostasis, inflammation and innate and adaptive immune responses. The anti-haemostatic properties of tick saliva have been described by many studies, but few show that tick infestations or its anti-haemostatic components exert systemic effects in vivo. In the present study, we extended these observations and show that, compared with normal skin, bovine hosts that are genetically susceptible to tick infestations present an increase in the clotting time of blood collected from the immediate vicinity of haemorrhagic feeding pools in skin infested with different developmental stages of Rhipicepahlus microplus; conversely, we determined that clotting time of tick-infested skin from genetically resistant bovines was shorter than that of normal skin. Coagulation and inflammation have many components in common and we determined that in resistant bovines, eosinophils and basophils, which are known to contain tissue factor, are recruited in greater numbers to the inflammatory site of tick bites than in susceptible hosts. Finally, we correlated the observed differences in clotting times with the expression profiles of transcripts for putative anti-haemostatic proteins in different developmental stages of R. microplus fed on genetically susceptible and resistant hosts: we determined that transcripts coding for proteins similar to these molecules are overrepresented in salivary glands from nymphs and males fed on susceptible bovines. Our data indicate that ticks are able to modulate their host`s local haemostatic reactions. In the resistant phenotype, larger amounts of inflammatory cells are recruited and expression of anti-coagulant molecules is decreased tick salivary glands, features that can hamper the tick`s blood meal. (C) 2010 Elsevier Inc. All rights reserved.

Lp95, a novel leptospiral protein that binds extracellular matrix components and activates e-selectin on endothelial cells

Objectives: The study of a predicted outer membrane leptospiral protein encoded by the gene LIC12690 in mediating the adhesion process. Methods: The gene was cloned and expressed in Escherichia coli BL21 (SI) strain by using the expression vector pAE. The recombinant protein tagged with N-terminal hexahistidine was purified by metal-charged chromatography and used to assess its ability to activate human umbilical vein endothelial cells (HUVECs). Results: The recombinant leptospiral protein of 95 kDa, named Lp95, activated E-selectin in a dose-dependent fashion but not the intercellular adhesion molecule 1 (ICAM-1). In addition, we show that pathogenic and non-pathogenic Leptospira are both capable to stimulate endothelium E-selectin and ICAM-1, but the pathogenic L. interrogans serovar Copenhageni strain promotes a statistically significant higher activation than the non-pathogenic L. biflexa serovar Patoc (P < 0.01). The Lp95 was identified in vivo in the renal tubules of animal during experimental infection with L. interrogans. The whole Lp95 as well as its fragments, the C-terminal containing the domain of unknown function (DUF), the N-terminal and the central overlap regions bind laminin and fibronectin ECM molecules, being the binding stronger with the DUF containing fragment. Conclusion: This is the first leptospiral protein capable to mediate the adhesion to ECM components and the activation of HUVECS, thus suggesting its participation in the pathogenesis of Leptospira. (C) 2009 The British Infection Society. Published by Elsevier Ltd. All rights reserved.