Search

Clinical applications of precision oncology require accurate tests that can distinguish true cancer-specific mutations from errors introduced at each step of next-generation sequencing (NGS). To date, no bulk sequencing study has addressed the effects of cross-site reproducibility, nor the biological, technical and computational factors that influence variant identification. Here we report a systematic interrogation of somatic mutations in paired tumor–normal cell lines to identify factors affecting detection reproducibility and accuracy at six different centers. Using whole-genome sequencing (WGS) and whole-exome sequencing (WES), we evaluated the reproducibility of different sample types with varying input amount and tumor purity, and multiple library construction protocols, followed by processing with nine bioinformatics pipelines. We found that read coverage and callers affected both WGS and WES reproducibility, but WES performance was influenced by insert fragment size, genomic copy content and the global imbalance score (GIV; G > T/C > A). Finally, taking into account library preparation protocol, tumor content, read coverage and bioinformatics processes concomitantly, we recommend actionable practices to improve the reproducibility and accuracy of NGS experiments for cancer mutation detection.

Viruses with large double-stranded DNA genomes appear to have captured the majority of their genes from the hosts at different stages of evolution. The origin of many virus genes is readily detected through highly significant sequence similarity with cellular homologs. This is the case, in particular, for virus enzymes, such as DNA and RNA polymerases or nucleotide kinases, that retain their catalytic activity after capture by an ancestral virus. However, a large fraction of virus genes have no readily detectable cellular homologs so that their origin remains enigmatic. We sought to explore potential origins of proteins of unknown provenance encoded in the genomes of orthopoxviruses, a thoroughly studied virus genus which includes major human pathogens. To this end, we used AlphaFold2, to predict the structures of all 214 proteins encoded by orthopoxviruses. Among the proteins of unknown provenance, structure prediction yielded a clear indication of origin for 14, along with validating several inferences previously made by sequence analysis. The major trend that emerges from these findings is the exaptation of enzymes from cellular organisms for non-enzymatic, structural roles in virus reproduction which is accompanied by disruption of catalytic sites and overall drastic divergence which precludes detection of homology at the sequence level. Among the 16 orthopoxvirus proteins found to be inactivated enzyme derivatives, are the poxvirus replication processivity factor A20, an inactivated derivative of bacterial NAD-dependent DNA ligase; major core protein A3, an inactivated deubiquitinase; F11, an inactivated prolyl hydroxylase; and more similar cases. However, for nearly one third of the orthopoxvirus virion proteins, no significantly similar structures were identified, suggesting exaptation with subsequent major structural rearrangement, yielding novel protein folds.

Clinical applications of precision oncology require accurate tests that can distinguish cancer-specific mutations from errors introduced at each step of next generation sequencing (NGS). For NGS to successfully improve patient lives, discriminating between true mutations and artifacts is crucial. To date, no study has addressed the effects of cross site reproducibility together with the potentially influential interactions between biological, technical, and computational factors on the accurate identification of variants. Here we systematically interrogated somatic mutations in paired tumor-normal cell lines to identify factors affecting detection reproducibility and accuracy. Different types of samples with varying input amount and tumor purity were processed using multiple library construction protocols. Whole-genome (WGS) and whole-exome sequencing (WES) were carried out at six sequencing centers followed by processing with nine bioinformatics pipelines to evaluate reproducibility. We identified artifacts of C>A mutations in WES due to sample and library processing and highlighted limitations of bioinformatics tools for artifact detection and removal. By examining the interactions and effects of NGS platforms, library preparation protocol, tumor content, read coverage, and bioinformatics processes concomitantly, here we recommend actionable practices for cancer mutation detection experiments and analysis using NGS technologies.

Fixing Broken DNA Some physiological processes, such as immunoglobulin class switching and telomere attrition, result in double-stranded DNA breaks. The DNA damage repair protein, 53BP1, prevents nucleolytic processing of these breaks, but the proteins it partners with to do this are unknown (see the Perspective by Lukas and Luka s ). Di Virgilio et al. (p. 711 , published online 10 January), using mass spectroscopy–based methods, and Zimmermann et al. (p. 700 , published online 10 January), using a telomere-based assay, identify Rif1 as a 53BP1 phosphorylation- and DNA damage–dependent interaction partner. Mice with a B cell–specific deletion in Rif1 showed impaired immunoglobulin class switching. Rif1-deficient cells exhibited extensive 5′-3′ resection at DNA ends, with enhanced genetic instability. Thus, Rif1 partners with 53BP1 to promote the proper repair of double-stranded DNA breaks.

Deficiencies in factors that regulate the DNA damage response enhance the incidence of malignancy by destabilizing the genome. However, the precise influence of the DNA damage response on regulation of cancer-associated rearrangements is not well defined. Here we examine the genome-wide impact of tumor protein P53-binding protein 1 (53BP1) deficiency in lymphoma and translocation. While both activation-induced cytidine deaminase (AID) and 53BP1 have been associated with cancer in humans, neither AID overexpression nor loss of 53BP1 is sufficient to produce malignancy. However, the combination of 53BP1 deficiency and AID deregulation results in B cell lymphoma. Deep sequencing of the genome of 53BP1(-/-) cancer cells and translocation capture sequencing (TC-Seq) of primary 53BP1(-/-) B cells revealed that their chromosomal rearrangements differ from those found in wild-type cells in that they show increased DNA end resection. Moreover, loss of 53BP1 alters the translocatome by increasing rearrangements to intergenic regions.

