PrimerSuite: A High-Throughput Web-Based Primer Design Program for Multiplex Bisulfite PCR

Typing of Echinococcus multilocularis by Region-Specific Extraction and Next-Generation Sequencing of the mitogenome

Background

Infection by the fox tapeworm Echinococcus multilocularis may lead to a severe zoonosis in humans, alveolar echinococcosis, which may be fatal if left untreated. Typing is important to understand the epidemiology of this parasite, yet there is limited knowledge on the microdiversity of E. multilocularis on the local scale, since the typing resolution of established methods is restricted.

Methods

The mitogenome of E. multilocularis was used as the target regions to modify, apply and validate the Region-Specific Extraction (RSE) method in combination with Next-Generation Sequencing (NGS). Single Nucleotide Polymorphisms (SNPs) were detected in the mitochondrial DNA (mtDNA) and analysed bioinformatically. To validate the success and the accuracy of the RSE protocol, the mitogenomes of some E. multilocularis isolates were also analysed by the Whole-Genome Sequencing (WGS).

Results

With the chosen combination of methods, the entire mitogenome (~13 kb) of E. multilocularis could be captured and amplified. The read depth (median ≥ 156X) was sufficient to detect existing SNPs. The comparison of mitogenome sequences extracted by RSE with mitogenome sequences obtained by WGS showed that the accuracy of the RSE method was consistently comparable to direct Whole-Genome Sequencing.

Conclusion

The results demonstrate that the RSE method in combination with NGS is suitable to analyse the microdiversity of E. multilocularis at the whole mitogenome level. For the capture and sequencing of large (several kb) genomic regions of E. multilocularis and other applications, this method can be very helpful.

Digital PCR-Based Gene Expression Analysis Using a Highly Multiplexed Assay with Universal Detection Probes to Study Induced Pluripotent Stem Cell Differentiation into Cranial Neural Crest Cells

Induced pluripotent stem cells (iPSCs) have the potential to differentiate into any cell type, offering a valuable tool for research in developmental biology, regenerative medicine, and disease modeling. In this study, iPSCs were differentiated into cranial neural crest cells (CNCCs) over a 14-day period. RNA was extracted from these cells at day 0 (iPSCs), day 7, and day 14 to evaluate successful differentiation through the expression analysis of pluripotency and CNCC marker genes.A key focus was the conversion of existing qPCR assays into multiplexed RT-dPCR assays utilizing universal detection probes for precise gene expression analysis during the differentiation of induced pluripotent stem cells (iPSCs) into cranial neural crest cells (CNCCs). We aimed to leverage the superior precision, sensitivity, and multiplexing-degree of dPCR, particularly in quantifying low-abundance targets. We conducted a comparative analysis of the temporal expression patterns of crucial marker genes using both qPCR and dPCR.Our experiments revealed that the four five-plex dPCR assays could successfully detect and quantify the pluripotency and CNCC marker genes and evaluate CNCC differentiation. We observed the expected downregulation of pluripotency genes during differentiation. Conversely, the upregulation of CNCC markers validates the successful differentiation process. In conclusion, SYBR Green I gene expression qPCR assays can be readily converted into multiplex dPCR assays using universal detection probes.Overall, this work underscores the potential of dPCR as a valuable tool for molecular profiling in stem cell research, offering robust, precise, and efficient gene expression analysis. The findings suggest that while qPCR remains a reliable method for routine applications, dPCR provides particular advantages for high-precision, low-sample input studies, expanding the analytical toolbox for stem cell differentiation and gene expression research.

Multiplex digital profiling of DNA methylation heterogeneity for sensitive and cost-effective cancer detection in low-volume liquid biopsies

Molecular alterations in cancerous tissues exhibit intercellular genetic and epigenetic heterogeneity, complicating the performance of diagnostic assays, particularly for early cancer detection. Conventional liquid biopsy methods have limited sensitivity and/or ability to assess epigenetic heterogeneity of rare epiallelic variants cost-effectively. We report an approach, named REM-DREAMing (Ratiometric-Encoded Multiplex Discrimination of Rare EpiAlleles by Melt), which leverages a digital microfluidic platform that incorporates a ratiometric fluorescence multiplex detection scheme and precise digital high-resolution melt analysis to enable low-cost, parallelized analysis of heterogeneous methylation patterns on a molecule-by-molecule basis for the detection of cancer in liquid biopsies. We applied the platform to simultaneously assess intermolecular epigenetic heterogeneity in five methylation biomarkers for improved, blood-based screening for early-stage non-small cell lung cancer. In a cohort of 48 low-volume liquid biopsy specimens from patients with indeterminant pulmonary nodules, we show that assessment of intermolecular methylation density distributions can notably improve the performance of multigene methylation biomarker panels for the early detection of cancer.

Does Every Strain of Pseudomonas aeruginosa Attack the Same? Results of a Study of the Prevalence of Virulence Factors of Strains Obtained from Different Animal Species in Northeastern Poland

BACKGROUND

Pseudomonas aeruginosa is a pathogen that causes infections in animals and humans, with veterinary implications including ear infections in dogs, respiratory diseases in cats, and mastitis in ruminants. In humans, it causes severe hospital-acquired infections, particularly in immunosuppressed patients. This study aimed to identify and assess the prevalence of specific virulence factors in Pseudomonas aeruginosa isolates.

METHODS

We analyzed 98 Pseudomonas aeruginosa isolates from various animal samples (dogs, cats, ruminants, fowl) from northeastern Poland in 2019-2022 for virulence-related genes (toxA, exoU, exoT, exoS, lasB, plcN, plcH, pldA, aprA, gacA, algD, pelA, endA, and oprF) by PCR and assessed biofilm formation at 48 and 72 h. Genomic diversity was assessed by ERIC-PCR.

RESULTS

The obtained results showed that all strains harbored the pelA gene (100%), while the lowest prevalence was found for pldA (24%) and exoU (36%). Regardless of the animal species, strong biofilm forming ability was prevalent among the strains after both 48 h (75%) and 72 h (74%). We obtained as many as 87 different genotyping profiles, where the dominant one was profile ERIC-48, observed in four strains.

CONCLUSIONS

