High fidelity of whole-genome amplified DNA on high-density single nucleotide polymorphism arrays

Successful use of whole genome amplified DNA from multiple source types for high-density Illumina SNP microarrays

BACKGROUND

The recommended genomic DNA input requirements for whole genome single nucleotide polymorphism microarrays can limit the scope of molecular epidemiological studies. We performed a large-scale evaluation of whole genome amplified DNA as input into high-density, whole-genome Illumina® Infinium® SNP microarray.

RESULTS

Overall, 6622 DNA samples from 5970 individuals were obtained from three distinct biospecimen sources and genotyped using gDNA and/or wgaDNA inputs. When genotypes from the same individual were compared with standard, native gDNA input amount, we observed 99.94% mean concordance with wgaDNA input.

CONCLUSIONS

Our results demonstrate that carefully conducted studies with wgaDNA inputs can yield high-quality genotyping results. These findings should enable investigators to consider expansion of ongoing studies using high-density SNP microarrays, currently challenged by small amounts of available DNA.

Variant discovery in targeted resequencing using whole genome amplified DNA

Background

Next generation sequencing and advances in genomic enrichment technologies have enabled the discovery of the full spectrum of variants from common to rare alleles in the human population. The application of such technologies can be limited by the amount of DNA available. Whole genome amplification (WGA) can overcome such limitations. Here we investigate applicability of using WGA by comparing SNP and INDEL variant calls from a single genomic/WGA sample pair from two capture separate experiments: a 50 Mbp whole exome capture and a custom capture array of 4 Mbp region on chr12.

Results

Our results comparing variant calls derived from genomic and WGA DNA show that the majority of variant SNP and INDEL calls are common to both callsets, both at the site and genotype level and suggest that allele bias plays a minimal role when using WGA DNA in re-sequencing studies.

Conclusions

Although the results of this study are based on a limited sample size, they suggest that using WGA DNA allows the discovery of the vast majority of variants, and achieves high concordance metrics, when comparing to genomic DNA calls.

Inter-chromosomal level of genome organization and longevity-related phenotypes in humans

Studies focusing on unraveling the genetic origin of health span in humans assume that polygenic, aging-related phenotypes are inherited through Mendelian mechanisms of inheritance of individual genes. We use the Framingham Heart Study (FHS) data to examine whether non-Mendelian mechanisms of inheritance can drive linkage of loci on non-homologous chromosomes and whether such mechanisms can be relevant to longevity-related phenotypes. We report on genome-wide inter-chromosomal linkage disequilibrium (LD) and on chromosome-wide intra-chromosomal LD and show that these are real phenomena in the FHS data. Genetic analysis of inheritance in families based on Mendelian segregation reveals that the alleles of single nucleotide polymorphisms (SNPs) in LD at loci on non-homologous chromosomes are inherited as a complex resembling haplotypes of a genetic unit. This result implies that the inter-chromosomal LD is likely caused by non-random assortment of non-homologous chromosomes during meiosis. The risk allele haplotypes can be subject to dominant-negative selection primary through the mechanisms of non-Mendelian inheritance. They can go to extinction within two human generations. The set of SNPs in inter-chromosomal LD (N=68) is nearly threefold enriched, with high significance (p=1.6 × 10(-9)), on non-synonymous coding variants (N=28) compared to the entire qualified set of the studied SNPs. Genes for the tightly linked SNPs are involved in fundamental biological processes in an organism. Survival analyses show that the revealed non-genetic linkage is associated with heritable complex phenotype of premature death. Our results suggest the presence of inter-chromosomal level of functional organization in the human genome and highlight a challenging problem of genomics of human health and aging.

Sortilin receptor 1 predicts longitudinal cognitive change

The gene encoding sortilin receptor 1 (SORL1) has been associated with Alzheimer's disease risk. We examined 15 SORL1 variants and single nucleotide polymorphism (SNP) set risk scores in relation to longitudinal verbal, spatial, memory, and perceptual speed performance, testing for age trends and sex-specific effects. Altogether, 1609 individuals from 3 population-based Swedish twin studies were assessed up to 5 times across 16 years. Controlling for apolipoprotein E genotype (APOE), multiple simple and sex-moderated associations were observed for spatial, episodic memory, and verbal trajectories (p = 1.25E-03 to p = 4.83E-02). Five variants (rs11600875, rs753780, rs7105365, rs11820794, rs2070045) were associated across domains. Notably, in those homozygous for the rs2070045 risk allele, men demonstrated initially favorable performance but accelerating declines, and women showed overall lower performance. SNP set risk scores predicted spatial (Card Rotations, p = 5.92E-03) and episodic memory trajectories (Thurstone Picture Memory, p = 3.34E-02), where higher risk scores benefited men's versus women's performance up to age 75 but with accelerating declines. SORL1 is associated with cognitive aging, and might contribute differentially to change in men and women.

