Comprehensive Genome-wide Perturbations via CRISPR Adaptation Reveal Complex Genetics of Antibiotic Sensitivity.

Genome-wide CRISPR screens enable systematic interrogation of gene function. However, guide RNA libraries are costly to synthesize, and their limited diversity compromises the sensitivity of CRISPR screens. Using the Streptococcus pyogenes CRISPR-Cas adaptation machinery, we developed CRISPR adaptation-mediated library manufacturing (CALM), which turns bacterial cells into "factories" for generating hundreds of thousands of crRNAs covering 95% of all targetable genomic sites. With an average gene targeted by more than 100 distinct crRNAs, these highly comprehensive CRISPRi libraries produced varying degrees of transcriptional repression critical for uncovering novel antibiotic resistance determinants. Furthermore, by iterating CRISPR adaptation, we rapidly generated dual-crRNA libraries representing more than 100,000 dual-gene perturbations. The polarized nature of spacer adaptation revealed the historical contingency in the stepwise acquisition of genetic perturbations leading to increasing antibiotic resistance. CALM circumvents the expense, labor, and time required for synthesis and cloning of gRNAs, allowing generation of CRISPRi libraries in wild-type bacteria refractory to routine genetic manipulation.

Genome-wide analysis of the host intracellular network that regulates survival of Mycobacterium tuberculosis.

We performed a genome-wide siRNA screen to identify host factors that regulated pathogen load in human macrophages infected with a virulent strain of Mycobacterium tuberculosis. Iterative rounds of confirmation, followed by validation, identified 275 such molecules that were all found to functionally associate with each other through a dense network of interactions. This network then yielded to a molecular description of the host cell functional modules that were both engaged and perturbed by the pathogen. Importantly, a subscreen against a panel of field isolates revealed that the molecular composition of the host interface varied with both genotype and the phenotypic properties of the pathogen. An analysis of these differences, however, permitted identification of those host factors that were invariantly involved, regardless of the diversification in adaptive mechanisms employed by the pathogen. Interestingly, these factors were found to predominantly function through the regulation of autophagy.

High-Throughput Approaches to Pinpoint Function within the Noncoding Genome.

The clustered regularly interspaced short palindromic repeats (CRISPR)-Cas nuclease system is a powerful tool for genome editing, and its simple programmability has enabled high-throughput genetic and epigenetic studies. These high-throughput approaches offer investigators a toolkit for functional interrogation of not only protein-coding genes but also noncoding DNA. Historically, noncoding DNA has lacked the detailed characterization that has been applied to protein-coding genes in large part because there has not been a robust set of methodologies for perturbing these regions. Although the majority of high-throughput CRISPR screens have focused on the coding genome to date, an increasing number of CRISPR screens targeting noncoding genomic regions continue to emerge. Here, we review high-throughput CRISPR-based approaches to uncover and understand functional elements within the noncoding genome and discuss practical aspects of noncoding library design and screen analysis.

Assessment of DNA quality for whole genome library preparation.

In recent years, more sophisticated DNA technologies for genotyping have enabled considerable progress in various fields such as clinical genetics, archaeogenetics and forensic genetics. DNA samples previously rejected as too challenging to analyze due to low amounts of degraded DNA can now provide useful information. To increase the chances of success with the new methodologies, it is crucial to know the fragment size of the template DNA molecules, and whether the DNA in a sample is mostly single or double stranded. With this knowledge, an appropriate library preparation method can be chosen, and the DNA shearing parameters of the protocol can be adjusted to the DNA fragment size in the sample. In this study, we first developed and evaluated a user-friendly fluorometry-based protocol for estimation of DNA strandedness. We also evaluated different capillary electrophoresis methods for estimation of DNA fragmentation levels. Next, we applied the developed methodologies to a broad variety of DNA samples processed with different DNA extraction protocols. Our findings show that both the applied DNA extraction method and the sample type affect the DNA strandedness and fragmentation. The established protocols and the gained knowledge will be applicable for future sequencing-based high-density SNP genotyping in various fields.

