Transcriptome analysis of a parasitic clade V nematode:comparative analysis of potential molecular anthelmintic targets in Cylicostephanus goldi

Clade V nematodes comprise several parasitic species that include the cyathostomins, primary helminth pathogens of horses. Next generation transcriptome datasets are available for eight parasitic clade V nematodes, although no equine parasites are included in this group. Here, we report next generation transcriptome sequencing analysis for the common cyathostomin species, Cylicostephanus goldi. A cDNA library was generated from RNA extracted from 17 C. goldi male and female adult parasites. Following sequencing using a 454 GS FLX pyrosequencer, a total of 475,215 sequencing reads were generated, which were assembled into 26,910 contigs. Using Gene Ontology and Kyoto Encyclopedia of Genes and Genomes databases, 27% of the transcriptome was annotated. Further in-depth analysis was carried out by comparing the C. goldi dataset with the next generation transcriptomes and genomes of other clade V nematodes, with the Oesophagostomum dentatum transcriptome and the Haemonchus contortus genome showing the highest levels of sequence identity with the cyathostomin dataset (45%). The C. goldi transcriptome was mined for genes associated with anthelmintic mode of action and/or resistance. Sequences encoding proteins previously associated with the three major anthelmintic classes used in horses were identified, with the exception of the P-glycoprotein group. Targeted resequencing of the glutamate gated chloride channel α4 subunit (glc-3), one of the primary targets of the macrocyclic lactone anthelmintics, was performed for several cyathostomin species. We believe this study reports the first transcriptome dataset for an equine helminth parasite, providing the opportunity for in-depth analysis of these important parasites at the molecular level. Sequences encoding enzymes involved in key processes and genes associated with levamisole/pyrantel and macrocyclic lactone resistance, in particular the glutamate gated chloride channels, were identified. This novel data will inform cyathostomin biology and anthelmintic resistance studies in future.

Molecular pedigree reconstruction and estimation of evolutionary parameters in a wild Atlantic salmon river system with incomplete sampling: a power analysis

Background: Pedigree reconstruction using genetic analysis provides a useful means to estimate fundamental population biology parameters relating to population demography, trait heritability and individual fitness when combined with other sources of data. However, there remain limitations to pedigree reconstruction in wild populations, particularly in systems where parent-offspring relationships cannot be directly observed, there is incomplete sampling of individuals, or molecular parentage inference relies on low quality DNA from archived material. While much can still be inferred from incomplete or sparse pedigrees, it is crucial to evaluate the quality and power of available genetic information a priori to testing specific biological hypotheses. Here, we used microsatellite markers to reconstruct a multi-generation pedigree of wild Atlantic salmon (Salmo salar L.) using archived scale samples collected with a total trapping system within a river over a 10 year period. Using a simulation-based approach, we determined the optimal microsatellite marker number for accurate parentage assignment, and evaluated the power of the resulting partial pedigree to investigate important evolutionary and quantitative genetic characteristics of salmon in the system.

Results: We show that at least 20 microsatellites (ave. 12 alleles/locus) are required to maximise parentage assignment and to improve the power to estimate reproductive success and heritability in this study system. We also show that 1.5 fold differences can be detected between groups simulated to have differing reproductive success, and that it is possible to detect moderate heritability values for continuous traits (h(2) similar to 0.40) with more than 80% power when using 28 moderately to highly polymorphic markers.

Conclusion: The methodologies and work flow described provide a robust approach for evaluating archived samples for pedigree-based research, even where only a proportion of the total population is sampled. The results demonstrate the feasibility of pedigree-based studies to address challenging ecological and evolutionary questions in free-living populations, where genealogies can be traced only using molecular tools, and that significant increases in pedigree assignment power can be achieved by using higher numbers of markers.

Genome analysis of a simultaneously predatory and prey-independent, novel Bdellovibrio bacteriovorus from the River Tiber, supports in silico predictions of both ancient and recent lateral gene transfer from diverse bacteria

