Niche and local geography shape the pangenome of wastewater- and livestock-associated Enterobacteriaceae.

Escherichia coli and other Enterobacteriaceae are diverse species with "open" pangenomes, where genes move intra- and interspecies via horizontal gene transfer. However, most analyses focus on clinical isolates. The pangenome dynamics of natural populations remain understudied, despite their suggested role as reservoirs for antimicrobial resistance (AMR) genes. Here, we analyze near-complete genomes for 827 Enterobacteriaceae (553 Escherichia and 274 non-Escherichia spp.) with 2292 circularized plasmids in total, collected from 19 locations (livestock farms and wastewater treatment works in the United Kingdom) within a 30-km radius at three time points over a year. We find different dynamics for chromosomal and plasmid-borne genes. Plasmids have a higher burden of AMR genes and insertion sequences, and AMR-gene-carrying plasmids show evidence of being under stronger selective pressure. Environmental niche and local geography both play a role in shaping plasmid dynamics. Our results highlight the importance of local strategies for controlling the spread of AMR.

Genomic network analysis of environmental and livestock F-type plasmid populations.

F-type plasmids are diverse and of great clinical significance, often carrying genes conferring antimicrobial resistance (AMR) such as extended-spectrum ß-lactamases, particularly in Enterobacterales. Organising this plasmid diversity is challenging, and current knowledge is largely based on plasmids from clinical settings. Here, we present a network community analysis of a large survey of F-type plasmids from environmental (influent, effluent and upstream/downstream waterways surrounding wastewater treatment works) and livestock settings. We use a tractable and scalable methodology to examine the relationship between plasmid metadata and network communities. This reveals how niche (sampling compartment and host genera) partition and shape plasmid diversity. We also perform pangenome-style analyses on network communities. We show that such communities define unique combinations of core genes, with limited overlap. Building plasmid phylogenies based on alignments of these core genes, we demonstrate that plasmid accessory function is closely linked to core gene content. Taken together, our results suggest that stable F-type plasmid backbone structures can persist in environmental settings while allowing dramatic variation in accessory gene content that may be linked to niche adaptation. The association of F-type plasmids with AMR may reflect their suitability for rapid niche adaptation.

DNA Thermo-Protection Facilitates Whole-Genome Sequencing of Mycobacteria Direct from Clinical Samples.

Mycobacterium tuberculosis is the leading cause of death from bacterial infection. Improved rapid diagnosis and antimicrobial resistance determination, such as by whole-genome sequencing, are required. Our aim was to develop a simple, low-cost method of preparing DNA for sequencing direct from M. tuberculosis-positive clinical samples (without culture). Simultaneous sputum liquefaction, bacteria heat inactivation (99°C/30 min), and enrichment for mycobacteria DNA were achieved using an equal volume of thermo-protection buffer (4 M KCl, 0.05 M HEPES buffer, pH 7.5, 0.1% dithiothreitol [DTT]). The buffer emulated intracellular conditions found in hyperthermophiles, thus protecting DNA from rapid thermodegradation, which renders it a poor template for sequencing. Initial validation experiments employed mycobacteria DNA, either extracted or intracellular. Next, mock clinical samples (infection-negative human sputum spiked with 0 to 105Mycobacterium bovis BCG cells/ml) underwent liquefaction in thermo-protection buffer and heat inactivation. DNA was extracted and sequenced. Human DNA degraded faster than mycobacteria DNA, resulting in target enrichment. Four replicate experiments achieved M. tuberculosis detection at 101 BCG cells/ml, with 31 to 59 M. tuberculosis complex reads. Maximal genome coverage (>97% at 5× depth) occurred at 104 BCG cells/ml; >91% coverage (1× depth) occurred at 103 BCG cells/ml. Final validation employed M. tuberculosis-positive clinical samples (n = 20), revealing that initial sample volumes of ≥1 ml typically yielded higher mean depths of M. tuberculosis genome coverage, with an overall range of 0.55 to 81.02. A mean depth of 3 gave >96% 1-fold tuberculosis (TB) genome coverage (in 15/20 clinical samples). A mean depth of 15 achieved >99% 5-fold genome coverage (in 9/20 clinical samples). In summary, direct-from-sample sequencing of M. tuberculosis genomes was facilitated by a low-cost thermo-protection buffer.

