Evaluation of pre-analytical factors affecting plasma DNA analysis.

Pre-analytical factors can significantly affect circulating cell-free DNA (cfDNA) analysis. However, there are few robust methods to rapidly assess sample quality and the impact of pre-analytical processing. To address this gap and to evaluate effects of DNA extraction methods and blood collection tubes on cfDNA yield and fragment size, we developed a multiplexed droplet digital PCR (ddPCR) assay with 5 short and 4 long amplicons targeting single copy genomic loci. Using this assay, we compared 7 cfDNA extraction kits and found cfDNA yield and fragment size vary significantly. We also compared 3 blood collection protocols using plasma samples from 23 healthy volunteers (EDTA tubes processed within 1 hour and Cell-free DNA Blood Collection Tubes processed within 24 and 72 hours) and found no significant differences in cfDNA yield, fragment size and background noise between these protocols. In 219 clinical samples, cfDNA fragments were shorter in plasma samples processed immediately after venipuncture compared to archived samples, suggesting contribution of background DNA by lysed peripheral blood cells. In summary, we have described a multiplexed ddPCR assay to assess quality of cfDNA samples prior to downstream molecular analyses and we have evaluated potential sources of pre-analytical variation in cfDNA studies

Refined characterization of circulating tumor DNA through biological feature integration

AbstractCirculating tumor DNA (ctDNA) in blood plasma is present at very low concentrations compared to cell-free DNA (cfDNA) of non-tumor origin. To enhance ctDNA detection, recent studies have been focused on understanding the non-random fragmentation pattern of cfDNA. These studies have investigated fragment sizes, genomic position of fragment end points, and fragment end motifs. Although these features have been described and shown to be aberrant in cancer patients, there is a lack of understanding of how the individual and integrated analysis of these features enrich ctDNA fraction and enhance ctDNA detection. Using whole genome sequencing and copy number analysis of plasma samples from 5 high grade serious ovarian cancer patients, we observed that (1) ctDNA is enriched not only in fragments shorter than mono-nucleosomes (~ 167 bp), but also in those shorter than di-nucleosomes (~ 240–330 bp) (28–159% enrichment). (2) fragments that start and end at the border or within the nucleosome core are enriched in ctDNA (5–46% enrichment). (3) certain DNA motifs conserved in regions 10 bp up- and down- stream of fragment ends (i.e. cleavage sites) could be used to detect tumor-derived fragments (10–44% enrichment). We further show that the integrated analysis of these three features resulted in a higher enrichment of ctDNA when compared to using fragment size alone (additional 7–25% enrichment after fragment size selection). We believe these genome wide features, which are independent of genetic mutational changes, could allow new ways to analyze and interpret cfDNA data, as significant aberrations of these features from a healthy state could improve its utility as a diagnostic biomarker.</jats:p

The Genomic and Immune Landscapes of Lethal Metastatic Breast Cancer

TCR repertoire; Breast cancer; Clade mutationsRepertori TCR; Càncer de mama; Mutacions cladeRepertorio TCR; Cáncer de mama; Mutaciones cladoThe detailed molecular characterization of lethal cancers is a prerequisite to understanding resistance to therapy and escape from cancer immunoediting. We performed extensive multi-platform profiling of multi-regional metastases in autopsies from 10 patients with therapy-resistant breast cancer. The integrated genomic and immune landscapes show that metastases propagate and evolve as communities of clones, reveal their predicted neo-antigen landscapes, and show that they can accumulate HLA loss of heterozygosity (LOH). The data further identify variable tumor microenvironments and reveal, through analyses of T cell receptor repertoires, that adaptive immune responses appear to co-evolve with the metastatic genomes. These findings reveal in fine detail the landscapes of lethal metastatic breast cancer

Neuronally enriched microvesicle RNAs are differentially expressed in the serums of Parkinson’s patients