SummaryActivation-induced cytidine deaminase (AID) promotes chromosomal translocations by inducing DNA double-strand breaks (DSBs) at immunoglobulin (Ig) genes and oncogenes in the G1 phase. RPA is a single-stranded DNA (ssDNA)-binding protein that associates with resected DSBs in the S phase and facilitates the assembly of factors involved in homologous repair (HR), such as Rad51. Notably, RPA deposition also marks sites of AID-mediated damage, but its role in Ig gene recombination remains unclear. Here, we demonstrate that RPA associates asymmetrically with resected ssDNA in response to lesions created by AID, recombination-activating genes (RAG), or other nucleases. Small amounts of RPA are deposited at AID targets in G1 in an ATM-dependent manner. In contrast, recruitment in the S–G2/M phase is extensive, ATM independent, and associated with Rad51 accumulation. In the S–G2/M phase, RPA increases in nonhomologous-end-joining-deficient lymphocytes, where there is more extensive DNA-end resection. Thus, most RPA recruitment during class switch recombination represents salvage of unrepaired breaks by homology-based pathways during the S–G2/M phase of the cell cycle.

DNA in cells is predominantly B-form double helix. Though certain DNA sequences in vitro may fold into other structures, such as triplex, left-handed Z form, or quadruplex DNA, the stability and prevalence of these structures in vivo are not known. Here, using computational analysis of sequence motifs, RNA polymerase II binding data, and genome-wide potassium permanganate-dependent nuclease footprinting data, we map thousands of putative non-B DNA sites at high resolution in mouse B cells. Computational analysis associates these non-B DNAs with particular structures and indicates that they form at locations compatible with an involvement in gene regulation. Further analyses support the notion that non-B DNA structure formation influences the occupancy and positioning of nucleosomes in chromatin. These results suggest that non-B DNAs contribute to the control of a variety of critical cellular and organismal processes.

SummaryThe c-Myc HLH-bZIP protein has been implicated in physiological or pathological growth, proliferation, apoptosis, metabolism, and differentiation at the cellular, tissue, or organismal levels via regulation of numerous target genes. No principle yet unifies Myc action due partly to an incomplete inventory and functional accounting of Myc’s targets. To observe Myc target expression and function in a system where Myc is temporally and physiologically regulated, the transcriptomes and the genome-wide distributions of Myc, RNA polymerase II, and chromatin modifications were compared during lymphocyte activation and in ES cells as well. A remarkably simple rule emerged from this quantitative analysis: Myc is not an on-off specifier of gene activity, but is a nonlinear amplifier of expression, acting universally at active genes, except for immediate early genes that are strongly induced before Myc. This rule of Myc action explains the vast majority of Myc biology observed in literature.

Recurrent chromosomal translocations underlie both haematopoietic and solid tumours. Their origin has been ascribed to selection of random rearrangements, targeted DNA damage, or frequent nuclear interactions between translocation partners; however, the relative contribution of each of these elements has not been measured directly or on a large scale. Here we examine the role of nuclear architecture and frequency of DNA damage in the genesis of chromosomal translocations by measuring these parameters simultaneously in cultured mouse B lymphocytes. In the absence of recurrent DNA damage, translocations between Igh or Myc and all other genes are directly related to their contact frequency. Conversely, translocations associated with recurrent site-directed DNA damage are proportional to the rate of DNA break formation, as measured by replication protein A accumulation at the site of damage. Thus, non-targeted rearrangements reflect nuclear organization whereas DNA break formation governs the location and frequency of recurrent translocations, including those driving B-cell malignancies.

Chromosomal translocations require formation and joining of DNA double strand breaks (DSBs). These events disrupt the integrity of the genome and are involved in producing leukemias, lymphomas and sarcomas. Translocations are frequent, clonal and recurrent in mature B cell lymphomas, which bear a particularly high DNA damage burden by virtue of Activation Induced Cytidine Deaminase (AID) expression. Despite the ubiquity of genomic rearrangements, the forces that underlie their genesis are not well understood. Here, we provide a detailed description of a new method for studying these events, Translocation Capture Sequencing (TC-Seq). TC-Seq provides the means to document chromosomal rearrangements genome-wide in primary cells, and to discover recombination hotspots. Demonstrating its effectiveness, we successfully estimate the frequency of c-myc/IgH translocations in primary B cells, and identify hotspots of AID-mediated recombination. Furthermore, TC-Seq can be adapted to generate genome-wide rearrangement maps in any cell type and under any condition.

