Sophia Yancopoulos

Efficient sorting of genomic permutations by translocation, inversion and block interchange

MOTIVATION Finding genomic distance based on gene order is a classic problem in genome rearrangements. Efficient exact algorithms for genomic distances based on inversions and/or translocations have been found but are complicated by special cases, rare in simulations and empirical data. We seek a universal operation underlying a more inclusive set of evolutionary operations and yielding a tractable genomic distance with simple mathematical form. RESULTS We study a universal double-cut-and-join operation that accounts for inversions, translocations, fissions and fusions, but also produces circular intermediates which can be reabsorbed. The genomic distance, computable in linear time, is given by the number of breakpoints minus the number of cycles (b-c) in the comparison graph of the two genomes; the number of hurdles does not enter into it. Without changing the formula, we can replace generation and re-absorption of a circular intermediate by a generalized transposition, equivalent to a block interchange, with weight two. Our simple algorithm converts one multi-linear chromosome genome to another in the minimum distance.

Intraclonal complexity in chronic lymphocytic leukemia: fractions enriched in recently born/divided and older/quiescent cells.

The failure of chemotherapeutic regimens to eradicate cancers often results from the outgrowth of minor subclones with more dangerous genomic abnormalities or with self-renewing capacity. To explore such intratumor complexities in B-cell chronic lymphocytic leukemia (CLL), we measured B-cell kinetics in vivo by quantifying deuterium (2H)-labeled cells as an indicator of a cell that had divided. Separating CLL clones on the basis of reciprocal densities of chemokine (C-X-C motif) receptor 4 (CXCR4) and cluster designation 5 (CD5) revealed that the CXCR4dimCD5bright (proliferative) fraction contained more 2H-labeled DNA and hence divided cells than the CXCR4brightCD5dim (resting) fraction. This enrichment was confirmed by the relative expression of two cell cycle-associated molecules in the same fractions, Ki-67 and minichromosome maintenance protein 6 (MCM6). Comparisons of global gene expression between the CXCR4dimCD5bright and CXCR4brightCD5dim fractions indicated higher levels of pro-proliferation and antiapoptotic genes and genes involved in oxidative injury in the proliferative fraction. An extended immunophenotype was also defined, providing a wider range of surface molecules characteristic of each fraction. These intraclonal analyses suggest a model of CLL cell biology in which the leukemic clone contains a spectrum of cells from the proliferative fraction, enriched in recently divided robust cells that are lymphoid tissue emigrants, to the resting fraction enriched in older, less vital cells that need to immigrate to lymphoid tissue or die. The model also suggests several targets preferentially expressed in the two populations amenable for therapeutic attack. Finally, the study lays the groundwork for future analyses that might provide a more robust understanding of the development and clonal evolution of this currently incurable disease.

Identification of outcome-correlated cytokine clusters in chronic lymphocytic leukemia

Individual cytokines and groups of cytokines that might represent networks in chronic lymphocytic leukemia (CLL) were analyzed and their prognostic values determined. Serum levels of 23 cytokines were measured in 84 patients and 49 age-matched controls; 17 levels were significantly elevated in patients. Unsupervised hierarchical bicluster analysis identified 3 clusters (CLs) of highly correlated but differentially expressed cytokines: CL1 (CXCL9, CXCL10, CXCL11, CCL3, CCL4, CCL19, IL-5, IL-12, and IFNγ), CL2 (TNFα, IL-6, IL-8, and GM-CSF), and CL3 (IL-1β, IL-2, IL-4, IL-15, IL-17, and IFNα). Combination scores integrating expression of CL1/CL2 or CL1/CL3 strongly correlated (P < .005) with time-to-first-treatment and overall survival (OS), respectively. Patients with the worst course had high CL1 and low CL2 or CL3 levels. Multivariate analysis revealed that CL1/CL2 combination score and immunoglobulin heavy chain variable region mutation status were independent prognostic indicators for time-to-first-treatment, whereas CL1/CL3 combination score and immunoglobulin heavy chain variable region mutation status were independent markers for OS. Thus, we identified groups of cytokines differentially expressed in CLL that are independent prognostic indicators of aggressive disease and OS. These findings indicate the value of multicytokine analyses for prognosis and suggest therapeutic strategies in CLL aimed at reducing CL1 and increasing CL2/CL3 cytokines.

