Identity and nature of neural stem cells in the adult human subventricular zone.

The existence of neural stem cells (NSCs) in adult human brain neurogenic regions remains unresolved. To address this, we created a cell atlas of the adult human subventricular zone (SVZ) derived from fresh neurosurgical samples using single-cell transcriptomics. We discovered 2 adult radial glia (RG)-like populations, aRG1 and aRG2. aRG1 shared features with fetal early RG (eRG) and aRG2 were transcriptomically similar to fetal outer RG (oRG). We also captured early neuronal and oligodendrocytic NSC states. We found that the biological programs driven by their transcriptomes support their roles as early lineage NSCs. Finally, we show that these NSCs have the potential to transition between states and along lineage trajectories. These data reveal that multipotent NSCs reside in the adult human SVZ.

Invasive growth of brain metastases is linked to CHI3L1 release from pSTAT3-positive astrocytes.

BACKGROUND: Compared to minimally invasive brain metastases (MI BrM), highly invasive (HI) lesions form abundant contacts with cells in the peritumoral brain parenchyma and are associated with poor prognosis. Reactive astrocytes (RAs) labeled by phosphorylated STAT3 (pSTAT3) have recently emerged as a promising therapeutic target for BrM. Here, we explore whether BrM invasion pattern is influenced by pSTAT3+ RAs and may serve as a predictive biomarker for STAT3 inhibition. METHODS: We used immunohistochemistry to identify pSTAT3+ RAs in HI and MI human and patient-derived xenograft (PDX) BrM. Using PDX, syngeneic, and transgenic mouse models of HI and MI BrM, we assessed how pharmacological STAT3 inhibition or RA-specific STAT3 genetic ablation affected BrM growth in vivo. Cancer cell invasion was modeled in vitro using a brain slice-tumor co-culture assay. We performed single-cell RNA sequencing of human BrM and adjacent brain tissue. RESULTS: RAs expressing pSTAT3 are situated at the brain-tumor interface and drive BrM invasive growth. HI BrM invasion pattern was associated with delayed growth in the context of STAT3 inhibition or genetic ablation. We demonstrate that pSTAT3+ RAs secrete Chitinase 3-like-1 (CHI3L1), which is a known STAT3 transcriptional target. Furthermore, single-cell RNA sequencing identified CHI3L1-expressing RAs in human HI BrM. STAT3 activation, or recombinant CHI3L1 alone, induced cancer cell invasion into the brain parenchyma using a brain slice-tumor plug co-culture assay. CONCLUSIONS: Together, these data reveal that pSTAT3+ RA-derived CHI3L1 is associated with BrM invasion, implicating STAT3 and CHI3L1 as clinically relevant therapeutic targets for the treatment of HI BrM.

Infantile onset encephalomyopathy, retinopathy, optic atrophy, and mitochondrial DNA depletion associated with a novel pathogenic DHX16 variant.

We studied a patient with mitochondrial DNA depletion in skeletal muscle and a multiorgan phenotype, including fatal encephalomyopathy, retinopathy, optic atrophy, and sensorineural hearing loss. Instead of pathogenic variants in the mitochondrial maintenance genes, we identified previously unpublished variant in DHX16 gene, a de novo heterozygous c.1360C>T (p. Arg454Trp). Variants in DHX16 encoding for DEAH-box RNA helicase have previously been reported only in five patients with a phenotype called as neuromuscular oculoauditory syndrome including developmental delay, neuromuscular symptoms, and ocular or auditory defects with or without seizures. We performed functional studies on patient-derived fibroblasts and skeletal muscle revealing, that the DHX16 expression was decreased. Clinical features together with functional data suggest, that our patient's disease is associated with a novel pathogenic DHX16 variant, and mtDNA depletion could be a secondary manifestation of the disease.

Germline Missense Variants in CDC20 Result in Aberrant Mitotic Progression and Familial Cancer.

