Identification of germline susceptibility loci in ETV6-RUNX1-rearranged childhood acute lymphoblastic leukemia

Acute lymphoblastic leukemia (ALL) is a malignant disease of the white blood cells. The etiology of ALL is believed to be multifactorial and likely to involve an interplay of environmental and genetic variables. We performed a genome-wide association study of 355 750 single-nucleotide polymorphisms (SNPs) in 474 controls and 419 childhood ALL cases characterized by a t(12;21)(p13;q22) — the most common chromosomal translocation observed in childhood ALL — which leads to an ETV6–RUNX1 gene fusion. The eight most strongly associated SNPs were followed-up in 951 ETV6-RUNX1-positive cases and 3061 controls from Germany/Austria and Italy, respectively. We identified a novel, genome-wide significant risk locus at 3q28 (TP63, rs17505102, PCMH=8.94 × 10−9, OR=0.65). The separate analysis of the combined German/Austrian sample only, revealed additional genome-wide significant associations at 11q11 (OR8U8, rs1945213, P=9.14 × 10−11, OR=0.69) and 8p21.3 (near INTS10, rs920590, P=6.12 × 10−9, OR=1.36). These associations and another association at 11p11.2 (PTPRJ, rs3942852, P=4.95 × 10−7, OR=0.72) remained significant in the German/Austrian replication panel after correction for multiple testing. Our findings demonstrate that germline genetic variation can specifically contribute to the risk of ETV6–RUNX1-positive childhood ALL. The identification of TP63 and PTPRJ as susceptibility genes emphasize the role of the TP53 gene family and the importance of proteins regulating cellular processes in connection with tumorigenesis

Corrigendum: An ethnically relevant consensus Korean reference genome is a step towards personal reference genomes.

This corrects the article DOI: 10.1038/ncomms13637

KoVariome: Korean National Standard Reference Variome database of whole genomes with comprehensive SNV, indel, CNV, and SV analyses

High-coverage whole-genome sequencing data of a single ethnicity can provide a useful catalogue of population-specific genetic variations, and provides a critical resource that can be used to more accurately identify pathogenic genetic variants. We report a comprehensive analysis of the Korean population, and present the Korean National Standard Reference Variome (KoVariome). As a part of the Korean Personal Genome Project (KPGP), we constructed the KoVariome database using 5.5 terabases of whole genome sequence data from 50 healthy Korean individuals in order to characterize the benign ethnicity-relevant genetic variation present in the Korean population. In total, KoVariome includes 12.7M single-nucleotide variants (SNVs), 1.7M short insertions and deletions (indels), 4K structural variations (SVs), and 3.6K copy number variations (CNVs). Among them, 2.4M (19%) SNVs and 0.4M (24%) indels were identified as novel. We also discovered selective enrichment of 3.8M SNVs and 0.5M indels in Korean individuals, which were used to filter out 1,271 coding-SNVs not originally removed from the 1,000 Genomes Project when prioritizing disease-causing variants. KoVariome health records were used to identify novel disease-causing variants in the Korean population, demonstrating the value of high-quality ethnic variation databases for the accurate interpretation of individual genomes and the precise characterization of genetic variation

NCI60 Cancer Cell Line Panel Data and RNAi Analysis Help Identify EAF2 as a Modulator of Simvastatin and Lovastatin Response in HCT-116 Cells

Simvastatin and lovastatin are statins traditionally used for lowering serum cholesterol levels. However, there exists evidence indicating their potential chemotherapeutic characteristics in cancer. In this study, we used bioinformatic analysis of publicly available data in order to systematically identify the genes involved in resistance to cytotoxic effects of these two drugs in the NCI60 cell line panel. We used the pharmacological data available for all the NCI60 cell lines to classify simvastatin or lovastatin resistant and sensitive cell lines, respectively. Next, we performed whole-genome single marker case-control association tests for the lovastatin and simvastatin resistant and sensitive cells using their publicly available Affymetrix 125K SNP genomic data. The results were then evaluated using RNAi methodology. After correction of the p-values for multiple testing using False Discovery Rate, our results identified three genes (NRP1, COL13A1, MRPS31) and six genes (EAF2, ANK2, AKAP7, STEAP2, LPIN2, PARVB) associated with resistance to simvastatin and lovastatin, respectively. Functional validation using RNAi confirmed that silencing of EAF2 expression modulated the response of HCT-116 colon cancer cells to both statins. In summary, we have successfully utilized the publicly available data on the NCI60 cell lines to perform whole-genome association studies for simvastatin and lovastatin. Our results indicated genes involved in the cellular response to these statins and siRNA studies confirmed the role of the EAF2 in response to these drugs in HCT-116 colon cancer cells

