Large scale read classification for next generation sequencing

Next Generation Sequencing (NGS) has revolutionised molecular biology, resulting in an explosion of data sets and an increasing role in clinical practice. Such applications necessarily require rapid identification of the organism as a prelude to annotation and further analysis. NGS data consist of a substantial number of short sequence reads, given context through downstream assembly and annotation, a process requiring reads consistent with the assumed species or species group. Highly accurate results have been obtained for restricted sets using SVM classifiers, but such methods are difficult to parallelise and success depends on careful attention to feature selection. This work examines the problem at very large scale, using a mix of synthetic and real data with a view to determining the overall structure of the problem and the effectiveness of parallel ensembles of simpler classifiers (principally random forests) in addressing the challenges of large scale genomics.

Correlation between BRAF mutation and the clinicopathological parameters in papillary thyroid carcinoma with particular reference to follicular variant

Mutation of the BRAF gene is common in thyroid cancer. Follicular variant of papillary thyroid carcinoma is a variant of papillary thyroid carcinoma that has created continuous diagnostic controversies among pathologists. The aims of this study are to (1) investigate whether follicular variant of papillary thyroid carcinoma has a different pattern of BRAF mutation than conventional papillary thyroid carcinoma in a large cohort of patients with typical features of follicular variant of papillary thyroid carcinoma and (2) to study the relationship of clinicopathological features of papillary thyroid carcinomas with BRAF mutation. Tissue blocks from 76 patients with diagnostic features of papillary thyroid carcinomas (40 with conventional type and 36 with follicular variant) were included in the study. From these, DNA was extracted and BRAF V600E mutations were detected by polymerase chain reaction followed by restriction enzyme digestion and sequencing of exon 15. Analysis of the data indicated that BRAF V600E mutation is significantly more common in conventional papillary thyroid carcinoma (58% versus 31%, P = .022). Furthermore, the mutation was often noted in female patients (P = .017), in high-stage cancers (P = .034), and in tumors with mild lymphocytic thyroiditis (P = .006). We concluded that follicular variant of papillary thyroid carcinoma differs from conventional papillary thyroid carcinoma in the rate of BRAF mutation. The results of this study add further information indicating that mutations in BRAF play a role in thyroid cancer development and progression.

Offspring of mothers fed a high fat diet display hepatic cell cycle inhibition and associated changes in gene expression and DNA methylation

The association between an adverse early life environment and increased susceptibility to later-life metabolic disorders such as obesity, type 2 diabetes and cardiovascular disease is described by the developmental origins of health and disease hypothesis. Employing a rat model of maternal high fat (MHF) nutrition, we recently reported that offspring born to MHF mothers are small at birth and develop a postnatal phenotype that closely resembles that of the human metabolic syndrome. Livers of offspring born to MHF mothers also display a fatty phenotype reflecting hepatic steatosis and characteristics of non-alcoholic fatty liver disease. In the present study we hypothesised that a MHF diet leads to altered regulation of liver development in offspring; a derangement that may be detectable during early postnatal life. Livers were collected at postnatal days 2 (P2) and 27 (P27) from male offspring of control and MHF mothers (n = 8 per group). Cell cycle dynamics, measured by flow cytometry, revealed significant G0/G1 arrest in the livers of P2 offspring born to MHF mothers, associated with an increased expression of the hepatic cell cycle inhibitor Cdkn1a. In P2 livers, Cdkn1a was hypomethylated at specific CpG dinucleotides and first exon in offspring of MHF mothers and was shown to correlate with a demonstrable increase in mRNA expression levels. These modifications at P2 preceded observable reductions in liver weight and liver:brain weight ratio at P27, but there were no persistent changes in cell cycle dynamics or DNA methylation in MHF offspring at this time. Since Cdkn1a up-regulation has been associated with hepatocyte growth in pathologic states, our data may be suggestive of early hepatic dysfunction in neonates born to high fat fed mothers. It is likely that these offspring are predisposed to long-term hepatic dysfunction.

miRPlant : an integrated tool for identification of plant miRNA from RNA sequencing data

Background Small RNA sequencing is commonly used to identify novel miRNAs and to determine their expression levels in plants. There are several miRNA identification tools for animals such as miRDeep, miRDeep2 and miRDeep*. miRDeep-P was developed to identify plant miRNA using miRDeep’s probabilistic model of miRNA biogenesis, but it depends on several third party tools and lacks a user-friendly interface. The objective of our miRPlant program is to predict novel plant miRNA, while providing a user-friendly interface with improved accuracy of prediction. Result We have developed a user-friendly plant miRNA prediction tool called miRPlant. We show using 16 plant miRNA datasets from four different plant species that miRPlant has at least a 10% improvement in accuracy compared to miRDeep-P, which is the most popular plant miRNA prediction tool. Furthermore, miRPlant uses a Graphical User Interface for data input and output, and identified miRNA are shown with all RNAseq reads in a hairpin diagram. Conclusions We have developed miRPlant which extends miRDeep* to various plant species by adopting suitable strategies to identify hairpin excision regions and hairpin structure filtering for plants. miRPlant does not require any third party tools such as mapping or RNA secondary structure prediction tools. miRPlant is also the first plant miRNA prediction tool that dynamically plots miRNA hairpin structure with small reads for identified novel miRNAs. This feature will enable biologists to visualize novel pre-miRNA structure and the location of small RNA reads relative to the hairpin. Moreover, miRPlant can be easily used by biologists with limited bioinformatics skills.

