High-resolution identification of chromosomal abnormalities using oligonucleotide arrays containing 116,204 SNPs

The deleterious variants of N-acetylgalactosamine-6-sulfatase (GalN6S) enzyme trigger Morquio a syndrome by disrupting protein foldings

Lysosomal enzymes degrade cellular macromolecules, while their inactivation causes human hereditary metabolic disorders. Mucopolysaccharidosis IVA (MPS IVA; Moquio A syndrome) is one of the lysosomal storage disorders caused by a defective Galactosamine-6-sulfatase (GalN6S) enzyme. In several populations, disease incidence is elevated due to missense mutations brought on by non-synonymous allelic variation in the GalN6S enzyme. Here, we studied the effect of non-synonymous single nucleotide polymorphism (nsSNPs) on the structural dynamics of the GalN6S enzyme and its binding with N-acetylgalactosamine (GalNAc) using all-atom molecular dynamics simulation and an essential dynamics approach. Consequently, in this study, we have identified three functionally disruptive mutations in domain-I and domain-II, that is, S80L, R90W, and S162F, which presumably contribute to post-translational modifications. The study delineated that both domains work cooperatively, and alteration in domain II (S80L, R90W) leads to conformational changes in the catalytic site in domain-I, while mutation S162F mainly provokes higher residual flexibility of domain II. These results show that these mutations impair the hydrophobic core, implying that Morquio A syndrome is caused by misfolding of the GalN6S enzyme. The results also show the instability of the GalN6S-GalNAc complex upon substitution. Overall, the structural dynamics resulting from point mutations give the molecular rationale for Moquio A syndrome and, more importantly, the Mucopolysaccharidoses (MPS) family of diseases, re-establishing MPS IVA as a protein-folding disease.Communicated by Ramaswamy H. Sarma.

In Silico Characterization and Analysis of Clinically Significant Variants of Lipase-H (LIPH Gene) Protein Associated with Hypotrichosis

Hypotrichosis is an uncommon type of alopecia (hair loss) characterized by coarse scalp hair caused by the reduced or fully terminated activity of the Lipase-H (LIPH) enzyme. LIPH gene mutations contribute to the development of irregular or non-functional proteins. Because several cellular processes, including cell maturation and proliferation, are inhibited when this enzyme is inactive, the hair follicles become structurally unreliable, undeveloped, and immature. This results in brittle hair, as well as altered hair shaft development and structure. Because of these nsSNPs, the protein's structure and/or function may be altered. Given the difficulty in discovering functional SNPs in genes associated with disease, it is possible to assess potential functional SNPs before conducting broader population investigations. As a result, in our in silico analysis, we separated potentially hazardous nsSNPs of the LIPH gene from benign representatives using a variety of sequencing and architecture-based bioinformatics approaches. Using seven prediction algorithms, 9 out of a total of 215 nsSNPs were shown to be the most likely to cause harm. In order to distinguish between potentially harmful and benign nsSNPs of the LIPH gene, in our in silico investigation, we employed a range of sequence- and architecture-based bioinformatics techniques. Three nsSNPs (W108R, C246S, and H248N) were chosen as potentially harmful. The present findings will likely be helpful in future large population-based studies, as well as in drug discovery, particularly in the creation of personalized medicine, since this study provides an initial thorough investigation of the functional nsSNPs of LIPH.

ADGRV1 is implicated in myoclonic epilepsy

OBJECTIVE

To investigate the significance of variation in ADGRV1 (also known as GPR98, MASS1, and VLGR1), MEF2C, and other genes at the 5q14.3 chromosomal locus in myoclonic epilepsy.

METHODS

We studied the epilepsy phenotypes of 4 individuals with 5q14.3 deletion and found that all had myoclonic seizures. We then screened 6 contiguous genes at 5q14.3, MEF2C, CETN3, MBLAC2, POLR3G, LYSMD3, and ADGRV1, in a 95-patient cohort with epilepsy and myoclonic seizures. Of these genes, point mutations in MEF2C cause a phenotype involving seizures and intellectual disability. A role for ADGRV1 in epilepsy has been proposed previously, based on a recessive mutation in the Frings mouse model of audiogenic seizures, as well as a shared homologous region with another epilepsy gene, LGI1.

RESULTS

Six patients from the myoclonic epilepsy cohort had likely pathogenic ultra-rare ADGRV1 variants, and statistical analysis showed that ultra-rare variants were significantly overrepresented when compared to healthy population data from the Genome Aggregation Database. Of the remaining genes, no definite pathogenic variants were identified.

SIGNIFICANCE

Our data suggest that the ADGRV1 variation contributes to epilepsy with myoclonic seizures, although the inheritance pattern may be complex in many cases. In patients with 5q14.3 deletion and epilepsy, ADGRV1 haploinsufficiency likely contributes to seizure development. The latter is a shift from current thinking, as MEF2C haploinsufficiency has been considered the main cause of epilepsy in 5q14.3 deletion syndrome. In cases of 5q14.3 deletion and epilepsy, seizures likely occur due to haploinsufficiency of one or both of ADGRV1 and MEF2C.

A copy number variation genotyping method for aneuploidy detection in spontaneous abortion specimens

OBJECTIVE

Chromosomal abnormalities such as aneuploidy have been shown to be responsible for causing spontaneous abortion. Genetic evaluation of abortions is currently underperformed. Screening for aneuploidy in the products of conception can help determine the etiology. We designed a high-throughput ligation-dependent probe amplification (HLPA) assay to examine aneuploidy of 24 chromosomes in miscarriage tissues and aimed to validate the performance of this technique.

METHODS

We carried out aneuploidy screening in 98 fetal tissue samples collected from female subjects with singleton pregnancies who experienced spontaneous abortion. The mean maternal age was 31.6 years (range: 24-43), and the mean gestational age was 10.2 weeks (range: 4.6-14.1). HLPA was performed in parallel with array comparative genomic hybridization, which is the gold standard for aneuploidy detection in clinical practices. The results from the two platforms were compared.

RESULTS

Forty-nine out of ninety-eight samples were found to be aneuploid. HLPA showed concordance with array comparative genomic hybridization in diagnosing aneuploidy.

CONCLUSION

High-throughput ligation-dependent probe amplification is a rapid and accurate method for aneuploidy detection. It can be used as a cost-effective screening procedure in clinical spontaneous abortions. © 2016 John Wiley & Sons, Ltd.

Combined Analysis of SNP Array Data Identifies Novel CNV Candidates and Pathways in Ependymoma and Mesothelioma

