Whole-genome DNA Sequencing Research Articles

Patient-derived Enterococcus faecium with inflammatory genotypes promote colitis.

Dysbiosis of gut microbiota promotes colitis in ulcerative colitis (UC). Enterococcus faecium is an important constituent of dysbiotic microbiota. However, the mechanisms underlying E. faecium-induced colitis remain unclear. Overall, 23 E. faecium strains isolated from human feces and 3 commercial strains were inoculated into Il10-/- mice. Mouse colons were histologically evaluated and analyzed using real-time PCR analysis of cytokines. Genes in 26 E. faecium strains were identified by whole-genome shotgun sequencing of genomic DNA. The production of reactive oxygen species (ROS) from each strain was measured. An antioxidant, lipoic acid, was orally administered to the colitis mouse model. Inoculation of E. faecium derived from patients with UC resulted in colitis in Il10-/- mice. The genotypes of 26 strains were characterized by identifying 1893 known genes; clustering all the strains based on the genotypes showed two major groups-inflammatory and probiotic clusters. Additionally, linear discriminant analysis clarified that lipoic acid metabolism was a significantly abundant pathway in the probiotic cluster compared to the inflammatory cluster. Further, the production of ROS was greater in inflammatory than in probiotic strains. Administration of lipoic acid in E. faecium-inoculated mice ameliorated colitis. Enterococcus faecium strains in the inflammatory cluster promoted colitis with higher production of ROS than the strains in the probiotic cluster.

Genomic map of candidate human imprint control regions: the imprintome

