Genomic origin, fragmentomics, and transcriptional properties of long cell-free DNA molecules in human plasma.

Recent studies have revealed an unexplored population of long cell-free DNA (cfDNA) molecules in human plasma using long-read sequencing technologies. However, the biological properties of long cfDNA molecules (>500 bp) remain largely unknown. To this end, we have investigated the origins of long cfDNA molecules from different genomic elements. Analysis of plasma cfDNA using long-read sequencing reveals an uneven distribution of long molecules from across the genome. Long cfDNA molecules show overrepresentation in euchromatic regions of the genome, in sharp contrast to short DNA molecules. We observe a stronger relationship between the abundance of long molecules and mRNA gene expression levels, compared with short molecules (Pearson's r = 0.71 vs. -0.14). Moreover, long and short molecules show distinct fragmentation patterns surrounding CpG sites. Leveraging the cleavage preferences surrounding CpG sites, the combined cleavage ratios of long and short molecules can differentiate patients with hepatocellular carcinoma (HCC) from non-HCC subjects (AUC = 0.87). We also investigated knockout mice in which selected nuclease genes had been inactivated in comparison with wild-type mice. The proportion of long molecules originating from transcription start sites are lower in Dffb-deficient mice but higher in Dnase1l3-deficient mice compared with that of wild-type mice. This work thus provides new insights into the biological properties and potential clinical applications of long cfDNA molecules.

Universal Targeted Haplotyping by Droplet Digital PCR Sequencing and Its Applications in Noninvasive Prenatal Testing and Pharmacogenetics Analysis.

BACKGROUND: The analysis of haplotypes of variants is important for pharmacogenomics analysis and noninvasive prenatal testing for monogenic diseases. However, there is a lack of robust methods for targeted haplotyping. METHODS: We developed digital PCR haplotype sequencing (dHapSeq) for targeted haplotyping of variants, which is a method that compartmentalizes long DNA molecules into droplets. Within one droplet, 2 target regions are PCR amplified from one template molecule, and their amplicons are fused together. The fused products are then sequenced to determine the phase relationship of the single nucleotide polymorphism (SNP) alleles. The entire haplotype of 10s of SNPs can be deduced after the phase relationship of individual SNPs are determined in a pairwise manner. We applied dHapSeq to noninvasive prenatal testing in 4 families at risk for thalassemia and utilized it to detect NUDT15 diplotypes for predicting drug tolerance in pediatric acute lymphoblastic leukemia (72 cases and 506 controls). RESULTS: For SNPs within 40 kb, phase relation can be determined with 100% accuracy. In 7 trio families, the haplotyping results for 97 SNPs spanning 185 kb determined by dHapSeq were concordant with the results deduced from the genotypes of both parents and the fetus. In 4 thalassemia families, a 19.3-kb Southeast Asian deletion was successfully phased with 97 downstream SNPs, enabling noninvasive determination of fetal inheritance using relative haplotype dosage analysis. In the NUDT15 analysis, the variant status and phase of the variants were successfully determined in all cases and controls. CONCLUSIONS: The dHapSeq represents a robust and scalable haplotyping approach with numerous clinical and research applications.

Methylation-Associated Nucleosomal Patterns of Cell-Free DNA in Cancer Patients and Pregnant Women.

BACKGROUND: Cell-free DNA (cfDNA) analysis offers an attractive noninvasive means of detecting and monitoring diseases. cfDNA cleavage patterns within a short range (e.g., 11 nucleotides) have been reported to correlate with cytosine-phosphate-guanine (CpG) methylation, allowing fragmentomics-based methylation analysis (FRAGMA). Here, we adopted FRAGMA to the extended region harboring multiple nucleosomes, termed FRAGMAXR. METHODS: We profiled cfDNA nucleosomal patterns over the genomic regions from -800 to 800 bp surrounding differentially methylated CpG sites, harboring approximately 8 nucleosomes, referred to as CpG-associated cfDNA nucleosomal patterns. Such nucleosomal patterns were analyzed by FRAGMAXR in cancer patients and pregnant women. RESULTS: We identified distinct cfDNA nucleosomal patterns around differentially methylated CpG sites. Compared with subjects without cancer, patients with hepatocellular carcinoma (HCC) showed reduced amplitude of nucleosomal patterns, with a gradual decrease over tumor stages. Nucleosomal patterns associated with differentially methylated CpG sites could be used to train a machine learning model, resulting in the detection of HCC patients with an area under the receiver operating characteristic curve of 0.93. We further demonstrated the feasibility of multicancer detection using a dataset comprising lung, breast, and ovarian cancers. The tissue-of-origin analysis of plasma cfDNA from pregnant women and cancer patients revealed that the placental DNA and tumoral DNA contributions deduced by FRAGMAXR correlated well with values measured using genetic variants (Pearson r: 0.85 and 0.94, respectively). CONCLUSIONS: CpG-associated cfDNA nucleosomal patterns of cfDNA molecules are influenced by DNA methylation and might be useful for biomarker developments for cancer liquid biopsy and noninvasive prenatal testing.

Single-molecule sequencing reveals a large population of long cell-free DNA molecules in maternal plasma.

In the field of circulating cell-free DNA, most of the studies have focused on short DNA molecules (e.g., <500 bp). The existence of long cell-free DNA molecules has been poorly explored. In this study, we demonstrated that single-molecule real-time sequencing allowed us to detect and analyze a substantial proportion of long DNA molecules from both fetal and maternal sources in maternal plasma. Such molecules were beyond the size detection limits of short-read sequencing technologies. The proportions of long cell-free DNA molecules in maternal plasma over 500 bp were 15.5%, 19.8%, and 32.3% for the first, second, and third trimesters, respectively. The longest fetal-derived plasma DNA molecule observed was 23,635 bp. Long plasma DNA molecules demonstrated predominance of A or G 5' fragment ends. Pregnancies with preeclampsia demonstrated a reduction in long maternal plasma DNA molecules, reduced frequencies for selected 5' 4-mer end motifs ending with G or A, and increased frequencies for selected motifs ending with T or C. Finally, we have developed an approach that employs the analysis of methylation patterns of the series of CpG sites on a long DNA molecule for determining its tissue origin. This approach achieved an area under the curve of 0.88 in differentiating between fetal and maternal plasma DNA molecules, enabling the determination of maternal inheritance and recombination events in the fetal genome. This work opens up potential clinical utilities of long cell-free DNA analysis in maternal plasma including noninvasive prenatal testing of monogenic diseases and detection/monitoring of pregnancy-associated disorders such as preeclampsia.

Genome-wide detection of cytosine methylation by single molecule real-time sequencing.

5-Methylcytosine (5mC) is an important type of epigenetic modification. Bisulfite sequencing (BS-seq) has limitations, such as severe DNA degradation. Using single molecule real-time sequencing, we developed a methodology to directly examine 5mC. This approach holistically examined kinetic signals of a DNA polymerase (including interpulse duration and pulse width) and sequence context for every nucleotide within a measurement window, termed the holistic kinetic (HK) model. The measurement window of each analyzed double-stranded DNA molecule comprised 21 nucleotides with a cytosine in a CpG site in the center. We used amplified DNA (unmethylated) and M.SssI-treated DNA (methylated) (M.SssI being a CpG methyltransferase) to train a convolutional neural network. The area under the curve for differentiating methylation states using such samples was up to 0.97. The sensitivity and specificity for genome-wide 5mC detection at single-base resolution reached 90% and 94%, respectively. The HK model was then tested on human-mouse hybrid fragments in which each member of the hybrid had a different methylation status. The model was also tested on human genomic DNA molecules extracted from various biological samples, such as buffy coat, placental, and tumoral tissues. The overall methylation levels deduced by the HK model were well correlated with those by BS-seq (r = 0.99; P < 0.0001) and allowed the measurement of allele-specific methylation patterns in imprinted genes. Taken together, this methodology has provided a system for simultaneous genome-wide genetic and epigenetic analyses.

Detection and characterization of jagged ends of double-stranded DNA in plasma.

Cell-free DNA in plasma has been used for noninvasive prenatal testing and cancer liquid biopsy. The physical properties of cell-free DNA fragments in plasma, such as fragment sizes and ends, have attracted much recent interest, leading to the emerging field of cell-free DNA fragmentomics. However, one aspect of plasma DNA fragmentomics as to whether double-stranded plasma molecules might carry single-stranded ends, termed a jagged end in this study, remains underexplored. We have developed two approaches for investigating the presence of jagged ends in a plasma DNA pool. These approaches utilized DNA end repair to introduce differential methylation signals between the original sequence and the jagged ends, depending on whether unmethylated or methylated cytosines were used in the DNA end-repair procedure. The majority of plasma DNA molecules (87.8%) were found to bear jagged ends. The jaggedness varied according to plasma DNA fragment sizes and appeared to be in association with nucleosomal patterns. In the plasma of pregnant women, the jaggedness of fetal DNA molecules was higher than that of the maternal counterparts. The jaggedness of plasma DNA correlated with the fetal DNA fraction. Similarly, in the plasma of cancer patients, tumor-derived DNA molecules in patients with hepatocellular carcinoma showed an elevated jaggedness compared with nontumoral DNA. In mouse models, knocking out of the Dnase1 gene reduced jaggedness, whereas knocking out of the Dnase1l3 gene enhanced jaggedness. Hence, plasma DNA jagged ends represent an intrinsic property of plasma DNA and provide a link between nuclease activities and the fragmentation of plasma DNA.

Comparison of Single Molecule, Real-Time Sequencing and Nanopore Sequencing for Analysis of the Size, End-Motif, and Tissue-of-Origin of Long Cell-Free DNA in Plasma.

BACKGROUND: Recent studies using single molecule, real-time (SMRT) sequencing revealed a substantial population of analyzable long cell-free DNA (cfDNA) in plasma. Potential clinical utilities of such long cfDNA in pregnancy and cancer have been demonstrated. However, the performance of different long-read sequencing platforms for the analysis of long cfDNA remains unknown. METHODS: Size biases of SMRT sequencing by Pacific Biosciences (PacBio) and nanopore sequencing by Oxford Nanopore Technologies (ONT) were evaluated using artificial mixtures of sonicated human and mouse DNA of different sizes. cfDNA from plasma samples of pregnant women at different trimesters, hepatitis B carriers, and patients with hepatocellular carcinoma were sequenced with the 2 platforms. RESULTS: Both platforms showed biases to sequence longer (1500 bp vs 200 bp) DNA fragments, with PacBio showing a stronger bias (5-fold overrepresentation of long fragments vs 2-fold in ONT). Percentages of cfDNA fragments 500 bp were around 6-fold higher in PacBio compared with ONT. End motif profiles of cfDNA from PacBio and ONT were similar, yet exhibited platform-dependent patterns. Tissue-of-origin analysis based on single-molecule methylation patterns showed comparable performance on both platforms. CONCLUSIONS: SMRT sequencing generated data with higher percentages of long cfDNA compared with nanopore sequencing. Yet, a higher number of long cfDNA fragments eligible for the tissue-of-origin analysis could be obtained from nanopore sequencing due to its much higher throughput. When analyzing the size and end motif of cfDNA, one should be aware of the analytical characteristics and possible biases of the sequencing platforms being used.

Droplet digital PCR is a cost-effective method for analyzing long cell-free DNA in maternal plasma: Application in preeclampsia.

OBJECTIVE: Long cell-free DNA (cfDNA) can be found in the plasma of pregnant women and cancer patients. We investigated if droplet digital PCR (ddPCR) can analyze such molecules for diagnostic purposes using preeclampsia as a model. METHOD: Plasma samples from ten preeclamptic and sixteen normal pregnancies were analyzed. Two ddPCR assays targeting a single-copy gene, VCP, and one ddPCR assay targeting LINE-1 repetitive regions were used to measure the percentages of long cfDNA >533, 1001, and 170 bp, respectively. The LINE-1 assay was developed as guided by in silico PCR analyses to better differentiate preeclamptic and normal pregnancies. RESULTS: Preeclamptic patients had a significantly lower median percentage of long cfDNA than healthy pregnant controls, as determined by the LINE-1 170 bp assay (28.9% vs. 35.1%, p < 0.0001) and the VCP 533 bp assay (6.6% vs. 8.7%, p = 0.014). The LINE-1 assay provided a better differentiation than the VCP 533 bp assay (area under ROC curves, 0.94 vs. 0.79). CONCLUSION: ddPCR is a cost-effective approach for unlocking diagnostic information carried by long cfDNA in plasma and may have applications for the detection of preeclampsia. Further longitudinal studies with larger cohorts are required to assess the clinical utility of this test.

Identification and characterization of extrachromosomal circular DNA in maternal plasma.

We explored the presence of extrachromosomal circular DNA (eccDNA) in the plasma of pregnant women. Through sequencing following either restriction enzyme or Tn5 transposase treatment, we identified eccDNA molecules in the plasma of pregnant women. These eccDNA molecules showed bimodal size distributions peaking at ∼202 and ∼338 bp with distinct 10-bp periodicity observed throughout the size ranges within both peaks, suggestive of their nucleosomal origin. Also, the predominance of the 338-bp peak of eccDNA indicated that eccDNA had a larger size distribution than linear DNA in human plasma. Moreover, eccDNA of fetal origin were shorter than the maternal eccDNA. Genomic annotation of the overall population of eccDNA molecules revealed a preference of these molecules to be generated from 5'-untranslated regions (5'-UTRs), exonic regions, and CpG island regions. Two sets of trinucleotide repeat motifs flanking the junctional sites of eccDNA supported multiple possible models for eccDNA generation. This work highlights the topologic analysis of plasma DNA, which is an emerging direction for circulating nucleic acid research and applications.

