Early epigenetic downregulation of WNK2 kinase during pancreatic ductal adenocarcinoma development

Pancreatic ductal adenocarcinoma (PDAC) is usually incurable. Contrary to genetic mechanisms involved in PDAC pathogenesis, epigenetic alterations are ill defined. Here, we determine the contribution of epigenetically silenced genes to the development of PDAC. We analyzed enriched, highly methylated DNAs from PDACs, chronic pancreatitis (CP) and normal tissues using CpG island microarrays and identified WNK2 as a prominent candidate tumor suppressor gene being downregulated early in PDAC development. WNK2 was further investigated in tissue microarrays, methylation analysis of early pancreatic intraepithelial neoplasia (PanIN), mouse models for PDAC and pancreatitis, re-expression studies after demethylation, and cell growth assays using WNK2 overexpression. Demethylation assays confirmed the link between methylation and expression. WNK2 hypermethylation was higher in tumor than in surrounding inflamed tissues and was observed in PanIN lesions as well as in a PDAC mouse model. WNK2 mRNA and protein expressions were lower in PDAC and CP compared with normal tissues both in patients and mouse models. Overexpression of WNK2 led to reduced cell growth, and WNK2 expression in tissues correlated negatively with pERK1/2 expression, a downstream target of WNK2 responsible for cell proliferation. Downregulation of WNK2 by promoter hypermethylation occurs early in PDAC pathogenesis and may support tumor cell growth via the ERK-MAPK pathway.

Hypermethylation of TUSC5 genes in breast cancer tissue

AIM

Breast cancer (BC) is one of the most common forms of cancer amongst females. Early diagnosis, prognosis and therapy plays crucial role in the survival of patients with breast cancer. The study was aimed on identification of potential markers for early BC diagnostics by means of genome-wide comparative analysis of gene expression in cancer and normal tissue of breast.

METHODS

The analysis of gene expression in 15 invasive adenocarcinoma specimens and 15 normal breast tissue was conducted using the full-genome microarrays Sentrix HumanWD-6V3 BeadChip (Illumina). Methylation of TP53INK1 and TUSС5 promoters was interrogated by the combined bisulfite restriction analysis (COBRA).

RESULTS

Analysis of gene expression in the samples of breast adenocarcinoma revealed abnormal expression of more than 2,300 genes. While genes TFF1, S100P, ERBB2, TOP2A, CDF15, HOOK1, DNAJC12, CORO2A were up-regulated in cancer, decreased expression was found for genes TUSC5, SFRP1, PPPQR1B, NTRK4, TIMP4, BARD1, AKR1C2, TP53INK1 and others. Analysis of DNA methylation of TUSC5 by COBRA revealed higher levels of exon 1 methylation (11/12) in samples of breast cancer, whereas the gene was essentially unmethylated in matched normal appearing tissue of breast (2/12). TP53INK1 gene was methylated neither in cancer nor in normalcy.

CONCLUSION

A total of 149 genes exhibited the highest difference in expression in cancer versus normal appearing tissue of breast. Most prominent down-regulated candidates, TUSC5 and TP53INK1, were reported for the first time in breast cancer and may be considered as potential markers of the disease. Aberrant DNA hypermethylation of TUSC5 suggests epigenetic mechanism of cancer associated down-regulation.

Identification of immunohistochemical prognostic markers for survival after resection of pulmonary metastases from colorectal carcinoma

BACKGROUND

Although aggressive resection of pulmonary metastases prolongs the survival of patients with metastatic colorectal cancer, there is a need for predictive pathologic parameters to understand the key molecular events of metastatic progression. The aim of this study was to verify immunohistochemical markers in addition to established clinical parameters after surgery.

METHODS

