Mapping and sequencing of structural variation from eight human genomes.

Genetic variation among individual humans occurs on many different scales, ranging from gross alterations in the human karyotype to single nucleotide changes. Here we explore variation on an intermediate scale--particularly insertions, deletions and inversions affecting from a few thousand to a few million base pairs. We employed a clone-based method to interrogate this intermediate structural variation in eight individuals of diverse geographic ancestry. Our analysis provides a comprehensive overview of the normal pattern of structural variation present in these genomes, refining the location of 1,695 structural variants. We find that 50% were seen in more than one individual and that nearly half lay outside regions of the genome previously described as structurally variant. We discover 525 new insertion sequences that are not present in the human reference genome and show that many of these are variable in copy number between individuals. Complete sequencing of 261 structural variants reveals considerable locus complexity and provides insights into the different mutational processes that have shaped the human genome. These data provide the first high-resolution sequence map of human structural variation--a standard for genotyping platforms and a prelude to future individual genome sequencing projects.

Power to detect risk alleles using genome-wide tag SNP panels.

Advances in high-throughput genotyping and the International HapMap Project have enabled association studies at the whole-genome level. We have constructed whole-genome genotyping panels of over 550,000 (HumanHap550) and 650,000 (HumanHap650Y) SNP loci by choosing tag SNPs from all populations genotyped by the International HapMap Project. These panels also contain additional SNP content in regions that have historically been overrepresented in diseases, such as nonsynonymous sites, the MHC region, copy number variant regions and mitochondrial DNA. We estimate that the tag SNP loci in these panels cover the majority of all common variation in the genome as measured by coverage of both all common HapMap SNPs and an independent set of SNPs derived from complete resequencing of genes obtained from SeattleSNPs. We also estimate that, given a sample size of 1,000 cases and 1,000 controls, these panels have the power to detect single disease loci of moderate risk (lambda approximately 1.8-2.0). Relative risks as low as lambda approximately 1.1-1.3 can be detected using 10,000 cases and 10,000 controls depending on the sample population and disease model. If multiple loci are involved, the power increases significantly to detect at least one locus such that relative risks 20%-35% lower can be detected with 80% power if between two and four independent loci are involved. Although our SNP selection was based on HapMap data, which is a subset of all common SNPs, these panels effectively capture the majority of all common variation and provide high power to detect risk alleles that are not represented in the HapMap data.

Design of tag SNP whole genome genotyping arrays.

Whole genome association studies have recently been enabled by combining tag SNP information derived from the International HapMap project with novel whole genome genotyping array technologies. In particular, Infinium whole genome genotyping (WGG) technology now has the power to genotype over 1 million SNPs on a single array. Additionally, this assay provides access to virtually any SNP in the genome enabling selection of optimized SNP content . In this chapter, we provide an overview of the tag SNP-based selection strategy for Infinium whole-genome genotyping BeadChips, including the Human 1 M BeadChip. These advances in both SNP content and technology have enabled both large-scale whole-genome disease association (WGAS) and copy number variation (CNV) studies with the ultimate goal of identifying common genetic variants, disease-associated loci, proteins, and biomarkers.

Genomic profiling of 766 cancer-related genes in archived esophageal normal and carcinoma tissues.

We employed the BeadArraytrade mark technology to perform a genetic analysis in 33 formalin-fixed, paraffin-embedded (FFPE) human esophageal carcinomas, mostly squamous-cell-carcinoma (ESCC), and their adjacent normal tissues. A total of 1,432 single nucleotide polymorphisms (SNPs) derived from 766 cancer-related genes were genotyped with partially degraded genomic DNAs isolated from these samples. This directly targeted genomic profiling identified not only previously reported somatic gene amplifications (e.g., CCND1) and deletions (e.g., CDKN2A and CDKN2B) but also novel genomic aberrations. Among these novel targets, the most frequently deleted genomic regions were chromosome 3p (including tumor suppressor genes FANCD2 and CTNNB1) and chromosome 5 (including tumor suppressor gene APC). The most frequently amplified genomic region was chromosome 3q (containing DVL3, MLF1, ABCC5, BCL6, AGTR1 and known oncogenes TNK2, TNFSF10, FGF12). The chromosome 3p deletion and 3q amplification occurred coincidently in nearly all of the affected cases, suggesting a molecular mechanism for the generation of somatic chromosomal aberrations. We also detected significant differences in germline allele frequency between the esophageal cohort of our study and normal control samples from the International HapMap Project for 10 genes (CSF1, KIAA1804, IL2, PMS2, IRF7, FLT3, NTRK2, MAP3K9, ERBB2 and PRKAR1A), suggesting that they might play roles in esophageal cancer susceptibility and/or development. Taken together, our results demonstrated the utility of the BeadArray technology for high-throughput genetic analysis in FFPE tumor tissues and provided a detailed genetic profiling of cancer-related genes in human esophageal cancer.