Copy number variation is a class of structural genomic modifications that includes the gain and loss of a specific genomic region, which may include an entire gene. Many studies have used low-resolution techniques to identify regions that are frequently lost or amplified in cancer. Usually, researchers choose to use proprietary or non-open-source software to detect these regions because the graphical interface tends to be easier to use. In this study, we combined two different open-source packages into an innovative strategy to identify novel copy number variations and pathways associated with cancer. We used a mesothelioma and ependymoma published datasets to assess our tool. We detected previously described and novel copy number variations that are associated with cancer chemotherapy resistance. We also identified altered pathways associated with these diseases, like cell adhesion in patients with mesothelioma and negative regulation of glutamatergic synaptic transmission in ependymoma patients. In conclusion, we present a novel strategy using open-source software to identify copy number variations and altered pathways associated with cancer.

Application of array-based comparative genomic hybridization to pediatric neurologic diseases

PURPOSE

Array comparative genomic hybridization (array-CGH) is a technique used to analyze quantitative increase or decrease of chromosomes by competitive DNA hybridization of patients and controls. This study aimed to evaluate the benefits and yield of array-CGH in comparison with conventional karyotyping in pediatric neurology patients.

MATERIALS AND METHODS

We included 87 patients from the pediatric neurology clinic with at least one of the following features: developmental delay, mental retardation, dysmorphic face, or epilepsy. DNA extracted from patients and controls was hybridized on the Roche NimbleGen 135K oligonucleotide array and compared with G-band karyotyping. The results were analyzed with findings reported in recent publications and internet databases.

RESULTS

Chromosome imbalances, including 9 cases detected also by G-band karyotyping, were found in 28 patients (32.2%), and at least 19 of them seemed to be causally related to the abnormal phenotypes. Regarding each clinical symptom, 26.2% of 42 developmental delay patients, 44.4% of 18 mental retardation patients, 42.9% of 28 dysmorphic face patients, and 34.6% of 26 epilepsy patients showed abnormal array results.

CONCLUSION

Although there were relatively small number of tests in patients with pediatric neurologic disease, this study demonstrated that array-CGH is a very useful tool for clinical diagnosis of unknown genome abnormalities performed in pediatric neurology clinics.

Hybridization and amplification rate correction for affymetrix SNP arrays

BACKGROUND

Copy number variation (CNV) is essential to understand the pathology of many complex diseases at the DNA level. Affymetrix SNP arrays, which are widely used for CNV studies, significantly depend on accurate copy number (CN) estimation. Nevertheless, CN estimation may be biased by several factors, including cross-hybridization and training sample batch, as well as genomic waves of intensities induced by sequence-dependent hybridization rate and amplification efficiency. Since many available algorithms only address one or two of the three factors, a high false discovery rate (FDR) often results when identifying CNV. Therefore, we have developed a new CNV detection pipeline which is based on hybridization and amplification rate correction (CNVhac).

METHODS

CNVhac first estimates the allelic concentrations (ACs) of target sequences by using the sample independent parameters trained through physicochemical hybridization law. Then the raw CN is estimated by taking the ratio of AC to the corresponding average AC from a reference sample set for one specific site. Finally, a hidden Markov model (HMM) segmentation process is implemented to detect CNV regions.

RESULTS

Based on public HapMap data, the results show that CNVhac effectively smoothes the genomic waves and facilitates more accurate raw CN estimates compared to other methods. Moreover, CNVhac alleviates, to a certain extent, the sample dependence of inference and makes CNV calling with appreciable low FDRs.

CONCLUSION

CNVhac is an effective approach to address the common difficulties in SNP array analysis, and the working principles of CNVhac can be easily extended to other platforms.

Amplification ratio control system for copy number variation genotyping

We describe a generic design for ratiometric analysis suitable for determination of copy number variation (CNV) class of a gene. Following two initial sequence-specific PCR priming cycles, both ends of both amplicons (one test and one reference) in a duplex reaction, are all primed by the same universal primer (UP). Following each amplification denaturation step, the UP target and its reverse complement (UP′) in each strand form a hairpin. The bases immediately beyond the 3′-end of the UP and 5′ of UP′ are chosen such as not to base pair in the hairpin (otherwise priming is ablated). This hairpin creates a single constant environment for priming events and chaperones free 3′-ends of amplicon strands. The resultant ‘amplification ratio control system’ (ARCS) permits ratiometric representation of amplicons relative to the original template into PCR plateau phase. These advantages circumvent the need for real-time PCR for quantitation. Choice of different %(G+C) content for the target and reference amplicons allows liquid phase thermal melt discrimination and quantitation of amplicons. The design is generic, simple to set up and economical. Comparisons with real-time PCR and other techniques are made and CNV assays demonstrated for haptoglobin duplicon and ‘chemokine (C-C motif) ligand 3-like 1’ gene.

Genome rearrangements detected by SNP microarrays in individuals with intellectual disability referred with possible Williams syndrome

Background

Intellectual disability (ID) affects 2–3% of the population and may occur with or without multiple congenital anomalies (MCA) or other medical conditions. Established genetic syndromes and visible chromosome abnormalities account for a substantial percentage of ID diagnoses, although for ∼50% the molecular etiology is unknown. Individuals with features suggestive of various syndromes but lacking their associated genetic anomalies pose a formidable clinical challenge. With the advent of microarray techniques, submicroscopic genome alterations not associated with known syndromes are emerging as a significant cause of ID and MCA.

Methodology/Principal Findings

High-density SNP microarrays were used to determine genome wide copy number in 42 individuals: 7 with confirmed alterations in the WS region but atypical clinical phenotypes, 31 with ID and/or MCA, and 4 controls. One individual from the first group had the most telomeric gene in the WS critical region deleted along with 2 Mb of flanking sequence. A second person had the classic WS deletion and a rearrangement on chromosome 5p within the Cri du Chat syndrome (OMIM:123450) region. Six individuals from the ID/MCA group had large rearrangements (3 deletions, 3 duplications), one of whom had a large inversion associated with a deletion that was not detected by the SNP arrays.

Conclusions/Significance

Combining SNP microarray analyses and qPCR allowed us to clone and sequence 21 deletion breakpoints in individuals with atypical deletions in the WS region and/or ID or MCA. Comparison of these breakpoints to databases of genomic variation revealed that 52% occurred in regions harboring structural variants in the general population. For two probands the genomic alterations were flanked by segmental duplications, which frequently mediate recurrent genome rearrangements; these may represent new genomic disorders. While SNP arrays and related technologies can identify potentially pathogenic deletions and duplications, obtaining sequence information from the breakpoints frequently provides additional information.

A maternally inherited chromosome 18q22.1 deletion in a male with late-presenting diaphragmatic hernia and microphthalmia-evaluation of DSEL as a candidate gene for the diaphragmatic defect

