Citation: Yonghua Zhuang, Kristen Wade, Laura M. Saba, Katerina Kechris. Development of a tissue augmented Bayesian model for expression quantitative trait loci analysis[J]. Mathematical Biosciences and Engineering, 2020, 17(1): 122-143. doi: 10.3934/mbe.2020007
[1] | Ignacio Alvarez-Castro, Jarad Niemi . Fully Bayesian analysis of allele-specific RNA-seq data. Mathematical Biosciences and Engineering, 2019, 16(6): 7751-7770. doi: 10.3934/mbe.2019389 |
[2] | Xue Wang, Fang Yang, Hongyuan Liu, Qingxuan Shi . Error correction of semantic mathematical expressions based on bayesian algorithm. Mathematical Biosciences and Engineering, 2022, 19(6): 5428-5445. doi: 10.3934/mbe.2022255 |
[3] | Katrine O. Bangsgaard, Morten Andersen, James G. Heaf, Johnny T. Ottesen . Bayesian parameter estimation for phosphate dynamics during hemodialysis. Mathematical Biosciences and Engineering, 2023, 20(3): 4455-4492. doi: 10.3934/mbe.2023207 |
[4] | Dao-Peng Chen, Fangyuan Zhang, Shili Lin . AIJ: joint test for simultaneous detection of imprinting and non-imprinting allelic expression imbalance. Mathematical Biosciences and Engineering, 2020, 17(1): 366-386. doi: 10.3934/mbe.2020020 |
[5] | Saskya Mary Soemartojo, Titin Siswantining, Yoel Fernando, Devvi Sarwinda, Herley Shaori Al-Ash, Sarah Syarofina, Noval Saputra . Iterative bicluster-based Bayesian principal component analysis and least squares for missing-value imputation in microarray and RNA-sequencing data. Mathematical Biosciences and Engineering, 2022, 19(9): 8741-8759. doi: 10.3934/mbe.2022405 |
[6] | Samuel Bronstein, Stefan Engblom, Robin Marin . Bayesian inference in epidemics: linear noise analysis. Mathematical Biosciences and Engineering, 2023, 20(2): 4128-4152. doi: 10.3934/mbe.2023193 |
[7] | M. Nagy, M. H. Abu-Moussa, Adel Fahad Alrasheedi, A. Rabie . Expected Bayesian estimation for exponential model based on simple step stress with Type-I hybrid censored data. Mathematical Biosciences and Engineering, 2022, 19(10): 9773-9791. doi: 10.3934/mbe.2022455 |
[8] | Elena Izquierdo-Kulich, Margarita Amigó de Quesada, Carlos Manuel Pérez-Amor, José Manuel Nieto-Villar . Morphogenesis and aggressiveness of cervix carcinoma. Mathematical Biosciences and Engineering, 2011, 8(4): 987-997. doi: 10.3934/mbe.2011.8.987 |
[9] | Olaiya Mathilde Adéoti, Aliou Diop, Romain Glèlè Kakaï . Bayesian inference and impact of parameter prior specification in flexible multilevel nonlinear models in the context of infectious disease modeling. Mathematical Biosciences and Engineering, 2025, 22(4): 897-919. doi: 10.3934/mbe.2025032 |
[10] | Jun Xu, Bei Wang, Zhengtao Liu, Mingchun Lai, Mangli Zhang, Shusen Zheng . miR-223-3p regulating the occurrence and development of liver cancer cells by targeting FAT1 gene. Mathematical Biosciences and Engineering, 2020, 17(2): 1534-1547. doi: 10.3934/mbe.2020079 |
PSMD9 (26S Proteasome Non-ATPase Regulatory Subunit 9; also known as bridge-1, p27 or Rpn4) is a 25kDa protein encoded by the PSMD9 gene on chromosome 12q24.31-q24.32. Also known in humans as p27 (not to be confused with p27KIP1), the protein is highly conserved in model organisms such as rodents (in which it is often referred to as bridge-1), yeast (nas2) and nematodes (PSMD9). PSMD9 is ubiquitously expressed in mammalian cell types and appears to be at particularly high levels in lymphatic, endocrine, kidney, and skin epidermis tissues [1].
PSMD9 was first characterised though its interaction with the proteasomal chaperones PSMC6 and PSMC3, with which it forms a complex that modulates the assembly of the 26S proteasome from the constituent 20S catalytic and 19S (PA700) regulatory subunits (Figure 1) [2]. In addition to this proteasomal role, it is also apparent that PSMD9 has a separate role as a transcriptional co-regulator although this role was first described under the alternative name of bridge-1. PSMD9 regulates gene expression of specific target genes by binding of its PDZ domain to the transcription factors PDX-1, E12 and E47 [3,4]. A comprehensive 3D structure of PSMD9 is yet to be elucidated but analyses of the secondary structure have revealed the presence of a unique PDZ domain, a motif associated with specific protein-protein interactions, and it is assumed that this motif underpins both functional roles. The regulatory impact of PSMD9 on both proteasomal activity and transcription remains unclear. The ubiquitin/proteasome pathway is the main non-lysosomal route for intracellular protein degradation and is responsible for the turnover of > 80% of cellular proteins [5]. Despite its role in proteasome assembly, it is unknown whether changes in PSMD9 expression impact on overall proteasomal activity, although it has been reported that PSMD9 mediates the specific proteasomal degradation of IκBα leading to regulation of NF-κB signalling [6]. Similarly, the range of transcriptional targets for the transcriptional co-activator role of PSMD9 is unknown, with experimentally confirmed targets limited to the insulin gene at this time [4]. PSMD9 has, however, been linked with a wide range of human diseases and may, in some contexts, present a novel therapeutic target.
Banz et al. were the first to comment on the potential relevance of PSMD9 in cancer [7]. They demonstrated that treatment of breast cancer cell lines with the TGFβ family ligand Activin A stimulated PSMD9 expression together with increased levels of the signal transduction proteins Smad2, Smad3, and Smad4. PSMD9 knock-down was also associated with decreases in the levels of these Smads, implicating PSMD9 in the transmission of Activin A signalling through the regulation of Smad proteins. Activin A has roles in various potentially cancer-related pathways including those controlling inflammation, immunity, wound repair [8] and glucose metabolism [9], and has been implicated in pro- and anti-carcinogenic roles in several tissue types [10,11,12]. The team concluded that PSMD9 is a potential regulator of cell cycle progression. Further to this, we have recently published that low expression of PSMD9 within breast cancer tissues significantly associates with reduced incidence of local recurrence in patients receiving adjuvant radiotherapy but not in those treated without [13]. There is currently no standardized methodology by which to assess the radio-sensitivity of a tumour before treatment, although potential molecular markers have been reported, including p53 [14], Bcl-2 [15] and survivin [16]. Notably, we used a clinically routine immunohistochemistry methodology for our assessment, therefore identifying PSMD9 as a clinically-translatable biomarker of radiotherapy response. We also showed that in vitro depletion of PSMD9 in breast cancer cell lines acts to sensitize cells to radiotherapy, indicating that PSMD9 may have a direct role in cellular response to radiotherapy and may represent a novel target for radio-sensitizing drugs. As yet, it is unclear whether this influence of PSMD9 relates to its role in the proteasome or as a transcriptional co-regulator, however data showing the proteasome inhibitor Bortezomib to induce similar radio-sensitization [17] imply it may be the former.
Genome-wide linkage studies have identified a type 2 diabetes (T2D) locus called “non-insulin-dependent diabetes mellitus” (NIDDM) 2 [18,19]. The PSMD9 gene lies within the NIDDM2 locus and SNPs within the PSMD9 gene have been linked to the onset of T2D [20,21]. In a fine-mapping study Gragnoli and Cronsell performed direct sequencing in 237 Italian T2D patients and uncovered four rare PSMD9 mutations/variants that were associated with the disease [21]. Significant linkage of the PSMD9 SNPs rs74421874, rs3825172 and rs1043307/rs2514259 was later identified in late onset T2D [20] as well as early-onset T2D/maturity-onset diabetes of the young (MODY) type 3 [22] in Italian families. These PSMD9 SNPs have also been significantly linked to particular pathological features of this disease including T2D-neuropathy [23], microvascular [24] and macrovascular pathology [25] and hypertension [26]. The impact of these SNPs on PSMD9 function or expression has not yet been determined.
Thomas et al demonstrated that PSMD9 interacted with the basic helix-loop-helix (BHLH) transcription factors E12 and E47 via its PDZ domain and increased activation of enhancers of glucose-response within the insulin promoter. Functional investigations showed that increased expression of PSMD9 promoted pancreatic β-cell survival in vivo [4,27,28] whereas decreased expression of endogenous PSMD9 in rat insulinoma (INS-1) cells diminished insulin promoter activity suggesting that reduced expression of this protein may be linked with diabetes development [4]. In contrast, Volinic et al. found that increased expression of pancreatic PSMD9 led to insulin deficiency and diabetes in a murine model, and proposed a biphasic pattern of insulin promoter activation to explain their findings [29]. They showed that increasing levels of exogenous PSMD9 resulted in a shift from transcriptional activation to repression of the insulin promoter, and suggested that this may be the result of an incremental dosage effect in which increased levels first promote but then uncouple the optimal stoichiometric assembly of transcriptional regulatory complexes [29]. This suggested that individuals with reduced or augmented PSMD9 expression levels could be at an increased risk of diabetes and that modulation of this molecule may provide an effective therapeutic strategy.
It is now established that the inflammatory response is a major contributing pathway in both the development and progression of T2D. Chronic inflammation is triggered by the metabolic stress associated with hyperglycaemia in pancreatic islets that may lead to defective insulin secretion or response [30]. A study by Liu et al. into the role of the proteasome in diabetes demonstrated that the production of oxidative species due to hyperglycaemic conditions in diabetic models led to an increase in proteasome function in a PA700 dependant manner and that this caused subsequent activation of NF-κB [31], providing a further pointer for a mechanism of PSMD9’s influence as both a component of the PA700 complex and a known regulator of NF-κB.
Major depressive disorder (MDD) and generalised anxiety disorder (GAD) are the most prevalent mental health disorders in the UK, with approximately 1 in 5 people affected by these disorders in their lifetime. Familial connections have been indicated in these conditions, with heritability values of 30-40% [32]. Similarly, the rarer condition of Schizophrenia (SCZ), which affects approximately 1 in 100 individuals during their lifetimes [33], has a remarkably high heritability value of up to 80% [34], suggesting strong genetic contributions to the pathogenesis of these mental health illnesses. In addition, frequent connections have been made between T2D and mental health issues, with each an established comorbidity of the other [35,36]. Due to this, Gragnoli et al. further investigated the PSMD9 SNPs from their previous T2D studies and found that rs74421874, rs3825172 and rs1043307/rs2514259 were also significantly associated with both GAD [37] and MDD [38] development in T2D backgrounds. In support of this, a number of other studies have also linked PSMD9 SNPs to mental health disorders [39,40,41]. A number of PSMD9 SNPs has also been demonstrated to associate with individuals with clinically diagnosed cases of SCZ [41] and MDD [40]. Furthermore, PSMD9 rs1043307 has been correlated with a more positive response to tricyclic antidepressants [39]. This study predicted a link between PSMD9 and BHLH transcription factor binding pathway, again suggesting a possible transcriptional co-activator role for the protein, this time in the context of mental health. These data suggest that SNPs in PSMD9 may be used to identify individuals with increased risk of developing MDD or GAD in T2D patients or SCZ in the general population.
A definitive understanding of the molecular pathogenesis of these mental health disorders remains frustratingly evasive. Current theory focuses on dysregulation of various important proteins within the brain including stress hormones, cytokines and monoamines which may lead to decreased neuronal activity or inappropriate activation of inflammatory pathways [42,43,44,45]. In particular, a decrease in the expression of several specific proteins has been associated with SCZ with the most prominent example being the protein dysbindin [46,47,48,49]. Han et al. recently demonstrated that PSMD9 and dysbindin could be co-immunoprecipitated from mouse neuronal tissue [50], a direct interaction that supports a role for PSMD9 analogous to that previously discussed with respect to the degradation of Iκβα, in which PSMD9 may act to enhance delivery of dysbindin to the proteasome for degradation. Similarly, as part of their study into antidepressant treatment response, Wong et al. proposed that the PSMD9 SNP rs1043307 may be associated with a change in antigen processing or degradation that may lead to an alteration in the immune response of those individuals who responded more positively to anti-depressant treatments in comparison to those who do not [39]. Gragnoli et al. added to this with the suggestion that PSMD9, as a component of the 26S proteasome, may be involved in the immune-tolerance checking process [38] and dysregulation of this pathway through changes in proteasome-associated PSMD9 may lead to an autoimmune response that could trigger these mental health conditions.
Polycystic ovary syndrome (PCOS) is the most common endocrinopathy in women of reproductive age. Polycystic ovaries are generally enlarged and contain multiple follicles that do not undergo cyclic luteal differentiation, increasing infertility rates. Previous to their described work in breast cancer, Banz et el. examined PSMD9 function in ovary cells [51]. They performed in vitro experiments treating granulosa cells with activin A, and showed the same trends as those in cancer cells, with an increase in PSMD9 levels following activin A treatment and a decrease in Smad levels following PSMD9 knock down. The team suggested that PSMD9 may perform a regulatory role in granulosa cell proliferation and differentiation. Interestingly, PCOS demonstrates frequent correlations with T2D [52] and in many women raised levels of testosterone are indicative of the condition, which is reportedly triggered by hyperinsulinemia following a reduction in insulin sensitivity [53]. Therefore, through its relationship with the insulin signalling and activin A pathways, it is plausible that PSMD9 may become a potentially interesting new target in the treatment of PCOS and associated infertility, once data in humans with these disorders have been established.
PSMD9 has been associated with a group of chronic neurodegenerative diseases broadly termed synucleinopathies, which include Parkinson’s disease and multiple system atrophy [54]. These diseases are characterised by the presence of amyloid fibrils containing α-synuclein (α-Syn) in neuronal cells [55,56,57]. Over-expression of α-Syn in the dopaminergic neurons of the model organism C. elegans resulted in an up-regulation of PSMD9 expression [58]. The authors conclude that proteasome-associated PSMD9 is likely to be a component of the α-Syn degradation pathway that may have an important role in all synucleinopathies. Indeed, proteasomal disruptions have been strongly implicated in the formation of α-Syn aggregates [59,60] and although the mechanisms remain unclear, depletion of the 26S proteasome in transgenic mice has led to the development of Parkinson’s disease-like symptoms [61]. Further to this, aggregation of α-Syn is attributed to misfolding of the protein and Sangith et al. speculated that PSMD9 may have the capacity to perform quality control of protein misfolding. They suggest that proteasome associated PSMD9 located on the endoplasmic reticulum would have the ability to bind directly to the C-terminal motifs of exiting proteins via its PDZ domain, offering a way to circumvent the ubiquitination pathway [62]. Interestingly, although α-Syn has traditionally been studied within the brain, this protein is also present in pancreatic beta cells where it is proposed to have a regulatory role in insulin secretion [63,64] and might well be involved in the development of T2D. Epidemiological data suggests that T2D may be a risk factor in the development of PD [65] and may be associated with a more severe form of the disease [66] and, although conclusive data are limited, it has been suggested that 50-80% of PD patients may have abnormal glucose tolerance [67,68]. Through its association with metabolism and proteostasis PSMD9 may therefore provide a novel avenue of research for ageing related neurodegenerative diseases.
PSMD9 is a component of the 26S proteasome; a mediator of protein turnover, and has been shown to act as a transcriptional co-activator. PSMD9 has been associated with the pathogenic pathways of several important human disorders. Commonalities in the pathways regulating these disorders include changes in insulin secretion, alterations in protein turnover and immune dysregulation. Here, we suggest that targeting PSMD9 might be a novel therapy for the treatment of a range of diseases.
All authors declare no conflicts of interest in this paper.
[1] | A. C. Nica and E. T. Dermitzakis, Expression quantitative trait loci: present and future, Philos. Trans. R. Soc. Lond. B Biol. Sci., 368 (2013), 20120362. |
[2] | L. A. Hindorff, P. Sethupathy, H. A. Junkins, et al., Potential etiologic and functional implications of genome-wide association loci for human diseases and traits, Proc. Natl. Acad. Sci. U S A, 106 (2009), 9362-9367. |
[3] |
B. Hrdlickova, R. C. de Almeida, Z. Borek, et al., Genetic variation in the non-coding genome: Involvement of micro-rnas and long non-coding rnas in disease, Biochim. Biophys. Acta, 1842 (2014), 1910-1922. doi: 10.1016/j.bbadis.2014.03.011
![]() |
[4] | I. Ricaño-Ponce and C. Wijmenga, Mapping of immune-mediated disease genes, Annu. Rev. Genomics Hum. Genet., 14 (2013), 325-353. |
[5] | R. C. Jansen and J. P. Nap, Genetical genomics: the added value from segregation, Trends Genet., 17 (2001), 388-391. |
[6] | L. J. Carithers, K. Ardlie, M. Barcus, et al., A novel approach to high-quality postmortem tissue procurement: The gtex project, Biopreserv. Biobank., 13 (2015), 311-319, |
[7] | W. Cookson, L. Liang, G. Abecasis, et al., Mapping complex disease traits with global gene expression, Nat. Rev. Genet., 10 (2009), 184-194. |
[8] | A. C. Nica and E. T. Dermitzakis, Using gene expression to investigate the genetic basis of complex disorders, Hum. Mol. Genet., 17 (2008), R129-134. |
[9] | M. V. Rockman and L. Kruglyak, Genetics of global gene expression, Nat. Rev. Genet., 7 (2006), 862-872. |
[10] | F. A. Cubillos, V. Coustham and O. Loudet, Lessons from eqtl mapping studies: non-coding regions and their role behind natural phenotypic variation in plants, Curr. Opin. Plant. Biol., 15 (2012), 192-198. |
[11] | H. B. Fraser, A. M. Moses and E. E. Schadt, Evidence for widespread adaptive evolution of gene expression in budding yeast, Proc. Natl. Acad. Sci. U S A, 107 (2010), 2977-2982. |
[12] | A. L. Dixon, L. Liang, M. F. Moffatt, et al., A genome-wide association study of global gene expression, Nat. Genet., 39 (2007), 1202-1207. |
[13] | H. H. H. Göring, J. E. Curran, M. P. Johnson, et al., Discovery of expression qtls using large-scale transcriptional profiling in human lymphocytes, Nat. Genet., 39 (2007), 1208-1216. |
[14] | E. E. Schadt, C. Molony, E. Chudin, et al., Mapping the genetic architecture of gene expression in human liver, PLoS Biol., 6 (2008), e107. |
[15] |
A. Gerrits, Y. Li, B. M. Tesson, et al., Expression quantitative trait loci are highly sensitive to cellular differentiation state, PLoS Genet., 5 (2009), e1000692. doi: 10.1371/journal.pgen.1000692
![]() |
[16] | G. K. Chen and J. S. Witte, Enriching the analysis of genomewide association studies with hierarchical modeling, Am. J. Hum. Genet., 81 (2007), 397-404. |
[17] | X. Zhang, S. Huang, W. Sun, et al., Rapid and robust resampling-based multiple-testing correction with application in a genome-wide expression quantitative trait loci study, Genetics, 190 (2012), 1511-1520. |
[18] | M. P. Scott-Boyer, G. C. Imholte, A. Tayeb, et al., An integrated hierarchical bayesian model for multivariate eqtl mapping, Stat. Appl. Genet. Mol. Biol., 11 (2012), 10.1515/1544-6115.1760. |
[19] | O. Stegle, L. Parts, R. Durbin, et al., A bayesian framework to account for complex non-genetic factors in gene expression levels greatly increases power in eqtl studies, PLoS Comput. Biol., 6 (2010), e1000770. |
[20] | M. Stephens and D. J. Balding, Bayesian statistical methods for genetic association studies, Nat. Rev. Genet., 10 (2009), 681-690. |
[21] | J.-B. Veyrieras, S. Kudaravalli, S. Y. Kim, et al., High-resolution mapping of expression-qtls yields insight into human gene regulation, PLoS Genet., 4 (2008), e1000214. |
[22] | M. Banterle, L. Bottolo, S. Richardson, et al., Sparse variable and covariance selection for highdimensional seemingly unrelated bayesian regression, bioRxiv, 467019. |
[23] | G. C. Imholte, M.-P. Scott-Boyer, A. Labbe, et al., ibmq: a r/bioconductor package for integrated bayesian modeling of eqtl data, Bioinformatics, 29 (2013), 2797-2798. |
[24] | D. Duong, L. Gai, S. Snir, et al., Applying meta-analysis to genotype-tissue expression data from multiple tissues to identify eqtls and increase the number of egenes, Bioinformatics, 33 (2017), i67-i74. |
[25] | J. H. Sul, B. Han, C. Ye, et al., Effectively identifying eqtls from multiple tissues by combining mixed model and meta-analytic approaches, PLoS Genet., 9 (2013), e1003491. |
[26] | T. Flutre, X. Wen, J. Pritchard, et al., A statistical framework for joint eqtl analysis in multiple tissues, PLoS Genet., 9 (2013), e1003486. |
[27] | G. Li, A. A. Shabalin, I. Rusyn, et al., An empirical bayes approach for multiple tissue eqtl analysis, Biostatistics, 19 (2018), 391-406. |
[28] | A. Das, M. Morley, C. S. Moravec, et al., Bayesian integration of genetics and epigenetics detects causal regulatory snps underlying expression variability, Nat. Commun., 6 (2015), 8555. |
[29] | E. J. Chesler, L. Lu, J. Wang, et al., Webqtl: rapid exploratory analysis of gene expression and genetic networks for brain and behavior, Nat. Neurosci., 7 (2004), 485-486. |
[30] | J. Wang, R. W. Williams and K. F. Manly, Webqtl: web-based complex trait analysis, Neuroinformatics, 1 (2003), 299-308. |
[31] |
T. J. Phillips, M. Huson, C. Gwiazdon, et al., Effects of acute and repeated ethanol exposures on the locomotor activity of bxd recombinant inbred mice, Alcohol. Clin. Exp. Res., 19 (1995), 269-278. doi: 10.1111/j.1530-0277.1995.tb01502.x
![]() |
[32] |
B. Tabakoff, L. Saba, K. Kechris, et al., The genomic determinants of alcohol preference in mice, Mamm. Genome., 19 (2008), 352-365. doi: 10.1007/s00335-008-9115-z
![]() |
[33] |
B. J. Bennett, C. R. Farber, L. Orozco, et al., A high-resolution association mapping panel for the dissection of complex traits in mice, Genome. Res., 20 (2010), 281-290. doi: 10.1101/gr.099234.109
![]() |
[34] | R. Alberts, L. Lu, R. W. Williams, et al., Genome-wide analysis of the mouse lung transcriptome reveals novel molecular gene interaction networks and cell-specific expression signatures, Respir. Res., 12 (2011), 61. |
[35] |
C. Blauwendraat, M. Francescatto, J. R. Gibbs, et al., Comprehensive promoter level expression quantitative trait loci analysis of the human frontal lobe, Genome Med., 8 (2016), 65. doi: 10.1186/s13073-016-0320-1
![]() |
[36] |
J. A. Webster, J. R. Gibbs, J. Clarke, et al., Genetic control of human brain transcript expression in alzheimer disease, Am. J. Hum. Genet., 84 (2009), 445-458. doi: 10.1016/j.ajhg.2009.03.011
![]() |
[37] |
S. Lagarrigue, L. Martin, F. Hormozdiari, et al., Analysis of allele-specific expression in mouse liver by rna-seq: a comparison with cis-eqtl identified using genetic linkage, Genetics, 195 (2013), 1157-1166. doi: 10.1534/genetics.113.153882
![]() |
[38] | A. Gelman, J. B. Carlin, H. S. Stern, et al., Bayesian data analysis, vol. 2, Chapman & Hall/CRC Boca Raton, FL, USA, 2014. |
[39] | P. D. Hoff, A first course in Bayesian statistical methods, vol. 580, Springer, 2009. |
[40] | E. Lesaffre and A. B. Lawson, Bayesian biostatistics, John Wiley & Sons, 2012. |
[41] | X. Robin, N. Turck, A. Hainard, et al., proc: an open-source package for r and s+ to analyze and compare roc curves, BMC Bioinformatics, 12 (2011), 77. |
[42] | E. R. DeLong, D. M. DeLong and D. L. Clarke-Pearson, Comparing the areas under two or more correlated receiver operating characteristic curves: a nonparametric approach, Biometrics, 44 (1988), 837-845. |
[43] | S. A. Stouffer, E. A. Suchman, L. C. DeVinney, et al., The american soldier: Adjustment during army life, Princeton University Press, Vol. 1. |
[44] | L. T., On the combination of independent tests, Magyar Tud Akad Mat Kutato Int Közl. |
[45] | M. C. Whitlock, Combining probability from independent tests: the weighted z-method is superior to fisher's approach, J. Evol. Biol., 18 (2005), 1368-1373. |
[46] | D. Bates, M. Mächler, B. Bolker, et al., Fitting linear mixed-effects models using lme4, J. Stat. Software, 67 (2015), 1-48. |
[47] | R. V. Lenth, Least-squares means: The R package lsmeans, J. Stat. Software, 69 (2016), 1-33. |
[48] | RStudio Team, RStudio: Integrated Development Environment for R, RStudio, Inc., Boston, MA, 2015. |
[49] | R Core Team, R: A Language and Environment for Statistical Computing, R Foundation for Statistical Computing, Vienna, Austria, 2015. |
[50] | A. A. Shabalin, Matrix eqtl: ultra fast eqtl analysis via large matrix operations, Bioinformatics, 28 (2012), 1353-1358. |
[51] | H. Wickham, ggplot2: Elegant Graphics for Data Analysis, Springer-Verlag New York, 2009. |
[52] | R. C. Team, D. Wuertz, T. Setz, et al., fBasics: Rmetrics-Markets and Basic Statistics, 2014, R package version 3011.87. |
[53] | D. B. Dahl, xtable: Export Tables to LaTeX or HTML, 2016, R package version 1.8-2. |
[54] | S. Durinck, Y. Moreau, A. Kasprzyk, et al., Biomart and bioconductor: a powerful link between biological databases and microarray data analysis, Bioinformatics, 21 (2005), 3439-3440. |
[55] | H. Wickham, The split-apply-combine strategy for data analysis, J. Stat. Software, 40 (2011), 1-29. |
[56] | M. Dowle, A. Srinivasan, T. Short, et al., data.table: Extension of Data.frame, 2015, R package version 1.9.6. |
[57] | S. Anders and W. Huber, Differential expression analysis for sequence count data, Genome biol., 11 (2010), R106. |
[58] | C. W. Law, Y. Chen, W. Shi, et al., voom: Precision weights unlock linear model analysis tools for rna-seq read counts, Genome Biol., 15 (2014), R29. |
[59] | W. Sun, A statistical framework for eqtl mapping using rna-seq data, Biometrics, 68 (2012), 1-11. |
[60] | L. Bottolo and S. Richardson, Evolutionary stochastic search for bayesian model exploration, Bayesian Anal., 5 (2010), 583-618. |
[61] | N. A. Walter, S. K. McWeeney, S. T. Peters, et al., Snps matter: impact on detection of differential expression, Nature Methods, 4 (2007), 679. |
1. | Adrien Rousseau, Anne Bertolotti, Regulation of proteasome assembly and activity in health and disease, 2018, 19, 1471-0072, 697, 10.1038/s41580-018-0040-z | |
2. | Peter Larsson, Daniella Pettersson, Hanna Engqvist, Elisabeth Werner Rönnerman, Eva Forssell-Aronsson, Anikó Kovács, Per Karlsson, Khalil Helou, Toshima Z. Parris, Pan-cancer analysis of genomic and transcriptomic data reveals the prognostic relevance of human proteasome genes in different cancer types, 2022, 22, 1471-2407, 10.1186/s12885-022-10079-4 |