Reference Grade Characterization of Polymorphisms in Full-Length HLA Class I and II Genes With Short-Read Sequencing on the ION PGM System and Long-Reads Generated by Single Molecule, Real-Time Sequencing on the PacBio Platform.

Although NGS technologies fuel advances in high-throughput HLA genotyping methods for identification and classification of HLA genes to assist with precision medicine efforts in disease and transplantation, the efficiency of these methods are impeded by the absence of adequately-characterized high-frequency HLA allele reference sequence databases for the highly polymorphic HLA gene system. Here, we report on producing a comprehensive collection of full-length HLA allele sequences for eight classical HLA loci found in the Japanese population. We augmented the second-generation short read data generated by the Ion Torrent technology with long amplicon spanning consensus reads delivered by the third-generation SMRT sequencing method to create reference grade high-quality sequences of HLA class I and II gene alleles resolved at the genomic coding and non-coding level. Forty-six DNAs were obtained from a reference set used previously to establish the HLA allele frequency data in Japanese subjects. The samples included alleles with a collective allele frequency in the Japanese population of more than 99.2%. The HLA loci were independently amplified by long-range PCR using previously designed HLA-locus specific primers and subsequently sequenced using SMRT and Ion PGM sequencers. The mapped long and short-reads were used to produce a reference library of consensus HLA allelic sequences with the help of the reference-aware software tool LAA for SMRT Sequencing. A total of 253 distinct alleles were determined for 46 healthy subjects. Of them, 137 were novel alleles: 101 SNVs and/or indels and 36 extended alleles at a partial or full-length level. Comparing the HLA sequences from the perspective of nucleotide diversity revealed that HLA-DRB1 was the most divergent among the eight HLA genes, and that the HLA-DPB1 gene sequences diverged into two distinct groups, DP2 and DP5, with evidence of independent polymorphisms generated in exon 2. We also identified two specific intronic variations in HLA-DRB1 that might be involved in rheumatoid arthritis. In conclusion, full-length HLA allele sequencing by third-generation and second-generation technologies has provided polymorphic gene reference sequences at a genomic allelic resolution including allelic variations assigned up to the field-4 level for a stronger foundation in precision medicine and HLA-related disease and transplantation studies.

HLA Typing for the Next Generation.

Allele-level resolution data at primary HLA typing is the ideal for most histocompatibility testing laboratories. Many high-throughput molecular HLA typing approaches are unable to determine the phase of observed DNA sequence polymorphisms, leading to ambiguous results. The use of higher resolution methods is often restricted due to cost and time limitations. Here we report on the feasibility of using Pacific Biosciences' Single Molecule Real-Time (SMRT) DNA sequencing technology for high-resolution and high-throughput HLA typing. Seven DNA samples were typed for HLA-A, -B and -C. The results showed that SMRT DNA sequencing technology was able to generate sequences that spanned entire HLA Class I genes that allowed for accurate allele calling. Eight novel genomic HLA class I sequences were identified, four were novel alleles, three were confirmed as genomic sequence extensions and one corrected an existing genomic reference sequence. This method has the potential to revolutionize the field of HLA typing. The clinical impact of achieving this level of resolution HLA typing data is likely to considerable, particularly in applications such as organ and blood stem cell transplantation where matching donors and recipients for their HLA is of utmost importance.

Reconstructing complex regions of genomes using long-read sequencing technology.

Obtaining high-quality sequence continuity of complex regions of recent segmental duplication remains one of the major challenges of finishing genome assemblies. In the human and mouse genomes, this was achieved by targeting large-insert clones using costly and laborious capillary-based sequencing approaches. Sanger shotgun sequencing of clone inserts, however, has now been largely abandoned, leaving most of these regions unresolved in newer genome assemblies generated primarily by next-generation sequencing hybrid approaches. Here we show that it is possible to resolve regions that are complex in a genome-wide context but simple in isolation for a fraction of the time and cost of traditional methods using long-read single molecule, real-time (SMRT) sequencing and assembly technology from Pacific Biosciences (PacBio). We sequenced and assembled BAC clones corresponding to a 1.3-Mbp complex region of chromosome 17q21.31, demonstrating 99.994% identity to Sanger assemblies of the same clones. We targeted 44 differences using Illumina sequencing and find that PacBio and Sanger assemblies share a comparable number of validated variants, albeit with different sequence context biases. Finally, we targeted a poorly assembled 766-kbp duplicated region of the chimpanzee genome and resolved the structure and organization for a fraction of the cost and time of traditional finishing approaches. Our data suggest a straightforward path for upgrading genomes to a higher quality finished state.

