Genetic Heterogeneity of Induced Pluripotent Stem Cells: Results from 24 Clones Derived from a Single C57BL/6 Mouse

Induced pluripotent stem cells (iPSCs) have tremendous potential as a tool for disease modeling, drug testing, and other applications. Since the generation of iPSCs “captures” the genetic history of the individual cell that was reprogrammed, iPSC clones (even those derived from the same individual) would be expected to demonstrate genetic heterogeneity. To assess the degree of genetic heterogeneity, and to determine whether some cells are more genetically “fit” for reprogramming, we performed exome sequencing on 24 mouse iPSC clones derived from skin fibroblasts obtained from two different sites of the same 8-week-old C57BL/6J male mouse. While no differences in the coding regions were detected in the two parental fibroblast pools, each clone had a unique genetic signature with a wide range of heterogeneity observed among the individual clones: a total of 383 iPSC variants were validated for the 24 clones (mean 16.0/clone, range 0–45). Since these variants were all present in the vast majority of the cells in each clone (variant allele frequencies of 40–60% for heterozygous variants), they most likely preexisted in the individual cells that were reprogrammed, rather than being acquired during reprogramming or cell passaging. We then tested whether this genetic heterogeneity had functional consequences for hematopoietic development by generating hematopoietic progenitors in vitro and enumerating colony forming units (CFUs). While there was a range of hematopoietic potentials among the 24 clones, only one clone failed to differentiate into hematopoietic cells; however, it was able to form a teratoma, proving its pluripotent nature. Further, no specific association was found between the mutational spectrum and the hematopoietic potential of each iPSC clone. These data clearly highlight the genetic heterogeneity present within individual fibroblasts that is captured by iPSC generation, and suggest that most of the changes are random, and functionally benign.     

Rapid interactome profiling by massive sequencing

We have developed a high-throughput protein expression and interaction analysis platform that combines cDNA phage display library selection and massive gene sequencing using the 454 platform. A phage display library of open reading frame (ORF) fragments was created from mRNA derived from different tissues. This was used to study the interaction network of the enzyme transglutaminase 2 (TG2), a multifunctional enzyme involved in the regulation of cell growth, differentiation and apoptosis, associated with many different pathologies. After two rounds of panning with TG2 we assayed the frequency of ORFs within the selected phage population using 454 sequencing. Ranking and analysis of more than 120 000 sequences allowed us to identify several potential interactors, which were subsequently confirmed in functional assays. Within the identified clones, three had been previously described as interacting proteins (fibronectin, SMOC1 and GSTO2), while all the others were new. When compared with standard systems, such as microtiter enzyme-linked immunosorbant assay, the method described here is dramatically faster and yields far more information about the interaction under study, allowing better characterization of complex systems. For example, in the case of fibronectin, it was possible to identify the specific domains involved in the interaction.     

Massively parallel tag sequencing reveals the complexity of anaerobic marine protistan communities

Background  

Recent advances in sequencing strategies make possible unprecedented depth and scale of sampling for molecular detection of microbial diversity. Two major paradigm-shifting discoveries include the detection of bacterial diversity that is one to two orders of magnitude greater than previous estimates, and the discovery of an exciting 'rare biosphere' of molecular signatures ('species') of poorly understood ecological significance. We applied a high-throughput parallel tag sequencing (454 sequencing) protocol adopted for eukaryotes to investigate protistan community complexity in two contrasting anoxic marine ecosystems (Framvaren Fjord, Norway; Cariaco deep-sea basin, Venezuela). Both sampling sites have previously been scrutinized for protistan diversity by traditional clone library construction and Sanger sequencing. By comparing these clone library data with 454 amplicon library data, we assess the efficiency of high-throughput tag sequencing strategies. We here present a novel, highly conservative bioinformatic analysis pipeline for the processing of large tag sequence data sets.     

Non-gridded library: a new approach for BAC (bacterial artificial chromosome) exploitation in hexaploid wheat (Triticum aestivum)

The feasibility of exploiting non-gridded bacterial artificial chromosome (BAC) libraries and some major factors affecting the efficiency of handling such libraries were studied in hexaploid wheat. Even for a bacterial culture containing only 55% recombinants, some 2000 BAC clones with inserts ranging from 45 to 245 kb could be pooled. The pooled BAC clones could be amplified by culturing for up to 6 h without losing any target clones. These results imply that even for hexaploid wheat, which has an extremely large genome, some 250 pools are sufficient for a BAC library that should satisfy many research objectives. This non-gridded strategy would dramatically reduce the cost and make robotic equipment non-essential in exploiting BAC technology. To construct a representative library and to minimise clone competition, thawing and re-freezing ligation mixtures and bacterial cultures should be avoided in BAC library construction and application.     

