Christina Sison

Chromatin stretch enhancer states drive cell-specific gene regulation and harbor human disease risk variants

Significance Using high-throughput experiments, we determined the functional epigenomic landscape in pancreatic islet cells. Computational integration of these data along with similar data from the ENCODE project revealed the presence of large gene control elements across diverse cell types that we refer to as “stretch enhancers.” Stretch enhancers are cell type specific and are associated with increased expression of genes involved in cell-specific processes. We find that genetic variations associated with common disease are highly enriched in stretch enhancers; notably, stretch enhancers specific to pancreatic islets harbor variants linked to type 2 diabetes and related traits. We propose that stretch enhancers form as pluripotent cells differentiate into committed lineages, to program important cell-specific gene expression. Chromatin-based functional genomic analyses and genomewide association studies (GWASs) together implicate enhancers as critical elements influencing gene expression and risk for common diseases. Here, we performed systematic chromatin and transcriptome profiling in human pancreatic islets. Integrated analysis of islet data with those from nine cell types identified specific and significant enrichment of type 2 diabetes and related quantitative trait GWAS variants in islet enhancers. Our integrated chromatin maps reveal that most enhancers are short (median = 0.8 kb). Each cell type also contains a substantial number of more extended (≥3 kb) enhancers. Interestingly, these stretch enhancers are often tissue-specific and overlap locus control regions, suggesting that they are important chromatin regulatory beacons. Indeed, we show that (i) tissue specificity of enhancers and nearby gene expression increase with enhancer length; (ii) neighborhoods containing stretch enhancers are enriched for important cell type–specific genes; and (iii) GWAS variants associated with traits relevant to a particular cell type are more enriched in stretch enhancers compared with short enhancers. Reporter constructs containing stretch enhancer sequences exhibited tissue-specific activity in cell culture experiments and in transgenic mice. These results suggest that stretch enhancers are critical chromatin elements for coordinating cell type–specific regulatory programs and that sequence variation in stretch enhancers affects risk of major common human diseases.

Maturation and Diversity of the VRC01-Antibody Lineage over 15 Years of Chronic HIV-1 Infection

HIV-1-neutralizing antibodies develop in most HIV-1-infected individuals, although highly effective antibodies are generally observed only after years of chronic infection. Here, we characterize the rate of maturation and extent of diversity for the lineage that produced the broadly neutralizing antibody VRC01 through longitudinal sampling of peripheral B cell transcripts over 15 years and co-crystal structures of lineage members. Next-generation sequencing identified VRC01-lineage transcripts, which encompassed diverse antibodies organized into distinct phylogenetic clades. Prevalent clades maintained characteristic features of antigen recognition, though each evolved binding loops and disulfides that formed distinct recognition surfaces. Over the course of the study period, VRC01-lineage clades showed continuous evolution, with rates of ∼2 substitutions per 100 nucleotides per year, comparable to that of HIV-1 evolution. This high rate of antibody evolution provides a mechanism by which antibody lineages can achieve extraordinary diversity and, over years of chronic infection, develop effective HIV-1 neutralization.

Mining the antibodyome for HIV-1–neutralizing antibodies with next-generation sequencing and phylogenetic pairing of heavy/light chains

Next-generation sequencing of antibody transcripts from HIV-1–infected individuals with broadly neutralizing antibodies could provide an efficient means for identifying somatic variants and characterizing their lineages. Here, we used 454 pyrosequencing and identity/divergence grid sampling to analyze heavy- and light-chain sequences from donor N152, the source of the broadly neutralizing antibody 10E8. We identified variants with up to 28% difference in amino acid sequence. Heavy- and light-chain phylogenetic trees of identified 10E8 variants displayed similar architectures, and 10E8 variants reconstituted from matched and unmatched phylogenetic branches displayed significantly lower autoreactivity when matched. To test the generality of phylogenetic pairing, we analyzed donor International AIDS Vaccine Initiative 84, the source of antibodies PGT141–145. Heavy- and light-chain phylogenetic trees of PGT141–145 somatic variants also displayed remarkably similar architectures; in this case, branch pairings could be anchored by known PGT141–145 antibodies. Altogether, our findings suggest that phylogenetic matching of heavy and light chains can provide a means to approximate natural pairings.

