The role of HERC2 and RNF8 ubiquitin E3 ligases in the promotion of translesion DNA synthesis in the chicken DT40 cell line

Implications of ubiquitination and the maintenance of replication fork stability in cancer therapy

DNA replication forks are subject to intricate surveillance and strict regulation by sophisticated cellular machinery. Such close regulation is necessary to ensure the accurate duplication of genetic information and to tackle the diverse endogenous and exogenous stresses that impede this process. Stalled replication forks are vulnerable to collapse, which is a major cause of genomic instability and carcinogenesis. Replication stress responses, which are organized via a series of coordinated molecular events, stabilize stalled replication forks and carry out fork reversal and restoration. DNA damage tolerance and repair pathways such as homologous recombination and Fanconi anemia also contribute to replication fork stabilization. The signaling network that mediates the transduction and interplay of these pathways is regulated by a series of post-translational modifications, including ubiquitination, which affects the activity, stability, and interactome of substrates. In particular, the ubiquitination of replication protein A and proliferating cell nuclear antigen at stalled replication forks promotes the recruitment of downstream regulators. In this review, we describe the ubiquitination-mediated signaling cascades that regulate replication fork progression and stabilization. In addition, we discuss the targeting of replication fork stability and ubiquitination system components as a potential therapeutic approach for the treatment of cancer.

NEDD4 E3 ubiquitin ligases: promising biomarkers and therapeutic targets for cancer

Accumulating evidence has demonstrated that NEDD4 E3 ubiquitin ligase family plays a pivotal oncogenic role in a variety of malignancies via mediating ubiquitin dependent degradation processes. Moreover, aberrant expression of NEDD4 E3 ubiquitin ligases is often indicative of cancer progression and correlated with poor prognosis. In this review, we are going to address association of expression of NEDD4 E3 ubiquitin ligases with cancers, the signaling pathways and the molecular mechanisms by which the NEDD4 E3 ubiquitin ligases regulate oncogenesis and progression, and the therapies targeting the NEDD4 E3 ubiquitin ligases. This review provides the systematic and comprehensive summary of the latest research status of E3 ubiquitin ligases in the NEDD4 subfamily, and proposes that NEDD4 family E3 ubiquitin ligases are promising anti-cancer drug targets, aiming to provide research direction for clinical targeting of NEDD4 E3 ubiquitin ligase therapy.

DNA replication: Mechanisms and therapeutic interventions for diseases

Accurate and integral cellular DNA replication is modulated by multiple replication-associated proteins, which is fundamental to preserve genome stability. Furthermore, replication proteins cooperate with multiple DNA damage factors to deal with replication stress through mechanisms beyond their role in replication. Cancer cells with chronic replication stress exhibit aberrant DNA replication and DNA damage response, providing an exploitable therapeutic target in tumors. Numerous evidence has indicated that posttranslational modifications (PTMs) of replication proteins present distinct functions in DNA replication and respond to replication stress. In addition, abundant replication proteins are involved in tumorigenesis and development, which act as diagnostic and prognostic biomarkers in some tumors, implying these proteins act as therapeutic targets in clinical. Replication-target cancer therapy emerges as the times require. In this context, we outline the current investigation of the DNA replication mechanism, and simultaneously enumerate the aberrant expression of replication proteins as hallmark for various diseases, revealing their therapeutic potential for target therapy. Meanwhile, we also discuss current observations that the novel PTM of replication proteins in response to replication stress, which seems to be a promising strategy to eliminate diseases.

The regulatory roles of the E3 ubiquitin ligase NEDD4 family in DNA damage response

E3 ubiquitin ligases, an important part of ubiquitin proteasome system, catalyze the covalent binding of ubiquitin to target substrates, which plays a role in protein ubiquitination and regulates different biological process. DNA damage response (DDR) is induced in response to DNA damage to maintain genome integrity and stability, and this process has crucial significance to a series of cell activities such as differentiation, apoptosis, cell cycle. The NEDD4 family, belonging to HECT E3 ubiquitin ligases, is reported as regulators that participate in the DDR process by recognizing different substrates. In this review, we summarize recent researches on NEDD4 family members in the DDR and discuss the roles of NEDD4 family members in the cascade reactions induced by DNA damage. This review may contribute to the further study of pathophysiology for certain diseases and pharmacology for targeted drugs.

Ube2N is present and functions within listeria Actin-rich structures and lamellipodia: A localization and pharmacological inhibition study

