Closing the gap: Solving complex medically relevant genes at scale.

Comprehending the mechanism behind human diseases with an established heritable component represents the forefront of personalized medicine. Nevertheless, numerous medically important genes are inaccurately represented in short-read sequencing data analysis due to their complexity and repetitiveness or the so-called 'dark regions' of the human genome. The advent of PacBio as a long-read platform has provided new insights, yet HiFi whole-genome sequencing (WGS) cost remains frequently prohibitive. We introduce a targeted sequencing and analysis framework, Twist Alliance Dark Genes Panel (TADGP), designed to offer phased variants across 389 medically important yet complex autosomal genes. We highlight TADGP accuracy across eleven control samples and compare it to WGS. This demonstrates that TADGP achieves variant calling accuracy comparable to HiFi-WGS data, but at a fraction of the cost. Thus, enabling scalability and broad applicability for studying rare diseases or complementing previously sequenced samples to gain insights into these complex genes. TADGP revealed several candidate variants across all cases and provided insight into LPA diversity when tested on samples from rare disease and cardiovascular disease cohorts. In both cohorts, we identified novel variants affecting individual disease-associated genes (e.g., IKZF1, KCNE1). Nevertheless, the annotation of the variants across these 389 medically important genes remains challenging due to their underrepresentation in ClinVar and gnomAD. Consequently, we also offer an annotation resource to enhance the evaluation and prioritization of these variants. Overall, we can demonstrate that TADGP offers a cost-efficient and scalable approach to routinely assess the dark regions of the human genome with clinical relevance.

A transcription factor ZmGLK36 confers broad resistance to maize rough dwarf disease in cereal crops.

Maize rough dwarf disease (MRDD), caused by maize rough dwarf virus (MRDV) or rice black-streaked dwarf virus (RBSDV), seriously threatens worldwide production of all major cereal crops, including maize, rice, wheat and barley. Here we report fine mapping and cloning of a previously reported major quantitative trait locus (QTL) (qMrdd2) for RBSDV resistance in maize. Subsequently, we show that qMrdd2 encodes a G2-like transcription factor named ZmGLK36 that promotes resistance to RBSDV by enhancing jasmonic acid (JA) biosynthesis and JA-mediated defence response. We identify a 26-bp indel located in the 5' UTR of ZmGLK36 that contributes to differential expression and resistance to RBSDV in maize inbred lines. Moreover, we show that ZmDBF2, an AP2/EREBP family transcription factor, directly binds to the 26-bp indel and represses ZmGLK36 expression. We further demonstrate that ZmGLK36 plays a conserved role in conferring resistance to RBSDV in rice and wheat using transgenic or marker-assisted breeding approaches. Our results provide insights into the molecular mechanisms of RBSDV resistance and effective strategies to breed RBSDV-resistant cereal crops.

Functional Genomics of Gastrointestinal Escherichia coli Isolated from Patients with Cancer and Diarrhea.

We describe the epidemiology and clinical characteristics of 29 patients with cancer and diarrhea in whom Enteroaggregative Escherichia coli (EAEC) was initially identified by GI BioFire panel multiplex. E. coli strains were successfully isolated from fecal cultures in 14 of 29 patients. Six of the 14 strains were identified as EAEC and 8 belonged to other diverse E. coli groups of unknown pathogenesis. We investigated these strains by their adherence to human intestinal organoids, cytotoxic responses, antibiotic resistance profile, full sequencing of their genomes, and annotation of their functional virulome. Interestingly, we discovered novel and enhanced adherence and aggregative patterns for several diarrheagenic pathotypes that were not previously seen when co-cultured with immortalized cell lines. EAEC isolates displayed exceptional adherence and aggregation to human colonoids compared not only to diverse GI E. coli , but also compared to prototype strains of other diarrheagenic E. coli . Some of the diverse E. coli strains that could not be classified as a conventional pathotype also showed an enhanced aggregative and cytotoxic response. Notably, we found a high carriage rate of antibiotic resistance genes in both EAEC strains and diverse GI E. coli isolates and observed a positive correlation between adherence to colonoids and the number of metal acquisition genes carried in both EAEC and the diverse E. coli strains. This work indicates that E. coli from cancer patients constitute strains of remarkable pathotypic and genomic divergence, including strains of unknown disease etiology with unique virulomes. Future studies will allow for the opportunity to re-define E. coli pathotypes with greater diagnostic accuracy and into more clinically relevant groupings.