BACKGROUND: Evolution equipped Bdellovibrio bacteriovorus predatory bacteria to invade other bacteria, digesting and replicating, sealed within them thus preventing nutrient-sharing with organisms in the surrounding environment. Bdellovibrio were previously described as "obligate predators" because only by mutations, often in gene bd0108, are 1 in ~1x10(7) of predatory lab strains of Bdellovibrio converted to prey-independent growth. A previous genomic analysis of B. bacteriovorus strain HD100 suggested that predatory consumption of prey DNA by lytic enzymes made Bdellovibrio less likely than other bacteria to acquire DNA by lateral gene transfer (LGT). However the Doolittle and Pan groups predicted, in silico, both ancient and recent lateral gene transfer into the B. bacteriovorus HD100 genome.

RESULTS: To test these predictions, we isolated a predatory bacterium from the River Tiber- a good potential source of LGT as it is rich in diverse bacteria and organic pollutants- by enrichment culturing with E. coli prey cells. The isolate was identified as B. bacteriovorus and named as strain Tiberius. Unusually, this Tiberius strain showed simultaneous prey-independent growth on organic nutrients and predatory growth on live prey. Despite the prey-independent growth, the homolog of bd0108 did not have typical prey-independent-type mutations. The dual growth mode may reflect the high carbon content of the river, and gives B. bacteriovorus Tiberius extended non-predatory contact with the other bacteria present. The HD100 and Tiberius genomes were extensively syntenic despite their different cultured-terrestrial/freshly-isolated aquatic histories; but there were significant differences in gene content indicative of genomic flux and LGT. Gene content comparisons support previously published in silico predictions for LGT in strain HD100 with substantial conservation of genes predicted to have ancient LGT origins but little conservation of AT-rich genes predicted to be recently acquired.

CONCLUSIONS: The natural niche and dual predatory, and prey-independent growth of the B. bacteriovorus Tiberius strain afforded it extensive non-predatory contact with other marine and freshwater bacteria from which LGT is evident in its genome. Thus despite their arsenal of DNA-lytic enzymes; Bdellovibrio are not always predatory in natural niches and their genomes are shaped by acquiring whole genes from other bacteria.

IGFBP7 promoter methylation and gene expression analysis in prostate cancer

PURPOSE: IGFBP7 belongs to a family of insulin-like growth factor-1 regulatory binding proteins. IGFBP7 hypermethylation is associated with its down-regulation in various carcinomas. In prostate cancer IGFBP7 down-regulation has been widely reported but to our knowledge the mechanisms behind this event are unknown. We performed a denaturing high performance liquid chromatography screening and validation strategy to profile the methylation status of IGFBP7 in prostate cancer.

MATERIALS AND METHODS: We combined denaturing high performance liquid chromatography and bisulfite sequencing to examine IGFBP7 methylation in a panel of prostate cancer cell lines. Quantitative methylation specific polymerase chain reaction was used to determine methylation levels in prostate tissue specimens of primary prostate cancer, histologically benign prostate adjacent to tumor, high grade prostatic intraepithelial neoplasia and benign prostatic hyperplasia. IGFBP7 gene expression was measured by quantitative methylation specific polymerase chain reaction in cell lines and tissue specimens.

RESULTS: IGFBP7 was methylated in the 4 prostate cancer cell lines DU145, LNCaP, PC-3 and 22RV1. Quantitative methylation specific polymerase chain reaction analysis revealed that promoter methylation was associated with decreased IGFBP7 expression. Quantitative methylation specific polymerase chain reaction showed that IGFBP7 methylation was more frequently detected in prostate cancer (60% (31/52)) and high grade prostatic intraepithelial neoplasia (40% (6/15)) samples compared to histologically benign prostate adjacent to tumor (10%) and benign prostatic hyperplasia (0%) samples.

CONCLUSIONS: To our knowledge this is the first report of aberrant IGFBP7 promoter hypermethylation and concurrent IGFBP7 gene silencing in prostate cancer cell lines. Results demonstrate that CpG methylation of IGFBP7 may represent a novel biomarker of prostate cancer and pre-invasive neoplasms. Thus, future examination of IGFBP7 methylation and expression in a larger patient cohort, including bodily fluids, is justified to further evaluate its role in a diagnostic and prognostic setting.

Recognition of O6MeG lesions by MGMT and mismatch repair proficiency may be a prerequisite for low-dose radiation hypersensitivity

