STATISTICAL METHODS FOR THE ANALYSIS OF CANCER GENOME SEQUENCING DATA

The purpose of cancer genome sequencing studies is to determine the nature and types of alterations present in a typical cancer and to discover genes mutated at high frequencies. In this article we discuss statistical methods for the analysis of data generated in these studies. We place special emphasis on a two-stage study design introduced by Sjoblom et al.[1]. In this context, we describe statistical methods for constructing scores that can be used to prioritize candidate genes for further investigation and to assess the statistical signicance of the candidates thus identfied

Modulation of enhancer looping and differential gene targeting by Epstein-Barr virus transcription factors directs cellular reprogramming

Epstein-Barr virus (EBV) epigenetically reprogrammes B-lymphocytes to drive immortalization and facilitate viral persistence. Host-cell transcription is perturbed principally through the actions of EBV EBNA 2, 3A, 3B and 3C, with cellular genes deregulated by specific combinations of these EBNAs through unknown mechanisms. Comparing human genome binding by these viral transcription factors, we discovered that 25% of binding sites were shared by EBNA 2 and the EBNA 3s and were located predominantly in enhancers. Moreover, 80% of potential EBNA 3A, 3B or 3C target genes were also targeted by EBNA 2, implicating extensive interplay between EBNA 2 and 3 proteins in cellular reprogramming. Investigating shared enhancer sites neighbouring two new targets (WEE1 and CTBP2) we discovered that EBNA 3 proteins repress transcription by modulating enhancer-promoter loop formation to establish repressive chromatin hubs or prevent assembly of active hubs. Re-ChIP analysis revealed that EBNA 2 and 3 proteins do not bind simultaneously at shared sites but compete for binding thereby modulating enhancer-promoter interactions. At an EBNA 3-only intergenic enhancer site between ADAM28 and ADAMDEC1 EBNA 3C was also able to independently direct epigenetic repression of both genes through enhancer-promoter looping. Significantly, studying shared or unique EBNA 3 binding sites at WEE1, CTBP2, ITGAL (LFA-1 alpha chain), BCL2L11 (Bim) and the ADAMs, we also discovered that different sets of EBNA 3 proteins bind regulatory elements in a gene and cell-type specific manner. Binding profiles correlated with the effects of individual EBNA 3 proteins on the expression of these genes, providing a molecular basis for the targeting of different sets of cellular genes by the EBNA 3s. Our results therefore highlight the influence of the genomic and cellular context in determining the specificity of gene deregulation by EBV and provide a paradigm for host-cell reprogramming through modulation of enhancer-promoter interactions by viral transcription factors

Transcription factor 7-like 2 (TCF7L2) variant is associated with familial breast cancer risk: a case-control study

BACKGROUND: The transcription factor 7-like 2 (TCF7L2) is a critical component of the Wnt/β-catenin pathway. Aberrant TCF7L2 expression modifies Wnt signaling and mediates oncogenic effects through the upregulation of c-MYC and cyclin D. Genetic alterations in TCF7L2 may therefore affect cancer risk. Recently, TCF7L2 variants, including the microsatellite marker DG10S478 and the nearly perfectly linked SNP rs12233372, were identified to associate with type 2 diabetes. METHODS: We investigated the effect of the TCF7L2 rs12255372 variant on familial breast cancer (BC) risk by means of TaqMan allelic discrimination, analyzing BRCA1/2 mutation-negative index patients of 592 German BC families and 735 control individuals. RESULTS: The T allele of rs12255372 showed an association with borderline significance (OR = 1.19, 95% C.I. = 1.01-1.42, P = 0.04), and the Cochran-Armitage test for trend revealed an allele dose-dependent association of rs12255372 with BC risk (P(trend )= 0.04). CONCLUSION: Our results suggest a possible influence of TCF7L2 rs12255372 on the risk of familial BC

Molecular Mechanics of the α-Actinin Rod Domain: Bending, Torsional, and Extensional Behavior

