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

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

Capsule carbohydrate structure determines virulence in Acinetobacter baumannii.

Acinetobacter baumannii is a highly antibiotic-resistant bacterial pathogen for which novel therapeutic approaches are needed. Unfortunately, the drivers of virulence in A. baumannii remain uncertain. By comparing genomes among a panel of A. baumannii strains we identified a specific gene variation in the capsule locus that correlated with altered virulence. While less virulent strains possessed the intact gene gtr6, a hypervirulent clinical isolate contained a spontaneous transposon insertion in the same gene, resulting in the loss of a branchpoint in capsular carbohydrate structure. By constructing isogenic gtr6 mutants, we confirmed that gtr6-disrupted strains were protected from phagocytosis in vitro and displayed higher bacterial burden and lethality in vivo. Gtr6+ strains were phagocytized more readily and caused lower bacterial burden and no clinical illness in vivo. We found that the CR3 receptor mediated phagocytosis of gtr6+, but not gtr6-, strains in a complement-dependent manner. Furthermore, hypovirulent gtr6+ strains demonstrated increased virulence in vivo when CR3 function was abrogated. In summary, loss-of-function in a single capsule assembly gene dramatically altered virulence by inhibiting complement deposition and recognition by phagocytes across multiple A. baumannii strains. Thus, capsular structure can determine virulence among A. baumannii strains by altering bacterial interactions with host complement-mediated opsonophagocytosis.

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

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

Abnormal SCID Newborn Screening and Spontaneous Recovery Associated with a Novel Haploinsufficiency IKZF1 Mutation.

PURPOSE: IKAROS, encoded by IKZF1, is a member of the IKAROS family of zinc-finger transcription factors playing critical roles in lymphocyte development, differentiation, and tumor suppression. Several studies demonstrated that IKZF1 mutations affecting DNA binding or homo-/hetero-dimerization are mostly associated with common variable immunodeficiency, combined immunodeficiency, or hematologic manifestations. Herein we report a likely de novo, nonsense IKZF1 mutation (p.C182*) in a baby with low T cell receptor excision circles (TREC) identified by newborn screening testing for severe combined immunodeficiency. The patient also presented a profound B cell deficiency at birth. METHODS: Genetic, functional, immunologic, and clinical outcome data associated with this patient and her mutation were evaluated. RESULTS: Mutant p.C182* was detected in the cytoplasm of the patient's primary cells, in contrast to wild type (WT) IKAROS protein, only detected in the nucleus. Functional in vitro assessments revealed that p.C182* was less stable than WT IKAROS protein and failed to bind to its target DNA binding sequence and dimerize with WT IKAROS protein, resulting in impaired pericentromeric targeting and transcriptional repression by means of haploinsufficiency. During follow-up, while a spontaneous recovery of TREC and T cells was observed, B cells improved but not to sustained normal ranges. CONCLUSIONS: Patients with IKAROS-associated diseases can present with SCID-like TREC values through newborn screening testing. IKZF1 mutations should be added to the low TREC differential, although spontaneous recovery has to be considered.

CRISPR/Cas9 gene editing for the creation of an MGAT1-deficient CHO cell line to control HIV-1 vaccine glycosylation.

Over the last decade, multiple broadly neutralizing monoclonal antibodies (bN-mAbs) to the HIV-1 envelope protein (Env) gp120 have been described. Many of these recognize epitopes consisting of both amino acid and glycan residues. Moreover, the glycans required for binding of these bN-mAbs are early intermediates in the N-linked glycosylation pathway. This type of glycosylation substantially alters the mass and net charge of Envs compared to molecules with the same amino acid sequence but possessing mature, complex (sialic acid-containing) carbohydrates. Since cell lines suitable for biopharmaceutical production that limit N-linked glycosylation to mannose-5 (Man5) or earlier intermediates are not readily available, the production of vaccine immunogens displaying these glycan-dependent epitopes has been challenging. Here, we report the development of a stable suspension-adapted Chinese hamster ovary (CHO) cell line that limits glycosylation to Man5 and earlier intermediates. This cell line was created using the clustered regularly interspaced short palindromic repeat (CRISPR)/CRISPR-associated protein 9 (Cas9) gene editing system and contains a mutation that inactivates the gene encoding Mannosyl (Alpha-1,3-)-Glycoprotein Beta-1,2-N-Acetylglucosaminyltransferase (MGAT1). Monomeric gp120s produced in the MGAT1- CHO cell line exhibit improved binding to prototypic glycan-dependent bN-mAbs directed to the V1/V2 domain (e.g., PG9) and the V3 stem (e.g., PGT128 and 10-1074) while preserving the structure of the important glycan-independent epitopes (e.g., VRC01). The ability of the MGAT1- CHO cell line to limit glycosylation to early intermediates in the N-linked glycosylation pathway without impairing the doubling time or ability to grow at high cell densities suggests that it will be a useful substrate for the biopharmaceutical production of HIV-1 vaccine immunogens.