Low-dose hyper-radiosensitivity (HRS) is the phenomenon whereby cells exposed to radiation doses of less than approximately 0.5 Gy exhibit increased cell killing relative to that predicted from back-extrapolating high-dose survival data using a linear-quadratic model. While the exact mechanism remains to be elucidated, the involvement of several molecular repair pathways has been documented. These processes in turn are also associated with the response of cells to O6-methylguanine (O6MeG) lesions. We propose a model in which the level of low-dose cell killing is determined by the efficiency of both pre-replicative repair by the DNA repair enzyme O6-methylguanine methyltransferase (MGMT) and post-replicative repair by the DNA mismatch repair (MMR) system. We therefore hypothesized that the response of cells to low doses of radiation is dependent on the expression status of MGMT and MMR proteins. MMR (MSH2, MSH6, MLH1, PMS1, PMS2) and MGMT protein expression signatures were determined in a panel of normal (PWR1E, RWPE1) and malignant (22RV1, DU145, PC3) prostate cell lines and correlated with clonogenic survival and cell cycle analysis. PC3 and RWPE1 cells (HRS positive) were associated with MGMT and MMR proficiency, whereas HRS negative cell lines lacked expression of at least one (MGMT or MMR) protein. MGMT inactivation had no significant effect on cell survival. These results indicate a possible role for MMR-dependent processing of damage produced by low doses of radiation.

Variation in DNA repair genes XRCC3, XRCC4, XRCC5 and susceptibility to myeloma

Cytogenetic analysis in myeloma reveals marked chromosomal instability. Both widespread genomic alterations and evidence of aberrant class switch recombination, the physiological process that regulates maturation of the antibody response, implicate the DNA repair pathway in disease pathogenesis. We therefore assessed 27 SNPs in three genes (XRCC3, XRCC4 and XRCC5) central to DNA repair in patients with myeloma and controls from the EpiLymph study and from an Irish hospital registry (n = 306 cases, 263 controls). For the haplotype-tagging SNP (htSNP) rs963248 in XRCC4, Allele A was significantly more frequent in cases than in controls (86.4 versus 80.8%; odds ratio 1.51; 95% confidence interval 1.10-2.08; P = 0.0133), as was the AA genotype (74 versus 65%) (P = 0.026). Haplotype analysis was performed using Unphased for rs963248 in combination with additional SNPs in XRCC4. The strongest evidence of association came from the A-T haplotype from rs963248-rs2891980 (P = 0.008). For XRCC5, the genotype GG from rs1051685 was detected in 10 cases from different national populations but in only one control (P = 0.015). This SNP is located in the 3'-UTR of XRCC5. Overall, these data provide support for the hypothesis that common variation in the genes encoding DNA repair proteins contributes to susceptibility to myeloma.

The emerging roles of DNA methylation in the clinical management of prostate cancer

Aberrant DNA methylation is one of the hallmarks of carcinogenesis and has been recognized in cancer cells for more than 20 years. The role of DNA methylation in malignant transformation of the prostate has been intensely studied, from its contribution to the early stages of tumour development to the advanced stages of androgen independence. The most significant advances have involved the discovery of numerous targets such as GSTP1, Ras-association domain family 1A (RASSF1A) and retinoic acid receptor beta2 (RARbeta2) that become inactivated through promoter hypermethylation during the course of disease initiation and progression. This has provided the basis for translational research into methylation biomarkers for early detection and prognosis of prostate cancer. Investigations into the causes of these methylation events have yielded little definitive data. Aberrant hypomethylation and how it impacts upon prostate cancer has been less well studied. Herein we discuss the major developments in the fields of prostate cancer and DNA methylation, and how this epigenetic modification can be harnessed to address some of the key issues impeding the successful clinical management of prostate cancer.

Expression of the NF-kappaB targets BCL2 and BIRC5/Survivin characterizes small B-cell and aggressive B-cell lymphomas, respectively