Inhibition of transcription leads to rewiring of locus-specific chromatin proteomes.

Transcription of a chromatin template involves the concerted interaction of many different proteins and protein complexes. Analyses of specific factors showed that these interactions change during stress and upon developmental switches. However, how the binding of multiple factors at any given locus is coordinated has been technically challenging to investigate. Here we used Epi-Decoder in yeast to systematically decode, at one transcribed locus, the chromatin binding changes of hundreds of proteins in parallel upon perturbation of transcription. By taking advantage of improved Epi-Decoder libraries, we observed broad rewiring of local chromatin proteomes following chemical inhibition of RNA polymerase. Rapid reduction of RNA polymerase II binding was accompanied by reduced binding of many other core transcription proteins and gain of chromatin remodelers. In quiescent cells, where strong transcriptional repression is induced by physiological signals, eviction of the core transcriptional machinery was accompanied by the appearance of quiescent cell-specific repressors and rewiring of the interactions of protein-folding factors and metabolic enzymes. These results show that Epi-Decoder provides a powerful strategy for capturing the temporal binding dynamics of multiple chromatin proteins under varying conditions and cell states. The systematic and comprehensive delineation of dynamic local chromatin proteomes will greatly aid in uncovering protein-protein relationships and protein functions at the chromatin template.

Development of novel microsatellite markers for population differentiation and detection of natural selection in wild populations of butter catfish, Ompok bimaculatus (Bloch, 1794).

BACKGROUND: Butter catfish (Ompok bimaculatus) is a preferred species in South East Asia, with huge aquaculture potential. However, there is limited information about genetic stock composition due to insufficient markers. The goal of this study was to develop de novo microsatellite markers. METHODS AND RESULTS: For sequencing, genomic SMRT bell libraries (1.5 Kbp size) were prepared for O. bimaculatus. A total of 114 SSR containing sequences were used for primer designing. Polymorphic loci were validated by genotyping 83 individuals from four distant riverine populations, viz., Brahmaputra, Bichiya, Gomti and Kaveri. A total of 30 microsatellite loci were polymorphic, of which five were found to be associated with functional genes and eight (four positive and four negative) loci were found to be under selection pressure. A total of 115 alleles were detected in all loci and PIC ranged from 0.539 to 0.927 and pair-wise FST values from 0.1267 to 0.26002 (p < 0.001), with an overall FST value of 0.17047, indicating the presence of population sub-structure. Cross-species transferability of 29 loci (96.67%) was successful in congener species, Ompok pabda. CONCLUSION: The novel SSR markers developed in this study would facilitate stock characterization of natural populations, to be used in future selection breeding programs and planning conservation strategies in these species. Identified non-neutral markers will give insights into the effect of local adaptation on genetic differentiation in the natural population of this species.

Microsatellite DNA markers for the commercially important fish Hypophthalmus donascimientoi (Pisces: Siluriformes: Pimelodidae) in the Amazon basin: isolation, characterization and amplification in congeneric species.

