Robust Performance of SARS-CoV-2 Whole-Genome Sequencing from Wastewater with a Nonselective Virus Concentration Method.

The sequencing of human virus genomes from wastewater samples is an efficient method for tracking viral transmission and evolution at the community level. However, this requires the recovery of viral nucleic acids of high quality. We developed a reusable tangential-flow filtration system to concentrate and purify viruses from wastewater for genome sequencing. A pilot study was conducted with 94 wastewater samples from four local sewersheds, from which viral nucleic acids were extracted, and the whole genome of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) was sequenced using the ARTIC V4.0 primers. Our method yielded a high probability (0.9) of recovering complete or near-complete SARS-CoV-2 genomes (>90% coverage at 10× depth) from wastewater when the COVID-19 incidence rate exceeded 33 cases per 100 000 people. The relative abundances of sequenced SARS-CoV-2 variants followed the trends observed from patient-derived samples. We also identified SARS-CoV-2 lineages in wastewater that were underrepresented or not present in the clinical whole-genome sequencing data. The developed tangential-flow filtration system can be easily adopted for the sequencing of other viruses in wastewater, particularly those at low concentrations.

Key transcription factors influence the epigenetic landscape to regulate retinal cell differentiation.

How the diverse neural cell types emerge from multipotent neural progenitor cells during central nervous system development remains poorly understood. Recent scRNA-seq studies have delineated the developmental trajectories of individual neural cell types in many neural systems including the neural retina. Further understanding of the formation of neural cell diversity requires knowledge about how the epigenetic landscape shifts along individual cell lineages and how key transcription factors regulate these changes. In this study, we dissect the changes in the epigenetic landscape during early retinal cell differentiation by scATAC-seq and identify globally the enhancers, enriched motifs, and potential interacting transcription factors underlying the cell state/type specific gene expression in individual lineages. Using CUT&Tag, we further identify the enhancers bound directly by four key transcription factors, Otx2, Atoh7, Pou4f2 and Isl1, including those dependent on Atoh7, and uncover the sequential and combinatorial interactions of these factors with the epigenetic landscape to control gene expression along individual retinal cell lineages such as retinal ganglion cells (RGCs). Our results reveal a general paradigm in which transcription factors collaborate and compete to regulate the emergence of distinct retinal cell types such as RGCs from multipotent retinal progenitor cells (RPCs).

Core needle biopsy wash as a tool for acquiring additional diagnostic material for laboratory testing.

AIMS: To determine if a simple prewash step added to the processing workflow of tissue procurement by a core needle biopsy device will recover enough cells to expand the laboratory testing armamentarium. METHODS: Tissue was obtained from unfixed resection specimens using a core needle device and washed in a buffered solution before fixation. This creates a liquid aliquot from which dislodged cells can be kept and separated from the tissue specimen, the latter of which can then undergo traditional formalin-fixed, paraffin-embedded processing. RESULTS: Cells dislodged from the tissue during the biopsy procedure are recoverable, are representative of the tissue section and of sufficient quantities for additional laboratory testing. CONCLUSIONS: The core needle biopsy wash is an under-recognised and underutilised approach to extending the diagnostic capabilities of the limited amount of targeted material obtained during this common procedure. The ability to recover supplemental amounts of diagnostic material yields great potential as a substrate for a multitude of current and developing laboratory assays.

Integrative Genomic Analysis of Pediatric Myeloid-Related Acute Leukemias Identifies Novel Subtypes and Prognostic Indicators.

Genomic characterization of pediatric patients with acute myeloid leukemia (AML) has led to the discovery of somatic mutations with prognostic implications. Although gene-expression profiling can differentiate subsets of pediatric AML, its clinical utility in risk stratification remains limited. Here, we evaluate gene expression, pathogenic somatic mutations, and outcome in a cohort of 435 pediatric patients with a spectrum of pediatric myeloid-related acute leukemias for biological subtype discovery. This analysis revealed 63 patients with varying immunophenotypes that span a T-lineage and myeloid continuum designated as acute myeloid/T-lymphoblastic leukemia (AMTL). Within AMTL, two patient subgroups distinguished by FLT3-ITD and PRC2 mutations have different outcomes, demonstrating the impact of mutational composition on survival. Across the cohort, variability in outcomes of patients within isomutational subsets is influenced by transcriptional identity and the presence of a stem cell-like gene-expression signature. Integration of gene expression and somatic mutations leads to improved risk stratification. SIGNIFICANCE: Immunophenotype and somatic mutations play a significant role in treatment approach and risk stratification of acute leukemia. We conducted an integrated genomic analysis of pediatric myeloid malignancies and found that a combination of genetic and transcriptional readouts was superior to immunophenotype and genomic mutations in identifying biological subtypes and predicting outcomes. This article is highlighted in the In This Issue feature, p. 549.

