Comprehensive genome methylation and whole genome expression analysis in penile carcinoma: Uncovering new molecular markers

Background: Penile carcinoma (PeCa) is frequently associated with high morbidity rates. Unlikely of the vast majority of tumors, there is no molecular markers described that are able to assist in diagnosis and prognosis or with potential to be therapeutic targets in PeCa. Patients and methods: DNA methylation status (244K Human DNA Methylation Microarray platform, Agilent Technologies) and large-scale expression analysis (4x44K Whole Human Genome Microarray, Agilent Technologies) were performed in 35 and 37 PeCa, respectively. Quantitative bisulfite pyrosequencing (qBP) and RT-qPCR were used to validate the findings in 93 samples. HPV status was assessed using the Linear Array HPV Genotyping kit (Roche Molecular Diagnostics, CA, USA). Results: Methylome analysis revealed 171 hypermethylated and 449 hypomethylated CpGs sites and the transcriptome profiling showed 2986 down- and 2817 over-expressed genes. HPV positivity was found in 32.7% of the cases, mainly the HPV16. The integrative analysis in 32 PeCa revealed a panel of 96 genes with inverse correlation between methylation and gene expression levels. The CpG hypermetlylation and gene downexpression, was confirmed for TWIST1, RSOP2, SOX3, SOX17, CD133, OTX2, HOXA3 and MEIS. In addition, BIRC5, DNMT1 and DNMT3B presented low levels of methylation and overexpression. The comparison of the results with clinical findings revealed that LIN28A, NKX2.2, NKX2.3, LHX5, BDNF, FOXA1 and CDX2 were associated with poor prognosis features. Conclusion: Putative prognostic markers were detected revealing that DNA methylation modulates the expression of several genes in PeCa. These data may prove instrumental for biomarker discovery in clinics and molecular epidemiology of PeCa.

Gene expression analysis of laryngeal squamous cell carcinoma

Laryngeal squamous cell carcinoma (LSCC) is one of the most common malignancies of the head and neck tumors Zhang et al., 2013 [1]). Previous studies have associated its occurrence with social activities, such as tobacco and alcohol consumption (Hashibe et al., 2007a [2]; Hashibe et al., 2007b [3]; Shangina et al., 2006 [4]). Here, we performed a genome-wide gene expression profiling in thirty-one patients positively diagnosed for LSCC, in order to investigate new targets involved in tumorigenesis.

Xylella fastidiosa: An in vivo system to study possible survival strategies within citrus xylem vessels based on global gene expression analysis

Xylella fastidiosa inhabits the plant xylem, a nutrient-poor environment, so that mechanisms to sense and respond to adverse environmental conditions are extremely important for bacterial survival in the plant host. Although the complete genome sequences of different Xylella strains have been determined, little is known about stress responses and gene regulation in these organisms. In this work, a DNA microarray was constructed containing 2,600 ORFs identified in the genome sequencing project of Xylella fastidiosa 9a5c strain, and used to check global gene expression differences in the bacteria when it is infecting a symptomatic and a tolerant citrus tree. Different patterns of expression were found in each variety, suggesting that bacteria are responding differentially according to each plant xylem environment. The global gene expression profile was determined and several genes related to bacterial survival in stressed conditions were found to be differentially expressed between varieties, suggesting the involvement of different strategies for adaptation to the environment. The expression pattern of some genes related to the heat shock response, toxin and detoxification processes, adaptation to atypical conditions, repair systems as well as some regulatory genes are discussed in this paper. DNA microarray proved to be a powerful technique for global transcriptome analyses. This is one of the first studies of Xylella fastidiosa gene expression in vivo which helped to increase insight into stress responses and possible bacterial survival mechanisms in the nutrient-poor environment of xylem vessels.

Resistance to octreotide LAR in acromegalic patients with high SSTR2 expression: analysis of AIP expression

We present here the clinical and molecular data of two patients with acromegaly treated with octreotide LAR after non-curative surgery, and who presented different responses to therapy. Somatostatin receptor type 2 and 5 (SSTR2 and SSTR5), and aryl hydrocarbon receptor-interacting protein (AIP) expression levels were analyzed by qPCR. In both cases, high SSTR2 and low SSTR5 expression levels were detected; however, only one of the patients achieved disease control after octreotide LAR therapy. When we analyzed AIP expression levels of both cases, the patient whose disease was controlled after therapy exhibited AIP expression levels that were two times higher than the patient whose disease was still active. These two cases illustrate that, although the currently available somatostatin analogs bind preferentially to SSTR2, some patients are not responsive to therapy despite high expression of this receptor. This difference could be explained by differences in post-receptor signaling pathways, including the recently described involvement of AIP. Arq Bras Endocrinol Metab. 2012;56(8):501-6