BACKGROUND: The genus Hypophthalmus comprises six species (H. edentatus, H. marginatus, H. fimbriatus, and H. oremaculatus), and the recently described: H. donascimientoi and H. celiae. The popular name for Hypophthalmus spp. in Brazil is mapará, this name refers to the six species. This group of fish has commercial importance for the states of Amazonas and Pará and, for this reason, requires studies to identify fish stocks. One approach is to use molecular markers, which have been very useful in studies with identification and population analysis of fish. Microsatellite molecular markers (SSRs) are one of the most informative markers for this purpose. There is little populations study of Hypophthalmus using SSRs, and there are less than six loci for the species Hypophthalmus marginatus available in the literature. With the construction of a genomic library of H. donascimientoi, we aimed to isolate and characterize SSRS markers and evaluate the extent of interspecific amplification. METHODS AND RESULTS: A genomic library was constructed with regions enriched of microsatellite for Hypophthalmus donascimientoi. A total of 126 contigs with 42 SSRs were used to design flanking primers for 39 microsatellites. Fifteen loci were characterized in three locations of the Solimões/Amazonas Rivers. The number of alleles ranged from one to 17 with a total of 126 alleles. The mean observed heterozygosity (HO) and expected heterozygosity (HE) were 0.721-0.692, respectively (S.d. HO 0.061 and HE 0.060). Two loci showed significant deviation in the HWE. The PIC ranged from 0.375 to 0.908. Such identified, 12 highly informative loci, and two moderately informative loci. Among the fifteen loci characterized, seven were successfully amplified in four other species of the genus. CONCLUSIONS: The microsatellite showed promise for estimating the genetic variability of H. donascimientoi and can be used as an efficient tool in population analyses of this species and in congeneric species analyzed.

Protocol variations in run-on transcription dataset preparation produce detectable signatures in sequencing libraries.

BACKGROUND: A variety of protocols exist for producing whole genome run-on transcription datasets. However, little is known about how differences between these protocols affect the signal within the resulting libraries. RESULTS: Using run-on transcription datasets generated from the same biological system, we show that a variety of GRO- and PRO-seq preparation methods leave identifiable signatures within each library. Specifically we show that the library preparation method results in differences in quality control metrics, as well as differences in the signal distribution at the 5 ' end of transcribed regions. These shifts lead to disparities in eRNA identification, but do not impact analyses aimed at inferring the key regulators involved in changes to transcription. CONCLUSIONS: Run-on sequencing protocol variations result in technical signatures that can be used to identify both the enrichment and library preparation method of a particular data set. These technical signatures are batch effects that limit detailed comparisons of pausing ratios and eRNAs identified across protocols. However, these batch effects have only limited impact on our ability to infer which regulators underlie the observed transcriptional changes.

Validation of a random Vibrio parahaemolyticus genomic library by selection of quinolone resistance in a heterologous host.

Vibrio parahaemolyticus is a shellfish-borne pathogen that is a highly prevalent causative agent of inflammatory gastroenteritis in humans. Genomic libraries have proven useful for the identification of novel gene functions in many bacterial species. In this study we prepared a library containing 40 kb fragments of randomly sheared V. parahaemolyticus genomic DNA and introduced this into Escherichia coli HB101 using a commercially available low copy cosmid system. In order to estimate coverage and suitability of the library and potentially identify novel antimicrobial resistance determinants, we screened for the acquisition of resistance to the fluoroquinolone norfloxacin - a phenotype exhibited by V. parahaemolyticus but not the heterologous E. coli host. Upon selection on solid medium containing norfloxacin, 0.52% of the library population was resistant, consistent with the selection of a single resistance locus. End-sequencing identified six distinct insert fragments. All clones displayed fourfold increased norfloxacin MIC compared with E. coli HB101 carrying an empty vector. The common locus contained within resistant clones included qnr, a previously described quinolone resistance gene. These results indicate that the library was unbiased, of sufficient coverage and that heterologous expression was possible. While we hope that this library proves useful for identifying the genetic determinants of complex phenotypes such as those related to virulence, not all norfloxacin resistance genes were detected in our screen. As such, we discuss the benefits and limitations of this approach for identifying the genetic basis of uncharacterized bacterial phenotypes.

Deep mutational scanning reveals the structural basis for α-synuclein activity.

Defining the biologically active structures of proteins in their cellular environments remains challenging for proteins with multiple conformations and functions, where only a minor conformer might be associated with a given function. Here, we use deep mutational scanning to probe the structure and dynamics of α-synuclein, a protein known to adopt disordered, helical and amyloid conformations. We examined the effects of 2,600 single-residue substitutions on the ability of intracellularly expressed α-synuclein to slow the growth of yeast. Computational analysis of the data showed that the conformation responsible for this phenotype is a long, uninterrupted, amphiphilic helix with increasing dynamics toward the C terminus. Deep mutational scanning can therefore determine biologically active conformations in cellular environments, even for a highly dynamic multi-conformational protein.

