Differential regulation of LRRC37A2 in gastric cancer by DNA methylation.

Gastric cancer (GC) is one of the leading types of fatal cancer worldwide. Epigenetic manipulation of cancer cells is a useful tool to better understand gene expression regulatory mechanisms and contributes to the discovery of novel biomarkers. Our research group recently reported a list of 83 genes that are potentially modulated by DNA methylation in GC cell lines. Herein, we further explored the regulation of one of these genes, LRRC37A2, in clinical samples. LRRC37A2 expression was evaluated by RT-qPCR, and DNA methylation was studied using next-generation bisulphite sequencing in 36 GC and paired adjacent nonneoplastic tissue samples. We showed that both reduced LRRC37A2 mRNA levels and increased LRRC37A2 exon methylation were associated with undifferentiated and poorly differentiated tumours. Moreover, LRRC37A2 gene expression and methylation levels were inversely correlated at the +45 exon CpG site. We suggest that DNA hypermethylation may contribute to reducing LRRC37A2 expression in undifferentiated and poorly differentiated GC. Therefore, our results show how some genes may be useful to stratify patients who are more likely to benefit from epigenetic therapy.Abbreviations: AR: androgen receptor; 5-AZAdC: 5-aza-2'-deoxycytidine; B2M: beta-2-microglobulin; GAPDH: glyceraldehyde-3-phosphate dehydrogenase; GC: gastric cancer; GLM: general linear model; LRRC37A2: leucine-rich repeat containing 37 member A2; SD: standard deviation; TFII-I: general transcription factor II-I; TSS: transcription start site; XBP1: X-box binding protein 1.

Multicenter Performance Evaluation of Elecsys Anti-HBc II, Anti-HCV II, HIV combi PT, HBsAg II, and Syphilis Immunoassays.

BACKGROUND: The WHO recommends mandatory serological testing of blood donors for hepatitis B virus, hepatitis C virus (HCV), human immunodeficiency virus (HIV), and syphilis. We evaluated the performance of Elecsys® infectious disease immunoassays against commercially available comparator assays. METHODS: Prospective, routine, anonymized patient or donor samples (n = 8,821) were analyzed at three German sites using Elecsys antihepatitis B core antigen (Anti-HBc II), Anti-HCV II, HIV combi PT, hepatitis B surface antigen (HBsAg II), and Syphilis immunoassays (cobas e 411 analyzer) versus ARCHITECT comparator assays. RESULTS: The Elecsys immunoassays demonstrated comparable sensitivity (≤ 1.54% difference) and equivalent specificity (≤ 0.63% difference) to the respective ARCHITECT comparator assays. Overall sensitivity for the Elecsys and ARCHITECT infectious disease panels was 99.78% vs. 99.40%, respectively, and overall specificity was 99.74% vs. 99.80%, respectively. CONCLUSIONS: The Elecsys infectious disease immunoassays demonstrated high sensitivity and specificity, which were similar to comparator assays, supporting their suitability for routine laboratory practice.

Evidence of two deeply divergent co-existing mitochondrial genomes in the Tuatara reveals an extremely complex genomic organization.

Animal mitochondrial genomic polymorphism occurs as low-level mitochondrial heteroplasmy and deeply divergent co-existing molecules. The latter is rare, known only in bivalvian mollusks. Here we show two deeply divergent co-existing mt-genomes in a vertebrate through genomic sequencing of the Tuatara (Sphenodon punctatus), the sole-representative of an ancient reptilian Order. The two molecules, revealed using a combination of short-read and long-read sequencing technologies, differ by 10.4% nucleotide divergence. A single long-read covers an entire mt-molecule for both strands. Phylogenetic analyses suggest a 7-8 million-year divergence between genomes. Contrary to earlier reports, all 37 genes typical of animal mitochondria, with drastic gene rearrangements, are confirmed for both mt-genomes. Also unique to vertebrates, concerted evolution drives three near-identical putative Control Region non-coding blocks. Evidence of positive selection at sites linked to metabolically important transmembrane regions of encoded proteins suggests these two mt-genomes may confer an adaptive advantage for an unusually cold-tolerant reptile.