The vaccinia virus E6R gene (VACVWR062) is conserved in all members of the poxvirus family and encodes a protein associated with the mature virion. We confirmed this association and provided evidence for an internal location. An inducible mutant that conditionally expresses E6 was constructed. In the absence of inducer, plaque formation and virus production were severely inhibited in several cell lines, whereas some replication occurred in others. This difference could be due to variation in the stringency of repression, since we could not isolate a stable deletion mutant even in the more “permissive” cells. Under non-permissive conditions, viral late proteins were synthesized but processing of core proteins was inefficient, indicative of an assembly block. Transmission electron microscopy of sections of cells infected with the mutant in the absence of inducer revealed morphogenetic defects with crescents and empty immature virions adjacent to dense inclusions of viroplasm. Mature virions were infrequent and cores appeared to have lucent centers.

Rapid increases in DNA sequencing capabilities have led to a vast increase in the data generated from prokaryotic genomic studies, which has been a boon to scientists studying micro-organism evolution and to those who wish to understand the biological underpinnings of microbial systems. The NCBI Protein Clusters Database (ProtClustDB) has been created to efficiently maintain and keep the deluge of data up to date. ProtClustDB contains both curated and uncurated clusters of proteins grouped by sequence similarity. The May 2008 release contains a total of 285 386 clusters derived from over 1.7 million proteins encoded by 3806 nt sequences from the RefSeq collection of complete chromosomes and plasmids from four major groups: prokaryotes, bacteriophages and the mitochondrial and chloroplast organelles. There are 7180 clusters containing 376 513 proteins with curated gene and protein functional annotation. PubMed identifiers and external cross references are collected for all clusters and provide additional information resources. A suite of web tools is available to explore more detailed information, such as multiple alignments, phylogenetic trees and genomic neighborhoods. ProtClustDB provides an efficient method to aggregate gene and protein annotation for researchers and is available at http://www.ncbi.nlm.nih.gov/sites/entrez?db=proteinclusters.

Previous studies indicated that the vaccinia virus D10 protein, which is conserved in all sequenced poxviruses, participates in the rapid turnover of host and viral mRNAs. D10 contains a motif present in the family of Nudix/MutT enzymes, a subset of which has been shown to enhance mRNA turnover in eukaryotic cells through cleavage of the 5′ cap (m 7 GpppNm-). Here, we demonstrate that a purified recombinant D10 fusion protein possesses an intrinsic activity that liberates m 7 GDP from capped RNA substrates. Furthermore, point mutations in the Nudix/MutT motif abolished decapping activity. D10 has a strong affinity for capped RNA substrates ( K m ≈ 3 nm). RNAs of 24–309 nt were decapped to comparable extents, whereas the cap of a 12-nt RNA was uncleaved. At large molar ratios relative to capped RNA substrate, competitor m 7 GpppG, m 7 GTP, or m 7 GDP inhibited decapping, whereas even higher concentrations of unmethylated analogs did not. High concentrations of uncapped RNA were also inhibitory, suggesting that D10 recognizes its substrate through interaction with both cap and RNA moieties. Thus far, poxviruses represent the only virus family shown to encode a Nudix hydrolase-decapping enzyme. Although it may seem self-destructive for a virus to encode a decapping and a capping enzyme, accelerated mRNA turnover helps eliminate competing host mRNAs and allows stage-specific synthesis of viral proteins.

We provide the initial characterization of the product of the vaccinia virus L3L open reading frame (VACWR090), which is conserved in all sequenced members of the poxvirus family. The predicted polypeptide contains no motifs or other features that provided a clue to the role of the L3 protein, and no functional information was available regarding a homolog discovered in Plasmodium falciparum . The L3 protein was expressed following viral DNA replication, a finding consistent with a putative late promoter sequence, and was packaged as a non-membrane protein in mature virus particles. A recombinant virus, in which the L3L gene was regulated by the Escherichia coli lac operator/repressor system, had a conditional lethal phenotype. The virus replicated in the presence of inducer, but in its absence, the yields of infectious virus were reduced by 99%. When cells were infected without inducer, however, no defect in gene expression or morphogenesis was noted. Virus particles lacking L3, which assembled in the absence of inducer, were indistinguishable from wild-type virions with regard to morphology, major structural proteins, and DNA content but were noninfectious. L3-deficient virions were able to bind and penetrate cells but produced extremely small amounts of viral early mRNA. A defect in transcription was demonstrated by in vitro studies with permeabilized virions, but soluble extracts of L3-deficient virions showed normal levels of template-dependent transcriptional activity, indicating that only transcription of the packaged genome is impaired.

