Reference-unbiased copy number variant analysis using CGH microarrays

Comparative genomic hybridization (CGH) microarrays have been used to determine copy number variations (CNVs) and their effects on complex diseases. Detection of absolute CNVs independent of genomic variants of an arbitrary reference sample has been a critical issue in CGH array experiments. Whole genome analysis using massively parallel sequencing with multiple ultra-high resolution CGH arrays provides an opportunity to catalog highly accurate genomic variants of the reference DNA (NA10851). Using information on variants, we developed a new method, the CGH array reference-free algorithm (CARA), which can determine reference-unbiased absolute CNVs from any CGH array platform. The algorithm enables the removal and rescue of false positive and false negative CNVs, respectively, which appear due to the effects of genomic variants of the reference sample in raw CGH array experiments. We found that the CARA remarkably enhanced the accuracy of CGH array in determining absolute CNVs. Our method thus provides a new approach to interpret CGH array data for personalized medicine.     

Bisulfite-modified target DNA array for aberrant methylation analysis

Aberrant DNA methylation of CpG islands is among the earliest and most frequent alterations in cancer. It is of great importance to develop simple and high-throughput methods of methylation analysis for earlier cancer diagnosis or the detection of recurrence. In this study, bisulfite-modified target DNA arrays were prepared on positively charged nylon membrane with two different procedures: fixing PCR products and fixing genomic DNA. First, a bisulfite PCR product array was prepared through fixing PCR products amplified in bisulfite sequencing primers from the bisulfite-modified genomic DNA of different clinical samples on membrane. Furthermore, bisulfite-modified genomic DNA of the different samples was directly fixed on membrane to fabricate bisulfite genomic DNA arrays. The two kinds of arrays were hybridized by probes labeled with digoxigenin, and the hybridization signals were obtained through chemiluminescent detection. The methylation statuses of the IGFBP7 gene for breast tumor and normal tissue samples and for normal human blood cell samples were detected successfully by the two procedures. It was shown that the methods are reliable and sensitive and that they have high potential in screening molecular methylation markers from a large number of clinical samples.     

A powerful tool for genome analysis in maize: development and evaluation of the high density 600 k SNP genotyping array

Background

High density genotyping data are indispensable for genomic analyses of complex traits in animal and crop species. Maize is one of the most important crop plants worldwide, however a high density SNP genotyping array for analysis of its large and highly dynamic genome was not available so far.

Results

We developed a high density maize SNP array composed of 616,201 variants (SNPs and small indels). Initially, 57 M variants were discovered by sequencing 30 representative temperate maize lines and then stringently filtered for sequence quality scores and predicted conversion performance on the array resulting in the selection of 1.2 M polymorphic variants assayed on two screening arrays. To identify high-confidence variants, 285 DNA samples from a broad genetic diversity panel of worldwide maize lines including the samples used for sequencing, important founder lines for European maize breeding, hybrids, and proprietary samples with European, US, semi-tropical, and tropical origin were used for experimental validation. We selected 616 k variants according to their performance during validation, support of genotype calls through sequencing data, and physical distribution for further analysis and for the design of the commercially available Affymetrix® Axiom® Maize Genotyping Array. This array is composed of 609,442 SNPs and 6,759 indels. Among these are 116,224 variants in coding regions and 45,655 SNPs of the Illumina® MaizeSNP50 BeadChip for study comparison. In a subset of 45,974 variants, apart from the target SNP additional off-target variants are detected, which show only a minor bias towards intermediate allele frequencies. We performed principal coordinate and admixture analyses to determine the ability of the array to detect and resolve population structure and investigated the extent of LD within a worldwide validation panel.

Conclusions

The high density Affymetrix® Axiom® Maize Genotyping Array is optimized for European and American temperate maize and was developed based on a diverse sample panel by applying stringent quality filter criteria to ensure its suitability for a broad range of applications. With 600 k variants it is the largest currently publically available genotyping array in crop species.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-823) contains supplementary material, which is available to authorized users.     

Infinium Assay for Large-scale SNP Genotyping Applications