Integrated analysis of epigenomic and genomic changes by DNA methylation dependent mechanisms provides potential novel biomarkers for prostate cancer

Epigenetic silencing mediated by CpG methylation is a common feature of many cancers. Characterizing aberrant DNA methylation changes associated with tumor progression may identify potential prognostic markers for prostate cancer (PCa). We treated two PCa cell lines, 22Rv1 and DU-145 with the demethylating agent 5-Aza 2’–deoxycitidine (DAC) and global methylation status was analyzed by performing methylation-sensitive restriction enzyme based differential methylation hybridization strategy followed by genome-wide CpG methylation array profiling. In addition, we examined gene expression changes using a custom microarray. Gene Set Enrichment Analysis (GSEA) identified the most significantly dysregulated pathways. In addition, we assessed methylation status of candidate genes that showed reduced CpG methylation and increased gene expression after DAC treatment, in Gleason score (GS) 8 vs. GS6 patients using three independent cohorts of patients; the publically available The Cancer Genome Atlas (TCGA) dataset, and two separate patient cohorts. Our analysis, by integrating methylation and gene expression in PCa cell lines, combined with patient tumor data, identified novel potential biomarkers for PCa patients. These markers may help elucidate the pathogenesis of PCa and represent potential prognostic markers for PCa patients.

DNA adduction by the potent carcinogen aflatoxin B1 : mechanistic studies

Aflatoxin B1, a potently carcinogenic fungal metabolite, is converted to the biologically active form by chemical oxidation using dimethyldioxirane and enzymatically by cytochrome P450 mixed-function oxidases. Both processes give rise to mixtures of the exo- and endo-8,9-epoxides. Methanolysis studies reveal exclusive trans opening of both epoxides under neutral conditions in CH3OH and CH3OH/H2O mixtures; an SN2 mechanism is postulated. Under acidic conditions, the exo isomer gives mixtures of trans and cis solvolysis products, suggesting that the reaction is, at least in part, SN1; the endo isomer gives only the trans product. The exo isomer reacts with DNA by attack of the nitrogen atom at the 7 position of guanine on C8 of the epoxide to give the trans adduct; the endo epoxide fails to form an adduct at this or any other site in DNA. The exo isomer is strongly mutagenic in a base-pair reversion assay employing Salmonella typhimurium; the endo isomer is essentially nonmutagenic. Aflatoxin B1 and its derivatives intercalate in DNA. These results are consistent with a mechanism in which intercalation of the exo epoxide optimally orients the epoxide for an SN2 reaction with guanine but intercalation of the endo isomer places the epoxide in an orientation which precludes reaction. Thus, while the exo epoxide is a potent mutagen, the endo epoxide fails to react with DNA.

Determination of urinary thymidine glycol using affinity chromatography, HPLC and post-column reaction detection : a biomarker of oxidative DNA damage upon kidney transplantation

Reactive oxygen species are generated during ischaemia-reperfusion of tissue. Oxidation of thymidine by hydroxyl radicals (HO) leads to the formation of 5,6-dihydroxy-5,6-dihydrothymidine (thymidine glycol). Thymidine glycol is excreted in urine and can be used as biomarker of oxidative DNA damage. Time dependent changes in urinary excretion rates of thymidine glycol were determined in six patients after kidney transplantation and in six healthy controls. A new analytical method was developed involving affinity chromatography and subsequent reverse-phase high-performance liquid chromatography (RP-HPLC) with a post-column chemical reaction detector and endpoint fluorescence detection. The detection limit of this fluorimetric assay was 1.6 ng thymidine glycol per ml urine, which corresponds to about half of the physiological excretion level in healthy control persons. After kidney transplantation the urinary excretion rate of thymidine glycol increased gradually reaching a maximum around 48 h. The excretion rate remained elevated until the end of the observation period of 10 days. Severe proteinuria with an excretion rate of up to 7.2 g of total protein per mmol creatinine was also observed immediately after transplantation and declined within the first 24 h of allograft function (0.35 + 0.26 g/mmol creatinine). The protein excretion pattern, based on separation of urinary proteins on sodium dodecyl sulphate-polyacrylamide gel electrophorosis (SDS-PAGE), as well as excretion of individual biomarker proteins, indicated nonselective glomerular and tubular damage. The increased excretion of thymidine glycol after kidney transplantation may be explained by ischaemia-reperfusion induced oxidative DNA damage of the transplanted kidney.

DNA binding study of isophorone in rats and mice

