Quality control parameters on a large dataset of regionally dissected human control brains for whole genome expression studies

We are building an open-access database of regional human brain expression designed to allow the genome-wide assessment of genetic variability on expression. Array and RNA sequencing technologies make assessment of genome-wide expression possible. Human brain tissue is a challenging source for this work because it can only be obtained several and variable hours post-mortem and after varying agonal states. These variables alter RNA integrity in a complex manner. In this report, we assess the effect of post-mortem delay, agonal state and age on gene expression, and the utility of pH and RNA integrity number as predictors of gene expression as measured on 1266 Affymetrix Exon Arrays. We assessed the accuracy of the array data using QuantiGene, as an independent non-PCR-based method. These quality control parameters will allow database users to assess data accuracy. We report that within the parameters of this study post-mortem delay, agonal state and age have little impact on array quality, array data are robust to variable RNA integrity, and brain pH has only a small effect on array performance. QuantiGene gave very similar expression profiles as array data. This study is the first step in our initiative to make human, regional brain expression freely available.     

Screening the human exome: a comparison of whole genome and whole transcriptome sequencing

Background  

There is considerable interest in the development of methods to efficiently identify all coding variants present in large sample sets of humans. There are three approaches possible: whole-genome sequencing, whole-exome sequencing using exon capture methods, and RNA-Seq. While whole-genome sequencing is the most complete, it remains sufficiently expensive that cost effective alternatives are important.     

Exploring whole genome amplification as a DNA recovery tool for molecular genetic studies.

The whole genome amplification (WGA) protocol evaluated during this study, GenomiPhi DNA amplification kit, is a novel method that is not based on polymerase chain reaction but rather relies on the highly processive and high fidelity Phi29 DNA polymerase to replicate linear genomic DNA by multiple strand displacement amplification. As little as 1 ng of genomic DNA template is sufficient to produce microgram quantities of high molecular weight DNA. The question explored during this study is whether such a WGA method is appropriate to reliably replenish and even recover depleted DNA samples that can be used for downstream genetic analysis. A series of human DNA samples was tested in our laboratory and validated using such analytical methods as gene-specific polymerase chain reaction, direct sequencing, microsatellite marker analysis, and single nucleotide polymorphism allelic discrimination using TaqMan and Pyrosequencing chemistries. Although degraded genomic DNA is not a good template for Phi29 WGA, this method is a powerful tool to replenish depleted DNA stocks and to increase the amount of sample for which biological tissue availability is scarce. The testing performed during the validation phase of the study indicates no discernable difference between WGA samples and the original DNA templates. Thus, GenomiPhi WGA can be used to increase precious or depleted DNA stocks, thereby extending the life of a family-based linkage analysis project and increasing statistical power.     

Viral clearance studies on new and used chromatography resins: critical review of a large dataset.

Viral clearance studies for na?ve and maximally cycled chromatographic resins used for cGMP recombinant protein production are reviewed for three products, comprising 10 different chromatographic steps, including affinity, ion exchange, immobilized metal ion affinity, and hydrophobic interaction modes. Thirty-two separate studies were conducted (over 90 runs in total). No consistent reductions in model virus clearance were observed with used resins. The results address the reproducibility of virus clearance studies conducted by different scientists over several years at multiple contract labs. The log reduction values (LRVs) are typically within 0.5 LRVs for new and used resin, but varied as much as 2 LRVs for resins showing no functional deterioration. This relatively large difference is not believed to reflect resin changes, but highlights the challenges encountered in modeling column clearance. Production column performance and cleaning efficacy are demonstrated for these steps by trending mock runs, impurity removal and product recovery. No deterioration in cGMP column performance is seen over the established resin lifetimes, confirming that the resin regeneration and sanitization procedures restore the resins to a suitable initial state without damage. It is proposed that for some chromatography steps, the combination of lab-scale cycling studies confirming consistent performance throughout the resin lifetime and monitoring of cGMP manufacturing preclude the need for virus clearance studies on maximally cycled resin.     

Simulating association studies: a data-based resampling method for candidate regions or whole genome scans

