Pamela J. Thomas

Comparative analyses of multi-species sequences from targeted genomic regions

The systematic comparison of genomic sequences from different organisms represents a central focus of contemporary genome analysis. Comparative analyses of vertebrate sequences can identify coding and conserved non-coding regions, including regulatory elements, and provide insight into the forces that have rendered modern-day genomes. As a complement to whole-genome sequencing efforts, we are sequencing and comparing targeted genomic regions in multiple, evolutionarily diverse vertebrates. Here we report the generation and analysis of over 12 megabases (Mb) of sequence from 12 species, all derived from the genomic region orthologous to a segment of about 1.8 Mb on human chromosome 7 containing ten genes, including the gene mutated in cystic fibrosis. These sequences show conservation reflecting both functional constraints and the neutral mutational events that shaped this genomic region. In particular, we identify substantial numbers of conserved non-coding segments beyond those previously identified experimentally, most of which are not detectable by pair-wise sequence comparisons alone. Analysis of transposable element insertions highlights the variation in genome dynamics among these species and confirms the placement of rodents as a sister group to the primates.

Tracking a Hospital Outbreak of Carbapenem-Resistant Klebsiella pneumoniae with Whole-Genome Sequencing

Tracking a hospital outbreak of carbapenem-resistant Klebsiella pneumoniae with whole-genome sequencing revealed its origin and probable modes of transmission. A Detective Story Some infections are largely a thing of the past—plague, syphilis. The unfortunate result of these antibiotic-driven successes is the emergence of drug-resistant pathogens. And, ironically enough, hospitals are at the center of the problem. An example of this occurred in 2011 at the Clinical Center of the U.S. National Institutes of Health (NIH), in which an outbreak of drug-resistant Klebsiella pneumoniae infected 18 patients, causing the death of 6 of them. Using a combination of whole-genome sequencing and patient tracking, Snitkin and his colleagues examined how the bacteria was spreading through the hospital. The results outline a complicated path of transmission within the hospital that defied standard containment methods, yielding lessons for the future. A patient known to be infected with a drug-resistant form of K. pneumoniae was admitted to the NIH Clinical Center on 13 June 2011. Enhanced isolation procedures were immediately implemented, and no spread of the bacteria was seen for the month she was in the hospital. Although all seemed well, a few weeks later on August 5th, a second infected patient was discovered, followed by a series of other patients with infection or colonization—about 1 a week to a total of 18 by the end of 2011. Six people ultimately died as a result of the bacteria. The outbreak was finally contained by rigorous control procedures. A careful survey of the bed locations of each patient did not shed much light on how the bacteria traveled on its deadly path: The first patient did not even come into contact with any of the others. So the authors performed whole-genome sequencing on all of the bacteria that were found, determining the most likely evolutionary relationships among them by comparing the variations at single nucleotides that arise as bacteria grow. Combining this evolutionary information with the physical tracking of the patients pointed to the most likely transmission scenario. The authors concluded that all of the K. pneumoniae cases likely originated with the index patient, from at least two different sites on her body, rather than by independently introduced bacteria. There were at least three different initial transmission events. Particularly disturbing was the fact that one of the infections could be linked to contamination of a ventilator that had been cleaned by thorough methods. Sophisticated deployment of whole-genome sequencing revealed the weaknesses in this medical who-done-it, informing improvements in hospital preventive measures. If applied rapidly, such analysis can even expose the causes of nosocomial infections in real time. The Gram-negative bacteria Klebsiella pneumoniae is a major cause of nosocomial infections, primarily among immunocompromised patients. The emergence of strains resistant to carbapenems has left few treatment options, making infection containment critical. In 2011, the U.S. National Institutes of Health Clinical Center experienced an outbreak of carbapenem-resistant K. pneumoniae that affected 18 patients, 11 of whom died. Whole-genome sequencing was performed on K. pneumoniae isolates to gain insight into why the outbreak progressed despite early implementation of infection control procedures. Integrated genomic and epidemiological analysis traced the outbreak to three independent transmissions from a single patient who was discharged 3 weeks before the next case became clinically apparent. Additional genomic comparisons provided evidence for unexpected transmission routes, with subsequent mining of epidemiological data pointing to possible explanations for these transmissions. Our analysis demonstrates that integration of genomic and epidemiological data can yield actionable insights and facilitate the control of nosocomial transmission.

A selective jumonji H3K27 demethylase inhibitor modulates the proinflammatory macrophage response

