EMMA 2--a MAGE-compliant system for the collaborative analysis and integration of microarray data

Background: Understanding transcriptional regulation by genome-wide microarray studies can contribute to unravel complex relationships between genes. Attempts to standardize the annotation of microarray data include the Minimum Information About a Microarray Experiment (MIAME) recommendations, the MAGE-ML format for data interchange, and the use of controlled vocabularies or ontologies. The existing software systems for microarray data analysis implement the mentioned standards only partially and are often hard to use and extend. Integration of genomic annotation data and other sources of external knowledge using open standards is therefore a key requirement for future integrated analysis systems. Results: The EMMA 2 software has been designed to resolve shortcomings with respect to full MAGE-ML and ontology support and makes use of modern data integration techniques. We present a software system that features comprehensive data analysis functions for spotted arrays, and for the most common synthesized oligo arrays such as Agilent, Affymetrix and NimbleGen. The system is based on the full MAGE object model. Analysis functionality is based on R and Bioconductor packages and can make use of a compute cluster for distributed services. Conclusion: Our model-driven approach for automatically implementing a full MAGE object model provides high flexibility and compatibility. Data integration via SOAP-based web-services is advantageous in a distributed client-server environment as the collaborative analysis of microarray data is gaining more and more relevance in international research consortia. The adequacy of the EMMA 2 software design and implementation has been proven by its application in many distributed functional genomics projects. Its scalability makes the current architecture suited for extensions towards future transcriptomics methods based on high-throughput sequencing approaches which have much higher computational requirements than microarrays.

ArrayMining: a modular web-application for microarray analysis combining ensemble and consensus methods with cross-study normalization

Background: Statistical analysis of DNA microarray data provides a valuable diagnostic tool for the investigation of genetic components of diseases. To take advantage of the multitude of available data sets and analysis methods, it is desirable to combine both different algorithms and data from different studies. Applying ensemble learning, consensus clustering and cross-study normalization methods for this purpose in an almost fully automated process and linking different analysis modules together under a single interface would simplify many microarray analysis tasks. Results: We present ArrayMining.net, a web-application for microarray analysis that provides easy access to a wide choice of feature selection, clustering, prediction, gene set analysis and cross-study normalization methods. In contrast to other microarray-related web-tools, multiple algorithms and data sets for an analysis task can be combined using ensemble feature selection, ensemble prediction, consensus clustering and cross-platform data integration. By interlinking different analysis tools in a modular fashion, new exploratory routes become available, e.g. ensemble sample classification using features obtained from a gene set analysis and data from multiple studies. The analysis is further simplified by automatic parameter selection mechanisms and linkage to web tools and databases for functional annotation and literature mining. Conclusion: ArrayMining.net is a free web-application for microarray analysis combining a broad choice of algorithms based on ensemble and consensus methods, using automatic parameter selection and integration with annotation databases.

mRNA-Seq and microarray development for the Grooved carpet shell clam, Ruditapes decussatus: a functional approach to unravel host -parasite interaction

Background The Grooved Carpet shell clam Ruditapes decussatus is the autochthonous European clam and the most appreciated from a gastronomic and economic point of view. The production is in decline due to several factors such as Perkinsiosis and habitat invasion and competition by the introduced exotic species, the manila clam Ruditapes philippinarum. After we sequenced R. decussatus transcriptome we have designed an oligo microarray capable of contributing to provide some clues on molecular response of the clam to Perkinsiosis. Results A database consisting of 41,119 unique transcripts was constructed, of which 12,479 (30.3%) were annotated by similarity. An oligo-DNA microarray platform was then designed and applied to profile gene expression in R. decussatus heavily infected by Perkinsus olseni. Functional annotation of differentially expressed genes between those two conditionswas performed by gene set enrichment analysis. As expected, microarrays unveil genes related with stress/infectious agents such as hydrolases, proteases and others. The extensive role of innate immune system was also analyzed and effect of parasitosis upon expression of important molecules such as lectins reviewed. Conclusions This study represents a first attempt to characterize Ruditapes decussatus transcriptome, an important marine resource for the European aquaculture. The trancriptome sequencing and consequent annotation will increase the available tools and resources for this specie, introducing the possibility of high throughput experiments such as microarrays analysis. In this specific case microarray approach was used to unveil some important aspects of host-parasite interaction between the Carpet shell clam and Perkinsus, two non-model species, highlighting some genes associated with this interaction. Ample information was obtained to identify biological processes significantly enriched among differentially expressed genes in Perkinsus infected versus non-infected gills. An overview on the genes related with the immune system on R. decussatus transcriptome is also reported.

