Bayesian model accounting for within-class biological variability in Serial Analysis of Gene Expression (SAGE)

Abstract Background An important challenge for transcript counting methods such as Serial Analysis of Gene Expression (SAGE), "Digital Northern" or Massively Parallel Signature Sequencing (MPSS), is to carry out statistical analyses that account for the within-class variability, i.e., variability due to the intrinsic biological differences among sampled individuals of the same class, and not only variability due to technical sampling error. Results We introduce a Bayesian model that accounts for the within-class variability by means of mixture distribution. We show that the previously available approaches of aggregation in pools ("pseudo-libraries") and the Beta-Binomial model, are particular cases of the mixture model. We illustrate our method with a brain tumor vs. normal comparison using SAGE data from public databases. We show examples of tags regarded as differentially expressed with high significance if the within-class variability is ignored, but clearly not so significant if one accounts for it. Conclusion Using available information about biological replicates, one can transform a list of candidate transcripts showing differential expression to a more reliable one. Our method is freely available, under GPL/GNU copyleft, through a user friendly web-based on-line tool or as R language scripts at supplemental web-site.

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

Abstract 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/THY1. 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.

Estradiol induces transcriptional and posttranscriptional modifications in versican expression in the mouse uterus

We have previously shown the differential expression of versican in the mouse uterus under ovarian hormone influence. We also demonstrated there is not a direct correlation between mRNA levels and protein expression, suggesting posttranscriptional events, such as alteration in mRNA stability. This posttranscriptional effect may result in the elongation and stabilization of transcripts poly(A) tail. Thus, the aim of this study was to analyze whether estradiol (E2) regulates versican mRNA stability and expression in a dose-related and time-dependent manner. For this purpose female mice were ovariectomized and treated with a single injection of 0.1 or 10 μg E2. To block transcription a group of females received a single injection of alpha-amanitin before hormone administration. Uterine tissues were collected 30 min, 1, 3, 6, 12 and 24 h after treatments and processed for quantitative real time PCR (qPCR), RACE-PAT Assay and immunohistochemistry. qPCR showed that versican mRNA levels are higher than control from 3 to 24 h after E2 administration, whereas after transcription inhibition versican mRNA unexpectedly increases within 3 h, which can be explained when transcriptional blockers alter the degradation rate of the transcript, resulting in the superinduction of this mRNA. Accordingly, analysis of versican transcript poly(A) tail evidenced a longer product 3 h after treatment, but not after 12 h. Versican immunoreaction becomes conspicuous in the superficial stroma only 3 h after E2 injection, whereas the whole stroma is immunoreactive from 6 h onward. These results demonstrate that E2 modulates versican at the transcriptional and posttranscriptional levels in a time-dependent manner.

A novel phase variation mechanism in the meningococcus driven by a ligand-responsive repressor and differential spacing of distal promoter elements

Phase variable expression, mediated by high frequency reversible changes in the length of simple sequence repeats, facilitates adaptation of bacterial populations to changing environments and is frequently important in bacterial virulence. Here we elucidate a novel phase variable mechanism for NadA expression, an adhesin and invasin of Neisseria meningitidis. The NadR repressor protein binds to operators flanking the phase variable tract of the nadA promoter gene and contributes to the differential expression levels of phase variant promoters with different numbers of repeats, likely due to different spacing between operators. It is shown that IHF binds between these operators, and may permit looping of the promoter, allowing interaction of NadR at operators located distally or overlapping the promoter. The 4-hydroxyphenylacetic acid, a metabolite of aromatic amino acid catabolism that is secreted in saliva, induces nadA expression by inhibiting the DNA binding activity of the NadR repressor. When induced, only minor differences are evident between NadR-independent transcription levels of promoter phase variants, which are likely due to differential RNA polymerase contacts leading to altered promoter activity. These results suggest that NadA expression is under both stochastic and tight environmental-sensing regulatory control, and both regulations are mediated by the NadR repressor that and may be induced during colonization of the oropharynx where it plays a major role in the successful adhesion and invasion of the mucosa. Hence, simple sequence repeats in promoter regions may be a strategy used by host-adapted bacterial pathogens to randomly switch between expression states that may nonetheless still be induced by appropriate niche-specific signals.