Genotyping variants in the human genome has proven to be an efficient method to identify genetic associations with phenotypes. The distribution of variants within families or populations can facilitate identification of the genetic factors of disease. Illumina''s panel of genotyping BeadChips allows investigators to genotype thousands or millions of single nucleotide polymorphisms (SNPs) or to analyze other genomic variants, such as copy number, across a large number of DNA samples. These SNPs can be spread throughout the genome or targeted in specific regions in order to maximize potential discovery. The Infinium assay has been optimized to yield high-quality, accurate results quickly. With proper setup, a single technician can process from a few hundred to over a thousand DNA samples per week, depending on the type of array. This assay guides users through every step, starting with genomic DNA and ending with the scanning of the array. Using propriety reagents, samples are amplified, fragmented, precipitated, resuspended, hybridized to the chip, extended by a single base, stained, and scanned on either an iScan or Hi Scan high-resolution optical imaging system. One overnight step is required to amplify the DNA. The DNA is denatured and isothermally amplified by whole-genome amplification; therefore, no PCR is required. Samples are hybridized to the arrays during a second overnight step. By the third day, the samples are ready to be scanned and analyzed. Amplified DNA may be stockpiled in large quantities, allowing bead arrays to be processed every day of the week, thereby maximizing throughput.     

Characterising chromosome rearrangements: recent technical advances in molecular cytogenetics

Genomic rearrangements can result in losses, amplifications, translocations and inversions of DNA fragments thereby modifying genome architecture, and potentially having clinical consequences. Many genomic disorders caused by structural variation have initially been uncovered by early cytogenetic methods. The last decade has seen significant progression in molecular cytogenetic techniques, allowing rapid and precise detection of structural rearrangements on a whole-genome scale. The high resolution attainable with these recently developed techniques has also uncovered the role of structural variants in normal genetic variation alongside single-nucleotide polymorphisms (SNPs). We describe how array-based comparative genomic hybridisation, SNP arrays, array painting and next-generation sequencing analytical methods (read depth, read pair and split read) allow the extensive characterisation of chromosome rearrangements in human genomes.     

Oligonucleotide array for detection of common severe determinants of alpha thalassemia

A simple and an efficient oligonucleotide array was developed to identify common severe determinants of alpha (alpha) thalassemia. A total of 14 probes were designed to detect the most frequently three deletions (-alpha(3.7), -alpha(4.2), -(SEA)) and two non-deletions (alpha(Quong Sze), alpha(Constant Spring)). PCR products were amplified from human genomic DNA and allowed to hybridize with the oligonucleotide array. Hybridization was detected by fluorescence scanning, and alpha globin genotypes were assigned by quantitative analysis of the hybridization results. The efficiency and specificity of identifying alpha globin genotypes using the oligonucleotide arrays was evaluated by blinded analysis of 690 samples from unrelated individuals. The oligonucleotide array method described in this paper provides unambiguous detection of complex combinations of heterozygous, compound heterozygous and homozygous alpha thalassemia genotypes. The experimental results demonstrate that this methodological approach may be applied for screening and for hemological diagnosis in population at large.     

Rapid detection of gene mutations responsible for non-syndromic aortic aneurysm and dissection using two different methods: resequencing microarray technology and next-generation sequencing

Aortic aneurysm and/or dissection (AAD) is a life-threatening condition, and several syndromes are known to be related to AAD. In this study, two new technologies, resequencing array technology (ResAT) and next-generation sequencing (NGS), were used to analyze eight genes associated with syndromic AAD in 70 patients with non-syndromic AAD. Eighteen sequence variants were detected using both ResAT and NGS. In addition one of these sequence variants was detected by ResAT only and two additional variants by NGS only. Three of the 18 variants are likely to be pathogenic (in 4.3% of AAD patients and in 8.6% of a subset of patients with thoracic AAD), highlighting the importance of genetic analysis in non-syndromic AAD. ResAT and NGS similarly detected most, but not all, of the variants. Resequencing array technology was a rapid and efficient method for detecting most nucleotide substitutions, but was unable to detect short insertions/deletions, and it is impractical to update custom arrays frequently. Next-generation sequencing was able to detect almost all types of mutation, but requires improved informatics methods.     