Identification and CRISPR/Cas9 Inactivation of the C1s Protease Responsible for Proteolysis of Recombinant Proteins Produced in CHO Cells.

Proteolysis associated with recombinant protein expression in Chinese Hamster Ovary (CHO) cells has hindered the development of biologics including HIV vaccines. When expressed in CHO cells, the recombinant HIV envelope protein, gp120, undergoes proteolytic clipping by a serine protease at a key epitope recognized by neutralizing antibodies. The problem is particularly acute for envelope proteins from clade B viruses that represent the major genetic subtype circulating in much of the developed world, including the US and Europe. In this paper, we have identified complement Component 1's (C1s), a serine protease from the complement cascade, as the protease responsible for the proteolysis of gp120 in CHO cells. CRISPR/Cas9 knockout of the C1s protease in a CHO cell line was shown to eliminate the proteolytic activity against the recombinantly expressed gp120. In addition, the C1s-/- MGAT1- CHO cell line, with the C1s protease and the MGAT1 glycosyltransferase knocked out, enabled the production of unclipped gp120 from a clade B isolate (BaL-rgp120) and enriched for mannose-5 glycans on gp120 that are required for the binding of multiple broadly neutralizing monoclonal antibodies (bN-mAbs). The availability of this technology will allow for the scale-up and testing of multiple vaccine concepts in regions of the world where clade B viruses are in circulation. Furthermore, the proteolysis issues caused by the C1s protease suggests a broader need for a C1s-deficient CHO cell line to express other recombinant proteins that are susceptible to serine protease activity in CHO cells. Similarly, the workflow described here to identify and knockout C1s in a CHO cell line can be applied to remedy the proteolysis of biologics by other CHO proteases.

Assessment of Insertion Sequence Mobilization as an Adaptive Response to Oxidative Stress in Acinetobacter baumannii Using IS-seq.

Insertion sequence (IS) elements are found throughout bacterial genomes and contribute to genome variation by interrupting genes or altering gene expression. Few of the more than 30 IS elements described in Acinetobacter baumannii have been characterized for transposition activity or expression effects. A targeted sequencing method, IS-seq, was developed to efficiently map the locations of new insertion events in A. baumannii genomes and was used to identify novel IS sites following growth in the presence of hydrogen peroxide, which causes oxidative stress. Serial subculture in the presence of subinhibitory concentrations of hydrogen peroxide led to rapid selection of cells carrying an ISAba1 element upstream of the catalase-peroxidase gene katG Several additional sites for the elements ISAba1, ISAba13, ISAba25, ISAba26, and ISAba125 were found at low abundance after serial subculture, indicating that each element is active and contributes to genetic variation that may be subject to selection. Following hydrogen peroxide exposure, rapid changes in gene expression were observed in genes related to iron homeostasis. The IS insertions adjacent to katG resulted in more than 20-fold overexpression of the gene and increased hydrogen peroxide tolerance.IMPORTANCE Insertion sequences (IS) contribute to genomic and phenotypic variation in many bacterial species, but little is known about how transposition rates vary among elements or how selective pressure influences this process. A new method for identifying new insertion locations that arise under experimental growth conditions in the genome, termed IS-seq, was developed and tested with cells grown in the presence of hydrogen peroxide, which causes oxidative stress. Gene expression changes in response to hydrogen peroxide exposure are similar to those observed in other species and include genes that control free iron concentrations. New IS insertions adjacent to a gene encoding a catalase enzyme confirm that IS elements can rapidly contribute to adaptive variation in the presence of selection.

Chemical Genetic Interaction Profiling Reveals Determinants of Intrinsic Antibiotic Resistance in Mycobacterium tuberculosis.

Chemotherapy for tuberculosis (TB) is lengthy and could benefit from synergistic adjuvant therapeutics that enhance current and novel drug regimens. To identify genetic determinants of intrinsic antibiotic susceptibility in Mycobacterium tuberculosis, we applied a chemical genetic interaction (CGI) profiling approach. We screened a saturated transposon mutant library and identified mutants that exhibit altered fitness in the presence of partially inhibitory concentrations of rifampin, ethambutol, isoniazid, vancomycin, and meropenem, antibiotics with diverse mechanisms of action. This screen identified the M. tuberculosis cell envelope to be a major determinant of antibiotic susceptibility but did not yield mutants whose increase in susceptibility was due to transposon insertions in genes encoding efflux pumps. Intrinsic antibiotic resistance determinants affecting resistance to multiple antibiotics included the peptidoglycan-arabinogalactan ligase Lcp1, the mycolic acid synthase MmaA4, the protein translocase SecA2, the mannosyltransferase PimE, the cell envelope-associated protease CaeA/Hip1, and FecB, a putative iron dicitrate-binding protein. Characterization of a deletion mutant confirmed FecB to be involved in the intrinsic resistance to every antibiotic analyzed. In contrast to its predicted function, FecB was dispensable for growth in low-iron medium and instead functioned as a critical mediator of envelope integrity.

