Deep whole-genome sequencing of 3 cancer cell lines on 2 sequencing platforms

Rare genetic variation in fibronectin 1 (FN1) protects against APOEε4 in Alzheimer's disease

The risk of developing Alzheimer's disease (AD) significantly increases in individuals carrying the APOEε4 allele. Elderly cognitively healthy individuals with APOEε4 also exist, suggesting the presence of cellular mechanisms that counteract the pathological effects of APOEε4; however, these mechanisms are unknown. We hypothesized that APOEε4 carriers without dementia might carry genetic variations that could protect them from developing APOEε4-mediated AD pathology. To test this, we leveraged whole-genome sequencing (WGS) data in the National Institute on Aging Alzheimer's Disease Family Based Study (NIA-AD FBS), Washington Heights/Inwood Columbia Aging Project (WHICAP), and Estudio Familiar de Influencia Genetica en Alzheimer (EFIGA) cohorts and identified potentially protective variants segregating exclusively among unaffected APOEε4 carriers. In homozygous unaffected carriers above 70 years old, we identified 510 rare coding variants. Pathway analysis of the genes harboring these variants showed significant enrichment in extracellular matrix (ECM)-related processes, suggesting protective effects of functional modifications in ECM proteins. We prioritized two genes that were highly represented in the ECM-related gene ontology terms, (FN1) and collagen type VI alpha 2 chain (COL6A2) and are known to be expressed at the blood-brain barrier (BBB), for postmortem validation and in vivo functional studies. An independent analysis in a large cohort of 7185 APOEε4 homozygous carriers found that rs140926439 variant in FN1 was protective of AD (OR = 0.29; 95% CI [0.11, 0.78], P = 0.014) and delayed age at onset of disease by 3.37 years (95% CI [0.42, 6.32], P = 0.025). The FN1 and COL6A2 protein levels were increased at the BBB in APOEε4 carriers with AD. Brain expression of cognitively unaffected homozygous APOEε4 carriers had significantly lower FN1 deposition and less reactive gliosis compared to homozygous APOEε4 carriers with AD, suggesting that FN1 might be a downstream driver of APOEε4-mediated AD-related pathology and cognitive decline. To validate our findings, we used zebrafish models with loss-of-function (LOF) mutations in fn1b-the ortholog for human FN1. We found that fibronectin LOF reduced gliosis, enhanced gliovascular remodeling, and potentiated the microglial response, suggesting that pathological accumulation of FN1 could impair toxic protein clearance, which is ameliorated with FN1 LOF. Our study suggests that vascular deposition of FN1 is related to the pathogenicity of APOEε4, and LOF variants in FN1 may reduce APOEε4-related AD risk, providing novel clues to potential therapeutic interventions targeting the ECM to mitigate AD risk.

Childhood cancer mutagenesis caused by transposase-derived PGBD5

Genomic rearrangements are a hallmark of most childhood tumors, including medulloblastoma, one of the most common brain tumors in children, but their causes remain largely unknown. Here, we show that PiggyBac transposable element derived 5 (Pgbd5) promotes tumor development in multiple developmentally accurate mouse models of Sonic Hedgehog (SHH) medulloblastoma. Most Pgbd5-deficient mice do not develop tumors, while maintaining normal cerebellar development. Ectopic activation of SHH signaling is sufficient to enforce cerebellar granule cell progenitor-like cell states, which exhibit Pgbd5-dependent expression of distinct DNA repair and neurodevelopmental factors. Mouse medulloblastomas expressing Pgbd5 have increased numbers of somatic structural DNA rearrangements, some of which carry PGBD5-specific sequences at their breakpoints. Similar sequence breakpoints recurrently affect somatic DNA rearrangements of known tumor suppressors and oncogenes in medulloblastomas in 329 children. This identifies PGBD5 as a medulloblastoma mutator and provides a genetic mechanism for the generation of oncogenic DNA rearrangements in childhood cancer.

Detecting Somatic Mutations Without Matched Normal Samples Using Long Reads

