Accurate detection of subclonal single nucleotide variants in whole genome amplified and pooled cancer samples using HaloPlex target enrichment

Background

Target enrichment and resequencing is a widely used approach for identification of cancer genes and genetic variants associated with diseases. Although cost effective compared to whole genome sequencing, analysis of many samples constitutes a significant cost, which could be reduced by pooling samples before capture. Another limitation to the number of cancer samples that can be analyzed is often the amount of available tumor DNA. We evaluated the performance of whole genome amplified DNA and the power to detect subclonal somatic single nucleotide variants in non-indexed pools of cancer samples using the HaloPlex technology for target enrichment and next generation sequencing.

Results

We captured a set of 1528 putative somatic single nucleotide variants and germline SNPs, which were identified by whole genome sequencing, with the HaloPlex technology and sequenced to a depth of 792–1752. We found that the allele fractions of the analyzed variants are well preserved during whole genome amplification and that capture specificity or variant calling is not affected. We detected a large majority of the known single nucleotide variants present uniquely in one sample with allele fractions as low as 0.1 in non-indexed pools of up to ten samples. We also identified and experimentally validated six novel variants in the samples included in the pools.

Conclusion

Our work demonstrates that whole genome amplified DNA can be used for target enrichment equally well as genomic DNA and that accurate variant detection is possible in non-indexed pools of cancer samples. These findings show that analysis of a large number of samples is feasible at low cost, even when only small amounts of DNA is available, and thereby significantly increases the chances of indentifying recurrent mutations in cancer samples.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-14-856) contains supplementary material, which is available to authorized users.     

Deep mtDNA subdivision within Linnean species in an endemic radiation of tiger beetles from New Zealand (genus Neocicindela)

The invertebrate fauna of New Zealand is of great interest as a geologically tractable model for the study of species diversification, but direct comparisons with closely related lineages elsewhere are lacking. Integrating population-level analyses with studies of taxonomy and clade diversification, we performed mtDNA analysis on Neocicindela (Cicindelidae, tiger beetles) for a broad sample of populations from 11 of 12 known species and 161 specimens (three loci, 1883 nucleotides), revealing 123 distinct haplotypes. Phylogenetic reconstruction recovered two main lineages, each composed of 5-6 Linnean species whose origin was dated to 6.66 and 7.26 Mya, while the Neocicindela stem group was placed at 10.82 ± 0.48 Mya. Species delimitation implementing a character-based (diagnostic) species concept recognized 19 species-level groups that were in general agreement with Linnean species but split some of these into mostly allopatric subgroups. Tree-based methods of species delimitation using a mixed Yule-coalescence model were inconclusive, and recognized 32-51 entities (including singletons), splitting existing species into up to 8 partially sympatric groups. These findings were different from patterns in the Australian sister genus Rivacindela, where character-based and tree-based methods were previously shown to produce highly congruent groupings. In Neocicindela, the pattern of mtDNA variation was characterized by high intra-population and intra-species haplotype divergence, the coexistence of divergent haplotypes in sympatry, and a poor correlation of genetic and geographic distance. These observations combined suggest a scenario of phylogeographic divergence and secondary contact driven by orogenetic and climatic changes of the Pleistocene/Pliocene. The complex evolutionary history of most species of Neocicindela due to the relative instability of the New Zealand biota resulted in populations of mixed ancestry but not in a general loss of genetic variation.     

Tissue-specific Salmonella Typhimurium gene expression during persistence in pigs

Salmonellosis caused by Salmonella Typhimurium is one of the most important bacterial zoonotic diseases. The bacterium persists in pigs resulting in asymptomatic 'carrier pigs', generating a major source for Salmonella contamination of pork. Until now, very little is known concerning the mechanisms used by Salmonella Typhimurium during persistence in pigs. Using in vivo expression technology (IVET), a promoter-trap method based on ΔpurA attenuation of the parent strain, we identified 37 Salmonella Typhimurium genes that were expressed 3 weeks post oral inoculation in the tonsils, ileum and ileocaecal lymph nodes of pigs. Several genes were expressed in all three analyzed organs, while other genes were only expressed in one or two organs. Subsequently, the identified IVET transformants were pooled and reintroduced in pigs to detect tissue-specific gene expression patterns. We found that efp and rpoZ were specifically expressed in the ileocaecal lymph nodes during Salmonella peristence in pigs. Furthermore, we compared the persistence ability of substitution mutants for the IVET-identified genes sifB and STM4067 to that of the wild type in a mixed infection model. The ΔSTM4067::kanR was significantly attenuated in the ileum contents, caecum and caecum contents and faeces of pigs 3 weeks post inoculation, while deletion of the SPI-2 effector gene sifB did not affect Salmonella Typhimurium persistence. Although our list of identified genes is not exhaustive, we found that efp and rpoZ were specifically expressed in the ileocaecal lymph nodes of pigs and we identified STM4067 as a factor involved in Salmonella persistence in pigs. To our knowledge, our study is the first to identify Salmonella Typhimurium genes expressed during persistence in pigs.     

