Neighboring Single Nucleotide Polymorphisms Research Articles

Genome-Wide Detection of Structural Variations Reveals New Regions Associated with Domestication in Small Ruminants.

During domestication processes, changes in selective pressures induce multiple phenotypical, physiological, and behavioral changes in target species. The rise of next-generation sequencing has provided a chance to study the genetics bases of these changes, most of the time based on single nucleotide polymorphisms (SNPs). However, several studies have highlighted the impact of structural variations (SVs) on individual fitness, particularly in domestic species. We aimed at unraveling the role of SVs during the domestication and later improvement of small ruminants by analyzing whole-genome sequences of 40 domestic sheep and 11 of their close wild relatives (Ovis orientalis), and 40 goats and 18 of their close wild relatives (Capra aegagrus). Using a combination of detection tools, we called 45,796 SVs in Ovis and 15,047 SVs in Capra genomes, including insertions, deletions, inversions, copy number variations, and chromosomal translocations. Most of these SVs were previously unreported in small ruminants. 69 and 45 SVs in sheep and goats, respectively, were in genomic regions with neighboring SNPs highly differentiated between wilds and domestics (i.e., putatively related to domestication). Among them, 25 and 20 SVs were close to or overlapping with genes related to physiological and morpho-anatomical traits linked with productivity (e.g., size, meat or milk quality, wool color), reproduction, or immunity. Finally, several of the SVs differentiated between wilds and domestics would not have been detected by screening only the differentiation of SNPs surrounding them, highlighting the complementarity of SVs and SNPs based approaches to detect signatures of selection.

Development of a multiplex methylation-sensitive restriction enzyme-based SNP typing system for deconvolution of semen-containing mixtures.

The identification of mixed stains has always been a difficult problem in personal identification in the forensic field. In recent years, tissue-specific methylation sites have proven to be very stable biomarkers for distinguishing tissue origin. However, it is still challenging to perform tissue source identification and individual identification simultaneously. In this study, we developed a method that uses tissue-specific methylation markers combined with single-nucleotide polymorphism (SNP) markers to detect semen from mixed biofluids and to identify individuals simultaneously. Semen-specific CpG markers were chosen from the literature and further validated utilizing methylation-sensitive restriction endonuclease (MSRE) combined with PCR technology. The neighboring SNP markers were searched in the flanking sequence of the target CpG within 400bp, and SNP typing was then carried out through a single-base extension reaction followed by capillary electrophoresis. Eventually, a method of MSRE combined with SNaPshot that could detect 12 compound CpG-SNP markers was developed. Using this system, 10ng of total DNA and DNA mixture with semen content up to 25% could be typed successfully. Moreover, the cumulative discrimination power of the system in the northern Chinese Han population is 0.9998. This study provides a valuable strategy for forensic practice to perform tissue origin and individual identification from mixed stains simultaneously.

Interweaving Between Genetic and Epigenetic Studies on Childhood Asthma.

The cause and underlying mechanisms that contribute to asthma pathogenesis are not well known. Both genome- and epigenome-wide association studies have identified genes associated with asthma risk. It is unknown to what extent genes identified in these two types of studies overlap. Based on existing literature and the DisGeNET database, we extracted overlapping genes identified in genetic and epigenetic studies of childhood asthma. Through analyses of variance, we assessed whether DNA methylation (DNAm) at 5′-C-phosphate-G-3′ (CpGs) on the overlapping genes was associated with neighboring single-nucleotide polymorphisms (SNPs) within 1M base pairs (bps) and with low linkage disequilibrium (r2 < 0.2) in the childhood asthma-related genes. In total, 285 genes from genetic studies and 226 genes from epigenetic studies were shown to be associated with asthma risk, of which six overlap. Of the six genes, 79 CpGs and 8229 unique neighboring SNPs (1M bps) were included in methylation quantitative loci (methQTL) assessment analyses. We tested the association of DNAm at each of the 79 CpG sites with its neighboring SNPs. After adjusting for multiple testing by controlling the false discovery rate to 0.05 when testing methQTL for each CpG site, we found statistically significant associations in three genes with their neighboring SNPs and identified 34 unique methQTLs. The rather limited overlap in genes between genetic and epigenetic studies on asthma and the absence of methQTL in some of the overlapping genes highlight a need to jointly, rather than independently, examine genetic and epigenetic effects on asthma risk to improve our understanding of the underlying mechanisms of asthma.