High consistency between replicate 454 pyrosequencing analyses of ectomycorrhizal plant root samples

In this methodological study, we compare 454 sequencing and a conventional cloning and Sanger sequencing approach in their ability to characterize fungal communities PCR amplified from four root systems of the ectomycorrhizal plant Bistorta vivipara. To examine variation introduced by stochastic processes during the laboratory work, we replicated all analyses using two independently obtained DNA extractions from the same root systems. The ITS1 region was used as DNA barcode and the sequences were clustered into OTUs as proxies for species using single linkage clustering (BLASTClust) and 97% sequence similarity cut-off. A relatively low overlap in fungal OTUs was observed between the 454 and the clone library datasets — even among the most abundant OTUs. In a non-metric multidimensional scaling analysis, the samples grouped more according to methodology compared to plant. Some OTUs frequently detected by 454, most notably those OTUs with taxonomic affinity to Glomales, were not detected in the Sanger dataset. Likewise, a few OTUs, including Cenococcum sp., only appeared in the clone libraries. Surprisingly, we observed a significant relationship between GC/AT content of the OTUs and their proportional abundances in the 454 versus the clone library datasets. Reassuringly, a very good consistency in OTU recovery was observed between replicate runs of both sequencing methods. This indicates that stochastic processes had little impact when applying the same sequencing technique on replicate samples.     

Mitochondrial DNA Variant Discovery and Evaluation in Human Cardiomyopathies through Next-Generation Sequencing

Mutations in mitochondrial DNA (mtDNA) may cause maternally-inherited cardiomyopathy and heart failure. In homoplasmy all mtDNA copies contain the mutation. In heteroplasmy there is a mixture of normal and mutant copies of mtDNA. The clinical phenotype of an affected individual depends on the type of genetic defect and the ratios of mutant and normal mtDNA in affected tissues. We aimed at determining the sensitivity of next-generation sequencing compared to Sanger sequencing for mutation detection in patients with mitochondrial cardiomyopathy. We studied 18 patients with mitochondrial cardiomyopathy and two with suspected mitochondrial disease. We “shotgun” sequenced PCR-amplified mtDNA and multiplexed using a single run on Roche''s 454 Genome Sequencer. By mapping to the reference sequence, we obtained 1,300× average coverage per case and identified high-confidence variants. By comparing these to >400 mtDNA substitution variants detected by Sanger, we found 98% concordance in variant detection. Simulation studies showed that >95% of the homoplasmic variants were detected at a minimum sequence coverage of 20× while heteroplasmic variants required >200× coverage. Several Sanger “misses” were detected by 454 sequencing. These included the novel heteroplasmic 7501T>C in tRNA serine 1 in a patient with sudden cardiac death. These results support a potential role of next-generation sequencing in the discovery of novel mtDNA variants with heteroplasmy below the level reliably detected with Sanger sequencing. We hope that this will assist in the identification of mtDNA mutations and key genetic determinants for cardiomyopathy and mitochondrial disease.     

Use of Subtractive Hybridization To Design Habitat-Based Oligonucleotide Probes for Investigation of Natural Bacterial Communities

We describe a rapid oligonucleotide probe design strategy based on subtractive hybridization which yields probes for 16S rRNA or rRNA genes of individual members of microbial communities that are specific within the context of those communities. This strategy circumvents the need to sequence many similar or identical clones of dominant members of a community. Radioactively labeled subfragments of a cloned 16S rRNA gene sequence for which a probe is required (target) were hybridized with biotinylated total 16S ribosomal DNA (rDNA) amplified from the microbial community, and the hybrids formed were subsequently discarded. The remaining enriched fragments were used to screen a library consisting of cloned subfragments of the target sequence by colony hybridization in order to identify the variable regions of the 16S rRNA gene with the required specificity. The sequencing of random clones in one 16S rDNA library demonstrated that only those clones with 100% sequence identity with the probe fragment were detected by it. Moreover, sequencing of other, randomly selected, probe-positive clones revealed 100% sequence identity with the probe. Probes developed in this way tended to correspond to more variable regions of the 16S rRNA if the target sequences were similar to the sequences of other clones in the library and to less variable regions if the target sequences were phylogenetically isolated within the clone library. Although the absolute specificity of the latter probes, as assessed by comparison with available database sequences, was lower than the absolute specificity of the probes from the more variable regions, they were specific within the context of the environmental samples from which they were derived.     