The jumonji (JMJ) family of histone demethylases are Fe2+- and α-ketoglutarate-dependent oxygenases that are essential components of regulatory transcriptional chromatin complexes. These enzymes demethylate lysine residues in histones in a methylation-state and sequence-specific context. Considerable effort has been devoted to gaining a mechanistic understanding of the roles of histone lysine demethylases in eukaryotic transcription, genome integrity and epigenetic inheritance, as well as in development, physiology and disease. However, because of the absence of any selective inhibitors, the relevance of the demethylase activity of JMJ enzymes in regulating cellular responses remains poorly understood. Here we present a structure-guided small-molecule and chemoproteomics approach to elucidating the functional role of the H3K27me3-specific demethylase subfamily (KDM6 subfamily members JMJD3 and UTX). The liganded structures of human and mouse JMJD3 provide novel insight into the specificity determinants for cofactor, substrate and inhibitor recognition by the KDM6 subfamily of demethylases. We exploited these structural features to generate the first small-molecule catalytic site inhibitor that is selective for the H3K27me3-specific JMJ subfamily. We demonstrate that this inhibitor binds in a novel manner and reduces lipopolysaccharide-induced proinflammatory cytokine production by human primary macrophages, a process that depends on both JMJD3 and UTX. Our results resolve the ambiguity associated with the catalytic function of H3K27-specific JMJs in regulating disease-relevant inflammatory responses and provide encouragement for designing small-molecule inhibitors to allow selective pharmacological intervention across the JMJ family.

Transcription-associated mutational asymmetry in mammalian evolution

Although mutation is commonly thought of as a random process, evolutionary studies show that different types of nucleotide substitution occur with widely varying rates that presumably reflect biases intrinsic to mutation and repair mechanisms. A strand asymmetry, the occurrence of particular substitution types at higher rates than their complementary types, that is associated with DNA replication has been found in bacteria and mitochondria. A strand asymmetry that is associated with transcription and attributable to higher rates of cytosine deamination on the coding strand has been observed in enterobacteria. Here, we describe a qualitatively different transcription-associated strand asymmetry in mammals, which may be a byproduct of transcription-coupled repair in germline cells. This mutational asymmetry has acted over long periods of time to produce a compositional asymmetry, an excess of G+T over A+C on the coding strand, in most genes. The mutational and compositional asymmetries can be used to detect the orientations and approximate extents of transcribed regions.

Single-molecule sequencing to track plasmid diversity of hospital-associated carbapenemase-producing Enterobacteriaceae

Single-molecule sequencing of bacteria at the NIH Clinical Center documents diverse plasmids encoding antibiotic resistance and their transfer between microbes. How Antibiotic Resistance Spreads Among Bacteria Antibiotic-resistant microbes are spreading at an alarming rate in health care facilities throughout the world. Conlan et al. use a new DNA sequencing method to take a close look at one way in which antibiotic resistance spreads. With single-molecule sequencing, the authors completely characterized individual plasmids, the circular bits of DNA that carry the genes for antibiotic resistance in bacteria. They focused on resistance to the carbapenems, a class of antibiotics that is often used for infections that do not respond to more conventional antimicrobial agents. By using this approach in their microbial surveillance program at the NIH Clinical Center, the authors found evidence that plasmids carrying carbapenemase genes moved from one microbial species to another within the hospital environment. They also used the technique to test hypotheses about patient-to-patient transmission and to characterize a previously undescribed carbapenemase-encoding plasmid carried by diverse bacterial species that could cause dangerous clinical infections. Public health officials have raised concerns that plasmid transfer between Enterobacteriaceae species may spread resistance to carbapenems, an antibiotic class of last resort, thereby rendering common health care–associated infections nearly impossible to treat. To determine the diversity of carbapenemase-encoding plasmids and assess their mobility among bacterial species, we performed comprehensive surveillance and genomic sequencing of carbapenem-resistant Enterobacteriaceae in the National Institutes of Health (NIH) Clinical Center patient population and hospital environment. We isolated a repertoire of carbapenemase-encoding Enterobacteriaceae, including multiple strains of Klebsiella pneumoniae, Klebsiella oxytoca, Escherichia coli, Enterobacter cloacae, Citrobacter freundii, and Pantoea species. Long-read genome sequencing with full end-to-end assembly revealed that these organisms carry the carbapenem resistance genes on a wide array of plasmids. K. pneumoniae and E. cloacae isolated simultaneously from a single patient harbored two different carbapenemase-encoding plasmids, indicating that plasmid transfer between organisms was unlikely within this patient. We did, however, find evidence of horizontal transfer of carbapenemase-encoding plasmids between K. pneumoniae, E. cloacae, and C. freundii in the hospital environment. Our data, including full plasmid identification, challenge assumptions about horizontal gene transfer events within patients and identify possible connections between patients and the hospital environment. In addition, we identified a new carbapenemase-encoding plasmid of potentially high clinical impact carried by K. pneumoniae, E. coli, E. cloacae, and Pantoea species, in unrelated patients and in the hospital environment.