Replicability analysis in genome-wide association studies via Cartesian hidden Markov models

BackgroundReplicability analysis which aims to detect replicated signals attracts more and more attentions in modern scientific applications. For example, in genome-wide association studies (GWAS), it would be of convincing to detect an association which can be replicated in more than one study. Since the neighboring single nucleotide polymorphisms (SNPs) often exhibit high correlation, it is desirable to exploit the dependency information among adjacent SNPs properly in replicability analysis. In this paper, we propose a novel multiple testing procedure based on the Cartesian hidden Markov model (CHMM), called repLIS procedure, for replicability analysis across two studies, which can characterize the local dependence structure among adjacent SNPs via a four-state Markov chain.ResultsTheoretical results show that the repLIS procedure can control the false discovery rate (FDR) at the nominal level α and is shown to be optimal in the sense that it has the smallest false non-discovery rate (FNR) among all α-level multiple testing procedures. We carry out simulation studies to compare our repLIS procedure with the existing methods, including the Benjamini-Hochberg (BH) procedure and the empirical Bayes approach, called repfdr. Finally, we apply our repLIS procedure and repfdr procedure in the replicability analyses of psychiatric disorders data sets collected by Psychiatric Genomics Consortium (PGC) and Wellcome Trust Case Control Consortium (WTCCC). Both the simulation studies and real data analysis show that the repLIS procedure is valid and achieves a higher efficiency compared with its competitors.ConclusionsIn replicability analysis, our repLIS procedure controls the FDR at the pre-specified level α and can achieve more efficiency by exploiting the dependency information among adjacent SNPs.

Complex polymorphisms in endocytosis genes suggest alpha-cyclodextrin as a treatment for breast cancer.

Most breast cancer deaths are caused by metastasis and treatment options beyond radiation and cytotoxic drugs, which have severe side effects, and hormonal treatments, which are or become ineffective for many patients, are urgently needed. This study reanalyzed existing data from three genome-wide association studies (GWAS) using a novel computational biostatistics approach (muGWAS), which had been validated in studies of 600–2000 subjects in epilepsy and autism. MuGWAS jointly analyzes several neighboring single nucleotide polymorphisms while incorporating knowledge about genetics of heritable diseases into the statistical method and about GWAS into the rules for determining adaptive genome-wide significance. Results from three independent GWAS of 1000–2000 subjects each, which were made available under the National Institute of Health’s “Up For A Challenge” (U4C) project, not only confirmed cell-cycle control and receptor/AKT signaling, but, for the first time in breast cancer GWAS, also consistently identified many genes involved in endo-/exocytosis (EEC), most of which had already been observed in functional and expression studies of breast cancer. In particular, the findings include genes that translocate (ATP8A1, ATP8B1, ANO4, ABCA1) and metabolize (AGPAT3, AGPAT4, DGKQ, LPPR1) phospholipids entering the phosphatidylinositol cycle, which controls EEC. These novel findings suggest scavenging phospholipids as a novel intervention to control local spread of cancer, packaging of exosomes (which prepare distant microenvironment for organ-specific metastases), and endocytosis of β1 integrins (which are required for spread of metastatic phenotype and mesenchymal migration of tumor cells). Beta-cyclodextrins (βCD) have already been shown to be effective in in vitro and animal studies of breast cancer, but exhibits cholesterol-related ototoxicity. The smaller alpha-cyclodextrins (αCD) also scavenges phospholipids, but cannot fit cholesterol. An in-vitro study presented here confirms hydroxypropyl (HP)-αCD to be twice as effective as HPβCD against migration of human cells of both receptor negative and estrogen-receptor positive breast cancer. If the previous successful animal studies with βCDs are replicated with the safer and more effective αCDs, clinical trials of adjuvant treatment with αCDs are warranted. Ultimately, all breast cancer are expected to benefit from treatment with HPαCD, but women with triple-negative breast cancer (TNBC) will benefit most, because they have fewer treatment options and their cancer advances more aggressively.

Abstract 2261: Synergistic effect and VEGF/HSP70-hom haplotype analysis: Relationship to prostate cancer risk and clinical outcome