ABSTRACT Imprinted genes – critical for growth, metabolism, and neuronal function – are expressed from one parental allele. Parent-of-origin-dependent CpG methylation regulates this expression at imprint control regions (ICRs). Since ICRs are established before tissue specification, these methylation marks are similar across cell types. Thus, they are attractive for investigating the developmental origins of adult diseases using accessible tissues, but remain unknown. We determined genome-wide candidate ICRs in humans by performing whole-genome bisulphite sequencing (WGBS) of DNA derived from the three germ layers and from gametes. We identified 1,488 hemi-methylated candidate ICRs, including 19 of 25 previously characterized ICRs (https://humanicr.org/). Gamete methylation approached 0% or 100% in 332 ICRs (178 paternally and 154 maternally methylated), supporting parent-of-origin-specific methylation, and 65% were in well-described CTCF-binding or DNaseI hypersensitive regions. This draft of the human imprintome will allow for the systematic determination of the role of early-acquired imprinting dysregulation in the pathogenesis of human diseases and developmental and behavioural disorders.

Genomic dynamics of brown trout populations released to a novel environment.

Population translocations occur for a variety of reasons, from displacement due to climate change to human‐induced transfers. Such actions have adverse effects on genetic variation and understanding their microevolutionary consequences requires monitoring. Here, we return to an experimental release of brown trout (Salmo trutta) in order to monitor the genomic effects of population translocations. In 1979, fish from each of two genetically (F ST = 0.16) and ecologically separate populations were simultaneously released, at one point in time, to a lake system previously void of brown trout. Here, whole‐genome sequencing of pooled DNA (Pool‐seq) is used to characterize diversity within and divergence between the introduced populations and fish inhabiting two lakes downstream of the release sites, sampled 30 years later (c. 5 generations). Present results suggest that while extensive hybridization has occurred, the two introduced populations are unequally represented in the lakes downstream of the release sites. One population, which is ecologically resident in its original habitat, mainly contributes to the lake closest to the release site. The other population, migratory in its natal habitat, is genetically more represented in the lake further downstream. Genomic regions putatively under directional selection in the new habitat are identified, where allele frequencies in both established populations are more similar to the introduced population stemming from a resident population than the migratory one. Results suggest that the microevolutionary consequences of population translocations, for example, hybridization and adaptation, can be rapid and that Pool‐seq can be used as an initial tool to monitor genome‐wide effects.

Sulphate repression of ssuD-dependent alkanesulphonate-sulphur assimilation in Escherichia coli.

Escherichia coli cells utilize alkanesulphonates including taurine as the sulphur source. We previously reported that when E. coli cells carrying a double deletion in tauD and cysN were inoculated into a taurine-containing minimal medium, they started to grow only after long-term incubation (Nishikawa et al. 2018, Microbiology 164: 1446-1456). We show here that cells that can induce ssuD-dependent alkanesulphonate-sulphur assimilation (SASSA) are essentially rare, but suppressors that can induce SASSA appear during long-term incubation. Mutant cells carrying ΔtauD and ΔcysN, ΔcysC or ΔcysH generated suppressor cells that can induce SASSA at a frequency of about 10-6 in a population. Whereas ΔtauD ΔcysN cells without prior SASSA did not express ssuD even when necessary, the cells with prior SASSA properly expressed ssuD. Whole-genome DNA sequencing of a clone isolated from ΔtauD ΔcysN cells with prior SASSA revealed that the influx of sulphate or thiosulphate may be related to the regulation of SASSA. To clarify whether sulphate or thiosulphate affects the induction of SASSA, the effect of mutations in sbp and cysP, which are responsible for sulphate and thiosulphate uptake with different preferences for substrates, was examined. Only the ΔtauD ΔcysN Δsbp mutant did not show repression of SASSA when no sulphate was added to the medium. When the concentration of the sulphate added was over 10 μM, the Δsbp mutant showed repression of SASSA. Therefore, it was considered that the influx of extracellular sulphate resulted in repression of SASSA.

Risk assessment with low-pass whole-genome sequencing of cell-free DNA before CD19 CAR T-cell therapy for large B-cell lymphoma

AbstractPatients with relapsed or refractory large B-cell lymphomas (rrLBCL) can achieve long-term remission after CD19 chimeric antigen receptor T-cell therapy (CART19). However, more than half of recipients will experience treatment failure. Thus, approaches are needed to identify high-risk patients who may benefit from alternative or consolidative therapy. We evaluated low-pass whole-genome sequencing (lpWGS) of cell-free DNA (cfDNA) before CART19 as a new approach for risk stratification. We performed lpWGS on pretreatment plasma samples from 122 patients at time of leukapheresis who received standard-of-care CART19 for rrLBCL to define DNA copy number alterations (CNAs). In multivariable selection, high focal CNA score (FCS) denoting genomic instability was the most significant pretreatment variable associated with inferior 3-month complete response rates (28% vs 56%, P = .0029), progression-free survival (PFS; P = .0007; hazard ratio, 2.11), and overall survival (OS; P = .0026; hazard ratio, 2.10). We identified 34 unique focal CNAs in 108 (89%) patients; of these, deletion 10q23.3 leading to loss of FAS death receptor was the most highly associated with poor outcomes, leading to inferior PFS (P < .0001; hazard ratio, 3.49) and OS (P = .0027; hazard ratio, 2.68). By combining FCS with traditional markers of increased tumor bulk (elevated lactate dehydrogenase and >1 extranodal site), we built a simple risk model that could reliably risk stratify patients. Thus, lpWGS of cfDNA is a minimally invasive assay that could rapidly identify high-risk patients and may guide patient selection for and targeted therapies to evaluate in future clinical trials.

A suspension technique for efficient large-scale cancer organoid culturing and perturbation screens

Organoid cell culture methodologies are enabling the generation of cell models from healthy and diseased tissue. Patient-derived cancer organoids that recapitulate the genetic and histopathological diversity of patient tumours are being systematically generated, providing an opportunity to investigate new cancer biology and therapeutic approaches. The use of organoid cultures for many applications, including genetic and chemical perturbation screens, is limited due to the technical demands and cost associated with their handling and propagation. Here we report and benchmark a suspension culture technique for cancer organoids which allows for the expansion of models to tens of millions of cells with increased efficiency in comparison to standard organoid culturing protocols. Using whole-genome DNA and RNA sequencing analyses, as well as medium-throughput drug sensitivity testing and genome-wide CRISPR-Cas9 screening, we demonstrate that cancer organoids grown as a suspension culture are genetically and phenotypically similar to their counterparts grown in standard conditions. This culture technique simplifies organoid cell culture and extends the range of organoid applications, including for routine use in large-scale perturbation screens.

Interpretable prediction of necrotizing enterocolitis from machine learning analysis of premature infant stool microbiota

BackgroundNecrotizing enterocolitis (NEC) is a common, potentially catastrophic intestinal disease among very low birthweight premature infants. Affecting up to 15% of neonates born weighing less than 1500 g, NEC causes sudden-onset, progressive intestinal inflammation and necrosis, which can lead to significant bowel loss, multi-organ injury, or death. No unifying cause of NEC has been identified, nor is there any reliable biomarker that indicates an individual patient’s risk of the disease. Without a way to predict NEC in advance, the current medical strategy involves close clinical monitoring in an effort to treat babies with NEC as quickly as possible before irrecoverable intestinal damage occurs. In this report, we describe a novel machine learning application for generating dynamic, individualized NEC risk scores based on intestinal microbiota data, which can be determined from sequencing bacterial DNA from otherwise discarded infant stool. A central insight that differentiates our work from past efforts was the recognition that disease prediction from stool microbiota represents a specific subtype of machine learning problem known as multiple instance learning (MIL).ResultsWe used a neural network-based MIL architecture, which we tested on independent datasets from two cohorts encompassing 3595 stool samples from 261 at-risk infants. Our report also introduces a new concept called the “growing bag” analysis, which applies MIL over time, allowing incorporation of past data into each new risk calculation. This approach allowed early, accurate NEC prediction, with a mean sensitivity of 86% and specificity of 90%. True-positive NEC predictions occurred an average of 8 days before disease onset. We also demonstrate that an attention-gated mechanism incorporated into our MIL algorithm permits interpretation of NEC risk, identifying several bacterial taxa that past work has associated with NEC, and potentially pointing the way toward new hypotheses about NEC pathogenesis. Our system is flexible, accepting microbiota data generated from targeted 16S or “shotgun” whole-genome DNA sequencing. It performs well in the setting of common, potentially confounding preterm neonatal clinical events such as perinatal cardiopulmonary depression, antibiotic administration, feeding disruptions, or transitions between breast feeding and formula.ConclusionsWe have developed and validated a robust MIL-based system for NEC prediction from harmlessly collected premature infant stool. While this system was developed for NEC prediction, our MIL approach may also be applicable to other diseases characterized by changes in the human microbiota.

Comprehensive Molecular Characterization Reveals Genomic and Transcriptomic Subtypes of Metastatic Urothelial Carcinoma

Recent molecular characterization of primary urothelial carcinoma (UC) may guide future clinical decision-making. For metastatic UC (mUC), a comprehensive molecular characterization is still lacking. We analyzed whole-genome DNA and RNA sequencing data for fresh-frozen metastatic tumor biopsies from 116 mUC patients who were scheduled for palliative systemic treatment within the context of a clinical trial (NCT01855477 and NCT02925234). Hierarchical clustering for mutational signatures revealed two major genomic subtypes: GenS1 (67%), which was APOBEC-driven; and GenS2 (24%), which had a high fraction of de novo mutational signatures related to reactive oxygen species and is putatively clock-like. Significantly mutated genes (SMGs) did not differ between the genomic subtypes. Transcriptomic analysis revealed five mUC subtypes: luminal-a and luminal-b (40%), stroma-rich (24%), basal/squamous (23%), and a nonspecified subtype (12%). These subtypes differed regarding expression of key genes, SMGs, oncogenic pathway activity, and immune cell infiltration. We integrated the genomic and transcriptomic data to propose potential therapeutic options by transcriptomic subtype and for individual patients. This in-depth analysis of a large cohort of patients with mUC may serve as a reference for subtype-oriented and patient-specific research on the etiology of mUC and for novel drug development. Patient summaryWe carried out an in-depth analysis of the molecular and genetic features of metastatic cancer involving the cells that line the urinary tract. We showed that this is a heterogeneous disease with different molecular subtypes and we identified possible targets for therapy for each subtype.

Detection of PKD1 and PKD2 Somatic Variants in Autosomal Dominant Polycystic Kidney Cyst Epithelial Cells by Whole-Genome Sequencing.

Autosomal dominant polycystic kidney disease (ADPKD) is a genetic disorder characterized by the development of multiple cysts in the kidneys. It is often caused by pathogenic mutations in PKD1 and PKD2 genes that encode polycystin proteins. Although the molecular mechanisms for cystogenesis are not established, concurrent inactivating germline and somatic mutations in PKD1 and PKD2 have been previously observed in renal tubular epithelium (RTE). To further investigate the cellular recessive mechanism of cystogenesis in RTE, we conducted whole-genome DNA sequencing analysis to identify germline variants and somatic alterations in RTE of 90 unique kidney cysts obtained during nephrectomy from 24 unrelated participants. Kidney cysts were overall genomically stable, with low burdens of somatic short mutations or large-scale structural alterations. Pathogenic somatic "second hit" alterations disrupting PKD1 or PKD2 were identified in 93% of the cysts. Of these, 77% of cysts acquired short mutations in PKD1 or PKD2 ; specifically, 60% resulted in protein truncations (nonsense, frameshift, or splice site) and 17% caused non-truncating mutations (missense, in-frame insertions, or deletions). Another 18% of cysts acquired somatic chromosomal loss of heterozygosity (LOH) events encompassing PKD1 or PKD2 ranging from 2.6 to 81.3 Mb. 14% of these cysts harbored copy number neutral LOH events, while the other 3% had hemizygous chromosomal deletions. LOH events frequently occurred at chromosomal fragile sites, or in regions comprising chromosome microdeletion diseases/syndromes. Almost all somatic "second hit" alterations occurred at the same germline mutated PKD1/2 gene. These findings further support a cellular recessive mechanism for cystogenesis in ADPKD primarily caused by inactivating germline and somatic variants of PKD1 or PKD2 genes in kidney cyst epithelium.

Genomic Instability of Circulating Tumor DNA as a Prognostic Marker for Pancreatic Cancer Survival: A Prospective Cohort Study.

Simple SummaryThis prospective cohort study showed that circulating tumor DNA-genomic instability (ctDNA-GI) I-scores, which was defined as the natural logarithm of the sum of LOESS-normalized Z-scores of sequenced reads in 1 Mb bins, are prognostic of the outcome of either localized or metastatic pancreatic adenocarcinoma. At baseline, 24.1% of patients had high genomic instability with I-score. Multivariable analyses demonstrated I-score was a significant factor for progression-free survival and overall survival.Genomic instability of circulating tumor DNA (ctDNA) as a prognostic biomarker has not been evaluated in pancreatic cancer. We investigated the role of the genomic instability index of ctDNA in pancreatic ductal adenocarcinoma (PDAC). We prospectively enrolled 315 patients newly diagnosed with resectable (n = 110), locally advanced (n = 78), and metastatic (n = 127) PDAC from March 2015 through January 2020. Low-depth whole-genome cell-free DNA sequencing identified genome-wide copy number alterations using instability score (I-score) to reflect genome-wide instability. Plasma cell-free and matched tumor tissue DNA from 15 patients with resectable pancreatic cancer was sequenced to assess the concordance of chromosomal copy number alteration profiles. Associations of I-score with clinical factors or survival were assessed. Seventy-six patients had high genomic instability with I-score > 7.3 in pre-treatment ctDNA; proportions of high I-score were 5.5%, 5.1%, and 52% in resectable, locally advanced, and metastatic stages, respectively. Correlation coefficients between Z-scores of plasma and tissue DNA at segment resolution were high (r2 = 0.82). Univariable analysis showed the association of I-score with progression-free survival in each stage. Multivariable analyses demonstrated that clinical stage-adjusted I-scores were significant factors for progression-free and overall survival. In these patients, ctDNA genomic I-scores provided prognostic information relevant to progression-free survival in each clinical stage.

Isolation of High-Molecular-Weight (HMW ) DNA from Fusarium oxysporum for Long-Read Sequencing.

In many instances Fusarium oxysporum genomes are complex and challenging to assemble mainly due to the increased number of repetitive elements and variable numbers of supernumerary chromosomes, which are primarily associated with pathogenicity. Hence, to obtain the accurate F. oxysporum genome assembly and high-resolution sequence information, protocols for versatile, reliable, and high recovery of high-quality DNA for diverse sequencing platforms are instrumental. Here, we describe two protocols for the isolation of DNA from isolates of F. oxysporum. One is a quick and easy method for high-throughput extraction of DNA to rapidly screen diverse isolates by marker-assisted PCR analysis. Another is to harvest high-quality and high-molecular-weight DNA for whole-genome sequencing. In addition, we also include a library preparation protocol optimized for the third-generation sequencing technology using the portable MinION device to obtain long-read sequences. These protocols can be potentially further applied for all Fusarium spp. and other fungal pathogens, including Verticillium.

First detection of porcine circovirus type 3 in Ukraine

To date, there is no information regarding the occurrence of porcine circovirus type 3 (PCV-3) in pigs in Ukraine. Aim. The aim of this work was to study the probable occurrence of the little-studied PCV-3 in pigs with different health status in Dnipropetrovsk, Donetsk, Kyiv, and Kharkiv regions of Ukraine. Methods. Blood, semen, liver, spleen, lung samples and nasal swabs of sows and boars of different ages and with different health status, belonging to farms from Dnipro, Donetsk, Kyiv, and Kharkiv regions of Ukraine, were used for the study. PCV-3 genomic material was detected by the standard polymerase chain reaction using specific primers, flanking a fragment of the rep gene of the virus with the length of 418 bp. To visualize the amplicons, horizontal gel electrophoresis was used and ethidium bromide staining after electrophoresis, followed by photographing the gels using Image Lab 5.2.1 software. Results. DNA of PCV-3 was found in two liver samples and four nasal swabs in two different farms, obtained from clinically healthy pigs, which suggests the possibility of the circulation of this infectious agent at the subclinical level of infection at the farm under investigation. No PCV-3 coinfection with the causative agents of porcine reproductive and respiratory syndrome (PRRS), Aujeszky’s disease, PCV-2, and mycoplasmas was found at this farm. Conclusions. Porcine circovirus type 3 (PCV-3) – (a little-studied causative agent of swine disease) was detected in 6 out of 61 samples, originating from two farms in the Kyiv and Kharkiv regions, obtained from clinically healthy animals) for the first time in Ukraine. This indicates possible circulation of the pathogen among pig farms in Ukraine and demonstrates the need to create and implement a target risk analysis, an extensive survey, as well as to develop control measures of the disease spreading (both organizational and technical preventive). Molecular genetic surveying and subsequent monitoring of PCV-3 among domestic and wild animals, which can cross the borders, will give a possibility to determine the risks of its spreading and related economic and epidemiological consequences. The whole-genome DNA sequencing of the detected virus isolates is planned to determine the relation of Ukrainian strains of the virus to other strains circulating in Europe and other parts of the world. Better understanding the risks, epidemiology and pathology, associated with this new virus for the Ukrainian pig breeding industry, will help to prevent and control its further spread and harmful effects.

