SNP discovery using Next Generation Transcriptomic Sequencing in Atlantic herring (Clupea harengus)

The introduction of Next Generation Sequencing (NGS) has revolutionised population genetics, providing studies of non-model species with unprecedented genomic coverage, allowing evolutionary biologists to address questions previously far beyond the reach of available resources. Furthermore, the simple mutation model of Single Nucleotide Polymorphisms (SNPs) permits cost-effective high-throughput genotyping in thousands of individuals simultaneously. Genomic resources are scarce for the Atlantic herring (Clupea harengus), a small pelagic species that sustains high revenue fisheries. This paper details the development of 578 SNPs using a combined NGS and high-throughput genotyping approach. Eight individuals covering the species distribution in the eastern Atlantic were bar-coded and multiplexed into a single cDNA library and sequenced using the 454 GS FLX platform. SNP discovery was performed by de novo sequence clustering and contig assembly, followed by the mapping of reads against consensus contig sequences. Selection of candidate SNPs for genotyping was conducted using an in silico approach. SNP validation and genotyping were performed simultaneously using an Illumina 1,536 GoldenGate assay. Although the conversion rate of candidate SNPs in the genotyping assay cannot be predicted in advance, this approach has the potential to maximise cost and time efficiencies by avoiding expensive and time-consuming laboratory stages of SNP validation. Additionally, the in silico approach leads to lower ascertainment bias in the resulting SNP panel as marker selection is based only on the ability to design primers and the predicted presence of intron-exon boundaries. Consequently SNPs with a wider spectrum of minor allele frequencies (MAFs) will be genotyped in the final panel. The genomic resources presented here represent a valuable multi-purpose resource for developing informative marker panels for population discrimination, microarray development and for population genomic studies in the wild.

Novel tools for conservation genomics: comparing two high-throughput approaches for SNP discovery in the transcriptome of the European hake

The growing accessibility to genomic resources using next-generation sequencing (NGS) technologies has revolutionized the application of molecular genetic tools to ecology and evolutionary studies in non-model organisms. Here we present the case study of the European hake (Merluccius merluccius), one of the most important demersal resources of European fisheries. Two sequencing platforms, the Roche 454 FLX (454) and the Illumina Genome Analyzer (GAII), were used for Single Nucleotide Polymorphisms (SNPs) discovery in the hake muscle transcriptome. De novo transcriptome assembly into unique contigs, annotation, and in silico SNP detection were carried out in parallel for 454 and GAII sequence data. High-throughput genotyping using the Illumina GoldenGate assay was performed for validating 1,536 putative SNPs. Validation results were analysed to compare the performances of 454 and GAII methods and to evaluate the role of several variables (e.g. sequencing depth, intron-exon structure, sequence quality and annotation). Despite well-known differences in sequence length and throughput, the two approaches showed similar assay conversion rates (approximately 43%) and percentages of polymorphic loci (67.5% and 63.3% for GAII and 454, respectively). Both NGS platforms therefore demonstrated to be suitable for large scale identification of SNPs in transcribed regions of non-model species, although the lack of a reference genome profoundly affects the genotyping success rate. The overall efficiency, however, can be improved using strict quality and filtering criteria for SNP selection (sequence quality, intron-exon structure, target region score).

Efficient Genotyping of KRAS Mutant Non-Small Cell Lung Cancer Using a Multiplexed Droplet Digital PCR Approach.

Droplet digital PCR (ddPCR) can be used to detect low frequency mutations in oncogene-driven lung cancer. The range of KRAS point mutations observed in NSCLC necessitates a multiplex approach to efficient mutation detection in circulating DNA. Here we report the design and optimisation of three discriminatory ddPCR multiplex assays investigating nine different KRAS mutations using PrimePCR™ ddPCR™ Mutation Assays and the Bio-Rad QX100 system. Together these mutations account for 95% of the nucleotide changes found in KRAS in human cancer. Multiplex reactions were optimised on genomic DNA extracted from KRAS mutant cell lines and tested on DNA extracted from fixed tumour tissue from a cohort of lung cancer patients without prior knowledge of the specific KRAS genotype. The multiplex ddPCR assays had a limit of detection of better than 1 mutant KRAS molecule in 2,000 wild-type KRAS molecules, which compared favourably with a limit of detection of 1 in 50 for next generation sequencing and 1 in 10 for Sanger sequencing. Multiplex ddPCR assays thus provide a highly efficient methodology to identify KRAS mutations in lung adenocarcinoma.

Establishing an EGFR mutation screening service for non-small cell lung cancer - sample quality criteria and candidate histological predictors.