rpb2 is a reliable reference gene for quantitative gene expression analysis in the dermatophyte Trichophyton rubrum

The selection of reference genes used for data normalization to quantify gene expression by real-time PCR amplifications (qRT-PCR) is crucial for the accuracy of this technique. In spite of this, little information regarding such genes for qRT-PCR is available for gene expression analyses in pathogenic fungi. Thus, we investigated the suitability of eight candidate reference genes in isolates of the human dermatophyte Trichophyton rubrum subjected to several environmental challenges, such as drug exposure, interaction with human nail and skin, and heat stress. The stability of these genes was determined by geNorm, NormFinder and Best-Keeper programs. The gene with the most stable expression in the majority of the conditions tested was rpb2 (DNA-dependent RNA polymerase II), which was validated in three T. rubrum strains. Moreover, the combination of rpb2 and chs1 (chitin synthase) genes provided for the most reliable qRT-PCR data normalization in T. rubrum under a broad range of biological conditions. To the best of our knowledge this is the first report on the selection of reference genes for qRT-PCR data normalization in dermatophytes and the results of these studies should permit further analysis of gene expression under several experimental conditions, with improved accuracy and reliability.

CLONING AND EXPRESSION ANALYSIS OF ALLOGRAFT INFLAMMATORY FACTOR TYPE 1 IN COELOMOCYTES OF ANTARCTIC SEA URCHIN (STERECHINUS NEUMAYERI)

We have cloned and characterized for the first time an allograft inflammatory factor 1 (Sn-AIF-1) from the Antarctic sea urchin. We report the cloning of Sn-AIF-1 cDNA and the characterization of its expression in coelomocytes after a bacterial challenge. The cDNA Sn-AIF-1 has a size of 608 bp and encodes a polypeptide of 151 aa. The deduced amino acid sequence has a putative size of 17.430 Da, an isoelectric point of 4.92, and shows 2 elongation factor handlike motifs that normally bind calcium ions. BLAST analysis revealed close matches with other known AIF-1. The deduced amino acid sequence of Sn-AIF-1 showed high homology with AIF-1 in vertebrates such as fish, mice, and humans; and in the case of invertebrates, the major degree of identity (55%) was with a predicted sequence of the purple sea urchin AIF-1, and 52% corresponded to a sponge. Expression of Sn-AIF-1 mRNA was analyzed by qPCR. Sn-AIF-1 mRNA expression was measured from coelomocytes after a bacterial challenge using RT-PCR and revealed that the gene was upregulated after 24 h. Sn-AIF-1 could participate in the inflammatory response, particularly in the activation of coelomocytes and their survival.

Expression analysis of stem cell-related genes reveal OCT4 as a predictor of poor clinical outcome in medulloblastoma

Aberrant expression of stem cell-related genes in tumors may confer more primitive and aggressive traits affecting clinical outcome. Here, we investigated expression and prognostic value of the neural stem cell marker CD133, as well as of the pluripotency genes LIN28 and OCT4 in 37 samples of pediatric medulloblastoma, the most common and challenging type of embryonal tumor. While most medulloblastoma samples expressed CD133 and LIN28, OCT4 expression was found to be more sporadic, with detectable levels occurring in 48% of tumors. Expression levels of OCT4, but not CD133 or LIN28, were significantly correlated with shorter survival (P <= 0.0001). Median survival time of patients with tumors hyperexpressing OCT4 and tumors displaying low/undetectable OCT4 expression were 6 and 153 months, respectively. More importantly, when patients were clinically stratified according to their risk of tumor recurrence, positive OCT4 expression in primary tumor specimens could discriminate patients classified as average risk but which further deceased within 5 years of diagnosis (median survival time of 28 months), a poor clinical outcome typical of high risk patients. Our findings reveal a previously unknown prognostic value for OCT4 expression status in medulloblastoma, which might be used as a further indicator of poor survival and aid postoperative treatment selection, with a particular potential benefit for clinically average risk patients.