Discovery of Molecular DNA Methylation-Based Biomarkers through Genome-Wide Analysis of Response Patterns to BCG for Bladder Cancer.

BACKGROUND: Bacillus Calmette-Guérin (BCG) immunotherapy, the standard adjuvant intravesical therapy for some intermediate and most high-risk non-muscle invasive bladder cancers (NMIBCs), suffers from a heterogenous response rate. Molecular markers to help guide responses are scarce and currently not used in the clinical setting. METHODS: To identify novel biomarkers and pathways involved in response to BCG immunotherapy, we performed a genome-wide DNA methylation analysis of NMIBCs before BCG therapy. Genome-wide DNA methylation profiles of DNA isolated from tumors of 26 BCG responders and 27 failures were obtained using the Infinium MethylationEPIC BeadChip. RESULTS: Distinct DNA methylation patterns were found by genome-wide analysis in the two groups. Differentially methylated CpG sites were predominantly located in gene promoters and gene bodies associated with bacterial invasion of epithelial cells, chemokine signaling, endocytosis, and focal adhesion. In total, 40 genomic regions with a significant difference in methylation between responders and failures were detected. The differential methylation state of six of these regions, localized in the promoters of the genes GPR158, KLF8, C12orf42, WDR44, FLT1, and CHST11, were internally validated by bisulfite-sequencing. GPR158 promoter hypermethylation was the best predictor of BCG failure with an AUC of 0.809 (p-value < 0.001). CONCLUSIONS: Tumors from BCG responders and BCG failures harbor distinct DNA methylation profiles. Differentially methylated DNA regions were detected in genes related to pathways involved in bacterial invasion of cells or focal adhesion. We identified candidate DNA methylation biomarkers that may help to predict patient prognosis after external validation in larger, well-designed cohorts.

The impact of DNA demethylation on the upregulation of the NRN1 and TNFAIP3 genes associated with advanced gastric cancer.

Gastric cancer (GC) is the third leading cause of cancer-related death worldwide. Very few therapeutic options are currently available in this neoplasia. The use of 5-Aza-2'-deoxycytidine (5-AZAdC) was approved for the treatment of myelodysplastic syndromes, and this drug can treat solid tumours at low doses. Epigenetic manipulation of GC cell lines is a useful tool to better understand gene expression regulatory mechanisms for clinical applications. Therefore, we compared the gene expression profile of 5-AZAdC-treated and untreated GC cell lines by a microarray assay. Among the genes identified in this analysis, we selected NRN1 and TNFAIP3 to be evaluated for gene expression by RT-qPCR and DNA methylation by bisulfite DNA next-generation sequencing in 43 and 52 pairs of GC and adjacent non-neoplastic tissue samples, respectively. We identified 83 candidate genes modulated by DNA methylation in GC cell lines. Increased expression of NRN1 and TNFAIP3 was associated with advanced tumours (P < 0.05). We showed that increased NRN1 and TNFAIP3 expression seems to be regulated by DNA demethylation in GC samples: inverse correlations between the mRNA and DNA methylation levels in the promoter of NRN1 (P < 0.05) and the intron of TNFAIP3 (P < 0.05) were detected. Reduced NRN1 promoter methylation was associated with III/IV TNM stage tumours (P = 0.03) and the presence of Helicobacter pylori infection (P = 0.02). The identification of demethylated activated genes in GC may be useful in clinical practice, stratifying patients who are less likely to benefit from 5-AZAdC-based therapies. KEY MESSAGES: Higher expression of NRN1 and TNFAIP3 is associated with advanced gastric cancer (GC). NRN1 promoter hypomethylation contributes to gene upregulation in advanced GC. TNFAIP3 intronic-specific CpG site demethylation contributes to gene upregulation in GC. These findings may be useful to stratify GC patients who are less likely to benefit from DNA demethylating-based therapies.