Using an Affymetrix GeneChip(R) Human Mapping 100K Set to study a patient with a late-presenting, right-sided diaphragmatic hernia and microphthalmia, we found a maternally inherited deletion that was 2.7 Mb in size at chromosome 18q22.1. Mapping of this deletion using fluorescence in situ hybridization revealed three deleted genes-CDH19, DSEL, and TXNDC10, and one gene that contained the deletion breakpoint, CCDC102B. We selected DSEL for further study in 125 patients with diaphragmatic hernias, as it is involved in the synthesis of decorin, a protein that is required for normal collagen formation and that is upregulated during myogenesis. We found p.Met14Ile in an unrelated patient with a late-presenting, anterior diaphragmatic hernia. In the murine diaphragm, Dsel was only weakly expressed at the time of diaphragm closure and its expression in C2C12 myoblast cells did not change significantly during myoblast differentiation, thus reducing the likelihood that the gene is involved in myogenesis of the diaphragm. Although it is possible that the 18q22.1 deletion and haploinsufficiency for DSEL contributed to the diaphragmatic defect in the patient, a definite role for DSEL and decorin in the formation of the collagen-containing, central tendon of the diaphragm has not yet been established.

Advances in genome studies: The PAG 2010 conference

An overview is provided of the advances in plant, animal and human genome studies by summarizing the contents of seven plenary lectures presented at the Plant and Animal Genome (PAG) meeting in January 2010. The area of biology covered was wide and reflected the nature of this fast moving science.

Statistics and bioinformatics in nutritional sciences: analysis of complex data in the era of systems biology

Over the past 2 decades, there have been revolutionary developments in life science technologies characterized by high throughput, high efficiency, and rapid computation. Nutritionists now have the advanced methodologies for the analysis of DNA, RNA, protein, low-molecular-weight metabolites, as well as access to bioinformatics databases. Statistics, which can be defined as the process of making scientific inferences from data that contain variability, has historically played an integral role in advancing nutritional sciences. Currently, in the era of systems biology, statistics has become an increasingly important tool to quantitatively analyze information about biological macromolecules. This article describes general terms used in statistical analysis of large, complex experimental data. These terms include experimental design, power analysis, sample size calculation, and experimental errors (Type I and II errors) for nutritional studies at population, tissue, cellular, and molecular levels. In addition, we highlighted various sources of experimental variations in studies involving microarray gene expression, real-time polymerase chain reaction, proteomics, and other bioinformatics technologies. Moreover, we provided guidelines for nutritionists and other biomedical scientists to plan and conduct studies and to analyze the complex data. Appropriate statistical analyses are expected to make an important contribution to solving major nutrition-associated problems in humans and animals (including obesity, diabetes, cardiovascular disease, cancer, ageing, and intrauterine growth retardation).

A comprehensive profile of DNA copy number variations in a Korean population: identification of copy number invariant regions among Koreans

To examine copy number variations among the Korean population, we compared individual genomes with the Korean reference genome assembly using the publicly available Korean HapMap SNP 50 k chip data from 90 individuals. Korean individuals exhibited 123 copy number variation regions (CNVRs) covering 27.2 mb, equivalent to 1.0% of the genome in the copy number variation (CNV) analysis using the combined criteria of P value (P<0.01) and standard deviation of copy numbers (SD>or= 0.25) among study subjects. In contrast, when compared to the Affymetrix reference genome assembly from multiple ethnic groups, considerably more CNVRs (n=643) were detected in larger proportions (5.0%) of the genome covering 135.1 mb even by more stringent criteria (P<0.001 and SD>or=0.25), reflecting ethnic diversity of structural variations between Korean and other populations. Some CNVRs were validated by the quantitative multiplex PCR of short fluorescent fragment (QMPSF) method, and then copy number invariant regions were detected among the study subjects. These copy number invariant regions would be used as good internal controls for further CNV studies. Lastly, we demonstrated that the CNV information could stratify even a single ethnic population with a proper reference genome assembly from multiple heterogeneous populations.

Genomic landscape of meningiomas

Meningiomas are one of the most common adult brain tumors. For most patients, surgical excision is curative. However, up to 20% recur. Currently, the molecular determinants predicting recurrence and malignant transformation are lacking. We performed retrospective global genetic and genomic analysis of 85 meningioma samples of various grades. Copy number alterations were assessed by 100K single-nucleotide polymorphism arrays and correlated with gene expression, proliferation indices and clinical outcome. In addition to chromosome 22q loss, which was detected in the majority of clinical samples, chromosome 6q and 14q loss was significantly more common in recurrent tumors and was associated with anaplastic histology. Five "classes" of meningiomas were detected by gene expression analysis that correlated with copy number alterations, recurrent status and malignant histology. These classes more accurately identified recurrent tumors relative to Ki-67 index and extent of surgical resection, and highlight substantial expression heterogeneity between meningiomas. These data offer the most complete description of the genomic landscape of meningiomas, and provide broad genomic information that may be used to further stratify meningioma patients into prognostic risk groups.

Detection of pathogenic copy number variants in children with idiopathic intellectual disability using 500 K SNP array genomic hybridization

BACKGROUND

Array genomic hybridization is being used clinically to detect pathogenic copy number variants in children with intellectual disability and other birth defects. However, there is no agreement regarding the kind of array, the distribution of probes across the genome, or the resolution that is most appropriate for clinical use.

RESULTS

We performed 500 K Affymetrix GeneChip array genomic hybridization in 100 idiopathic intellectual disability trios, each comprised of a child with intellectual disability of unknown cause and both unaffected parents. We found pathogenic genomic imbalance in 16 of these 100 individuals with idiopathic intellectual disability. In comparison, we had found pathogenic genomic imbalance in 11 of 100 children with idiopathic intellectual disability in a previous cohort who had been studied by 100 K GeneChip array genomic hybridization. Among 54 intellectual disability trios selected from the previous cohort who were re-tested with 500 K GeneChip array genomic hybridization, we identified all 10 previously-detected pathogenic genomic alterations and at least one additional pathogenic copy number variant that had not been detected with 100 K GeneChip array genomic hybridization. Many benign copy number variants, including one that was de novo, were also detected with 500 K array genomic hybridization, but it was possible to distinguish the benign and pathogenic copy number variants with confidence in all but 3 (1.9%) of the 154 intellectual disability trios studied.

CONCLUSION

Affymetrix GeneChip 500 K array genomic hybridization detected pathogenic genomic imbalance in 10 of 10 patients with idiopathic developmental disability in whom 100 K GeneChip array genomic hybridization had found genomic imbalance, 1 of 44 patients in whom 100 K GeneChip array genomic hybridization had found no abnormality, and 16 of 100 patients who had not previously been tested. Effective clinical interpretation of these studies requires considerable skill and experience.

Hybridization modeling of oligonucleotide SNP arrays for accurate DNA copy number estimation