A polygenic score for acute vaso-occlusive pain in pediatric sickle cell disease.

Individuals with monogenic disorders can experience variable phenotypes that are influenced by genetic variation. To investigate this in sickle cell disease (SCD), we performed whole-genome sequencing (WGS) of 722 individuals with hemoglobin HbSS or HbSß0-thalassemia from Baylor College of Medicine and from the St. Jude Children's Research Hospital Sickle Cell Clinical Research and Intervention Program (SCCRIP) longitudinal cohort study. We developed pipelines to identify genetic variants that modulate sickle hemoglobin polymerization in red blood cells and combined these with pain-associated variants to build a polygenic score (PGS) for acute vaso-occlusive pain (VOP). Overall, we interrogated the α-thalassemia deletion -α3.7 and 133 candidate single-nucleotide polymorphisms (SNPs) across 66 genes for associations with VOP in 327 SCCRIP participants followed longitudinally over 6 years. Twenty-one SNPs in 9 loci were associated with VOP, including 3 (BCL11A, MYB, and the ß-like globin gene cluster) that regulate erythrocyte fetal hemoglobin (HbF) levels and 6 (COMT, TBC1D1, KCNJ6, FAAH, NR3C1, and IL1A) that were associated previously with various pain syndromes. An unweighted PGS integrating all 21 SNPs was associated with the VOP event rate (estimate, 0.35; standard error, 0.04; P = 5.9 × 10-14) and VOP event occurrence (estimate, 0.42; standard error, 0.06; P = 4.1 × 10-13). These associations were stronger than those of any single locus. Our findings provide insights into the genetic modulation of VOP in children with SCD. More generally, we demonstrate the utility of WGS for investigating genetic contributions to the variable expression of SCD-associated morbidities.

Generation of mouse-human chimeric embryos.

Naive human pluripotent stem cells (hPSCs) can be used to generate mature human cells of all three germ layers in mouse-human chimeric embryos. Here, we describe a protocol for generating mouse-human chimeric embryos by injecting naive hPSCs converted from the primed state. Primed hPSCs are treated with a mammalian target of rapamycin inhibitor (Torin1) for 3 h and dissociated to single cells, which are plated on mouse embryonic fibroblasts in 2iLI medium, a condition essentially the same for culturing mouse embryonic stem cells. After 3-4 d, bright, dome-shaped colonies with mouse embryonic stem cell morphology are passaged in 2iLI medium. Established naive hPSCs are injected into mouse blastocysts, which produce E17.5 mouse embryos containing 0.1-4.0% human cells as quantified by next-generation sequencing of 18S ribosomal DNA amplicons. The protocol is suitable for studying the development of hPSCs in mouse embryos and may facilitate the generation of human cells, tissues and organs in animals.

Single cell transcriptomics reveals lineage trajectory of retinal ganglion cells in wild-type and Atoh7-null retinas.

Atoh7 has been believed to be essential for establishing the retinal ganglion cell (RGC) lineage, and Pou4f2 and Isl1 are known to regulate RGC specification and differentiation. Here we report our further study of the roles of these transcription factors. Using bulk RNA-seq, we identify genes regulated by the three transcription factors, which expand our understanding of the scope of downstream events. Using scRNA-seq on wild-type and mutant retinal cells, we reveal a transitional cell state of retinal progenitor cells (RPCs) co-marked by Atoh7 and other genes for different lineages and shared by all early retinal lineages. We further discover the unexpected emergence of the RGC lineage in the absence of Atoh7. We conclude that competence of RPCs for different retinal fates is defined by lineage-specific genes co-expressed in the transitional state and that Atoh7 defines the RGC competence and collaborates with other factors to shepherd transitional RPCs to the RGC lineage.