A Novel Germline Compound Heterozygous Mutation of BRCA2 Gene Associated With Familial Peripheral Neuroblastic Tumors in Two Siblings.

ObjectivesTo investigate the genetic variants that are responsible for peripheral neuroblastic tumors (PNTs) oncogenesis in one family case.Materials and MethodsOne family was recruited, including the healthy parents, sister affected by neuroblastoma (NB), and brother who suffered from ganglioneuroma (GN). Whole-genome sequencing (WGS) of germline DNA from all the family members and RNA-seq of tumor RNA from the siblings were performed. Mutants were validated by Sanger sequencing and co-IP was performed to assess the impact of the mutant on chemosensitivity in the SH-SY5Y cell line.ResultsA novel compound heterozygous mutation of BRCA2 was locked as the cause of carcinogenesis. One allele was BRCA2-S871X (stop-gain) from the siblings’ mother, the other was BRCA2-N372H (missense) from their father. This novel compound heterozygous mutations of the BRCA2 gene associated with PNTs by disordering DNA damage and response (DDR) signal pathway. Moreover, chemosensitivity was reduced in the NB cell line due to the BRCA2-N372H mutant.ConclusionIn summary, these results revealed a novel germline compound heterozygous mutation of the BRCA2 gene associated with familial PNTs.

