Belt and braces: Two escape ways to maintain the cassette reservoir of large chromosomal integrons.

Integrons are adaptive devices that capture, stockpile, shuffle and express gene cassettes thereby sampling combinatorial phenotypic diversity. Some integrons called sedentary chromosomal integrons (SCIs) can be massive structures containing hundreds of cassettes. Since most of these cassettes are non-expressed, it is not clear how they remain stable over long evolutionary timescales. Recently, it was found that the experimental inversion of the SCI of Vibrio cholerae led to a dramatic increase of the cassette excision rate associated with a fitness defect. Here, we question the evolutionary sustainability of this apparently counter selected genetic context. Through experimental evolution, we find that the integrase is rapidly inactivated and that the inverted SCI can recover its original orientation by homologous recombination between two insertion sequences (ISs) present in the array. These two outcomes of SCI inversion restore the normal growth and prevent the loss of cassettes, enabling SCIs to retain their roles as reservoirs of functions. These results illustrate a nice interplay between gene orientation, genome rearrangement, bacterial fitness and demonstrate how integrons can benefit from their embedded ISs.

Chlamydia trachomatis development requires both host glycolysis and oxidative phosphorylation but has only minor effects on these pathways.

The obligate intracellular bacteria Chlamydia trachomatis obtain all nutrients from the cytoplasm of their epithelial host cells and stimulate glucose uptake by these cells. They even hijack host ATP, exerting a strong metabolic pressure on their host at the peak of the proliferative stage of their developmental cycle. However, it is largely unknown whether infection modulates the metabolism of the host cell. Also, the reliance of the bacteria on host metabolism might change during their progression through their biphasic developmental cycle. Herein, using primary epithelial cells and 2 cell lines of nontumoral origin, we showed that between the 2 main ATP-producing pathways of the host, oxidative phosphorylation (OxPhos) remained stable and glycolysis was slightly increased. Inhibition of either pathway strongly reduced bacterial proliferation, implicating that optimal bacterial growth required both pathways to function at full capacity. While we found C. trachomatis displayed some degree of energetic autonomy in the synthesis of proteins expressed at the onset of infection, functional host glycolysis was necessary for the establishment of early inclusions, whereas OxPhos contributed less. These observations correlated with the relative contributions of the pathways in maintaining ATP levels in epithelial cells, with glycolysis contributing the most. Altogether, this work highlights the dependence of C. trachomatis on both host glycolysis and OxPhos for efficient bacterial replication. However, ATP consumption appears at equilibrium with the normal production capacity of the host and the bacteria, so that no major shift between these pathways is required to meet bacterial needs.

Serpentine: a flexible 2D binning method for differential Hi-C analysis.

MOTIVATION: Hi-C contact maps reflect the relative contact frequencies between pairs of genomic loci, quantified through deep sequencing. Differential analyses of these maps enable downstream biological interpretations. However, the multi-fractal nature of the chromatin polymer inside the cellular envelope results in contact frequency values spanning several orders of magnitude: contacts between loci pairs separated by large genomic distances are much sparser than closer pairs. The same is true for poorly covered regions, such as repeated sequences. Both distant and poorly covered regions translate into low signal-to-noise ratios. There is no clear consensus to address this limitation. RESULTS: We present Serpentine, a fast, flexible procedure operating on raw data, which considers the contacts in each region of a contact map. Binning is performed only when necessary on noisy regions, preserving informative ones. This results in high-quality, low-noise contact maps that can be conveniently visualized for rigorous comparative analyses. AVAILABILITY AND IMPLEMENTATION: Serpentine is available on the PyPI repository and https://github.com/koszullab/serpentine; documentation and tutorials are provided at https://serpentine.readthedocs.io/en/latest/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Immune response profiling of patients with spondyloarthritis reveals signalling networks mediating TNF-blocker function in vivo.