Utility of a panviral microarray for detection of swine respiratory viruses in clinical samples

Several factors have recently converged, elevating the need for highly parallel diagnostic platforms that have the ability to detect many known, novel, and emerging pathogenic agents simultaneously. Panviral DNA microarrays represent the most robust approach for massively parallel viral surveillance and detection. The Virochip is a panviral DNA microarray that is capable of detecting all known viruses, as well as novel viruses related to known viral families, in a single assay and has been used to successfully identify known and novel viral agents in clinical human specimens. However, the usefulness and the sensitivity of the Virochip platform have not been tested on a set of clinical veterinary specimens with the high degree of genetic variance that is frequently observed with swine virus field isolates. In this report, we investigate the utility and sensitivity of the Virochip to positively detect swine viruses in both cell culture-derived samples and clinical swine samples. The Virochip successfully detected porcine reproductive and respiratory syndrome virus (PRRSV) in serum containing 6.10 × 10(2) viral copies per microliter and influenza A virus in lung lavage fluid containing 2.08 × 10(6) viral copies per microliter. The Virochip also successfully detected porcine circovirus type 2 (PCV2) in serum containing 2.50 × 10(8) viral copies per microliter and porcine respiratory coronavirus (PRCV) in turbinate tissue homogenate. Collectively, the data in this report demonstrate that the Virochip can successfully detect pathogenic viruses frequently found in swine in a variety of solid and liquid specimens, such as turbinate tissue homogenate and lung lavage fluid, as well as antemortem samples, such as serum.

Analysis of the genome content of Lactococcus garvieae by genomic interspecies microarray hybridization

BACKGROUND Lactococcus garvieae is a bacterial pathogen that affects different animal species in addition to humans. Despite the widespread distribution and emerging clinical significance of L. garvieae in both veterinary and human medicine, there is almost a complete lack of knowledge about the genetic content of this microorganism. In the present study, the genomic content of L. garvieae CECT 4531 was analysed using bioinformatics tools and microarray-based comparative genomic hybridization (CGH) experiments. Lactococcus lactis subsp. lactis IL1403 and Streptococcus pneumoniae TIGR4 were used as reference microorganisms. RESULTS The combination and integration of in silico analyses and in vitro CGH experiments, performed in comparison with the reference microorganisms, allowed establishment of an inter-species hybridization framework with a detection threshold based on a sequence similarity of >or= 70%. With this threshold value, 267 genes were identified as having an analogue in L. garvieae, most of which (n = 258) have been documented for the first time in this pathogen. Most of the genes are related to ribosomal, sugar metabolism or energy conversion systems. Some of the identified genes, such as als and mycA, could be involved in the pathogenesis of L. garvieae infections. CONCLUSIONS In this study, we identified 267 genes that were potentially present in L. garvieae CECT 4531. Some of the identified genes could be involved in the pathogenesis of L. garvieae infections. These results provide the first insight into the genome content of L. garvieae.

Discovery of stable and variable differences in the Mycobacterium avium subsp. paratuberculosis type I, II, and III genomes by pan-genome microarray analysis

Mycobacterium avium subsp. paratuberculosis is an important animal pathogen widely disseminated in the environment that has also been associated with Crohn's disease in humans. Three M. avium subsp. paratuberculosis genomotypes are recognized, but genomic differences have not been fully described. To further investigate these potential differences, a 60-mer oligonucleotide microarray (designated the MAPAC array), based on the combined genomes of M. avium subsp. paratuberculosis (strain K-10) and Mycobacterium avium subsp. hominissuis (strain 104), was designed and validated. By use of a test panel of defined M. avium subsp. paratuberculosis strains, the MAPAC array was able to identify a set of large sequence polymorphisms (LSPs) diagnostic for each of the three major M. avium subsp. paratuberculosis types. M. avium subsp. paratuberculosis type II strains contained a smaller genomic complement than M. avium subsp. paratuberculosis type I and M. avium subsp. paratuberculosis type III genomotypes, which included a set of genomic regions also found in M. avium subsp. hominissuis 104. Specific PCRs for genes within LSPs that differentiated M. avium subsp. paratuberculosis types were devised and shown to accurately screen a panel (n = 78) of M. avium subsp. paratuberculosis strains. Analysis of insertion/deletion region INDEL12 showed deletion events causing a reduction in the complement of mycobacterial cell entry genes in M. avium subsp. paratuberculosis type II strains and significantly altering the coding of a major immunologic protein (MPT64) associated with persistence and granuloma formation. Analysis of MAPAC data also identified signal variations in several genomic regions, termed variable genomic islands (vGIs), suggestive of transient duplication/deletion events. vGIs contained significantly low GC% and were immediately flanked by insertion sequences, integrases, or short inverted repeat sequences. Quantitative PCR demonstrated that variation in vGI signals could be associated with colony growth rate and morphology.

