Matthew S. Hestand

Polymerase specific error rates and profiles identified by single molecule sequencing.

DNA polymerases have an innate error rate which is polymerase and DNA context specific. Historically the mutational rate and profiles have been measured using a variety of methods, each with their own technical limitations. Here we used the unique properties of single molecule sequencing to evaluate the mutational rate and profiles of six DNA polymerases at the sequence level. In addition to accurately determining mutations in double strands, single molecule sequencing also captures direction specific transversions and transitions through the analysis of heteroduplexes. Not only did the error rates vary, but also the direction specific transitions differed among polymerases.

GBSX: a toolkit for experimental design and demultiplexing genotyping by sequencing experiments

BackgroundMassive parallel sequencing is a powerful tool for variant discovery and genotyping. To reduce costs, sequencing of restriction enzyme based reduced representation libraries can be utilized. This technology is generally referred to as Genotyping By Sequencing (GBS). To deal with GBS experimental design and initial processing specific bioinformatic tools are needed.ResultsGBSX is a package that assists in selecting the appropriate enzyme and the design of compatible in-line barcodes. Post sequencing, it performs optimized demultiplexing using these barcodes to create fastq files per barcode which can easily be plugged into existing variant analysis pipelines. Here we demonstrate the usability of the GBSX toolkit and demonstrate improved in-line barcode demultiplexing and trimming performance compared to existing tools.ConclusionsGBSX provides an easy to use suite of tools for designing and demultiplexing of GBS experiments.

A Distinct Class of Chromoanagenesis Events Characterized by Focal Copy Number Gains.

Chromoanagenesis is the process by which a single catastrophic event creates complex rearrangements confined to a single or a few chromosomes. It is usually characterized by the presence of multiple deletions and/or duplications, as well as by copy neutral rearrangements. In contrast, an array CGH screen of patients with developmental anomalies revealed three patients in which a single chromosome carries from 8 to 11 large copy number gains confined to a single chromosome or chromosomal arm, but the absence of deletions. Subsequent fluorescence in situ hybiridization and massive parallel sequencing revealed the duplicons to be clustered together in distinct locations across the altered chromosomes. Breakpoint junction sequences showed both microhomology and non‐templated insertions of up to 40 bp. Hence, these patients each demonstrate a single altered chromosome of clustered insertional duplications, no deletions, and breakpoint junction sequences showing microhomology and/or non‐templated insertions. These observations are difficult to reconcile with current mechanistic descriptions of chromothripsis and chromoanasynthesis. Therefore, we hypothesize those rearrangements to be of a mechanistically different origin. In addition, we suggest that large untemplated insertional sequences observed at breakpoints are driven by a non‐canonical non‐homologous end joining mechanism.

Single molecule real-time (SMRT) sequencing comes of age: applications and utilities for medical diagnostics

Abstract Short read massive parallel sequencing has emerged as a standard diagnostic tool in the medical setting. However, short read technologies have inherent limitations such as GC bias, difficulties mapping to repetitive elements, trouble discriminating paralogous sequences, and difficulties in phasing alleles. Long read single molecule sequencers resolve these obstacles. Moreover, they offer higher consensus accuracies and can detect epigenetic modifications from native DNA. The first commercially available long read single molecule platform was the RS system based on PacBios single molecule real-time (SMRT) sequencing technology, which has since evolved into their RSII and Sequel systems. Here we capsulize how SMRT sequencing is revolutionizing constitutional, reproductive, cancer, microbial and viral genetic testing.

Detecting AGG Interruptions in Male and Female FMR1 Premutation Carriers by Single-Molecule Sequencing

The FMR1 gene contains an unstable CGG repeat in its 5′ untranslated region. Premutation alleles range between 55 and 200 repeat units and confer a risk for developing fragile X‐associated tremor/ataxia syndrome or fragile X‐associated primary ovarian insufficiency. Furthermore, the premutation allele often expands to a full mutation during female germline transmission giving rise to the fragile X syndrome. The risk for a premutation to expand depends mainly on the number of CGG units and the presence of AGG interruptions in the CGG repeat. Unfortunately, the detection of AGG interruptions is hampered by technical difficulties. Here, we demonstrate that single‐molecule sequencing enables the determination of not only the repeat size, but also the complete repeat sequence including AGG interruptions in male and female alleles with repeats ranging from 45 to 100 CGG units. We envision this method will facilitate research and diagnostic analysis of the FMR1 repeat expansion.

Variant discovery and breakpoint region prediction for studying the human 22q11.2 deletion using BAC clone and whole genome sequencing analysis