No correlation was found between presence or absence of determined genes and the nature of infection. Similarly, no correlation was found between biofilm-forming genes and biofilm strength. The high genetic diversity indicates challenges for effective prevention, emphasizing the need for ongoing monitoring and research.

A novel technique for estimating age and demography of long-lived seabirds (genus Pterodroma) using an epigenetic clock for Gould's petrel (Pterodroma leucoptera)

Understanding the demography of wildlife populations is a key component for ecological research, and where necessary, supporting the conservation and management of long-lived animals. However, many animals lack phenological changes with which to determine individual age; therefore, gathering this fundamental information presents difficulties. More so for species that are rare, highly mobile, migratory and those that reside in inaccessible habitats. Until recently, the primary method to measure demography is through labour intensive mark-recapture approaches, necessitating decades of effort for long-lived species. Gadfly petrels (genus: Pterodroma) are one such taxa that are overrepresented with threatened and declining species, and for which numerous aspects of their ecology present challenges for research, monitoring and recovery efforts. To overcome some of these challenges, we developed the first DNA methylation (DNAm) demography technique to estimate the age of petrels, using the epigenetic clock of Gould's petrels (Pterodroma leucoptera). We collected reference blood samples from known-aged Gould's petrels at a long-term monitored population on Cabbage Tree Island, Australia. Epigenetic ages were successfully estimated for 121 individuals ranging in age from zero (fledgling) to 30 years of age, showing a mean error of 2.24 ± 0.17 years between the estimated and real age across the population. This is the first development of an epigenetic clock using multiplex PCR sequencing in a bird. This method enables demography to be measured with relative accuracy in a single sampling trip. This technique can provide information for emerging demographic risks that can mask declines in long-lived seabird populations and be applied to other Pterodroma populations.

Single-molecule epiallelic profiling of DNA derived from routinely collected Pap specimens for noninvasive detection of ovarian cancer

Recent advances in molecular analyses of ovarian cancer have revealed a wealth of promising tumour-specific biomarkers, including protein, DNA mutations and methylation; however, reliably detecting such alterations at satisfactorily high sensitivity and specificity through low-cost methods remains challenging, especially in early-stage diseases. Here we present PapDREAM, a new approach that enables detection of rare, ovarian-cancer-specific aberrations of DNA methylation from routinely-collected cervical Pap specimens. The PapDREAM approach employs a microfluidic platform that performs highly parallelized digital high-resolution melt to analyze locus-specific DNA methylation patterns on a molecule-by-molecule basis at or near single CpG-site resolution at a fraction (< 1/10th) of the cost of next-generation sequencing techniques. We demonstrate the feasibility of the platform by assessing intermolecular heterogeneity of DNA methylation in a panel of methylation biomarker loci using DNA derived from Pap specimens obtained from a cohort of 43 women, including 18 cases with ovarian cancer and 25 cancer-free controls. PapDREAM leverages systematic multidimensional bioinformatic analyses of locus-specific methylation heterogeneity to improve upon Pap-specimen-based detection of ovarian cancer, demonstrating a clinical sensitivity of 50% at 99% specificity in detecting ovarian cancer cases with an area under the receiver operator curve of 0.90. We then establish a logistic regression model that could be used to identify high-risk patients for subsequent clinical follow-up and monitoring. The results of this study support the utility of PapDREAM as a simple, low-cost screening method with the potential to integrate with existing clinical workflows for early detection of ovarian cancer. KEY POINTS: We present a microfluidic platform for detection and analysis of rare, heterogeneously methylated DNA within Pap specimens towards detection of ovarian cancer. The platform achieves high sensitivity (fractions <0.00005%) at a suitably low cost (∼$25) for routine screening applications. Furthermore, it provides molecule-by-molecule quantitative analysis to facilitate further study on the effect of heterogeneous methylation on cancer development.

guidedNOMe-seq quantifies chromatin states at single allele resolution for hundreds of custom regions in parallel

Since the introduction of next generation sequencing technologies, the field of epigenomics has evolved rapidly. However, most commonly used assays are enrichment-based methods and thus only semi-quantitative. Nucleosome occupancy and methylome sequencing (NOMe-seq) allows for quantitative inference of chromatin states with single locus resolution, but this requires high sequencing depth and is therefore prohibitively expensive to routinely apply to organisms with large genomes. To overcome this limitation, we introduce guidedNOMe-seq, where we combine NOMe profiling with large scale sgRNA synthesis and Cas9-mediated region-of-interest (ROI) liberation. To facilitate quantitative comparisons between multiple samples, we additionally develop an R package to standardize differential analysis of any type of NOMe-seq data. We extensively benchmark guidedNOMe-seq in a proof-of-concept study, dissecting the interplay of ChAHP and CTCF on chromatin. In summary we present a cost-effective, scalable, and customizable target enrichment extension to the existing NOMe-seq protocol allowing genome-scale quantification of nucleosome occupancy and transcription factor binding at single allele resolution.

Multiplexed electrochemical liposomes applied to the detection of nucleic acids for Influenza A, Influenza B and SARS-CoV-2

Multiplexing is a relevant strategy for biosensors to improve accuracy and decision-making due to the increased amount of simultaneously obtained information. Liposomes offer unique benefits for label-based multiplexing since a variety of different marker molecules can be encapsulated, leading to intrinsic signal amplification and enabling a variety of detection formats. We successfully developed an electrochemical (EC) liposome-based platform technology for the simultaneous detection of at least three analytes by studying parameters to ensure specific and sensitive bioassay performance. Influenza A and B and SARS-CoV-2 sequences served as model system in a standard sandwich hybridization assay. Studies included encapsulants, probe distribution on liposomes and capture beads, assay setup and interferences between liposomes to also ensure a generalization of the platform. Ruthenium hexamine(III), potassium hexacyanoferrate(II) and m-carboxy luminol, when encapsulated separately into a liposome, provided desirable long-term stability of at least 12 months and no cross-signals between liposomes. Through the optimization process, low limits of detections of 1.6 nmol L-1, 125 pmol L-1 and 130 pmol L-1, respectively, were achieved in a multiplexed assay setup, which were similar to singleplex assays. Non-specific interactions were limited to 25.1%, 7.6% and 7.5%, respectively, through sequential liposome incubations and singleplex capture bead designs. Here, ruthenium hexamine liposomes had only mediocre performance so that low overall signal strength translated into higher LODs and worse specificity. A different marker such as ferroin may be an option in the future. The identification of further electrochemical markers will provide new opportunities for liposomes to function as multiplex, orthogonal or internal standard labels in electrochemical bioassays.