Molecular Characterization Reveals Genomic and Transcriptomic Subtypes of Metastatic Urothelial Carcinoma

Background: Metastatic urothelial carcinoma (mUC) is a lethal cancer for which few therapeutic options exist. Methods: To define the molecular landscape of mUC and identify targets for therapy, we performed whole-genome DNA sequencing on fresh-frozen metastatic tumor biopsies of 116 mUC patients, and mRNA sequencing on 90 matched biopsies. Results: Hierarchical clustering based on mutational signatures revealed two major genomic subtypes. The most prevalent subtype (67%) consisted of tumors with high APOBEC mutagenesis. Transcriptomic analysis revealed five mRNA-based subtypes: two luminal subtypes (40%), a stroma-rich (24%), basal/squamous (23%), and non-specified subtype (12%). The transcriptomic subtypes were different regarding driver gene alterations, gene amplifications, pathway activity, and immune cell infiltration. By integrating the genomic and transcriptomic data, potential therapeutic options per transcriptomic subtype and individual patient were proposed. Conclusions: This study serves as a reference for subtype-oriented and patient-specific research on the etiology of mUC, and for novel drug development.

Transient genomic instability drives tumorigenesis through accelerated clonal evolution

Abnormal numerical and structural chromosome content is frequently found in human cancer. To test the role of aneuploidy in tumor initiation and progression, we generated mice with random aneuploidies by transient induction of polo-like kinase 4 (Plk4), a master regulator of centrosome number. Short-term chromosome instability (CIN) from transient Plk4 induction resulted in formation of aggressive T-cell lymphomas in mice with heterozygous inactivation of one p53 allele and accelerated tumor development in the absence of p53. Transient CIN increased the frequency of lymphoma-initiating cells with a specific karyotype profile, including trisomy of chromosomes 4, 5, 14, and 15 occurring early in tumorigenesis. Tumor development in mice with chronic CIN induced by an independent mechanism (through inactivation of the spindle assembly checkpoint) gradually trended toward a similar karyotypic profile, as determined by single-cell whole-genome DNA sequencing. Overall, we show how transient CIN generates cells with random aneuploidies from which ones that acquire a karyotype with specific chromosome gains are sufficient to drive cancer formation, and that distinct CIN mechanisms can lead to similar karyotypic cancer-causing outcomes.

