Significance Associated with Phenotype Score Aids in Variant Prioritization for Exome Sequencing Analysis

Several in silico annotation-based methods have been developed to prioritize variants in exome sequencing analysis. This study introduced a novel metric Significance Associated with Phenotypes (SAP) score, which generates a statistical score by comparing an individual's observed phenotypes against existing gene-phenotype associations. To evaluate the SAP score, a retrospective analysis was performed on 219 exomes. Among them, 82 family-based and 35 singleton exomes had at least one disease-causing variant that explained the patient's clinical features. SAP scores were calculated, and the rank of the disease-causing variant was compared with a known method, Exomiser. Using the SAP score, the known causative variant was ranked in the top 10 retained variants for 94% (77 of 82) of the family-based exomes and in first place for 73% of these cases. For singleton exomes, the SAP score analysis ranked the known pathogenic variants within the top 10 for 80% (28 of 35) of cases. The SAP score, which is independent of detected variants, demonstrates comparable performance with Exomiser, which considers both phenotype and variant-level evidence simultaneously. Among 102 cases with negative results or variants of uncertain significance, SAP score analysis revealed two cases with a potential new diagnosis based on rank. The SAP score, a phenotypic quantitative metric, can be used in conjunction with standard variant filtration and annotation to enhance variant prioritization in exome analysis.

Mitochondrial DNA haplogroup, genetic ancestry, and susceptibility to Ewing sarcoma

Based on current studies, the incidence of Ewing sarcoma (ES) varies significantly by race and ethnicity, with the disease being most common in patients of European ancestry. However, race/ethnicity has generally been self-reported rather than formally evaluated at a population level using DNA evidence. Additionally, mitochondrial dysfunction is a hallmark of ES, yet there have been no reported studies of mitochondrial genetics in ES. Thus, we evaluated both the mitochondrial and nuclear ancestries of 420 pediatric ES patients in the United States using whole-genome sequencing. We found that the mitochondrial DNA (mtDNA) genomes of only six (1.4 %) patients belonged to African L haplogroups, while those of 90 % of the patients belonged to macrohaplogroup R, which includes haplogroup H, the most common maternal lineage in Europe. Compared to the general US population, European haplogroups were significantly enriched in ES patients (p < 2.2e-16) and the African haplogroups are significantly impoverished (p < 4.6e-16). Using the ancestry informative markers defined in a National Genographic study, the vast majority of patients exhibited significant nuclear ancestry originating from the Mediterranean, Northern Europe, and Southwest Asia, including all six patients with African L mtDNAs. Very few had primarily African nuclear ancestry. This is the first genomic epidemiology study to simultaneously interrogate the mitochondrial and nuclear ancestries of ES patients. While supporting previous findings of enriched European ancestry in ES patients, these results also suggest alternative hypotheses for the significant contribution of mitochondrial ancestry in ES patients, as well as the protective role of African ancestry.

High Prevalence of SARS-CoV-2 Genetic Variation and D614G Mutation in Pediatric Patients With COVID-19

BACKGROUND

The full spectrum of the disease phenotype and viral genotype of coronavirus disease 2019 (COVID-19) have yet to be thoroughly explored in children. Here, we analyze the relationships between viral genetic variants and clinical characteristics in children.

METHODS

Whole-genome sequencing was performed on respiratory specimens collected for all SARS-CoV-2-positive children (n = 141) between March 13 and June 16, 2020. Viral genetic variations across the SARS-CoV-2 genome were identified and investigated to evaluate genomic correlates of disease severity.

RESULTS

Higher viral load was detected in symptomatic patients (P = .0007) and in children <5 years old (P = .0004). Genomic analysis revealed a mean pairwise difference of 10.8 single nucleotide variants (SNVs), and the majority (55.4%) of SNVs led to an amino acid change in the viral proteins. The D614G mutation in the spike protein was present in 99.3% of the isolates. The calculated viral mutational rate of 22.2 substitutions/year contrasts the 13.5 substitutions/year observed in California isolates without the D614G mutation. Phylogenetic clade 20C was associated with severe cases of COVID-19 (odds ratio, 6.95; P = .0467). Epidemiological investigation revealed major representation of 3 of 5 major Nextstrain clades (20A, 20B, and 20C) consistent with multiple introductions of SARS-CoV-2 in Southern California.