A Putative Otu Domain-containing Protein 1 Deubiquitinating Enzyme Is Differentially Expressed In Thyroid Cancer And Identifies Less-aggressive Tumours.

This study aimed to identify novel biomarkers for thyroid carcinoma diagnosis and prognosis. We have constructed a human single-chain variable fragment (scFv) antibody library that was selected against tumour thyroid cells using the BRASIL method (biopanning and rapid analysis of selective interactive ligands) and phage display technology. One highly reactive clone, scFv-C1, with specific binding to papillary thyroid tumour proteins was confirmed by ELISA, which was further tested against a tissue microarray that comprised of 229 thyroid tissues, including: 110 carcinomas (38 papillary thyroid carcinomas (PTCs), 42 follicular carcinomas, 30 follicular variants of PTC), 18 normal thyroid tissues, 49 nodular goitres (NG) and 52 follicular adenomas. The scFv-C1 was able to distinguish carcinomas from benign lesions (P=0.0001) and reacted preferentially against T1 and T2 tumour stages (P=0.0108). We have further identified an OTU domain-containing protein 1, DUBA-7 deubiquitinating enzyme as the scFv-binding antigen using two-dimensional polyacrylamide gel electrophoresis and mass spectrometry. The strategy of screening and identifying a cell-surface-binding antibody against thyroid tissues was highly effective and resulted in a useful biomarker that recognises malignancy among thyroid nodules and may help identify lower-risk cases that can benefit from less-aggressive management.

Ki-1/57 And Cgi-55 Ectopic Expression Impact Cellular Pathways Involved In Proliferation And Stress Response Regulation.

Ki-1/57 (HABP4) and CGI-55 (SERBP1) are regulatory proteins and paralogs with 40.7% amino acid sequence identity and 67.4% similarity. Functionally, they have been implicated in the regulation of gene expression on both the transcriptional and mRNA metabolism levels. A link with tumorigenesis is suggested, since both paralogs show altered expression levels in tumor cells and the Ki-1/57 gene is found in a region of chromosome 9q that represents a haplotype for familiar colon cancer. However, the target genes regulated by Ki-1/57 and CGI-55 are unknown. Here, we analyzed the alterations of the global transcriptome profile after Ki-1/57 or CGI-55 overexpression in HEK293T cells by DNA microchip technology. We were able to identify 363 or 190 down-regulated and 50 or 27 up-regulated genes for Ki-1/57 and CGI-55, respectively, of which 20 were shared between both proteins. Expression levels of selected genes were confirmed by qRT-PCR both after protein overexpression and siRNA knockdown. The majority of the genes with altered expression were associated to proliferation, apoptosis and cell cycle control processes, prompting us to further explore these contexts experimentally. We observed that overexpression of Ki-1/57 or CGI-55 results in reduced cell proliferation, mainly due to a G1 phase arrest, whereas siRNA knockdown of CGI-55 caused an increase in proliferation. In the case of Ki-1/57 overexpression, we found protection from apoptosis after treatment with the ER-stress inducer thapsigargin. Together, our data give important new insights that may help to explain these proteins putative involvement in tumorigenic events.

Analysis Of The Contribution Of Immunologically-detectable Her2, Steroid Receptors And Of The Triple-negative" Tumor Status To Disease-free And Overall Survival Of Women With Epithelial Ovarian Cancer."

We assessed associations between steroid receptors including: estrogen-alpha, estrogen-beta, androgen receptor, progesterone receptor, the HER2 status and triple-negative epithelial ovarian cancer (ERα-/PR-/HER2-; TNEOC) status and survival in women with epithelial ovarian cancer. The study included 152 women with primary epithelial ovarian cancer. The status of steroid receptor and HER2 was determined by immunohistochemistry. Disease-free and overall survival were calculated and compared with steroid receptor and HER2 status as well as clinicopathological features using the Cox Proportional Hazards model. A mean follow-up period of 43.6 months (interquartile range=41.4 months) was achieved where 44% of patients had serous tumor, followed by mucinous (23%), endometrioid (9%), mixed (9%), undifferentiated (8.5%) and clear cell tumors (5.3%). ER-alpha staining was associated with grade II-III tumors. Progesterone receptor staining was positively associated with a Body Mass Index≥25. Androgen receptor positivity was higher in serous tumors. In stand-alone analysis of receptor contribution to survival, estrogen-alpha positivity was associated with greater disease-free survival. However, there was no significant association between steroid receptor expression, HER2 status, or TNEOC status, and overall survival. Although estrogen-alpha, androgen receptor, progesterone receptor and the HER2 status were associated with key clinical features of the women and pathological characteristics of the tumors, these associations were not implicated in survival. Interestingly, women with TNEOC seem to fare the same way as their counterparts with non-TNEOC.