INTRODUCTION: EGFR screening requires good quality tissue, sensitivity and turn-around time (TAT). We report our experience of routine screening, describing sample type, TAT, specimen quality (cellularity and DNA yield), histopathological description, mutation result and clinical outcome. METHODS: Non-small cell lung cancer (NSCLC) sections were screened for EGFR mutations (M+) in exons 18-21. Clinical, pathological and screening outcome data were collected for year 1 of testing. Screening outcome alone was collected for year 2. RESULTS: In year 1, 152 samples were tested, most (72%) were diagnostic. TAT was 4.9 days (95%confidence interval (CI)=4.5-5.5). EGFR-M+ prevalence was 11% and higher (20%) among never-smoking women with adenocarcinomas (ADCs), but 30% of mutations occurred in current/ex-smoking men. EGFR-M+ tumours were non-mucinous ADCs and 100% thyroid transcription factor (TTF1+). No mutations were detected in poorly differentiated NSCLC-not otherwise specified (NOS). There was a trend for improved overall survival (OS) among EGFR-M+ versus EGFR-M- patients (median OS=78 versus 17 months). In year 1, test failure rate was 19%, and associated with scant cellularity and low DNA concentrations. However 75% of samples with poor cellularity but representative of tumour were informative and mutation prevalence was 9%. In year 2, 755 samples were tested; mutation prevalence was 13% and test failure only 5.4%. Although samples with low DNA concentration (2.2 ng/μL), the mutation rate was 9.2%. CONCLUSION: Routine epidermal growth factor receptor (EGFR) screening using diagnostic samples is fast and feasible even on samples with poor cellularity and DNA content. Mutations tend to occur in better-differentiated non-mucinous TTF1+ ADCs. Whether these histological criteria may be useful to select patients for EGFR testing merits further investigation.

Frequency of cytokine gene promoter polymorphisms in the Northern Ireland Cystic Fibrosis population

It has been postulated that cytokine allele frequencies are gender and perhaps geographically-specific. Cytokine release is crucial in the regulation of the type and magnitude of the immune response. This study observed no differences in the frequency of cytokine promoter polymorphisms associated with variant levels of expression in patients with CIF and a non-CF population of Northern Ireland. (c) 2007 European Cystic Fibrosis Society. Published by Elsevier B.V. All rights reserved.

Serotonin Transporter Gene Polymorphism and Myocardial Infarction - Etude Cas-te'moins de L'Infarctus du Myocarde (ECTIM)

Background— Depression is a risk factor for myocardial infarction (MI). Selective serotonin reuptake inhibitors reduce this risk. The site of action is the serotonin transporter (SLC6A4), which is expressed in brain and blood cells. A functional polymorphism in the promoter region of the SLC6A4 gene has been described. This polymorphism may be associated with the risk of MI. Methods and Results— The SLC6A4 polymorphism has been investigated by polymerase chain reaction in 671 male patients with MI and in 688 controls from the Etude Cas-Témoins de l’Infarctus du Myocarde (ECTIM) multicentric study. Percentages for LL, LS, and SS genotypes were 35.5%, 45.4%, and 19.1%, respectively, for cases versus 28.1%, 49.1%, and 22.8%, respectively, for controls. S allele frequency was 41.8% and 47.4% for cases and controls, respectively. After adjustment for age and center by using multivariable logistic regression, the odds ratio for MI associated with the LL genotype was 1.40 (95% CI 1.11 to 1.76, P=0.0047). Conclusions— The LL genotype of the SLC6A4 polymorphism is associated with a higher risk of MI. This could be attributable to the effect of the polymorphism on serotonin-mediated platelet activation or smooth muscle cell proliferation or on other risk factors, such as depression or response to stress

A novel diagnostic test detects a low frequency of the hemicentin Gln5345Arg variant among Northern Irish age related macular degeneration patients.

Purpose: Age related macular degeneration (AMD) is a common cause of severe vision loss. Identification of genes involved in AMD will facilitate early detection and ultimately help to identify pathways for treatment for this disorder. The A16,263G mutation in the HEMICENTIN-1 gene produces a non-conservative substitution of arginine for glutamine at codon 5345 which has been implicated in familial AMD. The aim of this study is to develop a rapid diagnostic assay for the detection of this mutation and to evaluate its frequency in a sample of AMD patients. Methods: A primer probe set was designed from exon 104 of the HEMICENTIN-1 gene to differentiate between mutant and wild type alleles. A region spanning the mutation was amplified by PCR using a LightCycler (Roche Diagnostic). The mutation was then detected by melt curve analysis of the hybrid formed between the PCR product and a specific fluorescent probe. The frequency of the mutation within the Northern Ireland population was evaluated by assaying 508 affected AMD patients, 25 possibly affected and 163 controls. Results: This assay clearly discriminates between the A16,263G mutant and wild type HEMICENTIN-1 alleles. The wild type sequence has a single base mismatch with the probe which decreases the stability of the hybrid, resulting in a lower TM (TM=51.27 °C) than that observed for the perfectly matched mutant allele (TM=59.9 °C). The mutant allele was detected in only one of the 696 subjects, an affected AMD patient. Conclusions: We describe a rapid assay for the genotyping of the Gln5345Arg mutation using real-time fluorescence PCR to facilitate rapid processing of samples through combined amplification and detection steps. These characteristics are suitable for a clinical setting where high throughput diagnostic procedures are required. The frequency of this mutation within the Northern Ireland population has been estimated at 0.2%, concurring with previous findings that this mutation is a rare variant associated with AMD. A rapid diagnostic assay will facilitate a reliable and convenient evaluation of the frequency of the Gln5345Arg mutation and its association with AMD within other populations.

