Reanalysis of genomic data in rare disease: current practice and attitudes among Australian clinical and laboratory genetics services.

Reanalyzing stored genomic data over time is highly effective in increasing diagnostic yield in rare disease. Automation holds the promise of delivering the benefits of reanalysis at scale. Our study aimed to understand current reanalysis practices among Australian clinical and laboratory genetics services and explore attitudes towards large-scale automated re-analysis. We collected audit data regarding testing and reanalysis volumes, policies and procedures from all Australian diagnostic laboratories providing rare disease genomic testing. A genetic health professionals' survey explored current practices, barriers to reanalysis, preferences and attitudes towards automation. Between 2018 and 2021, Australian diagnostic laboratories performed over 25,000 new genomic tests and 950 reanalyses, predominantly in response to clinician requests. Laboratory and clinical genetic health professionals (N = 134) identified workforce capacity as the principal barrier to reanalysis. No specific laboratory or clinical guidelines for genomic data reanalysis or policies were identified nationally. Perceptions of acceptability and feasibility of automating reanalysis were positive, with professionals emphasizing clinical and workflow benefits. In conclusion, there is a large and rapidly growing unmet need for reanalysis of existing genomic data. Beyond developing scalable automated reanalysis pipelines, leadership and policy are needed to successfully transform service delivery models and maximize clinical benefit.

Shariant platform: Enabling evidence sharing across Australian clinical genetic-testing laboratories to support variant interpretation.

Sharing genomic variant interpretations across laboratories promotes consistency in variant assertions. A landscape analysis of Australian clinical genetic-testing laboratories in 2017 identified that, despite the national-accreditation-body recommendations encouraging laboratories to submit genotypic data to clinical databases, fewer than 300 variants had been shared to the ClinVar public database. Consultations with Australian laboratories identified resource constraints limiting routine application of manual processes, consent issues, and differences in interpretation systems as barriers to sharing. This information was used to define key needs and solutions required to enable national sharing of variant interpretations. The Shariant platform, using both the GRCh37 and GRCh38 genome builds, was developed to enable ongoing sharing of variant interpretations and associated evidence between Australian clinical genetic-testing laboratories. Where possible, two-way automated sharing was implemented so that disruption to laboratory workflows would be minimized. Terms of use were developed through consultation and currently restrict access to Australian clinical genetic-testing laboratories. Shariant was designed to store and compare structured evidence, to promote and record resolution of inter-laboratory classification discrepancies, and to streamline the submission of variant assertions to ClinVar. As of December 2021, more than 14,000 largely prospectively curated variant records from 11 participating laboratories have been shared. Discrepant classifications have been identified for 11% (28/260) of variants submitted by more than one laboratory. We have demonstrated that co-design with clinical laboratories is vital to developing and implementing a national variant-interpretation sharing effort. This approach has improved inter-laboratory concordance and enabled opportunities to standardize interpretation practices.

Molecular epidemiology of enteroviruses associated with hand, foot, and mouth disease in South India from 2015 to 2017.

Hand, foot, and mouth disease (HFMD) is a common childhood infection caused by human enteroviruses and is clinically characterised by fever with vesicular rash on the hands, feet, and mouth. While enterovirus A71 (EV-A71) and coxsackievirus A16 (CVA16) were the major etiological agents of HFMD in India earlier, the data on recently circulating enteroviruses associated with HFMD are sparse. Here, we describe the molecular epidemiology of enteroviruses associated with HFMD in South India from 2015 to 2017. We used archived enterovirus real-time reverse transcription (RT) PCR-positive vesicle swab and/or throat swab specimens from clinically suspected HFMD cases collected from four secondary-care hospitals in South India between July 2015 and December 2017. PCR amplification and sequencing were done based on the 5'VP1, 3'VP1, VP2, or 5´NCR regions to identify enterovirus types. Genetic diversity among enteroviruses was inferred by phylogenetic analysis. Of the 107 enterovirus RNA real-time RT-PCR-positive HFMD cases, 69 (64%) were typed as CVA6, 16 (15%) were CVA16, and one (1%) was CVA10, whereas in 21 (20%) cases, the virus was not typeable by any of the methods used in the study. The majority of HFMD cases (89, 83%) were in children less than five years old, while 11 (10.3%) were in adults. 5'VP1 yielded the maximum number of enteroviruses genotyped, and phylogenetic analysis showed that the CVA6 strains belonged to subclade D3, while the subclades of CVA16 and CVA10 were B1c and D, respectively. The predominant etiological agent of HFMD in South India during 2015-2017 was CVA6, followed by CVA16 and CVA10.