DNA sequencing of tumours to identify somatic mutations has become a critical tool to guide the type of treatment given to cancer patients. The gold standard for mutation calling is comparing sequencing data from the tumour to a matched normal sample to avoid mis-classifying inherited SNPs as mutations. This procedure works extremely well, but in certain situations only a tumour sample is available. While approaches have been developed to find mutations without a matched normal, they have limited accuracy or require specific types of input data (e.g. ultra-deep sequencing). Here we explore the application of single molecule long read sequencing to calling somatic mutations without matched normal samples. We develop a simple theoretical framework to show how haplotype phasing is an important source of information for determining whether a variant is a somatic mutation. We then use simulations to assess the range of experimental parameters (tumour purity, sequencing depth) where this approach is effective. These ideas are developed into a prototype somatic mutation caller, smrest, and its use is demonstrated on two highly mutated cancer cell lines. Finally, we argue that this approach has potential to measure clinically important biomarkers that are based on the genome-wide distribution of mutations: tumour mutation burden and mutation signatures.

Rare genetic variation in Fibronectin 1 ( FN1 ) protects against APOEe4 in Alzheimer's disease

The risk of developing Alzheimer's disease (AD) significantly increases in individuals carrying the APOEε4 allele. Elderly cognitively healthy individuals with APOEε4 also exist, suggesting the presence of cellular mechanisms that counteract the pathological effects of APOEε4 ; however, these mechanisms are unknown. We hypothesized that APOEε4 carriers without dementia might carry genetic variations that could protect them from developing APOEε4- mediated AD pathology. To test this, we leveraged whole genome sequencing (WGS) data in National Institute on Aging Alzheimer's Disease Family Based Study (NIA-AD FBS), Washington Heights/Inwood Columbia Aging Project (WHICAP), and Estudio Familiar de Influencia Genetica en Alzheimer (EFIGA) cohorts and identified potentially protective variants segregating exclusively among unaffected APOEε4 carriers. In homozygous unaffected carriers above 70 years old, we identified 510 rare coding variants. Pathway analysis of the genes harboring these variants showed significant enrichment in extracellular matrix (ECM)-related processes, suggesting protective effects of functional modifications in ECM proteins. We prioritized two genes that were highly represented in the ECM-related gene ontology terms, (FN1) and collagen type VI alpha 2 chain ( COL6A2 ) and are known to be expressed at the blood-brain barrier (BBB), for postmortem validation and in vivo functional studies. The FN1 and COL6A2 protein levels were increased at the BBB in APOEε4 carriers with AD. Brain expression of cognitively unaffected homozygous APOEε4 carriers had significantly lower FN1 deposition and less reactive gliosis compared to homozygous APOEε4 carriers with AD, suggesting that FN1 might be a downstream driver of APOEε4 -mediated AD-related pathology and cognitive decline. To validate our findings, we used zebrafish models with loss-of-function (LOF) mutations in fn1b - the ortholog for human FN1 . We found that fibronectin LOF reduced gliosis, enhanced gliovascular remodeling and potentiated the microglial response, suggesting that pathological accumulation of FN1 could impair toxic protein clearance, which is ameliorated with FN1 LOF. Our study suggests vascular deposition of FN1 is related to the pathogenicity of APOEε4 , LOF variants in FN1 may reduce APOEε4 -related AD risk, providing novel clues to potential therapeutic interventions targeting the ECM to mitigate AD risk.

Identification of Rare Damaging Missense and Loss of Function Variants in GWAS Loci Using Genome Sequencing Data from Two Cohorts of Familial Late-Onset Alzheimer's Disease

PURPOSE

Few rare pathogenic variants have been identified in the 70+ genetic loci from genome wide association studies of late-onset Alzheimer's disease (AD), limiting research on underlying mechanisms, risk assessment, and genetic counseling.

METHODS

Using genome sequencing data from 197 families in The National Institute on Aging Alzheimer's Disease Family Based Study (AD-FBS), and 214 families in The Estudio Familiar de la Influencia Genética en Alzheimer (EFIGA), we characterized rare coding variants predicted to highly damaging missense or loss of function variants (LoF) within known GWAS loci.

RESULTS

Eight coding and one LoF variant segregated in 10 (5.1%) AD-FBS families and 16 coding and two LoF variants segregated in 18 (8.4%) EFIGA families. ABCA7 and AKAP9 contained the most damaging variants. In 51 (25.9%) of the AD-FBS and in 26 (12.1%) of the EFIGA families, APOE-ε4 was the only variant segregating with familial AD (fAD). Neither APOE-ε4 nor missense or LoF variants were found in 44.1% of the AD-FBS and 62.1% of the EFIGA families.

