Associations of differentially expressed genes with psychoneurological symptoms in patients with head and neck cancer: A longitudinal study.

OBJECTIVE: Patients with head and neck cancer (HNC) experience psychoneurological symptoms (PNS, i.e., depression, fatigue, sleep disturbance, pain, and cognitive dysfunction) during intensity-modulated radiotherapy (IMRT) that negatively impact their functional status, quality of life, and overall survival. The underlying mechanisms for PNS are still not fully understood. This study aimed to examine differentially expressed genes and pathways related to PNS for patients undergoing IMRT (i.e., before, end of, 6 months, and 12 months after IMRT). METHODS: Participants included 142 patients with HNC (mean age 58.9 ± 10.3 years, 72.5% male, 83.1% White). Total RNA extracted from blood leukocytes were used for genome-wide gene expression assays. Linear mixed effects model was used to examine the association between PNS and gene expression across time. Gene Ontology (GO) enrichment analysis was employed to identify pathways related to PNS. RESULTS: A total of 1352 genes (162 upregulated, 1190 downregulated) were significantly associated with PNS across time (false discovery rate (FDR) < 0.05). Among these genes, 112 GO terms were identified (FDR < 0.05). The top 20 GO terms among the significant upregulated genes were related to immune and inflammatory responses, while the top 20 GO terms among the significant downregulated genes were associated with telomere maintenance. CONCLUSION: This study is the first to identify genes and pathways linked to immune and inflammatory responses and telomere maintenance that are associated with PNS in patients with HNC receiving IMRT. Inflammation and aging markers may be candidate biomarkers for PNS. Understanding biological markers may produce targets for novel interventions.

Live-cell invasive phenotyping uncovers the ALK2/BMP6 iron homeostasis pathway as a therapeutic vulnerability in LKB1-mutant lung cancer.

The acquisition of invasive properties is a prerequisite for tumor progression and metastasis. Molecular subtypes of KRAS-driven lung cancer exhibit distinct modes of invasion that likely contribute to unique growth properties and therapeutic susceptibilities. Despite this, pre-clinical discovery strategies designed to exploit invasive phenotypes are lacking. To address this, we designed an experimental system to screen for targetable signaling pathways linked to active early invasion phenotypes in the two most prominent molecular subtypes, TP53 and LKB1, of KRAS-driven lung adenocarcinoma (LUAD). By combining live-cell imaging of human bronchial epithelial cells in a 3D invasion matrix with RNA transcriptome profiling, we identified the LKB1-specific upregulation of bone morphogenetic protein 6 (BMP6). Examination of early-stage lung cancer patients confirmed upregulation of BMP6 in LKB1-mutant lung tumors. At the molecular level, we find that the canonical iron regulatory hormone Hepcidin is induced via BMP6 signaling upon LKB1 loss, where intact LKB1 kinase activity is necessary to maintain signaling homeostasis. Furthermore, pre-clinical studies in a novel Kras/Lkb1-mutant syngeneic mouse model show that potent growth suppression was achieved by inhibiting the ALK2/BMP6 signaling axis with single agents that are currently in clinical trials. We show that alterations in the iron homeostasis pathway are accompanied by simultaneous upregulation of ferroptosis protection proteins. Thus, LKB1 is sufficient to regulate both the 'gas' and 'breaks' to finely tune iron-regulated tumor progression.

Chromatin regulators in the TBX1 network confer risk for conotruncal heart defects in 22q11.2DS.

Congenital heart disease (CHD) affecting the conotruncal region of the heart, occurs in 40-50% of patients with 22q11.2 deletion syndrome (22q11.2DS). This syndrome is a rare disorder with relative genetic homogeneity that can facilitate identification of genetic modifiers. Haploinsufficiency of TBX1, encoding a T-box transcription factor, is one of the main genes responsible for the etiology of the syndrome. We suggest that genetic modifiers of conotruncal defects in patients with 22q11.2DS may be in the TBX1 gene network. To identify genetic modifiers, we analyzed rare, predicted damaging variants in whole genome sequence of 456 cases with conotruncal defects and 537 controls, with 22q11.2DS. We then performed gene set approaches and identified chromatin regulatory genes as modifiers. Chromatin genes with recurrent damaging variants include EP400, KAT6A, KMT2C, KMT2D, NSD1, CHD7 and PHF21A. In total, we identified 37 chromatin regulatory genes, that may increase risk for conotruncal heart defects in 8.5% of 22q11.2DS cases. Many of these genes were identified as risk factors for sporadic CHD in the general population. These genes are co-expressed in cardiac progenitor cells with TBX1, suggesting that they may be in the same genetic network. The genes KAT6A, KMT2C, CHD7 and EZH2, have been previously shown to genetically interact with TBX1 in mouse models. Our findings indicate that disturbance of chromatin regulatory genes impact the TBX1 gene network serving as genetic modifiers of 22q11.2DS and sporadic CHD, suggesting that there are some shared mechanisms involving the TBX1 gene network in the etiology of CHD.

