A HIDDEN MARKOV MODEL FOR JOINT ESTIMATION OF GENOTYPE AND COPY NUMBER IN HIGH-THROUGHPUT SNP CHIPS

Amplifications and deletions of chromosomal DNA, as well as copy-neutral loss of heterozygosity have been associated with diseases processes. High-throughput single nucleotide polymorphism (SNP) arrays are useful for making genome-wide estimates of copy number and genotype calls. Because neighboring SNPs in high throughput SNP arrays are likely to have dependent copy number and genotype due to the underlying haplotype structure and linkage disequilibrium, hidden Markov models (HMM) may be useful for improving genotype calls and copy number estimates that do not incorporate information from nearby SNPs. We improve previous approaches that utilize a HMM framework for inference in high throughput SNP arrays by integrating copy number, genotype calls, and the corresponding confidence scores when available. Using simulated data, we demonstrate how confidence scores control smoothing in a probabilistic framework. Software for fitting HMMs to SNP array data is available in the R package ICE.

A noncoding melanophilin gene (MLPH) SNP at the splice donor of exon 1 represents a candidate causal mutation for coat color dilution in dogs

Coat color dilution in several breeds of dog is characterized by a specific pigmentation phenotype and sometimes accompanied by hair loss and recurrent skin inflammation, the so-called color dilution alopecia or black hair follicular dysplasia. Coat color dilution (d) is inherited as a Mendelian autosomal recessive trait. In a previous study, MLPH polymorphisms showed perfect cosegregation with the dilute phenotype within breeds. However, different dilute haplotypes were found in different breeds, and no single polymorphism was identified in the coding sequence that was likely to be causative for the dilute phenotype. We resequenced the 5'-region of the canine MLPH gene and identified a strong candidate single nucleotide polymorphism within the nontranslated exon 1, which showed perfect association to the dilute phenotype in 65 dilute dogs from 7 different breeds. The A/G polymorphism is located at the last nucleotide of exon 1 and the mutant A-allele is predicted to reduce splicing efficiency 8-fold. An MLPH mRNA expression study using quantitative reverse transcriptase-polymerase chain reaction confirmed that dd animals had only about approximately 25% of the MLPH transcript compared with DD animals. These results provide preliminary evidence that the reported regulatory MLPH mutation might represent a causal mutation for coat color dilution in dogs.

A polymorphism within the equine CRISP3 gene is associated with stallion fertility in Hanoverian warmblood horses

Fertility of stallions is of high economic importance, especially for large breeding organisations and studs. Breeding schemes with respect to fertility traits and selection of stallions at an early stage may be improved by including molecular genetic markers associated with traits. The genes coding for equine cysteine-rich secretory proteins (CRISPs) are promising candidate genes because previous studies have shown that CRISPs play a role in the fertilising ability of male animals. We have previously characterised the three equine CRISP genes and identified a non-synonymous polymorphism in the CRISP1 gene. In this study, we report one non-synonymous polymorphism in the CRISP2 gene and four non-synonymous polymorphisms in the CRISP3 gene. All six CRISP polymorphisms were genotyped in 107 Hanoverian breeding stallions. Insemination records of stallions were used to analyse the association between CRISP polymorphisms and fertility traits. Three statistical models were used to evaluate the influence of single mutations, genotypes and haplotypes of the polymorphisms. The CRISP3 AJ459965:c.+622G>A SNP leading to the amino acid substitution E208K was significantly associated with the fertility of stallions. Stallions heterozygous for the CRISP3 c.+622G>A SNP had lower fertility than homozygous stallions (P = 0.0234). The pregnancy rate per cycle in these stallions was estimated to be approximately 7% lower than in stallions homozygous at this position.

Polymorphism and chromosomal location of the MC4R (melanocortin-4 receptor) gene in the dog and red fox

