Technical Considerations for Reduced Representation Bisulfite Sequencing with Multiplexed Libraries

Reduced representation bisulfite sequencing (RRBS), which couples bisulfite conversion and next generation sequencing, is an innovative method that specifically enriches genomic regions with a high density of potential methylation sites and enables investigation of DNA methylation at single-nucleotide resolution. Recent advances in the Illumina DNA sample preparation protocol and sequencing technology have vastly improved sequencing throughput capacity. Although the new Illumina technology is now widely used, the unique challenges associated with multiplexed RRBS libraries on this platform have not been previously described. We have made modifications to the RRBS library preparation protocol to sequence multiplexed libraries on a single flow cell lane of the Illumina HiSeq 2000. Furthermore, our analysis incorporates a bioinformatics pipeline specifically designed to process bisulfite-converted sequencing reads and evaluate the output and quality of the sequencing data generated from the multiplexed libraries. We obtained an average of 42 million paired-end reads per sample for each flow-cell lane, with a high unique mapping efficiency to the reference human genome. Here we provide a roadmap of modifications, strategies, and trouble shooting approaches we implemented to optimize sequencing of multiplexed libraries on an a RRBS background

Methodology Report Technical Considerations for Reduced Representation Bisulfite Sequencing with Multiplexed Libraries

Reduced representation bisulfite sequencing (RRBS), which couples bisulfite conversion and next generation sequencing, is an innovative method that specifically enriches genomic regions with a high density of potential methylation sites and enables investigation of DNA methylation at single-nucleotide resolution. Recent advances in the Illumina DNA sample preparation protocol and sequencing technology have vastly improved sequencing throughput capacity. Although the new Illumina technology is now widely used, the unique challenges associated with multiplexed RRBS libraries on this platform have not been previously described. We have made modifications to the RRBS library preparation protocol to sequence multiplexed libraries on a single flow cell lane of the Illumina HiSeq 2000. Furthermore, our analysis incorporates a bioinformatics pipeline specifically designed to process bisulfite-converted sequencing reads and evaluate the output and quality of the sequencing data generated from the multiplexed libraries. We obtained an average of 42 million paired-end reads per sample for each flow-cell lane, with a high unique mapping efficiency to the reference human genome. Here we provide a roadmap of modifications, strategies, and trouble shooting approaches we implemented to optimize sequencing of multiplexed libraries on an a RRBS background

Assessment of a size-based method for enriching circulating tumour cells in colorectal cancer

Circulating tumour cells (CTC) from solid tumours are a prerequisite for metastasis. Isolating CTCs and understanding their biology is essential for developing new clinical tests and precision oncology. Currently, CellSearch is the only FDA (U.S. Food and Drug Administration)-approved method for CTC enrichment but possesses several drawbacks owing to a reliance on the epithelial cell adhesion molecule (EpCAM) and a resource-intensive nature. Addressing these shortcomings, we optimised an existing size-based method, MetaCell, to enrich CTCs from blood of colorectal cancer (CRC) patients. We evaluated the ability of MetaCell to enrich CTCs by spiking blood with CRC cell lines and assessing the cell recovery rates and WBC depletion via immunostaining and gene expression. We then applied MetaCell to samples from 17 CRC patients and seven controls. Recovery rates were \u3e85% in cell lines, with \u3e95% depletion in WBCs. MetaCell yielded CTCs and CTC clusters in 52.9% and 23.5% of the patients, respectively, without false positives in control patients. CTCs and cluster detection did not correlate with histopathological parameters. Overall, we demonstrated that the MetaCell platform enriched CRC cells with high recovery rates and high purity. Our pilot study also demonstrated the ability of MetaCell to detect CTCs in CRC patients

Comparison of alignment software for genome-wide bisulphite sequence data

Recent advances in next generation sequencing (NGS) technology now provide the opportunity to rapidly interrogate the methylation status of the genome. However, there are challenges in handling and interpretation of the methylation sequence data because of its large volume and the consequences of bisulphite modification. We sequenced reduced representation human genomes on the Illumina platform and efficiently mapped and visualized the data with different pipelines and software packages. We examined three pipelines for aligning bisulphite converted sequencing reads and compared their performance. We also comment on pre-processing and quality control of Illumina data. This comparison highlights differences in methods for NGS data processing and provides guidance to advance sequence-based methylation data analysis for molecular biologists

Proteomic Analysis of Aortae from Human Lipoprotein(a) Transgenic Mice Shows an Early Metabolic Response Independent of Atherosclerosis

Background: Elevated low density lipoprotein (LDL) and lipoprotein(a) are independent risk factors for the development of atherosclerosis. Using a proteomic approach we aimed to determine early changes in arterial protein expression in transgenic mice containing both human LDL and lipoprotein(a) in circulation. Methods and Results: Plasma lipid analyses showed the lipoprotein(a) transgenic mice had significantly higher lipid levels than wildtype, including a much increased LDL and high density lipoprotein (HDL) cholesterol. Analysis of aortae from lipoprotein(a) mice showed lipoprotein(a) accumulation but no lipid accumulation or foam cells, leaving the arteries essentially atherosclerosis free. Using two-dimensional gel electrophoresis and mass spectrometry, we identified 34 arterial proteins with significantly altered abundance (P,0.05) in lipoprotein(a) transgenic mice compared to wildtype including 17 that showed a $2 fold difference. Some proteins of interest showed a similarly altered abundance at the transcript level. These changes collectively indicated an initial metabolic response that included a down regulation in energy, redox and lipid metabolism proteins and changes in structural proteins at a stage when atherosclerosis had not yet developed. Conclusions: Our study shows that human LDL and lipoprotein(a) promote changes in the expression of a unique set o