A Whole-Genome SNP Association Study of NCI60 Cell Line Panel Indicates a Role of Ca2+ Signaling in Selenium Resistance

Epidemiological studies have suggested an association between selenium intake and protection from a variety of cancer. Considering this clinical importance of selenium, we aimed to identify the genes associated with resistance to selenium treatment. We have applied a previous methodology developed by our group, which is based on the genetic and pharmacological data publicly available for the NCI60 cancer cell line panel. In short, we have categorized the NCI60 cell lines as selenium resistant and sensitive based on their growth inhibition (GI50) data. Then, we have utilized the Affymetrix 125K SNP chip data available and carried out a genome-wide case-control association study for the selenium sensitive and resistant NCI60 cell lines. Our results showed statistically significant association of four SNPs in 5q33–34, 10q11.2, 10q22.3 and 14q13.1 with selenium resistance. These SNPs were located in introns of the genes encoding for a kinase-scaffolding protein (AKAP6), a membrane protein (SGCD), a channel protein (KCNMA1), and a protein kinase (PRKG1). The knock-down of KCNMA1 by siRNA showed increased sensitivity to selenium in both LNCaP and PC3 cell lines. Furthermore, SNP-SNP interaction (epistasis) analysis indicated the interactions of the SNPs in AKAP6 with SGCD as well as SNPs in AKAP6 with KCNMA1 with each other, assuming additive genetic model. These genes were also all involved in the Ca2+ signaling, which has a direct role in induction of apoptosis and induction of apoptosis in tumor cells is consistent with the chemopreventive action of selenium. Once our findings are further validated, this knowledge can be translated into clinics where individuals who can benefit from the chemopreventive characteristics of the selenium supplementation will be easily identified using a simple DNA analysis

Recessive mutations in SPTBN2 implicate β-III spectrin in both cognitive and motor development

β-III spectrin is present in the brain and is known to be important in the function of the cerebellum. Heterozygous mutations in SPTBN2, the gene encoding β-III spectrin, cause Spinocerebellar Ataxia Type 5 (SCA5), an adult-onset, slowly progressive, autosomal-dominant pure cerebellar ataxia. SCA5 is sometimes known as "Lincoln ataxia," because the largest known family is descended from relatives of the United States President Abraham Lincoln. Using targeted capture and next-generation sequencing, we identified a homozygous stop codon in SPTBN2 in a consanguineous family in which childhood developmental ataxia co-segregates with cognitive impairment. The cognitive impairment could result from mutations in a second gene, but further analysis using whole-genome sequencing combined with SNP array analysis did not reveal any evidence of other mutations. We also examined a mouse knockout of β-III spectrin in which ataxia and progressive degeneration of cerebellar Purkinje cells has been previously reported and found morphological abnormalities in neurons from prefrontal cortex and deficits in object recognition tasks, consistent with the human cognitive phenotype. These data provide the first evidence that β-III spectrin plays an important role in cortical brain development and cognition, in addition to its function in the cerebellum; and we conclude that cognitive impairment is an integral part of this novel recessive ataxic syndrome, Spectrin-associated Autosomal Recessive Cerebellar Ataxia type 1 (SPARCA1). In addition, the identification of SPARCA1 and normal heterozygous carriers of the stop codon in SPTBN2 provides insights into the mechanism of molecular dominance in SCA5 and demonstrates that the cell-specific repertoire of spectrin subunits underlies a novel group of disorders, the neuronal spectrinopathies, which includes SCA5, SPARCA1, and a form of West syndrome

Epistasis between COMT and MTHFR in Maternal-Fetal Dyads Increases Risk for Preeclampsia