Gene expression profiling identifies pathways involved in seed maturation of Jatropha curcas.

BACKGROUND: Jatropha curcas, a tropical shrub, is a promising biofuel crop, which produces seeds with high content of oil and protein. To better understand the maturation process of J. curcas seeds and to improve its agronomic performance, a two-step approach was performed in six different maturation stages of seeds: 1) generation of the entire transcriptome of J. curcas seeds using 454-Roche sequencing of a cDNA library, 2) comparison of transcriptional expression levels using a custom Agilent 8x60K oligonucleotide microarray. RESULTS: A total of 793,875 high-quality reads were assembled into 19,382 unique full-length contigs, of which 13,507 could be annotated with Gene Ontology (GO) terms. Microarray data analysis identified 9111 probes (out of 57,842 probes), which were differentially expressed between the six maturation stages. The expression results were validated for 75 selected transcripts based on expression levels, predicted function, pathway, and length. Result from cluster analyses showed that transcripts associated with fatty acid, flavonoid, and phenylpropanoid biosynthesis were over-represented in the early stages, while those of lipid storage were over-represented in the late stages. Expression analyses of different maturation stages of J. curcas seed showed that most changes in transcript abundance occurred between the two last stages, suggesting that the timing of metabolic pathways during seed maturation in J. curcas occurs in late stages. The co-expression results showed that the hubs (CB5-D, CDR1, TT8, DFR, HVA22) with the highest number of edges, associated with fatty acid and flavonoid biosynthesis, are showing a decrease in their expression during seed maturation. Furthermore, seed development and hormone pathways are significantly well connected. CONCLUSION: The obtained results revealed differentially expressed sequences (DESs) regulating important pathways related to seed maturation, which could contribute to the understanding of the complex regulatory network during seed maturation with the focus on lipid, flavonoid and phenylpropanoid biosynthesis. This study provides detailed information on transcriptional changes during J. curcas seed maturation and provides a starting point for a genomic survey of seed quality traits. The results highlighted specific genes and processes relevant to the molecular mechanisms involved in Jatropha seed maturation. These data can also be utilized regarding other Euphorbiaceae species.

PRIMEval: Optimization and screening of multiplex oligonucleotide assays.

The development of multiplex polymerase chain reaction and microarray assays is challenging due to primer dimer formation, unspecific hybridization events, the generation of unspecific by-products, primer depletion, and thus lower amplification efficiencies. We have developed a software workflow with three underlying algorithms that differ in their use case and specificity, allowing the complete in silico evaluation of such assays on user-derived data sets. We experimentally evaluated the method for the prediction of oligonucleotide hybridization events including resulting products and probes, self-dimers, cross-dimers and hairpins at different experimental conditions. The developed method allows explaining the observed artefacts through in silico WGS data and thermodynamic predictions. PRIMEval is available publicly at https://primeval.ait.ac.at.

Autologous induced pluripotent stem cell-derived four-organ-chip.

Microphysiological systems play a pivotal role in progressing toward a global paradigm shift in drug development. Here, we designed a four-organ-chip interconnecting miniaturized human intestine, liver, brain and kidney equivalents. All four organ models were predifferentiated from induced pluripotent stem cells from the same healthy donor and integrated into the microphysiological system. The coculture of the four autologous tissue models in one common medium deprived of tissue specific growth factors was successful over 14-days. Although there were no added growth factors present in the coculture medium, the intestine, liver and neuronal model maintained defined marker expression. Only the renal model was overgrown by coexisting cells and did not further differentiate. This model platform will pave the way for autologous coculture cross-talk assays, disease induction and subsequent drug testing.

EPIC-TABSAT: analysis tool for targeted bisulfite sequencing experiments and array-based methylation studies.