CDC20 is a coactivator of the anaphase promoting complex/cyclosome (APC/C) and is essential for mitotic progression. APC/CCDC20 is inhibited by the spindle assembly checkpoint (SAC), which prevents premature separation of sister chromatids and aneuploidy in daughter cells. Although overexpression of CDC20 is common in many cancers, oncogenic mutations have never been identified in humans. Using whole-exome sequencing, we identified heterozygous missense CDC20 variants (L151R and N331K) that segregate with ovarian germ cell tumors in two families. Functional characterization showed these mutants retain APC/C activation activity but have impaired binding to BUBR1, a component of the SAC. Expression of L151R and N331K variants promoted mitotic slippage in HeLa cells and primary skin fibroblasts derived from carriers. Generation of mice carrying the N331K variant using CRISPR-Cas9 showed that, although homozygous N331K mice were nonviable, heterozygotes displayed accelerated oncogenicity of Myc-driven cancers. These findings highlight an unappreciated role for CDC20 variants as tumor-promoting genes. SIGNIFICANCE: Two germline CDC20 missense variants that segregate with cancer in two families compromise the spindle assembly checkpoint and lead to aberrant mitotic progression, which could predispose cells to transformation. See related commentary by Villarroya-Beltri and Malumbres, p. 3432.

Glioblastoma scRNA-seq shows treatment-induced, immune-dependent increase in mesenchymal cancer cells and structural variants in distal neural stem cells.

BACKGROUND: Glioblastoma is a treatment-resistant brain cancer. Its hierarchical cellular nature and its tumor microenvironment (TME) before, during, and after treatments remain unresolved. METHODS: Here, we used single-cell RNA sequencing to analyze new and recurrent glioblastoma and the nearby subventricular zone (SVZ). RESULTS: We found 4 glioblastoma neural lineages are present in new and recurrent glioblastoma with an enrichment of the cancer mesenchymal lineage, immune cells, and reactive astrocytes in early recurrences. Cancer lineages were hierarchically organized around cycling oligodendrocytic and astrocytic progenitors that are transcriptomically similar but distinct to SVZ neural stem cells (NSCs). Furthermore, NSCs from the SVZ of patients with glioblastoma harbored glioblastoma chromosomal anomalies. Lastly, mesenchymal cancer cells and TME reactive astrocytes shared similar gene signatures which were induced by radiotherapy in a myeloid-dependent fashion in vivo. CONCLUSION: These data reveal the dynamic, immune-dependent nature of glioblastoma's response to treatments and identify distant NSCs as likely cells of origin.

PRAME protein expression in DICER1-related tumours.

DICER1 syndrome is an autosomal dominant tumour predisposition syndrome usually affecting persons under 30 years of age. Many of the associated benign and malignant lesions occur almost exclusively in DICER1 syndrome. One such tumour, pituitary blastoma (pitB), overexpresses PRAME 500x above control levels. PRAME (PReferentially expressed Antigen in MElanoma) is expressed in malignancies that are not DICER1-related (e.g. melanoma). To address whether PRAME expression is part of the DICER1 phenotype, or simply a feature of pitB, a series of 75 DICER1-mutated specimens and 33 non-mutated specimens was surveyed using immunohistochemistry for PRAME, together with EZH2, which complexes with PRAME. In DICER1-mutated specimens, positive staining for PRAME was only seen in malignant tumours; 7 of 11 histological types and 34/62 individual tumours were positive, while non-tumourous lesions were always negative. Pleuropulmonary blastoma (PPB) showed a continuum in staining, with type I lesions being PRAME negative (n = 7) but all type II and type III lesions PRAME positive (n = 7). Similarly, cystic nephroma (CN) was negative (n = 8), with anaplastic sarcoma of the kidney being positive (n = 2). However, one atypical CN with mesenchymal cell proliferation was PRAME-positive. Embryonal rhabdomyosarcoma (RMS) with DICER1 pathogenic variants (PVs) was positive for PRAME (5/6), but the same tumour type without DICER1 PVs was also positive (9/15). Staining for EZH2 corresponded to that seen with PRAME, validating the latter. This study leads us to conclude that (1) PRAME expression occurs in two-thirds of DICER1-related malignancies; (2) PRAME may be a marker for the progression that certain DICER1-related lesions are thought to undergo, such as PPB and CN; and (3) PRAME expression in some tumours, such as RMS, appears to be an intrinsic feature of the tumour, rather than specifically related to DICER1 PVs. Therapy directed against PRAME may offer novel treatment options in patients with the DICER1 syndrome.

A functionally impaired missense variant identified in French Canadian families implicates FANCI as a candidate ovarian cancer-predisposing gene.

