Single Nucleus RNA Sequencing of Remnant Kidney Biopsies and Urine Cell RNA Sequencing Reveal Cell Specific Markers of Covid-19 Acute Kidney Injury.

Acute kidney injury (AKI) in COVID-19 patients is associated with high mortality and morbidity. Critically ill COVID-19 patients are at twice the risk of in-hospital mortality compared to non-COVID AKI patients. We know little about the cell-specific mechanism in the kidney that contributes to worse clinical outcomes in these patients. New generation single cell technologies have the potential to provide insights into physiological states and molecular mechanisms in COVID-AKI. One of the key limitations is that these patients are severely ill posing significant risks in procuring additional biopsy tissue. We recently generated single nucleus RNA-sequencing data using COVID-AKI patient biopsy tissue as part of the human kidney atlas. Here we describe this approach in detail and report deeper comparative analysis of snRNAseq of 4 COVID-AKI, 4 reference, and 6 non-COVID-AKI biopsies. We also generated and analyzed urine transcriptomics data to find overlapping COVID-AKI-enriched genes and their corresponding cell types in the kidney from snRNA-seq data. We identified all major and minor cell types and states by using by using less than a few cubic millimeters of leftover tissue after pathological workup in our approach. Differential expression analysis of COVID-AKI biopsies showed pathways enriched in viral response, WNT signaling, kidney development, and cytokines in several nephron epithelial cells. COVID-AKI profiles showed a much higher proportion of altered TAL cells than non-COVID AKI and the reference samples. In addition to kidney injury and fibrosis markers indicating robust remodeling we found that, 17 genes overlap between urine cell COVID-AKI transcriptome and the snRNA-seq data from COVID-AKI biopsies. A key feature was that several of the distal nephron and collecting system cell types express these markers. Some of these markers have been previously observed in COVID-19 studies suggesting a common mechanism of injury and potentially the kidney as one of the sources of soluble factors with a potential role in disease progression. Translational Statement: The manuscript describes innovation, application and discovery that impact clinical care in kidney disease. First, the approach to maximize use of remnant frozen clinical biopsies to inform on clinically relevant molecular features can augment existing pathological workflow for any frozen tissue without much change in the protocol. Second, this approach is transformational in medical crises such as pandemics where mechanistic insights are needed to evaluate organ injury, targets for drug therapy and diagnostic and prognostic markers. Third, the cell type specific and soluble markers identified and validated can be used for diagnoses or prognoses in AKI due to different etiologies and in multiorgan injury.

Evaluation of Plasma Biomarkers to Predict Major Adverse Kidney Events in Hospitalized Patients With COVID-19.

RATIONALE & OBJECTIVE: Patients hospitalized with COVID-19 are at increased risk for major adverse kidney events (MAKE). We sought to identify plasma biomarkers predictive of MAKE in patients hospitalized with COVID-19. STUDY DESIGN: Prospective cohort study. SETTING & PARTICIPANTS: A total of 576 patients hospitalized with COVID-19 between March 2020 and January 2021 across 3 academic medical centers. EXPOSURE: Twenty-six plasma biomarkers of injury, inflammation, and repair from first available blood samples collected during hospitalization. OUTCOME: MAKE, defined as KDIGO stage 3 acute kidney injury (AKI), dialysis-requiring AKI, or mortality up to 60 days. ANALYTICAL APPROACH: Cox proportional hazards regression to associate biomarker level with MAKE. We additionally applied the least absolute shrinkage and selection operator (LASSO) and random forest regression for prediction modeling and estimated model discrimination with time-varying C index. RESULTS: The median length of stay for COVID-19 hospitalization was 9 (IQR, 5-16) days. In total, 95 patients (16%) experienced MAKE. Each 1 SD increase in soluble tumor necrosis factor receptor 1 (sTNFR1) and sTNFR2 was significantly associated with an increased risk of MAKE (adjusted HR [AHR], 2.30 [95% CI, 1.86-2.85], and AHR, 2.26 [95% CI, 1.73-2.95], respectively). The C index of sTNFR1 alone was 0.80 (95% CI, 0.78-0.84), and the C index of sTNFR2 was 0.81 (95% CI, 0.77-0.84). LASSO and random forest regression modeling using all biomarkers yielded C indexes of 0.86 (95% CI, 0.83-0.89) and 0.84 (95% CI, 0.78-0.91), respectively. LIMITATIONS: No control group of hospitalized patients without COVID-19. CONCLUSIONS: We found that sTNFR1 and sTNFR2 are independently associated with MAKE in patients hospitalized with COVID-19 and can both also serve as predictors for adverse kidney outcomes. PLAIN-LANGUAGE SUMMARY: Patients hospitalized with COVID-19 are at increased risk for long-term adverse health outcomes, but not all patients suffer long-term kidney dysfunction. Identification of patients with COVID-19 who are at high risk for adverse kidney events may have important implications in terms of nephrology follow-up and patient counseling. In this study, we found that the plasma biomarkers soluble tumor necrosis factor receptor 1 (sTNFR1) and sTNFR2 measured in hospitalized patients with COVID-19 were associated with a greater risk of adverse kidney outcomes. Along with clinical variables previously shown to predict adverse kidney events in patients with COVID-19, both sTNFR1 and sTNFR2 are also strong predictors of adverse kidney outcomes.