OBJECTIVES: Antitumour necrosis factor (TNF) therapy has revolutionised treatment of several chronic inflammatory diseases, including spondyloarthritis (SpA). However, TNF inhibitors (TNFi) are not effective in all patients and the biological basis for treatment failure remains unknown. We have analysed induced immune responses to define the mechanism of action of TNF blockers in SpA and to identify immunological correlates of responsiveness to TNFi. METHODS: Immune responses to microbial and pathway-specific stimuli were analysed in peripheral blood samples from 80 patients with axial SpA before and after TNFi treatment, using highly standardised whole-blood stimulation assays. Cytokines and chemokines were measured in a Clinical Laboratory Improvement Amendments (CLIA)-certified laboratory, and gene expression was monitored using nCounter assays. RESULTS: Anti-TNF therapy induced profound changes in patients' innate immune responses. TNFi action was selective, and had only minor effects on Th1/Th17 immunity. Modular transcriptional repertoire analysis identified prostaglandin E2 synthesis and signalling, leucocyte recirculation, macrophage polarisation, dectin and interleukin (IL)-1 signalling, as well as the nuclear factor kappa B (NF-kB) transcription factor family as key pathways targeted by TNF blockers in vivo. Analysis of induced immune responses before treatment initiation revealed that expression of molecules associated with leucocyte adhesion and invasion, chemotaxis and IL-1 signalling are correlated with therapeutic responses to anti-TNF. CONCLUSIONS: We show that TNFi target multiple immune cell pathways that cooperate to resolve inflammation. We propose that immune response profiling provides new insight into the biology of TNF-blocker action in patients and can identify signalling pathways associated with therapeutic responses to biological therapies.

A Parent-of-Origin Effect Impacts the Phenotype in Low Penetrance Retinoblastoma Families Segregating the c.1981C>T/p.Arg661Trp Mutation of RB1.

Retinoblastoma (Rb), the most common pediatric intraocular neoplasm, results from inactivation of both alleles of the RB1 tumor suppressor gene. The second allele is most commonly lost, as demonstrated by loss of heterozygosity studies. RB1 germline carriers usually develop bilateral tumors, but some Rb families display low penetrance and variable expressivity. In order to decipher the underlying mechanisms, 23 unrelated low penetrance pedigrees segregating the common c.1981C>T/p.Arg661Trp mutation and other low penetrance mutations were studied. In families segregating the c.1981C>T mutation, we demonstrated, for the first time, a correlation between the gender of the transmitting carrier and penetrance, as evidenced by Fisher's exact test: the probability of being unaffected is 90.3% and 32.5% when the mutation is inherited from the mother and the father, respectively (p-value = 7.10(-7). Interestingly, a similar correlation was observed in families segregating other low penetrance alleles. Consequently, we investigated the putative involvement of an imprinted, modifier gene in low penetrance Rb. We first ruled out a MED4-driven mechanism by MED4 methylation and expression analyses. We then focused on the differentially methylated CpG85 island located in intron 2 of RB1 and showing parent-of-origin-specific DNA methylation. This differential methylation promotes expression of the maternal c.1981C>T allele. We propose that the maternally inherited c.1981C>T/p.Arg661Trp allele retains sufficient tumor suppressor activity to prevent retinoblastoma development. In contrast, when the mutation is paternally transmitted, the low residual activity would mimic a null mutation and subsequently lead to retinoblastoma. This implies that the c.1981C>T mutation is not deleterious per se but needs to be destabilized in order to reach pRb haploinsufficiency and initiate tumorigenesis. We suggest that this phenomenon might be a general mechanism to explain phenotypic differences in low penetrance Rb families.

Functional Assessment of Genetic Variants with Outcomes Adapted to Clinical Decision-Making.

Understanding the medical effect of an ever-growing number of human variants detected is a long term challenge in genetic counseling. Functional assays, based on in vitro or in vivo evaluations of the variant effects, provide essential information, but they require robust statistical validation, as well as adapted outputs, to be implemented in the clinical decision-making process. Here, we assessed 25 pathogenic and 15 neutral missense variants of the BRCA1 breast/ovarian cancer susceptibility gene in four BRCA1 functional assays. Next, we developed a novel approach that refines the variant ranking in these functional assays. Lastly, we developed a computational system that provides a probabilistic classification of variants, adapted to clinical interpretation. Using this system, the best functional assay exhibits a variant classification accuracy estimated at 93%. Additional theoretical simulations highlight the benefit of this ready-to-use system in the classification of variants after functional assessment, which should facilitate the consideration of functional evidences in the decision-making process after genetic testing. Finally, we demonstrate the versatility of the system with the classification of siRNAs tested for human cell growth inhibition in high throughput screening.

Yeast cells reveal the misfolding and the cellular mislocalization of the human BRCA1 protein.