Developmental constraint shaped genome evolution and erythrocyte loss in Antarctic fishes following paleoclimate change.

In the frigid, oxygen-rich Southern Ocean (SO), Antarctic icefishes (Channichthyidae; Notothenioidei) evolved the ability to survive without producing erythrocytes and hemoglobin, the oxygen-transport system of virtually all vertebrates. Here, we integrate paleoclimate records with an extensive phylogenomic dataset of notothenioid fishes to understand the evolution of trait loss associated with climate change. In contrast to buoyancy adaptations in this clade, we find relaxed selection on the genetic regions controlling erythropoiesis evolved only after sustained cooling in the SO. This pattern is seen not only within icefishes but also occurred independently in other high-latitude notothenioids. We show that one species of the red-blooded dragonfish clade evolved a spherocytic anemia that phenocopies human patients with this disease via orthologous mutations. The genomic imprint of SO climate change is biased toward erythrocyte-associated conserved noncoding elements (CNEs) rather than to coding regions, which are largely preserved through pleiotropy. The drift in CNEs is specifically enriched near genes that are preferentially expressed late in erythropoiesis. Furthermore, we find that the hematopoietic marrow of icefish species retained proerythroblasts, which indicates that early erythroid development remains intact. Our results provide a framework for understanding the interactions between development and the genome in shaping the response of species to climate change.

Comparative genomic analysis of a Shiga toxin-producing Escherichia coli (STEC) O145:H25 associated with a severe pediatric case of hemolytic uremic syndrome in Davidson County, Tennessee, US.

BACKGROUND: Shiga toxin-producing E. coli (STECs) are foodborne pathogens associated with bloody diarrhea and hemolytic uremic syndrome (HUS). Although the STEC O157 serogroup accounts for the highest number of infections, HUS-related complications and deaths, the STEC non-O157, as a group, accounts for a larger proportion of STEC infections and lower HUS cases. There is limited information available on how to recognize non-O157 serotypes associated with severe disease. The objectives of this study were to describe a patient with STEC non-O157 infection complicated with HUS and to conduct a comparative whole genome sequence (WGS) analysis among the patient's STEC clinical isolate and STEC O157 and non-O157 strains. RESULTS: The STEC O145:H25 strain EN1I-0044-2 was isolated from a pediatric patient with diarrhea, HUS and severe neurologic and cardiorespiratory complications, who was enrolled in a previously reported case-control study of acute gastroenteritis conducted in Davidson County, Tennessee in 2013. The strain EN1I-0044-2 genome sequence contained a chromosome and three plasmids. Two of the plasmids were similar to those present in O145:H25 strains whereas the third unique plasmid EN1I-0044-2_03 shared no similarity with other STEC plasmids, and it carried 23 genes of unknown function. Strain EN1I-0044-2, compared with O145:H25 and O157 serogroup strains shared chromosome- and plasmid-encoded virulence factors, including Shiga toxin, LEE type III secretion system, LEE effectors, SFP fimbriae, and additional toxins and colonization factors. CONCLUSIONS: A STEC O145:H25 strain EN1I-0044-2 was isolated from a pediatric patient with severe disease, including HUS, in Davidson County, TN. Phylogenetic and comparison WGS analysis provided evidence that strain EN1I-0044-2 closely resembles O145:H25, and confirmed an independent evolutionary path of STEC O145:H25 and O145:H28 serotypes. The strain EN1I-0044-2 virulence make up was similar to other O145:H25 and O157 serogroups. It carried stx2 and the LEE pathogenicity island, and additional colonization factors and enterotoxin genes. A unique feature of strain EN1I-0044-2 was the presence of plasmid pEN1I-0044-2_03 carrying genes with functions to be determined. Further studies will be necessary to elucidate the role that newly acquired genes by O145:H25 strains play in pathogenesis, and to determine if they may serve as genetic markers of severe disease.

Pathogenic genetic variations of C. acnes are associated with clinically relevant orthopedic shoulder infections.

