De novo missense variants in exon 9 of SEPHS1 cause a neurodevelopmental condition with developmental delay, poor growth, hypotonia, and dysmorphic features.

Selenophosphate synthetase (SEPHS) plays an essential role in selenium metabolism. Two mammalian SEPHS paralogues, SEPHS1 and SEPHS2, share high sequence identity and structural homology with SEPHS. Here, we report nine individuals from eight families with developmental delay, growth and feeding problems, hypotonia, and dysmorphic features, all with heterozygous missense variants in SEPHS1. Eight of these individuals had a recurrent variant at amino acid position 371 of SEPHS1 (p.Arg371Trp, p.Arg371Gln, and p.Arg371Gly); seven of these variants were known to be de novo. Structural modeling and biochemical assays were used to understand the effect of these variants on SEPHS1 function. We found that a variant at residue Trp352 results in local structural changes of the C-terminal region of SEPHS1 that decrease the overall thermal stability of the enzyme. In contrast, variants of a solvent-exposed residue Arg371 do not impact enzyme stability and folding but could modulate direct protein-protein interactions of SEPSH1 with cellular factors in promoting cell proliferation and development. In neuronal SH-SY5Y cells, we assessed the impact of SEPHS1 variants on cell proliferation and ROS production and investigated the mRNA expression levels of genes encoding stress-related selenoproteins. Our findings provided evidence that the identified SEPHS1 variants enhance cell proliferation by modulating ROS homeostasis. Our study supports the hypothesis that SEPHS1 plays a critical role during human development and provides a basis for further investigation into the molecular mechanisms employed by SEPHS1. Furthermore, our data suggest that variants in SEPHS1 are associated with a neurodevelopmental disorder.

Evaluation of the diagnostic accuracy of exome sequencing and its impact on diagnostic thinking for patients with rare disease in a publicly funded health care system: A prospective cohort study.

PURPOSE: To evaluate the diagnostic utility of publicly funded clinical exome sequencing (ES) for patients with suspected rare genetic diseases. METHODS: We prospectively enrolled 297 probands who met eligibility criteria and received ES across 5 sites in Ontario, Canada, and extracted data from medical records and clinician surveys. Using the Fryback and Thornbury Efficacy Framework, we assessed diagnostic accuracy by examining laboratory interpretation of results and assessed diagnostic thinking by examining the clinical interpretation of results and whether clinical-molecular diagnoses would have been achieved via alternative hypothetical molecular tests. RESULTS: Laboratories reported 105 molecular diagnoses and 165 uncertain results in known and novel genes. Of these, clinicians interpreted 102 of 105 (97%) molecular diagnoses and 6 of 165 (4%) uncertain results as clinical-molecular diagnoses. The 108 clinical-molecular diagnoses were in 104 families (35% diagnostic yield). Each eligibility criteria resulted in diagnostic yields of 30% to 40%, and higher yields were achieved when >2 eligibility criteria were met (up to 45%). Hypothetical tests would have identified 61% of clinical-molecular diagnoses. CONCLUSION: We demonstrate robustness in eligibility criteria and high clinical validity of laboratory results from ES testing. The importance of ES was highlighted by the potential 40% of patients that would have gone undiagnosed without this test.

Landscape heterogeneity provides co-benefits to predator and prey.

Predator populations are imperiled globally, due in part to changing habitat and trophic interactions. Theoretical and laboratory studies suggest that heterogeneous landscapes containing prey refuges acting as source habitats can benefit both predator and prey populations, although the importance of heterogeneity in natural systems is uncertain. Here, we tested the hypothesis that landscape heterogeneity mediates predator-prey interactions between the California spotted owl (Strix occidentalis occidentalis)-a mature forest species-and one of its principal prey, the dusky-footed woodrat (Neotoma fuscipes)-a younger forest species-to the benefit of both. We did so by combining estimates of woodrat density and survival from live trapping and very high frequency tracking with direct observations of prey deliveries to dependent young by owls in both heterogeneous and homogeneous home ranges. Woodrat abundance was ~2.5 times higher in owl home ranges (14.12 km2 ) featuring greater heterogeneity in vegetation types (1805.0 ± 50.2 SE) compared to those dominated by mature forest (727.3 ± 51.9 SE), in large part because of high densities in young forests appearing to act as sources promoting woodrat densities in nearby mature forests. Woodrat mortality rates were low across vegetation types and did not differ between heterogeneous and homogeneous home ranges, yet all observed predation by owls occurred within mature forests, suggesting young forests may act as woodrat refuges. Owls exhibited a type 1 functional response, consuming ~2.5 times more woodrats in heterogeneous (31.1/month ± 5.2 SE) versus homogeneous (12.7/month ± 3.7 SE) home ranges. While consumption of smaller-bodied alternative prey partially compensated for lower woodrat consumption in homogeneous home ranges, owls nevertheless consumed 30% more biomass in heterogeneous home ranges-approximately equivalent to the energetic needs of producing one additional offspring. Thus, a mosaic of vegetation types including young forest patches increased woodrat abundance and availability that, in turn, provided energetic and potentially reproductive benefits to mature forest-associated spotted owls. More broadly, our findings provide strong empirical evidence that heterogeneous landscapes containing prey refuges can benefit both predator and prey populations. As anthropogenic activities continue to homogenize landscapes globally, promoting heterogeneous systems with prey refuges may benefit imperiled predators.