Delineation of the proximal 3q microdeletion syndrome.

Interstitial deletions of the proximal long arm of chromosome 3 are very rare and a defined clinical phenotype is not established yet. We report on the clinical, cytogenetic and molecular findings of a 20-month-old Hispanic male with a 2.5 Mb de novo deletion on q13.11q13.12. Up to now, this is the smallest deletion reported among patients with the proximal 3q microdeletion syndrome. The patient has distinct facial features including brachycephaly, broad and prominent forehead, flat nasal bridge, prominent ears, anteverted nose, tetralogy of Fallot, bilateral cryptorchidism, and peripheral skeletal abnormalities. To further delineate the proximal 3q deletion syndrome, the phenotype of our patient was compared with 10 other patients previously described. We found that ALCAM and CBLB are the only genes deleted in our patient and based on previously published data, we propose that the CBLB gene is responsible for the craniofacial phenotype in patients with deletions of proximal 3q region.

Mosaic tetrasomy 12p with triplication of 12p detected by array-based comparative genomic hybridization of peripheral blood DNA.

A patient whose dysmorphism at birth was not diagnostic for Pallister-Killian syndrome (PKS) was found to have mosaic tetrasomy 12p by an array-based comparative genomic hybridization of peripheral blood DNA. He was determined to be mosaic for 46,XY,trp(12)(p11.2 --> p13) in cultured skin fibroblasts. His appearance was typical for PKS at 4 months of age.

Microarray optimizations: increasing spot accuracy and automated identification of true microarray signals.

In this paper, fluorescent microarray images and various analysis techniques are described to improve the microarray data acquisition processes. Signal intensities produced by rarely expressed genes are initially correctly detected, but they are often lost in corrections for background, log or ratio. Our analyses indicate that a simple correlation between the mean and median signal intensities may be the best way to eliminate inaccurate microarray signals. Unlike traditional quality control methods, the low intensity signals are retained and inaccurate signals are eliminated in this mean and median correlation. With larger amounts of microarray data being generated, it becomes increasingly more difficult to analyze data on a visual basis. Our method allows for the automatic quantitative determination of accurate and reliable signals, which can then be used for normalization. We found that a mean to median correlation of 85% or higher not only retains more data than current methods, but the retained data is more accurate than traditional thresholds or common spot flagging algorithms. We have also found that by using pin microtapping and microvibrations, we can control spot quality independent from initial PCR volume.

Deletion at fragile sites is a common and early event in Barrett's esophagus.

Barrett's esophagus (BE) is a premalignant intermediate to esophageal adenocarcinoma, which develops in the context of chronic inflammation and exposure to bile and acid. We asked whether there might be common genomic alterations that could be identified as potential clinical biomarker(s) for BE by whole genome profiling. We detected copy number alterations and/or loss of heterozygosity at 56 fragile sites in 20 patients with premalignant BE. Chromosomal fragile sites are particularly sensitive to DNA breaks and are frequent sites of rearrangement or loss in many human cancers. Seventy-eight percent of all genomic alterations detected by array-CGH were associated with fragile sites. Copy number losses in early BE were observed at particularly high frequency at FRA3B (81%), FRA9A/C (71.4%), FRA5E (52.4%), and FRA 4D (52.4%), and at lower frequencies in other fragile sites, including FRA1K (42.9%), FRAXC (42.9%), FRA 12B (33.3%), and FRA16D (33.3%). Due to the consistency of the region of copy number loss, we were able to verify these results by quantitative PCR, which detected the loss of FRA3B and FRA16D, in 83% and 40% of early molecular stage BE patients, respectively. Loss of heterozygosity in these cases was confirmed through pyrosequencing at FRA3B and FRA16D (75% and 70%, respectively). Deletion and genomic instability at FRA3B and other fragile sites could thus be a biomarker of genetic damage in BE patients and a potential biomarker of cancer risk.