Methotrexate-related central neurotoxicity: clinical characteristics, risk factors and genome-wide association study in children treated for acute lymphoblastic leukemia.

Symptomatic methotrexate-related central neurotoxicity (MTX neurotoxicity) is a severe toxicity experienced during acute lymphoblastic leukemia (ALL) therapy with potential long-term neurologic complications. Risk factors and long-term outcomes require further study. We conducted a systematic, retrospective review of 1,251 consecutive Australian children enrolled on Berlin-Frankfurt-Münster or Children's Oncology Group-based protocols between 1998-2013. Clinical risk predictors for MTX neurotoxicity were analyzed using regression. A genome-wide association study (GWAS) was performed on 48 cases and 537 controls. The incidence of MTX neurotoxicity was 7.6% (n=95 of 1,251), at a median of 4 months from ALL diagnosis and 8 days after intravenous or intrathecal MTX. Grade 3 elevation of serum aspartate aminotransferase (P=0.005, odds ratio 2.31 [range, 1.28-4.16]) in induction/consolidation was associated with MTX neurotoxicity, after accounting for the only established risk factor, age ≥10 years. Cumulative incidence of CNS relapse was increased in children where intrathecal MTX was omitted following symptomatic MTX neurotoxicity (n=48) compared to where intrathecal MTX was continued throughout therapy (n=1,174) (P=0.047). Five-year central nervous system relapse-free survival was 89.2 4.6% when intrathecal MTX was ceased compared to 95.4 0.6% when intrathecal MTX was continued. Recurrence of MTX neurotoxicity was low (12.9%) for patients whose intrathecal MTX was continued after their first episode. The GWAS identified single-nucletide polymorphism associated with MTX neurotoxicity near genes regulating neuronal growth, neuronal differentiation and cytoskeletal organization (P<1x10-6). In conclusion, increased serum aspartate aminotransferase and age ≥10 years at diagnosis were independent risk factors for MTX neurotoxicity. Our data do not support cessation of intrathecal MTX after a first MTX neurotoxicity event.

The relation between age of retirement and frailty in later life? A cross-sectional study in Flemish older adults.

BACKGROUND: Policymakers in several European countries, concerned about the sustainability of their pension system, have raised the statutory retirement age. While several studies investigated the effect of retirement on health, the relationship between retirement and frailty is neglected. Notwithstanding, frailty is associated with adverse outcomes. OBJECTIVE: The aim of this study was to examine the relationship between age of retirement and frailty in later life. METHODS: Data of the Belgian Ageing Studies, a cross-sectional research project was used. The present study includes N=12659 participants (>60y) in 83 Flemish municipalities. To address reverse causality, only participants not retired because of health-related reasons were included. The Comprehensive Frailty Assessment Instrument, a multidimensional frailty scale with four domains (physical, psychological, social and environmental) was used to operationalize frailty. Univariate general linear regression analyses (GLM) were performed for scores on the total frailty scale and the four subdomains separately. The analysis was done for men and women separately, since both groups have different labor trajectories. RESULTS: The present study found a negative association between age of retirement and physical frailty for both men and women in later life, and total frailty for men, although the differences were small. No evidence was found for a relation between age of retirement and the other subdomains of frailty. CONCLUSIONS: The results suggest that age of retirement is not a clinically relevant predictor for frailty in later life. Differences within and between subpopulations (e.g., profession) can shed a new light on this relation.

Identification of a Locus Near ULK1 Associated With Progression-Free Survival in Ovarian Cancer.