The RACK1 signaling scaffold protein selectively interacts with Yersinia pseudotuberculosis virulence function

Many gram-negative bacteria use type III secretion systems to translocate effector proteins into host cells. These effectors interfere with cellular functions in a highly regulated manner resulting in effects that are beneficial for the bacteria. The pathogen Yersinia can resist phagocytosis by eukaryotic cells by translocating Yop effectors into the target cell cytoplasm. This is called antiphagocytosis, and constitutes an important virulence feature of this pathogen since it allows survival in immune cell rich lymphoid organs. We show here that the virulence protein YopK has a role in orchestrating effector translocation necessary for productive antiphagocytosis. We present data showing that YopK influences Yop effector translocation by modulating the ratio of the pore-forming proteins YopB and YopD in the target cell membrane. Further, we show that YopK that can interact with the translocators, is exposed inside target cells and binds to the eukaryotic signaling protein RACK1. This protein is engaged upon Y. pseudotuberculosis-mediated β1-integrin activation and localizes to phagocytic cups. Cells with downregulated RACK1 levels are protected from antiphagocytosis. This resistance is not due to altered levels of translocated antiphagocytic effectors, and cells with reduced levels of RACK1 are still sensitive to the later occurring cytotoxic effect caused by the Yop effectors. Further, a yopK mutant unable to bind RACK1 shows an avirulent phenotype during mouse infection, suggesting that RACK1 targeting by YopK is a requirement for virulence. Together, our data imply that the local event of Yersinia-mediated antiphagocytosis involves a step where YopK, by binding RACK1, ensures an immediate specific spatial delivery of antiphagocytic effectors leading to productive inhibition of phagocytosis.     

Tryptophan orientations in membrane-bound gramicidin and melittin-a comparative linear dichroism study on transmembrane and surface-bound peptides

In the search for methods to study structure and function of membrane-associated proteins and peptides flow linear dichroism, LD, spectroscopy has emerged as a promising technique. Using shear-aligned lipid vesicles, conformations and binding geometries of membrane-bound bio-macromolecules can be assessed. Here we investigate anchoring properties and specific orientations of tryptophan relative to the peptide backbone and to the membrane normal for the model peptides gramicidin and melittin. We have monitored the conformational change associated with the refolding of non-channel gramicidin into its channel form, and quantitatively determined the average orientations of its tryptophan transition moments, suggesting that these residues adopt a well-defined orientation at the membrane interface. An important conclusion regards the structural variation of gramicidin between these two distinct transmembrane forms. Whilst circular dichroism (CD) spectra, as has been reported before, vary strongly between the two forms suggesting their structures might be quite different, the LD results clearly evidence both the peptide backbone orientation and tryptophan side-chain positioning to be very similar. The latter are oriented in accord with what is expected from their role to anchor peptide termini to the membrane surface. The variations in CD could be due to, the in LD observed, minor shifts in mutual orientation and distance between neighbouring tryptophans sensitively determining their exciton interactions. Our data dispute that the non-channel form of membrane-bound gramicidin would be any of the intertwined forms often observed in crystal as the positioning of tryptophans along the peptide axis would not be compatible with the strong interfacial positioning observed here. The general role of tryptophans as interfacial anchors is further assessed for melittin whose conformation shows considerable angular spread, consistent with a carpet model of its mechanism for induced membrane leakage, and a predominantly surface-aligned membrane orientation governed by amphipathic interactions.     

Complete (1)H and (13)C NMR chemical shift assignments of mono-, di-, and trisaccharides as basis for NMR chemical shift predictions of polysaccharides using the computer program casper