HIV, asymptomatic STI, and the rectal mucosal immune environment among young men who have sex with men.

Young men who have sex with men (YMSM) are disproportionately affected by HIV and bacterial sexually transmitted infections (STI) including gonorrhea, chlamydia, and syphilis; yet research into the immunologic effects of these infections is typically pursued in siloes. Here, we employed a syndemic approach to understand potential interactions of these infections on the rectal mucosal immune environment among YMSM. We enrolled YMSM aged 18-29 years with and without HIV and/or asymptomatic bacterial STI and collected blood, rectal secretions, and rectal tissue biopsies. YMSM with HIV were on suppressive antiretroviral therapy (ART) with preserved blood CD4 cell counts. We defined 7 innate and 19 adaptive immune cell subsets by flow cytometry, the rectal mucosal transcriptome by RNAseq, and the rectal mucosal microbiome by 16S rRNA sequencing and examined the effects of HIV and STI and their interactions. We measured tissue HIV RNA viral loads among YMSM with HIV and HIV replication in rectal explant challenge experiments among YMSM without HIV. HIV, but not asymptomatic STI, was associated with profound alterations in the cellular composition of the rectal mucosa. We did not detect a difference in the microbiome composition associated with HIV, but asymptomatic bacterial STI was associated with a higher probability of presence of potentially pathogenic taxa. When examining the rectal mucosal transcriptome, there was evidence of statistical interaction; asymptomatic bacterial STI was associated with upregulation of numerous inflammatory genes and enrichment for immune response pathways among YMSM with HIV, but not YMSM without HIV. Asymptomatic bacterial STI was not associated with differences in tissue HIV RNA viral loads or in HIV replication in explant challenge experiments. Our results suggest that asymptomatic bacterial STI may contribute to inflammation particularly among YMSM with HIV, and that future research should examine potential harms and interventions to reduce the health impact of these syndemic infections.

Loss of the endocytic tumor suppressor HD-PTP phenocopies LKB1 and promotes RAS-driven oncogenesis.

Oncogenic RAS mutations drive aggressive cancers that are difficult to treat in the clinic, and while direct inhibition of the most common KRAS variant in lung adenocarcinoma (G12C) is undergoing clinical evaluation, a wide spectrum of oncogenic RAS variants together make up a large percentage of untargetable lung and GI cancers. Here we report that loss-of-function alterations (mutations and deep deletions) in the gene that encodes HD-PTP (PTPN23) occur in up to 14% of lung cancers in the ORIEN Avatar lung cancer cohort, associate with adenosquamous histology, and occur alongside an altered spectrum of KRAS alleles. Furthermore, we show that in publicly available early-stage NSCLC studies loss of HD-PTP is mutually exclusive with loss of LKB1, which suggests they restrict a common oncogenic pathway in early lung tumorigenesis. In support of this, knockdown of HD-PTP in RAS-transformed lung cancer cells is sufficient to promote FAK-dependent invasion. Lastly, knockdown of the Drosophila homolog of HD-PTP (dHD-PTP/Myopic) synergizes to promote RAS-dependent neoplastic progression. Our findings highlight a novel tumor suppressor that can restrict RAS-driven lung cancer oncogenesis and identify a targetable pathway for personalized therapeutic approaches for adenosquamous lung cancer.

Outbreak of severe acute respiratory coronavirus virus 2 (SARS-CoV-2) in hospitalized hemodialysis patients: An epidemiologic and genomic investigation.

We performed an epidemiological investigation and genome sequencing of severe acute respiratory coronavirus virus 2 (SARS-CoV-2) to define the source and scope of an outbreak in a cluster of hospitalized patients. Lack of appropriate respiratory hygiene led to SARS-CoV-2 transmission to patients and healthcare workers during a single hemodialysis session, highlighting the importance of infection prevention precautions.

