Affymetrix Genome-wide Single Nucleotide Polymorphism Research Articles

A map of copy number variations in the Tunisian population: a valuable tool for medical genomics in North Africa

Copy number variation (CNV) is considered as the most frequent type of structural variation in the human genome. Some CNVs can act on human phenotype diversity, encompassing rare Mendelian diseases and genomic disorders. The North African populations remain underrepresented in public genetic databases in terms of single-nucleotide variants as well as for larger genomic mutations. In this study, we present the first CNV map for a North African population using the Affymetrix Genome-Wide SNP (single-nucleotide polymorphism) array 6.0 array genotyping intensity data to call CNVs in 102 Tunisian healthy individuals. Two softwares, PennCNV and Birdsuite, were used to call CNVs in order to provide reliable data. Subsequent bioinformatic analyses were performed to explore their features and patterns. The CNV map of the Tunisian population includes 1083 CNVs spanning 61.443 Mb of the genome. The CNV length ranged from 1.017 kb to 2.074 Mb with an average of 56.734 kb. Deletions represent 57.43% of the identified CNVs, while duplications and the mixed loci are less represented. One hundred and three genes disrupted by CNVs are reported to cause 155 Mendelian diseases/phenotypes. Drug response genes were also reported to be affected by CNVs. Data on genes overlapped by deletions and duplications segments and the sequence properties in and around them also provided insights into the functional and health impacts of CNVs. These findings represent valuable clues to genetic diversity and personalized medicine in the Tunisian population as well as in the ethnically similar populations from North Africa.

Genome-Wide Association Meta-Analysis Reveals Novel Juvenile Idiopathic Arthritis Susceptibility Loci.

Juvenile idiopathic arthritis (JIA) is the most common childhood rheumatic disease and has a strong genomic component. To date, JIA genetic association studies have had limited sample sizes, usedheterogeneous patient populations, or included only candidate regions. The aim of this study was to identify new associations between JIA patients with oligoarticular disease and those with IgM rheumatoid factor (RF)-negative polyarticular disease, which are clinically similar and the most prevalent JIA disease subtypes. Three cohorts comprising 2,751 patients with oligoarticular or RF-negative polyarticular JIA were genotyped using the Affymetrix Genome-Wide SNP Array 6.0 or the Illumina HumanCoreExome-12+ Array. Overall, 15,886 local and out-of-study controls, typed on these platforms or the Illumina HumanOmni2.5, were used forassociation analyses. High-quality single-nucleotide polymorphisms (SNPs) were used for imputation to 1000 Genomes prior to SNP association analysis. Meta-analysis showed evidence of association (P < 1 × 10-6 ) at 9 regions: PRR9_LOR (P = 5.12 × 10-8 ), ILDR1_CD86 (P = 6.73 × 10-8 ), WDFY4 (P = 1.79 × 10-7 ), PTH1R (P = 1.87 × 10-7 ), RNF215 (P = 3.09 × 10-7 ), AHI1_LINC00271 (P = 3.48 × 10-7 ), JAK1 (P = 4.18 × 10-7 ), LINC00951 (P = 5.80 × 10-7 ), and HBP1 (P = 7.29 × 10-7 ). Of these, PRR9_LOR, ILDR1_CD86, RNF215, LINC00951, and HBP1 were shown, for the first time, to be autoimmune disease susceptibility loci. Furthermore, associated SNPs included cis expression quantitative trait loci for WDFY4, CCDC12, MTP18, SF3A1, AHI1, COG5, HBP1, and GPR22. This study provides evidence of both unique JIA risk loci and risk loci overlapping between JIA and other autoimmune diseases. These newly associated SNPs are shown to influence gene expression, and their bounding regions tie into molecular pathways of immunologic relevance. Thus, they likely represent regions that contribute to the pathology of oligoarticular JIA and RF-negative polyarticular JIA.

Genome-wide association study of sporadic brain arteriovenous malformations

