Jagged Ends of Urinary Cell-Free DNA: Characterization and Feasibility Assessment in Bladder Cancer Detection.

BACKGROUND: Double-stranded DNA in plasma is known to carry single-stranded ends, called jagged ends. Plasma DNA jagged ends are biomarkers for pathophysiologic states such as pregnancy and cancer. It remains unknown whether urinary cell-free DNA (cfDNA) molecules have jagged ends. METHODS: Jagged ends of cfDNA were detected by incorporating unmethylated cytosines during a DNA end-repair process, followed by bisulfite sequencing. Incorporation of unmethylated cytosines during the repair of the jagged ends lowered the apparent methylation levels measured by bisulfite sequencing and were used to calculate a jagged end index. This approach is called jagged end analysis by sequencing. RESULTS: The jagged end index of urinary cfDNA was higher than that of plasma DNA. The jagged end index profile of plasma DNA displayed several strongly oscillating major peaks at intervals of approximately 165 bp (i.e., nucleosome size) and weakly oscillating minor peaks with periodicities of approximately 10 bp. In contrast, the urinary DNA jagged end index profile showed weakly oscillating major peaks but strongly oscillating minor peaks. The jagged end index was generally higher in nucleosomal linker DNA regions. Patients with bladder cancer (n = 46) had lower jagged end indexed of urinary DNA than participants without bladder cancer (n = 39). The area under the curve for differentiating between patients with and without bladder cancer was 0.83. CONCLUSIONS: Jagged ends represent a property of urinary cfDNA. The generation of jagged ends might be related to nucleosomal structures, with enrichment in linker DNA regions. Jagged ends of urinary DNA could potentially serve as a new biomarker for bladder cancer detection.

Mendelian randomization analyses suggest a role for cholesterol in the development of endometrial cancer.

Blood lipids have been associated with the development of a range of cancers, including breast, lung and colorectal cancer. For endometrial cancer, observational studies have reported inconsistent associations between blood lipids and cancer risk. To reduce biases from unmeasured confounding, we performed a bidirectional, two-sample Mendelian randomization analysis to investigate the relationship between levels of three blood lipids (low-density lipoprotein [LDL] and high-density lipoprotein [HDL] cholesterol, and triglycerides) and endometrial cancer risk. Genetic variants associated with each of these blood lipid levels (P < 5 × 10-8 ) were identified as instrumental variables, and assessed using genome-wide association study data from the Endometrial Cancer Association Consortium (12 906 cases and 108 979 controls) and the Global Lipids Genetic Consortium (n = 188 578). Mendelian randomization analyses found genetically raised LDL cholesterol levels to be associated with lower risks of endometrial cancer of all histologies combined, and of endometrioid and non-endometrioid subtypes. Conversely, higher genetically predicted HDL cholesterol levels were associated with increased risk of non-endometrioid endometrial cancer. After accounting for the potential confounding role of obesity (as measured by genetic variants associated with body mass index), the association between genetically predicted increased LDL cholesterol levels and lower endometrial cancer risk remained significant, especially for non-endometrioid endometrial cancer. There was no evidence to support a role for triglycerides in endometrial cancer development. Our study supports a role for LDL and HDL cholesterol in the development of non-endometrioid endometrial cancer. Further studies are required to understand the mechanisms underlying these findings.

Comprehensive characterization and resolution of discrepant spectrophotometric bilirubin results in patients on eltrombopag therapy.

Background Eltrombopag is a thrombopoietin receptor agonist used for the treatment of thrombocytopenic conditions. It can cause pH-dependent discoloration of plasma/serum. Eltrombopag is potentially hepatotoxic. It can affect the assessment of hyperbilirubinemia because of its (i) absorbance at ~450 nm (bilirubin), (ii) absorbance at ~550 nm (diazo-bilirubin) and (iii) it can cause yellowish discoloration of the eyes at normal circulating bilirubin levels. Methods We collected 66 samples from patients on a range of eltrombopag dosages up to 150 mg daily. Bilirubin was measured using multiple routine spectrophotometric analyzers, the Doumas reference method and high-performance liquid chromatography (HPLC). Plasma/serum eltrombopag concentrations were determined using liquid chromatography tandem mass spectrometry (LC-MS/MS). Spike-in and admixture experiments delineated the effects of eltrombopag and its metabolites. Results Forty-nine of 52 samples from patients on ≥50 mg daily eltrombopag therapy showed significantly discrepant inter-analyzer total bilirubin results, a difference up to 64 µmol/L (3.7 mg/dL). In one sample, total bilirubin varied from 8 to 65 µmol/L (0.4-3.8 mg/dL) by different routine analyzers, with direct bilirubin ≤4 µmol/L (0.2 mg/dL). There was a positive correlation between total bilirubin difference and plasma eltrombopag concentration (r = 0.679), and spike-in experiments demonstrated that Beckman AU and Doumas reference methods were susceptible to positive interference. HPLC can quantify bilirubin after separating eltrombopag, and results suggest different analyzers are affected to varying degrees by eltrombopag and its metabolites. Conclusions Eltrombopag and its metabolites can cause positive interference to the spectrophotometric measurements of total bilirubin. Accurate measurements of total bilirubin may improve our understanding of the prevalence of hyperbilirubinemia in patients on eltrombopag therapy.