DNA methylation is one of the major epigenetic modifications and has frequently demonstrated its suitability as diagnostic and prognostic biomarker. In addition to chip and sequencing based epigenome wide methylation profiling methods, targeted bisulfite sequencing (TBS) has been established as a cost-effective approach for routine diagnostics and target validation applications. Yet, an easy-to-use tool for the analysis of TBS data in combination with array-based methylation results has been missing. Consequently, we have developed EPIC-TABSAT, a user-friendly web-based application for the analysis of targeted sequencing data that additionally allows the integration of array-based methylation results. The tool can handle multiple targets as well as multiple sequencing files in parallel and covers the complete data analysis workflow from calculation of quality metrics to methylation calling and interactive result presentation. The graphical user interface offers an unprecedented way to interpret TBS data alone or in combination with array-based methylation studies. Together with the computation of target-specific epialleles it is useful in validation, research, and routine diagnostic environments. EPIC-TABSAT is freely accessible to all users at https://tabsat.ait.ac.at/.

The Pattern and Distribution of Induced Mutations in J. curcas Using Reduced Representation Sequencing.

Mutagenesis in combination with Genotyping by Sequencing (GBS) is a powerful tool for introducing variation, studying gene function and identifying causal mutations underlying phenotypes of interest in crop plant genomes. About 400 million paired-end reads were obtained from 82 ethylmethane sulfonate (EMS) induced mutants and 14 wild-type accessions of Jatropha curcas for the detection of Single Nucleotide Polymorphisms (SNPs) and Insertion/Deletions (InDels) by two different approaches (nGBS and ddGBS) on an Illumina HiSeq 2000 sequencer. Using bioinformatics analyses, 1,452 induced SNPs and InDels were identified in coding regions, which were distributed across 995 genes. The predominantly observed mutations were G/C to A/T transitions (64%), while transversions were observed at a lower frequency (36%). Regarding the effect of mutations on gene function, 18% of the mutations were located in intergenic regions. In fact, mutants with the highest number of heterozygous SNPs were found in samples treated with 0.8% EMS for 3 h. Reconstruction of the metabolic pathways showed that in total 16 SNPs were located in six KEGG pathways by nGBS and two pathways by ddGBS. The most highly represented pathways were ether-lipid metabolism and glycerophospholipid metabolism, followed by starch and sucrose metabolism by nGBS and triterpenoid biosynthesis as well as steroid biosynthesis by ddGBS. Furthermore, high genome methylation was observed in J. curcas, which might help to understand the plasticity of the Jatropha genome in response to environmental factors. At last, the results showed that continuously vegetatively propagated tissue is a fast, efficient and accurate method to dissolve chimeras, especially for long-lived plants like J. curcas. Obtained data showed that allelic variations and in silico analyses of gene functions (gene function prediction), which control important traits, could be identified in mutant populations using nGBS and ddGBS. However, the handling of GBS data is more difficult and more challenging than the traditional TILLING strategy in mutated plants, since the Jatropha genome sequence is incomplete, which makes alignment and variant analysis of target sequence reads challenging to perform and interpret. Therefore, providing a complete Jatropha reference genome sequence with high quality should be a priority for any breeding program.

Oli2go: an automated multiplex oligonucleotide design tool.

The success of widely used oligonucleotide-based experiments, ranging from PCR to microarray, strongly depends on an accurate design. The design process involves a number of steps, which use specific parameters to produce high quality oligonucleotides. Oli2go is an efficient, user friendly, fully automated multiplex oligonucleotide design tool, which performs primer and different hybridization probe designs as well as specificity and cross dimer checks in a single run. The main improvement to existing oligonucleotide design web-tools is that oli2go combines multiple steps in an all-in-one solution, where other web applications only accomplish parts of the whole design workflow. Especially, the oli2go specificity check is not only performed against a single species (e.g. mouse), but against bacteria, viruses, fungi, invertebrates, plants, protozoa, archaea and sequences from whole genome shotgun sequence projects and environmental samples, at once. This allows the design of highly specific oligonucleotides in multiplex applications, which is further assured by performing dimer checks not only on the primers themselves, but in an all-against-all fashion. The software is freely accessible to all users at http://oli2go.ait.ac.at/.