BackgroundCirculating small RNAs (smRNAs) originate from diverse tissues and organs. Previous studies investigating smRNAs as potential biomarkers for Parkinson’s disease (PD) have yielded inconsistent results. We investigated whether smRNA profiles from neuronally-enriched serum exosomes and microvesicles are altered in PD patients and discriminate PD subjects from controls.MethodsDemographic, clinical, and serum samples were obtained from 60 PD subjects and 40 age- and sex-matched controls. Exosomes and microvesicles were extracted and isolated using a validated neuronal membrane marker (CD171). Sequencing and bioinformatics analyses were used to identify differentially expressed smRNAs in PD and control samples. SmRNAs also were tested for association with clinical metrics. Logistic regression and random forest classification models evaluated the discriminative value of the smRNAs.ResultsIn serum CD171 enriched exosomes and microvesicles, a panel of 29 smRNAs was expressed differentially between PD and controls (false discovery rate (FDR) &lt; 0.05). Among the smRNAs, 23 were upregulated and 6 were downregulated in PD patients. Pathway analysis revealed links to cellular proliferation regulation and signaling. Least absolute shrinkage and selection operator adjusted for the multicollinearity of these smRNAs and association tests to clinical parameters via linear regression did not yield significant results. Univariate logistic regression models showed that four smRNAs achieved an AUC ≥ 0.74 to discriminate PD subjects from controls. The random forest model had an AUC of 0.942 for the 29 smRNA panel.ConclusionCD171-enriched exosomes and microvesicles contain the differential expression of smRNAs between PD and controls. Future studies are warranted to follow up on the findings and understand the scientific and clinical relevance

Refined characterization of circulating tumor DNA through biological feature integration

Circulating tumor DNA (ctDNA) in blood plasma is present at very low concentrations compared to cell-free DNA (cfDNA) of non-tumor origin. To enhance ctDNA detection, recent studies have been focused on understanding the non-random fragmentation pattern of cfDNA. These studies have investigated fragment sizes, genomic position of fragment end points, and fragment end motifs. Although these features have been described and shown to be aberrant in cancer patients, there is a lack of understanding of how the individual and integrated analysis of these features enrich ctDNA fraction and enhance ctDNA detection. Using whole genome sequencing and copy number analysis of plasma samples from 5 high grade serious ovarian cancer patients, we observed that (1) ctDNA is enriched not only in fragments shorter than mono-nucleosomes (~ 167 bp), but also in those shorter than di-nucleosomes (~ 240-330 bp) (28-159% enrichment). (2) fragments that start and end at the border or within the nucleosome core are enriched in ctDNA (5-46% enrichment). (3) certain DNA motifs conserved in regions 10 bp up- and down- stream of fragment ends (i.e. cleavage sites) could be used to detect tumor-derived fragments (10-44% enrichment). We further show that the integrated analysis of these three features resulted in a higher enrichment of ctDNA when compared to using fragment size alone (additional 7-25% enrichment after fragment size selection). We believe these genome wide features, which are independent of genetic mutational changes, could allow new ways to analyze and interpret cfDNA data, as significant aberrations of these features from a healthy state could improve its utility as a diagnostic biomarker

Evaluation of pre-analytical factors affecting plasma DNA analysis.

Pre-analytical factors can significantly affect circulating cell-free DNA (cfDNA) analysis. However, there are few robust methods to rapidly assess sample quality and the impact of pre-analytical processing. To address this gap and to evaluate effects of DNA extraction methods and blood collection tubes on cfDNA yield and fragment size, we developed a multiplexed droplet digital PCR (ddPCR) assay with 5 short and 4 long amplicons targeting single copy genomic loci. Using this assay, we compared 7 cfDNA extraction kits and found cfDNA yield and fragment size vary significantly. We also compared 3 blood collection protocols using plasma samples from 23 healthy volunteers (EDTA tubes processed within 1 hour and Cell-free DNA Blood Collection Tubes processed within 24 and 72 hours) and found no significant differences in cfDNA yield, fragment size and background noise between these protocols. In 219 clinical samples, cfDNA fragments were shorter in plasma samples processed immediately after venipuncture compared to archived samples, suggesting contribution of background DNA by lysed peripheral blood cells. In summary, we have described a multiplexed ddPCR assay to assess quality of cfDNA samples prior to downstream molecular analyses and we have evaluated potential sources of pre-analytical variation in cfDNA studies

Genetic imputation of kidney transcriptome, proteome and multi-omics illuminates new blood pressure and hypertension targets