Evaluating the Genetic Background Effect on Dissecting the Genetic Basis of Kernel Traits in Reciprocal Maize Introgression Lines.

Maize yield is mostly determined by its grain size. Although numerous quantitative trait loci (QTL) have been identified for kernel-related traits, the application of these QTL in breeding programs has been strongly hindered because the populations used for QTL mapping are often different from breeding populations. However, the effect of genetic background on the efficiency of QTL and the accuracy of trait genomic prediction has not been fully studied. Here, we used a set of reciprocal introgression lines (ILs) derived from 417F × 517F to evaluate how genetic background affects the detection of QTLassociated with kernel shape traits. A total of 51 QTL for kernel size were identified by chromosome segment lines (CSL) and genome-wide association studies (GWAS) methods. These were subsequently clustered into 13 common QTL based on their physical position, including 7 genetic-background-independent and 6 genetic-background-dependent QTL, respectively. Additionally, different digenic epistatic marker pairs were identified in the 417F and 517F ILs. Therefore, our results demonstrated that genetic background strongly affected not only the kernel size QTL mapping via CSL and GWAS but also the genomic prediction accuracy and epistatic detection, thereby enhancing our understanding of how genetic background affects the genetic dissection of grain size-related traits.

Exome-wide assessment of isolated biliary atresia: A report from the National Birth Defects Prevention Study using child-parent trios and a case-control design to identify novel rare variants.

The etiology of biliary atresia (BA) is unknown, but recent studies suggest a role for rare protein-altering variants (PAVs). Exome sequencing data from the National Birth Defects Prevention Study on 54 child-parent trios, one child-mother duo, and 1513 parents of children with other birth defects were analyzed. Most (91%) cases were isolated BA. We performed (1) a trio-based analysis to identify rare de novo, homozygous, and compound heterozygous PAVs and (2) a case-control analysis using a sequence kernel-based association test to identify genes enriched with rare PAVs. While we replicated previous findings on PKD1L1, our results do not suggest that recurrent de novo PAVs play important roles in BA susceptibility. In fact, our finding in NOTCH2, a disease gene associated with Alagille syndrome, highlights the difficulty in BA diagnosis. Notably, IFRD2 has been implicated in other gastrointestinal conditions and warrants additional study. Overall, our findings strengthen the hypothesis that the etiology of BA is complex.

Structural variation across 138,134 samples in the TOPMed consortium.

Ever larger Structural Variant (SV) catalogs highlighting the diversity within and between populations help researchers better understand the links between SVs and disease. The identification of SVs from DNA sequence data is non-trivial and requires a balance between comprehensiveness and precision. Here we present a catalog of 355,667 SVs (59.34% novel) across autosomes and the X chromosome (50bp+) from 138,134 individuals in the diverse TOPMed consortium. We describe our methodologies for SV inference resulting in high variant quality and >90% allele concordance compared to long-read de-novo assemblies of well-characterized control samples. We demonstrate utility through significant associations between SVs and important various cardio-metabolic and hematologic traits. We have identified 690 SV hotspots and deserts and those that potentially impact the regulation of medically relevant genes. This catalog characterizes SVs across multiple populations and will serve as a valuable tool to understand the impact of SV on disease development and progression.

Structural variation across 138,134 samples in the TOPMed consortium.

Ever larger Structural Variant (SV) catalogs highlighting the diversity within and between populations help researchers better understand the links between SVs and disease. The identification of SVs from DNA sequence data is non-trivial and requires a balance between comprehensiveness and precision. Here we present a catalog of 355,667 SVs (59.34% novel) across autosomes and the X chromosome (50bp+) from 138,134 individuals in the diverse TOPMed consortium. We describe our methodologies for SV inference resulting in high variant quality and >90% allele concordance compared to long-read de-novo assemblies of well-characterized control samples. We demonstrate utility through significant associations between SVs and important various cardio-metabolic and hemotologic traits. We have identified 690 SV hotspots and deserts and those that potentially impact the regulation of medically relevant genes. This catalog characterizes SVs across multiple populations and will serve as a valuable tool to understand the impact of SV on disease development and progression.

Long read sequencing and expression studies of AHDC1 deletions in Xia-Gibbs syndrome reveal a novel genetic regulatory mechanism.