BackgroundThe pathogenesis of sporadic brain arteriovenous malformations (BAVMs) remains unknown, but studies suggest a genetic component. We estimated the heritability of sporadic BAVM and performed a genome-wide association study (GWAS)...

Common Variation at 6q25.3 (TULP4) Influences Risk for Arterial Thrombosis in Myeloproliferative Neoplasms

Inherent tendency for thrombosis is a major complication in myeloproliferative neoplasms (MPN). The molecular basis of thrombosis in MPN is not well understood, however, genetic factors have been proposed to play a role. To agnostically investigate the role of common germline variation in MPN thrombophilia, we performed genome-wide association studies (GWAS) in MPN patient cohorts characterized for arterial thrombosis (AT) events. In our discovery cohort (n=383) from Vienna, Austria, 18% of patients have suffered from AT after MPN diagnosis, and 33% of patients showed records for AT events at any time. Of this discovery cohort, a subset of patients (n=302) selected independently from thrombotic status was genotyped on the Affymetrix Genome-Wide SNP 6.0 array platform. After assessment of case-control setup and genotyping quality as standardly implemented in GWAS, confidently genotyped single nucleotide polymorphisms (SNPs) were tested for allelic association with occurrence of AT after MPN diagnosis. We observed an association signal beyond genome-wide statistical significance (P <5x10-8) at chromosome 6q25.3, tagged by rs6455579 and six additional closely related tag-SNPs. The 6q25.3 risk haplotype is common with a minor allele frequency (MAF) of ~8% in the Vienna MPN cohort. The minor allele rs6455579_C correlated with increased risk for AT. We next genotyped the full discovery cohort from Vienna (n=383) and a replication cohort from Pavia, Italy (n=505) using qPCR-based SNP genotyping for rs6455579 and computed odds ratios (ORs) to estimate the effect sizes of the association. In the Vienna cohort, rs6455579 genotype was significantly associated with both "AT after diagnosis" (P =1.90x10-8, ORhet=6.93) and "AT at any time" (P =6.96x10-4, ORhet=3.07). In the Pavia cohort, an association with "AT after diagnosis" could not be observed (P =0.14), however, the association of rs6455579 with "AT at any time" could be reproduced at formal statistical significance (P =7.92x10-3, ORhet=2.55). Moreover, Kaplan-Meier statistics on the Pavia cohort revealed a significant difference in AT-free survival after diagnosis upon rs6455579 genotype (log-rank-test P =0.009; rs6455579 homozygous major vs. heterozygous), underpinning the relevance of disease duration for the impact of the haplotype on MPN-related AT. To gain insight into the physiological mechanisms behind the 6q25.3 risk haplotype, we next evaluated the discovery cohort for correlations of the haplotype with clinical parameters other than AT. We could observe a significant trend for increased white blood cell (WBC) count in CALR mutated (n=90, P =2.50x10-3) but not in JAK2 mutated (n=182, P =0.52) ET and PMF patients carrying the risk allele. CALR mutated patients have been previously reported to exhibit significantly prolonged thrombosis-free survival as compared to JAK2 mutated patients. Thus, our observation indicates that the 6q25.3 germline risk haplotype might impact on WBC count most strongly in a subgroup of patients considered low risk for thrombosis based on the somatic mutational status. The 6q25.3 core haplotype (R2 >0.3) spans ~300 kilobases, covering intergenic sequence as well as promoter and 5' exons of the TULP4 gene. To test for the possibility of a rare coding variant in TULP4 or other more distant genes underlying the association through long-range linkage disequilibrium, we used genotype imputation in conjunction with the 1000 genomes reference panel to infer genotypes on all untyped variants (MAF>1%) in a 5 megabase region centered on the core haplotype. We did not detect any coding variants reflecting the association, and the tag-SNPs from the initial GWAS remained the most strongly associated variants. Causative non-coding genetic variation identified in GWAS is thought to exert its function through differential regulation of specific target genes. Therefore, we evaluated a potential influence of the 6q25.3 risk haplotype on TULP4 gene expression in peripheral blood (The Cancer Genome Atlas (TCGA) LAML dataset, RNA-Seq on 173 acute myeloid leukemia patients). Indeed we detected significantly decreased TULP4 expression in risk haplotype carriers (P =0.029), providing indirect evidence for reduced TULP4 transcript levels impacting on elevated risk for AT in MPN. Further studies will be required to functionally assess the potential role of TULP4 in MPN-related AT. DisclosuresGisslinger:Novartis: Honoraria, Research Funding, Speakers Bureau; Geron: Consultancy; Sanofi Aventis: Consultancy; Janssen Cilag: Honoraria, Speakers Bureau; Celgene: Consultancy, Honoraria, Research Funding, Speakers Bureau; AOP ORPHAN: Consultancy, Honoraria, Research Funding, Speakers Bureau. Kralovics:AOP Orphan: Research Funding; Qiagen: Membership on an entity's Board of Directors or advisory committees.