Velo-cardio-facial syndrome/DiGeorge syndrome/22q11.2 deletion syndrome (22q11.2DS) is caused by meiotic non-allelic homologous recombination events between flanking low copy repeats termed LCR22A and LCR22D, resulting in a 3 million base pair (Mb) deletion. Due to their complex structure, large size and high sequence identity, genetic variation within LCR22s among different individuals has not been well characterized. In this study, we sequenced 13 BAC clones derived from LCR22A/D and aligned them with 15 previously available BAC sequences to create a new genetic variation map. The thousands of variants identified by this analysis were not uniformly distributed in the two LCR22s. Moreover, shared single nucleotide variants between LCR22A and LCR22D were enriched in the Breakpoint Cluster Region pseudogene (BCRP) block, suggesting the existence of a possible recombination hotspot there. Interestingly, breakpoints for atypical 22q11.2 rearrangements have previously been located to BCRPs To further explore this finding, we carried out in-depth analyses of whole genome sequence (WGS) data from two unrelated probands harbouring a de novo 3Mb 22q11.2 deletion and their normal parents. By focusing primarily on WGS reads uniquely mapped to LCR22A, using the variation map from our BAC analysis to help resolve allele ambiguity, and by performing PCR analysis, we infer that the deletion breakpoints were most likely located near or within the BCRP module. In summary, we found a high degree of sequence variation in LCR22A and LCR22D and a potential recombination breakpoint near or within the BCRP block, providing a starting point for future breakpoint mapping using additional trios.

A catalog of hemizygous variation in 127 22q11 deletion patients

The 22q11.2 deletion syndrome is the most common microdeletion disorder, with wide phenotypic variability. To investigate variation within the non-deleted allele we performed targeted resequencing of the 22q11.2 region for 127 patients, identifying multiple deletion sizes, including two deletions with atypical breakpoints. We cataloged ~12,000 hemizygous variant positions, of which 84% were previously annotated. Within the coding regions 95 non-synonymous variants, three stop gains, and two frameshift insertions were identified, some of which we speculate could contribute to atypical phenotypes. We also catalog tolerability of 22q11 gene mutations based on related autosomal recessive disorders in man, embryonic lethality in mice, cross-species conservation and observations that some genes harbor more or less variants than expected. This extensive catalog of hemizygous variants will serve as a blueprint for future experiments to correlate 22q11DS variation with phenotype.

Nested Inversion Polymorphisms Predispose Chromosome 22q11.2 to Meiotic Rearrangements.

Inversion polymorphisms between low-copy repeats (LCRs) might predispose chromosomes to meiotic non-allelic homologous recombination (NAHR) events and thus lead to genomic disorders. However, for the 22q11.2 deletion syndrome (22q11.2DS), the most common genomic disorder, no such inversions have been uncovered as of yet. Using fiber-FISH, we demonstrate that parents transmitting the de novo 3 Mb LCR22A-D 22q11.2 deletion, the reciprocal duplication, and the smaller 1.5 Mb LCR22A-B 22q11.2 deletion carry inversions of LCR22B-D or LCR22C-D. Hence, the inversions predispose chromosome 22q11.2 to meiotic rearrangements and increase the individual risk for transmitting rearrangements. Interestingly, the inversions are nested or flanking rather than coinciding with the deletion or duplication sizes. This finding raises the possibility that inversions are a prerequisite not only for 22q11.2 rearrangements but also for all NAHR-mediated genomic disorders.

EquCab3, an Updated Reference Genome for the Domestic Horse

EquCab2, a high-quality reference genome for the domestic horse, was released in 2007. Since then, it has served as the foundation for nearly all genomic work done in equids. Recent advances in genomic sequencing technology and computational assembly methods have allowed scientists to improve reference assemblies of large animal and plant genomes in terms of contiguity and composition. In 2014, the equine genomics research community began a project to improve the reference sequence for the horse, building upon the solid foundation of EquCab2 and incorporating new short-read data, long-read data, and proximity ligation data. The result, EquCab3, is presented here. The count of non-N bases in the incorporated chromosomes is improved from 2.33Gb in EquCab2 to 2.41Gb from EquCab3. Contiguity has also been improved nearly 40-fold with a contig N50 of 4.5Mb and scaffold contiguity enhanced to where all but one of the 32 chromosomes is comprised of a single scaffold.

Deletion size analysis of 1680 22q11.2DS subjects identifies a new recombination hotspot on chromosome 22q11.2

Recurrent, de novo, meiotic non-allelic homologous recombination events between low copy repeats, termed LCR22s, leads to the 22q11.2 deletion syndrome (22q11.2DS; velo-cardio-facial syndrome/DiGeorge syndrome). Although most 22q11.2DS patients have a similar sized 3 million base pair (Mb), LCR22A-D deletion, some have nested LCR22A-B or LCR22A-C deletions. Our goal is to identify additional recurrent 22q11.2 deletions associated with 22q11.2DS, serving as recombination hotspots for meiotic chromosomal rearrangements. Here, using data from Affymetrix 6.0 microarrays on 1680 22q11.2DS subjects, we identified what appeared to be a nested proximal 22q11.2 deletion in 38 (2.3%) of them. Using molecular and haplotype analyses from 14 subjects and their parent(s) with available DNA, we found essentially three types of scenarios to explain this observation. In eight subjects, the proximal breakpoints occurred in a small sized 12 kb LCR distal to LCR22A, referred to LCR22A+, resulting in LCR22A+-B or LCR22A+-D deletions. Six of these eight subjects had a nested 22q11.2 deletion that occurred during meiosis in a parent carrying a benign 0.2 Mb duplication of the LCR22A-LCR22A+ region with a breakpoint in LCR22A+. Another six had a typical de novo LCR22A-D deletion on one allele and inherited the LCR22A-A+ duplication from the other parent thus appearing on microarrays to have a nested deletion. LCR22A+ maps to an evolutionary breakpoint between mice and humans and appears to serve as a local hotspot for chromosome rearrangements on 22q11.2.