Xia-Gibbs syndrome (XGS; MIM# 615829) is a rare mendelian disorder characterized by Development Delay (DD), intellectual disability (ID), and hypotonia. Individuals with XGS typically harbor de novo protein-truncating mutations in the AT-Hook DNA binding motif containing 1 (AHDC1) gene, although some missense mutations can also cause XGS. Large de novo heterozygous deletions that encompass the AHDC1 gene have also been ascribed as diagnostic for the disorder, without substantial evidence to support their pathogenicity. We analyzed 19 individuals with large contiguous deletions involving AHDC1, along with other genes. One individual bore the smallest known contiguous AHDC1 deletion (∼350 Kb), encompassing eight other genes within chr1p36.11 (Feline Gardner-Rasheed, IFI6, FAM76A, STX12, PPP1R8, THEMIS2, RPA2, SMPDL3B) and terminating within the first intron of AHDC1. The breakpoint junctions and phase of the deletion were identified using both short and long read sequencing (Oxford Nanopore). Quantification of RNA expression patterns in whole blood revealed that AHDC1 exhibited a mono-allelic expression pattern with no deficiency in overall AHDC1 expression levels, in contrast to the other deleted genes, which exhibited a 50% reduction in mRNA expression. These results suggest that AHDC1 expression in this individual is compensated by a novel regulatory mechanism and advances understanding of mutational and regulatory mechanisms in neurodevelopmental disorders.

Construction of a new chromosome-scale, long-read reference genome assembly for the Syrian hamster, Mesocricetus auratus.

BACKGROUND: The Syrian hamster (Mesocricetus auratus) has been suggested as a useful mammalian model for a variety of diseases and infections, including infection with respiratory viruses such as SARS-CoV-2. The MesAur1.0 genome assembly was generated in 2013 using whole-genome shotgun sequencing with short-read sequence data. Current more advanced sequencing technologies and assembly methods now permit the generation of near-complete genome assemblies with higher quality and greater continuity. FINDINGS: Here, we report an improved assembly of the M. auratus genome (BCM_Maur_2.0) using Oxford Nanopore Technologies long-read sequencing to produce a chromosome-scale assembly. The total length of the new assembly is 2.46 Gb, similar to the 2.50-Gb length of a previous assembly of this genome, MesAur1.0. BCM_Maur_2.0 exhibits significantly improved continuity, with a scaffold N50 that is 6.7 times greater than MesAur1.0. Furthermore, 21,616 protein-coding genes and 10,459 noncoding genes are annotated in BCM_Maur_2.0 compared to 20,495 protein-coding genes and 4,168 noncoding genes in MesAur1.0. This new assembly also improves the unresolved regions as measured by nucleotide ambiguities, where ∼17.11% of bases in MesAur1.0 were unresolved compared to BCM_Maur_2.0, in which the number of unresolved bases is reduced to 3.00%. CONCLUSIONS: Access to a more complete reference genome with improved accuracy and continuity will facilitate more detailed, comprehensive, and meaningful research results for a wide variety of future studies using Syrian hamsters as models.

Overexpression of a Zea mays Brassinosteroid-Signaling Kinase Gene ZmBSK1 Confers Salt Stress Tolerance in Maize.

Salinity has become a crucial environmental factor seriously restricting maize (Zea mays L.) growth, development and productivity. However, how plants respond to salt stress is still poorly understood. In this study, we report that a maize brassinosteroid-signaling kinase gene ZmBSK1 plays a significant role in salt stress response. Expression pattern analysis revealed that the transcript level of ZmBSK1 was upregulated by NaCl treatment both in maize leaves, roots, and stems. Phenotypic and physiological analysis showed that overexpression of ZmBSK1 in maize improved salt tolerance by reducing the malondialdehyde (MDA) content, the percentage of electrolyte leakage, O2 - and H2O2 accumulation under salt stress, relying on the increases of antioxidant defense enzyme activities and proline content. qRT-PCR analysis showed that overexpression of ZmBSK1 also positively modulated the expression levels of reactive oxygen species (ROS)-scavenging and proline biosynthesis-related genes under salt stress. Moreover, immunoprecipitation-mass spectrometry (IP-MS) assay and firefly luciferase complementation imaging (LCI) assay showed that ZmBSK1 could associate with heat shock protein ZmHSP8 and 14-3-3-like protein ZmGF14-6, and their gene expression levels could be significantly induced by NaCl treatment in different maize tissues. Our findings unravel the new function of ZmBSK1 in salt stress response, which provides the theoretical bases for the improvement of maize salt resistance.

Fully resolved assembly of Cryptosporidium parvum.

BACKGROUND: Cryptosporidium parvum is an apicomplexan parasite commonly found across many host species with a global infection prevalence in human populations of 7.6%. Understanding its diversity and genomic makeup can help in fighting established infections and prohibiting further transmission. The basis of every genomic study is a high-quality reference genome that has continuity and completeness, thus enabling comprehensive comparative studies. FINDINGS: Here, we provide a highly accurate and complete reference genome of Cryptosporidium parvum. The assembly is based on Oxford Nanopore reads and was improved using Illumina reads for error correction. We also outline how to evaluate and choose from different assembly methods based on 2 main approaches that can be applied to other Cryptosporidium species. The assembly encompasses 8 chromosomes and includes 13 telomeres that were resolved. Overall, the assembly shows a high completion rate with 98.4% single-copy BUSCO genes. CONCLUSIONS: This high-quality reference genome of a zoonotic IIaA17G2R1 C. parvum subtype isolate provides the basis for subsequent comparative genomic studies across the Cryptosporidium clade. This will enable improved understanding of diversity, functional, and association studies.

Identification and Fine-Mapping of a Novel QTL, qMrdd2, That Confers Resistance to Maize Rough Dwarf Disease.

Maize rough dwarf disease (MRDD) is caused by a virus and seriously affects maize quality and yield worldwide. MRDD can be most effectively controlled with disease-resistant hybrids of corn. Here, MRDD-resistant (Qi319) and -susceptible (Ye478) parental inbred maize lines and their 314 recombinant inbred lines (RILs) that were derived from a cross between them were evaluated across three environments. A stable resistance QTL, qMrdd2, was identified and mapped using best linear unbiased prediction (BLUP) values to a 0.55-Mb region between the markers MK807 and MK811 on chromosome 2 (B73 RefGen_v3) and was found to explain 8.6 to 11.0% of the total phenotypic variance in MRDD resistance. We validated the effect of qMrdd2 using a chromosome segment substitution line (CSSL) that was derived from a cross between maize inbred Qi319 as the MRDD resistance donor and Ye478 as the recipient. Disease severity index of the CSSL haplotype II harboring qMrdd2 was significantly lower than that of the susceptible parent Ye478. Subsequently, we fine-mapped qMrdd2 to a 315-kb region flanked by the markers RD81 and RD87, thus testing recombinant-derived progeny using selfed backcrossed families. In this study, we identified a novel QTL for MRDD resistance by combining the RIL and CSSL populations, thus providing important genetic information that can be used for breeding MRDD-resistant varieties of maize.

AHDC1 missense mutations in Xia-Gibbs syndrome.

Xia-Gibbs syndrome (XGS; MIM: 615829) is a phenotypically heterogeneous neurodevelopmental disorder (NDD) caused by newly arising mutations in the AT-Hook DNA-Binding Motif-Containing 1 (AHDC1) gene that are predicted to lead to truncated AHDC1 protein synthesis. More than 270 individuals have been diagnosed with XGS worldwide. Despite the absence of an independent assay for AHDC1 protein function to corroborate potential functional consequences of rare variant genetic findings, there are also reports of individuals with XGS-like trait manifestations who have de novo missense AHDC1 mutations and who have been provided a molecular diagnosis of the disorder. To investigate a potential contribution of missense mutations to XGS, we mapped the missense mutations from 10 such individuals to the AHDC1 conserved protein domain structure and detailed the observed phenotypes. Five newly identified individuals were ascertained from a local XGS Registry, and an additional five were taken from external reports or databases, including one publication. Where clinical data were available, individuals with missense mutations all displayed phenotypes consistent with those observed in individuals with AHDC1 truncating mutations, including delayed motor milestones, intellectual disability (ID), hypotonia, and speech delay. A subset of the 10 reported missense mutations cluster in two regions of the AHDC1 protein with known conserved domains, likely representing functional motifs. Variants outside the clustered regions score lower for computational prediction of their likely damaging effects. Overall, de novo missense variants in AHDC1 are likely diagnostic of XGS when in silico analysis of their position relative to conserved regions is considered together with disease trait manifestations.

Transmission event of SARS-CoV-2 delta variant reveals multiple vaccine breakthrough infections.

BACKGROUND: This study aims to identify the causative strain of SARS-CoV-2 in a cluster of vaccine breakthroughs. Vaccine breakthrough by a highly transmissible SARS-CoV-2 strain is a risk to global public health. METHODS: Nasopharyngeal swabs from suspected vaccine breakthrough cases were tested for SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2) by qPCR (quantitative polymerase chain reaction) for Wuhan-Hu1 and alpha variant. Positive samples were then sequenced by Swift Normalase Amplicon Panels to determine the causal variant. GATK (genome analysis toolkit) variants were filtered with allele fraction ≥80 and min read depth 30x. RESULTS: Viral sequencing revealed an infection cluster of 6 vaccinated patients infected with the delta (B.1.617.2) SARS-CoV-2 variant. With no history of vaccine breakthrough, this suggests the delta variant may possess immune evasion in patients that received the Pfizer BNT162b2, Moderna mRNA-1273, and Covaxin BBV152. CONCLUSIONS: Delta variant may pose the highest risk out of any currently circulating SARS-CoV-2 variants, with previously described increased transmissibility over alpha variant and now, possible vaccine breakthrough. FUNDING: Parts of this work was supported by the National Institute of Allergy and Infectious Diseases (1U19AI144297) and Baylor College of Medicine internal funding.