The actin cytoskeleton forms much of the structure needed for the intracellular motility of an assortment of microbes as well as entire cells. The co-factor to the ubiquitin conjugating enzyme Ube2N (Ube2V1) has been implicated in both cancer cell metastasis and lysine-63 ubiquitylation of β actin. As this protein complexes with Ube2N, we sought to investigate whether Ube2N itself was involved in actin-based events occurring during the Listeria monocytogenes infections as well as within motile whole cells. Through examination of L. monocytogenes actin clouds, comet tails and membrane protrusions as well as lamellipodia in migrating cells, we show that Ube2N is recruited to actin-rich structures. When pharmacologically inhibited we demonstrate that Ube2N is crucial for the function of actin-rich structures when associated with the plasma membrane.

A large Canadian cohort provides insights into the genetic architecture of human hair colour

Hair colour is a polygenic phenotype that results from differences in the amount and ratio of melanins located in the hair bulb. Genome-wide association studies (GWAS) have identified many loci involved in the pigmentation pathway affecting hair colour. However, most of the associated loci overlap non-protein coding regions and many of the molecular mechanisms underlying pigmentation variation are still not understood. Here, we conduct GWAS meta-analyses of hair colour in a Canadian cohort of 12,741 individuals of European ancestry. By performing fine-mapping analyses we identify candidate causal variants in pigmentation loci associated with blonde, red and brown hair colour. Additionally, we observe colocalization of several GWAS hits with expression and methylation quantitative trait loci (QTLs) of cultured melanocytes. Finally, transcriptome-wide association studies (TWAS) further nominate the expression of EDNRB and CDK10 as significantly associated with hair colour. Our results provide insights on the mechanisms regulating pigmentation biology in humans.

Post-Translational Modifications of PCNA in Control of DNA Synthesis and DNA Damage Tolerance-the Implications in Carcinogenesis

The faithful DNA replication is a critical event for cell survival and inheritance. However, exogenous or endogenous sources of damage challenge the accurate synthesis of DNA, which causes DNA lesions. The DNA lesions are obstacles for replication fork progression. However, the prolonged replication fork stalling leads to replication fork collapse, which may cause DNA double-strand breaks (DSB). In order to maintain genomic stability, eukaryotic cells evolve translesion synthesis (TLS) and template switching (TS) to resolve the replication stalling. Proliferating cell nuclear antigen (PCNA) trimer acts as a slide clamp and encircles DNA to orchestrate DNA synthesis and DNA damage tolerance (DDT). The post-translational modifications (PTMs) of PCNA regulate these functions to ensure the appropriate initiation and termination of replication and DDT. The aberrant regulation of PCNA PTMs will result in DSB, which causes mutagenesis and poor response to chemotherapy. Here, we review the roles of the PCNA PTMs in DNA duplication and DDT. We propose that clarifying the regulation of PCNA PTMs may provide insights into understanding the development of cancers.

Genetic evidence for the involvement of mismatch repair proteins, PMS2 and MLH3, in a late step of homologous recombination

Homologous recombination (HR) repairs DNA double-strand breaks using intact homologous sequences as template DNA. Broken DNA and intact homologous sequences form joint molecules (JMs), including Holliday junctions (HJs), as HR intermediates. HJs are resolved to form crossover and noncrossover products. A mismatch repair factor, MLH3 endonuclease, produces the majority of crossovers during meiotic HR, but it remains elusive whether mismatch repair factors promote HR in nonmeiotic cells. We disrupted genes encoding the MLH3 and PMS2 endonucleases in the human B cell line, TK6, generating null MLH3^-/- and PMS2^-/- mutant cells. We also inserted point mutations into the endonuclease motif of MLH3 and PMS2 genes, generating endonuclease death MLH3^DN/DN and PMS2^EK/EK cells. MLH3^-/- and MLH3^DN/DN cells showed a very similar phenotype, a 2.5-fold decrease in the frequency of heteroallelic HR-dependent repair of restriction enzyme-induced double-strand breaks. PMS2^-/- and PMS2^EK/EK cells showed a phenotype very similar to that of the MLH3 mutants. These data indicate that MLH3 and PMS2 promote HR as an endonuclease. The MLH3^DN/DN and PMS2^EK/EK mutations had an additive effect on the heteroallelic HR. MLH3^DN/DN/PMS2^EK/EK cells showed normal kinetics of γ-irradiation-induced Rad51 foci but a significant delay in the resolution of Rad51 foci and a 3-fold decrease in the number of cisplatin-induced sister chromatid exchanges. The ectopic expression of the Gen1 HJ re-solvase partially reversed the defective heteroallelic HR of MLH3^DN/DN/PMS2^EK/EK cells. Taken together, we propose that MLH3 and PMS2 promote HR as endonucleases, most likely by processing JMs in mammalian somatic cells.