Gene expression analysis in ‘Candidatus Phytoplasma mali’-resistant and -susceptible Malus genotypes

Apple proliferation (AP) disease is the most important graft-transmissible and vector-borne disease of apple in Europe. ‘Candidatus Phytoplasma mali’ (Ca. P. mali) is the causal agent of AP. Apple (Malus x domestica) and other Malus species are the only known woody hosts. In European apple orchards, the cultivars are mainly grafted on one rootstock, M. x domestica cv. M9. M9 like all other M. x domestica cultivars is susceptible to ‘Ca. P. mali’. Resistance to AP was found in the wild genotype Malus sieboldii (MS) and in MS-derived hybrids but they were characterised by poor agronomic value. The breeding of a new rootstock carrying the resistant and the agronomic traits was the major aim of a project of which this work is a part. The objective was to shed light into the unknown resistance mechanism. The plant-phytoplasma interaction was studied by analysing differences between the ‘Ca. P. mali’-resistant and -susceptible genotypes related to constitutively expressed genes or to induced genes during infection. The cDNA-Amplified Fragment Length Polymorphism (cDNA-AFLP) technique was employed in both approaches. Differences related to constitutively expressed genes were identified between two ‘Ca. P. mali’-resistant hybrid genotypes (4551 and H0909) and the ‘Ca. P. mali’-susceptible M9. 232 cDNA-AFLP bands present in the two resistant genotypes but absent in the susceptible one were isolated but several different products associated to each band were found. Therefore, two different macroarray hybridisation experiments were performed with the cDNA-AFLP fragments yielding 40 sequences encoding for genes of unknown function or a wide array of functions including plant defence. In the second approach, individuation and analysis of the induced genes was carried out exploiting an in vitro system in which healthy and ‘Ca. P. mali’-infected micropropagated plants were maintained under controlled conditions. Infection trials using in vitro grafting of ‘Ca. P. mali’ showed that the resistance phenotype could be reproduced in this system. In addition, ex vitro plants were generated as an independent control of the genes differentially expressed in the in vitro plants. The cDNA-AFLP analysis in in vitro plants yielded 63 bands characterised by over-expression in the infected state of both the H0909 and MS genotypes. The major part (37 %) of the associated sequences showed homology with products of unknown function. The other genes were involved in plant defence, energy transport/oxidative stress response, protein metabolism and cellular growth. Real-time qPCR analysis was employed to validate the differential expression of the genes individuated in the cDNA-AFLP analysis. Since no internal controls were available for the study of the gene expression in Malus, an analysis on housekeeping genes was performed. The most stably expressed genes were the elongation factor-1 α (EF1) and the eukaryotic translation initiation factor 4-A (eIF4A). Twelve out of 20 genes investigated through qPCR were significantly differentially expressed in at least one genotype either in in vitro plants or in ex vitro plants. Overall, about 20% of the genes confirmed their cDNA-AFLP expression pattern in M. sieboldii or H0909. On the contrary, 30 % of the genes showed down-regulation or were not differentially expressed. For the remaining 50 % of the genes a contrasting behaviour was observed. The qPCR data could be interpreted as follows: the phytoplasma infection unbalance photosynthetic activity and photorespiration down-regulating genes involved in photosynthesis and in the electron transfer chain. As result, and in contrast to M. x domestica genotypes, an up-regulation of genes of the general response against pathogens was found in MS. These genes involved the pathway of H2O2 and the production of secondary metabolites leading to the hypothesis that a response based on the accumulation of H2O2 in MS would be at the base of its resistance. This resembles a phenomenon known as “recovery” where the spontaneous remission of the symptoms is observed in old susceptible plants but occurring in a stochastic way while the resistance in MS is an inducible but stable feature. As additional product of this work three cDNA-AFLP-derived markers were developed which showed independent distribution among the seedlings of two breeding progenies and were associated to a genomic region characteristic of MS. These markers will contribute to the development of molecular markers for the resistance as well as to map the resistance on the Malus genome.

