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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

LGP2 binds to PACT to regulate RIG-I- and MDA5-mediated antiviral responses.

The retinoic acid-inducible gene I (RIG-I)-like receptors (RLRs) RIG-I, MDA5, and LGP2 stimulate inflammatory and antiviral responses by sensing nonself RNA molecules produced during viral replication. Here, we investigated how LGP2 regulates the RIG-I- and MDA5-dependent induction of type I interferon (IFN) signaling and showed that LGP2 interacted with different components of the RNA-silencing machinery. We identified a direct protein-protein interaction between LGP2 and the IFN-inducible, double-stranded RNA binding protein PACT. The LGP2-PACT interaction was mediated by the regulatory C-terminal domain of LGP2 and was necessary for inhibiting RIG-I-dependent responses and for amplifying MDA5-dependent responses. We described a point mutation within LGP2 that disrupted the LGP2-PACT interaction and led to the loss of LGP2-mediated regulation of RIG-I and MDA5 signaling. These results suggest a model in which the LGP2-PACT interaction regulates the inflammatory responses mediated by RIG-I and MDA5 and enables the cellular RNA-silencing machinery to coordinate with the innate immune response.

Transcription-dependent regulation of replication dynamics modulates genome stability.

Common fragile sites (CFSs) are loci that are hypersensitive to replication stress and hotspots for chromosomal rearrangements in cancers. CFSs replicate late in S phase, are cell-type specific and nest in large genes. The relative impact of transcription-replication conflicts versus a low density in initiation events on fragility is currently debated. Here we addressed the relationships between transcription, replication, and instability by manipulating the transcription of endogenous large genes in chicken and human cells. We found that inducing low transcription with a weak promoter destabilized large genes, whereas stimulating their transcription with strong promoters alleviated instability. Notably, strong promoters triggered a switch to an earlier replication timing, supporting a model in which high transcription levels give cells more time to complete replication before mitosis. Transcription could therefore contribute to maintaining genome integrity, challenging the dominant view that it is exclusively a threat.

The Loss of Expression of a Single Type 3 Effector (CT622) Strongly Reduces Chlamydia trachomatis Infectivity and Growth.

Invasion of epithelial cells by the obligate intracellular bacterium Chlamydia trachomatis results in its enclosure inside a membrane-bound compartment termed an inclusion. The bacterium quickly begins manipulating interactions between host intracellular trafficking and the inclusion interface, diverging from the endocytic pathway and escaping lysosomal fusion. We have identified a previously uncharacterized protein, CT622, unique to the Chlamydiaceae, in the absence of which most bacteria failed to establish a successful infection. CT622 is abundant in the infectious form of the bacteria, in which it associates with CT635, a putative novel chaperone protein. We show that CT622 is translocated into the host cytoplasm via type three secretion throughout the developmental cycle of the bacteria. Two separate domains of roughly equal size have been identified within CT622 and a 1.9 Å crystal structure of the C-terminal domain has been determined. Genetic disruption of ct622 expression resulted in a strong bacterial growth defect, which was due to deficiencies in proliferation and in the generation of infectious bacteria. Our results converge to identify CT622 as a secreted protein that plays multiple and crucial roles in the initiation and support of the C. trachomatis growth cycle. They reveal that genetic disruption of a single effector can deeply affect bacterial fitness.

Contribution of germline deleterious variants in the RAD51 paralogs to breast and ovarian cancers.

RAD51 paralogs (RAD51B, RAD51C, RAD51D, XRCC2, and XRCC3) have recently been involved in breast and ovarian cancer predisposition: RAD51B, RAD51C, and RAD51D in ovarian cancer, RAD51B and XRCC2 in breast cancer. The aim of this study was to estimate the contribution of deleterious variants in the five RAD51 paralogs to breast and ovarian cancers. The five RAD51 paralog genes were analyzed by next-generation sequencing technologies in germline DNA from 2649 consecutive patients diagnosed with breast and/or ovarian cancer. Twenty-one different deleterious variants were identified in the RAD51 paralogs in 30 patients: RAD51B (n = 4), RAD51C (n = 12), RAD51D (n = 7), XRCC2 (n = 2), and XRCC3 (n = 5). The overall deleterious variant rate was 1.13% (95% confidence interval (CI): 0.72-1.55%) (30/2649), including 15 variants in breast cancer only cases (15/2063; 0.73% (95% CI: 0.34-1.11%)) and 15 variants in cases with at least one ovarian cancer (15/570; 2.63% (95% CI: 1.24-4.02%)). This study is the first evaluation of the five RAD51 paralogs in breast and ovarian cancer predisposition and it demonstrates that deleterious variants can be present in breast cancer only cases. Moreover, this is the first time that XRCC3 deleterious variants have been identified in breast and ovarian cancer cases.

Heritability of the melatonin synthesis variability in autism spectrum disorders.

Autism Spectrum Disorders (ASD) are heterogeneous neurodevelopmental disorders with a complex genetic architecture. They are characterized by impaired social communication, stereotyped behaviors and restricted interests and are frequently associated with comorbidities such as intellectual disability, epilepsy and severe sleep disorders. Hyperserotonemia and low melatonin levels are among the most replicated endophenotypes reported in ASD, but their genetic causes remain largely unknown. Based on the biochemical profile of 717 individuals including 213 children with ASD, 128 unaffected siblings and 376 parents and other relatives, we estimated the heritability of whole-blood serotonin, platelet N-acetylserotonin (NAS) and plasma melatonin levels, as well as the two enzymes arylalkylamine N-acetyltransferase (AANAT) and acetylserotonin O-methyltransferase (ASMT) activities measured in platelets. Overall, heritability was higher for NAS (0.72 ± 0.091) and ASMT (0.59 ± 0.097) compared with serotonin (0.31 ± 0.078), AANAT (0.34 ± 0.077) and melatonin (0.22 ± 0.071). Bivariate analyses showed high phenotypic and genetic correlations between traits of the second step of the metabolic pathway (NAS, ASMT and melatonin) indicating the contribution of shared genetic factors. A better knowledge of the heritability of the melatonin synthesis variability constitutes an important step to identify the factors that perturb this pathway in individuals with ASD.

Mosaicism and prenatal diagnosis options: insights from retinoblastoma.

In sporadic cases, a post-zygotic mutational event signifies a somatic mosaicism in the affected child only, which implies that these mutations affect only a portion of the body. Therefore siblings do not need follow-up. On the other hand, a pre-zygotic mutation transmitted by an unaffected mosaic parent implies recurrent risks in offspring. To better estimate the contribution of pre- and post-zygotic events, we analysed 124 consecutive bilateral retinoblastoma probands, carrying a heterozygous RB1 pathogenic variant and their unaffected, non-carrier parents. In order to evaluate somatic mosaicism in blood, the deleterious RB1 pathogenic variant identified in the proband, was searched for in the unaffected parents, using targeted deep sequencing. Observed recurrences, which represent an estimation of germline and somatic mosaicisms, were recorded and computed in the sibships. Deep sequencing revealed one mosaic-unaffected parent out of 124 tested couples, which provides an estimation of the maximal risk of recurrence, due to parental mosaicism, at 0.4%. Follow-up in the sibships showed one recurrence, providing a maximal recurrence risk, due to parental mosaicism, at 0.8%. Two different statistical strategies led to close estimates (0.4 and 0.8% risks) which appeared 266-533-fold higher, as compared with the general population. These recurrence estimates could be considered when counselling couples with retinoblastoma or diseases with a high de novo mutation rate.

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.

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.