HECT E3 ubiquitin ligases - emerging insights into their biological roles and disease relevance

Homologous to E6AP C-terminus (HECT) E3 ubiquitin ligases play a critical role in various cellular pathways, including but not limited to protein trafficking, subcellular localization, innate immune response, viral infections, DNA damage responses and apoptosis. To date, 28 HECT E3 ubiquitin ligases have been identified in humans, and recent studies have begun to reveal how these enzymes control various cellular pathways by catalyzing the post-translational attachment of ubiquitin to their respective substrates. New studies have identified substrates and/or interactors with different members of the HECT E3 ubiquitin ligase family, particularly for E6AP and members of the neuronal precursor cell-expressed developmentally downregulated 4 (NEDD4) family. However, there still remains many unanswered questions about the specific roles that each of the HECT E3 ubiquitin ligases have in maintaining cellular homeostasis. The present Review discusses our current understanding on the biological roles of the HECT E3 ubiquitin ligases in the cell and how they contribute to disease development. Expanded investigations on the molecular basis for how and why the HECT E3 ubiquitin ligases recognize and regulate their intracellular substrates will help to clarify the biochemical mechanisms employed by these important enzymes in ubiquitin biology.

Spatiotemporal regulation of PCNA ubiquitination in damage tolerance pathways

DNA is constantly exposed to a wide variety of exogenous and endogenous agents, and most DNA lesions inhibit DNA synthesis. To cope with such problems during replication, cells have molecular mechanisms to resume DNA synthesis in the presence of DNA lesions, which are known as DNA damage tolerance (DDT) pathways. The concept of ubiquitination-mediated regulation of DDT pathways in eukaryotes was established via genetic studies in the yeast Saccharomyces cerevisiae, in which two branches of the DDT pathway are regulated via ubiquitination of proliferating cell nuclear antigen (PCNA): translesion DNA synthesis (TLS) and homology-dependent repair (HDR), which are stimulated by mono- and polyubiquitination of PCNA, respectively. Over the subsequent nearly two decades, significant progress has been made in understanding the mechanisms that regulate DDT pathways in other eukaryotes. Importantly, TLS is intrinsically error-prone because of the miscoding nature of most damaged nucleotides and inaccurate replication of undamaged templates by TLS polymerases (pols), whereas HDR is theoretically error-free because the DNA synthesis is thought to be predominantly performed by pol δ, an accurate replicative DNA pol, using the undamaged sister chromatid as its template. Thus, the regulation of the choice between the TLS and HDR pathways is critical to determine the appropriate biological outcomes caused by DNA damage. In this review, we summarize our current understanding of the species-specific regulatory mechanisms of PCNA ubiquitination and how cells choose between TLS and HDR. We then provide a hypothetical model for the spatiotemporal regulation of DDT pathways in human cells.

CtIP-BRCA1 complex and MRE11 maintain replication forks in the presence of chain terminating nucleoside analogs

Chain-terminating nucleoside analogs (CTNAs), which cannot be extended by DNA polymerases, are widely used as antivirals or anti-cancer agents, and can induce cell death. Processing of blocked DNA ends, like camptothecin-induced trapped-topoisomerase I, can be mediated by TDP1, BRCA1, CtIP and MRE11. Here, we investigated whether the CtIP-BRCA1 complex and MRE11 also contribute to cellular tolerance to CTNAs, including 2',3'-dideoxycytidine (ddC), cytarabine (ara-C) and zidovudine (Azidothymidine, AZT). We show that BRCA1-/-, CtIPS332A/-/- and nuclease-dead MRE11D20A/- mutants display increased sensitivity to CTNAs, accumulate more DNA damage (chromosomal breaks, γ-H2AX and neutral comets) when treated with CTNAs and exhibit significant delays in replication fork progression during exposure to CTNAs. Moreover, BRCA1-/-, CtIPS332A/-/- and nuclease-dead MRE11D20A/- mutants failed to resume DNA replication in response to CTNAs, whereas control and CtIP+/-/- cells experienced extensive recovery of DNA replication. In summary, we provide clear evidence that MRE11 and the collaborative action of BRCA1 and CtIP play a critical role in the nuclease-dependent removal of incorporated ddC from replicating genomic DNA. We propose that BRCA1-CTIP and MRE11 prepare nascent DNA ends, blocked from synthesis by CTNAs, for further repair.

