Bacterial droplet-based single-cell RNA-seq reveals antibiotic-associated heterogeneous cellular states.

We introduce BacDrop, a highly scalable technology for bacterial single-cell RNA sequencing that has overcome many challenges hindering the development of scRNA-seq in bacteria. BacDrop can be applied to thousands to millions of cells from both gram-negative and gram-positive species. It features universal ribosomal RNA depletion and combinatorial barcodes that enable multiplexing and massively parallel sequencing. We applied BacDrop to study Klebsiella pneumoniae clinical isolates and to elucidate their heterogeneous responses to antibiotic stress. In an unperturbed population presumed to be homogeneous, we found within-population heterogeneity largely driven by the expression of mobile genetic elements that promote the evolution of antibiotic resistance. Under antibiotic perturbation, BacDrop revealed transcriptionally distinct subpopulations associated with different phenotypic outcomes including antibiotic persistence. BacDrop thus can capture cellular states that cannot be detected by bulk RNA-seq, which will unlock new microbiological insights into bacterial responses to perturbations and larger bacterial communities such as the microbiome.

An Aqueous Analog MAC Machine.

Using ions in aqueous milieu for signal processing, like in biological circuits, may potentially lead to a bioinspired information processing platform. Studies, however, have focused on individual ionic diodes and transistors rather than circuits comprising many such devices. Here a 16 × 16 array of new ionic transistors is developed in an aqueous quinone solution. Each transistor features a concentric ring electrode pair with a disk electrode at the center. The electrochemistry of these electrodes in the solution provides the basis for the transistor operation. The ring pair electrochemically tunes the local electrolytic concentration to modulate the disk's Faradaic reaction rate. Thus, the disk current as a Faradaic reaction to the disk voltage is gated by the ring pair. The 16 × 16 array of these transistors performs analog multiply-accumulate (MAC) operations, a computing modality hotly pursued for low-power artificial neural networks. This exploits the transistor's operating regime where the disk current is a multiplication of the disk voltage and a weight parameter tuned by the ring pair gating. Such disk currents from multiple transistors are summated in a global reference electrode to complete a MAC task. This ionic circuit demonstrating analog computing is a step toward sophisticated aqueous ionics.

CMOS electrochemical pH localizer-imager.

pH controls a large repertoire of chemical and biochemical processes in water. Densely arrayed pH microenvironments would parallelize these processes, enabling their high-throughput studies and applications. However, pH localization, let alone its arrayed realization, remains challenging because of fast diffusion of protons in water. Here, we demonstrate arrayed localizations of picoliter-scale aqueous acids, using a 256-electrochemical cell array defined on and operated by a complementary metal oxide semiconductor (CMOS)-integrated circuit. Each cell, comprising a concentric pair of cathode and anode with their current injections controlled with a sub-nanoampere resolution by the CMOS electronics, creates a local pH environment, or a pH "voxel," via confined electrochemistry. The system also monitors the spatiotemporal pH profile across the array in real time for precision pH control. We highlight the utility of this CMOS pH localizer-imager for high-throughput tasks by parallelizing pH-gated molecular state encoding and pH-regulated enzymatic DNA elongation at any selected set of cells.

High-Throughput Regulatory Part Prototyping and Analysis by Cell-Free Protein Synthesis and Droplet Microfluidics.

Engineering regulatory parts for improved performance in genetic programs has played a pivotal role in the development of the synthetic biology cell programming toolbox. Here, we report the development of a novel high-throughput platform for regulatory part prototyping and analysis that leverages the advantages of engineered DNA libraries, cell-free protein synthesis (CFPS), high-throughput emulsion droplet microfluidics, standard flow sorting adapted to screen droplet reactions, and next-generation sequencing (NGS). With this integrated platform, we screened the activity of millions of genetic parts within hours, followed by NGS retrieval of the improved designs. This in vitro platform is particularly valuable for engineering regulatory parts of nonmodel organisms, where in vivo high-throughput screening methods are not readily available. The platform can be extended to multipart screening of complete genetic programs to optimize yield and stability.

Increasing cell-free gene expression yields from linear templates in Escherichia coli and Vibrio natriegens extracts by using DNA-binding proteins.

In crude extract-based cell-free protein synthesis (CFPS), DNA templates are transcribed and translated into functional proteins. Although linear expression templates (LETs) are less laborious and expensive to generate, plasmid templates are often desired over polymerase chain reaction-generated LETs due to increased stability and protection against exonucleases present in the extract of the reaction. Here we demonstrate that addition of a double stranded DNA-binding protein to the CFPS reaction, termed single-chain Cro protein (scCro), achieves terminal protection of LETs. This CroP-LET (scCro-based protection of LET) method effectively increases superfolder green fluorescent protein (sfGFP) expression levels from LETs in Escherichia coli CFPS reactions by sixfold. Our yields are comparable to other strategies that provide chemical and enzymatic DNA stabilization in E. coli CFPS. Notably, we also report that the CroP-LET method successfully enhanced yields in CFPS platforms derived from nonmodel organisms. Our results show that CroP-LET increased sfGFP yields by 18-fold in the Vibrio natriegens CFPS platform. With the fast-expanding applications of CFPS platforms, this method provides a practical and generalizable solution to protect linear expression DNA templates.