CONCLUSIONS

Genomic evaluation demonstrated greater than expected genetic diversity, presence of the D614G mutation, increased mutation rate, and evidence of multiple introductions of SARS-CoV-2 into Southern California. Our findings suggest a possible association of phylogenetic clade 20C with severe disease, but small sample size precludes a definitive conclusion. Our study warrants larger and multi-institutional genomic evaluation and has implications for infection control practices.

Increased viral variants in children and young adults with impaired humoral immunity and persistent SARS-CoV-2 infection: A consecutive case series

BACKGROUND

There is increasing concern that persistent infection of SARS-CoV-2 within immunocompromised hosts could serve as a reservoir for mutation accumulation and subsequent emergence of novel strains with the potential to evade immune responses.

METHODS

We describe three patients with acute lymphoblastic leukemia who were persistently positive for SARS-CoV-2 by real-time polymerase chain reaction. Viral viability from longitudinally-collected specimens was assessed. Whole-genome sequencing and serological studies were performed to measure viral evolution and evidence of immune escape.

FINDINGS

We found compelling evidence of ongoing replication and infectivity for up to 162 days from initial positive by subgenomic RNA, single-stranded RNA, and viral culture analysis. Our results reveal a broad spectrum of infectivity, host immune responses, and accumulation of mutations, some with the potential for immune escape.

INTERPRETATION

Our results highlight the potential need to reassess infection control precautions in the management and care of immunocompromised patients. Routine surveillance of mutations and evaluation of their potential impact on viral transmission and immune escape should be considered.

Persistent SARS-CoV-2 infection and increasing viral variants in children and young adults with impaired humoral immunity

Background

There is increasing concern that persistent infection of SARS-CoV-2 within immunocompromised hosts could serve as a reservoir for mutation accumulation and subsequent emergence of novel strains with the potential to evade immune responses.

Methods

We describe three patients with acute lymphoblastic leukemia who were persistently positive for SARS-CoV-2 by real-time polymerase chain reaction. Viral viability from longitudinally-collected specimens was assessed. Whole-genome sequencing and serological studies were performed to measure viral evolution and evidence of immune escape.

Findings

We found compelling evidence of ongoing replication and infectivity for up to 162 days from initial positive by subgenomic RNA, single-stranded RNA, and viral culture analysis. Our results reveal a broad spectrum of infectivity, host immune responses, and accumulation of mutations, some with the potential for immune escape.

Interpretation

Our results highlight the need to reassess infection control precautions in the management and care of immunocompromised patients. Routine surveillance of mutations and evaluation of their potential impact on viral transmission and immune escape should be considered.

Funding

The work was partially funded by The Saban Research Institute at Children's Hospital Los Angeles intramural support for COVID-19 Directed Research (X.G. and J.D.B.), the Johns Hopkins Center of Excellence in Influenza Research and Surveillance HHSN272201400007C (A.P.), NIH/NIAID R01AI127877 (S.D.B.), NIH/NIAID R01AI130398 (S.D.B.), NIH 1U54CA260517 (S.D.B.), an endowment to S.D.B. from the Crown Family Foundation, an Early Postdoc.Mobility Fellowship Stipend to O.F.W. from the Swiss National Science Foundation (SNSF), and a Coulter COVID-19 Rapid Response Award to S.D.B. L.G. is a SHARE Research Fellow in Pediatric Hematology-Oncology.

Early pandemic molecular diversity of SARS-CoV-2 in children

Background

In the US, community circulation of the SARS-CoV-2 virus likely began in February 2020 after mostly travel-related cases. Children's Hospital of Philadelphia began testing on 3/9/2020 for pediatric and adult patients, and for all admitted patients on 4/1/2020, allowing an early glimpse into the local molecular epidemiology of the virus.

Methods

We obtained 169 SARS-CoV-2 samples (83 from patients <21 years old) from March through May and produced whole genome sequences. We used genotyping tools to track variants over time and to test for possible genotype associated clinical presentations and outcomes in children.

