Epigenetic analysis of cell-free DNA by fragmentomic profiling.

Cell-free DNA (cfDNA) fragmentation patterns contain important molecular information linked to tissues of origin. We explored the possibility of using fragmentation patterns to predict cytosine-phosphate-guanine (CpG) methylation of cfDNA, obviating the use of bisulfite treatment and associated risks of DNA degradation. This study investigated the cfDNA cleavage profile surrounding a CpG (i.e., within an 11-nucleotide [nt] window) to analyze cfDNA methylation. The cfDNA cleavage proportion across positions within the window appeared nonrandom and exhibited correlation with methylation status. The mean cleavage proportion was ∼twofold higher at the cytosine of methylated CpGs than unmethylated ones in healthy controls. In contrast, the mean cleavage proportion rapidly decreased at the 1-nt position immediately preceding methylated CpGs. Such differential cleavages resulted in a characteristic change in relative presentations of CGN and NCG motifs at 5' ends, where N represented any nucleotide. CGN/NCG motif ratios were correlated with methylation levels at tissue-specific methylated CpGs (e.g., placenta or liver) (Pearson's absolute r > 0.86). cfDNA cleavage profiles were thus informative for cfDNA methylation and tissue-of-origin analyses. Using CG-containing end motifs, we achieved an area under a receiver operating characteristic curve (AUC) of 0.98 in differentiating patients with and without hepatocellular carcinoma and enhanced the positive predictive value of nasopharyngeal carcinoma screening (from 19.6 to 26.8%). Furthermore, we elucidated the feasibility of using cfDNA cleavage patterns to deduce CpG methylation at single CpG resolution using a deep learning algorithm and achieved an AUC of 0.93. FRAGmentomics-based Methylation Analysis (FRAGMA) presents many possibilities for noninvasive prenatal, cancer, and organ transplantation assessment.

Single-Molecule Sequencing Enables Long Cell-Free DNA Detection and Direct Methylation Analysis for Cancer Patients.

BACKGROUND: Analysis of circulating tumor DNA has become increasingly important as a tool for cancer care. However, the focus of previous studies has been on short fragments of DNA. Also, bisulfite sequencing, a conventional approach for methylation analysis, causes DNA degradation, which is not ideal for the assessment of long DNA properties and methylation patterns. This study attempted to overcome such obstacles by single-molecule sequencing. METHODS: Single-molecule real-time (SMRT) sequencing was used to sequence plasma DNA. We performed fragment size and direct methylation analysis for each molecule. A methylation score concerning single-molecule methylation patterns was used for cancer detection. RESULTS: A substantial proportion of plasma DNA was longer than 1 kb with a median of 16% in hepatocellular carcinoma (HCC) patients, hepatitis B virus carriers, and healthy individuals. The longest plasma DNA molecule in the HCC patients was 39.8 kb. Tumoral cell-free DNA (cfDNA) was generally shorter than nontumoral cfDNA. The longest tumoral cfDNA was 13.6 kb. Tumoral cfDNA had lower methylation levels compared with nontumoral cfDNA (median: 59.3% vs 76.9%). We developed and analyzed a metric reflecting single-molecule methylation patterns associated with cancer, named the HCC methylation score. HCC patients displayed significantly higher HCC methylation scores than those without HCC. Interestingly, compared to using short cfDNA (area under the receiver operating characteristic [ROC] curve, AUC: 0.75), the use of long cfDNA molecules greatly enhanced the discriminatory power (AUC: 0.91). CONCLUSIONS: A previously unidentified long cfDNA population was revealed in cancer patients. The presence and direct methylation analysis of these molecules open new possibilities for cancer liquid biopsy.

N,O-Bis(trimethylsilyl)trifluoroacetamide as an Effective Interface Film Additive on Lithium Anodes.

Lithium anodes have attracted much attention because of their high energy density, but the existence of lithium dendrites tremendously limits their practical application. Herein, it is creatively proposed to employ N,O-bis(trimethylsilyl)trifluoroacetamide (BSTFA) as an electrolyte additive to stabilize the solid electrolyte interface. BSTFA is reduced on the lithium anode surface prior to other components to form a passivation layer composed of LiF, Li3N, and SiOx, which not only significantly prevents the continuous consumption of the electrolyte and reduces side reactions but also effectively promotes the uniform deposition of lithium ions with fast Li+ transmission, thereby solving the problem of lithium dendrites. Electrochemical results indicate that BSTFA can obviously reduce polarization in a Li||Li battery at a current density of 1 mA cm-2. Besides, an excellent cycling performance (107 mA h g-1) and Coulombic efficiency (99%) can be obtained for a Li||LiNi0.6Co0.2Mn0.2O2 (NCM622) battery with 0.5 wt % BSTFA at 2 C after 200 cycles, even at a high NCM622 loading of 6 mg cm-2.

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.