Registry in a tube: multiplexed pools of retrievable parts for genetic design space exploration.

Genetic designs can consist of dozens of genes and hundreds of genetic parts. After evaluating a design, it is desirable to implement changes without the cost and burden of starting the construction process from scratch. Here, we report a two-step process where a large design space is divided into deep pools of composite parts, from which individuals are retrieved and assembled to build a final construct. The pools are built via multiplexed assembly and sequenced using next-generation sequencing. Each pool consists of ∼20 Mb of up to 5000 unique and sequence-verified composite parts that are barcoded for retrieval by PCR. This approach is applied to a 16-gene nitrogen fixation pathway, which is broken into pools containing a total of 55 848 composite parts (71.0 Mb). The pools encompass an enormous design space (1043 possible 23 kb constructs), from which an algorithm-guided 192-member 4.5 Mb library is built. Next, all 1030 possible genetic circuits based on 10 repressors (NOR/NOT gates) are encoded in pools where each repressor is fused to all permutations of input promoters. These demonstrate that multiplexing can be applied to encompass entire design spaces from which individuals can be accessed and evaluated.

Functional optimization of gene clusters by combinatorial design and assembly.

Large microbial gene clusters encode useful functions, including energy utilization and natural product biosynthesis, but genetic manipulation of such systems is slow, difficult and complicated by complex regulation. We exploit the modularity of a refactored Klebsiella oxytoca nitrogen fixation (nif) gene cluster (16 genes, 103 parts) to build genetic permutations that could not be achieved by starting from the wild-type cluster. Constraint-based combinatorial design and DNA assembly are used to build libraries of radically different cluster architectures by varying part choice, gene order, gene orientation and operon occupancy. We construct 84 variants of the nifUSVWZM operon, 145 variants of the nifHDKY operon, 155 variants of the nifHDKYENJ operon and 122 variants of the complete 16-gene pathway. The performance and behavior of these variants are characterized by nitrogenase assay and strand-specific RNA sequencing (RNA-seq), and the results are incorporated into subsequent design cycles. We have produced a fully synthetic cluster that recovers 57% of wild-type activity. Our approach allows the performance of genetic parts to be quantified simultaneously in hundreds of genetic contexts. This parallelized design-build-test-learn cycle, which can access previously unattainable regions of genetic space, should provide a useful, fast tool for genetic optimization and hypothesis testing.

Staphylococcal enterotoxin P predicts bacteremia in hospitalized patients colonized with methicillin-resistant Staphylococcus aureus.

BACKGROUND: Methicillin-resistant Staphylococcus aureus (MRSA) colonization predicts later infection, with both host and pathogen determinants of invasive disease. METHODS: This nested case-control study evaluates predictors of MRSA bacteremia in an 8-intensive care unit (ICU) prospective adult cohort from 1 September 2003 through 30 April 2005 with active MRSA surveillance and collection of ICU, post-ICU, and readmission MRSA isolates. We selected MRSA carriers who did (cases) and those who did not (controls) develop MRSA bacteremia. Generating assembled genome sequences, we evaluated 30 MRSA genes potentially associated with virulence and invasion. Using multivariable Cox proportional hazards regression, we assessed the association of these genes with MRSA bacteremia, controlling for host risk factors. RESULTS: We collected 1578 MRSA isolates from 520 patients. We analyzed host and pathogen factors for 33 cases and 121 controls. Predictors of MRSA bacteremia included a diagnosis of cancer, presence of a central venous catheter, hyperglycemia (glucose level, >200 mg/dL), and infection with a MRSA strain carrying the gene for staphylococcal enterotoxin P (sep). Receipt of an anti-MRSA medication had a significant protective effect. CONCLUSIONS: In an analysis controlling for host factors, colonization with MRSA carrying sep increased the risk of MRSA bacteremia. Identification of risk-adjusted genetic determinants of virulence may help to improve prediction of invasive disease and suggest new targets for therapeutic intervention.

The genomic basis of adaptive evolution in threespine sticklebacks.

