Deep sequencing reveals clonal evolution patterns and mutation events associated with relapse in B-cell lymphomas

Background

Molecular mechanisms associated with frequent relapse of diffuse large B-cell lymphoma (DLBCL) are poorly defined. It is especially unclear how primary tumor clonal heterogeneity contributes to relapse. Here, we explore unique features of B-cell lymphomas - VDJ recombination and somatic hypermutation - to address this question.

Results

We performed high-throughput sequencing of rearranged VDJ junctions in 14 pairs of matched diagnosis-relapse tumors, among which 7 pairs were further characterized by exome sequencing. We identify two distinctive modes of clonal evolution of DLBCL relapse: an early-divergent mode in which clonally related diagnosis and relapse tumors diverged early and developed in parallel; and a late-divergent mode in which relapse tumors developed directly from diagnosis tumors with minor divergence. By examining mutation patterns in the context of phylogenetic information provided by VDJ junctions, we identified mutations in epigenetic modifiers such as KMT2D as potential early driving events in lymphomagenesis and immune escape alterations as relapse-associated events.

Conclusions

Altogether, our study for the first time provides important evidence that DLBCL relapse may result from multiple, distinct tumor evolutionary mechanisms, providing rationale for therapies for each mechanism. Moreover, this study highlights the urgent need to understand the driving roles of epigenetic modifier mutations in lymphomagenesis, and immune surveillance factor genetic lesions in relapse.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0432-0) contains supplementary material, which is available to authorized users.     

CIP2A oncoprotein controls cell growth and autophagy through mTORC1 activation

mTORC1 (mammalian target of rapamycin complex 1) integrates information regarding availability of nutrients and energy to coordinate protein synthesis and autophagy. Using ribonucleic acid interference screens for autophagy-regulating phosphatases in human breast cancer cells, we identify CIP2A (cancerous inhibitor of PP2A [protein phosphatase 2A]) as a key modulator of mTORC1 and autophagy. CIP2A associates with mTORC1 and acts as an allosteric inhibitor of mTORC1-associated PP2A, thereby enhancing mTORC1-dependent growth signaling and inhibiting autophagy. This regulatory circuit is reversed by ubiquitination and p62/SQSTM1-dependent autophagic degradation of CIP2A and subsequent inhibition of mTORC1 activity. Consistent with CIP2A’s reported ability to protect c-Myc against proteasome-mediated degradation, autophagic degradation of CIP2A upon mTORC1 inhibition leads to destabilization of c-Myc. These data characterize CIP2A as a distinct regulator of mTORC1 and reveals mTORC1-dependent control of CIP2A degradation as a mechanism that links mTORC1 activity with c-Myc stability to coordinate cellular metabolism, growth, and proliferation.     

Prevalence of class I–III BRAF mutations among 114,662 cancer patients in a large genomic database

BRAF mutations are relatively common in many cancers, particularly melanoma, colorectal cancer, and thyroid cancer and to a lesser extent in lung cancer. These mutations can be targeted by BRAF and MEK inhibitors, which exhibit good clinical activity. There are conflicting reports of the various relative rates of BRAF Class I mutations (V600 locus), defined as those that exhibit extremely strong kinase activity by stimulating monomeric activation of BRAF, Class II, define as non-V600 mutations that activate BRAF to signal as a RAS-independent dimer, and Class III mutations, defined as “kinase-dead” with low kinase activity as compared to wild type BRAF. Prospective studies have largely focused on patients with tumors harboring Class I BRAF mutations (limited to the V600 locus) where response rates up to 70% with BRAF plus MEK inhibition have been demonstrated. We report on the relative prevalence of various types of BRAF mutations across human cancers in a cohort of 114,662 patients that received comprehensive genomic profiling using next-generation sequencing. Of these patients, 4517 (3.9%) a pathogenic or presumed pathogenic BRAF mutation (3.9%). Of these, 1271 were seen in melanoma, representing 39.7% of all melanomas sequenced, representing the highest rate in all tumors. Class I (V600) mutations were seen overall in 2841 patients (62.1% of BRAF mutations, 2.4% of total cancers). Class II mutations were seen in 746 tumors (16.5% of BRAF mutant, 0.7% of total), and Class III mutations were seen in 801 tumors (17.7% of BRAF, 0.7% of total). Knowledge of the relative prevalence of these types of mutations can aid in the development of agents that might better address non-V600 mutations in cancer.Impact statementThese data represent the largest aggregation of BRAF mutations within a single clinical database to our knowledge. The relative proportions of both BRAF V600 mutations and non-V600 mutations are informative in all cancers and by malignancy, and can serve as a definitive gold-standard for BRAF mutation cancer incidence by malignancy. The rate of BRAF mutation in human cancer in a real-world large database is lower than previously reported likely representing testing more broadly across tumor types. The relative percentages of Class II and Class III BRAF mutations are higher than previously reported, representing almost 35% of BRAF mutations in cancer. These findings provide support for the development of effective treatments for non-V600 BRAF mutations in cancer.     