BACKGROUND: Familial ovarian cancer (OC) cases not harbouring pathogenic variants in either of the BRCA1 and BRCA2 OC-predisposing genes, which function in homologous recombination (HR) of DNA, could involve pathogenic variants in other DNA repair pathway genes. METHODS: Whole exome sequencing was used to identify rare variants in HR genes in a BRCA1 and BRCA2 pathogenic variant negative OC family of French Canadian (FC) ancestry, a population exhibiting genetic drift. OC cases and cancer-free individuals from FC and non-FC populations were investigated for carrier frequency of FANCI c.1813C>T; p.L605F, the top-ranking candidate. Gene and protein expression were investigated in cancer cell lines and tissue microarrays, respectively. RESULTS: In FC subjects, c.1813C>T was more common in familial (7.1%, 3/42) than sporadic (1.6%, 7/439) OC cases (P = 0.048). Carriers were detected in 2.5% (74/2950) of cancer-free females though female/male carriers were more likely to have a first-degree relative with OC (121/5249, 2.3%; Spearman correlation = 0.037; P = 0.011), suggesting a role in risk. Many of the cancer-free females had host factors known to reduce risk to OC which could influence cancer risk in this population. There was an increased carrier frequency of FANCI c.1813C>T in BRCA1 and BRCA2 pathogenic variant negative OC families, when including the discovery family, compared to cancer-free females (3/23, 13%; OR = 5.8; 95%CI = 1.7-19; P = 0.005). In non-FC subjects, 10 candidate FANCI variants were identified in 4.1% (21/516) of Australian OC cases negative for pathogenic variants in BRCA1 and BRCA2, including 10 carriers of FANCI c.1813C>T. Candidate variants were significantly more common in familial OC than in sporadic OC (P = 0.04). Localization of FANCD2, part of the FANCI-FANCD2 (ID2) binding complex in the Fanconi anaemia (FA) pathway, to sites of induced DNA damage was severely impeded in cells expressing the p.L605F isoform. This isoform was expressed at a reduced level, destabilized by DNA damaging agent treatment in both HeLa and OC cell lines, and exhibited sensitivity to cisplatin but not to a poly (ADP-ribose) polymerase inhibitor. By tissue microarray analyses, FANCI protein was consistently expressed in fallopian tube epithelial cells and only expressed at low-to-moderate levels in 88% (83/94) of OC samples. CONCLUSIONS: This is the first study to describe candidate OC variants in FANCI, a member of the ID2 complex of the FA DNA repair pathway. Our data suggest that pathogenic FANCI variants may modify OC risk in cancer families.

Molecular characterization of DICER1-mutated pituitary blastoma.

Pituitary blastoma (PitB) has recently been identified as a rare and potentially lethal pediatric intracranial tumor. All cases that have been studied molecularly possess at least one DICER1 pathogenic variant. Here, we characterized nine pituitary samples, including three fresh frozen PitBs, three normal fetal pituitary glands and three normal postnatal pituitary glands using small-RNA-Seq, RNA-Seq, methylation profiling, whole genome sequencing and Nanostring® miRNA analyses; an extended series of 21 pituitary samples was used for validation purposes. These analyses demonstrated that DICER1 RNase IIIb hotspot mutations in PitBs induced improper processing of miRNA precursors, resulting in aberrant 5p-derived miRNA products and a skewed distribution of miRNAs favoring mature 3p over 5p miRNAs. This led to dysregulation of hundreds of 5p and 3p miRNAs and concomitant dysregulation of numerous mRNA targets. Gene expression analysis revealed PRAME as the most significantly upregulated gene (500-fold increase). PRAME is a member of the Retinoic Acid Receptor (RAR) signaling pathway and in PitBs, the RAR, WNT and NOTCH pathways are dysregulated. Cancer Hallmarks analysis showed that PI3K pathway is activated in the tumors. Whole genome sequencing demonstrated a quiet genome with very few somatic alterations. The comparison of methylation profiles to publicly available data from ~ 3000 other central nervous system tumors revealed that PitBs have a distinct methylation profile compared to all other tumors, including pituitary adenomas. In conclusion, this comprehensive characterization of DICER1-related PitB revealed key molecular underpinnings of PitB and identified pathways that could potentially be exploited in the treatment of this tumor.

Invasive growth associated with cold-inducible RNA-binding protein expression drives recurrence of surgically resected brain metastases.

