Finite-size effects in transcript sequencing count distribution: its power-law correction necessarily precedes downstream normalization and comparative analysis.

BACKGROUND: Though earlier works on modelling transcript abundance from vertebrates to lower eukaroytes have specifically singled out the Zip's law, the observed distributions often deviate from a single power-law slope. In hindsight, while power-laws of critical phenomena are derived asymptotically under the conditions of infinite observations, real world observations are finite where the finite-size effects will set in to force a power-law distribution into an exponential decay and consequently, manifests as a curvature (i.e., varying exponent values) in a log-log plot. If transcript abundance is truly power-law distributed, the varying exponent signifies changing mathematical moments (e.g., mean, variance) and creates heteroskedasticity which compromises statistical rigor in analysis. The impact of this deviation from the asymptotic power-law on sequencing count data has never truly been examined and quantified. RESULTS: The anecdotal description of transcript abundance being almost Zipf's law-like distributed can be conceptualized as the imperfect mathematical rendition of the Pareto power-law distribution when subjected to the finite-size effects in the real world; This is regardless of the advancement in sequencing technology since sampling is finite in practice. Our conceptualization agrees well with our empirical analysis of two modern day NGS (Next-generation sequencing) datasets: an in-house generated dilution miRNA study of two gastric cancer cell lines (NUGC3 and AGS) and a publicly available spike-in miRNA data; Firstly, the finite-size effects causes the deviations of sequencing count data from Zipf's law and issues of reproducibility in sequencing experiments. Secondly, it manifests as heteroskedasticity among experimental replicates to bring about statistical woes. Surprisingly, a straightforward power-law correction that restores the distribution distortion to a single exponent value can dramatically reduce data heteroskedasticity to invoke an instant increase in signal-to-noise ratio by 50% and the statistical/detection sensitivity by as high as 30% regardless of the downstream mapping and normalization methods. Most importantly, the power-law correction improves concordance in significant calls among different normalization methods of a data series averagely by 22%. When presented with a higher sequence depth (4 times difference), the improvement in concordance is asymmetrical (32% for the higher sequencing depth instance versus 13% for the lower instance) and demonstrates that the simple power-law correction can increase significant detection with higher sequencing depths. Finally, the correction dramatically enhances the statistical conclusions and eludes the metastasis potential of the NUGC3 cell line against AGS of our dilution analysis. CONCLUSIONS: The finite-size effects due to undersampling generally plagues transcript count data with reproducibility issues but can be minimized through a simple power-law correction of the count distribution. This distribution correction has direct implication on the biological interpretation of the study and the rigor of the scientific findings. REVIEWERS: This article was reviewed by Oliviero Carugo, Thomas Dandekar and Sandor Pongor.

Feasibility of low-throughput next generation sequencing for germline DNA screening.

BACKGROUND: Next generation sequencing (NGS) promises many benefits for clinical diagnostics. However, current barriers to its adoption include suboptimal amenability for low clinical throughputs and uncertainty over data accuracy and analytical procedures. We assessed the feasibility and performance of low-throughput NGS for detecting germline mutations for Lynch syndrome (LS). METHODS: Sequencing depth, time, and cost of 6 formats on the MiSeq and Personal Genome Machine platforms at 1-12 samples/run were calculated. Analytical performance was assessed from 3 runs of 3 DNA samples annotated for 7500 nucleotides by BeadChip arrays. The clinical performance of low-throughput NGS and 9 analytical processes were assessed through blinded analysis of DNA samples from 12 LS cases confirmed by Sanger sequencing, and 3 control cases. RESULTS: The feasibility analysis revealed different formats were optimal at different throughputs. Detection was reproducible for 2619/2635 (99.39%) replicate variants, and sensitivity and specificity to array annotation were 99.42% and 99.99% respectively. Eleven of 16 inconsistently detected variants could be specifically identified by having allele frequencies ≤ 0.15, strand biases >-35, or genotype quality scores ≤ 80. Positive selection for variants in the Human Genome Mutation Database (colorectal cancer, nonpolyposis) and variants with ≤ 5% frequency in the Asian population gave the best clinical performance (92% sensitivity, 67% specificity). CONCLUSIONS: Low-throughput NGS can be a cost-efficient and reliable approach for screening germline variants; however, its clinical utility is subject to the quality of annotation of clinically relevant variants.

Association of maternal and nutrient supply line factors with DNA methylation at the imprinted IGF2/H19 locus in multiple tissues of newborn twins.