CONCLUSIONS

Although rare variants were found in both family groups, many families had no gene variant segregating within the family, indicating that the genetic basis for AD has yet to be fully defined.

Evolution of structural rearrangements in prostate cancer intracranial metastases

Intracranial metastases in prostate cancer are uncommon but clinically aggressive. A detailed molecular characterization of prostate cancer intracranial metastases would improve our understanding of their pathogenesis and the search for new treatment strategies. We evaluated the clinical and molecular characteristics of 36 patients with metastatic prostate cancer to either the dura or brain parenchyma. We performed whole genome sequencing (WGS) of 10 intracranial prostate cancer metastases, as well as WGS of primary prostate tumors from men who later developed metastatic disease (n = 6) and nonbrain prostate cancer metastases (n = 36). This first whole genome sequencing study of prostate intracranial metastases led to several new insights. First, there was a higher diversity of complex structural alterations in prostate cancer intracranial metastases compared to primary tumor tissues. Chromothripsis and chromoplexy events seemed to dominate, yet there were few enrichments of specific categories of structural variants compared with non-brain metastases. Second, aberrations involving the AR gene, including AR enhancer gain were observed in 7/10 (70%) of intracranial metastases, as well as recurrent loss of function aberrations involving TP53 in 8/10 (80%), RB1 in 2/10 (20%), BRCA2 in 2/10 (20%), and activation of the PI3K/AKT/PTEN pathway in 8/10 (80%). These alterations were frequently present in tumor tissues from other sites of disease obtained concurrently or sequentially from the same individuals. Third, clonality analysis points to genomic factors and evolutionary bottlenecks that contribute to metastatic spread in patients with prostate cancer. These results describe the aggressive molecular features underlying intracranial metastasis that may inform future diagnostic and treatment approaches.

Precise characterization of somatic complex structural variations from tumor/control paired long-read sequencing data with nanomonsv

We present our novel software, nanomonsv, for detecting somatic structural variations (SVs) using tumor and matched control long-read sequencing data with a single-base resolution. The current version of nanomonsv includes two detection modules, Canonical SV module, and Single breakend SV module. Using tumor/control paired long-read sequencing data from three cancer and their matched lymphoblastoid lines, we demonstrate that Canonical SV module can identify somatic SVs that can be captured by short-read technologies with higher precision and recall than existing methods. In addition, we have developed a workflow to classify mobile element insertions while elucidating their in-depth properties, such as 5' truncations, internal inversions, as well as source sites for 3' transductions. Furthermore, Single breakend SV module enables the detection of complex SVs that can only be identified by long-reads, such as SVs involving highly-repetitive centromeric sequences, and LINE1- and virus-mediated rearrangements. In summary, our approaches applied to cancer long-read sequencing data can reveal various features of somatic SVs and will lead to a better understanding of mutational processes and functional consequences of somatic SVs.

Mapinsights: deep exploration of quality issues and error profiles in high-throughput sequence data

High-throughput sequencing (HTS) has revolutionized science by enabling super-fast detection of genomic variants at base-pair resolution. Consequently, it poses the challenging problem of identification of technical artifacts, i.e. hidden non-random error patterns. Understanding the properties of sequencing artifacts holds the key in separating true variants from false positives. Here, we develop Mapinsights, a toolkit that performs quality control (QC) analysis of sequence alignment files, capable of detecting outliers based on sequencing artifacts of HTS data at a deeper resolution compared with existing methods. Mapinsights performs a cluster analysis based on novel and existing QC features derived from the sequence alignment for outlier detection. We applied Mapinsights on community standard open-source datasets and identified various quality issues including technical errors related to sequencing cycles, sequencing chemistry, sequencing libraries and across various orthogonal sequencing platforms. Mapinsights also enables identification of anomalies related to sequencing depth. A logistic regression-based model built on the features of Mapinsights shows high accuracy in detecting 'low-confidence' variant sites. Quantitative estimates and probabilistic arguments provided by Mapinsights can be utilized in identifying errors, bias and outlier samples, and also aid in improving the authenticity of variant calls.