MOTIVATION: Reductions in genotyping costs have heightened interest in performing whole genome association scans and in the fine mapping of candidate regions. Improvements in study design and analytic techniques will require the simulation of datasets with realistic patterns of linkage disequilibrium and allele frequencies for typed SNPs. METHODS: We describe a general approach to simulate genotyped datasets for standard case-control or affected child trio data, by resampling from existing phased datasets. The approach allows for considerable flexibility in disease models, potentially involving a large number of interacting loci. The method is most applicable for diseases caused by common variants that have not been under strong selection, a class specifically targeted by the International HapMap project. RESULTS: Using the three population Phase I/II HapMap data as a testbed for our approach, we have implemented the approach in HAP-SAMPLE, a web-based simulation tool.     

An ileal Crohn's disease gene signature based on whole human genome expression profiles of disease unaffected ileal mucosal biopsies

Previous genome-wide expression studies have highlighted distinct gene expression patterns in inflammatory bowel disease (IBD) compared to control samples, but the interpretation of these studies has been limited by sample heterogeneity with respect to disease phenotype, disease activity, and anatomic sites. To further improve molecular classification of inflammatory bowel disease phenotypes we focused on a single anatomic site, the disease unaffected proximal ileal margin of resected ileum, and three phenotypes that were unlikely to overlap: ileal Crohn's disease (ileal CD), ulcerative colitis (UC), and control patients without IBD. Whole human genome (Agilent) expression profiling was conducted on two independent sets of disease-unaffected ileal samples collected from the proximal margin of resected ileum. Set 1 (47 ileal CD, 27 UC, and 25 Control non-IBD patients) was used as the training set and Set 2 was subsequently collected as an independent test set (10 ileal CD, 10 UC, and 10 control non-IBD patients). We compared the 17 gene signatures selected by four different feature-selection methods to distinguish ileal CD phenotype with non-CD phenotype. The four methods yielded different but overlapping solutions that were highly discriminating. All four of these methods selected FOLH1 as a common feature. This gene is an established biomarker for prostate cancer, but has not previously been associated with Crohn's disease. Immunohistochemical staining confirmed increased expression of FOLH1 in the ileal epithelium. These results provide evidence for convergent molecular abnormalities in the macroscopically disease unaffected proximal margin of resected ileum from ileal CD subjects.     

A comparative whole genome analysis of Helicobacter pylori from a human dense South Asian setting

Helicobacter pylori, a Gram-negative bacterium, is associated with a wide range of gastric diseases such as gastritis, duodenal ulcer, and gastric cancer. The prevalence of H pylori and risk of disease vary in different parts of the world based on the prevailing bacterial lineage. Here, we present a contextual and comparative genomics analysis of 20 clinical isolates of H pylori from patients in Bangladesh. Despite a uniform host ethnicity (Bengali), isolates were classified as being part of the HpAsia2 (50%) or HpEurope (50%) population. Out of twenty isolates, eighteen isolates were cagA positive, with two HpEurope isolates being cagA negative, three EPIYA motif patterns (AB, AB-C, and ABC-C) were observed among the cagA-positive isolates. Three vacA genotypes were observed with the s1m1i1dic1 genotype observed in 75% of isolates; the s1m2i1d1c2 and s2m2i2d2c2 genotypes were found to be 15% and 10% of isolates, respectively. The non-virulent genotypes s2m2i2d2c2 was only observed in HpEurope population isolates. Genotypic analysis of oipA gene, present in all isolates, revealed five different patterns of the CT repeat; all HpAsia2 isolates were in “ON” while 20% of HpEurope isolates were genotypically “OFF.” The three blood group antigen binding adhesins encoded genes (bab genes) examined and we observed that the most common genotype was (babA/babB/-) found in eight isolates, notably six were HpAsia2 isolates. The babA gene was found in all HpAsia2 isolates but present in only half of the HpEurope isolates. In silico antibiotic susceptibility analysis revealed that 40% of the strains were multi-drug resistant. Mutations associated with resistance to metronidazole, fluoroquinolone, and clarithromycin were detected 90%, 45%, and 5%, respectively, in H pylori strain. In conclusion, it is evident that two populations of H pylori with similar antibiotic profiles are predominant in Bangladesh, and it appears that genotypically the HpAisa2 isolates are potentially more virulent than the HpEurope isolates.     

Optimized PCR fragments for heteroduplex analysis of the whole human mitochondrial genome with denaturing HPLC