Orientation-aware plasma cell-free DNA fragmentation analysis in open chromatin regions informs tissue of origin.

Cell-free DNA (cfDNA) in human plasma is a class of biomarkers with many current and potential future diagnostic applications. Recent studies have shown that cfDNA molecules are not randomly fragmented and possess information related to their tissues of origin. Pathologies causing death of cells from particular tissues result in perturbations in the relative distribution of DNA from the affected tissues. Such tissue-of-origin analysis is particularly useful in the development of liquid biopsies for cancer. It is therefore of value to accurately determine the relative contributions of the tissues to the plasma DNA pool in a simultaneous manner. In this work, we report that in open chromatin regions, cfDNA molecules show characteristic fragmentation patterns reflected by sequencing coverage imbalance and differentially phased fragment end signals. The latter refers to differences in the read densities of sequences corresponding to the orientation of the upstream and downstream ends of cfDNA molecules in relation to the reference genome. Such cfDNA fragmentation patterns preferentially occur in tissue-specific open chromatin regions where the corresponding tissues contributed DNA into the plasma. Quantitative analyses of such signals allow measurement of the relative contributions of various tissues toward the plasma DNA pool. These findings were validated by plasma DNA sequencing data obtained from pregnant women, organ transplantation recipients, and cancer patients. Orientation-aware plasma DNA fragmentation analysis therefore has potential diagnostic applications in noninvasive prenatal testing, organ transplantation monitoring, and cancer liquid biopsy.

