Genome sequencing detects a wide range of clinically relevant copy-number variants and other genomic alterations.

PURPOSE: Copy-number variants (CNVs) and other non-single nucleotide variant/indel variant types contribute an important proportion of diagnoses in individuals with suspected genetic disease. This study describes the range of such variants detected by genome sequencing (GS). METHODS: For a pediatric cohort of 1032 participants undergoing clinical GS, we characterize the CNVs and other non-single nucleotide variant/indel variant types that were reported, including aneuploidies, mobile element insertions, and uniparental disomies, and we describe the bioinformatic pipeline used to detect these variants. RESULTS: Together, these genetic alterations accounted for 15.8% of reported variants. Notably, 67.9% of these were deletions, 32.9% of which overlapped a single gene, and many deletions were reported together with a second variant in the same gene in cases of recessive disease. A retrospective medical record review in a subset of this cohort revealed that up to 6 additional genetic tests were ordered in 68% (26/38) of cases, some of which failed to report the CNVs/rare variants reported on GS. CONCLUSION: GS detected a broad range of reported variant types, including CNVs ranging in size from 1 Kb to 46 Mb.

Rapid Whole-Genomic Sequencing and a Targeted Neonatal Gene Panel in Infants With a Suspected Genetic Disorder.

Importance: Genomic testing in infancy guides medical decisions and can improve health outcomes. However, it is unclear whether genomic sequencing or a targeted neonatal gene-sequencing test provides comparable molecular diagnostic yields and times to return of results. Objective: To compare outcomes of genomic sequencing with those of a targeted neonatal gene-sequencing test. Design, Setting, and Participants: The Genomic Medicine for Ill Neonates and Infants (GEMINI) study was a prospective, comparative, multicenter study of 400 hospitalized infants younger than 1 year of age (proband) and their parents, when available, suspected of having a genetic disorder. The study was conducted at 6 US hospitals from June 2019 to November 2021. Exposure: Enrolled participants underwent simultaneous testing with genomic sequencing and a targeted neonatal gene-sequencing test. Each laboratory performed an independent interpretation of variants guided by knowledge of the patient's phenotype and returned results to the clinical care team. Change in clinical management, therapies offered, and redirection of care was provided to families based on genetic findings from either platform. Main Outcomes and Measures: Primary end points were molecular diagnostic yield (participants with ≥1 pathogenic variant or variant of unknown significance), time to return of results, and clinical utility (changes in patient care). Results: A molecular diagnostic variant was identified in 51% of participants (n = 204; 297 variants identified with 134 being novel). Molecular diagnostic yield of genomic sequencing was 49% (95% CI, 44%-54%) vs 27% (95% CI, 23%-32%) with the targeted gene-sequencing test. Genomic sequencing did not report 19 variants found by the targeted neonatal gene-sequencing test; the targeted gene-sequencing test did not report 164 variants identified by genomic sequencing as diagnostic. Variants unidentified by the targeted genomic-sequencing test included structural variants longer than 1 kilobase (25.1%) and genes excluded from the test (24.6%) (McNemar odds ratio, 8.6 [95% CI, 5.4-14.7]). Variant interpretation by laboratories differed by 43%. Median time to return of results was 6.1 days for genomic sequencing and 4.2 days for the targeted genomic-sequencing test; for urgent cases (n = 107) the time was 3.3 days for genomic sequencing and 4.0 days for the targeted gene-sequencing test. Changes in clinical care affected 19% of participants, and 76% of clinicians viewed genomic testing as useful or very useful in clinical decision-making, irrespective of a diagnosis. Conclusions and Relevance: The molecular diagnostic yield for genomic sequencing was higher than a targeted neonatal gene-sequencing test, but the time to return of routine results was slower. Interlaboratory variant interpretation contributes to differences in molecular diagnostic yield and may have important consequences for clinical management.

Rapid Whole Genome Sequencing for Diagnosis of Single Locus Genetic Diseases in Critically Ill Children.