Results

Our analysis uncovered 13 major lineages that changed in relative abundance as cases peaked in mid-April in Philadelphia. We detected at least 6 introductions of distinct viral variants into the population. As a group, children had more diverse virus genotypes than the adults tested. No strong differences in clinical variables were associated with genotypes.

Conclusions

Whole genome analysis revealed unexpected diversity, and distinct circulating viral variants within the initial peak of cases in Philadelphia. Most introductions appeared to be local from nearby states. Although limited by sample size, we found no evidence that different genotypes had different clinical impacts in children in this study.

Abstract S08-01: Highly sensitive and full-genome interrogation of SARS-CoV-2 using multiplexed PCR enrichment followed by next-generation sequencing

Many detection methods have been used or reported for the diagnosis and/or surveillance of COVID-19. Among them, reverse transcription polymerase chain reaction (RT-PCR) is the most commonly used because of its high sensitivity, typically claiming detection of about 5 copies of viruses. However, it has been reported that only 47-59% of the positive cases were identified by some RT-PCR methods, probably due to low viral load, timing of sampling, degradation of virus RNA in the sampling process, or possible mutations spanning the primer binding sites. Therefore, alternative and highly sensitive methods are imperative. With the goal of improving sensitivity and accommodating various application settings, we developed a multiplex-PCR-based method comprising 343 pairs of specific primers and demonstrated its efficiency at detecting SARS-CoV-2 at low copy numbers. The assay produced clean characteristic target peaks of defined sizes, which allowed for direct identification of positives by electrophoresis. We further amplified the entire SARS-CoV-2 genome from 8 to half a million viral copies purified from 13 COVID-19 positive specimens and detected mutations through next-generation sequencing. Finally, we developed a multiplex-PCR-based metagenomic method in parallel that required modest sequencing depth for uncovering SARS-CoV-2 mutational diversity and potentially novel or emerging isolates.

Citation Format: Chenyu Li, David N. Debruyne, Julia Spencer, Vidushi Kapoor, Lily Y. Liu, Bo Zhou, Utsav Pandey, Moiz Bootwalla, Dejerianne Ostrow, Dennis T. Maglinte, David Ruble, Alex Ryutov, Lishuang Shen, Lucie Lee, Rounak Feigelman, Grayson Burdon, Jeffrey Liu, Alejandra Oliva, Adam Borcherding, Hongdong Tan, Alexander E. Urban, Xiaowu Gai, Jennifer Dien Bard, Guoying Liu, Zhitong Liu. Highly sensitive and full-genome interrogation of SARS-CoV-2 using multiplexed PCR enrichment followed by next-generation sequencing [abstract]. In: Proceedings of the AACR Virtual Meeting: COVID-19 and Cancer; 2020 Jul 20-22. Philadelphia (PA): AACR; Clin Cancer Res 2020;26(18_Suppl):Abstract nr S08-01.

Pan-cancer analysis of whole genomes

Cancer is driven by genetic change, and the advent of massively parallel sequencing has enabled systematic documentation of this variation at the whole-genome scale1-3. Here we report the integrative analysis of 2,658 whole-cancer genomes and their matching normal tissues across 38 tumour types from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). We describe the generation of the PCAWG resource, facilitated by international data sharing using compute clouds. On average, cancer genomes contained 4-5 driver mutations when combining coding and non-coding genomic elements; however, in around 5% of cases no drivers were identified, suggesting that cancer driver discovery is not yet complete. Chromothripsis, in which many clustered structural variants arise in a single catastrophic event, is frequently an early event in tumour evolution; in acral melanoma, for example, these events precede most somatic point mutations and affect several cancer-associated genes simultaneously. Cancers with abnormal telomere maintenance often originate from tissues with low replicative activity and show several mechanisms of preventing telomere attrition to critical levels. Common and rare germline variants affect patterns of somatic mutation, including point mutations, structural variants and somatic retrotransposition. A collection of papers from the PCAWG Consortium describes non-coding mutations that drive cancer beyond those in the TERT promoter4; identifies new signatures of mutational processes that cause base substitutions, small insertions and deletions and structural variation5,6; analyses timings and patterns of tumour evolution7; describes the diverse transcriptional consequences of somatic mutation on splicing, expression levels, fusion genes and promoter activity8,9; and evaluates a range of more-specialized features of cancer genomes8,10-18.