Staphylococcus epidermidis pan-genome sequence analysis reveals diversity of skin commensal and hospital infection-associated isolates

BackgroundWhile Staphylococcus epidermidis is commonly isolated from healthy human skin, it is also the most frequent cause of nosocomial infections on indwelling medical devices. Despite its importance, few genome sequences existed and the most frequent hospital-associated lineage, ST2, had not been fully sequenced.ResultsWe cultivated 71 commensal S. epidermidis isolates from 15 skin sites and compared them with 28 nosocomial isolates from venous catheters and blood cultures. We produced 21 commensal and 9 nosocomial draft genomes, and annotated and compared their gene content, phylogenetic relatedness and biochemical functions. The commensal strains had an open pan-genome with 80% core genes and 20% variable genes. The variable genome was characterized by an overabundance of transposable elements, transcription factors and transporters. Biochemical diversity, as assayed by antibiotic resistance and in vitro biofilm formation, demonstrated the varied phenotypic consequences of this genomic diversity. The nosocomial isolates exhibited both large-scale rearrangements and single-nucleotide variation. We showed that S. epidermidis genomes separate into two phylogenetic groups, one consisting only of commensals. The formate dehydrogenase gene, present only in commensals, is a discriminatory marker between the two groups.ConclusionsCommensal skin S. epidermidis have an open pan-genome and show considerable diversity between isolates, even when derived from a single individual or body site. For ST2, the most common nosocomial lineage, we detect variation between three independent isolates sequenced. Finally, phylogenetic analyses revealed a previously unrecognized group of S. epidermidis strains characterized by reduced virulence and formate dehydrogenase, which we propose as a clinical molecular marker.

Adaptive Evolution of Foundation Kinetochore Proteins in Primates

Rapid evolution is a hallmark of centromeric DNA in eukaryotic genomes. Yet, the centromere itself has a conserved functional role that is mediated by the kinetochore protein complex. To broaden our understanding about both the DNA and proteins that interact at the functional centromere, we sought to gain a detailed view of the evolutionary events that have shaped the primate kinetochore. Specifically, we performed comparative mapping and sequencing of the genomic regions encompassing the genes encoding three foundation kinetochore proteins: Centromere Proteins A, B, and C (CENP-A, CENP-B, and CENP-C). A histone H3 variant, CENP-A provides the foundation of the centromere-specific nucleosome. Comparative sequence analyses of the CENP-A gene in 14 primate species revealed encoded amino-acid residues within both the histone-fold domain and the N-terminal tail that are under strong positive selection. Similar comparative analyses of CENP-C, another foundation protein essential for centromere function, identified amino-acid residues throughout the protein under positive selection in the primate lineage, including several in the centromere localization and DNA-binding regions. Perhaps surprisingly, the gene encoding CENP-B, a kinetochore protein that binds specifically to alpha-satellite DNA, was not found to be associated with signatures of positive selection. These findings point to important and distinct evolutionary forces operating on the DNA and proteins of the primate centromere.

Plasmid Dynamics in KPC-Positive Klebsiella pneumoniae during Long-Term Patient Colonization

