Search

The human genome is packaged within a three-dimensional (3D) nucleus and organized into structural units known as compartments, topologically associating domains (TADs), and loops. TAD boundaries, separating adjacent TADs, have been found to be well conserved across mammalian species and more evolutionarily constrained than TADs themselves. Recent studies show that structural variants (SVs) can modify 3D genomes through the disruption of TADs, which play an essential role in insulating genes from outside regulatory elements' aberrant regulation. However, how SV affects the 3D genome structure and their association among different aspects of gene regulation and candidate cis -regulatory elements (cCREs) have rarely been studied systematically. Here, we assess the impact of SVs intersecting with TAD boundaries by developing an integrative Hi-C analysis pipeline, which enables the generation of an in-depth catalog of TADs and TAD boundaries in human lymphoblastoid cell lines (LCLs) to fill the gap of limited resources. Our catalog contains 18,865 TADs, including 4596 sub-TADs, with 185 SVs (TAD-SVs) that alter chromatin architecture. By leveraging the ENCODE registry of cCREs in humans, we determine that 34 of 185 TAD-SVs intersect with cCREs and observe significant enrichment of TAD-SVs within cCREs. This study provides a database of TADs and TAD-SVs in the human genome that will facilitate future investigations of the impact of SVs on chromatin structure and gene regulation in health and disease., (© 2024 Li et al.; Published by Cold Spring Harbor Laboratory Press.)

Fewer than half of individuals with a suspected Mendelian or monogenic condition receive a precise molecular diagnosis after comprehensive clinical genetic testing. Improvements in data quality and costs have heightened interest in using long-read sequencing (LRS) to streamline clinical genomic testing, but the absence of control data sets for variant filtering and prioritization has made tertiary analysis of LRS data challenging. To address this, the 1000 Genomes Project (1KGP) Oxford Nanopore Technologies Sequencing Consortium aims to generate LRS data from at least 800 of the 1KGP samples. Our goal is to use LRS to identify a broader spectrum of variation so we may improve our understanding of normal patterns of human variation. Here, we present data from analysis of the first 100 samples, representing all 5 superpopulations and 19 subpopulations. These samples, sequenced to an average depth of coverage of 37× and sequence read N50 of 54 kbp, have high concordance with previous studies for identifying single nucleotide and indel variants outside of homopolymer regions. Using multiple structural variant (SV) callers, we identify an average of 24,543 high-confidence SVs per genome, including shared and private SVs likely to disrupt gene function as well as pathogenic expansions within disease-associated repeats that were not detected using short reads. Evaluation of methylation signatures revealed expected patterns at known imprinted loci, samples with skewed X-inactivation patterns, and novel differentially methylated regions. All raw sequencing data, processed data, and summary statistics are publicly available, providing a valuable resource for the clinical genetics community to discover pathogenic SVs., (© 2024 Gustafson et al.; Published by Cold Spring Harbor Laboratory Press.)

Diverse sets of complete human genomes are required to construct a pangenome reference and to understand the extent of complex structural variation. Here, we sequence 65 diverse human genomes and build 130 haplotype-resolved assemblies (130 Mbp median continuity), closing 92% of all previous assembly gaps 1,2 and reaching telomere-to-telomere (T2T) status for 39% of the chromosomes. We highlight complete sequence continuity of complex loci, including the major histocompatibility complex (MHC), SMN1 / SMN2 , NBPF8 , and AMY1/AMY2 , and fully resolve 1,852 complex structural variants (SVs). In addition, we completely assemble and validate 1,246 human centromeres. We find up to 30-fold variation in α-satellite high-order repeat (HOR) array length and characterize the pattern of mobile element insertions into α-satellite HOR arrays. While most centromeres predict a single site of kinetochore attachment, epigenetic analysis suggests the presence of two hypomethylated regions for 7% of centromeres. Combining our data with the draft pangenome reference 1 significantly enhances genotyping accuracy from short-read data, enabling whole-genome inference 3 to a median quality value (QV) of 45. Using this approach, 26,115 SVs per sample are detected, substantially increasing the number of SVs now amenable to downstream disease association studies., Competing Interests: Competing Interests E.E.E. is a scientific advisory board member of Variant Bio, Inc. C. Lee is a scientific advisory board member of Nabsys and Genome Insight. S.K. has received travel funds to speak at events hosted by Oxford Nanopore Technologies. The following authors have previously disclosed a patent application (No. EP19169090) relevant to Strand-seq: J.O.K., T.M., and D.P. The other authors declare no competing interests.

The mechanisms underlying the selective regional vulnerability to neurodegeneration in Huntington's disease (HD) have not been fully defined. To explore the role of astrocytes in this phenomenon, we used single-nucleus and bulk RNAseq, lipidomics, HTT gene CAG repeat-length measurements, and multiplexed immunofluorescence on HD and control post-mortem brains. We identified genes that correlated with CAG repeat length, which were enriched in astrocyte genes, and lipidomic signatures that implicated poly-unsaturated fatty acids in sensitizing neurons to cell death. Because astrocytes play essential roles in lipid metabolism, we explored the heterogeneity of astrocytic states in both protoplasmic and fibrous-like (CD44+) astrocytes. Significantly, one protoplasmic astrocyte state showed high levels of metallothioneins and was correlated with the selective vulnerability of distinct striatal neuronal populations. When modeled in vitro, this state improved the viability of HD-patient-derived spiny projection neurons. Our findings uncover key roles of astrocytic states in protecting against neurodegeneration in HD., (© 2024. The Author(s).)

Autism spectrum disorder (ASD) is a complex neurodevelopmental condition with a strong genetic component in which rare variants contribute significantly to risk. We performed whole genome and/or exome sequencing (WGS and WES) and SNP-array analysis to identify both rare sequence and copy number variants (SNVs and CNVs) in 435 individuals from 116 ASD families. We identified 37 rare potentially damaging de novo SNVs (pdSNVs) in the cases (n = 144). Interestingly, two of them (one stop-gain and one missense variant) occurred in the same gene, BRSK2. Moreover, the identification of 8 severe de novo pdSNVs in genes not previously implicated in ASD (AGPAT3, IRX5, MGAT5B, RAB8B, RAP1A, RASAL2, SLC9A1, YME1L1) highlighted promising candidates. Potentially damaging CNVs (pdCNVs) provided support to the involvement of inherited variants in PHF3, NEGR1, TIAM1 and HOMER1 in neurodevelopmental disorders (NDD), although mostly acting as susceptibility factors with incomplete penetrance. Interpretation of identified pdSNVs/pdCNVs according to the ACMG guidelines led to a molecular diagnosis in 19/144 cases, although this figure represents a lower limit and is expected to increase thanks to further clarification of the role of likely pathogenic variants in ASD/NDD candidate genes not yet established. In conclusion, our study highlights promising ASD candidate genes and contributes to characterize the allelic diversity, mode of inheritance and phenotypic impact of de novo and inherited risk variants in ASD/NDD genes., (© 2024. The Author(s).)