Genome reference and sequence variation in the large repetitive central exon of human MUC5AC.

Despite modern sequencing efforts, the difficulty in assembly of highly repetitive sequences has prevented resolution of human genome gaps, including some in the coding regions of genes with important biological functions. One such gene, MUC5AC, encodes a large, secreted mucin, which is one of the two major secreted mucins in human airways. The MUC5AC region contains a gap in the human genome reference (hg19) across the large, highly repetitive, and complex central exon. This exon is predicted to contain imperfect tandem repeat sequences and multiple conserved cysteine-rich (CysD) domains. To resolve the MUC5AC genomic gap, we used high-fidelity long PCR followed by single molecule real-time (SMRT) sequencing. This technology yielded long sequence reads and robust coverage that allowed for de novo sequence assembly spanning the entire repetitive region. Furthermore, we used SMRT sequencing of PCR amplicons covering the central exon to identify genetic variation in four individuals. The results demonstrated the presence of segmental duplications of CysD domains, insertions/deletions (indels) of tandem repeats, and single nucleotide variants. Additional studies demonstrated that one of the identified tandem repeat insertions is tagged by nonexonic single nucleotide polymorphisms. Taken together, these data illustrate the successful utility of SMRT sequencing long reads for de novo assembly of large repetitive sequences to fill the gaps in the human genome. Characterization of the MUC5AC gene and the sequence variation in the central exon will facilitate genetic and functional studies for this critical airway mucin.

An Inv(16)(p13.3q24.3)-encoded CBFA2T3-GLIS2 fusion protein defines an aggressive subtype of pediatric acute megakaryoblastic leukemia.

To define the mutation spectrum in non-Down syndrome acute megakaryoblastic leukemia (non-DS-AMKL), we performed transcriptome sequencing on diagnostic blasts from 14 pediatric patients and validated our findings in a recurrency/validation cohort consisting of 34 pediatric and 28 adult AMKL samples. Our analysis identified a cryptic chromosome 16 inversion (inv(16)(p13.3q24.3)) in 27% of pediatric cases, which encodes a CBFA2T3-GLIS2 fusion protein. Expression of CBFA2T3-GLIS2 in Drosophila and murine hematopoietic cells induced bone morphogenic protein (BMP) signaling and resulted in a marked increase in the self-renewal capacity of hematopoietic progenitors. These data suggest that expression of CBFA2T3-GLIS2 directly contributes to leukemogenesis.

Polymorphism discovery in high-throughput resequenced microarray-enriched human genomic loci.

Identifying genetic variants and mutations that underlie human diseases requires development of robust, cost-effective tools for routine resequencing of regions of interest in the human genome. Here, we demonstrate that coupling Applied Biosystems SOLiD system-sequencing platform with microarray capture of targeted regions provides an efficient and robust method for high-coverage resequencing and polymorphism discovery in human protein-coding exons.

ALLPATHS 2: small genomes assembled accurately and with high continuity from short paired reads.

We demonstrate that genome sequences approaching finished quality can be generated from short paired reads. Using 36 base (fragment) and 26 base (jumping) reads from five microbial genomes of varied GC composition and sizes up to 40 Mb, ALLPATHS2 generated assemblies with long, accurate contigs and scaffolds. Velvet and EULER-SR were less accurate. For example, for Escherichia coli, the fraction of 10-kb stretches that were perfect was 99.8% (ALLPATHS2), 68.7% (Velvet), and 42.1% (EULER-SR).

Sequence and structural variation in a human genome uncovered by short-read, massively parallel ligation sequencing using two-base encoding.