Characteristics of Fetal Extrachromosomal Circular DNA in Maternal Plasma: Methylation Status and Clearance.

BACKGROUND: Although the characterization of cell-free extrachromosomal circular DNA (eccDNA) has gained much research interest, the methylation status of these molecules is yet to be elucidated. We set out to compare the methylation densities of plasma eccDNA of maternal and fetal origins, and between small and large molecules. The clearance of fetal eccDNA from maternal circulation was also investigated. METHODS: We developed a sequencing protocol for eccDNA methylation analysis using tagmentation and enzymatic conversion approaches. A restriction enzyme-based approach was applied to verify the tagmentation results. The efficiency of cell-free fetal eccDNA clearance was investigated by fetal eccDNA fraction evaluations at various postpartum time points. RESULTS: The methylation densities of fetal eccDNA (median: 56.3%; range: 40.5-67.6%) were lower than the maternal eccDNA (median: 66.7%; range: 56.5-75.7%) (P = 0.02, paired t-test). In addition, eccDNA molecules from the smaller peak cluster (180-230 bp) were of lower methylation levels than those from the larger peak cluster (300-450 bp). Both of these findings were confirmed using the restriction enzyme approach. We also observed comparable methylation densities between linear and eccDNA of both maternal and fetal origins. The average half-lives of fetal linear and eccDNA in the maternal blood were 30.2 and 29.7 min, respectively. CONCLUSIONS: We found that fetal eccDNA in plasma was relatively hypomethylated compared to the maternal eccDNA. The methylation densities of eccDNA were positively correlated with their sizes. In addition, fetal eccDNA was found to be rapidly cleared from the maternal blood after delivery, similar to fetal linear DNA.