Antibody testing for COVID-19: A report from the National COVID Scientific Advisory Panel.

Background: The COVID-19 pandemic caused >1 million infections during January-March 2020. There is an urgent need for reliable antibody detection approaches to support diagnosis, vaccine development, safe release of individuals from quarantine, and population lock-down exit strategies. We set out to evaluate the performance of ELISA and lateral flow immunoassay (LFIA) devices. Methods: We tested plasma for COVID (severe acute respiratory syndrome coronavirus 2; SARS-CoV-2) IgM and IgG antibodies by ELISA and using nine different LFIA devices. We used a panel of plasma samples from individuals who have had confirmed COVID infection based on a PCR result (n=40), and pre-pandemic negative control samples banked in the UK prior to December-2019 (n=142). Results: ELISA detected IgM or IgG in 34/40 individuals with a confirmed history of COVID infection (sensitivity 85%, 95%CI 70-94%), vs. 0/50 pre-pandemic controls (specificity 100% [95%CI 93-100%]). IgG levels were detected in 31/31 COVID-positive individuals tested ≥10 days after symptom onset (sensitivity 100%, 95%CI 89-100%). IgG titres rose during the 3 weeks post symptom onset and began to fall by 8 weeks, but remained above the detection threshold. Point estimates for the sensitivity of LFIA devices ranged from 55-70% versus RT-PCR and 65-85% versus ELISA, with specificity 95-100% and 93-100% respectively. Within the limits of the study size, the performance of most LFIA devices was similar. Conclusions: Currently available commercial LFIA devices do not perform sufficiently well for individual patient applications. However, ELISA can be calibrated to be specific for detecting and quantifying SARS-CoV-2 IgM and IgG and is highly sensitive for IgG from 10 days following first symptoms.

Comparison of long-read sequencing technologies in the hybrid assembly of complex bacterial genomes.

Illumina sequencing allows rapid, cheap and accurate whole genome bacterial analyses, but short reads (<300 bp) do not usually enable complete genome assembly. Long-read sequencing greatly assists with resolving complex bacterial genomes, particularly when combined with short-read Illumina data (hybrid assembly). However, it is not clear how different long-read sequencing methods affect hybrid assembly accuracy. Relative automation of the assembly process is also crucial to facilitating high-throughput complete bacterial genome reconstruction, avoiding multiple bespoke filtering and data manipulation steps. In this study, we compared hybrid assemblies for 20 bacterial isolates, including two reference strains, using Illumina sequencing and long reads from either Oxford Nanopore Technologies (ONT) or SMRT Pacific Biosciences (PacBio) sequencing platforms. We chose isolates from the family Enterobacteriaceae, as these frequently have highly plastic, repetitive genetic structures, and complete genome reconstruction for these species is relevant for a precise understanding of the epidemiology of antimicrobial resistance. We de novo assembled genomes using the hybrid assembler Unicycler and compared different read processing strategies, as well as comparing to long-read-only assembly with Flye followed by short-read polishing with Pilon. Hybrid assembly with either PacBio or ONT reads facilitated high-quality genome reconstruction, and was superior to the long-read assembly and polishing approach evaluated with respect to accuracy and completeness. Combining ONT and Illumina reads fully resolved most genomes without additional manual steps, and at a lower consumables cost per isolate in our setting. Automated hybrid assembly is a powerful tool for complete and accurate bacterial genome assembly.

The impact of sequencing depth on the inferred taxonomic composition and AMR gene content of metagenomic samples.