Genome-wide recombination drives diversification of epidemic strains of Acinetobacter baumannii

Acinetobacter baumannii is an emerging human pathogen and a significant cause of nosocomial infections among hospital patients worldwide. The enormous increase in multidrug resistance among hospital isolates and the recent emergence of pan-drug–resistant strains underscores the urgency to understand how A. baumannii evolves in hospital environments. To this end, we undertook a genomic study of a polyclonal outbreak of multidrug-resistant A. baumannii at the research-based National Institutes of Health Clinical Center. Comparing the complete genome sequences of the three dominant outbreak strain types enabled us to conclude that, despite all belonging to the same epidemic lineage, the three strains diverged before their arrival at the National Institutes of Health. The simultaneous presence of three divergent strains from this lineage supports its increasing prevalence in international hospitals and suggests an ongoing adaptation to the hospital environment. Further genomic comparisons uncovered that much of the diversification that occurred since the divergence of the three outbreak strains was mediated by homologous recombination across 20% of their genomes. Inspection of recombinant regions revealed that several regions were associated with either the loss or swapping out of genes encoding proteins that are exposed to the cell surface or that synthesize cell-surface molecules. Extending our analysis to a larger set of international clinical isolates revealed a previously unappreciated ability of A. baumannii to vary surface molecules through horizontal gene transfer, with subsequent intraspecies dissemination by homologous recombination. These findings have immediate implications in surveillance, prevention, and treatment of A. baumannii infections.

Circadian changes in long noncoding RNAs in the pineal gland.

Long noncoding RNAs (lncRNAs) play a broad range of biological roles, including regulation of expression of genes and chromosomes. Here, we present evidence that lncRNAs are involved in vertebrate circadian biology. Differential night/day expression of 112 lncRNAs (0.3 to >50 kb) occurs in the rat pineal gland, which is the source of melatonin, the hormone of the night. Approximately one-half of these changes reflect nocturnal increases. Studies of eight lncRNAs with 2- to >100-fold daily rhythms indicate that, in most cases, the change results from neural stimulation from the central circadian oscillator in the suprachiasmatic nucleus (doubling time = 0.5–1.3 h). Light exposure at night rapidly reverses (halving time = 9–32 min) levels of some of these lncRNAs. Organ culture studies indicate that expression of these lncRNAs is regulated by norepinephrine acting through cAMP. These findings point to a dynamic role of lncRNAs in the circadian system.

Analysis of immunoglobulin transcripts and hypermutation following SHIV(AD8) infection and protein-plus-adjuvant immunization.

Developing predictive animal models to assess how candidate vaccines and infection influence the ontogenies of Envelope (Env)-specific antibodies is critical for the development of an HIV vaccine. Here we use two nonhuman primate models to compare the roles of antigen persistence, diversity and innate immunity. We perform longitudinal analyses of HIV Env-specific B-cell receptor responses to SHIVAD8 infection and Env protein vaccination with eight different adjuvants. A subset of the SHIVAD8-infected animals with higher viral loads and greater Env diversity show increased neutralization associated with increasing somatic hypermutation (SHM) levels over time. The use of adjuvants results in increased ELISA titres but does not affect the mean SHM levels or CDR H3 lengths. Our study shows how the ontogeny of Env-specific B cells can be tracked, and provides insights into the requirements for developing neutralizing antibodies that should facilitate translation to human vaccine studies.

Addressing Bias in Small RNA Library Preparation for Sequencing: A New Protocol Recovers MicroRNAs that Evade Capture by Current Methods

Recent advances in sequencing technology have helped unveil the unexpected complexity and diversity of small RNAs. A critical step in small RNA library preparation for sequencing is the ligation of adapter sequences to both the 5′ and 3′ ends of small RNAs. Studies have shown that adapter ligation introduces a significant but widely unappreciated bias in the results of high-throughput small RNA sequencing. We show that due to this bias the two widely used Illumina library preparation protocols produce strikingly different microRNA (miRNA) expression profiles in the same batch of cells. There are 102 highly expressed miRNAs that are >5-fold differentially detected and some miRNAs, such as miR-24-3p, are over 30-fold differentially detected. While some level of bias in library preparation is not surprising, the apparent massive differential bias between these two widely used adapter sets is not well appreciated. In an attempt to mitigate this bias, the new Bioo Scientific NEXTflex V2 protocol utilizes a pool of adapters with random nucleotides at the ligation boundary. We show that this protocol is able to detect robustly several miRNAs that evade capture by the Illumina-based methods. While these analyses do not indicate a definitive gold standard for small RNA library preparation, the results of the NEXTflex protocol do correlate best with RT-qPCR. As increasingly more laboratories seek to study small RNAs, researchers should be aware of the extent to which the results may differ with different protocols, and should make an informed decision about the protocol that best fits their study.

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.