Mammalian and Malaria Parasite Cyclase-associated Proteins Catalyze Nucleotide Exchange on G-actin through a Conserved Mechanism

Cyclase-associated proteins (CAPs) are among the most highly conserved regulators of actin dynamics, being present in organisms from mammals to apicomplexan parasites. Yeast, plant, and mammalian CAPs are large multidomain proteins, which catalyze nucleotide exchange on actin monomers from ADP to ATP and recycle actin monomers from actin-depolymerizing factor (ADF)/cofilin for new rounds of filament assembly. However, the mechanism by which CAPs promote nucleotide exchange is not known. Furthermore, how apicomplexan CAPs, which lack many domains present in yeast and mammalian CAPs, contribute to actin dynamics is not understood. We show that, like yeast Srv2/CAP, mouse CAP1 interacts with ADF/cofilin and ADP-G-actin through its N-terminal α-helical and C-terminal β-strand domains, respectively. However, in the variation to yeast Srv2/CAP, mouse CAP1 has two adjacent profilin-binding sites, and it interacts with ATP-actin monomers with high affinity through its WH2 domain. Importantly, we revealed that the C-terminal β-sheet domain of mouse CAP1 is essential and sufficient for catalyzing nucleotide exchange on actin monomers, although the adjacent WH2 domain is not required for this function. Supporting these data, we show that the malaria parasite Plasmodium falciparum CAP, which is entirely composed of the β-sheet domain, efficiently promotes nucleotide exchange on actin monomers. Collectively, this study provides evidence that catalyzing nucleotide exchange on actin monomers via the β-sheet domain is the most highly conserved function of CAPs from mammals to apicomplexan parasites. Other functions, including interactions with profilin and ADF/cofilin, evolved in more complex organisms to adjust the specific role of CAPs in actin dynamics.     

The selective regulation of alpha Vbeta 1 integrin expression is based on the hierarchical formation of alpha V-containing heterodimers

The integrin beta1 subunit can form a heterodimer with 12 different alpha subunits. According to the present model, the expression level of any alphabeta complex is regulated by the availability of the specific alpha subunit, whereas beta1 subunit is constantly present in a large excess. The expression of several heterodimers containing the alphaV subunit seems to be regulated by an identical mechanism. The fact that many cells express alphaVbeta1 heterodimer, and that this fibronectin/vitronectin receptor may be selectively regulated, compromises the present model of the regulation of beta1 and alphaV integrins. We have tried to solve this problem by assuming that distinct alphabeta heterodimers are formed with different tendency. To test the hypothesis, we analyzed WM-266-4 melanoma cells transfected with a cDNA construct coding for an intracellular single-chain anti-alphaV integrin antibody. We could see 70-80% reduction in the cell surface expression of alphaV subunit. However, the only one of the alphaV integrins reduced on the cell surface was alphaVbeta1. This suggests that the cell surface expression level of alphaVbeta1 is dependent on the number of alphaV subunits available after the formation of other alphaV-containing heterodimers. Thus, there seems to be a hierarchy in the complex formation between alphaV and its different beta-partners. These observations explain how alphaVbeta1 can be specifically regulated without concomitant changes in the expression of other alphaV or beta1 integrins.     

Mating success in lekking males: a meta-analysis

Traits that are correlated with mating success are likely tobe subject to sexual selection. In lekking species, a male'smating success can be estimated as the number of females thathe copulates with. Earlier reviews of sexual selection in lekkingspecies have been inconclusive, suggesting that different traitsmay be important in different species. To obtain a more completeunderstanding of the outcome of sexual selection in this matingsystem, we performed a meta-analysis in which we combined theresults from different studies across a wide variety of taxa.Our aim was to synthesize available information about correlatesof male mating success in lekking species. We found that behavioraltraits such as male display activit aggression rate, and lekattendance were positively correlated with male mating success.Further, territory position was negatively correlated with malemating success, such that males with territories close to thegeometric center of the leks had higher mating success thanother males. The size of "extravagant" traits, such as birdstails and ungulate antlers, and age were positively correlatedwith male mating success. Male morphology (measure of body size)and territory size showed small effects on male mating success.Our results confirm some of the suggestions put forward by earlierreviews but add more rigor to the condusions drawn. Part ofthe variation across studies still remain unaccounted for. Furtherstudies are needed to perform proper meta-analyses that cantake factors like phylogeny and sexual dimorphism into account.     

Engagement of overexpressed Her2 with GEP100 induces autonomous invasive activities and provides a biomarker for metastases of lung adenocarcinoma