Genetic testing in ambulatory cardiology clinics reveals high rate of findings with clinical management implications.

PURPOSE: Cardiovascular disease (CVD) is the leading cause of death in adults in the United States, yet the benefits of genetic testing are not universally accepted. METHODS: We developed the "HeartCare" panel of genes associated with CVD, evaluating high-penetrance Mendelian conditions, coronary artery disease (CAD) polygenic risk, LPA gene polymorphisms, and specific pharmacogenetic (PGx) variants. We enrolled 709 individuals from cardiology clinics at Baylor College of Medicine, and samples were analyzed in a CAP/CLIA-certified laboratory. Results were returned to the ordering physician and uploaded to the electronic medical record. RESULTS: Notably, 32% of patients had a genetic finding with clinical management implications, even after excluding PGx results, including 9% who were molecularly diagnosed with a Mendelian condition. Among surveyed physicians, 84% reported medical management changes based on these results, including specialist referrals, cardiac tests, and medication changes. LPA polymorphisms and high polygenic risk of CAD were found in 20% and 9% of patients, respectively, leading to diet, lifestyle, and other changes. Warfarin and simvastatin pharmacogenetic variants were present in roughly half of the cohort. CONCLUSION: Our results support the use of genetic information in routine cardiovascular health management and provide a roadmap for accompanying research.