Affymetrix SNP arrays have been widely used for single-nucleotide polymorphism (SNP) genotype calling and DNA copy number variation inference. Although numerous methods have achieved high accuracy in these fields, most studies have paid little attention to the modeling of hybridization of probes to off-target allele sequences, which can affect the accuracy greatly. In this study, we address this issue and demonstrate that hybridization with mismatch nucleotides (HWMMN) occurs in all SNP probe-sets and has a critical effect on the estimation of allelic concentrations (ACs). We study sequence binding through binding free energy and then binding affinity, and develop a probe intensity composite representation (PICR) model. The PICR model allows the estimation of ACs at a given SNP through statistical regression. Furthermore, we demonstrate with cell-line data of known true copy numbers that the PICR model can achieve reasonable accuracy in copy number estimation at a single SNP locus, by using the ratio of the estimated AC of each sample to that of the reference sample, and can reveal subtle genotype structure of SNPs at abnormal loci. We also demonstrate with HapMap data that the PICR model yields accurate SNP genotype calls consistently across samples, laboratories and even across array platforms.

A genome-wide survey of copy number variations in Han Chinese residing in Taiwan

Copy number variation (CNV) is a form of DNA sequence variation in the human genome. CNVs can affect expression of nearby and distant genes, and some of them might cause certain phenotypic differences. CNVs vary slightly in location and frequency among different populations. Because currently-available CNV information from Asian population was limited to fewer small-scale studies with only dozens of subjects, a high-resolution CNV survey was conducted using a large number of Han Chinese in this study. The Illumina HumanMap550K single-nucleotide polymorphism array was used to identify CNVs from 813 unrelated Han Chinese residing in Taiwan. A total of 365 CNV regions were identified in this population, and the average size of the CNV regions was 235 kb (covering a total of 2.86% of the human genome), and 67 (18.4%) were newly-discovered CNV regions. Two hundred and seventy-nine CNV regions (76%) were verified from 304 randomly-selected samples by Affymetrix 500K GeneChip and qPCR experiments. These regions contain 1029 genes, some of which are associated with diseases. Consistent with previous studies, most CNVs were rare structural variations in the human genome, and only 64 regions (17.5%) had a CNV allele frequency greater than 1%. Our discovery of 67 new CNV regions indicates that previous CNV coverage of the human genome is incomplete and there is diversity among different ethnic populations. The comprehensive knowledge of CNVs in the human genome is very important and useful in further genetic studies.

Progress in searching for the febrile seizure susceptibility genes

Febrile seizures (FS) represent the most common form of childhood seizures. They affect 2-5% of infants in the Caucasian population and are even more common in the Japanese population, affecting 6-9% of infants. Some familial FS are associated with a wide variety of afebrile seizures. Generalized epilepsy with febrile seizures plus (GEFS+) is a familial epilepsy syndrome with a spectrum of phenotypes including FS, atypical FS (FS+) and afebrile seizures. A significant genetic component exists for susceptibility to FS and GEFS+: extensive genetic studies have shown that at least nine loci are responsible for FS. Furthermore, mutations in the voltage-gated sodium channel subunit genes (SCN1A, SCN2A and SCN1B) and the GABA(A) receptor subunit genes (GABRG2 and GABRD) have been identified in GEFS+. However, the causative genes have not been identified in most patients with FS or GEFS+. Common forms of FS are genetically complex disorders believed to be influenced by variations in several susceptibility genes. Recently, several association studies on FS have been reported, but the results vary among different groups and no consistent or convincing FS susceptibility gene has emerged. Herein, we review the genetic data reported in FS, including the linkage analysis, association studies, and genetic abnormalities found in the FS-related disorders such as GEFS+ and severe myoclonic epilepsy in infancy.

A large-scale survey of genetic copy number variations among Han Chinese residing in Taiwan

BACKGROUND

Copy number variations (CNVs) have recently been recognized as important structural variations in the human genome. CNVs can affect gene expression and thus may contribute to phenotypic differences. The copy number inferring tool (CNIT) is an effective hidden Markov model-based algorithm for estimating allele-specific copy number and predicting chromosomal alterations from single nucleotide polymorphism microarrays. The CNIT algorithm, which was constructed using data from 270 HapMap multi-ethnic individuals, was applied to identify CNVs from 300 unrelated Han Chinese individuals in Taiwan.

RESULTS

Using stringent selection criteria, 230 regions with variable copy numbers were identified in the Han Chinese population; 133 (57.83%) had been reported previously, 64 displayed greater than 1% CNV allele frequency. The average size of the CNV regions was 322 kb (ranging from 1.48 kb to 5.68 Mb) and covered a total of 2.47% of the human genome. A total of 196 of the CNV regions were simple deletions and 27 were simple amplifications. There were 449 genes and 5 microRNAs within these CNV regions; some of these genes are known to be associated with diseases.

CONCLUSION

The identified CNVs are characteristic of the Han Chinese population and should be considered when genetic studies are conducted. The CNV distribution in the human genome is still poorly characterized, and there is much diversity among different ethnic populations.

A novel procedure for genotyping of single nucleotide polymorphisms in trisomy with genomic DNA and the invader assay

Individuals with trisomy 21 display complex phenotypes with differing degrees of severity. Numerous reliable methods have been established to diagnose the initial trisomy in these patients, but the identification and characterization of the genetic basis of the phenotypic variation in individuals with trisomy remains challenging. To date, methods that can accurately determine genotypes in trisomic DNA samples are expensive, require specialized equipment and complicated analyses. Here we report proof-of-concept results for an Invader® assay-based genotyping procedure that can determine SNP genotypes in trisomic genomic DNA samples in a simple and cost-effective manner. The procedure requires only two experimental steps: a real-time measurement of the fluorescent Invader® signal and analysis with a specifically designed clustering algorithm. The approach was tested using genomic DNA samples from 23 individuals with trisomy 21, and results were compared to genotypes previously determined with pyrosequencing. Additional assays for 15 SNPs were tested in a set of 21 DNA samples to assess assay performance. Our method successfully identified the correct SNP genotypes for the trisomic genomic DNA samples tested, and thus provides an alternative to determine SNP genotypes in trisomic DNA samples for subsequent association studies in patients with Down syndrome and other trisomies.

Ploidy status and copy number aberrations in primary glioblastomas defined by integrated analysis of allelic ratios, signal ratios and loss of heterozygosity using 500K SNP Mapping Arrays

Background

Genomic hybridization platforms, including BAC-CGH and genotyping arrays, have been used to estimate chromosome copy number (CN) in tumor samples by detecting the relative strength of genomic signal. The methods rely on the assumption that the predominant chromosomal background of the samples is diploid, an assumption that is frequently incorrect for tumor samples. In addition to generally greater resolution, an advantage of genotyping arrays over CGH arrays is the ability to detect signals from individual alleles, allowing estimation of loss-of-heterozygosity (LOH) and allelic ratios to enhance the interpretation of copy number alterations. Copy number events associated with LOH potentially have the same genetic consequences as deletions.

Results