From our subset of patients undergoing resection of pulmonary metastases from metastatic colorectal carcinoma, we analyzed 39 patients (23 men and 16 women) between 2003 and 2007. Only patients who met the criteria for a potentially curative operation were included. All patients were analyzed with regard to age and sex, primary tumor location, stage of the primary tumor, history of hepatic metastases, number of pulmonary metastases, pre-thoracotomy carcinoembryonic (CEA) serum antigen level, and the presence of thoracic lymph node metastasis. Furthermore, we immunohistochemically investigated the expression of vascular endothelial growth factor (VEGF)-D, FBJ murine osteosarcoma viral oncogene homolog B (FOS-B), and melanoma antigen (MAGE)-A in the surgical specimens of pulmonary metastatic lesions.

RESULTS

The overall 3-year survival was 50.6 %. A significantly longer survival was observed with multivariate analysis in patients with a pre-thoracotomy serum carcinoembryonic antigen level of no more than 4.2 ng/mL ( P = 0.001), and Dukes stage A or B primary tumor ( P = 0.001). A significantly longer recurrence-free survival was observed with multivariate analysis in patients without thoracic lymph node involvement compared to patients with pulmonary and/or mediastinal lymph node metastases ( P = 0.006). The stage of the primary tumor remained significant ( P = 0.029), and FOS-B expression in tumor cells showed a trend towards favorable recurrence-free survival after pulmonary metastasectomy ( P = 0.059). No statistically significant difference was found in the overall survival rate or recurrence-free survival rate of patients with expression of VEGF-D or MAGE-A antigen in pulmonary metastatic tumor cells.

CONCLUSIONS

Our results suggest that in addition to clinically prognostic factors, FOS-B expression has a debatable impact on patient survival. We conclude that the evaluation of molecular and clinical prognostic parameters at the time of pulmonary metastasectomy offers a greater understanding of the metastatic process and provides important information for patient selection.

Complete genome sequence and comparative analysis of the metabolically versatile Pseudomonas putida KT2440

Pseudomonas putida is a metabolically versatile saprophytic soil bacterium that has been certified as a biosafety host for the cloning of foreign genes. The bacterium also has considerable potential for biotechnological applications. Sequence analysis of the 6.18 Mb genome of strain KT2440 reveals diverse transport and metabolic systems. Although there is a high level of genome conservation with the pathogenic Pseudomonad Pseudomonas aeruginosa (85% of the predicted coding regions are shared), key virulence factors including exotoxin A and type III secretion systems are absent. Analysis of the genome gives insight into the non-pathogenic nature of P. putida and points to potential new applications in agriculture, biocatalysis, bioremediation and bioplastic production.

Contig selection in physical mapping

In physical mapping, one orders a set of genetic landmarks or a library of cloned fragments of DNA according to their position in the genome. Our approach to physical mapping divides the problem into smaller and easier subproblems by partitioning the probe set into independent parts (probe contigs). For this purpose we introduce a new distance function between probes, the averaged rank distance (ARD) derived from bootstrap resampling of the raw data. The ARD measures the pairwise distances of probes within a contig and smoothes the distances of probes across different contigs. It shows distinct jumps at contig borders. This makes it appropriate for contig selection by clustering. We have designed a physical mapping algorithm that makes use of these observations and seems to be particularly well suited to the delineation of reliable contigs. We evaluated our method on data sets from two physical mapping projects. On data from the recently sequenced bacterium Xylella fastidiosa, the probe contig set produced by the new method was evaluated using the probe order derived from the sequence information. Our approach yielded a basically correct contig set. On this data we also compared our method to an approach which uses the number of supporting clones to determine contigs. Our map is much more accurate. In comparison to a physical map of Pasteurella haemolytica that was computed using simulated annealing, the newly computed map is considerably cleaner. The results of our method have already proven helpful for the design of experiments aimed at further improving the quality of a map.

Identification and characterization of a SET/NAP protein encoded by a brain-specific gene, MB20

