Comparing genotyping algorithms for Illumina's Infinium whole-genome SNP BeadChips

Background Illumina's Infinium SNP BeadChips are extensively used in both small and large-scale genetic studies. A fundamental step in any analysis is the processing of raw allele A and allele B intensities from each SNP into genotype calls (AA, AB, BB). Various algorithms which make use of different statistical models are available for this task. We compare four methods (GenCall, Illuminus, GenoSNP and CRLMM) on data where the true genotypes are known in advance and data from a recently published genome-wide association study. Results In general, differences in accuracy are relatively small between the methods evaluated, although CRLMM and GenoSNP were found to consistently outperform GenCall. The performance of Illuminus is heavily dependent on sample size, with lower no call rates and improved accuracy as the number of samples available increases. For X chromosome SNPs, methods with sex-dependent models (Illuminus, CRLMM) perform better than methods which ignore gender information (GenCall, GenoSNP). We observe that CRLMM and GenoSNP are more accurate at calling SNPs with low minor allele frequency than GenCall or Illuminus. The sample quality metrics from each of the four methods were found to have a high level of agreement at flagging samples with unusual signal characteristics. Conclusions CRLMM, GenoSNP and GenCall can be applied with confidence in studies of any size, as their performance was shown to be invariant to the number of samples available. Illuminus on the other hand requires a larger number of samples to achieve comparable levels of accuracy and its use in smaller studies (50 or fewer individuals) is not recommended.

The waddlia genome: A window into chlamydial biology

Growing evidence suggests that a novel member of the Chlamydiales order, Waddlia chondrophila, is a potential agent of miscarriage in humans and abortion in ruminants. Due to the lack of genetic tools to manipulate chlamydia, genomic analysis is proving to be the most incisive tool in stimulating investigations into the biology of these obligate intracellular bacteria. 454/Roche and Solexa/Illumina technologies were thus used to sequence and assemble de novo the full genome of the first representative of the Waddliaceae family, W. chondrophila. The bacteria possesses a 2′116′312bp chromosome and a 15′593 bp low-copy number plasmid that might integrate into the bacterial chromosome. The Waddlia genome displays numerous repeated sequences indicating different genome dynamics from classical chlamydia which almost completely lack repetitive elements. Moreover, W. chondrophila exhibits many virulence factors also present in classical chlamydia, including a functional type III secretion system, but also a large complement of specific factors for resistance to host or environmental stresses. Large families of outer membrane proteins were identified indicating that these highly immunogenic proteins are not Chlamydiaceae specific and might have been present in their last common ancestor. Enhanced metabolic capability for the synthesis of nucleotides, amino acids, lipids and other co-factors suggests that the common ancestor of the modern Chlamydiales may have been less dependent on their eukaryotic host. The fine-detailed analysis of biosynthetic pathways brings us closer to possibly developing a synthetic medium to grow W. chondrophila, a critical step in the development of genetic tools. As a whole, the availability of the W. chondrophila genome opens new possibilities in Chlamydiales research, providing new insights into the evolution of members of the order Chlamydiales and the biology of the Waddliaceae.

A systematic review of environmental factors and obesogenic dietary intakes among adults : are we getting closer to understanding obesogenic environments?

This study examined whether physical, social, cultural and economical environmental factors are associated with obesogenic dietary behaviours and overweight/obesity among adults. Literature searches of databases (i.e. PubMed, CSA Illumina, Web of Science, PsychInfo) identified studies examining environmental factors and the consumption of energy, fat, fibre, fruit, vegetables, sugar-sweetened drinks, meal patterns and weight status. Twenty-eight studies were in-scope, the majority (n= 16) were conducted in the USA. Weight status was consistently associated with the food environment; greater accessibility to supermarkets or less access to takeaway outlets were associated with a lower BMI or prevalence of overweight/obesity. However, obesogenic dietary behaviours did not mirror these associations; mixed associations were found between the environment and obesogenic dietary behaviours. Living in a socioeconomically-deprived area was the only environmental factor consistently associated with a number of obesogenic dietary behaviours. Associations between the environment and weight status are more consistent than that seen between the environment and dietary behaviours. The environment may play an important role in the development of overweight/obesity, however the dietary mechanisms that contribute to this remain unclear and the physical activity environment may also play an important role in weight gain, overweight and obesity.

Next-generation sequencing of cervical DNA detects human papillomavirus types not detected by commercial kits