We describe the genome sequencing of an anonymous individual of African origin using a novel ligation-based sequencing assay that enables a unique form of error correction that improves the raw accuracy of the aligned reads to >99.9%, allowing us to accurately call SNPs with as few as two reads per allele. We collected several billion mate-paired reads yielding approximately 18x haploid coverage of aligned sequence and close to 300x clone coverage. Over 98% of the reference genome is covered with at least one uniquely placed read, and 99.65% is spanned by at least one uniquely placed mate-paired clone. We identify over 3.8 million SNPs, 19% of which are novel. Mate-paired data are used to physically resolve haplotype phases of nearly two-thirds of the genotypes obtained and produce phased segments of up to 215 kb. We detect 226,529 intra-read indels, 5590 indels between mate-paired reads, 91 inversions, and four gene fusions. We use a novel approach for detecting indels between mate-paired reads that are smaller than the standard deviation of the insert size of the library and discover deletions in common with those detected with our intra-read approach. Dozens of mutations previously described in OMIM and hundreds of nonsynonymous single-nucleotide and structural variants in genes previously implicated in disease are identified in this individual. There is more genetic variation in the human genome still to be uncovered, and we provide guidance for future surveys in populations and cancer biopsies.

Preparation of genome-wide DNA fragment libraries using bisulfite in polyacrylamide gel electrophoresis slices with formamide denaturation and quality control for massively parallel sequencing by oligonucleotide ligation and detection.

Bisulfite sequencing is widely used for analysis of DNA methylation status (i.e., 5-methylcytosine [5mC] vs. cytosine [C]) in CpG-rich or other loci in genomic DNA (gDNA). Such methods typically involve reaction of gDNA with bisulfite followed by polymerase chain reaction (PCR) amplification of specific regions of interest that, overall, converts C-->T (thymine) and 5mC-->C and then capillary sequencing to measure C versus T composition at CpG sites. Massively parallel sequencing by oligonucleotide ligation and detection (SOLiD) has recently enabled relatively low-cost whole genome sequencing, and it would be highly desirable to apply such massively parallel sequencing to bisulfite-converted whole genomes to determine DNA methylation status of an entire genome, which has heretofore not been reported. As an initial step toward achieving this goal, we have extended our ongoing interest in improving bisulfite conversion sample preparation to include a human genome-wide fragment library for SOliD. The current article features novel use of formamide denaturant during bisulfite conversion of a suitably constructed library directly in a band slice from polyacryamide gel electrophoresis (PAGE). To validate this new protocol for 5mC-protected fragment library conversion, which we refer to as Bis-PAGE, capillary-based size analysis and Sanger sequencing were carried out for individual amplicons derived from single-molecule PCR (smPCR) of randomly selected library fragments. smPCR/Capillary Sanger sequencing of approximately 200 amplicons unambiguously demonstrated greater than 99% C-->T conversion. All of these approximately 200 Sanger sequences were analyzed with a previously published web-accessible bioinformatics tool (methBLAST) for mapping to human chromosomes, the results of which indicated random distribution of analyzed fragments across all chromosomes. Although these particular Bis-PAGE conversion and quality control methods were exemplified in the context of a fragment library for SOLiD, the concepts can be generalized to include other genome-wide library constructions intended for DNA methylation analysis by alternative high-throughput or massively parallelized methods that are currently available.

Rapid whole-genome mutational profiling using next-generation sequencing technologies.

Forward genetic mutational studies, adaptive evolution, and phenotypic screening are powerful tools for creating new variant organisms with desirable traits. However, mutations generated in the process cannot be easily identified with traditional genetic tools. We show that new high-throughput, massively parallel sequencing technologies can completely and accurately characterize a mutant genome relative to a previously sequenced parental (reference) strain. We studied a mutant strain of Pichia stipitis, a yeast capable of converting xylose to ethanol. This unusually efficient mutant strain was developed through repeated rounds of chemical mutagenesis, strain selection, transformation, and genetic manipulation over a period of seven years. We resequenced this strain on three different sequencing platforms. Surprisingly, we found fewer than a dozen mutations in open reading frames. All three sequencing technologies were able to identify each single nucleotide mutation given at least 10-15-fold nominal sequence coverage. Our results show that detecting mutations in evolved and engineered organisms is rapid and cost-effective at the whole-genome level using new sequencing technologies. Identification of specific mutations in strains with altered phenotypes will add insight into specific gene functions and guide further metabolic engineering efforts.