We examined the population dynamics of human immunodeficiency virus type 1 pro variants during the evolution of resistance to the protease inhibitor ritonavir (RTV) in vivo. pro variants were followed in subjects who had added RTV to their previously failed reverse transcriptase inhibitor therapy using a heteroduplex tracking assay designed to detect common resistance-associated mutations. In most cases the initial variant appeared rapidly within 2 to 3 months followed by one or more subsequent population turnovers. Some of the subsequent transitions between variants were rapid, and some were prolonged with the coexistence of multiple variants. In several cases variants without resistance mutations persisted despite the emergence of new variants with an increasing number of resistance-associated mutations. Based on the rate of turnover of pro variants in the RTV-treated subjects we estimated that the mean fitness of newly emerging variants was increased 1.2-fold (range, 1.02 to 1.8) relative to their predecessors. A subset of pro genes was introduced into infectious molecular clones. The corresponding viruses displayed impaired replication capacity and reduced susceptibility to RTV. A subset of these clones also showed increased susceptibility to two nonnucleoside reverse transcriptase inhibitors and the protease inhibitor saquinavir. Finally, a significant correlation between the reduced replication capacity and reduced processing at the gag NC-p1 processing site was noted. Our results reveal a complexity of patterns in the evolution of resistance to a protease inhibitor. In addition, these results suggest that selection for resistance to one protease inhibitor can have pleiotropic effects that can affect fitness and susceptibility to other drugs.

The initiation of drug therapy or the addition of a new drug to preexisting therapy can have a significant impact on human immunodeficiency virus type 1 (HIV-1) populations within a person. Drug therapy directed at reverse transcriptase and protease can result in dramatic decreases in virus load, causing a contraction in the virus population that represents a potential genetic bottleneck as a subset of virus with genomes carrying resistance mutations repopulate the host. While this bottleneck exerts an effect directly on the region that is being targeted by the drugs, it also affects other regions of the viral genome. We have applied heteroduplex tracking assays (HTA) specific to variable regions 1 and 2 (V1/V2) and variable region 3 (V3) of the HIV-1 env gene to analyze the effect of a genetic bottleneck created by the selection of resistance to ritonavir, a protease inhibitor. Subjects were classified into groups on the basis of the extent of the initial drop in virus load and the duration of virus load reduction. Subjects with a strong initial drop in virus load exhibited a loss of heterogeneity in the env region at virus load rebound; in contrast, subjects with a weak initial drop in virus load exhibited little to no loss of heterogeneity at virus load rebound in either region of env examined. The duration of virus load reduction also affected env populations. Subjects that had prolonged reductions exhibited slower population diversification and the appearance of new V1/V2 species after rebound. The longer reduction of virus load in these subjects may have allowed for improved immune system function, which in turn could have selected for new escape mutants. Subjects with rapid rebound quickly reequilibrated the entry env variants back into the resistant population. When the pro gene developed further resistance mutations subsequent to virus load rebound, no changes were observed in V1/V2 or V3 populations, suggesting that the high virus loads allowed the env populations to reequilibrate rapidly. The rapid equilibration of env variants during pro gene sequence transitions at high virus load suggests that recombination is active in defining the HIV-1 sequence population. Conversely, part of the success of suppressive antiviral therapy may be to limit the potential for evolution through recombination, which requires dually infected cells.

The vaccinia virus A11R gene has orthologs in all known poxvirus genomes, and the A11 protein has been previously reported to interact with the putative DNA packaging protein A32 in a yeast two-hybrid screen. Using antisera raised against A11 peptides, we show that the A11 protein was (i) expressed at late times with an apparent mass of 40 kDa, (ii) not incorporated into virus particles, (iii) phosphorylated independently of the viral F10 kinase, (iv) coimmunoprecipitated with A32, and (v) localized to the viral factory. To determine the role of the A11 protein and test whether it is indeed involved in DNA packaging, we constructed a recombinant vaccinia virus with an inducible A11R gene. This recombinant was dependent on inducer for single-cycle growth and plaque formation. In the absence of inducer, viral late proteins were produced at normal levels, but proteolytic processing and other posttranslational modifications of some proteins were inhibited, suggesting a block in virus particle assembly. Consistent with this observation, electron microscopy of cells infected in the absence of inducer showed virus factories with abnormal electron-dense viroplasms and intermediate density regions associated with membranes and containing the D13 protein. However, no viral membrane crescents, immature virions, or mature virions were produced. The requirement for nonvirion protein A11 in order to make normal viral membranes was an unexpected and exciting finding, since neither the origin of these membranes nor their mechanism of formation in the cytoplasm of infected cells is understood.

50 years ago, Vincent Allfrey and colleagues discovered that lymphocyte activation triggers massive acetylation of chromatin. However, the molecular mechanisms driving epigenetic accessibility are still unknown. We here show that stimulated lymphocytes decondense chromatin by three differentially regulated steps. First, chromatin is repositioned away from the nuclear periphery in response to global acetylation. Second, histone nanodomain clusters decompact into mononucleosome fibers through a mechanism that requires Myc and continual energy input. Single-molecule imaging shows that this step lowers transcription factor residence time and non-specific collisions during sampling for DNA targets. Third, chromatin interactions shift from long range to predominantly short range, and CTCF-mediated loops and contact domains double in numbers. This architectural change facilitates cognate promoter-enhancer contacts and also requires Myc and continual ATP production. Our results thus define the nature and transcriptional impact of chromatin decondensation and reveal an unexpected role for Myc in the establishment of nuclear topology in mammalian cells.