Genetic mechanisms of blood pressure (BP) regulation remain poorly defined. Using kidney-specific epigenomic annotations and 3D genome information we generated and validated gene expression prediction models for the purpose of transcriptome-wide association studies in 700 human kidneys. We identified 889 kidney genes associated with BP of which 399 were prioritised as contributors to BP regulation. Imputation of kidney proteome and microRNAome uncovered 97 renal proteins and 11 miRNAs associated with BP. Integration with plasma proteomics and metabolomics illuminated circulating levels of myo-inositol, 4-guanidinobutanoate and angiotensinogen as downstream effectors of several kidney BP genes (SLC5A11, AGMAT, AGT, respectively). We showed that genetically determined reduction in renal expression may mimic the effects of rare loss-of-function variants on kidney mRNA/protein and lead to an increase in BP (e.g., ENPEP). We demonstrated a strong correlation (r = 0.81) in expression of protein-coding genes between cells harvested from urine and the kidney highlighting a diagnostic potential of urinary cell transcriptomics. We uncovered adenylyl cyclase activators as a repurposing opportunity for hypertension and illustrated examples of BP-elevating effects of anticancer drugs (e.g. tubulin polymerisation inhibitors). Collectively, our studies provide new biological insights into genetic regulation of BP with potential to drive clinical translation in hypertension.</p

Genome-wide analysis of aberrant position and sequence of plasma DNA fragment ends in patients with cancer.

Genome-wide fragmentation patterns in cell-free DNA (cfDNA) in plasma are strongly influenced by cellular origin due to variation in chromatin accessibility across cell types. Such differences between healthy and cancer cells provide the opportunity for development of novel cancer diagnostics. Here, we investigated whether analysis of cfDNA fragment end positions and their surrounding DNA sequences reveals the presence of tumor-derived DNA in blood. We performed genome-wide analysis of cfDNA from 521 samples and analyzed sequencing data from an additional 2147 samples, including healthy individuals and patients with 11 different cancer types. We developed a metric based on genome-wide differences in fragment positioning, weighted by fragment length and GC content [information-weighted fraction of aberrant fragments (iwFAF)]. We observed that iwFAF strongly correlated with tumor fraction, was higher for DNA fragments carrying somatic mutations, and was higher within genomic regions affected by copy number amplifications. We also calculated sample-level means of nucleotide frequencies observed at genomic positions spanning fragment ends. Using a combination of iwFAF and nine nucleotide frequencies from three positions surrounding fragment ends, we developed a machine learning model to differentiate healthy individuals from patients with cancer. We observed an area under the receiver operative characteristic curve (AUC) of 0.91 for detection of cancer at any stage and an AUC of 0.87 for detection of stage I cancer. Our findings remained robust with as few as 1 million fragments analyzed per sample, demonstrating that analysis of fragment ends can become a cost-effective and accessible approach for cancer detection and monitoring

Summary statistics of UK Biobank blood pressure genome-wide association studies (GWAS) using 337,422 unrelated white European individuals

Three blood pressure traits were analysed: systolic blood pressure (SBP), diastolic blood pressure (DBP) and pulse pressure (PP; the difference between SBP and DBP). Mean SBP and DBP values from automated values were calculated. After calculating blood pressure values, SBP and DBP were adjusted for medication use by adding 15 and 10 mm Hg to their values, respectively, for individuals reported to be taking blood pressure–lowering medication.For the UK Biobank genome-wide association studies (GWAS), we performed linear mixed model (LMM) association testing under an additive genetic model of the three continuous, medication-adjusted blood pressure traits (SBP, DBP, PP) for all measured and imputed genetic variants (Data Field-22828) with minor allele frequency (MAF) &gt;=1% and imputation score&gt;=0.3 in dosage format using the BOLT-LMM (v2.4.1) software. Covariates were age, age2, sex, BMI, genotyping array and 10PCs. Genomic inflation was not applied to the GWAS summary statistics.Sample QC was described below:We included up to 337,422 individuals from UK Biobank for the purpose of this project. We followed UK Biobank sample-based quality control criteria (Nature 2018;562:203-209); excluded were samples/individuals based on the following criteria: (i) outliers in heterozygosity and missingness, (ii) self-reported gender not consistent with genetic data inferred gender (ii) sample call rate (computed using probesets internal to Affymetrix