Deep phenotyping of the neuroimaging and skeletal features in KBG syndrome: a study of 53 patients and review of the literature.

BACKGROUND: KBG syndrome is caused by haploinsufficiency of ANKRD11 and is characterised by macrodontia of upper central incisors, distinctive facial features, short stature, skeletal anomalies, developmental delay, brain malformations and seizures. The central nervous system (CNS) and skeletal features remain poorly defined. METHODS: CNS and/or skeletal imaging were collected from molecularly confirmed individuals with KBG syndrome through an international network. We evaluated the original imaging and compared our results with data in the literature. RESULTS: We identified 53 individuals, 44 with CNS and 40 with skeletal imaging. Common CNS findings included incomplete hippocampal inversion and posterior fossa malformations; these were significantly more common than previously reported (63.4% and 65.9% vs 1.1% and 24.7%, respectively). Additional features included patulous internal auditory canal, never described before in KBG syndrome, and the recurrence of ventriculomegaly, encephalic cysts, empty sella and low-lying conus medullaris. We found no correlation between these structural anomalies and epilepsy or intellectual disability. Prevalent skeletal findings comprised abnormalities of the spine including scoliosis, coccygeal anomalies and cervical ribs. Hand X-rays revealed frequent abnormalities of carpal bone morphology and maturation, including a greater delay in ossification compared with metacarpal/phalanx bones. CONCLUSION: This cohort enabled us to describe the prevalence of very heterogeneous neuroradiological and skeletal anomalies in KBG syndrome. Knowledge of the spectrum of such anomalies will aid diagnostic accuracy, improve patient care and provide a reference for future research on the effects of ANKRD11 variants in skeletal and brain development.

Macrocephaly and developmental delay caused by missense variants in RAB5C.

Rab GTPases are important regulators of intracellular vesicular trafficking. RAB5C is a member of the Rab GTPase family that plays an important role in the endocytic pathway, membrane protein recycling and signaling. Here we report on 12 individuals with nine different heterozygous de novo variants in RAB5C. All but one patient with missense variants (n = 9) exhibited macrocephaly, combined with mild-to-moderate developmental delay. Patients with loss of function variants (n = 2) had an apparently more severe clinical phenotype with refractory epilepsy and intellectual disability but a normal head circumference. Four missense variants were investigated experimentally. In vitro biochemical studies revealed that all four variants were damaging, resulting in increased nucleotide exchange rate, attenuated responsivity to guanine exchange factors and heterogeneous effects on interactions with effector proteins. Studies in C. elegans confirmed that all four variants were damaging in vivo and showed defects in endocytic pathway function. The variant heterozygotes displayed phenotypes that were not observed in null heterozygotes, with two shown to be through a dominant negative mechanism. Expression of the human RAB5C variants in zebrafish embryos resulted in defective development, further underscoring the damaging effects of the RAB5C variants. Our combined bioinformatic, in vitro and in vivo experimental studies and clinical data support the association of RAB5C missense variants with a neurodevelopmental disorder characterized by macrocephaly and mild-to-moderate developmental delay through disruption of the endocytic pathway.

Protocol for a Prospective, Observational Cost-effectiveness Analysis of Returning Secondary Findings of Genome Sequencing for Unexplained Suspected Genetic Conditions.