Enzymatic C-H activation of aromatic compounds through CO2 fixation.

The direct C-H carboxylation of aromatic compounds is an attractive route to the corresponding carboxylic acids, but remains challenging under mild conditions. It has been proposed that the first step in anaerobic microbial degradation of recalcitrant aromatic compounds is a UbiD-mediated carboxylation. In this study, we use the UbiD enzyme ferulic acid decarboxylase (Fdc) in combination with a carboxylic acid reductase to create aromatic degradation-inspired cascade reactions, leading to efficient functionalization of styrene through CO2 fixation. We reveal that rational structure-guided laboratory evolution can expand the substrate scope of Fdc, resulting in activity on a range of mono- and bicyclic aromatic compounds through a single mutation. Selected variants demonstrated 150-fold improvement in the conversion of coumarillic acid to benzofuran + CO2 and unlocked reactivity towards naphthoic acid. Our data demonstrate that UbiD-mediated C-H activation is a versatile tool for the transformation of aryl/alkene compounds and CO2 into commodity chemicals.

Pyruvate-responsive genetic circuits for dynamic control of central metabolism.

Dynamic regulation is a promising strategy for fine-tuning metabolic fluxes in microbial cell factories. However, few of these synthetic regulatory systems have been developed for central carbon metabolites. Here we created a set of programmable and bifunctional pyruvate-responsive genetic circuits for dynamic dual control (activation and inhibition) of central metabolism in Bacillus subtilis. We used these genetic circuits to design a feedback loop control system that relies on the intracellular concentration of pyruvate to fine-tune the target metabolic modules, leading to the glucaric acid titer increasing from 207 to 527 mg l-1. The designed logic gate-based circuits were enabled by the characterization of a new antisense transcription mechanism in B. subtilis. In addition, a further increase to 802 mg l-1 was achieved by blocking the formation of by-products. Here, the constructed pyruvate-responsive genetic circuits are presented as effective tools for the dynamic control of central metabolism of microbial cell factories.

Development of highly polymorphic microsatellite markers for Campanula glomerata L. (Campanulaceae).

BACKGROUND: Campanula glomerata L. (Campanulaceae) is a self-incompatible, insect-pollinated herb occurring in calcareous grasslands, and is declining and (critically) endangered in many parts of its European distribution range. It often exists as small and isolated populations. A recovery plan of C. glomerata has been implemented in southern Belgium, involving plant translocations. METHODS AND RESULTS: We developed microsatellite markers using an enriched genomic library and characterized 16 loci in 111 individuals from eight populations. These 16 loci were highly polymorphic, with 11 to 31 alleles per locus for a total of 329 alleles, and expected heterozygosity (He) ranging from 0.470 to 0.938. CONCLUSIONS: These highly polymorphic loci constitute a promising tool for detailed genetic analyses: assigning individuals to distinct multilocus genotypes will allow quantifying pollen dispersal, clonal propagation and sexual recruitment and identifying admixed seed progeny and their pollen donors. Evaluating the genetic status of existing populations and a genetic monitoring of the translocated populations will contribute to optimize success in restoring viable and evolutionary resilient populations.

Microfluidic automated plasmid library enrichment for biosynthetic gene cluster discovery.

Microbial biosynthetic gene clusters are a valuable source of bioactive molecules. However, because they typically represent a small fraction of genomic material in most metagenomic samples, it remains challenging to deeply sequence them. We present an approach to isolate and sequence gene clusters in metagenomic samples using microfluidic automated plasmid library enrichment. Our approach provides deep coverage of the target gene cluster, facilitating reassembly. We demonstrate the approach by isolating and sequencing type I polyketide synthase gene clusters from an Antarctic soil metagenome. Our method promotes the discovery of functional-related genes and biosynthetic pathways.