BACKGROUND: Shotgun metagenomics is increasingly used to characterise microbial communities, particularly for the investigation of antimicrobial resistance (AMR) in different animal and environmental contexts. There are many different approaches for inferring the taxonomic composition and AMR gene content of complex community samples from shotgun metagenomic data, but there has been little work establishing the optimum sequencing depth, data processing and analysis methods for these samples. In this study we used shotgun metagenomics and sequencing of cultured isolates from the same samples to address these issues. We sampled three potential environmental AMR gene reservoirs (pig caeca, river sediment, effluent) and sequenced samples with shotgun metagenomics at high depth (~ 200 million reads per sample). Alongside this, we cultured single-colony isolates of Enterobacteriaceae from the same samples and used hybrid sequencing (short- and long-reads) to create high-quality assemblies for comparison to the metagenomic data. To automate data processing, we developed an open-source software pipeline, 'ResPipe'. RESULTS: Taxonomic profiling was much more stable to sequencing depth than AMR gene content. 1 million reads per sample was sufficient to achieve < 1% dissimilarity to the full taxonomic composition. However, at least 80 million reads per sample were required to recover the full richness of different AMR gene families present in the sample, and additional allelic diversity of AMR genes was still being discovered in effluent at 200 million reads per sample. Normalising the number of reads mapping to AMR genes using gene length and an exogenous spike of Thermus thermophilus DNA substantially changed the estimated gene abundance distributions. While the majority of genomic content from cultured isolates from effluent was recoverable using shotgun metagenomics, this was not the case for pig caeca or river sediment. CONCLUSIONS: Sequencing depth and profiling method can critically affect the profiling of polymicrobial animal and environmental samples with shotgun metagenomics. Both sequencing of cultured isolates and shotgun metagenomics can recover substantial diversity that is not identified using the other methods. Particular consideration is required when inferring AMR gene content or presence by mapping metagenomic reads to a database. ResPipe, the open-source software pipeline we have developed, is freely available ( https://gitlab.com/hsgweon/ResPipe ).

Automated detection of bacterial growth on 96-well plates for high-throughput drug susceptibility testing of Mycobacterium tuberculosis.

M. tuberculosis grows slowly and is challenging to work with experimentally compared with many other bacteria. Although microtitre plates have the potential to enable high-throughput phenotypic testing of M. tuberculosis, they can be difficult to read and interpret. Here we present a software package, the Automated Mycobacterial Growth Detection Algorithm (AMyGDA), that measures how much M. tuberculosis is growing in each well of a 96-well microtitre plate. The plate used here has serial dilutions of 14 anti-tuberculosis drugs, thereby permitting the MICs to be elucidated. The three participating laboratories each inoculated 38 96-well plates with 15 known M. tuberculosis strains (including the standard H37Rv reference strain) and, after 2 weeks' incubation, measured the MICs for all 14 drugs on each plate and took a photograph. By analysing the images, we demonstrate that AMyGDA is reproducible, and that the MICs measured are comparable to those measured by a laboratory scientist. The AMyGDA software will be used by the Comprehensive Resistance Prediction for Tuberculosis: an International Consortium (CRyPTIC) to measure the drug susceptibility profile of a large number (>30000) of samples of M. tuberculosis from patients over the next few years.

Validating a 14-Drug Microtiter Plate Containing Bedaquiline and Delamanid for Large-Scale Research Susceptibility Testing of Mycobacterium tuberculosis.

The UKMYC5 plate is a 96-well microtiter plate designed by the CRyPTIC Consortium (Comprehensive Resistance Prediction for Tuberculosis: an International Consortium) to enable the measurement of MICs of 14 different antituberculosis (anti-TB) compounds for >30,000 clinical Mycobacterium tuberculosis isolates. Unlike the MYCOTB plate, on which the UKMYC5 plate is based, the UKMYC5 plate includes two new (bedaquiline and delamanid) and two repurposed (clofazimine and linezolid) compounds. UKMYC5 plates were tested by seven laboratories on four continents by use of a panel of 19 external quality assessment (EQA) strains, including H37Rv. To assess the optimal combination of reading method and incubation time, MICs were measured from each plate by two readers, using three methods (mirrored box, microscope, and Vizion digital viewing system), after 7, 10, 14, and 21 days of incubation. In addition, all EQA strains were subjected to whole-genome sequencing and phenotypically characterized by the 7H10/7H11 agar proportion method (APM) and by use of MGIT960 mycobacterial growth indicator tubes. We concluded that the UKMYC5 plate is optimally read using the Vizion system after 14 days of incubation, achieving an interreader agreement of 97.9% and intra- and interlaboratory reproducibility rates of 95.6% and 93.1%, respectively. The mirrored box had a similar reproducibility. Strains classified as resistant by APM, MGIT960, or the presence of mutations known to confer resistance consistently showed elevated MICs compared to those for strains classified as susceptible. Finally, the UKMYC5 plate records intermediate MICs for one strain for which the APM measured MICs close to the applied critical concentration, providing early evidence that the UKMYC5 plate can quantitatively measure the magnitude of resistance to anti-TB compounds that is due to specific genetic variation.

Structural basis for LMO2-driven recruitment of the SCL:E47bHLH heterodimer to hematopoietic-specific transcriptional targets.

Cell fate is governed by combinatorial actions of transcriptional regulators assembling into multiprotein complexes. However, the molecular details of how these complexes form are poorly understood. One such complex, which contains the basic-helix-loop-helix heterodimer SCL:E47 and bridging proteins LMO2:LDB1, critically regulates hematopoiesis and induces T cell leukemia. Here, we report the crystal structure of (SCL:E47)bHLH:LMO2:LDB1LID bound to DNA, providing a molecular account of the network of interactions assembling this complex. This reveals an unexpected role for LMO2. Upon binding to SCL, LMO2 induces new hydrogen bonds in SCL:E47, thereby strengthening heterodimer formation. This imposes a rotation movement onto E47 that weakens the heterodimer:DNA interaction, shifting the main DNA-binding activity onto additional protein partners. Along with biochemical analyses, this illustrates, at an atomic level, how hematopoietic-specific SCL sequesters ubiquitous E47 and associated cofactors and supports SCL's reported DNA-binding-independent functions. Importantly, this work will drive the design of small molecules inhibiting leukemogenic processes.

Complementary functions of the flippase ATP8B1 and the floppase ABCB4 in maintaining canalicular membrane integrity.

BACKGROUND & AIMS: Progressive familial intrahepatic cholestasis can be caused by mutations in ABCB4 or ATP8B1; each encodes a protein that translocates phospholipids, but in opposite directions. ABCB4 flops phosphatidylcholine from the inner to the outer leaflet, where it is extracted by bile salts. ATP8B1, in complex with the accessory protein CDC50A, flips phosphatidylserine in the reverse direction. Abcb4(-/-) mice lack biliary secretion of phosphatidylcholine, whereas Atp8b1-deficient mice have increased excretion of phosphatidylserine into bile. Each system is thought to have a role protecting the canalicular membrane from bile salts. METHODS: To investigate the relationship between the mechanisms of ABCB4 and ATP8B1, we expressed the transporters separately and together in cultured cells and studied viability and phospholipid transport. We also created mice with disruptions in ABCB4 and ATP8B1 (double knockouts) and studied bile formation and hepatic damage in mice fed bile salts. RESULTS: Overexpression of ABCB4 was toxic to HEK293T cells; the toxicity was counteracted by coexpression of the ATP8B1-CDC50A complex. In Atp8b1-deficient mice, bile salts induced extraction of phosphatidylserine and ectoenzymes from the canalicular membrane; this process was not observed in the double-knockout mice. CONCLUSIONS: ATP8B1 is required for hepatocyte function, particularly in the presence of ABCB4. This is most likely because the phosphatidylserine flippase complex of ATP8B1-CDC50A counteracts the destabilization of the membrane that occurs when ABCB4 flops phosphatidylcholine. Lipid asymmetry is therefore important for the integrity of the canalicular membrane; ABCB4 and ATP8B1 cooperate to protect hepatocytes from bile salts.

Structure of the leukemia oncogene LMO2: implications for the assembly of a hematopoietic transcription factor complex.

The LIM only protein 2 (LMO2) is a key regulator of hematopoietic stem cell development whose ectopic expression in T cells leads to the onset of acute lymphoblastic leukemia. Through its LIM domains, LMO2 is thought to function as the scaffold for a DNA-binding transcription regulator complex, including the basic helix-loop-helix proteins SCL/TAL1 and E47, the zinc finger protein GATA-1, and LIM-domain interacting protein LDB1. To understand the role of LMO2 in the formation of this complex and ultimately to dissect its function in normal and aberrant hematopoiesis, we solved the crystal structure of LMO2 in complex with the LID domain of LDB1 at 2.4 Å resolution. We observe a largely unstructured LMO2 kept in register by the LID binding both LIM domains. Comparison of independently determined crystal structures of LMO2 reveals large movements around a conserved hinge between the LIM domains. We demonstrate that such conformational flexibility is necessary for binding of LMO2 to its partner protein SCL/TAL1 in vitro and for the function of this complex in vivo. These results, together with molecular docking and analysis of evolutionarily conserved residues, yield the first structural model of the DNA-binding complex containing LMO2, LDB1, SCL/TAL1, and GATA-1.

The Human Scavenger Receptor CD36: glycosylation status and its role in trafficking and function.

Human CD36 is a class B scavenger receptor expressed in a variety of cell types such as macrophage and adipocytes. This plasma membrane glycoprotein has a wide range of ligands including oxidized low density lipoprotein and long chain fatty acids which involves the receptor in diseases such as atherosclerosis and insulin resistance. CD36 is heavily modified post-translationally by N-linked glycosylation, and 10 putative glycosylation sites situated in the large extracellular loop of the protein have been identified; however, their utilization and role in the folding and function of the protein have not been characterized. Using mass spectrometry on purified and peptide N-glycosidase F-deglycosylated CD36 and also by comparing the electrophoretic mobility of different glycosylation site mutants, we have determined that 9 of the 10 sites can be modified by glycosylation. Flow cytometric analysis of the different glycosylation mutants expressed in mammalian cells established that glycosylation is necessary for trafficking to the plasma membrane. Minimally glycosylated mutants that supported trafficking were identified and indicated the importance of carboxyl-terminal sites Asn-247, Asn-321, and Asn-417. However, unlike SRBI, no individual site was found to be essential for proper trafficking of CD36. Surprisingly, these minimally glycosylated mutants appear to be predominantly core-glycosylated, indicating that mature glycosylation is not necessary for surface expression in mammalian cells. The data also show that neither the nature nor the pattern of glycosylation is relevant to binding of modified low density lipoprotein.

Cd36, a class B scavenger receptor, functions as a monomer to bind acetylated and oxidized low-density lipoproteins.

Cd36 is a small-molecular-weight integral membrane protein expressed in a diverse, but select, range of cell types. It has an equally diverse range of ligands and physiological functions, which has implicated Cd36 in a number of diseases including insulin resistance, diabetes, and, most notably, atherosclerosis. The protein is reported to reside in detergent-resistant microdomains within the plasma membrane and to form homo- and hetero-intermolecular interactions. These data suggest that this class B scavenger receptor may gain functionality for ligand binding, and/or ligand internalization, by formation of protein complexes at the cell surface. Here, we have overexpressed Cd36 in insect cells, purified the recombinant protein to homogeneity, and analyzed its stability and solubility in a variety of nonionic and zwitterionic detergents. Octylglucoside conferred the greatest degree of stability, and by analytical ultracentrifugation we show that the protein is monomeric. A solid-phase ligand-binding assay demonstrated that the purified monomeric protein retains high affinity for acetylated and oxidized low-density lipoproteins. Therefore, no accessory proteins are required for interaction with ligand, and binding is a property of the monomeric fold of the protein. Thus, the highly purified and functional Cd36 should be suitable for crystallization in octylglucoside, and the in vitro ligand-binding assay represents a promising screen for identification of bioactive molecules targeting atherogenesis at the level of ligand binding.