OpenPBTA: The Open Pediatric Brain Tumor Atlas

Pediatric brain and spinal cancers are collectively the leading disease-related cause of death in children; thus, we urgently need curative therapeutic strategies for these tumors. To accelerate such discoveries, the Children's Brain Tumor Network (CBTN) and Pacific Pediatric Neuro-Oncology Consortium (PNOC) created a systematic process for tumor biobanking, model generation, and sequencing with immediate access to harmonized data. We leverage these data to establish OpenPBTA, an open collaborative project with over 40 scalable analysis modules that genomically characterize 1,074 pediatric brain tumors. Transcriptomic classification reveals universal TP53 dysregulation in mismatch repair-deficient hypermutant high-grade gliomas and TP53 loss as a significant marker for poor overall survival in ependymomas and H3 K28-mutant diffuse midline gliomas. Already being actively applied to other pediatric cancers and PNOC molecular tumor board decision-making, OpenPBTA is an invaluable resource to the pediatric oncology community.

Progranulin mutations in clinical and neuropathological Alzheimer's disease

INTRODUCTION

Progranulin (GRN) mutations occur in frontotemporal lobar degeneration (FTLD) and in Alzheimer's disease (AD), often with TDP-43 pathology.

METHODS

We determined the frequency of rs5848 and rare, pathogenic GRN mutations in two autopsy and one family cohort. We compared Braak stage, β-amyloid load, hyperphosphorylated tau (PHFtau) tangle density and TDP-43 pathology in GRN carriers and non-carriers.

RESULTS

Pathogenic GRN mutations were more frequent in all cohorts compared to the Genome Aggregation Database (gnomAD), but there was no evidence for association with AD. Pathogenic GRN carriers had significantly higher PHFtau tangle density adjusting for age, sex and APOE ε4 genotype. AD patients with rs5848 had higher frequencies of hippocampal sclerosis and TDP-43 deposits. Twenty-two rare, pathogenic GRN variants were observed in the family cohort.

DISCUSSION

GRN mutations in clinical and neuropathological AD increase the burden of tau-related brain pathology but show no specific association with β-amyloid load or AD.

Implementation of Whole-Genome and Transcriptome Sequencing Into Clinical Cancer Care

PURPOSE

The combination of whole-genome and transcriptome sequencing (WGTS) is expected to transform diagnosis and treatment for patients with cancer. WGTS is a comprehensive precision diagnostic test that is starting to replace the standard of care for oncology molecular testing in health care systems around the world; however, the implementation and widescale adoption of this best-in-class testing is lacking.

METHODS

Here, we address the barriers in integrating WGTS for cancer diagnostics and treatment selection and answer questions regarding utility in different cancer types, cost-effectiveness and affordability, and other practical considerations for WGTS implementation.

RESULTS

We review the current studies implementing WGTS in health care systems and provide a synopsis of the clinical evidence and insights into practical considerations for WGTS implementation. We reflect on regulatory, costs, reimbursement, and incidental findings aspects of this test.

CONCLUSION

WGTS is an appropriate comprehensive clinical test for many tumor types and can replace multiple, cascade testing approaches currently performed. Decreasing sequencing cost, increasing number of clinically relevant aberrations and discovery of more complex biomarkers of treatment response, should pave the way for health care systems and laboratories in implementing WGTS into clinical practice, to transform diagnosis and treatment for patients with cancer.

TLsub: A transfer learning based enhancement to accurately detect mutations with wide-spectrum sub-clonal proportion

Mutation detecting is a routine work for sequencing data analysis and the trading of existing tools often involves the combinations of signals on a set of overlapped sequencing reads. However, the subclonal mutations, which are reported to contribute to tumor recurrence and metastasis, are sometimes eliminated by existing signals. When the clonal proportion decreases, signals often present ambiguous, while complicated interactions among signals break the IID assumption for most of the machine learning models. Although the mutation callers could lower the thresholds, false positives are significantly introduced. The main aim here was to detect the subclonal mutations with high specificity from the scenario of ambiguous sample purities or clonal proportions. We proposed a novel machine learning approach for filtering false positive calls to accurately detect mutations with wide spectrum subclonal proportion. We have carried out a series of experiments on both simulated and real datasets, and compared to several state-of-art approaches, including freebayes, MuTect2, Sentieon and SiNVICT. The results demonstrated that the proposed method adapts well to different diluted sequencing signals and can significantly reduce the false positive when detecting subclonal mutations. The codes have been uploaded at https://github.com/TrinaZ/TL-fpFilter for academic usage only.