Frequency of hepatitis B surface antigen variants (HBsAg) in hepatitis B virus genotype B and C infected East- and Southeast Asian patients: Detection by the Elecsys® HBsAg II assay.

BACKGROUND: To avoid false negative results, hepatitis B surface antigen (HBsAg) assays need to detect samples with mutations in the immunodominant 'a' determinant region, which vary by ethnographic region. OBJECTIVE: We evaluated the prevalence and type of HBsAg mutations in a hepatitis B virus (HBV)-infected East- and Southeast Asian population, and the diagnostic performance of the Elecsys® HBsAg II Qualitative assay. STUDY DESIGN: We analyzed 898 samples from patients with HBV infection from four sites (China [Beijing and Guangzhou], Korea and Vietnam). HBsAg mutations were detected and sequenced using highly sensitive ultra-deep sequencing and compared between the first (amino acids 124-137) and second (amino acids 139-147) loops of the 'a' determinant region using the Elecsys® HBsAg II Qualitative assay. RESULTS: Overall, 237 distinct amino acid mutations in the major hydrophilic region were identified; mutations were present in 660 of 898 HBV-infected patient samples (73.5%). Within the pool of 237 distinct mutations, the majority of the amino acid mutations were found in HBV genotype C (64.8%). We identified 25 previously unknown distinct mutations, mostly prevalent in genotype C-infected Korean patients (n = 18) followed by Chinese (n = 12) patients. All 898 samples were correctly identified by the Elecsys® HBsAg II Qualitative assay. CONCLUSIONS: We observed 237 distinct (including 25 novel) mutations, demonstrating the complexity of HBsAg variants in HBV-infected East- and Southeast Asian patients. The Elecsys® HBsAg II Qualitative assay can reliably detect HBV-positive samples and is suitable for routine diagnostic use in East and Southeast Asia.

Substantial variation in the hepatitis B surface antigen (HBsAg) in hepatitis B virus (HBV)-positive patients from South Africa: Reliable detection of HBV by the Elecsys HBsAg II assay.

BACKGROUND: It is essential that hepatitis B surface antigen (HBsAg) diagnostic assays reliably detect genetic diversity in the major hydrophilic region (MHR) of HBsAg to avoid false-negative results. Mutations in this domain display marked ethno-geographic variation and may lead to failure to diagnose hepatitis B virus (HBV) infection. OBJECTIVES: Evaluate diagnostic performance of the Elecsys® HBsAg II Qualitative assay in a cohort of South African HBV-positive blood donors. STUDY DESIGN: A total of 179 South African HBsAg- and HBV DNA > 100 IU/mL-positive blood donor samples were included. Samples were sequenced for genetic variation in HBsAg MHR using next-generation ultra-deep sequencing. HBsAg seropositivity was determined using the Roche Elecsys HBsAg II Qualitative assay. Mutation rates were compared between the first (amino acids 124-137) and second (amino acids 139-147) loops of the immunodominant MHR 'a' determinant region. Frequency of occult HBV infection-associated Y100C mutations was also determined. RESULTS: We observed a total of 279 MHR mutations (117 variants) in 102 (57%) samples, of which 91 were located in the 'a' determinant region. The major vaccine-induced escape mutation G145R was observed in two samples. All occult HBV infection-associated Y100C and common diagnostic and vaccine-escape-associated P120T, G145R, K122R, M133L, M133T, Q129H, G130N, and T126S mutations were reliably detected by the assay, which consistently detected the presence of HBsAg in all 179 samples including samples with 11 novel mutations. CONCLUSIONS: Despite substantial variation in HBsAg MHR, the Elecsys HBsAg II Qualitative assay robustly detects HBV infection in this South African cohort.