BACKGROUND: Many loci have been found to be associated with risk of epithelial ovarian cancer (EOC). However, although there is considerable variation in progression-free survival (PFS), no loci have been found to be associated with outcome at genome-wide levels of significance. METHODS: We carried out a genome-wide association study (GWAS) of PFS in 2,352 women with EOC who had undergone cytoreductive surgery and standard carboplatin/paclitaxel chemotherapy. RESULTS: We found seven SNPs at 12q24.33 associated with PFS (P < 5 × 10-8), the top SNP being rs10794418 (HR = 1.24; 95% CI, 1.15-1.34; P = 1.47 × 10-8). High expression of a nearby gene, ULK1, is associated with shorter PFS in EOC, and with poor prognosis in other cancers. SNP rs10794418 is also associated with expression of ULK1 in ovarian tumors, with the allele associated with shorter PFS being associated with higher expression, and chromatin interactions were detected between the ULK1 promoter and associated SNPs in serous and endometrioid EOC cell lines. ULK1 knockout ovarian cancer cell lines showed significantly increased sensitivity to carboplatin in vitro. CONCLUSIONS: The locus at 12q24.33 represents one of the first genome-wide significant loci for survival for any cancer. ULK1 is a plausible candidate for the target of this association. IMPACT: This finding provides insight into genetic markers associated with EOC outcome and potential treatment options.See related commentary by Peres and Monteiro, p. 1604.

Genetic analysis of Enterovirus D68 associated with pneumonia in children from South India.

EV-D68 is an emerging enterovirus infection associated with severe acute respiratory illness (SARI), acute flaccid myelitis (AFM) and acute flaccid paralysis (AFP). While EV-D68 outbreaks and sporadic cases are reported globally, a single case has been reported from India. The present study aims to investigate the molecular epidemiology and clinical characteristics of EV-D68-associated SARI cases from South India. We screened influenza-negative archived throat swab specimens from Influenza-Like Illness (ILI) and SARI cases (n=959; 2016 to 2018 period) for enteroviruses by pan-enterovirus real-time RT-PCR. Thirteen samples positive for enteroviruses were typed by PCR and sequencing based on VPI, VP2 and/or 5'NCR regions. One EV-D68 RNA sample was subjected to next-generation sequencing for whole genome characterisation. Among 13 enterovirus cases, four were ECHO-11, three EV-D68, two CV-A16 and one each EV-71, CV-B1, CV-B2 and CV-A9. All three cases of EV-D68 infection were reported in children below 2 years of age from Kerala state of South India during June and July 2017. The patients developed pneumonia without any neurological complications. Sequencing based on VPI and 5'NCR regions showed that EV-D68 strains belong to the novel subclade B3. The EV-D68 complete genome identified with two unique amino acid substitutions in VP1 (T-246-I) and 3D (K-344-R) regions. This study reiterates the EV-D68 novel subclade B3 circulation in India and indicates the urgent need for structured EV-D68 surveillance in the country to describe the epidemiology.

Investigation of current models of care for genetic heart disease in Australia: A national clinical audit.

BACKGROUND: This sub-study of the Australian Genomics Cardiovascular Genetic Disorders Flagship sought to conduct the first nation-wide audit in Australia to establish the current practices across cardiac genetics clinics. METHOD: An audit of records of patients with a suspected genetic heart disease (cardiomyopathy, primary arrhythmia, autosomal dominant congenital heart disease) who had a cardiac genetics consultation between 1st January 2016 and 31 July 2018 and were offered a diagnostic genetic test. RESULTS: This audit included 536 records at multidisciplinary cardiac genetics clinics from 11 public tertiary hospitals across five Australian states. Most genetic consultations occurred in a clinic setting (90%), followed by inpatient (6%) and Telehealth (4%). Queensland had the highest proportion of Telehealth consultations (9% of state total). Sixty-six percent of patients had a clinical diagnosis of a cardiomyopathy, 28% a primary arrhythmia, and 0.7% congenital heart disease. The reason for diagnosis was most commonly as a result of investigations of symptoms (73%). Most patients were referred by a cardiologist (85%), followed by a general practitioner (9%) and most genetic tests were funded by the state Genetic Health Service (73%). Nationally, 29% of genetic tests identified a pathogenic or likely pathogenic gene variant; 32% of cardiomyopathies, 26% of primary arrhythmia syndromes, and 25% of congenital heart disease. CONCLUSION: We provide important information describing the current models of care for genetic heart diseases throughout Australia. These baseline data will inform the implementation and impact of whole genome sequencing in the Australian healthcare landscape.

Validation of replacement questions for slowness and weakness to assess the Fried Phenotype: a cross-sectional study.