Impact of antibiotics on conjugational resistance gene transfer in Staphylococcus aureus in sewage.

Ohlsen K, Ternes T, Werner G, Wallner U, Löffler D, Ziebuhr W, Witte W, Hacker J. Institute for Molecular Biology of Infectious Diseases, The University of Würzburg, Röntgenring 11, D-97070 Würzburg, Germany. knut.ohlsen@mail.uni-wuerzburg.de The growing rate of microbial pathogens becoming resistant to standard antibiotics is an important threat to public health. In order to assess the role of antibiotics in the environment on the spread of resistance factors, the impact of subinhibitory concentrations of antibiotics in sewage on gene transfer was investigated using conjugative gentamicin resistance (aacA-aphD) plasmids of Staphylococcus aureus. Furthermore, the concentration of antibiotics in hospital sewage was measured by high-performance liquid chromatography (HPLC)-electrospray tandem mass spectrometry. Several antibiotics were found to be present in sewage, e.g. ciprofloxacin up to 0.051 mgl(-1) and erythromycin up to 0.027 mgl(-1). Resistance plasmid transfer occurred both on solidified (dewatered) sewage and in liquid sewage in a bioreactor with a frequency of 1.1x10(-5)-5.0x10(-8). However, low-level concentrations of antibiotics measured in sewage are below concentrations that can increase plasmid transfer frequencies of gentamicin resistance plasmids of staphylococci.

Cyclooxygenase-2 and inducible nitric oxide synthase gene polymorphisms and risk of reflux esophagitis, Barrett's Esophagus, and esophageal adenocarcinoma

The incidence of esophageal adenocarcinoma has increased in recent years, and Barrett's esophagus is a recognized risk factor. Gastroesophageal reflux of acid and/or bile is linked to these conditions and to reflux esophagitis. Inflammatory disorders can lead to carcinogenesis through activation of "prosurvival genes," including cyclooxygenase-2 (COX-2) and inducible nitric oxide synthase (iNOS). Increased expression of these enzymes has been found in esophageal adenocarcinoma, Barrett's esophagus, and reflux esophagitis. Polymorphic variants in COX-2 and iNOS genes may be modifiers of risk of these conditions. In a population-based case-control study, we examined associations of the COX-2 8473 T>C and iNOS Ser 608 Leu (C>T) polymorphisms with risk of esophageal adenocarcinoma, Barrett's esophagus, and reflux esophagitis. Genomic DNA was extracted from blood samples collected from cases of esophageal adenocarcinoma (n = 210), Barrett's esophagus (n = 212), and reflux esophagitis (n = 230) and normal population controls frequency matched for age and sex (n = 248). Polymorphisms were genotyped using TaqMan allelic discrimination assays. Odds ratios and 95% confidence intervals were obtained from logistic regression models adjusted for potential confounding factors. The presence of at least one COX-2 8473 C allele was associated with a significantly increased risk of esophageal adenocarcinoma (adjusted odds ratio, 1.58; 95% confidence interval, 1.04-2.40). There was no significant association between this polymorphism and risk of Barrett's esophagus or reflux esophagitis or between the iNOS Ser 608 Leu polymorphism and risk of these esophageal conditions. Our study suggests that the COX-2 8473 C allele is a potential genetic marker for susceptibility to esophageal adenocarcinoma.

The role of human demographic history in determining the distribution and frequency of transferase-deficient galactosaemia mutations

Classical or transferase-deficient galactosaemia is an inherited metabolic disorder caused by mutation in the human Galactose-1-phosphate uridyl transferase (GALT) gene. Of some 170 causative mutations reported, fewer than 10% are observed in more than one geographic region or ethnic group. To better understand the population history of the common GALT mutations, we have established a haplotyping system for the GALT locus incorporating eight single nucleotide polymorphisms and three short tandem repeat markers. We analysed haplotypes associated with the three most frequent GALT gene mutations, Q188R, K285N and Duarte-2 (D2), and estimated their age. Haplotype diversity, in conjunction with measures of genetic diversity and of linkage disequilibrium, indicated that Q188R and K285N are European mutations. The Q188R mutation arose in central Europe within the last 20 000 years, with its observed east-west cline of increasing relative allele frequency possibly being due to population expansion during the re-colonization of Europe by Homo sapiens in the Mesolithic age. K285N was found to be a younger mutation that originated in Eastern Europe and is probably more geographically restricted as it arose after all major European population expansions. The D2 variant was found to be an ancient mutation that originated before the expansion of Homo sapiens out of Africa. Heredity (2010) 104, 148-154; doi:10.1038/hdy.2009.84; published online 29 July 2009