Insertion-deletions burden in copy number polymorphisms of the Tibetan population.

BACKGROUND:Many studies have been conducted to identify either insertions-deletions (inDels) or copy number variations (CNVs) in humans, but few studies have been conducted to identify both of these forms coexisting in the same region.AIMS AND OBJECTIVES:To map the functionally significant sites within human genes that are likely to influence human traits and diseases.MATERIALS AND METHODS:In this report, we describe an inDel map in the 1051 Tibetan CNV regions obtained through CNV genotyping using Affymetrix Genome-wide single nucleotide polymorphism 6.0 chip. InDel polymorphisms in these copy number polymorphism regions were identified with a computational approach using the 2500 deoxyribonucleic acid sequences obtained from the 1000 Genome Project.RESULTS:The study identified a total of 95935 inDels that range from 1 bp to several bps in length which were found scattered across regulatory regions, exons and in introns of genes underlying the CNVs. A study on the distribution of inDels revealed that the majority of inDels were found in coding regions of the genome than the noncoding, while within the genes, inDels in intron regions were more followed by exonic regions and finally the regulatory regions.CONCLUSION:Study of inDels in CNV regions contribute to the enhanced understanding of the role played by the two variations and their collective influence on the genome. Further, a collection of these inDel genetic markers will aid in genetic mapping, further understanding of the phenotypic variability, identification of disease genes and in detecting novel CNVs.

Genome-Wide Association Study (GWAS) Of Venous Thromboembolism (VTE) In African-Americans From The Electronic Medical Records & Genomics (eMERGE) Networkm