A microfluidic-based filtration system to enrich for bone marrow disseminated tumor cells from breast cancer patients.

Disseminated tumors cells (DTCs) present in the bone marrow (BM) are believed to be the progenitors of distant metastatic spread, a major cause of mortality in breast cancer patients. To better understand the behavior and therapeutic vulnerabilities of these rare cell populations, unbiased methods for selective cell enrichment are required. In this study, we have evaluated a microfluidic-based filtration system (ParsortixR, Angle PLC), previously demonstrated for use in circulating tumor cell (CTC) capture, to capture BM DTCs. Performance using BM samples was also compared directly to enrichment of CTCs in the peripheral blood (PB) from both metastatic and non-metastatic breast cancer patients. Although the non-specific capture of BM immune cells was significant, the device could routinely achieve significant cytoreduction of BM and PB WBCs and at least 1,000-fold enrichment of DTCs, based on labeled tumor cell spike-in experiments. Detection of previously characterized DTC-associated gene expression biomarkers was greatly enhanced by the enrichment method, as demonstrated by droplet digital PCR assay. Cells eluted from the device were viable and suitable for single cell RNA sequencing experiments. DTCs in enriched BM samples comprised up to 5% of the total cell population, allowing for effective single cell and population-based transcriptional profiling of these rare cells. Use of the Parsortix instrument will be an effective approach to enrich for rare BM DTCs in order to better understand their diverse molecular phenotypes and develop approaches to eradicate these cells to prevent distant disease development in breast cancer patients.

Whole exome sequencing reveals the maintained polyclonal nature from primary to metastatic malignant peripheral nerve sheath tumor in two patients with NF1.

BACKGROUND: Malignant peripheral nerve sheath tumors (MPNSTs) are aggressive soft tissue sarcomas with high metastatic rates and poor overall patient survival. There are currently no effective therapies, underscoring the pressing need to define the molecular etiologies that underlie MPNST progression. The aim of this study was to examine clonal progression and identify the molecular events critical for MPNST spread. METHODS: In two patients with temporally and spatially distinct metastatic lesions, we employed whole exome sequencing (WES) to elucidate the genetic events of clonal progression, thus identifying the molecular events critical for MPNST spread. RESULTS: First, we demonstrated shared clonal origins for the metastatic lesions relative to the primary tumors, which were maintained throughout the course of MPNST progression, supporting the conclusion that cancer cells with metastatic potential already exist in the primary neoplasm. Second, we discovered TRIM23, a member of the Tripartite Motif family of proteins, as a regulator of MPNST lung metastatic spread in vivo. CONCLUSIONS: The ability to track the genomic evolution from primary to metastatic MPNST offers new insights into the sequence of genetic events required for tumor progression and has identified TRIM23 as a novel target for future study in this rare cancer.

Tumor mutation burden and checkpoint immunotherapy markers in primary and metastatic synovial sarcoma.