Primary open-angle glaucoma (POAG), the leading cause of irreversible blindness worldwide, disproportionately affects individuals of African ancestry. We conducted a genome-wide association study (GWAS) for POAG in 11,275 individuals of African ancestry (6,003 cases; 5,272 controls). We detected 46 risk loci associated with POAG at genome-wide significance. Replication and post-GWAS analyses, including functionally informed fine-mapping, multiple trait co-localization, and in silico validation, implicated two previously undescribed variants (rs1666698 mapping to DBF4P2; rs34957764 mapping to ROCK1P1) and one previously associated variant (rs11824032 mapping to ARHGEF12) as likely causal. For individuals of African ancestry, a polygenic risk score (PRS) for POAG from our mega-analysis (African ancestry individuals) outperformed a PRS from summary statistics of a much larger GWAS derived from European ancestry individuals. This study quantifies the genetic architecture similarities and differences between African and non-African ancestry populations for this blinding disease., Competing Interests: Declaration of interests J.M.O. is a member of the scientific advisory board of Life Biosciences and a paid consultant of Atheneum Partners, Cerner Enviza (includes Kantar Health), and Calico. A.G.R. holds intellectual property for the use of gene therapy to treat glaucoma. E.M.-E. is a scientific advisor for Avisi and a paid consultant of Aerie Pharmaceuticals, Allergan, Eyenovia, and Thea Pharma. J.L. receives instrument support from Carl Zeiss Meditech, Inc., and Heidelberg Engineering, GmBH; receives research support from Novartis, Inc.; and is a paid consultant at Thea, Inc., Alcon Laboratories, Inc., Johnson & Johnson, Inc., Abbvie, Inc., Carl Zeiss Meditech, Inc., Genetech, Inc., and ONL Therapeutics, Inc., (Copyright © 2023 Elsevier Inc. All rights reserved.)

New Jersey was among the first states impacted by the COVID-19 pandemic, with one of the highest overall death rates in the nation. Nevertheless, relatively few reports have been published focusing specifically on New Jersey. Here we report on molecular, clinical, and epidemiologic observations, from the largest healthcare network in the state, in a cohort of vaccinated and unvaccinated individuals with laboratory-confirmed SARS-CoV-2 infection. We conducted molecular surveillance of SARS-CoV-2-positive nasopharyngeal swabs collected in nine hospitals from December 2020 through June 2022, using both whole genome sequencing (WGS) and a real-time RT-PCR screening assay targeting spike protein mutations found in variants of concern (VOCs) within our region. De-identified clinical data were obtained retrospectively, including demographics, COVID-19 vaccination status, ICU admission, ventilator support, mortality, and medical history. Statistical analyses were performed to identify associations between SARS-CoV-2 variants, vaccination status, clinical outcomes, and medical risk factors. A total of 5007 SARS-CoV-2-positive nasopharyngeal swabs were successfully screened and/or sequenced. Variant screening identified three predominant VOCs, including Alpha ( n = 714), Delta ( n = 1877), and Omicron ( n = 1802). Omicron isolates were further sub-typed as BA.1 ( n = 899), BA.2 ( n = 853), or BA.4/BA.5 ( n = 50); the remaining 614 isolates were classified as "Other". Approximately 31.5% (1577/5007) of the samples were associated with vaccine breakthrough infections, which increased in frequency following the emergence of Delta and Omicron. Severe clinical outcomes included ICU admission (336/5007 = 6.7%), ventilator support (236/5007 = 4.7%), and mortality (430/5007 = 8.6%), with increasing age being the most significant contributor to each ( p < 0.001). Unvaccinated individuals accounted for 79.7% (268/336) of ICU admissions, 78.3% (185/236) of ventilator cases, and 74.4% (320/430) of deaths. Highly significant ( p < 0.001) increases in mortality were observed in individuals with cardiovascular disease, hypertension, cancer, diabetes, and hyperlipidemia, but not with obesity, thyroid disease, or respiratory disease. Significant differences ( p < 0.001) in clinical outcomes were also noted between SARS-CoV-2 variants, including Delta, Omicron BA.1, and Omicron BA.2. Vaccination was associated with significantly improved clinical outcomes in our study, despite an increase in breakthrough infections associated with waning immunity, greater antigenic variability, or both. Underlying comorbidities contributed significantly to mortality in both vaccinated and unvaccinated individuals, with increasing risk based on the total number of comorbidities. Real-time RT-PCR-based screening facilitated timely identification of predominant variants using a minimal number of spike protein mutations, with faster turnaround time and reduced cost compared to WGS. Continued evolution of SARS-CoV-2 variants will likely require ongoing surveillance for new VOCs, with real-time assessment of clinical impact.

Several studies have identified rare and common genetic variants associated with severe COVID-19, but no study has reported genetic determinants as predisposition factors for neurological complications. In this report, we identified rare/unique structural variants (SVs) implicated in neurological functions in two individuals with neurological manifestations of COVID-19. This report highlights the possible genetic link to the neurological symptoms with COVID-19 and calls for a collective effort to study these cohorts for a possible genetic linkage.