Single-cell expression analysis of BMP15 and GDF9 in mature oocytes and BMPR2 in cumulus cells of women with polycystic ovary syndrome undergoing controlled ovarian hyperstimulation

To detect expression of bone morphogenetic protein 15 (BMP15) and growth differentiation factor 9 (GDF9) in oocytes, and their receptor type 2 receptor for BMPs (BMPR2) in cumulus cells in women with polycystic ovary syndrome (PCOS) undergoing in vitro fertilization (IVF), and determine if BMPR2, BMP15, and GDF9 expression correlate with hyperandrogenism in FF of PCOS patients. Prospective case-control study. Eighteen MII-oocytes and their respective cumulus cells were obtained from 18 patients with PCOS, and 48 MII-oocytes and cumulus cells (CCs) from 35 controls, both subjected to controlled ovarian hyperstimulation (COH), and follicular fluid (FF) was collected from small (10-14 mm) and large (> 18 mm) follicles. RNeasy Micro Kit (Qiagen(A (R))) was used for RNA extraction and gene expression was quantified in each oocyte individually and in microdissected cumulus cells from cumulus-oocyte complexes retrieved from preovulatory follicles using qRT-PCR. Chemiluminescence and RIA assays were used for hormone assays. BMP15 and GDF9 expression per oocyte was higher among women with PCOS than the control group. A positive correlation was found between BMPR2 transcripts and hyperandrogenism in FF of PCOS patients. Progesterone values in FF were lower in the PCOS group. We inferred that BMP15 and GDF9 transcript levels increase in mature PCOS oocytes after COH, and might inhibit the progesterone secretion by follicular cells in PCOS follicles, preventing premature luteinization in cumulus cells. BMPR2 expression in PCOS cumulus cells might be regulated by androgens.

Expression Analysis of Wound Healing Genes in Human Periapical Granulomas of Progressive and Stable Nature

Introduction: Wound healing process involves the activation of extracellular matrix components, remodeling enzymes, cellular adhesion molecules, growth factors, cytokines and chemokines genes. However, the molecular patterns underlying the healing process periapical environment remain unclear. Here we hypothesized that endodontic infection might result in an imbalance in the expression of wound healing genes involved in the pathogenesis of periapical lesions. Furthermore, we suggest that differential expression of wound healing markers in active and latent granulomas could account for different clinical outcomes for such lesions. Methods: Study samples consisted of 93 periapical granulomas collected after endodontic surgeries and 24 healthy periodontal ligament tissues collected from premolars extracted for orthodontic purposes as control samples. Of these, 10 periapical granulomas and 5 healthy periapical tissues were used for expression analysis of 84 wound healing genes by using a pathway-specific real-time polymerase chain reaction array. The remaining 83 granulomas and all 24 control specimens were used to validate the obtained array data by real-time polymerase chain reaction. Observed variations in expression of wound healing genes were analyzed according to the classification of periapical granulomas as active/progressive versus inactive/stable (as determined by receptor activator for nuclear factor kappa B ligand/osteoprotegerin expression ratio). Results: We observed a marked increase of 5-fold or greater in SERPINE1, TIMP1, COL1A1, COL5A1, VTN, CTGF, FGF7, TGFB1, TNF, CXCL11, ITGA4, and ITGA5 genes in the periapical granulomas when compared with control samples. SERPINE1, TIMP1, COL1A1, TGFB1, and ITGA4 mRNA expression was significantly higher in inactive compared with active periapical granulomas (P < .001), whereas TNF and CXCL11 mRNA expression was higher in active lesions (P < .001). Conclusions: The identification of novel gene targets that curb the progression status of periapical lesions might contribute to a more accurate diagnosis and lead to treatment modalities more conducive to endodontic success. (J Endod 2012;38:185-190)

Identification and expression analysis of microRNAs and targets in the biofuel crop sugarcane