Background Human papillomavirus (HPV) is the aetiological agent for cervical cancer and genital warts. Concurrent HPV and HIV infection in the South African population is high. HIV positive (+) women are often infected with multiple, rare and undetermined HPV types. Data on HPV incidence and genotype distribution are based on commercial HPV detection kits, but these kits may not detect all HPV types in HIV + women. The objectives of this study were to (i) identify the HPV types not detected by commercial genotyping kits present in a cervical specimen from an HIV positive South African woman using next generation sequencing, and (ii) determine if these types were prevalent in a cohort of HIV-infected South African women. Methods Total DNA was isolated from 109 cervical specimens from South African HIV + women. A specimen within this cohort representing a complex multiple HPV infection, with 12 HPV genotypes detected by the Roche Linear Array HPV genotyping (LA) kit, was selected for next generation sequencing analysis. All HPV types present in this cervical specimen were identified by Illumina sequencing of the extracted DNA following rolling circle amplification. The prevalence of the HPV types identified by sequencing, but not included in the Roche LA, was then determined in the 109 HIV positive South African women by type-specific PCR. Results Illumina sequencing identified a total of 16 HPV genotypes in the selected specimen, with four genotypes (HPV-30, 74, 86 and 90) not included in the commercial kit. The prevalence's of HPV-30, 74, 86 and 90 in 109 HIV positive South African women were found to be 14.6 %, 12.8 %, 4.6 % and 8.3 % respectively. Conclusions Our results indicate that there are HPV types, with substantial prevalence, in HIV positive women not being detected in molecular epidemiology studies using commercial kits. The significance of these types in relation to cervical disease remains to be investigated.

Identification of 23 new prostate cancer susceptibility loci using the iCOGS custom genotyping

Prostate cancer is the most frequently diagnosed cancer in males in developed countries. To identify common prostate cancer susceptibility alleles, we genotyped 211,155 SNPs on a custom Illumina array (iCOGS) in blood DNA from 25,074 prostate cancer cases and 24,272 controls from the international PRACTICAL Consortium. Twenty-three new prostate cancer susceptibility loci were identified at genome-wide significance (P < 5 × 10−8). More than 70 prostate cancer susceptibility loci, explaining ~30% of the familial risk for this disease, have now been identified. On the basis of combined risks conferred by the new and previously known risk loci, the top 1% of the risk distribution has a 4.7-fold higher risk than the average of the population being profiled. These results will facilitate population risk stratification for clinical studies.

Fine-mapping identifies multiple prostate cancer risk loci on 5p15 some associating with TERT expression

Associations between single nucleotide polymorphisms (SNPs) at 5p15 and multiple cancer types have been reported. We have previously shown evidence for a strong association between prostate cancer (PrCa) risk and rs2242652 at 5p15, intronic in the telomerase reverse transcriptase (TERT) gene that encodes TERT. To comprehensively evaluate the association between genetic variation across this region and PrCa, we performed a fine-mapping analysis by genotyping 134 SNPs using a custom Illumina iSelect array or Sequenom MassArray iPlex, followed by imputation of 1094 SNPs in 22 301 PrCa cases and 22 320 controls in The PRACTICAL consortium. Multiple stepwise logistic regression analysis identified four signals in the promoter or intronic regions of TERT that independently associated with PrCa risk. Gene expression analysis of normal prostate tissue showed evidence that SNPs within one of these regions also associated with TERT expression, providing a potential mechanism for predisposition to disease.

Transcriptome analysis of neutrophils after endurance exercise reveals novel signaling mechanisms in the immune response to physiological stress

Neutrophils serve as an intriguing model for the study of innate immune cellular activity induced by physiological stress. We measured changes in the transcriptome of circulating neutrophils following an experimental exercise trial (EXTRI) consisting of 1 h of intense cycling immediately followed by 1 h of intense running. Blood samples were taken at baseline, 3 h, 48 h, and 96 h post-EXTRI from eight healthy, endurance-trained, male subjects. RNA was extracted from isolated neutrophils. Differential gene expression was evaluated using Illumina microarrays and validated with quantitative PCR. Gene set enrichment analysis identified enriched molecular signatures chosen from the Molecular Signatures Database. Blood concentrations of muscle damage indexes, neutrophils, interleukin (IL)-6 and IL-10 were increased (P < 0.05) 3 h post-EXTRI. Upregulated groups of functionally related genes 3 h post-EXTRI included gene sets associated with the recognition of tissue damage, the IL-1 receptor, and Toll-like receptor (TLR) pathways (familywise error rate, P value < 0.05). The core enrichment for these pathways included TLRs, low-affinity immunoglobulin receptors, S100 calcium binding protein A12, and negative regulators of innate immunity, e.g., IL-1 receptor antagonist, and IL-1 receptor associated kinase-3. Plasma myoglobin changes correlated with neutrophil TLR4 gene expression (r = 0.74; P < 0.05). Neutrophils had returned to their nonactivated state 48 h post-EXTRI, indicating that their initial proinflammatory response was transient and rapidly counterregulated. This study provides novel insight into the signaling mechanisms underlying the neutrophil responses to endurance exercise, suggesting that their transcriptional activity was particularly induced by damage-associated molecule patterns, hypothetically originating from the leakage of muscle components into the circulation.