PINCER: improved CRISPR/Cas9 screening by efficient cleavage at conserved residues.

CRISPR/Cas9 functional genomic screens have emerged as essential tools in drug target discovery. However, the sensitivity of available genome-wide CRISPR libraries is impaired by guides which inefficiently abrogate gene function. While Cas9 cleavage efficiency optimization and essential domain targeting have been developed as independent guide design rationales, no library has yet combined these into a single cohesive strategy to knock out gene function. Here, in a massive reanalysis of CRISPR tiling data using the most comprehensive feature database assembled, we determine which features of guides and their targets best predict activity and how to best combine them into a single guide design algorithm. We present the ProteIN ConsERvation (PINCER) genome-wide CRISPR library, which for the first time combines enzymatic efficiency optimization with conserved length protein region targeting, and also incorporates domains, coding sequence position, U6 termination (TTT), restriction sites, polymorphisms and specificity. Finally, we demonstrate superior performance of the PINCER library compared to alternative genome-wide CRISPR libraries in head-to-head validation. PINCER is available for individual gene knockout and genome-wide screening for both the human and mouse genomes.

Chromatin dynamics enable transcriptional rhythms in the cnidarian Nematostella vectensis.

In animals, circadian rhythms are driven by oscillations in transcription, translation, and proteasomal degradation of highly conserved genes, resulting in diel cycles in the expression of numerous clock-regulated genes. Transcription is largely regulated through the binding of transcription factors to cis-regulatory elements within accessible regions of the chromatin. Chromatin remodeling is linked to circadian regulation in mammals, but it is unknown whether cycles in chromatin accessibility are a general feature of clock-regulated genes throughout evolution. To assess this, we applied an ATAC-seq approach using Nematostella vectensis, grown under two separate light regimes (light:dark (LD) and constant darkness (DD)). Based on previously identified N. vectensis circadian genes, our results show the coupling of chromatin accessibility and circadian transcription rhythmicity under LD conditions. Out of 180 known circadian genes, we were able to list 139 gene promoters that were highly accessible compared to common promoters. Furthermore, under LD conditions, we identified 259 active enhancers as opposed to 333 active enhancers under DD conditions, with 171 enhancers shared between the two treatments. The development of a highly reproducible ATAC-seq protocol integrated with published RNA-seq and ChIP-seq databases revealed the enrichment of transcription factor binding sites (such as C/EBP, homeobox, and MYB), which have not been previously associated with circadian signaling in cnidarians. These results provide new insight into the regulation of cnidarian circadian machinery. Broadly speaking, this supports the notion that the association between chromatin remodeling and circadian regulation arose early in animal evolution as reflected in this non-bilaterian lineage.

Flexible chromosome painting based on multiplex PCR of oligonucleotides and its application for comparative chromosome analyses in Cucumis.

Chromosome painting is a powerful technique for chromosome and genome studies. We developed a flexible chromosome painting technique based on multiplex PCR of a synthetic oligonucleotide (oligo) library in cucumber (Cucumis sativus L., 2n = 14). Each oligo in the library was associated with a universal as well as nested specific primers for amplification, which allow the generation of different probes from the same oligo library. We were also able to generate double-stranded labelled oligos, which produced much stronger signals than single-stranded labelled oligos, by amplification using fluorophore-conjugated primer pairs. Oligos covering cucumber chromosome 1 (Chr1) and chromosome 4 (Chr4) consisting of eight segments were synthesized in one library. Different oligo probes generated from the library painted the corresponding chromosomes/segments unambiguously, especially on pachytene chromosomes. This technique was then applied to study the homoeologous relationships among cucumber, C. hystrix and C. melo chromosomes based on cross-species chromosome painting using Chr4 probes. We demonstrated that the probe was feasible to detect interspecies chromosome homoeologous relationships and chromosomal rearrangement events. Based on its advantages and great convenience, we anticipate that this flexible oligo-painting technique has great potential for the studies of the structure, organization, and evolution of chromosomes in any species with a sequenced genome.