Background The incidence of VTE in African-Americans (AAs) is similar to or higher than in Americans of European ancestry. However, the carrier frequencies of inherited thrombophilias common in whites (i.e., Factor V Leiden, Prothrombin G20210A) are very low in AAs, suggesting that other inherited thrombophilias may be associated with VTE in AAs. Objective To identify potentially novel non-hypothesis-driven single nucleotide polymorphisms (SNPs) associated with VTE in AAs. Methods We used the resources of the eMERGE Network to perform a GWAS of VTE in AAs. The eMERGE Network (funded by NHGRI) is comprised of nine sites each with DNA biobanks linked to electronic health records (EHRs). Approximately 39,206 unique DNA samples have been genotyped using either Affymetrix or Illumina genome-wide SNP arrays. Led by the Coordinating Center and eMERGE genomics workgroup, an imputation pipeline was developed for merging genomic data across the different SNP arrays used by the eMERGE sites, to maximize sample size and the power to detect associations. The imputation was performed using the 1000 Genomes Cosmopolitan reference panel which includes 1092 individuals and over 36 million SNPs. From our previously-identified cohort of Olmsted County, MN residents with incident or recurrent VTE, 1996-2005, we derived and validated an EHR-driven phenotype extraction algorithm that leveraged structured data based on ICD-9-CM codes and unstructured data from clinical notes via natural language processing to identify VTE cases and controls with 100% and 94% positive and negative predictive values, respectively. We tested for an association between each SNP and VTE among AAs using unconditional logistic regression, adjusting for age, sex and eMERGE site. Results Among 294 AA VTE cases and 3,661 AA controls (total n=3,955; females, n=2,512), the Factor V Leiden (F5 rs6025) was not analyzed due to a very low minor allele frequency (MAF=0.0036). The prothrombin G20210A (F2 rs1799963) and ABO blood type O (ABO rs8176719) SNPs were not genotyped in any of the arrays and could not be imputed. Among SNPs with an imputation score &gt;0.8, the most significant SNPs associated with VTE were ITPR3 (inositol 1,4,5-triphosphate receptor type 3) rs2229637 (OR=1.65; p=3.61E-07; MAF=0.19) and CLEC7A (C-type lectin domain family 7, member A) rs59819090 (OR=2.16; p=1.06E-06; MAF=0.056). ITPR3 SNPs have been associated with coronary artery aneurysm in Kawasaki disease, type 1 diabetes mellitus and other autoimmune disorders. CLEC7A rs59819090 encodes for a serine82 to leucine nonsynonymous amino acid change in dectin-1. Dectin-1 is a transmembrane innate immune pattern recognition receptor on myeloid cells that upon binding it’s agonist, β-glucans, stimulates cytokine production and the respiratory burst, a prerequisite for formation of neutrophil extracelluar traps (NETs). NETs have been associated with VTE (van Montfoort, et al. Ateriosclero Thromb Vasc Biol 2013;33:147-51). Dectin-1 serine82 is not evolutionarily conserved and the serine82 to leucine amino acid change is predicted to be tolerated, with a moderate effect on protein function. Conclusions ITPR3 rs2229637 and CLAC7A rs59819090 may be associated with VTE in African-Americans. These observations require replication and functional studies toward understanding how they may lie on the causal pathway to VTE. Disclosures: No relevant conflicts of interest to declare.

Genome-wide association analysis of juvenile idiopathic arthritis identifies a new susceptibility locus at chromosomal region 3q13.

In a genome-wide association study of Caucasian patients with juvenile idiopathic arthritis (JIA), we have previously described findings limited to autoimmunity loci shared by JIA and other diseases. The present study was undertaken to identify novel JIA-predisposing loci using genome-wide approaches. The discovery cohort consisted of Caucasian JIA cases (n = 814) and local controls (n = 658) genotyped on the Affymetrix Genome-Wide SNP 6.0 Array, along with 2,400 out-of-study controls. In a replication study, we genotyped 10 single-nucleotide polymorphisms (SNPs) in 1,744 cases and 7,010 controls from the US and Europe. Analysis within the discovery cohort provided evidence of associations at 3q13 within C3orf1 and near CD80 (rs4688011) (odds ratio [OR] 1.37, P = 1.88 × 10(-6) ) and at 10q21 near JMJD1C (rs647989 [OR 1.59, P = 6.1 × 10(-8) ], rs12411988 [OR 1.57, P = 1.16 × 10(-7) ], and rs10995450 [OR 1.31, P = 6.74 × 10(-5) ]). Meta-analysis provided further evidence of association for these 4 SNPs (P = 3.6 × 10(-7) for rs4688011, P = 4.33 × 10(-5) for rs6479891, P = 2.71 × 10(-5) for rs12411988, and P = 5.39 × 10(-5) for rs10995450). Gene expression data on 68 JIA cases and 23 local controls showed cis expression quantitative trait locus associations for C3orf1 SNP rs4688011 (P = 0.024 or P = 0.034, depending on the probe set) and JMJD1C SNPs rs6479891 and rs12411988 (P = 0.01 or P = 0.04, depending on the probe set and P = 0.008, respectively). Using a variance component liability model, it was estimated that common SNP variation accounts for approximately one-third of JIA susceptibility. Genetic association results and correlated gene expression findings provide evidence of JIA association at 3q13 and suggest novel genes as plausible candidates in disease pathology.