Upon admission to intensive care units (ICU), the differential diagnosis of almost all infants with diseases of unclear etiology includes single locus genetic diseases. Rapid whole genome sequencing (rWGS), including sample preparation, short-read sequencing-by-synthesis, informatics pipelining, and semiautomated interpretation, can now identify nucleotide and structural variants associated with most genetic diseases with robust analytic and diagnostic performance in as little as 13.5 h. Early diagnosis of genetic diseases transforms medical and surgical management of infants in ICUs, minimizing both the duration of empiric treatment and the delay to start of specific treatment. Both positive and negative rWGS tests have clinical utility and can improve outcomes. Since first described 10 years ago, rWGS has evolved considerably. Here we describe our current methods for routine diagnostic testing for genetic diseases by rWGS in as little as 18 h.

A genome sequencing system for universal newborn screening, diagnosis, and precision medicine for severe genetic diseases.

Newborn screening (NBS) dramatically improves outcomes in severe childhood disorders by treatment before symptom onset. In many genetic diseases, however, outcomes remain poor because NBS has lagged behind drug development. Rapid whole-genome sequencing (rWGS) is attractive for comprehensive NBS because it concomitantly examines almost all genetic diseases and is gaining acceptance for genetic disease diagnosis in ill newborns. We describe prototypic methods for scalable, parentally consented, feedback-informed NBS and diagnosis of genetic diseases by rWGS and virtual, acute management guidance (NBS-rWGS). Using established criteria and the Delphi method, we reviewed 457 genetic diseases for NBS-rWGS, retaining 388 (85%) with effective treatments. Simulated NBS-rWGS in 454,707 UK Biobank subjects with 29,865 pathogenic or likely pathogenic variants associated with 388 disorders had a true negative rate (specificity) of 99.7% following root cause analysis. In 2,208 critically ill children with suspected genetic disorders and 2,168 of their parents, simulated NBS-rWGS for 388 disorders identified 104 (87%) of 119 diagnoses previously made by rWGS and 15 findings not previously reported (NBS-rWGS negative predictive value 99.6%, true positive rate [sensitivity] 88.8%). Retrospective NBS-rWGS diagnosed 15 children with disorders that had been undetected by conventional NBS. In 43 of the 104 children, had NBS-rWGS-based interventions been started on day of life 5, the Delphi consensus was that symptoms could have been avoided completely in seven critically ill children, mostly in 21, and partially in 13. We invite groups worldwide to refine these NBS-rWGS conditions and join us to prospectively examine clinical utility and cost effectiveness.

An automated 13.5 hour system for scalable diagnosis and acute management guidance for genetic diseases.