Bradycardia-to-delivery interval and fetal outcomes in umbilical cord prolapse.

INTRODUCTION: Umbilical cord prolapse is a major obstetric emergency associated with significant perinatal complications. However, there is no consensus on the optimal decision-to-delivery interval, as many previous studies have shown poor correlation between the interval and umbilical cord arterial blood gas or perinatal outcomes. We aim to investigate whether bradycardia-to-delivery or decision-to-delivery interval was related to poor cord arterial pH or adverse perinatal outcome in umbilical cord prolapse. MATERIAL AND METHODS: This was a retrospective study conducted at a university tertiary obstetric unit in Hong Kong. All women with singleton pregnancy complicated by cord prolapse during labor between 1995 and 2018 were included. Women were categorized into three groups. Group 1: persistent bradycardia; Group 2: any type of decelerations without bradycardia; and Group 3: normal fetal heart rate. The main outcome was cord arterial blood gas results of the newborns in different groups. Maternal demographic data and perinatal outcomes were reviewed. Correlation analysis between cord arterial blood gas result and time intervals including bradycardia-to-delivery, deceleration-to-delivery, and decision-to-delivery were performed for the different groups with Spearman test. RESULTS: There were 34, 30, and 50 women in Groups 1, 2, and 3, respectively. Cord arterial pH and base excess did not correlate with decision-to-delivery interval in any of the groups, but they were inversely correlated with bradycardia-to-delivery interval in Group 1 (Spearman's ρ = -.349; P = .043 and Spearman's ρ = -.558; P = .001, respectively). The cord arterial pH drops at 0.009 per minute with bradycardia-to-delivery interval in Group 1 (95% CI 0.0180-0.0003). The risk of significant acidosis (pH < 7) was 80% when bradycardia-to-delivery interval was >20 minutes, and 17.2% when the interval was <20 minutes. CONCLUSIONS: There is significant correlation between bradycardia-to-delivery interval and cord arterial pH in umbilical cord prolapse with fetal bradycardia but not in cases with decelerations or normal heart rate. The drop of cord arterial pH is rapid and urgent delivery is essential in such situations.