A new member of the NAP/SET gene family, named MB20, was isolated from a mouse brain cDNA library by virtue of its CAG trinucleotide repetitive sequence and a brain-specific gene expression pattern. The complementary DNA sequence predicted an open reading frame of 545 amino acids, with four copies of an 11-amino-acid direct repeat. The consensus sequence for these repeats, PKE-P--K-EE, is present in the largest subunit of murine neurofilament (NF-H). The MB20 protein sequence is homologous to nucleosome assembly proteins of several species, and its C-terminus is homologous to SET proteins. Immunoblot analysis revealed that MB20 protein is expressed in the brain. Transient transfection and immunofluorescence microscopy demonstrated that MB20 is distributed in the cytoplasm as well as in the nucleus. Deletion of the N-terminal end imparts the complete localization of MB20 protein to the nucleus. The ability of MB20 to bind histone proteins was analyzed by sucrose gradient sedimentation and by retention of histone proteins by immobilized MB20 protein. On the basis of its expression pattern, predicted sequence, and protein properties, we propose that MB20 plays a unique role in modulating nucleosome structure and gene expression during brain development.

Application of bootstrap techniques to physical mapping

Ordering genetic markers or clones from a genomic library into a physical map is a central problem in genetics. In the presence of errors, there is no efficient algorithm known that solves this problem. Based on a standard heuristic algorithm for it, we present a method to construct a confidence neighborhood for a computed solution. We compute a confidence value for putative local solutions derived from bootstrap replicates of the original solution. In the reliable parts, the confidence neighborhood and the computed solution tend to coincide. In regions that are ill-defined by the data, the neighborhood contains additional reasonable alternatives. This offers the possibility of designing further experiments for the badly defined regions to improve the quality of the physical map. We analyze our approach by a simulation study and by application to a dataset of the genome of the bacterium Xylella fastidiosa.

Progress in Arabidopsis genome sequencing and functional genomics

Arabidopsis thaliana has a relatively small genome of approximately 130 Mb containing about 10% repetitive DNA. Genome sequencing studies reveal a gene-rich genome, predicted to contain approximately 25000 genes spaced on average every 4.5 kb. Between 10 to 20% of the predicted genes occur as clusters of related genes, indicating that local sequence duplication and subsequent divergence generates a significant proportion of gene families. In addition to gene families, repetitive sequences comprise individual and small clusters of two to three retroelements and other classes of smaller repeats. The clustering of highly repetitive elements is a striking feature of the A. thaliana genome emerging from sequence and other analyses.

Sequence and analysis of chromosome 4 of the plant Arabidopsis thaliana

The higher plant Arabidopsis thaliana (Arabidopsis) is an important model for identifying plant genes and determining their function. To assist biological investigations and to define chromosome structure, a coordinated effort to sequence the Arabidopsis genome was initiated in late 1996. Here we report one of the first milestones of this project, the sequence of chromosome 4. Analysis of 17.38 megabases of unique sequence, representing about 17% of the genome, reveals 3,744 protein coding genes, 81 transfer RNAs and numerous repeat elements. Heterochromatic regions surrounding the putative centromere, which has not yet been completely sequenced, are characterized by an increased frequency of a variety of repeats, new repeats, reduced recombination, lowered gene density and lowered gene expression. Roughly 60% of the predicted protein-coding genes have been functionally characterized on the basis of their homology to known genes. Many genes encode predicted proteins that are homologous to human and Caenorhabditis elegans proteins.

Computational aspects of expression data

Several experimental techniques are available nowadays to study the spectrum of genes expressed in a cell at a specific moment. Typically, such methods generate large amounts of expression data that may be hard to interpret. Here we review computational questions and approaches resulting from the various experimental techniques.

Mutation of a gene for a Drosophila kinesin-like protein, Klp38B, leads to failure of cytokinesis