NDM-5 and OXA-181 Beta-Lactamases, a Significant Threat Continues To Spread in the Americas.

Among Gram-negative bacteria, carbapenem-resistant infections pose a serious and life-threatening challenge. Here, the CRACKLE network reports a sentinel detection and characterization of a carbapenem-resistant Klebsiella pneumoniae ST147 isolate harboring blaNDM-5 and blaOXA-181 from a young man who underwent abdominal surgery in India. blaNDM-5 was located on an IncFII plasmid of ≈90 kb, whereas blaOXA-181 was chromosomally encoded. Resistome and genome analysis demonstrated multiple copies of the transposable element IS26 and a "hot-spot region" in the IncFII plasmid.

Nosocomial Outbreak of Extensively Drug-Resistant Acinetobacter baumannii Isolates Containing blaOXA-237 Carried on a Plasmid.

Carbapenem antibiotics are among the mainstays for treating infections caused by Acinetobacter baumannii, especially in the Northwest United States, where carbapenem-resistant A. baumannii remains relatively rare. However, between June 2012 and October 2014, an outbreak of carbapenem-resistant A. baumannii occurred in 16 patients from five health care facilities in the state of Oregon. All isolates were defined as extensively drug resistant. Multilocus sequence typing revealed that the isolates belonged to sequence type 2 (international clone 2 [IC2]) and were >95% similar as determined by repetitive-sequence-based PCR analysis. Multiplex PCR revealed the presence of a blaOXA carbapenemase gene, later identified as blaOXA-237 Whole-genome sequencing of all isolates revealed a well-supported separate branch within a global A. baumannii phylogeny. Pacific Biosciences (PacBio) SMRT sequencing was also performed on one isolate to gain insight into the genetic location of the carbapenem resistance gene. We discovered that blaOXA-237, flanked on either side by ISAba1 elements in opposite orientations, was carried on a 15,198-bp plasmid designated pORAB01-3 and was present in all 16 isolates. The plasmid also contained genes encoding a TonB-dependent receptor, septicolysin, a type IV secretory pathway (VirD4 component, TraG/TraD family) ATPase, an integrase, a RepB family plasmid DNA replication initiator protein, an alpha/beta hydrolase, and a BrnT/BrnA type II toxin-antitoxin system. This is the first reported outbreak in the northwestern United States associated with this carbapenemase. Particularly worrisome is that blaOXA-237 was carried on a plasmid and found in the most prominent worldwide clonal group IC2, potentially giving pORAB01-3 great capacity for future widespread dissemination.

Initial Assessment of the Molecular Epidemiology of blaNDM-1 in Colombia.

We report complete genome sequences of four blaNDM-1-harboring Gram-negative multidrug-resistant (MDR) isolates from Colombia. The blaNDM-1 genes were located on 193-kb Inc FIA, 178-kb Inc A/C2, and 47-kb (unknown Inc type) plasmids. Multilocus sequence typing (MLST) revealed that these isolates belong to sequence type 10 (ST10) (Escherichia coli), ST392 (Klebsiella pneumoniae), and ST322 and ST464 (Acinetobacter baumannii and Acinetobacter nosocomialis, respectively). Our analysis identified that the Inc A/C2 plasmid in E. coli contained a novel complex transposon (Tn125 and Tn5393 with three copies of blaNDM-1) and a recombination "hot spot" for the acquisition of new resistance determinants.

Host fate is rapidly determined by innate effector-microbial interactions during Acinetobacter baumannii bacteremia.

BACKGROUND: Acinetobacter baumannii is one of the most antibiotic-resistant pathogens. Defining mechanisms driving pathogenesis is critical to enable new therapeutic approaches. METHODS: We studied virulence differences across a diverse panel of A. baumannii clinical isolates during murine bacteremia to elucidate host-microbe interactions that drive outcome. RESULTS: We identified hypervirulent strains that were lethal at low intravenous inocula and achieved very high early, and persistent, blood bacterial densities. Virulent strains were nonlethal at low inocula but lethal at 2.5-fold higher inocula. Finally, relatively avirulent (hypovirulent) strains were nonlethal at 20-fold higher inocula and were efficiently cleared by early time points. In vivo virulence correlated with in vitro resistance to complement and macrophage uptake. Depletion of complement, macrophages, and neutrophils each independently increased bacterial density of the hypovirulent strain but insufficiently to change lethality. However, disruption of all 3 effector mechanisms enabled early bacterial densities similar to hypervirulent strains, rendering infection 100% fatal. CONCLUSIONS: The lethality of A. baumannii strains depends on distinct stages. Strains resistant to early innate effectors are able to establish very high early bacterial blood density, and subsequent sustained bacteremia leads to Toll-like receptor 4-mediated hyperinflammation and lethality. These results have important implications for translational efforts to develop therapies that modulate host-microbe interactions.

