Modulation of the proteostasis network promotes tumor resistance to oncogenic KRAS inhibitors.

Despite substantial advances in targeting mutant KRAS, tumor resistance to KRAS inhibitors (KRASi) remains a major barrier to progress. Here, we report proteostasis reprogramming as a key convergence point of multiple KRASi-resistance mechanisms. Inactivation of oncogenic KRAS down-regulated both the heat shock response and the inositol-requiring enzyme 1α (IRE1α) branch of the unfolded protein response, causing severe proteostasis disturbances. However, IRE1α was selectively reactivated in an ER stress-independent manner in acquired KRASi-resistant tumors, restoring proteostasis. Oncogenic KRAS promoted IRE1α protein stability through extracellular signal-regulated kinase (ERK)-dependent phosphorylation of IRE1α, leading to IRE1α disassociation from 3-hydroxy-3-methylglutaryl reductase degradation (HRD1) E3-ligase. In KRASi-resistant tumors, both reactivated ERK and hyperactivated AKT restored IRE1α phosphorylation and stability. Suppression of IRE1α overcame resistance to KRASi. This study reveals a druggable mechanism that leads to proteostasis reprogramming and facilitates KRASi resistance.

Alternative polyadenylation transcriptome-wide association study identifies APA-linked susceptibility genes in brain disorders.