PURPOSE: When screening large populations, performance-based measures can be difficult to conduct because they are time consuming and costly, and require well-trained assessors. The aim of the present study is to validate a set of questions replacing the performance-based measures slowness and weakness as part of the Fried frailty phenotype (FRIED-P). METHODS: A cross-sectional study was conducted among community-dwelling older adults (≥ 60 years) in three Flemish municipalities. The Fried Phenotype (FRIED-P) was used to measure physical frailty. The two performance-based measures of the Fried Phenotype (slowness and weakness) were also measured by means of six substituting questions (FRIED-Q). These questions were validated through sensitivity, specificity, Cohen's kappa value, observed agreement, correlation analysis, and the area under the curve (AUC, ROC curve). RESULTS: 196 older adults participated. According to the FRIED-P, 19.5% of them were frail, 56.9% were pre-frail and 23.6% were non-frail. For slowness, the observed sensitivity was 47.0%, the specificity was 96.5% and the AUC was 0.717. For weakness, the sensitivity was 46.2%, the specificity was 83.7%, and the AUC was 0.649. The overall Spearman correlation between the FRIED-P and the FRIED-Q was r = 0.721 with an observed agreement of 76.6% (weighted linear kappa value = 0.663, quadratic kappa value = 0.738). CONCLUSIONS: The concordance between the FRIED-P and FRIED-Q was substantial, characterized by a very high specificity, but a moderate sensitivity. This alternative operationalization of the Fried Phenotype-i.e., including six replacement questions instead of two performance-based tests-can be considered to apply as screening tool to screen physical frailty in large populations.

Feasibility of Ultra-Rapid Exome Sequencing in Critically Ill Infants and Children With Suspected Monogenic Conditions in the Australian Public Health Care System.

Importance: Widespread adoption of rapid genomic testing in pediatric critical care requires robust clinical and laboratory pathways that provide equitable and consistent service across health care systems. Objective: To prospectively evaluate the performance of a multicenter network for ultra-rapid genomic diagnosis in a public health care system. Design, Setting, and Participants: Descriptive feasibility study of critically ill pediatric patients with suspected monogenic conditions treated at 12 Australian hospitals between March 2018 and February 2019, with data collected to May 2019. A formal implementation strategy emphasizing communication and feedback, standardized processes, coordination, distributed leadership, and collective learning was used to facilitate adoption. Exposures: Ultra-rapid exome sequencing. Main Outcomes and Measures: The primary outcome was time from sample receipt to ultra-rapid exome sequencing report. The secondary outcomes were the molecular diagnostic yield, the change in clinical management after the ultra-rapid exome sequencing report, the time from hospital admission to the laboratory report, and the proportion of laboratory reports returned prior to death or hospital discharge. Results: The study population included 108 patients with a median age of 28 days (range, 0 days to 17 years); 34% were female; and 57% were from neonatal intensive care units, 33% were from pediatric intensive care units, and 9% were from other hospital wards. The mean time from sample receipt to ultra-rapid exome sequencing report was 3.3 days (95% CI, 3.2-3.5 days) and the median time was 3 days (range, 2-7 days). The mean time from hospital admission to ultra-rapid exome sequencing report was 17.5 days (95% CI, 14.6-21.1 days) and 93 reports (86%) were issued prior to death or hospital discharge. A molecular diagnosis was established in 55 patients (51%). Eleven diagnoses (20%) resulted from using the following approaches to augment standard exome sequencing analysis: mitochondrial genome sequencing analysis, exome sequencing-based copy number analysis, use of international databases to identify novel gene-disease associations, and additional phenotyping and RNA analysis. In 42 of 55 patients (76%) with a molecular diagnosis and 6 of 53 patients (11%) without a molecular diagnosis, the ultra-rapid exome sequencing result was considered as having influenced clinical management. Targeted treatments were initiated in 12 patients (11%), treatment was redirected toward palliative care in 14 patients (13%), and surveillance for specific complications was initiated in 19 patients (18%). Conclusions and Relevance: This study suggests feasibility of ultra-rapid genomic testing in critically ill pediatric patients with suspected monogenic conditions in the Australian public health care system. However, further research is needed to understand the clinical value of such testing, and the generalizability of the findings to other health care settings.

Concordances and differences between a unidimensional and multidimensional assessment of frailty: a cross-sectional study.