Exome sequencing in children with clinically suspected maturity-onset diabetes of the young.

OBJECTIVE: Commercial gene panels identify pathogenic variants in as low as 27% of patients suspected to have MODY, suggesting the role of yet unidentified pathogenic variants. We sought to identify novel gene variants associated with MODY. RESEARCH DESIGN AND METHODS: We recruited 10 children with a clinical suspicion of MODY but non-diagnostic commercial MODY gene panels. We performed exome sequencing (ES) in them and their parents. RESULTS: Mean age at diabetes diagnosis was 10 (± 3.8) years. Six were females; 4 were non-Hispanic white, 5 Hispanic, and 1 Asian. Our variant prioritization analysis identified a pathogenic, de novo variant in INS (c.94G > A, p.Gly32Ser), confirmed by Sanger sequencing, in a proband who was previously diagnosed with "autoantibody-negative type 1 diabetes (T1D)" at 3 y/o. This rare variant, absent in the general population (gnomAD database), has been reported previously in neonatal diabetes. We also identified a frameshift deletion (c.2650delC, p.Gln884AsnfsTer57) in RFX6 in a child with a previous diagnosis of "autoantibody-negative T1D" at 12 y/o. The variant was inherited from the mother, who was diagnosed with "thin type 2 diabetes" at 25 y/o. Heterozygous protein-truncating variants in RFX6 gene have been recently reported in individuals with MODY. CONCLUSIONS: We diagnosed two patients with MODY using ES in children initially classified as "T1D". One has a likely pathogenic novel gene variant not previously associated with MODY. We demonstrate the clinical utility of ES in patients with clinical suspicion of MODY.