Mutations of STAT3 underlie the autosomal dominant form of hyperimmunoglobulin E syndrome (HIES). STAT3 has critical roles in immune cells and thus, hematopoietic stem cell transplantation (HSCT), might be a reasonable therapeutic strategy in this disease. However, STAT3 also has critical functions in nonhematopoietic cells and dissecting the protean roles of STAT3 is limited by the lethality associated with germline deletion of Stat3. Thus, predicting the efficacy of HSCT for HIES is difficult. To begin to dissect the importance of STAT3 in hematopoietic and nonhematopoietic cells as it relates to HIES, we generated a mouse model of this disease. We found that these transgenic mice recapitulate multiple aspects of HIES, including elevated serum IgE and failure to generate Th17 cells. We found that these mice were susceptible to bacterial infection that was partially corrected by HSCT using wild-type bone marrow, emphasizing the role played by the epithelium in the pathophysiology of HIES.

SummaryA key finding of the ENCODE project is that the enhancer landscape of mammalian cells undergoes marked alterations during ontogeny. However, the nature and extent of these changes are unclear. As part of the NIH Mouse Regulome Project, we here combined DNaseI hypersensitivity, ChIP-seq, and ChIA-PET technologies to map the promoter-enhancer interactomes of pluripotent ES cells and differentiated B lymphocytes. We confirm that enhancer usage varies widely across tissues. Unexpectedly, we find that this feature extends to broadly transcribed genes, including Myc and Pim1 cell-cycle regulators, which associate with an entirely different set of enhancers in ES and B cells. By means of high-resolution CpG methylomes, genome editing, and digital footprinting, we show that these enhancers recruit lineage-determining factors. Furthermore, we demonstrate that the turning on and off of enhancers during development correlates with promoter activity. We propose that organisms rely on a dynamic enhancer landscape to control basic cellular functions in a tissue-specific manner.

Rapidly evolving entities, such as viruses, can undergo complex genetic changes in the face of strong selective pressure. We have developed a modified heteroduplex tracking assay (HTA) capable of detecting the presence of single, specific mutations or sets of linked mutations. The initial application of this approach, termed multiple-site-specific (MSS) HTA, was directed toward the detection of mutations in the HIV-1 pro gene at positions 46, 48, 54, 82, 84, and 90, which are associated with resistance to multiple protease inhibitors. We demonstrate that MSS HTA is sensitive and largely specific to all targeted mutations. The assay allows the accurate and reproducible quantitation of viral subpopulations comprising 3% or more of the total population. Furthermore, we used MSS HTA in longitudinal studies of pro gene evolution in vitro and in vivo . In the examples shown here, populations turned over rapidly and more than one population was present frequently. To demonstrate the versatility of MSS HTA, we also constructed a probe sensitive to changes at positions 181 and 184 of the RT coding domain. Changes at these positions are involved in resistance to nevirapine and 2′,3′-dideoxy-3′-thiacytidine (3TC), respectively. This assay easily detected the evolution of resistance to 3TC. MSS HTA provides a rapid and sensitive approach for detecting the presence of and quantifying complex mixtures of distinct genotypes, including genetically linked mutations, and, as one example, represents a useful tool for following the evolution of drug resistance during failure of HIV-1 antiviral therapy.

SUMMARY The ‘‘CTCF code’’ hypothesis posits that CTCF pleiotropic functions are driven by recognition of diverse sequences through combinatorial use of its 11 zinc fingers (ZFs). This model, however, is supported by in vitro binding studies of a limited number of sequences. To study CTCF multivalency in vivo, we define ZF binding requirements at � 50,000 genomic sites in primary lymphocytes. We find that CTCF reads sequence diversity through ZF clustering. ZFs 4‐7 anchor CTCF to � 80% of targets containing the core motif. Nonconserved flanking sequences are recognized by ZFs 1‐2 and ZFs 8‐11 clusters, which also stabilize CTCF broadly. Alternatively, ZFs 9‐11 associate with a second phylogenetically conserved upstream motif at � 15% of its sites. Individually, ZFs increase overall binding and chromatin residence time. Unexpectedly, we also uncovered a conserved downstream DNA motif that destabilizes CTCF occupancy. Thus, CTCF associates with a wide array of DNA modules via combinatorial clustering of its 11 ZFs.

Human Burkitt lymphomas are divided into two main clinical variants: the endemic form, affecting African children infected with malaria and the Epstein-Barr virus, and the sporadic form, distributed across the rest of the world. However, whereas sporadic translocations decapitate Myc from 5′ proximal regulatory elements, most endemic events occur hundreds of kilobases away from Myc . The origin of these rearrangements and how they deregulate oncogenes at such distances remain unclear. We here recapitulate endemic Burkitt lymphoma-like translocations in plasmacytomas from uracil N -glycosylase and activation-induced cytidine deaminase-deficient mice. Mapping of translocation breakpoints using an acetylated histone H3 lysine 9 chromatin immunoprecipitation sequencing approach reveals Igh fusions up to ∼350 kb upstream of Myc or the related oncogene Mycn . A comprehensive analysis of epigenetic marks, PolII recruitment, and transcription in tumor cells demonstrates that the 3′ Igh enhancer (Eα) vastly remodels ∼450 kb of chromatin into translocated sequences, leading to significant polymerase occupancy and constitutive oncogene expression. We show that this long-range epigenetic reprogramming is directly proportional to the physical interaction of Eα with translocated sites. Our studies thus uncover the extent of epigenetic remodeling by Ig 3′ enhancers and provide a rationale for the long-range deregulation of translocated oncogenes in endemic Burkitt lymphomas. The data also shed light on the origin of endemic-like chromosomal rearrangements.