BACKGROUND: Sixty percent of surgically resected brain metastases (BrM) recur within 1 year. These recurrences have long been thought to result from the dispersion of cancer cells during surgery. We tested the alternative hypothesis that invasion of cancer cells into the adjacent brain plays a significant role in local recurrence and shortened overall survival. METHODS: We determined the invasion pattern of 164 surgically resected BrM and correlated with local recurrence and overall survival. We performed single-cell RNA sequencing (scRNAseq) of >15,000 cells from BrM and adjacent brain tissue. Validation of targets was performed with a novel cohort of BrM patient-derived xenografts (PDX) and patient tissues. RESULTS: We demonstrate that invasion of metastatic cancer cells into the adjacent brain is associated with local recurrence and shortened overall survival. scRNAseq of paired tumor and adjacent brain samples confirmed the existence of invasive cancer cells in the tumor-adjacent brain. Analysis of these cells identified cold-inducible RNA-binding protein (CIRBP) overexpression in invasive cancer cells compared to cancer cells located within the metastases. Applying PDX models that recapitulate the invasion pattern observed in patients, we show that CIRBP is overexpressed in highly invasive BrM and is required for efficient invasive growth in the brain. CONCLUSIONS: These data demonstrate peritumoral invasion as a driver of treatment failure in BrM that is functionally mediated by CIRBP. These findings improve our understanding of the biology underlying postoperative treatment failure and lay the groundwork for rational clinical trial development based upon invasion pattern in surgically resected BrM.

Author Correction: Single-cell RNA-seq reveals that glioblastoma recapitulates a normal neurodevelopmental hierarchy.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Single-cell RNA-seq reveals that glioblastoma recapitulates a normal neurodevelopmental hierarchy.

Cancer stem cells are critical for cancer initiation, development, and treatment resistance. Our understanding of these processes, and how they relate to glioblastoma heterogeneity, is limited. To overcome these limitations, we performed single-cell RNA sequencing on 53586 adult glioblastoma cells and 22637 normal human fetal brain cells, and compared the lineage hierarchy of the developing human brain to the transcriptome of cancer cells. We find a conserved neural tri-lineage cancer hierarchy centered around glial progenitor-like cells. We also find that this progenitor population contains the majority of the cancer's cycling cells, and, using RNA velocity, is often the originator of the other cell types. Finally, we show that this hierarchal map can be used to identify therapeutic targets specific to progenitor cancer stem cells. Our analyses show that normal brain development reconciles glioblastoma development, suggests a possible origin for glioblastoma hierarchy, and helps to identify cancer stem cell-specific targets.

Exome sequencing study of partial agenesis of the corpus callosum in men with developmental delay, epilepsy, and microcephaly.

BACKGROUND: This study reports the genetic features of four Caucasian males from the Saguenay-Lac-St-Jean region affected by partial agenesis of the corpus callosum (ACC) with hypotonia, epilepsy, developmental delay, microcephaly, hypoplasia, and autistic behavior. METHODS: We performed whole exome sequencing (WES) to identify new genes involved in this pathological phenotype. The regions of interest were subsequently sequenced for family members. RESULTS: Single-nucleotide variations (SNVs) and insertions or deletions were detected in genes potentially implicated in brain defects observed in these patients. One patient did not have mutations in genes related to ACC, but carried a de novo pathogenic mutation in Mucolipin-1 (MCOLN1) and was diagnosed with mucolipidosis type IV. Among the other probands, missense SNVs were observed in DCLK2 (Doublecortin Like Kinase 2), HERC2 (HECT And RLD Domain Containing E3 Ubiquitin Protein Ligase 2), and KCNH3 (Potassium channel, voltage-gated, subfamily H, member 3). One patient also carried a non-frameshift insertion in CACNA1A (Cav2.1(P/Q-type) calcium channels). CONCLUSION: Although no common genetic defect was observed in this study, we provide evidence for new avenues of investigation for ACC, such as molecular pathways involving HERC2, CACNA1A, KCNH3, and more importantly DCLK2. We also allowed to diagnose an individual with mucolipidosis type IV.

DGCR8 microprocessor defect characterizes familial multinodular goiter with schwannomatosis.