From Samples to Germline and Somatic Sequence Variation: A Focus on Next-Generation Sequencing in Melanoma Research

Next-generation sequencing (NGS) applications have flourished in the last decade, permitting the identification of cancer driver genes and profoundly expanding the possibilities of genomic studies of cancer, including melanoma. Here we aimed to present a technical review across many of the methodological approaches brought by the use of NGS applications with a focus on assessing germline and somatic sequence variation. We provide cautionary notes and discuss key technical details involved in library preparation, the most common problems with the samples, and guidance to circumvent them. We also provide an overview of the sequence-based methods for cancer genomics, exposing the pros and cons of targeted sequencing vs. exome or whole-genome sequencing (WGS), the fundamentals of the most common commercial platforms, and a comparison of throughputs and key applications. Details of the steps and the main software involved in the bioinformatics processing of the sequencing results, from preprocessing to variant prioritization and filtering, are also provided in the context of the full spectrum of genetic variation (SNVs, indels, CNVs, structural variation, and gene fusions). Finally, we put the emphasis on selected bioinformatic pipelines behind (a) short-read WGS identification of small germline and somatic variants, (b) detection of gene fusions from transcriptomes, and (c) de novo assembly of genomes from long-read WGS data. Overall, we provide comprehensive guidance across the main methodological procedures involved in obtaining sequencing results for the most common short- and long-read NGS platforms, highlighting key applications in melanoma research.

A multi-platform reference for somatic structural variation detection

Accurate detection of somatic structural variation (SV) in cancer genomes remains a challenging problem. This is in part due to the lack of high-quality, gold-standard datasets that enable the benchmarking of experimental approaches and bioinformatic analysis pipelines. Here, we performed somatic SV analysis of the paired melanoma and normal lymphoblastoid COLO829 cell lines using four different sequencing technologies. Based on the evidence from multiple technologies combined with extensive experimental validation, we compiled a comprehensive set of carefully curated and validated somatic SVs, comprising all SV types. We demonstrate the utility of this resource by determining the SV detection performance as a function of tumor purity and sequence depth, highlighting the importance of assessing these parameters in cancer genomics projects. The truth somatic SV dataset as well as the underlying raw multi-platform sequencing data are freely available and are an important resource for community somatic benchmarking efforts.

Whole-genome sequencing illuminates multifaceted targets of selection to humic substances in Eurasian perch

Extreme environments are inhospitable to the majority of species, but some organisms are able to survive in such hostile conditions due to evolutionary adaptations. For example, modern bony fishes have colonized various aquatic environments, including perpetually dark, hypoxic, hypersaline and toxic habitats. Eurasian perch (Perca fluviatilis) is among the few fish species of northern latitudes that is able to live in very acidic humic lakes. Such lakes represent almost "nocturnal" environments; they contain high levels of dissolved organic matter, which in addition to creating a challenging visual environment, also affects a large number of other habitat parameters and biotic interactions. To reveal the genomic targets of humic-associated selection, we performed whole-genome sequencing of perch originating from 16 humic and 16 clear-water lakes in northern Europe. We identified over 800,000 SNPs, of which >10,000 were identified as potential candidates under selection (associated with >3,000 genes) using multiple outlier approaches. Our findings suggest that adaptation to the humic environment may involve hundreds of regions scattered across the genome. Putative signals of adaptation were detected in genes and gene families with diverse functions, including organism development and ion transportation. The observed excess of variants under selection in regulatory regions highlights the importance of adaptive evolution via regulatory elements, rather than via protein sequence modification. Our study demonstrates the power of whole-genome analysis to illuminate multifaceted nature of humic adaptation and provides the foundation for further investigation of causal mutations underlying phenotypic traits of ecological and evolutionary importance.

OUP accepted manuscript