Leaf-, Panel- And Latex-expressed Sequenced Tags From The Rubber Tree (hevea Brasiliensis) Under Cold-stressed And Suboptimal Growing Conditions: The Development Of Gene-targeted Functional Markers For Stress Response.

Hevea brasiliensis is a native species of the Amazon Basin of South America and the primary source of natural rubber worldwide. Due to the occurrence of South American Leaf Blight disease in this area, rubber plantations have been extended to suboptimal regions. Rubber tree breeding is time-consuming and expensive, but molecular markers can serve as a tool for early evaluation, thus reducing time and costs. In this work, we constructed six different cDNA libraries with the aim of developing gene-targeted molecular markers for the rubber tree. A total of 8,263 reads were assembled, generating 5,025 unigenes that were analyzed; 912 expressed sequence tags (ESTs) represented new transcripts, and two sequences were highly up-regulated by cold stress. These unigenes were scanned for microsatellite (SSR) regions and single nucleotide polymorphisms (SNPs). In total, 169 novel EST-SSR markers were developed; 138 loci were polymorphic in the rubber tree, and 98 % presented transferability to six other Hevea species. Locus duplication was observed in H. brasiliensis and other species. Additionally, 43 SNP markers in 13 sequences that showed similarity to proteins involved in stress response, latex biosynthesis and developmental processes were characterized. cDNA libraries are a rich source of SSR and SNP markers and enable the identification of new transcripts. The new markers developed here will be a valuable resource for linkage mapping, QTL identification and other studies in the rubber tree and can also be used to evaluate the genetic variability of other Hevea species, which are valuable assets in rubber tree breeding.

[the New Classification Of Breast Cancers: Finding The Luminal A].

To compare the distributions of patients with clinical-pathological subtypes of luminal B-like breast cancer according to the 2011 and 2013 St. Gallen International Breast Cancer Conference Expert Panel. We studied 142 women with breast cancer who were positive to estrogen receptor and had been treated in São Paulo state, southeast Brazil. The expression of the following receptors was assessed by immunohistochemistry: estrogen, progesterone (PR) and Ki-67. The expression of HER-2 was measured by fluorescent in situ hybridization analysis in tissue microarray. There were 29 cases of luminal A breast cancers according to the 2011 St. Gallen International Breast Cancer Conference Expert Panel that were classified as luminal B-like in the 2013 version. Among the 65 luminal B-like breast cancer cases, 29 (45%) were previous luminal A tumors, 15 cases (20%) had a Ki-67 >14% and were at least 20% PR positive and 21 cases (35%) had Ki-67 >14% and more than 20% were PR positive. The 2013 St. Gallen consensus updated the definition of intrinsic molecular subtypes and increased the number of patients classified as having luminal B-like breast cancer in our series, for whom the use of cytotoxic drugs will probably be proposed with additional treatment cost.

Down-regulation Of Slc8a1 As A Putative Apoptosis-evasion Mechanism By Modulation Of Calcium Levels In Penile Carcinoma.

The SLC8A1 gene, which encodes the Na(+)/Ca(2+) exchanger, plays a key role in calcium homeostasis. Our previous gene expression oligoarray data revealed SLC8A1 underexpression in penile carcinoma (PeCa). The aim of this study was to investigate whether the dysregulation of SLC8A1 expression is associated with apoptosis and cell proliferation in PeCa, via modulation of calcium concentration. The underlying mechanisms of SLC8A1 underexpression were also explored, focusing on copy number alteration and microRNA. Transcript levels of SLC8A1 gene and miR-223 were evaluated by quantitative PCR, comparing PeCa samples with normal glans tissues. SLC8A1 copy number was evaluated by microarray-based comparative genomic hybridization (array-CGH). Caspase-3 and Ki-67 immunostaining, as well as calcium distribution by Laser Ablation Imaging Inductively Coupled Plasma Mass Spectrometry [LA(i)-ICP-MS], were investigated in both normal and tumor samples. Confirming our previous data, SLC8A1 underexpression was detected in PeCa samples (P=0.001) and was not associated with gene copy number loss. In contrast, overexpression of miR-223 (P=0.002) was inversely correlated with SLC8A1 (P=0.015, r=-0.426), its putative repressor. In addition, SLC8A1 underexpression was associated with decreased calcium distribution, high Ki-67 and low caspase-3 immunoexpression in PeCa when compared with normal tissues. Down-regulation of the SLC8A1 gene, most likely mediated by its regulator miR-223, can lead to reduced calcium levels in PeCa and, consequently, to suppression of apoptosis and increased tumor cell proliferation. These data suggest that the miR-223-NCX1-calcium-signaling axis may represent a potential therapeutic approach in PeCa.