Stem cell transcriptome profiling via massive-scale mRNA sequencing.

We developed a massive-scale RNA sequencing protocol, short quantitative random RNA libraries or SQRL, to survey the complexity, dynamics and sequence content of transcriptomes in a near-complete fashion. This method generates directional, random-primed, linear cDNA libraries that are optimized for next-generation short-tag sequencing. We surveyed the poly(A)(+) transcriptomes of undifferentiated mouse embryonic stem cells (ESCs) and embryoid bodies (EBs) at an unprecedented depth (10 Gb), using the Applied Biosystems SOLiD technology. These libraries capture the genomic landscape of expression, state-specific expression, single-nucleotide polymorphisms (SNPs), the transcriptional activity of repeat elements, and both known and new alternative splicing events. We investigated the impact of transcriptional complexity on current models of key signaling pathways controlling ESC pluripotency and differentiation, highlighting how SQRL can be used to characterize transcriptome content and dynamics in a quantitative and reproducible manner, and suggesting that our understanding of transcriptional complexity is far from complete.

A high-resolution, nucleosome position map of C. elegans reveals a lack of universal sequence-dictated positioning.

Using the massively parallel technique of sequencing by oligonucleotide ligation and detection (SOLiD; Applied Biosystems), we have assessed the in vivo positions of more than 44 million putative nucleosome cores in the multicellular genetic model organism Caenorhabditis elegans. These analyses provide a global view of the chromatin architecture of a multicellular animal at extremely high density and resolution. While we observe some degree of reproducible positioning throughout the genome in our mixed stage population of animals, we note that the major chromatin feature in the worm is a diversity of allowed nucleosome positions at the vast majority of individual loci. While absolute positioning of nucleosomes can vary substantially, relative positioning of nucleosomes (in a repeated array structure likely to be maintained at least in part by steric constraints) appears to be a significant property of chromatin structure. The high density of nucleosomal reads enabled a substantial extension of previous analysis describing the usage of individual oligonucleotide sequences along the span of the nucleosome core and linker. We release this data set, via the UCSC Genome Browser, as a resource for the high-resolution analysis of chromatin conformation and DNA accessibility at individual loci within the C. elegans genome.

Analysis of the Otd-dependent transcriptome supports the evolutionary conservation of CRX/OTX/OTD functions in flies and vertebrates.

Homeobox transcription factors of the vertebrate CRX/OTX family play critical roles in photoreceptor neurons, the rostral brain and circadian processes. In mouse, the three related proteins, CRX, OTX1, and OTX2, fulfill these functions. In Drosophila, the single founding member of this gene family, called orthodenticle (otd), is required during embryonic brain and photoreceptor neuron development. We have used global gene expression analysis in late pupal heads to better characterize the post-embryonic functions of Otd in Drosophila. We have identified 61 genes that are differentially expressed between wild type and a viable eye-specific otd mutant allele. Among them, about one-third represent potentially direct targets of Otd based on their association with evolutionarily conserved Otd-binding sequences. The spectrum of biological functions associated with these gene targets establishes Otd as a critical regulator of photoreceptor morphology and phototransduction, as well as suggests its involvement in circadian processes. Together with the well-documented role of otd in embryonic patterning, this evidence shows that vertebrate and fly genes contribute to analogous biological processes, notwithstanding the significant divergence of the underlying genetic pathways. Our findings underscore the common evolutionary history of photoperception-based functions in vertebrates and invertebrates and support the view that a complex nervous system was already present in the last common ancestor of all bilateria.

Direct control of neurogenesis by selector factors in the fly eye: regulation of atonal by Ey and So.

During eye development, the selector factors of the Eyeless/Pax6 or Retinal Determination (RD) network control specification of organ-type whereas the bHLH-type proneural factor Atonal drives neurogenesis. Although significant progress has been made in dissecting the acquisition of ;eye identity' at the transcriptional level, the molecular mechanisms underlying the progression from neuronal progenitor to differentiating neuron remain unclear. A recently proposed model for the integration of organ specification and neurogenesis hypothesizes that atonal expression in the eye is RD-network-independent and that Eyeless works in parallel or downstream of atonal to modify the neurogenetic program. We show here that distinct cis-regulatory elements control atonal expression specifically in the eye and that the RD factors Eyeless and Sine oculis function as direct regulators. We find that these transcription factors interact in vitro and provide indirect evidence that this interaction may be required in vivo. The subordination of neurogenesis to the RD pathway in the eye provides a direct mechanism for the coordination of neurogenesis and tissue specification during sensory organ formation.

