Search

Calling cards technology using self-reporting transposons enables the identification of DNA-protein interactions through RNA sequencing. Although immensely powerful, current implementations of calling cards in bulk experiments on populations of cells are technically cumbersome and require many replicates to identify independent insertions into the same genomic locus. Here, we have drastically reduced the cost and labor requirements of calling card experiments in bulk populations of cells by introducing a DNA barcode into the calling card itself. An additional barcode incorporated during reverse transcription enables simultaneous transcriptome measurement in a facile and affordable protocol. We demonstrate that barcoded self-reporting transposons recover in vitro binding sites for four basic helix-loop-helix transcription factors with important roles in cell fate specification: ASCL1, MYOD1, NEUROD2, and NGN1. Further, simultaneous calling cards and transcriptional profiling during transcription factor overexpression identified both binding sites and gene expression changes for two of these factors. Lastly, we demonstrated barcoded calling cards can record binding in vivo in the mouse brain. In sum, RNA-based identification of transcription factor binding sites and gene expression through barcoded self-reporting transposon calling cards and transcriptomes is an efficient and powerful method to infer gene regulatory networks in a population of cells.

Transcription factors (TFs) enact precise regulation of gene expression through site-specific, genome-wide binding. Common methods for TF-occupancy profiling, such as chromatin immunoprecipitation, are limited by requirement of TF-specific antibodies and provide only end-point snapshots of TF binding. Alternatively, TF-tagging techniques, in which a TF is fused to a DNA-modifying enzyme that marks TF-binding events across the genome as they occur, do not require TF-specific antibodies and offer the potential for unique applications, such as recording of TF occupancy over time and cell type specificity through conditional expression of the TF–enzyme fusion. Here, we create a viral toolkit for one such method, calling cards, and demonstrate that these reagents can be delivered to the live mouse brain and used to report TF occupancy. Further, we establish a Cre-dependent calling cards system and, in proof-of-principle experiments, show utility in defining cell type-specific TF profiles and recording and integrating TF-binding events across time. This versatile approach will enable unique studies of TF-mediated gene regulation in live animal models.

Differential epidermal growth factor receptor (EGFR) phosphorylation is thought to couple receptor activation to distinct signaling pathways. However, the molecular mechanisms responsible for biased signaling are unresolved due to a lack of insight into the phosphorylation patterns of full-length EGFR. We extended a single-molecule pull-down technique previously used to study protein-protein interactions to allow for robust measurement of receptor phosphorylation. We found that EGFR is predominantly phosphorylated at multiple sites, yet phosphorylation at specific tyrosines is variable and only a subset of receptors share phosphorylation at the same site, even with saturating ligand concentrations. We found distinct populations of receptors as soon as 1 min after ligand stimulation, indicating early diversification of function. To understand this heterogeneity, we developed a mathematical model. The model predicted that variations in phosphorylation are dependent on the abundances of signaling partners, while phosphorylation levels are dependent on dimer lifetimes. The predictions were confirmed in studies of cell lines with different expression levels of signaling partners, and in experiments comparing low- and high-affinity ligands and oncogenic EGFR mutants. These results reveal how ligand-regulated receptor dimerization dynamics and adaptor protein concentrations play critical roles in EGFR signaling.

A purchased lot of the antibiotic hygromycin B was found to be contaminated with a novel bacterial species, which we designate Pseudomonas hygromyciniae. Characteristics of P. hygromyciniae include its ability to use a variety of compounds as carbon sources, its pathogenicity towards lettuce and Galleria mellonella, and its ability to inhibit the growth of an E. coli strain. P. hygromyciniae is unlikely to be a human pathogen, as it did not survive at 37 °C and was not cytotoxic towards a mammalian cell line. The P. hygromyciniae strain harbors a novel 250 kb megaplasmid which confers resistance to hygromycin B and contains numerous other genes predicted to encode replication and conjugation machinery. These findings indicate that commercially manufactured antibiotics represent another extreme environment that may support the growth of novel bacterial species.IMPORTANCEMicrobial ecologists have surveyed numerous natural and manmade environments in search of new microbial species. In some instances, these microbes are discovered in harsh conditions, such as deep-sea vents, and their discovery leads to better understanding of how microbes adapt to their environment. Here, we have discovered a new species of bacteria from an extreme manmade environment: a lyophilized, commercially available bottle of the antibiotic hygromycin B.

Significance Consistent sex differences in incidence and outcome have been reported in numerous cancers including brain tumors. GBM, the most common and aggressive primary brain tumor, occurs with higher incidence and shorter survival in males compared to females. Brd4 is essential for regulating transcriptome-wide gene expression and specifying cell identity, including that of GBM. We report that sex-biased Brd4 activity drives sex differences in GBM and renders male and female tumor cells differentially sensitive to BET inhibitors. The observed sex differences in BETi treatment strongly indicate that sex differences in disease biology translate into sex differences in therapeutic responses. This has critical implications for clinical use of BET inhibitors further affirming the importance of inclusion of sex as a biological variable., Sex can be an important determinant of cancer phenotype, and exploring sex-biased tumor biology holds promise for identifying novel therapeutic targets and new approaches to cancer treatment. In an established isogenic murine model of glioblastoma (GBM), we discovered correlated transcriptome-wide sex differences in gene expression, H3K27ac marks, large Brd4-bound enhancer usage, and Brd4 localization to Myc and p53 genomic binding sites. These sex-biased gene expression patterns were also evident in human glioblastoma stem cells (GSCs). These observations led us to hypothesize that Brd4-bound enhancers might underlie sex differences in stem cell function and tumorigenicity in GBM. We found that male and female GBM cells exhibited sex-specific responses to pharmacological or genetic inhibition of Brd4. Brd4 knockdown or pharmacologic inhibition decreased male GBM cell clonogenicity and in vivo tumorigenesis while increasing both in female GBM cells. These results were validated in male and female patient-derived GBM cell lines. Furthermore, analysis of the Cancer Therapeutic Response Portal of human GBM samples segregated by sex revealed that male GBM cells are significantly more sensitive to BET (bromodomain and extraterminal) inhibitors than are female cells. Thus, Brd4 activity is revealed to drive sex differences in stem cell and tumorigenic phenotypes, which can be abrogated by sex-specific responses to BET inhibition. This has important implications for the clinical evaluation and use of BET inhibitors.