The melanocortin-4 receptor (MC4R) is expressed in the hypothalamus and regulates energy intake and body weight. In silico screening of the canine chromosome 1 sequence and a comparison with the porcine MC4R sequence by BLAST were performed. The nucleotide sequence of the whole coding region and 3'- and 5'-flanking regions of the dog (1214 bp) and red fox (1177 bp) MC4R gene was established and high conservation of the nucleotide sequences was revealed (99%). Five sets of PCR primers were designed and a search for polymorphism was performed by the SSCP technique in a group of 31 dogs representing nineteen breeds and 35 farm red foxes. Sequencing of DNA fragments, representing the identified SSCP patterns, revealed three single nucleotide polymorphisms (including a missense one) in dogs and four silent SNPs in red foxes. An average SNP frequency was approx. 1/400 bp in the dog and 1/300 bp in the red fox. We mapped the MC4R gene by FISH to the canine chromosome 1 (CFA1q1.1) and to the red fox chromosome 5 (VVU5p1.2).

The HCP5 single-nucleotide polymorphism: a simple screening tool for prediction of hypersensitivity reaction to abacavir

The HLA-B 5701 allele is predictive of hypersensitivity reaction to abacavir, a response herein termed "ABC-HSR." This study of 1,103 individuals infected with human immunodeficiency virus assessed the usefulness of genotyping a HCP5 single-nucleotide polymorphism (SNP), rs2395029, in relation to ABC-HSR. In populations with European ancestry, rs2395029 is in linkage disequilibrium with HLA-B 5701. The HCP5 SNP was present in all 98 HLA-B 5701-positive individuals and was absent in 999 of 1005 HLA-B 5701-negative individuals. rs2395029 was overrepresented in 25 individuals with clinically likely ABC-HSR, compared with its frequency in 175 ABC-tolerant individuals (80% vs. 2%, respectively; P < .0001). Therefore, HCP5 genotyping could serve as a simple screening tool for ABC-HSR, particularly in settings where sequence-based HLA typing is not available.

Concept, design and implementation of a cardiovascular gene-centric 50 k SNP array for large-scale genomic association studies.

A wealth of genetic associations for cardiovascular and metabolic phenotypes in humans has been accumulating over the last decade, in particular a large number of loci derived from recent genome wide association studies (GWAS). True complex disease-associated loci often exert modest effects, so their delineation currently requires integration of diverse phenotypic data from large studies to ensure robust meta-analyses. We have designed a gene-centric 50 K single nucleotide polymorphism (SNP) array to assess potentially relevant loci across a range of cardiovascular, metabolic and inflammatory syndromes. The array utilizes a "cosmopolitan" tagging approach to capture the genetic diversity across approximately 2,000 loci in populations represented in the HapMap and SeattleSNPs projects. The array content is informed by GWAS of vascular and inflammatory disease, expression quantitative trait loci implicated in atherosclerosis, pathway based approaches and comprehensive literature searching. The custom flexibility of the array platform facilitated interrogation of loci at differing stringencies, according to a gene prioritization strategy that allows saturation of high priority loci with a greater density of markers than the existing GWAS tools, particularly in African HapMap samples. We also demonstrate that the IBC array can be used to complement GWAS, increasing coverage in high priority CVD-related loci across all major HapMap populations. DNA from over 200,000 extensively phenotyped individuals will be genotyped with this array with a significant portion of the generated data being released into the academic domain facilitating in silico replication attempts, analyses of rare variants and cross-cohort meta-analyses in diverse populations. These datasets will also facilitate more robust secondary analyses, such as explorations with alternative genetic models, epistasis and gene-environment interactions.

Tracking a tuberculosis outbreak over 21 years: strain-specific single nucleotide polymorphism-typing combined with targeted whole genome sequencing.

BACKGROUND  Whole genome sequencing (WGS) is increasingly used in molecular-epidemiological investigations of bacterial pathogens, despite cost- and time-intensive analyses. We combined strain-specific single nucleotide polymorphism (SNP)-typing and targeted WGS to investigate a tuberculosis cluster spanning 21 years in Bern, Switzerland. METHODS  Based on genome sequences of three historical outbreak Mycobacterium tuberculosis isolates, we developed a strain-specific SNP-typing assay to identify further cases. We screened 1,642 patient isolates, and performed WGS on all identified cluster isolates. We extracted SNPs to construct genomic networks. Clinical and social data were retrospectively collected. RESULTS  We identified 68 patients associated with the outbreak strain. Most were diagnosed in 1991-1995, but cases were observed until 2011. Two thirds belonged to the homeless and substance abuser milieu. Targeted WGS revealed 133 variable SNP positions among outbreak isolates. Genomic network analyses suggested a single origin of the outbreak, with subsequent division into three sub-clusters. Isolates from patients with confirmed epidemiological links differed by 0-11 SNPs. CONCLUSIONS  Strain-specific SNP-genotyping allowed rapid and inexpensive identification of M. tuberculosis outbreak isolates in a population-based strain collection. Subsequent targeted WGS provided detailed insights into transmission dynamics. This combined approach could be applied to track bacterial pathogens in real-time and at high resolution.