Many surgeons continue to face the clinical dilemma of interpreting a positive aspiration or unexpected positive Cutibacterium acnes (C. acnes) culture. There are factors that complicate the interpretation of positive cultures including variations in both frequency of false positive cultures and virulence properties. As indices of virulence, hemolytic strains, from previously confirmed clinically infected shoulders, were compared with non-hemolytic isolates determined to be contaminants, by RNA-sequencing (RNA-Seq). Six C. acnes isolates from patients who underwent revision total shoulder arthroplasty (TSA) were identified based on previously described infection criteria. Three C. acnes isolates from each group underwent RNA-Seq. Differential gene expression analysis, principal component analysis (PCA), and heatmap analysis were used to determine the gene variation and patterning between the definite infection and probable contaminant isolates. Differential gene expression analysis identified genes that were differentially expressed between the isolates classified as definite infection and isolates classified as probable contaminants. PCA using a 500 gene subset of identified genes was able to find combinations of these genes that separated out the definite infection and probable contaminants isolates. The heatmap demonstrated similar gene expression in the three Definite Infections isolates, and significantly different expression when compared with the probable contaminant isolates. Clinical significance: C. acnes revision TSA isolates classified as definite infection and probable contaminant demonstrated a similar gene expression pattern to each respective group and different gene expression pattern when compared between groups. These findings indicate distinct differences in C. acnes strains associated with clinically relevant orthopedic TSA infections.

Endothelialization of arterial vascular grafts by circulating monocytes.

Recently our group demonstrated that acellular tissue engineered vessels (A-TEVs) comprised of small intestinal submucosa (SIS) immobilized with heparin and vascular endothelial growth factor (VEGF) could be implanted into the arterial system of a pre-clinical ovine animal model, where they endothelialized within one month and remained patent. Here we report that immobilized VEGF captures blood circulating monocytes (MC) with high specificity under a range of shear stresses. Adherent MC differentiate into a mixed endothelial (EC) and macrophage (Mφ) phenotype and further develop into mature EC that align in the direction of flow and produce nitric oxide under high shear stress. In-vivo, newly recruited cells on the vascular lumen express MC markers and at later times they co-express MC and EC-specific proteins and maintain graft patency. This novel finding indicates that the highly prevalent circulating MC contribute directly to the endothelialization of acellular vascular grafts under the right chemical and biomechanical cues.

MYCN amplification and ATRX mutations are incompatible in neuroblastoma.

Aggressive cancers often have activating mutations in growth-controlling oncogenes and inactivating mutations in tumor-suppressor genes. In neuroblastoma, amplification of the MYCN oncogene and inactivation of the ATRX tumor-suppressor gene correlate with high-risk disease and poor prognosis. Here we show that ATRX mutations and MYCN amplification are mutually exclusive across all ages and stages in neuroblastoma. Using human cell lines and mouse models, we found that elevated MYCN expression and ATRX mutations are incompatible. Elevated MYCN levels promote metabolic reprogramming, mitochondrial dysfunction, reactive-oxygen species generation, and DNA-replicative stress. The combination of replicative stress caused by defects in the ATRX-histone chaperone complex, and that induced by MYCN-mediated metabolic reprogramming, leads to synthetic lethality. Therefore, ATRX and MYCN represent an unusual example, where inactivation of a tumor-suppressor gene and activation of an oncogene are incompatible. This synthetic lethality may eventually be exploited to improve outcomes for patients with high-risk neuroblastoma.

Global Genome Conformational Programming during Neuronal Development Is Associated with CTCF and Nuclear FGFR1-The Genome Archipelago Model.

During the development of mouse embryonic stem cells (ESC) to neuronal committed cells (NCC), coordinated changes in the expression of 2851 genes take place, mediated by the nuclear form of FGFR1. In this paper, widespread differences are demonstrated in the ESC and NCC inter- and intra-chromosomal interactions, chromatin looping, the formation of CTCF- and nFGFR1-linked Topologically Associating Domains (TADs) on a genome-wide scale and in exemplary HoxA-D loci. The analysis centered on HoxA cluster shows that blocking FGFR1 disrupts the loop formation. FGFR1 binding and genome locales are predictive of the genome interactions; likewise, chromatin interactions along with nFGFR1 binding are predictive of the genome function and correlate with genome regulatory attributes and gene expression. This study advances a topologically integrated genome archipelago model that undergoes structural transformations through the formation of nFGFR1-associated TADs. The makeover of the TAD islands serves to recruit distinct ontogenic programs during the development of the ESC to NCC.