Expression profiling of stem cell-related genes in neoadjuvant-treated gastric cancer: a NOTCH2, GSK3B and β-catenin gene signature predicts survival

Cancer stem cell (CSC) based gene expression signatures are associated with prognosis in various tumour types and CSCs are suggested to be particularly drug resistant. The aim of our study was first, to determine the prognostic significance of CSC-related gene expression in residual tumour cells of neoadjuvant-treated gastric cancer (GC) patients. Second, we wished to examine, whether expression alterations between pre- and post-therapeutic tumour samples exist, consistent with an enrichment of drug resistant tumour cells. The expression of 44 genes was analysed in 63 formalin-fixed, paraffin embedded tumour specimens with partial tumour regression (10-50% residual tumour) after neoadjuvant chemotherapy by quantitative real time PCR low-density arrays. A signature of combined GSK3B(high), β-catenin (CTNNB1)(high) and NOTCH2(low) expression was strongly correlated with better patient survival (p<0.001). A prognostic relevance of these genes was also found analysing publically available gene expression data. The expression of 9 genes was compared between pre-therapeutic biopsies and post-therapeutic resected specimens. A significant post-therapeutic increase in NOTCH2, LGR5 and POU5F1 expression was found in tumours with different tumour regression grades. No significant alterations were observed for GSK3B and CTNNB1. Immunohistochemical analysis demonstrated a chemotherapy-associated increase in the intensity of NOTCH2 staining, but not in the percentage of NOTCH2. Taken together, the GSK3B, CTNNB1 and NOTCH2 expression signature is a novel, promising prognostic parameter for GC. The results of the differential expression analysis indicate a prominent role for NOTCH2 and chemotherapy resistance in GC, which seems to be related to an effect of the drugs on NOTCH2 expression rather than to an enrichment of NOTCH2 expressing tumour cells.

Expression of Toll-like receptor 2 in duodenal biopsies from dogs with inflammatory bowel disease is associated with severity of disease

There is growing evidence that aberrant innate immune responses towards the bacterial flora of the gut play a role in the pathogenesis of canine inflammatory bowel disease (IBD). Toll-like receptors (TLR) play an important role as primary sensors of invading pathogens and have gained significant attention in human IBD as differential expression and polymorphisms of certain TLR have been shown to occur in ulcerative colitis (UC) and Crohn's disease (CD). The aim of the current study was to evaluate the expression of two TLR important for recognition of commensals in the gut. TLR2 and TLR4 mRNA expression in duodenal biopsies from dogs with IBD was measured and correlated with clinical and histological disease severity. Endoscopic duodenal biopsies from 20 clinical cases and 7 healthy control dogs were used to extract mRNA. TLR2 and TLR4 mRNA expression was assessed using quantitative real-time PCR. TLR2 mRNA expression was significantly increased in the IBD dogs compared to controls, whereas TLR4 mRNA expression was similar in IBD and control cases. In addition, TLR2 mRNA expression was mildly correlated with clinical severity of disease, however, there was no correlation between TLR2 expression and histological severity of disease.

Aberrant low expression level of bovine beta-lactoglobulin is associated with a C to A transversion in the BLG promoter region