Cell therapy using induced pluripotent stem cell (iPSC) derivatives may result in abnormal tissue generation because the cells undergo numerous cycles of mitosis before clinical application, potentially increasing the accumulation of genetic abnormalities. Therefore, genetic tests may predict abnormal tissue formation after transplantation. Here, we administered iPSC derivatives with or without single-nucleotide variants (SNVs) and deletions in cancer-related genes with various genomic copy number variant (CNV) profiles into immunodeficient mice and examined the relationships between mutations and abnormal tissue formation after transplantation. No positive correlations were found between the presence of SNVs/deletions and the formation of abnormal tissues; the overall predictivity was 29%. However, a copy number higher than 3 was correlated, with an overall predictivity of 86%. Furthermore, we found CNV hotspots at 14q32.33 and 17q12 loci. Thus, CNV analysis may predict abnormal tissue formation after transplantation of iPSC derivatives and reduce the number of tumorigenicity tests.

Validation of HER2 Status in Whole Genome Sequencing Data of Breast Cancers with the Ploidy-Corrected Copy Number Approach

BACKGROUND AND OBJECTIVE

Human epidermal growth factor receptor 2 (HER2) protein overexpression is one of the most significant biomarkers for breast cancer diagnostics, treatment prediction, and prognostics. The high accessibility of HER2 inhibitors in routine clinical practice directly translates into the diagnostic need for precise and robust marker identification. Even though multigene next-generation sequencing methodologies have slowly taken over the field of single-biomarker molecular tests, the copy number alterations such as amplification of the HER2-coding ERBB2 gene are hard to validate on next-generation sequencing platforms as they are characterized by chromosomal structural heterogeneity, polysomy, and genomic context of ploidy. In our study, we tested the approach of using whole genome sequencing instead of next-generation sequencing panels to determine HER2 status in the clinical set-up.

METHODS

We used a large dataset of 876 patients with breast cancer whole genomes with curated clinical data and an additional set of 551 patients' external genomic data. We used the decision-tree-based algorithm for optimization of the diagnostic tool for HER2 status assessment by whole genome sequencing.

RESULTS

The most efficient approach to assess HER2 status in whole genome sequencing data was the ploidy-corrected copy number, utilizing ERBB2 copy number and mean tumor ploidy. The classifier achieved sensitivity of 91.18% and specificity of 98.69% on the internal validation dataset and 89.86% and 96.06% on the external data, which is similar to other next-generation sequencing methods, currently tested in the clinic.

CONCLUSIONS

We provide evidence that the HER2 status may be reliably determined by whole genome sequencing and is applicable across different laboratory protocols and pipelines. We suggest using the ploidy-corrected copy number for diagnostic purposes.

New York's Polyethnic-1000: a regional initiative to understand how diverse ancestries influence the risk, progression, and treatment of cancers

Research consortia can help to repair deficiencies in knowledge about the influence of inherited genetic diversity on disease. The New York Genome Center (NYGC) recently established Polyethnic-1000 (P-1000), a multi-institutional collaboration to study hereditary factors affecting several types of cancer. Here, we describe its rationale, organization, development, current activities, and prospects.

Batch effects in population genomic studies with low-coverage whole genome sequencing data: Causes, detection and mitigation

Over the past few decades, there has been an explosion in the amount of publicly available sequencing data. This opens new opportunities for combining data sets to achieve unprecedented sample sizes, spatial coverage or temporal replication in population genomic studies. However, a common concern is that nonbiological differences between data sets may generate patterns of variation in the data that can confound real biological patterns, a problem known as batch effects. In this paper, we compare two batches of low-coverage whole genome sequencing (lcWGS) data generated from the same populations of Atlantic cod (Gadus morhua). First, we show that with a "batch-effect-naive" bioinformatic pipeline, batch effects systematically biased our genetic diversity estimates, population structure inference and selection scans. We then demonstrate that these batch effects resulted from multiple technical differences between our data sets, including the sequencing chemistry (four-channel vs. two-channel), sequencing run, read type (single-end vs. paired-end), read length (125 vs. 150 bp), DNA degradation level (degraded vs. well preserved) and sequencing depth (0.8× vs. 0.3× on average). Lastly, we illustrate that a set of simple bioinformatic strategies (such as different read trimming and single nucleotide polymorphism filtering) can be used to detect batch effects in our data and substantially mitigate their impact. We conclude that combining data sets remains a powerful approach as long as batch effects are explicitly accounted for. We focus on lcWGS data in this paper, which may be particularly vulnerable to certain causes of batch effects, but many of our conclusions also apply to other sequencing strategies.