Alternative polyadenylation (APA) plays an essential role in brain development; however, current transcriptome-wide association studies (TWAS) largely overlook APA in nominating susceptibility genes. Here, we performed a 3' untranslated region (3'UTR) APA TWAS (3'aTWAS) for 11 brain disorders by combining their genome-wide association studies data with 17,300 RNA-seq samples across 2,937 individuals. We identified 354 3'aTWAS-significant genes, including known APA-linked risk genes, such as SNCA in Parkinson's disease. Among these 354 genes, ~57% are not significant in traditional expression- and splicing-TWAS studies, since APA may regulate the translation, localization and protein-protein interaction of the target genes independent of mRNA level expression or splicing. Furthermore, we discovered ATXN3 as a 3'aTWAS-significant gene for amyotrophic lateral sclerosis, and its modulation substantially impacted pathological hallmarks of amyotrophic lateral sclerosis in vitro. Together, 3'aTWAS is a powerful strategy to nominate important APA-linked brain disorder susceptibility genes, most of which are largely overlooked by conventional expression and splicing analyses.

STAT5 confers lactogenic properties in breast tumorigenesis and restricts metastatic potential.

Signal transducer and activator of transcription 5 (STAT5) promotes cell survival and instigates breast tumor formation, and in the normal breast it also drives alveolar differentiation and lactogenesis. However, whether STAT5 drives a differentiated phenotype in breast tumorigenesis and therefore impacts cancer spread and metastasis is unclear. We found in two genetically engineered mouse models of breast cancer that constitutively activated Stat5a (Stat5aca) caused precancerous mammary epithelial cells to become lactogenic and evolve into tumors with diminished potential to metastasize. We also showed that STAT5aca reduced the migratory and invasive ability of human breast cancer cell lines in vitro. Furthermore, we demonstrated that STAT5aca overexpression in human breast cancer cells lowered their metastatic burden in xenografted mice. Moreover, RPPA, Western blotting, and studies of ChIPseq data identified several EMT drivers regulated by STAT5. In addition, bioinformatic studies detected a correlation between STAT5 activity and better prognosis of breast cancer patients. Together, we conclude that STAT5 activation during mammary tumorigenesis specifies a tumor phenotype of lactogenic differentiation, suppresses EMT, and diminishes potential for subsequent metastasis.

Exaggerated false positives by popular differential expression methods when analyzing human population samples.

When identifying differentially expressed genes between two conditions using human population RNA-seq samples, we found a phenomenon by permutation analysis: two popular bioinformatics methods, DESeq2 and edgeR, have unexpectedly high false discovery rates. Expanding the analysis to limma-voom, NOISeq, dearseq, and Wilcoxon rank-sum test, we found that FDR control is often failed except for the Wilcoxon rank-sum test. Particularly, the actual FDRs of DESeq2 and edgeR sometimes exceed 20% when the target FDR is 5%. Based on these results, for population-level RNA-seq studies with large sample sizes, we recommend the Wilcoxon rank-sum test.

3'aQTL-atlas: an atlas of 3'UTR alternative polyadenylation quantitative trait loci across human normal tissues.

Genome-wide association studies (GWAS) have identified thousands of non-coding single-nucleotide polymorphisms (SNPs) associated with human traits and diseases. However, functional interpretation of these SNPs remains a significant challenge. Our recent study established the concept of 3' untranslated region (3'UTR) alternative polyadenylation (APA) quantitative trait loci (3'aQTLs), which can be used to interpret ∼16.1% of GWAS SNPs and are distinct from gene expression QTLs and splicing QTLs. Despite the growing interest in 3'aQTLs, there is no comprehensive database for users to search and visualize them across human normal tissues. In the 3'aQTL-atlas (https://wlcb.oit.uci.edu/3aQTLatlas), we provide a comprehensive list of 3'aQTLs containing ∼1.49 million SNPs associated with APA of target genes, based on 15,201 RNA-seq samples across 49 human Genotype-Tissue Expression (GTEx v8) tissues isolated from 838 individuals. The 3'aQTL-atlas provides a ∼2-fold increase in sample size compared with our published study. It also includes 3'aQTL searches by Gene/SNP across tissues, a 3'aQTL genome browser, 3'aQTL boxplots, and GWAS-3'aQTL colocalization event visualization. The 3'aQTL-atlas aims to establish APA as an emerging molecular phenotype to explain a large fraction of GWAS risk SNPs, leading to significant novel insights into the genetic basis of APA and APA-linked susceptibility genes in human traits and diseases.

Notch-induced endoplasmic reticulum-associated degradation governs mouse thymocyte ß-selection.

Signals from the pre-T cell receptor and Notch coordinately instruct ß-selection of CD4-CD8-double negative (DN) thymocytes to generate αß T cells in the thymus. However, how these signals ensure a high-fidelity proteome and safeguard the clonal diversification of the pre-selection TCR repertoire given the considerable translational activity imposed by ß-selection is largely unknown. Here, we identify the endoplasmic reticulum (ER)-associated degradation (ERAD) machinery as a critical proteostasis checkpoint during ß-selection. Expression of the SEL1L-HRD1 complex, the most conserved branch of ERAD, is directly regulated by the transcriptional activity of the Notch intracellular domain. Deletion of Sel1l impaired DN3 to DN4 thymocyte transition and severely impaired mouse αß T cell development. Mechanistically, Sel1l deficiency induced unresolved ER stress that triggered thymocyte apoptosis through the PERK pathway. Accordingly, genetically inactivating PERK rescued T cell development from Sel1l-deficient thymocytes. In contrast, IRE1α/XBP1 pathway was induced as a compensatory adaptation to alleviate Sel1l-deficiency-induced ER stress. Dual loss of Sel1l and Xbp1 markedly exacerbated the thymic defect. Our study reveals a critical developmental signal controlled proteostasis mechanism that enforces T cell development to ensure a healthy adaptive immunity.

An atlas of alternative polyadenylation quantitative trait loci contributing to complex trait and disease heritability.

Genome-wide association studies have identified thousands of noncoding variants associated with human traits and diseases. However, the functional interpretation of these variants is a major challenge. Here, we constructed a multi-tissue atlas of human 3'UTR alternative polyadenylation (APA) quantitative trait loci (3'aQTLs), containing approximately 0.4 million common genetic variants associated with the APA of target genes, identified in 46 tissues isolated from 467 individuals (Genotype-Tissue Expression Project). Mechanistically, 3'aQTLs can alter poly(A) motifs, RNA secondary structure and RNA-binding protein-binding sites, leading to thousands of APA changes. Our CRISPR-based experiments indicate that such 3'aQTLs can alter APA regulation. Furthermore, we demonstrate that mapping 3'aQTLs can identify APA regulators, such as La-related protein 4. Finally, 3'aQTLs are colocalized with approximately 16.1% of trait-associated variants and are largely distinct from other QTLs, such as expression QTLs. Together, our findings show that 3'aQTLs contribute substantially to the molecular mechanisms underlying human complex traits and diseases.