A Rapid and Inexpensive PCR Test for Mastitis Diagnosis Based on NGS Data

Mastitis is a common mammary gland disease of dairy cattle caused by a wide range of organisms including bacteria, fungi and algae. Mastitis contributes to economic losses of dairy farms due to reduced yield and poor quality of milk. Since the correct identification of pathogens responsible for the development of mastitis is crucial to the success of treatment, it is necessary to develop a quick and accurate test to distinguish the main pathogens causing this disease. In this paper, we describe the development of a test based on the multiplex polymerase chain reaction (PCR) method allowing for the identification of Streptococcus agalactiae, Streptococcus dysgalactiae, Streptococcus uberis and Staphylococcus aureus. When creating our test, we relied on the results from new generation sequencing (NGS) for accurate determination of species affiliation. The multiplex PCR test was verified on 100 strains including veterinary samples, ATCC and Polish Collection of Microorganisms (PCM) reference strains. The obtained results indicate that this test is accurate and displays high specificity. It may serve as a valuable molecular tool for the detection of major mastitis pathogens.

MethMarkerDB: a comprehensive cancer DNA methylation biomarker database

DNA methylation plays a crucial role in tumorigenesis and tumor progression, sparking substantial interest in the clinical applications of cancer DNA methylation biomarkers. Cancer-related whole-genome bisulfite sequencing (WGBS) data offers a promising approach to precisely identify these biomarkers with differentially methylated regions (DMRs). However, currently there is no dedicated resource for cancer DNA methylation biomarkers with WGBS data. Here, we developed a comprehensive cancer DNA methylation biomarker database (MethMarkerDB, https://methmarkerdb.hzau.edu.cn/), which integrated 658 WGBS datasets, incorporating 724 curated DNA methylation biomarker genes from 1425 PubMed published articles. Based on WGBS data, we documented 5.4 million DMRs from 13 common types of cancer as candidate DNA methylation biomarkers. We provided search and annotation functions for these DMRs with different resources, such as enhancers and SNPs, and developed diagnostic and prognostic models for further biomarker evaluation. With the database, we not only identified known DNA methylation biomarkers, but also identified 781 hypermethylated and 5245 hypomethylated pan-cancer DMRs, corresponding to 693 and 2172 genes, respectively. These novel potential pan-cancer DNA methylation biomarkers hold significant clinical translational value. We hope that MethMarkerDB will help identify novel cancer DNA methylation biomarkers and propel the clinical application of these biomarkers.

Amplicon-Based Bisulfite Conversion-NGS DNA Methylation Analysis Protocol

DNA methylation is an important epigenetic modification that regulates chromatin structure and the cell-type-specific expression of genes. The association of aberrant DNA methylation with many diseases, as well as the increasing interest in modifying the methylation mark in a directed manner at genomic sites using epigenome editing for research and therapeutic purposes, increases the need for easy and efficient DNA methylation analysis methods. The standard approach to analyze DNA methylation with a single-cytosine resolution is bisulfite conversion of DNA followed by next-generation sequencing (NGS). In this chapter, we describe a robust, powerful, and cost-efficient protocol for the amplification of target regions from bisulfite-converted DNA, followed by a second PCR step to generate libraries for Illumina NGS. In the two consecutive PCR steps, first, barcodes are added to individual amplicons, and in the second PCR, indices and Illumina adapters are added to the samples. Finally, we describe a detailed bioinformatics approach to extract DNA methylation levels of the target regions from the sequencing data. Combining barcodes with indices enables a high level of multiplexing allowing to sequence multiple pooled samples in the same sequencing run. Therefore, this method is a robust, accurate, quantitative, and cheap approach for the readout of DNA methylation patterns at defined genomic regions.

Development of a Multiplex PCR Assay for Efficient Detection of Two Potential Probiotic Strains Using Whole Genome-Based Primers

Probiotics are microorganisms that exert strain-specific health-promoting effects on the host. Τhey are employed in the production of functional dairy or non-dairy food products; still, their detection in these complex matrices is a challenging task. Several culture-dependent and culture-independent methods have been developed in this direction; however, they present low discrimination at the strain level. Here, we developed a multiplex PCR assay for the detection of two potential probiotic lactic acid bacteria (LAB) strains, Lactiplantibacillus plantarum L125 and Lp. pentosus L33, in monocultures and yogurt samples. Unique genomic regions were identified via comparative genomic analysis and were used to produce strain-specific primers. Then, primer sets were selected that produced distinct electrophoretic DNA banding patterns in multiplex PCR for each target strain. This method was further implemented for the detection of the two strains in yogurt samples, highlighting its biotechnological applicability. Moreover, it can be applied with appropriate modifications to detect any bacterial strain with available WGS.

Contribution of Spontaneous Mutations to Quantitative and Molecular Variation at the Highly Repetitive rDNA Locus in Yeast

The ribosomal DNA array in Saccharomyces cerevisiae consists of many tandem repeats whose copy number is believed to be functionally important but highly labile. Regulatory mechanisms have evolved to maintain copy number by directed mutation, but how spontaneous variation at this locus is generated and selected has not been well characterized. We applied a mutation accumulation approach to quantify the impacts of mutation and selection on this unique genomic feature across hundreds of mutant strains. We find that mutational variance for this trait is relatively high, and that unselected mutations elsewhere in the genome can disrupt copy number maintenance. In consequence, copy number generally declines gradually, consistent with a previously proposed model of rDNA maintenance where a downward mutational bias is normally compensated by mechanisms that increase copy number when it is low. This pattern holds across ploidy levels and strains in the standard lab environment but differs under some stressful conditions. We identify several alleles, gene categories, and genomic features that likely affect copy number, including aneuploidy for chromosome XII. Copy number change is associated with reduced growth in diploids, consistent with stabilizing selection. Levels of standing variation in copy number are well predicted by a balance between mutation and stabilizing selection, suggesting this trait is not subject to strong diversifying selection in the wild. The rate and spectrum of point mutations within the rDNA locus itself are distinct from the rest of the genome and predictive of polymorphism locations. Our findings help differentiate the roles of mutation and selection and indicate that spontaneous mutation patterns shape several aspects of ribosomal DNA evolution.

Comprehensive methylome sequencing reveals prognostic epigenetic biomarkers for prostate cancer mortality

BACKGROUND

Prostate cancer is a clinically heterogeneous disease with a subset of patients rapidly progressing to lethal-metastatic prostate cancer. Current clinicopathological measures are imperfect predictors of disease progression. Epigenetic changes are amongst the earliest molecular changes in tumourigenesis. To find new prognostic biomarkers to enable earlier intervention and improved outcomes, we performed methylome sequencing of DNA from patients with localised prostate cancer and long-term clinical follow-up.

METHODS

We used whole-genome bisulphite sequencing (WGBS) to comprehensively map and compare DNA methylation of radical prostatectomy tissue between patients with lethal disease (n = 7) and non-lethal (n = 8) disease (median follow-up 19.5 years). Validation of differentially methylated regions (DMRs) was performed in an independent cohort (n = 185, median follow-up 15 years) using targeted multiplex bisulphite sequencing of candidate regions. Survival was assessed via univariable and multivariable analyses including clinicopathological measures (log-rank and Cox regression models).

RESULTS

WGBS data analysis identified cancer-specific methylation patterns including CpG island hypermethylation, and hypomethylation of repetitive elements, with increasing disease risk. We identified 1420 DMRs associated with prostate cancer-specific mortality (PCSM), which showed enrichment for gene sets downregulated in prostate cancer and de novo methylated in cancer. Through comparison with public prostate cancer datasets, we refined the DMRs to develop an 18-gene prognostic panel. Applying this panel to an independent cohort, we found significant associations between PCSM and hypermethylation at EPHB3, PARP6, TBX1, MARCH6 and a regulatory element within CACNA2D4. Strikingly in a multivariable model, inclusion of CACNA2D4 methylation was a better predictor of PCSM versus grade alone (Harrell's C-index: 0.779 vs. 0.684).

CONCLUSIONS

Our study provides detailed methylome maps of non-lethal and lethal prostate cancer and identifies novel genic regions that distinguish these patient groups. Inclusion of our DNA methylation biomarkers with existing clinicopathological measures improves prognostic models of prostate cancer mortality, and holds promise for clinical application.

Modelling clinical DNA fragmentation in the development of universal PCR-based assays for bisulfite-converted, formalin-fixed and cell-free DNA sample analysis

In fragmented DNA, PCR-based methods quantify the number of intact regions at a specific amplicon length. However, the relationship between the population of DNA fragments within a sample and the likelihood they will amplify has not been fully described. To address this, we have derived a mathematical equation that relates the distribution profile of a stochastically fragmented DNA sample to the probability that a DNA fragment within that sample can be amplified by any PCR assay of arbitrary length. Two panels of multiplex PCR assays for quantifying fragmented DNA were then developed: a four-plex panel that can be applied to any human DNA sample and used to estimate the percentage of regions that are intact at any length; and a two-plex panel optimized for quantifying circulating cell-free DNA (cfDNA). For these assays, regions of the human genome least affected by copy number aberration were identified and selected; within these copy-neutral regions, each PCR assay was designed to amplify both genomic and bisulfite-converted DNA; and all assays were validated for use in both conventional qPCR and droplet-digital PCR. Finally, using the cfDNA-optimized assays we find evidence of universally conserved nucleosome positioning among individuals.

MS-HRM protocol: a simple and low-cost approach for technical validation of next-generation methylation sequencing data

DNA methylation is a fundamental epigenetic process and have a critical role in many biological processes. The study of DNA methylation at a large scale of genomic levels is widely conducted by several techniques that are next-generation sequencing (NGS)-based methods. Methylome data revealed by DNA methylation next-generation sequencing (mNGS), should be always verified by another technique which they usually have a high cost. In this study, we offered a low-cost approach to corroborate the mNGS data. In this regard, mNGS was performed on 6 colorectal cancer (case group) and 6 healthy individual colon tissue (control group) samples. An R-script detected differentially methylated regions (DMRs), was further validated by high resolution melting (MS-HRM) analysis. After analyzing the data, the algorithm found 194 DMRs. Two locations with the highest level of methylation difference were verified by MS-HRM, which their results were in accordance with the mNGS. Therefore, in the present study, we suggested MS-HRM as a simple, accurate and low-cost method, useful for confirming methylation sequencing results.

Technical considerations in PCR-based assay design for diagnostic DNA methylation cancer biomarkers

Background

DNA methylation biomarkers for early detection, risk stratification and treatment response in cancer have been of great interest over the past decades. Nevertheless, clinical implementation of these biomarkers is limited, as only < 1% of the identified biomarkers is translated into a clinical or commercial setting. Technical factors such as a suboptimal genomic location of the assay and inefficient primer or probe design have been emphasized as important pitfalls in biomarker research. Here, we use eleven diagnostic DNA methylation biomarkers for colorectal cancer (ALX4, APC, CDKN2A, MGMT, MLH1, NDRG4, SDC2, SFRP1, SFRP2, TFPI1 and VIM), previously described in a systematic literature search, to evaluate these pitfalls.

Results

To assess the genomic assay location, the optimal genomic locations according to TCGA data were extracted and compared to the genomic locations used in the published assays for all eleven biomarkers. In addition, all primers and probes were technically evaluated according to several criteria, based on literature and expert opinion. Both assay location and assay design quality varied widely among studies.

Conclusions

Large variation in both assay location and design hinders the development of future DNA methylation biomarkers as well as inter-study comparability.

Novel method utilizing bisulfite conversion with dual amplification-refractory mutation system polymerase chain reaction to detect circulating pancreatic β-cell cfDNA

Aims/Introduction

Several research groups have reported methods for quantifying pancreatic beta cell (β‐cell) injury by measuring β‐cell‐specific CpG unmethylation of the insulin gene in circulation using digital droplet PCR or next‐generation sequencing. However, these methods have certain disadvantages, such as the need to consider the background signal owing to the small number of target CpG sites and the need for unique equipment.

Materials and Methods

We established a novel method for detecting four CpG unmethylations of the insulin gene using two‐step amplification refractory mutation system PCR. We applied it to type 1 diabetes (T1D) patients with a wide range of disease durations and to healthy adults.

Results

The assay showed high linearity and could detect a single copy of unmethylated insulin DNA in experiments using methylated and unmethylated plasmid DNA. The unmethylated insulin DNA level in the type 1 diabetes group, whose β‐cell mass was considerably reduced, was similar to that of healthy adults. An inverse correlation was observed between copy number and disease duration in patients with unmethylated insulin DNA‐positive type 1 diabetes.

Conclusions

We developed a novel method for detecting unmethylated insulin DNA in circulation that can be performed using a conventional real‐time PCR system. This method would be useful for analyzing dynamic profiles of β‐cells in human disease such as type 1 diabetes.

End-point RT-PCR based on a conservation landscape for SARS-COV-2 detection

End-point RT-PCR is a suitable alternative diagnostic technique since it is cheaper than RT-qPCR tests and can be implemented on a massive scale in low- and middle-income countries. In this work, a bioinformatic approach to guide the design of PCR primers was developed, and an alternative diagnostic test based on end-point PCR was designed. End-point PCR primers were designed through conservation analysis based on kmer frequency in SARS-CoV-2 and human respiratory pathogen genomes. Highly conserved regions were identified for primer design, and the resulting PCR primers were used to amplify 871 nasopharyngeal human samples with a previous RT-qPCR based SARS-CoV-2 diagnosis. The diagnostic test showed high accuracy in identifying SARS-CoV-2-positive samples including B.1.1.7, P.1, B.1.427/B.1.429 and B.1.617.2/ AY samples with a detection limit of 7.2 viral copies/µL. In addition, this test could discern SARS-CoV-2 infection from other viral infections with COVID-19-like symptomatology. The designed end-point PCR diagnostic test to detect SARS-CoV-2 is a suitable alternative to RT-qPCR. Since the proposed bioinformatic approach can be easily applied in thousands of viral genomes and over highly divergent strains, it can be used as a PCR design tool as new SARS-CoV-2 variants emerge. Therefore, this end-point PCR test could be employed in epidemiological surveillance to detect new SARS-CoV-2 variants as they emerge and propagate.

Epigenome - A mediator for host-microbiome crosstalk

The interaction between the gut and its eventual trillions of microbe inhabitants during microbial colonization, represents a critical time period for establishing the overall health and wellbeing of an individual. The gut microbiome represents a diverse community of microbes that are critical for many physiological roles of the host including host metabolism. These processes are controlled by a fine-tuned chemical cross talk between the host and microbiota. Although the exact mechanisms behind this cross talk remains elusive, microbiota induced epigenetic mechanisms like DNA methylation and histone modifications may be key. This review presents our perspective on the epigenome as a mediator for host-microbiota cross talk, as well as methodology to study epigenetics, the role of dysbiosis in disease, and how the gut microbiome-host axis may be used in personal medicine.

Age-related DNA methylation analysis for forensic age estimation using post-mortem blood samples from Japanese individuals

As one of external visible characteristics (EVCs) in forensic phenotyping, age estimation is essential to providing additional information about a sample donor. With the development of epigenetics, age-related DNA methylation may be used as a reliable predictor of age estimation. With the aim of building a feasible age estimation model for Japanese individuals, 53 CpG sites distributed between 11 candidate genes were selected from previous studies. The DNA methylation level of each target CpG site was identified and measured on a massive parallel platform (synthesis by sequencing, Illumina, California, United States) from 60 forensic blood samples during the initial training phase. Multiple linear regression and quantile regression analyses were later performed to build linear and quantile age estimation models, respectively. Four CpG sites on four genes- ASPA, ELOVL2, ITGA2B, and PDE4C -, were found to be highly correlated with chronological age in DNA samples from Japanese individuals (|R| > 0.75). Subsequently, an independent validation dataset (n = 30) was used to verify and evaluate the performance of the two models. Comparison of mean absolute deviation (MAD) with other indicators showed that both models provide accurate age predictions (MAD: linear = 6.493 years; quantile = 6.243 years). The quantile model, however, can provide the changeable prediction intervals that grow wider with increasing age, and this tendency is consistent with the natural aging process in humans. Hence, the quantile model is recommended in this study.

Application of fragment analysis based on methylation status mobility difference to identify vaginal secretions

Identifying vaginal secretions attaching or adhering to a suspect's belongings would be beneficial for reconstructing the events that have taken place during a sexual assault. The present study describes a novel approach to identify vaginal secretions by fragment analysis using capillary electrophoresis, based on the mobility differences of PCR amplicons from bisulfite-treated DNA depending on methylation status. We targeted three genome regions including each of three vaginal secretion-specific methylated CpG sites reported previously: cg25416153, cg09765089, and cg14991487. In all three genome regions, the amplicon peaks for methylated genomic DNA (gDNA) sequences were only detected in vaginal samples, whereas samples of other body fluids (blood, saliva, semen, and deposit on skin surface) only showed amplicon peaks for unmethylated gDNA sequences. In vaginal secretions, the methylation ratio of each of the three targeted regions between samples was variable, while the ratios at the three regions in each sample were similar. Furthermore, commercial vaginal epithelial cells were completely methylated at the three regions. Therefore, vaginal secretion-specific methylation may derive from vaginal epithelial cells present in the sample. In forensic cases with a limited amount of DNA, the reproducibility of a detected peak using the present method is not high due to degradation of DNA by bisulfite treatment and subsequent stochastic PCR bias. However, it was possible to detect peaks from methylated DNA sequences by performing PCR and capillary electrophoresis in triplicate after bisulfite treatment, even when bisulfite treatment was performed using 0.5 ng of gDNA from vaginal secretions. In addition, the level of methylation at each targeted region was found to be stable in vaginal secretions stored for 1 year at room temperature. Therefore, we conclude that detection of the visual peak from vaginal secretion-specific methylated DNA sequence is useful to prove the presence of vaginal secretions. This approach has the potential to analyze multiple marker regions simultaneously, and may provide a new multiplex assay to identify various body fluids.

Hairpin-Bisulfite PCR

Ligation of a hairpin oligonucleotide to genomic DNA prior to bisulfite conversion and PCR amplification physically links the two complementary DNA strands. This additional step in the conversion procedure overcomes the limitations of conventional bisulfite sequencing where information of the cytosine methylation status is only obtained from one of the two strands of an individual DNA molecule. Sequences derived from hairpin bisulfite PCR products reveal the dynamics of this epigenetic memory system on both strands of individual DNA molecules. The chapter describes a reliable step-by-step procedure to generate hairpin-linked DNA. It also provides a guide for efficient bisulfite conversion that is suitable for both conventional and hairpin bisulfite sequencing approaches.

Development of the VISAGE enhanced tool and statistical models for epigenetic age estimation in blood, buccal cells and bones

DNA methylation is known as a biomarker for age with applications in forensics. Here we describe the VISAGE (VISible Attributes through GEnomics) Consortium's enhanced tool for epigenetic age estimation in somatic tissues. The tool is based on eight DNA methylation markers (44 CpGs), bisulfite multiplex PCR followed by sequencing on the MiSeq FGx platform, and three statistical prediction models for blood, buccal cells and bones. The model for blood is based on six CpGs from ELOVL2, MIR29B2CHG, KLF14, FHL2, TRIM59 and PDE4C, and predicts age with a mean absolute error (MAE) of 3.2 years, while the model for buccal cells includes five CpGs from PDE4C, MIR29B2CHG, ELOVL2, KLF14 and EDARADD and predicts age with MAE of 3.7 years, and the model for bones has six CpGs from ELOVL2, KLF14, PDE4C and ASPA and predicts age with MAE of 3.4 years. The VISAGE enhanced tool for age estimation in somatic tissues enables reliable collection of DNA methylation data from small amounts of DNA using a sensitive multiplex MPS assay that provides accurate estimation of age in blood, buccal swabs, and bones using the statistical model tailored to each tissue.

MethPanel: a parallel pipeline and interactive analysis tool for multiplex bisulphite PCR sequencing to assess DNA methylation biomarker panels for disease detection

Summary

DNA methylation patterns in a cell are associated with gene expression and the phenotype of a cell, including disease states. Bisulphite PCR sequencing is commonly used to assess the methylation profile of genomic regions between different cells. Here we have developed MethPanel, a computational pipeline with an interactive graphical interface to rapidly analyse multiplex bisulphite PCR sequencing data. MethPanel comprises a complete analysis workflow from genomic alignment to DNA methylation calling and supports an unlimited number of PCR amplicons and input samples. MethPanel offers important and unique features, such as calculation of an epipolymorphism score and bisulphite PCR bias correction capabilities, and is designed so that the methylation data from all samples can be processed in parallel. The outputs are automatically forwarded to a shinyApp for convenient display, visualization and remotely sharing data with collaborators and clinicians.

Availabilityand implementation

MethPanel is freely available at https://github.com/thinhong/MethPanel.

Supplementary information

Supplementary data are available at Bioinformatics online.

A DNA methylation age predictor for zebrafish

Changes in DNA methylation at specific CpG sites have been used to build predictive models to estimate animal age, predominantly in mammals. Little testing for this effect has been conducted in other vertebrate groups, such as bony fish, the largest vertebrate class. The development of most age-predictive models has relied on a genome-wide sequencing method to obtain a DNA methylation level, which makes it costly to deploy as an assay to estimate age in many samples. Here, we have generated a reduced representation bisulfite sequencing data set of caudal fin tissue from a model fish species, zebrafish (Danio rerio), aged from 11.9-60.1 weeks. We identified changes in methylation at specific CpG sites that correlated strongly with increasing age. Using an optimised unique set of 26 CpG sites we developed a multiplex PCR assay that predicts age with an average median absolute error rate of 3.2 weeks in zebrafish between 10.9-78.1 weeks of age. We also demonstrate the use of multiplex PCR as an efficient quantitative approach to measure DNA methylation for the use of age estimation. This study highlights the potential further use of DNA methylation as an age estimation method in non-mammalian vertebrate species.

AmpliconDesign - an interactive web server for the design of high-throughput targeted DNA methylation assays

Targeted analysis of DNA methylation patterns based on bisulfite-treated genomic DNA (BT-DNA) is considered as a gold-standard for epigenetic biomarker development. Existing software tools facilitate primer design, primer quality control or visualization of primer localization. However, high-throughput design of primers for BT-DNA amplification is hampered by limits in throughput and functionality of existing tools, requiring users to repeatedly perform specific tasks manually. Consequently, the design of PCR primers for BT-DNA remains a tedious and time-consuming process. To bridge this gap, we developed AmpliconDesign, a webserver providing a scalable and user-friendly platform for the design and analysis of targeted DNA methylation studies based on BT-DNA, e.g. deep amplicon bisulfite sequencing (ampBS-seq) or EpiTYPER MassArray. Core functionality of the web server includes high-throughput primer design and binding site validation based on in silico bisulfite-converted DNA sequences, prediction of fragmentation patterns for EpiTYPER MassArray, an interactive quality control as well as a streamlined analysis workflow for ampBS-seq.

SNP barcodes provide higher resolution than microsatellite markers to measure Plasmodium vivax population genetics

BACKGROUND

Genomic surveillance of malaria parasite populations has the potential to inform control strategies and to monitor the impact of interventions. Barcodes comprising large numbers of single nucleotide polymorphism (SNP) markers are accurate and efficient genotyping tools, however may need to be tailored to specific malaria transmission settings, since 'universal' barcodes can lack resolution at the local scale. A SNP barcode was developed that captures the diversity and structure of Plasmodium vivax populations of Papua New Guinea (PNG) for research and surveillance.

METHODS

Using 20 high-quality P. vivax genome sequences from PNG, a total of 178 evenly spaced neutral SNPs were selected for development of an amplicon sequencing assay combining a series of multiplex PCRs and sequencing on the Illumina MiSeq platform. For initial testing, 20 SNPs were amplified in a small number of mono- and polyclonal P. vivax infections. The full barcode was then validated by genotyping and population genetic analyses of 94 P. vivax isolates collected between 2012 and 2014 from four distinct catchment areas on the highly endemic north coast of PNG. Diversity and population structure determined from the SNP barcode data was then benchmarked against that of ten microsatellite markers used in previous population genetics studies.

RESULTS

From a total of 28,934,460 reads generated from the MiSeq Illumina run, 87% mapped to the PvSalI reference genome with deep coverage (median = 563, range 56-7586) per locus across genotyped samples. Of 178 SNPs assayed, 146 produced high-quality genotypes (minimum coverage = 56X) in more than 85% of P. vivax isolates. No amplification bias was introduced due to either polyclonal infection or whole genome amplification (WGA) of samples before genotyping. Compared to the microsatellite panels, the SNP barcode revealed greater variability in genetic diversity between populations and geographical population structure. The SNP barcode also enabled assignment of genotypes according to their geographic origins with a significant association between genetic distance and geographic distance at the sub-provincial level.

CONCLUSIONS

High-throughput SNP barcoding can be used to map variation of malaria transmission dynamics at sub-national resolution. The low cost per sample and genotyping strategy makes the transfer of this technology to field settings highly feasible.

Comprehensive evaluation of targeted multiplex bisulphite PCR sequencing for validation of DNA methylation biomarker panels

BACKGROUND

DNA methylation is a well-studied epigenetic mark that is frequently altered in diseases such as cancer, where specific changes are known to reflect the type and severity of the disease. Therefore, there is a growing interest in assessing the clinical utility of DNA methylation as a biomarker for diagnosing disease and guiding treatment. The development of an accurate loci-specific methylation assay, suitable for use on low-input clinical material, is crucial for advancing DNA methylation biomarkers into a clinical setting. A targeted multiplex bisulphite PCR sequencing approach meets these needs by allowing multiple DNA methylated regions to be interrogated simultaneously in one experiment on limited clinical material.

RESULTS

Here, we provide an updated protocol and recommendations for multiplex bisulphite PCR sequencing (MBPS) assays for target DNA methylation analysis. We describe additional steps to improve performance and reliability: (1) pre-sequencing PCR optimisation which includes assessing the optimal PCR cycling temperature and primer concentration and (2) post-sequencing PCR optimisation to achieve uniform coverage of each amplicon. We use a gradient of methylated controls to demonstrate how PCR bias can be assessed and corrected. Methylated controls also allow assessment of the sensitivity of methylation detection for each amplicon. Here, we show that the MBPS assay can amplify as little as 0.625 ng starting DNA and can detect methylation differences of 1% with a sequencing coverage of 1000 reads. Furthermore, the multiplex bisulphite PCR assay can comprehensively interrogate multiple regions on 1-5 ng of formalin-fixed paraffin-embedded DNA or circulating cell-free DNA.

CONCLUSIONS

The MBPS assay is a valuable approach for assessing methylated DNA regions in clinical samples with limited material. The optimisation and additional quality control steps described here improve the performance and reliability of this method, advancing it towards potential clinical applications in biomarker studies.

Integrative Analysis of Methylation and Transcriptional Profiles to Reveal the Genetic Stability of Cashmere Traits in the Tβ4 Overexpression of Cashmere Goats

Simple Summary

Cashmere goats have double coats consisting of non-medullated fine inner hairs or cashmere fibers produced by secondary hair follicles (SHFs) and guard hairs produced by primary hair follicles (PHFs). Cashmere is an important economic product worldwide and the world market for cashmere is increasing while the current production of cashmere is limited. Thymosin β4 (Tβ4), a 4.9-kDa protein, contains 43 amino acids. Here, we produced Tβ4 overexpression (Tβ4-OE) offspring using two methods. The somatic cell nuclear transfer (SCNT) goats had increased hair follicle development and higher cashmere yields than wild type (WT) and natural mating (NM) goats. Taken together, our results showed that DNA methylation affected the expression of differentially expressed genes (DEGs) between generations and the genetic stability of cashmere traits.

Abstract

DNA methylation alteration is frequently observed in exogenous gene silencing and may play important roles in the genetic stability of traits. Cashmere is derived from the secondary hair follicles (SHFs) of cashmere goats, which are morphogenetically distinct from primary hair follicles (PHFs). Here, in light of having initially produced 15 Tβ4 overexpression (Tβ4-OE) cashmere goats which had more SHFs than the wild type (WT) goats, and produced more cashmere, we produced Tβ4-OE offsprings both via somatic cell nuclear transfer (SCNT) and via natural mating (NM). However, the desired trait exhibited lower fixation in the line-bred offspring compared to the SCNT offspring. Integrative analysis of methylation and transcriptional profiles showed that this might be due to the influence of methylation on the expression of differentially expressed genes (DEGs) between generations, which was mutually consistent with the results of the functional and pathway enrichment analysis of differentially methylated regions (DMRs) and DEGs. Overall, our study systematically describes the DNA methylation characteristics between generations of cashmere goats and provides a basis for improving genetic stability.

Endometrial genome-wide DNA methylation patterns of Guanzhong dairy goats at days 5 and 15 of the gestation period

Uterine receptivity for the embryo is established and maintained through a series of precise cellular and molecular events, such as DNA methylation. There have been no studies to elucidate entire genome DNA methylation changes associated with embryo receptivity development of the endometrium (RE). In the present study, there was development of a complete genome-wide DNA methylome maps of the RE using whole-genome bisulphite sequencing and bioinformatics analysis. As many as 163.06 Gb of sequencing data averaging 81.53 Gb per sample were obtained for genome bisulphite sequencing of endometrium samples. There were distinct genome-wide DNA methylation patterns in pre-receptive endometrium (PE; Day 5 of gestation) and RE (Day 15 of gestation). There were as many as 16,467 differentially methylated regions (DMRs); 21,391 DMRs were less methylated in RE samples compared with PE samples (P-values ≤ 0.05 and |log2 (fold change)| ≥ 2). Compared with PE samples, methylation ratios of IGF2BP2, ACOX2, PTGDS, VEGFB and PTGDR2 genes were markedly less in RE samples (P-value ≤ 0.05 and |log2 (fold change)| ≥ 2). Conversely, in RE samples there was a markedly greater methylation ratio of IGFBP3 and IGF1R genes. The results of KEGG analysis indicated that these genes were involved in the signalling pathways for insulin, mitogen-activated protein kinase, gonadotropin-releasing hormone, vascular endothelial growth factor and progesterone-mediated oocyte maturation, which participated in differential regulation of goat endometrial development during receptive and prereceptive phases. The results of previous and the present study indicate resulting proteins of IGF2BP2, PTGDS, VEGFB, PGR, IGFBP3 and IGF1R gene expression may have important functions in regulating endometrial receptivity for the embryo.

PrimerROC: accurate condition-independent dimer prediction using ROC analysis

To-date systematic testing and comparison of the accuracy of available primer-dimer prediction software has never been conducted, due in part to a lack of tools able to measure the efficacy of Gibbs free energy (ΔG) calculations at predicting dimer formation in PCR. To address this we have developed a novel online tool called PrimerROC (www.primer-dimer.com/roc/), which uses epidemiologically-based Receiver Operating Characteristic (ROC) curves to assess dimer prediction accuracy. Moreover, by integrating PrimerROC with our PrimerDimer prediction software we can determine a ΔG-based dimer-free threshold above which dimer formation is predicted unlikely to occur. Notably, PrimerROC determines this cut-off without any additional information such as salt concentration or annealing temperature, meaning that our PrimerROC method is an assay and condition independent prediction tool. To demonstrate the broad utility of PrimerROC we assessed the performance of seven publically available primer design and dimer analysis tools using a dataset of over 300 primer pairs. We found that our PrimerROC/PrimerDimer software consistently outperforms these other tools and can achieve predictive accuracies greater than 92%. To illustrate its predictive power this method was used in multiplex PCR design to successfully generate four resequencing assays containing up to 126 primers with no observable primer-primer amplification artefacts.

Genome-wide differences in DNA methylation changes in caprine ovaries between oestrous and dioestrous phases

Background

DNA methylation plays a vital role in reproduction. Entire genome DNA methylation changes during the oestrous phase (ES) and dioestrous phase (DS) in the ovaries of Guanzhong dairy goats were investigated using bisulphite sequencing to understand the molecular biological mechanisms of these goats’ oestrous cycle.

Results

We discovered distinct genome-wide DNA methylation patterns in ES and DS ovaries. A total of 26,910 differentially methylated regions were upregulated and 21,453 differentially methylated regions were downregulated in the ES samples compared with the DS samples (P-values ≤0.05 and fold change of methylation ratios ≥2). Differentially methylated region analysis showed hypomethylation in the gene body regions and hypermethylation in the joining region between upstream regions and gene bodies. The methylation ratios of the STAR, FGF2, FGF12, BMP5 and SMAD6 genes in the ES samples were lower than those of the DS samples (P-values ≤0.05 and fold change of methylation ratios ≥2). Conversely, the methylation ratios of the EGFR, TGFBR2, IGF2BP1 and MMD2 genes increased in the ES samples compared with the DS samples. In addition, 223 differentially methylated genes were found in the GnRH signalling pathway (KO04912), ovarian steroidogenesis pathway (KO04913), oestrogen signalling pathway (KO04915), oxytocin signalling pathway (KO04921), insulin secretion pathway (KO04911) and MAPK signalling pathway (KO04010).

Conclusions

This study is the first large-scale comparison of the high-resolution DNA methylation landscapes of oestrous and dioestrous ovaries from dairy goats. Previous studies and our investigations have shown that the NR5A2, STAR, FGF2 and BMP5 genes might have potential application value in regulating caprine oestrus.

Electronic supplementary material

The online version of this article (10.1186/s40104-018-0301-x) contains supplementary material, which is available to authorized users.

MethyMer: Design of combinations of specific primers for bisulfite sequencing of complete CpG islands

We present MethyMer, a Python-based tool aimed at selecting primers for amplification of complete CpG islands. These regions are difficult in terms of selecting appropriate primers because of their low-complexity, high GC content. Moreover, bisulfite treatment, in fact, leads to the reduction of the 4-letter alphabet (ATGC) to 3-letter one (ATG, except for methylated cytosines), and this also reduces region complexity and increases mispriming potential. MethyMer has a flexible scoring system, which optimizes the balance between various characteristics such as nucleotide composition, thermodynamic features (melting temperature, dimers [Formula: see text]G, etc.), the presence of CpG sites and polyN tracts, and primer specificity, which is assessed with aligning primers to the bisulfite-treated genome using bowtie (up to three mismatches are allowed). Users are able to customize desired or limit ranges of various parameters as well as penalties for non-desired values. Moreover, MethyMer allows picking up the optimal combination of PCR primer pairs to perform the amplification of a large genomic locus, e.g. CpG island or other hard-to-study region, with minimal overlap of the individual amplicons. MethyMer incorporates ENCODE genome annotation records (promoter/enhancer/insulator), The Cancer Genome Atlas (TCGA) CpG methylation data derived with Illumina Infinium 450K microarrays, and records on correlations between TCGA RNA-Seq and CpG methylation data for 20 cancer types. These databases are included in the MethyMer release. Our tool is available at https://sourceforge.net/projects/methymer/ .

Evaluation of Different Oligonucleotide Base Substitutions at CpG Binding sites in Multiplex Bisulfite-PCR sequencing

Multiplex bisulfite-PCR sequencing is a convenient and scalable method for the quantitative determination of the methylation state of target DNA regions. A challenge of this application is the presence of CpGs in the same region where primers are being placed. A common solution to the presence of CpGs within a primer-binding region is to substitute a base degeneracy at the cytosine position. However, the efficacy of different substitutions and the extent to which bias towards methylated or unmethylated templates may occur has never been evaluated in bisulfite multiplex sequencing applications. In response, we examined the performance of four different primer substitutions at the cytosine position of CpG's contained within the PCR primers. In this study, deoxyinosine-, 5-nitroindole-, mixed-base primers and primers with an abasic site were evaluated across a series of methylated controls. Primers that contained mixed- or deoxyinosine- base modifications performed most robustly. Mixed-base primers were further selected to determine the conditions that induce bias towards methylated templates. This identified an optimized set of conditions where the methylated state of bisulfite DNA templates can be accurately assessed using mixed-base primers, and expands the scope of bisulfite resequencing assays when working with challenging templates.

A multiplex microplatform for the detection of multiple DNA methylation events using gold–DNA affinity

We report a new multiplexed strategy for the electrochemical detection of regional DNA methylation across multiple regions.