PURPOSE: Although costly, genome-wide sequencing (GWS) detects an extensive range of variants, enhancing our ability to diagnose and assess risk for an increasing number of diseases. In addition to detecting variants related to the indication for testing, GWS can detect secondary variants in BRCA1, BRCA2, and other genes for which early intervention may improve health. As the list of secondary findings grows, there is increased demand for surveillance and management by multiple specialists, adding pressure to constrained health care budgets. Secondary finding testing is actively debated because some consider it opportunistic screening for future health risks that may not manifest. Given the economic implications of secondary finding testing and follow-up and its unproven clinical utility, the objective is to assess the incremental cost-effectiveness of secondary finding ascertainment per case detected and per unit of improved clinical utility in families of children with unexplained suspected genetic conditions undergoing clinical GWS. METHODS: Those undergoing trio genome or exome sequencing are eligible for the study. Positive secondary finding index cases will be matched to negative controls (1:2) based on age group, primary result(s) type, and clinical indication. During the 2-year study, 71 cases and 142 matched controls are expected. Health service use will be collected in patients and 1 adult family member every 6 months. The per-child and per-dyad total cost will be determined by multiplying use of each resource by a corresponding unit price and summing all cost items. Costs will be estimated from the public and societal payer perspectives. The mean cost per child and per dyad for secondary finding-positive and secondary finding-negative groups will be compared statistically. If important demographic differences are observed between groups, ordinary least-squares regression, log transformation, or other nonparametric technique will be used to compare adjusted mean costs. The ratio of the difference in mean cost to the secondary finding yield will be used to estimate incremental cost-effectiveness. In secondary analyses, effectiveness will be estimated using the number of clinical management changes due to secondary findings or the Clinician-Reported Genetic Testing Utility Index (C-GUIDE) score, a validated measure of clinical utility. Sensitivity analysis will be undertaken to assess the robustness of the findings to variation in key parameters. IMPLICATIONS: This study generates key evidence to inform clinical practice and funding allocation related to secondary finding testing. The inclusion of family members and a new measure of clinical utility represent important advancements in economic evaluation in genomics.

Canakinumab in addition to phosphate-binding and phosphaturia-inducing therapy were effective in achieving remission in a child with a large familial calcinotic tumour.

We describe the clinical evolution of a patient with tumoral calcinosis due to a pathogenic variant in the GALNT3 gene presented with a large mass overlying her left hip associated complicated by inflammatory flares. Therapy (sevelamer, acetazolamide, and probenecid) was unsuccessful in preventing tumour surgeries, therefore, interleukin-1ß monoclonal antibody therapy was added; this was successful in the prevention of tumour re-growth. This case highlights the importance of assessing and treating the inflammatory aspect of calcinotic tumour.

Multi-trophic occupancy modeling connects temporal dynamics of woodpeckers and beetle sign following fire.

In conifer forests of western North America, wildlife populations can change rapidly in the decade following wildfire as trees die and animals respond to concomitant resource pulses that occur across multiple trophic levels. In particular, black-backed woodpeckers (Picoides arcticus) show predictable temporal increases then declines following fire; this trajectory is widely believed to be a response to the woodpeckers' main prey, woodboring beetle larvae of the families Buprestidae and Cerambycidae, but we lack understanding of how abundances of these predators and prey may be associated in time or space. Here, we pair woodpecker surveys over 10 years with surveys of woodboring beetle sign and activity, collected at 128 survey plots across 22 recent fires, to ask whether accumulated beetle sign indicates current or past black-backed woodpecker occurrence, and whether that relationship is mediated by the number of years since fire. We test this relationship using an integrative multi-trophic occupancy model. Our results demonstrate that woodboring beetle sign is a positive indicator of woodpecker presence 1-3 years following fire, an uninformative indicator from 4-6 years after fire, and a negative indicator beginning 7 years following fire. Woodboring beetle activity, itself, is temporally variable and dependent on tree species composition, with beetle sign generally accumulating over time, particularly in stands with diverse tree communities, but decreasing over time in Pinus-dominated stands where faster bark decay rates lead to brief pulses of beetle activity followed by rapid degradation of tree substrate and accumulated beetle sign. Altogether, the strong connections of woodpecker occurrence to beetle activity support prior hypotheses of how multi-trophic interactions govern rapid temporal dynamics of primary and secondary consumers in burned forests. While our results indicate that beetle sign is, at best, a rapidly shifting and potentially misleading measure of woodpecker occurrence, the better we understand the interacting mechanisms underlying temporally dynamic systems, the more successfully we will be able to predict the outcomes of management actions.