BACKGROUND: Many instruments to identify frail older people have been developed. One of the consequences is that the prevalence rates of frailty vary widely dependent on the instrument selected. The aims of this study were 1) to examine the concordances and differences between a unidimensional and multidimensional assessment of frailty, 2) to assess to what extent the characteristics of a 'frail sample' differ depending on the selected frailty measurement because 'being frail' is used in many studies as an inclusion criterion. METHOD: A cross-sectional study was conducted among 196 community-dwelling older adults (≥60 years), which were selected from the census records. Unidimensional frailty was operationalized according to the Fried Phenotype (FP) and multidimensional frailty was measured with the Comprehensive Frailty Assessment Instrument (CFAI). The concordances and differences were examined by prevalence, correlations, observed agreement and Kappa values. Differences between sample characteristics (e.g., age, physical activity, life satisfaction) were investigated with ANOVA and Kruskall-Wallis test. RESULTS: The mean age was 72.74 (SD 8.04) and 48.98% was male. According to the FP 23.59% was not-frail, 56.92% pre-frail and 19.49% frail. According to the CFAI, 44.33% was no-to-low frail, 37.63% was mild frail and 18.04% was high frail. The correlation between FP and the CFAI was r = 0.46 and the observed agreement was 52.85%. The Kappa value was κ = 0.35 (quadratic κ = 0.45). In total, 11.92% of the participants were frail according to both measurements, 7.77% was solely frail according to the FP and 6.21% was solely frail according to the CFAI. The 'frail sample respondents' according to the FP had higher levels of life satisfaction and net income, but performed less physical activities in comparison to high frail people according to the CFAI. CONCLUSION: The present study shows that the FP and CFAI partly measure the same 'frailty-construct', although differences were found for instance in the prevalence of frailty and the composition of the 'frail participants'. Since 'being frail' is an inclusion criterion in many studies, researchers must be aware that the choice of the frailty measurement has an impact on both the estimates of frailty prevalence and the characteristics of the selected sample.

The Social Environment's Relationship With Frailty: Evidence From Existing Studies.

Increasingly, policymakers assume that informal networks will provide care for frail older people. While the literature has mainly discussed the role of the family, broader social networks are also considered to be important. However, these social networks can diminish in later life. This systematic review investigates whether the social environment increases the risk of frailty or helps to prevent it. Findings from 15 original studies were classified using five different factors, which denoted five dimensions of the social environment: (a) social networks, (b) social support, (c) social participation, (d) subjective neighborhood experience, and (e) socioeconomic neighborhood characteristics. The discussion highlights that the social environment and frailty are indeed related, and how the neighborhood dimensions and social participation had more consistent results than social support and social networks. Conclusively, recommendations are formulated to contemplate all dimensions of the social environment for further research examining frailty and community care.

Assessment of moderate coffee consumption and risk of epithelial ovarian cancer: a Mendelian randomization study.

Background: Coffee consumption has been shown to be associated with various health outcomes in observational studies. However, evidence for its association with epithelial ovarian cancer (EOC) is inconsistent and it is unclear whether these associations are causal. Methods: We used single nucleotide polymorphisms associated with (i) coffee and (ii) caffeine consumption to perform Mendelian randomization (MR) on EOC risk. We conducted a two-sample MR using genetic data on 44 062 individuals of European ancestry from the Ovarian Cancer Association Consortium (OCAC), and combined instrumental variable estimates using a Wald-type ratio estimator. Results: For all EOC cases, the causal odds ratio (COR) for genetically predicted consumption of one additional cup of coffee per day was 0.92 [95% confidence interval (CI): 0.79, 1.06]. The COR was 0.90 (95% CI: 0.73, 1.10) for high-grade serous EOC. The COR for genetically predicted consumption of an additional 80 mg caffeine was 1.01 (95% CI: 0.92, 1.11) for all EOC cases and 0.90 (95% CI: 0.73, 1.10) for high-grade serous cases. Conclusions: We found no evidence indicative of a strong association between EOC risk and genetically predicted coffee or caffeine levels. However, our estimates were not statistically inconsistent with earlier observational studies and we were unable to rule out small protective associations.

Analyses of germline variants associated with ovarian cancer survival identify functional candidates at the 1q22 and 19p12 outcome loci.