Birt-Hogg-Dubé (BHD) syndrome is an autosomal dominant disorder characterized by cutaneous fibrofolliculomas, pulmonary cysts, and kidney malignancies. Affected individuals carry germ line mutations in folliculin (FLCN), a tumor suppressor gene that becomes biallelically inactivated in kidney tumors by second-hit mutations. Similar to other factors implicated in kidney cancer, FLCN has been shown to modulate activation of mammalian target of rapamycin (mTOR). However, its precise in vivo function is largely unknown because germ line deletion of Flcn results in early embryonic lethality in animal models. Here, we describe mice deficient in the newly characterized folliculin-interacting protein 1 (Fnip1). In contrast to Flcn, Fnip1-/- mice develop normally, are not susceptible to kidney neoplasia, but display a striking pro-B cell block that is entirely independent of mTOR activity.We show that this developmental arrest results from rapid caspase-induced pre-B cell death, and that a Bcl2 transgene reconstitutes mature B-cell populations, respectively.We also demonstrate that conditional deletion of Flcn recapitulates the pro-B cell arrest of Fnip1-/- mice. Our studies thus demonstrate that the FLCN-FNIP complex deregulated in BHD syndrome is absolutely required for B-cell differentiation, and that it functions through both mTOR-dependent and independent pathways., This work was supported in part by the Intramural Research Program of NIAMS and NCI, Center for Cancer Research, National Institutes of Health (NIH). This project was funded in part with federal funds from the NCI, NIH, under contract no. HHSN261200800001E.

The cytidine deaminase AID hypermutates immunoglobulin genes but can also target oncogenes, leading to tumorigenesis. The extent of AID's promiscuity and its predilection for immunoglobulin genes are unknown. We report here that AID interacted broadly with promoter-proximal sequences associated with stalled polymerases and chromatin-activating marks. In contrast, genomic occupancy of replication protein A (RPA), an AID cofactor, was restricted to immunoglobulin genes. The recruitment of RPA to the immunoglobulin loci was facilitated by phosphorylation of AID at Ser38 and Thr140. We propose that stalled polymerases recruit AID, thereby resulting in low frequencies of hypermutation across the B cell genome. Efficient hypermutation and switch recombination required AID phosphorylation and correlated with recruitment of RPA. Our findings provide a rationale for the oncogenic role of AID in B cell malignancy.

Background There is an increasing number of complete and incomplete virus genome sequences available in public databases. This large body of sequence data harbors information about epidemiology, phylogeny, and virulence. Several specialized databases, such as the NCBI Influenza Virus Resource or the Los Alamos HIV database, offer sophisticated query interfaces along with integrated exploratory data analysis tools for individual virus species to facilitate extracting this information. Thus far, there has not been a comprehensive database for dengue virus, a significant public health threat. Results We have created an integrated web resource for dengue virus. The technology developed for the NCBI Influenza Virus Resource has been extended to process non-segmented dengue virus genomes. In order to allow efficient processing of the dengue genome, which is large in comparison with individual influenza segments, we developed an offline pre-alignment procedure which generates a multiple sequence alignment of all dengue sequences. The pre-calculated alignment is then used to rapidly create alignments of sequence subsets in response to user queries. This improvement in technology will also facilitate the incorporation of additional virus species in the future. The set of virus-specific databases at NCBI, which will be referred to as Virus Variation Resources (VVR), allow users to build complex queries against virus-specific databases and then apply exploratory data analysis tools to the results. The metadata is automatically collected where possible, and extended with data extracted from the literature. Conclusion The NCBI Dengue Virus Resource integrates dengue sequence information with relevant metadata (sample collection time and location, disease severity, serotype, sequenced genome region) and facilitates retrieval and preliminary analysis of dengue sequences using integrated web analysis and visualization tools.

Immunization with recombinant proteins may provide a safer alternative to live vaccinia virus for prophylaxis of poxvirus infections. Although antibody protects against vaccinia virus infection, the mechanism is not understood and the selection of immunogens is daunting as there are dozens of surface proteins and two infectious forms known as the mature virion (MV) and the enveloped virion (EV). Our previous studies showed that mice immunized with soluble forms of EV membrane proteins A33 and B5 and MV membrane protein L1 or passively immunized with antibodies to these proteins survived an intranasal challenge with vaccinia virus. The present study compared MV protein A27, which has a role in virus attachment to glycosaminoglycans on the cell surface, to L1 with respect to immunogenicity and protection. Although mice developed similar levels of neutralizing antibody after immunizations with A27 or L1, A27-immunized mice exhibited more severe disease upon an intranasal challenge with vaccinia virus. In addition, mice immunized with A27 and A33 were not as well protected as mice receiving L1 and A33. Polyclonal rabbit anti-A27 and anti-L1 IgG had equivalent MV-neutralizing activities when measured by the prevention of infection of human or mouse cells or cells deficient in glycosaminoglycans or by adding antibody prior to or after virus adsorption. Nevertheless, the passive administration of antibody to A27 was poorly protective compared to the antibody to L1. These studies raise questions regarding the basis for antibody protection against poxvirus disease and highlight the importance of animal models for the early evaluation of vaccine candidates.