OncoKids: A Comprehensive Next-Generation Sequencing Panel for Pediatric Malignancies

The OncoKids panel is an amplification-based next-generation sequencing assay designed to detect diagnostic, prognostic, and therapeutic markers across the spectrum of pediatric malignancies, including leukemias, sarcomas, brain tumors, and embryonal tumors. This panel uses low input amounts of DNA (20 ng) and RNA (20 ng) and is compatible with formalin-fixed, paraffin-embedded and frozen tissue, bone marrow, and peripheral blood. The DNA content of this panel covers the full coding regions of 44 cancer predisposition loci, tumor suppressor genes, and oncogenes; hotspots for mutations in 82 genes; and amplification events in 24 genes. The RNA content includes 1421 targeted gene fusions. We describe the validation of this panel by using a large cohort of 192 unique clinical samples that included a wide range of tumor types and alterations. Robust performance was observed for analytical sensitivity, reproducibility, and limit of detection studies. The results from this study support the use of OncoKids for routine clinical testing of a wide variety of pediatric malignancies.

Clinical utility of the low-density Infinium QC genotyping Array in a genomics-based diagnostics laboratory

Background

With 15,949 markers, the low-density Infinium QC Array-24 BeadChip enables linkage analysis, HLA haplotyping, fingerprinting, ethnicity determination, mitochondrial genome variations, blood groups and pharmacogenomics. It represents an attractive independent QC option for NGS-based diagnostic laboratories, and provides cost-efficient means for determining gender, ethnic ancestry, and sample kinships, that are important for data interpretation of NGS-based genetic tests.

Methods

We evaluated accuracy and reproducibility of Infinium QC genotyping calls by comparing them with genotyping data of the same samples from other genotyping platforms, whole genome/exome sequencing. Accuracy and robustness of determining gender, provenance, and kinships were assessed.

Results

Concordance of genotype calls between Infinium QC and other platforms was above 99%. Here we show that the chip’s ancestry informative markers are sufficient for ethnicity determination at continental and sometimes subcontinental levels, with assignment accuracy varying with the coverage for a particular region and ethnic groups. Mean accuracies of provenance prediction at a regional level were varied from 81% for Asia, to 89% for Americas, 86% for Africa, 97% for Oceania, 98% for Europe, and 100% for India. Mean accuracy of ethnicity assignment predictions was 63%. Pairwise concordances of AFR samples with the samples from any other super populations were the lowest (0.39–0.43), while the concordances within the same population were relatively high (0.55–0.61). For all populations except African, cross-population comparisons were similar in their concordance ranges to the range of within-population concordances (0.54–0.57). Gender determination was correct in all tested cases.

Conclusions

Our results indicate that the Infinium QC Array-24 chip is suitable for cost-efficient, independent QC assaying in the settings of an NGS-based molecular diagnostic laboratory; hence, we recommend its integration into the standard laboratory workflow. Low-density chips can provide sample-specific measures for variant call accuracy, prevent sample mix-ups, validate self-reported ethnicities, and detect consanguineous cases. Integration of low-density chips into QC procedures aids proper interpretation of candidate sequence variants. To enhance utility of this low-density chip, we recommend expansion of ADME and mitochondrial markers. Inexpensive Infinium-like low-density human chips have a potential to become a “Swiss army knife” among genotyping assays suitable for many applications requiring high-throughput assays.

Electronic supplementary material

The online version of this article (10.1186/s12920-017-0297-7) contains supplementary material, which is available to authorized users.

Inherited germline ATRX mutation in two brothers with ATR-X syndrome and osteosarcoma

