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.

Heat and chemical stress modulate the expression of the alpha-RYR gene in broiler chickens

The biological cause of Pork Stress syndrome, which leads to PSE (pale, soft, exudative) meat, is excessive release of Ca(2+) ions, which is promoted by a genetic mutation in the ryanodine receptors (RyR) located in the sarcoplasmic reticulum of the skeletal muscle cells. We examined the relationship between the formation of PSE meat under halothane treatment and heat stress exposure in chicken alpha RYR hot spot fragments. Four test groups were compared: 1) birds slaughtered without any treatment, i.e., the control group (C); 2) birds slaughtered immediately after halothane treatment (H); 3) birds slaughtered immediately after heat stress treatment (HS), and 4) birds exposed to halothane and to heat stress (H+HS), before slaughtering. Breast muscle mRNA was extracted, amplified by RT-PCR, and sequenced. PSE meat was evaluated using color determination (L*value). The most common alteration was deletion of a single nucleotide, which generated a premature stop codon, resulting in the production of truncated proteins. The highest incidence of nonsense transcripts came with exposure to halothane; 80% of these abnormal transcripts were detected in H and H+HS groups. As a consequence, the incidence of abnormal meat was highest in the H+HS group (66%). In HS, H, and C groups, PSE meat developed in 60, 50, and 33% of the samples, respectively. Thus, halothane apparently modulates alpha RYR gene expression in this region, and synergically with exposure to heat stress, causes Avian Stress syndrome, resulting in PSE meat in broiler chickens.

Chicken skeletal muscle-associated macroarray for gene discovery

Macro- and microarrays are well-established technologies to determine gene functions through repeated measurements of transcript abundance. We constructed a chicken skeletal muscle-associated array based on a muscle-specific EST database, which was used to generate a tissue expression dataset of similar to 4500 chicken genes across 5 adult tissues (skeletal muscle, heart, liver, brain, and skin). Only a small number of ESTs were sufficiently well characterized by BLAST searches to determine their probable cellular functions. Evidence of a particular tissue-characteristic expression can be considered an indication that the transcript is likely to be functionally significant. The skeletal muscle macroarray platform was first used to search for evidence of tissue-specific expression, focusing on the biological function of genes/transcripts, since gene expression profiles generated across tissues were found to be reliable and consistent. Hierarchical clustering analysis revealed consistent clustering among genes assigned to 'developmental growth', such as the ontology genes and germ layers. Accuracy of the expression data was supported by comparing information from known transcripts and tissue from which the transcript was derived with macroarray data. Hybridization assays resulted in consistent tissue expression profile, which will be useful to dissect tissue-regulatory networks and to predict functions of novel genes identified after extensive sequencing of the genomes of model organisms. Screening our skeletal-muscle platform using 5 chicken adult tissues allowed us identifying 43 'tissue-specific' transcripts, and 112 co-expressed uncharacterized transcripts with 62 putative motifs. This platform also represents an important tool for functional investigation of novel genes; to determine expression pattern according to developmental stages; to evaluate differences in muscular growth potential between chicken lines, and to identify tissue-specific genes.

Effects of trans-10, cis-12 conjugated linoleic acid on gene expression and lipid metabolism of adipose tissue of growing pigs

The objective of the present study was to determine the effects of trans-10, cis-12 conjugated linoleic acid (CLA) in adipose tissue explant cultures of growing pigs on the following responses: lipogenesis (measured as rate of C-14-labeled glucose incorporation over a subsequent 2-h incubation in the presence or absence of insulin), lipolysis (release of non-esterified fatty acid over a 2-h incubation in the presence or absence of isoproterenol), activities of lipogenic enzymes, and mRNA abundance of fatty acid synthase (FAS). Adipose tissue explants from nine growing pigs (78 +/- 3 kg) were cultured in 199 medium with insulin, dexamethasone and antibiotics for 4, 12, 24, and 48 h. The treatments were 1) control: 100 mu M polyvinyl alcohol (PVA); 2) pGH: 100 ng/mL porcine growth hormone (pGH) plus 100 mu M PVA; 3) CLA200: 200 mu M trans-10, cis-12 CLA; 4) CLA50: 50 mu M trans-10, cis-12 CLA, and 5) LA: 200 mu M linoleic acid. Fatty acids were added along with PVA (2: 1), respectively, for 24 h. Explants were collected after each culture period and assayed for lipogenesis. Transcripts of FAS mRNA were quantified by real-time RT-PCR after 24 and 48 h. Lipolysis and activities of FAS, glucose 6-phosphate dehydrogenase, 6-phosphogluconate dehydrogenase, and NADP-malate dehydrogenase were determined after 48 h. As expected, glucose incorporation was decreased (P < 0.05) in response to pGH treatment (positive control). LA had no effect on any parameter evaluated. Treatment with trans-10, cis-12 CLA decreased FAS activity (P < 0.05), but NADPH-generating enzymes were unaffected by treatments. Consistent with reduction in FAS activity, both lipid synthesis and FAS mRNA abundance were reduced with chronic CLA treatment, pGH increased baseline and stimulated lipolysis (P < 0.05) after 48 h of culture, while CLA treatment had no effect on non-esterified fatty acid release. Results of this study showed that trans-10, cis-12 CLA alters lipogenesis but has no effect on lipolysis in cultures of pig adipose tissue.