Thermotolerance in the pathogen Cryptococcus neoformans is linked to antigen masking via mRNA decay-dependent reprogramming.

A common feature shared by systemic fungal pathogens of environmental origin, such as Cryptococcus neoformans, is their ability to adapt to mammalian core body temperature. In C. neoformans, this adaptation is accompanied by Ccr4-mediated decay of ribosomal protein mRNAs. Here we use the related, but thermo-intolerant species Cryptococcus amylolentus to demonstrate that this response contributes to host-temperature adaptation and pathogenicity of cryptococci. In a C. neoformans ccr4Δ mutant, stabilized ribosomal protein mRNAs are retained in the translating pool, and stress-induced transcriptomic changes are reduced in comparison with the wild type strain, likely due to ineffective translation of transcription factors. In addition, the mutant displays increased exposure of cell wall glucans, and recognition by Dectin-1 results in increased phagocytosis by lung macrophages, linking mRNA decay to adaptation and immune evasion.

Correction to: H3.3 K27M depletion increases differentiation and extends latency of diffuse intrinsic pontine glioma growth in vivo.

The original article can be found online.

H3.3 K27M depletion increases differentiation and extends latency of diffuse intrinsic pontine glioma growth in vivo.

Histone H3 K27M mutation is the defining molecular feature of the devastating pediatric brain tumor, diffuse intrinsic pontine glioma (DIPG). The prevalence of histone H3 K27M mutations indicates a critical role in DIPGs, but the contribution of the mutation to disease pathogenesis remains unclear. We show that knockdown of this mutation in DIPG xenografts restores K27M-dependent loss of H3K27me3 and delays tumor growth. Comparisons of matched DIPG xenografts with and without K27M knockdown allowed identification of mutation-specific effects on the transcriptome and epigenome. The resulting transcriptional changes recapitulate expression signatures from K27M primary DIPG tumors and are strongly enriched for genes associated with nervous system development. Integrated analysis of ChIP-seq and expression data showed that genes upregulated by the mutation are overrepresented in apparently bivalent promoters. Many of these targets are associated with more immature differentiation states. Expression profiles indicate K27M knockdown decreases proliferation and increases differentiation within lineages represented in DIPG. These data suggest that K27M-mediated loss of H3K27me3 directly regulates a subset of genes by releasing poised promoters, and contributes to tumor phenotype and growth by limiting differentiation. The delayed tumor growth associated with knockdown of H3 K27M provides evidence that this highly recurrent mutation is a relevant therapeutic target.

The Clonal Evolution of Metastatic Osteosarcoma as Shaped by Cisplatin Treatment.

To investigate the genomic evolution of metastatic pediatric osteosarcoma, we performed whole-genome and targeted deep sequencing on 14 osteosarcoma metastases and two primary tumors from four patients (two to eight samples per patient). All four patients harbored ancestral (truncal) somatic variants resulting in TP53 inactivation and cell-cycle aberrations, followed by divergence into relapse-specific lineages exhibiting a cisplatin-induced mutation signature. In three of the four patients, the cisplatin signature accounted for >40% of mutations detected in the metastatic samples. Mutations potentially acquired during cisplatin treatment included NF1 missense mutations of uncertain significance in two patients and a KIT G565R activating mutation in one patient. Three of four patients demonstrated widespread ploidy differences between samples from the sample patient. Single-cell seeding of metastasis was detected in most metastatic samples. Cross-seeding between metastatic sites was observed in one patient, whereas in another patient a minor clone from the primary tumor seeded both metastases analyzed. These results reveal extensive clonal heterogeneity in metastatic osteosarcoma, much of which is likely cisplatin-induced. IMPLICATIONS: The extent and consequences of chemotherapy-induced damage in pediatric cancers is unknown. We found that cisplatin treatment can potentially double the mutational burden in osteosarcoma, which has implications for optimizing therapy for recurrent, chemotherapy-resistant disease.