Replication Dynamics for Six Gram-Negative Bacterial Species during Bloodstream Infection.

ABSTRACTBloodstream infections (BSI) are a major public health burden due to high mortality rates and the cost of treatment. The impact of BSI is further compounded by a rise in antibiotic resistance among Gram-negative species associated with these infections. Escherichia coli, Serratia marcescens, Klebsiella pneumoniae, Enterobacter hormaechei, Citrobacter freundii, and Acinetobacter baumannii are all common causes of BSI, which can be recapitulated in a murine model. The objective of this study was to characterize infection kinetics and bacterial replication rates during bacteremia for these six pathogens to gain a better understanding of bacterial physiology during infection. Temporal observations of bacterial burdens of the tested species demonstrated varied abilities to establish colonization in the spleen, liver, or kidney. K. pneumoniae and S. marcescens expanded rapidly in the liver and kidney, respectively. Other organisms, such as C. freundii and E. hormaechei, were steadily cleared from all three target organs throughout the infection. In situ replication rates measured by whole-genome sequencing of bacterial DNA recovered from murine spleens demonstrated that each species was capable of sustained replication at 24 h postinfection, and several species demonstrated <60-min generation times. The relatively short generation times observed in the spleen were in contrast to an overall decrease in bacterial burden for some species, suggesting that the rate of immune-mediated clearance exceeded replication. Furthermore, bacterial generation times measured in the murine spleen approximated those measured during growth in human serum cultures. Together, these findings provide insight into the infection kinetics of six medically important species during bacteremia.