Increased LIS1 expression affects human and mouse brain development.

Deletions of the PAFAH1B1 gene (encoding LIS1) in 17p13.3 result in isolated lissencephaly sequence, and extended deletions including the YWHAE gene (encoding 14-3-3epsilon) cause Miller-Dieker syndrome. We identified seven unrelated individuals with submicroscopic duplication in 17p13.3 involving the PAFAH1B1 and/or YWHAE genes, and using a 'reverse genomics' approach, characterized the clinical consequences of these duplications. Increased PAFAH1B1 dosage causes mild brain structural abnormalities, moderate to severe developmental delay and failure to thrive. Duplication of YWHAE and surrounding genes increases the risk for macrosomia, mild developmental delay and pervasive developmental disorder, and results in shared facial dysmorphologies. Transgenic mice conditionally overexpressing LIS1 in the developing brain showed a decrease in brain size, an increase in apoptotic cells and a distorted cellular organization in the ventricular zone, including reduced cellular polarity but preserved cortical cell layer identity. Collectively, our results show that an increase in LIS1 expression in the developing brain results in brain abnormalities in mice and humans.

Single nucleotide polymorphism-based genome-wide chromosome copy change, loss of heterozygosity, and aneuploidy in Barrett's esophagus neoplastic progression.

Chromosome copy gain, loss, and loss of heterozygosity (LOH) involving most chromosomes have been reported in many cancers; however, less is known about chromosome instability in premalignant conditions. 17p LOH and DNA content abnormalities have been previously reported to predict progression from Barrett's esophagus (BE) to esophageal adenocarcinoma (EA). Here, we evaluated genome-wide chromosomal instability in multiple stages of BE and EA in whole biopsies. Forty-two patients were selected to represent different stages of progression from BE to EA. Whole BE or EA biopsies were minced, and aliquots were processed for flow cytometry and genotyped with a paired constitutive control for each patient using 33,423 single nucleotide polymorphisms (SNP). Copy gains, losses, and LOH increased in frequency and size between early- and late-stage BE (P < 0.001), with SNP abnormalities increasing from <2% to >30% in early and late stages, respectively. A set of statistically significant events was unique to either early or late, or both, stages, including previously reported and novel abnormalities. The total number of SNP alterations was highly correlated with DNA content aneuploidy and was sensitive and specific to identify patients with concurrent EA (empirical receiver operating characteristic area under the curve = 0.91). With the exception of 9p LOH, most copy gains, losses, and LOH detected in early stages of BE were smaller than those detected in later stages, and few chromosomal events were common in all stages of progression. Measures of chromosomal instability can be quantified in whole biopsies using SNP-based genotyping and have potential to be an integrated platform for cancer risk stratification in BE.

Speech delay and autism spectrum behaviors are frequently associated with duplication of the 7q11.23 Williams-Beuren syndrome region.

