Targeted isolation of cloned genomic regions by recombineering for haplotype phasing and isogenic targeting

Universal Targeted Haplotyping by Droplet Digital PCR Sequencing and Its Applications in Noninvasive Prenatal Testing and Pharmacogenetics Analysis

BACKGROUND

The analysis of haplotypes of variants is important for pharmacogenomics analysis and noninvasive prenatal testing for monogenic diseases. However, there is a lack of robust methods for targeted haplotyping.

METHODS

We developed digital PCR haplotype sequencing (dHapSeq) for targeted haplotyping of variants, which is a method that compartmentalizes long DNA molecules into droplets. Within one droplet, 2 target regions are PCR amplified from one template molecule, and their amplicons are fused together. The fused products are then sequenced to determine the phase relationship of the single nucleotide polymorphism (SNP) alleles. The entire haplotype of 10s of SNPs can be deduced after the phase relationship of individual SNPs are determined in a pairwise manner. We applied dHapSeq to noninvasive prenatal testing in 4 families at risk for thalassemia and utilized it to detect NUDT15 diplotypes for predicting drug tolerance in pediatric acute lymphoblastic leukemia (72 cases and 506 controls).

RESULTS

For SNPs within 40 kb, phase relation can be determined with 100% accuracy. In 7 trio families, the haplotyping results for 97 SNPs spanning 185 kb determined by dHapSeq were concordant with the results deduced from the genotypes of both parents and the fetus. In 4 thalassemia families, a 19.3-kb Southeast Asian deletion was successfully phased with 97 downstream SNPs, enabling noninvasive determination of fetal inheritance using relative haplotype dosage analysis. In the NUDT15 analysis, the variant status and phase of the variants were successfully determined in all cases and controls.

CONCLUSIONS

The dHapSeq represents a robust and scalable haplotyping approach with numerous clinical and research applications.

Reassembly of the Biosynthetic Gene Cluster Enables High Epothilone Yield in Engineered Schlegelella brevitalea

Epothilones, as a new class of microtubule-stabilizing anticancer drugs, exhibit strong bioactivity against taxane-resistant cells and show clinical activity for the treatment of advanced breast cancer. Additionally, they also show great potential for a central nervous system injury and Alzheimer's disease. However, due to the long fermentation period of the original producer and challenges of genetic engineering of nonribosomal peptide/polyketide (NRP/PK) megasynthase genes, the application of epothilones is severely limited. Here, we addressed these problems by reassembling a novel 56-kb epothilone biosynthetic gene cluster, optimizing the promoter of each gene based on RNA-seq profiling, and completing precursor synthetic pathways in engineered Schlegella brevitalea. Furthermore, we debottlenecked the cell autolysis by optimizing culture conditions. Finally, the yield of epothilones in shake flasks was improved to 82 mg/L in six-day fermentation. Overall, we not only constructed epothilone overproducers for further drug development but also provided a rational strategy for high-level NRP/PK compound production.

ExoCET: exonuclease in vitro assembly combined with RecET recombination for highly efficient direct DNA cloning from complex genomes

The exponentially increasing volumes of DNA sequence data highlight the need for new DNA cloning methods to explore the new information. Here, we describe ‘ExoCET’ (Exonuclease Combined with RecET recombination) to directly clone any chosen region from bacterial and mammalian genomes with nucleotide precision into operational plasmids. ExoCET combines in vitro exonuclease and annealing with the remarkable capacity of full length RecET homologous recombination (HR) to retrieve specified regions from genomic DNA preparations. Using T4 polymerase (T4pol) as the in vitro exonuclease for ExoCET, we directly cloned large regions (>50 kb) from bacterial and mammalian genomes, including DNA isolated from blood. Employing RecET HR or Cas9 cleavage in vitro, the directly cloned region can be chosen with nucleotide precision to position, for example, a gene into an expression vector without the need for further subcloning. In addition to its utility for bioprospecting in bacterial genomes, ExoCET presents straightforward access to mammalian genomes for various applications such as region-specific DNA sequencing that retains haplotype phasing, the rapid construction of optimal, haplotypic, isogenic targeting constructs or a new way to genotype that presents advantages over Southern blotting or polymerase chain reaction. The direct cloning capacities of ExoCET present new freedoms in recombinant DNA technology.

Revealing complete complex KIR haplotypes phased by long-read sequencing technology

The killer cell immunoglobulin-like receptor (KIR) region of human chromosome 19 contains up to 16 genes for natural killer (NK) cell receptors that recognize human leukocyte antigen (HLA)/peptide complexes and other ligands. The KIR proteins fulfill functional roles in infections, pregnancy, autoimmune diseases and transplantation. However, their characterization remains a constant challenge. Not only are the genes highly homologous due to their recent evolution by tandem duplications, but the region is structurally dynamic due to frequent transposon-mediated recombination. A sequencing approach that precisely captures the complexity of KIR haplotypes for functional annotation is desirable. We present a unique approach to haplotype the KIR loci using single-molecule, real-time (SMRT) sequencing. Using this method, we have-for the first time-comprehensively sequenced and phased sixteen KIR haplotypes from eight individuals without imputation. The information revealed four novel haplotype structures, a novel gene-fusion allele, novel and confirmed insertion/deletion events, a homozygous individual, and overall diversity for the structural haplotypes and their alleles. These KIR haplotypes augment our existing knowledge by providing high-quality references, evolutionary informers, and source material for imputation. The haplotype sequences and gene annotations provide alternative loci for the KIR region in the human genome reference GrCh38.p8.