Noninvasive Detection of Bladder Cancer by Shallow-Depth Genome-Wide Bisulfite Sequencing of Urinary Cell-Free DNA for Methylation and Copy Number Profiling.

BACKGROUND: The current diagnosis and monitoring of bladder cancer are heavily reliant on cystoscopy, an invasive and costly procedure. Previous efforts in urine-based detection of bladder cancer focused on targeted approaches that are predicated on the tumor expressing specific aberrations. We aimed to noninvasively detect bladder cancer by the genome-wide assessment of methylomic and copy number aberrations (CNAs). We also investigated the size of tumor cell-free (cf)DNA fragments. METHODS: Shallow-depth paired-end genome-wide bisulfite sequencing of urinary cfDNA was done for 46 bladder cancer patients and 39 cancer-free controls with hematuria. We assessed (a) proportional contribution from different tissues by methylation deconvolution, (b) global hypomethylation, (c) CNA, and (d) cfDNA size profile. RESULTS: Methylomic and copy number approaches were synergistically combined to detect bladder cancer with a sensitivity of 93.5% (84.2% for low-grade nonmuscle-invasive disease) and a specificity of 95.8%. The prevalence of methylomic and CNAs reflected disease stage and tumor size. Sampling over multiple time points could assess residual disease and changes in tumor load. Muscle-invasive bladder cancer was associated with a higher proportion of long cfDNA, as well as longer cfDNA fragments originating from genomic regions enriched for tumor DNA. CONCLUSIONS: Bladder cancer can be detected noninvasively in urinary cfDNA by methylomic and copy number analysis without previous knowledge or assumptions of specific aberrations. Such analysis could be used as a liquid biopsy to aid diagnosis and for potential longitudinal monitoring of tumor load. Further understanding of the differential size and fragmentation of cfDNA could improve the detection of bladder cancer.

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.

A Targeted Gene Panel That Covers Coding, Non-coding and Short Tandem Repeat Regions Improves the Diagnosis of Patients With Neurodegenerative Diseases.

Genetic testing for neurodegenerative diseases (NDs) is highly challenging because of genetic heterogeneity and overlapping manifestations. Targeted-gene panels (TGPs), coupled with next-generation sequencing (NGS), can facilitate the profiling of a large repertoire of ND-related genes. Due to the technical limitations inherent in NGS and TGPs, short tandem repeat (STR) variations are often ignored. However, STR expansions are known to cause such NDs as Huntington's disease and spinocerebellar ataxias type 3 (SCA3). Here, we studied the clinical utility of a custom-made TGP that targets 199 NDs and 311 ND-associated genes on 118 undiagnosed patients. At least one known or likely pathogenic variation was found in 54 patients; 27 patients demonstrated clinical profiles that matched the variants; and 16 patients whose original diagnosis were refined. A high concordance of variant calling were observed when comparing the results from TGP and whole-exome sequencing of four patients. Our in-house STR detection algorithm has reached a specificity of 0.88 and a sensitivity of 0.82 in our SCA3 cohort. This study also uncovered a trove of novel and recurrent variants that may enrich the repertoire of ND-related genetic markers. We propose that a combined comprehensive TGPs-bioinformatics pipeline can improve the clinical diagnosis of NDs.

Preferred end coordinates and somatic variants as signatures of circulating tumor DNA associated with hepatocellular carcinoma.