Beta-lactoglobulin (beta-LG) is the major whey protein in cow's milk. It is well established that the predominant 2 genetic variants, beta-LG A and B, are differentially expressed. Extensive investigation of the genetic variation in the promoter region of the BLG gene revealed the existence of specific haplotypes associated with the A and B variants, respectively. However, the genetic basis for the differential expression of BLG A and B alleles is still elusive. We have previously reported a quantitative beta-LG B variant, characterized by a very low beta-LG protein expression level. Here, we report that the corresponding BLG allele (BLG B*) shows a correspondingly low mRNA expression level. Comparative DNA sequencing of 7,670 bp of the BLG B* allele and the established BLG B allele revealed a unique difference of a C to A transversion at position 215 bp upstream of the translation initiation site (g.-215C>A). This mutation segregated perfectly with the differential phenotypic expression in a paternal half-sib family and could be confirmed in 2 independent cases. The sequence of the BLG B allele in the region of the mutation is highly conserved among 4 related ruminant species. The site of the mutation corresponds to a putative consensus-binding sequence for the transcription factors c-Rel and Elk-1 as predicted by searching the TRANSFAC database. The beta-LG B* site might be relevant in the natural production of milk of low beta-LG content.

Reduced atrial connexin43 expression after pediatric heart surgery

Myocardial dysfunction and arrhythmias may be induced by congenital heart defects, but also be the result of heart surgery with cardiopulmonary bypass (CPB), potentially caused by differential expression of connexin40 (Cx40) and connexin43 (Cx43). In 16 pediatric patients undergoing corrective heart surgery, connexin mRNA expression was studied in volume overloaded (VO group, n=8) and not overloaded (NO group, n=8) right atrial myocardium, excised before and after CPB. Additionally, in eight of these patients ventricular specimens were investigated. The atrial Cx43 expression decreased during CPB, which was restricted to the VO group (p=0.008). In contrast, atrial Cx40 mRNA did not change during CPB. In ventricular myocardium compared to atrial mRNA levels, Cx40 was lower (p=0.006) and Cx43 higher (p=0.017) expressed, without significant change during CPB. This study revealed a significant influence of CPB and the underlying heart defect on Cx43 expression.

Non-invasive gene-expression-based detection of well-developed collateral function in individuals with and without coronary artery disease

BACKGROUND: In patients with coronary artery disease (CAD), a well grown collateral circulation has been shown to be important. The aim of this prospective study using peripheral blood monocytes was to identify marker genes for an extensively grown coronary collateral circulation. METHODS: Collateral flow index (CFI) was obtained invasively by angioplasty pressure sensor guidewire in 160 individuals (110 patients with CAD, and 50 individuals without CAD). RNA was extracted from monocytes followed by microarray-based gene-expression analysis. 76 selected genes were analysed by real-time polymerase chain reaction (PCR). A receiver operating characteristics analysis based on differential gene expression was then performed to separate individuals with poor (CFI<0.21) and well-developed collaterals (CFI>or=0.21) Thereafter, the influence of the chemokine MCP-1 on the expression of six selected genes was tested by PCR. RESULTS: The expression of 203 genes significantly correlated with CFI (p = 0.000002-0.00267) in patients with CAD and 56 genes in individuals without CAD (p = 00079-0.0430). Biological pathway analysis revealed 76 of those genes belonging to four different pathways: angiogenesis, integrin-, platelet-derived growth factor-, and transforming growth factor beta-signalling. Three genes in each subgroup differentiated with high specificity among individuals with low and high CFI (>or=0.21). Two out of these genes showed pronounced differential expression between the two groups after cell stimulation with MCP-1. CONCLUSIONS: Genetic factors play a role in the formation and the preformation of the coronary collateral circulation. Gene expression analysis in peripheral blood monocytes can be used for non-invasive differentiation between individuals with poorly and with well grown collaterals. MCP-1 can influence the arteriogenic potential of monocytes.

Gene expression in cultured endometrium from women with different outcomes following IVF