Marine stickleback fish have colonized and adapted to thousands of streams and lakes formed since the last ice age, providing an exceptional opportunity to characterize genomic mechanisms underlying repeated ecological adaptation in nature. Here we develop a high-quality reference genome assembly for threespine sticklebacks. By sequencing the genomes of twenty additional individuals from a global set of marine and freshwater populations, we identify a genome-wide set of loci that are consistently associated with marine-freshwater divergence. Our results indicate that reuse of globally shared standing genetic variation, including chromosomal inversions, has an important role in repeated evolution of distinct marine and freshwater sticklebacks, and in the maintenance of divergent ecotypes during early stages of reproductive isolation. Both coding and regulatory changes occur in the set of loci underlying marine-freshwater evolution, but regulatory changes appear to predominate in this well known example of repeated adaptive evolution in nature.

A high-resolution map of human evolutionary constraint using 29 mammals.

The comparison of related genomes has emerged as a powerful lens for genome interpretation. Here we report the sequencing and comparative analysis of 29 eutherian genomes. We confirm that at least 5.5% of the human genome has undergone purifying selection, and locate constrained elements covering ∼4.2% of the genome. We use evolutionary signatures and comparisons with experimental data sets to suggest candidate functions for ∼60% of constrained bases. These elements reveal a small number of new coding exons, candidate stop codon readthrough events and over 10,000 regions of overlapping synonymous constraint within protein-coding exons. We find 220 candidate RNA structural families, and nearly a million elements overlapping potential promoter, enhancer and insulator regions. We report specific amino acid residues that have undergone positive selection, 280,000 non-coding elements exapted from mobile elements and more than 1,000 primate- and human-accelerated elements. Overlap with disease-associated variants indicates that our findings will be relevant for studies of human biology, health and disease.

A scalable, fully automated process for construction of sequence-ready human exome targeted capture libraries.

Genome targeting methods enable cost-effective capture of specific subsets of the genome for sequencing. We present here an automated, highly scalable method for carrying out the Solution Hybrid Selection capture approach that provides a dramatic increase in scale and throughput of sequence-ready libraries produced. Significant process improvements and a series of in-process quality control checkpoints are also added. These process improvements can also be used in a manual version of the protocol.

High-quality draft assemblies of mammalian genomes from massively parallel sequence data.

Massively parallel DNA sequencing technologies are revolutionizing genomics by making it possible to generate billions of relatively short (~100-base) sequence reads at very low cost. Whereas such data can be readily used for a wide range of biomedical applications, it has proven difficult to use them to generate high-quality de novo genome assemblies of large, repeat-rich vertebrate genomes. To date, the genome assemblies generated from such data have fallen far short of those obtained with the older (but much more expensive) capillary-based sequencing approach. Here, we report the development of an algorithm for genome assembly, ALLPATHS-LG, and its application to massively parallel DNA sequence data from the human and mouse genomes, generated on the Illumina platform. The resulting draft genome assemblies have good accuracy, short-range contiguity, long-range connectivity, and coverage of the genome. In particular, the base accuracy is high (≥99.95%) and the scaffold sizes (N50 size = 11.5 Mb for human and 7.2 Mb for mouse) approach those obtained with capillary-based sequencing. The combination of improved sequencing technology and improved computational methods should now make it possible to increase dramatically the de novo sequencing of large genomes. The ALLPATHS-LG program is available at http://www.broadinstitute.org/science/programs/genome-biology/crd.

A scalable, fully automated process for construction of sequence-ready barcoded libraries for 454.

We present an automated, high throughput library construction process for 454 technology. Sample handling errors and cross-contamination are minimized via end-to-end barcoding of plasticware, along with molecular DNA barcoding of constructs. Automation-friendly magnetic bead-based size selection and cleanup steps have been devised, eliminating major bottlenecks and significant sources of error. Using this methodology, one technician can create 96 sequence-ready 454 libraries in 2 days, a dramatic improvement over the standard method.

Human chromosome 11 DNA sequence and analysis including novel gene identification.

Chromosome 11, although average in size, is one of the most gene- and disease-rich chromosomes in the human genome. Initial gene annotation indicates an average gene density of 11.6 genes per megabase, including 1,524 protein-coding genes, some of which were identified using novel methods, and 765 pseudogenes. One-quarter of the protein-coding genes shows overlap with other genes. Of the 856 olfactory receptor genes in the human genome, more than 40% are located in 28 single- and multi-gene clusters along this chromosome. Out of the 171 disorders currently attributed to the chromosome, 86 remain for which the underlying molecular basis is not yet known, including several mendelian traits, cancer and susceptibility loci. The high-quality data presented here--nearly 134.5 million base pairs representing 99.8% coverage of the euchromatic sequence--provide scientists with a solid foundation for understanding the genetic basis of these disorders and other biological phenomena.

DNA sequence and analysis of human chromosome 8.