We have utilized allelic ratios to detect patterns that are indicative of higher ploidy levels. An integrated analysis using allelic ratios, total signal and LOH indicates that many or most of the chromosomes from 24 glioblastoma tumors are in fact aneuploid. Some putative whole-chromosome losses actually represent trisomy, and many apparent sub-chromosomal losses are in fact relative losses against a triploid or tetraploid background.

Conclusion

These results suggest a re-interpretation of previous findings based only on total signal ratios. One interesting observation is that many single or multiple-copy deletions occur at common putative tumor suppressor sites subsequent to chromosomal duplication; these losses do not necessarily result in LOH, but nonetheless occur in conspicuous patterns. The 500 K Mapping array was also capable of detecting many sub-mega base losses and gains that were overlooked by CGH-BAC arrays, and was superior to CGH-BAC arrays in resolving regions of complex CN variation.

Genome-wide copy number analysis using copy number inferring tool (CNIT) and DNA pooling

Copy number variation (CNV) has become an important genomic structure element in the human population, and some CNVs are related to specific traits and diseases. Moreover, analysis of human genomes has been potentiated by the use of high-resolution microarrays that assess single nucleotide polymorphisms (SNPs). Although many programs have been designed to analyze data from Affymetrix SNP microarrays, they all have high false-positive rates (FPRs) in copy number (CN) analyses. Copy number analysis tool (CNAT) 4.0 is a recently developed program that offers improved CN estimation, but small amplifications and deletions are lost when using the smoothing procedure. Here, we propose a copy number inferring tool (CNIT) algorithm for the 100K SNP microarray to investigate CNVs at 29.6-kb resolution. CNIT estimated SNP allelic and total CN with reliable P values based on intensity data. In addition, the hidden Markov model (HMM) method was applied to predict regions having altered CN by considering contiguous SNPs. Based on a CN analysis of 23 unrelated Taiwanese and 30 HapMap Centre d'Etude du Polymorphisme Humain (CEPH) trios, CNIT showed higher accuracy and power than other programs. The FPRs and false-negative rates (FNRs) of CNIT were 0.1% and 0.16%, respectively. CNIT also showed better sensitivity for detecting small amplifications and deletions. Furthermore, DNA pooling of 10 and 30 normal unrelated individuals were applied to the 100K SNP microarray, respectively, and 12 common CN-variable regions were identified, suggesting that DNA pooling can be applied to discover common CNVs.

A mutation in the canalicular phospholipid transporter gene, ABCB4, is associated with cholestasis, ductopenia, and cirrhosis in adults

Cholestatic liver disease (CLD) is a major cause of progressive liver damage and liver failure. Several forms of biliary cirrhosis are caused by mutations in specific genes. We sought to identify a genetic defect in a family with CLD impossible to assign to a distinct pathogenic entity. Clinical and histopathological characterization of the family members, microarray-based single-nucleotide polymorphism genotyping, and analysis of candidate genes were performed. Among six of 11 siblings severely affected by idiopathic CLD in a family from a population isolate in Transylvania, three died of cirrhosis (aged 5, 7, and 43 years) and three had adult-onset disease with small duct cholangiopathy, including ductopenia. Others were mildly affected and experienced intrahepatic cholestasis of pregnancy, miscarriages, or stillbirth. Pedigree studies revealed distant parental consanguinity. Genome-wide linkage analysis and autozygosity mapping yielded a single maximal lod-score of 3.88 on chromosome 7q21.1-7q22, excluding other genomic loci. Sequencing of ABCB4 at this locus revealed a novel missense mutation c.2362C>T (p.Arg788Trp) which cosegregated with severity of disease. Bile from a mutation homozygote showed a reduced phosphatidylcholine/bile acid ratio, consistent with reduced ABCB4 phosphatidylcholine transport activity.

CONCLUSION

We show that a missense mutation in ABCB4 is a cause for ductopenic CLD in adulthood. Allelic status correlated with severity of liver disease ranging from intrahepatic cholestasis of pregnancy through fibrosis to cirrhosis and death in childhood and adulthood. Mutational analysis of ABCB4 should be generally considered in all patients with cholestatic liver disease of unknown etiology regardless of age and onset of disease.

Identification of copy number alterations and its association with pathological features in clear cell and papillary RCC

We report and characterize the copy number alterations (CNAs) in 35 clear cell and 12 papillary renal cell carcinomas (RCC) using Affymetrix 100K SNP arrays. Novel gain and loss regions are identified in both subtypes. In addition, statistically significant CNA are detected and associated with the pathological features: VHL mutation status, tumor grades, and sarcomatoid component in clear cell RCC and in types 1 and 2 of papillary RCC. Florescence in situ hybridization confirmed the copy number gain in the transforming growth factor, beta-induced gene (TGFBI), which is a possible oncogene for clear cell RCC.

Clinical and molecular delineation of the 17q21.31 microdeletion syndrome

BACKGROUND

The chromosome 17q21.31 microdeletion syndrome is a novel genomic disorder that has originally been identified using high resolution genome analyses in patients with unexplained mental retardation.

AIM

We report the molecular and/or clinical characterisation of 22 individuals with the 17q21.31 microdeletion syndrome.

RESULTS

We estimate the prevalence of the syndrome to be 1 in 16,000 and show that it is highly underdiagnosed. Extensive clinical examination reveals that developmental delay, hypotonia, facial dysmorphisms including a long face, a tubular or pear-shaped nose and a bulbous nasal tip, and a friendly/amiable behaviour are the most characteristic features. Other clinically important features include epilepsy, heart defects and kidney/urologic anomalies. Using high resolution oligonucleotide arrays we narrow the 17q21.31 critical region to a 424 kb genomic segment (chr17: 41046729-41470954, hg17) encompassing at least six genes, among which is the gene encoding microtubule associated protein tau (MAPT). Mutation screening of MAPT in 122 individuals with a phenotype suggestive of 17q21.31 deletion carriers, but who do not carry the recurrent deletion, failed to identify any disease associated variants. In five deletion carriers we identify a <500 bp rearrangement hotspot at the proximal breakpoint contained within an L2 LINE motif and show that in every case examined the parent originating the deletion carries a common 900 kb 17q21.31 inversion polymorphism, indicating that this inversion is a necessary factor for deletion to occur (p<10(-5)).

CONCLUSION

Our data establish the 17q21.31 microdeletion syndrome as a clinically and molecularly well recognisable genomic disorder.

MLPA and MAPH: sensitive detection of deletions and duplications

The detection of quantitative changes in genomic DNA, i.e., deletions and duplications or so called Copy Number Variants (CNV), is an important element of a complete mutation screening strategy. However, because of practical difficulties, screening for quantitative changes in genomic DNA is often ignored. Hitherto, the techniques available were technically challenging and laborious and thus too costly to be applied on a routine basis. The development of MAPH (Multiplex Amplifiable Probe Hybridization) and more recently MLPA (Multiplex Ligation-dependent Probe Amplification) have revived interest in the detection of deletions and duplications, primarily due to the simplicity and flexibility of these two approaches. Compared to previous technologies, e.g., Southern blotting, fluorescence in situ hybridization (FISH), quantitative PCR (qPCR), and breakpoint PCR, they have some clear advantages, including high resolution, high throughput, amenability to multiplexing, and simplicity.