Circulating tumor-derived cell-free DNA (ctDNA) analysis offers an attractive noninvasive means for detection and monitoring of cancers. Evidence for the presence of cancer is dependent on the ability to detect features in the peripheral circulation that are deemed as cancer-associated. We explored approaches to improve the chance of detecting the presence of cancer based on sequence information present on ctDNA molecules. We developed an approach to detect the total pool of somatic mutations. We then investigated if there existed a class of ctDNA signature in the form of preferred plasma DNA end coordinates. Cell-free DNA fragmentation is a nonrandom process. Using plasma samples obtained from liver transplant recipients, we showed that liver contributed cell-free DNA molecules ended more frequently at certain genomic coordinates than the nonliver-derived molecules. The abundance of plasma DNA molecules with these liver-associated ends correlated with the liver DNA fractions in the plasma samples. Studying the DNA end characteristics in plasma of patients with hepatocellular carcinoma and chronic hepatitis B, we showed that there were millions of tumor-associated plasma DNA end coordinates in the genome. Abundance of plasma DNA molecules with tumor-associated DNA ends correlated with the tumor DNA fractions even in plasma samples of hepatocellular carcinoma patients that were subjected to shallow-depth sequencing analysis. Plasma DNA end coordinates may therefore serve as hallmarks of ctDNA that could be sampled readily and, hence, may improve the cost-effectiveness of liquid biopsy assessment.

Identification of nine new susceptibility loci for endometrial cancer.

Endometrial cancer is the most commonly diagnosed cancer of the female reproductive tract in developed countries. Through genome-wide association studies (GWAS), we have previously identified eight risk loci for endometrial cancer. Here, we present an expanded meta-analysis of 12,906 endometrial cancer cases and 108,979 controls (including new genotype data for 5624 cases) and identify nine novel genome-wide significant loci, including a locus on 12q24.12 previously identified by meta-GWAS of endometrial and colorectal cancer. At five loci, expression quantitative trait locus (eQTL) analyses identify candidate causal genes; risk alleles at two of these loci associate with decreased expression of genes, which encode negative regulators of oncogenic signal transduction proteins (SH2B3 (12q24.12) and NF1 (17q11.2)). In summary, this study has doubled the number of known endometrial cancer risk loci and revealed candidate causal genes for future study.

Genetic overlap between endometriosis and endometrial cancer: evidence from cross-disease genetic correlation and GWAS meta-analyses.

Epidemiological, biological, and molecular data suggest links between endometriosis and endometrial cancer, with recent epidemiological studies providing evidence for an association between a previous diagnosis of endometriosis and risk of endometrial cancer. We used genetic data as an alternative approach to investigate shared biological etiology of these two diseases. Genetic correlation analysis of summary level statistics from genomewide association studies (GWAS) using LD Score regression revealed moderate but significant genetic correlation (rg  = 0.23, P = 9.3 × 10-3 ), and SNP effect concordance analysis provided evidence for significant SNP pleiotropy (P = 6.0 × 10-3 ) and concordance in effect direction (P = 2.0 × 10-3 ) between the two diseases. Cross-disease GWAS meta-analysis highlighted 13 distinct loci associated at P ≤ 10-5 with both endometriosis and endometrial cancer, with one locus (SNP rs2475335) located within PTPRD associated at a genomewide significant level (P = 4.9 × 10-8 , OR = 1.11, 95% CI = 1.07-1.15). PTPRD acts in the STAT3 pathway, which has been implicated in both endometriosis and endometrial cancer. This study demonstrates the value of cross-disease genetic analysis to support epidemiological observations and to identify biological pathways of relevance to multiple diseases.

Genomewide bisulfite sequencing reveals the origin and time-dependent fragmentation of urinary cfDNA.

Urinary cell-free (cf) DNA holds great potential as a completely noninvasive form of liquid biopsy. Knowledge of the composition of cfDNA by tissue of origin is useful for guiding its clinical uses. We conducted a global survey of urinary cfDNA composition using genomewide bisulfite sequencing. While previous studies focused on detecting cfDNA from a single source at a time, genomewide tissue specific methylation signatures allow us to simultaneously deduce the proportional contribution from each contributing tissue. The proportional contributions derived from methylation deconvolution are highly correlated with those calculated using allograft-derived donor-specific genetic markers in the urine of hematopoetic stem cell and renal transplant recipients. We found a large variation of proportional contributions from different tissues. We then assessed if cfDNA undergoes time-dependent fragmentation in urine by conducting in vitro incubation experiments. In vitro incubation at 37°C showed that urinary cfDNA concentration decreased under first order kinetics with a half-life of 2.6 to 5.1h. This is reflected in parallel by a decrease in the proportion of long fragments and increase in amplitude of 10bp periodicity seen in the cfDNA size profile. This global survey of urinary cfDNA has deepened our understanding of the composition, degradation and variation of cfDNA in the urinary tract and has laid a foundation for the use of genomewide urinary cfDNA sequencing as a molecular diagnostics tool.