Validation of reference genes for gene expression studies in the honey bee, Apis mellifera, by quantitative real-time RT-PCR

For obtaining accurate and reliable gene expression results it is essential that quantitative real-time RT-PCR (qRT-PCR) data are normalized with appropriate reference genes. The current exponential increase in postgenomic studies on the honey bee, Apis mellifera, makes the standardization of qRT-PCR results an important task for ongoing community efforts. For this aim we selected four candidate reference genes (actin, ribosomal protein 49, elongation factor 1-alpha, tbp-association factor) and used three software-based approaches (geNorm, BestKeeper and NormFinder) to evaluate the suitability of these genes as endogenous controls. Their expression was examined during honey bee development, in different tissues, and after juvenile hormone exposure. Furthermore, the importance of choosing an appropriate reference gene was investigated for two developmentally regulated target genes. The results led us to consider all four candidate genes as suitable genes for normalization in A. mellifera. However, each condition evaluated in this study revealed a specific set of genes as the most appropriated ones.

Structural basis of GC-1 selectivity for thyroid hormone receptor isoforms

Background: Thyroid receptors, TRa and TR beta, are involved in important physiological functions such as metabolism, cholesterol level and heart activities. Whereas metabolism increase and cholesterol level lowering could be achieved by TR beta isoform activation, TRa activation affects heart rates. Therefore, beta-selective thyromimetics have been developed as promising drug-candidates for treatment of obesity and elevated cholesterol level. GC-1 [ 3,5-dimethyl-4-(4'-hydroxy- 3'-isopropylbenzyl)-phenoxy acetic acid] has ability to lower LDL cholesterol with 600-to 1400-fold more potency and approximately two-to threefold more efficacy than atorvastatin (Lipitor(C)) in studies in rats, mice and monkeys. Results: To investigate GC-1 specificity, we solved crystal structures and performed molecular dynamics simulations of both isoforms complexed with GC-1. Crystal structures reveal that, in TRa Arg228 is observed in multiple conformations, an effect triggered by the differences in the interactions between GC-1 and Ser277 or the corresponding asparagine (Asn331) of TR beta. The corresponding Arg282 of TR beta is observed in only one single stable conformation, interacting effectively with the ligand. Molecular dynamics support this model: our simulations show that the multiple conformations can be observed for the Arg228 in TR alpha, in which the ligand interacts either strongly with the ligand or with the Ser277 residue. In contrast, a single stable Arg282 conformation is observed for TR beta, in which it strongly interacts with both GC-1 and the Asn331. Conclusion: Our analysis suggests that the key factors for GC-1 selectivity are the presence of an oxyacetic acid ester oxygen and the absence of the amino group relative to T(3). These results shed light into the beta-selectivity of GC-1 and may assist the development of new compounds with potential as drug candidates to the treatment of hypercholesterolemia and obesity.

Gene expression down-regulation in CD90(+) prostate tumor-associated stromal cells involves potential organ-specific genes

Background: The prostate stroma is a key mediator of epithelial differentiation and development, and potentially plays a role in the initiation and progression of prostate cancer. The tumor-associated stroma is marked by increased expression of CD90/THYI. Isolation and characterization of these stromal cells could provide valuable insight into the biology of the tumor microenvironment. Methods: Prostate CD90(+) stromal fibromuscular cells from tumor specimens were isolated by cell-sorting and analyzed by DNA microarray. Dataset analysis was used to compare gene expression between histologically normal and tumor-associated stromal cells. For comparison, stromal cells were also isolated and analyzed from the urinary bladder. Results: The tumor-associated stromal cells were found to have decreased expression of genes involved in smooth muscle differentiation, and those detected in prostate but not bladder. Other differential expression between the stromal cell types included that of the CXC-chemokine genes. Conclusion: CD90(+) prostate tumor-associated stromal cells differed from their normal counterpart in expression of multiple genes, some of which are potentially involved in organ development.

Gene Expression Profiles in Radiation Workers Occupationally Exposed to Ionizing Radiation

