Parallel mapping of genotypes to phenotypes contributing to overall biological fitness

Laboratory selection is a powerful approach for engineering new traits in metabolic engineering applications. This approach is limited because determining the genetic basis of improved strains can be difficult using conventional methods. We have recently reported a new method that enables the measurement of fitness for all clones contained within comprehensive genomic libraries, thus enabling the genome-scale mapping of fitness altering genes. Here, we demonstrate a strategy for relating these measurements to the individual phenotypes selected for in a particular environment. We first provide a mathematical framework for decomposing fitness into selectable phenotypes. We then employed this framework to predict that single-batch selections would enrich primarily for library clones with increased growth rate, serial-batch would enrich for a broad collection of clones enhanced via a combination of increased growth rate and/or reduced lag times, and that overlap among selected clones would be minimal. We used the SCalar Analysis of Library Enrichments (SCALEs) method to test these predictions. We mapped all genomic regions for which increased copy number conferred a selective advantage to Escherichia coli when cultured via single- or serial-batch in the presence of 1-naphthol. We identified a surprisingly large collection (163 total) of tolerance regions, including all previously identified solvent tolerance genes in E. coli. We show that the majority of the identified regions were unique to the different selection strategies examined and that such differences were indeed due to differences among enriched clones in growth rate and lag times over the solvent concentrations examined. The combination of a framework for decomposing overall fitness into selectable phenotypes along with a genome-scale method for mapping genes to such phenotypes lays the groundwork for improving the rational design of laboratory selections.     

Engineering improved ethanol production in Escherichia coli with a genome-wide approach

A key challenge to the commercial production of commodity chemical and fuels is the toxicity of such molecules to the microbial host. While a number of studies have attempted to engineer improved tolerance for such compounds, the majority of these studies have been performed in wild-type strains and culturing conditions that differ considerably from production conditions. Here we applied the multiscalar analysis of library enrichments (SCALEs) method and performed a growth selection in an ethanol production system to quantitatively map in parallel all genes in the genome onto ethanol tolerance and production. In order to perform the selection in an ethanol-producing system, we used a previously engineered Escherichia coli ethanol production strain (LW06; ATCC BAA-2466) (Woodruff et al., in press), as the host strain for the multiscalar genomic library analysis (>10⁶ clones for each library of 1, 2, or 4 kb overlapping genomic fragments). By testing individually selected clones, we confirmed that growth selections enriched for clones with both improved ethanol tolerance and production phenotypes. We performed combinatorial testing of the top genes identified (uspC, otsA, otsB) to investigate their ability to confer improved ethanol tolerance or ethanol production. We determined that overexpression of otsA was required for improved tolerance and productivity phenotypes, with the best performing strains showing up to 75% improvement relative to the parent production strain.     

SCALEs: multiscale analysis of library enrichment

We report a genome-wide, multiscale approach to simultaneously measure the effect that the increased copy of each gene and/or operon has on a desired trait or phenotype. The method involves (i) growth selections on a mixture of several different plasmid-based genomic libraries of defined insert sizes or SCALEs, (ii) microarray studies of enriched plasmid DNA, and a (iii) mathematical multiscale analysis that precisely identifies the relevant genetic elements. This approach allows for identification of all single open reading frames and larger multigene fragments within a genomic library that alter the expression of a given phenotype. We have demonstrated this method in Escherichia coli by monitoring, in parallel, a population of >10(6) genomic library clones of different insert sizes, throughout continuous selections over a period of 100 generations.     

Hybridization selections, using immobilized driver DNA, to enrich for clones unique to a library

We have developed a technique which allowed us to isolate sex-specific repeats of chickens (Gallus g. domesticus) from a genomic library originally containing one sex-specific repeat clone per 300 clones. Using this plus/minus selection method, we were able to enrich a sub-library in which the sex-specific repetitive clones made up over half of the population (170-fold increase in representation). This enrichment technique used hybridization kinetics to collect clones which are bound (plus selection) or not bound (minus selection) to their respective driver DNAs immobilized on nitrocellulose paper. Plus selections enriched for repetitive sequences and removed from the library most unique sequences as well as vectors containing no inserted sequences. These plus-selected sub-libraries were enriched for repetitive clones, yet still contained sequences representing as little as 10(-4) of the genome. Minus selection removed repetitive sequences shared between the driver and the library. When a plus-selected sub-library was minus selected, the doubly selected sub-library was enriched in repetitive sequences present in the library but not the minus driver.     

A genomics approach to improve the analysis and design of strain selections