BACKGROUNDDICER1 is the only miRNA biogenesis component associated with an inherited tumor syndrome, featuring multinodular goiter (MNG) and rare pediatric-onset lesions. Other susceptibility genes for familial forms of MNG likely exist.METHODSWhole-exome sequencing of a kindred with early-onset MNG and schwannomatosis was followed by investigation of germline pathogenic variants that fully segregated with the disease. Genome-wide analyses were performed on 13 tissue samples from familial and nonfamilial DGCR8-E518K-positive tumors, including MNG, schwannomas, papillary thyroid cancers (PTCs), and Wilms tumors. miRNA profiles of 4 tissue types were compared, and sequencing of miRNA, pre-miRNA, and mRNA was performed in a subset of 9 schwannomas, 4 of which harbor DGCR8-E518K.RESULTSWe identified c.1552G>A;p.E518K in DGCR8, a microprocessor component located in 22q, in the kindred. The variant identified is a somatic hotspot in Wilms tumors and has been identified in 2 PTCs. Copy number loss of chromosome 22q, leading to loss of heterozygosity at the DGCR8 locus, was found in all 13 samples harboring c.1552G>A;p.E518K. miRNA profiling of PTCs, MNG, schwannomas, and Wilms tumors revealed a common profile among E518K hemizygous tumors. In vitro cleavage demonstrated improper processing of pre-miRNA by DGCR8-E518K. MicroRNA and RNA profiling show that this variant disrupts precursor microRNA production, impacting populations of canonical microRNAs and mirtrons.CONCLUSIONWe identified DGCR8 as the cause of an unreported autosomal dominant mendelian tumor susceptibility syndrome: familial multinodular goiter with schwannomatosis.FUNDINGCanadian Institutes of Health Research, Compute Canada, Alex's Lemonade Stand Foundation, the Mia Neri Foundation for Childhood Cancer, Cassa di Sovvenzioni e Risparmio fra il Personale della Banca d'Italia, and the KinderKrebsInitiative Buchholz/Holm-Seppensen.

Duplication 2p16 is associated with perisylvian polymicrogyria.

Polymicrogyria (PMG) is a heterogeneous brain malformation that may result from prenatal vascular disruption or infection, or from numerous genetic causes that still remain difficult to identify. We identified three unrelated patients with polymicrogyria and duplications of chromosome 2p, defined the smallest region of overlap, and performed gene pathway analysis using Cytoscape. The smallest region of overlap in all three children involved 2p16.1-p16.3. All three children have bilateral perisylvian polymicrogyria (BPP), intrauterine and postnatal growth deficiency, similar dysmorphic features, and poor feeding. Two of the three children had documented intellectual disability. Gene pathway analysis suggested a number of developmentally relevant genes and gene clusters that were over-represented in the critical region. We narrowed a rare locus for polymicrogyria to a region of 2p16.1-p16.3 that contains 33-34 genes, 23 of which are expressed in cerebral cortex during human fetal development. Using pathway analysis, we showed that several of the duplicated genes contribute to neurodevelopmental pathways including morphogen, cytokine, hormonal and growth factor signaling, regulation of cell cycle progression, cell morphogenesis, axonal guidance, and neuronal migration. These findings strengthen the evidence for a novel locus associated with polymicrogyria on 2p16.1-p16.3, and comprise the first step in defining the underlying genetic etiology.

Ovarian small cell carcinoma in one of a pair of monozygous twins.

One of a pair of monozygous twins was diagnosed and died of small cell carcinoma of the ovary of hypercalcemic type (SCCOHT) at the age of 30 years. Her sister remained unaffected and was very concerned about her risk for developing SCCOHT. By performing comprehensive molecular analysis using whole exome sequencing (WES) approach, we showed that the deceased twin's tumour has bi-allelic somatic genetic defects (a pathogenic frameshift deletion in SMARCA4 and LOH on chr19p). Results of WES of constitutional DNA from her unaffected sister were confirmatory. Based on our findings, we concluded that the living twin is not at risk for SCCOHT and does not need to consider preventive oophorectomy.

Publisher Correction: A PRDX1 mutant allele causes a MMACHC secondary epimutation in cblC patients.

The original version of this Article contained an error in the title, which was incorrectly given as 'APRDX1 mutant allele causes a MMACHC secondary epimutation in cblC patients'. This has now been corrected in both the PDF and HTML versions of the Article to read 'A PRDX1 mutant allele causes a MMACHC secondary epimutation in cblC patients'.