HERCing: Structural and Functional Relevance of the Large HERC Ubiquitin Ligases

Homologous to the E6AP carboxyl terminus (HECT) and regulator of chromosome condensation 1 (RCC1)-like domain-containing proteins (HERCs) belong to the superfamily of ubiquitin ligases. HERC proteins are divided into two subfamilies, Large and Small HERCs. Despite their similarities in terms of both structure and domains, these subfamilies are evolutionarily very distant and result from a convergence phenomenon rather than from a common origin. Large HERC genes, HERC1 and HERC2, are present in most metazoan taxa. They encode very large proteins (approximately 5,000 amino acid residues in a single polypeptide chain) that contain more than one RCC1-like domain as a structural characteristic. Accumulating evidences show that these unusually large proteins play key roles in a wide range of cellular functions which include neurodevelopment, DNA damage repair, and cell proliferation. To better understand the origin, evolution, and function of the Large HERC family, this minireview provides with an integrated overview of their structure and function and details their physiological implications. This study also highlights and discusses how dysregulation of these proteins is associated with severe human diseases such as neurological disorders and cancer.

PDIP38/PolDIP2 controls the DNA damage tolerance pathways by increasing the relative usage of translesion DNA synthesis over template switching

Replicative DNA polymerases are frequently stalled at damaged template strands. Stalled replication forks are restored by the DNA damage tolerance (DDT) pathways, error-prone translesion DNA synthesis (TLS) to cope with excessive DNA damage, and error-free template switching (TS) by homologous DNA recombination. PDIP38 (Pol-delta interacting protein of 38 kDa), also called Pol δ-interacting protein 2 (PolDIP2), physically associates with TLS DNA polymerases, polymerase η (Polη), Polλ, and PrimPol, and activates them in vitro. It remains unclear whether PDIP38 promotes TLS in vivo, since no method allows for measuring individual TLS events in mammalian cells. We disrupted the PDIP38 gene, generating PDIP38^-/- cells from the chicken DT40 and human TK6 B cell lines. These PDIP38^-/- cells did not show a significant sensitivity to either UV or H₂O₂, a phenotype not seen in any TLS-polymerase-deficient DT40 or TK6 mutants. DT40 provides a unique opportunity of examining individual TLS and TS events by the nucleotide sequence analysis of the immunoglobulin variable (Ig V) gene as the cells continuously diversify Ig V by TLS (non-templated Ig V hypermutation) and TS (Ig gene conversion) during in vitro culture. PDIP38^-/- cells showed a shift in Ig V diversification from TLS to TS. We measured the relative usage of TLS and TS in TK6 cells at a chemically synthesized UV damage (CPD) integrated into genomic DNA. The loss of PDIP38 also caused an increase in the relative usage of TS. The number of UV-induced sister chromatid exchanges, TS events associated with crossover, was increased a few times in PDIP38^-/- human and chicken cells. Collectively, the loss of PDIP38 consistently causes a shift in DDT from TLS to TS without enhancing cellular sensitivity to DNA damage. We propose that PDIP38 controls the relative usage of TLS and TS increasing usage of TLS without changing the overall capability of DDT.

DNA Damage Tolerance Mechanisms Revealed from the Analysis of Immunoglobulin V Gene Diversification in Avian DT40 Cells

DNA replication is an essential biochemical reaction in dividing cells that frequently stalls at damaged sites. Homologous/homeologous recombination (HR)-mediated template switch and translesion DNA synthesis (TLS)-mediated bypass processes release arrested DNA replication forks. These mechanisms are pivotal for replication fork maintenance and play critical roles in DNA damage tolerance (DDT) and gap-filling. The avian DT40 B lymphocyte cell line provides an opportunity to examine HR-mediated template switch and TLS triggered by abasic sites by sequencing the constitutively diversifying immunoglobulin light-chain variable gene (IgV). During IgV diversification, activation-induced deaminase (AID) converts dC to dU, which in turn is excised by uracil DNA glycosylase and yields abasic sites within a defined window of around 500 base pairs. These abasic sites can induce gene conversion with a set of homeologous upstream pseudogenes via the HR-mediated template switch, resulting in templated mutagenesis, or can be bypassed directly by TLS, resulting in non-templated somatic hypermutation at dC/dG base pairs. In this review, we discuss recent works unveiling IgV diversification mechanisms in avian DT40 cells, which shed light on DDT mode usage in vertebrate cells and tolerance of abasic sites.