A Common Variant at the 14q32 Endometrial Cancer Risk Locus Activates AKT1 through YY1 Binding.

A recent meta-analysis of multiple genome-wide association and follow-up endometrial cancer case-control datasets identified a novel genetic risk locus for this disease at chromosome 14q32.33. To prioritize the functional SNP(s) and target gene(s) at this locus, we employed an in silico fine-mapping approach using genotyped and imputed SNP data for 6,608 endometrial cancer cases and 37,925 controls of European ancestry. Association and functional analyses provide evidence that the best candidate causal SNP is rs2494737. Multiple experimental analyses show that SNP rs2494737 maps to a silencer element located within AKT1, a member of the PI3K/AKT/MTOR intracellular signaling pathway activated in endometrial tumors. The rs2494737 risk A allele creates a YY1 transcription factor-binding site and abrogates the silencer activity in luciferase assays, an effect mimicked by transfection of YY1 siRNA. Our findings suggest YY1 is a positive regulator of AKT1, mediating the stimulatory effects of rs2494737 increasing endometrial cancer risk. Identification of an endometrial cancer risk allele within a member of the PI3K/AKT signaling pathway, more commonly activated in tumors by somatic alterations, raises the possibility that well tolerated inhibitors targeting this pathway could be candidates for evaluation as chemopreventive agents in individuals at high risk of developing endometrial cancer.

Meta-analysis of genome-wide association studies identifies common susceptibility polymorphisms for colorectal and endometrial cancer near SH2B3 and TSHZ1.

High-risk mutations in several genes predispose to both colorectal cancer (CRC) and endometrial cancer (EC). We therefore hypothesised that some lower-risk genetic variants might also predispose to both CRC and EC. Using CRC and EC genome-wide association series, totalling 13,265 cancer cases and 40,245 controls, we found that the protective allele [G] at one previously-identified CRC polymorphism, rs2736100 near TERT, was associated with EC risk (odds ratio (OR) = 1.08, P = 0.000167); this polymorphism influences the risk of several other cancers. A further CRC polymorphism near TERC also showed evidence of association with EC (OR = 0.92; P = 0.03). Overall, however, there was no good evidence that the set of CRC polymorphisms was associated with EC risk, and neither of two previously-reported EC polymorphisms was associated with CRC risk. A combined analysis revealed one genome-wide significant polymorphism, rs3184504, on chromosome 12q24 (OR = 1.10, P = 7.23 × 10(-9)) with shared effects on CRC and EC risk. This polymorphism, a missense variant in the gene SH2B3, is also associated with haematological and autoimmune disorders, suggesting that it influences cancer risk through the immune response. Another polymorphism, rs12970291 near gene TSHZ1, was associated with both CRC and EC (OR = 1.26, P = 4.82 × 10(-8)), with the alleles showing opposite effects on the risks of the two cancers.

Common colorectal cancer risk alleles contribute to the multiple colorectal adenoma phenotype, but do not influence colonic polyposis in FAP.

The presence of multiple (5-100) colorectal adenomas suggests an inherited predisposition, but the genetic aetiology of this phenotype is undetermined if patients test negative for Mendelian polyposis syndromes such as familial adenomatous polyposis (FAP) and MUTYH-associated polyposis (MAP). We investigated whether 18 common colorectal cancer (CRC) predisposition single-nucleotide polymorphisms (SNPs) could help to explain some cases with multiple adenomas who phenocopied FAP or MAP, but had no pathogenic APC or MUTYH variant. No multiple adenoma case had an outlying number of CRC SNP risk alleles, but multiple adenoma patients did have a significantly higher number of risk alleles than population controls (P=5.7 × 10(-7)). The association was stronger in those with ≥10 adenomas. The CRC SNPs accounted for 4.3% of the variation in multiple adenoma risk, with three SNPs (rs6983267, rs10795668, rs3802842) explaining 3.0% of the variation. In FAP patients, the CRC risk score did not differ significantly from the controls, as we expected given the overwhelming effect of pathogenic germline APC variants on the phenotype of these cases. More unexpectedly, we found no evidence that the CRC SNPs act as modifier genes for the number of colorectal adenomas in FAP patients. In conclusion, common colorectal tumour risk alleles contribute to the development of multiple adenomas in patients without pathogenic germline APC or MUTYH variants. This phenotype may have 'polygenic' or monogenic origins. The risk of CRC in relatives of multiple adenoma cases is probably much lower for cases with polygenic disease, and this should be taken into account when counselling such patients.