Embryonic lethal genetic variants and chromosomally normal pregnancy loss

OBJECTIVE

To examine whether rare damaging genetic variants are associated with chromosomally normal pregnancy loss and estimate the magnitude of the association.

DESIGN

Case-control.

SETTING

Cases were derived from a consecutive series of karyotyped losses at one New Jersey hospital. Controls were derived from the National Database for Autism Research.

PATIENT(S)

Cases comprised 19 chromosomally normal loss conceptus-parent trios. Controls comprised 547 unaffected siblings of autism case-parent trios.

INTERVENTION(S)

None.

MAIN OUTCOME MEASURE(S)

The rate of damaging variants in the exome (loss of function and missense-damaging) and the proportions of probands with at least one such variant among cases vs. controls.

RESULTS

The proportions of probands with at least one rare damaging variant were 36.8% among cases and 22.9% among controls (odds ratio, 2.0; 99% confidence interval, 0.5-7.3). No case had a variant in a known fetal anomaly gene. The proportion with variants in possibly embryonic lethal genes increased in case probands (odds ratio, 14.5; 99% confidence interval, 1.5-89.7); variants occurred in BAZ1A, FBN2, and TIMP2.

CONCLUSION(S)

Rare genetic variants in the conceptus may be a cause of chromosomally normal pregnancy loss. A larger sample is needed to estimate the magnitude of the association with precision and identify relevant biologic pathways.

Prospects and challenges of cancer systems medicine: from genes to disease networks

It is becoming evident that holistic perspectives toward cancer are crucial in deciphering the overwhelming complexity of tumors. Single-layer analysis of genome-wide data has greatly contributed to our understanding of cellular systems and their perturbations. However, fundamental gaps in our knowledge persist and hamper the design of effective interventions. It is becoming more apparent than ever, that cancer should not only be viewed as a disease of the genome but as a disease of the cellular system. Integrative multilayer approaches are emerging as vigorous assets in our endeavors to achieve systemic views on cancer biology. Herein, we provide a comprehensive review of the approaches, methods and technologies that can serve to achieve systemic perspectives of cancer. We start with genome-wide single-layer approaches of omics analyses of cellular systems and move on to multilayer integrative approaches in which in-depth descriptions of proteogenomics and network-based data analysis are provided. Proteogenomics is a remarkable example of how the integration of multiple levels of information can reduce our blind spots and increase the accuracy and reliability of our interpretations and network-based data analysis is a major approach for data interpretation and a robust scaffold for data integration and modeling. Overall, this review aims to increase cross-field awareness of the approaches and challenges regarding the omics-based study of cancer and to facilitate the necessary shift toward holistic approaches.

Somatic variant analysis of linked-reads sequencing data with Lancet

SUMMARY

We present a new version of the popular somatic variant caller, Lancet, that supports the analysis of linked-reads sequencing data. By seamlessly integrating barcodes and haplotype read assignments within the colored De Bruijn graph local-assembly framework, Lancet computes a barcode-aware coverage and identifies variants that disagree with the local haplotype structure.

AVAILABILITY AND IMPLEMENTATION

Lancet is implemented in C++ and available for academic and non-commercial research purposes as an open-source package at https://github.com/nygenome/lancet.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Phage display: an ideal platform for coupling protein to nucleic acid

Display technology, especially phage display technology, has been widely applied in many fields. The theoretical core of display technology is the physical linkage between the protein/peptide on the surface of a phage and the coding DNA sequence inside the same phage. Starting from phage-displayed peptide/protein/antibody libraries and taking advantage of the ever-growing power of next-generation sequencing (NGS) for DNA sequencing/decoding, rich protein-related information can easily be obtained in a high-throughput way. Based on this information, many scientific and clinical questions can be readily addressed. In the past few years, aided by the development of NGS, droplet technology, and massive oligonucleotide synthesis, we have witnessed and continue to witness large advances of phage display technology, in both technology development and application. The aim of this review is to summarize and discuss these recent advances.