Overexpression of Her2/ErbB2/Neu in cancer is often correlated with recurrent distant metastasis, although the mechanism still remains largely elusive. We have previously shown that EGFR, when tyrosine-phosphorylated, binds to GEP100/BRAG2 to activate Arf6, which induces cancer invasion and metastasis. We now show that overexpressed Her2 in lung adenocarcinoma cells also employs GEP100. Like EGFR-GEP100 binding, this association is primarily mediated by the pleckstrin homology (PH) domain of GEP100 and Tyr1139/Tyr1196 of Her2. Tyr1139/Tyr1196 are autonomously phosphorylated, when Her2 is overexpressed. Accordingly, invasive activities mediated by the Her2-GEP100 pathway are not dependent on external factors. Blocking Her2-GEP100 binding, as well as its signaling pathway all inhibit cancer invasive activities. Moreover, our clinical study indicates that co-overexpression of Her2 with GEP100 in primary lung adenocarcinomas of patients is correlated with the presence of their node-metastasis with a statistical significance. Since the GEP100 PH domain interacts with both Her2 and EGFR, targeting this domain may provide novel cancer therapeutics.     

The Schizosaccharomyces pombe Pzh1 protein phosphatase regulates Na+ ion influx in a Trk1-independent fashion.

We have previously shown that fission yeast encodes a PPZ-like phosphatase, designated Pzhl, which is an important determinant of cation homeostasis. pzh1 delta mutants display increased tolerance to Na+ ions, but they are hypersensitive to KC1 [Balcells, L., Gómez, N., Casamayor, A., Clotet, J. & Ari?o, J. (1997) Eur. J. Biochem. 250, 476-483]. We have immunodetected Pzh1 in yeast extracts and found that this phosphatase is largely associated with particulate fractions. Cells defective in Pzh1 do not show altered efflux of Na+ or Li+ ions, but they accumulate these cations more slowly than wild-type cells. K+ ion content of pzh1 delta cells is about twice that of wild-type cells, and this can be explained by decreased efflux of K+. Therefore, Pzh1 may regulate both Na+ influx and K+ efflux in fission yeast. To test the possible relationship between K+ uptake, Na+ tolerance and Pzh1 function, we deleted the trk1+ gene, which encodes a putative high-affinity transporter of K+ ions. trkl delta mutants grew well even at relatively low concentrations of KCl and did not show significantly altered content or influx of K+ ions. However, they showed a Na(+)-sensitive phenotype which was greatly intensified by deletion of the sod2+ gene (which encodes the major determinant for efflux of Na+ ions), and clearly ameliorated by deletion of the pzh1 phosphatase, as well as by moderate concentrations of KCl in the medium. These results suggest that Trk1 does not mediate the effect of Pzh1 on NaCl tolerance and that fission yeast contains efficient systems, other than Trk1, for uptake of K+ ions.     

Independent hybrid populations of <Emphasis Type="Italic">Formica polyctena X</Emphasis><Emphasis Type="Italic">rufa</Emphasis> wood ants (Hymenoptera: Formicidae) abound under conditions of forest fragmentation

Combined genetic and morphological data indicate frequent hybridisation between the wood ants Formica polyctena Förster 1850 and F. rufa Linnaeus 1761 in Central Europe. The genetic and morphological traits give a concordant picture of hybridisation with a strong correlation between the genotypic admixture proportions at 19 microsatellite loci and the first vectors of a principal component analysis (P < 0.001) and of a 3-class discriminant analysis (P < 0.001) of 15 quantitative morphological characters. This integrative approach enabled a grouping into F. polyctena, the hybrid and rufa. Genetic differentiation between the hybrid and F. rufa is significantly larger than between the hybrid and polyctena, indicating gene flow mainly between the latter entities. A suggested gene flow bias towards F. polyctena agrees with differential queen acceptance and mating behaviour. Both genetic and phenotypic colony parameters indicate predominance of monogyny in F. rufa but of polygyny in polyctena and the hybrid. Hybrids are intermediate between the parental species in body size, diagnostic morphological characters, monogyny frequency, size of nest population, nest diameter and infestation rate with epizootic fungi. The three entities respond differently to woodland fragmentation. Hybrids are significantly more abundant in forests with a coherent area <300 ha than in woodland above this size. Regions with high hybrid frequency in Germany—the Eastern Oberlausitz (23%) and the Baltic Sea islands Darss, Hiddensee and Rügen (28%)—are characterised by a fragmented woodland structure whereas regions with low hybrid frequency—Brandenburg and the lower Erzgebirge (3.4%)—have clearly larger and more coherent forest systems. Data from other European countries indicate habitat fragmentation to be a facilitating factor but no essential precondition for interspecific hybridisation in these ants. Hybrids are hypothesised to have selective advantage in fragmented systems because of combining the main reproductive and dispersal strategies of the parental species.