Size-tagged preferred ends in maternal plasma DNA shed light on the production mechanism and show utility in noninvasive prenatal testing.

Cell-free DNA in human plasma is nonrandomly fragmented and reflects genomewide nucleosomal organization. Previous studies had demonstrated tissue-specific preferred end sites in plasma DNA of pregnant women. In this study, we performed integrative analysis of preferred end sites with the size characteristics of plasma DNA fragments. We mined the preferred end sites in short and long plasma DNA molecules separately and found that these "size-tagged" ends showed improved accuracy in fetal DNA fraction estimation and enhanced noninvasive fetal trisomy 21 testing. Further analysis revealed that the fetal and maternal preferred ends were generated from different locations within the nucleosomal structure. Hence, fetal DNA was frequently cut within the nucleosome core while maternal DNA was mostly cut within the linker region. We further demonstrated that the nucleosome accessibility in placental cells was higher than that for white blood cells, which might explain the difference in the cutting positions and the shortness of fetal DNA in maternal plasma. Interestingly, short and long size-tagged ends were also observable in the plasma of nonpregnant healthy subjects and demonstrated size differences similar to those in the pregnant samples. Because the nonpregnant samples did not contain fetal DNA, the data suggested that the interrelationship of preferred DNA ends, chromatin accessibility, and plasma DNA size profile is likely a general one, extending beyond the context of pregnancy. Plasma DNA fragment end patterns have thus shed light on production mechanisms and show utility in future developments in plasma DNA-based noninvasive molecular diagnostics.