Metagenomic Shotgun Sequencing Analysis of Canalicular Concretions in Lacrimal Canaliculitis Cases.

Lacrimal canaliculitis is a rare infection of the lacrimal canaliculi with canalicular concretions formed by aggregation of organisms. Metagenomic shotgun sequencing analysis using next-generation sequencing has been used to detect pathogens directly from clinical samples. Using this technology, we report cases of successful pathogen detection of canalicular concretions in lacrimal canaliculitis cases. We investigated patients with primary lacrimal canaliculitis examined in the eye clinics of four hospitals from February 2015 to July 2017. Eighteen canalicular concretion specimens collected from 18 eyes of 17 patients were analyzed by shotgun metagenomics sequencing using the MiSeq platform (Illumina). Taxonomic classification was performed using the GenBank NT database. The canalicular concretion diversity was characterized using the Shannon diversity index. This study included 18 eyes (17 patients, 77.1 ± 6.1 years): 82.4% were women with lacrimal canaliculitis; canalicular concretions were obtained from 12 eyes using lacrimal endoscopy and six eyes using canaliculotomy with curettage. Sequencing analysis detected bacteria in all samples (Shannon diversity index, 0.05-1.47). The following genera of anaerobic bacteria (>1% abundance) were identified: Actinomyces spp. in 15 eyes, Propionibacterium spp., Parvimonas spp. in 11 eyes, Prevotella spp. in 9 eyes, Fusobacterium spp. in 6 eyes, Selenomonas spp. in 5 eyes, Aggregatibacter spp. in 3 eyes, facultative and aerobic bacteria such as Streptococcus spp. in 13 eyes, Campylobacter spp. in 6 eyes, and Haemophilus spp. in 3 eyes. The most common combinations were Actinomyces spp. and Streptococcus spp. and Parvinomonas spp. and Streptococcus spp., found in 10 cases. Pathogens were identified successfully using metagenomic shotgun sequencing analysis in patients with canalicular concretions. Canalicular concretions are polymicrobial with anaerobic and facultative, aerobic bacteria.

BtsT/ BtsS is involved in glyoxylate transport in E. coli and its mutations facilitated glyoxylate utilization.

Glyoxylate is an important chemical and is also an intermediate involved in metabolic pathways of living microorganisms. However, it cannot be rapidly utilized by many microbes. We observed a very long lag phase (up to 120 h) when E. coli is growing in a mineral medium supplemented with 50 mM glyoxylate. To better understand this strange growth pattern on glyoxylate and accelerate glyoxylate utilization, a random genomic library of E. coli was transformed into E. coli BW25113, and mutants that showed significantly shortened lag phase on glyoxylate were obtained. Interestingly, mutations in BtsT/BtsS, a two component system that is involved in pyruvate transport, were found to be a common feature in all mutants retrieved. We further demonstrated, through reverse engineering, that the mutations in BtsT/BtsS can indeed increase glyoxylate uptake. Growth experiments with different concentration of glyoxylate also showed the higher the concentration of glyoxylate, the shorter the lag phase. These new findings thus increased our understanding on microbial utilization of glyoxylate.

Genome-wide C-SWAT library for high-throughput yeast genome tagging.

Here we describe a C-SWAT library for high-throughput tagging of Saccharomyces cerevisiae open reading frames (ORFs). In 5,661 strains, we inserted an acceptor module after each ORF that can be efficiently replaced with tags or regulatory elements. We validated the library with targeted sequencing and tagged the proteome with bright fluorescent proteins to quantify the effect of heterologous transcription terminators on protein expression and to localize previously undetected proteins.