Comparative genomic hybridization arrays in clinical pathology: progress and challenges

Array-based comparative genomic hybridization (array CGH) genome scanning is a powerful method for the global detection of gains and losses of genetic material in both congenital and neoplastic disorders. When used as a clinical diagnostic test, array CGH combines the whole genome perspective of traditional G-banded cytogenetics with the targeted identification of cryptic chromosomal abnormalities characteristic of fluorescence in situ hybridization (FISH). However, the presence of structural variants in the human genome can complicate analysis of patient samples, and array CGH does not provide morphologic information about chromosome structure, balanced translocations, or the actual chromosomal location of segmental duplications. Identification of such anomalies has significant diagnostic and prognostic implications for the patient. We therefore propose that array CGH should be used as a guide to the presence of genomic structural rearrangements in germline and tumor genomes that can then be further characterized by FISH or G-banding, depending on the clinical scenario. In this article, we share some of our experiences with diagnostic array CGH and discuss recent progress and challenges involved with the integration of array CGH into clinical laboratory medicine.     

Regeneration of comparative genomic hybridization oligonucleotide microarrays with dimethylurea

Reuse of materials in DNA hybridization-based methods has been known since the advent of Southern membranes. Array-based comparative genomic hybridization is essentially Southern hybridization with multiple probes immobilized on a solid surface. We show that comparative genomic hybridization microarrays fabricated with maskless array synthesizer technology can be used up to four times with the application of 1,3-dimethylurea as an array-stripping agent. We reproducibly detected chromosomal aberrations (0.6-22.4Mb in size) in four hybridization rounds using regenerated microarray slides. We also demonstrated that regenerated arrays can detect smaller alterations (16-200kbp), such as common copy number variants, as well as complex aberration profiles in tumor DNA.     

Array Comparative Genomic Hybridization (Array CGH) for Detection of Genomic Copy Number Variants

Array CGH for the detection of genomic copy number variants has replaced G-banded karyotype analysis. This paper describes the technology and its application in a clinical diagnostic service laboratory. DNA extracted from a patient’s sample (blood, saliva or other tissue types) is labeled with a fluorochrome (either cyanine 5 or cyanine 3). A reference DNA sample is labeled with the opposite fluorochrome. There follows a cleanup step to remove unincorporated nucleotides before the labeled DNAs are mixed and resuspended in a hybridization buffer and applied to an array comprising ~60,000 oligonucleotide probes from loci across the genome, with high probe density in clinically important areas. Following hybridization, the arrays are washed, then scanned and the resulting images are analyzed to measure the red and green fluorescence for each probe. Software is used to assess the quality of each probe measurement, calculate the ratio of red to green fluorescence and detect potential copy number variants.     

Suspension array technology: evolution of the flat-array paradigm.

Suspension arrays of microspheres analyzed using flow cytometry offer a new approach to multiplexed assays for large-scale screening applications. By optically encoding micron-sized polymer particles, suspension microarrays can be created to enable highly multiplexed analysis of complex samples. Each element in the array is comprised of a subpopulation of particles with distinct optical properties and each array element bears a different surface receptor. Nucleic acids, proteins, lipids or carbohydrates can serve as receptors to support the analysis of a wide range of biomolecular assemblies, and applications in genomic and proteomic research are being developed. Coupled with recent innovations for rapid serial analysis of samples, molecular analysis with microsphere arrays holds significant potential as a general analysis platform for both research and clinical applications.     

Copy number imbalances detected with a BAC-based array comparative genomic hybridization platform in congenital diaphragmatic hernia fetuses

