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.

Optimization of sample preparation for next-generation sequencing with Illumina Genome Analyzer II

The growing interest for sequencing with higher throughput in the last decade has led to the development of new sequencing applications. This thesis concentrates on optimizing DNA library preparation for Illumina Genome Analyzer II sequencer. The library preparation steps that were optimized include fragmentation, PCR purification and quantification. DNA fragmentation was performed with focused sonication in different concentrations and durations. Two column based PCR purification method, gel matrix method and magnetic bead based method were compared. Quantitative PCR and gel electrophoresis in a chip were compared for DNA quantification. The magnetic bead purification was found to be the most efficient and flexible purification method. The fragmentation protocol was changed to produce longer fragments to be compatible with longer sequencing reads. Quantitative PCR correlates better with the cluster number and should thus be considered to be the default quantification method for sequencing. As a result of this study more data have been acquired from sequencing with lower costs and troubleshooting has become easier as qualification steps have been added to the protocol. New sequencing instruments and applications will create a demand for further optimizations in future.

Development and characterization of twenty-two polymorphic microsatellite markers for the leafcutter ant, Acromyrmex lundii, utilizing Illumina sequencing

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Development of twenty-one polymorphic microsatellite markers for the fungus-growing ant, Mycocepurus goeldii (Formicidae: Attini), using Illumina paired-end genomic sequencing

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Water buffalo genome characterization by the Illumina BovineHD BeadChip

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Analise microbiológica da flora intestinal do mosquito Aedes Aegypti através do sequenciamento do gene 16S rRNA na plataforma Illumina MiSeq

Dengue é uma arbovirose que afeta cerca de 100 milhões de pessoas anualmente, em mais de 100 países situados nas regiões tropicais e subtropicais. Foi considerada a doença viral que mais cresceu no ultimo ano, repercutindo em impactos sociais e econômicos nas regiões endêmicas devido às altas taxas de morbidade e mortalidade desencadeadas pela infecção. O principal vetor da dengue é o mosquito Aedes aegypti, presente em toda a faixa tropical e subtropical. Por apresentar hematofagia antropofílica, rápido desenvolvimento e características comportamentais especificas, é um excelente transmissor do vírus dengue. Medidas de controle da disseminação da dengue são restritas à eliminação do mosquito vetor, e um tratamento específico ainda não foi desenvolvido, bem como a criação de uma vacina que previna simultaneamente a infecção pelos quatro sorotipos do arbovírus. Uma característica que determina a disseminação de doenças é a alta competência vetorial de seus mosquitos transmissores, que tem sido associada à composição da microbiota intestinal do inseto. As bactérias presente no intestino do mosquito exercem funções relacionadas a sua nutrição, desenvolvimento e reprodução, e são também um importante fator na eliminação de patógenos, por interferirem diretamente na atividade viral, ou indiretamente a partir da ativação das vias antivirais pelos micro-organismos. Dessa forma, este trabalho visa estudar a diversidade microbiana intestinal do mosquito Aedes aegypti em diferentes estágios de vida, através de sequenciamento de última geração com a plataforma MiSeq Illumina

Comparative clinical sample preparation of DNA and RNA viral nucleic acids for a commercial deep sequencing system (Illumina MiSeq (R))

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Generation and analysis of a mouse intestinal metatranscriptome through Illumina based RNA-sequencing

With the advent of high through-put sequencing (HTS), the emerging science of metagenomics is transforming our understanding of the relationships of microbial communities with their environments. While metagenomics aims to catalogue the genes present in a sample through assessing which genes are actively expressed, metatranscriptomics can provide a mechanistic understanding of community inter-relationships. To achieve these goals, several challenges need to be addressed from sample preparation to sequence processing, statistical analysis and functional annotation. Here we use an inbred non-obese diabetic (NOD) mouse model in which germ-free animals were colonized with a defined mixture of eight commensal bacteria, to explore methods of RNA extraction and to develop a pipeline for the generation and analysis of metatranscriptomic data. Applying the Illumina HTS platform, we sequenced 12 NOD cecal samples prepared using multiple RNA-extraction protocols. The absence of a complete set of reference genomes necessitated a peptide-based search strategy. Up to 16% of sequence reads could be matched to a known bacterial gene. Phylogenetic analysis of the mapped ORFs revealed a distribution consistent with ribosomal RNA, the majority from Bacteroides or Clostridium species. To place these HTS data within a systems context, we mapped the relative abundance of corresponding Escherichia coli homologs onto metabolic and protein-protein interaction networks. These maps identified bacterial processes with components that were well-represented in the datasets. In summary this study highlights the potential of exploiting the economy of HTS platforms for metatranscriptomics.

Isolation and Characterization of Microsatellite DNA Markers in the Deep-Sea Amphipod Paralicella tenuipes by Illumina MiSeq Sequencing

Funding This work was supported by the HADEEP projects, funded by the Nippon Foundation, Japan (2009765188), the Natural Environmental Research Council, UK (NE/E007171/1) and the Total Foundation, France. We acknowledge additional support from the Marine Alliance for Science and Technology for Scotland (MASTS) funded by the Scottish Funding Council (Ref: HR09011) and contributing institutions. We also acknowledge support from the Leverhulme Trust to SBP. Additional sea time was supported by NIWA’s ‘Impact of Resource Use on Vulnerable Deep-Sea Communities’ project (CO1_0906)

Mapping Protein-DNA Interactions Using ChIP-exo and Illumina-Based Sequencing

Chromatin immunoprecipitation (ChIP) provides a means of enriching DNA associated with transcription factors, histone modifications, and indeed any other proteins for which suitably characterized antibodies are available. Over the years, sequence detection has progressed from quantitative real-time PCR and Southern blotting to microarrays (ChIP-chip) and now high-throughput sequencing (ChIP-seq). This progression has vastly increased the sequence coverage and data volumes generated. This in turn has enabled informaticians to predict the identity of multi-protein complexes on DNA based on the overrepresentation of sequence motifs in DNA enriched by ChIP with a single antibody against a single protein. In the course of the development of high-throughput sequencing, little has changed in the ChIP methodology until recently. In the last three years, a number of modifications have been made to the ChIP protocol with the goal of enhancing the sensitivity of the method and further reducing the levels of nonspecific background sequences in ChIPped samples. In this chapter, we provide a brief commentary on these methodological changes and describe a detailed ChIP-exo method able to generate narrower peaks and greater peak coverage from ChIPped material.

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.