Sequencing platform shifts provide opportunities but pose challenges for combining genomic data sets

Technological advances in DNA sequencing over the last decade now permit the production and curation of large genomic data sets in an increasing number of nonmodel species. Additionally, these new data provide the opportunity for combining data sets, resulting in larger studies with a broader taxonomic range. Whilst the development of new sequencing platforms has been beneficial, resulting in a higher throughput of data at a lower per-base cost, shifts in sequencing technology can also pose challenges for those wishing to combine new sequencing data with data sequenced on older platforms. Here, we outline the types of studies where the use of curated data might be beneficial, and highlight potential biases that might be introduced by combining data from different sequencing platforms. As an example of the challenges associated with combining data across sequencing platforms, we focus on the impact of the shift in Illumina's base calling technology from a four-channel system to a two-channel system. We caution that when data are combined from these two systems, erroneous guanine base calls that result from the two-channel chemistry can make their way through a bioinformatic pipeline, eventually leading to inaccurate and potentially misleading conclusions. We also suggest solutions for dealing with such potential artefacts, which make samples sequenced on different sequencing platforms appear more differentiated from one another than they really are. Finally, we stress the importance of archiving tissue samples and the associated sequences for the continued reproducibility and reusability of sequencing data in the face of ever-changing sequencing platform technology.

Cross-Site Concordance Evaluation of Tumor DNA and RNA Sequencing Platforms for the CIMAC-CIDC Network

PURPOSE

Whole-exome (WES) and RNA sequencing (RNA-seq) are key components of cancer immunogenomic analyses. To evaluate the consistency of tumor WES and RNA-seq profiling platforms across different centers, the Cancer Immune Monitoring and Analysis Centers (CIMAC) and the Cancer Immunologic Data Commons (CIDC) conducted a systematic harmonization study.

EXPERIMENTAL DESIGN

DNA and RNA were centrally extracted from fresh frozen and formalin-fixed paraffin-embedded non-small cell lung carcinoma tumors and distributed to three centers for WES and RNA-seq profiling. In addition, two 10-plex HapMap cell line pools with known mutations were used to evaluate the accuracy of the WES platforms.

RESULTS

The WES platforms achieved high precision (> 0.98) and recall (> 0.87) on the HapMap pools when evaluated on loci using > 50× common coverage. Nonsynonymous mutations clustered by tumor sample, achieving an index of specific agreement above 0.67 among replicates, centers, and sample processing. A DV200 > 24% for RNA, as a putative presequencing RNA quality control (QC) metric, was found to be a reliable threshold for generating consistent expression readouts in RNA-seq and NanoString data. MedTIN > 30 was likewise assessed as a reliable RNA-seq QC metric, above which samples from the same tumor across replicates, centers, and sample processing runs could be robustly clustered and HLA typing, immune infiltration, and immune repertoire inference could be performed.

CONCLUSIONS

The CIMAC collaborating laboratory platforms effectively generated consistent WES and RNA-seq data and enable robust cross-trial comparisons and meta-analyses of highly complex immuno-oncology biomarker data across the NCI CIMAC-CIDC Network.

Shiny-SoSV: A web-based performance calculator for somatic structural variant detection

Somatic structural variants are an important contributor to cancer development and evolution. Accurate detection of these complex variants from whole genome sequencing data is influenced by a multitude of parameters. However, there are currently no tools for guiding study design nor are there applications that could predict the performance of somatic structural variant detection. To address this gap, we developed Shiny-SoSV, a user-friendly web-based calculator for determining the impact of common variables on the sensitivity, precision and F1 score of somatic structural variant detection, including choice of variant detection tool, sequencing depth of coverage, variant allele fraction, and variant breakpoint resolution. Using simulation studies, we determined singular and combinatoric effects of these variables, modelled the results using a generalised additive model, allowing structural variant detection performance to be predicted for any combination of predictors. Shiny-SoSV provides an interactive and visual platform for users to easily compare individual and combined impact of different parameters. It predicts the performance of a proposed study design, on somatic structural variant detection, prior to the commencement of benchwork. Shiny-SoSV is freely available at https://hcpcg.shinyapps.io/Shiny-SoSV with accompanying user’s guide and example use-cases.