Following isophorone exposure, in a 2-year study with F344 rats and B6C3F1 mice performed under the National Toxicology Program (NTP), an elevated incidence of tumors was observed in male rats (kidney tumors) and male mice (liver tumors). Female rats and mice showed no elevation of tumor rates by isophorone (NTP 1986).

DNA methylation at the novel CpG sites in the promoter of MED15/PCQAP gene as a biomarker for head and neck cancers

Head and neck cancers (HNCs) represent a significant and ever-growing burden to the modern society, mainly due to the lack of early diagnostic methods. A significant number of HNCs is often associated with drinking, smoking, chewing beetle nut, and human papilloma virus (HPV) infections. We have analyzed DNA methylation patterns in tumor and normal tissue samples collected from head and neck squamous cell carcinoma (HNSCC) patients who were smokers. We have identified novel methylation sites in the promoter of the mediator complex subunit 15 (MED15/PCQAP) gene (encoing a co-factor important for regulation of transcription initiation for promoters of many genes), hypermethylated specifically in tumor cells. Two clusters of CpG dinucleotides methylated in tumors, but not in normal tissue from the same patients, were identified. These CpG methylation events in saliva samples were further validated in a separate cohort of HNSCC patients (who developed cancer due to smoking or HPV infections) and healthy controls using methylation-specific PCR (MSP). We used saliva as a biological medium because of its non-invasive nature, close proximity to the tumors, easiness and it is an economically viable option for large-scale screening studies. The methylation levels for the two identified CpG clusters were significantly different between the saliva samples collected from healthy controls and HNSCC individuals (Welch's t-test returning P, 0.05 and Mann-Whitney test P, 0.01 for both). The developed MSP assays also provided a good discriminative ability with AUC values of 0.70 (P, 0.01) and 0.63 (P, 0.05). The identified novel CpG methylation sites may serve as potential non-invasive biomarkers for detecting HNSCC. © the authors.

Epigenetics and migraine; complex mitochondrial interactions contributing to disease susceptibility

Migraine is a common neurological disorder classified by the World Health Organisation (WHO) as one of the top twenty most debilitating diseases in the developed world. Current therapies are only effective for a proportion of sufferers and new therapeutic targets are desperately needed to alleviate this burden. Recently the role of epigenetics in the development of many complex diseases including migraine has become an emerging topic. By understanding the importance of acetylation, methylation and other epigenetic modifications, it then follows that this modification process is a potential target to manipulate epigenetic status with the goal of treating disease. Bisulphite sequencing and methylated DNA immunoprecipitation have been used to demonstrate the presence of methylated cytosines in the human D-loop of mitochondrial DNA (mtDNA), proving that the mitochondrial genome is methylated. For the first time, it has been shown that there is a difference in mtDNA epigenetic status between healthy controls and those with disease, especially for neurodegenerative and age related conditions. Given co-morbidities with migraine and the suggestive link between mitochondrial dysfunction and the lowered threshold for triggering a migraine attack, mitochondrial methylation may be a new avenue to pursue. Creative thinking and new approaches are needed to solve complex problems and a systems biology approach, where multiple layers of information are integrated is becoming more important in complex disease modelling.

Methods for extracting genomic DNA from whole blood samples : current perspectives

Deoxyribonucleic acid (DNA) extraction has considerably evolved since it was initially performed back in 1869. It is the first step required for many of the available downstream applications used in the field of molecular biology. Whole blood samples are one of the main sources used to obtain DNA, and there are many different protocols available to perform nucleic acid extraction on such samples. These methods vary from very basic manual protocols to more sophisticated methods included in automated DNA extraction protocols. Based on the wide range of available options, it would be ideal to determine the ones that perform best in terms of cost-effectiveness and time efficiency. We have reviewed DNA extraction history and the most commonly used methods for DNA extraction from whole blood samples, highlighting their individual advantages and disadvantages. We also searched current scientific literature to find studies comparing different nucleic acid extraction methods, to determine the best available choice. Based on our research, we have determined that there is not enough scientific evidence to support one particular DNA extraction method from whole blood samples. Choosing a suitable method is still a process that requires consideration of many different factors, and more research is needed to validate choices made at facilities around the world.

Detection of DNA mutations using novel SERS (Surface-Enhanced Raman Spectroscopy) diagnostic platform

This article describes the detection of DNA mutations using novel Au-Ag coated GaN substrate as SERS (surface-enhanced Raman spectroscopy) diagnostic platform. Oligonucleotide sequences corresponding to the BCR-ABL (breakpoint cluster region-Abelson) gene responsible for development of chronic myelogenous leukemia were used as a model system to demonstrate the discrimination between the wild type and Met244Val mutations. The thiolated ssDNA (single-strand DNA) was immobilized on the SERS-active surface and then hybridized to a labeled target sequence from solution. An intense SERS signal of the reporter molecule MGITC was detected from the complementary target due to formation of double helix. The SERS signal was either not observed, or decreased dramatically for a negative control sample consisting of labeled DNA that was not complementary to the DNA probe. The results indicate that our SERS substrate offers an opportunity for the development of novel diagnostic assays.