Denaturing high pressure liquid chromatography (dHPLC) is an efficient method for discovery of unknown mutations by heteroduplex analysis of PCR fragments. For comprehensive mutation scanning of the whole 16.569 bp human mitochondrial genome, we developed a set of 67 primer pairs defining overlapping PCR fragments that are well suited for heteroduplex analysis. The aim of our optimization efforts was to ensure that point mutations are detectable at every nucleotide position of each amplicon. Some GC-rich regions of mitochondrial DNA (mtDNA) were found to have unfavourable melting profiles in all possible amplicons, therefore requiring GC-clamps at the end of one or both oligonucleotide PCR primers. Following detection of a heteroduplex pattern by dHPLC, our primers can also be employed for DNA sequencing to identify the underlying mutation. In case of heteroplasmic mutations with a low proportion of mutant mtDNA, a fragment collector is useful to recover the heteroduplex peak, which contains mutant and wildtype DNA molecules in a 1:1 ratio.     

Quality control and integration of genotypes from two calling pipelines for whole genome sequence data in the Alzheimer's disease sequencing project

The Alzheimer's Disease Sequencing Project (ADSP) performed whole genome sequencing (WGS) of 584 subjects from 111 multiplex families at three sequencing centers. Genotype calling of single nucleotide variants (SNVs) and insertion-deletion variants (indels) was performed centrally using GATK-HaplotypeCaller and Atlas V2. The ADSP Quality Control (QC) Working Group applied QC protocols to project-level variant call format files (VCFs) from each pipeline, and developed and implemented a novel protocol, termed “consensus calling,” to combine genotype calls from both pipelines into a single high-quality set. QC was applied to autosomal bi-allelic SNVs and indels, and included pipeline-recommended QC filters, variant-level QC, and sample-level QC. Low-quality variants or genotypes were excluded, and sample outliers were noted. Quality was assessed by examining Mendelian inconsistencies (MIs) among 67 parent-offspring pairs, and MIs were used to establish additional genotype-specific filters for GATK calls. After QC, 578 subjects remained. Pipeline-specific QC excluded ~12.0% of GATK and 14.5% of Atlas SNVs. Between pipelines, ~91% of SNV genotypes across all QCed variants were concordant; 4.23% and 4.56% of genotypes were exclusive to Atlas or GATK, respectively; the remaining ~0.01% of discordant genotypes were excluded. For indels, variant-level QC excluded ~36.8% of GATK and 35.3% of Atlas indels. Between pipelines, ~55.6% of indel genotypes were concordant; while 10.3% and 28.3% were exclusive to Atlas or GATK, respectively; and ~0.29% of discordant genotypes were. The final WGS consensus dataset contains 27,896,774 SNVs and 3,133,926 indels and is publicly available.     

A system for assessment of monokine gene expression using human whole blood

Monocyte derived cytokines (monokines) are important mediators in inflammatory diseases and cancer. Control of monokine expression is also a major therapeutic target in autoimmune inflammation. Whole blood cultures permit examination of monokine expression under conditions which emulate the in-vivo environment whilst avoiding many of the artefacts associated with monocyte separation and culture. Here we describe a system for measuring interleukin-1β, interleukin-1α, interleukin-6 and tumour necrosis factor-α mRNA in stimulated human whole blood ex-vivo, which can be applied to specimens from treated patients. Oligodeoxyribonucleotide probes are designed to allow standardisation of hybridisation and washing procedures. Washing and reprobing of membranes in appropriate sequence permits measurement of each monokine mRNA and mRNA for glyceraldehyde-3-phosphate dehydrogenase in only 7 ml of lipopolysaccharide-stimulated human blood. The method has been used successfully in studies of dexamethasone and methotrexate action on lipopolysaccharide stimulated IL-β gene expression.     

A whole genome sequence association study for puberty in a large Duroc × Erhualian F2 population