Genome-scale analysis of anti-metabolite directed strain engineering

Classic strain engineering methods have previously been limited by the low-throughput of conventional sequencing technology. Here, we applied a new genomics technology, scalar analysis of library enrichments (SCALEs), to measure >3 million Escherichia coli genomic library clone enrichment patterns resulting from growth selections employing three aspartic-acid anti-metabolites. Our objective was to assess the extent to which access to genome-scale enrichment patterns would provide strain-engineering insights not reasonably accessible through the use of conventional sequencing. We determined that the SCALEs method identified a surprisingly large range of anti-metabolite tolerance regions (423, 865, or 909 regions for each of the three anti-metabolites) when compared to the number of regions (1-3 regions) indicated by conventional sequencing. Genome-scale methods uniquely enable the calculation of clone fitness values by providing concentration data for all clones within a genomic library before and after a period of selection. We observed that clone fitness values differ substantially from clone concentration values and that this is due to differences in overall clone fitness distributions for each selection. Finally, we show that many of the clones of highest fitness overlapped across all selections, suggesting that inhibition of aspartate metabolism, as opposed to specific inhibited enzymes, dominated each selection. Our follow up studies confirmed our observed growth phenotypes and showed that intracellular amino-acid levels were also altered in several of the identified clones. These results demonstrate that genome-scale methods, such as SCALEs, can be used to dramatically improve understanding of classic strain engineering approaches.     

An Anchored Framework BAC Map of Mouse Chromosome 11 Assembled Using Multiplex Oligonucleotide Hybridization

Despite abundant library resources for many organisms, physical mapping of these organisms has been seriously limited due to lack of efficient library screening techniques. We have developed a highly efficient strategy for large-scale screening of genomic libraries based on multiplex oligonucleotide hybridization on high-density genomic filters. We have applied this strategy to generate a bacterial artificial chromosome (BAC) anchored map of mouse chromosome 11. Using the MIT mouse SSLP data, 320 pairs of oligonucleotide probes were designed with an “overgo” computer program that selects new primer sequences that avoid the microsatellite repeat. BACs identified by these probes are automatically anchored to the chromosome. Ninety-two percent of the probes identified positive clones from a 5.9-fold coverage mouse BAC library with an average of 7 positive clones per marker. An average of 4.2 clones was confirmed for 204 markers by PCR. Our data show that a large number of clones can be efficiently isolated from a large genomic library using this strategy with minimal effort. This strategy will have wide application for large-scale mapping and sequencing of human and other large genomes.     

Recombinant Environmental Libraries Provide Access to Microbial Diversity for Drug Discovery from Natural Products

To further explore possible avenues for accessing microbial biodiversity for drug discovery from natural products, we constructed and screened a 5,000-clone “shotgun” environmental DNA library by using an Escherichia coli-Streptomyces lividans shuttle cosmid vector and DNA inserts from microbes derived directly (without cultivation) from soil. The library was analyzed by several means to assess diversity, genetic content, and expression of heterologous genes in both expression hosts. We found that the phylogenetic content of the DNA library was extremely diverse, representing mostly microorganisms that have not been described previously. The library was screened by PCR for sequences similar to parts of type I polyketide synthase genes and tested for the expression of new molecules by screening of live colonies and cell extracts. The results revealed new polyketide synthase genes in at least eight clones. In addition, at least five additional clones were confirmed by high-pressure liquid chromatography analysis and/or biological activity to produce heterologous molecules. These data reinforce the idea that exploiting previously unknown or uncultivated microorganisms for the discovery of novel natural products has potential value and, most importantly, suggest a strategy for developing this technology into a realistic and effective drug discovery tool.     

From micrograms to picograms: quantitative PCR reduces the material demands of high-throughput sequencing