Plasma DNA End-Motif Profiling as a Fragmentomic Marker in Cancer, Pregnancy, and Transplantation.

Plasma DNA fragmentomics is an emerging area of research covering plasma DNA sizes, end points, and nucleosome footprints. In the present study, we found a significant increase in the diversity of plasma DNA end motifs in patients with hepatocellular carcinoma (HCC). Compared with patients without HCC, patients with HCC showed a preferential pattern of 4-mer end motifs. In particular, the abundance of plasma DNA motif CCCA was much lower in patients with HCC than in subjects without HCC. The aberrant end motifs were also observed in patients with other cancer types, including colorectal cancer, lung cancer, nasopharyngeal carcinoma, and head and neck squamous cell carcinoma. We further observed that the profile of plasma DNA end motifs originating from the same organ, such as the liver, placenta, and hematopoietic cells, generally clustered together. The profile of end motifs may therefore serve as a class of biomarkers for liquid biopsy in oncology, noninvasive prenatal testing, and transplantation monitoring. SIGNIFICANCE: Plasma DNA molecules originating from the liver, HCC and other cancers, placenta, and hematopoietic cells each harbor a set of characteristic plasma DNA end motifs. Such markers carry tissue-of-origin information and represent a new class of biomarkers in the nascent field of fragmentomics.This article is highlighted in the In This Issue feature, p. 627.

Methylation analysis of plasma DNA informs etiologies of Epstein-Barr virus-associated diseases.

Epstein-Barr virus (EBV) is associated with a number of diseases, including malignancies. Currently, it is not known whether patients with different EBV-associated diseases have different methylation profiles of circulating EBV DNA. Through whole-genome methylation analysis of plasma samples from patients with nasopharyngeal carcinoma (NPC), EBV-associated lymphoma and infectious mononucleosis, we demonstrate that EBV DNA methylation profiles exhibit a disease-associated pattern. This observation implies a significant potential for the development of methylation analysis of plasma EBV DNA for NPC diagnostics. We further analyse the plasma EBV DNA methylome of NPC and non-NPC subjects from a prospective screening cohort. Plasma EBV DNA fragments demonstrate differential methylation patterns between NPC and non-NPC subjects. Combining such differential methylation patterns with the fractional concentration (count) and size of plasma EBV DNA, population screening of NPC is performed with an improved positive predictive value of 35.1%, compared to a count- and size-based only protocol.

A Peripheral Blood Gene Expression Signature to Diagnose Subclinical Acute Rejection.

BACKGROUND: In kidney transplant recipients, surveillance biopsies can reveal, despite stable graft function, histologic features of acute rejection and borderline changes that are associated with undesirable graft outcomes. Noninvasive biomarkers of subclinical acute rejection are needed to avoid the risks and costs associated with repeated biopsies. METHODS: We examined subclinical histologic and functional changes in kidney transplant recipients from the prospective Genomics of Chronic Allograft Rejection (GoCAR) study who underwent surveillance biopsies over 2 years, identifying those with subclinical or borderline acute cellular rejection (ACR) at 3 months (ACR-3) post-transplant. We performed RNA sequencing on whole blood collected from 88 individuals at the time of 3-month surveillance biopsy to identify transcripts associated with ACR-3, developed a novel sequencing-based targeted expression assay, and validated this gene signature in an independent cohort. RESULTS: Study participants with ACR-3 had significantly higher risk than those without ACR-3 of subsequent clinical acute rejection at 12 and 24 months, faster decline in graft function, and decreased graft survival in adjusted Cox analysis. We identified a 17-gene signature in peripheral blood that accurately diagnosed ACR-3, and validated it using microarray expression profiles of blood samples from 65 transplant recipients in the GoCAR cohort and three public microarray datasets. In an independent cohort of 110 transplant recipients, tests of the targeted expression assay on the basis of the 17-gene set showed that it identified individuals at higher risk of ongoing acute rejection and future graft loss. CONCLUSIONS: Our targeted expression assay enabled noninvasive diagnosis of subclinical acute rejection and inflammation in the graft and may represent a useful tool to risk-stratify kidney transplant recipients.