Congenital diaphragmatic hernia (CDH) is a phenotypically and genetically heterogeneous disorder, with a complex inheritance pattern. Structural abnormalities of almost all chromosomes have been described in association with CDH. We made a molecular analysis through array comparative genomic hybridization (array CGH) of a group of fetuses with prenatal ultrasound diagnosis of CDH and normal G-banded karyotypes. A whole genome BAC-array CGH, composed of approximately 5000 BAC clones, was carried out on blood samples from fetuses with prenatal ultrasound diagnosis of CDH and a normal karyotype (500-band level). All potential cytogenetic alterations detected on the arrays were reported. The array CGH analysis showed copy number gains and losses in 10 of 12 cases. Eighty-five clones showed genomic imbalances, and 29 clones displayed described copy number variations. We identified a recurrent gain in 17q12 in two of 12 cases, which has not been previously described. Our results may contribute to determining the effectiveness and applicability of array CGH for prenatal diagnosis purposes, and also to elucidate the submicroscopic genomic instability of CDH fetuses.     

Genome-wide mapping of copy number variation in humans: comparative analysis of high resolution array platforms

Accurate and efficient genome-wide detection of copy number variants (CNVs) is essential for understanding human genomic variation, genome-wide CNV association type studies, cytogenetics research and diagnostics, and independent validation of CNVs identified from sequencing based technologies. Numerous, array-based platforms for CNV detection exist utilizing array Comparative Genome Hybridization (aCGH), Single Nucleotide Polymorphism (SNP) genotyping or both. We have quantitatively assessed the abilities of twelve leading genome-wide CNV detection platforms to accurately detect Gold Standard sets of CNVs in the genome of HapMap CEU sample NA12878, and found significant differences in performance. The technologies analyzed were the NimbleGen 4.2 M, 2.1 M and 3×720 K Whole Genome and CNV focused arrays, the Agilent 1×1 M CGH and High Resolution and 2×400 K CNV and SNP+CGH arrays, the Illumina Human Omni1Quad array and the Affymetrix SNP 6.0 array. The Gold Standards used were a 1000 Genomes Project sequencing-based set of 3997 validated CNVs and an ultra high-resolution aCGH-based set of 756 validated CNVs. We found that sensitivity, total number, size range and breakpoint resolution of CNV calls were highest for CNV focused arrays. Our results are important for cost effective CNV detection and validation for both basic and clinical applications.     

Molekulare Karyotypisierung in der genetischen Diagnostik

Molecular karyotyping by array comparative genomic hybridization (CGH) and single nucleotide polymorphism (SNP) arrays allows for a high resolution scan of the entire genome. It detects gains or losses (copy number variants) that might be the underlying cause of a genetic disorder. This technique is mainly applied to cases with syndromal or non-syndromal impaired (intellectual) development and is used to characterize genetic aberrations of tumor samples. Furthermore, molecular karyotyping might be useful for resolving prenatal cases with abnormal ultrasound findings. The purpose of this article is to explain the basic techniques, their limitations and strengths and to provide an outlook on future prospects.     

Frozen protein arrays: a new method for arraying and detecting recombinant and native tissue proteins

DNA microarrays are powerful tools for high throughput analysis of gene expression; however, they do not measure protein expression. Current methods for producing protein arrays require sophisticated equipment or extensive protein modification. We developed a low overhead, customizable assay platform called frozen protein arrays that can detect native proteins in protein lysates. Frozen protein arrays were formed from a block of frozen histologic embedding compound containing an array of wells. The wells were filled with samples, which freeze and bond to the block. Cryosections were cut and transferred to nitrocellulose-coated slides. The reproducibility, linearity, and sensitivity was confirmed using frozen protein arrays filled with prostate specific antigen. Frozen protein arrays could detect native tissue proteins. The alpha1 subunit of NaK-ATPase was detected in rat kidneys with a coefficient of variation of 4.3-6.6%. Frozen protein array analysis indicated that the protein abundance decreased by 48.7% following renal ischemia, similar to the 40% decrease by Western blotting. We conclude that frozen protein arrays are a low cost, moderate size platform for arraying samples including protein lysates. Production of many identical frozen protein arrays is easy, inexpensive, and requires only small sample volumes. The method is gentle on proteins as they remain frozen during production.     