A Fosmid Pool-Based Next Generation Sequencing Approach to Haplotype-Resolve Whole Genomes

Haplotype resolution of human genomes is essential to describe and interpret genetic variation and its impact on biology and disease. Our approach to haplotyping relies on converting genomic DNA into a fosmid library, which represents the entire diploid genome as a collection of haploid DNA clones of ~40 kb in size. These can be partitioned into pools such that the probability that the same pool contains both parental haplotypes is reduced to ~1 %. This is the key principle of this method, allowing entire pools of fosmids to be massively parallel sequenced, yielding haploid sequence output. Here, we present a detailed protocol for fosmid pool-based next generation sequencing to haplotype-resolve whole genomes including the following steps: (1) generation of high molecular weight DNA fragments of ~40 kb in size from genomic DNA; (2) fosmid cloning and partitioning into 96-well plates; (3) barcoded sequencing library preparation from fosmid pools for next generation sequencing; and (4) computational analysis of fosmid sequences and assembly into contiguous haploid sequences.This method can be used in combination with, but also without, whole genome shotgun sequencing to extensively resolve heterozygous SNPs and structural variants within genomic regions, resulting in haploid contigs of several hundred kb up to several Mb. This method has a broad range of applications including population and ancestry genetics, the clinical interpretation of mutations in personal genomes, the analysis of cancer genomes and highly complex disease gene regions such as MHC. Moreover, haplotype-resolved genome sequencing allows description and interpretation of the diploid nature of genome biology, for example through the analysis of haploid gene forms and allele-specific phenomena. Application of this method has enabled the production of most of the molecular haplotype-resolved genomes reported to date.

Signatures of Crested Ibis MHC Revealed by Recombination Screening and Short-Reads Assembly Strategy

Whole-genome shotgun (WGS) sequencing has become a routine method in genome research over the past decade. However, the assembly of highly polymorphic regions in WGS projects remains a challenge, especially for large genomes. Employing BAC library constructing, PCR screening and Sanger sequencing, traditional strategy is laborious and expensive, which hampers research on polymorphic genomic regions. As one of the most highly polymorphic regions, the major histocompatibility complex (MHC) plays a central role in the adaptive immunity of all jawed vertebrates. In this study, we introduced an efficient procedure based on recombination screening and short-reads assembly. With this procedure, we constructed a high quality 488-kb region of crested ibis MHC that consists of 3 superscaffolds and contains 50 genes. Our sequence showed comparable quality (97.29% identity) to traditional Sanger assembly, while the workload was reduced almost 7 times. Comparative study revealed distinctive features of crested ibis by exhibiting the COL11A2-BLA-BLB-BRD2 cluster and presenting both ADPRH and odorant receptor (OR) gene in the MHC region. Furthermore, the conservation of the BF-TAP1-TAP2 structure in crested ibis and other vertebrate lineages is interesting in light of the hypothesis that coevolution of functionally related genes in the primordial MHC is responsible for the appearance of the antigen presentation pathways at the birth of the adaptive immune system.

RecET direct cloning and Redαβ recombineering of biosynthetic gene clusters, large operons or single genes for heterologous expression

Full-length RecE and RecT from Rac prophage mediate highly efficient linear-linear homologous recombination that can be used to clone large DNA regions directly from genomic DNA into expression vectors, bypassing library construction and screening. Homologous recombination mediated by Redαβ from lambda phage has been widely used for recombinant DNA engineering. Here we present a protocol for direct cloning and engineering of biosynthetic gene clusters, large operons or single genes from genomic DNA using one Escherichia coli host that harbors both RecET and Redαβ systems. The pipeline uses standardized cassettes for horizontal gene transfer options, as well as vectors with different replication origins configured to minimize recombineering background through the use of selectively replicating templates or CcdB counterselection. These optimized reagents and protocols facilitate fast acquisition of transgenes from genomic DNA preparations, which are ready for heterologous expression within 1 week.

Haplotype-resolved genome sequencing: experimental methods and applications

Human genomes are diploid and, for their complete description and interpretation, it is necessary not only to discover the variation they contain but also to arrange it onto chromosomal haplotypes. Although whole-genome sequencing is becoming increasingly routine, nearly all such individual genomes are mostly unresolved with respect to haplotype, particularly for rare alleles, which remain poorly resolved by inferential methods. Here, we review emerging technologies for experimentally resolving (that is, 'phasing') haplotypes across individual whole-genome sequences. We also discuss computational methods relevant to their implementation, metrics for assessing their accuracy and completeness, and the relevance of haplotype information to applications of genome sequencing in research and clinical medicine.

Lambda red mediated gap repair utilizes a novel replicative intermediate in Escherichia coli

