Glioblastoma gene network reconstruction and ontology analysis by online bioinformatics tools.

Glioblastoma is the most aggressive type of brain tumors resistant to a number of antitumor drugs. The problem of therapy and drug treatment course is complicated by extremely high heterogeneity in the benign cell populations, the random arrangement of tumor cells, and polymorphism of their nuclei. The pathogenesis of gliomas needs to be studied using modern cellular technologies, genome- and transcriptome-wide technologies of high-throughput sequencing, analysis of gene expression on microarrays, and methods of modern bioinformatics to find new therapy targets. Functional annotation of genes related to the disease could be retrieved based on genetic databases and cross-validated by integrating complementary experimental data. Gene network reconstruction for a set of genes (proteins) proved to be effective approach to study mechanisms underlying disease progression. We used online bioinformatics tools for annotation of gene list for glioma, reconstruction of gene network and comparative analysis of gene ontology categories. The available tools and the databases for glioblastoma gene analysis are discussed together with the recent progress in this field.

Prognostic Significance of BIRC5/Survivin in Breast Cancer: Results from Three Independent Cohorts.

Breast cancer (BC) histological and molecular classifications significantly improved the treatment strategy and prognosis. Inhibitor of apoptosis BIRC5/survivin is often overexpressed in cancers, however, indications of its importance in BC are inconsistent. We integrate BIRC5 protein and mRNA measures with clinical associates and long-term outcome in three independent cohorts Protein levels of BIRC5 were measured in primary lysates of 845 patients of the West Swedish BC cohort (VGR-BC) and linked to 5- and 27-years survival. The results were externally validated in transcriptomic data from METABRIC and SCAN-B cohorts. Survival analysis showed that high levels of BIRC5 were consistently associated with a poor probability of 5-year overall survival. High BIRC5 in VGR-BC contributed negatively to the disease-specific survival at 5 and 27 years. Subsets with different status by ER (estrogen receptor) expression and presence of nodal metastasis supported independent association of high BIRC5 with poor prognosis in all cohorts. In METABRIC and SCAN-B cohorts, high levels of BIRC5 mRNA were associated with the basal-like and luminal B molecular BC subtypes and with increasing histologic grade. BIRC5 is a sensitive survival marker that acts independent of ER and nodal status, and its levels need to be considered when making treatment decisions.

Trans-Ethnic Mapping of BANK1 Identifies Two Independent SLE-Risk Linkage Groups Enriched for Co-Transcriptional Splicing Marks.

BANK1 is a susceptibility gene for several systemic autoimmune diseases in several populations. Using the genome-wide association study (GWAS) data from Europeans (EUR) and African Americans (AA), we performed an extensive fine mapping of ankyrin repeats 1 (BANK1). To increase the SNP density, we used imputation followed by univariate and conditional analysis, combined with a haplotypic and expression quantitative trait locus (eQTL) analysis. The data from Europeans showed that the associated region was restricted to a minimal and dependent set of SNPs covering introns two and three, and exon two. In AA, the signal found in the Europeans was split into two independent effects. All of the major risk associated SNPs were eQTLs, and the risks were associated with an increased BANK1 gene expression. Functional annotation analysis revealed the enrichment of repressive B cell epigenomic marks (EZH2 and H3K27me3) and a strong enrichment of splice junctions. Furthermore, one eQTL located in intron two, rs13106926, was found within the binding site for RUNX3, a transcriptional activator. These results connect the local genome topography, chromatin structure, and the regulatory landscape of BANK1 with co-transcriptional splicing of exon two. Our data defines a minimal set of risk associated eQTLs predicted to be involved in the expression of BANK1 modulated through epigenetic regulation and splicing. These findings allow us to suggest that the increased expression of BANK1 will have an impact on B-cell mediated disease pathways.

Altered expression of multiple genes involved in retinoic acid biosynthesis in human colorectal cancer.

All-trans-retinoic acid (atRA), the oxidized form of vitamin A (retinol), regulates a wide variety of biological processes, such as cell proliferation and differentiation. Multiple alcohol, retinol and retinaldehyde dehydrogenases (ADHs, RDHs, RALDHs) as well as aldo-keto reductases (AKRs) catalyze atRA production. The reduced atRA biosynthesis has been observed in several human tumors, including colorectal cancer. However, subsets of atRA-synthesizing enzymes have not been determined in colorectal tumors. We investigated the expression patterns of genes involved in atRA biosynthesis in normal human colorectal tissues, primary carcinomas and cancer cell lines by RT-PCR. These genes were identified using transcriptomic data analysis (expressed sequence tags, RNA-sequencing, microarrays). Our results indicate that each step of the atRA biosynthesis pathway is dysregulated in colorectal cancer. Frequent and significant decreases in the mRNA levels of the ADH1B, ADH1C, RDHL, RDH5 and AKR1B10 genes were observed in a majority of colorectal carcinomas. The expression levels of the RALDH1 gene were reduced, and the expression levels of the cytochrome CYP26A1 gene increased. The human colon cancer cell lines showed a similar pattern of changes in the mRNA levels of these genes. A dramatic reduction in the expression of genes encoding the predominant retinol-oxidizing enzymes could impair atRA production. The most abundant of these genes, ADH1B and ADH1C, display decreased expression during progression from adenoma to early and more advanced stage of colorectal carcinomas. The diminished atRA biosynthesis may lead to alteration of cell growth and differentiation in the colon and rectum, thus contributing to the progression of colorectal cancer.