Current efforts to recover the Neandertal and mammoth genomes by 454 DNA sequencing demonstrate the sensitivity of this technology. However, routine 454 sequencing applications still require microgram quantities of initial material. This is due to a lack of effective methods for quantifying 454 sequencing libraries, necessitating expensive and labour-intensive procedures when sequencing ancient DNA and other poor DNA samples. Here we report a 454 sequencing library quantification method based on quantitative PCR that effectively eliminates these limitations. We estimated both the molecule numbers and the fragment size distributions in sequencing libraries derived from Neandertal DNA extracts, SAGE ditags and bonobo genomic DNA, obtaining optimal sequencing yields without performing any titration runs. Using this method, 454 sequencing can routinely be performed from as little as 50 pg of initial material without titration runs, thereby drastically reducing costs while increasing the scope of sample throughput and protocol development on the 454 platform. The method should also apply to Illumina/Solexa and ABI/SOLiD sequencing, and should therefore help to widen the accessibility of all three platforms.     

Rapid Identification of Sequences for Orphan Enzymes to Power Accurate Protein Annotation

The power of genome sequencing depends on the ability to understand what those genes and their proteins products actually do. The automated methods used to assign functions to putative proteins in newly sequenced organisms are limited by the size of our library of proteins with both known function and sequence. Unfortunately this library grows slowly, lagging well behind the rapid increase in novel protein sequences produced by modern genome sequencing methods. One potential source for rapidly expanding this functional library is the “back catalog” of enzymology – “orphan enzymes,” those enzymes that have been characterized and yet lack any associated sequence. There are hundreds of orphan enzymes in the Enzyme Commission (EC) database alone. In this study, we demonstrate how this orphan enzyme “back catalog” is a fertile source for rapidly advancing the state of protein annotation. Starting from three orphan enzyme samples, we applied mass-spectrometry based analysis and computational methods (including sequence similarity networks, sequence and structural alignments, and operon context analysis) to rapidly identify the specific sequence for each orphan while avoiding the most time- and labor-intensive aspects of typical sequence identifications. We then used these three new sequences to more accurately predict the catalytic function of 385 previously uncharacterized or misannotated proteins. We expect that this kind of rapid sequence identification could be efficiently applied on a larger scale to make enzymology’s “back catalog” another powerful tool to drive accurate genome annotation.     

A low-cost,high-throughput microfluidic nano-culture platform for functional metagenomics

Functional metagenomics is an attractive culture-independent approach for functional screening of diverse microbiomes to identify known and novel genes. Since functional screening can involve sifting through tens of thousands of metagenomic library clones, an easy high-throughput screening approach is desirable. Here, we demonstrate a proof-of-concept application of a low-cost, high-throughput droplet based microfluidic assay to the selection of antibiotic resistance genes from a soil metagenomic library. Metagenomic library members encapsulated in nanoliter volume water-in-oil droplets were printed on glass slides robotically, and cell growth in individual drops in the presence of ampicillin was imaged and quantified to identify ampicillin-resistant clones. From the hits, true positives were confirmed by sequencing and functional validation. The ease of liquid handling, ease of set-up, low cost, and robust workflow makes the droplet-based nano-culture platform a promising candidate for screening and selection assays for functional metagenomic libraries.     

Engineering Bacteria to Search for Specific Concentrations of Molecules by a Systematic Synthetic Biology Design Method

Bacteria navigate environments full of various chemicals to seek favorable places for survival by controlling the flagella’s rotation using a complicated signal transduction pathway. By influencing the pathway, bacteria can be engineered to search for specific molecules, which has great potential for application to biomedicine and bioremediation. In this study, genetic circuits were constructed to make bacteria search for a specific molecule at particular concentrations in their environment through a synthetic biology method. In addition, by replacing the “brake component” in the synthetic circuit with some specific sensitivities, the bacteria can be engineered to locate areas containing specific concentrations of the molecule. Measured by the swarm assay qualitatively and microfluidic techniques quantitatively, the characteristics of each “brake component” were identified and represented by a mathematical model. Furthermore, we established another mathematical model to anticipate the characteristics of the “brake component”. Based on this model, an abundant component library can be established to provide adequate component selection for different searching conditions without identifying all components individually. Finally, a systematic design procedure was proposed. Following this systematic procedure, one can design a genetic circuit for bacteria to rapidly search for and locate different concentrations of particular molecules by selecting the most adequate “brake component” in the library. Moreover, following simple procedures, one can also establish an exclusive component library suitable for other cultivated environments, promoter systems, or bacterial strains.     

Successful development of microsatellite markers in a challenging species: the horizontal borer Austroplatypus incompertus (Coleoptera: Curculionidae)