Incorporating pyrodiversity into wildlife habitat assessments for rapid post-fire management: A woodpecker case study.

Spatial and temporal variation in fire characteristics-termed pyrodiversity-are increasingly recognized as important factors that structure wildlife communities in fire-prone ecosystems, yet there have been few attempts to incorporate pyrodiversity or post-fire habitat dynamics into predictive models of animal distributions and abundance to support post-fire management. We use the black-backed woodpecker-a species associated with burned forests-as a case study to demonstrate a pathway for incorporating pyrodiversity into wildlife habitat assessments for adaptive management. Employing monitoring data (2009-2019) from post-fire forests in California, we developed three competing occupancy models describing different hypotheses for habitat associations: (1) a static model representing an existing management tool, (2) a temporal model accounting for years since fire, and (3) a temporal-landscape model which additionally incorporates emerging evidence from field studies about the influence of pyrodiversity. Evaluating predictive ability, we found superior support for the temporal-landscape model, which showed a positive relationship between occupancy and pyrodiversity and interactions between habitat associations and years since fire. We incorporated the new temporal-landscape model into an RShiny application to make this decision-support tool accessible to decision-makers.

Vaginal Metastases of Wilms' Tumor in a Pediatric Patient: A Rare Case.

BACKGROUND: Wilms' tumor is the second most common pediatric abdominal cancer; however, it rarely involves the female reproductive tract. There are few cases reported in the literature describing uterine, ovarian, cervical, and vaginal involvement. CASE: We report the case of a 7-year-old girl presenting with a large renal mass with retroperitoneal nodal and lung metastases; she was diagnosed with stage 4 favorable histology Wilms' tumor. She was treated with surgery, chemotherapy, and radiation. She presented with vaginal bleeding 10 months after completing treatment; biopsy of a vaginal mass confirmed recurrence, and this was sent for molecular profiling, which did not identify an inherited cancer predisposition or targetable mutation. She was again treated with chemotherapy; examination redemonstrated a small vaginal mass, but re-biopsy of the lesion was negative for malignancy. Due to high risk of local relapse, ongoing chemotherapy and pelvic radiation ensued. End-of-treatment imaging and vaginoscopy showed no residual disease. SUMMARY AND CONCLUSION: Vaginal metastases of Wilms' tumor are very rare; this is the second reported case in the literature. Pediatric clinicians should have a strong suspicion for vaginal metastases in cancer patients presenting with vaginal bleeding, especially when their pubertal development does not suggest that bleeding would be secondary to menarche. Long-term gynecologic care for these patients is paramount to reduce morbidity from chemotherapy and pelvic radiation. Fertility preservation counselling should be made early, through referral to a specialist.

Rhabdomyosarcoma as the first presentation in Neurofibromatosis Type 1: case series and review of the literature.

Neurofibromatosis Type 1 (NF1) is a neurocutaneous syndrome characterized by multiple café-au-lait macules, neurofibromas, and predisposition to malignancies, including rhabdomyosarcomas (RMS). Somatic NF1 mutations occur in RMS and other cancers, and ∼1% of patients with RMS have NF1. We describe three patients who presented prior to one year of age with RMS and were subsequently diagnosed with NF1. Compared to sporadic RMS, patients with this cancer predisposition syndrome are diagnosed younger, genitourinary sites are more common, and tumors are almost exclusively the embryonal subtype. Genomic sequencing of the tumor was initiated in one patient, and we identified a second sequence variant in NF1. The identification of molecular drivers in tumors is changing the nature of pediatric oncology by informing therapeutics targeted to specific molecular pathways and selecting patients who are likely to harbor germline variants in cancer predisposition genes who would benefit from a Medical Genetics assessment.

The implementation of an enhanced clinical model to improve the diagnostic yield of exome sequencing for patients with a rare genetic disease: A Canadian experience.

The introduction of clinical exome sequencing (ES) has provided a unique opportunity to decrease the diagnostic odyssey for patients living with a rare genetic disease (RGD). ES has been shown to provide a diagnosis in 29%-57% of patients with a suspected RGD, with as many as 70% remaining undiagnosed. There is a need to advance the clinical model of care by more formally integrating approaches that were previously considered research into an enhanced diagnostic workflow. We developed an Exome Clinic, which set out to evaluate a workflow for improving the diagnostic yield of ES for patients with an undiagnosed RGD. Here, we report the outcomes of 47 families who underwent clinical ES in the first year of the clinic. The diagnostic yield from clinical ES was 40% (19/47). Families who remained undiagnosed after ES had the opportunity for follow-up studies that included phenotyping and candidate variant segregation in relatives, genomic matchmaking, and ES reanalysis. This enhanced diagnostic workflow increased the diagnostic yield to 55% (26/47), predominantly through the resolution of variants and genes of uncertain significance. We advocate that this approach be integrated into mainstream clinical practice and highlight the importance of a coordinated translational approach for patients with RGD.