Ultra-deep sequencing reveals high prevalence and broad structural diversity of hepatitis B surface antigen mutations in a global population.

The diversity of the hepatitis B surface antigen (HBsAg) has a significant impact on the performance of diagnostic screening tests and the clinical outcome of hepatitis B infection. Neutralizing or diagnostic antibodies against the HBsAg are directed towards its highly conserved major hydrophilic region (MHR), in particular towards its "a" determinant subdomain. Here, we explored, on a global scale, the genetic diversity of the HBsAg MHR in a large, multi-ethnic cohort of randomly selected subjects with HBV infection from four continents. A total of 1553 HBsAg positive blood samples of subjects originating from 20 different countries across Africa, America, Asia and central Europe were characterized for amino acid variation in the MHR. Using highly sensitive ultra-deep sequencing, we found 72.8% of the successfully sequenced subjects (n = 1391) demonstrated amino acid sequence variation in the HBsAg MHR. This indicates that the global variation frequency in the HBsAg MHR is threefold higher than previously reported. The majority of the amino acid mutations were found in the HBV genotypes B (28.9%) and C (25.4%). Collectively, we identified 345 distinct amino acid mutations in the MHR. Among these, we report 62 previously unknown mutations, which extends the worldwide pool of currently known HBsAg MHR mutations by 22%. Importantly, topological analysis identified the "a" determinant upstream flanking region as the structurally most diverse subdomain of the HBsAg MHR. The highest prevalence of "a" determinant region mutations was observed in subjects from Asia, followed by the African, American and European cohorts, respectively. Finally, we found that more than half (59.3%) of all HBV subjects investigated carried multiple MHR mutations. Together, this worldwide ultra-deep sequencing based genotyping study reveals that the global prevalence and structural complexity of variation in the hepatitis B surface antigen have, to date, been significantly underappreciated.

DaMold: A data-mining platform for variant annotation and visualization in molecular diagnostics research.

Next-generation sequencing (NGS) has become a powerful and efficient tool for routine mutation screening in clinical research. As each NGS test yields hundreds of variants, the current challenge is to meaningfully interpret the data and select potential candidates. Analyzing each variant while manually investigating several relevant databases to collect specific information is a cumbersome and time-consuming process, and it requires expertise and familiarity with these databases. Thus, a tool that can seamlessly annotate variants with clinically relevant databases under one common interface would be of great help for variant annotation, cross-referencing, and visualization. This tool would allow variants to be processed in an automated and high-throughput manner and facilitate the investigation of variants in several genome browsers. Several analysis tools are available for raw sequencing-read processing and variant identification, but an automated variant filtering, annotation, cross-referencing, and visualization tool is still lacking. To fulfill these requirements, we developed DaMold, a Web-based, user-friendly tool that can filter and annotate variants and can access and compile information from 37 resources. It is easy to use, provides flexible input options, and accepts variants from NGS and Sanger sequencing as well as hotspots in VCF and BED formats. DaMold is available as an online application at http://damold.platomics.com/index.html, and as a Docker container and virtual machine at https://sourceforge.net/projects/damold/.

MSP-HTPrimer: a high-throughput primer design tool to improve assay design for DNA methylation analysis in epigenetics.