We report a family in which two brothers had an undiagnosed genetic disorder comprised of dysmorphic features, microcephaly, severe intellectual disability (non-verbal), mild anemia, and cryptorchidism. Both developed osteosarcoma. Trio exome sequencing (using blood samples from the younger brother and both parents) was performed and a nonsense NM_000489.4:c.7156C>T (p.Arg2386*) mutation in the ATRX gene was identified in the proband (hemizygous) and in the mother's peripheral blood DNA (heterozygous). The mother is healthy, does not exhibit any clinical manifestations of ATR-X syndrome and there was no family history of cancer. The same hemizygous pathogenic variant was confirmed in the affected older brother's skin tissue by subsequent Sanger sequencing. Chromosomal microarray studies of both brothers' osteosarcomas revealed complex copy number alterations consistent with the clinical diagnosis of osteosarcoma. Recently, somatic mutations in the ATRX gene have been observed as recurrent alterations in both osteosarcoma and brain tumors. However, it is unclear if there is any association between osteosarcoma and germline ATRX mutations, specifically in patients with constitutional ATR-X syndrome. This is the first report of osteosarcoma diagnosed in two males with ATR-X syndrome, suggesting a potential increased risk for cancer in patients with this disorder.

A Rapid and Sensitive Next-Generation Sequencing Method to Detect RB1 Mutations Improves Care for Retinoblastoma Patients and Their Families

Retinoblastoma is a childhood eye malignancy that can lead to the loss of vision, eye(s), and sometimes life. The tumors are initiated by inactivating mutations in both alleles of the tumor-suppressor gene, RB1, or, rarely, by MYCN amplification. Timely identification of a germline RB1 mutation in blood samples or either somatic RB1 mutation or MYCN amplification in tumors is important for effective care and management of retinoblastoma patients and their families. However, current procedures to thoroughly test RB1 mutations are complicated and lengthy. Herein, we report a next-generation sequencing-based method capable of detecting point mutations, small indels, and large deletions or duplications across the entire RB1 gene and amplification of MYCN gene on a single platform. From DNA extraction to clinical interpretation requires only 3 days, enabling early molecular diagnosis of retinoblastoma and optimal treatment outcomes. This method can also detect low-level mosaic mutations in blood samples that can be missed by routine Sanger sequencing. In addition, it can differentiate between RB1 mutation- and MYCN amplification-driven retinoblastomas. This rapid, comprehensive, and sensitive method for detecting RB1 mutations and MYCN amplification can readily identify RB1 mutation carriers and thus improve the management and genetic counseling for retinoblastoma patients and their families.

Comprehensive, Integrative Genomic Analysis of Diffuse Lower-Grade Gliomas

BACKGROUND

Diffuse low-grade and intermediate-grade gliomas (which together make up the lower-grade gliomas, World Health Organization grades II and III) have highly variable clinical behavior that is not adequately predicted on the basis of histologic class. Some are indolent; others quickly progress to glioblastoma. The uncertainty is compounded by interobserver variability in histologic diagnosis. Mutations in IDH, TP53, and ATRX and codeletion of chromosome arms 1p and 19q (1p/19q codeletion) have been implicated as clinically relevant markers of lower-grade gliomas.

METHODS

We performed genomewide analyses of 293 lower-grade gliomas from adults, incorporating exome sequence, DNA copy number, DNA methylation, messenger RNA expression, microRNA expression, and targeted protein expression. These data were integrated and tested for correlation with clinical outcomes.

RESULTS

Unsupervised clustering of mutations and data from RNA, DNA-copy-number, and DNA-methylation platforms uncovered concordant classification of three robust, nonoverlapping, prognostically significant subtypes of lower-grade glioma that were captured more accurately by IDH, 1p/19q, and TP53 status than by histologic class. Patients who had lower-grade gliomas with an IDH mutation and 1p/19q codeletion had the most favorable clinical outcomes. Their gliomas harbored mutations in CIC, FUBP1, NOTCH1, and the TERT promoter. Nearly all lower-grade gliomas with IDH mutations and no 1p/19q codeletion had mutations in TP53 (94%) and ATRX inactivation (86%). The large majority of lower-grade gliomas without an IDH mutation had genomic aberrations and clinical behavior strikingly similar to those found in primary glioblastoma.

CONCLUSIONS