Recombinant proteins are being evaluated as smallpox and monkeypox vaccines because of their perceived safety compared to live vaccinia virus. Previously, we demonstrated that three or more injections of a Ribi-type adjuvant with a combination of three proteins from the outer membranes of intracellular (L1 protein) and extracellular (A33 and B5 proteins) forms of vaccinia virus protected mice against a lethal intranasal challenge with vaccinia virus. Here, we compared several adjuvants and found that QS-21 and to a lesser extent alum+CpG oligodeoxynucleotides accelerated and enhanced neutralizing antibody responses to a mixture of L1 and A33 proteins, provided the highest ratio of IgG2a to IgG1 isotype response, and protected mice against disease and death after only two immunizations 3 weeks apart. In addition, monkeys immunized with recombinant vaccinia virus proteins and QS-21 developed neutralizing antibody to monkeypox virus and had reduced virus load, skin lesions, and morbidity compared to the non-immunized group following monkeypox virus challenge.

The protein content of vaccinia virus mature virions, purified by rate zonal and isopycnic centrifugations and solubilized by SDS or a solution of urea and thiourea, was determined by the accurate mass and time tag technology which uses both tandem mass spectrometry and Fourier transform-ion cyclotron resonance mass spectrometry to detect tryptic peptides separated by high-resolution liquid chromatography. Eighty vaccinia virus-encoded proteins representing 37% of the 218 genes annotated in the complete genome sequence were detected in at least three analyses. Ten proteins accounted for approximately 80% of the virion mass. Thirteen identified proteins were not previously reported as components of virions. On the other hand, 8 previously described virion proteins were not detected here, presumably due to technical reasons including small size and hydrophobicity. In addition to vaccinia virus-encoded proteins, 24 host proteins omitting isoforms were detected. The most abundant of these were cytoskeletal proteins, heat shock proteins and proteins involved in translation.

Protease inhibitors represent some of the most potent agents available for therapeutic strategies designed to inhibit human immunodeficiency virus type 1 (HIV-1) replication. Under certain circumstances the virus develops resistance to the inhibitor, thereby negating the benefits of this therapy. We have carried out selections for high-level resistance to each of three protease inhibitors (indinavir, ritonavir, and saquinavir) in cell culture. Mutations accumulated over most of the course of the increasing selective pressure. There was significant overlap in the identity of the mutations selected with the different inhibitors, and this gave rise to high levels of cross-resistance. Virus particles from the resistant variants all showed defects in processing at the NC/p1 protease cleavage site in Gag. Selections with pairs of inhibitors yielded similar patterns of resistance mutations. A virus that could replicate at near-toxic levels of the three protease inhibitors combined was selected. Theprosequence of this virus was similar to that of the viruses that had been selected for high-level resistance to each of the drugs singly. Finally, a molecular clone carrying the eight most common resistance mutations seen in these selections was characterized. The sequence of this virus was relatively stable during selection for revertants in spite of displaying poor processing at the NC/p1 site and having significantly reduced fitness. These results reveal patterns of drug resistance that extend to near the limits of attainable selective pressure with these inhibitors and confirm the patterns of cross-resistance for these three inhibitors and the attenuation of virion protein processing and fitness that accompanies high-level resistance.

The evolution of human immunodeficiency virus type 1 (HIV-1) strains with reduced susceptibility to protease inhibitors (PIs) is a major cause of PI treatment failure. A subset of subjects failing a therapy regimen containing the PI nelfinavir developed mutations at position 88 in the protease region. The N88S mutation occurring in some of these subjects induces amprenavir hypersusceptibility and a reduction of fitness and replication capacity. Here we demonstrate that substitutions L63P and V77I in protease, in combination, partially compensate for the loss of fitness, loss of replication capacity, loss of specific infectivity, and aberrant Gag processing induced by the N88S mutation. In addition, these mutations partially ablate amprenavir hypersusceptibility. Addition of mutation M46L to a strain harboring mutations L63P, V77I, and N88S resulted in a reduction of fitness and infectivity without changing Gag-processing efficiency, while amprenavir hypersusceptibility was further diminished. The ratio of reverse transcriptase activity to p24 protein was reduced in this strain compared to that in the other variants, suggesting that the M46L effect on fitness occurred through a mechanism different from a Gag-processing defect. We utilized these mutant strains to undertake a systematic comparison of indirect, single, cycle-based measures of fitness with direct, replication-based fitness assays and demonstrated that both yield consistent results. However, we observed that the magnitude of the fitness loss for one of the mutants varied depending on the assay used. Fil: Resch, W.. University of North Carolina; Estados Unidos Fil: Ziermann, R.. ViroLogic ; Estados Unidos Fil: Parkin, N.. ViroLogic ; Estados Unidos Fil: Gamarnik, Andrea Vanesa. ViroLogic ; Estados Unidos. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Parque Centenario. Instituto de Investigaciones Bioquímicas de Buenos Aires. Fundación Instituto Leloir. Instituto de Investigaciones Bioquímicas de Buenos Aires; Argentina Fil: Swanstrom, R.. University of North Carolina; Estados Unidos