Ionizing radiation OR) imposes risks to human health and the environment. IR at low doses and low (lose rates has the potency to initiate carcinogenesis. Genotoxic environmental agents such as IR trigger a cascade of signal transduction pathways for cellular protection. In this study, using cDNA microarray technique, we monitored the gene expression profiles in lymphocytes derived from radiation-ex posed individuals (radiation workers). Physical dosimetry records on these patients indicated that the absorbed dose ranged from 0.696 to 39.088 mSv. Gene expression analysis revealed statistically significant transcriptional changes in a total of 78 genes (21 up-regulated and 57 clown-regulated) involved in several biological processes such as ubiquitin cycle (UHRF2 and PIAS1), DNA repair (LIG3, XPA, ERCC5, RAD52, DCLRE1C), cell cycle regulation/proliferation (RHOA, CABLES2, TGFB2, IL16), and stress response (GSTP1, PPP2R5A, DUSP22). Some of the genes that showed altered expression profiles in this study call be used as biomarkers for monitoring the chronic low level exposure in humans. Additionally, alterations in gene expression patterns observed in chronically exposed radiation workers reinforces the need for defining the effective radiation dose that causes immediate genetic damage as well as the long-term effects on genomic instability, including cancer.

Unexpected Diversity of Cellular Immune Responses against Nef and Vif in HIV-1-Infected Patients Who Spontaneously Control Viral Replication

Background: HIV-1-infected individuals who spontaneously control viral replication represent an example of successful containment of the AIDS virus. Understanding the anti-viral immune responses in these individuals may help in vaccine design. However, immune responses against HIV-1 are normally analyzed using HIV-1 consensus B 15-mers that overlap by 11 amino acids. Unfortunately, this method may underestimate the real breadth of the cellular immune responses against the autologous sequence of the infecting virus. Methodology and Principal Findings: Here we compared cellular immune responses against nef and vif-encoded consensus B 15-mer peptides to responses against HLA class I-predicted minimal optimal epitopes from consensus B and autologous sequences in six patients who have controlled HIV-1 replication. Interestingly, our analysis revealed that three of our patients had broader cellular immune responses against HLA class I-predicted minimal optimal epitopes from either autologous viruses or from the HIV-1 consensus B sequence, when compared to responses against the 15-mer HIV-1 type B consensus peptides. Conclusion and Significance: This suggests that the cellular immune responses against HIV-1 in controller patients may be broader than we had previously anticipated.

The frequency of CD127(low) expressing CD4(+)CD25(high) T regulatory cells is inversely correlated with human T lymphotrophic virus type-1 (HTLV-1) proviral load in HTLV-1-infection and HTLV-1-associated myelopathy/tropical spastic paraparesis

Background: CD4(+)CD25(high) regulatory T (T(Reg)) cells modulate antigen-specific T cell responses, and can suppress anti-viral immunity. In HTLV-1 infection, a selective decrease in the function of T(Reg) cell mediated HTLV-1-tax inhibition of FOXP3 expression has been described. The purpose of this study was to assess the frequency and phenotype of T(Reg) cells in HTLV-1 asymptomatic carriers and in HTLV-1-associated neurological disease (HAM/TSP) patients, and to correlate with measures of T cell activation. Results: We were able to confirm that HTLV-1 drives activation, spontaneous IFN gamma production, and proliferation of CD4+ T cells. We also observed a significantly lower proportion of CTLA-4(+) T(Reg) cells (CD4(+)CD25(high) T cells) in subjects with HAM/TSP patients compared to healthy controls. Ki-67 expression was negatively correlated to the frequency of CTLA-4(+) T(Reg) cells in HAM/TSP only, although Ki-67 expression was inversely correlated with the percentage of CD127(low) T(Reg) cells in healthy control subjects. Finally, the proportion of CD127(low) T(Reg) cells correlated inversely with HTLV-1 proviral load. Conclusion: Taken together, the results suggest that T(Reg) cells may be subverted in HAM/TSP patients, which could explain the marked cellular activation, spontaneous cytokine production, and proliferation of CD4(+) T cells, in particular those expressing the CD25(high)CD127(low) phenotype. T(Reg) cells represent a potential target for therapeutic intervention for patients with HTLV-1-related neurological diseases.

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.

CYP2C19 and ABCB1 gene polymorphisms are differently distributed according to ethnicity in the Brazilian general population

Background: Recent studies have reported the clinical importance of CYP2C19 and ABCB1 polymorphisms in an individualized approach to clopidogrel treatment. The aims of this study were to evaluate the frequencies of CYP2C19 and ABCB1 polymorphisms and to identify the clopidogrel-predicted metabolic phenotypes according to ethnic groups in a sample of individuals representative of a highly admixtured population. Methods: One hundred and eighty-three Amerindians and 1,029 subjects of the general population of 4 regions of the country were included. Genotypes for the ABCB1c.C3435T (rs1045642), CYP2C19*2 (rs4244285), CYP2C19*3 (rs4986893), CYP2C19*4 (rs28399504), CYP2C19*5 (rs56337013), and CYP2C19*17 (rs12248560) polymorphisms were detected by polymerase chain reaction followed by high resolution melting analysis. The CYP2C19*3, CYP2C19*4 and CYP2C19*5 variants were genotyped in a subsample of subjects (300 samples randomly selected). Results: The CYP2C19*3 and CYP2C19*5 variant alleles were not detected and the CYP2C19*4 variant allele presented a frequency of 0.3%. The allelic frequencies for the ABCB1c.C3435T, CYP2C19*2 and CYP2C19*17 polymorphisms were differently distributed according to ethnicity: Amerindian (51.4%, 10.4%, 15.8%); Caucasian descent (43.2%, 16.9%, 18.0%); Mulatto (35.9%, 16.5%, 21.3%); and African descent (32.8%, 20.2%, 26.3%) individuals, respectively. As a result, self-referred ethnicity was able to predict significantly different clopidogrel-predicted metabolic phenotypes prevalence even for a highly admixtured population. Conclusion: Our findings indicate the existence of inter-ethnic differences in the ABCB1 and CYP2C19 variant allele frequencies in the Brazilian general population plus Amerindians. This information could help in stratifying individuals from this population regarding clopidogrel-predicted metabolic phenotypes and design more cost-effective programs towards individualization of clopidogrel therapy.

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.

Collagen V-induced nasal tolerance downregulates pulmonary collagen mRNA gene and TGF-beta expression in experimental systemic sclerosis

Background: The purpose of this study was to evaluate collagen deposition, mRNA collagen synthesis and TGFbeta expression in the lung tissue in an experimental model of scleroderma after collagen V-induced nasal tolerance. Methods: Female New Zealand rabbits (N = 12) were immunized with 1 mg/ml of collagen V in Freund's adjuvant (IM). After 150 days, six immunized animals were tolerated by nasal administration of collagen V ( 25 mu g/day) (IM-TOL) daily for 60 days. The collagen content was determined by morphometry, and mRNA expressions of types I, III and V collagen were determined by Real-time PCR. The TGF-beta expression was evaluated by immunostaining and quantified by point counting methods. To statistic analysis ANOVA with Bonferroni test were employed for multiple comparison when appropriate and the level of significance was determined to be p < 0.05. Results: IM-TOL, when compared to IM, showed significant reduction in total collagen content around the vessels (0.371 +/- 0.118 vs. 0.874 +/- 0.282, p < 0.001), bronchioles (0.294 +/- 0.139 vs. 0.646 +/- 0.172, p < 0.001) and in the septal interstitium (0.027 +/- 0.014 vs. 0.067 +/- 0.039, p = 0.026). The lung tissue of IM-TOL, when compared to IM, showed decreased immunostaining of types I, III and V collagen, reduced mRNA expression of types I (0.10 +/- 0.07 vs. 1.0 +/- 0.528, p = 0.002) and V (1.12 +/- 0.42 vs. 4.74 +/- 2.25, p = 0.009) collagen, in addition to decreased TGF-beta expression ( p < 0.0001). Conclusions: Collagen V-induced nasal tolerance in the experimental model of SSc regulated the pulmonary remodeling process, inhibiting collagen deposition and collagen I and V mRNA synthesis. Additionally, it decreased TGF-beta expression, suggesting a promising therapeutic option for scleroderma treatment.

Differential gene expression profiling in mucus glands of honey bee (Apis mellifera) drones during sexual maturation

The mating sign that each drone leaves when mating with a queen essentially consists of mucus gland proteins. We employed a Representational Difference Analysis (RDA) methodology to identify genes that are differentially expressed in mucus glands during sexual maturation of drones. The RDA library for mucus glands of newly emerged drones was more complex than that of 8 day-old drones, with matches to 20 predicted genes. Another 26 reads matched to the Apis genome but not to any predicted gene. Since these ESTs were located within ORFs they may represent novel honey bee genes, possibly fast evolving mucus gland proteins. In the RDA library for mucus glands of 8 day-old drones, most reads corresponded to a capsid protein of deformed wing virus, indicating high viral loads in these glands. The expression of two genes encoding venom allergens, acid phosphatase-1 and hyaluronidase, in drone mucus glands argues for their homology with the female venom glands, both associated with the reproductive system.