Clinical cancer genomic profiling by three-platform sequencing of whole genome, whole exome and transcriptome.

To evaluate the potential of an integrated clinical test to detect diverse classes of somatic and germline mutations relevant to pediatric oncology, we performed three-platform whole-genome (WGS), whole exome (WES) and transcriptome (RNA-Seq) sequencing of tumors and normal tissue from 78 pediatric cancer patients in a CLIA-certified, CAP-accredited laboratory. Our analysis pipeline achieves high accuracy by cross-validating variants between sequencing types, thereby removing the need for confirmatory testing, and facilitates comprehensive reporting in a clinically-relevant timeframe. Three-platform sequencing has a positive predictive value of 97-99, 99, and 91% for somatic SNVs, indels and structural variations, respectively, based on independent experimental verification of 15,225 variants. We report 240 pathogenic variants across all cases, including 84 of 86 known from previous diagnostic testing (98% sensitivity). Combined WES and RNA-Seq, the current standard for precision oncology, achieved only 78% sensitivity. These results emphasize the critical need for incorporating WGS in pediatric oncology testing.

Identification of Therapeutic Targets in Rhabdomyosarcoma through Integrated Genomic, Epigenomic, and Proteomic Analyses.

Personalized cancer therapy targeting somatic mutations in patient tumors is increasingly being incorporated into practice. Other therapeutic vulnerabilities resulting from changes in gene expression due to tumor specific epigenetic perturbations are progressively being recognized. These genomic and epigenomic changes are ultimately manifest in the tumor proteome and phosphoproteome. We integrated transcriptomic, epigenomic, and proteomic/phosphoproteomic data to elucidate the cellular origins and therapeutic vulnerabilities of rhabdomyosarcoma (RMS). We discovered that alveolar RMS occurs further along the developmental program than embryonal RMS. We also identified deregulation of the RAS/MEK/ERK/CDK4/6, G2/M, and unfolded protein response pathways through our integrated analysis. Comprehensive preclinical testing revealed that targeting the WEE1 kinase in the G2/M pathway is the most effective approach in vivo for high-risk RMS.

Genetic Risk for Subsequent Neoplasms Among Long-Term Survivors of Childhood Cancer.

Purpose Childhood cancer survivors are at increased risk of subsequent neoplasms (SNs), but the germline genetic contribution is largely unknown. We assessed the contribution of pathogenic/likely pathogenic (P/LP) mutations in cancer predisposition genes to their SN risk. Patients and Methods Whole-genome sequencing (30-fold) was performed on samples from childhood cancer survivors who were ≥ 5 years since initial cancer diagnosis and participants in the St Jude Lifetime Cohort Study, a retrospective hospital-based study with prospective clinical follow-up. Germline mutations in 60 genes known to be associated with autosomal dominant cancer predisposition syndromes with moderate to high penetrance were classified by their pathogenicity according to the American College of Medical Genetics and Genomics guidelines. Relative rates (RRs) and 95% CIs of SN occurrence by mutation status were estimated using multivariable piecewise exponential regression stratified by radiation exposure. Results Participants were 3,006 survivors (53% male; median age, 35.8 years [range, 7.1 to 69.8 years]; 56% received radiotherapy), 1,120 SNs were diagnosed among 439 survivors (14.6%), and 175 P/LP mutations were identified in 5.8% (95% CI, 5.0% to 6.7%) of survivors. Mutations were associated with significantly increased rates of breast cancer (RR, 13.9; 95% CI, 6.0 to 32.2) and sarcoma (RR, 10.6; 95% CI, 4.3 to 26.3) among irradiated survivors and with increased rates of developing any SN (RR, 4.7; 95% CI, 2.4 to 9.3), breast cancer (RR, 7.7; 95% CI, 2.4 to 24.4), nonmelanoma skin cancer (RR, 11.0; 95% CI, 2.9 to 41.4), and two or more histologically distinct SNs (RR, 18.6; 95% CI, 3.5 to 99.3) among nonirradiated survivors. Conclusion The findings support referral of all survivors for genetic counseling for potential clinical genetic testing, which should be prioritized for nonirradiated survivors with any SN and for those with breast cancer or sarcoma in the field of prior irradiation.