Nuclear factor kappa B (NF-kappaB) activation has been proposed as a cardinal feature of tumourigenesis, although the precise mechanism, frequency, relevance, and extent of NF-kappaB activation in lymphomas remain to be fully elucidated. In this study, expression profiling and tissue microarray studies of 209 and 323 non-Hodgkin's lymphomas (NHLs) respectively, including the most frequent sub-types of NHL, were employed to generate a hypothesis concerning the most common NF-kappaB targets in NHL. These analyses showed that NF-kappaB activation is a common phenomenon in NHL, resulting in the expression of distinct sets of NF-kappaB target genes, depending on the cell context. BCL2 and BIRC5/Survivin were identified as key NF-kappaB targets and their expression distinguished small and aggressive B-cell lymphomas, respectively. Interestingly, in the vast majority of B-cell lymphomas, the expression of these markers was mutually exclusive. A set of genes was identified whose expression correlates either with BIRC5/Survivin or with BCL2. BIRC5/Survivin expression, in contrast to BCL2, was associated with a signature of cell proliferation (overexpression of cell cycle control, DNA repair, and polymerase genes), which may contribute to the aggressive phenotype and poor prognosis of these lymphomas. Strikingly, mantle cell lymphoma and chronic lymphocytic leukaemia expressed highly elevated levels of BCL2 protein and mRNA, higher than that observed in reactive mantle zone cells or even in follicular lymphomas, where BCL2 expression is deregulated through the t(14;18) translocation. In parallel with this observation, BIRC5/Survivin expression was higher in Burkitt's lymphoma and diffuse large B-cell lymphoma than in non-tumoural germinal centre cells. In vitro studies confirmed that NF-kappaB activation contributes to the expression of both markers. In cell lines representing aggressive lymphomas, NF-kappaB inhibition resulted in a decrease in BIRC5/Survivin expression. Meanwhile, in chronic lymphocytic leukaemia (CLL)-derived lymphocytes, NF-kappaB inhibition resulted in a marked decrease in BCL2 expression.

Automated tumor analysis for molecular profiling in lung cancer

The discovery and clinical application of molecular biomarkers in solid tumors, increasingly relies on nucleic acid extraction from FFPE tissue sections and subsequent molecular profiling. This in turn requires the pathological review of haematoxylin & eosin (H&E) stained slides, to ensure sample quality, tumor DNA sufficiency by visually estimating the percentage tumor nuclei and tumor annotation for manual macrodissection. In this study on NSCLC, we demonstrate considerable variation in tumor nuclei percentage between pathologists, potentially undermining the precision of NSCLC molecular evaluation and emphasising the need for quantitative tumor evaluation. We subsequently describe the development and validation of a system called TissueMark for automated tumor annotation and percentage tumor nuclei measurement in NSCLC using computerized image analysis. Evaluation of 245 NSCLC slides showed precise automated tumor annotation of cases using Tissuemark, strong concordance with manually drawn boundaries and identical EGFR mutational status, following manual macrodissection from the image analysis generated tumor boundaries. Automated analysis of cell counts for % tumor measurements by Tissuemark showed reduced variability and significant correlation (p < 0.001) with benchmark tumor cell counts. This study demonstrates a robust image analysis technology that can facilitate the automated quantitative analysis of tissue samples for molecular profiling in discovery and diagnostics.

Mapping protein-DNA interactions using ChIP-sequencing

Chromatin immunoprecipitation (ChIP) allows enrichment of genomic regions which are associated with specific transcription factors, histone modifications, and indeed any other epitopes which are present on chromatin. The original ChIP methods used site-specific PCR and Southern blotting to confirm which regions of the genome were enriched, on a candidate basis. The combination of ChIP with genomic tiling arrays (ChIP-chip) allowed a more unbiased approach to map ChIP-enriched sites. However, limitations of microarray probe design and probe number have a detrimental impact on the coverage, resolution, sensitivity, and cost of whole-genome tiling microarray sets for higher eukaryotes with large genomes. The combination of ChIP with high-throughput sequencing technology has allowed more comprehensive surveys of genome occupancy, greater resolution, and lower cost for whole genome coverage. Herein, we provide a comparison of high-throughput sequencing platforms and a survey of ChIP-seq analysis tools, discuss experimental design, and describe a detailed ChIP-seq method.Chromatin immunoprecipitation (ChIP) allows enrichment of genomic regions which are associated with specific transcription factors, histone modifications, and indeed any other epitopes which are present on chromatin. The original ChIP methods used site-specific PCR and Southern blotting to confirm which regions of the genome were enriched, on a candidate basis. The combination of ChIP with genomic tiling arrays (ChIP-chip) allowed a more unbiased approach to map ChIP-enriched sites. However, limitations of microarray probe design and probe number have a detrimental impact on the coverage, resolution, sensitivity, and cost of whole-genome tiling microarray sets for higher eukaryotes with large genomes. The combination of ChIP with high-throughput sequencing technology has allowed more comprehensive surveys of genome occupancy, greater resolution, and lower cost for whole genome coverage. Herein, we provide a comparison of high-throughput sequencing platforms and a survey of ChIP-seq analysis tools, discuss experimental design, and describe a detailed ChIP-seq method.