Strain engineering has been traditionally centered on the use of mutation, selection, and screening to develop improved strains. Although mutational and screening methods are well-characterized, selection remains poorly understood. We hypothesized that we could use a genome-wide method for assessing laboratory selections to design selections with enhanced sensitivity (true positives) and specificity (true negatives) towards a single desired phenotype. To test this hypothesis, we first applied multi-SCale Analysis of Library Enrichments (SCALEs) to identify genes conferring increased fitness in continuous flow selections with increasing levels of 3-hydroxypropionic acid (3-HP). We found that this selection not only enriched for 3-HP tolerance phenotypes but also for wall adherence phenotypes (41% false positives). Using this genome-wide data, we designed a serial-batch selection with a decreasing 3-HP gradient. Further examination by ROC analysis confirmed that the serial-batch approach resulted in significantly increased sensitivity (46%) and specificity (10%) for our desired phenotype (3-HP tolerance).     

Rapid and efficient screening of a Representational Difference Analysis library using reverse Southern hybridisation: identification of genetic differences between Haemophilus parasuis isolates

Representational Difference Analysis (RDA) is an established technique used for isolation of specific genetic differences between or within bacterial species. This method was used to investigate the genetic basis of serovar-specificity and the relationship between serovar and virulence in Haemophilus parasuis. An RDA clone library of 96 isolates was constructed using H. parasuis strains H425(P) (serovar 12) and HS1967 (serovar 4). To screen such a large clone library to determine which clones are strain-specific would typically involved separately labelling each clone for use in Southern hybridisation against genomic DNA from each of the strains. In this study, a novel application of reverse Southern hybridisation was used to screen the RDA library: genomic DNA from each strain was labelled and used to probe the library to identify strain-specific clones. This novel approach represents a significant improvement in methodology that is rapid and efficient.     

CpG_MPs: identification of CpG methylation patterns of genomic regions from high-throughput bisulfite sequencing data

High-throughput bisulfite sequencing is widely used to measure cytosine methylation at single-base resolution in eukaryotes. It permits systems-level analysis of genomic methylation patterns associated with gene expression and chromatin structure. However, methods for large-scale identification of methylation patterns from bisulfite sequencing are lacking. We developed a comprehensive tool, CpG_MPs, for identification and analysis of the methylation patterns of genomic regions from bisulfite sequencing data. CpG_MPs first normalizes bisulfite sequencing reads into methylation level of CpGs. Then it identifies unmethylated and methylated regions using the methylation status of neighboring CpGs by hotspot extension algorithm without knowledge of pre-defined regions. Furthermore, the conservatively and differentially methylated regions across paired or multiple samples (cells or tissues) are identified by combining a combinatorial algorithm with Shannon entropy. CpG_MPs identified large amounts of genomic regions with different methylation patterns across five human bisulfite sequencing data during cellular differentiation. Different sequence features and significantly cell-specific methylation patterns were observed. These potentially functional regions form candidate regions for functional analysis of DNA methylation during cellular differentiation. CpG_MPs is the first user-friendly tool for identifying methylation patterns of genomic regions from bisulfite sequencing data, permitting further investigation of the biological functions of genome-scale methylation patterns.     

High throughput direct end sequencing of BAC clones

Libraries constructed in bacterial artificial chromosome (BAC) vectors have become the choice for clone sets in high throughput genomic sequencing projects primarily because of their high stability. BAC libraries have been proposed as a source for minimally over-lapping clones for sequencing large genomic regions, and the use of BAC end sequences (i.e. sequences adjoining the insert sites) has been proposed as a primary means for selecting minimally overlapping clones for sequencing large genomic regions. For this strategy to be effective, high throughput methods for BAC end sequencing of all the clones in deep coverage BAC libraries needed to be developed. Here we describe a low cost, efficient, 96 well procedure for BAC end sequencing. These methods allow us to generate BAC end sequences from human and Arabidoposis libraries with an average read length of >450 bases and with a single pass sequencing average accuracy of >98%. Application of BAC end sequences in genomic sequen-cing is discussed.     

Metagenomic profiling: microarray analysis of an environmental genomic library