The International Human Genome Sequencing Consortium (IHGSC) recently completed a sequence of the human genome. As part of this project, we have focused on chromosome 8. Although some chromosomes exhibit extreme characteristics in terms of length, gene content, repeat content and fraction segmentally duplicated, chromosome 8 is distinctly typical in character, being very close to the genome median in each of these aspects. This work describes a finished sequence and gene catalogue for the chromosome, which represents just over 5% of the euchromatic human genome. A unique feature of the chromosome is a vast region of approximately 15 megabases on distal 8p that appears to have a strikingly high mutation rate, which has accelerated in the hominids relative to other sequenced mammals. This fast-evolving region contains a number of genes related to innate immunity and the nervous system, including loci that appear to be under positive selection--these include the major defensin (DEF) gene cluster and MCPH1, a gene that may have contributed to the evolution of expanded brain size in the great apes. The data from chromosome 8 should allow a better understanding of both normal and disease biology and genome evolution.

DNA sequence and analysis of human chromosome 18.

Chromosome 18 appears to have the lowest gene density of any human chromosome and is one of only three chromosomes for which trisomic individuals survive to term. There are also a number of genetic disorders stemming from chromosome 18 trisomy and aneuploidy. Here we report the finished sequence and gene annotation of human chromosome 18, which will allow a better understanding of the normal and disease biology of this chromosome. Despite the low density of protein-coding genes on chromosome 18, we find that the proportion of non-protein-coding sequences evolutionarily conserved among mammals is close to the genome-wide average. Extending this analysis to the entire human genome, we find that the density of conserved non-protein-coding sequences is largely uncorrelated with gene density. This has important implications for the nature and roles of non-protein-coding sequence elements.

The genome sequence of the rice blast fungus Magnaporthe grisea.

Magnaporthe grisea is the most destructive pathogen of rice worldwide and the principal model organism for elucidating the molecular basis of fungal disease of plants. Here, we report the draft sequence of the M. grisea genome. Analysis of the gene set provides an insight into the adaptations required by a fungus to cause disease. The genome encodes a large and diverse set of secreted proteins, including those defined by unusual carbohydrate-binding domains. This fungus also possesses an expanded family of G-protein-coupled receptors, several new virulence-associated genes and large suites of enzymes involved in secondary metabolism. Consistent with a role in fungal pathogenesis, the expression of several of these genes is upregulated during the early stages of infection-related development. The M. grisea genome has been subject to invasion and proliferation of active transposable elements, reflecting the clonal nature of this fungus imposed by widespread rice cultivation.

Effectiveness of health care worker training on the oral health of elderly residents of nursing homes.

OBJECTIVES: The aim of this study was to evaluate the effect of a staff training programme on mouth care on the oral health of elderly residents of long-term care institutions. METHODS: Seventy-eight residents of five long-stay institutions were enrolled and underwent a baseline oral health assessment. Staff caring for residents in three of the sites received intensive training in mouth care. This comprised lecture and video material complemented by clinical demonstrations. The oral health of residents at all five sites was reassessed at periods of 3 and 9 months. Staff caring for residents in the remaining two institutions were then provided with mouth care training and all patients were reassessed at 18 months. Statistical analyses were undertaken to examine for significant changes in selected oral health parameters after training, within each group. RESULTS: Oral mucosal disease and oral dryness were common at baseline. The staff training was well received. Following staff training, there was a significant reduction in the number of residents left to undertake their own oral care. There were significant improvements in denture hygiene and a reduction in the number of residents wearing dentures overnight. The prevalence of oral mucosal disease dropped, with significant reductions in angular cheilitis and denture stomatitis. CONCLUSION: This education programme was effective in changing oral health care procedures within long-stay institutions for the elderly, with measurable improvements in oral health of the residents.

The complete genome and proteome of Mycoplasma mobile.

Although often considered "minimal" organisms, mycoplasmas show a wide range of diversity with respect to host environment, phenotypic traits, and pathogenicity. Here we report the complete genomic sequence and proteogenomic map for the piscine mycoplasma Mycoplasma mobile, noted for its robust gliding motility. For the first time, proteomic data are used in the primary annotation of a new genome, providing validation of expression for many of the predicted proteins. Several novel features were discovered including a long repeating unit of DNA of approximately 2435 bp present in five complete copies that are shown to code for nearly identical yet uniquely expressed proteins. M. mobile has among the lowest DNA GC contents (24.9%) and most reduced set of tRNAs of any organism yet reported (28). Numerous instances of tandem duplication as well as lateral gene transfer are evident in the genome. The multiple available complete genome sequences for other motile and immotile mycoplasmas enabled us to use comparative genomic and phylogenetic methods to suggest several candidate genes that might be involved in motility. The results of these analyses leave open the possibility that gliding motility might have arisen independently more than once in the mycoplasma lineage.