Abstract Background MicroRNAs (miRNAs) are small regulatory RNAs, some of which are conserved in diverse plant genomes. Therefore, computational identification and further experimental validation of miRNAs from non-model organisms is both feasible and instrumental for addressing miRNA-based gene regulation and evolution. Sugarcane (Saccharum spp.) is an important biofuel crop with publicly available expressed sequence tag and genomic survey sequence databases, but little is known about miRNAs and their targets in this highly polyploid species. Results In this study, we have computationally identified 19 distinct sugarcane miRNA precursors, of which several are highly similar with their sorghum homologs at both nucleotide and secondary structure levels. The accumulation pattern of mature miRNAs varies in organs/tissues from the commercial sugarcane hybrid as well as in its corresponding founder species S. officinarum and S. spontaneum. Using sugarcane MIR827 as a query, we found a novel MIR827 precursor in the sorghum genome. Based on our computational tool, a total of 46 potential targets were identified for the 19 sugarcane miRNAs. Several targets for highly conserved miRNAs are transcription factors that play important roles in plant development. Conversely, target genes of lineage-specific miRNAs seem to play roles in diverse physiological processes, such as SsCBP1. SsCBP1 was experimentally confirmed to be a target for the monocot-specific miR528. Our findings support the notion that the regulation of SsCBP1 by miR528 is shared at least within graminaceous monocots, and this miRNA-based post-transcriptional regulation evolved exclusively within the monocots lineage after the divergence from eudicots. Conclusions Using publicly available nucleotide databases, 19 sugarcane miRNA precursors and one new sorghum miRNA precursor were identified and classified into 14 families. Comparative analyses between sugarcane and sorghum suggest that these two species retain homologous miRNAs and targets in their genomes. Such conservation may help to clarify specific aspects of miRNA regulation and evolution in the polyploid sugarcane. Finally, our dataset provides a framework for future studies on sugarcane RNAi-dependent regulatory mechanisms.

Multivariate gene expression analysis reveals functional connectivity changes between normal/tumoral prostates

Abstract Background Prostate cancer is a leading cause of death in the male population, therefore, a comprehensive study about the genes and the molecular networks involved in the tumoral prostate process becomes necessary. In order to understand the biological process behind potential biomarkers, we have analyzed a set of 57 cDNA microarrays containing ~25,000 genes. Results Principal Component Analysis (PCA) combined with the Maximum-entropy Linear Discriminant Analysis (MLDA) were applied in order to identify genes with the most discriminative information between normal and tumoral prostatic tissues. Data analysis was carried out using three different approaches, namely: (i) differences in gene expression levels between normal and tumoral conditions from an univariate point of view; (ii) in a multivariate fashion using MLDA; and (iii) with a dependence network approach. Our results show that malignant transformation in the prostatic tissue is more related to functional connectivity changes in their dependence networks than to differential gene expression. The MYLK, KLK2, KLK3, HAN11, LTF, CSRP1 and TGM4 genes presented significant changes in their functional connectivity between normal and tumoral conditions and were also classified as the top seven most informative genes for the prostate cancer genesis process by our discriminant analysis. Moreover, among the identified genes we found classically known biomarkers and genes which are closely related to tumoral prostate, such as KLK3 and KLK2 and several other potential ones. Conclusion We have demonstrated that changes in functional connectivity may be implicit in the biological process which renders some genes more informative to discriminate between normal and tumoral conditions. Using the proposed method, namely, MLDA, in order to analyze the multivariate characteristic of genes, it was possible to capture the changes in dependence networks which are related to cell transformation.

Semantic integration of gene expression analysis tools and data sources using software connectors

Abstract Background The study and analysis of gene expression measurements is the primary focus of functional genomics. Once expression data is available, biologists are faced with the task of extracting (new) knowledge associated to the underlying biological phenomenon. Most often, in order to perform this task, biologists execute a number of analysis activities on the available gene expression dataset rather than a single analysis activity. The integration of heteregeneous tools and data sources to create an integrated analysis environment represents a challenging and error-prone task. Semantic integration enables the assignment of unambiguous meanings to data shared among different applications in an integrated environment, allowing the exchange of data in a semantically consistent and meaningful way. This work aims at developing an ontology-based methodology for the semantic integration of gene expression analysis tools and data sources. The proposed methodology relies on software connectors to support not only the access to heterogeneous data sources but also the definition of transformation rules on exchanged data. Results We have studied the different challenges involved in the integration of computer systems and the role software connectors play in this task. We have also studied a number of gene expression technologies, analysis tools and related ontologies in order to devise basic integration scenarios and propose a reference ontology for the gene expression domain. Then, we have defined a number of activities and associated guidelines to prescribe how the development of connectors should be carried out. Finally, we have applied the proposed methodology in the construction of three different integration scenarios involving the use of different tools for the analysis of different types of gene expression data. Conclusions The proposed methodology facilitates the development of connectors capable of semantically integrating different gene expression analysis tools and data sources. The methodology can be used in the development of connectors supporting both simple and nontrivial processing requirements, thus assuring accurate data exchange and information interpretation from exchanged data.