Exome Sequencing and Prediction of Long-Term Kidney Allograft Function.

Current strategies to improve graft outcome following kidney transplantation consider information at the human leukocyte antigen (HLA) loci. Cell surface antigens, in addition to HLA, may serve as the stimuli as well as the targets for the anti-allograft immune response and influence long-term graft outcomes. We therefore performed exome sequencing of DNA from kidney graft recipients and their living donors and estimated all possible cell surface antigens mismatches for a given donor/recipient pair by computing the number of amino acid mismatches in trans-membrane proteins. We designated this tally as the allogenomics mismatch score (AMS). We examined the association between the AMS and post-transplant estimated glomerular filtration rate (eGFR) using mixed models, considering transplants from three independent cohorts (a total of 53 donor-recipient pairs, 106 exomes, and 239 eGFR measurements). We found that the AMS has a significant effect on eGFR (mixed model, effect size across the entire range of the score: -19.4 [-37.7, -1.1], P = 0.0042, χ2 = 8.1919, d.f. = 1) that is independent of the HLA-A, B, DR matching, donor age, and time post-transplantation. The AMS effect is consistent across the three independent cohorts studied and similar to the strong effect size of donor age. Taken together, these results show that the AMS, a novel tool to quantify amino acid mismatches in trans-membrane proteins in individual donor/recipient pair, is a strong, robust predictor of long-term graft function in kidney transplant recipients.

[Analysis of the deafness gene screening results from newborns in Shijiazhuang].

OBJECTIVE: To build information repository of the carrying rate of neonatal deafness gene in Shijiazhuang. METHOD: Blood samples were collected from the heel in 3-days neonates. Mutations of the deafness related genes were detected by the method of fluorescent PCR. Neonates received the detection of 6 mutation sites from 3 genes, including GJB2 (235delC, 299-300delAT), SLC26A4 (IVS7-2A> G, 2168A> G), mitochondrial DNA12S rRNA(1494C>T,1555A>G). RESULT: There were 384 neonates who carried mutations among 421 subjects and the carrying rate was 4.08%, 158 (1.68%) newborns carried heterozygous mutations and 1 (0.01%) case carried homogeneous mutation of GJB2 (235 delC), 55 (0.58%) neonates carried heterozygous mutations of GJB2 (299-300delAT); 133 (1.41%) neonates carried heterozygous mutations and 1 (0.01%) homogeneous of SLC26A4(IVS7-2A>G),19 (0.20%) newborns carried heterozygous mutations of SLC26A4 (2168A>G). The numbers of neonates who carried homogeneous and heterogeneous mutation of mitochondrial 12S rRNA gene were 14 and 3 with carring rates of 0.15% and 0.03%. Two newborns were found to carry more than one mutation. One carried 235delC, IVS7-2A>G and 1555A>G and another carried 235delC and IVS7-2A>G. CONCLUSION: The main mutational patterns were 235delC from GJB2 gene and IVS7-2A>G from SLC26A4 gene in Shijiazhuang newborns. The carrying rate information repository of neonatal deafness gene has been built preliminarily.

Analysis of the internal transcribed spacer 2 (ITS2) region of scuticociliates and related taxa (Ciliophora, Oligohymenophorea) to infer their evolution and phylogeny.

The ITS2 (ITS--internal transcribed spacer) region of the rDNA in 11 representative scuticociliates and two ambiguously related genera was analyzed. In common with other eukaryotes, the putative ITS2 folding pattern consists of a closed loop with four helices supported by minimum free energy and compensatory base changes (CBCs), although two of these helices are variable and sometimes absent. Three topologies were obtained on the basis of traditional primary sequence analysis, "string" strategy of secondary structure and analysis of the combined data. It was found that the secondary structure information could help to improve alignment and utilize appropriately phylogenetic strategies. The proposed phylogenies, though differing between sequence- and structure-based results, provide consistent support for high-level clades: the systematically questionable genera Dexiotrichides and Cardiostomatella always cluster together in a clade basal to the scuticociliates s.s., whereas Pleuronema branches from other uronematids at a deep level, and is hence a divergent taxon. Within the well-supported monophyletic philasterids, a sister relationship exists between Orchitophrya and Mesanophrys, while Uronema shows a close relationship with the group including Paranophrys and Parauronema. The positions of Metanophrys, Pseudocohnilembus and Anophryoides among the philasterids remain poorly resolved. Our findings firmly support the proposed evolutionary scenario inferred previously both from morphological and molecular data.