The Vacuolar H+-ATPase B1 Subunit Polymorphism p.E161K Associates with Impaired Urinary Acidification in Recurrent Stone Formers

Mutations in the vacuolar–type H+-ATPase B1 subunit gene ATP6V1B1 cause autosomal–recessive distal renal tubular acidosis (dRTA). We previously identified a single-nucleotide polymorphism (SNP) in the human B1 subunit (c.481G.A; p.E161K) that causes greatly diminished pump function in vitro. To investigate the effect of this SNP on urinary acidification, we conducted a genotype-phenotype analysis of recurrent stone formers in theDallas and Bern kidney stone registries. Of 555 patients examined, 32 (5.8%) were heterozygous for the p.E161K SNP, and the remaining 523 (94.2%) carried two wild–type alleles. After adjustment for sex, age, body mass index, and dietary acid and alkali intake, p.E161K SNP carriers had a nonsignificant tendency to higher urinary pH on a random diet (6.31 versus 6.09; P=0.09). Under an instructed low–Ca and low–Na diet, urinary pH was higher in p.E161K SNP carriers (6.56 versus 6.01; P,0.01). Kidney stones of p.E161K carriers were more likely to contain calcium phosphate than stones of wild-type patients. In acute NH4Cl loading, p.E161K carriers displayed a higher trough urinary pH (5.34 versus 4.89; P=0.01) than wild-type patients. Overall, 14.6% of wild-type patients and 52.4% of p.E161K carriers were unable to acidify their urine below pH 5.3 and thus, can be considered to have incomplete dRTA. In summary, our data indicate that recurrent stone formers with the vacuolar H+-ATPase B1 subunit p.E161K SNP exhibit a urinary acidification deficit with an increased prevalence of calcium phosphate– containing kidney stones. The burden of E161K heterozygosity may be a forme fruste of dRTA.

Dopaminergic denervation severity depends on COMT Val158Met polymorphism in Parkinson's disease.

BACKGROUND Catecholamine-O-methyl-tranferase (COMT) initiates dopamine degradation. Its activity is mainly determined by a single nucleotide polymorphism in the COMT gene (Val158Met, rs4680) separating high (Val/Val, COMT(HH)), intermediate (Val/Met, COMT(HL)) and low metabolizers (Met/Met, COMT(LL)). We investigated dopaminergic denervation in the striatum in PD patients according to COMT rs4680 genotype. METHODS Patients with idiopathic PD were assessed for motor severity (UPDRS-III rating scale in OFF-state), dopaminergic denervation using [123I]-FP-CIT SPECT imaging, and genotyped for the COMT rs4680 enzyme. [123I]-FP-CIT binding potential (BP) for each voxel was defined by the ratio of tracer-binding in the region of interest (striatum, caudate nucleus and putamen) to that in a region of non-specific activity. Genotyping was performed using TaqMan(®) SNP genotyping assay. We used a regression model to evaluate the effect of COMT genotype on the BP in the striatum and its sub-regions. RESULTS Genotype distribution was: 11 (27.5%) COMT(HH), 26 (65%) COMT(HL) and 3 (7.5%) COMT(LL). There were no significant differences in disease severity, treatments, or motor scores between genotypes. When adjusted to clinical severity, gender and age, low and intermediate metabolizers showed significantly higher rates of striatal denervation (COMT(HL+LL) BP = 1.32 ± 0.04) than high metabolizers (COMT(HH), BP = 1.6 ± 0.08; F(1.34) = 9.0, p = 0.005). Striatal sub-regions showed similar results. BP and UPDRS-III motor scores (r = 0.44, p = 0.04) (p < 0.001) were highly correlated. There was a gender effect, but no gender-genotype interaction. CONCLUSIONS Striatal denervation differs according to COMT-Val158Met polymorphism. COMT activity may play a role as a compensatory mechanism in PD motor symptoms.