Examining the neutralizing capacity of monoclonal antibodies (MAbs) used to treat COVID-19, as well as antibodies recovered from unvaccinated, previously vaccinated, and infected individuals, against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern (VOCs) remains critical to study. Here, we report on a SARS-CoV-2 nosocomial outbreak caused by the SARS-CoV-2 R.1 variant harboring the E484K mutation in a 281-bed psychiatric facility in New Jersey among unvaccinated inpatients and health care professionals (HCPs). A total of 81 inpatients and HCPs tested positive for SARS-Cov-2 by reverse transcription (RT)-PCR from 29 October 9 to 30 November 2020. The R.1 variant exhibits partial or complete resistance to two MAbs in clinical use, as well as 2 receptor binding domain MAbs and 4 N-terminal domain (NTD) MAbs. NTD MAbs against pseudovirus harboring single characteristic R.1 mutations highlight the role of S255F in loss of activity. Additionally, we note dampened neutralization capacity by plasma from individuals with previous SARS-CoV-2 infection or sera from vaccinated individuals. The relative resistance of the R.1 variant is likely lower than that of B.1.351 and closer to that of P.1 and B.1.526. The R.1 lineage has been reported in 47 states in the United States and 40 countries. Although high proportions exhibited symptoms (26% and 61% among patients and HCPs, respectively) and relative antibody resistance, we detected only 10 R.1 variants from over 2,900 samples (~0.34%) collected from January to October 2021. Among 3 vaccinated individuals previously infected with R.1, we observed robust neutralizing antibody responses against SARS-CoV-2 wild type and VOCs. IMPORTANCE The neutralizing capacities of monoclonal antibodies used to treat COVID-19 and of those recovered from previously infected and vaccinated individuals against SARS-CoV-2 variants of concern (VOCs) remain important questions. We report on a nosocomial outbreak caused by a SARS-CoV-2 R.1 variant harboring an E484K mutation among 81 unvaccinated inpatients and health care professionals. We note high attack rates with symptoms in nearly 50% of infected individuals, in sharp contrast to an unrelated institutional outbreak caused by the R.1 variant among a vaccinated population. We found little evidence of significant community spillover. This variant exhibits partial or complete resistance to two monoclonal antibodies in clinical use and dampened the neutralization capacity of convalescent-phase plasma from individuals with previous infection or sera from vaccinated individuals. Among three vaccinated individuals previously infected with R.1, we observed robust neutralizing antibody responses against SARS-CoV-2 wild type and VOCs. These findings underscore the importance of vaccination for prevention of symptomatic COVID-19 disease.

The 1000 Genomes Project (1kGP) is the largest fully open resource of whole-genome sequencing (WGS) data consented for public distribution without access or use restrictions. The final, phase 3 release of the 1kGP included 2,504 unrelated samples from 26 populations and was based primarily on low-coverage WGS. Here, we present a high-coverage 3,202-sample WGS 1kGP resource, which now includes 602 complete trios, sequenced to a depth of 30X using Illumina. We performed single-nucleotide variant (SNV) and short insertion and deletion (INDEL) discovery and generated a comprehensive set of structural variants (SVs) by integrating multiple analytic methods through a machine learning model. We show gains in sensitivity and precision of variant calls compared to phase 3, especially among rare SNVs as well as INDELs and SVs spanning frequency spectrum. We also generated an improved reference imputation panel, making variants discovered here accessible for association studies., Competing Interests: Declaration of interests E.E.E. is a scientific advisory board (SAB) member of Variant Bio, Inc. P.F. is an SAB member of Fabric Genomics, Inc., and Eagle Genomics, Ltd., (Copyright © 2022 The Authors. Published by Elsevier Inc. All rights reserved.)