Partner specificity is essential for proper function of the SIX-type homeodomain proteins Sine oculis and Optix during fly eye development.

The development of the Drosophila visual system utilizes two members of the highly conserved Six-Homeobox family of transcription factor, Sine oculis and Optix. Although in vitro studies have detected differences in DNA-binding and interactions with some co-factors, questions remain as to what extent the activity for these two transcriptional regulators is redundant or specific in vivo. In this work, we show that the SoD mutation within the Six domain does not abolish DNA-protein interactions, but alters co-factor binding specificity to resemble that of Optix. A mutation in the same region of Optix alters its activity in vivo. We propose that the dominant mutant phenotype is primarily due to an alteration in binding properties of the Sine oculis protein and that distinct partner interactions is one important mechanism in determining significant functional differences between these highly conserved factors during eye development.

von Willebrand's disease and psychotic disorders: co-segregation and genetic associations.

OBJECTIVES: To evaluate co-segregation and genetic associations between von Willebrand's disease (vWD) and psychotic disorders. METHODS: The study was initiated following ascertainment of a nuclear family in which four members were diagnosed with vWD and psychotic/mood disorders. As co-segregation was uncertain in the extended pedigree, we also investigated population-based linkage and association using polymorphisms of vWF, the gene conferring susceptibility to vWD. Three common vWF polymorphisms were investigated among 194 patients with psychotic disorders (bipolar I disorder, BD I; schizoaffective disorder, SZA and schizophrenia, SZ) and their parents. The cases were also compared with unrelated population-based controls (n = 183). RESULTS: The transmission disequilibrium test and related analyses suggested nominally significant transmission distortion of one allele and related haplotypes to the probands from their parents. The most significant results were obtained among patients with BD I, and similar trends were also evident in the SZ sample. Comparisons between the cases and population-based controls did not reveal associations, though marginally significant case-control differences were obtained in the BD I sample. CONCLUSIONS: These studies are consistent with association and linkage between vWF and BD I. However, given the relatively small sample, stochastic variation is an alternative explanation.

Coordinating proliferation and tissue specification to promote regional identity in the Drosophila head.

The Decapentaplegic and Notch signaling pathways are thought to direct regional specification in the Drosophila eye-antennal epithelium by controlling the expression of selector genes for the eye (Eyeless/Pax6, Eyes absent) and/or antenna (Distal-less). Here, we investigate the function of these signaling pathways in this process. We find that organ primordia formation is indeed controlled at the level of Decapentaplegic expression but critical steps in regional specification occur earlier than previously proposed. Contrary to previous findings, Notch does not specify eye field identity by promoting Eyeless expression but it influences eye primordium formation through its control of proliferation. Our analysis of Notch function reveals an important connection between proliferation, field size, and regional specification. We propose that field size modulates the interaction between the Decapentaplegic and Wingless pathways, thereby linking proliferation and patterning in eye primordium development.

Linkage and association between serotonin 2A receptor gene polymorphisms and bipolar I disorder.

Several inconsistent associations between bipolar I disorder (BD1) and polymorphisms of the genes encoding the serotonin 2A receptor (HTR2A) have been published. We conducted the Transmission Disequilibrium Test (TDT) and case-control comparisons involving nine single nucleotide polymorphisms at the serotonin 2A receptor gene (four SNPs of HTR2A exons and five flanking SNPs). Comparison of BD1 cases (n = 93) with a group of unrelated population based controls (n = 92) revealed associations with SNPs on exons 2 and 3 (516C/T and 1354C/T, respectively), consistent with haplotype-based differences. Analysis of the cases and their available parents using the TDT suggested significant linkage and associations with 1354C/T, as well as haplotypes bearing this SNP. Our results support an etiological role for HTR2A in BD1. In view of the relatively small sample, replicate studies using large samples are needed.