Sex differences in DNA methylation and expression in zebrafish brain: a test of an extended ‘male sex drive’ hypothesis

The sex drive hypothesis predicts that stronger selection on male traits has resulted in masculinization of the genome. Here we test whether such masculinizing effects can be detected at the level of the transcriptome and methylome in the adult zebrafish brain. Although methylation is globally similar, we identified 914 specific differentially methylated CpGs (DMCs) between males and females (435 were hypermethylated and 479 were hypomethylated in males compared to females). These DMCs were prevalent in gene body, intergenic regions and CpG island shores. We also discovered 15 distinct CpG clusters with striking sex-specific DNA methylation differences. In contrast, at transcriptome level, more female-biased genes than male-biased genes were expressed, giving little support for the male sex drive hypothesis. Our study provides genome-wide methylome and transcriptome assessment and sheds light on sex-specific epigenetic patterns and in zebrafish for the first time

Comparison of Global DNA Methylation Patterns in Human Melanoma Tissues and Their Derivative Cell Lines

DNA methylation is a heritable epigenetic mark that is fundamental to mammalian development. Aberrant DNA methylation is an epigenetic hallmark of cancer cells. Cell lines are a critical in vitro model and very widely used to unravel mechanisms of cancer cell biology. However, limited data are available to assess whether DNA methylation patterns in tissues are retained when cell lines are established. Here, we provide the first genome-scale sequencing-based methylation map of metastatic melanoma tumour tissues and their derivative cell lines. We show that DNA methylation profiles are globally conserved in vitro compared to the tumour tissue of origin. However, we identify sites that are consistently hypermethylated in cell lines compared to their tumour tissue of origin. The genes associated with these common differentially methylated regions are involved in cell metabolism, cell cycle and apoptosis and are also strongly enriched for the H3K27me3 histone mark and PRC2 complex-related genes. Our data indicate that although global methylation patterns are similar between tissues and cell lines, there are site-specific epigenomic differences that could potentially impact gene expression. Our work provides a valuable resource for identifying false positives due to cell culture and for better interpretation of cancer epigenetics studies in the future

Single-Cell DNA Methylation Analysis in Cancer

Morphological, transcriptomic, and genomic defects are well-explored parameters of cancer biology. In more recent years, the impact of epigenetic influences, such as DNA methylation, is becoming more appreciated. Aberrant DNA methylation has been implicated in many types of cancers, influencing cell type, state, transcriptional regulation, and genomic stability to name a few. Traditionally, large populations of cells from the tissue of interest are coalesced for analysis, producing averaged methylome data. Considering the inherent heterogeneity of cancer, analysing populations of cells as a whole denies the ability to discover novel aberrant methylation patterns, identify subpopulations, and trace cell lineages. Due to recent advancements in technology, it is now possible to obtain methylome data from single cells. This has both research and clinical implications, ranging from the identification of biomarkers to improved diagnostic tools. As with all emerging technologies, distinct experimental, bioinformatic, and practical challenges present themselves. This review begins with exploring the potential impact of single-cell sequencing on understanding cancer biology and how it could eventually benefit a clinical setting. Following this, the techniques and experimental approaches which made this technology possible are explored. Finally, the present challenges currently associated with single-cell DNA methylation sequencing are described

Recent Discoveries in Epigenetic Modifications of Polycystic Kidney Disease

Autosomal Dominant Polycystic Kidney Disease (ADPKD) is a heritable renal disease that results in end-stage kidney disease, due to the uncontrolled bilateral growth of cysts throughout the kidneys. While it is known that a mutation within a PKD-causing gene is required for the development of ADPKD, the underlying mechanism(s) causing cystogenesis and progression of the disease are not well understood. Limited therapeutic options are currently available to slow the rate of cystic growth. Epigenetic modifications, including DNA methylation, are known to be altered in neoplasia, and several FDA-approved therapeutics target these disease-specific changes. As there are many similarities between ADPKD and neoplasia, we (and others) have postulated that ADPKD kidneys contain alterations to their epigenetic landscape that could be exploited for future therapeutic discovery. Here we summarise the current understanding of epigenetic changes that are associated with ADPKD, with a particular focus on the burgeoning field of ADPKD-specific alterations in DNA methylation

Proteomic Analysis of Liver from Human Lipoprotein(a) Transgenic Mice Shows an Oxidative Stress and Lipid Export Response

Background. Mouse models of hypercholesterolaemia have been used to identify arterial proteins involved in atherosclerosis. As the liver is extremely sensitive to dyslipidemia, one might expect major changes in the abundance of liver proteins in these models even before atherosclerosis develops. Methods. Lipid levels were measured and a proteomic approach was used to quantify proteins in the livers of mice with an elevated low-density lipoprotein (LDL) and the presence of lipoprotein(a) [Lp(a)] but no atherosclerosis. Results. The livers of Lp(a) mice showed an increased triglyceride but reduced phospholipid and oxidised lipid content. Two-dimensional gel electrophoresis and mass spectrometry analysis identified 24 liver proteins with significantly increased abundance in Lp(a) mice (P<0.05). A bioinformatic analysis of the 24 proteins showed the major effect was that of an enhanced antioxidant and lipid efflux response with significant increases in antioxidant (Park7, Gpx1, Prdx6, and Sod1) and lipid metabolism proteins (Fabp4, Acaa2, apoA4, and ApoA1). Interestingly, human liver cells treated with Lp(a) showed significant increases in Gpx1 and Prdx6 but not Sod1 or Park7. Conclusions. The presence of human LDL and Lp(a) in mice promotes an enhanced flux of lipids into the liver which elicits an antioxidant and lipid export response before the onset of atherosclerosis. The antioxidant response can be reproduced in human liver cells treated with Lp(a)