The integration of genomewide data from multiple platforms delineated three molecular classes of lower-grade gliomas that were more concordant with IDH, 1p/19q, and TP53 status than with histologic class. Lower-grade gliomas with an IDH mutation either had 1p/19q codeletion or carried a TP53 mutation. Most lower-grade gliomas without an IDH mutation were molecularly and clinically similar to glioblastoma. (Funded by the National Institutes of Health.).

Epigenetic analysis of KSHV latent and lytic genomes

Epigenetic modifications of the herpesviral genome play a key role in the transcriptional control of latent and lytic genes during a productive viral lifecycle. In this study, we describe for the first time a comprehensive genome-wide ChIP-on-Chip analysis of the chromatin associated with the Kaposi's sarcoma-associated herpesvirus (KSHV) genome during latency and lytic reactivation. Depending on the gene expression class, different combinations of activating [acetylated H3 (AcH3) and H3K4me3] and repressive [H3K9me3 and H3K27me3] histone modifications are associated with the viral latent genome, which changes upon reactivation in a manner that is correlated with their expression. Specifically, both the activating marks co-localize on the KSHV latent genome, as do the repressive marks. However, the activating and repressive histone modifications are mutually exclusive of each other on the bulk of the latent KSHV genome. The genomic region encoding the IE genes ORF50 and ORF48 possesses the features of a bivalent chromatin structure characterized by the concomitant presence of the activating H3K4me3 and the repressive H3K27me3 marks during latency, which rapidly changes upon reactivation with increasing AcH3 and H3K4me3 marks and decreasing H3K27me3. Furthermore, EZH2, the H3K27me3 histone methyltransferase of the Polycomb group proteins (PcG), colocalizes with the H3K27me3 mark on the entire KSHV genome during latency, whereas RTA-mediated reactivation induces EZH2 dissociation from the genomic regions encoding IE and E genes concurrent with decreasing H3K27me3 level and increasing IE/E lytic gene expression. Moreover, either the inhibition of EZH2 expression by a small molecule inhibitor DZNep and RNAi knockdown, or the expression of H3K27me3-specific histone demethylases apparently induced the KSHV lytic gene expression cascade. These data indicate that histone modifications associated with the KSHV latent genome are involved in the regulation of latency and ultimately in the control of the temporal and sequential expression of the lytic gene cascade. In addition, the PcG proteins play a critical role in the control of KSHV latency by maintaining a reversible heterochromatin on the KSHV lytic genes. Thus, the regulation of the spatial and temporal association of the PcG proteins with the KSHV genome may be crucial for propagating the KSHV lifecycle.

KSHV is a ubiquitous herpesvirus that establishes a life-long persistent infection in humans and is associated with Kaposi's sarcoma and several lymphoid malignancies. During latency, the KSHV genome persists as a multicopy circular DNA assembled into nucleosomal structures. While viral latency is characterized by restricted viral gene expression, reactivation induces the lytic replication program and the expression of viral genes in defined sequential and temporal order. Posttranslational modifications of the viral chromatin structure have been implicated to regulate viral gene expressions but the underlying gene regulatory mechanisms are still elusive. Here, we demonstrate that the latent and lytic chromatins of KSHV are associated with a distinctive pattern of activating and repressive histone modifications whose distribution changes upon reactivation in an organized manner in correlation with the temporally ordered expression of viral lytic genes. Furthermore, we demonstrate that the evolutionarily conserved Polycomb group proteins, that maintain the repression of genes involved in hematopoiesis, X-chromosome inactivation, cell proliferation and stem cell differentiation, also play a critical role in the regulation of KSHV latency by maintaining a repressive chromatin structure. Thus, the epigenetic program of KSHV is at the crux of restricting latent gene expression and the orderly expression of lytic genes.

Small bowel cancer. Radiologic diagnosis

The diagnosis of small bowel cancer before it has metastasized continues to be a clinical and radiologic dilemma. This article presents an overview of the role of radiology in the diagnosis of small bowel cancer, describes techniques for optimal radiologic examination, and presents a pictorial essay of the radiologic features and differential diagnoses of these tumors.