Comprehensive Somatic Copy Number Analysis Using Aqueous Humor Liquid Biopsy for Retinoblastoma

Simple SummaryAqueous humor (AH) liquid biopsy is an enriched source of cell-free circulating tumor-derived DNA for retinoblastoma (RB). The use of this AH liquid biopsy allows for genomic analysis of eyes in the absence of tumor tissue. Development of this platform was critical because direct tumor biopsy is prohibited in RB due to risk of extraocular tumor spread. In this retrospective study, we provide comprehensive, whole-genome analysis of the somatic copy number alterations (SCNAs) in 68 eyes of 64 RB patients. We show that the prevalence of specific SCNAs differ between eyes that required immediate enucleation (surgical removal) and eyes that were attempted to be saved but subsequently failed treatment, requiring secondary enucleation. Increases in chromosomal instability, or higher number of broad genomic alterations, predict higher risk clinical and biomarker features in these eyes. Prospective analyses are needed to further determine the clinical relevance and application of these findings.Aqueous humor (AH) liquid biopsy has been established as a surrogate tumor biopsy for retinoblastoma (RB). Previous AH studies have focused on highly recurrent RB somatic copy number alterations (SCNAs) including gain of 1q, 2p, 6p, and loss of 13q and 16q. In this retrospective study, we provide a comprehensive, whole-genome analysis of RB SCNAs and evaluate associated clinical features for 68 eyes of 64 RB patients from whom AH was obtained between December 2014 and October 2020. Shallow whole-genome sequencing of AH cell-free DNA was performed to assess for SCNAs. The prevalence of specific non-highly recurrent SCNAs, such as 20q gain and 8p loss, differed between primarily and secondarily enucleated eyes. Increases in chromosomal instability predict more advanced seeding morphology (p = 0.015); later age of diagnosis (p < 0.0001); greater odds of an endophytic tumor growth pattern (without retinal detachment; p = 0.047); tumor heights >10 mm (p = 0.09); and containing 6p gain, a biomarker of poor ocular prognosis (p = 0.004). The AH liquid biopsy platform is a high-yield method of whole-genome RB SCNA analysis, and SCNAs are associated with numerous clinical findings in RB eyes. Prospective analyses are encouraged to further elucidate the clinical relevance of specific SCNAs in RB.