α-Actinin is an actin crosslinking molecule that can serve as a scaffold and maintain dynamic actin filament networks. As a crosslinker in the stressed cytoskeleton, α-actinin can retain conformation, function, and strength. α-Actinin has an actin binding domain and a calmodulin homology domain separated by a long rod domain. Using molecular dynamics and normal mode analysis, we suggest that the α-actinin rod domain has flexible terminal regions which can twist and extend under mechanical stress, yet has a highly rigid interior region stabilized by aromatic packing within each spectrin repeat, by electrostatic interactions between the spectrin repeats, and by strong salt bridges between its two anti-parallel monomers. By exploring the natural vibrations of the α-actinin rod domain and by conducting bending molecular dynamics simulations we also predict that bending of the rod domain is possible with minimal force. We introduce computational methods for analyzing the torsional strain of molecules using rotating constraints. Molecular dynamics extension of the α-actinin rod is also performed, demonstrating transduction of the unfolding forces across salt bridges to the associated monomer of the α-actinin rod domain

Knowledge driven decomposition of tumor expression profiles

<p>Abstract</p> <p>Background</p> <p>Tumors have been hypothesized to be the result of a mixture of oncogenic events, some of which will be reflected in the gene expression of the tumor. Based on this hypothesis a variety of data-driven methods have been employed to decompose tumor expression profiles into component profiles, hypothetically linked to these events. Interpretation of the resulting data-driven components is often done by post-hoc comparison to, for instance, functional groupings of genes into gene sets. None of the data-driven methods allow the incorporation of that type of knowledge directly into the decomposition.</p> <p>Results</p> <p>We present a linear model which uses knowledge driven, pre-defined components to perform the decomposition. We solve this decomposition model in a constrained linear least squares fashion. From a variety of options, a lasso-based solution to the model performs best in linking single gene perturbation data to mouse data. Moreover, we show the decomposition of expression profiles from human breast cancer samples into single gene perturbation profiles and gene sets that are linked to the hallmarks of cancer. For these breast cancer samples we were able to discern several links between clinical parameters, and the decomposition weights, providing new insights into the biology of these tumors. Lastly, we show that the order in which the Lasso regularization shrinks the weights, unveils consensus patterns within clinical subgroups of the breast cancer samples.</p> <p>Conclusion</p> <p>The proposed lasso-based constrained least squares decomposition provides a stable and relevant relation between samples and knowledge-based components, and is thus a viable alternative to data-driven methods. In addition, the consensus order of component importance within clinical subgroups provides a better molecular characterization of the subtypes.</p

Electrostatic phase separation: a review

The current understanding and developments in the electrostatic phase separation are reviewed. The literature covers predominantly two immiscible and inter-dispersed liquids following the last review on the topic some 15 years. Electrocoalescence kinetics and governing parameters, such as the applied field, liquid properties, drop shape and flow, are considered. The unfavorable effects, such as chain formation and partial coalescence, are discussed in detail. Moreover, the prospects of microfluidics platforms, non-uniform fields, coalescence on the dielectric surfaces to enhance the electrocoalescence rate are also considered. In addition to the electrocoalescence in water-in-oil emulsions the research in oil-in-oil coalescence is also discussed. Finally the studies in electrocoalescer development and commercial devices are also surveyed. The analysis of the literature reveals that the use of pulsed DC and AC electric fields is preferred over constant DC fields for efficient coalescence; but the selection of the optimum field frequency a priori is still not possible and requires further research. Some recent studies have helped to clarify important aspects of the process such as partial coalescence and drop–drop non-coalescence. On the other hand, some key phenomena such as thin film breakup and chain formation are still unclear. Some designs of inline electrocoalescers have recently been proposed; however with limited success: the inadequate knowledge of the underlying physics still prevents this technology from leaving the realm of empiricism and fully developing in one based on rigorous scientific methodology

ER-Alpha-cDNA As Part of a Bicistronic Transcript Gives Rise to High Frequency, Long Term, Receptor Expressing Cell Clones