DNA copy number analysis in gastrointestinal stromal tumors using gene expression microarrays

We report a method, Expression-Microarray Copy Number Analysis (ECNA) for the detection of copy number changes using Affymetrix Human Genome U133 Plus 2.0 arrays, starting with as little as 5 ng input genomic DNA. An analytical approach was developed using DNA isolated from cell lines containing various X-chromosome numbers, and validated with DNA from cell lines with defined deletions and amplifications in other chromosomal locations. We applied this method to examine the copy number changes in DNA from 5 frozen gastrointestinal stromal tumors (GIST). We detected known copy number aberrations consistent with previously published results using conventional or BAC-array CGH, as well as novel changes in GIST tumors. These changes were concordant with results from Affymetrix 100K human SNP mapping arrays. Gene expression data for these GIST samples had previously been generated on U133A arrays, allowing us to explore correlations between chromosomal copy number and RNA expression levels. One of the novel aberrations identified in the GIST samples, a previously unreported gain on 1q21.1 containing the PEX11B gene, was confirmed in this study by FISH and was also shown to have significant differences in expression pattern when compared to a control sample. In summary, we have demonstrated the use of gene expression microarrays for the detection of genomic copy number aberrations in tumor samples. This method may be used to study copy number changes in other species for which RNA expression arrays are available, e.g. other mammals, plants, etc., and for which SNPs have not yet been mapped.

Resolution of trisomic mosaicism in prenatal diagnosis: estimated performance of a 50K SNP microarray

OBJECTIVE

To evaluate the ability of a DNA single nucleotide polymorphism (SNP) microarray to detect chromosome mosaicism for trisomy in prenatal samples in order to compare this with conventional cytogenetics.

METHOD

We created a dilution series of mock mosaic samples, by mixing measured amounts of fibroblast cells containing trisomy 8 from a male with aliquots of cells with a normal female karyotype. DNAs were extracted from these mosaic mixtures, then analysed on the Affymetrix 50K Xba SNP chip. Duplicate aliquots of each mosaic sample were probed using interphase FISH, with centromeric probes for chromosomes X, Y and 8, to estimate independently the proportion of male trisomy 8 in each sample. Data from the arrays were analysed using publicly available analysis tools. Statistical calculations were then performed using a Student's t-test to determine if there was a significant difference between the copy numbers of each chromosome.

RESULTS

These experiments using the Affymetrix 50K Xba SNP microarray showed mosaicism to be obvious at 20% and with additional statistical calculations, the lower limit for detection is about 10%.

CONCLUSION

The SNP microarray platform tested can detect mosaicism for trisomy in prenatal samples at levels comparable with conventional cytogenetic techniques in routine use.

Representational oligonucleotide microarray analysis (ROMA) and comparison of binning and change-point methods of analysis: application to detection of del22q11.2 (DiGeorge) syndrome

DiGeorge (del22q11.2) syndrome is estimated to occur in 1:4,000 births, is the most common contiguous-gene deletion syndrome in humans, and is caused by autosomal dominant deletions in the 22q11.2 DiGeorge syndrome critical region (DGCR). Multiple microarray methods have been developed recently for analyzing such copy number changes, but data analysis and accurate deletion detection remains challenging. Clinical use of these microarray methods would have many advantages, particularly when the possibility of a chromosomal disorder cannot be determined simply on the basis of history and physical examination data alone. We investigated the use of the microarray technique, representational oligonucleotide microarray analysis (ROMA), in the detection of del22q11.2 syndrome. Genomic DNA was isolated from three well-characterized cell lines with 22q11.2 DGCR deletions and from the blood of a patient suspected of having del22q11.2 syndrome, and analyzed using both the binning and change-point model algorithms. Though the 22q11.2 deletion was easily identified with either method, change-point models provide clearer identification of deleted regions, with the potential for fewer false-positive results. For circumstances in which a clear, a priori, copy-number change hypothesis is not present, such as in many clinical samples, change-point methods of analysis may be easier to interpret.

Visualization of uniparental inheritance, Mendelian inconsistencies, deletions, and parent of origin effects in single nucleotide polymorphism trio data with SNPtrio

A variety of alterations occur in chromosomal DNA, many of which can be detected using high density single nucleotide polymorphism (SNP) microarrays. These include deletions and duplications (assessed by observing changes in copy number) and regions of homozygosity. The analysis of SNP data from trios can provide an additional category of information about the nature and origin of inheritance patterns, including uniparental disomy (UPD), loss of transmitted allele (LTA), and nonparental relationship. The main purpose of SNPtrio is to locate regions of uniparental inheritance (UPI) and Mendelian inconsistency (MI), identify the type (paternal vs. maternal, iso- vs. hetero-), and assess the associated statistical probability of occurrence by chance. SNPtrio's schema permits the identification of hemizygous or homozygous deletions as well as UPD. We validated the performance of SNPtrio on three platforms (Affymetrix 10 K and 100 K arrays and Illumina 550 K arrays) using SNP data obtained from DNA samples of patients known to have UPD and diagnosed with Prader-Willi syndrome, Angelman syndrome, Beckwith-Wiedemann syndrome, pseudohypoparathyroidism, and a complex chromosome 2 abnormality. We further validated SNPtrio using DNA from patients previously shown to have microdeletions that were verified by fluorescence in situ hybridization (FISH). SNPtrio successfully identified previously known UPD and deletion regions, and generated associated probability values. SNPtrio analysis of trisomy 21 (Down syndrome) cases and their parents permitted identification of the parent of origin of the extra chromosomal copy. SNPtrio is freely accessible at http://pevsnerlab.kennedykrieger.org/SNPtrio.htm (Last accessed: 20 June 2007).

Side effects of genome structural changes

The first extensive catalog of structural human variation was recently released. It showed that large stretches of genomic DNA that vary considerably in copy number were extremely abundant. Thus it is conceivable that they play a major role in functional variation. Consistently, genomic insertions and deletions were shown to contribute to phenotypic differences by modifying not only the expression levels of genes within the aneuploid segments but also of normal copy-number neighboring genes. In this report, we review the possible mechanisms behind this latter effect.

Assessment of algorithms for high throughput detection of genomic copy number variation in oligonucleotide microarray data

Background

Genomic deletions and duplications are important in the pathogenesis of diseases, such as cancer and mental retardation, and have recently been shown to occur frequently in unaffected individuals as polymorphisms. Affymetrix GeneChip whole genome sampling analysis (WGSA) combined with 100 K single nucleotide polymorphism (SNP) genotyping arrays is one of several microarray-based approaches that are now being used to detect such structural genomic changes. The popularity of this technology and its associated open source data format have resulted in the development of an increasing number of software packages for the analysis of copy number changes using these SNP arrays.