While many genetic diseases have effective treatments, they frequently progress rapidly to severe morbidity or mortality if those treatments are not implemented immediately. Since front-line physicians frequently lack familiarity with these diseases, timely molecular diagnosis may not improve outcomes. Herein we describe Genome-to-Treatment, an automated, virtual system for genetic disease diagnosis and acute management guidance. Diagnosis is achieved in 13.5 h by expedited whole genome sequencing, with superior analytic performance for structural and copy number variants. An expert panel adjudicated the indications, contraindications, efficacy, and evidence-of-efficacy of 9911 drug, device, dietary, and surgical interventions for 563 severe, childhood, genetic diseases. The 421 (75%) diseases and 1527 (15%) effective interventions retained are integrated with 13 genetic disease information resources and appended to diagnostic reports ( https://gtrx.radygenomiclab.com ). This system provided correct diagnoses in four retrospectively and two prospectively tested infants. The Genome-to-Treatment system facilitates optimal outcomes in children with rapidly progressive genetic diseases.

Best practices for the interpretation and reporting of clinical whole genome sequencing.

Whole genome sequencing (WGS) shows promise as a first-tier diagnostic test for patients with rare genetic disorders. However, standards addressing the definition and deployment practice of a best-in-class test are lacking. To address these gaps, the Medical Genome Initiative, a consortium of leading health care and research organizations in the US and Canada, was formed to expand access to high quality clinical WGS by convening experts and publishing best practices. Here, we present best practice recommendations for the interpretation and reporting of clinical diagnostic WGS, including discussion of challenges and emerging approaches that will be critical to harness the full potential of this comprehensive test.

Expanding the phenotypic and molecular spectrum of NFS1-related disorders that cause functional deficiencies in mitochondrial and cytosolic iron-sulfur cluster containing enzymes.

Iron-sulfur cluster proteins are involved in critical functions for gene expression regulation and mitochondrial bioenergetics including the oxidative phosphorylation system. The c.215G>A p.(Arg72Gln) variant in NFS1 has been previously reported to cause infantile mitochondrial complex II and III deficiency. We describe three additional unrelated patients with the same missense variant. Two infants with the same homozygous variant presented with hypotonia, weakness and lactic acidosis, and one patient with compound heterozygous p.(Arg72Gln) and p.(Arg412His) variants presented as a young adult with gastrointestinal symptoms and fatigue. Skeletal muscle biopsy from patients 1 and 3 showed abnormal mitochondrial morphology, and functional analyses demonstrated decreased activity in respiratory chain complex II and variably in complexes I and III. We found decreased mitochondrial and cytosolic aconitase activities but only mildly affected lipoylation of pyruvate dehydrogenase and 2-oxoglutarate dehydrogenase enzymes. Our studies expand the phenotypic spectrum and provide further evidence for the pathogenicity and functional sequelae of NFS1-related disorders with disturbances in both mitochondrial and cytosolic iron-sulfur cluster containing enzymes.

Uncertain, Not Unimportant: Callosal Dysgenesis and Variants of Uncertain Significance in ROBO1.

Congenital anomalies affect 3% to 5% of births and remain the leading cause of infant death in the United States. As whole exome and genome sequencing are increasingly used to diagnose underlying genetic disease, the patient's clinical presentation remains the most important context for interpreting sequencing results, including frequently reported variants of uncertain significance (VUS). Classification of a variant as VUS acknowledges limits on evidence to establish whether a variant can be classified as pathogenic or benign according to the American College of Medical Genetics guidelines. Importantly, the VUS designation reflects limits on the breadth of evidence linking the genetic variant to a disease. However, available evidence, although limited, may be surprisingly relevant in an individual patient's case. Accordingly, a VUS result should be approached neither as nondiagnostic genetic result nor as automatically "uncertain" in its potential to guide clinical decision-making. In this article, we discuss a case of an infant born at 29 weeks 4 days without a corpus callosum, whose whole genome sequencing yielded compound heterozygous variants both classified as VUS in ROBO1, a gene encoding for a receptor involved in a canonical signaling mechanism that guides axons across midline. Approaching the VUS result as potentially pathogenic, we found the infant ultimately had pituitary dysfunction and renal anomalies consistent with other reported ROBO1 variants and basic science literature. Accordingly, we highlight resources for variant interpretation available to clinicians to evaluate VUS results, particularly as they inform the diagnosis of individually rare but collectively common rare diseases.

Project Baby Bear: Rapid precision care incorporating rWGS in 5 California children's hospitals demonstrates improved clinical outcomes and reduced costs of care.

Genetic disorders are a leading contributor to mortality in neonatal and pediatric intensive care units (ICUs). Rapid whole-genome sequencing (rWGS)-based rapid precision medicine (RPM) is an intervention that has demonstrated improved clinical outcomes and reduced costs of care. However, the feasibility of broad clinical deployment has not been established. The objective of this study was to implement RPM based on rWGS and evaluate the clinical and economic impact of this implementation as a first line diagnostic test in the California Medicaid (Medi-Cal) program. Project Baby Bear was a payor funded, prospective, real-world quality improvement project in the regional ICUs of five tertiary care children's hospitals. Participation was limited to acutely ill Medi-Cal beneficiaries who were admitted November 2018 to May 2020, were <1 year old and within one week of hospitalization, or had just developed an abnormal response to therapy. The whole cohort received RPM. There were two prespecified primary outcomes-changes in medical care reported by physicians and changes in the cost of care. The majority of infants were from underserved populations. Of 184 infants enrolled, 74 (40%) received a diagnosis by rWGS that explained their admission in a median time of 3 days. In 58 (32%) affected individuals, rWGS led to changes in medical care. Testing and precision medicine cost $1.7 million and led to $2.2-2.9 million cost savings. rWGS-based RPM had clinical utility and reduced net health care expenditures for infants in regional ICUs. rWGS should be considered early in ICU admission when the underlying etiology is unclear.

Expanding the genotypic spectrum of ACTG2-related visceral myopathy.

Visceral myopathies (VMs) encompass a spectrum of disorders characterized by chronic disruption of gastrointestinal function, with or without urinary system involvement. Pathogenic missense variation in smooth muscle γ-actin gene (ACTG2) is associated with autosomal dominant VM. Whole-genome sequencing of an infant presenting with chronic intestinal pseudo-obstruction revealed a homozygous 187 bp (c.589_613 + 163del188) deletion spanning the exon 6-intron 6 boundary within ACTG2 The patient's clinical course was marked by prolonged hospitalizations, multiple surgeries, and intermittent total parenteral nutrition dependence. This case supports the emerging understanding of allelic heterogeneity in ACTG2-related VM, in which both biallelic and monoallelic variants in ACTG2 are associated with gastrointestinal dysfunction of similar severity and overlapped clinical presentation. Moreover, it illustrates the clinical utility of rapid whole-genome sequencing, which can comprehensively and precisely detect different types of genomic variants including small deletions, leading to guidance of clinical care decisions.

Hemimegalencephaly and intractable seizures associated with the NPRL3 gene variant in a newborn: A case report.

Hemimegalencephaly (HME) is a rare hamartomatous congenital malformation of the brain characterized by dysplastic overgrowth of either one of the cerebral hemispheres. HME is associated with early onset seizures, abnormal neurological findings, and with subsequent cognitive and behavioral disabilities. Seizures associated with HME are often refractory to antiepileptic medications. Hemispherectomy is usually necessary to provide effective seizure control. The exact etiology of HME is not fully understood, but involves a disturbance in early brain development and likely involves genes responsible for patterning and symmetry of the brain. We present a female newborn who had refractory seizures due to HME. Whole genome sequencing revealed a novel, likely pathogenic, maternally inherited, 3Kb deletion encompassing exon 5 of the NPRL3 gene (chr16:161898-164745x1). The NPRL3 gene encodes for a nitrogen permease regulator 3-like protein, a subunit of the GATOR complex, which regulates the mTOR signaling pathway. A trial of mTOR inhibitor drug, Sirolimus, did not improve her seizure control. Functional hemispherectomy at 3 months of age resulted in total abatement of clinical seizures.

ECHS1 disease in two unrelated families of Samoan descent: Common variant - rare disorder.

Mutations in the short-chain enoyl-CoA hydratase (SCEH) gene, ECHS1, cause a rare autosomal recessive disorder of valine catabolism. Patients usually present with developmental delay, regression, dystonia, feeding difficulties, and abnormal MRI with bilateral basal ganglia involvement. We present clinical, biochemical, molecular, and functional data for four affected patients from two unrelated families of Samoan descent with identical novel compound heterozygous mutations. Family 1 has three affected boys while Family 2 has an affected daughter, all with clinical and MRI findings of Leigh syndrome and intermittent episodes of acidosis and ketosis. WES identified a single heterozygous variant in ECHS1 at position c.832G > A (p.Ala278Thr). However, western blot revealed significantly reduced ECHS1 protein for all affected family members. Decreased SCEH activity in fibroblasts and a mild increase in marker metabolites in urine further supported ECHS1 as the underlying gene defect. Additional investigations at the DNA (aCGH, WGS) and RNA (qPCR, RT-PCR, RNA-Seq, RNA-Array) level identified a silent, common variant at position c.489G > A (p.Pro163=) as the second mutation. This substitution, present at high frequency in the Samoan population, is associated with decreased levels of normally spliced mRNA. To our understanding, this is the first report of a novel, hypomorphic allele c.489G > A (p.Pro163=), associated with SCEH deficiency.

Clinical utility of ultra-rapid whole-genome sequencing in an infant with atypical presentation of WT1-associated nephrotic syndrome type 4.

Relatively little is known about phenotypic variability in nonsyndromic nephropathy associated with the gene encoding the WT1 transcription factor. We report a 12-mo-old female who presented with vomiting, diarrhea, and fatigue in the setting of renal failure and malignant hypertension. Trio ultra-rapid whole-genome sequencing identified a novel, likely pathogenic, de novo missense variant (c.485T > A, p.Val162Asp) in WT1 in 46 h, consistent with a diagnosis of nephrotic syndrome type 4 (NPHS4; OMIM 256370). This disorder typically presents with nephrotic syndrome (gross proteinuria, hypoalbuminemia, and edema). Rapid diagnosis had an immediate impact on her clinical management in the pediatric intensive care unit. Diagnostic renal biopsy was avoided, and placement of permanent dialysis access, a gastrostomy tube, and bilateral nephrectomy were accelerated. This report expands the presenting phenotype of nonsyndromic nephrotic syndrome and/or renal failure due to heterozygous variants in WT1 (NPHS4). It also highlights the relationship between time to genomic diagnosis and clinical utility in critically ill infants.

Genotype-phenotype correlation at codon 1740 of SETD2.

The SET domain containing 2, histone lysine methyltransferase encoded by SETD2 is a dual-function methyltransferase for histones and microtubules and plays an important role for transcriptional regulation, genomic stability, and cytoskeletal functions. Specifically, SETD2 is associated with trimethylation of histone H3 at lysine 36 (H3K36me3) and methylation of α-tubulin at lysine 40. Heterozygous loss of function and missense variants have previously been described with Luscan-Lumish syndrome (LLS), which is characterized by overgrowth, neurodevelopmental features, and absence of overt congenital anomalies. We have identified 15 individuals with de novo variants in codon 1740 of SETD2 whose features differ from those with LLS. Group 1 consists of 12 individuals with heterozygous variant c.5218C>T p.(Arg1740Trp) and Group 2 consists of 3 individuals with heterozygous variant c.5219G>A p.(Arg1740Gln). The phenotype of Group 1 includes microcephaly, profound intellectual disability, congenital anomalies affecting several organ systems, and similar facial features. Individuals in Group 2 had moderate to severe intellectual disability, low normal head circumference, and absence of additional major congenital anomalies. While LLS is likely due to loss of function of SETD2, the clinical features seen in individuals with variants affecting codon 1740 are more severe suggesting an alternative mechanism, such as gain of function, effects on epigenetic regulation, or posttranslational modification of the cytoskeleton. Our report is a prime example of different mutations in the same gene causing diverging phenotypes and the features observed in Group 1 suggest a new clinically recognizable syndrome uniquely associated with the heterozygous variant c.5218C>T p.(Arg1740Trp) in SETD2.

A Randomized, Controlled Trial of the Analytic and Diagnostic Performance of Singleton and Trio, Rapid Genome and Exome Sequencing in Ill Infants.

The second Newborn Sequencing in Genomic Medicine and Public Health study was a randomized, controlled trial of the effectiveness of rapid whole-genome or -exome sequencing (rWGS or rWES, respectively) in seriously ill infants with diseases of unknown etiology. Here we report comparisons of analytic and diagnostic performance. Of 1,248 ill inpatient infants, 578 (46%) had diseases of unknown etiology. 213 infants (37% of those eligible) were enrolled within 96 h of admission. 24 infants (11%) were very ill and received ultra-rapid whole-genome sequencing (urWGS). The remaining infants were randomized, 95 to rWES and 94 to rWGS. The analytic performance of rWGS was superior to rWES, including variants likely to affect protein function, and ClinVar pathogenic/likely pathogenic variants (p < 0.0001). The diagnostic performance of rWGS and rWES were similar (18 diagnoses in 94 infants [19%] versus 19 diagnoses in 95 infants [20%], respectively), as was time to result (median 11.0 versus 11.2 days, respectively). However, the proportion diagnosed by urWGS (11 of 24 [46%]) was higher than rWES/rWGS (p = 0.004) and time to result was less (median 4.6 days, p < 0.0001). The incremental diagnostic yield of reflexing to trio after negative proband analysis was 0.7% (1 of 147). In conclusion, rapid genomic sequencing can be performed as a first-tier diagnostic test in inpatient infants. urWGS had the shortest time to result, which was important in unstable infants, and those in whom a genetic diagnosis was likely to impact immediate management. Further comparison of urWGS and rWES is warranted because genomic technologies and knowledge of variant pathogenicity are evolving rapidly.

Diagnosis of genetic diseases in seriously ill children by rapid whole-genome sequencing and automated phenotyping and interpretation.

By informing timely targeted treatments, rapid whole-genome sequencing can improve the outcomes of seriously ill children with genetic diseases, particularly infants in neonatal and pediatric intensive care units (ICUs). The need for highly qualified professionals to decipher results, however, precludes widespread implementation. We describe a platform for population-scale, provisional diagnosis of genetic diseases with automated phenotyping and interpretation. Genome sequencing was expedited by bead-based genome library preparation directly from blood samples and sequencing of paired 100-nt reads in 15.5 hours. Clinical natural language processing (CNLP) automatically extracted children's deep phenomes from electronic health records with 80% precision and 93% recall. In 101 children with 105 genetic diseases, a mean of 4.3 CNLP-extracted phenotypic features matched the expected phenotypic features of those diseases, compared with a match of 0.9 phenotypic features used in manual interpretation. We automated provisional diagnosis by combining the ranking of the similarity of a patient's CNLP phenome with respect to the expected phenotypic features of all genetic diseases, together with the ranking of the pathogenicity of all of the patient's genomic variants. Automated, retrospective diagnoses concurred well with expert manual interpretation (97% recall and 99% precision in 95 children with 97 genetic diseases). Prospectively, our platform correctly diagnosed three of seven seriously ill ICU infants (100% precision and recall) with a mean time saving of 22:19 hours. In each case, the diagnosis affected treatment. Genome sequencing with automated phenotyping and interpretation in a median of 20:10 hours may increase adoption in ICUs and, thereby, timely implementation of precise treatments.

Measurement of Elevated Concentrations of Urine Keratan Sulfate by UPLC-MSMS in Lysosomal Storage Disorders (LSDs): Comparison of Urine Keratan Sulfate Levels in MPS IVA Versus Other LSDs.

Keratan sulfate (KS) is commonly elevated in urine samples from patients with mucopolysaccharidosis type IVA (MPS IVA) and is considered pathognomonic for the condition. Recently, a new method has been described by Martell et al. to detect and measure urinary KS utilizing LC-MS/MS. As a part of the validation of this method in our laboratory, we studied the sensitivity and specificity of elevated urine KS levels using 25 samples from 15 MPS IVA patients, and 138 samples from 102 patients with other lysosomal storage disorders, including MPS I (n = 9), MPS II (n = 13), MPS III (n = 23), MPS VI (n = 7), beta-galactosidase deficiency (n = 7), mucolipidosis (ML) type II, II/III and III (n = 51), alpha-mannosidosis (n = 11), fucosidosis (n = 4), sialidosis (n = 5), Pompe disease (n = 3), aspartylglucosaminuria (n = 4), and galactosialidosis (n = 1). As expected, urine KS values were significantly higher (fivefold average increase) than age-matched controls in all MPS IVA patients. Urine KS levels were also significantly elevated (threefold to fourfold increase) in patients with GM-1 gangliosidosis, MPS IVB, ML II and ML II/III, and fucosidosis. Urine KS was also elevated to a smaller degree (1.1-fold to 1.7-fold average increase) in patients with MPS I, MPS II, and ML III. These findings suggest that while the UPLC-MS/MS urine KS method is 100% sensitive for the detection of patients with MPS IVA, elevated urine KS is not specific for this condition. Therefore, caution is advised when interpreting urinary keratan sulfate results.

Contribution of FKBP5 genetic variation to gemcitabine treatment and survival in pancreatic adenocarcinoma.

PURPOSE: FKBP51, (FKBP5), is a negative regulator of Akt. Variability in FKBP5 expression level is a major factor contributing to variation in response to chemotherapeutic agents including gemcitabine, a first line treatment for pancreatic cancer. Genetic variation in FKBP5 could influence its function and, ultimately, treatment response of pancreatic cancer. EXPERIMENTAL DESIGN: We set out to comprehensively study the role of genetic variation in FKBP5 identified by Next Generation DNA resequencing on response to gemcitabine treatment of pancreatic cancer by utilizing both tumor and germline DNA samples from 43 pancreatic cancer patients, including 19 paired normal-tumor samples. Next, genotype-phenotype association studies were performed with overall survival as well as with FKBP5 gene expression in tumor using the same samples in which resequencing had been performed, followed by functional genomics studies. RESULTS: In-depth resequencing identified 404 FKBP5 single nucleotide polymorphisms (SNPs) in normal and tumor DNA. SNPs with the strongest associations with survival or FKBP5 expression were subjected to functional genomic study. Electromobility shift assay showed that the rs73748206 "A(T)" SNP altered DNA-protein binding patterns, consistent with significantly increased reporter gene activity, possibly through its increased binding to Glucocorticoid Receptor (GR). The effect of rs73748206 was confirmed on the basis of its association with FKBP5 expression by affecting the binding to GR in lymphoblastoid cell lines derived from the same patients for whom DNA was used for resequencing. CONCLUSION: This comprehensive FKBP5 resequencing study provides insights into the role of genetic variation in variation of gemcitabine response.