Targeted DNA Methylation Analysis Facilitates Leukocyte Counts in Dried Blood Samples.

BACKGROUND: Cell-type specific DNA methylation (DNAm) can be employed to determine the numbers of leukocyte subsets in blood. In contrast to conventional methods for leukocyte counts, which are based on cellular morphology or surface marker protein expression, the cellular deconvolution based on DNAm levels is applicable for frozen or dried blood. Here, we further enhanced targeted DNAm assays for leukocyte counts in clinical application. METHODS: DNAm profiles of 40 different studies were compiled to identify CG dinucleotides (CpGs) with cell-type specific DNAm using a computational framework, CimpleG. DNAm levels at these CpGs were then measured with digital droplet PCR in venous blood from 160 healthy donors and 150 patients with various hematological disorders. Deconvolution was further validated with venous blood (n = 75) and capillary blood (n = 31) that was dried on Whatman paper or on Mitra microsampling devices. RESULTS: In venous blood, automated cell counting or flow cytometry correlated well with epigenetic estimates of relative leukocyte counts for granulocytes (r = 0.95), lymphocytes (r = 0.97), monocytes (r = 0.82), CD4 T cells (r = 0.84), CD8 T cells (r = 0.94), B cells (r = 0.96), and NK cells (r = 0.72). Similar correlations and precisions were achieved for dried blood samples. Spike-in with a reference plasmid enabled accurate epigenetic estimation of absolute leukocyte counts from dried blood samples, correlating with conventional venous (r = 0.86) and capillary (r = 0.80) blood measurements. CONCLUSIONS: The advanced selection of cell-type specific CpGs and utilization of digital droplet PCR analysis provided accurate epigenetic blood counts. Analysis of dried blood facilitates self-sampling with a finger prick, thereby enabling easier accessibility to testing.

First-time application of droplet digital PCR for methylation testing of the 11p15.5 imprinting regions.

BACKGROUND: Beckwith-Wiedemann syndrome and Silver-Russel syndrome are two imprinting disorders caused by opposite molecular alterations in 11p15.5. With the current diagnostic tests, their molecular diagnosis is challenging due to molecular heterogeneity and mosaic occurrence of the most frequent alterations. As the determination of precise (epi)genotype of patients is relevant as the basis for a personalized treatment, different approaches are needed to increase the sensitivity of diagnostic testing of imprinting disorders. METHODS: We established methylation-specific droplet digital PCR approaches (MS-ddPCR) for the two imprinting centers in 11p15.5, and analyzed patients with paternal uniparental disomy of chromosome 11p15.5 (upd(11)pat) and other imprinting defects in the region. The results were compared to those from MS-MLPA (multiplex ligation-dependent probe amplification) and MS-pyrosequencing. RESULTS: MS-ddPCR confirmed the molecular alterations in all patients and the results matched well with MS-MLPA. The results of MS-pyrosequencing varied between different runs, whereas MS-ddPCR results were reproducible. CONCLUSION: We show for the first time that MS-ddPCR is a reliable and easy applicable method for determination of MS-associated changes in imprinting disorders. It is therefore an additional tool for multimethod diagnostics of imprinting disorders suitable to improve the diagnostic yield.

A systematically derived overview of the non-ubiquitous pathways and genes that define the molecular and genetic signature of the healthy trabecular meshwork.