Genomic and transcriptomic analyses of colistin-resistant clinical isolates of Klebsiella pneumoniae reveal multiple pathways of resistance.

The emergence of multidrug-resistant (MDR) Klebsiella pneumoniae has resulted in a more frequent reliance on treatment using colistin. However, resistance to colistin (Col(r)) is increasingly reported from clinical settings. The genetic mechanisms that lead to Col(r) in K. pneumoniae are not fully characterized. Using a combination of genome sequencing and transcriptional profiling by RNA sequencing (RNA-Seq) analysis, distinct genetic mechanisms were found among nine Col(r) clinical isolates. Col(r) was related to mutations in three different genes in K. pneumoniae strains, with distinct impacts on gene expression. Upregulation of the pmrH operon encoding 4-amino-4-deoxy-L-arabinose (Ara4N) modification of lipid A was found in all Col(r) strains. Alteration of the mgrB gene was observed in six strains. One strain had a mutation in phoQ. Common among these seven strains was elevated expression of phoPQ and unaltered expression of pmrCAB, which is involved in phosphoethanolamine addition to lipopolysaccharide (LPS). In two strains, separate mutations were found in a previously uncharacterized histidine kinase gene that is part of a two-component regulatory system (TCRS) now designated crrAB. In these strains, expression of pmrCAB, crrAB, and an adjacent glycosyltransferase gene, but not that of phoPQ, was elevated. Complementation with the wild-type allele restored colistin susceptibility in both strains. The crrAB genes are present in most K. pneumoniae genomes, but not in Escherichia coli. Additional upregulated genes in all strains include those involved in cation transport and maintenance of membrane integrity. Because the crrAB genes are present in only some strains, Col(r) mechanisms may be dependent on the genetic background.

Population structure of KPC-producing Klebsiella pneumoniae isolates from midwestern U.S. hospitals.

Genome sequencing of carbapenem-resistant Klebsiella pneumoniae isolates from regional U.S. hospitals was used to characterize strain diversity and the bla(KPC) genetic context. A phylogeny based on core single-nucleotide variants (SNVs) supports a division of sequence type 258 (ST258) into two distinct groups. The primary differences between the groups are in the capsular polysaccharide locus (cps) and their plasmid contents. A strict association between clade and KPC variant was found. The bla(KPC) gene was found on variants of two plasmid backbones. This study indicates that highly similar K. pneumoniae subpopulations coexist within the same hospitals over time.

Spliceosome malfunction causes neurodevelopmental disorders with overlapping features.

Pre-mRNA splicing is a highly coordinated process. While its dysregulation has been linked to neurological deficits, our understanding of the underlying molecular and cellular mechanisms remains limited. We implicated pathogenic variants in U2AF2 and PRPF19, encoding spliceosome subunits in neurodevelopmental disorders (NDDs), by identifying 46 unrelated individuals with 23 de novo U2AF2 missense variants (including 7 recurrent variants in 30 individuals) and 6 individuals with de novo PRPF19 variants. Eight U2AF2 variants dysregulated splicing of a model substrate. Neuritogenesis was reduced in human neurons differentiated from human pluripotent stem cells carrying two U2AF2 hyper-recurrent variants. Neural loss of function (LoF) of the Drosophila orthologs U2af50 and Prp19 led to lethality, abnormal mushroom body (MB) patterning, and social deficits, which were differentially rescued by wild-type and mutant U2AF2 or PRPF19. Transcriptome profiling revealed splicing substrates or effectors (including Rbfox1, a third splicing factor), which rescued MB defects in U2af50-deficient flies. Upon reanalysis of negative clinical exomes followed by data sharing, we further identified 6 patients with NDD who carried RBFOX1 missense variants which, by in vitro testing, showed LoF. Our study implicates 3 splicing factors as NDD-causative genes and establishes a genetic network with hierarchy underlying human brain development and function.

Report of two cases of Schaaf-Yang syndrome: Same genotype and different phenotype.

We report two, genotypically identical but phenotypically distinct cases of Schaaf-Yang syndrome and propose the early use of Genome Sequencing in patients with nonspecific presentations to facilitate the early diagnosis of children with rare genetic diseases and improve overall health care outcomes.

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

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