CFTR F508del (c.1521&#95;1523delCTT, p.Phe508delPhe) is the most common pathogenic allele underlying cystic fibrosis (CF), and its frequency varies in a geographic cline across Europe. We hypothesized that genetic variation associated with this cline is overrepresented in a large cohort (N > 5,000) of persons with CF who underwent whole-genome sequencing and that this pattern could result in spurious associations between variants correlated with both the F508del genotype and CF-related outcomes. Using principal-component (PC) analyses, we showed that variation in the CFTR region disproportionately contributes to a PC explaining a relatively high proportion of genetic variance. Variation near CFTR was correlated with population structure among persons with CF, and this correlation was driven by a subset of the sample inferred to have European ancestry. We performed genome-wide association studies comparing persons with CF with one versus two copies of the F508del allele; this allowed us to identify genetic variation associated with the F508del allele and to determine that standard PC-adjustment strategies eliminated the significant association signals. Our results suggest that PC adjustment can adequately prevent spurious associations between genetic variants and CF-related traits and are therefore effective tools to control for population structure even when population structure is confounded with disease severity and a common pathogenic variant., Competing Interests: M.B. is the editor-in-chief and J.X.C. (member of the Cystic Fibrosis Genome Project) is the deputy editor of HGG Advances. The authors declare no other competing interests., (© 2022 The Author(s).)

Studies of de novo mutation (DNM) have typically excluded some of the most repetitive and complex regions of the genome because these regions cannot be unambiguously mapped with short-read sequencing data. To better understand the genome-wide pattern of DNM, we generated long-read sequence data from an autism parent-child quad with an affected female where no pathogenic variant had been discovered in short-read Illumina sequence data. We deeply sequenced all four individuals by using three sequencing platforms (Illumina, Oxford Nanopore, and Pacific Biosciences) and three complementary technologies (Strand-seq, optical mapping, and 10X Genomics). Using long-read sequencing, we initially discovered and validated 171 DNMs across two children-a 20% increase in the number of de novo single-nucleotide variants (SNVs) and indels when compared to short-read callsets. The number of DNMs further increased by 5% when considering a more complete human reference (T2T-CHM13) because of the recovery of events in regions absent from GRCh38 (e.g., three DNMs in heterochromatic satellites). In total, we validated 195 de novo germline mutations and 23 potential post-zygotic mosaic mutations across both children; the overall true substitution rate based on this integrated callset is at least 1.41 × 10 -8 substitutions per nucleotide per generation. We also identified six de novo insertions and deletions in tandem repeats, two of which represent structural variants. We demonstrate that long-read sequencing and assembly, especially when combined with a more complete reference genome, increases the number of DNMs by >25% compared to previous studies, providing a more complete catalog of DNM compared to short-read data alone., Competing Interests: Declaration of interests E.E.E. is a scientific advisory board (SAB) member of Variant Bio, Inc., (Copyright © 2022 American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

Typical genotyping workflows map reads to a reference genome before identifying genetic variants. Generating such alignments introduces reference biases and comes with substantial computational burden. Furthermore, short-read lengths limit the ability to characterize repetitive genomic regions, which are particularly challenging for fast k-mer-based genotypers. In the present study, we propose a new algorithm, PanGenie, that leverages a haplotype-resolved pangenome reference together with k-mer counts from short-read sequencing data to genotype a wide spectrum of genetic variation-a process we refer to as genome inference. Compared with mapping-based approaches, PanGenie is more than 4 times faster at 30-fold coverage and achieves better genotype concordances for almost all variant types and coverages tested. Improvements are especially pronounced for large insertions (≥50 bp) and variants in repetitive regions, enabling the inclusion of these classes of variants in genome-wide association studies. PanGenie efficiently leverages the increasing amount of haplotype-resolved assemblies to unravel the functional impact of previously inaccessible variants while being faster compared with alignment-based workflows., (© 2022. The Author(s).)

To identify genetic determinants of airway dysfunction, we performed a transcriptome-wide association study for asthma by combining RNA-seq data from the nasal airway epithelium of 681 children, with UK Biobank genetic association data. Our airway analysis identified 95 asthma genes, 58 of which were not identified by transcriptome-wide association analyses using other asthma-relevant tissues. Among these genes were MUC5AC, an airway mucin, and FOXA3, a transcriptional driver of mucus metaplasia. Muco-ciliary epithelial cultures from genotyped donors revealed that the MUC5AC risk variant increases MUC5AC protein secretion and mucus secretory cell frequency. Airway transcriptome-wide association analyses for mucus production and chronic cough also identified MUC5AC. These cis-expression variants were associated with trans effects on expression; the MUC5AC variant was associated with upregulation of non-inflammatory mucus secretory network genes, while the FOXA3 variant was associated with upregulation of type-2 inflammation-induced mucus-metaplasia pathway genes. Our results reveal genetic mechanisms of airway mucus pathobiology., (© 2022. The Author(s).)

Impressive global efforts have identified both rare and common gene variants associated with severe COVID-19 using sequencing technologies. However, these studies lack the sensitivity to accurately detect several classes of variants, especially large structural variants (SVs), which account for a substantial proportion of genetic diversity including clinically relevant variation. We performed optical genome mapping on 52 severely ill COVID-19 patients to identify rare/unique SVs as decisive predisposition factors associated with COVID-19. We identified 7 SVs involving genes implicated in two key host-viral interaction pathways: innate immunity and inflammatory response, and viral replication and spread in nine patients, of which SVs in STK26 and DPP4 genes are the most intriguing candidates. This study is the first to systematically assess the potential role of SVs in the pathogenesis of COVID-19 severity and highlights the need to evaluate SVs along with sequencing variants to comprehensively associate genomic information with interindividual variability in COVID-19 phenotypes., Competing Interests: J.L., A.H., and A.C. are salaried employees of Bionano Genomics Inc. R.K. has received honoraria, travel funding, and research support from Illumina, Asuragen, QIAGEN, Perkin Elmer Inc, Bionano Genomics, and BMS. A.H.B. has received funding from the NIH, MDA (USA), AFM Telethon, Alexion Pharmaceuticals Inc., Audentes Therapeutics Inc., Dynacure SAS, and Pfizer Inc. He has consulted and received compensation or honoraria from Asklepios BioPharmaceutical Inc, Audentes Therapeutics, Biogen, F. Hoffman-La Roche AG, GLG Inc, Guidepoint Global, and Kate Therapeutics, and holds equity in Ballard Biologics and Kate Therapeutics., (© 2022 The Author(s).)

The molecular mechanisms underlying the clinical manifestations of coronavirus disease 2019 (COVID-19), and what distinguishes them from common seasonal influenza virus and other lung injury states such as acute respiratory distress syndrome, remain poorly understood. To address these challenges, we combine transcriptional profiling of 646 clinical nasopharyngeal swabs and 39 patient autopsy tissues to define body-wide transcriptome changes in response to COVID-19. We then match these data with spatial protein and expression profiling across 357 tissue sections from 16 representative patient lung samples and identify tissue-compartment-specific damage wrought by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, evident as a function of varying viral loads during the clinical course of infection and tissue-type-specific expression states. Overall, our findings reveal a systemic disruption of canonical cellular and transcriptional pathways across all tissues, which can inform subsequent studies to combat the mortality of COVID-19 and to better understand the molecular dynamics of lethal SARS-CoV-2 and other respiratory infections., Competing Interests: O.E. is scientific adviser and equity holder in Freenome, Owkin, Volastra Therapeutics, and OneThree Biotech. R.E.S. is on the scientific advisory board of Miromatrix, Inc., and is a consultant and speaker for Alnylam, Inc. L.S. is a scientific co-founder and paid consultant. C.M. and E.E.A. are consultants for Onegevity Health. C.E.M. is a co-founder of Biotia and Onegevity Health and an advisor to Nanostring. T.H., S.W., Y.K., and J.R. are employees of Nanostring, Inc. All other authors declare no competing interests., (© 2022 The Author(s).)

Emergence of SARS-CoV-2 with high transmission and immune evasion potential, the so-called variants of concern (VOC), is a major concern. We describe the early genomic epidemiology of SARS-CoV-2 recovered from vaccinated health care professionals (HCP). Our postvaccination COVID-19 symptoms-based surveillance program among HCPs in a 17-hospital network identified all vaccinated HCPs who tested positive for COVID-19 after routine screening or after self-reporting. From 1 January 2021 to 30 April 2021, 23,687 HCPs received either mRNA-1273 or BNT162b2 mRNA vaccine. All available postvaccination SARS-CoV-2 samples and a random collection from nonvaccinated patients during the similar time frame were subjected to VOC screening and whole-genome sequencing (WGS). Sixty-two percent (23,697/37,500) of HCPs received at least one vaccine dose, with 60% (22,458) fully vaccinated. We detected 138 (0.58%, 138/23,697) COVID-19 cases, 105 among partially vaccinated and 33 (0.15%, 33/22,458) among fully vaccinated. Five partially vaccinated required hospitalization, four with supplemental oxygen. VOC screening from 16 fully vaccinated HCPs identified 6 (38%) harboring N501Y and 1 (6%) with E484K polymorphisms; percentage of concurrent nonvaccinated samples was 37% (523/1,404) and 20% (284/1,394), respectively. There was an upward trend from January to April for E484K/Q (3% to 26%) and N501Y (1% to 49%). WGS analysis from vaccinated and nonvaccinated individuals indicated highly congruent phylogenies. We did not detect an increased frequency of any receptor-binding domain (RBD)/N-terminal domain (NTD) polymorphism between groups ( P  > 0.05). Our results support robust protection by vaccination, particularly among recipients of both doses. Despite VOCs accounting for over 40% of SARS-CoV-2 from fully vaccinated individuals, the genomic diversity appears to proportionally represent VOCs among nonvaccinated populations. IMPORTANCE A number of highly effective vaccines have been developed and deployed to combat the COVID-19 pandemic. The emergence and epidemiological dominance of SARS-CoV-2 mutants with high transmission potential and immune evasion properties, the so-called variants of concern (VOC), continue to be a major concern. Whether these VOCs alter the efficacy of the administered vaccines is of great concern and a critical question to study. We describe the initial genomic epidemiology of SARS-CoV-2 recovered from partial/fully vaccinated health care professionals and probe specifically for VOC enrichment. Our findings support the high level of protection provided by full vaccination despite a steep increase in the prevalence of polymorphisms associated with increased transmission potential (N501Y) and immune evasion (E484K) in the nonvaccinated population. Thus, we do not find evidence of VOC enrichment among vaccinated groups. Overall, the genomic diversity of SARS-CoV-2 recovered postvaccination appears to proportionally represent the observed viral diversity within the community.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern (VOCs), harboring spike protein N-terminal domain (NTD) or receptor-binding domain (RBD) mutations, exhibit reduced in vitro susceptibility to convalescent-phase serum, commercial antibody cocktails, and vaccine neutralization and have been associated with reinfections. The accumulation of these mutations could be the consequence of intrahost viral evolution due to prolonged infection in immunocompromised hosts. In this study, we document the microevolution of SARS-CoV-2 recovered from sequential tracheal aspirates from an immunosuppressed patient on steroids and convalescent plasma therapy and identify the emergence of multiple NTD and RBD mutations. SARS-CoV-2 genomes from the first swab (day 0) and from three tracheal aspirates (days 7, 21, and 27) were compared at the sequence level. We identified a mixed viral population with five different S protein mutations (141 to 144 deletion, 243 to 244 deletion, E484K, Q493K, and Q493R) at the NTD or RBD region from the second tracheal aspirate sample (day 21) and a predominance of the S protein 141 to 144 LGVY deletion and E484K mutant on day 27. The neutralizing antibodies against various S protein lentiviral pseudovirus mutants, as well as the anti-SARS-CoV-2 total Ig and IgG, showed "U" shape dynamics, in support of the endogenous development of neutralizing antibodies. The patient's compromised immune status, the antirejection regiment, convalescent plasma treatment, and the development of neutralizing antibodies may have resulted in unique selective pressures on the intrahost genomic evolution, and this observation supports the hypotheses that VOCs can independently arise and that immunocompromised patients on convalescent plasma therapy are potential breeding grounds for immune escape mutants. IMPORTANCE Over a year of the COVID-19 pandemic, distinct severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) lineages have arisen in multiple geographic areas around the world. SARS-CoV-2 variants of concern (VOCs), i.e., B.1.1.7 (alpha), B.1.351 (beta), P.1 (gamma), and B.1.617.2 (delta), harboring mutations and/or deletions in spike protein N-terminal domain (NTD) or receptor-binding domain (RBD) regions showed evidence of increased transmissibility and disease severity and possible reduced vaccine efficacy. In this study, we report the emergence of five different NTD and RBD mutations in an uncommon SARS-CoV-2 B.1.369 lineage from an immunosuppressed patient undergoing steroid and convalescent plasma therapy. The observation highlighted that VOCs can independently arise in immunocompromised populations undergoing anti-SARS-CoV-2 therapy, and enhanced measures will be required to reduce the transmission.

Autism is a highly heritable complex disorder in which de novo mutation (DNM) variation contributes significantly to risk. Using whole-genome sequencing data from 3,474 families, we investigate another source of large-effect risk variation, ultra-rare variants. We report and replicate a transmission disequilibrium of private, likely gene-disruptive (LGD) variants in probands but find that 95% of this burden resides outside of known DNM-enriched genes. This variant class more strongly affects multiplex family probands and supports a multi-hit model for autism. Candidate genes with private LGD variants preferentially transmitted to probands converge on the E3 ubiquitin-protein ligase complex, intracellular transport and Erb signaling protein networks. We estimate that these variants are approximately 2.5 generations old and significantly younger than other variants of similar type and frequency in siblings. Overall, private LGD variants are under strong purifying selection and appear to act on a distinct set of genes not yet associated with autism., (© 2021. The Author(s), under exclusive licence to Springer Nature America, Inc.)

Summary: We present a new version of the popular somatic variant caller, Lancet, that supports the analysis of linked-reads sequencing data. By seamlessly integrating barcodes and haplotype read assignments within the colored De Bruijn graph local-assembly framework, Lancet computes a barcode-aware coverage and identifies variants that disagree with the local haplotype structure., Availability and Implementation: Lancet is implemented in C++ and available for academic and non-commercial research purposes as an open-source package at https://github.com/nygenome/lancet., Supplementary Information: Supplementary data are available at Bioinformatics online., (© The Author(s) 2020. Published by Oxford University Press. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Background: Previous research in autism and other neurodevelopmental disorders (NDDs) has indicated an important contribution of protein-coding (coding) de novo variants (DNVs) within specific genes. The role of de novo noncoding variation has been observable as a general increase in genetic burden but has yet to be resolved to individual functional elements. In this study, we assessed whole-genome sequencing data in 2671 families with autism (discovery cohort of 516 families, replication cohort of 2155 families). We focused on DNVs in enhancers with characterized in vivo activity in the brain and identified an excess of DNVs in an enhancer named hs737., Results: We adapted the fitDNM statistical model to work in noncoding regions and tested enhancers for excess of DNVs in families with autism. We found only one enhancer (hs737) with nominal significance in the discovery (p = 0.0172), replication (p = 2.5 × 10 -3 ), and combined dataset (p = 1.1 × 10 -4 ). Each individual with a DNV in hs737 had shared phenotypes including being male, intact cognitive function, and hypotonia or motor delay. Our in vitro assessment of the DNVs showed they all reduce enhancer activity in a neuronal cell line. By epigenomic analyses, we found that hs737 is brain-specific and targets the transcription factor gene EBF3 in human fetal brain. EBF3 is genome-wide significant for coding DNVs in NDDs (missense p = 8.12 × 10 -35 , loss-of-function p = 2.26 × 10 -13 ) and is widely expressed in the body. Through characterization of promoters bound by EBF3 in neuronal cells, we saw enrichment for binding to NDD genes (p = 7.43 × 10 -6 , OR = 1.87) involved in gene regulation. Individuals with coding DNVs have greater phenotypic severity (hypotonia, ataxia, and delayed development syndrome [HADDS]) in comparison to individuals with noncoding DNVs that have autism and hypotonia., Conclusions: In this study, we identify DNVs in the hs737 enhancer in individuals with autism. Through multiple approaches, we find hs737 targets the gene EBF3 that is genome-wide significant in NDDs. By assessment of noncoding variation and the genes they affect, we are beginning to understand their impact on gene regulatory networks in NDDs., (© 2021. The Author(s).)

Long-read and strand-specific sequencing technologies together facilitate the de novo assembly of high-quality haplotype-resolved human genomes without parent-child trio data. We present 64 assembled haplotypes from 32 diverse human genomes. These highly contiguous haplotype assemblies (average minimum contig length needed to cover 50% of the genome: 26 million base pairs) integrate all forms of genetic variation, even across complex loci. We identified 107,590 structural variants (SVs), of which 68% were not discovered with short-read sequencing, and 278 SV hotspots (spanning megabases of gene-rich sequence). We characterized 130 of the most active mobile element source elements and found that 63% of all SVs arise through homology-mediated mechanisms. This resource enables reliable graph-based genotyping from short reads of up to 50,340 SVs, resulting in the identification of 1526 expression quantitative trait loci as well as SV candidates for adaptive selection within the human population., (Copyright © 2021 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

RTK/RAS/RAF pathway alterations (RPAs) are a hallmark of lung adenocarcinoma (LUAD). In this study, we use whole-genome sequencing (WGS) of 85 cases found to be RPA(-) by previous studies from The Cancer Genome Atlas (TCGA) to characterize the minority of LUADs lacking apparent alterations in this pathway. We show that WGS analysis uncovers RPA(+) in 28 (33%) of the 85 samples. Among the remaining 57 cases, we observe focal deletions targeting the promoter or transcription start site of STK11 (n = 7) or KEAP1 (n = 3), and promoter mutations associated with the increased expression of ILF2 (n = 6). We also identify complex structural variations associated with high-level copy number amplifications. Moreover, an enrichment of focal deletions is found in TP53 mutant cases. Our results indicate that RPA(-) cases demonstrate tumor suppressor deletions and genome instability, but lack unique or recurrent genetic lesions compensating for the lack of RPAs. Larger WGS studies of RPA(-) cases are required to understand this important LUAD subset., Competing Interests: Declaration of interests M.M. is listed as the inventor on the patent for the EGFR mutation analysis for lung cancer diagnosis licensed to LabCorp; holds research support from Bayer, Janssen, and Ono; and serves as scientific advisory board chair for OrigiMed. P.W.L. serves on the scientific advisory boards of AnchorDx and Progenity. D.J.K. has research support from Revolution Medicines and Genentech and is a consultant to AADi., (Copyright © 2021 The Authors. Published by Elsevier Inc. All rights reserved.)

The Trans-Omics for Precision Medicine (TOPMed) programme seeks to elucidate the genetic architecture and biology of heart, lung, blood and sleep disorders, with the ultimate goal of improving diagnosis, treatment and prevention of these diseases. The initial phases of the programme focused on whole-genome sequencing of individuals with rich phenotypic data and diverse backgrounds. Here we describe the TOPMed goals and design as well as the available resources and early insights obtained from the sequence data. The resources include a variant browser, a genotype imputation server, and genomic and phenotypic data that are available through dbGaP (Database of Genotypes and Phenotypes) 1 . In the first 53,831 TOPMed samples, we detected more than 400 million single-nucleotide and insertion or deletion variants after alignment with the reference genome. Additional previously undescribed variants were detected through assembly of unmapped reads and customized analysis in highly variable loci. Among the more than 400 million detected variants, 97% have frequencies of less than 1% and 46% are singletons that are present in only one individual (53% among unrelated individuals). These rare variants provide insights into mutational processes and recent human evolutionary history. The extensive catalogue of genetic variation in TOPMed studies provides unique opportunities for exploring the contributions of rare and noncoding sequence variants to phenotypic variation. Furthermore, combining TOPMed haplotypes with modern imputation methods improves the power and reach of genome-wide association studies to include variants down to a frequency of approximately 0.01%.

Purpose: POAG is the leading cause of irreversible blindness in African Americans. In this study, we quantitatively assess the association of autosomal ancestry with POAG risk in a large cohort of self-identified African Americans., Methods: Subjects recruited to the Primary Open-Angle African American Glaucoma Genetics (POAAGG) study were classified as glaucoma cases or controls by fellowship-trained glaucoma specialists. POAAGG subjects were genotyped using the MEGA Ex array (discovery cohort, n = 3830; replication cohort, n = 2135). Population structure was interrogated using principal component analysis in the context of the 1000 Genomes Project superpopulations., Results: The majority of POAAGG samples lie on an axis between African and European superpopulations, with great variation in admixture. Cases had a significantly lower mean value of the ancestral component q0 than controls for both cohorts (P = 6.14-4; P = 3-6), consistent with higher degree of African ancestry. Among POAG cases, higher African ancestry was also associated with thinner central corneal thickness (P = 2-4). Admixture mapping showed that local genetic ancestry was not a significant risk factor for POAG. A polygenic risk score, comprised of 23 glaucoma-associated single nucleotide polymorphisms from the NHGRI-EBI genome-wide association study catalog, was significant in both cohorts (P < 0.001), suggesting that both known POAG single nucleotide polymorphisms and an omnigenic ancestry effect influence POAG risk., Conclusions: In sum, the POAAGG study population is very admixed, with a higher degree of African ancestry associated with an increased POAG risk. Further analyses should consider social and environmental factors as possible confounding factors for disease predisposition.

Coronavirus disease 2019 (COVID-19) is caused by SARS-CoV-2, an emerging virus that utilizes host proteins ACE2 and TMPRSS2 as entry factors. Understanding the factors affecting the pattern and levels of expression of these genes is important for deeper understanding of SARS-CoV-2 tropism and pathogenesis. Here we explore the role of genetics and co-expression networks in regulating these genes in the airway, through the analysis of nasal airway transcriptome data from 695 children. We identify expression quantitative trait loci for both ACE2 and TMPRSS2, that vary in frequency across world populations. We find TMPRSS2 is part of a mucus secretory network, highly upregulated by type 2 (T2) inflammation through the action of interleukin-13, and that the interferon response to respiratory viruses highly upregulates ACE2 expression. IL-13 and virus infection mediated effects on ACE2 expression were also observed at the protein level in the airway epithelium. Finally, we define airway responses to common coronavirus infections in children, finding that these infections generate host responses similar to other viral species, including upregulation of IL6 and ACE2. Our results reveal possible mechanisms influencing SARS-CoV-2 infectivity and COVID-19 clinical outcomes.

Rationale: Puerto Ricans have the highest childhood asthma prevalence in the United States (23.6%); however, the etiology is uncertain. Objectives: In this study, we sought to uncover the genetic architecture of lung function in Puerto Rican youth with and without asthma who were recruited from the island ( n  = 836). Methods: We used admixture-mapping and whole-genome sequencing data to discover genomic regions associated with lung function. Functional roles of the prioritized candidate SNPs were examined with chromatin immunoprecipitation sequencing, RNA sequencing, and expression quantitative trait loci data. Measurements and Main Results: We discovered a genomic region at 1q32 that was significantly associated with a 0.12-L decrease in the lung volume of exhaled air (95% confidence interval, -0.17 to -0.07; P  = 6.62 × 10 -8 ) with each allele of African ancestry. Within this region, two SNPs were expression quantitative trait loci of TMEM9 in nasal airway epithelial cells and MROH3P in esophagus mucosa. The minor alleles of these SNPs were associated with significantly decreased lung function and decreased TMEM9 gene expression. Another admixture-mapping peak was observed on chromosome 5q35.1, indicating that each Native American ancestry allele was associated with a 0.15-L increase in lung function (95% confidence interval, 0.08-0.21; P  = 5.03 × 10 -6 ). The region-based association tests identified four suggestive windows that harbored candidate rare variants associated with lung function. Conclusions: We identified common and rare genetic variants that may play a critical role in lung function among Puerto Rican youth. We independently validated an inflammatory pathway that could potentially be used to develop more targeted treatments and interventions for patients with asthma.

Cancer genomes often harbor hundreds of somatic DNA rearrangement junctions, many of which cannot be easily classified into simple (e.g., deletion) or complex (e.g., chromothripsis) structural variant classes. Applying a novel genome graph computational paradigm to analyze the topology of junction copy number (JCN) across 2,778 tumor whole-genome sequences, we uncovered three novel complex rearrangement phenomena: pyrgo, rigma, and tyfonas. Pyrgo are "towers" of low-JCN duplications associated with early-replicating regions, superenhancers, and breast or ovarian cancers. Rigma comprise "chasms" of low-JCN deletions enriched in late-replicating fragile sites and gastrointestinal carcinomas. Tyfonas are "typhoons" of high-JCN junctions and fold-back inversions associated with expressed protein-coding fusions, breakend hypermutation, and acral, but not cutaneous, melanomas. Clustering of tumors according to genome graph-derived features identified subgroups associated with DNA repair defects and poor prognosis., Competing Interests: Declaration of Interests J.S.R.-F. reports receiving personal/consultancy fees from VolitionRx, Paige.AI, Goldman Sachs, REPARE Therapeutics, GRAIL, Ventana Medical Systems, Roche, Genentech, and InviCRO outside of the scope of the submitted work., (Copyright © 2020 Elsevier Inc. All rights reserved.)

A key goal of whole-genome sequencing for studies of human genetics is to interrogate all forms of variation, including single-nucleotide variants, small insertion or deletion (indel) variants and structural variants. However, tools and resources for the study of structural variants have lagged behind those for smaller variants. Here we used a scalable pipeline 1 to map and characterize structural variants in 17,795 deeply sequenced human genomes. We publicly release site-frequency data to create the largest, to our knowledge, whole-genome-sequencing-based structural variant resource so far. On average, individuals carry 2.9 rare structural variants that alter coding regions; these variants affect the dosage or structure of 4.2 genes and account for 4.0-11.2% of rare high-impact coding alleles. Using a computational model, we estimate that structural variants account for 17.2% of rare alleles genome-wide, with predicted deleterious effects that are equivalent to loss-of-function coding alleles; approximately 90% of such structural variants are noncoding deletions (mean 19.1 per genome). We report 158,991 ultra-rare structural variants and show that 2% of individuals carry ultra-rare megabase-scale structural variants, nearly half of which are balanced or complex rearrangements. Finally, we infer the dosage sensitivity of genes and noncoding elements, and reveal trends that relate to element class and conservation. This work will help to guide the analysis and interpretation of structural variants in the era of whole-genome sequencing.

Baseline lung function, quantified as forced expiratory volume in the first second of exhalation (FEV 1 ), is a standard diagnostic criterion used by clinicians to identify and classify lung diseases. Using whole-genome sequencing data from the National Heart, Lung, and Blood Institute Trans-Omics for Precision Medicine project, we identified a novel genetic association with FEV 1 on chromosome 12 in 867 African American children with asthma ( P = 1.26 × 10 -8 , β = 0.302). Conditional analysis within 1 Mb of the tag signal (rs73429450) yielded one major and two other weaker independent signals within this peak. We explored statistical and functional evidence for all variants in linkage disequilibrium with the three independent signals and yielded nine variants as the most likely candidates responsible for the association with FEV 1 Hi-C data and expression QTL analysis demonstrated that these variants physically interacted with KITLG (KIT ligand, also known as SCF ), and their minor alleles were associated with increased expression of the KITLG gene in nasal epithelial cells. Gene-by-air-pollution interaction analysis found that the candidate variant rs58475486 interacted with past-year ambient sulfur dioxide exposure ( P = 0.003, β = 0.32). This study identified a novel protective genetic association with FEV 1 , possibly mediated through KITLG , in African American children with asthma. This is the first study that has identified a genetic association between lung function and KITLG , which has established a role in orchestrating allergic inflammation in asthma., (Copyright © 2020 by the Genetics Society of America.)

To test the performance of a new sequencing platform, develop an updated somatic calling pipeline and establish a reference for future benchmarking experiments, we performed whole-genome sequencing of 3 common cancer cell lines (COLO-829, HCC-1143 and HCC-1187) along with their matched normal cell lines to great sequencing depths (up to 278x coverage) on both Illumina HiSeqX and NovaSeq sequencing instruments. Somatic calling was generally consistent between the two platforms despite minor differences at the read level. We designed and implemented a novel pipeline for the analysis of tumor-normal samples, using multiple variant callers. We show that coupled with a high-confidence filtering strategy, the use of combination of tools improves the accuracy of somatic variant calling. We also demonstrate the utility of the dataset by creating an artificial purity ladder to evaluate the somatic pipeline and benchmark methods for estimating purity and ploidy from tumor-normal pairs. The data and results of the pipeline are made accessible to the cancer genomics community.

Purpose: To maximize the discovery of potentially pathogenic variants to better understand the diagnostic utility of genome sequencing (GS) and to assess how the presence of multiple risk events might affect the phenotypic severity in autism spectrum disorders (ASD)., Methods: GS was applied to 180 simplex and multiplex ASD families (578 individuals, 213 patients) with exome sequencing and array comparative genomic hybridization further applied to a subset for validation and cross-platform comparisons., Results: We found that 40.8% of patients carried variants with evidence of disease risk, including a de novo frameshift variant in NR4A2 and two de novo missense variants in SYNCRIP, while 21.1% carried clinically relevant pathogenic or likely pathogenic variants. Patients with more than one risk variant (9.9%) were more severely affected with respect to cognitive ability compared with patients with a single or no-risk variant. We observed no instance among the 27 multiplex families where a pathogenic or likely pathogenic variant was transmitted to all affected members in the family., Conclusion: The study demonstrates the diagnostic utility of GS, especially for multiple risk variants that contribute to the phenotypic severity, shows the genetic heterogeneity in multiplex families, and provides evidence for new genes for follow up.

Background: Prompted by the revolution in high-throughput sequencing and its potential impact for treating cancer patients, we initiated a clinical research study to compare the ability of different sequencing assays and analysis methods to analyze glioblastoma tumors and generate real-time potential treatment options for physicians., Methods: A consortium of seven institutions in New York City enrolled 30 patients with glioblastoma and performed tumor whole genome sequencing (WGS) and RNA sequencing (RNA-seq; collectively WGS/RNA-seq); 20 of these patients were also analyzed with independent targeted panel sequencing. We also compared results of expert manual annotations with those from an automated annotation system, Watson Genomic Analysis (WGA), to assess the reliability and time required to identify potentially relevant pharmacologic interventions., Results: WGS/RNAseq identified more potentially actionable clinical results than targeted panels in 90% of cases, with an average of 16-fold more unique potentially actionable variants identified per individual; 84 clinically actionable calls were made using WGS/RNA-seq that were not identified by panels. Expert annotation and WGA had good agreement on identifying variants [mean sensitivity = 0.71, SD = 0.18 and positive predictive value (PPV) = 0.80, SD = 0.20] and drug targets when the same variants were called (mean sensitivity = 0.74, SD = 0.34 and PPV = 0.79, SD = 0.23) across patients. Clinicians used the information to modify their treatment plan 10% of the time., Conclusion: These results present the first comprehensive comparison of technical and machine augmented analysis of targeted panel and WGS/RNA-seq to identify potential cancer treatments.

Hundreds of thousands of human whole genome sequencing (WGS) datasets will be generated over the next few years. These data are more valuable in aggregate: joint analysis of genomes from many sources increases sample size and statistical power. A central challenge for joint analysis is that different WGS data processing pipelines cause substantial differences in variant calling in combined datasets, necessitating computationally expensive reprocessing. This approach is no longer tenable given the scale of current studies and data volumes. Here, we define WGS data processing standards that allow different groups to produce functionally equivalent (FE) results, yet still innovate on data processing pipelines. We present initial FE pipelines developed at five genome centers and show that they yield similar variant calling results and produce significantly less variability than sequencing replicates. This work alleviates a key technical bottleneck for genome aggregation and helps lay the foundation for community-wide human genetics studies.

DNA methylation patterns in the genome both reflect and help to mediate transcriptional regulatory processes. The digital nature of DNA methylation, present or absent on each allele, makes this assay capable of quantifying events in subpopulations of cells, whereas genome-wide chromatin studies lack the same quantitative capacity. Testing DNA methylation throughout the genome is possible using whole-genome bisulfite sequencing (WGBS), but the high costs associated with the assay have made it impractical for studies involving more than limited numbers of samples. We have optimized a new transposase-based library preparation assay for the Illumina HiSeq X platform suitable for limited amounts of DNA and providing a major cost reduction for WGBS. By incorporating methylated cytosines during fragment end repair, we reveal an end-repair artifact affecting 1%-2% of reads that we can remove analytically. We show that the use of a high (G + C) content spike-in performs better than PhiX in terms of bisulfite sequencing quality. As expected, the loci with transposase-accessible chromatin are DNA hypomethylated and enriched in flanking regions by post-translational modifications of histones usually associated with positive effects on gene expression. Using these transposase-accessible loci to represent the cis -regulatory loci in the genome, we compared the representation of these loci between WGBS and other genome-wide DNA methylation assays, showing WGBS to outperform substantially all of the alternatives. We conclude that it is now technologically and financially feasible to perform WGBS in larger numbers of samples with greater accuracy than previously possible., (© 2018 Suzuki et al.; Published by Cold Spring Harbor Laboratory Press.)

High-throughput sequencing is a revolutionary technology for the analysis of metagenomic samples. However, querying large volumes of reads against comprehensive DNA/RNA databases in a sensitive manner can be compute-intensive. Here, we present taxMaps, a highly efficient, sensitive, and fully scalable taxonomic classification tool. Using a combination of simulated and real metagenomics data sets, we demonstrate that taxMaps is more sensitive and more precise than widely used taxonomic classifiers and is capable of delivering classification accuracy comparable to that of BLASTN, but at up to three orders of magnitude less computational cost., (© 2018 Corvelo et al.; Published by Cold Spring Harbor Laboratory Press.)

Reliable detection of somatic variations is of critical importance in cancer research. Here we present Lancet, an accurate and sensitive somatic variant caller, which detects SNVs and indels by jointly analyzing reads from tumor and matched normal samples using colored de Bruijn graphs. We demonstrate, through extensive experimental comparison on synthetic and real whole-genome sequencing datasets, that Lancet has better accuracy, especially for indel detection, than widely used somatic callers, such as MuTect, MuTect2, LoFreq, Strelka, and Strelka2. Lancet features a reliable variant scoring system, which is essential for variant prioritization, and detects low-frequency mutations without sacrificing the sensitivity to call longer insertions and deletions empowered by the local-assembly engine. In addition to genome-wide analysis, Lancet allows inspection of somatic variants in graph space, which augments the traditional read alignment visualization to help confirm a variant of interest. Lancet is available as an open-source program at https://github.com/nygenome/lancet., Competing Interests: The authors declare no competing interests.

To identify novel genes associated with ALS, we undertook two lines of investigation. We carried out a genome-wide association study comparing 20,806 ALS cases and 59,804 controls. Independently, we performed a rare variant burden analysis comparing 1,138 index familial ALS cases and 19,494 controls. Through both approaches, we identified kinesin family member 5A (KIF5A) as a novel gene associated with ALS. Interestingly, mutations predominantly in the N-terminal motor domain of KIF5A are causative for two neurodegenerative diseases: hereditary spastic paraplegia (SPG10) and Charcot-Marie-Tooth type 2 (CMT2). In contrast, ALS-associated mutations are primarily located at the C-terminal cargo-binding tail domain and patients harboring loss-of-function mutations displayed an extended survival relative to typical ALS cases. Taken together, these results broaden the phenotype spectrum resulting from mutations in KIF5A and strengthen the role of cytoskeletal defects in the pathogenesis of ALS., (Copyright © 2018 Elsevier Inc. All rights reserved.)