Introduction {#S1}
============

G-protein coupled receptors (GPCRs) represent the largest class of cell surface receptors and when activated by extracellular ligands, they transmit signals involved in cell growth, proliferation, survival, differentiation and motility ([@R1]). GPCRs initiate signal transduction cascades by coupling to and activating heterotrimeric G-proteins. Depending on the coupling specificity, a given GPCR may couple to one or more different Gα proteins, which are grouped into four main families: Gα~12~ (including Gα~13~), Gα~s~, Gα~i~, and Gα~q~ (including Gα~11~). Each family of G-proteins activates distinct sets of second messenger and kinase signaling cascades ([@R1]).

Heterotrimeric G-proteins are major cellular signaling hubs, and they are activated upon binding to GPCRs or non-receptor guanine nucleotide exchange factors (GEFs) ([@R2]). Due to their important roles in regulating cell survival, proliferation and movement, it is not surprising that tumors often harbor mutations in or exhibit aberrant expression of G-proteins and their GEFs ([@R3]). For example, activating mutations in *GNA11* and *GNAQ*, encoding Gα~11~ and Gα~q~ respectively, were recently discovered to be the key oncogenic drivers in uveal melanoma ([@R4]--[@R6]). Gα~s~ is the most frequently mutated G-protein in human cancers and activating mutations in the gene, *GNAS,* have been found in a variety of neoplasms including pituitary, thyroid, pancreatic, biliary tract, colon and small intestine and a variety of other tumors ([@R3]). Furthermore, constitutively active mutants of genes encoding Gα~i~, Gα~o~, Gα~q~ Gα~12~, and Gα~13~ were found to induce cellular transformation in experimental systems (reviewed in ([@R1], [@R7])).

Despite the transforming capacity of constitutive Gα~12~ and Gα~13~ activity in experimental systems and numerous implications of this G-protein family and downstream targets in cancer metastasis ([@R8]--[@R13]), activating mutations in the *GNA13* and *GNA12* genes in patient tumor samples have not been described. However, recent large-scale sequencing efforts have revealed the presence of *GNA13* mutations in Burkitt's lymphoma and Diffuse Large B cell Lymphoma (DLBCL) ([@R14], [@R15]). Interestingly, recent studies in mouse models demonstrated that conditional B cell deficiency in Gα~13~ or the Gα~13~ --coupled sphingosine 1 phosphate receptor 2 (S1P2) result in DLBCL-like phenotypes ([@R16], [@R17]). Based on the analysis of deposited sequencing data from tumors in the Catalog of Somatic Mutations in Cancer (COSMIC), mutations in *GNA13* in human Burkitt's lymphoma and DLBCL are highly statistically significant over background cancer mutation rate, with p-value and q-value scores of 0 ([@R3]), suggesting these mutations are likely not random. However, unlike the activating GTPase domain mutations found in other G-proteins in cancers, including Gα~q~ and Gα~s~, the mutations in *GNA13* are distributed throughout the gene. Furthermore, we identified additional mutations in the gene of the major downstream effector of Gα~13~ signaling, RhoA.

In this study, we characterize the mutations identified in *GNA13* and *RHOA* in Burkitt's lymphoma and DLBCL tumor samples to determine how these mutations affect protein function and signaling capacity. We also evaluated the effects of mutations and wild type (WT) Gα~13~ expression on tumor growth and progression in xenograft models. Overall, our results support a tumor suppressive role for the Gα~13~/RhoA axis in Burkitt's lymphoma and DLBCL. Our data also extend recent findings supporting the presence of disruptive *RHOA* mutations in peripheral T cell lymphomas, suggesting that disruption of RhoA function may have a broad impact in multiple haematological malignancies ([@R18]--[@R22]).

Results {#S2}
=======

Mutations in GNA13 and RHOA are frequent in Burkitt's Lymphoma and DLBCL tumors {#S3}
-------------------------------------------------------------------------------

Data from genome-wide sequencing analyses collected through the Catalog of Somatic Mutations in Cancer (COSMIC v72) database reveal the presence of *GNA13* mutations in nearly 2% of all haematopoietic and lymphoid malignancies ([Figure 1A](#F1){ref-type="fig"}). Previous statistical analyses of these mutations indicated p-value and q-value (for false discovery rate) scores approaching 0, suggesting they are unlikely to be random, but rather could have important driver mutation functions ([@R3]). Of the haematopoietic and lymphoid malignancies evaluated in COSMIC, most of the *GNA13* mutations are present in B cell lymphomas, primarily Diffuse Large B cell lymphoma (DLBCL) and Burkitt's Lymphoma, for which mutations are harbored in approximately 10% of patient tumor samples ([Figure 1A](#F1){ref-type="fig"}). Mutations in *GNA13* found in Burkitt's lymphoma and DLBCL appeared likely to result in loss of function because nearly 22% (5/23) of the DLBCL mutations (17% of overall in both lymphomas) were nonsense, resulting in a premature STOP codon and all other mutations were non-synonymous ([Figure 1A](#F1){ref-type="fig"}). When mapped onto the crystal structure of *GNA13*, the mutations clustered into multiple regions including the binding interface for the Gβ/γ subunits, the helical domain, conformational switches 1 and 2 (important for nucleotide binding as well as interaction with regulators and effectors), conformational switch 3, the GPCR interface, and directly in the nucleotide pocket ([Figures 1A](#F1){ref-type="fig"}, [Supplemental Table 1](#SD1){ref-type="supplementary-material"}). An interesting in-frame deletion was observed in the loop connecting the GTPase and the helical domain that can dramatically impair the relative mobility of the two domains. Mutations are therefore not localized into one region and so they could exert a broad influence on the protein function or stability.

RhoA is the main downstream effector of Gα~13~ signaling ([@R23]--[@R25]), therefore we searched the COSMIC database to determine if there were mutations in *RHOA* in Burkitt's Lymphoma and DLBCL tumors as well. Although less frequent than the *GNA13* mutations, non-redundant *RHOA* mutations were observed in nearly 6% of Burkitt's Lymphoma and almost 1% of DLBCL tumors ([Figure 1B](#F1){ref-type="fig"}), as these *RHOA* mutations were present in separate tumor samples from those with *GNA13* mutations. Taken together, the Gα~13~/RhoA axis is disrupted in more than 10% of these B cell lymphomas. Interestingly, there was a recurrent R5Q mutation, observed in 7/16 *RHOA* mutations identified in these tumors. Structural analysis of RhoA suggests that R5 is important for GEF binding. Several of the other mutations in *RHOA* also localized to GEF binding residues, indicating that these mutations could influence RhoA activation due to changes in RhoA-GEF interactions ([Figure 1B](#F1){ref-type="fig"}, [Supplemental Table 2](#SD1){ref-type="supplementary-material"}).

Characterization of GNA13 mutations indicates loss of function {#S4}
--------------------------------------------------------------

In order to determine the consequences of the *GNA13* mutations on Gα~13~ protein function, we prepared and characterized 9 of the mutations from the COSMIC data repository, representing a variety of the structural and functional domains. Plasmids encoding empty vector, Gα~13~ WT, Gα~13~ mutants based on Burkitt's lymphoma and DLBCL tumor samples and an established constitutively active Q226L mutant were transfected into HEK293 cells and Gα~13~ protein expression was assessed by western blot ([Figure 2A](#F2){ref-type="fig"}). Although most of the Gα~13~ mutants were expressed at similar levels to WT, there were a few including L197Q and F245S that had reduced protein expression. Because the Gα~13~ antibody recognizes an epitope in the C-terminus of Gα~13~, disrupted antibody recognition is unlikely to fully explain the relatively low protein expression observed in the L197Q and F245S mutants. Rather, it could be that these mutants affect protein stability, thus reducing expression.

To evaluate the functional consequences of these mutants, we tested their activity in a serum response factor response element (SRF-RE) luciferase assay, which is used to selectively monitor gene transcription downstream of Gα~13~-RhoA signaling. The SRF-RE luciferase assay demonstrated that all 9 Gα~13~ mutants had significantly reduced basal transcriptional activity relative to WT and constitutively active (Q226L) Gα~13~ ([Figure 2B](#F2){ref-type="fig"}). Due to promiscuity of most Gα~13~-coupled GPCRs, which often also couple to other G proteins and/or have ligands that target other GPCRs, a synthetic biology approach was employed to determine the effects of the Gα~13~ mutants on ligand-activated Gα~13~ signaling events. The synthetic GPCR (SyR), also known as DREADD ([@R26]), is activated by an otherwise biologically inert ligand, clozapine N-oxide (CNO). Because no Gα~13~-SyR is currently available, an SyR that couples to Gα~i~ was used in combination with a chimeric G-protein that has Gα~i~-coupling specificity (based on the last 5 amino acids at the C-terminus) but retains Gα~13~ signaling response. Additionally, in order to ensure that the responses observed were due to Gα~13~ activity and not influenced by potential Gα~i~ signaling, a C-I mutation was made in the last 5 amino acids of the Gα~13-i5~ chimeras to make them pertussis toxin (PTX)-insensitive and the assay was performed in the presence of PTX. The SRF-RE assay demonstrated that some of the Gα~13~ mutants were completely dead whereas others still retained some ligand-activated response ([Figure 2C](#F2){ref-type="fig"}). Nevertheless, overall SRF-RE activity was lower in all the mutants compared with WT (and basal activity was again lower in all the mutants). A recently developed TGF-α shedding assay was used as another measurement of Gα~13~ activity ([@R27]). Gα~13~ and Gα~q~ activity are known to induce the shedding of tethered TGF-α, and so we could monitor the release of AP-TGF-α into media as a measurement of Gα~13~ signaling ([Figure 3A](#F3){ref-type="fig"}). Similar to the SRF-RE results, the spontaneous (or basal) activities of all the Gα~13~ mutants were significantly lower than WT or the Q226L constitutively active mutant ([Figure 3B](#F3){ref-type="fig"}, [Supplemental Table 3](#SD1){ref-type="supplementary-material"}), whereas a few of the mutants still retained some ligand-induced signaling response as determined using the synthetic receptor system, SyR-Gi with the Gα~13-i5~ chimeras, ([Figure 3C](#F3){ref-type="fig"}, [Supplemental Table 3](#SD1){ref-type="supplementary-material"}) and using another receptor that couples to Gα~13~, the dopamine D2 receptor ([Supplemental Figure 1](#SD2){ref-type="supplementary-material"}, [Supplemental Table 3](#SD1){ref-type="supplementary-material"}).

GNA13 and RHOA mutations impair RhoA activity {#S5}
---------------------------------------------

RhoA is the main downstream effector of Gα~13~ signaling, thus we sought to analyze the effects of the *GNA13* mutations on RhoA activity and to determine the effect of the recurrent RhoA R5Q mutation. We monitored active GTP-bound RhoA using a pull-down assay for three of the Gα~13~ mutants and controls. In line with the previous experiments, the Gα~13~ mutants had considerably reduced levels of active RhoA compared to WT Gα~13~, constitutively active Gα~13~ Q226L and the thrombin-stimulated positive control ([Figure 4A](#F4){ref-type="fig"}). Next, we evaluated the activity of the R5Q RhoA mutant that had been identified in several DLBCL and Burkitt's lymphoma tumors. The R5Q mutant exhibited reduced activity compared with WT RhoA and the constitutively active Q63L RhoA mutant in both an active RhoA pull down assay ([Figure 4B](#F4){ref-type="fig"}) and an SRF-RE luciferase assay ([Figure 4C](#F4){ref-type="fig"}). However, the RhoA R5Q mutant appeared to retain a little activity compared with the dominant negative T19N RhoA mutant ([Figure 4C](#F4){ref-type="fig"}). Formation of actin stress fibers is often used as another indication of RhoA activity, so we performed a stress fiber formation assay using LifeActin-GFP in MDCK cells. Consistent with the SRF and RhoA pull down results, we observed significant actin stress fiber formation in WT and Q63L constitutively active RhoA-transfected MDCK cells, but not with the R5Q mutant or the dominant negative T19N mutant ([Figure 4D](#F4){ref-type="fig"}). Overall, these data indicate that the *GNA13* and *RHOA* mutations observed in DLBCL and Burkitt's lymphoma tumors result in loss of function of RhoA activity and downstream signaling events.

Restoration of Gα~13~ WT suppresses tumor growth {#S6}
------------------------------------------------

The COSMIC sequencing data indicated the presence of a *GNA13* mutation at L184R in the frequently used Raji Burkitt's lymphoma cell line. Therefore, in order to determine if mutation and thus loss of Gα~13~ activity provides the cells with a growth or survival advantage, we expressed Gα~13~ WT in the cells. First we verified the COSMIC sequencing data in the Raji cell line and confirmed the presence of a heterozygous L184R mutation ([Figure 5A](#F5){ref-type="fig"}). Then, we transduced Raji cells using MSCV puro IRES GFP retrovirus with Gα~13~ WT or vector control. After cell sorting, we were able to achieve a highly GFP-enriched population of cells and verified the expression of Gα~13~ ([Figure 5B](#F5){ref-type="fig"}). B cell lines are notorious for their difficulties in terms of feasibility of genetic manipulations, which prevented us from deploying similar or complementary strategies to express Gα~13~ WT in additional relevant cellular systems. Regarding Gα~13~ function, it was previously reported that Gα~13~ activity has suppressive effects on Akt phosphorylation ([@R17]). Thus, we measured basal Akt phosphorylation at S473 and found that Gα~13~ WT expression in the Raji cells reduced phosphorylated Akt ([Figure 5C](#F5){ref-type="fig"}). Using *in vitro* viability and proliferation assays, we did not observe any differences between the control and Gα~13~ WT-expressing Raji cells (data not shown). In a methocellulose clonogenic assay, WT Gα~13~--expressing cells showed only a slight reduction in the number of colonies; however, the difference in cell morphology was striking ([Figures 5D--E](#F5){ref-type="fig"}). While the majority of the control Raji cells formed spheroid colonies in the methocellulose, the Raji cells expressing Gα~13~ WT tended to exhibit a more flattened morphology ([Figure 5E](#F5){ref-type="fig"}). Next we tested the Raji control cells and Raji Gα~13~ WT-transduced cells in a tumor xenograft model using NSG (NOD SCID gamma) mice. In line with the role of a tumor suppressor, the restoration of Gα~13~ WT to the Raji cells resulted in significantly delayed and impaired tumor growth and smaller tumor masses and volumes ([Figures 5F--H](#F5){ref-type="fig"}). Histological evaluation of the cells indicated that the Raji cells with Gα~13~ WT expressed were considerably smaller and more necrotic than the control Raji cells ([Figures 5I--J](#F5){ref-type="fig"}). Additionally, the Raji cells expressing Gα~13~ WT were significantly less proliferative, based on Ki67 staining and analysis of tumor sections ([Figure 5K](#F5){ref-type="fig"}). These data indicate that restoration and/or overexpression of Gα~13~ WT influences the survival and/or proliferation of the Raji cells within the tumor microenvironment, which may be the mechanism by which it suppresses tumor growth *in vivo*.

Discussion {#S7}
==========

Burkitt's lymphoma and DLBCL are aggressive and prevalent B cell malignancies ([@R28]). Although some genetic alterations have been identified and associated with malignancy of these B cell lymphomas, much is still unknown about the majority of the genetic mutations and how these alterations may affect the development and progression of disease ([@R14], [@R29]). Previous genome wide sequencing studies unveiled frequent mutations in *GNA13*, the gene encoding the G-protein, Gα~13~ ([@R14], [@R15]), and these mutations were found to be highly statistically significant over background cancer mutation rates ([@R3]). Although activity of Gα~13~ has previously been linked to cellular transformation in fibroblasts as well as tumor cell invasion in several cancers including breast cancer and prostate cancer ([@R7], [@R30]--[@R32]), our characterization of the *GNA13* mutations present in human Burkitt's lymphoma and DLBCL tumors indicates that these mutations rather result in loss of function and/or reduced protein stability. Furthermore, we demonstrated that expression of wild-type Gα~13~ to Burkitt's lymphoma cells with mutated *GNA13* results in reduced tumor growth in xenograft models and altered cell morphology and growth in a methocellulose colony assay.

Recent elegant genetic mouse models support our loss-of-function characterization of *GNA13* mutations present in DLBCL and Burkitt's lymphoma, as mice deficient in lymphoid lineage for Gα~13~ or the Gα~13~-coupled receptor, sphingosine 1 phosphate receptor 2 (S1P2), exhibited an increase in total T cell and B cell numbers and resulted in an outgrowth of cells classified as diffuse large B cell lymphoma ([@R16], [@R17], [@R33]). Specifically, these studies identified a migration defect and an increase in phosphorylated Akt in germinal center (GC) B cells caused by S1P2 or Gα~12~/Gα~13~ deficiency. Thus, signaling through the Gα~12~/Gα~13~ axis via GPCRs, including S1P2 and the recently characterized P2Ry8, is proposed to help confine B cells to the GC and limit B cell expansion ([@R17]). Loss of this axis would therefore allow cells to escape the GC and populate lymph nodes as is characteristic of diffuse large B cell lymphomas. However, in addition to the migration defects observed in these previous studies, our data indicate a direct role for Gα~13~ in growth suppression of B cells *in vivo* and suggest Gα~13~ may influence their survival, and/or differentiation state.

Previous sequencing of non-Hodgkin's lymphomas reported frequent *GNA13* mutations that were enriched in germinal center B cell (GCB) type DLBCL ([@R34]). Of interest, DLBCL have been recently classified in multiple distinct groups based on their histological and biological characteristics and morphological features ([@R35]). Thus, based on our findings it is possible that DLBCL harboring *RHOA* or *GNA13* mutations may exhibit particular characteristics. Indeed, further characterization of *RHOA* and *GNA13* mutations in relation to subgroup, stage and patient survival is warranted once more information becomes available.

To evaluate the mechanisms by which the *GNA13* mutations may alter signaling that could affect the growth and morphology of the B cells, we first investigated possible effects on the major downstream effector, RhoA. The constitutive/basal activity of RhoA was found to be significantly reduced in all 9 of the Gα~13~ mutants tested compared with WT Gα~13~. Although a few of the mutants still had some ligand-induced signaling responses, the overall levels of activity were reduced compared with WT. Furthermore, these ligand-induced signaling assays were performed with receptor overexpression systems, so it is possible that these ligand-responses may not be as pronounced in more physiological systems.

Because the Gα~13~ mutants critically impacted RhoA activity, we investigated whether independent mutations in *RHOA* were present in DLBCL and Burkitt's lymphoma tumors by searching the COSMIC database. Indeed, mutations in *RHOA* were identified in several patient tumors that were independent from those with *GNA13* mutations. Intriguingly, a recurrent glutamine substitution at R5 of RhoA was observed in several patient tumors, aligned with the recent description of *RHOA* mutations, including a recurrent R5Q mutation, in 8.5% of pediatric Burkitt's Lymphoma cases studied ([@R36]), Correlating with the loss of function data observed with the Gα~13~ mutations, our data indicate that the activity of R5Q is significantly impaired compared with WT RhoA in active RhoA pull down, SRF-RE luciferase and actin stress fiber formation assays. Based on structural modeling and known interactions, R5 appears to be directly involved in the interaction with RhoGEFs (particularly RhoGEFs 11, 12, and DBS). Although no mutations were found in the RhoGEFs most frequently activated by Gα13, p115 (ARHGEF1), LARG (ARHGEF12), or PDZ RhoGEF (ARHGEF11), this is not surprising due to the redundancy of the RhoGEFs and the opportunities for compensation ([@R37], [@R38]). Interestingly, inactivating or dominant negative mutations in RhoA were recently identified in 50--70% of peripheral T cell lymphomas ([@R18]--[@R22]), suggesting that the Gα~13~/RhoA pathway may have tumor suppressive functions in other haematological malignancies as well.

Although Gα~13~ and RhoA have been associated with cellular transformation and characterized as growth promoting in fibroblasts and some epithelial cancer models, this axis appears to have the opposite effect in B cell lymphomas. This sort of duality in function based on cellular context is not unheard of and has been observed with other key pathways, including JNK and transforming growth factor beta signaling pathways ([@R39], [@R40]). Although the precise mechanism by which Gα~13~/RhoA axis may inhibit cell growth in B cells requires further exploration, the reduction in phosphorylated Akt (Ser473) associated with Gα~13~ activity could in part affect the growth and survival of B cells ([@R16], [@R17]). Additionally, we observed distinct morphological changes in the cells with mutant Gα~13~ compared with WT-Gα~13~ both in the methocelluolose clonogenic assay and in the histology analysis from tumor xenografts, suggesting other factors may play a role. Nevertheless, the reasons why Gα~13~ and RhoA may exhibit pro-tumorigenic effects in some cancers while having tumor suppressive functions in other cancers requires further investigation. In this regard, once better tools become available for manipulation of B cells and B cell lines, it will be possible to determine how disruption of *GNA13* or *RHOA* in normal human B cells and other Burkitt's lymphoma and DLBCL lines influences their growth and morphology, thereby helping to elucidate fully the repertoire of tumor suppressive signals induced by Gα~13~/RhoA in B cells.

Characterization of tumor mutations provides valuable information for how to precisely identify causes of disease and best treatment options. Our data suggest that the mutations identified in the Gα~13~ signaling pathway in Burkitt's lymphoma and DLBCL tumors result in loss of function, and that restoration and/or activation of this signaling pathway may help reduce tumor growth and progression. Overall, the characterization of these mutations may enable more precise treatment in the \>10% of patients with mutations in *GNA13* and *RHOA* in DLBCL and Burkitt's lymphoma.

Materials and Methods {#S8}
=====================

Identification of interaction interfaces for mapping and structural characterization of mutations {#S9}
-------------------------------------------------------------------------------------------------

Structures of RhoA protein and its complexes were retrieved from the PDB. Per-residue contact strengths of RhoA with each of the available co-crystallized binding partners were measured as described in ([@R41], [@R42]) for residue side-chains only and averaged over multiple structures of homologous complexes. Based on this analysis, residues were classified as a part of GAP interface, GEF interface, or nucleotide binding pocket. Conformational switches I and II were found to be extensively involved in binding of both GAPs and GEFs. The analysis for Gα~13~ was conducted similarly; however, because there are no structures of trimeric α/β/γ complexes for Gα~13~ and because the N-terminal helix is absent from all available structures, this part of the protein was modeled and Gβ/γ interface was mapped by homology with bovine Gαt (transducin, PDB 1got), rat Gαi1 (PDB 1gp2), and mouse Gαq (PDB 3ah8).

DNA constructs and site directed mutagenesis {#S10}
--------------------------------------------

MSCV puro IRES GFP (PIG) and pBABE puro retroviral vectors were obtained from Addgene ([@R43], [@R44]). The SyR-Gi plasmid (also known as Gi-DREADD) was originally provided by Dr. Bryan Roth and was subcloned into the pBABE retroviral vector. Gα~13~ and N-terminally tagged AU5-RhoA were cloned into our in house pCEFL vector and mutations as indicated were prepared by site directed mutagenesis using the QuikChange Lightning Kit (Agilent, Santa Clara, CA). Gα~13-i5~ chimeras were prepared by PCR substitution of the last 5 amino acids of Gα~13~ with the Gαi1 sequence including a C-I mutation to make the plasmid pertussis toxin insensitive (the last 5 amino acids were as follows: DIGLF). Sequences were confirmed by DNA sequence analysis (NIDCR shared resource facility). A plasmid for alkaline phosphatase-fused transforming growth factor-α (AP-TGFα) ([@R45]) was kindly provided by Prof. Shigeki Higashiyama, Ehime University, Japan. A human dopamine D2 receptor (D2R, long isoform) was inserted into the pCAGGS expression plasmid as previously described ([@R27]).

Genomic DNA isolation and sequencing {#S11}
------------------------------------

Raji cells were obtained from ATCC (ATCC CCL-86). Genomic DNA was isolated from 5×10^6^ Raji cells using the QiaAmp DNA isolation kit (Qiagen) and PCR was performed on the genomic DNA in the region surrounding the 184 amino acid position using AccuPrime Pfx Supermix (Life Technologies, Grand Island, NY). PCR primer details are included in [supplemental information](#SD1){ref-type="supplementary-material"}.

Retrovirus production and infection {#S12}
-----------------------------------

Retrovirus from MSCV PIG empty vector and expressing Gα~13~ WT was produced by transfection of HEK293t/17 cells in a poly-lysine coated 6-well plate with 1 ug of MSCV plasmid DNA, 0.67 μg of gag/pol and 0.33 μg of VSV-G using Turbofect transfection reagent (Life Technologies). Virus for pBABE SyR-Gi was similarly prepared by transfecting 1 μg of pBABE SyR-Gi, 0.67 μg of gag/pol and 0.33 μg of VSV-G. Viral supernatants were collected at 48 h and 72 h post-transfection and filtered with 0.45 μm PVDF-filters onto target cells. For HEK293 cells transduced with pBABE SyR-Gi, 6 μg/ml of polybrene was added with the viral supernatants; cells were then selected with 1 μg/ml puromycin (Life Technologies) to generate a stable line. For transduction of Raji cells with MSCV virus, Raji cells were seeded onto retronectin-coated plates (Takara, Madison, WI) and centrifuged with viral supernatants at 2000 rpm for 2 hr and then incubated overnight at 37°C/5% CO~2~. MSCV-transduced Raji cells were allowed to recover before FACS sorting for GFP expression using a BD FACSAria III Cell Sorter (NIDCR Shared Resource Facility). Raji cells were cultured in RPMI media (Sigma, St. Louis, MO) supplemented with 10% FBS.

SRF-RE luciferase assay {#S13}
-----------------------

A serum response factor response element (SRF-RE) luciferase assay was performed by seeding HEK293 cells in a poly-lysine coated 24-well plate and culturing in DMEM supplemented with 10% FBS for 24 h. Cells were then co-transfected with 100 ng Gα~13~ plasmid DNA or control, 50 ng pSRF-RE-firefly luciferase reporter DNA, and 20 ng pRL-renilla luciferase using lipofectamine 2000 (Life Technologies) transfection reagent. For stimulated SRF response, HEK293 cells stably transduced with SyR-Gi were similarly transfected. The day after transfection, cells were serum starved overnight in the presence of 50 ng/ml pertussis toxin (List Biological Laboratories, Campbell, CA) and then stimulated for 6 h with CNO prior to harvesting the cells. For all SRF-luciferase assays, cells were lysed and luciferase activity was determined using Dual-Glo Luciferase Assay Kit (Promega, Madison, WI). Chemiluminescence was measured using a BioTek Synergy Neo plate reader and the SRF-RE activation was calculated as the ratio of firefly to renilla luciferase levels. The assays were performed three times in duplicates and results are presented as the mean +/− standard deviation.

Rho-GTP pull down assay and western blot {#S14}
----------------------------------------

Antibodies to RhoA (catalog \# 2117), pAkt (Ser473) (catalog \#4060), Akt (catalog \#9272), and α-tubulin (catalog \#3873) were obtained from Cell Signaling Technology. Gα~13~ antibody was obtained from Santa Cruz (catalog \# sc-26788). Cells for western blot lysates were lysed on ice in RIPA buffer (Sigma) containing protease and phosphatase inhibitors (Thermo Scientific, Pittsburgh PA) and clarified by centrifugation at 13,000 rpm for 10 min at 4°C. The relative Gα~13~ protein expression levels were quantified by densitometry using ImageJ and normalized to α-tubulin; data represents an average of three independent western blots (mean +/− standard deviation). GST-tagged Rhotekin RBD beads were purchased from Millipore for pull down assays to detect active RhoA. HEK293 cells were seeded in 6 cm plates and transfected with lipofectamine 2000 for 48 h with 4 μg of plasmid DNAs. Cells were serum starved overnight and harvested for pull downs with 350 μl of recommended lysis buffer. Lysates were clarified and 50 μl of lysate was saved for input and the remaining lysate was rotated in tubes with the GST-tagged Rhotekin RBD beads. A three minute thrombin (Sigma) stimulation (1 U/mL) was used as a positive control. Lysates were resolved on SDS-PAGE gels, transferred onto PVDF membranes (Millipore, Billerica, MA), probed with appropriate antibodies and developed with ECL and film (GE healthcare, Piscataway, NJ) or scanned on an Odyssey imager.

Tumor xenograft model and Tissue Analysis {#S15}
-----------------------------------------

Ten-week old female NOD-*scid* IL2Rgamma^null^ (NOD.Cg-*Prkdc^scid^ Il2rg^tm1Wjl^*/SzJ) mice were obtained from the Jackson Laboratory and randomly assigned for injection with tumor cells. Raji cells transduced with MSCV PIG control or MSCV PIG Gα~13~ WT were resuspended in phosphate buffered saline (PBS) and 5×10^6^ cells in a 200 μl volume were injected into the right and left flanks of five mice for each group (ten tumors in total for each group. One injection was excluded from the Raji control group due to an injection error). Tumors were monitored two to three times a week and mice were sacrificed all together once the largest tumors reached approximately 1 cm^3^ volume. Tumor measurements were performed by a blinded individual and the measurements were recorded by an un-blinded person. An initial pilot study was performed with 4 tumors in each group, which was used to determine the appropriate sample size. This study was then expanded and repeated two more times; data trends were the same in all three experiments. Necropsies were performed and tumors were weighed, measured and imaged and tumor sections were cut and fixed in z-fix and embedded in paraffin for histology and immunohistochemistry (IHC) analysis as previously described ([@R46]) and elaborated in [supplemental information](#SD1){ref-type="supplementary-material"}. Individual data points are shown in the scatter dot plots to indicate variance in the data and indicates the mean +/− standard error of the mean. These animal studies were carried out according to National Institutes of Health (NIH) approved protocols (ASP \# 13-695) in compliance with the NIH Guide for the Care and Use of Laboratory Mice.

TGFα shedding assay {#S16}
-------------------

TGFα shedding assay was performed as described previously ([@R27]) and as detailed in [supplemental information](#SD1){ref-type="supplementary-material"}. Briefly, HEK293A cells were seeded in 12-well plates and transfected using lipofectamine 2000 (Life Technologies) with a combination of the following plasmids: AP-TGFα-encoding plasmid, a GPCR-encoding plasmid and a Gα-encoding plasmid. After 24 h, cells were harvested with a trypsin-EDTA solution, washed with D-PBS, suspended in Hank's balanced salt solution containing 5 mM HEPES (pH 7.4), and seeded in a 96-well plate and incubated for 30 min. The seeded cells were treated with various concentrations of a GPCR ligand (10 μl per well) and incubated for 1 h. Conditioned media (CM, 80 μl/well) was transferred into a blank 96-well plate. Alkaline phosphatase (AP) solution was added to both the CM plate and the cell plate (80 μl/well). Absorbance at 405 nm of the two plates was measured before and after 1 h incubation at room temperature, using a SpectraMax 340PC384 microplate reader (Molecular Devices, Sunnyvale, CA). TGFα shedding data are shown as the mean +/− standard deviation.

Spontaneous activity of Gα~13~ was measured by a previously described method with a slight modification ([@R47]). Briefly, HEK293A cells were seeded in 96-well plates in Opti-MEM (80 μl/well), transfected with lipofectamine 2000 reagent (0.1 μl/well), the AP-TGFα-encoding plasmids (20 ng/well), Gα~13~-encoding plasmid (5 to 20 ng/well) and the pCEFL vector plasmid, which was used to adjust total volume of transfected plasmid. After adding the transfection solution, the cells were incubated for 24 h and then AP-TGFα release was measured as described above.

Actin stress fiber formation and confocal imaging {#S17}
-------------------------------------------------

Madin-Darby Canine Kidney (MDCK) cells were seeded onto poly-lysine coated glass cover slips (Fisher) in a 6-well plate and transfected with 1 ug LifeActin GFP and 1 ug of pCEFL RhoA plasmid or vector control using lipofectamine 2000. Cells for immunofluorescence (IF) imaging were washed with PBS and fixed with 2% formaldehyde -- PBS solution for 12 min at room temperature. Fixed cells were washed twice with PBS and then mounted onto glass slides (Thermo Scientific/Fisher) with Fluorosave (Millipore/Calbiochem) mounting solution. Confocal images were collected on a Zeiss LSM-700 laser scanning microscope with a 40× oil immersion lens using 488 nm excitation for the LifeActin GFP.

Methocellulose clonogenic assay {#S18}
-------------------------------

Methocult^™^ (H4434) was purchased from Stem Cell Technologies and a clonogenic assay was performed by seeding 2000 Raji cells in 1 mL of Methocult^™^ in triplicate in a 6 well plate. Colonies were allowed to grow for 14 days and then number of colonies per well and morphology of the colonies were analyzed. Results are shown as the mean +/− standard deviation. Statistics were calculated using GraphPad Prism 6 using a paired t-test.

Supplementary Material {#S19}
======================

This study was supported by the National Institute of Dental and Craniofacial Research intramural program at NIH. We thank Maria S. Degese for her help with immunohistochemistry. A.I. was funded by PRESTO from JST. J.A. was funded by CREST from JST. We thank Miho Morikawa for technical assistance with the TGFα shedding assay. RAD and MSL are supported by the Division of Intramural Research, NIAID, NIH. IK is supported by NIH grants R01 GM071872 and R01 AI118985.

**Conflict of Interest**

Authors have no conflicts of interest to declare.

![Sequence and structure localization of mutations in *GNA13* and *RHOA* that are observed in Burkitt's Lymphoma and DLBCL tumors. A) Table of the number and percentage of *GNA13* mutations in haematopoietic malignancies overall and in Burkitt's Lymphoma and DLBCL based on data from COSMIC v72 (top). Linear diagram of mutations along the *GNA13* gene and the functional and structural domains of the protein Gα~13~ (bottom). B) Table of the number and percentage of *RHOA* mutations in haematopoietic malignancies overall and in Burkitt's Lymphoma and DLBCL based on data from COSMIC v72 (top). Linear diagram of mutations along the *RHOA* gene and the functional and structural domains of the protein and the ribbon diagram of RhoA crystal structure (bottom). Structural and functional domains of both proteins are color-coded and indicated on separate lines above the linear diagram and mapped onto the 3D structures. The mutated residue positions are shown as stars on over the corresponding residues of the linear diagrams and as spheres in the 3D structure representation. The nucleotide is shown as yellow skin.](nihms731271f1){#F1}

![*GNA13* mutations in Burkitt's lymphoma and DLBCL tumor samples result in loss of downstream transcriptional activity based on SRF-RE luciferase assay. A) Representative western blot for Gα~13~ expression of cells transfected with Gα~13~ WT, mutants based on COSMIC data, constitutively active Q226L mutant, or vector control. Antibody detection of α-tubulin was used as a loading control (top). Densitometry quantification of Gα~13~ protein expression of the mutants relative to WT (bottom) as an average of three independent western blots. Error bars indicate standard deviation. B) SRF-RE luciferase assay of spontaneous activity of cells transfected with Gα~13~ WT, mutants, constitutively active Q226L mutant, or vector control. Error bars indicate standard deviation. C) SRF-RE luciferase activity of cells expressing SyR-Gi transfected with WT or mutant Gα~13-i5~ chimeras, in the presence (+) and absence (−) of CNO stimulation. Error bars indicate standard deviation.](nihms731271f2){#F2}

![*GNA13* mutations in Burkitt's lymphoma and DLBCL tumor samples result in loss of function based on a TGFα shedding assay. A) Schematic representation of the AP-TGFα shedding assay used to monitor activity of the Gα~13~ mutants. B) AP-TGFα shedding assay for spontaneous/basal activity of Gα~13~ WT and mutants. Error indicates standard deviation. C) AP-TGFα shedding assay dose response curves of CNO stimulation of cells expressing SyR-Gi transfected with WT or mutant Gα~13-i5~ chimeras.](nihms731271f3){#F3}

![Burkitt's lymphoma and DLBCL *GNA13* and *RHOA* mutations impair RhoA activity. A) GST-Rhotekin-RBD pull down assay to monitor RhoA activity in HEK293 cells transfected with vector control (V), WT, constitutively active Q226L or mutant Gα~13~ plasmids. Thrombin (Th) stimulation was used as a positive control. Input lysates were probed for Gα~13~ and for total RhoA and α-tubulin as loading controls. B) GST-Rhotekin-RBD pull down assay to monitor RhoA activity in HEK293 cells transfected with vector control (V) or RhoA WT, constitutively active Q63L, dominant negative T19N or R5Q mutant plasmids. C) SRF-RE luciferase assay detecting spontaneous activity of cells transfected with vector control, RhoA WT, constitutively active Q63L, dominant negative T19N and the R5Q mutant plasmids. Error bars indicate standard deviation. D) Immunofluorescence images of actin stress fibers in MDCK cells co-transfected with LifeActin-GFP and vector control or RhoA WT, constitutively active Q63L, dominant negative T19N or R5Q mutant plasmids.](nihms731271f4){#F4}

![Gα~13~ WT expression has tumor suppressive effects on Raji Burkitt's Lymphoma xenograft tumor growth. A) Sequence chromatogram of Raji cell DNA validating the presence of a heterozygous T→G mutation resulting in the L184R mutation. B) Immunofluorescence imaging for GFP expression of Raji cells transduced with MSCV puro IRES GFP (PIG) vector control or expressing Gα~13~ WT. C) Western blot of Raji cells transduced with MSCV puro IRES GFP (PIG) vector control or expressing Gα~13~ WT. Lysates were probed for Gα~13~ and pAkt S473 expression and for total Akt and α-tubulin as loading controls. D) Number of colonies formed in a methocellulose clonogenic assay of Raji control and Gα~13~ WT-expressing cells. Data represents the averages of three independent experiments performed in duplicate or triplicate. Error indicates standard deviation and a two-tailed t-test was used to determine statistical significance. E) Percentage of colonies exhibiting spheroid, flat or mixed morphologies (representative images on right) in Raji control and Gα~13~ WT-expressing cells. F) Tumor volume measurements (error bars represent standard error of the mean (SEM)) of Raji control (n=9) and Gα~13~ WT-expressing (G13 WT, n=10) tumor xenografts over time, starting with injections of the tumor cells at Day 0. G) Graphs of the mean (error bars represent SEM) of final tumor masses (left) and final tumor volumes (right) of Raji control and Gα~13~ WT (G13 WT) tumors. H) Representative image of Raji control and Gα~13~ WT tumors. A two-tailed t-test was used to determine statistical significance. I) H&E staining of tissue sections from Raji control and Gα~13~ WT tumors. H&E shows viable round atypical cells infiltrating skeletal muscle, aberrant mitotic figures (arrowhead) an apoptotic cells (star) in the Raji control image. H&E shows the Raji Gα~13~ WT tumor is highly necrotic with numerous ghost cells, eosinophilic cell bodies lacking nuclear structure (star), and chromatin dust (arrowhead) observed. J) Quantification of necrotic area in Raji control (n=5) and Raji Gα~13~ WT (n=5) tumors analyzed. Statistics were based on a Mann Whitney test. K) Quantification of Ki67 immunohistochemistry staining in Raji control (n=5) and Raji Gα~13~ WT (n=5) tumors. Statistics were based on a t-test.](nihms731271f5){#F5}