Most human genetic variation is classified as VUS - variants of uncertain significance. While advances in genome editing have allowed innovation in pooled screening platforms, many screens deal with relatively simple readouts (viability, fluorescence) and cannot identify the complex cellular phenotypes that underlie most human diseases. In this paper, we present a generalizable functional genomics platform that combines high-content imaging, machine learning, and microraft isolation in a new method termed “Raft-Seq”. We highlight the efficacy of our platform by showing its ability to distinguish pathogenic point mutations of the mitochondrial regulator MFN2, even when the cellular phenotype is subtle. We also show that our platform achieves its efficacy using multiple cellular features, which can be configured on-the-fly. Raft-Seq enables a new way to perform pooled screening on sets of mutations in biologically relevant cells, with the ability to physically capture any cell with a perturbed phenotype and expand it clonally, directly from the primary screen.Graphical AbstractHere, we address the need to evaluate the impact of numerous genetic variants. This manuscript depicts the methods of using machine learning on a biologically relevant phenotype to predict specific point mutations, followed by physically capturing those mutated cells.

BackgroundKlebsiella pneumoniae (K. pneumoniae) strains have been divided into two major categories: classical K. pneumonia, which are frequently multidrug-resistant and cause hospital-acquired infections in patients with impaired defenses, and hypervirulent K. pneumoniae, which cause severe disseminated community-acquired infections in immunologically normal hosts. Both types of infections may lead to bacteremia and are associated with significant morbidity and mortality. The relative burden of these two types of K. pneumoniae among bloodstream isolates within the United States is not well understood.MethodsWe examined consecutive K. pneumoniae isolates cultured from the blood of hospitalized patients at Northwestern Memorial Hospital (NMH) in Chicago, Illinois April 2015 – April 2017. Bloodstream isolates underwent whole genome sequencing, and clinically relevant data (sequence type [ST], capsule loci, virulence factors, and antimicrobial resistance genes) were inferred from the genome using the bioinformatic tools Kleborate and Kaptive. Patient demographic, comorbidity, and infection information, as well as phenotypic antimicrobial resistance of the isolate were extracted from the electronic medical record. Candidate hypervirulent isolates were tested in a murine model of pneumonia. Complete genome sequences of these candidate hypervirulent isolates were obtained and plasmid content was evaluated.ResultsK. pneumoniae bloodstream isolates (n=104) from NMH were highly diverse consisting of 75 distinct STs and 51 unique capsule loci. The majority of these isolates (n=58, 55.8%) were susceptible to all tested antibiotics except ampicillin, but 17 (16.3%) were multidrug-resistant. A total of 32 (30.8%) of these isolates were STs of known high-risk clones, including ST258 and ST45. In particular, 18 (17.3%) were resistant to ceftriaxone, 17 harbored extended-spectrum beta-lactamases, and 9 (8.7%) were resistant to meropenem. Four (3.8%) of the 104 isolates were hypervirulent K. pneumoniae, as evidenced by hypermucoviscous phenotypes, high levels of virulence in a murine model of pneumonia, and the presence of large plasmids similar to characterized hypervirulence plasmids. Of particular concern, several of these plasmids contained tra conjugation loci suggesting the potential for transmission. Two of these hypervirulent strains belonged to the well characterized ST23 lineage and one to the emerging ST66 lineage.ConclusionsWhile NMH K. pneumoniae bloodstream infections are caused by highly diverse strains, the most prevalent STs were those of multidrug-resistant high-risk clones. A small number of hypervirulent K. pneumoniae were observed in patients with no recent travel history, suggesting that these isolates are undergoing community spread in the United States.

Nanopore sequencing has been investigated as a rapid and cost-efficient option for HLA typing in recent years. Despite the lower raw read accuracy, encouraging typing accuracy has been reported, and long reads from the platform offer additional benefits of the improved phasing of distant variants. The newly released R10.3 flow cells are expected to provide higher read-level accuracy than previous chemistries. We examined the performance of R10.3 flow cells on the MinION device in HLA typing after enrichment of target genes by multiplexed PCR. We also aimed to mimic a 1-day workflow with 8-24 samples per sequencing run. A diverse collection of 102 unique samples were typed for HLA-A, -B, -C, -DPA1, -DPB1, -DQA1, -DQB1, -DRB1, -DRB3/4/5 loci. The concordance rates at 2-field and 3-field resolutions were 99.5% (1836 alleles) and 99.3% (1710 alleles). We also report important quality metrics from these sequencing runs. Continued research and independent validations are warranted to increase the robustness of nanopore-based HLA typing for broad clinical application.

Long INterspersed Elements-1 (L1s) constitute >17% of the human genome and still actively transpose in it. Characterizing L1 transposition across the genome is critical for understanding genome evolution and somatic mutations. However, to date, L1 insertion and fixation patterns have not been studied comprehensively. To fill this gap, we investigated three genome-wide data sets of L1s that integrated at different evolutionary times: 17,037 de novo L1s (from an L1 insertion cell-line experiment conducted in-house), and 1,212 polymorphic and 1,205 human-specific L1s (from public databases). We characterized 49 genomic features—proxying chromatin accessibility, transcriptional activity, replication, recombination, etc.—in the ±50 kb flanks of these elements. These features were contrasted between the three L1 data sets and L1-free regions using state-of-the-art Functional Data Analysis statistical methods, which treat high-resolution data as mathematical functions. Our results indicate that de novo, polymorphic, and human-specific L1s are surrounded by different genomic features acting at specific locations and scales. This led to an integrative model of L1 transposition, according to which L1s preferentially integrate into open-chromatin regions enriched in non-B DNA motifs, whereas they are fixed in regions largely free of purifying selection—depleted of genes and noncoding most conserved elements. Intriguingly, our results suggest that L1 insertions modify local genomic landscape by extending CpG methylation and increasing mononucleotide microsatellite density. Altogether, our findings substantially facilitate understanding of L1 integration and fixation preferences, pave the way for uncovering their role in aging and cancer, and inform their use as mutagenesis tools in genetic studies.

Enrichment of electrochemically-active microorganisms (EAM) to date has mostly relied on microbial fuel cells fed with wastewater. This study aims to enrich novel EAM by exposing tropical sediment, not frequently reported in the literature, to sustained anodic potentials. Voltamperometric techniques and electrochemical impedance spectroscopy, performed over a wide range of potentials, characterise extracellular electron transfer (EET) over time. Applied potential is found to affect biofilm electrochemical signature. Geobacter metallireducens is heavily enriched on the electrodes, as determined by metagenomic and metatranscriptomic analysis, in the first report of the species in a lactate-fed system. Two novel isolates are grown in pure culture from the enrichment, identified by 16S rRNA gene sequencing as Aeromonas and Enterobacter, respectively. The names proposed are Aeromonas sp. CL-1 and Enterobacter sp. EA-1. Both isolates are capable of EET on carbon felt and screen-printed carbon electrodes without the addition of exogenous redox mediators. Enterobacter sp. EA-1 can also perform mediated electron transfer using the soluble redox mediator 2-hydroxy-1,4-naphthoquinone (HNQ). Both isolates are able to use acetate and lactate as electron donors. This work outlines a comprehensive methodology for characterising novel EAM from unconventional inocula.