Results

We evaluated four publicly available software packages for high throughput copy number analysis using synthetic and empirical 100 K SNP array data sets, the latter obtained from 107 mental retardation (MR) patients and their unaffected parents and siblings. We evaluated the software with regards to overall suitability for high-throughput 100 K SNP array data analysis, as well as effectiveness of normalization, scaling with various reference sets and feature extraction, as well as true and false positive rates of genomic copy number variant (CNV) detection.

Conclusion

We observed considerable variation among the numbers and types of candidate CNVs detected by different analysis approaches, and found that multiple programs were needed to find all real aberrations in our test set. The frequency of false positive deletions was substantial, but could be greatly reduced by using the SNP genotype information to confirm loss of heterozygosity.

Copy-number variations measured by single-nucleotide-polymorphism oligonucleotide arrays in patients with mental retardation

Whole-genome analysis using high-density single-nucleotide-polymorphism oligonucleotide arrays allows identification of microdeletions, microduplications, and uniparental disomies. We studied 67 children with unexplained mental retardation with normal karyotypes, as assessed by G-banded chromosome analyses. Their DNAs were analyzed with Affymetrix 100K arrays. We detected 11 copy-number variations that most likely are causative of mental retardation, because they either arose de novo (9 cases) and/or overlapped with known microdeletions (2 cases). The eight deletions and three duplications varied in size from 200 kb to 7.5 Mb. Of the 11 copy-number variations, 5 were flanked by low-copy repeats. Two of those, on chromosomes 15q25.2 and Xp22.31, have not been described before and have a high probability of being causative of new deletion and duplication syndromes, respectively. In one patient, we found a deletion affecting only a single gene, MBD5, which codes for the methyl-CpG-binding domain protein 5. In addition to the 67 children, we investigated 4 mentally retarded children with apparent balanced translocations and detected four deletions at breakpoint regions ranging in size from 1.1 to 14 Mb.

Oligonucleotide arrays for high-resolution analysis of copy number alteration in mental retardation/multiple congenital anomalies

Genetic diseases arising from microdeletions and microduplications lead to copy number alterations of genomic regions containing one or more genes. Clinically, these rearrangements may be detected by routine cytogenetic testing, which may include karyotype analysis, subtelomeric analysis with fluorescence in situ hybridization, and/or fluorescence in situ hybridization directed at known chromosomal rearrangement-based disorders. The major limitations of these tests are low resolution and limited coverage of the genome. Array-based comparative genomic hybridization has recently become a widely used approach in the genome-wide analysis of copy number alterations in children with mental retardation and/or multiple congenital anomalies. Oligonucleotide-based arrays provide a genome-wide coverage at a much higher resolution than microarrays currently used in clinical diagnostics, greatly improving the rate of detection of submicroscopic copy number alterations in children with mental retardation and/or multiple congenital anomalies.

Application of a BRAF pyrosequencing assay for mutation detection and copy number analysis in malignant melanoma

Mutations in the BRAF gene are found in the majority of cutaneous malignant melanomas and subsets of other tumors. These mutations lead to constitutive activation of BRAF with increased downstream ERK (extracellular signal-regulated kinase) signaling; therefore, the development of RAF kinase inhibitors for targeted therapy is being actively pursued. A methodology that allows sensitive, cost-effective, high-throughput analysis of BRAF mutations will be needed to triage patients for specific molecular-based therapies. Pyrosequencing is a high-throughput, sequencing-by-synthesis method that is particularly useful for analysis of single nucleotide polymorphisms or hotspot mutations. Mutational analysis of BRAF is highly amenable to pyrosequencing because the majority of mutations in this gene localize to codons 600 and 601 and consist of single or dinucleotide substitutions. In this study, DNAs from a panel of melanocyte cell lines, melanoma cell lines, and melanoma tumors were used to validate a pyrosequencing assay to detect BRAF mutations. The assay demonstrates high accuracy and precision for detecting common and variant exon 15 BRAF mutations. Further, comparison of pyrosequencing data with 100K single nucleotide polymorphism microarray data allows characterization of BRAF amplification events that may accompany BRAF mutation. Pyro-sequencing serves as an excellent platform for BRAF genotyping of tumors from patients entering clinical trial.

Clonality of the stem cell compartment during evolution of myelodysplastic syndromes and other bone marrow failure syndromes

Clonal hematopoiesis, observed in certain forms of marrow failure including aplastic anemia (AA), may be due to stem cell depletion. Alternatively, oligoclonality may be a result of recruitment of a preexisting defective clone, such as in paroxysmal nocturnal hemoglobinuria (PNH) or myelodysplastic syndromes (MDS). In PNH, exogenous permissive factors may be required for dominance of the abnormal clone, while in MDS, stem cells undergo transformation steps leading to a growth advantage. Stem or multipotent progenitor cell involvement in PNH is evidenced by long-term persistence of a clonal defect and its presence in all blood cells. In MDS, some clonal aberrations may have a 'founder-effect' and additional defects are secondary. Metaphase cytogenetics measures the proportion of clonal cells within dividing progenitor but not mature cells. Owing to low resolution, lesions can be found in only approximately 50% of MDS patients. This shortcoming may be overcome by application of newer technologies such as comparative genomic hybridization and SNP array-based karyotyping (SNP-A). SNP-A facilitates identification of cryptic lesions in bone marrow failure patients with normal or abnormal cytogenetics and allows for detection of loss of heterozygosity as a result of uniparental disomy, a lesion frequently found in MDS.

Chromosomal abnormality at 6p25.1-25.3 identifies a susceptibility locus for hypothalamic hamartoma associated with epilepsy

The pathogenesis of hypothalamic hamartoma (HH) associated with epilepsy is unknown. We have identified an individual with HH and refractory epilepsy exhibiting subtle dysmorphic features. High-resolution karyotype identified a duplication of the terminal end of 6p (6p25.1-25.3), confirmed by fluorescent in situ-hybridization (FISH). Copy number analysis with high-density (250K) single nucleotide polymorphism (SNP) genotyping microarrays characterized the abnormality as a series of amplified regions between 1.4 Mb and 10.2 Mb, with a small tandem deletion from 8.8 Mb to 9.7 Mb. There are 38 RefSeq genes within the duplicated regions, and no known coding sequences within the deletion. This unique patient helps identify 6p25.1-25.3 as a possible susceptibility locus for sporadic HH.

DNA copy number alterations in prostate cancers: a combined analysis of published CGH studies

BACKGROUND

Identifying genomic regions that are commonly deleted or gained in neoplastic cells is an important approach to identify tumor suppressor genes and oncogenes. Studies in the last two decades have identified a number of common DNA copy number alterations in prostate cancer. However, because of various sample sizes, diverse tumor types and sources, as well as a variety of detection methods with various sensitivities and resolutions, it is difficult to summarize and fully interpret the overall results.