DORGE: Discovery of Oncogenes and tumoR suppressor genes using Genetic and Epigenetic features.

Data-driven discovery of cancer driver genes, including tumor suppressor genes (TSGs) and oncogenes (OGs), is imperative for cancer prevention, diagnosis, and treatment. Although epigenetic alterations are important for tumor initiation and progression, most known driver genes were identified based on genetic alterations alone. Here, we developed an algorithm, DORGE (Discovery of Oncogenes and tumor suppressoR genes using Genetic and Epigenetic features), to identify TSGs and OGs by integrating comprehensive genetic and epigenetic data. DORGE identified histone modifications as strong predictors for TSGs, and it found missense mutations, super enhancers, and methylation differences as strong predictors for OGs. We extensively validated DORGE-predicted cancer driver genes using independent functional genomics data. We also found that DORGE-predicted dual-functional genes (both TSGs and OGs) are enriched at hubs in protein-protein interaction and drug-gene networks. Overall, our study has deepened the understanding of epigenetic mechanisms in tumorigenesis and revealed previously undetected cancer driver genes.

Protein quality control through endoplasmic reticulum-associated degradation maintains haematopoietic stem cell identity and niche interactions.

Stem cells need to be protected from genotoxic and proteotoxic stress to maintain a healthy pool throughout life1-3. Little is known about the proteostasis mechanism that safeguards stem cells. Here we report endoplasmic reticulum-associated degradation (ERAD) as a protein quality checkpoint that controls the haematopoietic stem cell (HSC)-niche interaction and determines the fate of HSCs. The SEL1L-HRD1 complex, the most conserved branch of ERAD4, is highly expressed in HSCs. Deletion of Sel1l led to niche displacement of HSCs and a complete loss of HSC identity, and allowed highly efficient donor-HSC engraftment without irradiation. Mechanistic studies identified MPL, the master regulator of HSC identity5, as a bona fide ERAD substrate that became aggregated in the endoplasmic reticulum following ERAD deficiency. Restoration of MPL signalling with an agonist partially rescued the number and reconstitution capacity of Sel1l-deficient HSCs. Our study defines ERAD as an essential proteostasis mechanism to safeguard a healthy stem cell pool by regulating the stem cell-niche interaction.

Adaptive endoplasmic reticulum stress signalling via IRE1α-XBP1 preserves self-renewal of haematopoietic and pre-leukaemic stem cells.

Over their lifetime, long-term haematopoietic stem cells (HSC) are exposed to a variety of stress conditions that they must endure. Many stresses, such as infection/inflammation, reactive oxygen species, nutritional deprivation and hypoxia, activate unfolded protein response signalling, which induces either adaptive changes to resolve the stress or apoptosis to clear the damaged cell. Whether unfolded-protein-response signalling plays any role in HSC regulation remains to be established. Here, we report that the adaptive signalling of the unfolded protein response, IRE1α-XBP1, protects HSCs from endoplasmic reticulum stress-induced apoptosis. IRE1α knockout leads to reduced reconstitution of HSCs. Furthermore, we show that oncogenic N-RasG12D activates IRE1α-XBP1, through MEK-GSK3ß, to promote HSC survival under endoplasmic reticulum stress. Inhibiting IRE1α-XBP1 abolished N-RasG12D-mediated survival under endoplasmic reticulum stress and diminished the competitive advantage of NrasG12D HSCs in transplant recipients. Our studies illuminate how the adaptive endoplasmic reticulum stress response is advantageous in sustaining self-renewal of HSCs and promoting pre-leukaemic clonal dominance.