Preeclampsia is a leading cause of perinatal morbidity and mortality. This disorder is thought to be multifactorial in origin, with multiple genes, environmental and social factors, contributing to disease. One proposed mechanism is placental hypoxia-driven imbalances in angiogenic and anti-angiogenic factors, causing endothelial cell dysfunction. Catechol-O-methyltransferase (Comt)-deficient pregnant mice have a preeclampsia phenotype that is reversed by exogenous 2-methoxyestradiol (2-ME), an estrogen metabolite generated by COMT. 2-ME inhibits Hypoxia Inducible Factor 1α, a transcription factor mediating hypoxic responses. COMT has been shown to interact with methylenetetrahydrofolate reductase (MTHFR), which modulates the availability of S-adenosylmethionine (SAM), a COMT cofactor. Variations in MTHFR have been associated with preeclampsia. By accounting for allelic variation in both genes, the role of COMT has been clarified. COMT allelic variation is linked to enzyme activity and four single nucleotide polymorphisms (SNPs) (rs6269, rs4633, rs4680, and rs4818) form haplotypes that characterize COMT activity. We tested for association between COMT haplotypes and the MTHFR 677 C→T polymorphism and preeclampsia risk in 1103 Chilean maternal-fetal dyads. The maternal ACCG COMT haplotype was associated with reduced risk for preeclampsia (P = 0.004), and that risk increased linearly from low to high activity haplotypes (P = 0.003). In fetal samples, we found that the fetal ATCA COMT haplotype and the fetal MTHFR minor “T” allele interact to increase preeclampsia risk (p = 0.022). We found a higher than expected number of patients with preeclampsia with both the fetal risk alleles alone (P = 0.052) and the fetal risk alleles in combination with a maternal balancing allele (P<0.001). This non-random distribution was not observed in controls (P = 0.341 and P = 0.219, respectively). Our findings demonstrate a role for both maternal and fetal COMT in preeclampsia and highlight the importance of including allelic variation in MTHFR

Patterns of Ancestry, Signatures of Natural Selection, and Genetic Association with Stature in Western African Pygmies

African Pygmy groups show a distinctive pattern of phenotypic variation, including short stature, which is thought to reflect past adaptation to a tropical environment. Here, we analyze Illumina 1M SNP array data in three Western Pygmy populations from Cameroon and three neighboring Bantu-speaking agricultural populations with whom they have admixed. We infer genome-wide ancestry, scan for signals of positive selection, and perform targeted genetic association with measured height variation. We identify multiple regions throughout the genome that may have played a role in adaptive evolution, many of which contain loci with roles in growth hormone, insulin, and insulin-like growth factor signaling pathways, as well as immunity and neuroendocrine signaling involved in reproduction and metabolism. The most striking results are found on chromosome 3, which harbors a cluster of selection and association signals between approximately 45 and 60 Mb. This region also includes the positional candidate genes DOCK3, which is known to be associated with height variation in Europeans, and CISH, a negative regulator of cytokine signaling known to inhibit growth hormone-stimulated STAT5 signaling. Finally, pathway analysis for genes near the strongest signals of association with height indicates enrichment for loci involved in insulin and insulin-like growth factor signaling

A second generation human haplotype map of over 3.1 million SNPs

We describe the Phase II HapMap, which characterizes over 3.1 million human single nucleotide polymorphisms (SNPs) genotyped in 270 individuals from four geographically diverse populations and includes 25-35% of common SNP variation in the populations surveyed. The map is estimated to capture untyped common variation with an average maximum r(2) of between 0.9 and 0.96 depending on population. We demonstrate that the current generation of commercial genome-wide genotyping products captures common Phase II SNPs with an average maximum r(2) of up to 0.8 in African and up to 0.95 in non-African populations, and that potential gains in power in association studies can be obtained through imputation. These data also reveal novel aspects of the structure of linkage disequilibrium. We show that 10-30% of pairs of individuals within a population share at least one region of extended genetic identity arising from recent ancestry and that up to 1% of all common variants are untaggable, primarily because they lie within recombination hotspots. We show that recombination rates vary systematically around genes and between genes of different function. Finally, we demonstrate increased differentiation at non-synonymous, compared to synonymous, SNPs, resulting from systematic differences in the strength or efficacy of natural selection between populations.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/62863/1/nature06258.pd