METHODS

We performed a combined analysis of all published comparative genomic hybridization (CGH) studies of prostate cancer and estimated the frequency of alterations across the genome for all tumors, as well as in advanced and localized tumors separately. A total of 41 studies examining 872 cancers were included in this study.

RESULTS

The frequency of deletions and gains were estimated in all tumors, as well as in advanced and localized tumors. Eight deleted and five gained regions were found in more than 10% of the prostate tumors. An additional six regions were commonly deleted and seven were commonly gained in advanced tumors. While 8p was the most common location of deletion, occurring in about a third of all tumors and about half of advanced tumors, 8q was the most commonly gained region, affecting about a quarter of all tumors and about half of all advanced tumors.

CONCLUSIONS

The large number of tumors examined in this combined analysis provides better estimates of the frequency of specific alterations in the prostate cancer cell genome, and offers important clues for prioritizing efforts to identify tumor suppressor genes and oncogenes in these altered regions.

Comprehensive cytogenetic characterization of an esthesioneuroblastoma

Esthesioneuroblastoma is a malignant neuroectodermal tumor originating from olfactory epithelial cells in the nasal vault. Due to the rarity of this tumor entity, cytogenetic data are very limited. Therefore, we performed comprehensive cytogenetic analyses of an esthesioneuroblastoma, Hyam's grade III-IV, using trypsin-Giemsa staining (GTG banding), multicolor fluorescence in situ hybridization (M-FISH), and locus-specific FISH complemented by molecular karyotyping using high-density single nucleotide polymorphism arrays. GTG banding of 25 metaphases revealed 54 structural intrachromosomal aberrations, predominantly located on 2q, 6q, 21q, and 22q, which were confirmed by FISH analysis. Interestingly, we found two novel, so far not described deletions, del(2)(q37) and del(21)(q22). Using GTG banding, locus-specific FISH, and M-FISH, we detected numeric changes of chromosomes 5, 17, 19, and 22, as well as trisomy 8 at low frequency. Applying SNP array karyotyping, we confirmed the chromosomal aberrations del(2)(q37.3), del(3)(q27.2), del(10)(q26.11), chromosomal imbalance on 17q, del(21)(q22), and revealed a number of so far unknown aberrations (gain of 2q14.3, 13q33.3, and 13q34). While the cytogenetically revealed low frequency mosaic del(6)(q22q24) was not visible using SNP array karyotyping, some of the smaller imbalances (SNP array data) could not have been detected by classic cytogenetic analysis. Therefore, our study supports the usefulness of applying complementary methods for cytogenetic analysis.

Genome-wide copy number profiling on high-density bacterial artificial chromosomes, single-nucleotide polymorphisms, and oligonucleotide microarrays: a platform comparison based on statistical power analysis

Recently, comparative genomic hybridization onto bacterial artificial chromosome (BAC) arrays (array-based comparative genomic hybridization) has proved to be successful for the detection of submicroscopic DNA copy-number variations in health and disease. Technological improvements to achieve a higher resolution have resulted in the generation of additional microarray platforms encompassing larger numbers of shorter DNA targets (oligonucleotides). Here, we present a novel method to estimate the ability of a microarray to detect genomic copy-number variations of different sizes and types (i.e. deletions or duplications). We applied our method, which is based on statistical power analysis, to four widely used high-density genomic microarray platforms. By doing so, we found that the high-density oligonucleotide platforms are superior to the BAC platform for the genome-wide detection of copy-number variations smaller than 1 Mb. The capacity to reliably detect single copy-number variations below 100 kb, however, appeared to be limited for all platforms tested. In addition, our analysis revealed an unexpected platform-dependent difference in sensitivity to detect a single copy-number loss and a single copy-number gain. These analyses provide a first objective insight into the true capacities and limitations of different genomic microarrays to detect and define DNA copy-number variations.

QuantiSNP: an Objective Bayes Hidden-Markov Model to detect and accurately map copy number variation using SNP genotyping data

Array-based technologies have been used to detect chromosomal copy number changes (aneuploidies) in the human genome. Recent studies identified numerous copy number variants (CNV) and some are common polymorphisms that may contribute to disease susceptibility. We developed, and experimentally validated, a novel computational framework (QuantiSNP) for detecting regions of copy number variation from BeadArray SNP genotyping data using an Objective Bayes Hidden-Markov Model (OB-HMM). Objective Bayes measures are used to set certain hyperparameters in the priors using a novel re-sampling framework to calibrate the model to a fixed Type I (false positive) error rate. Other parameters are set via maximum marginal likelihood to prior training data of known structure. QuantiSNP provides probabilistic quantification of state classifications and significantly improves the accuracy of segmental aneuploidy identification and mapping, relative to existing analytical tools (Beadstudio, Illumina), as demonstrated by validation of breakpoint boundaries. QuantiSNP identified both novel and validated CNVs. QuantiSNP was developed using BeadArray SNP data but it can be adapted to other platforms and we believe that the OB-HMM framework has widespread applicability in genomic research. In conclusion, QuantiSNP is a novel algorithm for high-resolution CNV/aneuploidy detection with application to clinical genetics, cancer and disease association studies.

Deletion mapping in Xp21 for patients with complex glycerol kinase deficiency using SNP mapping arrays

Infantile or complex glycerol kinase deficiency (cGKD) is a contiguous gene deletion syndrome caused by a loss of GK (MIM# 300474), along with its neighboring genes, Duchenne muscular dystrophy (DMD; MIM# 300377) and/or Nuclear Receptor Subfamily 0, Group B, Member 1 (NR0B1; MIM# 300473). Patients with cGKD present with glyceroluria and hyperglycerolemia in association with DMD and/or adrenal hypoplasia congenita (AHC). The purpose of these investigations was to determine whether the Affymetrix GeneChip Mapping Array (SNP chip) could be utilized to detect and map breakpoints in patients with cGKD. Genomic DNAs from several primary lymphoblastoid cell lines from patients with cGKD were analyzed on the Affymetrix platform. The Affymetrix SNP chip is a high-density oligonucleotide array that allows a standardized, parallel interrogation of thousands of SNPs across the entire genome (except for the Y chromosome). Analysis of the array features' hybridization intensities enabled clear delineation of the patient deletions with a high degree of confidence. Many of these patient deletions had been mapped by PCR and their breakpoints confirmed by sequencing. This study demonstrates the utility of the Affymetrix Mapping GeneChips for molecular cytogenetic analysis, beyond the SNP genotyping for which the arrays were initially designed. With one out of 160 live births (approximately 25,000 U.S. neonates annually) reported to have cytogenetic disorders, we envision a significant need for such a standardized platform to carry out rapid, high-throughput, genomic analyses for molecular cytogenetics applications.