We previously identified associations with ovarian cancer outcome at five genetic loci. To identify putatively causal genetic variants and target genes, we prioritized two ovarian outcome loci (1q22 and 19p12) for further study. Bioinformatic and functional genetic analyses indicated that MEF2D and ZNF100 are targets of candidate outcome variants at 1q22 and 19p12, respectively. At 19p12, the chromatin interaction of a putative regulatory element with the ZNF100 promoter region correlated with candidate outcome variants. At 1q22, putative regulatory elements enhanced MEF2D promoter activity and haplotypes containing candidate outcome variants modulated these effects. In a public dataset, MEF2D and ZNF100 expression were both associated with ovarian cancer progression-free or overall survival time. In an extended set of 6,162 epithelial ovarian cancer patients, we found that functional candidates at the 1q22 and 19p12 loci, as well as other regional variants, were nominally associated with patient outcome; however, no associations reached our threshold for statistical significance (p<1×10-5). Larger patient numbers will be needed to convincingly identify any true associations at these loci.

Genome-wide association study evaluating single-nucleotide polymorphisms and outcomes in patients with advanced stage serous ovarian or primary peritoneal cancer: An NRG Oncology/Gynecologic Oncology Group study.

OBJECTIVE: This study evaluated single nucleotide polymorphisms (SNPs) associated with progression free (PFS) and overall survival (OS) in patients with advanced stage serous EOC. METHODS: Patients enrolled in GOG-172 and 182 who provided specimens for translational research and consent were included. Germline DNA was evaluated with the Illumina's HumanOMNI1-Quad beadchips and scanned using Illumina's iScan optical imaging system. SNPs with allele frequency>0.05 and genotyping rate>0.98 were included. Analysis of SNPs for PFS and OS was done using Cox regression. Statistical significance was determined using Bonferroni corrected p-values with genomic control adjustment. RESULTS: The initial GWAS analysis included 1,124,677 markers in 396 patients. To obtain the final data set, quality control checks were performed and limited to serous tumors and self-identified Caucasian race. In total 636,555 SNPs and 289 patients passed all the filters. The pre-specified statistical level of significance was 7.855e-08. No SNPs met this criteria for PFS or OS, however, two SNPs were close to significance (rs10899426 p-2.144e-08) (rs6256 p-9.774e-07) for PFS and 2 different SNPs were identified (rs295315 p-7.536e-07; rs17693104 p-7.734e-07) which were close to significance for OS. CONCLUSIONS: Using the pre-specified level of significance of 1×10-08, we did not identify any SNPs of statistical significance for OS or PFS, however several were close. The SNP's identified in this GWAS study will require validation and these preliminary findings may lead to identification of novel pathways and biomarkers.

Corrigendum: Whole-genome landscape of pancreatic neuroendocrine tumours.

This corrects the article DOI: 10.1038/nature21063.

Whole-genome landscape of pancreatic neuroendocrine tumours.

The diagnosis of pancreatic neuroendocrine tumours (PanNETs) is increasing owing to more sensitive detection methods, and this increase is creating challenges for clinical management. We performed whole-genome sequencing of 102 primary PanNETs and defined the genomic events that characterize their pathogenesis. Here we describe the mutational signatures they harbour, including a deficiency in G:C > T:A base excision repair due to inactivation of MUTYH, which encodes a DNA glycosylase. Clinically sporadic PanNETs contain a larger-than-expected proportion of germline mutations, including previously unreported mutations in the DNA repair genes MUTYH, CHEK2 and BRCA2. Together with mutations in MEN1 and VHL, these mutations occur in 17% of patients. Somatic mutations, including point mutations and gene fusions, were commonly found in genes involved in four main pathways: chromatin remodelling, DNA damage repair, activation of mTOR signalling (including previously undescribed EWSR1 gene fusions), and telomere maintenance. In addition, our gene expression analyses identified a subgroup of tumours associated with hypoxia and HIF signalling.

Hypermutation In Pancreatic Cancer.

Pancreatic cancer is molecularly diverse, with few effective therapies. Increased mutation burden and defective DNA repair are associated with response to immune checkpoint inhibitors in several other cancer types. We interrogated 385 pancreatic cancer genomes to define hypermutation and its causes. Mutational signatures inferring defects in DNA repair were enriched in those with the highest mutation burdens. Mismatch repair deficiency was identified in 1% of tumors harboring different mechanisms of somatic inactivation of MLH1 and MSH2. Defining mutation load in individual pancreatic cancers and the optimal assay for patient selection may inform clinical trial design for immunotherapy in pancreatic cancer.