Topologic Analysis of Plasma Mitochondrial DNA Reveals the Coexistence of Both Linear and Circular Molecules.

BACKGROUND: Cellular mitochondrial DNA (mtDNA) is organized as circular, covalently closed and double-stranded DNA. Studies have demonstrated the presence of short mtDNA fragments in plasma. It is not known whether circular mtDNA might concurrently exist with linear mtDNA in plasma. METHODS: We elucidated the topology of plasma mtDNA using restriction enzyme BfaI cleavage signatures on mtDNA fragment ends to differentiate linear and circular mtDNA. mtDNA fragments with both ends carrying BfaI cleavage signatures were defined as circular-derived mtDNA, whereas those with no cleavage signature or with 1 cleavage signature were defined as linear-derived mtDNA. An independent assay using exonuclease V to remove linear DNA followed by restriction enzyme MspI digestion was used for confirming the conclusions based on BfaI cleavage analysis. We analyzed the presence of BfaI cleavage signatures on plasma DNA ends in nonhematopoietically and hematopoietically derived DNA molecules by sequencing plasma DNA of patients with liver transplantation and bone marrow transplantation. RESULTS: Both linear and circular mtDNA coexisted in plasma. In patients with liver transplantation, donor-derived (i.e., liver) mtDNA molecules were mainly linear (median fraction, 91%; range, 75%-97%), whereas recipient-derived (i.e., hematopoietic) mtDNA molecules were mainly circular (median fraction, 88%; range, 77%-93%). The proportion of linear mtDNA was well correlated with liver DNA contribution in the plasma DNA pool (r = 0.83; P value = 0.0008). Consistent data were obtained from a bone marrow transplantation recipient in whom the donor-derived (i.e., hematopoietic) mtDNA molecules were predominantly circular. CONCLUSIONS: Linear and circular mtDNA molecules coexist in plasma and may have different tissue origins.

Second generation noninvasive fetal genome analysis reveals de novo mutations, single-base parental inheritance, and preferred DNA ends.

Plasma DNA obtained from a pregnant woman was sequenced to a depth of 270× haploid genome coverage. Comparing the maternal plasma DNA sequencing data with the parental genomic DNA data and using a series of bioinformatics filters, fetal de novo mutations were detected at a sensitivity of 85% and a positive predictive value of 74%. These results represent a 169-fold improvement in the positive predictive value over previous attempts. Improvements in the interpretation of the sequence information of every base position in the genome allowed us to interrogate the maternal inheritance of the fetus for 618,271 of 656,676 (94.2%) heterozygous SNPs within the maternal genome. The fetal genotype at each of these sites was deduced individually, unlike previously, where the inheritance was determined for a collection of sites within a haplotype. These results represent a 90-fold enhancement in the resolution in determining the fetus's maternal inheritance. Selected genomic locations were more likely to be found at the ends of plasma DNA molecules. We found that a subset of such preferred ends exhibited selectivity for fetal- or maternal-derived DNA in maternal plasma. The ratio of the number of maternal plasma DNA molecules with fetal preferred ends to those with maternal preferred ends showed a correlation with the fetal DNA fraction. Finally, this second generation approach for noninvasive fetal whole-genome analysis was validated in a pregnancy diagnosed with cardiofaciocutaneous syndrome with maternal plasma DNA sequenced to 195× coverage. The causative de novo BRAF mutation was successfully detected through the maternal plasma DNA analysis.