Human immunodeficiency virus type 1 (HIV-1) only recently established an epidemic world-wide infection in the human population. The virus persists in the human host through active replication and is able to avoid clearance by the immune system. Active replication is an important component of the rapid evolutionary potential of HIV-1, a potential which manifests itself in the evolution of immune escape variants, drug resistant variants, and variants with the ability to use different cell surface coreceptors in conjunction with CD4. Multiple zoonotic introductions, compartmentalization of virus replication in the body, and genetic bottlenecks associated with sampling during transmission, antiretroviral therapy, and geographic and/or host population isolation further contribute to the range of sequences present in extant viruses. The sum of the history of all of these phenomena is reflected in HIV-1 sequence variability, and most of these phenomena are ongoing today. Here we review the use of HIV-1 sequence variability to explore its underlying biology.

The critical initial step in V(D)J recombination, binding of RAG1 and RAG2 to recombination signal sequences flanking antigen receptor V, D, and J gene segments, has not previously been characterized in vivo. Here, we demonstrate that RAG protein binding occurs in a highly focal manner to a small region of active chromatin encompassing Ig kappa and Tcr alpha J gene segments and Igh and Tcr beta J and J-proximal D gene segments. Formation of these small RAG-bound regions, which we refer to as recombination centers, occurs in a developmental stage- and lineage-specific manner. Each RAG protein is independently capable of specific binding within recombination centers. While RAG1 binding was detected only at regions containing recombination signal sequences, RAG2 binds at thousands of sites in the genome containing histone 3 trimethylated at lysine 4. We propose that recombination centers coordinate V(D)J recombination by providing discrete sites within which gene segments are captured for recombination.

We have assembled two sets of HIV-1 V3 sequences with defined epidemiologic relationships associated with experimentally determined coreceptor usage or MT-2 cell tropism. These data sets were used for three purposes. First, they were employed to test existing methods for predicting coreceptor usage and MT-2 cell tropism. Of these methods, the presence of one basic amino acid at position 11 or 25 proved to be most reliable for both phenotypic classifications, although its predictive power for the X4 phenotype was less than 50%. Second, we used the sequence sets to train neural networks to infer coreceptor usage from V3 genotype with better success than the best available motif-based method, and with a predictive power equal to that of the best motif-based method for MT-2 cell tropism. Third, we used the sequence sets to reexamine patterns of variability associated with the different phenotypes, and we showed that the phenotype-associated sequence patterns could be reproduced from large sets of V3 sequences using phenotypes predicted by the trained neural network.

SummaryThe RAG1 endonuclease, together with its cofactor RAG2, is essential for V(D)J recombination but is a potent threat to genome stability. The sources of RAG1 mis-targeting and the mechanisms that have evolved to suppress it are poorly understood. Here, we report that RAG1 associates with chromatin at thousands of active promoters and enhancers in the genome of developing lymphocytes. The mouse and human genomes appear to have responded by reducing the abundance of “cryptic” recombination signals near RAG1 binding sites. This depletion operates specifically on the RSS heptamer, whereas nonamers are enriched at RAG1 binding sites. Reversing this RAG-driven depletion of cleavage sites by insertion of strong recombination signals creates an ectopic hub of RAG-mediated V(D)J recombination and chromosomal translocations. Our findings delineate rules governing RAG binding in the genome, identify areas at risk of RAG-mediated damage, and highlight the evolutionary struggle to accommodate programmed DNA damage in developing lymphocytes.

Summary The antibody gene mutator activation-induced cytidine deaminase (AID) promiscuously damages oncogenes, leading to chromosomal translocations and tumorigenesis. Why nonimmunoglobulin loci are susceptible to AID activity is unknown. Here, we study AID-mediated lesions in the context of nuclear architecture and the B cell regulome. We show that AID targets are not randomly distributed across the genome but are predominantly grouped within super-enhancers and regulatory clusters. Unexpectedly, in these domains, AID deaminates active promoters and eRNA + enhancers interconnected in some instances over megabases of linear chromatin. Using genome editing, we demonstrate that 3D-linked targets cooperate to recruit AID-mediated breaks. Furthermore, a comparison of hypermutation in mouse B cells, AID-induced kataegis in human lymphomas, and translocations in MEFs reveals that AID damages different genes in different cell types. Yet, in all cases, the targets are predominantly associated with topological complex, highly transcribed super-enhancers, demonstrating that these compartments are key mediators of AID recruitment.