Synovial sarcoma (SS) is a soft-tissue malignancy that most often affects patients aged between 15 and 40 years, and the prognosis for patients with metastatic disease is generally poor. This study was performed to evaluate checkpoint blockade immunotherapy markers in SS, including tumor mutational burden (TMB), DNA mismatch repair (MMR) status, and PDL-1 (programmed cell death ligand 1), PD1 (programmed cell death 1), and CD8 expression by normal-tumor paired whole-exome sequencing (WES) and immunohistochemistry (IHC). Outcomes evaluated included event-free and overall survival. Twenty one (21) FISH (Fluorescence In Situ Hybridization)-confirmed SS cases (11 F, 10 M) were studied, with age ranging from 8 to 89 years at diagnosis and follow-up ranging from 1 to 16 years. TMB (n = 16) ranged from 0.83 to 212/Mb (median, 1.7). Only one case showed a high TMB of 212/Mb and missense variants of MMR genes in the primary tumor, while the other 15 cases had a low TMB of less than 5/Mb. IHC was performed on all 21 tumor samples for PD-L1, PD1, CD8, and MMR proteins. PD-L1 membranous staining was detected in 3 of 21 cases (14.3%), ranging from 1 to 5% for tumor proportion score and 1-10 for combined positive score. PD1 was detected in 15 of 21 cases (71.4%), ranging from 1 to 25/HPF (high power field) (median, 2). CD8 stain was seen in all cases, ranging from 2 to 60/HPF (median, 5). PD1 staining results correlated with CD8 staining results (P < 0.0001). No correlation of TMB or IHC markers was found with survival. No fixed pattern of TMB or IHCs between primary and metastatic tumors was observed; there was no correlation between TMB or IHCs and age, location, or diagnosis subtype. All of the cases tested showed retained expression of MMR proteins. The results show that for SS, a tumor with strong driver translocation, most cases have a low TMB, but occasionally a high TMB may be present, as observed in 1 of the 16 (6.25%) cases. The demonstration of a subgroup of SS cases with high TMB might explain the 10% response rate to checkpoint immunotherapy observed in clinical trials in patients with SS.

Tumor mutation burden, DNA mismatch repair status and checkpoint immunotherapy markers in primary and relapsed malignant rhabdoid tumors.

INTRODUCTION: Malignant rhabdoid tumor (MRT) is a rare, aggressive pediatric tumor of nuclear lineage. It is mainly characterized by germline or somatic SMARCB1 (INI1) driver mutations. To characterize the potential for immunotherapy in untreated and treated MRT, current study investigated tumor mutational burden (TMB) and other biomarkers in MRT. MATERIAL AND METHODS: Normal-tumor paired whole exome sequencing (WES) and/or immunohistochemistry (IHC) of DNA mismatch repair (MMR) proteins, PD-L1, PD-1 and CD8 were performed in 16 cases, some with both primary and relapsed tumor. RESULTS: Five cases subjected to WES demonstrated germline SMARCB1 (INI1) mutations. TMB was 0.7-1.07/Mb in 4 of the 5 primary untreated tumors, and 33.81/Mb in one case with pathogenic MMR, POLD, and POLE mutations. Ten cases tested for MMR status by IHC showed retained nuclear expression of the proteins. Eight of the 16 cases (8/16, 50%) showed membranous expression of PD-L1 in 10-70% of tumor cells (tumor proportion score, TPS). Nine cases (9/16, 56.3%) showed high (>2/HPF) tumor infiltrating lymphocytes with PD-1 staining ranging 10-60%, correlating with tumor PD-L1 staining (p < 0.0001). Between post-treatment metastatic tumors and the pre-treatment primary tumors, TMB was similar while PD-L1 TPS was similar or lower. CONCLUSION: MRT has a low TMB. Nonetheless, because a subset of MRT cases have a PD-L1 TPS greater than the cutoff for checkpoint therapy in other malignancies, the utility of immune checkpoint inhibitors should be studied in this patient population.

Exome sequencing identifies gene variants and networks associated with extreme respiratory outcomes following preterm birth.