The analysis of microsatellite loci has allowed significant advances in evolutionary biology and pest management. However, until very recently, the potential benefits have been compromised by the high costs of developing these neutral markers. High-throughput sequencing provides a solution to this problem. We describe the development of 13 microsatellite markers for the eusocial ambrosia beetle, Austroplatypus incompertus, a significant pest of forests in southeast Australia. The frequency of microsatellite repeats in the genome of A. incompertus was determined to be low, and previous attempts at microsatellite isolation using a traditional genomic library were problematic. Here, we utilised two protocols, microsatellite-enriched genomic library construction and high-throughput 454 sequencing and characterised 13 loci which were polymorphic among 32 individuals. Numbers of alleles per locus ranged from 2 to 17, and observed and expected heterozygosities from 0.344 to 0.767 and from 0.507 to 0.860, respectively. These microsatellites have the resolution required to analyse fine-scale colony and population genetic structure. Our work demonstrates the utility of next-generation 454 sequencing as a method for rapid and cost-effective acquisition of microsatellites where other techniques have failed, or for taxa where marker development has historically been both complicated and expensive.     

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.     

Intracellular Screen To Identify Metagenomic Clones That Induce or Inhibit a Quorum-Sensing Biosensor

The goal of this study was to design and evaluate a rapid screen to identify metagenomic clones that produce biologically active small molecules. We built metagenomic libraries with DNA from soil on the floodplain of the Tanana River in Alaska. We extracted DNA directly from the soil and cloned it into fosmid and bacterial artificial chromosome vectors, constructing eight metagenomic libraries that contain 53,000 clones with inserts ranging from 1 to 190 kb. To identify clones of interest, we designed a high throughput “intracellular” screen, designated METREX, in which metagenomic DNA is in a host cell containing a biosensor for compounds that induce bacterial quorum sensing. If the metagenomic clone produces a quorum-sensing inducer, the cell produces green fluorescent protein (GFP) and can be identified by fluorescence microscopy or captured by fluorescence-activated cell sorting. Our initial screen identified 11 clones that induce and two that inhibit expression of GFP. The intracellular screen detected quorum-sensing inducers among metagenomic clones that a traditional overlay screen would not. One inducing clone carries a LuxI homologue that directs the synthesis of an N-acyl homoserine lactone quorum-sensing signal molecule. The LuxI homologue has 62% amino acid sequence identity to its closest match in GenBank, AmfI from Pseudomonas fluorescens, and is on a 78-kb insert that contains 67 open reading frames. Another inducing clone carries a gene with homology to homocitrate synthase. Our results demonstrate the power of an intracellular screen to identify functionally active clones and biologically active small molecules in metagenomic libraries.     

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.     

Microbial Diversity in Maras Salterns, a Hypersaline Environment in the Peruvian Andes

Maras salterns are located 3,380 m above sea level in the Peruvian Andes. These salterns consist of more than 3,000 little ponds which are not interconnected and act as crystallizers where salt precipitates. These ponds are fed by hypersaline spring water rich in sodium and chloride. The microbiota inhabiting these salterns was examined by fluorescence in situ hybridization (FISH), 16S rRNA gene clone library analysis, and cultivation techniques. The total counts per milliliter in the ponds were around 2 × 10⁶ to 3 × 10⁶ cells/ml, while the spring water contained less than 100 cells/ml and did not yield any detectable FISH signal. The microbiota inhabiting the ponds was dominated (80 to 86% of the total counts) by Archaea, while Bacteria accounted for 10 to 13% of the 4′,6′-diamidino-2-phenylindole (DAPI) counts. A total of 239 16S rRNA gene clones were analyzed (132 Archaea clones and 107 Bacteria clones). According to the clone libraries, the archaeal assemblage was dominated by microorganisms related to the cosmopolitan square archaeon “Haloquadra walsbyi,” although a substantial number of the sequences in the libraries (31% of the 16S rRNA gene archaeal clones) were related to Halobacterium sp., which is not normally found in clone libraries from solar salterns. All the bacterial clones were closely related to each other and to the γ-proteobacterium “Pseudomonas halophila” DSM 3050. FISH analysis with a probe specific for this bacterial assemblage revealed that it accounted for 69 to 76% of the total bacterial counts detected with a Bacteria-specific probe. When pond water was used to inoculate solid media containing 25% total salts, both extremely halophilic Archaea and Bacteria were isolated. Archaeal isolates were not related to the isolates in clone libraries, although several bacterial isolates were very closely related to the “P. halophila” cluster found in the libraries. As observed for other hypersaline environments, extremely halophilic bacteria that had ecological relevance seemed to be easier to culture than their archaeal counterparts.