Noninvasive Detection of Hepatocellular Carcinoma with Circulating Tumor DNA Features and α-Fetoprotein

Liver cancer is the fifth-most common cancer worldwide, with the third-highest rate of cancer-related mortality. Hepatocellular carcinoma (HCC) is the leading pathologic subtype, contributing 85% to 90% of cases of primary liver cancer. Most HCC patients are diagnosed at an advanced stage at which treatment is not curative. This study assessed the performance of a newly developed blood-based assay that utilizes genomic features and protein markers for the early detection of HCC. Two cancer-associated hallmarks, copy-number aberrations (CNA) and fragment size (FS), were characterized by shallow whole-genome sequencing of cell-free DNA and utilized to differentiate cancer patients from healthy subjects. As a clinically implemented biomarker of HCC, plasma α-fetoprotein (AFP) was also used with the genomic surrogates to optimize the detection of HCCs. The sensitivity of AFP ≥20.0 μg/L in detecting HCC was 57.9%. The combined genomic classifier CNA + FS via cell-free DNA shallow whole-genome sequencing identified nearly half of AFP-negative HCC patients (43.8%). By integrating CNA, FS as well as AFP (HCCseek), 75.0% sensitivity was achieved at 98.0% specificity, resulting in 92.6% accuracy, with 58.6% sensitivity in stage I HCC. The quantitative output of HCCseek was correlated with the severity of the disease (tumor size, stage, and recurrence-free survival). In summary, this study describes an efficient, noninvasive, and cost-effective method to detect HCC.