Noninvasive reconstruction of placental methylome from maternal plasma DNA: Potential for prenatal testing and monitoring.

OBJECTIVE: During human pregnancy, the DNA methylation of placental tissue is highly relevant to the normal growth and development of the fetus; therefore, methylomic analysis of the placental tissue possesses high research and clinical value in prenatal testing and monitoring. Thus, our aim is to develop an approach for reconstruction of the placental methylome, which should be completely noninvasive and achieve high accuracy and resolution. RESULTS: We propose a novel size-based algorithm, FEtal MEthylome Reconstructor (FEMER), to noninvasively reconstruct the placental methylome by genomewide bisulfite sequencing and size-based analysis of maternal plasma DNA. By applying FEMER on a real clinical dataset, we demonstrate that FEMER achieves both high accuracy and resolution, thus provides a high-quality view of the placental methylome from maternal plasma DNA. FEtal MEthylome Reconstructor could also predict the DNA methylation profile of CpG islands with high accuracy, thus shows potential in monitoring of key genes involved in placental/fetal development. Source code and testing datasets for FEMER are available at http://sunlab.cpy.cuhk.edu.hk/FEMER/. CONCLUSION: FEtal MEthylome Reconstructor could enhance the noninvasive fetal/placental methylomic analysis and facilitate its application in prenatal testing and monitoring.

Plasma DNA tissue mapping by genome-wide methylation sequencing for noninvasive prenatal, cancer, and transplantation assessments.

Plasma consists of DNA released from multiple tissues within the body. Using genome-wide bisulfite sequencing of plasma DNA and deconvolution of the sequencing data with reference to methylation profiles of different tissues, we developed a general approach for studying the major tissue contributors to the circulating DNA pool. We tested this method in pregnant women, patients with hepatocellular carcinoma, and subjects following bone marrow and liver transplantation. In most subjects, white blood cells were the predominant contributors to the circulating DNA pool. The placental contributions in the plasma of pregnant women correlated with the proportional contributions as revealed by fetal-specific genetic markers. The graft-derived contributions to the plasma in the transplant recipients correlated with those determined using donor-specific genetic markers. Patients with hepatocellular carcinoma showed elevated plasma DNA contributions from the liver, which correlated with measurements made using tumor-associated copy number aberrations. In hepatocellular carcinoma patients and in pregnant women exhibiting copy number aberrations in plasma, comparison of methylation deconvolution results using genomic regions with different copy number status pinpointed the tissue type responsible for the aberrations. In a pregnant woman diagnosed as having follicular lymphoma during pregnancy, methylation deconvolution indicated a grossly elevated contribution from B cells into the plasma DNA pool and localized B cells as the origin of the copy number aberrations observed in plasma. This method may serve as a powerful tool for assessing a wide range of physiological and pathological conditions based on the identification of perturbed proportional contributions of different tissues into plasma.

Combined Count- and Size-Based Analysis of Maternal Plasma DNA for Noninvasive Prenatal Detection of Fetal Subchromosomal Aberrations Facilitates Elucidation of the Fetal and/or Maternal Origin of the Aberrations.