Genomic libraries derived from environmental DNA (metagenomic libraries) are useful for characterizing uncultured microorganisms. However, conventional library-screening techniques permit characterization of relatively few environmental clones. Here we describe a novel approach for characterization of a metagenomic library by hybridizing the library with DNA from a set of groundwater isolates, reference strains, and communities. A cosmid library derived from a microcosm of groundwater microorganisms was used to construct a microarray (COSMO) containing approximately 1-kb PCR products amplified from the inserts of 672 cosmids plus a set of 16S ribosomal DNA controls. COSMO was hybridized with Cy5-labeled genomic DNA from each bacterial strain, and the results were compared with the results for a common Cy3-labeled reference DNA sample consisting of a composite of genomic DNA from multiple species. The accuracy of the results was confirmed by the preferential hybridization of each strain to its corresponding rDNA probe. Cosmid clones were identified that hybridized specifically to each of 10 microcosm isolates, and other clones produced positive results with multiple related species, which is indicative of conserved genes. Many clones did not hybridize to any microcosm isolate; however, some of these clones hybridized to community genomic DNA, suggesting that they were derived from microbes that we failed to isolate in pure culture. Based on identification of genes by end sequencing of 17 such clones, DNA could be assigned to functions that have potential ecological importance, including hydrogen oxidation, nitrate reduction, and transposition. Metagenomic profiling offers an effective approach for rapidly characterizing many clones and identifying the clones corresponding to unidentified species of microorganisms.     

Metagenomic Profiling: Microarray Analysis of an Environmental Genomic Library

Genomic libraries derived from environmental DNA (metagenomic libraries) are useful for characterizing uncultured microorganisms. However, conventional library-screening techniques permit characterization of relatively few environmental clones. Here we describe a novel approach for characterization of a metagenomic library by hybridizing the library with DNA from a set of groundwater isolates, reference strains, and communities. A cosmid library derived from a microcosm of groundwater microorganisms was used to construct a microarray (COSMO) containing ~1-kb PCR products amplified from the inserts of 672 cosmids plus a set of 16S ribosomal DNA controls. COSMO was hybridized with Cy5-labeled genomic DNA from each bacterial strain, and the results were compared with the results for a common Cy3-labeled reference DNA sample consisting of a composite of genomic DNA from multiple species. The accuracy of the results was confirmed by the preferential hybridization of each strain to its corresponding rDNA probe. Cosmid clones were identified that hybridized specifically to each of 10 microcosm isolates, and other clones produced positive results with multiple related species, which is indicative of conserved genes. Many clones did not hybridize to any microcosm isolate; however, some of these clones hybridized to community genomic DNA, suggesting that they were derived from microbes that we failed to isolate in pure culture. Based on identification of genes by end sequencing of 17 such clones, DNA could be assigned to functions that have potential ecological importance, including hydrogen oxidation, nitrate reduction, and transposition. Metagenomic profiling offers an effective approach for rapidly characterizing many clones and identifying the clones corresponding to unidentified species of microorganisms.     

用克隆池PCR法筛选小麦-簇毛麦易位系6VS/6AL基因组TAC文库获得抗病基因候选克隆

用根据抗病基因保守区设计的一对简并性引物,从小麦－簇毛麦易位系6VS／6AL cDNA中PCR扩增获得一个具有抗病基因核苷酸结合位点（Nucleotide binding site,NBS）结构特点的DNA片段克隆N7。从小麦－簇毛麦易位系6VS／6AL基因组TAC（Transformation-competent artificial chromosome,TAC）文库的22块96孔板提取所有2112个克隆池（每个池含约1000个克隆）的质粒,再根据N7的核苷酸序列设计一对特异引物,用克隆池PCR（pooled PCR)法经分级筛选从文库中获得一个阳性克隆。以N7为探针,通过Southern杂交证实了该TAC克隆为真正含有抗病候选基因的克隆。研究结果表明克隆池PCR法对克隆数目巨大的基因组文库的筛选很有效。     

Construction and utility of 10-kb libraries for efficient clone-gap closure for rice genome sequencing

Rice is an important crop and a model system for monocot genomics, and is a target for whole genome sequencing by the International Rice Genome Sequencing Project (IRGSP). The IRGSP is using a clone by clone approach to sequence rice based on minimum tiles of BAC or PAC clones. For chromosomes 10 and 3 we are using an integrated physical map based on two fingerprinted and end-sequenced BAC libraries to identifying a minimum tiling path of clones. In this study we constructed and tested two rice genomic libraries with an average insert size of 10 kb (10-kb library) to support the gap closure and finishing phases of the rice genome sequencing project. The HaeIII library contains 166,752 clones covering approximately 4.6x rice genome equivalents with an average insert size of 10.5 kb. The Sau3AI library contains 138,960 clones covering 4.2x genome equivalents with an average insert size of 11.6 kb. Both libraries were gridded in duplicate onto 11 high-density filters in a 5 x 5 pattern to facilitate screening by hybridization. The libraries contain an unbiased coverage of the rice genome with less than 5% contamination by clones containing organelle DNA or no insert. An efficient method was developed, consisting of pooled overgo hybridization, the selection of 10-kb gap spanning clones using end sequences, transposon sequencing and utilization of in silico draft sequence, to close relatively small gaps between sequenced BAC clones. Using this method we were able to close a majority of the gaps (up to approximately 50 kb) identified during the finishing phase of chromosome-10 sequencing. This method represents a useful way to close clone gaps and thus to complete the entire rice genome.     