Novel Missense CNTNAP2 Variant Identified in Two Consanguineous Pakistani Families With Developmental Delay, Epilepsy, Intellectual Disability, and Aggressive Behavior.

We report the genetic analysis of two consanguineous pedigrees of Pakistani ancestry in which two siblings in each family exhibited developmental delay, epilepsy, intellectual disability and aggressive behavior. Whole-genome sequencing was performed in Family 1, and we identified ~80,000 variants located in regions of homozygosity. Of these, 615 variants had a minor allele frequency ≤ 0.001, and 21 variants had CADD scores ≥ 15. Four homozygous exonic variants were identified in both affected siblings: PDZD7 (c.1348_1350delGAG, p.Glu450del), ALG6 (c.1033G>C, p.Glu345Gln), RBM20 (c.1587C>G, p.Ser529Arg), and CNTNAP2 (c.785G>A, p.Gly228Arg). Sanger sequencing revealed co-segregation of the PDZD7, RBM20, and CNTNAP2 variants with disease in Family 1. Pathogenic variants in PDZD7 and RBM20 are associated with autosomal recessive non-syndromic hearing loss and autosomal dominant dilated cardiomyopathy, respectively, suggesting that these variants are unlikely likely to contribute to the clinical presentation. Gene panel analysis was performed on the two affected siblings in Family 2, and they were found to also be homozygous for the p.Gly228Arg CNTNAP2 variant. Together these families provide a LOD score 2.9 toward p.Gly228Arg CNTNAP2 being a completely penetrant recessive cause of this disease. The clinical presentation of the affected siblings in both families is also consistent with previous reports from individuals with homozygous CNTNAP2 variants where at least one allele was a nonsense variant, frameshift or small deletion. Our data suggests that homozygous CNTNAP2 missense variants can also contribute to disease, thereby expanding the genetic landscape of CNTNAP2 dysfunction.

Staff Scientist Perspectives on Onboarding and Professional Development: A Case Study.