Bridging clinical care and research in Ontario, Canada: Maximizing diagnoses from reanalysis of clinical exome sequencing data.

We examined the utility of clinical and research processes in the reanalysis of publicly-funded clinical exome sequencing data in Ontario, Canada. In partnership with eight sites, we recruited 287 families with suspected rare genetic diseases tested between 2014 and 2020. Data from seven laboratories was reanalyzed with the referring clinicians. Reanalysis of clinically relevant genes identified diagnoses in 4% (13/287); four were missed by clinical testing. Translational research methods, including analysis of novel candidate genes, identified candidates in 21% (61/287). Of these, 24 families have additional evidence through data sharing to support likely diagnoses (8% of cohort). This study indicates few diagnoses are missed by clinical laboratories, the incremental gain from reanalysis of clinically-relevant genes is modest, and the highest yield comes from validation of novel disease-gene associations. Future implementation of translational research methods, including continued reporting of compelling genes of uncertain significance by clinical laboratories, should be considered to maximize diagnoses.

Mega-disturbances cause rapid decline of mature conifer forest habitat in California.

Mature forests provide important wildlife habitat and support critical ecosystem functions globally. Within the dry conifer forests of the western United States, past management and fire exclusion have contributed to forest conditions that are susceptible to increasingly severe wildfire and drought. We evaluated declines in conifer forest cover in the southern Sierra Nevada of California during a decade of record disturbance by using spatially comprehensive forest structure estimates, wildfire perimeter data, and the eDaRT forest disturbance tracking algorithm. Primarily due to the combination of wildfires, drought, and drought-associated beetle epidemics, 30% of the region's conifer forest extent transitioned to nonforest vegetation during 2011-2020. In total, 50% of mature forest habitat and 85% of high density mature forests either transitioned to lower density forest or nonforest vegetation types. California spotted owl protected activity centers (PAC) experienced greater canopy cover decline (49% of 2011 cover) than non-PAC areas (42% decline). Areas with high initial canopy cover and without tall trees were most vulnerable to canopy cover declines, likely explaining the disproportionate declines of mature forest habitat and within PACs. Drought and beetle attack caused greater cumulative declines than areas where drought and wildfire mortality overlapped, and both types of natural disturbance far outpaced declines attributable to mechanical activities. Drought mortality that disproportionately affects large conifers is particularly problematic to mature forest specialist species reliant on large trees. However, patches of degraded forests within wildfire perimeters were larger with greater core area than those outside burned areas, and remnant forest habitats were more fragmented within burned perimeters than those affected by drought and beetle mortality alone. The percentage of mature forest that survived and potentially benefited from lower severity wildfire increased over time as the total extent of mature forest declined. These areas provide some opportunity for improved resilience to future disturbances, but strategic management interventions are likely also necessary to mitigate worsening mega-disturbances. Remaining dry mature forest habitat in California may be susceptible to complete loss in the coming decades without a rapid transition from a conservation paradigm that attempts to maintain static conditions to one that manages for sustainable disturbance dynamics.

Circulating tumor DNA monitoring for early recurrence detection in epithelial ovarian cancer.