ABSTRACT Carbapenem-resistant Klebsiella pneumoniae strains are formidable hospital pathogens that pose a serious threat to patients around the globe due to a rising incidence in health care facilities, high mortality rates associated with infection, and potential to spread antibiotic resistance to other bacterial species, such as Escherichia coli. Over 6 months in 2011, 17 patients at the National Institutes of Health (NIH) Clinical Center became colonized with a highly virulent, transmissible carbapenem-resistant strain of K. pneumoniae. Our real-time genomic sequencing tracked patient-to-patient routes of transmission and informed epidemiologists’ actions to monitor and control this outbreak. Two of these patients remained colonized with carbapenemase-producing organisms for at least 2 to 4 years, providing the opportunity to undertake a focused genomic study of long-term colonization with antibiotic-resistant bacteria. Whole-genome sequencing studies shed light on the underlying complex microbial colonization, including mixed or evolving bacterial populations and gain or loss of plasmids. Isolates from NIH patient 15 showed complex plasmid rearrangements, leaving the chromosome and the blaKPC-carrying plasmid intact but rearranging the two other plasmids of this outbreak strain. NIH patient 16 has shown continuous colonization with blaKPC-positive organisms across multiple time points spanning 2011 to 2015. Genomic studies defined a complex pattern of succession and plasmid transmission across two different K. pneumoniae sequence types and an E. coli isolate. These findings demonstrate the utility of genomic methods for understanding strain succession, genome plasticity, and long-term carriage of antibiotic-resistant organisms. IMPORTANCE In 2011, the NIH Clinical Center had a nosocomial outbreak involving 19 patients who became colonized or infected with blaKPC-positive Klebsiella pneumoniae. Patients who have intestinal colonization with blaKPC-positive K. pneumoniae are at risk for developing infections that are difficult or nearly impossible to treat with existing antibiotic options. Two of those patients remained colonized with blaKPC-positive Klebsiella pneumoniae for over a year, leading to the initiation of a detailed genomic analysis exploring mixed colonization, plasmid recombination, and plasmid diversification. Whole-genome sequence analysis identified a variety of changes, both subtle and large, in the blaKPC-positive organisms. Long-term colonization of patients with blaKPC-positive Klebsiella pneumoniae creates new opportunities for horizontal gene transfer of plasmids encoding antibiotic resistance genes and poses complications for the delivery of health care. In 2011, the NIH Clinical Center had a nosocomial outbreak involving 19 patients who became colonized or infected with blaKPC-positive Klebsiella pneumoniae. Patients who have intestinal colonization with blaKPC-positive K. pneumoniae are at risk for developing infections that are difficult or nearly impossible to treat with existing antibiotic options. Two of those patients remained colonized with blaKPC-positive Klebsiella pneumoniae for over a year, leading to the initiation of a detailed genomic analysis exploring mixed colonization, plasmid recombination, and plasmid diversification. Whole-genome sequence analysis identified a variety of changes, both subtle and large, in the blaKPC-positive organisms. Long-term colonization of patients with blaKPC-positive Klebsiella pneumoniae creates new opportunities for horizontal gene transfer of plasmids encoding antibiotic resistance genes and poses complications for the delivery of health care.

Evolution of a bitter taste receptor gene cluster in a New World sparrow.

Bitter taste perception likely evolved as a protective mechanism against the ingestion of harmful compounds in food. The evolution of the taste receptor type 2 (TAS2R) gene family, which encodes the chemoreceptors that are directly responsible for the detection of bitter compounds, has therefore been of considerable interest. Though TAS2R repertoires have been characterized for a number of species, to date the complement of TAS2Rs from just one bird, the chicken, which had a notably small number of TAS2Rs, has been established. Here, we used targeted mapping and genomic sequencing in the white-throated sparrow (Zonotrichia albicollis) and sample sequencing in other closely related birds to reconstruct the history of a TAS2R gene cluster physically linked to the break points of an evolutionary chromosomal rearrangement. In the white-throated sparrow, this TAS2R cluster encodes up to 18 functional bitter taste receptors and likely underwent a large expansion that predates and/or coincides with the radiation of the Emberizinae subfamily into the New World. In addition to signatures of gene birth-and-death evolution within this cluster, estimates of Ka/Ks for the songbird TAS2Rs were similar to those previously observed in mammals, including humans. Finally, comparison of the complete genomic sequence of the cluster from two common haplotypes in the white-throated sparrow revealed a number of nonsynonymous variants and differences in functional gene content within this species. These results suggest that interspecies and intraspecies genetic variability does exist in avian TAS2Rs and that these differences could contribute to variation in bitter taste perception in birds.

Haplotype-Based Genomic Sequencing of a Chromosomal Polymorphism in the White-Throated Sparrow (Zonotrichia albicollis)

Inversion polymorphisms have been linked to a variety of fundamental biological and evolutionary processes. Yet few studies have used large-scale genomic sequencing to directly compare the haplotypes associated with the standard and inverted chromosome arrangements. Here we describe the targeted genomic sequencing and comparison of haplotypes representing alternative arrangements of a common inversion polymorphism linked to a suite of phenotypes in the white-throated sparrow (Zonotrichia albicollis). More than 7.4 Mb of genomic sequence was generated and assembled from both the standard (ZAL2) and inverted (ZAL2(m)) arrangements. Sequencing of a pair of inversion breakpoints led to the identification of a ZAL2-specific segmental duplication, as well as evidence of breakpoint reusage. Comparison of the haplotype-based sequence assemblies revealed low genetic differentiation outside versus inside the inversion indicative of historical patterns of gene flow and suppressed recombination between ZAL2 and ZAL2(m). Finally, despite ZAL2(m) being maintained in a near constant state of heterozygosity, no signatures of genetic degeneration were detected on this chromosome. Overall, these results provide important insights into the genomic attributes of an inversion polymorphism linked to mate choice and variation in social behavior.