Gene-Specific Substitution Profiles Describe the Types and Frequencies of Amino Acid Changes during Antibody Somatic Hypermutation

Somatic hypermutation (SHM) plays a critical role in the maturation of antibodies, optimizing recognition initiated by recombination of V(D)J genes. Previous studies have shown that the propensity to mutate is modulated by the context of surrounding nucleotides and that SHM machinery generates biased substitutions. To investigate the intrinsic mutation frequency and substitution bias of SHMs at the amino acid level, we analyzed functional human antibody repertoires and developed mGSSP (method for gene-specific substitution profile), a method to construct amino acid substitution profiles from next-generation sequencing-determined B cell transcripts. We demonstrated that these gene-specific substitution profiles (GSSPs) are unique to each V gene and highly consistent between donors. We also showed that the GSSPs constructed from functional antibody repertoires are highly similar to those constructed from antibody sequences amplified from non-productively rearranged passenger alleles, which do not undergo functional selection. This suggests the types and frequencies, or mutational space, of a majority of amino acid changes sampled by the SHM machinery to be well captured by GSSPs. We further observed the rates of mutational exchange between some amino acids to be both asymmetric and context dependent and to correlate weakly with their biochemical properties. GSSPs provide an improved, position-dependent alternative to standard substitution matrices, and can be utilized to developing software for accurately modeling the SHM process. GSSPs can also be used for predicting the amino acid mutational space available for antigen-driven selection and for understanding factors modulating the maturation pathways of antibody lineages in a gene-specific context. The mGSSP method can be used to build, compare, and plot GSSPs1; we report the GSSPs constructed for 69 common human V genes (DOI: 10.6084/m9.figshare.3511083) and provide high-resolution logo plots for each (DOI: 10.6084/m9.figshare.3511085).

A Neutralizing Antibody Recognizing Primarily N-Linked Glycan Targets the Silent Face of the HIV Envelope

&NA; Virtually the entire surface of the HIV‐1‐envelope trimer is recognized by neutralizing antibodies, except for a highly glycosylated region at the center of the “silent face” on the gp120 subunit. From an HIV‐1‐infected donor, #74, we identified antibody VRC‐PG05, which neutralized 27% of HIV‐1 strains. The crystal structure of the antigen‐binding fragment of VRC‐PG05 in complex with gp120 revealed an epitope comprised primarily of N‐linked glycans from N262, N295, and N448 at the silent face center. Somatic hypermutation occurred preferentially at antibody residues that interacted with these glycans, suggesting somatic development of glycan recognition. Resistance to VRC‐PG05 in donor #74 involved shifting of glycan‐N448 to N446 or mutation of glycan‐proximal residue E293. HIV‐1 neutralization can thus be achieved at the silent face center by glycan‐recognizing antibody; along with other known epitopes, the VRC‐PG05 epitope completes coverage by neutralizing antibody of all major exposed regions of the prefusion closed trimer. Graphical Abstract Figure. No caption available. HighlightsIdentified and defined crystal structure of antibody VRC‐PG05 in complex with gp120VRC‐PG05 epitope is at the center of the glycosylated silent face of HIV‐1 gp120VRC‐PG05 utilizes both glycopeptide and glycan‐cluster mechanisms of recognitionVRC‐PG05 completes neutralizing antibody coverage of the prefusion‐closed Env trimer &NA; The center of the “silent face” on the HIV‐1 envelope is shielded by glycans and has been devoid of antibody recognition. Zhou et al. identify the antibody VRC‐PG05, which binds a glycan‐dominated epitope at the silent face center and completes antibody recognition of all major exposed regions of the envelope trimer.