Variants in the human potassium channel gene (KCNN3) are associated with migraine in a high risk genetic isolate

The calcium-activated potassium ion channel gene (KCNN3) is located in the vicinity of the familial hemiplegic migraine type 2 locus on chromosome 1q21.3. This gene is expressed in the central nervous system and plays a role in neural excitability. Previous association studies have provided some, although not conclusive, evidence for involvement of this gene in migraine susceptibility. To elucidate KCNN3 involvement in migraine, we performed gene-wide SNP genotyping in a high-risk genetic isolate from Norfolk Island, a population descended from a small number of eighteenth century Isle of Man ‘Bounty Mutineer’ and Tahitian founders. Phenotype information was available for 377 individuals who are related through the single, well-defined Norfolk pedigree (96 were affected: 64 MA, 32 MO). A total of 85 SNPs spanning the KCNN3 gene were genotyped in a sub-sample of 285 related individuals (76 affected), all core members of the extensive Norfolk Island ‘Bounty Mutineer’ genealogy. All genotyping was performed using the Illumina BeadArray platform. The analysis was performed using the statistical program SOLAR v4.0.6 assuming an additive model of allelic effect adjusted for the effects of age and sex. Haplotype analysis was undertaken using the program HAPLOVIEW v4.0. A total of four intronic SNPs in the KCNN3 gene displayed significant association (P < 0.05) with migraine. Two SNPs, rs73532286 and rs6426929, separated by approximately 0.1 kb, displayed complete LD (r 2 = 1.00, D′ = 1.00, D′ 95% CI = 0.96–1.00). In all cases, the minor allele led to a decrease in migraine risk (beta coefficient = 0.286–0.315), suggesting that common gene variants confer an increased risk of migraine in the Norfolk pedigree. This effect may be explained by founder effect in this genetic isolate. This study provides evidence for association of variants in the KCNN3 ion channel gene with migraine susceptibility in the Norfolk genetic isolate with the rarer allelic variants conferring a possible protective role. This the first comprehensive analysis of this potential candidate gene in migraine and also the first study that has utilised the unique Norfolk Island large pedigree isolate to implicate a specific migraine gene. Studies of additional variants in KCNN3 in the Norfolk pedigree are now required (e.g. polyglutamine variants) and further analyses in other population data sets are required to clarify the association of the KCNN3 gene and migraine risk in the general outbred population.

Saliva-derived DNA performs well in large-scale, high-density single-nucleotide polymorphism microarray studies

As of June 2009, 361 genome-wide association studies (GWAS) had been referenced by the HuGE database. GWAS require DNA from many thousands of individuals, relying on suitable DNA collections. We recently performed a multiple sclerosis (MS) GWAS where a substantial component of the cases (24%) had DNA derived from saliva. Genotyping was done on the Illumina genotyping platform using the Infinium Hap370CNV DUO microarray. Additionally, we genotyped 10 individuals in duplicate using both saliva- and blood-derived DNA. The performance of blood- versus saliva-derived DNA was compared using genotyping call rate, which reflects both the quantity and quality of genotyping per sample and the “GCScore,” an Illumina genotyping quality score, which is a measure of DNA quality. We also compared genotype calls and GCScores for the 10 sample pairs. Call rates were assessed for each sample individually. For the GWAS samples, we compared data according to source of DNA and center of origin. We observed high concordance in genotyping quality and quantity between the paired samples and minimal loss of quality and quantity of DNA in the saliva samples in the large GWAS sample, with the blood samples showing greater variation between centers of origin. This large data set highlights the usefulness of saliva DNA for genotyping, especially in high-density single-nucleotide polymorphism microarray studies such as GWAS.

Identification of TFG (TRK-fused gene) as a putative metastatic melanoma tumor suppressor gene