Across the United States, the number of staff scientists (master's- or doctoral-level professionals working in nonfaculty roles) has grown by 35% since 2010, and they play an increasingly important role in research efforts. However, few targeted resources are available, which potentially limits the effectiveness of this group. Launched in 2016, the staff scientist path at Emory has tripled in size over 4 y to 138 staff. The present case study evaluated the perceptions of staff scientists related to onboarding experiences and professional development needs, including those needs arising from coronavirus disease 2019 (COVID-19) impacts in the workplace. A survey of Emory staff scientists was conducted from May to June 2019 as part of a program evaluation initiative to assess perceptions of onboarding and professional development opportunities. Interviews with a subset of scientists informed the survey development and identified COVID-19-related impacts on daily work. Results indicated the need for targeted orientation resources specific to staff scientists, accurate and timely information and resources to support scientists' supervisors, and professional development for scientists in leadership and management-related skills. Remote work associated with COVID-19 accentuated the need for managerial skills, including team development in digital work environments. Findings from this case study can inform policies and practices at Emory and other institutions that employ a similar staff scientist model.

Sex-specific recombination patterns predict parent of origin for recurrent genomic disorders.

BACKGROUND: Structural rearrangements of the genome, which generally occur during meiosis and result in large-scale (> 1 kb) copy number variants (CNV; deletions or duplications ≥ 1 kb), underlie genomic disorders. Recurrent pathogenic CNVs harbor similar breakpoints in multiple unrelated individuals and are primarily formed via non-allelic homologous recombination (NAHR). Several pathogenic NAHR-mediated recurrent CNV loci demonstrate biases for parental origin of de novo CNVs. However, the mechanism underlying these biases is not well understood. METHODS: We performed a systematic, comprehensive literature search to curate parent of origin data for multiple pathogenic CNV loci. Using a regression framework, we assessed the relationship between parental CNV origin and the male to female recombination rate ratio. RESULTS: We demonstrate significant association between sex-specific differences in meiotic recombination and parental origin biases at these loci (p = 1.07 × 10-14). CONCLUSIONS: Our results suggest that parental origin of CNVs is largely influenced by sex-specific recombination rates and highlight the need to consider these differences when investigating mechanisms that cause structural variation.

Identifying susceptibility genes for primary ovarian insufficiency on the high-risk genetic background of a fragile X premutation.

OBJECTIVE: To identify modifying genes that explains the risk of fragile X-associated primary ovarian insufficiency (FXPOI). DESIGN: Gene-based, case/control association study, followed by a functional screen of highly ranked genes using a Drosophila model. SETTING: Participants were recruited from academic and clinical settings. PATIENT(S): Women with a premutation (PM) who experienced FXPOI at the age of 35 years or younger (n = 63) and women with a PM who experienced menopause at the age of 50 years or older (n = 51) provided clinical information and a deoxyribonucleic acid sample for whole genome sequencing. The functional screen was on the basis of Drosophila TRiP lines. INTERVENTION(S): Clinical information and a DNA sample were collected for whole genome sequencing. MAIN OUTCOME MEASURES: A polygenic risk score derived from common variants associated with natural age at menopause was calculated and associated with the risk of FXPOI. Genes associated with the risk of FXPOI were identified on the basis of the P-value from gene-based association test and an altered level of fecundity when knocked down in the Drosophila PM model. RESULTS: The polygenic risk score on the basis of common variants associated with natural age at menopause explained approximately 8% of the variance in the risk of FXPOI. Further, SUMO1 and KRR1 were identified as possible modifying genes associated with the risk of FXPOI on the basis of an untargeted gene analysis of rare variants. CONCLUSIONS: In addition to the large genetic effect of a PM on ovarian function, the additive effects of common variants associated with natural age at menopause and the effect of rare modifying variants appear to play a role in FXPOI risk.

Direct Immunodetection of Global A-to-I RNA Editing Activity with a Chemiluminescent Bioassay.

Adenosine-to-inosine (A-to-I) editing is a conserved eukaryotic RNA modification that contributes to development, immune response, and overall cellular function. Here, we utilize Endonuclease V (EndoV), which binds specifically to inosine in RNA, to develop an EndoV-linked immunosorbency assay (EndoVLISA) as a rapid, plate-based chemiluminescent method for measuring global A-to-I editing signatures in cellular RNA. We first optimize and validate our assay with chemically synthesized oligonucleotides. We then demonstrate rapid detection of inosine content in treated cell lines, demonstrating equivalent performance against current standard RNA-seq approaches. Lastly, we deploy our EndoVLISA for profiling differential A-to-I RNA editing signatures in normal and diseased human tissue, illustrating the utility of our platform as a diagnostic bioassay. Together, the EndoVLISA method is cost-effective, straightforward, and utilizes common laboratory equipment, offering a highly accessible new approach for studying A-to-I editing. Moreover, the multi-well plate format makes this the first assay amenable for direct high-throughput quantification of A-to-I editing for applications in disease detection and drug development.

Genetic contributors to risk of schizophrenia in the presence of a 22q11.2 deletion.

Schizophrenia occurs in about one in four individuals with 22q11.2 deletion syndrome (22q11.2DS). The aim of this International Brain and Behavior 22q11.2DS Consortium (IBBC) study was to identify genetic factors that contribute to schizophrenia, in addition to the ~20-fold increased risk conveyed by the 22q11.2 deletion. Using whole-genome sequencing data from 519 unrelated individuals with 22q11.2DS, we conducted genome-wide comparisons of common and rare variants between those with schizophrenia and those with no psychotic disorder at age ≥25 years. Available microarray data enabled direct comparison of polygenic risk for schizophrenia between 22q11.2DS and independent population samples with no 22q11.2 deletion, with and without schizophrenia (total n = 35,182). Polygenic risk for schizophrenia within 22q11.2DS was significantly greater for those with schizophrenia (padj = 6.73 × 10-6). Novel reciprocal case-control comparisons between the 22q11.2DS and population-based cohorts showed that polygenic risk score was significantly greater in individuals with psychotic illness, regardless of the presence of the 22q11.2 deletion. Within the 22q11.2DS cohort, results of gene-set analyses showed some support for rare variants affecting synaptic genes. No common or rare variants within the 22q11.2 deletion region were significantly associated with schizophrenia. These findings suggest that in addition to the deletion conferring a greatly increased risk to schizophrenia, the risk is higher when the 22q11.2 deletion and common polygenic risk factors that contribute to schizophrenia in the general population are both present.

Identifying genetic factors that contribute to the increased risk of congenital heart defects in infants with Down syndrome.

Atrioventricular septal defects (AVSD) are a severe congenital heart defect present in individuals with Down syndrome (DS) at a > 2000-fold increased prevalence compared to the general population. This study aimed to identify risk-associated genes and pathways and to examine a potential polygenic contribution to AVSD in DS. We analyzed a total cohort of 702 individuals with DS with or without AVSD, with genomic data from whole exome sequencing, whole genome sequencing, and/or array-based imputation. We utilized sequence kernel association testing and polygenic risk score (PRS) methods to examine rare and common variants. Our findings suggest that the Notch pathway, particularly NOTCH4, as well as genes involved in the ciliome including CEP290 may play a role in AVSD in DS. These pathways have also been implicated in DS-associated AVSD in prior studies. A polygenic component for AVSD in DS has not been examined previously. Using weights based on the largest genome-wide association study of congenital heart defects available (2594 cases and 5159 controls; all general population samples), we found PRS to be associated with AVSD with odds ratios ranging from 1.2 to 1.3 per standard deviation increase in PRS and corresponding liability r2 values of approximately 1%, suggesting at least a small polygenic contribution to DS-associated AVSD. Future studies with larger sample sizes will improve identification and quantification of genetic contributions to AVSD in DS.

Two New Aspergillus flavus Reference Genomes Reveal a Large Insertion Potentially Contributing to Isolate Stress Tolerance and Aflatoxin Production.

Efforts in genome sequencing in the Aspergillus genus have led to the development of quality reference genomes for several important species including A. nidulans, A. fumigatus, and A. oryzae However, less progress has been made for A. flavus As part of the effort of the USDA-ARS Annual Aflatoxin Workshop Fungal Genome Project, the isolate NRRL3357 was sequenced and resulted in a scaffold-level genome released in 2005. Our goal has been biologically driven, focusing on two areas: isolate variation in aflatoxin production and drought stress exacerbating aflatoxin production by A. flavus Therefore, we developed two reference pseudomolecule genome assemblies derived from chromosome arms for two isolates: AF13, a MAT1-2, highly stress tolerant, and highly aflatoxigenic isolate; and NRRL3357, a MAT1-1, less stress tolerant, and moderate aflatoxin producer in comparison to AF13. Here, we report these two reference-grade assemblies for these isolates through a combination of PacBio long-read sequencing and optical mapping, and coupled them with comparative, functional, and phylogenetic analyses. This analysis resulted in the identification of 153 and 45 unique genes in AF13 and NRRL3357, respectively. We also confirmed the presence of a unique 310 Kb insertion in AF13 containing 60 genes. Analysis of this insertion revealed the presence of a bZIP transcription factor, named atfC, which may contribute to isolate pathogenicity and stress tolerance. Phylogenomic analyses comparing these and other available assemblies also suggest that the species complex of A. flavus is polyphyletic.

Selective Enrichment of A-to-I Edited Transcripts from Cellular RNA Using Endonuclease V.

Creating accurate maps of A-to-I RNA editing activity is vital to improving our understanding of the biological role of this process and harnessing it as a signal for disease diagnosis. Current RNA sequencing techniques are susceptible to random sampling limitations due to the complexity of the transcriptome and require large amounts of RNA material, specialized instrumentation, and high read counts to accurately interrogate A-to-I editing sites. To address these challenges, we show that Escherichia coli Endonuclease V (eEndoV), an inosine-cleaving enzyme, can be repurposed to bind and isolate A-to-I edited transcripts from cellular RNA. While Mg2+ enables eEndoV to catalyze RNA cleavage, we show that similar levels of Ca2+ instead promote binding of inosine without cleavage and thus enable high affinity capture of inosine in RNA. We leverage this capability to demonstrate EndoVIPER-seq (Endonuclease V inosine precipitation enrichment sequencing) as a facile and effective method to enrich A-to-I edited transcripts prior to RNA-seq, producing significant increases in the coverage and detection of identified editing sites. We envision the use of this approach as a straightforward and cost-effective strategy to improve the epitranscriptomic informational density of RNA samples, facilitating a deeper understanding of the functional roles of A-to-I editing.

Complete Sequence of the 22q11.2 Allele in 1,053 Subjects with 22q11.2 Deletion Syndrome Reveals Modifiers of Conotruncal Heart Defects.

The 22q11.2 deletion syndrome (22q11.2DS) results from non-allelic homologous recombination between low-copy repeats termed LCR22. About 60%-70% of individuals with the typical 3 megabase (Mb) deletion from LCR22A-D have congenital heart disease, mostly of the conotruncal type (CTD), whereas others have normal cardiac anatomy. In this study, we tested whether variants in the hemizygous LCR22A-D region are associated with risk for CTDs on the basis of the sequence of the 22q11.2 region from 1,053 22q11.2DS individuals. We found a significant association (FDR p < 0.05) of the CTD subset with 62 common variants in a single linkage disequilibrium (LD) block in a 350 kb interval harboring CRKL. A total of 45 of the 62 variants were associated with increased risk for CTDs (odds ratio [OR) ranges: 1.64-4.75). Associations of four variants were replicated in a meta-analysis of three genome-wide association studies of CTDs in affected individuals without 22q11.2DS. One of the replicated variants, rs178252, is located in an open chromatin region and resides in the double-elite enhancer, GH22J020947, that is predicted to regulate CRKL (CRK-like proto-oncogene, cytoplasmic adaptor) expression. Approximately 23% of patients with nested LCR22C-D deletions have CTDs, and inactivation of Crkl in mice causes CTDs, thus implicating this gene as a modifier. Rs178252 and rs6004160 are expression quantitative trait loci (eQTLs) of CRKL. Furthermore, set-based tests identified an enhancer that is predicted to target CRKL and is significantly associated with CTD risk (GH22J020946, sequence kernal association test (SKAT) p = 7.21 × 10-5) in the 22q11.2DS cohort. These findings suggest that variance in CTD penetrance in the 22q11.2DS population can be explained in part by variants affecting CRKL expression.

Genome-Wide Association Study to Find Modifiers for Tetralogy of Fallot in the 22q11.2 Deletion Syndrome Identifies Variants in the GPR98 Locus on 5q14.3.

BACKGROUND: The 22q11.2 deletion syndrome (22q11.2DS; DiGeorge syndrome/velocardiofacial syndrome) occurs in 1 of 4000 live births, and 60% to 70% of affected individuals have congenital heart disease, ranging from mild to severe. In our cohort of 1472 subjects with 22q11.2DS, a total of 62% (n=906) have congenital heart disease and 36% (n=326) of these have tetralogy of Fallot (TOF), comprising the largest subset of severe congenital heart disease in the cohort. METHODS AND RESULTS: To identify common genetic variants associated with TOF in individuals with 22q11.2DS, we performed a genome-wide association study using Affymetrix 6.0 array and imputed genotype data. In our cohort, TOF was significantly associated with a genotyped single-nucleotide polymorphism (rs12519770, P=2.98×10-8) in an intron of the adhesion GPR98 (G-protein-coupled receptor V1) gene on chromosome 5q14.3. There was also suggestive evidence of association between TOF and several additional single-nucleotide polymorphisms in this region. Some genome-wide significant loci in introns or noncoding regions could affect regulation of genes nearby or at a distance. On the basis of this possibility, we examined existing Hi-C chromatin conformation data to identify genes that might be under shared transcriptional regulation within the region on 5q14.3. There are 6 genes in a topologically associated domain of chromatin with GPR98, including MEF2C (Myocyte-specific enhancer factor 2C). MEF2C is the only gene that is known to affect heart development in mammals and might be of interest with respect to 22q11.2DS. CONCLUSIONS: In conclusion, common variants may contribute to TOF in 22q11.2DS and may function in cardiac outflow tract development.

PEMapper and PECaller provide a simplified approach to whole-genome sequencing.

The analysis of human whole-genome sequencing data presents significant computational challenges. The sheer size of datasets places an enormous burden on computational, disk array, and network resources. Here, we present an integrated computational package, PEMapper/PECaller, that was designed specifically to minimize the burden on networks and disk arrays, create output files that are minimal in size, and run in a highly computationally efficient way, with the single goal of enabling whole-genome sequencing at scale. In addition to improved computational efficiency, we implement a statistical framework that allows for a base by base error model, allowing this package to perform as well or better than the widely used Genome Analysis Toolkit (GATK) in all key measures of performance on human whole-genome sequences.