Evaluation of whole-genome DNA methylation sequencing library preparation protocols

BackgroundWith rapidly dropping sequencing cost, the popularity of whole-genome DNA methylation sequencing has been on the rise. Multiple library preparation protocols currently exist. We have performed 22 whole-genome DNA methylation sequencing experiments on snap frozen human samples, and extensively benchmarked common library preparation protocols for whole-genome DNA methylation sequencing, including three traditional bisulfite-based protocols and a new enzyme-based protocol. In addition, different input DNA quantities were compared for two kits compatible with a reduced starting quantity. In addition, we also present bioinformatic analysis pipelines for sequencing data from each of these library types.ResultsAn assortment of metrics were collected for each kit, including raw read statistics, library quality and uniformity metrics, cytosine retention, and CpG beta value consistency between technical replicates. Overall, the NEBNext Enzymatic Methyl-seq and Swift Accel-NGS Methyl-Seq kits performed quantitatively better than the other two protocols. In addition, the NEB and Swift kits performed well at low-input amounts, validating their utility in applications where DNA is the limiting factor.ResultsThe NEBNext Enzymatic Methyl-seq kit appeared to be the best option for whole-genome DNA methylation sequencing of high-quality DNA, closely followed by the Swift kit, which potentially works better for degraded samples. Further, a general bioinformatic pipeline is applicable across the four protocols, with the exception of extra trimming needed for the Swift Biosciences’s Accel-NGS Methyl-Seq protocol to remove the Adaptase sequence.

From Sampling to Sequencing: A Liquid Biopsy Pre-Analytic Workflow to Maximize Multi-Layer Genomic Information from a Single Tube.

Simple SummaryLiquid biopsies seek to isolate tumor derived genetic material that circulates in blood plasma or cerebrospinal fluid. The less-invasive character of liquid biopsies combined with the option for serial analyses bears enormous potential for treatment monitoring or surveillance. We aimed to establish robust sampling protocols and pre-analytical workflows to allow for site independent multi-layer liquid biopsy testing. For an optimal usage of precious material, we explored sample stabilization in various conservation tubes and describe a protocol for the parallel isolation of cell-free DNA and RNA. Quantification and quality control steps were optimized for minimal sample use with both high sensitivity and reproducibility. We provide detailed step-by-step information on how to i) choose the best-suited protocol and ii) implement this in the liquid biopsy workflow. We believe that our study has potential to increase comparability of liquid biopsy approaches to bring these one step closer to routine clinical application.Liquid biopsies hold great promise for the management of cancer. Reliable liquid biopsy data depend on stable and reproducible pre-analytical protocols that comply with quality measures, irrespective of the sampling and processing site. We established a workflow for plasma preservation, followed by processing, cell-free nucleic acid isolation, quantification, and enrichment of potentially tumor-derived cell-free DNA and RNA. Employing the same input material for a direct comparison of different kits and protocols allowed us to formulate unbiased recommendations for sample collection, storage, and processing. The presented workflow integrates the stabilization in Norgen, PAX, or Streck tubes and subsequent parallel isolation of cell-free DNA and RNA with NucleoSnap and NucleoSpin. Qubit, Bioanalyzer, and TapeStation quantification and quality control steps were optimized for minimal sample use and high sensitivity and reproducibility. We show the efficiency of the proposed workflow by successful droplet digital PCR amplification of both cell-free DNA and RNA and by detection of tumor-specific alterations in low-coverage whole-genome sequencing and DNA methylation profiling of plasma-derived cell-free DNA. For the first time, we demonstrated successful parallel extraction of cell-free DNA and RNA from plasma samples. This workflow paves the road towards multi-layer genomic analysis from one single liquid biopsy sample.