Understanding the effect of an ever-growing number of human variants detected by genome sequencing is a medical challenge. The yeast Saccharomyces cerevisiae model has held attention for its capacity to monitor the functional impact of missense mutations found in human genes, including the BRCA1 breast and ovarian cancer susceptibility gene. When expressed in yeast, the wild-type full-length BRCA1 protein forms a single nuclear aggregate and induces a growth inhibition. Both events are modified by pathogenic mutations of BRCA1. However, the biological processes behind these events in yeast remain to be determined. Here, we show that the BRCA1 nuclear aggregation and the growth inhibition are sensitive to misfolding effects induced by missense mutations. Moreover, misfolding mutations impair the nuclear targeting of BRCA1 in yeast cells and in a human cell line. In conclusion, we establish a connection between misfolding and nuclear transport impairment, and we illustrate that yeast is a suitable model to decipher the effect of misfolding mutations.

Cell-type-specific replication initiation programs set fragility of the FRA3B fragile site.

Common fragile sites have long been identified by cytogeneticists as chromosomal regions prone to breakage upon replication stress. They are increasingly recognized to be preferential targets for oncogene-induced DNA damage in pre-neoplastic lesions and hotspots for chromosomal rearrangements in various cancers. Common fragile site instability was attributed to the fact that they contain sequences prone to form secondary structures that may impair replication fork movement, possibly leading to fork collapse resulting in DNA breaks. Here we show, in contrast to this view, that the fragility of FRA3B--the most active common fragile site in human lymphocytes--does not rely on fork slowing or stalling but on a paucity of initiation events. Indeed, in lymphoblastoid cells, but not in fibroblasts, initiation events are excluded from a FRA3B core extending approximately 700 kilobases, which forces forks coming from flanking regions to cover long distances in order to complete replication. We also show that origins of the flanking regions fire in mid-S phase, leaving the site incompletely replicated upon fork slowing. Notably, FRA3B instability is specific to cells showing this particular initiation pattern. The fact that both origin setting and replication timing are highly plastic in mammalian cells explains the tissue specificity of common fragile site instability we observed. Thus, we propose that common fragile sites correspond to the latest initiation-poor regions to complete replication in a given cell type. For historical reasons, common fragile sites have been essentially mapped in lymphocytes. Therefore, common fragile site contribution to chromosomal rearrangements in tumours should be reassessed after mapping fragile sites in the cell type from which each tumour originates.

BRCA1 Circos: a visualisation resource for functional analysis of missense variants.