Quality requirements for genomic DNA preparations and storage conditions for a high-density oligonucleotide microarray

High-density oligonucleotide microarrays are widely used in genome-wide association studies. The purpose of this study was to assess the influence of various factors during the preparation of DNA on genotype calling for the Affymetrix high-density oligonucleotide microarray 250K GeneChip. DNA was extracted from peripheral whole blood by solution-based and silica-membrane-based methods. Blood was stored at 4°C or 25°C for 4 or 24 h, followed by DNA extraction. To examine the effects of freeze-thaw cycles, blood and DNA were also subjected to 5 and 10 or 20 of freeze-thaw cycles, respectively. The suitability of variously DNA preparations for the array was assessed by the call rate resulting from genotyping. All DNA samples showed mean call rates of more than 0.99, which passed the quality criteria for genotyping (greater than 0.95). The results indicated that the solution-based method and the silica-membrane-based DNA extraction method could provide DNA of sufficient quality for genotyping. In addition, DNA quality suitable for high-density oligonucleotide microarrays is not strongly dependent on the preparation conditions under standard procedures.

Characterization of whole genome amplified (WGA) DNA for use in genotyping assay development

Background

Genotyping assays often require substantial amounts of DNA. To overcome the problem of limiting amounts of available DNA, Whole Genome Amplification (WGA) methods have been developed. The multiple displacement amplification (MDA) method using Φ29 polymerase has become the preferred choice due to its high processivity and low error rate. However, the uniformity and fidelity of the amplification process across the genome has not been extensively characterized.

Results

To assess amplification uniformity, we used array-based comparative genomic hybridization (aCGH) to evaluate DNA copy number variations (CNVs) in DNAs amplified by two MDA kits: GenomiPhi and REPLI-g. The Agilent Human CGH array containing nearly one million probes was used in this study together with DNAs from a normal subject and 2 cystic fibrosis (CF) patients. Each DNA sample was amplified 4 independent times and compared to its native unamplified DNA. Komogorov distances and Phi correlations showed a high consistency within each sample group. Less than 2% of the probes showed more than 2-fold CNV introduced by the amplification process. The two amplification kits, REPLI-g and GenomiPhi, generate very similar amplified DNA samples despite the differences between the unamplified and amplified DNA samples. The results from aCGH analysis indicated that there were no obvious CNVs in the CFTR gene region due to WGA when compared to unamplified DNA. This was confirmed by quantitative real-time PCR copy number assays at 10 locations within the CFTR gene. DNA sequencing analysis of a 2-kb region within the CFTR gene showed no mutations introduced by WGA.

Conclusion

The relatively high uniformity and consistency of the WGA process, coupled with the low replication error rate, suggests that WGA DNA may be suitable for accurate genotyping. Regions of the genome that were consistently under-amplified were found to contain higher than average GC content. Because of the consistent differences between the WGA DNA and the native unamplified DNA, characterization of the genomic region of interest, as described here, will be necessary to ensure the reliability of genotyping results from WGA DNA.

Copy number and loss of heterozygosity detected by SNP array of formalin-fixed tissues using whole-genome amplification

The requirement for large amounts of good quality DNA for whole-genome applications prohibits their use for small, laser capture micro-dissected (LCM), and/or rare clinical samples, which are also often formalin-fixed and paraffin-embedded (FFPE). Whole-genome amplification of DNA from these samples could, potentially, overcome these limitations. However, little is known about the artefacts introduced by amplification of FFPE-derived DNA with regard to genotyping, and subsequent copy number and loss of heterozygosity (LOH) analyses. Using a ligation adaptor amplification method, we present data from a total of 22 Affymetrix SNP 6.0 experiments, using matched paired amplified and non-amplified DNA from 10 LCM FFPE normal and dysplastic oral epithelial tissues, and an internal method control. An average of 76.5% of SNPs were called in both matched amplified and non-amplified DNA samples, and concordance was a promising 82.4%. Paired analysis for copy number, LOH, and both combined, showed that copy number changes were reduced in amplified DNA, but were 99.5% concordant when detected, amplifications were the changes most likely to be 'missed', only 30% of non-amplified LOH changes were identified in amplified pairs, and when copy number and LOH are combined ∼50% of gene changes detected in the unamplified DNA were also detected in the amplified DNA and within these changes, 86.5% were concordant for both copy number and LOH status. However, there are also changes introduced as ∼20% of changes in the amplified DNA are not detected in the non-amplified DNA. An integrative network biology approach revealed that changes in amplified DNA of dysplastic oral epithelium localize to topologically critical regions of the human protein-protein interaction network, suggesting their functional implication in the pathobiology of this disease. Taken together, our results support the use of amplification of FFPE-derived DNA, provided sufficient samples are used to increase power and compensate for increased error rates.