Oligonucleotide capture sequencing of the SARS-CoV-2 genome and subgenomic fragments from COVID-19 individuals.

The newly emerged and rapidly spreading SARS-CoV-2 causes coronavirus disease 2019 (COVID-19). To facilitate a deeper understanding of the viral biology we developed a capture sequencing methodology to generate SARS-CoV-2 genomic and transcriptome sequences from infected patients. We utilized an oligonucleotide probe-set representing the full-length genome to obtain both genomic and transcriptome (subgenomic open reading frames [ORFs]) sequences from 45 SARS-CoV-2 clinical samples with varying viral titers. For samples with higher viral loads (cycle threshold value under 33, based on the CDC qPCR assay) complete genomes were generated. Analysis of junction reads revealed regions of differential transcriptional activity among samples. Mixed allelic frequencies along the 20kb ORF1ab gene in one sample, suggested the presence of a defective viral RNA species subpopulation maintained in mixture with functional RNA in one sample. The associated workflow is straightforward, and hybridization-based capture offers an effective and scalable approach for sequencing SARS-CoV-2 from patient samples.

Transmission event of SARS-CoV-2 Delta variant reveals multiple vaccine breakthrough infections.

Importance: Vaccine breakthrough by an emergent SARS-CoV-2 variant poses a great risk to global public health. Objective: To determine the SARS-CoV-2 variant responsible for 6 cases of vaccine breakthrough. Design: Nasopharyngeal swabs from suspected vaccine breakthrough cases were tested for SARS-CoV-2 by qPCR for Wuhan-Hu1 and Alpha variant. Positive samples were then sequenced by Swift Normalase Amplicon Panels to determine the causal variant. Setting: Transmission event occurred at events surrounding a wedding outside of Houston, TX. Two patients from India, likely transmitted the Delta variant to other guests. Participants: Following a positive SARS-CoV-2 qPCR test at a third-party site, six fully vaccinated patients were investigated. Three males and three females ranged from 53 to 69 years old. One patient suffered from diabetes while three others were classified as overweight. No significant other comorbidities were identified. None of the patients had a history of failed vaccination.

Neptune: an environment for the delivery of genomic medicine.

PURPOSE: Genomic medicine holds great promise for improving health care, but integrating searchable and actionable genetic data into electronic health records (EHRs) remains a challenge. Here we describe Neptune, a system for managing the interaction between a clinical laboratory and an EHR system during the clinical reporting process. METHODS: We developed Neptune and applied it to two clinical sequencing projects that required report customization, variant reanalysis, and EHR integration. RESULTS: Neptune has been applied for the generation and delivery of over 15,000 clinical genomic reports. This work spans two clinical tests based on targeted gene panels that contain 68 and 153 genes respectively. These projects demanded customizable clinical reports that contained a variety of genetic data types including single-nucleotide variants (SNVs), copy-number variants (CNVs), pharmacogenomics, and polygenic risk scores. Two variant reanalysis activities were also supported, highlighting this important workflow. CONCLUSION: Methods are needed for delivering structured genetic data to EHRs. This need extends beyond developing data formats to providing infrastructure that manages the reporting process itself. Neptune was successfully applied on two high-throughput clinical sequencing projects to build and deliver clinical reports to EHR systems. The software is open source and available at https://gitlab.com/bcm-hgsc/neptune .

Sequencing of a central nervous system tumor demonstrates cancer transmission in an organ transplant.

Four organ transplant recipients from an organ donor diagnosed with anaplastic pleomorphic xanthoastrocytoma developed fatal malignancies for which the origin could not be confirmed by standard methods. We identified the somatic mutational profiles of the neoplasms using next-generation sequencing technologies and tracked the relationship between the different samples. The data were consistent with the presence of an aggressive clonal entity in the donor and the subsequent proliferation of descendent tumors in each recipient. Deleterious mutations in BRAF, PIK3CA, SDHC, DDR2, and FANCD2, and a chromosomal deletion spanning the CDKN2A/B genes, were shared between the recipients' lesions. In addition to demonstrating that DNA sequencing tracked a donor/recipient cancer transmission, this study established that the genetic profile of a donor tumor and its potential aggressive phenotype could have been determined before transplantation was considered. As the genetic correlates of tumor invasion and metastases become better known, adding genetic profiling by DNA sequencing to the data considered for transplant safety should be considered.