BACKGROUND: Bisulfite (BS) conversion-based and methylation-sensitive restriction enzyme (MSRE)-based PCR methods have been the most commonly used techniques for locus-specific DNA methylation analysis. However, both methods have advantages and limitations. Thus, an integrated approach would be extremely useful to quantify the DNA methylation status successfully with great sensitivity and specificity. Designing specific and optimized primers for target regions is the most critical and challenging step in obtaining the adequate DNA methylation results using PCR-based methods. Currently, no integrated, optimized, and high-throughput methylation-specific primer design software methods are available for both BS- and MSRE-based methods. Therefore an integrated, powerful, and easy-to-use methylation-specific primer design pipeline with great accuracy and success rate will be very useful. RESULTS: We have developed a new web-based pipeline, called MSP-HTPrimer, to design primers pairs for MSP, BSP, pyrosequencing, COBRA, and MSRE assays on both genomic strands. First, our pipeline converts all target sequences into bisulfite-treated templates for both forward and reverse strand and designs all possible primer pairs, followed by filtering for single nucleotide polymorphisms (SNPs) and known repeat regions. Next, each primer pairs are annotated with the upstream and downstream RefSeq genes, CpG island, and cut sites (for COBRA and MSRE). Finally, MSP-HTPrimer selects specific primers from both strands based on custom and user-defined hierarchical selection criteria. MSP-HTPrimer produces a primer pair summary output table in TXT and HTML format for display and UCSC custom tracks for resulting primer pairs in GTF format. CONCLUSIONS: MSP-HTPrimer is an integrated, web-based, and high-throughput pipeline and has no limitation on the number and size of target sequences and designs MSP, BSP, pyrosequencing, COBRA, and MSRE assays. It is the only pipeline, which automatically designs primers on both genomic strands to increase the success rate. It is a standalone web-based pipeline, which is fully configured within a virtual machine and thus can be readily used without any configuration. We have experimentally validated primer pairs designed by our pipeline and shown a very high success rate of primer pairs: out of 66 BSP primer pairs, 63 were successfully validated without any further optimization step and using the same qPCR conditions. The MSP-HTPrimer pipeline is freely available from http://sourceforge.net/p/msp-htprimer.

Erratum to: MSRE-HTPrimer: a high-throughput and genome-wide primer design pipeline optimized for epigenetic research.

[This corrects the article DOI: 10.1186/s13148-016-0190-9.].

Analysis and Visualization Tool for Targeted Amplicon Bisulfite Sequencing on Ion Torrent Sequencers.

Targeted sequencing of PCR amplicons generated from bisulfite deaminated DNA is a flexible, cost-effective way to study methylation of a sample at single CpG resolution and perform subsequent multi-target, multi-sample comparisons. Currently, no platform specific protocol, support, or analysis solution is provided to perform targeted bisulfite sequencing on a Personal Genome Machine (PGM). Here, we present a novel tool, called TABSAT, for analyzing targeted bisulfite sequencing data generated on Ion Torrent sequencers. The workflow starts with raw sequencing data, performs quality assessment, and uses a tailored version of Bismark to map the reads to a reference genome. The pipeline visualizes results as lollipop plots and is able to deduce specific methylation-patterns present in a sample. The obtained profiles are then summarized and compared between samples. In order to assess the performance of the targeted bisulfite sequencing workflow, 48 samples were used to generate 53 different Bisulfite-Sequencing PCR amplicons from each sample, resulting in 2,544 amplicon targets. We obtained a mean coverage of 282X using 1,196,822 aligned reads. Next, we compared the sequencing results of these targets to the methylation level of the corresponding sites on an Illumina 450k methylation chip. The calculated average Pearson correlation coefficient of 0.91 confirms the sequencing results with one of the industry-leading CpG methylation platforms and shows that targeted amplicon bisulfite sequencing provides an accurate and cost-efficient method for DNA methylation studies, e.g., to provide platform-independent confirmation of Illumina Infinium 450k methylation data. TABSAT offers a novel way to analyze data generated by Ion Torrent instruments and can also be used with data from the Illumina MiSeq platform. It can be easily accessed via the Platomics platform, which offers a web-based graphical user interface along with sample and parameter storage. TABSAT is freely available under a GNU General Public License version 3.0 (GPLv3) at https://github.com/tadkeys/tabsat/ and http://demo.platomics.com/.

MSRE-HTPrimer: a high-throughput and genome-wide primer design pipeline optimized for epigenetic research.