The microbial detection array combined with random Phi29-amplification used as a diagnostic tool for virus detection in clinical samples

A common technique used for sensitive and specific diagnostic virus detection in clinical samples is PCR that can identify one or several viruses in one assay. However, a diagnostic microarray containing probes for all human pathogens could replace hundreds of individual PCR-reactions and remove the need for a clear clinical hypothesis regarding a suspected pathogen. We have established such a diagnostic platform for random amplification and subsequent microarray identification of viral pathogens in clinical samples. We show that Phi29 polymerase-amplification of a diverse set of clinical samples generates enough viral material for successful identification by the Microbial Detection Array, demonstrating the potential of the microarray technique for broad-spectrum pathogen detection. We conclude that this method detects both DNA and RNA virus, present in the same sample, as well as differentiates between different virus subtypes. We propose this assay for diagnostic analysis of viruses in clinical samples.

STAT1 mutations in autosomal dominant chronic mucocutaneous candidiasis

BACKGROUND

Chronic mucocutaneous candidiasis (CMC) is characterized by susceptibility to candida infection of skin, nails, and mucous membranes. Patients with recessive CMC and autoimmunity have mutations in the autoimmune regulator AIRE. The cause of autosomal dominant CMC is unknown.

METHODS

We evaluated 14 patients from five families with autosomal dominant CMC. We incubated their peripheral-blood mononuclear cells with different combinations of stimuli to test the integrity of pathways that mediate immunity, which led to the selection of 100 genes that were most likely to contain the genetic defect. We used an array-based sequence-capture assay, followed by next-generation sequencing, to identify mutations.

RESULTS

The mononuclear cells from the affected patients were characterized by poor production of interferon-γ, interleukin-17, and interleukin-22, suggesting that the defect lay within the interleukin-12 receptor and interleukin-23 receptor signaling pathways. We identified heterozygous missense mutations in the DNA sequence encoding the coiled-coil (CC) domain of signal transducer and activator of transcription 1 (STAT1) in the patients. These mutations lead to defective responses in type 1 and type 17 helper T cells (Th1 and Th17). The interferon-γ receptor pathway was intact in these patients.

CONCLUSIONS

Mutations in the CC domain of STAT1 underlie autosomal dominant CMC and lead to defective Th1 and Th17 responses, which may explain the increased susceptibility to fungal infection. (Funded by the Netherlands Organization for Scientific Research and others.).

Laser capture microdissection (LCM) and whole genome amplification (WGA) of DNA from normal breast tissue --- optimization for genome wide array analyses

Background

Laser capture microdissection (LCM) can be applied to tissues where cells of interest are distinguishable from surrounding cell populations. Here, we have optimized LCM for fresh frozen normal breast tissue where large amounts of fat can cause problems during microdissection. Since the amount of DNA needed for genome wide analyses, such as single nucleotide polymorphism (SNP) arrays, is often greater than what can be obtained from the dissected tissue, we have compared three different whole genome amplification (WGA) kits for amplification of DNA from LCM material. In addition, the genome wide profiling methods commonly used today require extremely high DNA quality compared to PCR based techniques and DNA quality is thus critical for successful downstream analyses.

Findings

We found that by using FrameSlides without glass backing for LCM and treating the slides with acetone after staining, the problems caused by excessive fat could be significantly decreased. The amount of DNA obtained after extraction from LCM tissue was not sufficient for direct SNP array analysis in our material. However, the two WGA kits based on Phi29 polymerase technology (Repli-g^® (Qiagen) and GenomiPhi (GE Healthcare)) gave relatively long amplification products, and amplified DNA from Repli-g^® gave call rates in the subsequent SNP analysis close to those from non-amplified DNA. Furthermore, the quality of the input DNA for WGA was found to be essential for successful SNP array results and initial DNA fragmentation problems could be reduced by switching from a regular halogen lamp to a VIS-LED lamp during LCM.

Conclusions

LCM must be optimized to work satisfactorily in difficult tissues. We describe a work flow for fresh frozen normal breast tissue where fat is inclined to cause problems if sample treatment is not adapted to this tissue. We also show that the Phi29-based Repli-g^® WGA kit (Qiagen) is a feasible approach to amplify DNA of high quality prior to genome wide analyses such as SNP profiling.