Within the large group of Estrogen Receptor alpha (ERα)-negative breast cancer patients, there is a subgroup carrying the phenotype ERα−, PR−, and Her2−, named accordingly “Triple-Negative” (TN). Using cell lines derived from this TN group, we wished to establish cell clones, in which ERα is ectopically expressed, forming part of a synthetic lethality screening system. Initially, we generated cell transfectants expressing a mono-cistronic ERα transcription unit, adjacent to a separate dominant selectable marker transcription unit. However, the yield of ERα expressing colonies was rather low (5–12.5%), and only about half of these displayed stable ectopic ERα expression over time. Generation and maintenance of such cell clones under minimal exposure to the ERα ligand, did not improve yield or expression stability. Indeed, other groups have also reported grave difficulties in obtaining ectopic expression of ERα in ERα-deficient breast carcinoma cells. We therefore switched to transfecting these cell lines with pERα-IRES, a plasmid vector encoding a bicistronic translation mRNA template: ERα Open Reading Frame (ORF) being upstream followed by a dominant-positive selectable marker (hygroR) ORF, directed for translation from an Internal Ribosome Entry Site (IRES). Through usage of this bicistronic vector linkage system, it was possible to generate a very high yield of ERα expressing cell clones (50–100%). The stability over time of these clones was also somewhat improved, though variations between individual cell clones were evident. Our successful experience with ERα in this system may serve as a paradigm for other genes where ectopic expression meets similar hardships

C-Terminal Region of EBNA-2 Determines the Superior Transforming Ability of Type 1 Epstein-Barr Virus by Enhanced Gene Regulation of LMP-1 and CXCR7

Type 1 Epstein-Barr virus (EBV) strains immortalize B lymphocytes in vitro much more efficiently than type 2 EBV, a difference previously mapped to the EBNA-2 locus. Here we demonstrate that the greater transforming activity of type 1 EBV correlates with a stronger and more rapid induction of the viral oncogene LMP-1 and the cell gene CXCR7 (which are both required for proliferation of EBV-LCLs) during infection of primary B cells with recombinant viruses. Surprisingly, although the major sequence differences between type 1 and type 2 EBNA-2 lie in N-terminal parts of the protein, the superior ability of type 1 EBNA-2 to induce proliferation of EBV-infected lymphoblasts is mostly determined by the C-terminus of EBNA-2. Substitution of the C-terminus of type 1 EBNA-2 into the type 2 protein is sufficient to confer a type 1 growth phenotype and type 1 expression levels of LMP-1 and CXCR7 in an EREB2.5 cell growth assay. Within this region, the RG, CR7 and TAD domains are the minimum type 1 sequences required. Sequencing the C-terminus of EBNA-2 from additional EBV isolates showed high sequence identity within type 1 isolates or within type 2 isolates, indicating that the functional differences mapped are typical of EBV type sequences. The results indicate that the C-terminus of EBNA-2 accounts for the greater ability of type 1 EBV to promote B cell proliferation, through mechanisms that include higher induction of genes (LMP-1 and CXCR7) required for proliferation and survival of EBV-LCLs

LRCH Proteins: A Novel Family of Cytoskeletal Regulators

Background: Comparative genomics has revealed an unexpected level of conservation for gene products across the evolution of animal species. However, the molecular function of only a few proteins has been investigated experimentally, and the role of many animal proteins still remains unknown. Here we report the characterization of a novel family of evolutionary conserved proteins, which display specific features of cytoskeletal scaffolding proteins, referred to as LRCHs. Principal Findings: Taking advantage of the existence of a single LRCH gene in flies, dLRCH, we explored its function in cultured cells, and show that dLRCH act to stabilize the cell cortex during cell division. dLRCH depletion leads to ectopic cortical blebs and alters positioning of the mitotic spindle. We further examined the consequences of dLRCH deletion throughout development and adult life. Although dLRCH is not essential for cell division in vivo, flies lacking dLRCH display a reduced fertility and fitness, particularly when raised at extreme temperatures. Conclusion/Significance: These results support the idea that some cytoskeletal regulators are important to buffer environmental variations and ensure the proper execution of basic cellular processes, such as the control of cell shape