NHLRC2 variants identified in patients with fibrosis, neurodegeneration, and cerebral angiomatosis (FINCA): characterisation of a novel cerebropulmonary disease.

A novel multi-organ disease that is fatal in early childhood was identified in three patients from two non-consanguineous families. These children were born asymptomatic but at the age of 2 months they manifested progressive multi-organ symptoms resembling no previously known disease. The main clinical features included progressive cerebropulmonary symptoms, malabsorption, progressive growth failure, recurrent infections, chronic haemolytic anaemia and transient liver dysfunction. In the affected children, neuropathology revealed increased angiomatosis-like leptomeningeal, cortical and superficial white matter vascularisation and congestion, vacuolar degeneration and myelin loss in white matter, as well as neuronal degeneration. Interstitial fibrosis and previously undescribed granuloma-like lesions were observed in the lungs. Hepatomegaly, steatosis and collagen accumulation were detected in the liver. A whole-exome sequencing of the two unrelated families with the affected children revealed the transmission of two heterozygous variants in the NHL repeat-containing protein 2 (NHLRC2); an amino acid substitution p.Asp148Tyr and a frameshift 2-bp deletion p.Arg201GlyfsTer6. NHLRC2 is highly conserved and expressed in multiple organs and its function is unknown. It contains a thioredoxin-like domain; however, an insulin turbidity assay on human recombinant NHLRC2 showed no thioredoxin activity. In patient-derived fibroblasts, NHLRC2 levels were low, and only p.Asp148Tyr was expressed. Therefore, the allele with the frameshift deletion is likely non-functional. Development of the Nhlrc2 null mouse strain stalled before the morula stage. Morpholino knockdown of nhlrc2 in zebrafish embryos affected the integrity of cells in the midbrain region. This is the first description of a fatal, early-onset disease; we have named it FINCA disease based on the combination of pathological features that include fibrosis, neurodegeneration, and cerebral angiomatosis.

APRDX1 mutant allele causes a MMACHC secondary epimutation in cblC patients.

To date, epimutations reported in man have been somatic and erased in germlines. Here, we identify a cause of the autosomal recessive cblC class of inborn errors of vitamin B12 metabolism that we name "epi-cblC". The subjects are compound heterozygotes for a genetic mutation and for a promoter epimutation, detected in blood, fibroblasts, and sperm, at the MMACHC locus; 5-azacytidine restores the expression of MMACHC in fibroblasts. MMACHC is flanked by CCDC163P and PRDX1, which are in the opposite orientation. The epimutation is present in three generations and results from PRDX1 mutations that force antisense transcription of MMACHC thereby possibly generating a H3K36me3 mark. The silencing of PRDX1 transcription leads to partial hypomethylation of the epiallele and restores the expression of MMACHC. This example of epi-cblC demonstrates the need to search for compound epigenetic-genetic heterozygosity in patients with typical disease manifestation and genetic heterozygosity in disease-causing genes located in other gene trios.

Functionally Null RAD51D Missense Mutation Associates Strongly with Ovarian Carcinoma.

RAD51D is a key player in DNA repair by homologous recombination (HR), and RAD51D truncating variant carriers have an increased risk for ovarian cancer. However, the contribution of nontruncating RAD51D variants to cancer predisposition remains uncertain. Using deep sequencing and case-control genotyping studies, we show that in French Canadians, the missense RAD51D variant c.620C>T;p.S207L is highly prevalent and is associated with a significantly increased risk for ovarian high-grade serous carcinoma (HGSC; 3.8% cases vs. 0.2% controls). The frequency of the p.S207L variant did not significantly differ from that of controls in breast, endometrial, pancreas, or colorectal adenocarcinomas. Functionally, we show that this mutation impairs HR by disrupting the RAD51D-XRCC2 interaction and confers PARP inhibitor sensitivity. These results highlight the importance of a functional RAD51D-XRCC2 interaction to promote HR and prevent the development of HGSC. This study identifies c.620C>T;p.S207L as the first bona fide pathogenic RAD51D missense cancer susceptibility allele and supports the use of targeted PARP-inhibitor therapies in ovarian cancer patients carrying deleterious missense RAD51D variants. Cancer Res; 77(16); 4517-29. ©2017 AACR.