Tshr Intronic Polymorphisms (rs179247 And Rs12885526) And Their Role In The Susceptibility Of The Brazilian Population To Graves' Disease And Graves' Ophthalmopathy.

Intronic thyroid-stimulating hormone receptor polymorphisms have been associated with the risk for both Graves' disease and Graves' ophthalmopathy, but results have been inconsistent among different populations. We aimed to investigate the influence of thyroid-stimulating hormone receptor intronic polymorphisms in a large well-characterized population of GD patients. We studied 279 Graves' disease patients (231 females and 48 males, 39.80 ± 11.69 years old), including 144 with Graves' ophthalmopathy, matched to 296 healthy control individuals. Thyroid-stimulating hormone receptor genotypes of rs179247 and rs12885526 were determined by Real Time PCR TaqMan(®) SNP Genotyping. A multivariate analysis showed that the inheritance of the thyroid-stimulating hormone receptor AA genotype for rs179247 increased the risk for Graves' disease (OR = 2.821; 95 % CI 1.595-4.990; p = 0.0004), whereas the thyroid-stimulating hormone receptor GG genotype for rs12885526 increased the risk for Graves' ophthalmopathy (OR = 2.940; 95 % CI 1.320-6.548; p = 0.0083). Individuals with Graves' ophthalmopathy also presented lower mean thyrotropin receptor antibodies levels (96.3 ± 143.9 U/L) than individuals without Graves' ophthalmopathy (98.3 ± 201.9 U/L). We did not find any association between the investigated polymorphisms and patients clinical features or outcome. We demonstrate that thyroid-stimulating hormone receptor intronic polymorphisms are associated with the susceptibility to Graves' disease and Graves' ophthalmopathy in the Brazilian population, but do not appear to influence the disease course.

Polymorphisms Of Cell Cycle Control Genes Influence The Development Of Sporadic Medullary Thyroid Carcinoma.

The role of key cell cycle regulation genes such as, CDKN1B, CDKN2A, CDKN2B, and CDKN2C in sporadic medullary thyroid carcinoma (s-MTC) is still largely unknown. In order to evaluate the influence of inherited polymorphisms of these genes on the pathogenesis of s-MTC, we used TaqMan SNP genotyping to examine 45 s-MTC patients carefully matched with 98 controls. A multivariate logistic regression analysis demonstrated that CDKN1B and CDKN2A genes were related to s-MTC susceptibility. The rs2066827*GT+GG CDKN1B genotype was more frequent in s-MTC patients (62.22%) than in controls (40.21%), increasing the susceptibility to s-MTC (OR=2.47; 95% CI=1.048-5.833; P=0.038). By contrast, the rs11515*CG+GG of CDKN2A gene was more frequent in the controls (32.65%) than in patients (15.56%), reducing the risk for s-MTC (OR=0.174; 95% CI=0.048-0.627; P=0.0075). A stepwise regression analysis indicated that two genotypes together could explain 11% of the total s-MTC risk. In addition, a relationship was found between disease progression and the presence of alterations in the CDKN1A (rs1801270), CDKN2C (rs12885), and CDKN2B (rs1063192) genes. WT rs1801270 CDKN1A patients presented extrathyroidal tumor extension more frequently (92%) than polymorphic CDKN1A rs1801270 patients (50%; P=0.0376). Patients with the WT CDKN2C gene (rs12885) presented larger tumors (2.9±1.8 cm) than polymorphic patients (1.5±0.7 cm; P=0.0324). On the other hand, patients with the polymorphic CDKN2B gene (rs1063192) presented distant metastases (36.3%; P=0.0261). In summary, we demonstrated that CDKN1B and CDKN2A genes are associated with susceptibility, whereas the inherited genetic profile of CDKN1A, CDKN2B, and CDKN2C is associated with aggressive features of tumors. This study suggests that profiling cell cycle genes may help define the risk and characterize s-MTC aggressiveness.