High density SNP arrays can be used to identify DNA copy number changes in tumors such as homozygous deletions of tumor suppressor genes and focal amplifications of oncogenes. Illumina Human CNV370 Bead chip arrays were used to assess the genome for unbalanced chromosomal events occurring in 39 cell lines derived from stage III metastatic melanomas. A number of genes previously recognized to have an important role in the development and progression of melanoma were identified including homozygous deletions of CDKN2A (13 of 39 samples), CDKN2B (10 of 39), PTEN (3 of 39), PTPRD (3 of 39), TP53 (1 of 39), and amplifications of CCND1 (2 of 39), MITF (2 of 39), MDM2 (1 of 39), and NRAS (1 of 39). In addition, a number of focal homozygous deletions potentially targeting novel melanoma tumor suppressor genes were identified. Because of their likely functional significance for melanoma progression, FAS, CH25H, BMPR1A, ACTA2, and TFG were investigated in a larger cohort of melanomas through sequencing. Nonsynonymous mutations were identified in BMPR1A (1 of 43), ACTA2 (3 of 43), and TFG (5 of 103). A number of potentially important mutation events occurred in TFG including the identification of a mini mutation ‘‘hotspot’’ at amino acid residue 380 (P380S and P380L) and the presence of multiple mutations in two melanomas. Mutations in TFG may have important clinical relevance for current therapeutic strategies to treat metastatic melanoma.

Time course-dependent changes in the transcriptome of human skeletal muscle during recovery from endurance exercise: From inflammation to adaptive remodeling

Re-programming of gene expression is fundamental for skeletal muscle adaptations in response to endurance exercise. This study investigated the time-course dependent changes in the muscular transcriptome following an endurance exercise trial consisting of 1 h of intense cycling immediately followed by 1 h of intense running. Skeletal muscle samples were taken at baseline, 3 h, 48 h, and 96 h post-exercise from eight healthy, endurance-trained, male individuals. RNA was extracted from muscle. Differential gene expression was evaluated using Illumina microarrays and validated with qPCR. Gene set enrichment analysis identified enriched molecular signatures chosen from the Molecular Signatures Database. Three h post-exercise, 102 gene sets were up-regulated [family wise error rate (FWER), P < 0.05]; including groups of genes related with leukocyte migration, immune and chaperone activation, and cyclic AMP responsive element binding protein (CREB) 1-signaling. Forty-eight h post-exercise, among 19 enriched gene sets (FWER, P < 0.05), two gene sets related to actin cytoskeleton remodeling were up-regulated. Ninety-six h post-exercise, 83 gene sets were enriched (FWER, P < 0.05), 80 of which were up-regulated; including gene groups related to chemokine signaling, cell stress management, and extracellular matrix remodeling. These data provide comprehensive insights into the molecular pathways involved in acute stress, recovery, and adaptive muscular responses to endurance exercise. The novel 96 h post-exercise transcriptome indicates substantial transcriptional activity, potentially associated with the prolonged presence of leukocytes in the muscles. This suggests that muscular recovery, from a transcriptional perspective, is incomplete 96 h after endurance exercise involving muscle damage.

Genome-wide DNA methylation analysis of human brain tissue from schizophrenia patients

Recent studies suggest that genetic and environmental factors do not account for all the schizophrenia risk and epigenetics also plays a role in disease susceptibility. DNA methylation is a heritable epigenetic modification that can regulate gene expression. Genome-Wide DNA methylation analysis was performed on post-mortem human brain tissue from 24 patients with schizophrenia and 24 unaffected controls. DNA methylation was assessed at over 485 000 CpG sites using the Illumina Infinium Human Methylation450 Bead Chip. After adjusting for age and post-mortem interval (PMI), 4 641 probes corresponding to 2 929 unique genes were found to be differentially methylated. Of those genes, 1 291 were located in a CpG island and 817 were in a promoter region. These include NOS1, AKT1, DTNBP1, DNMT1, PPP3CC and SOX10 which have previously been associated with schizophrenia. More than 100 of these genes overlap with a previous DNA methylation study of peripheral blood from schizophrenia patients in which 27 000 CpG sites were analysed. Unsupervised clustering analysis of the top 3 000 most variable probes revealed two distinct groups with significantly more people with schizophrenia in cluster one compared to controls (p = 1.74x10-4). The first cluster was composed of 88% of patients with schizophrenia and only 12% controls while the second cluster was composed of 27% of patients with schizophrenia and 73% controls. These results strongly suggest that differential DNA methylation is important in schizophrenia etiology and add support for the use of DNA methylation profiles as a future prognostic indicator of schizophrenia.