PURPOSE: Williams-Beuren syndrome is among the most well-characterized microdeletion syndromes, caused by recurrent de novo microdeletions at 7q11.23 mediated by nonallelic homologous recombination between low copy repeats flanking this critical region. However, the clinical phenotype associated with reciprocal microduplication of this genomic region is less well described. We investigated the molecular, clinical, neurodevelopmental, and behavioral features of seven patients with dup(7)(q11.23), including two children who inherited the microduplication from one of their parents, to more fully characterize this emerging microduplication syndrome. METHODS: Patients were identified by array-based comparative genomic hybridization. Clinical examinations were performed on seven affected probands, and detailed cognitive and behavioral evaluations were carried out on four of the affected probands. RESULTS: Our findings confirm initial reports of speech delay seen in patients with dup(7)(q11.23) and further delineate and expand the phenotypic spectrum of this condition to include communication, social interactions, and repetitive interests that are often observed in individuals diagnosed with autism spectrum disorders. CONCLUSIONS: Array-based comparative genomic hybridization is a powerful means of detecting genomic imbalances and identifying molecular etiologies in the clinic setting, including genomic disorders such as Williams-Beuren syndrome and dup(7)(q11.23). We propose that dup(7)(q11.23) syndrome may be as frequent as Williams-Beuren syndrome and a previously unrecognized cause of language delay and behavioral abnormalities. Indeed, these individuals may first be referred for evaluation of autism, even if they do not ultimately meet diagnostic criteria for an autism spectrum disorder.

Schnurri transcription factors from Drosophila and vertebrates can mediate Bmp signaling through a phylogenetically conserved mechanism.

Bone Morphogenetic Proteins (Bmps) are secreted growth factors that play crucial roles in animal development across the phylogenetic spectrum. Bmp signaling results in the phosphorylation and nuclear translocation of Smads, downstream signal transducers that bind DNA. In Drosophila, the zinc finger protein Schnurri (Shn) plays a key role in signaling by the Bmp2/Bmp4 homolog Decapentaplegic (Dpp), by forming a Shn/Smad repression complex on defined promoter elements in the brinker (brk) gene. Brk is a transcriptional repressor that downregulates Dpp target genes. Thus, brk inhibition by Shn results in the upregulation of Dpp-responsive genes. We present evidence that vertebrate Shn homologs can also mediate Bmp responsiveness through a mechanism similar to Drosophila Shn. We find that a Bmp response element (BRE) from the Xenopus Vent2 promoter drives Dpp-dependent expression in Drosophila. However, in sharp contrast to its activating role in vertebrates, the frog BRE mediates repression in Drosophila. Remarkably, despite these opposite transcriptional polarities, sequence changes that abolish cis-element activity in Drosophila also affect BRE function in Xenopus. These similar cis requirements reflect conservation of trans-acting factors, as human Shn1 (hShn1; HIVEP1) can interact with Smad1/Smad4 and assemble an hShn1/Smad complex on the BRE. Furthermore, both Shn and hShn1 activate the BRE in Xenopus embryos, and both repress brk and rescue embryonic patterning defects in shn mutants. Our results suggest that vertebrate Shn proteins function in Bmp signal transduction, and that Shn proteins recruit coactivators and co-repressors in a context-dependent manner, rather than acting as dedicated activators or repressors.

High-resolution genomic profiling of chromosomal aberrations using Infinium whole-genome genotyping.

Array-CGH is a powerful tool for the detection of chromosomal aberrations. The introduction of high-density SNP genotyping technology to genomic profiling, termed SNP-CGH, represents a further advance, since simultaneous measurement of both signal intensity variations and changes in allelic composition makes it possible to detect both copy number changes and copy-neutral loss-of-heterozygosity (LOH) events. We demonstrate the utility of SNP-CGH with two Infinium whole-genome genotyping BeadChips, assaying 109,000 and 317,000 SNP loci, to detect chromosomal aberrations in samples bearing constitutional aberrations as well tumor samples at sub-100 kb effective resolution. Detected aberrations include homozygous deletions, hemizygous deletions, copy-neutral LOH, duplications, and amplifications. The statistical ability to detect common aberrations was modeled by analysis of an X chromosome titration model system, and sensitivity was modeled by titration of gDNA from a tumor cell with that of its paired normal cell line. Analysis was facilitated by using a genome browser that plots log ratios of normalized intensities and allelic ratios along the chromosomes. We developed two modes of SNP-CGH analysis, a single sample and a paired sample mode. The single sample mode computes log intensity ratios and allelic ratios by referencing to canonical genotype clusters generated from approximately 120 reference samples, whereas the paired sample mode uses a paired normal reference sample from the same individual. Finally, the two analysis modes are compared and contrasted for their utility in analyzing different types of input gDNA: low input amounts, fragmented gDNA, and Phi29 whole-genome pre-amplified DNA.