Impact of whole genome amplification on analysis of copy number variants

Large-scale copy number variants (CNVs) have recently been recognized to play a role in human genome variation and disease. Approaches for analysis of CNVs in small samples such as microdissected tissues can be confounded by limited amounts of material. To facilitate analyses of such samples, whole genome amplification (WGA) techniques were developed. In this study, we explored the impact of Phi29 multiple-strand displacement amplification on detection of CNVs using oligonucleotide arrays. We extracted DNA from fresh frozen lymph node samples and used this for amplification and analysis on the Affymetrix Mapping 500k SNP array platform. We demonstrated that the WGA procedure introduces hundreds of potentially confounding CNV artifacts that can obscure detection of bona fide variants. Our analysis indicates that many artifacts are reproducible, and may correlate with proximity to chromosome ends and GC content. Pair-wise comparison of amplified products considerably reduced the number of apparent artifacts and partially restored the ability to detect real CNVs. Our results suggest WGA material may be appropriate for copy number analysis when amplified samples are compared to similarly amplified samples and that only the CNVs with the greatest significance values detected by such comparisons are likely to be representative of the unamplified samples.     

Linkage disequilibrium among commonly genotyped SNP variants detected from bull sequence,

Genomic prediction utilizing causal variants could increase selection accuracy above that achieved with SNPs genotyped by currently available arrays used for genomic selection. A number of variants detected from sequencing influential sires are likely to be causal, but noticeable improvements in prediction accuracy using imputed sequence variant genotypes have not been reported. Improvement in accuracy of predicted breeding values may be limited by the accuracy of imputed sequence variants. Using genotypes of SNPs on a high‐density array and non‐synonymous SNPs detected in sequence from influential sires of a multibreed population, results of this examination suggest that linkage disequilibrium between non‐synonymous and array SNPs may be insufficient for accurate imputation from the array to sequence. In contrast to 75% of array SNPs being strongly correlated to another SNP on the array, less than 25% of the non‐synonymous SNPs were strongly correlated to an array SNP. When correlations between non‐synonymous and array SNPs were strong, distances between the SNPs were greater than separation that might be expected based on linkage disequilibrium decay. Consistently near‐perfect whole‐genome linkage disequilibrium between the full array and each non‐synonymous SNP within the sequenced bulls suggests that whole‐genome approaches to infer sequence variants might be more accurate than imputation based on local haplotypes. Opportunity for strong linkage disequilibrium between sequence and array SNPs may be limited by discrepancies in allele frequency distributions, so investigating alternate genotyping approaches and panels providing greater chances of frequency‐matched SNPs strongly correlated to sequence variants is also warranted. Genotypes used for this study are available from https://www.animalgenome.org/repository/pub/ ;USDA2017.0519/.     

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.     

Using expression arrays for copy number detection: an example from E. coli

Background  

The sequencing of many genomes and tiling arrays consisting of millions of DNA segments spanning entire genomes have made high-resolution copy number analysis possible. Microarray-based comparative genomic hybridization (array CGH) has enabled the high-resolution detection of DNA copy number aberrations. While many of the methods and algorithms developed for the analysis microarrays have focused on expression analysis, the same technology can be used to detect genetic alterations, using for example standard commercial Affymetrix arrays. Due to the nature of the resultant data, standard techniques for processing GeneChip expression experiments are inapplicable.     

Efficient oligonucleotide probe selection for pan-genomic tiling arrays

Background  

Array comparative genomic hybridization is a fast and cost-effective method for detecting, genotyping, and comparing the genomic sequence of unknown bacterial isolates. This method, as with all microarray applications, requires adequate coverage of probes targeting the regions of interest. An unbiased tiling of probes across the entire length of the genome is the most flexible design approach. However, such a whole-genome tiling requires that the genome sequence is known in advance. For the accurate analysis of uncharacterized bacteria, an array must query a fully representative set of sequences from the species' pan-genome. Prior microarrays have included only a single strain per array or the conserved sequences of gene families. These arrays omit potentially important genes and sequence variants from the pan-genome.