用于高通量测序的基因组靶序列捕获方法的建立

文章旨在建立一种基因组目标靶序列捕捉文库的方法,并结合第二代测序技术,以实现候选基因区段的深度测序。利用Agilent公司的eArray在线平台,对1250个基因的11824个外显子共2414977bp的基因组序列进行120个碱基长度的捕捉探针(钓饵)设计,并制备成SureSelect液相靶序列捕获试剂。选用2例人基因组DNA,超声打断后末端补平并磷酸化,连接SOLiD接头,回收150bp～200bp的DNA片段,与靶序列探针杂交捕获目标序列,油包水微乳滴PCR扩增后,磁珠分离富集,上SOLiD测序系统通过工作流程分析(WFA)进行文库质量的评价,或正式测序反应。结果显示对所包含的11147个基因外显子片段设计出并合成了46509个捕捉探针,制备成SureSelect试剂盒。探针可有效地捕捉并富集基因组DNA的目标靶片段,定量PCR显示富集效率可达29倍。WFA分析表明文库可以在SOLiD仪器进行正式测序。测序结果显示靶序列区域的测序数占有效总测序数的比例达到70%,覆盖率均在200×以上。结果表明本研究所建立的SureSelect基因组靶序列捕捉、富集建立测序文库的技术路线可行,可直接用于SOLiD测序仪的测序。     

胜利油田单12高温油藏非培养和富集培养样品的细菌组成特性

[目的]通过比较分析油藏样品的微生物群落结构特点,认识油藏微生物的生态功能.[方法]利用3种油藏微生物研究中常用的富集培养方法,对胜利油田单12区块S12-4油井产出水样品进行了选择性富集培养,运用构建16S rRNA基因文库的方法分析了富集样品和非培养样品的细菌多样性.[结果]通过16S rRNA基因序列比对发现,非培养样品、异养菌富集样品、烃降解菌富集样品和硫酸盐还原菌富集样品中的优势菌分别为Pseudomonas属,Thermotoga属,Thermaerobacter属和Thermotoga属的成员.多样性分析结果表明,非培养样品的微生物多样性最丰富,同时非培养样品和富集样品的微生物群落结构存在很大的差异,富集样品中的微生物包括优势菌在油藏原位环境中含量很低.[结论]细菌组成差异的比较结果,对油藏微生物的生态功能研究和微生物驱油潜力评估具有重要意义.     

Combining phage display with de novo protein sequencing for reverse engineering of monoclonal antibodies

The enormous diversity created by gene recombination and somatic hypermutation makes de novo protein sequencing of monoclonal antibodies a uniquely challenging problem. Modern mass spectrometry-based sequencing will rarely, if ever, provide a single unambiguous sequence for the variable domains. A more likely outcome is computation of an ensemble of highly similar sequences that can satisfy the experimental data. This outcome can result in the need for empirical testing of many candidate sequences, sometimes iteratively, to identity one which can replicate the activity of the parental antibody. Here we describe an improved approach to antibody protein sequencing by using phage display technology to generate a combinatorial library of sequences that satisfy the mass spectrometry data, and selecting for functional candidates that bind antigen. This approach was used to reverse engineer 2 commercially-obtained monoclonal antibodies against murine CD137. Proteomic data enabled us to assign the majority of the variable domain sequences, with the exception of 3–5% of the sequence located within or adjacent to complementarity-determining regions. To efficiently resolve the sequence in these regions, small phage-displayed libraries were generated and subjected to antigen binding selection. Following enrichment of antigen-binding clones, 2 clones were selected for each antibody and recombinantly expressed as antigen-binding fragments (Fabs). In both cases, the reverse-engineered Fabs exhibited identical antigen binding affinity, within error, as Fabs produced from the commercial IgGs. This combination of proteomic and protein engineering techniques provides a useful approach to simplifying the technically challenging process of reverse engineering monoclonal antibodies from protein material.     

Construction and characterization of a gridded cattle BAC library