A large proportion of gilts and sows are culled from reproduction populations because of anestrus and pubertal reproductive failure. Selecting early onset of puberty gilts has a favorable effect on sows’ reproductivity. However, age at puberty is hard to be routinely measured in commercial herds. With molecular genetic predictors, identifying individuals that have a propensity for early onset of puberty can be simplified. We previously performed genome scanning and a genome‐wide association study for puberty in an F2 resource population using 183 microsatellites and 62 125 SNPs respectively. The detection power and resolution of identified quantitative trait loci were very low. Herein, we re‐sequenced 19 founders of the F2 resource population in high coverage, and whole genome sequences of F2 individuals were imputed to perform an association study for reproductive traits. A total of 2339 SNPs associated with pubertal reproductive failure were identified in the region of 30.94–40.74 Mb on SSC7, with the top one, positioned at 33.36 Mb, explaining 16% of the phenotypic variances. We improved the magnitude of the P‐value by 10E+5 to 10E+7 using the whole genome sequence rather than using low/middle density markers as in previous studies, and we narrowed down the QTL confidence interval to 5.25 Mb. Combining the annotation of gene function, RAB23 and BAK1 were perceived as the most compelling candidate genes. The identified loci may be useful in culling sows failing to show estrus by marker‐assisted selection to increase reproductive efficiency of swine herds.     

An inducible recA expression Bacillus subtilis genome vector for stable manipulation of large DNA fragments

Background

The Bacillus subtilis genome (BGM) vector is a novel cloning system based on the natural competence that enables B. subtilis to import extracellular DNA fragments into the cell and incorporate the recombinogenic DNA into the genome vector by homologous recombination. The BGM vector system has several attractive properties, such as a megabase cloning capacity, stable propagation of cloned DNA inserts, and various modification strategies using RecA-mediated homologous recombination. However, the endogenous RecA activity may cause undesirable recombination, as has been observed in yeast artificial chromosome systems. In this study, we developed a novel BGM vector system of an inducible recA expression BGM vector (iREX), in which the expression of recA can be controlled by xylose in the medium.

Results

We constructed the iREX system by introducing the xylose-inducible recA expression cassette followed by the targeted deletion of the endogenous recA. Western blot analysis showed that the expression of recA was strictly controlled by xylose in the medium. In the absence of xylose, recA was not expressed in the iREX, and the RecA-mediated recombination reactions were greatly suppressed. By contrast, the addition of xylose successfully induced RecA expression, which enabled the iREX to exploit the same capacities of transformation and gene modifications observed with the conventional BGM vector. In addition, an evaluation of the stability of the cloned DNA insert demonstrated that the DNA fragments containing homologous sequences were more stably maintained in the iREX by suppressing undesirable homologous recombination.

Conclusions

We developed a novel BGM vector with inducible recA expression system, iREX, which enables us to manipulate large DNA fragments more stably than the conventional BGM vector by suppressing undesirable recombination. In addition, we demonstrate that the iREX can be applied to handling the DNA, which has several homologous sequences, such as multiple-reporter expression cassettes. Thus, the iREX expands the utility of the BGM vector as a platform for engineering large DNA fragments.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1425-4) contains supplementary material, which is available to authorized users.     

Construction and validation of a Sinorhizobium meliloti whole genome DNA microarray: genome-wide profiling of osmoadaptive gene expression

Based on the complete Sinorhizobium meliloti genome sequence we established DNA microarrays as a comprehensive tool for systematic genome-wide gene expression analysis in S. meliloti 1021. For these PCR fragment-based microarrays, called Sm6kPCR, a collection of probes for the 6207 predicted protein-coding genes consisting of 6046 gene-specific PCR fragments and 161 70 mer oligonucleotides was arrayed in high density on glass slides. To obtain these PCR fragments primer pairs were designed to amplify internal gene-specific DNA fragments of 80-350 bp. Additionally, these primers were characterized by a 5' extension that allowed for reamplification using standard primers after the first amplification employing the specific primers. In order to ascertain the quality of the Sm6kPCR microarrays and to validate gene expression studies in S. meliloti parallel hybridizations based on RNA samples obtained from cells cultured under identical conditions were performed. In addition, gene expression in S. meliloti in response to an osmotic upshift imposed by the addition of 0.38 M NaCl was monitored. 137 genes were identified showing significant changes in gene expression resulting from the osmotic upshift. From these genes 52 were induced and 85 genes were repressed. Among the genes displaying different RNA levels some functional groups could be identified that are particularly remarkable. Repression was observed for 8 genes related to motility and chemotaxis, 7 genes encoding amino acid biosynthesis enzymes and 15 genes involved in iron uptake whereas 14 genes involved in transport of small molecules and 4 genes related to polysaccharide biosynthesis were induced.