Abstract Prostate cancer (PCa) is a complex disorder resulting from the combined effects of multiple environmental and genetic factors. Our previous single-locus analysis showed that VEGF and HSP70-hom polymorphisms were significantly associated with PCa susceptibility and prognosis. Both genes encoding these proteins were located on chromosome 6p21, and combining the neighboring single nucleotide polymorphisms (SNPs) into haplotypes may increase the association with the disease. Three tagging polymorphisms, the HSP70-hom 2437 T/C, the VEGF-1154 G/A, and the VEGF-634 G/C SNPs were genotyped in 101 cases and 80 controls. For the combined analysis of VEGF and HSP70-hom, we found a positive gradient in the odds ratios (ORs) related to the number of high-risk genotypes with a 3.53-fold increase of prostate carcinoma risk (OR= 3.53; p= 0.015). Furthermore, the TAG and CAG haplotypes at positions HSP70-hom, VEGF-1154 and VEGF -634 exhibited a two-fold (OR= 0.46; p= 0.014) and a seven-fold (OR= 0.14; p= 0.00005) reduction in PCa risk, respectively. Regarding PCa prognosis, the TAG haplotype had a negative association with the aggressive phenotype as defined by the histopathological grade (OR= 0.28; p= 0.006). Our findings confirm the role of at-risk haplotype across the HSP70-hom/VEGF gene cluster in determining susceptibility to PCa. Note: This abstract was not presented at the meeting. Citation Format: Sana Sfar, Hamadi Saad, Faouzi Mosbah, Lotfi Chouchane. Synergistic effect and VEGF/HSP70-hom haplotype analysis: Relationship to prostate cancer risk and clinical outcome [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 2261. doi:10.1158/1538-7445.AM2017-2261

Evaluation of a blood-specific DNA methylated region and trial for allele-specific blood identification from mixed body fluid DNA.

The identification of blood samples obtained from crime scenes has been an important step in forensic investigation. Recently, a novel approach using the blood-specific methylated CpG site cg06379435 has been reported. In this study, we developed a real-time polymerase-chain-reaction-based method that can simply and rapidly quantitate the methylation ratio of cg06379435 and its neighboring CpGs and set the threshold ratios for blood identification by analyzing various body fluid samples. Blood identification using the thresholds was successfully performed in the analysis of a small amount (1ng) of DNA from blood and various aged blood samples, including 29-year-old stains. We also demonstrated a test for allele-specific blood identification from a mixed DNA sample by bisulfite sequencing analysis of these CpG sites and their neighboring single nucleotide polymorphism, rs7359943 (A/G), which is of relevance in cases where mixed samples are obtained from crime scenes. The stability of DNA methylation in aged samples and the usefulness of neighboring genetic information shown in this study suggest that DNA-methylation-based body fluid identification will play a major role in future forensic investigations.

CNTNAP2 Is Significantly Associated With Speech Sound Disorder in the Chinese Han Population.

Speech sound disorder is the most common communication disorder. Some investigations support the possibility that the CNTNAP2 gene might be involved in the pathogenesis of speech-related diseases. To investigate single-nucleotide polymorphisms in the CNTNAP2 gene, 300 unrelated speech sound disorder patients and 200 normal controls were included in the study. Five single-nucleotide polymorphisms were amplified and directly sequenced. Significant differences were found in the genotype (P = .0003) and allele (P = .0056) frequencies of rs2538976 between patients and controls. The excess frequency of the A allele in the patient group remained significant after Bonferroni correction (P = .0280). A significant haplotype association with rs2710102T/+rs17236239A/+2538976A/+2710117A (P = 4.10e-006) was identified. A neighboring single-nucleotide polymorphism, rs10608123, was found in complete linkage disequilibrium with rs2538976, and the genotypes exactly corresponded to each other. The authors propose that these CNTNAP2 variants increase the susceptibility to speech sound disorder. The single-nucleotide polymorphisms rs10608123 and rs2538976 may merge into one single-nucleotide polymorphism.

A nonhomogeneous hidden markov model for gene mapping based on next-generation sequencing data.

The analysis of polygenetic characteristics for mapping quantitative trait loci (QTL) remains an important challenge. QTL analysis requires two or more strains of organisms that differ substantially in the (poly-)genetic trait of interest, resulting in a heterozygous offspring. The offspring with the trait of interest is selected and subsequently screened for molecular markers such as single-nucleotide polymorphisms (SNPs) with next-generation sequencing. Gene mapping relies on the co-segregation between genes and/or markers. Genes and/or markers that are linked to a QTL influencing the trait will segregate more frequently with this locus. For each identified marker, observed mismatch frequencies between the reads of the offspring and the parental reference strains can be modeled by a multinomial distribution with the probabilities depending on the state of an underlying, unobserved Markov process. The states indicate whether the SNP is located in a (vicinity of a) QTL or not. Consequently, genomic loci associated with the QTL can be discovered by analyzing hidden states along the genome. The aforementioned hidden Markov model assumes that the identified SNPs are equally distributed along the chromosome and does not take the distance between neighboring SNPs into account. The distance between the neighboring SNPs could influence the chance of co-segregation between genes and markers. To address this issue, we propose a nonhomogeneous hidden Markov model with a transition matrix that depends on a set of distance-varying observed covariates. The application of the model is illustrated on the data from a study of ethanol tolerance in yeast.

Revisiting the role of MCL1 in tumorigenesis of solid cancer: gene expression correlates with antiproliferative phenotype in breast cancer cells and its functional regulatory variants are associated with reduced cancer susceptibility.

Compared to the well-defined anti-apoptotic role of myeloid cell leukemia sequence 1 (MCL1), its antiproliferative function in tumorigenesis is less studied. We had recently reported that regulatory variants of MCL1 contribute to enhanced promoter activity but reduced risk of lung cancer. We hypothesized that MCL1 expression may manifest antiproliferative phenotype and its functional variations may have etiological relevance for breast cancer. We manipulated MCL1 expression in MCF-7 cells and MDA231 with overexpression and knockdown, analyzed the effects on cell viability and cell cycling phase, and characterized the correlation with expression profiles of key regulators of cell cycle. We further genotyped the -190 insertion polymorphism and the neighboring single nucleotide polymorphisms (SNPs) in 745 breast cancer patients and 537 controls and analyzed their association with cancer risk. We confirmed that heightened expression of MCL1 resulted in decreased proliferation ability of breast cancer cells. We further observed that MCL1 overexpression in breast cancer cells resulted in cell cycle progression arresting in S phase and concomitant enhanced expression of p27, which could be rescued by p27 knockdown with co-transfection of small interfering RNA (siRNA). Furthermore, we found a significant reduction in breast cancer risk [odds ratio (OR) = 0.74; 95 % confidence interval (CI) = 0.59-0.93] associated with -190 insertion genotype; the expression-enhancing regulatory haplotype (OR 0.79; 95 % CI 0.66-0.95) and diplotype (OR 0.71; 95 % CI 0.57-0.89) were consistently associated with decreased cancer susceptibility. The study demonstrates that the expression-enhancing regulatory variants of MCL1 are protective modifiers of breast cancer risk, and reduced cell proliferation and arrested cell cycle progression partly mediated by p27 might be the underlying mechanism.

Association mapping of stem rust race TTKSK resistance in US barley breeding germplasm

Key message Loci conferring resistance to the highly virulent African stem rust race TTKSK were identified in advanced barley breeding germplasm and positioned to chromosomes 5H and 7H using an association mapping approach. African races of the stem rust pathogen (Puccinia graminis f. sp. tritici) are a serious threat to barley production worldwide because of their wide virulence. To discover and characterize resistance to African stem rust race TTKSK in US barley breeding germplasm, over 3,000 lines/cultivars were assessed for resistance at the seedling stage in the greenhouse and also the adult plant stage in the field in Kenya. Only 12 (0.3 %) and 64 (2.1 %) lines exhibited a resistance level comparable to the resistant control at the seedling and adult plant stage, respectively. To map quantitative trait loci (QTL) for resistance to race TTKSK, an association mapping approach was conducted, utilizing 3,072 single nucleotide polymorphism (SNP) markers. At the seedling stage, two neighboring SNP markers (0.8 cM apart) on chromosome 7H (11_21491 and 12_30528) were found significantly associated with resistance. The most significant one found was 12_30528; thus, the resistance QTL was named Rpg-qtl-7H-12_30528. At the adult plant stage, two SNP markers on chromosome 5H (11_11355 and 12_31427) were found significantly associated with resistance. This resistance QTL was named Rpg-qtl-5H-11_11355 for the most significant marker identified. Adult plant resistance is of paramount importance for stem rust. The marker associated with Rpg-qtl-5H-11_11355 for adult plant resistance explained only a small portion of the phenotypic variation (0.02); however, this QTL reduced disease severity up to 55.0 % under low disease pressure and up to 21.1 % under heavy disease pressure. SNP marker 11_11355 will be valuable for marker-assisted selection of adult plant stem rust resistance in barley breeding.Electronic supplementary materialThe online version of this article (doi:10.1007/s00122-014-2297-8) contains supplementary material, which is available to authorized users.

Polymorphic L1 retrotransposons are frequently in strong linkage disequilibrium with neighboring SNPs

L1 retrotransposons have been the major driver of structural variation of the human genome. L1 insertion polymorphism (LIP)-mediated genomic variation can alter the transcriptome and contribute to the divergence of human phenotypes. To assess this possibility, a genome-wide association study (GWAS) including LIPs is required. Toward this ultimate goal, the present study examined linkage disequilibrium between six LIPs and their neighboring single nucleotide polymorphisms (SNPs). Genomic PCR and sequencing of L1-plus and -minus alleles from different donors revealed that all six LIPs were in strong linkage disequilibrium with at least one SNP. In addition, comparison of syntenic regions containing the identified SNP nucleotides was performed among modern humans (L1-plus and -minus alleles), archaic humans and non-human primates, revealing two different evolutionary schemes that might have resulted in the observed strong SNP-LIP linkage disequilibria. This study provides an experimental framework and guidance for a future SNP-LIP integrative GWAS.

The Role of Memory-related Gene WWC1 (KIBRA) in Lifetime Posttraumatic Stress Disorder: Evidence from Two Independent Samples from African Conflict Regions

Posttraumatic stress disorder (PTSD) results from the formation of a strong memory for the sensory-perceptual and affective representations of traumatic experiences, which is detached from the corresponding autobiographical context information. Because WWC1, the gene encoding protein KIBRA, is associated with long-term memory performance, we hypothesized that common WWC1 alleles influence the risk for a lifetime diagnosis of PTSD. Traumatic load and diagnosis of current and lifetime PTSD were assessed in two independent African samples of survivors from conflict zones who had faced severe trauma (n = 392, Rwanda, and n = 399, Northern Uganda, respectively). Array-based single nucleotide polymorphism (SNP) genotyping was performed. The influence of WWC1 tagging SNPs and traumatic load on lifetime PTSD was estimated by means of logistic regression models with correction for multiple comparisons in the Rwandan sample. Replication analysis was performed in the independent Ugandan sample. An association of two neighboring SNPs in almost complete linkage disequilibrium, rs10038727 and rs4576167, with lifetime PTSD was discovered in the Rwandan sample. Although each traumatic event added to the probability of lifetime PTSD in a dose-dependent manner in both genotype groups, carriers of the minor allele of both SNPs displayed a diminished risk (p = .007, odds ratio = .29 [95% confidence interval = .15-.54]). This effect was confirmed in the independent Ugandan sample. This study reveals an association between two WWC1 SNPs and the likelihood of PTSD development, indicating that this memory-related gene might be involved in processes that occur in response to traumatic stress and influence the strengthening of fear memories.

Prediction of a deletion copy number variant by a dense SNP panel

BackgroundA newly recognized type of genetic variation, Copy Number Variation (CNV), is detected in mammalian genomes, e.g. the cattle genome. This form of variation can potentially cause phenotypic variation. Our objective was to determine whether dense SNP (single nucleotide polymorphisms) panels can capture the genetic variation due to a simple bi-allelic CNV, with the prospect of including the effect of such structural variations into genomic predictions.MethodsA deletion type CNV on bovine chromosome 6 was predicted from its neighboring SNP with a multiple regression model. Our dataset consisted of CNV genotypes of 1,682 cows, along with 100 surrounding SNP genotypes. A prediction model was fitted considering 10 to 100 surrounding SNP and the accuracy obtained directly from the model was confirmed by cross-validation.Results and conclusionsThe accuracy of prediction increased with an increasing number of SNP in the model and the predicted accuracies were similar to those obtained by cross-validation. A substantial increase in accuracy was observed when the number of SNP increased from 10 to 50 but thereafter the increase was smaller, reaching the highest accuracy (0.94) with 100 surrounding SNP. Thus, we conclude that the genotype of a deletion type CNV and its putative QTL effect can be predicted with a maximum accuracy of 0.94 from surrounding SNP. This high prediction accuracy suggests that genetic variation due to simple deletion CNV is well captured by dense SNP panels. Since genomic selection relies on the availability of a dense marker panel with markers in close linkage disequilibrium to the QTL in order to predict their genetic values, we also discuss opportunities for genomic selection to predict the effects of CNV by dense SNP panels, when CNV cause variation in quantitative traits.

Identifying disease-associated SNP clusters via contiguous outlier detection

Although genome-wide association studies (GWAS) have identified many disease-susceptibility single-nucleotide polymorphisms (SNPs), these findings can only explain a small portion of genetic contributions to complex diseases, which is known as the missing heritability. A possible explanation is that genetic variants with small effects have not been detected. The chance is < 8 that a causal SNP will be directly genotyped. The effects of its neighboring SNPs may be too weak to be detected due to the effect decay caused by imperfect linkage disequilibrium. Moreover, it is still challenging to detect a causal SNP with a small effect even if it has been directly genotyped. In order to increase the statistical power when detecting disease-associated SNPs with relatively small effects, we propose a method using neighborhood information. Since the disease-associated SNPs account for only a small fraction of the entire SNP set, we formulate this problem as Contiguous Outlier DEtection (CODE), which is a discrete optimization problem. In our formulation, we cast the disease-associated SNPs as outliers and further impose a spatial continuity constraint for outlier detection. We show that this optimization can be solved exactly using graph cuts. We also employ the stability selection strategy to control the false positive results caused by imperfect parameter tuning. We demonstrate its advantage in simulations and real experiments. In particular, the newly identified SNP clusters are replicable in two independent datasets. The software is available at: http://bioinformatics.ust.hk/CODE.zip. eeyu@ust.hk Supplementary data are available at Bioinformatics online.

Association weight matrix for the genetic dissection of puberty in beef cattle

We describe a systems biology approach for the genetic dissection of complex traits based on applying gene network theory to the results from genome-wide associations. The associations of single-nucleotide polymorphisms (SNP) that were individually associated with a primary phenotype of interest, age at puberty in our study, were explored across 22 related traits. Genomic regions were surveyed for genes harboring the selected SNP. As a result, an association weight matrix (AWM) was constructed with as many rows as genes and as many columns as traits. Each {i, j} cell value in the AWM corresponds to the z-score normalized additive effect of the ith gene (via its neighboring SNP) on the jth trait. Columnwise, the AWM recovered the genetic correlations estimated via pedigree-based restricted maximum-likelihood methods. Rowwise, a combination of hierarchical clustering, gene network, and pathway analyses identified genetic drivers that would have been missed by standard genome-wide association studies. Finally, the promoter regions of the AWM-predicted targets of three key transcription factors (TFs), estrogen-related receptor gamma (ESRRG), Pal3 motif, bound by a PPAR-gamma homodimer, IR3 sites (PPARG), and Prophet of Pit 1, PROP paired-like homeobox 1 (PROP1), were surveyed to identify binding sites corresponding to those TFs. Applied to our case, the AWM results recapitulate the known biology of puberty, captured experimentally validated binding sites, and identified candidate genes and gene-gene interactions for further investigation.

Synergistic effect and VEGF/HSP70-hom haplotype analysis: Relationship to prostate cancer risk and clinical outcome

Prostate cancer (PCa) is a complex disorder resulting from the combined effects of multiple environmental and genetic factors. Our previous single-locus analysis showed that VEGF and HSP70-hom polymorphisms were significantly associated with PCa susceptibility and prognosis. Both genes encoding these proteins were located on chromosome 6p21, and combining the neighboring single nucleotide polymorphisms (SNPs) into haplotypes may increase the association with the disease. Three tagging polymorphisms, the HSP70-hom 2437 T/C, the VEGF-1154 G/A, and the VEGF-634 G/C SNPs were genotyped in 101 cases and 80 controls. For the combined analysis of VEGF and HSP70-hom, we found a positive gradient in the odds ratios (ORs) related to the number of high-risk genotypes with a 3.53-fold increase of prostate carcinoma risk (OR = 3.53; p = 0.015). Furthermore, the TAG and CAG haplotypes at positions HSP70-hom, VEGF −1154 and VEGF −634 exhibited a two-fold (OR = 0.46; p = 0.014) and a seven-fold (OR = 0.14; p = 0.00005) reduction in PCa risk, respectively. Regarding PCa prognosis, the TAG haplotype had a negative association with the aggressive phenotype as defined by the histopathological grade (OR = 0.28; p = 0.006). Our findings confirm the role of at-risk haplotype across the HSP70-hom/VEGF gene cluster in determining susceptibility to PCa.

Simultaneous analysis of all single-nucleotide polymorphisms in genome-wide association study of rheumatoid arthritis

The availability of very large number of markers by modern technology makes genome-wide association studies very popular. The usual approach is to test single-nucleotide polymorphisms (SNPs) one at a time for association with disease status. However, it may not be possible to detect marginally significant effects by single-SNP analysis. Simultaneous analysis of SNPs enables detection of even those SNPs with small effect by evaluating the collective impact of several neighboring SNPs. Also, false-positive signals may be weakened by the presence of other neighboring SNPs included in the analysis. We analyzed the North American Rheumatoid Arthritis Consortium data of Genetic Analysis Workshop 16 using HLasso, a new method for simultaneous analysis of SNPs. The simultaneous analysis approach has excellent control of type I error, and many of the previously reported results of single-SNP analyses were confirmed by this approach.

First genome-wide association scan on neurophysiological endophenotypes points to trans-regulation effects on SLC2A3 in dyslexic children

Dyslexia is one of the most common learning disorders affecting about 5% of all school-aged children. It has been shown that event-related potential measurements reveal differences between dyslexic children and age-matched controls. This holds particularly true for mismatch negativity (MMN), which reflects automatic speech deviance processing and is altered in dyslexic children. We performed a whole-genome association analysis in 200 dyslexic children, focusing on MMN measurements. We identified rs4234898, a marker located on chromosome 4q32.1, to be significantly associated with the late MMN component. This association could be replicated in an independent second sample of 186 dyslexic children, reaching genome-wide significance in the combined sample (P = 5.14e-08). We also found an association between the late MMN component and a two-marker haplotype of rs4234898 and rs11100040, one of its neighboring single nucleotide polymorphisms (SNPs). In the combined sample, this marker combination withstands correction for multiple testing (P = 6.71e-08). Both SNPs lie in a region devoid of any protein-coding genes; however, they both show significant association with mRNA-expression levels of SLC2A3 on chromosome 12, the predominant facilitative glucose transporter in neurons. Our results suggest a possible trans-regulation effect on SLC2A3, which might lead to glucose deficits in dyslexic children and could explain their attenuated MMN in passive listening tasks.

Polyclonality of BRAF Mutations in Acquired Melanocytic Nevi

Melanocytic nevi are thought to be senescent clones of melanocytes that have acquired an oncogenic BRAF mutation. BRAF mutation is considered to be a crucial step in the initiation of melanocyte transformation. However, using immunomagnetic separation or laser-capture microdissection, we examined BRAF mutations in sets of approximately 50 single cells isolated from acquired melanocytic nevi from 13 patients and found a substantial number of nevus cells that contained wild-type BRAF mixed with nevus cells that contained BRAF(V600E). Furthermore, we simultaneously amplified BRAF exon 15 and a neighboring single nucleotide polymorphism (SNP), rs7801086, from nevus cell samples obtained from four patients who were heterozygous for this SNP. Subcloning and sequencing of the polymerase chain reaction products showed that both SNP alleles harbored the BRAF(V600E) mutation, indicating that the same BRAF(V600E) mutation originated from different cells. The polyclonality of BRAF mutations in acquired melanocytic nevi suggests that mutation of BRAF may not be an initial event in melanocyte transformation.