Characterization of structurally defined epitopes recognized by monoclonal antibodies produced by chronic lymphocytic leukemia B cells

Despite a wealth of information about the structure of surface membrane immunoglobulin (smIg) on chronic lymphocytic leukemia (CLL) cells, little is known about epitopes reacting with their binding sites. Probing phage-displayed peptide libraries, we identified and characterized mimetopes for Igs of 4 patients with IGHV mutated CLL (M-CLL) and 4 with IGHV unmutated CLL (U-CLL). Six of these mAbs were representatives of stereotyped B-cell receptors characteristic of CLL. We found that mimetic epitopes for U- and M-CLL Igs differed significantly. M-CLL-derived peptides exhibited better amino acid motifs, were more similar to each other, aligned more easily, and formed tighter clusters than U-CLL-derived peptides. Mono-, oligo-, and polyreactivity of peptides correlated with structural changes within antigen-binding sites of selecting M-CLL mAbs. Although M-CLL-isolated peptides and certain U-CLL mAbs bound more effectively to the selecting mAb, others were not as specific, reacting with M-CLL and U-CLL mAbs; these data suggest that in vivo structurally diverse epitopes could bind smIgs of distinct CLL clones, thereby altering survival and growth. Finally, an M-CLL-derived peptide inhibited, in a dose-dependent manner, binding of its homologous mAb to human B lymphocytes; therefore peptides that inhibit or alter the consequences of antigen-smIg interactions may represent therapeutic modalities in CLL.

Sorting Genomes with Insertions, Deletions and Duplications by DCJ

We extend the DCJ paradigm to perform genome rearrangments on pairs of genomes having unequal gene content and/or multiple copies by permitting genes in one genome which are completely or partially unmatched in the other. The existence of unmatched gene ends introduces new kinds of paths in the adjacency graph, since some paths can now terminate internal to a chromosome and not on telomeres. We introduce Oghost adjacenciesO to supply the missing gene ends in the genome not containing them. Ghosts enable us to close paths that were due to incomplete matching, just as null points enable us to close even paths terminating in telomeres. We define generalalized DCJ operations on the generalized adjacency graph, and give a prescription for calculating the DCJ distance for the expanded repertoire of operations which includes insertions, deletions and duplications.

Autoantigen can promote progression to a more aggressive TCL1 leukemia by selecting variants with enhanced B-cell receptor signaling

Significance These studies indicate that autoantigen-reactivity plays a role in the progression of a murine leukemia that models human chronic lymphocytic leukemia. This indication is consistent with the notion that chronic lymphocytic leukemia evolves by selection of normal B cells that bind autoantigen via the B-cell antigen receptor. (Auto)antigen engagement by the B-cell receptor (BCR) and possibly the sites where this occurs influence the outcome of chronic lymphocytic leukemia (CLL). To test if selection for autoreactivity leads to increased aggressiveness and if this selection plays out equally in primary and secondary tissues, we used T-cell leukemia (TCL)1 cells reactive with the autoantigen phosphatidylcholine (PtC). After repeated transfers of splenic lymphocytes from a single mouse with oligoclonal PtC-reactive cells, outgrowth of cells expressing a single IGHV-D-J rearrangement and superior PtC-binding and disease virulence occurred. In secondary tissues, increased PtC-binding correlated with enhanced BCR signaling and cell proliferation, whereas reduced signaling and division of cells from the same clone was documented in cells residing in the bone marrow, blood, and peritoneum, even though cells from the last site had highest surface membrane IgM density. Gene-expression analyses revealed reciprocal changes of genes involved in BCR-, CD40-, and PI3K-signaling between splenic and peritoneal cells. Our results suggest autoantigen-stimulated BCR signaling in secondary tissues promotes selection, expansion, and disease progression by activating pro-oncogenic signaling pathways, and that—outside secondary lymphoid tissues—clonal evolution is retarded by diminished BCR-signaling. This transferrable, antigenic-specific murine B-cell clone (TCL1-192) provides a platform to study the types and sites of antigen-BCR interactions and genetic alterations that result and may have relevance to patients.

Indole-3-carbinol improves survival in lupus-prone mice by inducing tandem B- and T-cell differentiation blockades.