Identification of gene networks associated with acute myeloid leukemia by comparative molecular methylation and expression profiling

Around 80% of acute myeloid leukemia (AML) patients achieve a complete remission, however many will relapse and ultimately die of their disease. The association between karyotype and prognosis has been studied extensively and identified patient cohorts as having favourable [e.g. t(8; 21), inv (16)/t(16; 16), t(15; 17)], intermediate [e.g. cytogenetically normal (NK-AML)] or adverse risk [e.g. complex karyotypes]. Previous studies have shown that gene expression profiling signatures can classify the sub-types of AML, although few reports have shown a similar feature by using methylation markers. The global methylation patterns in 19 diagnostic AML samples were investigated using the Methylated CpG Island Amplification Microarray (MCAM) method and CpG island microarrays containing 12,000 CpG sites. The first analysis, comparing favourable and intermediate cytogenetic risk groups, revealed significantly differentially methylated CpG sites (594 CpG islands) between the two subgroups. Mutations in the NPM1 gene occur at a high frequency (40%) within the NK-AML subgroup and are associated with a more favourable prognosis in these patients. A second analysis comparing the NPM1 mutant and wild-type research study subjects again identified distinct methylation profiles between these two subgroups. Network and pathway analysis revealed possible molecular mechanisms associated with the different risk and/or mutation sub-groups. This may result in a better classification of the risk groups, improved monitoring targets, or the identification of novel molecular therapies.

Thrombopoietin, flt3-ligand and c-kit-ligand modulate HOX gene expression in expanding cord blood CD133(+) cells

Haemopoietic stem/progenitor cell (HSPC) development is regulated by extrinsic and intrinsic stimuli. Extrinsic modulators include growth factors and cell adhesion molecules, whereas intrinsic regulation is achieved with many transcription factor families, of which the HOX gene products are known to be important in haemopoiesis. Umbilical cord blood CD133(+) HSPC proliferation potential was tested in liquid culture with 'TPOFLK' (thrombopoietin, flt-3 ligand and c-kit ligand, promoting HSPC survival and self-renewal), in comparison to 'K36EG' (c-kit-ligand, interleukins-3 and -6, erythropoietin and granulocyte colony-stimulating factor, inducing haemopoietic differentiation). TPOFLK induced a higher CD133(+) HSPC proliferation (up to 60-fold more, at week 8) and maintained a higher frequency of the primitive colony-forming cells than K36EG. Quantitative polymerase chain reaction analysis revealed opposite expression patterns for specific HOX genes in expanding cord blood CD133(+) HSPC. After 8 weeks in liquid culture, TPOFLK increased the expression of HOX B3, B4 and A9 (associated with uncommitted HSPC) and reduced the expression of HOX B8 and A10 (expressed in committed myeloid cells) when compared to K36EG. These results suggest that TPOFLK induces CD133(+) HSPC proliferation, self-renewal and maintenance, up-regulation of HOX B3, B4 and A9 and down-regulation of HOX B8 and A10 gene expression.

A GREM1 gene variant associates with diabetic nephropathy

Gremlin, a cell growth and differentiation factor, promotes the development of diabetic nephropathy in animal models, but whether GREM1 gene variants associate with diabetic nephropathy is unknown. We comprehensively screened the 5' upstream region (including the predicted promoter), all exons, intron-exon boundaries, complete untranslated regions, and the 3' region downstream of the GREM1 gene. We identified 31 unique variants, including 24 with a minor allele frequency exceeding 5%, and 9 haplotype-tagging single nucleotide polymorphisms (htSNPs). We selected one additional variant that we predicted to alter transcription factor binding. We genotyped 709 individuals with type 1 diabetes of whom 267 had nephropathy (cases) and 442 had no evidence of kidney disease (controls). Three individual SNPs significantly associated with nephropathy at the 5% level, and two remained significant after adjustment for multiple testing. Subsequently, we genotyped a replicate population comprising 597 cases and 502 controls: this population supported an association with one of the SNPs (rs1129456; P = 0.0003). Combined analysis, adjusted for recruitment center (n = 8), suggested that the T allele conferred greater odds of nephropathy (OR 1.69; 95% CI 1.36 to 2.11). In summary, the GREM1 variant rs1129456 associates with diabetic nephropathy, perhaps explaining some of the genetic susceptibility to this condition.