SUMOylation of PCNA by PIAS1 and PIAS4 promotes template switch in the chicken and human B cell lines

DNA damage tolerance (DDT) releases replication blockage caused by damaged nucleotides on template strands employing two alternative pathways, error-prone translesion DNA synthesis (TLS) and error-free template switch (TS). Lys164 of proliferating cell nuclear antigen (PCNA) is SUMOylated during the physiological cell cycle. To explore the role for SUMOylation of PCNA in DDT, we characterized chicken DT40 and human TK6 B cells deficient in the PIAS1 and PIAS4 small ubiquitin-like modifier (SUMO) E3 ligases. DT40 cells have a unique advantage in the phenotypic analysis of DDT as they continuously diversify their immunoglobulin (Ig) variable genes by TLS and TS [Ig gene conversion (GC)], both relieving replication blocks at abasic sites without accompanying by DNA breakage. Remarkably, PIAS1 ^-/- /PIAS4 ^-/- cells displayed a multifold decrease in SUMOylation of PCNA at Lys164 and over a 90% decrease in the rate of TS. Likewise, PIAS1 ^-/- /PIAS4 ^-/- TK6 cells showed a shift of DDT from TS to TLS at a chemosynthetic UV lesion inserted into the genomic DNA. The PCNA ^K164R/K164R mutation caused a ∼90% decrease in the rate of Ig GC and no additional impact on PIAS1 ^-/- /PIAS4 ^-/- cells. This epistatic relationship between the PCNA ^K164R/K164R and the PIAS1 ^-/- /PIAS4 ^-/- mutations suggests that PIAS1 and PIAS4 promote TS mainly through SUMOylation of PCNA at Lys164. This idea is further supported by the data that overexpression of a PCNA-SUMO1 chimeric protein restores defects in TS in PIAS1 ^-/- /PIAS4 ^-/- cells. In conclusion, SUMOylation of PCNA at Lys164 promoted by PIAS1 and PIAS4 ensures the error-free release of replication blockage during physiological DNA replication in metazoan cells.

Network Analysis of UBE3A/E6AP-Associated Proteins Provides Connections to Several Distinct Cellular Processes

Perturbations in activity and dosage of the UBE3A ubiquitin-ligase have been linked to Angelman syndrome and autism spectrum disorders. UBE3A was initially identified as the cellular protein hijacked by the human papillomavirus E6 protein to mediate the ubiquitylation of p53, a function critical to the oncogenic potential of these viruses. Although a number of substrates have been identified, the normal cellular functions and pathways affected by UBE3A are largely unknown. Previously, we showed that UBE3A associates with HERC2, NEURL4, and MAPK6/ERK3 in a high-molecular-weight complex of unknown function that we refer to as the HUN complex (HERC2, UBE3A, and NEURL4). In this study, the combination of two complementary proteomic approaches with a rigorous network analysis revealed cellular functions and pathways in which UBE3A and the HUN complex are involved. In addition to finding new UBE3A-associated proteins, such as MCM6, SUGT1, EIF3C, and ASPP2, network analysis revealed that UBE3A-associated proteins are connected to several fundamental cellular processes including translation, DNA replication, intracellular trafficking, and centrosome regulation. Our analysis suggests that UBE3A could be involved in the control and/or integration of these cellular processes, in some cases as a component of the HUN complex, and also provides evidence for crosstalk between the HUN complex and CAMKII interaction networks. This study contributes to a deeper understanding of the cellular functions of UBE3A and its potential role in pathways that may be affected in Angelman syndrome, UBE3A-associated autism spectrum disorders, and human papillomavirus-associated cancers.

Molecular Basis for K63-Linked Ubiquitination Processes in Double-Strand DNA Break Repair: A Focus on Kinetics and Dynamics