Identification of neural genes using Xenopus DNA microarrays.

To isolate novel genes regulating neural induction, we used a DNA microarray approach. As neural induction is thought to occur by means of the inhibition of bone morphogenetic protein (BMP) signaling, BMP signaling was inhibited in ectodermal cells by overexpression of a dominant-negative receptor. RNAs were isolated from control animal cap explants and from dominant-negative BMP receptor expressing animal caps and subjected to a microarray experiment using newly generated high-density Xenopus DNA microarray chips representing over 17,000 unigenes. We have identified 77 genes that are induced in animal caps after inhibition of BMP signaling, and all of these genes were subjected to whole-mount in situ hybridization analysis. Thirty-two genes showed specific expression in neural tissues. Of the 32, 14 genes have never been linked to neural induction. Two genes that are highly induced by BMP inhibition are inhibitors of Wnt signaling, suggesting that a key step in neural induction is to produce Wnt antagonists to promote anterior neural plate development. Our current analysis also proves that a microarray approach is useful in identifying novel candidate factors involved in neural induction and patterning.

A Xenopus DNA microarray approach to identify novel direct BMP target genes involved in early embryonic development.

Bone morphogenetic proteins (BMPs), a subgroup of the transforming growth factor-beta (TGF-beta) superfamily, were originally isolated from bone on the basis of their ability to induce ectopic bone development. Although BMPs are involved in a wide range of developmental and physiological functions, very few vertebrate target genes in this pathway have been identified. To identify target genes regulated by the BMP growth factor family in Xenopus, large-scale microarray analyses were conducted to discover genes directly activated by this factor in dissociated animal cap tissues treated with a combination of the protein synthesis inhibitor cycloheximide and BMP2. Consequent expression patterns and behaviors of the most highly induced genes were analyzed by in situ and reverse transcriptase-polymerase chain reaction analyses. Here, we describe two sets of the most highly induced direct BMP target genes identified using microarrays prepared from two different stages of early Xenopus development. A wide variety of genes are induced by BMP2, ranging from cell cycle controllers, enzymes, signal transduction cascade components, and components of the blood and vascular system. The finding reinforces the notion that BMP signals play important roles in a variety of biological processes.

Phylogenetic footprinting and genome scanning identify vertebrate BMP response elements and new target genes.

The complex gene regulatory networks governed by growth factor signaling are still poorly understood. In order to accelerate the rate of progress in uncovering these networks, we explored the usefulness of interspecies sequence comparison (phylogenetic footprinting) to identify conserved growth factor response elements. The promoter regions of two direct target genes of Bone Morphogenetic Protein (BMP) signaling in Xenopus, Xvent2 and XId3, were compared with the corresponding human and/or mouse counterparts to identify conserved sequences. A comparison between the Xenopus and human Vent2 promoter sequences revealed a highly conserved 21 bp sequence that overlaps the previously reported Xvent2 BMP response element (BRE). Reporter gene assays using Xenopus animal pole ectodermal explants (animal caps) revealed that this conserved 21 bp BRE is both necessary and sufficient for BMP responsiveness. We combine the same phylogenetic footprinting approach with luciferase assays to identify a highly conserved 49 bp BMP responsive region in the Xenopus Id3 promoter. GFP reporters containing multimers of either the Xvent2 or XId3 BREs appear to recapitulate endogenous BMP signaling activity in transgenic Xenopus embryos. Comparison of the Xvent2 and the XId3 BRE revealed core sequence features that are both necessary and sufficient for BMP responsiveness: a Smad binding element (SBE) and a GC-rich element resembling an OAZ binding site. Based on these findings, we have implemented genome scanning to identify over 100 additional putative target genes containing 2 or more BRE-like sequences which are conserved between human and mouse. RT-PCR and in situ analyses revealed that this in silico approach can effectively be used to identify potential BMP target genes.