The lambda phage Red recombination system can mediate efficient homologous recombination in Escherichia coli, which is the basis of the DNA engineering technique termed recombineering. Red mediated insertion of DNA requires DNA replication, involves a single-stranded DNA intermediate and is more efficient on the lagging strand of the replication fork. Lagging strand recombination has also been postulated to explain the Red mediated repair of gapped plasmids by an Okazaki fragment gap filling model. Here, we demonstrate that gap repair involves a different strand independent mechanism. Gap repair assays examining the strand asymmetry of recombination did not show a lagging strand bias. Directly testing an ssDNA plasmid showed lagging strand recombination is possible but dsDNA plasmids did not employ this mechanism. Insertional recombination combined with gap repair also did not demonstrate preferential lagging strand bias, supporting a different gap repair mechanism. The predominant recombination route involved concerted insertion and subcloning though other routes also operated at lower frequencies. Simultaneous insertion of DNA resulted in modification of both strands and was unaffected by mutations to DNA polymerase I, responsible for Okazaki fragment maturation. The lower efficiency of an alternate Red mediated ends-in recombination pathway and the apparent lack of a Holliday junction intermediate suggested that gap repair does not involve a different Red recombination pathway. Our results may be explained by a novel replicative intermediate in gap repair that does not involve a replication fork. We exploited these observations by developing a new recombineering application based on concerted insertion and gap repair, termed SPI (subcloning plus insertion). SPI selected against empty vector background and selected for correct gap repair recombinants. We used SPI to simultaneously insert up to four different gene cassettes in a single recombineering reaction. Consequently, our findings have important implications for the understanding of E. coli replication and Red recombination.

Direct cloning of isogenic murine DNA in yeast and relevance of isogenicity for targeting in embryonic stem cells

Efficient gene targeting in embryonic stem cells requires that modifying DNA sequences are identical to those in the targeted chromosomal locus. Yet, there is a paucity of isogenic genomic clones for human cell lines and PCR amplification cannot be used in many mutation-sensitive applications. Here, we describe a novel method for the direct cloning of genomic DNA into a targeting vector, pRTVIR, using oligonucleotide-directed homologous recombination in yeast. We demonstrate the applicability of the method by constructing functional targeting vectors for mammalian genes Uhrf1 and Gfap. Whereas the isogenic targeting of the gene Uhrf1 showed a substantial increase in targeting efficiency compared to non-isogenic DNA in mouse E14 cells, E14-derived DNA performed better than the isogenic DNA in JM8 cells for both Uhrf1 and Gfap. Analysis of 70 C57BL/6-derived targeting vectors electroporated in JM8 and E14 cell lines in parallel showed a clear dependence on isogenicity for targeting, but for three genes isogenic DNA was found to be inhibitory. In summary, this study provides a straightforward methodological approach for the direct generation of isogenic gene targeting vectors.

Determination of haplotypes at structurally complex regions using emulsion haplotype fusion PCR

Background

Genotyping and massively-parallel sequencing projects result in a vast amount of diploid data that is only rarely resolved into its constituent haplotypes. It is nevertheless this phased information that is transmitted from one generation to the next and is most directly associated with biological function and the genetic causes of biological effects. Despite progress made in genome-wide sequencing and phasing algorithms and methods, problems assembling (and reconstructing linear haplotypes in) regions of repetitive DNA and structural variation remain. These dynamic and structurally complex regions are often poorly understood from a sequence point of view. Regions such as these that are highly similar in their sequence tend to be collapsed onto the genome assembly. This is turn means downstream determination of the true sequence haplotype in these regions poses a particular challenge. For structurally complex regions, a more focussed approach to assembling haplotypes may be required.

Results

In order to investigate reconstruction of spatial information at structurally complex regions, we have used an emulsion haplotype fusion PCR approach to reproducibly link sequences of up to 1kb in length to allow phasing of multiple variants from neighbouring loci, using allele-specific PCR and sequencing to detect the phase. By using emulsion systems linking flanking regions to amplicons within the CNV, this led to the reconstruction of a 59kb haplotype across the DEFA1A3 CNV in HapMap individuals.

Conclusion

This study has demonstrated a novel use for emulsion haplotype fusion PCR in addressing the issue of reconstructing structural haplotypes at multiallelic copy variable regions, using the DEFA1A3 locus as an example.

A genome-scale resource for in vivo tag-based protein function exploration in C. elegans

Understanding the in vivo dynamics of protein localization and their physical interactions is important for many problems in biology. To enable systematic protein function interrogation in a multicellular context, we built a genome-scale transgenic platform for in vivo expression of fluorescent- and affinity-tagged proteins in Caenorhabditis elegans under endogenous cis regulatory control. The platform combines computer-assisted transgene design, massively parallel DNA engineering, and next-generation sequencing to generate a resource of 14,637 genomic DNA transgenes, which covers 73% of the proteome. The multipurpose tag used allows any protein of interest to be localized in vivo or affinity purified using standard tag-based assays. We illustrate the utility of the resource by systematic chromatin immunopurification and automated 4D imaging, which produced detailed DNA binding and cell/tissue distribution maps for key transcription factor proteins.