BACKGROUND: Noninvasive prenatal detection of fetal subchromosomal copy number aberrations (CNAs) can be achieved through massively parallel sequencing of maternal plasma DNA. However, when a mother herself is a carrier of a CNA, one cannot discern if her fetus has inherited the CNA. In addition, false-positive results would become more prevalent when more subchromosomal regions are analyzed. METHODS: We used a strategy that combined count- and size-based analyses of maternal plasma DNA for the detection of fetal subchromosomal CNAs in 7 target regions for 10 test cases. RESULTS: For the 5 cases in which CNAs were present only in the fetus, the size-based approach confirmed the aberrations detected by the count-based approach. For the 5 cases in which the mother herself carried an aberration, we successfully deduced that 3 of the fetuses had inherited the aberrations and that the other 2 fetuses had not inherited the aberrations. No false positives were observed in this cohort. CONCLUSIONS: Combined count- and size-based analysis of maternal plasma DNA permits the noninvasive elucidation of whether a fetus has inherited a CNA from its mother who herself is a carrier of the CNA. This strategy has the potential to improve the diagnostic specificity of noninvasive prenatal testing.

Universal Haplotype-Based Noninvasive Prenatal Testing for Single Gene Diseases.

BACKGROUND: Researchers have developed approaches for the noninvasive prenatal testing of single gene diseases. One approach that allows for the noninvasive assessment of both maternally and paternally inherited mutations involves the analysis of single nucleotide polymorphisms (SNPs) in maternal plasma DNA with reference to parental haplotype information. In the past, parental haplotypes were resolved by complex experimental methods or inferential approaches, such as through the analysis of DNA from other affected family members. Recently, microfluidics-based linked-read sequencing technology has become available and allows the direct haplotype phasing of the whole genome rapidly. We explored the feasibility of applying this direct haplotyping technology in noninvasive prenatal testing. METHODS: We first resolved the haplotypes of parental genomes with the use of linked-read sequencing technology. Then, we identified SNPs within and flanking the genes of interest in maternal plasma DNA by targeted sequencing. Finally, we applied relative haplotype dosage analysis to deduce the mutation inheritance status of the fetus. RESULTS: Haplotype phasing and relative haplotype dosage analysis of 12 out of 13 families were successfully achieved. The mutational status of these 12 fetuses was correctly classified. CONCLUSIONS: High-throughput linked-read sequencing followed by maternal plasma-based relative haplotype dosage analysis represents a streamlined approach for noninvasive prenatal testing of inherited single gene diseases. The approach bypasses the need for mutation-specific assays and is not dependent on the availability of DNA from other affected family members. Thus, the approach is universally applicable to pregnancies at risk for the inheritance of a single gene disease.

Size-based molecular diagnostics using plasma DNA for noninvasive prenatal testing.

Noninvasive prenatal testing using fetal DNA in maternal plasma is an actively researched area. The current generation of tests using massively parallel sequencing is based on counting plasma DNA sequences originating from different genomic regions. In this study, we explored a different approach that is based on the use of DNA fragment size as a diagnostic parameter. This approach is dependent on the fact that circulating fetal DNA molecules are generally shorter than the corresponding maternal DNA molecules. First, we performed plasma DNA size analysis using paired-end massively parallel sequencing and microchip-based capillary electrophoresis. We demonstrated that the fetal DNA fraction in maternal plasma could be deduced from the overall size distribution of maternal plasma DNA. The fetal DNA fraction is a critical parameter affecting the accuracy of noninvasive prenatal testing using maternal plasma DNA. Second, we showed that fetal chromosomal aneuploidy could be detected by observing an aberrant proportion of short fragments from an aneuploid chromosome in the paired-end sequencing data. Using this approach, we detected fetal trisomy 21 and trisomy 18 with 100% sensitivity (T21: 36/36; T18: 27/27) and 100% specificity (non-T21: 88/88; non-T18: 97/97). For trisomy 13, the sensitivity and specificity were 95.2% (20/21) and 99% (102/103), respectively. For monosomy X, the sensitivity and specificity were both 100% (10/10 and 8/8). Thus, this study establishes the principle of size-based molecular diagnostics using plasma DNA. This approach has potential applications beyond noninvasive prenatal testing to areas such as oncology and transplantation monitoring.