BACKGROUND: Methylation-sensitive restriction enzymes-polymerase chain reaction (MSRE-PCR) has been used in epigenetic research to identify genome-wide and gene-specific DNA methylation. Currently, epigenome-wide discovery studies provide many candidate regions for which the MSREqPCR approach can be very effective to confirm the findings. MSREqPCR provides high multiplexing capabilities also when starting with limited amount of DNA-like cfDNA to validate many targets in a time- and cost-effective manner. Multiplex design is challenging and cumbersome to define specific primers in an effective manner, and no suitable software tools are freely available for high-throughput primer design in a time-effective manner and to automatically annotate the resulting primers with known SNPs, CpG, repeats, and RefSeq genes. Therefore a robust, powerful, high-throughput, optimized, and methylation-specific primer design tool with great accuracy will be very useful. RESULTS: We have developed a novel pipeline, called MSRE-HTPrimer, to design MSRE-PCR and genomic PCR primers pairs in a very efficient manner and with high success rate. First, our pipeline designs all possible PCR primer pairs and oligos, followed by filtering for SNPs loci and repeat regions. Next, each primer pair is annotated with the number of cut sites in primers and amplicons, upstream and downstream genes, and CpG islands loci. Finally, MSRE-HTPrimer selects resulting primer pairs for all target sequences based on a custom quality matrix defined by the user. MSRE-HTPrimer produces a table for all resulting primer pairs as well as a custom track in GTF file format for each target sequence to visualize it in UCSC genome browser. CONCLUSIONS: MSRE-HTPrimer, based on Primer3, is a high-throughput pipeline and has no limitation on the number and size of target sequences for primer design and provides full flexibility to customize it for specific requirements. It is a standalone web-based pipeline, which is fully configured within a virtual machine and thus can be readily used without any configuration. We have experimentally validated primer pairs designed by our pipeline and shown a very high success rate of primer pairs: out of 190 primer pairs, 71 % could be successfully validated. The MSRE-HTPrimer software is freely available from http://sourceforge.net/p/msrehtprimer/wiki/Virtual_Machine/ as a virtual machine.

ClinQC: a tool for quality control and cleaning of Sanger and NGS data in clinical research.

BACKGROUND: Traditional Sanger sequencing has been used as a gold standard method for genetic testing in clinic to perform single gene test, which has been a cumbersome and expensive method to test several genes in heterogeneous disease such as cancer. With the advent of Next Generation Sequencing technologies, which produce data on unprecedented speed in a cost effective manner have overcome the limitation of Sanger sequencing. Therefore, for the efficient and affordable genetic testing, Next Generation Sequencing has been used as a complementary method with Sanger sequencing for disease causing mutation identification and confirmation in clinical research. However, in order to identify the potential disease causing mutations with great sensitivity and specificity it is essential to ensure high quality sequencing data. Therefore, integrated software tools are lacking which can analyze Sanger and NGS data together and eliminate platform specific sequencing errors, low quality reads and support the analysis of several sample/patients data set in a single run. RESULTS: We have developed ClinQC, a flexible and user-friendly pipeline for format conversion, quality control, trimming and filtering of raw sequencing data generated from Sanger sequencing and three NGS sequencing platforms including Illumina, 454 and Ion Torrent. First, ClinQC convert input read files from their native formats to a common FASTQ format and remove adapters, and PCR primers. Next, it split bar-coded samples, filter duplicates, contamination and low quality sequences and generates a QC report. ClinQC output high quality reads in FASTQ format with Sanger quality encoding, which can be directly used in down-stream analysis. It can analyze hundreds of sample/patients data in a single run and generate unified output files for both Sanger and NGS sequencing data. Our tool is expected to be very useful for quality control and format conversion of Sanger and NGS data to facilitate improved downstream analysis and mutation screening. CONCLUSIONS: ClinQC is a powerful and easy to handle pipeline for quality control and trimming in clinical research. ClinQC is written in Python with multiprocessing capability, run on all major operating systems and is available at https://sourceforge.net/projects/clinqc.