Expression analysis and in silico characterization of intronic long noncoding RNAs in renal cell carcinoma: emerging functional associations

Abstract Background Intronic and intergenic long noncoding RNAs (lncRNAs) are emerging gene expression regulators. The molecular pathogenesis of renal cell carcinoma (RCC) is still poorly understood, and in particular, limited studies are available for intronic lncRNAs expressed in RCC Methods Microarray experiments were performed with custom-designed arrays enriched with probes for lncRNAs mapping to intronic genomic regions. Samples from 18 primary RCC tumors and 11 nontumor adjacent matched tissues were analyzed. Meta-analyses were performed with microarray expression data from three additional human tissues (normal liver, prostate tumor and kidney nontumor samples), and with large-scale public data for epigenetic regulatory marks and for evolutionarily conserved sequences. Results A signature of 29 intronic lncRNAs differentially expressed between RCC and nontumor samples was obtained (false discovery rate (FDR) <5%). A signature of 26 intronic lncRNAs significantly correlated with the RCC five-year patient survival outcome was identified (FDR <5%, p-value ≤0.01). We identified 4303 intronic antisense lncRNAs expressed in RCC, of which 22% were significantly (p <0.05) cis correlated with the expression of the mRNA in the same locus across RCC and three other human tissues. Gene Ontology (GO) analysis of those loci pointed to 'regulation of biological processes’ as the main enriched category. A module map analysis of the protein-coding genes significantly (p <0.05) trans correlated with the 20% most abundant lncRNAs, identified 51 enriched GO terms (p <0.05). We determined that 60% of the expressed lncRNAs are evolutionarily conserved. At the genomic loci containing the intronic RCC-expressed lncRNAs, a strong association (p <0.001) was found between their transcription start sites and genomic marks such as CpG islands, RNA Pol II binding and histones methylation and acetylation. Conclusion Intronic antisense lncRNAs are widely expressed in RCC tumors. Some of them are significantly altered in RCC in comparison with nontumor samples. The majority of these lncRNAs is evolutionarily conserved and possibly modulated by epigenetic modifications. Our data suggest that these RCC lncRNAs may contribute to the complex network of regulatory RNAs playing a role in renal cell malignant transformation.

MicroRNAs expression analysis in high grade gliomas

Several biomarkers had been proposed as useful parameters to better define the prognosis or to delineate new target therapy strategies for glioblastoma (GBM) patients. MicroRNAs could represent interesting molecules, for their role in tumorigenesis and cancer progression and for their specific tissue expression. Although many studies have tried to identify a specific microRNAs signature for glioblastoma, by now an exhaustive GBM microRNAs profile is far to be well defined. In this work we set up a real-time qPCR, based on LNA primers, to investigate the expression of 19 microRNAs in brain tumors, focusing our attention on GBMs. MiRNAs expression in 30 GBM paired FFPE-Fresh/Frozen samples was firstly analyzed. The good correlation obtained comparing miRNAs results confirmed the feasibility of performing miRNAs analysis starting from FFPE tissues. This leads to many advantages, as a good disposal of archival tumor and normal brain specimens and the possibility to verify the percentage of tumor cells in the analyzed sample. In the second part we compared 3 non-neoplastic brain references to use as control in miRNAs analysis. Normal adjacent the tumor, epileptic specimens and a commercial total RNA were analyzed for miRNAs expression and results showed that different non-neoplastic controls could lead to important discrepancies in GBM miRNAs profiles. Analyzing 50 FFPE GBMs using all 3 non-neoplastic references, we defined a putative GBM miRNAs signature: mir-10b, miR-21 and miR-27a resulted upregulated, while miR-7, miR-9, miR-26a, miR-31, miR-101, miR-137, miR-222 and miR-330 were downregulated. Comparing miRNAs expression among GBM group and gliomas of grade I, II and III, we obtained 3 miRNAs (miR-10b, mir-34a and miR-101) showing a different regulation status between high grade and low grade gliomas. Intriguingly, miR-10b was upregulated in high grade and significantly downregulated in low grade gliomas, suggesting that could be a candidate for a GBM target therapy.

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.