Polymorphisms in the glial glutamate transporter SLC1A2 are associated with essential tremor

Sporadic, genetically complex essential tremor (ET) is one of the most common movement disorders and may lead to severe impairment of the quality of life. Despite high heritability, the genetic determinants of ET are largely unknown. We performed the second genome-wide association study (GWAS) for ET to elucidate genetic risk factors of ET. Using the Affymetrix Genome-Wide SNP Array 6.0 (1000K) we conducted a two-stage GWAS in a total of 990 subjects and 1,537 control subjects from Europe to identify genetic variants associated with ET. We discovered association of an intronic variant of the main glial glutamate transporter (SLC1A2) gene with ET in the first-stage sample (rs3794087, p = 6.95 × 10(-5), odds ratio [OR] = 1.46). We verified the association of rs3794087 with ET in a second-stage sample (p = 1.25 × 10(-3), OR = 1.38). In the subgroup analysis of patients classified as definite ET, rs3794087 obtained genome-wide significance (p = 3.44 × 10(-10), OR = 1.59) in the combined first- and second-stage sample. Genetic fine mapping using nonsynonymous single nucleotide polymorphisms (SNPs) and SNPs in high linkage disequilibrium with rs3794087 did not reveal any SNP with a stronger association with ET than rs3794087. We identified SLC1A2 encoding the major glial high-affinity glutamate reuptake transporter in the brain as a potential ET susceptibility gene. Acute and chronic glutamatergic overexcitation is implied in the pathogenesis of ET. SLC1A2 is therefore a good functional candidate gene for ET.

Effect of Combining Multiple CNV Defining Algorithms on the Reliability of CNV Calls from SNP Genotyping Data

In addition to single-nucleotide polymorphisms (SNP), copy number variation (CNV) is a major component of human genetic diversity. Among many whole-genome analysis platforms, SNP arrays have been commonly used for genomewide CNV discovery. Recently, a number of CNV defining algorithms from SNP genotyping data have been developed; however, due to the fundamental limitation of SNP genotyping data for the measurement of signal intensity, there are still concerns regarding the possibility of false discovery or low sensitivity for detecting CNVs. In this study, we aimed to verify the effect of combining multiple CNV calling algorithms and set up the most reliable pipeline for CNV calling with Affymetrix Genomewide SNP 5.0 data. For this purpose, we selected the 3 most commonly used algorithms for CNV segmentation from SNP genotyping data, PennCNV, QuantiSNP; and BirdSuite. After defining the CNV loci using the 3 different algorithms, we assessed how many of them overlapped with each other, and we also validated the CNVs by genomic quantitative PCR. Through this analysis, we proposed that for reliable CNV-based genomewide association study using SNP array data, CNV calls must be performed with at least 3 different algorithms and that the CNVs consistently called from more than 2 algorithms must be used for association analysis, because they are more reliable than the CNVs called from a single algorithm. Our result will be helpful to set up the CNV analysis protocols for Affymetrix Genomewide SNP 5.0 genotyping data.

DNA microarray as a tool in establishing genetic relatedness—Current status and future prospects

In the past decades, microarray technology has definitely put an edge to the field of genetic research. Our aim was to determine whether single nucleotide polymorphism (SNP) microarrays could be used as a tool in establishing genetic relationships where current molecular genetic methods are not sufficient. We used the Genechip, Affymetrix GenomeWide SNP Array 6.0, which detects more than 900,000 SNP markers dispersed throughout the human genome. The intention was to find a good selection of SNP markers that could be used for statistical evaluation of relatedness in a forensic setting. We conducted pairwise comparisons in the R-package FEST as well as pedigree comparisons in Merlin. Our methods were applied on two separate families, where relationships as distant as 3rd cousins were known. In addition, a question about a possible common ancestry between the two families was tested. Relationships as distant as 2nd cousins could be readily distinguished both from unrelated and other, genetically, closer relationships. This was achieved with a selection of 5774 markers, where each pair of markers was separated by a genetic distance of at least 0.5cM (centiMorgan). When considering 3rd cousins, and more distant relationships, the number of markers needs to be extended, consequently decreasing the genetic distance between the markers. However, inclusion of a too large number of markers presents new challenges and our results imply that the use of too dense sets of markers always yields the highest probability for the genetically closest relationship hypothesis. Simulations confirm that this is most probably caused by the fact that the computational model assumes linkage equilibrium between markers, a problem that will be further evaluated. Our results do however suggest that SNP-data derived from microarrays are well suited for kinship determination provided linkage disequilibrium is properly accounted for.