The computer program casper uses (1)H and (13)C NMR chemical shift data of mono- to trisaccharides for the prediction of chemical shifts of oligo- and polysaccharides. In order to improve the quality of these predictions the (1)H and (13)C, as well as (31)P when applicable, NMR chemical shifts of 30 mono-, di-, and trisaccharides were assigned. The reducing sugars gave two distinct sets of NMR resonances due to the α- and β-anomeric forms. In total 35 (1)H and (13)C NMR chemical shift data sets were obtained from the oligosaccharides. One- and two-dimensional NMR experiments were used for the chemical shift assignments and special techniques were employed in some cases such as 2D (1)H,(13)C-HSQC Hadamard Transform methodology which was acquired approximately 45 times faster than a regular t(1) incremented (1)H,(13)C-HSQC experiment and a 1D (1)H,(1)H-CSSF-TOCSY experiment which was able to distinguish spin-systems in which the target protons were only 3.3Hz apart. The (1)H NMR chemical shifts were subsequently refined using total line-shape analysis with the PERCH NMR software. The acquired NMR data were then utilized in the casper program (http://www.casper.organ.su.se/casper/) for NMR chemical shift predictions of the O-antigen polysaccharides from Klebsiella O5, Shigella flexneri serotype X, and Salmonella arizonae O62. The data were compared to experimental data of the polysaccharides from the two former strains and the lipopolysaccharide of the latter strain showing excellent agreement between predicted and experimental (1)H and (13)C NMR chemical shifts.     

Gene Expression Signature of Normal Cell-of-Origin Predicts Ovarian Tumor Outcomes

The potential role of the cell-of-origin in determining the tumor phenotype has been raised, but not adequately examined. We hypothesized that distinct cells-of-origin may play a role in determining ovarian tumor phenotype and outcome. Here we describe a new cell culture medium for in vitro culture of paired normal human ovarian (OV) and fallopian tube (FT) epithelial cells from donors without cancer. While these cells have been cultured individually for short periods of time, to our knowledge this is the first long-term culture of both cell types from the same donors. Through analysis of the gene expression profiles of the cultured OV/FT cells we identified a normal cell-of-origin gene signature that classified primary ovarian cancers into OV-like and FT-like subgroups; this classification correlated with significant differences in clinical outcomes. The identification of a prognostically significant gene expression signature derived solely from normal untransformed cells is consistent with the hypothesis that the normal cell-of-origin may be a source of ovarian tumor heterogeneity and the associated differences in tumor outcome.     

Human Herpesvirus 6A Partially Suppresses Functional Properties of DC without Viral Replication

Human herpesvirus 6A (HHV-6A) is a common virus with a worldwide distribution that has been associated with multiple sclerosis. Whether HHV-6A can replicate in dendritic cells (DC) and how the infection might modulate the functional properties of the cell are currently not well known and need further investigations. Here, we show that a non-productive infection of HHV-6A in DC leads to the up-regulation of HLA-ABC, via autocrine IFN-α signaling, as well as the up-regulation of HLA-DR and CD86. However, HHV-6A exposure reduces IL-8 secretion by DC and their capacity to stimulate allogenic T cell proliferation. The ability to suppress DC functions important for activation of innate and adaptive immune responses might be one successful strategy by which HHV-6A avoids the induction of appropriate host defense mechanisms, and thus facilitating persistent infection.     

Blood Flow Changes Coincide with Cellular Rearrangements during Blood Vessel Pruning in Zebrafish Embryos

After the initial formation of a highly branched vascular plexus, blood vessel pruning generates a hierarchically structured network with improved flow characteristics. We report here on the cellular events that occur during the pruning of a defined blood vessel in the eye of developing zebrafish embryos. Time-lapse imaging reveals that the connection of a new blood vessel sprout with a previously perfused multicellular endothelial tube leads to the formation of a branched, Y-shaped structure. Subsequently, endothelial cells in parts of the previously perfused branch rearrange from a multicellular into a unicellular tube, followed by blood vessel detachment. This process is accompanied by endothelial cell death. Finally, we show that differences in blood flow between neighboring vessels are important for the completion of the pruning process. Our data suggest that flow induced changes in tubular architecture ensure proper blood vessel pruning.