A common polymorphism associated with antibiotic-induced cardiac arrhythmia

Drug-induced long QT syndrome (LQTS) is a prevalent disorder of uncertain etiology that predisposes to sudden death. KCNE2 encodes MinK-related peptide 1 (MiRP1), a subunit of the cardiac potassium channel IKr that has been associated previously with inherited LQTS. Here, we examine KCNE2 in 98 patients with drug-induced LQTS, identifying three individuals with sporadic mutations and a patient with sulfamethoxazole-associated LQTS who carried a single-nucleotide polymorphism (SNP) found in ≈1.6% of the general population. While mutant channels showed diminished potassium flux at baseline and wild-type drug sensitivity, channels with the SNP were normal at baseline but inhibited by sulfamethoxazole at therapeutic levels that did not affect wild-type channels. We conclude that allelic variants of MiRP1 contribute to a significant fraction of cases of drug-induced LQTS through multiple mechanisms and that common sequence variations that increase the risk of life-threatening drug reactions can be clinically silent before drug exposure.

Alkaline-mediated differential interaction (AMDI): A simple automatable single-nucleotide polymorphism assay

The key requirements for high-throughput single-nucleotide polymorphism (SNP) typing of DNA samples in large-scale disease case-control studies are automatability, simplicity, and robustness, coupled with minimal cost. In this paper we describe a fluorescence technique for the detection of SNPs that have been amplified by using the amplification refractory mutation system (ARMS)-PCR procedure. Its performance was evaluated using 32 sequence-specific primer mixes to assign the HLA-DRB alleles to 80 lymphoblastoid cell line DNAs chosen from our database for their diversity. All had been typed previously by alternative methods, either direct sequencing or gel electrophoresis. We believe the detection system that we call AMDI (alkaline-mediated differential interaction) satisfies the above criteria and is suitable for general high-throughput SNP typing.

Epidermal Growth Factor Gene (EGF) Polymorphism and Risk of Melanocytic Neoplasiay

A common single nucleotide polymorphism (SNP) in the 5' untranslated region (5'UTR) of the epidermal growth factor (EGF) gene modulates the level of transcription of this gene and hence is associated with serum levels of EGF. This variant may be associated with melanoma risk, but conflicting findings have been reported. An Australian melanoma case-control sample was typed for the EGF+61A>G transversion (rs4444903). The sample comprised 753 melanoma cases from 738 families stratified by family history of melanoma and 2387 controls from 645 unselected twin families. Ancestry of the cases and controls was recorded, and the twins had undergone skin examination to assess total body nevus count, degree of freckling and pigmentation phenotype. SNP genotyping was carried out via primer extension followed by matrix-assisted laser desorption time of flight (MALDI-TOF) mass spectroscopy. The EGIF+61 SNP was not found to be significantly associated with melanoma status or with development of nevi or freckles. Among melanoma cases, however, G homozygotes had thicker tumors (p=0.05), in keeping with two previous studies. The EGF polymorphism does not appear to predispose to melanoma or nevus development, but its significant association with tumor thickness implies that it may be a useful marker of prognosis.

Development of an oligonucleotide-based SNP detection method on lateral flow strips using hexapet tages