Single nucleotide polymorphism microarray analysis of genetic alterations in cancer

The identification of structural genetic alterations, including DNA amplifications, deletions, and loss of heterozygosity (LOH), using single nucleotide polymorphism (SNP) microarrays has provided important insights into the pathogenesis of a number of hematologic malignancies. Currently available SNP arrays comprise over a million SNP and copy number oligonucleotide probes that interrogate the genome at sub-kilobase resolution. The accurate detection of DNA copy number abnormalities and LOH is critically dependent on the use of high-quality DNA, the use of matched reference samples wherever possible, optimal normalization of raw microarray data, and computational algorithms to detect copy number alterations sensitively and robustly. This chapter provides methods and guidelines for preparing samples, processing and analyzing data, and validation of novel lesions. Specific examples are provided for Affymetrix SNP arrays in acute lymphoblastic leukemia.

Technical note: High fidelity of whole-genome amplified sheep (Ovis aries) deoxyribonucleic acid using a high-density single nucleotide polymorphism array-based genotyping platform

Advances in high-throughput genotyping technologies have afforded researchers the opportunity to study ever-increasing numbers of SNP in animal genomes. However, many studies encounter difficulties in obtaining sufficient quantities of high-quality DNA for such analyses, particularly when the source biological material is limited or degraded. The recent development of in vitro whole-genome amplification approaches has permitted researchers to circumvent these challenges by increasing the amount of usable DNA in normally small-quantity samples. Here, we assess the performance of whole-genome amplification products generated from ovine genomic DNA using a high-throughput SNP genotyping platform, the newly developed Illumina ovineSNP50 BeadChip. Our results demonstrate a high genotype call rate for conventional genomic DNA and whole-genome amplified genomic DNA. The data also reveal an exceptionally high concordance rate ( > or = 99%) between the genotypes generated from whole-genome amplified products and their conventional genomic DNA counterparts. This study supports the use of whole-genome amplification as a viable solution for the analysis of high-density SNP genotypic data using compromised or limited starting material.

Novel generating protective single nucleotide polymorphism barcode for breast cancer using particle swarm optimization

BACKGROUND

High-throughput single nucleotide polymorphism (SNP) genotyping generates a huge amount of SNP data in genome-wide association studies. Simultaneous analyses for multiple SNP interactions associated with many diseases and cancers are essential; however, these analyses are still computationally challenging.

METHODS

In this study, we propose an odds ratio-based binary particle swarm optimization (OR-BPSO) method to evaluate the risk of breast cancer.

RESULTS

BPSO provides the combinational SNPs with their corresponding genotype, called SNP barcodes, with the maximal difference of occurrence between the control and breast cancer groups. A specific SNP barcode with an optimized fitness value was identified among seven SNP combinations within the space of one minute. The identified SNP barcodes with the best performance between control and breast cancer groups were found to be control-dominant, suggesting that these SNP barcodes may prove protective against breast cancer. After statistical analysis, these control-dominant SNP barcodes were processed for odds ratio analysis for quantitative measurement with regard to the risk of breast cancer.

CONCLUSION

This study proposes an effective high-speed method to analyze the SNP-SNP interactions for breast cancer association study.

Fine-scaled human genetic structure revealed by SNP microarrays

We report an analysis of more than 240,000 loci genotyped using the Affymetrix SNP microarray in 554 individuals from 27 worldwide populations in Africa, Asia, and Europe. To provide a more extensive and complete sampling of human genetic variation, we have included caste and tribal samples from two states in South India, Daghestanis from eastern Europe, and the Iban from Malaysia. Consistent with observations made by Charles Darwin, our results highlight shared variation among human populations and demonstrate that much genetic variation is geographically continuous. At the same time, principal components analyses reveal discernible genetic differentiation among almost all identified populations in our sample, and in most cases, individuals can be clearly assigned to defined populations on the basis of SNP genotypes. All individuals are accurately classified into continental groups using a model-based clustering algorithm, but between closely related populations, genetic and self-classifications conflict for some individuals. The 250K data permitted high-level resolution of genetic variation among Indian caste and tribal populations and between highland and lowland Daghestani populations. In particular, upper-caste individuals from Tamil Nadu and Andhra Pradesh form one defined group, lower-caste individuals from these two states form another, and the tribal Irula samples form a third. Our results emphasize the correlation of genetic and geographic distances and highlight other elements, including social factors that have contributed to population structure.