The overwhelming success of exome- and genome-wide association studies in discovering thousands of disease-associated genes necessitates developing novel high-throughput functional genomics approaches to elucidate the molecular mechanisms of these genes. Here, we have coupled multiplexed repression of neurodevelopmental disease–associated genes to single-cell transcriptional profiling in differentiating human neurons to rapidly assay the functions of multiple genes in a disease-relevant context, assess potentially convergent mechanisms, and prioritize genes for specific functional assays. For a set of 13 autism spectrum disorder (ASD)–associated genes, we show that this approach generated important mechanistic insights, revealing two functionally convergent modules of ASD genes: one that delays neuron differentiation and one that accelerates it. Five genes that delay neuron differentiation (ADNP, ARID1B, ASH1L, CHD2, and DYRK1A) mechanistically converge, as they all dysregulate genes involved in cell-cycle control and progenitor cell proliferation. Live-cell imaging after individual ASD-gene repression validated this functional module, confirming that these genes reduce neural progenitor cell proliferation and neurite growth. Finally, these functionally convergent ASD gene modules predicted shared clinical phenotypes among individuals with mutations in these genes. Altogether, these results show the utility of a novel and simple approach for the rapid functional elucidation of neurodevelopmental disease-associated genes.

A high-confidence map of the direct, functional targets of each transcription factor (TF) requires convergent evidence from independent sources. Two significant sources of evidence are TF binding locations and the transcriptional responses to direct TF perturbations. Systematic data sets of both types exist for yeast and human, but they rarely converge on a common set of direct, functional targets for a TF. Even the few genes that are both bound and responsive may not be direct functional targets. Our analysis shows that when there are many nonfunctional binding sites and many indirect targets, nonfunctional sites are expected to occur in the cis-regulatory DNA of indirect targets by chance. To address this problem, we introduce dual threshold optimization (DTO), a new method for setting significance thresholds on binding and perturbation-response data, and show that it improves convergence. It also enables comparison of binding data to perturbation-response data that have been processed by network inference algorithms, which further improves convergence. The combination of dual threshold optimization and network inference greatly expands the high-confidence TF network map in both yeast and human. Next, we analyze a comprehensive new data set measuring the transcriptional response shortly after inducing overexpression of a yeast TF. We also present a new yeast binding location data set obtained by transposon calling cards and compare it to recent ChIP-exo data. These new data sets improve convergence and expand the high-confidence network synergistically.

BackgroundA transcription-factor (TF) network map indicates the direct, functional targets of each TF -- the genes it regulates by binding to their cis-regulatory DNA. Data on the genomic binding locations of each TF and the transcriptional responses to perturbations of its activity, such as overexpressing it, could support TF network mapping. Systematic data sets of both types exist for yeast and for human K562 and HEK293 cells.ResultsIn previous data, most TF binding sites appear to be non-functional, so one cannot take the genes in whose promoters a TF binds as its direct, functional (DF) targets. Taking the genes that are both bound by a TF and responsive to a perturbation of it as its DF targets (intersection algorithm) is also not safe, as we show by deriving a new lower bound on the expected false discovery rate of the intersection algorithm. When there are many non-functional binding sites and many indirect targets, non-functional sites are expected to occur in the cis-regulatory DNA of indirect targets by chance. Dual threshold optimization, a new method for setting significance thresholds on binding and response data, improves the intersection algorithm, as does post-processing perturbation-response data with NetProphet 2.0. A comprehensive new data set measuring the transcriptional response shortly after inducing overexpression of a TF also helps, as does transposon calling cards, a new method for identifying TF binding locations.ConclusionsThe combination of dual threshold optimization and NetProphet greatly expands the high-confidence TF network map in both yeast and human. In yeast, measuring the response shortly after inducing TF overexpression and measuring binding locations by using transposon calling cards improve the network synergistically.

Transcription factors (TFs) play a central role in the regulation of gene expression, controlling everything from cell fate decisions to activity dependent gene expression. However, widely-used methods for TF profiling in vivo (e.g. ChIP-seq) yield only an aggregated picture of TF binding across all cell types present within the harvested tissue; thus, it is challenging or impossible to determine how the same TF might bind different portions of the genome in different cell types, or even to identify its binding events at all in rare cell types in a complex tissue such as the brain. Here we present a versatile methodology, FLEX Calling Cards, for the mapping of TF occupancy in specific cell types from heterogenous tissues. In this method, the TF of interest is fused to a hyperactive piggyBac transposase (hypPB), and this bipartite gene is delivered, along with donor transposons, to mouse tissue via a Cre-dependent adeno-associated virus (AAV). The fusion protein is expressed in Cre-expressing cells where it inserts transposon “Calling Cards” near to TF binding sites. These transposons permanently mark TF binding events and can be mapped using high-throughput sequencing. Alternatively, unfused hypPB interacts with and records the binding of the super enhancer (SE)-associated bromodomain protein, Brd4. To demonstrate the FLEX Calling Card method, we first show that donor transposon and transposase constructs can be efficiently delivered to the postnatal day 1 (P1) mouse brain with AAV and that insertion profiles report TF occupancy. Then, using a Cre-dependent hypPB virus, we show utility of this tool in defining cell type-specific TF profiles in multiple cell types of the brain. This approach will enable important cell type-specific studies of TF-mediated gene regulation in the brain and will provide valuable insights into brain development, homeostasis, and disease.

Membrane proteins represent a significant proportion of total bacterial proteins and perform vital cellular functions ranging from exchanging metabolites and genetic material, secretion and sorting, sensing signal molecules, and cell division. Many of these functions are carried out at distinct foci on the bacterial membrane, and this subcellular localization can be coordinated by a number of factors, including lipid microdomains, protein-protein interactions, and membrane curvature. Elucidating the mechanisms behind focal protein localization in bacteria informs not only protein structure-function correlation, but also how to disrupt the protein function to limit virulence. Here we review recent advances describing a functional role for subcellular localization of membrane proteins involved in genetic transfer, secretion and sorting, cell division and growth, and signaling.