Detection of point mutations or single nucleotide polymorphisms (SNPs) is important in relation to disease susceptibility or detection in pathogens of mutations determining drug resistance or host range. There is an emergent need for rapid detection methods amenable to point-of-care applications. The purpose of this study was to reduce to practice a novel method for SNP detection and to demonstrate that this technology can be used downstream of nucleic acid amplification. The authors used a model system to develop an oligonucleotide-based SNP detection system on nitrocellulose lateral flow strips. To optimize the assay they used cloned sequences of the herpes simplex virus-1 (HSV-1) DNA polymerase gene into which they introduced a point mutation. The assay system uses chimeric polymerase chain reaction (PCR) primers that incorporate hexameric repeat tags ("hexapet tags"). The chimeric sequences allow capture of amplified products to predefined positions on a lateral flow strip. These "hexapet" sequences have minimal cross-reactivity and allow specific hybridization-based capture of the PCR products at room temperature onto lateral flow strips that have been striped with complementary hexapet tags. The allele-specific amplification was carried out with both mutant and wild-type primer sets present in the PCR mix ("competitive" format). The resulting PCR products carried a hexapet tag that corresponded with either a wild-type or mutant sequence. The lateral flow strips are dropped into the PCR reaction tube, and mutant sequence and wild-type sequences diffuse along the strip and are captured at the corresponding position on the strip. A red line indicative of a positive reaction is visible after 1 minute. Unlike other systems that require separate reactions and strips for each target sequence, this system allows multiplex PCR reactions and multiplex detection on a single strip or other suitable substrates. Unambiguous visual discrimination of a point mutation under room temperature hybridization conditions was achieved with this model system in 10 minutes after PCR. The authors have developed a capture-based hybridization method for the detection and discrimination of HSV-1 DNA polymerase genes that contain a single nucleotide change. It has been demonstrated that the hexapet oligonucleotides can be adapted for hybridization on the lateral flow strip platform for discrimination of SNPs. This is the first step in demonstrating SNP detection on lateral flow using the hexapet oligonucleotide capture system. It is anticipated that this novel system can be widely used in point-of-care settings.

The promoter polymorphism -232C/G of the PCK1 gene is associated with type 2 diabetes in a UK-resident South Asian population

Background - The PCK1 gene, encoding cytosolic phosphoenolpyruvate carboxykinase (PEPCKC), has previously been implicated as a candidate gene for type 2 diabetes (T2D) susceptibility. Rodent models demonstrate that over-expression of Pck1 can result in T2D development and a single nucleotide polymorphism (SNP) in the promoter region of human PCK1 (-232C/G) has exhibited significant association with the disease in several cohorts. Within the UK-resident South Asian population, T2D is 4 to 6 times more common than in indigenous white Caucasians. Despite this, few studies have reported on the genetic susceptibility to T2D in this ethnic group and none of these has investigated the possible effect of PCK1 variants. We therefore aimed to investigate the association between common variants of the PCK1 gene and T2D in a UK-resident South Asian population of Punjabi ancestry, originating predominantly from the Mirpur area of Azad Kashmir, Pakistan. Methods - We used TaqMan assays to genotype five tagSNPs covering the PCK1 gene, including the -232C/G variant, in 903 subjects with T2D and 471 normoglycaemic controls. Results - Of the variants studied, only the minor allele (G) of the -232C/G SNP demonstrated a significant association with T2D, displaying an OR of 1.21 (95% CI: 1.03 - 1.42, p = 0.019). Conclusion - This study is the first to investigate the association between variants of the PCK1 gene and T2D in South Asians. Our results suggest that the -232C/G promoter polymorphism confers susceptibility to T2D in this ethnic group.

Integration of global SNP-based mapping and expression arrays reveals key regions, mechanisms, and genes important in the pathogenesis of multiple myeloma.

Multiple myeloma is characterized by genomic alterations frequently involving gains and losses of chromosomes. Single nucleotide polymorphism (SNP)-based mapping arrays allow the identification of copy number changes at the sub-megabase level and the identification of loss of heterozygosity (LOH) due to monosomy and uniparental disomy (UPD). We have found that SNP-based mapping array data and fluorescence in situ hybridization (FISH) copy number data correlated well, making the technique robust as a tool to investigate myeloma genomics. The most frequently identified alterations are located at 1p, 1q, 6q, 8p, 13, and 16q. LOH is found in these large regions and also in smaller regions throughout the genome with a median size of 1 Mb. We have identified that UPD is prevalent in myeloma and occurs through a number of mechanisms including mitotic nondisjunction and mitotic recombination. For the first time in myeloma, integration of mapping and expression data has allowed us to reduce the complexity of standard gene expression data and identify candidate genes important in both the transition from normal to monoclonal gammopathy of unknown significance (MGUS) to myeloma and in different subgroups within myeloma. We have documented these genes, providing a focus for further studies to identify and characterize those that are key in the pathogenesis of myeloma.