Epigenetic events are crucial for early development, but can be influenced by environmental factors, potentially programming the genome for later adverse health outcomes. The insulin-like growth factor 2 (IGF2)/H19 locus is crucial for prenatal growth and the epigenetic state at this locus is environmentally labile. Recent studies have implicated maternal factors, including folate intake and smoking, in the regulation of DNA methylation at this locus, although data are often conflicting in the direction and magnitude of effect. Most studies have focused on single tissues and on one or two differentially-methylated regions (DMRs) regulating IGF2/H19 expression. In this study, we investigated the relationship between multiple shared and non-shared gestational/maternal factors and DNA methylation at four IGF2/H19 DMRs in five newborn cell types from 67 pairs of monozygotic and 49 pairs of dizygotic twins. Data on maternal and non-shared supply line factors were collected during the second and third trimesters of pregnancy and DNA methylation was measured via mass spectrometry using Sequenom MassArray EpiTyper analysis. Our exploratory approach showed that the site of umbilical cord insertion into the placenta in monochorionic twins has the strongest positive association with methylation in all IGF2/H19 DMRs (p<0.05). Further, evidence for tissue- and locus-specific effects were observed, emphasizing that responsiveness to environmental exposures in utero cannot be generalized across genes and tissues, potentially accounting for the lack of consistency in previous findings. Such complexity in responsiveness to environmental exposures in utero has implications for all epigenetic studies investigating the developmental origins of health and disease.

A distinct DNA methylation signature defines pediatric pre-B cell acute lymphoblastic leukemia.

Pre-B cell acute lymphoblastic leukemia (ALL) is the most prevalent childhood malignancy and remains one of the highest causes of childhood mortality. Despite this, the mechanisms leading to disease remain poorly understood. We asked if recurrent aberrant DNA methylation plays a role in childhood ALL and have defined a genome-scale DNA methylation profile associated with the ETV6-RUNX1 subtype of pediatric ALL. Archival bone marrow smears from 19 children collected at diagnosis and remission were used to derive a disease specific DNA methylation profile. The gene signature was confirmed in an independent cohort of 86 patients. A further 163 patients were analyzed for DNA methylation of a three gene signature. We found that the DNA methylation signature at diagnosis was unique from remission. Fifteen loci were sufficient to discriminate leukemia from disease-free samples and purified CD34+ cells. DNA methylation of these loci was recurrent irrespective of cytogenetic subtype of pre-B cell ALL. We show that recurrent aberrant genomic methylation is a common feature of pre-B ALL, suggesting a shared pathway for disease development. By revealing new DNA methylation markers associated with disease, this study has identified putative targets for development of novel epigenetic-based therapies.

Placenta-specific methylation of the vitamin D 24-hydroxylase gene: implications for feedback autoregulation of active vitamin D levels at the fetomaternal interface.

Plasma concentrations of biologically active vitamin D (1,25-(OH)(2)D) are tightly controlled via feedback regulation of renal 1alpha-hydroxylase (CYP27B1; positive) and 24-hydroxylase (CYP24A1; catabolic) enzymes. In pregnancy, this regulation is uncoupled, and 1,25-(OH)(2)D levels are significantly elevated, suggesting a role in pregnancy progression. Epigenetic regulation of CYP27B1 and CYP24A1 has previously been described in cell and animal models, and despite emerging evidence for a critical role of epigenetics in placentation generally, little is known about the regulation of enzymes modulating vitamin D homeostasis at the fetomaternal interface. In this study, we investigated the methylation status of genes regulating vitamin D bioavailability and activity in the placenta. No methylation of the VDR (vitamin D receptor) and CYP27B1 genes was found in any placental tissues. In contrast, the CYP24A1 gene is methylated in human placenta, purified cytotrophoblasts, and primary and cultured chorionic villus sampling tissue. No methylation was detected in any somatic human tissue tested. Methylation was also evident in marmoset and mouse placental tissue. All three genes were hypermethylated in choriocarcinoma cell lines, highlighting the role of vitamin D deregulation in this cancer. Gene expression analysis confirmed a reduced capacity for CYP24A1 induction with promoter methylation in primary cells and in vitro reporter analysis demonstrated that promoter methylation directly down-regulates basal promoter activity and abolishes vitamin D-mediated feedback activation. This study strongly suggests that epigenetic decoupling of vitamin D feedback catabolism plays an important role in maximizing active vitamin D bioavailability at the fetomaternal interface.