OBJECTIVE: Epithelial ovarian cancer (EOC) is the most lethal gynecologic malignancy. We examined the utility of circulating tumor DNA (ctDNA) as a prognostic biomarker for EOC by assessing its relationship with patient outcome and CA-125, pre-surgically and during post-treatment surveillance. METHODS: Plasma samples were collected from patients with stage I-IV EOC. Cohort A included patients with pre-surgical samples (N = 44, median follow-up: 2.7 years), cohort B and C included: patients with serially collected post-surgically (N = 12) and, during surveillance (N = 13), respectively (median follow-up: 2 years). Plasma samples were analyzed using a tumor-informed, personalized multiplex-PCR NGS assay; ctDNA status and CA-125 levels were correlated with clinical features and outcomes. RESULTS: Genomic profiling was performed on the entire cohort and was consistent with that seen in TCGA. In cohort A, ctDNA-positivity was observed in 73% (32/44) of presurgical samples and was higher in high nuclear grade disease. In cohort B and C, ctDNA was only detected in patients who relapsed (100% sensitivity and specificity) and preceded radiological findings by an average of 10 months. The presence of ctDNA at a single timepoint after completion of surgery +/- adjuvant chemotherapy and serially during surveillance was a strong predictor of relapse (HR:17.6, p = 0.001 and p < 0.0001, respectively), while CA-125 positivity was not (p = 0.113 and p = 0.056). CONCLUSIONS: The presence of ctDNA post-surgically is highly prognostic of reduced recurrence-free survival. CtDNA outperformed CA-125 in identifying patients at highest risk of recurrence. These results suggest that monitoring ctDNA could be beneficial in clinical decision-making for EOC patients.

Comparing genome sequencing technologies to improve rare disease diagnostics: a protocol for the evaluation of a pilot project, Genome-wide Sequencing Ontario.

BACKGROUND: Genome-wide sequencing has emerged as a promising strategy for the timely diagnosis of rare diseases, but it is not yet available as a clinical test performed in Canadian diagnostic laboratories. We describe the protocol for evaluating a 2-year pilot project, Genome-wide Sequencing Ontario, to offer high-quality clinical genome-wide sequencing in Ontario, Canada. METHODS: The Genome-wide Sequencing Ontario protocol was codesigned by the Ontario Ministry of Health, the Hospital for Sick Children in Toronto and the Children's Hospital of Eastern Ontario in Ottawa. Enrolment of a prospective cohort of patients began on Apr. 1, 2021. Eligible cases with blood samples available for the index case and both parents (i.e., trios) are randomized to receive exome sequencing or genome sequencing. We will collect patient-level data and ascertain costs associated with the laboratory workflow for exome sequencing and genome sequencing. We will compare point estimates for the diagnostic utility and timeliness of exome sequencing and genome sequencing, and we will determine an incremental cost-effectiveness ratio (expressed as the incremental cost of genome sequencing versus exome sequencing per additional patient with a causal variant detected). INTERPRETATION: Findings from this work will provide robust evidence for the diagnostic utility, cost-effectiveness and timeliness of exome sequencing and genome sequencing, and will be disseminated via academic publications and policy briefs. Findings will inform provincial and cross-provincial policy related to the long-term organization, delivery and reimbursement of clinical-grade genome diagnostics for rare disease.

BESPOKE IO protocol: a multicentre, prospective observational study evaluating the utility of ctDNA in guiding immunotherapy in patients with advanced solid tumours.

INTRODUCTION: Immunotherapy (IO) has transformed the treatment paradigm for a wide variety of solid tumours. However, assessment of response can be challenging with conventional radiological imaging (eg, iRECIST), which do not precisely capture the unique response patterns of tumours treated with IO. Emerging data suggest that circulating tumour DNA (ctDNA) can aid in response assessment in patients with solid tumours receiving IO. The short half-life of ctDNA puts it in a unique position for early treatment response monitoring. The BESPOKE IO study is designed to investigate the clinical utility of serial ctDNA testing to assess treatment response using a tumour-informed, bespoke ctDNA assay (Signatera) and to determine its impact on clinical decision-making with respect to continuation/discontinuation, or escalation/de-escalation of immunotherapy in patients with advanced solid tumours. METHODS AND ANALYSIS: The BESPOKE IO is a multicentre, prospective, observational study with a goal to enroll over 1500 patients with solid tumours receiving IO in up to 100 US sites. Patients will be followed for up to 2 years with serial ctDNA analysis, timed with every other treatment cycle. The primary endpoint is to determine the percentage of patients who will have their treatment regimen changed as guided by post-treatment bespoke ctDNA results along with standard response assessment tools. The major secondary endpoints include progression-free survival, overall survival and overall response rate based on the ctDNA dynamics. ETHICS AND DISSEMINATION: The BESPOKE IO study was approved by the WCG Institutional Review Board (Natera-20-043-NCP BESPOKE Study of ctDNA Guided Immunotherapy (BESPOKE IO)) on 22 February 2021. Data protection and privacy regulations will be strictly observed in the capturing, forwarding, processing and storing patients' data. Natera will approve the publication of any study results in accordance with the site-specific contract. TRIAL REGISTRATION NUMBER: NCT04761783.