Impact of recombination on polymorphism of genes encoding Kunitz-type protease inhibitors in the genus Solanum.

BACKGROUND: The group of Kunitz-type protease inhibitors (KPI) from potato is encoded by a polymorphic family of multiple allelic and non-allelic genes. The previous explanations of the KPI variability were based on the hypothesis of random mutagenesis as a key factor of KPI polymorphism. RESULTS: KPI-A genes from the genomes of Solanum tuberosum cv. Istrinskii and the wild species Solanum palustre were amplified by PCR with subsequent cloning in plasmids. True KPI sequences were derived from comparison of the cloned copies. "Hot spots" of recombination in KPI genes were independently identified by DnaSP 4.0 and TOPALi v2.5 software. The KPI-A sequence from potato cv. Istrinskii was found to be 100% identical to the gene from Solanum nigrum. This fact illustrates a high degree of similarity of KPI genes in the genus Solanum. Pairwise comparison of KPI A and B genes unambiguously showed a non-uniform extent of polymorphism at different nt positions. Moreover, the occurrence of substitutions was not random along the strand. Taken together, these facts contradict the traditional hypothesis of random mutagenesis as a principal source of KPI gene polymorphism. The experimentally found mosaic structure of KPI genes in both plants studied is consistent with the hypothesis suggesting recombination of ancestral genes. The same mechanism was proposed earlier for other resistance-conferring genes in the nightshade family (Solanaceae). Based on the data obtained, we searched for potential motifs of site-specific binding with plant DNA recombinases. During this work, we analyzed the sequencing data reported by the Potato Genome Sequencing Consortium (PGSC), 2011 and found considerable inconsistence of their data concerning the number, location, and orientation of KPI genes of groups A and B. CONCLUSIONS: The key role of recombination rather than random point mutagenesis in KPI polymorphism was demonstrated for the first time.

Activation of the hTERT expression in squamous cell cervical carcinoma is not associated with gene amplification.

The hTERT gene encodes the telomerase catalytic subunit that plays a key role in cancer cell immortalization. Earlier, hTERT amplification was detected in squamous cell cervical carcinomas (SCC), however possible relations between elevated hTERT mRNA level and gene amplification was not studied. Here, we compared the hTERT expression and copy number in the same tumors by quantitative real-time PCR. The hTERT DNA copy number was virtually unchanged in all 33 studied tumors, when compared to normal tissues. This result was confirmed using two reference genes beta-actin and beta-D-glucuronidase. Nevertheless, the activation of hTERT expression was found in 80% of cases (37/46, p<0.001). There was no correlation between the degree of mRNA increase and the tumor size and/or presence of metastases. No hTERT gene expression was observed in 20% of cases (9/46), while the control GADPH expression was unchanged. The detected elevation of the hTERT mRNA level was found using primers specific to functionally active full-length isoform of mRNA. Similar results were obtained with SCC cell lines carrying human papilloma virus (HPV) genomes. We conclude that frequent activation of hTERT expression in SCC is not associated with gene amplification.

Invariant amino acids essential for decoding function of polypeptide release factor eRF1.

In eukaryotic ribosome, the N domain of polypeptide release factor eRF1 is involved in decoding stop signals in mRNAs. However, structure of the decoding site remains obscure. Here, we specifically altered the stop codon recognition pattern of human eRF1 by point mutagenesis of the invariant Glu55 and Tyr125 residues in the N domain. The 3D structure of generated eRF1 mutants was not destabilized as demonstrated by calorimetric measurements and calculated free energy perturbations. In mutants, the UAG response was most profoundly and selectively affected. Surprisingly, Glu55Arg mutant completely retained its release activity. Substitution of the aromatic ring in position 125 reduced response toward all stop codons. This result demonstrates the critical importance of Tyr125 for maintenance of the intact structure of the eRF1 decoding site. The results also suggest that Tyr125 is implicated in recognition of the 3d stop codon position and probably forms an H-bond with Glu55. The data point to a pivotal role played by the YxCxxxF motif (positions 125-131) in purine discrimination of the stop codons. We speculate that eRF1 decoding site is formed by a 3D network of amino acids side chains.

Common and specific amino acid residues in the prokaryotic polypeptide release factors RF1 and RF2: possible functional implications.

Termination of protein synthesis is promoted in ribosomes by proper stop codon discrimination by class 1 polypeptide release factors (RFs). A large set of prokaryotic RFs differing in stop codon specificity, RF1 for UAG and UAA, and RF2 for UGA and UAA, was analyzed by means of a recently developed computational method allowing identification of the specificity-determining positions (SDPs) in families composed of proteins with similar but not identical function. Fifteen SDPs were identified within the RF1/2 superdomain II/IV known to be implicated in stop codon decoding. Three of these SDPs had particularly high scores. Five residues invariant for RF1 and RF2 [invariant amino acid residues (IRs)] were spatially clustered with the highest-scoring SDPs that in turn were located in two zones within the SDP/IR area. Zone 1 (domain II) included PxT and SPF motifs identified earlier by others as 'discriminator tripeptides'. We suggest that IRs in this zone take part in the recognition of U, the first base of all stop codons. Zone 2 (domain IV) possessed two SDPs with the highest scores not identified earlier. Presumably, they also take part in stop codon binding and discrimination. Elucidation of potential functional role(s) of the newly identified SDP/IR zones requires further experiments.