PURPOSE: The trabecular meshwork (TM) is situated in the most frontal part of the eye and is thought to play an important role in the regulation of the eye pressure. However, this tissue is rather difficult to harvest for research. The purpose of this study is therefore to integrate the existing gene expression data of the healthy TM to increase sample size and identify its signature genes and pathways. This provides a robust reference for the study of molecular disease processes and supports the selection of candidate target genes for new treatments. METHODS: A systematic search identified microarray data of healthy TM tissue. After quality control, datasets of low quality and deviating samples were excluded. Remaining individuals were jointly normalized and integrated into one database. The average gene expression of each tested gene over all individuals was calculated. The 25% genes with the highest average expression were identified as the most active genes in the healthy TM and used as input for pathway and network analysis. Additionally, ubiquitous pathways and genes were identified and excluded from the results. Lastly, we identified genes which are likely to be TM-specific. RESULTS: The gene expression data of 44 individuals, obtained from 18 datasets, were jointly normalized. Ubiquitous genes (n = 688) and ubiquitous pathways (n = 73) were identified and excluded. Following, 1882 genes and 211 pathways were identified as the signature genes and pathways of the healthy TM. Pathway analysis revealed multiple molecular processes of which some were already known to be active in the TM, for example extracellular matrix and elastic fiber formation. Forty-six candidate TM-specific genes were identified. These consist mainly of pseudogenes or novel transcripts of which the function is unknown. CONCLUSIONS: In this comprehensive meta-analysis we identified non-ubiquitous genes and pathways that form the signature of the functioning of the healthy TM. Additionally, 46 candidate TM-specific genes were identified. This method can also be used for other tissues that are difficult to obtain for study.

Toward Clinical Application of Leukocyte Counts Based on Targeted DNA Methylation Analysis.

BACKGROUND: Differential leukocyte counts are usually measured based on cellular morphology or surface marker expression. It has recently been shown that leukocyte counts can also be determined by cell-type-specific DNA methylation (DNAm). Such epigenetic leukocyte counting is applicable to small blood volumes and even frozen material, but for clinical translation, the method needs to be further refined and validated. METHODS: We further optimized and validated targeted DNAm assays for leukocyte deconvolution using 332 venous and 122 capillary blood samples from healthy donors. In addition, we tested 36 samples from ring trials and venous blood from 266 patients diagnosed with different hematological diseases. Deconvolution of cell types was determined with various models using DNAm values obtained by pyrosequencing or digital droplet PCR (ddPCR). RESULTS: Relative leukocyte quantification correlated with conventional blood counts for granulocytes, lymphocytes, B cells, T cells (CD4 or CD8), natural killer cells, and monocytes with pyrosequencing (r = 0.84; r = 0.82; r = 0.58; r = 0.50; r = 0.70; r = 0.61; and r = 0.59, respectively) and ddPCR measurements (r = 0.65; r = 0.79; r = 0.56; r = 0.57; r = 0.75; r = 0.49; and r = 0.46, respectively). In some patients, particularly with hematopoietic malignancies, we observed outliers in epigenetic leukocyte counts, which could be discerned if relative proportions of leukocyte subsets did not sum up to 100%. Furthermore, absolute quantification was obtained by spiking blood samples with a reference plasmid of known copy number. CONCLUSIONS: Targeted DNAm analysis by pyrosequencing or ddPCR is a valid alternative to quantify leukocyte subsets, but some assays require further optimization.

Small RNA Sequencing of Aqueous Humor and Plasma in Patients With Primary Open-Angle Glaucoma.

Purpose: Identify differentially expressed microRNAs (miRNAs) in aqueous humor (AH) and blood of primary open-angle glaucoma (POAG) patients by using small RNA sequencing. These may provide insight into POAG pathophysiology or serve as diagnostic biomarker. Methods: AH and plasma of nine POAG patients and 10 cataract control patients were small RNA sequenced on Illumina NovaSeq 6000. Identification of gene transcripts targeted by differentially expressed miRNAs was done with miRWalk and MirPath. These targets were used for pathway analysis and Gene Ontology enrichment. Diagnostic potential was evaluated by receiver operating characteristics analysis. Results: We identified 715 miRNAs in plasma and 62 miRNAs in AH. Plasma miRNA profile did not differ between POAG and control. In contrast, in AH, seven miRNAs were differentially expressed. Hsa-miR-30a-3p, hsa-miR-143-3p, hsa-miR-211-5p, and hsa-miR-221-3p were upregulated, whereas hsa-miR-92a-3p, hsa-miR-451a, and hsa-miR-486-5p were downregulated in POAG. Compared to previous studies, hsa-mir-143-3p, hsa-miR-211-5p, and hsa-miR-221-3p were reported previously, strengthening their involvement in POAG whereas hsa-miR-30a-3p, hsa-miR-92a-3p, and hsa-miR-486-5p are implicated in POAG for the first time. Identified gene transcripts were involved in several pathways, some implicated in glaucoma before (e.g., TGF-ß and neurotrophin signaling), whereas others are new (e.g., prolactin and apelin signaling). In respect to diagnostics, AH concentration of hsa-mir-143-3p had an area under the curve (AUC) of 0.889. Combined with hsa-miR-221-3p, AUC improved to 0.96. Conclusions: Small RNA sequencing identified seven differentially expressed miRNAs in AH of POAG patients. The differentially expressed miRNAs may be useful as POAG biomarkers or could become targets for new therapeutic strategies.