Indole-3-carbinol (I3C), derived from cruciferous vegetables, alters estrogen metabolism. Since lupus is estrogen dependent, we reasoned that I3C might be effective in SLE. I3C significantly thwarted disease progression and prolonged survival in (NZBxNZW) F1 mice. Immunofluorescent and serologic analyses in treated animals indicated a transient blockade in B-cell maturation with increased immature B cells, decreased mature B cells, and a significant reduction of certain autoantibodies. Subsequently, a delay in T-cell maturation occurred in the treated group, manifested by significantly increased naive T cells, decreased mature and memory T cells, and decreased CD4:CD8 T-cell ratios. T cells from the I3C cohort, stimulated in vitro with various mitogens, exhibited enhanced responsiveness. Con A-stimulated T cells from I3C-treated mice produced Th1 cytokines, whereas those from control animals produced Th2 cytokines. Our studies suggest immunological mechanisms by which I3C ameliorates SLE in mice and provide a rationale for its use as an adjunctive therapy for human lupus.

Surface expression of Bcl-2 in chronic lymphocytic leukemia and other B-cell leukemias and lymphomas without a breakpoint t(14;18).

Since its discovery in follicular lymphoma cells at the breakpoint t(14;18), Bcl-2 has been studied extensively in many basic and clinical science settings. Bcl-2 can locate as an integral mitochondrial membrane component, where its primary role is to block apoptosis by maintaining membrane integrity. Here we show that Bcl-2 also can position on the outer cell surface membrane of B cells from patients with chronic lymphocytic leukemia (B-CLL) and certain other leukemias that do not classically possess the chromosomal breakpoint t(14;18).Although low levels of Bcl-2 can be detected on the surface membrane of apparently healthy leukemic and normal B cells, expression of Bcl-2 correlates best with spontaneous or induced apoptosis. Notably, upon induction of apoptosis, B-CLL cells were much more efficient in upregulating surface Bcl-2 than normal B cells. It is not clear if this surface membrane expression is a passive consequence of the apoptotic process or an active attempt by the B cell to abort cell death by stabilizing the plasma membrane.

Genome Rearrangement by the Double Cut and Join Operation

The Double Cut and Join is an operation acting locally at four chromosomal positions without regard to chromosomal context. This chapter discusses its application and the resulting menu of operations for genomes consisting of arbitrary numbers of circular chromosomes, as well as for a general mix of linear and circular chromosomes. In the general case the menu includes: inversion, translocation, transposition, formation and absorption of circular intermediates, conversion between linear and circular chromosomes, block interchange, fission, and fusion. This chapter discusses the well-known edge graph and its dual, the adjacency graph, recently introduced by Bergeron et al. Step-by-step procedures are given for constructing and manipulating these graphs. Simple algorithms are given in the adjacency graph for computing the minimal DCJ distance between two genomes and finding a minimal sorting; and use of an online tool (Mauve) to generate synteny blocks and apply DCJ is described.

The Emperor Has No Caps! A Comparison of DCJ and Algebraic Distances

In this chapter we investigate the DCJ and algebraic distances and how they are found. We introduce a new graphical method to determine the permutation cycles which embody the composition permutation for the genome transformation in the algebraic method. This graphical method helps tie the two approaches together. In the usual approaches, the two methods differ only in the distance component due to the even paths in the adjacency graph of Bergeron, Mixtacki, and Stoye involving operations changing type and number of chromosomes, such as fission, fusion, altering chromosome type from circular to linear, and vice versa. Discussing each distance individually, we compare their underlying assumptions. Both methods resort to cycles to determine the distance, but the basic DCJ uses “caps” to close paths. Without caps the algebraic distance differs from the standard DCJ for even paths. However, if caps and null chromosomes are added, the weighting schemes agree. A convention which can be done in multiple ways is the method of path closure. We discuss implementation of the original closure rule to arrive at the usual weighting scheme for the DCJ. Instead, by a new alternative closure rule which we introduce, the distance diverts to the algebraic distance. Finally, we note that although the Bergeron, Mixtacki, and Stoye DCJ approach via the adjacency graph does away with “fictitious” caps and nulls, vestiges of fictitious operations may remain, as the resulting weighting scheme is equivalent to that of the basic DCJ.