Mutations in a gene (Klp38B) encoding a novel kinesin-like protein in Drosophila melanogaster lead to the formation of polyploid cells in the larval central nervous system and in the follicle cells of adult egg chambers. Some homozygous mutants survive to adulthood and also exhibit morphological defects indicative of abnormal cell cycle progression, including rough eyes, missing bristles, and abnormal abdominal cuticles. In larval brains, there is no accumulation of mitotic cells and the frequency of anaphase figures is comparable to wild type, suggesting that nuclear division is not affected. Such brains contain polyploid cells with metaphase and anaphase chromosomes associated with bipolar spindles. Such spindles have a number of unseparated centrosomes at their poles reflecting the degree of polyploidy of the cell. Follicle cells frequently contain two nuclei of roughly equal size. Taken together, we conclude that these Klp38B mutations lead to a failure of cytokinesis resulting in polyploidy, and discuss whether or not this is a direct effect of the mutation.

Expression of a Drosophila GATA transcription factor in multiple tissues in the developing embryos. Identification of homozygous lethal mutants with P-element insertion at the promoter region

GATA transcription factors are DNA-binding proteins that recognize the core consensus sequence, WGATAR. Previous studies indicated that GATA factors play ann important role in the development of tissue-specific functions in vertebrates. Here we report the identification of a new Drosophila melanogaster GATA factor, dGATAc, which displays a distinct expression pattern in embryos. The local concentration of dGATAc transcripts varies at different stages, being most prominent in the procephalic region at stages 6-10 and in the posterior spiracles, the gut, and the central nervous system at stages 11-13. On the basis of its predicted sequence, DNA-binding assays were performed to confirm that the dGATAc gene encodes a zinc finger protein that can bind the GATA consensus motif with predicted specificity. Two independent mutants carrying a P-element insertion at the dGATAc gene promoter region were identified that are homozygous lethal at the embryonic stage. Using a genetic scheme, it was demonstrated that the lack of dGATAc function can block normal embryonic development. Our results suggest that the dGATAc protein is a tissue-specific transcription factor that is vital to the development of multiple organ systems in D. melanogaster.

Cloning and expression of an evolutionary conserved single-domain angiotensin converting enzyme from Drosophila melanogaster

Mammalian somatic angiotensin converting enzyme (EC 3.4.15.1, ACE) consists of two highly homologous (N- and C-) domains encoded by a duplicated gene. We have identified an apparent single-domain (67 kDa) insect angiotensin converting enzyme (AnCE) in embryos of Drosophila melanogaster which converts angiotensin I to angiotensin II (Km, 365 microM), removes Phe-Arg from the C terminus of bradykinin (Km, 22 microM), and is inhibited by ACE inhibitors, captopril (IC50 = 1.1 x 10(-9) M) and trandolaprilat (IC50 = 1.6 x 10(-8) M). We also report the cloning and expression of a Drosophila AnCE cDNA which codes for a single-domain 615-amino acid protein with a predicted 17-amino acid signal peptide and regions with high levels of homology to both the N- and C-domains of mammalian somatic ACE, especially around the active site consensus sequence. Northern analysis identified a single 2.1-kilobase mRNA in Drosophila embryos, and Southern analysis of Drosophila genomic DNA indicates that the insect gene is not duplicated. When expressed in COS-7 cells, the AnCE protein is a secreted enzyme, which converts angiotensin I to angiotensin II and is inhibited by captopril (IC50 = 5.6 x 10(-9) M) and trandolaprilat (IC50 = 2 x 10(-8) M). The evolutionary significance of these results is discussed.

Algorithms and software tools for ordering clone libraries: application to the mapping of the genome of Schizosaccharomyces pombe

A complete set of software tools to aid the physical mapping of a genome has been developed and successfully applied to the genomic mapping of the fission yeast Schizosaccharomyces pombe. Two approaches were used for ordering single-copy hybridisation probes: one was based on the simulated annealing algorithm to order all probes, and another on inferring the minimum-spanning subset of the probes using a heuristic filtering procedure. Both algorithms produced almost identical maps, with minor differences in the order of repetitive probes and those having identical hybridisation patterns. A separate algorithm fitted the clones to the established probe order. Approaches for handling experimental noise and repetitive elements are discussed. In addition to these programs and the database management software, tools for visualizing and editing the data are described. The issues of combining the information from different libraries are addressed. Also, ways of handling multiple-copy probes and non-hybridisation data are discussed.