BACKGROUND: Inactivating germline mutations in the tumour suppressor gene BRCA1 are associated with a significantly increased risk of developing breast and ovarian cancer. A large number (>1500) of unique BRCA1 variants have been identified in the population and can be classified as pathogenic, non-pathogenic or as variants of unknown significance (VUS). Many VUS are rare missense variants leading to single amino acid changes. Their impact on protein function cannot be directly inferred from sequence information, precluding assessment of their pathogenicity. Thus, functional assays are critical to assess the impact of these VUS on protein activity. BRCA1 is a multifunctional protein and different assays have been used to assess the impact of variants on different biochemical activities and biological processes. METHODS AND RESULTS: To facilitate VUS analysis, we have developed a visualisation resource that compiles and displays functional data on all documented BRCA1 missense variants. BRCA1 Circos is a web-based visualisation tool based on the freely available Circos software package. The BRCA1 Circos web tool (http://research.nhgri.nih.gov/bic/circos/) aggregates data from all published BRCA1 missense variants for functional studies, harmonises their results and presents various functionalities to search and interpret individual-level functional information for each BRCA1 missense variant. CONCLUSIONS: This research visualisation tool will serve as a quick one-stop publically available reference for all the BRCA1 missense variants that have been functionally assessed. It will facilitate meta-analysis of functional data and improve assessment of pathogenicity of VUS.

Integron cassettes integrate into bacterial genomes via widespread non-classical attG sites.

Integrons are genetic elements involved in bacterial adaptation which capture, shuffle and express genes encoding adaptive functions embedded in cassettes. These events are governed by the integron integrase through site-specific recombination between attC and attI integron sites. Using computational and molecular genetic approaches, here we demonstrate that the integrase also catalyses cassette integration into bacterial genomes outside of its known att sites. Once integrated, these cassettes can be expressed if located near bacterial promoters and can be excised at the integration point or outside, inducing chromosomal modifications in the latter case. Analysis of more than 5 × 105 independent integration events revealed a very large genomic integration landscape. We identified consensus recombination sequences, named attG sites, which differ greatly in sequence and structure from classical att sites. These results unveil an alternative route for dissemination of adaptive functions in bacteria and expand the role of integrons in bacterial evolution.

A guide for functional analysis of BRCA1 variants of uncertain significance.

Germline mutations in the tumor suppressor gene BRCA1 confer an estimated lifetime risk of 56-80% for breast cancer and 15-60% for ovarian cancer. Since the mid 1990s when BRCA1 was identified, genetic testing has revealed over 1,500 unique germline variants. However, for a significant number of these variants, the effect on protein function is unknown making it difficult to infer the consequences on risks of breast and ovarian cancers. Thus, many individuals undergoing genetic testing for BRCA1 mutations receive test results reporting a variant of uncertain clinical significance (VUS), leading to issues in risk assessment, counseling, and preventive care. Here, we describe functional assays for BRCA1 to directly or indirectly assess the impact of a variant on protein conformation or function and how these results can be used to complement genetic data to classify a VUS as to its clinical significance. Importantly, these methods may provide a framework for genome-wide pathogenicity assignment.

Nucleotide supply, not local histone acetylation, sets replication origin usage in transcribed regions.

In eukaryotes, only a fraction of replication origins fire at each S phase. Local histone acetylation was proposed to control firing efficiency of origins, but conflicting results were obtained. We report that local histone acetylation does not reflect origin efficiencies along the adenosine monophosphate deaminase 2 locus in mammalian fibroblasts. Reciprocally, modulation of origin efficiency does not affect acetylation. However, treatment with a deacetylase inhibitor changes the initiation pattern. We demonstrate that this treatment alters pyrimidine biosynthesis and decreases fork speed, which recruits latent origins. Our findings reconcile results that seemed inconsistent and reveal an unsuspected effect of deacetylase inhibitors on replication dynamics.

Genome-wide mapping of individual replication fork velocities using nanopore sequencing.

Little is known about replication fork velocity variations along eukaryotic genomes, since reference techniques to determine fork speed either provide no sequence information or suffer from low throughput. Here we present NanoForkSpeed, a nanopore sequencing-based method to map and extract the velocity of individual forks detected as tracks of the thymidine analogue bromodeoxyuridine incorporated during a brief pulse-labelling of asynchronously growing cells. NanoForkSpeed retrieves previous Saccharomyces cerevisiae mean fork speed estimates (≈2 kb/min) in the BT1 strain exhibiting highly efficient bromodeoxyuridine incorporation and wild-type growth, and precisely quantifies speed changes in cells with altered replisome progression or exposed to hydroxyurea. The positioning of >125,000 fork velocities provides a genome-wide map of fork progression based on individual fork rates, showing a uniform fork speed across yeast chromosomes except for a marked slowdown at known pausing sites.

High-affinity autoreactive plasma cells disseminate through multiple organs in patients with immune thrombocytopenic purpura.

The major therapeutic goal for immune thrombocytopenic purpura (ITP) is to restore normal platelet counts using drugs to promote platelet production or by interfering with mechanisms responsible for platelet destruction. Eighty percent of patients with ITP possess anti-integrin αIIbß3 IgG autoantibodies that cause platelet opsonization and phagocytosis. The spleen is considered the primary site of autoantibody production by autoreactive B cells and platelet destruction. The immediate failure in approximately 50% of patients to recover a normal platelet count after anti-CD20 rituximab-mediated B cell depletion and splenectomy suggests that autoreactive, rituximab-resistant, IgG-secreting B cells (IgG-SCs) reside in other anatomical compartments. We analyzed more than 3,300 single IgG-SCs from spleen, bone marrow, and/or blood of 27 patients with ITP, revealing high interindividual variability in affinity for αIIbß3, with variations over 3 logs. IgG-SC dissemination and range of affinities were, however, similar for each patient. Longitudinal analysis of autoreactive IgG-SCs upon treatment with the anti-CD38 mAb daratumumab demonstrated variable outcomes, from complete remission to failure with persistence of high-affinity anti-αIIbß3 IgG-SCs in the bone marrow. This study demonstrates the existence and dissemination of high-affinity autoreactive plasma cells in multiple anatomical compartments of patients with ITP that may cause the failure of current therapies.

Assessment of human Nter and Cter BRCA1 mutations using growth and localization assays in yeast.

A large number of missense mutations have been identified within the tumor suppressor gene BRCA1. Most of them, called "variants of unknown significance" (VUS), cannot be classified as pathogenic or neutral by genetic methods, which complicates their cancer risk assessment. Functional assays have been developed to circumvent this uncertainty. They aim to determine how VUS impact the BRCA1 protein structure or function, thereby giving an indication of their potential to cause cancer. So far, three relevant assays have been designed in yeast and used on large sets of variants. However, they are limited to variants mapped in restricted domains of BRCA1. One of them, the small colony phenotype (SCP) assay, monitors the BRCA1-dependent growth of yeast colonies that increases with pathogenic but not neutral mutations positioned in the Cter region. Here, we extend this assay to the Nter part of BRCA1. We also designed a new assay, called the "yeast localization phenotype (YLP) assay," based on the accumulation of BRCA1 in a single inclusion body in the yeast nucleus. This phenotype is altered by variants positioned both in the Nter and Cter regions. Together, these assays provide new perspectives for the functional assessment of BRCA1 mutations in yeast.

MetaHiC phage-bacteria infection network reveals active cycling phages of the healthy human gut.

Bacteriophages play important roles in regulating the intestinal human microbiota composition, dynamics, and homeostasis, and characterizing their bacterial hosts is needed to understand their impact. We applied a metagenomic Hi-C approach on 10 healthy human gut samples to unveil a large infection network encompassing more than 6000 interactions bridging a metagenomic assembled genomes (MAGs) and a phage sequence, allowing to study in situ phage-host ratio. Whereas three-quarters of these sequences likely correspond to dormant prophages, 5% exhibit a much higher coverage than their associated MAG, representing potentially actively replicating phages. We detected 17 sequences of members of the crAss-like phage family, whose hosts diversity remained until recently relatively elusive. For each of them, a unique bacterial host was identified, all belonging to different genus of Bacteroidetes. Therefore, metaHiC deciphers infection network of microbial population with a high specificity paving the way to dynamic analysis of mobile genetic elements in complex ecosystems.

The anti-IgE mAb omalizumab induces adverse reactions by engaging Fcγ receptors.

Omalizumab is an anti-IgE monoclonal antibody (mAb) approved for the treatment of severe asthma and chronic spontaneous urticaria. Use of omalizumab is associated with reported side effects ranging from local skin inflammation at the injection site to systemic anaphylaxis. To date, the mechanisms through which omalizumab induces adverse reactions are still unknown. Here, we demonstrated that immune complexes formed between omalizumab and IgE can induce both skin inflammation and anaphylaxis through engagement of IgG receptors (FcγRs) in FcγR-humanized mice. We further developed an Fc-engineered mutant version of omalizumab, and demonstrated that this mAb is equally potent as omalizumab at blocking IgE-mediated allergic reactions, but does not induce FcγR-dependent adverse reactions. Overall, our data indicate that omalizumab can induce skin inflammation and anaphylaxis by engaging FcγRs, and demonstrate that Fc-engineered versions of the mAb could be used to reduce such adverse reactions.

High-throughput single-cell activity-based screening and sequencing of antibodies using droplet microfluidics.

Mining the antibody repertoire of plasma cells and plasmablasts could enable the discovery of useful antibodies for therapeutic or research purposes1. We present a method for high-throughput, single-cell screening of IgG-secreting primary cells to characterize antibody binding to soluble and membrane-bound antigens. CelliGO is a droplet microfluidics system that combines high-throughput screening for IgG activity, using fluorescence-based in-droplet single-cell bioassays2, with sequencing of paired antibody V genes, using in-droplet single-cell barcoded reverse transcription. We analyzed IgG repertoire diversity, clonal expansion and somatic hypermutation in cells from mice immunized with a vaccine target, a multifunctional enzyme or a membrane-bound cancer target. Immunization with these antigens yielded 100-1,000 IgG sequences per mouse. We generated 77 recombinant antibodies from the identified sequences and found that 93% recognized the soluble antigen and 14% the membrane antigen. The platform also allowed recovery of ~450-900 IgG sequences from ~2,200 IgG-secreting activated human memory B cells, activated ex vivo, demonstrating its versatility.

Author Correction: High-throughput single-cell activity-based screening and sequencing of antibodies using droplet microfluidics.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.