BACKGROUND: Previous studies have identified genetic variants associated with bronchopulmonary dysplasia (BPD) in extremely preterm infants. However, findings with genome-wide significance have been rare, and not replicated. We hypothesized that whole exome sequencing (WES) of premature subjects with extremely divergent phenotypic outcomes could facilitate the identification of genetic variants or gene networks contributing disease risk. RESULTS: The Prematurity and Respiratory Outcomes Program (PROP) recruited a cohort of > 765 extremely preterm infants for the identification of markers of respiratory morbidity. We completed WES on 146 PROP subjects (85 affected, 61 unaffected) representing extreme phenotypes of early respiratory morbidity. We tested for association between disease status and individual common variants, screened for rare variants exclusive to either affected or unaffected subjects, and tested the combined association of variants across gene loci. Pathway analysis was performed and disease-related expression patterns were assessed. Marginal association with BPD was observed for numerous common and rare variants. We identified 345 genes with variants unique to BPD-affected preterm subjects, and 292 genes with variants unique to our unaffected preterm subjects. Of these unique variants, 28 (19 in the affected cohort and 9 in unaffected cohort) replicate a prior WES study of BPD-associated variants. Pathway analysis of sets of variants, informed by disease-related gene expression, implicated protein kinase A, MAPK and Neuregulin/epidermal growth factor receptor signaling. CONCLUSIONS: We identified novel genes and associated pathways that may play an important role in susceptibility/resilience for the development of lung disease in preterm infants.

Clinical genomic profiling identifies TYK2 mutation and overexpression in patients with neurofibromatosis type 1-associated malignant peripheral nerve sheath tumors.

BACKGROUND: Malignant peripheral nerve sheath tumors (MPNSTs) are aggressive sarcomas that arise at an estimated frequency of 8% to 13% in individuals with neurofibromatosis type 1 (NF1). Compared with their sporadic counterparts, NF1-associated MPNSTs (NF1-MPNSTs) develop in young adults, frequently recur (approximately 50% of cases), and carry a dismal prognosis. As such, most individuals affected with NF1-MPNSTs die within 5 years of diagnosis, despite surgical resection combined with radiotherapy and chemotherapy. METHODS: Clinical genomic profiling was performed using 1000 ng of DNA from 7 cases of NF1-MPNST, and bioinformatic analyses were conducted to identify genes with actionable mutations. RESULTS: A total of 3 women and 4 men with NF1-MPNST were identified (median age, 38 years). Nonsynonymous mutations were discovered in 4 genes (neurofibromatosis type 1 [NF1], ROS proto-oncogene 1 [ROS1], tumor protein p53 [TP53], and tyrosine kinase 2 [TYK2]), which in addition were mutated in other MPNST cases in this sample set. Consistent with their occurrence in individuals with NF1, all tumors had at least 1 mutation in the NF1 gene. Whereas TP53 gene mutations are frequently observed in other cancers, ROS1 mutations are common in melanoma (15%-35%), another neural crest-derived malignancy. In contrast, TYK2 mutations are uncommon in other malignancies (<7%). In the current series, recurrent TYK2 mutations were identified in 2 cases of NF1-MPNST (30% of cases), whereas TYK2 protein overexpression was observed in 60% of MPNST cases using an independently generated tissue microarray, regardless of NF1 status. CONCLUSIONS: Clinical genomic analysis of the current series of NF1-MPNST cases found that TYK2 is a new gene mutated in MPNST. Future work will focus on examining the utility of TYK2 expression as a biomarker and therapeutic target for these cancers. Cancer 2017;123:1194-1201. © 2016 American Cancer Society.

Genome-wide significance testing of variation from single case exomes.

Standard techniques from genetic epidemiology are ill-suited to formally assess the significance of variants identified from a single case. We developed a statistical inference framework for identifying unusual functional variation from a single exome or genome, what we refer to as the 'n-of-one' problem. Using this approach we assessed our ability to identify the causal genotypes in over 5 million simulated cases of Mendelian disease, identifying 39% of disease genotypes as the most damaging unit in a typical exome background. We applied our approach to 129 n-of-one families from the Undiagnosed Diseases Program, nominating 60% of 30 disease genes determined to be diagnostic by a standard clinical workup. Our method can currently produce well-calibrated P values when applied to single genomes, can facilitate integration of multiple data types for n-of-one analyses, and, with further work, could become a widely used epidemiological method like linkage analysis or genome-wide association analysis.