AURKA overexpression accompanies dysregulation of DNA-damage response genes in invasive urothelial cell carcinoma

Invasive urothelial cell carcinoma (UCC) is characterized by increased chromosomal instability and follows an aggressive clinical course in contrast to non-invasive disease. To identify molecular processes that confer and maintain an aggressive malignant phenotype, we used a high-throughput genome-wide approach to interrogate a cohort of high and low clinical risk UCC tumors. Differential expression analyses highlighted cohesive dysregulation of critical genes involved in the G(2)/M checkpoint in aggressive UCC. Hierarchical clustering based on DNA Damage Response (DDR) genes separated tumors according to a pre-defined clinical risk phenotype. Using array-comparative genomic hybridization, we confirmed that the DDR was disrupted in tumors displaying high genomic instability. We identified DNA copy number gains at 20q13.2-q13.3 (AURKA locus) and determined that overexpression of AURKA accompanied dysregulation of DDR genes in high risk tumors. We postulated that DDR-deficient UCC tumors are advantaged by a selective pressure for AURKA associated override of M phase barriers and confirmed this in an independent tissue microarray series. This mechanism that enables cancer cells to maintain an aggressive phenotype forms a rationale for targeting AURKA as a therapeutic strategy in advanced stage UCC.

DNA microarray profiling of genes differentially regulated by the histone deacetylase inhibitors vorinostat and LBH589 in colon cancer cell lines

BACKGROUND: Despite the significant progress made in colon cancer chemotherapy, advanced disease remains largely incurable and novel efficacious chemotherapies are urgently needed. Histone deacetylase inhibitors (HDACi) represent a novel class of agents which have demonstrated promising preclinical activity and are undergoing clinical evaluation in colon cancer. The goal of this study was to identify genes in colon cancer cells that are differentially regulated by two clinically advanced hydroxamic acid HDACi, vorinostat and LBH589 to provide rationale for novel drug combination partners and identify a core set of HDACi-regulated genes.

METHODS: HCT116 and HT29 colon cancer cells were treated with LBH589 or vorinostat and growth inhibition, acetylation status and apoptosis were analyzed in response to treatment using MTS, Western blotting and flow cytometric analyses. In addition, gene expression was analyzed using the Illumina Human-6 V2 BeadChip array and Ingenuity Pathway Analysis.

RESULTS: Treatment with either vorinostat or LBH589 rapidly induced histone acetylation, cell cycle arrest and inhibited the growth of both HCT116 and HT29 cells. Bioinformatic analysis of the microarray profiling revealed significant similarity in the genes altered in expression following treatment with the two HDACi tested within each cell line. However, analysis of genes that were altered in expression in the HCT116 and HT29 cells revealed cell-line-specific responses to HDACi treatment. In addition a core cassette of 11 genes modulated by both vorinostat and LBH589 were identified in both colon cancer cell lines analyzed.

CONCLUSION: This study identified HDACi-induced alterations in critical genes involved in nucleotide metabolism, angiogenesis, mitosis and cell survival which may represent potential intervention points for novel therapeutic combinations in colon cancer. This information will assist in the identification of novel pathways and targets that are modulated by HDACi, providing much-needed information on HDACi mechanism of action and providing rationale for novel drug combination partners. We identified a core signature of 11 genes which were modulated by both vorinostat and LBH589 in a similar manner in both cell lines. These core genes will assist in the development and validation of a common gene set which may represent a molecular signature of HDAC inhibition in colon cancer.

The Cyclin D1 (CCND1) A870G polymorphism predicts clinical outcome to lapatinib and capecitabine in HER2-positive metastatic breast cancer

BACKGROUND: Lapatinib plus capecitabine emerged as an efficacious therapy in metastatic breast cancer (mBC). We aimed to identify germline single-nucleotide polymorphisms (SNPs) in genes involved in capecitabine catabolism and human epidermal receptor signaling that were associated with clinical outcome to assist in selecting patients likely to benefit from this combination.