Comparison of the Affymetrix SNP Array 5.0 and Oligoarray Platforms for Defining CNV

Together with single nucleotide polymorphism (SNP), copy number variations (CNV) are recognized to be the major component of human genetic diversity and used as a genetic marker in many disease association studies. Affymetrix Genome-wide SNP 5.0 is one of the commonly used SNP array platforms for SNP-GWAS as well as CNV analysis. However, there has been no report that validated the accuracy and reproducibility of CNVs identified by Affymetrix SNP array 5.0. In this study, we compared the characteristics of CNVs from the same set of genomic DNAs detected by three different array platforms; Affymetrix SNP array 5.0, Agilent 2X244K CNV array and NimbleGen 2.1M CNV array. In our analysis, Affymetrix SNP array 5.0 seems to detect CNVs in a reliable manner, which can be applied for association studies. However, for the purpose of defining CNVs in detail, Affymetrix Genome-wide SNP 5.0 might be relatively less ideal than NimbleGen 2.1M CNV array and Agilent 2X244K CNV array, which outperform Affymetrix array for defining the small-sized single copy variants. This result will help researchers to select a suitable array platform for CNV analysis.

Upregulation of FOXM1 induces genomic instability in human epidermal keratinocytes

BackgroundThe human cell cycle transcription factor FOXM1 is known to play a key role in regulating timely mitotic progression and accurate chromosomal segregation during cell division. Deregulation of FOXM1 has been linked to a majority of human cancers. We previously showed that FOXM1 was upregulated in basal cell carcinoma and recently reported that upregulation of FOXM1 precedes malignancy in a number of solid human cancer types including oral, oesophagus, lung, breast, kidney, bladder and uterus. This indicates that upregulation of FOXM1 may be an early molecular signal required for aberrant cell cycle and cancer initiation.ResultsThe present study investigated the putative early mechanism of UVB and FOXM1 in skin cancer initiation. We have demonstrated that UVB dose-dependently increased FOXM1 protein levels through protein stabilisation and accumulation rather than de novo mRNA expression in human epidermal keratinocytes. FOXM1 upregulation in primary human keratinocytes triggered pro-apoptotic/DNA-damage checkpoint response genes such as p21, p38 MAPK, p53 and PARP, however, without causing significant cell cycle arrest or cell death. Using a high-resolution Affymetrix genome-wide single nucleotide polymorphism (SNP) mapping technique, we provided the evidence that FOXM1 upregulation in epidermal keratinocytes is sufficient to induce genomic instability, in the form of loss of heterozygosity (LOH) and copy number variations (CNV). FOXM1-induced genomic instability was significantly enhanced and accumulated with increasing cell passage and this instability was increased even further upon exposure to UVB resulting in whole chromosomal gain (7p21.3-7q36.3) and segmental LOH (6q25.1-6q25.3).ConclusionWe hypothesise that prolonged and repeated UVB exposure selects for skin cells bearing stable FOXM1 protein causes aberrant cell cycle checkpoint thereby allowing ectopic cell cycle entry and subsequent genomic instability. The aberrant upregulation of FOXM1 serves as a 'first hit' where cells acquire genomic instability which in turn predisposes cells to a 'second hit' whereby DNA-damage checkpoint response (eg. p53 or p16) is abolished to allow damaged cells to proliferate and accumulate genetic aberrations/mutations required for cancer initiation.