Cells are exposed to thousands of DNA damage events on a daily basis. This damage must be repaired to preserve genetic information and prevent development of disease. The most deleterious damage is a double-strand break (DSB), which is detected and repaired by mechanisms known as non-homologous end-joining (NHEJ) and homologous recombination (HR), which are components of the DNA damage response system. NHEJ is an error-prone first line of defense, whereas HR invokes error-free repair and is the focus of this review. The functions of the protein components of HR-driven DNA repair are regulated by the coordinated action of post-translational modifications including lysine acetylation, phosphorylation, ubiquitination, and SUMOylation. The latter two mechanisms are fundamental for recognition of DSBs and reorganizing chromatin to facilitate repair. We focus on the structures and molecular mechanisms for the protein components underlying synthesis, recognition, and cleavage of K63-linked ubiquitin chains, which are abundant at damage sites and obligatory for DSB repair. The forward flux of the K63-linked ubiquitination cascade is driven by the combined activity of E1 enzyme, the heterodimeric E2 Mms2-Ubc13, and its cognate E3 ligases RNF8 and RNF168, which is balanced through the binding and cleavage of chains by the deubiquitinase BRCC36, and the proteasome, and through the binding of chains by recognition modules on repair proteins such as RAP80. We highlight a number of aspects regarding our current understanding for the role of kinetics and dynamics in determining the function of the enzymes and chain recognition modules that drive K63 ubiquitination.

HERC4 Is Overexpressed in Hepatocellular Carcinoma and Contributes to the Proliferation and Migration of Hepatocellular Carcinoma Cells

Recently, more and more evidences unveiled that ubiquitin-proteasome system (UPS) makes an important contribution to the occurrence and development of cancer. HERC4 is one identified Ubiqutin ligase E3, a member of UPS. Although some studies showed that HERC4 abnormally expresses in many cancer cells, till now, nothing has been reported about the function of HERC4 in the development of hepatoma carcinoma. To this end, in this study, we studied the function of HERC4 for the first time in hepatoma carcinoma cells. We detected the expression of HERC4 in tumor and normal tissues, and in hepatoma carcinoma cell lines by using qRT-PCR, Western blot, immunohistochemistry, and immunofluorescence. The data showed that tumor tissues expressed higher HERC4 than normal ones. HERC4 was expressed, although to a different extent, in hepatoma carcinoma cell lines. Colony formation assay, CCK-8 assay, EdU assay, wound healing assay, and FACS indicated that HERC4 plays a role in cell proliferative and migration ability. HERC4 overexpression increases the proliferative and migration ability and reduces apoptosis of hepatoma carcinoma cells; in contrast, knockdown of HERC4 decreases the proliferative and migration ability and increases the apoptosis rate of hepatoma carcinoma cells. Taken together, our findings showed that HERC4 has an effect on the occurrence and development of hepatoma carcinoma by promoting hepatoma carcinoma cell proliferation and migration, and by reducing cell apoptosis, further providing another therapeutic target for the intervention of related diseases.

Selection of Reference Genes for qRT-PCR Analysis of Gene Expression in Stipa grandis during Environmental Stresses

Stipa grandis P. Smirn. is a dominant plant species in the typical steppe of the Xilingole Plateau of Inner Mongolia. Selection of suitable reference genes for the quantitative real-time reverse transcription polymerase chain reaction (qRT-PCR) is important for gene expression analysis and research into the molecular mechanisms underlying the stress responses of S. grandis. In the present study, 15 candidate reference genes (EF1 beta, ACT, GAPDH, SamDC, CUL4, CAP, SNF2, SKIP1, SKIP5, SKIP11, UBC2, UBC15, UBC17, UCH, and HERC2) were evaluated for their stability as potential reference genes for qRT-PCR under different stresses. Four algorithms were used: GeNorm, NormFinder, BestKeeper, and RefFinder. The results showed that the most stable reference genes were different under different stress conditions: EF1beta and UBC15 during drought and salt stresses; ACT and GAPDH under heat stress; SKIP5 and UBC17 under cold stress; UBC15 and HERC2 under high pH stress; UBC2 and UBC15 under wounding stress; EF1beta and UBC17 under jasmonic acid treatment; UBC15 and CUL4 under abscisic acid treatment; and HERC2 and UBC17 under salicylic acid treatment. EF1beta and HERC2 were the most suitable genes for the global analysis of all samples. Furthermore, six target genes, SgPOD, SgPAL, SgLEA, SgLOX, SgHSP90 and SgPR1, were selected to validate the most and least stable reference genes under different treatments. Our results provide guidelines for reference gene selection for more accurate qRT-PCR quantification and will promote studies of gene expression in S. grandis subjected to environmental stress.