Estradiol and progesterone are crucial for the acquisition of receptivity and the change in transcriptional activity of target genes in the implantation window. The aim of this study was to differentiate the regulation of genes in the endometrium of patients with recurrent implantation failure (IF) versus those who became pregnant after in vitro fertilization (IVF) treatment. Moreover, the effect of embryo-derived factors on endometrial transcriptional activity was studied. Nine women with known IVF outcome (IF, M, miscarriage, OP, ongoing pregnancy) and undergoing hysteroscopy with endometrial biopsy were enrolled. Biopsies were taken during the midluteal phase. After culture in the presence of embryo-conditioned IVF media, total RNA was extracted and submitted to reverse transcription, target cDNA synthesis, biotin labelling, fragmentation and hybridization using the Affymetrix Human Genome U133A 2.0 Chip. Differential expression of selected genes was re-analysed by quantitative PCR, in which the results were calculated as threshold cycle differences between the groups and normalized to Glyceraldehyde phosphate dehydrogenase and beta-actin. Differences were seen for several genes from endometrial tissue between the IF and the pregnancy groups, and when comparing OP with M, 1875 up- and 1807 down-regulated genes were returned. Real-time PCR analysis confirmed up-regulation for somatostatin, PLAP-2, mucin 4 and CD163, and down-regulation of glycodelin, IL-24, CD69, leukaemia inhibitory factor and prolactin receptor between Op and M. When the different embryo-conditioned media were compared, no significant differential regulation could be demonstrated. Although microarray profiling may currently not be sensitive enough for studying the effects of embryo-derived factors on the endometrium, the observed differences in gene expression between M and OP suggest that it will become an interesting tool for the identification of fertility-relevant markers produced by the endometrium.

Gender modulates the expression of calcium-regulating proteins in pediatric atrial myocardium

A differential expression of sarcoplasmic reticulum calcium-ATPase (SERCA2a) and phospholamban (PLB) characterizes the remodeling process in heart failure and atrial arrhythmias in adult patients. Gender is known to modulate the course and prognosis of different forms of heart disease. We hypothesized that gender plays a role in molecular changes of myocardial calcium regulating components already in childhood. Moreover, we studied the influence of volume overloaded (VO) on SERCA2a and PLB in pediatric patients. Quantitative reverse transcription-polymerase chain reaction was used to measure mRNA expression of SERCA2a and PLB in atrial myocardium from 30 pediatric patients (12 girls, 18 boys). Eighteen patients had VO right atria, and 12 patients had not-overloaded atria (NO). Protein expression was studied by Western blot. In the entire population, SERCA2a and PLB expression was not different between girls and boys. If hemodynamic overload was taken into account, SERCA2a mRNA expression was significantly reduced in the VO group compared with the NO group (P = 0.021). The VO versus NO difference was restricted to boys, which corresponds to a highly significant interaction of gender versus VO status (P = 0.002). The PLB to SERCA2a protein ratio was significantly lower in girls (P = 0.028). The decrease in SERCA2a mRNA expression in VO atrial myocardium and the PLB to SERCA2a ratio of protein expression was modulated by gender in this pediatric population. To our knowledge, this study is the first to show the impact of gender on the differential expression of calcium-regulating components in pediatric cardiac patients.

Distinct microRNA Expression Profile in Prostate Cancer Patients with Early Clinical Failure and the Impact of let-7 as Prognostic Marker in High-Risk Prostate Cancer