A bovine genomic large-insert bacterial artificial chromosome (BAC) library has been constructed from leukocytes of a Holstein-Friesian male. Size fractionated DpnII-digested genomic DNA was ligated to the dephosphorylated BamH1 ends of a pBACe3.6 vector. Approximately 8.3 x 10(4) individual BAC clones were picked into 384-well plates. Two-hundred and sixty-seven randomly chosen clones were characterized by pulsed-field gel electrophoresis (PFGE). The average insert size was 104 kb with a frequency of clones without inserts of 5.5%. Thirty-four BAC clones were mapped by fluorescence in situ hybridization (FISH) to cattle chromosomes. Three showed signals at more than one location, one of them on the centromeric regions of all autosomes, indicating that the clone contains centromeric repeats. A subset of these BAC clones was used for the development of sequence tagged sites. Both subcloning and direct sequencing of the BACs were used for generating sequence tagged site information. The clones from the library were gridded onto high-density membranes, and PCR superpools were produced from the same set of clones. Membranes and superpools are available through the Resource Centre of the German Human Genome Project in Berlin (http:// www.rzpd.de).     

Enhanced CHO Clone Screening: Application of Targeted Locus Amplification and Next‐Generation Sequencing Technologies for Cell Line Development

Early analytical clone screening is important during Chinese hamster ovary (CHO) cell line development of biotherapeutic proteins to select a clonally derived cell line with most favorable stability and product quality. Sensitive sequence confirmation methods using mass spectrometry have limitations in throughput and turnaround time. Next‐generation sequencing (NGS) technologies emerged as alternatives for CHO clone analytics. We report an efficient NGS workflow applying the targeted locus amplification (TLA) strategy for genomic screening of antibody expressing CHO clones. In contrast to previously reported RNA sequencing approaches, TLA allows for targeted sequencing of genomic integrated transgenic DNA without prior locus information, robust detection of single‐nucleotide variants (SNVs) and transgenic rearrangements. During clone selection, TLA/NGS revealed CHO clones with high‐level SNVs within the antibody gene and we report in another case the utility of TLA/NGS to identify rearrangements at transgenic DNA level. We also determined detection limits for SNVs calling and the potential to identify clone contaminations by TLA/NGS. TLA/NGS also allows to identify genetically identical clones. In summary, we demonstrate that TLA/NGS is a robust screening method useful for routine clone analytics during cell line development with the potential to process up to 24 CHO clones in less than 7 workdays.     

Construction of a human chromosome 3 specific NotI linking library using a novel cloning procedure.

Two new diphasmid vectors (lambda SK17 and SK22) and a novel procedure to construct linking libraries are described. A partial filling-in reaction provides counter-selection against false linking clones in the library, and obviates the need for supF selection. The diphasmid vectors, in combination with the novel selection procedure, have been used to construct a chromosome 3 specific NotI linking library from a human chromosome 3/mouse microcell hybrid cell line. The application of the new vectors and the strong biochemical and biological selections resulted in a library of 60,000 NotI linking clones. As practically all of them are real NotI linking clones (no false recombinants) the library represents approximately 3,000 human recombinants (equal to 10-15 genomic equivalents of chromosome 3). Previously published methods for construction of linking libraries are compared with the procedure described in the present paper. The advantages of the new vectors and the novel protocol are discussed.     

Construction and Identification of Bacterial Artificial Chromosome Library for 0-613-2R in Upland Cotton

A bacterial artificial chromosome （BAC） library containing a large genomlc DNA insert is an important tool for genome physical mapping, map-based cloning, and genome sequencing. To Isolate genes via a map-based cloning strategy and to perform physical mapping of the cotton genome, a high-quality BAC library containing large cotton DNA Inserts Is needed. We have developed a BAC library of the restoring line 0-613-2R for Isolating the fertility restorer （Rf1） gene and genomic research in cotton （Gossypium hirsutum L.）. The BAC library contains 97 825 clones stored In 255 pieces of a 384-well mlcrotiter plate. Random samples of BACs digested with the Notl enzyme Indicated that the average Insert size Is approximately 130 kb, with a range of 80-275 kb, and 95.7% of the BAC clones in the library have an average insert size larger than 100 kb. Based on a cotton genome size of 2 250 Mb, library coverage is 5.7 × haploid genome equivalents. Four clones were selected randomly from the library to determine the stability of the BAC clones. There were no different fingerprints for 0 and 100 generations of each clone digested with Notl and Hlndiii enzymes. Thus, the atabiiity of a single BAC clone can be sustained at iesat for 100 generations. Eight simple sequence repeat （SSR） markers flanking the Rf; gene were chosen to screen the BAC library by pool using PCR method and 25 positive clones were identified with 3.1 positive clones per SSR marker.