Progesterone activates multiple innate immune pathways inChlamydia trachomatis-Infected endocervical cells

Problem Susceptibility to Chlamydia trachomatis infection is increased by oral con- traceptives and modulated by sex hormones. We therefore sought to determine the effects of female sex hormones on the innate immune response to C. trachomatis infection. Method of study ECC-1 endometrial cells, pre-treated with oestradiol or progesterone, were infected with C. trachomatis and the host transcriptome analysed by Illumina Sentrix HumanRef-8 microarray. Primary endocervical epithe- lial cells, prepared at either the proliferative or secretory phase of the menstrual cycle, were infected with C. trachomatis and cytokine gene expression determined by quantitative RT-PCR analysis. Results Chlamydia trachomatis yield from progesterone-primed ECC-1 cells was significantly reduced compared with oestradiol-treated cells. Genes upregulated in progesterone-treated and Chlamydia-infected cells only included multiple CC and CXC chemokines, IL-17C, IL-29, IL-32, TNF-a, DEFB4B, LCN2, S100A7-9, ITGAM, NOD2, JAK1, IL-6ST, type I and II interferon receptors, numerous interferon-stimulated genes and STAT6. CXCL10, CXCL11, CX3CL1 and IL-17C, which were also upregu- lated in infected secretory-stage primary cells, and there was a trend towards higher levels of immune mediators in infected secretory-phase compared with proliferative-phase cells. Conclusion Progesterone treatment primes multiple innate immune pathways in hormone-responsive epithelial cells that could potentially increase resis- tance to chlamydial infection.

Genome-wide SNP detection, validation, and development of an 8K SNP array for apple

As high-throughput genetic marker screening systems are essential for a range of genetics studies and plant breeding applications, the International RosBREED SNP Consortium (IRSC) has utilized the Illumina Infinium® II system to develop a medium- to high-throughput SNP screening tool for genome-wide evaluation of allelic variation in apple (Malus×domestica) breeding germplasm. For genome-wide SNP discovery, 27 apple cultivars were chosen to represent worldwide breeding germplasm and re-sequenced at low coverage with the Illumina Genome Analyzer II. Following alignment of these sequences to the whole genome sequence of 'Golden Delicious', SNPs were identified using SoapSNP. A total of 2,113,120 SNPs were detected, corresponding to one SNP to every 288 bp of the genome. The Illumina GoldenGate® assay was then used to validate a subset of 144 SNPs with a range of characteristics, using a set of 160 apple accessions. This validation assay enabled fine-tuning of the final subset of SNPs for the Illumina Infinium® II system. The set of stringent filtering criteria developed allowed choice of a set of SNPs that not only exhibited an even distribution across the apple genome and a range of minor allele frequencies to ensure utility across germplasm, but also were located in putative exonic regions to maximize genotyping success rate. A total of 7867 apple SNPs was established for the IRSC apple 8K SNP array v1, of which 5554 were polymorphic after evaluation in segregating families and a germplasm collection. This publicly available genomics resource will provide an unprecedented resolution of SNP haplotypes, which will enable marker-locus-trait association discovery, description of the genetic architecture of quantitative traits, investigation of genetic variation (neutral and functional), and genomic selection in apple.

Fine-mapping the HOXB region detects common variants tagging a rare coding allele : evidence for synthetic association in prostate cancer

The HOXB13 gene has been implicated in prostate cancer (PrCa) susceptibility. We performed a high resolution fine-mapping analysis to comprehensively evaluate the association between common genetic variation across the HOXB genetic locus at 17q21 and PrCa risk. This involved genotyping 700 SNPs using a custom Illumina iSelect array (iCOGS) followed by imputation of 3195 SNPs in 20,440 PrCa cases and 21,469 controls in The PRACTICAL consortium. We identified a cluster of highly correlated common variants situated within or closely upstream of HOXB13 that were significantly associated with PrCa risk, described by rs117576373 (OR 1.30, P = 2.62×10(-14)). Additional genotyping, conditional regression and haplotype analyses indicated that the newly identified common variants tag a rare, partially correlated coding variant in the HOXB13 gene (G84E, rs138213197), which has been identified recently as a moderate penetrance PrCa susceptibility allele. The potential for GWAS associations detected through common SNPs to be driven by rare causal variants with higher relative risks has long been proposed; however, to our knowledge this is the first experimental evidence for this phenomenon of synthetic association contributing to cancer susceptibility.