Background The identification of additional prognostic markers to improve risk stratification and to avoid overtreatment is one of the most urgent clinical needs in prostate cancer (PCa). MicroRNAs, being important regulators of gene expression, are promising biomarkers in various cancer entities, though the impact as prognostic predictors in PCa is poorly understood. The aim of this study was to identify specific miRNAs as potential prognostic markers in high-risk PCa and to validate their clinical impact. Methodology and Principal Findings We performed miRNA-microarray analysis in a high-risk PCa study group selected by their clinical outcome (clinical progression free survival (CPFS) vs. clinical failure (CF)). We identified seven candidate miRNAs (let-7a/b/c, miR-515-3p/5p, -181b, -146b, and -361) that showed differential expression between both groups. Further qRT-PCR analysis revealed down-regulation of members of the let-7 family in the majority of a large, well-characterized high-risk PCa cohort (n = 98). Expression of let-7a/b/and -c was correlated to clinical outcome parameters of this group. While let-7a showed no association or correlation with clinical relevant data, let-7b and let-7c were associated with CF in PCa patients and functioned partially as independent prognostic marker. Validation of the data using an independent high-risk study cohort revealed that let-7b, but not let-7c, has impact as an independent prognostic marker for BCR and CF. Furthermore, we identified HMGA1, a non-histone protein, as a new target of let-7b and found correlation of let-7b down-regulation with HMGA1 over-expression in primary PCa samples. Conclusion Our findings define a distinct miRNA expression profile in PCa cases with early CF and identified let-7b as prognostic biomarker in high-risk PCa. This study highlights the importance of let-7b as tumor suppressor miRNA in high-risk PCa and presents a basis to improve individual therapy for high-risk PCa patients.

Plakoglobin as a regulator of desmocollin gene expression.

Desmosomes are cell adhesion junctions required for the normal development and maintenance of mammalian tissues and organs such as the skin, skin appendages, and the heart. The goal of this study was to investigate how desmocollins (DSCs), transmembrane components of desmosomes, are regulated at the transcriptional level. We hypothesized that differential expression of the Dsc2 and Dsc3 genes is a prerequisite for normal development of skin appendages. We demonstrate that plakoglobin (Pg) in conjunction with lymphoid enhancer-binding factor 1 (Lef-1) differentially regulates the proximal promoters of these two genes. Specifically, we found that Lef-1 acts as a switch activating Dsc2 and repressing Dsc3 in the presence of Pg. Interestingly, we also determined that NF-κB pathway components, the downstream effectors of the ectodysplasin-A (EDA)/ ectodysplasin-A receptor (EDAR)/NF-κB signaling cascade, can activate Dsc2 expression. We hypothesize that Lef-1 and EDA/EDAR/NF-κB signaling contribute to a shift in Dsc isoform expression from Dsc3 to Dsc2 in placode keratinocytes. It is tempting to speculate that this shift is required for the invasive growth of placode keratinocytes into the dermis, a crucial step in skin appendage formation.

Independent polled mutations leading to complex gene expression differences in cattle

The molecular regulation of horn growth in ruminants is still poorly understood. To investigate this process, we collected 1019 hornless (polled) animals from different cattle breeds. High-density SNP genotyping confirmed the presence of two different polled associated haplotypes in Simmental and Holstein cattle co-localized on BTA 1. We refined the critical region of the Simmental polled mutation to 212 kb and identified an overlapping region of 932 kb containing the Holstein polled mutation. Subsequently, whole genome sequencing of polled Simmental and Holstein cows was used to determine polled associated genomic variants. By genotyping larger cohorts of animals with known horn status we found a single perfectly associated insertion/deletion variant in Simmental and other beef cattle confirming the recently published possible Celtic polled mutation. We identified a total of 182 sequence variants as candidate mutations for polledness in Holstein cattle, including an 80 kb genomic duplication and three SNPs reported before. For the first time we showed that hornless cattle with scurs are obligate heterozygous for one of the polled mutations. This is in contrast to published complex inheritance models for the bovine scurs phenotype. Studying differential expression of the annotated genes and loci within the mapped region on BTA 1 revealed a locus (LOC100848215), known in cow and buffalo only, which is higher expressed in fetal tissue of wildtype horn buds compared to tissue of polled fetuses. This implicates that the presence of this long noncoding RNA is a prerequisite for horn bud formation. In addition, both transcripts associated with polledness in goat and sheep (FOXL2 and RXFP2), show an overexpression in horn buds confirming their importance during horn development in cattle.