PATIENTS AND METHODS: DNA was extracted from 240 of 399 patients enrolled in EGF100151 clinical trial (NCT00078572; clinicaltrials.gov) and SNPs were successfully evaluated in 234 patients. The associations between SNPs and clinical outcome were analyzed using Fisher's exact test, Kaplan-Meier curves, log-rank tests, likelihood ratio test within logistic or Cox regression model, as appropriate.

RESULTS: There were significant interactions between CCND1 A870G and clinical outcome. Patients carrying the A-allele were more likely to benefit from lapatinib plus capecitabine versus capecitabine when compared with patients harboring G/G (P = 0.022, 0.024 and 0.04, respectively). In patients with the A-allele, the response rate (RR) was significantly higher with lapatinib plus capecitabine (35%) compared with capecitabine (11%; P = 0.001) but not between treatments in patients with G/G (RR = 24% and 32%, respectively; P = 0.85). Time to tumor progression (TTP) was longer in patients with the A-allele treated with lapatinib plus capecitabine compared with capecitabine (median TTP = 7.9 and 3.4 months; P < 0.001), but not in patients with G/G (median TTP = 6.1 and 6.6 months; P = 0.92).

CONCLUSION: Our findings suggest that CCND1A870G may be useful in predicting clinical outcome in HER2-positive mBC patients treated with lapatinib plus capecitabine.

Association of common gene variants in the WNT/β-catenin pathway with colon cancer recurrence

Wnt/β-catenin signaling has a central role in the development and progression of most colon cancers (CCs). Germline variants in Wnt/β-catenin pathway genes may result in altered gene function and/or activity, thereby causing inter-individual differences in relation to tumor recurrence capacity and chemoresistance. We investigated germline polymorphisms in a comprehensive panel of Wnt/β-catenin pathway genes to predict time to tumor recurrence (TTR) in patients with stage III and high-risk stage II CC. A total of 234 patients treated with 5-fluorouracil-based chemotherapy were included in this study. Whole-blood samples were analyzed for putative functional germline polymorphisms in SFRP3, SFRP4, DKK2, DKK3, Axin2, APC, TCF7L2, WNT5B, CXXC4, NOTCH2 and GLI1 genes by PCR-based restriction fragment-length polymorphism or direct DNA sequencing. Polymorphisms with statistical significance were validated in an independent study cohort. The minor allele of WNT5B rs2010851 T>G was significantly associated with a shorter TTR (10.7 vs 4.9 years; hazard ratio: 2.48; 95% CI, 0.96-6.38; P=0.04) in high-risk stage II CC patients. This result remained significant in multivariate Cox's regression analysis. This study shows that the WNT5B germline variant rs2010851 was significantly identified as a stage-dependent prognostic marker for CC patients after 5-fluorouracil-based adjuvant therapy.

Gender-specific profiling in SCN1A polymorphisms and time-to-recurrence in patients with stage II/III colorectal cancer treated with adjuvant 5-fluoruracil chemotherapy

This study was designed to analyze the gender-related association between SCN1A polymorphisms (voltage-gated sodium channels; α-subunit) and time-to-recurrence (TTR) in patients with colorectal cancer (CRC) treated with 5-fluoruracil (5-FU)-based adjuvant chemotherapy. We enrolled from a prospective database patients with stage II and III CRC treated with adjuvant 5-FU-based chemotherapy. Genotypes for SCN1A rs3812718 and rs229877 were determined by direct DNA sequencing. One hundred twenty-seven males and 107 females were included in the study. In the univariate and multivariate analysis, the shortest TTR was associated with female patients carrying the rs3812718-TT genotype (hazard ratio (HR): 2.26 (95% confidence interval (CI): 0.89, 5.70), P=0.039) but with male patients carrying the rs3812718-CC genotype (HR: 0.49 (95% CI: 0.18, 1.38), P=0.048). For rs229877 the CT genotype was associated with a trend for shorter TTR in both gender populations. The study validated gender-dependent association between genomic SCN1A rs3812718 polymorphism and TTR in CRC patients treated with adjuvant 5-FU-based chemotherapy. This study confirms that voltage-gated Na+ channels may be a potential therapeutic target and a useful predictive biomarker before 5-FU infusion.