Plasma GDF-15 concentration is not elevated in open-angle glaucoma.

AIM: Recently, the level of growth differentiation factor 15 (GDF-15) in blood, was proposed as biomarker to detect mitochondrial dysfunction. In the current study, we evaluate this biomarker in open-angle glaucoma (OAG), as there is increasing evidence that mitochondrial dysfunction plays a role in the pathophysiology of this disease. METHODS: Plasma GDF-15 concentrations were measured with ELISA in 200 OAG patients and 61 age-matched controls (cataract without glaucoma). The OAG patient group consisted of high tension glaucoma (HTG; n = 162) and normal tension glaucoma (NTG; n = 38). Groups were compared using the Kruskal-Wallis nonparametric test with Dunn's multiple comparison post-hoc correction. GDF-15 concentration was corrected for confounders identified with forward linear regression models. RESULTS: Before correcting for confounders, median plasma GDF-15 levels was significantly lower in the combined OAG group (p = 0.04), but not when analysing HTG and NTG patients separately. Forward linear regression analysis showed that age, gender, smoking and systemic hypertension were significant confounders affecting GDF-15 levels. After correction for these confounders, GDF-15 levels in OAG patients were no longer significantly different from controls. Subgroup analysis of the glaucoma patients did not show a correlation between disease severity and plasma GDF-15, but did reveal that for NTG patients, intake of dietary supplements, which potentially improve mitochondrial function, correlated with lower plasma GDF-15. CONCLUSION: The present study suggests that plasma GDF-15 is not suited as biomarker of mitochondrial dysfunction in OAG patients.

Mapping mRNA Expression of Glaucoma Genes in the Healthy Mouse Eye.

Purpose/Aim: Many genes have been associated with primary open-angle glaucoma (POAG). Knowing exactly where they are expressed in the eye helps to unravel POAG pathology and to select optimal targets for intervention. We investigated whether RNA in situ hybridization (RNA-ISH) is a convenient technique to obtain detailed pan-ocular expression data of these genes. We tested this for four diverse candidate POAG genes, selected because of unclear ocular distribution (F5 and Dusp1) and relevance for potential new therapies (Tnf, Tgfßr3). Optn, a POAG gene with well-known ocular expression pattern served as control. Methods: We made a list of candidate glaucoma genes reported in genetic studies. A table of their ocular expression at the tissue level was compiled using publicly available microarray data (the ocular tissue database). To add cellular detail we performed RNA-ISH for Optn, Tnf, Tgfßr3, F5, and Dusp1 on eyes of healthy, 2-month-old, pigmented, and albino mice. Results: Expression of the Optn control matched with published immunohistochemistry data. Ocular expression of Tnf was generally low, with patches of higher Tnf expression, superficially in the corneal epithelium. F5 had a restricted expression pattern with high expression in the nonpigmented ciliary body epithelium and moderate expression in the peripapillary region. Tgfßr3 and Dusp1 showed ubiquitous expression. Conclusions: RNA-ISH is a suitable technique to determine the ocular expression pattern of POAG genes, adding meaningful cellular detail to existing microarray expression data. For instance, the high expression of F5 in the nonpigmented ciliary body epithelium suggests a role of this gene in aqueous humor dynamics and intraocular pressure. In addition, the ubiquitous expression of Tgfßr3 has implications for designing TGF-ß-related glaucoma therapies, with respect to side effects. Creating pan-ocular expression maps of POAG genes with RNA-ISH will help to identify POAG pathways in specific cell types and to select targets for drug development.