Sex can be an important determinant of cancer phenotype and exploring sex-specific tumor biology holds promise for identifying novel therapeutic targets and new approaches to cancer treatment. In an established model of glioblastoma, we discovered transcriptome-wide sexual dimorphism in gene-expression that was concordant with sex differences in H3K27ac marks, large Brd4-bound enhancer usage, and Brd4 co-localization with Myc and p53. The sex-specific enhancer usage drove sex differences in stem cell function and tumorigenicity. Moreover, male and female GBM cells exhibited opposing responses to pharmacological or genetic inhibition of Brd4. Brd4 knockdown or inhibition decreased male GBM cell clonogenicity and in vivo tumorigenesis, while increasing both in female GBM cells. These results were validated in male and female patient-derived GBM cell lines. Thus, for the first time, Brd4 activity is revealed to drive a sexually dimorphic stem cell and tumorigenic phenotype, resulting in diametrically opposite responses to BET inhibition in male and female glioblastoma cells. This has critical implications for the clinical evaluation and use of BET inhibitors. Citation Format: Najla Kfoury-Beaumont, Zongtai Qi, Michael Wilkinson, Lauren Broestl, Kristopher Berrett, Arnav Moudgil, Sumithra Sankararaman, Xuhua Chen, Jay Gertz, Robi Mitra, Joshua B. Rubin. Brd4-bound enhancers drive critical sex differences in glioblastoma [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 3658.

Cellular heterogeneity confounds in situ assays of transcription factor (TF) binding. Single-cell RNA sequencing (scRNA-seq) deconvolves cell types from gene expression, but no technology links cell identity to TF binding sites (TFBS) in those cell types. We present self-reporting transposons (SRTs) and use them in single-cell calling cards (scCC), a novel assay for simultaneously measuring gene expression and mapping TFBS in single cells. The genomic locations of SRTs are recovered from mRNA, and SRTs deposited by exogenous, TF-transposase fusions can be used to map TFBS. We then present scCC, which map SRTs from scRNA-seq libraries, simultaneously identifying cell types and TFBS in those same cells. We benchmark multiple TFs with this technique. Next, we use scCC to discover BRD4-mediated cell-state transitions in K562 cells. Finally, we map BRD4 binding sites in the mouse cortex at single-cell resolution, establishing a new method for studying TF biology in situ.

We report a tool, Calling Cards Reporter Arrays (CCRA), that measures transcription factor (TF) binding and the consequences on gene expression for hundreds of synthetic promoters in yeast. Using Cbf1p and MAX, we demonstrate that the CCRA method is able to detect small changes in binding free energy with a sensitivity comparable to in vitro methods, enabling the measurement of energy landscapes in vivo. We then demonstrate the quantitative analysis of cooperative interactions by measuring Cbf1p binding at synthetic promoters with multiple sites. We find that the cooperativity between Cbf1p dimers varies sinusoidally with a period of 10.65 bp and energetic cost of 1.37 KBT for sites that are positioned ‘out of phase’. Finally, we characterize the binding and expression of a group of TFs, Tye7p, Gcr1p and Gcr2p, that act together as a ‘TF collective’, an important but poorly characterized model of TF cooperativity. We demonstrate that Tye7p often binds promoters without its recognition site because it is recruited by other collective members, whereas these other members require their recognition sites, suggesting a hierarchy where these factors recruit Tye7p but not vice versa. Our experiments establish CCRA as a useful tool for quantitative investigations into TF binding and function.

Background BAP1 is a histone deubiquitinase that acts as a tumor and metastasis suppressor associated with disease progression in human cancer. We have used the “Calling Card System” of transposase-directed transposon insertion mapping to identify the genomic targets of BAP1 in uveal melanoma (UM). This system was developed to identify the genomic loci visited by transcription factors that bind directly to DNA; our study is the first use of the system with a chromatin-remodeling factor that binds to histones but does not interact directly with DNA. Methods The transposase piggyBac (PBase) was fused to BAP1 and expressed in OCM-1A UM cells. The insertion of transposons near BAP1 binding sites in UM cells were identified by genomic sequencing. We also examined RNA expression in the same OCM-1A UM cells after BAP1 depletion to identify BAP1 binding sites associated with BAP1-responsive genes. Sets of significant genes were analyzed for common pathways, transcription factor binding sites, and ability to identify molecular tumor classes. Results We found a strong correlation between multiple calling-card transposon insertions targeted by BAP1-PBase and BAP1-responsive expression of adjacent genes. BAP1-bound genomic loci showed narrow distributions of insertions and were near transcription start sites, consistent with recruitment of BAP1 to these sites by specific DNA-binding proteins. Sequence consensus analysis of BAP1-bound sites showed enrichment of motifs specific for YY1, NRF1 and Ets transcription factors, which have been shown to interact with BAP1 in other cell types. Further, a subset of the BAP1 genomic target genes was able to discriminate aggressive tumors in published gene expression data from primary UM tumors. Conclusions The calling card methodology works equally well for chromatin regulatory factors that do not interact directly with DNA as for transcription factors. This technique has generated a new and expanded list of BAP1 targets in UM that provides important insight into metastasis pathways and identifies novel potential therapeutic targets. Electronic supplementary material The online version of this article (10.1186/s12920-018-0424-0) contains supplementary material, which is available to authorized users.

SUMMARYMolecular and behavioral responses to opioids are thought to be primarily mediated by neurons, although there is accumulating evidence that other cell types also play a role in drug addiction. To investigate cell-type-specific opioid responses, we performed single-cell RNA sequencing of the nucleus accumbens of mice following acute morphine treatment. Differential expression analysis uncovered robust morphine-dependent changes in gene expression in oligodendrocytes. We examined the expression of selected genes, including Cdkn1a and Sgk1, by FISH, confirming their induction by morphine in oligodendrocytes. Further analysis using RNAseq of FACS-purified oligodendrocytes revealed a large cohort of morphine-regulated genes. Importantly, the affected genes are enriched for roles in cellular pathways intimately linked to oligodendrocyte maturation and myelination, including the unfolded protein response. Altogether, our data shed light on a novel, morphine-dependent transcriptional response by oligodendrocytes that may contribute to the myelination defects observed in human opioid addicts.

Cre recombinase catalyzes the cleavage and religation of DNA at loxP sites. The enzyme is a homotetramer in its functional state, and the symmetry of the protein complex enforces a pseudo-palindromic symmetry upon the loxP sequence. The Cre-lox system is a powerful tool for many researchers. However, broader application of the system is limited by the fixed sequence preferences of Cre, which are determined by both the direct DNA contacts and the homotetrameric arrangement of the Cre monomers. As a first step toward achieving recombination at arbitrary asymmetric target sites, we have broken the symmetry of the Cre tetramer assembly. Using a combination of computational and rational protein design, we have engineered an alternative interface between Cre monomers that is functional yet incompatible with the wild-type interface. Wild-type and engineered interface halves can be mixed to create two distinct Cre mutants, neither of which are functional in isolation, but which can form an active heterotetramer when combined. When these distinct mutants possess different DNA specificities, control over complex assembly directly discourages recombination at unwanted half-site combinations, enhancing the specificity of asymmetric site recombination. The engineered Cre mutants exhibit this assembly pattern in a variety of contexts, including mammalian cells.

The Epidermal Growth Factor Receptor (EGFR/ErbB1/HER1) plays an important role in both physiological and cancer-related processes. To study the factors that influence receptor phosphorylation, we have coupled Single Molecule Pull-down (SiMPull) measurements with rule-based modeling of EGFR signaling. Using SiMPull, we quantified the phosphorylation state of thousands of individual receptors. These measurements enabled the first direct detection of multisite phosphorylation on full-length EGFR and revealed that the extent of phosphorylation varies by tyrosine site and is dependent on the relative abundance of signaling partners that limit access by tyrosine phosphatases. We also evaluated the impact of oncogenic mutations and ligands with varying affinity on phosphorylation kinetics. Simulations highlight the importance of dimer lifetimes on EGFR phosphorylation and signaling output.

The high-affinity receptor for IgE expressed on the surface of mast cells and basophils interacts with antigens, via bound IgE antibody, and triggers secretion of inflammatory mediators that contribute to allergic reactions. To understand how past inputs (memory) influence future inflammatory responses in mast cells, a microfluidic device was used to precisely control exposure of cells to alternating stimulatory and non-stimulatory inputs. We determined that the response to subsequent stimulation depends on the interval of signaling quiescence. For shorter intervals of signaling quiescence, the second response is blunted relative to the first response, whereas longer intervals of quiescence induce an enhanced second response. Through an iterative process of computational modeling and experimental tests, we found that these memory-like phenomena arise from a confluence of rapid, short-lived positive signals driven by the protein tyrosine kinase Syk; slow, long-lived negative signals driven by the lipid phosphatase Ship1; and slower degradation of Ship1 co-factors. This work advances our understanding of mast cell signaling and represents a generalizable approach for investigating the dynamics of signaling systems.

Sex can be an important determinant of cancer phenotype, and exploring sex-biased tumor biology holds promise for identifying novel therapeutic targets and new approaches to cancer treatment. In an established isogenic murine model of glioblastoma, we discovered correlated transcriptome-wide sex differences in gene expression, H3K27ac marks, large Brd4-bound enhancer usage, and Brd4 localization to Myc and p53 genomic binding sites. These sex-biased gene expression patterns were also evident in human glioblastoma stem cells (GSCs). These observations led us to hypothesize that Brd4-bound enhancers might underlie sex differences in stem cell function and tumorigenicity in GBM. We found that male and female GBM cells exhibited opposing responses to pharmacological or genetic inhibition of Brd4. Brd4 knockdown or pharmacologic inhibition decreased male GBM cell clonogenicity and in vivo tumorigenesis, while increasing both in female GBM cells. These results were validated in male and female patient-derived GBM cell lines. Furthermore, analysis of the Cancer Therapeutic Response Portal of human GBM samples segregated by sex revealed that male GBM cells are significantly more sensitive to BET inhibitors than are female cells. Thus, for the first time, Brd4 activity is revealed to drive a sex differences in stem cell and tumorigenic phenotype, resulting in diametrically opposite responses to BET inhibition in male and female GBM cells. This has important implications for the clinical evaluation and use of BET inhibitors.SignificanceConsistent sex differences in incidence and outcome have been reported in numerous cancers including brain tumors. GBM, the most common and aggressive primary brain tumor, occurs with higher incidence and shorter survival in males compared to females. Brd4 is essential for regulating transcriptome-wide gene expression and specifying cell identity, including that of GBM. We report that sex-biased Brd4 activity drive sex differences in GBM and render male and female tumor cells differentially sensitive to BET inhibitors. The observed sex differences in BETi treatment strongly indicate that sex differences in disease biology translate into sex differences in therapeutic responses. This has critical implications for clinical use of BET inhibitors further affirming the importance of inclusion of sex as a biological variable.

Oxford Nanopore Technologies’ MinION has expanded the current DNA sequencing toolkit by delivering long read lengths and extreme portability. The MinION has the potential to enable expedited point-of-care human leukocyte antigen (HLA) typing, an assay routinely used to assess the immunological compatibility between organ donors and recipients, but the platform’s high error rate makes it challenging to type alleles with clinical-grade accuracy. Here, we developed and validated Athlon, an algorithm that iteratively scores nanopore reads mapped to a hierarchical database of HLA alleles to arrive at a consensus diploid genotype; Athlon achieved a 100% accuracy in class I HLA typing at high resolution.

Objective To define the genetic landscape of amyotrophic lateral sclerosis (ALS) and assess the contribution of possible oligogenic inheritance, we aimed to comprehensively sequence 17 known ALS genes in 391 ALS patients from the United States. Methods Targeted pooled-sample sequencing was used to identify variants in 17 ALS genes. Fragment size analysis was used to define ATXN2 and C9ORF72 expansion sizes. Genotype–phenotype correlations were made with individual variants and total burden of variants. Rare variant associations for risk of ALS were investigated at both the single variant and gene level. Results A total of 64.3% of familial and 27.8% of sporadic subjects carried potentially pathogenic novel or rare coding variants identified by sequencing or an expanded repeat in C9ORF72 or ATXN2; 3.8% of subjects had variants in >1 ALS gene, and these individuals had disease onset 10 years earlier (p = 0.0046) than subjects with variants in a single gene. The number of potentially pathogenic coding variants did not influence disease duration or site of onset. Interpretation Rare and potentially pathogenic variants in known ALS genes are present in >25% of apparently sporadic and 64% of familial patients, significantly higher than previous reports using less comprehensive sequencing approaches. A significant number of subjects carried variants in >1 gene, which influenced the age of symptom onset and supports oligogenic inheritance as relevant to disease pathogenesis. ANN NEUROL 2015;77:100–113

Conservation and management of Sundarban mangrove forest is difficult chiefly due to inaccessibility and hostile condition. Remote sensing serves as an important tool to provide up-to date baseline information which is the primary requirement for the conservation planning of mangroves. In this study, supervised classification by maximum likelihood classifier (MLC) has been used to classify LANDSAT TM and LANDSAT ETM satellite data. This algorithm is used for computing likelihood of unknown measurement vector belonging to unknown classes based on Bayesian equation. Image spectra for various mangrove species were also generated from hyperspectral image. During field visits, GPS locations of five dominant mangrove species with appreciable distribution were taken and image spectra were generated for the same points from hyperion image. The result of this classification shows that, in 1999 total mangrove forest accounted for 55.01 % of the study area which has been reduced to 50.63 % in the year 2010. Avicennia sp. is found as most dominating species followed by Excoecaria sp. and Phoenix sp. but the aerial distribution of Avicennia sp., Bruguiera sp. and Ceriops sp. has reduced. In this classification technique the overall accuracy and Kappa value for 1999 and 2010 are 80 % and 0.77, 85.71 % and 0.81 respectively.

Currently, oncology testing includes molecular studies and cytogenetic analysis to detect genetic aberrations of clinical significance. Next-generation sequencing (NGS) allows rapid analysis of multiple genes for clinically actionable somatic variants. The WUCaMP assay uses targeted capture for NGS analysis of 25 cancer-associated genes to detect mutations at actionable loci. We present clinical validation of the assay and a detailed framework for design and validation of similar clinical assays. Deep sequencing of 78 tumor specimens (≥ 1000× average unique coverage across the capture region) achieved high sensitivity for detecting somatic variants at low allele fraction (AF). Validation revealed sensitivities and specificities of 100% for detection of single-nucleotide variants (SNVs) within coding regions, compared with SNP array sequence data (95% CI = 83.4-100.0 for sensitivity and 94.2-100.0 for specificity) or whole-genome sequencing (95% CI = 89.1-100.0 for sensitivity and 99.9-100.0 for specificity) of HapMap samples. Sensitivity for detecting variants at an observed 10% AF was 100% (95% CI = 93.2-100.0) in HapMap mixes. Analysis of 15 masked specimens harboring clinically reported variants yielded concordant calls for 13/13 variants at AF of ≥ 15%. The WUCaMP assay is a robust and sensitive method to detect somatic variants of clinical significance in molecular oncology laboratories, with reduced time and cost of genetic analysis allowing for strategic patient management.

Defensins found in mammals belong to mainly two subfamilies α- and β-defensins. Mammalian defensins are small molecules (18–45 residues) that are cysteine, arginine rich compounds. Antimicrobial activities of these peptides were shown against a wide variety of microbes including bacteria, fungi, viruses and protozoan parasites. To investigate the structure and activity relationship, amino acid substitutions that alter charge were introduced into synthetic defensin peptides by adding 2–2 Arg (RR) and Asp (DD) at both the terminal and tested their effects on HIV-1, E. coli, S. aureus, and P. aeruginosa. In the present study, we have chemically synthesized native defensin peptides and their variants with Arg (RR) and Asp (DD) amino acid residues at N- and C-termini. Later, we assayed their anti-HIV, anti-microbial activities, stability, cytotoxicity and hemolytic properties. We reported that anti-HIV and antimicrobial activities of native defensins is increased significantly by adding Arg (RR) residues at both the termini while the substitution of Arg (RR) with Asp (DD), eliminate anti-HIV and antimicrobial activity against all bacterial species tested. While other physical features i.e. stability, cell toxicity and hemolytic property were not affected by any of the changes in the sequence. The results suggest that the terminal residues in defensins are crucial functional elements that determine their microbicidal potency. The enhanced microbicidal activity observed for defensin peptides with Arg (RR) residues could be due to optimization of amphiphilicity of the structure, which could facilitate specific interactions with the microbial membranes.

A significant challenge of functional genomics is to develop methods for genome-scale acquisition and analysis of cell biological data. Here, we present an integrated method that combines genome-wide genetic perturbation of Saccharomyces cerevisiae with high-content screening to facilitate the genetic description of sub-cellular structures and compartment morphology. As proof of principle, we used a Rad52-GFP marker to examine DNA damage foci in ∼20 million single cells from ∼5,000 different mutant backgrounds in the context of selected genetic or chemical perturbations. Phenotypes were classified using a machine learning-based automated image analysis pipeline. 345 mutants were identified that had elevated numbers of DNA damage foci, almost half of which were identified only in sensitized backgrounds. Subsequent analysis of Vid22, a protein implicated in the DNA damage response, revealed that it acts together with the Sgs1 helicase at sites of DNA damage and preferentially binds G-quadruplex regions of the genome. This approach is extensible to numerous other cell biological markers and experimental systems.

Identifying the genomic targets of transcription factors is an important step in understanding the regulatory networks of gene transcription in yeast. We have developed a method that utilizes what we refer to as transposon “calling cards,” in which a transcription factor directs the Ty5 retrotransposase to insert transposons into the genome adjacent to where the transcription factor binds. This method is designed to be multiplexed with many barcoded transcription factors and has the potential to decrease the labor required for the study of large numbers of transcription factors.

Calling card analysis is a high-throughput method for identifying the genomic binding sites of multiple transcription factors in a single experiment in budding yeast. By tagging a DNA-binding protein with a targeting domain that directs the insertion of the Ty5 retrotransposon, the genomic binding sites for that transcription factor are marked. The transposition locations are then identified en masse by Illumina sequencing. The calling card protocol allows for simultaneous analysis of multiple transcription factors. By cloning barcodes into the Ty5 transposon, it is possible to pair a unique barcode with every transcription factor in the experiment. The method presented here uses expression of transcription factors from their native loci; however, it can also be altered to measure binding sites of transcription factors overexpressed from a plasmid.

A revelation of the genomic age has been the contributions of the mobile DNA segments called transposable elements to chromosome structure, function, and evolution in virtually all organisms. Substantial fractions of vertebrate genomes derive from transposable elements, being dominated by retroelements that move via RNA intermediates. Although many of these elements have been inactivated by mutation, several active retroelements remain. Vertebrate genomes also contain substantial quantities and a high diversity of cut-and-paste DNA transposons, but no active representative of this class has been identified in mammals. Here we show that a cut-and-paste element called piggyBat, which has recently invaded the genome of the little brown bat ( Myotis lucifugus ) and is a member of the piggyBac superfamily, is active in its native form in transposition assays in bat and human cultured cells, as well as in the yeast Saccharomyces cerevisiae . Our study suggests that some DNA transposons are still actively shaping some mammalian genomes and reveals an unprecedented opportunity to study the mechanism, regulation, and genomic impact of cut-and-paste transposition in a natural mammalian host.

Infectious Phenotype The pathogenic yeast Candida albicans needs to adopt a filamentous form to invade tissues. The distantly related yeast species Saccharomyces cerevisiae also takes on a filamentous form for nutrient foraging. Comparing genome-wide deletion libraries between the two species, Ryan et al. (p. 1353 ) identified genes involved in three different filamentous yeast phenotypes and found unique genes for each of these phenotypes. However, in addition, core genes, including a previously unknown conserved regulator, appear to have homologous roles in regulating filamentous growth in these distantly related yeast species.

The ability to chronicle transcription-factor binding events throughout the development of an organism would facilitate mapping of transcriptional networks that control cell-fate decisions. We describe a method for permanently recording protein–DNA interactions in mammalian cells. We endow transcription factors with the ability to deposit a transposon into the genome near to where they bind. The transposon becomes a “calling card” that the transcription factor leaves behind to record its visit to the genome. The locations of the calling cards can be determined by massively parallel DNA sequencing. We show that the transcription factor SP1 fused to the piggyBac transposase directs insertion of the piggyBac transposon near SP1 binding sites. The locations of transposon insertions are highly reproducible and agree with sites of SP1-binding determined by ChIP-seq. Genes bound by SP1 are more likely to be expressed in the HCT116 cell line we used, and SP1-bound CpG islands show a strong preference to be unmethylated. This method has the potential to trace transcription-factor binding throughout cellular and organismal development in a way that has heretofore not been possible.

Genome-wide association studies have identified common variation in the CHRNA5-CHRNA3-CHRNB4 and CHRNA6-CHRNB3 gene clusters that contribute to nicotine dependence. However, the role of rare variation in risk for nicotine dependence in these nicotinic receptor genes has not been studied. We undertook pooled sequencing of the coding regions and flanking sequence of the CHRNA5, CHRNA3, CHRNB4, CHRNA6 and CHRNB3 genes in African American and European American nicotine-dependent smokers and smokers without symptoms of dependence. Carrier status of individuals harboring rare missense variants at conserved sites in each of these genes was then compared in cases and controls to test for an association with nicotine dependence. Missense variants at conserved residues in CHRNB4 are associated with lower risk for nicotine dependence in African Americans and European Americans (AA P = 0.0025, odds-ratio (OR) = 0.31, 95% confidence-interval (CI) = 0.31-0.72; EA P = 0.023, OR = 0.69, 95% CI = 0.50-0.95). Furthermore, these individuals were found to smoke fewer cigarettes per day than non-carriers (AA P = 6.6 × 10(-5), EA P = 0.021). Given the possibility of stochastic differences in rare allele frequencies between groups replication of this association is necessary to confirm these findings. The functional effects of the two CHRNB4 variants contributing most to this association (T375I and T91I) and a missense variant in CHRNA3 (R37H) in strong linkage disequilibrium with T91I were examined in vitro. The minor allele of each polymorphism increased cellular response to nicotine (T375I P = 0.01, T91I P = 0.02, R37H P = 0.003), but the largest effect on in vitro receptor activity was seen in the presence of both CHRNB4 T91I and CHRNA3 R37H (P = 2 × 10(-6)).

Summary: Molecular and behavioral responses to opioids are thought to be primarily mediated by neurons, although there is accumulating evidence that other cell types play a prominent role in drug addiction. To investigate cell-type-specific opioid responses, we performed single-cell RNA sequencing (scRNA-seq) of the nucleus accumbens of mice following acute morphine treatment. Differential expression analysis uncovered unique morphine-dependent transcriptional responses by oligodendrocytes and astrocytes. We examined the expression of selected genes, including Cdkn1a and Sgk1, by FISH, confirming their induction by morphine in oligodendrocytes. Further analysis using RNA-seq of FACS-purified oligodendrocytes revealed a large cohort of morphine-regulated genes. The affected genes are enriched for roles in cellular pathways intimately linked to oligodendrocyte maturation and myelination, including the unfolded protein response. Altogether, our data illuminate the morphine-dependent transcriptional response by oligodendrocytes and offer mechanistic insights into myelination defects associated with opioid abuse. : Avey et al. use single-cell RNA-seq to uncover cell-type-specific responses to morphine in the brains of mice, revealing a unique and robust response by oligodendrocytes, which they further validate by multiple approaches. These analyses offer insights into the molecular mechanisms underlying oligodendrocyte dysfunction in the context of opioid abuse. Keywords: opioid, morphine, addiction, single-cell, RNA-seq, oligodendrocyte, UPR, glucocorticoid, nucleus accumbens, myelin

The expression of most genes is regulated by multiple transcription factors. The interactions between transcription factors produce complex patterns of gene expression that are not always obvious from the arrangement of cis-regulatory elements in a promoter. One critical element of promoters is the TATA box, the docking site for the RNA polymerase holoenzyme. Using a synthetic promoter system coupled to a thermodynamic model of combinatorial regulation, we analyze the effects of different strength TATA boxes on various aspects of combinatorial cis-regulation. The thermodynamic model explains 75% of the variance in gene expression in synthetic promoter libraries with different strength TATA boxes, suggesting that many of the salient aspects of cis-regulation are captured by the model. Our results demonstrate that the effect of changing the TATA box on gene expression is the same for all synthetic promoters regardless of the arrangement of cis-regulatory sites we studied. Our analysis also showed that in our synthetic system the strength of the RNA polymerase–TATA interaction does not alter the combinatorial interactions between transcription factors, or between transcription factors and RNA polymerase. Finally, we show that although stronger TATA boxes increase expression in a predictable fashion, stronger TATA boxes have very little effect on noise in our synthetic promoters, regardless of the arrangement of cis-regulatory sites. Our results support a modular model of promoter function, where cis-regulatory elements can be mixed and matched (programmed) with outcomes on expression that are predictable based on the rules of simple protein–protein and protein–DNA interactions.

SummaryThe human gut microbiota is a metabolic organ whose cellular composition is determined by a dynamic process of selection and competition. To identify microbial genes required for establishment of human symbionts in the gut, we developed an approach (insertion sequencing, or INSeq) based on a mutagenic transposon that allows capture of adjacent chromosomal DNA to define its genomic location. We used massively parallel sequencing to monitor the relative abundance of tens of thousands of transposon mutants of a saccharolytic human gut bacterium, Bacteroides thetaiotaomicron, as they established themselves in wild-type and immunodeficient gnotobiotic mice, in the presence or absence of other human gut commensals. In vivo selection transforms this population, revealing functions necessary for survival in the gut: we show how this selection is influenced by community composition and competition for nutrients (vitamin B12). INSeq provides a broadly applicable platform to explore microbial adaptation to the gut and other ecosystems.

Medical resequencing of candidate genes in individual patient samples is becoming increasingly important in the clinic and in clinical research. Medical resequencing requires the amplification and sequencing of many candidate genes in many patient samples. Here we introduce Nested Patch PCR, a novel method for highly multiplexed PCR that is very specific, can sensitively detect SNPs and mutations, and is easy to implement. This is the first method that couples multiplex PCR with sample-specific DNA barcodes and next-generation sequencing to enable highly multiplex mutation discovery in candidate genes for multiple samples in parallel. In our pilot study, we amplified exons from colon cancer and matched normal human genomic DNA. From each sample, we successfully amplified 96% (90 of 94) targeted exons from across the genome, totaling 21.6 kbp of sequence. Ninety percent of all sequencing reads were from targeted exons, demonstrating that Nested Patch PCR is highly specific. We found that the abundance of reads per exon was reproducible across samples. We reliably detected germline SNPs and discovered a colon tumor specific nonsense mutation in APC, a gene causally implicated in colorectal cancer. With Nested Patch PCR, candidate gene mutation discovery across multiple individual patient samples can now utilize the power of second-generation sequencing.

We present a protocol for a novel method for identifying the targets of DNA-binding proteins in the genome of the yeast Saccharomyces cerevisiae. This is accomplished by engineering a DNA-binding protein so that it leaves behind in the genome a permanent mark -- a 'calling card' -- that provides a record of that protein's visit to that region of the genome. The calling card is the yeast Ty5 retrotransposon, whose integrase interacts with the Sir4 protein. If Sir4 is fused to a DNA-binding protein, it recruits the Ty5 integrase, which directs insertion of a Ty5 calling card into the genome. The calling card along with the flanking genomic DNA is harvested by inverse PCR and its genomic location is determined by hybridization of the product to a DNA microarray. This method provides a straightforward alternative to the 'ChIP-chip' method for determining the targets of DNA-binding proteins. This protocol takes approximately 2 weeks to complete.

Identifying genomic targets of transcription factors is fundamental for understanding transcriptional regulatory networks. Current technology enables identification of all targets of a single transcription factor, but there is no realistic way to achieve the converse: identification of all proteins that bind to a promoter of interest. We have developed a method that promises to fill this void. It employs the yeast retrotransposon Ty5, whose integrase interacts with the Sir4 protein. A DNA-binding protein fused to Sir4 directs insertion of Ty5 into the genome near where it binds; the Ty5 becomes a “calling card” the DNA-binding protein leaves behind in the genome. We constructed customized calling cards for seven transcription factors of yeast by including in each Ty5 a unique DNA sequence that serves as a “molecular bar code.” Ty5 transposition was induced in a population of yeast cells, each expressing a different transcription factor–Sir4 fusion and its matched, bar-coded Ty5, and the calling cards deposited into selected regions of the genome were identified, revealing the transcription factors that visited that region of the genome. In each region we analyzed, we found calling cards for only the proteins known to bind there: In the GAL1–10 promoter we found only calling cards for Gal4; in the HIS4 promoter we found only Gcn4 calling cards; in the PHO5 promoter we found only Pho4 and Pho2 calling cards. We discuss how Ty5 calling cards might be implemented for mapping all targets of all transcription factors in a single experiment.

To better understand the complex role that alternative splicing plays in intracellular signaling, it is important to catalog the numerous splice variants involved in signal transduction. Therefore, we developed PASE (Prediction of Alternative Signaling Exons), a computational tool to identify novel alternative cassette exons that code for kinase phosphorylation or signaling protein-binding sites. We first applied PASE to the Caenorhabditis elegans genome. In this organism, our algorithm had an overall specificity of > or =76.4%, including 33 novel cassette exons that we experimentally verified. We then used PASE to analyze the human genome and made 804 predictions, of which 308 were found as alternative exons in the transcript database. We experimentally tested 384 of the remaining unobserved predictions and discovered 26 novel human exons for a total specificity of > or =41.5% in human. By using a test set of known alternatively spliced signaling exons, we determined that the sensitivity of PASE is approximately 70%. GO term analysis revealed that our exon predictions were found in the introns of known signal transduction genes more often than expected by chance, indicating PASE enriches for splice variants that function in signaling pathways. Overall, PASE was able to uncover 59 novel alternative cassette exons in C. elegans and humans through a genome-wide ab initio prediction method that enriches for exons involved in signaling.

Focal segmental glomerulosclerosis (FSGS) is a syndrome that involves kidney podocyte dysfunction and causes chronic kidney disease. Multiple factors including chemical toxicity, inflammation, and infection underlie FSGS; however, highly penetrant disease genes have been identified in a small fraction of patients with a family history of FSGS. Variants of apolipoprotein L1 (APOL1) have been linked to FSGS in African Americans with HIV or hypertension, supporting the proposal that genetic factors enhance FSGS susceptibility. Here, we used sequencing to investigate whether genetics plays a role in the majority of FSGS cases that are identified as primary or sporadic FSGS and have no known cause. Given the limited number of biopsy-proven cases with ethnically matched controls, we devised an analytic strategy to identify and rank potential candidate genes and used an animal model for validation. Nine candidate FSGS susceptibility genes were identified in our patient cohort, and three were validated using a high-throughput mouse method that we developed. Specifically, we introduced a podocyte-specific, doxycycline-inducible transactivator into a murine embryonic stem cell line with an FSGS-susceptible genetic background that allows shRNA-mediated targeting of candidate genes in the adult kidney. Our analysis supports a broader role for genetic susceptibility of both sporadic and familial cases of FSGS and provides a tool to rapidly evaluate candidate FSGS-associated genes.

We report a method for multilocus long-range haplotyping on human chromosome molecules in vitro based on the DNA polymerase colony (polony) technology. By immobilizing thousands of intact chromosome molecules within a polyacrylamide gel on a microscope slide and performing multiple amplifications from single molecules, we determined long-range haplotypes spanning a 153-Mb region of human chromosome 7 and found evidence of rare mitotic recombination events in human lymphocytes. Furthermore, the parallel nature of DNA polony technology allows efficient haplotyping on pooled DNAs from a population on one slide, with a throughput three orders of magnitudes higher than current molecular haplotyping methods. Linkage disequilibrium statistics established by our pooled DNA haplotyping method are more accurate than statistically inferred haplotypes. This haplotyping method is well suited for candidate gene-based association studies as well as for investigating the pattern of recombination in mammalian cells.

We describe a DNA sequencing technology in which a commonly available, inexpensive epifluorescence microscope is converted to rapid nonelectrophoretic DNA sequencing automation. We apply this technology to resequence an evolved strain of Escherichia coli at less than one error per million consensus bases. A cell-free, mate-paired library provided single DNA molecules that were amplified in parallel to 1-micrometer beads by emulsion polymerase chain reaction. Millions of beads were immobilized in a polyacrylamide gel and subjected to automated cycles of sequencing by ligation and four-color imaging. Cost per base was roughly one-ninth as much as that of conventional sequencing. Our protocols were implemented with off-the-shelf instrumentation and reagents.

One of the primary goals of functional genomics is to provide a quantitative understanding of gene function. However, the success of this enterprise is dependent on the accuracy and precision of the functional genomic data. A novel approach, digital analysis of gene expression (DAGE) described herein, is an accurate and precise technology for measuring digital gene expression on a relative or absolute scale by simply counting the number of transcripts of a gene being expressed at a given time. The result is a greatly improved technology sensitive enough for identifying and quantifying small (but biologically important and statistically relevant) changes in gene expression. Fourteen genes involved in galactose metabolism in Saccharomyces cerevisiae were analyzed for their expression levels in glucose and galactose minimal media. The quantitative expression results were characterized in terms of distributional and accuracy attributes; they were also in general agreement (in terms of direction of change) with corresponding results obtained using microarray technology. DAGE is likely to have profound implications in the field of functional genomics because the gene expression measurements are digital in nature and therefore more accurate than any other technologies.

Large-scale mutagenesis of target DNA sequences allows researchers to comprehensively assess the effects of single-nucleotide changes. Here we demonstrate the construction of a systematic allelic series (SAS) using massively parallel single-nucleotide mutagenesis with reversibly terminated deoxyinosine triphosphates (rtITP). We created a mutational library containing every possible single-nucleotide mutation surrounding the active site of the TEM-1 β-lactamase gene. When combined with high-throughput functional assays, SAS mutational libraries can expedite the functional assessment of genetic variation.