Next Article in Journal
Xeroderma Pigmentosum: Gene Variants and Splice Variants
Previous Article in Journal
Extracellular Vesicle-Derived DNA vs. CfDNA as a Biomarker for the Detection of Colon Cancer
Previous Article in Special Issue
Retinal Genomic Fabric Remodeling after Optic Nerve Injury
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

TWEAKing the Hippocampus: The Effects of TWEAK on the Genomic Fabric of the Hippocampus in a Neuropsychiatric Lupus Mouse Model

1
Personalized Genomics Laboratory, Center for Computational Systems Biology, Roy G Perry College of Engineering, Prairie View A&M University, Prairie View, TX 77446, USA
2
Dominick P. Purpura Department of Neuroscience, Albert Einstein College of Medicine, Bronx, NY 10461, USA
3
Department of Medicine (Rheumatology), Albert Einstein College of Medicine, Bronx, NY 10461, USA
4
Department of Pathology, New York Medical College, Valhalla, NY 10595, USA
5
Department of Microbiology & Immunology, Albert Einstein College of Medicine, Bronx, NY 10461, USA
6
Azrieli Faculty of Medicine, Bar-Ilan University, Zefat 1311502, Israel
7
Galilee Medical Center Research Institute, Nahariya 22100, Israel
*
Author to whom correspondence should be addressed.
Genes 2021, 12(8), 1172; https://doi.org/10.3390/genes12081172
Submission received: 25 June 2021 / Revised: 22 July 2021 / Accepted: 26 July 2021 / Published: 29 July 2021
(This article belongs to the Special Issue Genomic Fabric Remodeling in Neurological Disorders)

Abstract

:
Neuropsychiatric manifestations of systemic lupus erythematosus (SLE), specifically cognitive dysfunction and mood disorders, are widely prevalent in SLE patients, and yet poorly understood. TNF-like weak inducer of apoptosis (TWEAK) has previously been implicated in the pathogenesis of neuropsychiatric lupus (NPSLE), and we have recently shown its effects on the transcriptome of the cortex of the lupus-prone mice model MRL/lpr. As the hippocampus is thought to be an important focus of NPSLE processes, we explored the TWEAK-induced transcriptional changes that occur in the hippocampus, and isolated several genes (Dnajc28, Syne2, transthyretin) and pathways (PI3K-AKT, as well as chemokine-signaling and neurotransmission pathways) that are most differentially affected by TWEAK activation. While the functional roles of these genes and pathways within NPSLE need to be further investigated, an interesting link between neuroinflammation and neurodegeneration appears to emerge, which may prove to be a promising novel direction in NPSLE research.

1. Introduction

Neuropsychiatric Lupus (NPSLE) is one of the most prevalent manifestations of Sys-temic Lupus Erythematosus (SLE), occurring in up to 80% of lupus patients [1]. The 1999 American College of Rheumatology (ACR) ad hoc committee defined 19 clinical syndromes as manifestation of NPSLE, and those range from acute, overtly inflammatory presentations such as psychosis, transverse myelitis, and chorea to more subtle, non-specific symptoms such as headaches, mood disorders, and cognitive dysfunction [2]. Naturally, the more nebulous manifestations are significantly more prevalent (cognitive dysfunction and mood disorders range from 6.6% to 80%, while acute confusional state and cranial neuropathy affect 0.9–7% of patients [1]), but are more difficult to attribute directly to SLE, partly because they commonly occur regardless of systemic disease activity. Due to their apparent non-inflammatory presentation, cognitive dysfunction and mood disorders are thought to often be related to secondary causes, such as the patients being in a chronic disease state, neuro-affective medication use, or structural brain damage due to cerebrovascular disease, among others. Still, the fact that several lupus mice models pre-sent with similar neuropsychiatric clinical features [1]; advanced imaging studies of SLE patients demonstrate changes in neural networks that are common to patients compared with healthy controls [3]; in many patients, NPSLE manifestations occur early in the disease process, prior to the onset of chronic complications or treatment initiation [4,5], strongly support the notion that there are primary disease-related processes that lead to neurocognitive and psychiatric dysfunction in SLE. The non-specific nature of many of the NPSLE symptoms, and our current lack of understanding of the underlying processes make the diagnosis challenging, and directed treatment options scarce [6,7].
MRL/MpJ-Faslpr/lpr (MRL/lpr) lupus prone mice have a loss-of-function mutation in the Fas gene, superimposed on a complex MRL background. These mice manifest a systemic lupus-like phenotype, including anti-nuclear antibody formation, immune-complex mediated glomerular disease, and typical skin manifestations [8], as well as a range of cognitive and affective symptoms, with memory deficits and depression-like behavior [9]. Their neuropsychiatric behavior is evident as early as 8 weeks of age and peaks at a median of 16 weeks of age [10]. This protracted disease course together with systemic and neurologic presentation that is similar to human disease make this a popular model for NPSLE.
Tumor necrosis factor (TNF)-like weak inducer of apoptosis (TWEAK or TNFSF12), and its cognate receptor, Fn14 (TNFRSF12A), have been shown to play an important role in SLE in general [11,12,13] and NPSLE in particular. TWEAK, a secreted member of the TNF-ligand superfamily, and Fn14 are expressed in astrocytes, microglia, brain microvascular endothelial cells, and neurons; and their interaction activates pro-inflammatory cytokine production, among other effects [14,15,16]. MRL/lpr Fn14 knock-out mice (Fn14ko) display less depression and neurocognitive dysfunction than their background controls [17], and administration of TWEAK to the cerebrospinal fluid (CSF) of non-lupus mice induces neurocognitive behavioral changes [18]. Furthermore, NPSLE was associated with high TWEAK levels in the CSF of patients [19].
In a recent manuscript, we analyzed and compared the transcriptome of the cortex of the MRL/lpr lupus mouse model to the cortices of Fn14ko and non-lupus, MRL/MpJ (MRL/+), background control [20]. Utilizing a novel analysis method, incorporating both expression levels of the specific genes with their overall genetic importance to the model’s transcriptome, we identified several pathways that appear to have significant impact on the models’ differential phenotype. The MRL/lpr model displayed significant changes in neurotransmission processes, particularly in the dopaminergic pathway compared with Fn14ko and background MRL/+ controls. In addition, the Phosphoinositide 3-kinase (PI3K)-AKT intracellular pathway was identified as playing a role in the TWEAK-Fn14-mediated effects on NPSLE in the MRL/lpr mice. To further examine and validate the importance of these pathways, we proceeded to analyze the transcriptome of the cells in the hippocampus of the MRL/lpr mice compared with Fn14ko and background controls. The hippocampus has for long been implicated as a critical focus of neurocognitive and mood manifestation of NPSLE [21,22,23,24,25,26]. Thus, a careful analysis of its transcriptome may provide more specific and pertinent information regarding critical pathways in NPSLE, thereby identifying potential diagnostic markers, as well as promising treatment targets.

2. Materials and Methods

2.1. Animals

As described previously [20], MRL/lpr mice and Fn14ko (backcross generation #8) littermates were bred at Biogen Idec (Cambridge, MA, USA) and transferred to Albert Einstein College of Medicine (AECOM) at 8–10 weeks of age. Control MRL/+ mice were purchased from the Jackson Laboratory (Bar Harbor, ME, USA). Housing conditions were controlled, with a temperature of 21–23 °C and a 12:12 h light:dark cycle. All animal study protocols were approved by the institutional animal care and use committee (IACUC) at AECOM (protocol# 20170516).
There were 4 MRL/lpr, 4 Fn14ko and 4 MRL/+ mice used for this study. All mice were female and sacrificed at the diestrus phase of their hormonal cycle. At sacrifice, all were within one week of age (about 12 weeks old), and all were sacrificed within a 2-week time period. Following the sacrifice, the hippocampus was isolated, and immediately processed.

2.2. Transcriptomics

We used the experimental protocol and analyses presented in the previous report [20]. The hippocampus of each of the four mice from every group (MRL/+, MRL/lpr, and Fn14ko) was profiled separately. After reversed transcription in the presence of Cy3/Cy5 dUTP, total RNAs with different fluorescent labels from pairs of biological replicas were co-hybridized 17 h overnight at 65 °C with microarrays of 4 × 44 k Agilent 60-mer G2519F mouse chips.
The spots with corrupted pixels or with foreground fluorescence less than twice the background were eliminated from the analysis. Valid background subtracted foreground signals were normalized to the median and the results averaged separately for each group of spots probing redundantly the same gene.
Through the Genomic Fabric Paradigm (GFP) approach [27] we took full advantage of quantifying tens of thousands of genes at a time on four biological replicas. Thus, each quantified gene “i” in each region “B” (= cortex, hippocampus) and each phenotype “P” (= MRL/+, MRL/lpr, Fn14ko) was assigned the independent measures: average expression level (AVE), relative expression variability (REV) and expression correlation with each other gene (COR) according to the definitions (1)–(3):
A V E i ( B ; P ) = 1 R i k = 1 R i μ i , k ( B ; P ) = 1 R i k = 1 R i 1 4 ξ = 1 4 a i , k , ξ ( B ; P ) μ i , k ( B ; P ) , w h e r e :
B = cortex, hippocampus
R i = number of micro array spots probing redundantly gene i
a i , k , ξ ( B ; P ) = expression of gene “i” probed by spot “k” on biological replica “ ξ
a i , k , ξ ( B ; P ) = average expression of gene “i” probed by spot “k” on all biological replicas
R E V i ( B ; P ) = 1 2 r i χ 2 r i ; 0.975 + r i χ 2 r i ; 0.025 correction   coefficient 1 R i k = 1 R i s i k ( B ; P ) μ i k ( B ; P ) 2 pooled   CV × 100 %
χ 2 = chi-square score for r i degrees of freedom and α = 0.05
s i k = standard deviation of the expression of gene “i” probed by spot “k” on all biologicals replicas
r i = 4 R i 1 = number of degress of freedom
C O R i , j ( B ; P ) = k = 1 R i 1 4 ξ = 1 4 a i , k , ξ ( B ; P ) μ i , k ( B ; P ) 1 4 ξ = 1 4 a j , k , ξ ( B ; P ) μ j , k ( B ; P ) k = 1 R i 1 4 ξ = 1 4 a i , k , ξ ( B ; P ) μ i , k ( B ; P ) 2 k = 1 R i 1 4 ξ = 1 4 a j , k , ξ ( B ; P ) μ j , k ( B ; P ) 2   ( Pearson   correlation   cefficient )
Relative Expression Variability (REV) among biological replicas shows how much cellular homeostatic mechanisms control the transcript abundance against environmental slight random fluctuations [28], and the noise associated with the stochastic nature of the chemical reactions involved in the gene transcription. Genes critical for cell survival, proliferation, and integration in the multicellular structures are under strict control, while control of genes ensuring cell adaptation to environmental fluctuation is much more lenient.
Pearson pair-wise product-momentum correlation coefficient of expression levels (COR) with each other gene in that region and phenotype, or correlation with the same gene in other regions, reflects the Principle of Transcriptomic Stoichiometry [29], a generalization of Dalton’s Law of Multiple Proportions [30]. This principle states that genes networked in functional pathways are expressed in definite proportions, even under environmental fluctuations.
By combining REV and COR, we established gene hierarchy and identified Gene Master Regulator (GMR) in each region and each phenotype using Gene Commanding Height (GCH) [31,32,33]:
G C H i ( B ; P ) R E V ( B ; P ) R E V i ( B ; P ) exp 4 N j A L L , j i C O R i j ( B ; P ) 2 1
We compared AVEs of a gene in two phenotypes/regions and identified statistically significantly regulated/differentially expressed genes using the composite criterion of absolute fold change, exceeding the combined contributions of the expression variabilities, and the p-value of the heteroscedastic t-test being less than 0.05 (5). The expression ratio x was defined to clearly indicate the extent of the up- (positive ratio) or down- (negative ratio) regulation:
x i ( B 1 ; P 1   vs   B 2 ; P 2 ) > C U T i ( B 1 ; P 1   vs   B 2 ; P 2 ) = 1 + 1 100 2 R E V i ( B 1 ; P 1 ) 2 + R E V i ( B 2 ; P 2 ) 2 p i ( B 1 ; P 1   vs   B 2 ; P 2 ) < 0.05
where:
x i ( B 1 ; P 1   vs   B 2 ; P 2 ) μ i ( B 2 ; P 2 ) μ i ( B 1 ; P 1 ) ,   if   μ i ( B 2 ; P 2 ) μ i ( B 1 ; P 1 ) μ i ( B 1 ; P 1 ) μ i ( B 2 ; P 2 ) ,   if   μ i ( B 2 ; P 2 ) < μ i ( B 1 ; P 1 ) expression   ratio   ( negative   for   down - regulation )  
Traditionally, transcriptomic alterations are quantified by the percentages of up- and down-regulated genes. Not only is this method limited to only significantly regulated genes but it considers each affected gene as an equal +1 or −1 contributor.
For a more comprehensive characterization of the contribution of individual genes and functional pathways “Γ” to the expression difference between the compared phenotypes, we computed the Weighted Individual (Gene) Regulation (WIR) and the Weighted Pathway Regulation (WPR):
W I R i ( P 1   v s   P 2 ) A V E i ( P 2 ) x i ( P 1   v s   P 2 ) x i ( P 1   v s   P 2 ) regulation   sign x i ( P 1   v s   P 2 ) 1 absolute   net   fold - change 1 p i ( P 1   v s   P 2 ) confidence   of   the   regulation W P R Γ ( P 1   v s   P 2 ) 1 C a r d ( Γ ) i = 1 C a r d ( Γ ) W I R i Γ ( P 1   v s   P 2 ) 2
where:
P 1 = MRL / lpr ,   Fn 14 ko   P 2 =   MRL / + ,   Fn 14 ko   P 2 P 1 C a r d ( Γ ) = number of quantified genes in the partway Γ
As shown above, beyond the net-fold change, WIR takes into account the reference expression level of the gene and the statistical confidence (1 − p-value) of its regulation.

3. Results

Raw and normalized gene expression data were deposited and are publicly accessible at https://www.ncbi.nlm.nih.gov/geo/, (accessed on 25 March 2021) as GSE164140 (cortex) and GSE169486 (hippocampus). In total, we quantified 16,863 unigenes in each of all 24 profiled samples (2 regions × 3 phenotypes × 4 biological replicas). GFP approach turned expression data into: 101,178 average expression levels (AVE), 101,178 relative expression variabilities (REV), and 853,031,718 expression correlations (COR) among distinct genes in the same region and phenotype, and 50,679 between-regions correlations of the same genes. Thus, by fully exploiting the transcriptomic profiles, the workable experimental data was increased by 16,945 times of what is traditionally used in gene expression studies limited to only the average expression levels.

3.1. Independent Expression Characteristics of Individual Genes

For illustrative purposes, Figure 1 presents the average expression level (AVE), relative expression variability (REV) and correlation coefficient (COR) with Tnfrsf12a (Fn14) of the first 50 alphabetically ordered genes involved in the hippocampal PI3K-AKT signaling pathway of MRL/lpr, Fn14ko and MRL/+ mice. AVE indicates the expression level of each gene, and REV examines the genes’ degree of variability within each mouse model. It is assumed that genes that are critical for cell survival and function are highly preserved (low REV), and those that are meant to allow for adaptation would display higher REV. COR indicates the correlation of each gene’s expression with Tnfrsf12a expression. As Fn14ko mice do not express Tnfrsf12a, this model was not included in the COR analysis. Of note is the obvious independence of these three characteristics within each phenotype and the differences among the three models. Table 1 presents the 20 most expressed genes (highest AVE) in the hippocampus of each phenotype. Supplementary Tables S1 and S2 present the 20 most stably (low REV) and unstably (high REV) expressed genes in each phenotype and the corresponding values in the other two phenotypes.
As shown in Figure 1a,b, as well as Table 1, and in line with our previous cortex evaluation [20], Akt2 is highly expressed in the hippocampus of the MRL/lpr mice, and shows substantially increased REV in the MRL/lpr mice compared with both Fn14ko and MRL/+ controls. The apparent normalization of Akt2 expression in the Fn14ko model potentially points to Akt2 being directly related to TWEAK/Fn14 pathway activation. As discussed previously, the lack of correlation of Akt2 with Tnfrsf12a in the lupus mouse model (Figure 1c) may be due to a discrepancy between Tnfrsf12a gene expression (quantified here) and its activation.
Table 1 identifies Dnajc28 (DnaJ heat shock protein family 40 member C28) as the gene with the largest AVE in MRL/lpr compared with both Fn14ko and MRL/+. Dnajc28 was previously implicated in the pathogenesis of neurodegenerative diseases, such as Alzheimer’s Disease and Parkinson’s Disease [34,35,36], possibly as a protective protein that is expressed in high level in the setting of local injury or toxicity [37]. Its high level of expression in this case, therefore, is likely a compensatory mechanism to the TWEAK/Fn14-induced inflammation and its subsequent local damage.

3.2. Measures of Expression Regulation

Figure 2 presents the uniform contribution (reported as +1/−1, reflecting statistically significant up/down-regulated genes), expression ratio (or “fold-change”, negative for down-regulation) and weighted individual gene regulation (WIR) for the first 50 alphabetically ordered PI3K-AKT genes. The figure emphasizes the additive effect provided by the expression ratio (Figure 2b) and WIR (Figure 2c) measures, compared to the traditional uniform contribution analysis (Figure 2a). As discussed previously, the uniform contribution only identifies the genes that are significantly up-/down-regulated, without quantifying the level of regulation or its impact on the model’s transcriptome. Expression ratio discriminates the genes with respect to the magnitude of their regulation, while WIR (Figure 2c) weighs the net fold-change by the reference expression level and the statistical confidence of the regulation, thus providing a more comprehensive measure of transcriptomic impact.
Table 2 lists the 60 highest-contributor genes to the MRL/lpr model, compared with Fn14ko and MRL/+ controls, based on expression ratio and WIR values.
As demonstrated in Figure 2c and Table 2, Akt2 is among the most impactful genes within the PI3K-AKT pathway by an order of magnitude compared with other genes’ WIR. At the same time, as shown in Table 2, other genes also seem to play significant roles in regulating the genomic fabric of the different mice phenotypes through TWEAK/Fn14 activation (highlighted genes in Table 2 show similar effect of up- or down-regulation in both Fn14ko and MRL/+ controls compared with MRL/lpr indicating a likely role for the TWEAK/Fn14 pathway in these gene’s regulation). Of those, the most notable are Dnajc28, Syne2 (synaptic nuclear envelope 2) that are upregulated in the MRL/lpr model compared with both Fn14ko and MRL/+, as well as transthyretin (Ttr) that is downregulated in the lupus-prone mice compared with both controls. All three of these genes have been previously implicated in either CNS pathology [34,35,38,39,40,41,42] or autoimmunity [43], making them interesting targets of further study within the context of NPSLE.

3.3. Regulation of the Genomic Fabrics Responsible for Neurotransmission, Chemokine Signaling and PI3K-AKT Signaling Pathways

Figure 3a presents the percentages of up- and down-regulated genes within several relevant pathways, including PI3K-AKT (AKT), chemokine-signaling (CHS), and neurotransmission pathways. Figure 3b illustrates the weighted pathway regulation (WPR) scores of these pathways when comparing the three phenotypes. Notably, while there are no distinguishable differences between the models in the GABAergic (GAB), glutamatergic (GLU), and serotonergic (SER) pathways; the PI3K-AKT and chemokine-signaling, as well as the cholinergic (CHO) and dopaminergic (DOP) neurotransmission pathways are differentially regulated in the MRL/lpr mice compared with Fn14ko and MRL/+ controls. The low WPRs, when the 2 controls are compared with each other (Figure 3b; blue bars), indicate similar overall regulation of these pathways in the 2 control phenotypes; thus, pointing to a TWEAK/Fn14-mediated effect on those pathways that are differentially regulated.

3.4. Gene Hierarchy in the Hippocampus

The more important a gene is in preserving a particular phenotype, the more protected are its sequence and expression level by the cellular homeostatic mechanisms. In addition, critical genes play an important role in the regulation of major functional pathways, which can be evaluated through analyzing their expression coordination with the pathways’ genes. Combining the measure of expression control and expression coordination with other genes in the phenotype provides the Gene Commanding Height (GCH) score. We used GCH analysis as detailed in [20,31,32,33] to establish the gene hierarchy in the hippocampi of the three mouse models. Table 3 presents the top 15 genes with the highest GCH in each phenotype.
Of particular interest are the top GCH genes in the MRL/+ as these are important in the preservation of the healthy, non-lupus phenotype in this background control model. As expected, most of the genes encode for household proteins required for genetic material transcription and protection/repair regulation, intracellular structural and transportation mechanisms, and modulation of cell differentiation, proliferation, and signal transduction. Notably, one of the top GCH genes, Synaptotagmin XI (Syt11), plays an important role in regulating endocytosis and the vesicle-recycling process identified to be particularly significant in dopamine transmission, in addition to inhibiting cytokine secretion, such as interleukin-6 (IL6) and tumor necrosis factor (TNF), in macrophages and microglia [44]—both of these functions are potentially important in preventing the aberrant processes occurring in the lupus-prone brain.
The lack in overlap of the top GCH genes between the 3 phenotype is apparent. This is especially noteworthy when comparing the MRL/lpr model with the Fn14ko. As we discussed in our recent publication [20], while Fn14ko is thought to be physiologically similar to the MRL/lpr except for the knockout of 1 gene, transcriptomically, it appears that the differences between the models are more extensive and affect many more genes (in addition to the knocked down tweak) and functional pathways than would be expected.

3.5. Phenotype Dependence of the PI3K-AKT Pathway

The correlation expression levels of each of the AKT genes with their KEGG-established stimulators and inhibitors were analyzed (using Equation (3) in Materials and Methods section). Statistically significant (p < 0.05) COR of the AKT genes with their KEGG-determined stimulators and inhibitors [45] in the hippocampi of the three phenotypes are shown in Figure 4. According to the established and widely-used KEGG, the genes included in the groups labeled as PDK1, HSP90/Cd37, mTRC2, and TCL1 should be synergistically expressed (red lines) with all genes from the AKT group as they are thought to stimulate Akt gene expression within the pathway, while the genes labeled PP2A, CTMP and PHLPP should be antagonistically expressed (blue lines) with the Akt genes. As revealed by the coordination analysis in Figure 4, the genes are not uniformly correlated among the phenotypes, or consistently in-line with the KEGG-determined expected associations. The expression coordination is strongly dependent on the phenotype, contrary to the claimed universality of the KEGG-determined pathway. Interestingly, there are also substantial differences with expression coordination of the same genes between the hippocampus and cortex of the same phenotypes (Figure 7 in [20], as well as Figure 5 below).

3.6. Phenotype-Dependent Cortex-Hippocampus Synchronous Expression of Genes

Cognitive dysfunction in NPSLE has been frequently associated with hippocampal functional changes [21,22,23,24,25,26]. Previously, we presented data of gene expression analysis in the cortex of MRL/lpr, Fn14ko, and MRL/+ mice [20]. Figure 5 and Figure 6 illustrate the correlation between gene expression in the cortex, compared with the hippocampus of the 3 phenotypes, in the PI3K/AKT and neurotransmission pathways, respectively. These figures focus on the genes that are significantly expressed in-phase (expression levels are in the same direction, either enhanced or suppressed) and antiphase (gene expression is in opposite directions) among the two brain regions of each phenotype. As shown, there is enhanced in-phase gene expression between the regions in the MRL/lpr phenotype in about 10% of the analyzed genes in both pathways. However, the pattern of this association, both in the extent of general in-phase expression and which genes are synchronously expressed is phenotype-specific. Furthermore, when focusing on the asynchronously expressed genes between the 2 regions (blue lines), it is interesting that the lupus-prone MRL/lpr mice show the most asynchronous expression of genes in both analyzed pathways, while the MRL/+ controls show none.

4. Discussion

Analysis of the differential genomic expression and regulation of the hippocampi of lupus-prone mice highlights several important pathways that may play a role in inducing the NPSLE phenotype. More specifically, we focused on those pathways that appear to be TWEAK/Fn14-dependent, as TWEAK is an established key player in the pathogenesis of NPSLE [14,15,17,46], and an improved understanding of its downstream effects can provide important insight into the underlying pathologic processes, including potential targets for intervention. Our analysis highlights the importance of Akt2 in particular, and the PI3K-AKT pathway in general, in the genomic regulation of the MRL/lpr lupus-prone mouse. In addition to the PI3K-AKT pathway, other significant hippocampal pathways that seem to be associated with TWEAK/Fn14 activation are chemokine signaling, cholinergic, and dopaminergic neurotransmission pathways. In addition to these highlighted pathways, we also present genes that have significant impact on the MRL/lpr hippocampus transcriptome, such as Dnajc28, Syne2, and suppressed levels of Ttr, among others. Finally, we present evidence to differential pathway progression, or alterations to expected pathway sequences, as predicted by KEGG between the different mice phenotypes.
The hippocampus has long been implicated as a focus of NPSLE memory and learning impairment. In studies using advanced imaging techniques, the hippocampus is among the most consistently affected brain regions [3,47]. The specific localization of pathology to certain areas of the brain can be related to mechanical differences between the regions, such as increased regional vulnerability of the blood-brain barrier (BBB) increasing the local influx of inflammatory factors [48], as well as differences in local cell populations potentially making them selectively vulnerable to the NPSLE inflammatory drivers [49]. In addition, it has been demonstrated that different brain regions have variable cytokine profiles in the setting of NPSLE [50,51]. It was, therefore, important for us to specifically examine gene expression changes in the hippocampus, as it is likely more relevant to NPSLE neurocognitive changes than an evaluation that is not region-specific.
The neuropsychiatric manifestations of SLE have been shown to be triggered by an inflammatory process in a variety of contexts. Blood-brain-barrier (BBB) disruption allows infiltration of pathogenic antibodies to brains of mouse models, thereby causing neuropsychiatric manifestations [52]; tertiary lymph nodes at the site of the choroid plexus enable activated T- and B-cells to migrate to the brain parenchyma [53]; activated microglia cells are thought to play an active role in the local inflammatory process in this setting; abundance of pro-inflammatory cytokines is found in the CSF of NPSLE patients (reviewed in [1]). It is, therefore, not surprising that many of the NPSLE symptoms manifest in the setting of active SLE disease, and improve with immunosuppression. However, neurocognitive dysfunction in lupus patients has been a more elusive, less overtly inflammation-driven process. In human disease, it can often appear when overall disease activity is quiescent, and the symptoms often do not respond to immunosuppression. Previously, our group has shown in the MRL/lpr model that even with drastic attenuation of systemic inflammation, the neurocognitive behavioral phenotype persists, along with local cytokine production and neurodegeneration [54,55]. Accordingly, many of the most differentially regulated genes in the hippocampus of the lupus mouse model are ones that are related to neurodegenerative conditions, such as Alzheimer’s and Parkinson’s Disease, more than typical autoimmune, inflammation-related genes. At the same time, the prognosis of Alzheimer’s patients is associated with degree of systemic and local inflammation [56], and inflammatory processes, such as local cytokine production and T-cell infiltration, have been shown to play an important role in neurodegenerative conditions [57,58]. Thus, the interplay between inflammation and neurodegeneration is an important one to further explore.
Dnajc28 is a member of the Heat Shock Protein 40 (Hsp40) family. The members of this family of HSPs are thought to be molecular co-chaperones that bind to Hsp70 members, allowing them to interact with client proteins facilitating their proper folding, intracellular trafficking, and marking specific proteins for degradation [35]. A number of the Hsp40 family members have been implicated in familial forms of Parkinson’s Disease [35]. At the same time, several studies have shown neuroprotective effects of increased levels of extracellular HSPs, including Hsp40 and Hsp70, in several neurodegenerative disease models [59,60,61,62]. Thus, it is yet unclear whether the excess expression of Dnajc28 in the MRL/lpr model is damaging in and of itself, or whether its overexpression is a compensatory response to the stressed local environment. Similarly, Syne2 (spectrin repeat containing nuclear envelope protein 2), a member of the LINC (Linker of Nucleoskeleton and Cytoskeleton) complex that tethers the nuclear lamina to the cytoskeleton [63], has been previously identified to be associated with Alzheimer’s Disease (AD) and familial early-onset dementia [38,39]. Ttr is a systemic amyloid precursor that with abnormal folding due to genetic mutations or aging can lead to a form of systemic amyloidosis. In the setting of AD, it seems to have a protective effect by binding with the Aβ amyloid peptides and preventing fibril formation [42,64,65]. Interestingly, while historically Ttr was thought to be produced only by choroid plexus epithelial cells in the CNS, several groups have shown neuronal production of Ttr, particularly in the hippocampus and cortex, likely induced by Aβ precursor peptides as a local protective mechanism [42]. In our study, the increased expression of Ttr in the hippocampus of the MRL/lpr lupus model, compared with both the MRL/+ and Fn14ko controls, may indicate a compensatory, protective mechanism that was driven in the lupus mice due to stress or local inflammation (there is no evidence in the literature for increased production of amyloid precursor proteins in the context of SLE, or NPSLE [66], and so it is unclear what is driving the increased expression of Ttr here).
Similar to our previous findings in the MRL/lpr cortex, Akt2 appears to play a central role in the TWEAK/Fn14-induced effects of the hippocampus of the lupus model. Akt2 is one of 3 closely related serine/threonine-protein kinases (Akt1-3) that are key members of the PI3K-AKT pathway; regulating many essential processes, including cell proliferation and survival, growth, metabolism, and angiogenesis [67]. Studies of specific Akt1-3 null mice provided information regarding differences in roles and functions of the 3 isoforms. Akt3 is the most abundant in the brain, and plays a role in brain development and neurodegeneration [68]. Akt3-null mice have 25% smaller brain size [69], and increased susceptibility to demyelination in experimental autoimmune encephalitis (EAE), a widely used model for multiple sclerosis. Akt1 overexpression promotes enhanced myelination [70]. Conversely, Akt2-/- mice, lacking the Akt2 isoform that is known to be the most crucial in insulin-mediated glucose regulation [71], have normal brain size [69]. This lack of obvious brain manifestations in Akt2-null mice make our findings more curious. We show here that Akt2 is of the most differentially expressed genes in the hippocampus of MRL/lpr with significant WIR, indicating substantial effect on the MRL/lpr transcriptome and presumably phenotype. As Akt2 expression remains comparable to background control in the Fn14ko mice, the Akt2 overexpression in the MRL/lpr mice is likely mediated by TWEAK/Fn14 activation. TWEAK-Akt association has been shown in other systems, such as the heart [72], skeletal muscles [73], and tumors [74,75], including glioma, where Akt2 was specifically implicated in mediating TWEAK-induced cell survival [76]. Thus, it is conceivable that, while Akt2 does not play a major ongoing role in brain development and function, it can be activated in CNS inflammatory conditions such as NPSLE by TWEAK activity, possibly as a compensatory mechanism to improve cell survival in the context of inflammation-induced damage. Akt2-/- mice can be used to further explore whether Akt2 plays a direct pathogenic role contributing to the NPSLE-like phenotype, thereby clarifying whether Akt2 inhibition can be a viable and effective treatment in such scenarios.
In addition to the PI3K-AKT pathway, the chemokine-signaling pathway, as well as the cholinergic and dopaminergic neurotransmission pathways, were identified as important within the TWEAK/Fn14-mediated processes in the lupus-model hippocampus. Chemokine-signaling is an expected downstream event, especially in a cytokine-induced inflammatory context, as is the case in SLE in general, and in our experimental setting that focused on downstream effects of TWEAK/Fn14 activation. Highlighting the role of the cholinergic and dopaminergic neurotransmission pathways, as opposed to others, as significant in the pathogenesis of NPSLE is an important step in elucidating contributing mechanisms, as well as identifying targets of intervention for prevention and treatment.
Importantly, we also demonstrate that our current tools for genetic and pathway evaluations require further refinement. For example, KEGG-determined pathways widely relied upon for pathway identification and prediction studies, seem to be phenotype-dependent, and not universal across mice models. Furthermore, as we previously discussed in [20], gene knockout models, used in studies to isolate effects of one particular gene, display much more pervasive genetic and phenotypic differences from their controls. These observations should be taken into account when making predictions and trying to reach conclusions based on these methods. Of course, this presents an important limitation to our own study, and it should, therefore, be emphasized that this is an exploratory evaluation, and further confirmatory investigations need to be pursued prior to making definitive conclusions. These future studies should include protein analysis beyond gene expression, as there can be significant differences between a gene expression level and its actual translation. Another important limitation of our study is that we utilized only one mouse model: the MRL/lpr and its background control. Further studies need to be undertaken in other mice models to confirm our findings as relevant in NPSLE in general, and not just within the context of the MRL/lpr model.

5. Conclusions

SLE is a systemic autoimmune disease that drives inflammation in a myriad of organs, including the brain. While inflammation is a critical piece of the puzzle and probably its initial trigger, it is likely not the only driver of neuropsychiatric symptoms. Clinically, common and prevalent neuropsychiatric manifestations of SLE, such as cognitive dysfunction and mood disorders, often occur independent of disease activity (thus, not during times of increased systemic inflammation) and do not respond to immunosuppression. In this study, we focused on the effects of the pro-inflammatory cytokine TWEAK and its cognate receptor, Fn14, known to play an important role in NPSLE, on the transcriptome of the hippocampus of a lupus-prone mouse model. Notably, many of the most differentially regulated genes and pathways identified are those involved in neurodegenerative processes, as opposed to inflammatory ones, including Dnajc28, Ttr, and the PI3K-AKT pathway. Transitioning the focus of study to relevant neurodegenerative mechanisms, that potentially contribute to neuropsychiatric manifestations of SLE, may provide a clearer understanding of the underlying pathophysiology and enable the identification of effective treatment modalities that have, so far, remained elusive.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/genes12081172/s1, Table S1: The 20 most stably expressed genes (low relative expression variability, REV) in each phenotype; Table S2: The 20 most unstably expressed genes (high REV) in each phenotype.

Author Contributions

Conceptualization, D.A.I. and C.P.; methodology, D.A.I.; software, D.A.I.; validation, D.A.I. and N.S.; formal analysis, J.W., S.I. and N.S.; investigation, J.W. and S.I.; data curation, D.A.I.; writing—original draft preparation, D.A.I. and N.S..; writing—review and editing, D.A.I., C.P., N.S.; visualization, D.A.I.; supervision, C.P. and N.S.; project administration, C.P.; funding acquisition, C.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding. D.A.I. was supported by the Texas A&M University System Chancellor’s Research Initiative (CRI) funding for the Center for Computational Systems Biology at the Prairie View A&M University.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Raw and normalized gene expression data were deposited and are publicly accessible at https://www.ncbi.nlm.nih.gov/geo/, (accessed on 25 March 2021) as GSE164140 (cortex) and GSE169486 (hippocampus).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Schwartz, N.; Stock, A.D.; Putterman, C. Neuropsychiatric lupus: New mechanistic insights and future treatment directions. Nat. Rev. Rheumatol. 2019, 15, 137–152. [Google Scholar] [CrossRef] [PubMed]
  2. The American College of Rheumatology. Nomenclature and case definitions for neuropsychiatric lupus syndromes. Arthritis Rheum. 1999, 42, 599–608. [Google Scholar] [CrossRef]
  3. Mackay, M.; Vo, A.; Tang, C.C.; Small, M.; Anderson, E.W.; Ploran, E.J.; Storbeck, J.; Bascetta, B.; Kang, S.; Aranow, C.; et al. Metabolic and microstructural alterations in the SLE brain correlate with cognitive impairment. JCI Insight 2019, 4, e124002. [Google Scholar] [CrossRef]
  4. Steup-Beekman, G.M.; Zirkzee, E.J.; Cohen, D.; Gahrmann, B.M.; Emmer, B.J.; Steens, S.C.; Bollen, E.L.; van Buchem, M.A.; Huizinga, T.W. Neuropsychiatric manifestations in patients with systemic lupus erythematosus: Epidemiology and radiology pointing to an immune-mediated cause. Ann. Rheum. Dis. 2013, 72, ii76–ii79. [Google Scholar] [CrossRef]
  5. Bertsias, G.K.; Ioannidis, J.P.; Aringer, M.; Bollen, E.; Bombardieri, S.; Bruce, I.N.; Cervera, R.; Dalakas, M.; Doria, A.; Hanly, J.G.; et al. EULAR recommendations for the management of systemic lupus erythematosus with neuropsychiatric manifestations: Report of a task force of the EULAR standing committee for clinical affairs. Ann. Rheum. Dis. 2010, 69, 2074–2082. [Google Scholar] [CrossRef] [PubMed]
  6. Hanly, J.G.; Kozora, E.; Beyea, S.; Birnbaum, J. Nervous system disease in Systemic Lupus Erythematosus: Current status and future directions. Arthritis Rheumatol. 2018, 71, 33–42. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Kello, N.; Anderson, E.; Diamond, B. Cognitive dysfunction in Systemic Lupus Erythematosus: A case for initiating trials. Arthritis Rheumatol. 2019, 71, 1413–1425. [Google Scholar] [CrossRef] [PubMed]
  8. Du, Y.; Sanam, S.; Kate, K.; Mohan, C. Animal models of lupus and lupus nephritis. Curr. Pharm. Des. 2015, 21, 2320–2349. [Google Scholar] [CrossRef] [PubMed]
  9. Ballok, D.A. Neuroimmunopathology in a murine model of neuropsychiatric lupus. Brain Res. Rev. 2007, 54, 67–79. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Gulinello, M.; Putterman, C. The MRL/lpr mouse strain as a model for neuropsychiatric Systemic Lupus Erythematosus. J. Biomed. Biotechnol. 2011, 2011, 207504. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Doerner, J.L.; Wen, J.; Xia, Y.; Paz, K.B.; Schairer, D.; Wu, L.; Chalmers, S.A.; Izmirly, P.; Michaelson, J.S.; Burkly, L.C.; et al. TWEAK/Fn14 signaling involvement in the pathogenesis of cutaneous disease in the MRL/lpr model of spontaneous Lupus. J. Investig. Dermatol. 2015, 135, 1986–1995. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Michaelson, J.S.; Wisniacki, N.; Burkly, L.C.; Putterman, C. Role of TWEAK in lupus nephritis: A bench-to-bedside review. J. Autoimmun. 2012, 39, 130–142. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Zhao, Z.; Burkly, L.C.; Campbell, S.; Schwartz, N.; Molano, A.; Choudhury, A.; Eisenberg, R.A.; Michaelson, J.S.; Putterman, C. TWEAK/Fn14 interactions are instrumental in the pathogenesis of nephritis in the chronic graft-versus-host model of systemic lupus erythematosus. J. Immunol. 2007, 179, 7949–7958. [Google Scholar] [CrossRef] [Green Version]
  14. Desplat-Jego, S.; Varriale, S.; Creidy, R.; Terra, R.; Bernard, D.; Khrestchatisky, M.; Izui, S.; Chicheportiche, Y.; Boucraut, J. TWEAK is expressed by glial cells, induces astrocyte proliferation and increases EAE severity. J. Neuroimmunol. 2002, 133, 116–123. [Google Scholar] [CrossRef]
  15. Stock, A.D.; Wen, J.; Putterman, C. Neuropsychiatric Lupus, the Blood Brain Barrier, and the TWEAK/Fn14 Pathway. Front. Immunol. 2013, 4, 484. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Wen, J.; Doerner, J.; Weidenheim, K.; Xia, Y.; Stock, A.; Michaelson, J.S.; Baruch, K.; Deczkowska, A.; Gulinello, M.; Schwartz, M.; et al. TNF-like weak inducer of apoptosis promotes blood brain barrier disruption and increases neuronal cell death in MRL/lpr mice. J. Autoimmun. 2015, 60, 40–50. [Google Scholar] [CrossRef] [Green Version]
  17. Wen, J.; Xia, Y.; Stock, A.; Michaelson, J.S.; Burkly, L.C.; Gulinello, M.; Putterman, C. Neuropsychiatric disease in murine lupus is dependent on the TWEAK/Fn14 pathway. J. Autoimmun. 2013, 43, 44–54. [Google Scholar] [CrossRef] [Green Version]
  18. Wen, J.; Chen, C.H.; Stock, A.; Doerner, J.; Gulinello, M.; Putterman, C. Intracerebroventricular administration of TNF-like weak inducer of apoptosis induces depression-like behavior and cognitive dysfunction in non-autoimmune mice. Brain Behav. Immun. 2016, 54, 27–37. [Google Scholar] [CrossRef] [Green Version]
  19. Fragoso-Loyo, H.; Cabiedes, J.; Orozco-Narvaez, A.; Davila-Maldonado, L.; Atisha-Fregoso, Y.; Diamond, B.; Llorente, L.; Sanchez-Guerrero, J. Serum and cerebrospinal fluid autoantibodies in patients with neuropsychiatric lupus erythematosus. Implications for diagnosis and pathogenesis. PLoS ONE 2008, 3, e3347. [Google Scholar] [CrossRef] [Green Version]
  20. Iacobas, D.; Wen, J.; Iacobas, S.; Schwartz, N.; Putterman, C. Remodeling of neurotransmission, chemokine, and PI3K-AKT signaling genomic fabrics in neuropsychiatric Systemic Lupus Erythematosus. Genes 2021, 12, 251. [Google Scholar] [CrossRef]
  21. Appenzeller, S.; Carnevalle, A.D.; Li, L.M.; Costallat, L.T.; Cendes, F. Hippocampal atrophy in systemic lupus erythematosus. Ann. Rheum. Dis. 2006, 65, 1585–1589. [Google Scholar] [CrossRef] [Green Version]
  22. Ballok, D.A.; Woulfe, J.; Sur, M.; Cyr, M.; Sakic, B. Hippocampal damage in mouse and human forms of systemic autoimmune disease. Hippocampus 2004, 14, 649–661. [Google Scholar] [CrossRef] [Green Version]
  23. Kozora, E.; Brown, M.S.; Filley, C.M.; Zhang, L.; Miller, D.E.; West, S.G.; Pelzman, J.; Arciniegas, D.B. Memory impairment associated with neurometabolic abnormalities of the hippocampus in patients with non-neuropsychiatric systemic lupus erythematosus. Lupus 2011, 20, 598–606. [Google Scholar] [CrossRef]
  24. Lauvsnes, M.B.; Beyer, M.K.; Kvaloy, J.T.; Greve, O.J.; Appenzeller, S.; Kvivik, I.; Harboe, E.; Tjensvoll, A.B.; Goransson, L.G.; Omdal, R. Association of hippocampal atrophy with cerebrospinal fluid antibodies against the NR2 subtype of the N-methyl-D-aspartate receptor in patients with systemic lupus erythematosus and patients with primary Sjogren’s syndrome. Arthritis Rheumatol. 2014, 66, 3387–3394. [Google Scholar] [CrossRef]
  25. Liu, S.; Cheng, Y.; Zhao, Y.; Lai, A.; Lv, Z.; Xie, Z.; Upreti, B.; Wang, X.; Xu, X.; Luo, C.; et al. Hippocampal atrophy in Systemic Lupus Erythematosus patients without major neuropsychiatric manifestations. J. Immunol. Res. 2020, 2020, 2943848. [Google Scholar] [CrossRef] [PubMed]
  26. Zimmermann, N.; Correa, D.G.; Kubo, T.A.; Netto, T.M.; Pereira, D.B.; Fonseca, R.P.; Gasparetto, E.L. Global cognitive impairment in Systemic Lupus Erythematosus patients: A structural MRI study. Clin. Neuroradiol. 2017, 27, 23–29. [Google Scholar] [CrossRef] [PubMed]
  27. Iacobas, D.A.; Iacobas, S.; Stout, R.F.; Spray, D.C. Cellular environment remodels the genomic fabrics of functional pathways in astrocytes. Genes 2020, 11, 520. [Google Scholar] [CrossRef] [PubMed]
  28. Victorino, P.H.; Marra, C.; Iacobas, D.A.; Iacobas, S.; Spray, D.C.; Linden, R.; Adesse, D.; Petrs-Silva, H. Retinal genomic fabric remodeling after optic nerve injury. Genes 2021, 12, 403. [Google Scholar] [CrossRef]
  29. Iacobas, D.A.; Iacobas, S.; Lee, P.R.; Cohen, J.E.; Fields, R.D. Coordinated activity of transcriptional networks responding to the pattern of action potential firing in neurons. Genes 2019, 10, 754. [Google Scholar] [CrossRef] [Green Version]
  30. Law of Multiple Proportions. Available online: https://www.britannica.com/science/law-of-multiple-proportions (accessed on 24 April 2021).
  31. Iacobas, D.A. Biomarkers, master regulators and genomic fabric remodeling in a case of Papillary Thyroid Carcinoma. Genes 2020, 11, 1030. [Google Scholar] [CrossRef]
  32. Iacobas, S.; Ede, N.; Iacobas, D.A. The Gene Master Regulators (GMR) approach provides legitimate targets for personalized, time-sensitive cancer gene therapy. Genes 2019, 10, 560. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Iacobas, D.A.; Tuli, N.Y.; Iacobas, S.; Rasamny, J.K.; Moscatello, A.; Geliebter, J.; Tiwari, R.K. Gene master regulators of papillary and anaplastic thyroid cancers. Oncotarget 2018, 9, 2410–2424. [Google Scholar] [CrossRef] [Green Version]
  34. Zarouchlioti, C.; Parfitt, D.A.; Li, W.; Gittings, L.M.; Cheetham, M.E. DNAJ Proteins in neurodegeneration: Essential and protective factors. Philos. Trans. R. Soc. Lond. B Biol. Sci. 2018, 373, 20160534. [Google Scholar] [CrossRef] [Green Version]
  35. Hasegawa, T.; Yoshida, S.; Sugeno, N.; Kobayashi, J.; Aoki, M. DnaJ/Hsp40 family and Parkinson’s Disease. Front. Neurosci. 2017, 11, 743. [Google Scholar] [CrossRef]
  36. Park, S.K.; Arslan, F.; Kanneganti, V.; Barmada, S.J.; Purushothaman, P.; Verma, S.C.; Liebman, S.W. Overexpression of a conserved HSP40 chaperone reduces toxicity of several neurodegenerative disease proteins. Prion 2018, 12, 16–22. [Google Scholar] [CrossRef]
  37. Tanaka, S.; Kitagawa, K.; Ohtsuki, T.; Yagita, Y.; Takasawa, K.; Hori, M.; Matsumoto, M. Synergistic induction of HSP40 and HSC70 in the mouse hippocampal neurons after cerebral ischemia and ischemic tolerance in gerbil hippocampus. J. Neurosci. Res. 2002, 67, 37–47. [Google Scholar] [CrossRef] [PubMed]
  38. Guttula, S.V.; Allam, A.; Gumpeny, R.S. Analyzing microarray data of Alzheimer’s using cluster analysis to identify the biomarker genes. Int. J. Alzheimers Dis. 2012, 2012, 649456. [Google Scholar] [CrossRef] [Green Version]
  39. Alexander, J.; Kalev, O.; Mehrabian, S.; Traykov, L.; Raycheva, M.; Kanakis, D.; Drineas, P.; Lutz, M.I.; Strobel, T.; Penz, T.; et al. Familial early-onset dementia with complex neuropathologic phenotype and genomic background. Neurobiol. Aging 2016, 42, 199–204. [Google Scholar] [CrossRef] [PubMed]
  40. Alemi, M.; Silva, S.C.; Santana, I.; Cardoso, I. Transthyretin stability is critical in assisting β amyloid clearance—Relevance of transthyretin stabilization in Alzheimer’s disease. CNS Neurosci. Ther. 2017, 23, 605–619. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Giao, T.; Saavedra, J.; Cotrina, E.; Quintana, J.; Llop, J.; Arsequell, G.; Cardoso, I. Undiscovered roles for Transthyretin: From a transporter protein to a new therapeutic target for Alzheimer’s disease. Int. J. Mol. Sci. 2020, 21, 2075. [Google Scholar] [CrossRef] [Green Version]
  42. Li, X.; Masliah, E.; Reixach, N.; Buxbaum, J.N. Neuronal production of transthyretin in human and murine Alzheimer’s disease: Is it protective? J. Neurosci. 2011, 31, 12483–12490. [Google Scholar] [CrossRef] [Green Version]
  43. Zhao, J.; Li, Y.; Hu, Y.; Chen, C.; Zhou, Y.; Tao, Y.; Guo, M.; Qin, N.; Xu, L. MicroRNAs expression profile in CCR6(+) regulatory T cells. Peer J. 2014, 2, e575. [Google Scholar] [CrossRef] [Green Version]
  44. GeneCards. The Human Gene Database: Synaptotagmin 11. Available online: https://www.genecards.org/cgi-bin/carddisp.pl?gene=SYT11 (accessed on 31 May 2021).
  45. KEGG-Derived PI3K-Akt Signaling Pathway. Available online: https://www.genome.jp/kegg-bin/show_pathway?mmu04151 (accessed on 2 September 2020).
  46. Fragoso-Loyo, H.; Atisha-Fregoso, Y.; Nunez-Alvarez, C.A.; Llorente, L. Utility of TWEAK to assess neuropsychiatric disease activity in systemic lupus erhytematosus. Lupus 2016, 25, 364–369. [Google Scholar] [CrossRef]
  47. Mackay, M.; Tang, C.C.; Vo, A. Advanced neuroimaging in neuropsychiatric systemic lupus erythematosus. Curr. Opin. Neurol. 2020, 33, 353–361. [Google Scholar] [CrossRef]
  48. Chi, J.M.; Mackay, M.; Hoang, A.; Cheng, K.; Aranow, C.; Ivanidze, J.; Volpe, B.; Diamond, B.; Sanelli, P.C. Alterations in blood-brain barrier permeability in patients with Systemic Lupus Erythematosus. AJNR Am. J. Neuroradiol. 2019, 40, 470–477. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Leng, K.; Li, E.; Eser, R.; Piergies, A.; Sit, R.; Tan, M.; Neff, N.; Li, S.H.; Rodriguez, R.D.; Suemoto, C.K.; et al. Molecular characterization of selectively vulnerable neurons in Alzheimer’s disease. Nat. Neurosci. 2021, 24, 276–287. [Google Scholar] [CrossRef]
  50. Tomita, M.; Holman, B.J.; Santoro, T.J. Aberrant cytokine gene expression in the hippocampus in murine systemic lupus erythematosus. Neurosci. Lett. 2001, 302, 129–132. [Google Scholar] [CrossRef]
  51. Tomita, M.; Holman, B.J.; Williams, L.S.; Pang, K.C.; Santoro, T.J. Cerebellar dysfunction is associated with overexpression of proinflammatory cytokine genes in lupus. J. Neurosci. Res. 2001, 64, 26–33. [Google Scholar] [CrossRef] [PubMed]
  52. Kowal, C.; DeGiorgio, L.A.; Nakaoka, T.; Hetherington, H.; Huerta, P.T.; Diamond, B.; Volpe, B.T. Cognition and immunity; antibody impairs memory. Immunity 2004, 21, 179–188. [Google Scholar] [CrossRef] [Green Version]
  53. Stock, A.D.; Der, E.; Gelb, S.; Huang, M.; Weidenheim, K.; Ben-Zvi, A.; Putterman, C. Tertiary lymphoid structures in the choroid plexus in neuropsychiatric lupus. JCI Insight 2019, 4. [Google Scholar] [CrossRef] [Green Version]
  54. Stock, A.D.; Wen, J.; Doerner, J.; Herlitz, L.C.; Gulinello, M.; Putterman, C. Neuropsychiatric systemic lupus erythematosus persists despite attenuation of systemic disease in MRL/lpr mice. J. Neuroinflamm. 2015, 12, 205. [Google Scholar] [CrossRef] [Green Version]
  55. Wen, J.; Doerner, J.; Chalmers, S.; Stock, A.; Wang, H.; Gullinello, M.; Shlomchik, M.J.; Putterman, C. B cell and/or autoantibody deficiency do not prevent neuropsychiatric disease in murine systemic lupus erythematosus. J. Neuroinflamm. 2016, 13, 73. [Google Scholar] [CrossRef] [Green Version]
  56. Simen, A.A.; Bordner, K.A.; Martin, M.P.; Moy, L.A.; Barry, L.C. Cognitive dysfunction with aging and the role of inflammation. Ther. Adv. Chronic. Dis. 2011, 2, 175–195. [Google Scholar] [CrossRef] [Green Version]
  57. Togo, T.; Akiyama, H.; Iseki, E.; Kondo, H.; Ikeda, K.; Kato, M.; Oda, T.; Tsuchiya, K.; Kosaka, K. Occurrence of T cells in the brain of Alzheimer’s disease and other neurological diseases. J. Neuroimmunol. 2002, 124, 83–92. [Google Scholar] [CrossRef]
  58. Gate, D.; Saligrama, N.; Leventhal, O.; Yang, A.C.; Unger, M.S.; Middeldorp, J.; Chen, K.; Lehallier, B.; Channappa, D.; De Los Santos, M.B.; et al. Clonally expanded CD8 T cells patrol the cerebrospinal fluid in Alzheimer’s disease. Nature 2020, 577, 399–404. [Google Scholar] [CrossRef] [PubMed]
  59. Gao, X.; Carroni, M.; Nussbaum-Krammer, C.; Mogk, A.; Nillegoda, N.B.; Szlachcic, A.; Guilbride, D.L.; Saibil, H.R.; Mayer, M.P.; Bukau, B. Human Hsp70 disaggregase reverses Parkinson’s-linked α-synuclein amyloid fibrils. Mol. Cell 2015, 59, 781–793. [Google Scholar] [CrossRef] [Green Version]
  60. Kakimura, J.; Kitamura, Y.; Takata, K.; Umeki, M.; Suzuki, S.; Shibagaki, K.; Taniguchi, T.; Nomura, Y.; Gebicke-Haerter, P.J.; Smith, M.A.; et al. Microglial activation and amyloid-β clearance induced by exogenous heat-shock proteins. FASEB J. 2002, 16, 601–603. [Google Scholar] [CrossRef] [PubMed]
  61. Popiel, H.A.; Takeuchi, T.; Fujita, H.; Yamamoto, K.; Ito, C.; Yamane, H.; Muramatsu, S.; Toda, T.; Wada, K.; Nagai, Y. Hsp40 gene therapy exerts therapeutic effects on polyglutamine disease mice via a non-cell autonomous mechanism. PLoS ONE 2012, 7, e51069. [Google Scholar] [CrossRef] [PubMed]
  62. Takeuchi, T.; Suzuki, M.; Fujikake, N.; Popiel, H.A.; Kikuchi, H.; Futaki, S.; Wada, K.; Nagai, Y. Intercellular chaperone transmission via exosomes contributes to maintenance of protein homeostasis at the organismal level. Proc. Natl. Acad. Sci. USA 2015, 112, E2497–E2506. [Google Scholar] [CrossRef] [Green Version]
  63. GeneCards. The Human Gene Database: SYNE2. Available online: https://www.genecards.org/cgi-bin/carddisp.pl?gene=SYNE2 (accessed on 31 May 2021).
  64. Buxbaum, J.N.; Ye, Z.; Reixach, N.; Friske, L.; Levy, C.; Das, P.; Golde, T.; Masliah, E.; Roberts, A.R.; Bartfai, T. Transthyretin protects Alzheimer’s mice from the behavioral and biochemical effects of Abeta toxicity. Proc. Natl. Acad. Sci. USA 2008, 105, 2681–2686. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Schwarzman, A.L.; Gregori, L.; Vitek, M.P.; Lyubski, S.; Strittmatter, W.J.; Enghilde, J.J.; Bhasin, R.; Silverman, J.; Weisgraber, K.H.; Coyle, P.K.; et al. Transthyretin sequesters amyloid β protein and prevents amyloid formation. Proc. Natl. Acad. Sci. USA 1994, 91, 8368–8372. [Google Scholar] [CrossRef] [Green Version]
  66. Trysberg, E.; Hoglund, K.; Svenungsson, E.; Blennow, K.; Tarkowski, A. Decreased levels of soluble amyloid β-protein precursor and β-amyloid protein in cerebrospinal fluid of patients with systemic lupus erythematosus. Arthritis Res. Ther. 2004, 6, R129–R136. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Engelman, J.A.; Luo, J.; Cantley, L.C. The evolution of phosphatidylinositol 3-kinases as regulators of growth and metabolism. Nat. Rev. Genet. 2006, 7, 606–619. [Google Scholar] [CrossRef]
  68. Tschopp, O.; Yang, Z.Z.; Brodbeck, D.; Dummler, B.A.; Hemmings-Mieszczak, M.; Watanabe, T.; Michaelis, T.; Frahm, J.; Hemmings, B.A. Essential role of protein kinase B γ (PKB γ/Akt3) in postnatal brain development but not in glucose homeostasis. Development 2005, 132, 2943–2954. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Easton, R.M.; Cho, H.; Roovers, K.; Shineman, D.W.; Mizrahi, M.; Forman, M.S.; Lee, V.M.; Szabolcs, M.; de Jong, R.; Oltersdorf, T.; et al. Role for Akt3/protein kinase Bgamma in attainment of normal brain size. Mol. Cell. Biol. 2005, 25, 1869–1878. [Google Scholar] [CrossRef] [Green Version]
  70. Flores, A.I.; Narayanan, S.P.; Morse, E.N.; Shick, H.E.; Yin, X.; Kidd, G.; Avila, R.L.; Kirschner, D.A.; Macklin, W.B. Constitutively active Akt induces enhanced myelination in the CNS. J. Neurosci. 2008, 28, 7174–7183. [Google Scholar] [CrossRef] [Green Version]
  71. Cho, H.; Mu, J.; Kim, J.K.; Thorvaldsen, J.L.; Chu, Q.; Crenshaw, E.B., 3rd; Kaestner, K.H.; Bartolomei, M.S.; Shulman, G.I.; Birnbaum, M.J. Insulin resistance and a diabetes mellitus-like syndrome in mice lacking the protein kinase Akt2 (PKB β). Science 2001, 292, 1728–1731. [Google Scholar] [CrossRef] [PubMed]
  72. Yang, B.; Yan, P.; Gong, H.; Zuo, L.; Shi, Y.; Guo, J.; Guo, R.; Xie, J.; Li, B. TWEAK protects cardiomyocyte against apoptosis in a PI3K/AKT pathway dependent manner. Am. J. Transl. Res. 2016, 8, 3848–3860. [Google Scholar]
  73. Dogra, C.; Changotra, H.; Wedhas, N.; Qin, X.; Wergedal, J.E.; Kumar, A. TNF-related weak inducer of apoptosis (TWEAK) is a potent skeletal muscle-wasting cytokine. FASEB J. 2007, 21, 1857–1869. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Zhang, F.; Xu, M.; Yin, X.; Guo, H.; Zhang, B.; Wang, Y.; Xiao, J.; Zou, X.; Zhang, M.; Zhuge, Y. TWEAK promotes hepatic stellate cell migration through activating EGFR/Src and PI3K/AKT pathways. Cell Biol. Int. 2019, 44, 278–285. [Google Scholar] [CrossRef]
  75. Xu, R.D.; Feng, F.; Yu, X.S.; Liu, Z.D.; Lao, L.F. miR-149-5p inhibits cell growth by regulating TWEAK/Fn14/PI3K/AKT pathway and predicts favorable survival in human osteosarcoma. Int J. Immunopathol. Pharmacol. 2018, 32, 2058738418786656. [Google Scholar] [CrossRef] [PubMed]
  76. Fortin, S.P.; Ennis, M.J.; Savitch, B.A.; Carpentieri, D.; McDonough, W.S.; Winkles, J.A.; Loftus, J.C.; Kingsley, C.; Hostetter, G.; Tran, N.L. Tumor necrosis factor-like weak inducer of apoptosis stimulation of glioma cell survival is dependent on Akt2 function. Mol. Cancer Res. 2009, 7, 1871–1881. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. The independent characteristics of the first 50 alphabetically ordered hippocampal PI3K-AKT genes in the three phenotypes. (a) Average expression level (AVE); (b) Relative expression variability (REV); (c) Pair-wise product-momentum correlation coefficient (COR) of PI3K-AKT genes with Tnfrsf12a. Red/green arrows indicate the statistically significant (p < 0.05) correlations in the MRL/lpr and MRL/+ mice. With AVE = 262 and REV = 193, Akt2 is the most remarkable gene of the selection in the hippocampus of the MRL/lpr mice.
Figure 1. The independent characteristics of the first 50 alphabetically ordered hippocampal PI3K-AKT genes in the three phenotypes. (a) Average expression level (AVE); (b) Relative expression variability (REV); (c) Pair-wise product-momentum correlation coefficient (COR) of PI3K-AKT genes with Tnfrsf12a. Red/green arrows indicate the statistically significant (p < 0.05) correlations in the MRL/lpr and MRL/+ mice. With AVE = 262 and REV = 193, Akt2 is the most remarkable gene of the selection in the hippocampus of the MRL/lpr mice.
Genes 12 01172 g001
Figure 2. Measures of differential expression of the first 50 alphabetically ordered hippocampal PI3K-AKT genes in the MRL/lpr mouse with respect to the two controls. (a) Uniform contribution (+1 or −1) of the significantly regulated genes. (b) Fold-change (negative for down-regulation) of all genes within the selection. (c) Weighted Individual (gene) Regulation (WIR). In addition to the net fold-change, WIR considers also the reference expression level of the gene and the statistical confidence (1 − p-value) of its regulation. Taken together, it is a comprehensive measure of the relative contribution of the gene to the model’s transcriptome.
Figure 2. Measures of differential expression of the first 50 alphabetically ordered hippocampal PI3K-AKT genes in the MRL/lpr mouse with respect to the two controls. (a) Uniform contribution (+1 or −1) of the significantly regulated genes. (b) Fold-change (negative for down-regulation) of all genes within the selection. (c) Weighted Individual (gene) Regulation (WIR). In addition to the net fold-change, WIR considers also the reference expression level of the gene and the statistical confidence (1 − p-value) of its regulation. Taken together, it is a comprehensive measure of the relative contribution of the gene to the model’s transcriptome.
Genes 12 01172 g002
Figure 3. Overall transcriptomic differences in the analyzed pathways among the three phenotypes. (a) Percentages of up- and down-regulated genes in selected pathways for indicated comparisons. (b) Weighted Pathway Regulation (WPR) of selected pathways. AKT, PI3K-AKT pathway; CHS, chemokine signaling pathway; CHO, cholinergic neurotransmission; DOP, dopaminergic neurotransmission; GAB, GABAergic neurotransmission, GLU, glutamatergic neurotransmission; SER, serotonergic neurotransmission.
Figure 3. Overall transcriptomic differences in the analyzed pathways among the three phenotypes. (a) Percentages of up- and down-regulated genes in selected pathways for indicated comparisons. (b) Weighted Pathway Regulation (WPR) of selected pathways. AKT, PI3K-AKT pathway; CHS, chemokine signaling pathway; CHO, cholinergic neurotransmission; DOP, dopaminergic neurotransmission; GAB, GABAergic neurotransmission, GLU, glutamatergic neurotransmission; SER, serotonergic neurotransmission.
Genes 12 01172 g003
Figure 4. Phenotype-dependent transcriptomic network of AKT genes with their KEGG-derived activator and inhibitor genes in the hippocampus. Red/blue lines depict statistically significant (p < 0.05) expression synergism/antagonism between the linked genes. Red line indicates synergism, and blue line antagonism between the 2 genes.
Figure 4. Phenotype-dependent transcriptomic network of AKT genes with their KEGG-derived activator and inhibitor genes in the hippocampus. Red/blue lines depict statistically significant (p < 0.05) expression synergism/antagonism between the linked genes. Red line indicates synergism, and blue line antagonism between the 2 genes.
Genes 12 01172 g004
Figure 5. Significant in-phase and anti-phase expression of PI3K-AKT genes between the cortex and hippocampus in each phenotype. Red/blue lines indicate the genes that are significantly expressed in-phase and anti-phase between the two regions. Red numbers indicate percentages of in-phase expression, and blue numbers indicate percentages of anti-phase expression.
Figure 5. Significant in-phase and anti-phase expression of PI3K-AKT genes between the cortex and hippocampus in each phenotype. Red/blue lines indicate the genes that are significantly expressed in-phase and anti-phase between the two regions. Red numbers indicate percentages of in-phase expression, and blue numbers indicate percentages of anti-phase expression.
Genes 12 01172 g005
Figure 6. Significant in-phase and anti-phase expression of neurotransmission genes between the cortex and hippocampus in each phenotype. Red/blue lines indicate the genes that are significantly expressed in-phase and antiphase between the two regions. Red numbers indicate percentages of in-phase, and blue numbers of anti-phase expressions.
Figure 6. Significant in-phase and anti-phase expression of neurotransmission genes between the cortex and hippocampus in each phenotype. Red/blue lines indicate the genes that are significantly expressed in-phase and antiphase between the two regions. Red numbers indicate percentages of in-phase, and blue numbers of anti-phase expressions.
Genes 12 01172 g006
Table 1. The 20 most highly expressed genes in the hippocampus of MRL/lpr mice. Average values of gene expression (AVE) are expressed in the table as multiples of the median gene expression.
Table 1. The 20 most highly expressed genes in the hippocampus of MRL/lpr mice. Average values of gene expression (AVE) are expressed in the table as multiples of the median gene expression.
GeneDescriptionMRL/lprMRL/+Fn14ko
Dnajc28DnaJ (Hsp40) homolog, subfamily C, member 28334.060.170.09
Syne2synaptic nuclear envelope 2316.020.180.09
Rpp38ribonuclease P/MRP 38 subunit (human)314.631.031.31
Cox7bcytochrome c oxidase subunit VIIb311.8623.4031.39
Kctd10potassium channel tetramerisation domain containing 10307.630.320.09
Tle6transducin-like enhancer of split 6, homolog of Drosophila E307.410.370.54
Pcnx4pecanex homolog 4306.691.261.81
Dlx5distal-less homeobox 5302.190.480.14
Clec14aC-type lectin domain family 14, member a300.990.210.14
Coq6coenzyme Q6 homolog (yeast)290.922.953.18
CraddCASP2 and RIPK1 domain containing adaptor with death domain289.920.400.60
Ppp2caprotein phosphatase 2 (formerly 2A), catalytic subunit, α isoform289.2912.0819.16
Gps1G protein pathway suppressor 1282.0919.0221.20
Zfp455zinc finger protein 455276.700.200.12
Donsondownstream neighbor of SON267.230.831.36
Tmietransmembrane inner ear265.980.430.32
Akt2thymoma viral proto-oncogene 2261.944.097.64
Deaf1deformed epidermal autoregulatory factor 1 (Drosophila)259.891.350.82
Eif2s1eukaryotic translation initiation factor 2, subunit 1 α258.350.550.84
Dnajc14DnaJ (Hsp40) homolog, subfamily C, member 14256.520.770.63
MRL/MpJ-Faslpr/lpr lupus model (MRL/lpr); MRL/MpJ background control (MRL/+); MRL/lpr Fn14 knock-out (Fn14ko)
Table 2. The 60 genes that are the largest contributors to the MRL/lpr phenotype. Positive values indicate upregulation of the gene in the MRL/lpr compared with MRL/+ and Fn14ko, and negative values indicate its downregulation. Shaded genes are ones with a similar trend in both control phenotypes, indicating that TNF-like weak inducer of apoptosis (TWEAK)/Fn14 activation plays a role in the gene’s regulation. “x”, expression ratio/fold-change; “WIR”, Weighted Individual (gene) Regulation.
Table 2. The 60 genes that are the largest contributors to the MRL/lpr phenotype. Positive values indicate upregulation of the gene in the MRL/lpr compared with MRL/+ and Fn14ko, and negative values indicate its downregulation. Shaded genes are ones with a similar trend in both control phenotypes, indicating that TNF-like weak inducer of apoptosis (TWEAK)/Fn14 activation plays a role in the gene’s regulation. “x”, expression ratio/fold-change; “WIR”, Weighted Individual (gene) Regulation.
GeneDescriptionMRL/lpr vs. MRL/+MRL/lpr vs. Fn14ko
xWIRxWIR
Nrtnneurturin−17.12−1032.44−1.11−0.24
Dnajc28DnaJ (Hsp40) homolog, subfamily C, member 281979.55203.313817.12203.39
Syne2synaptic nuclear envelope 21714.27192.303672.93192.40
Rpp38ribonuclease P/MRP 38 subunit (human)304.48190.96241.04190.68
Tle6transducin-like enhancer of split 6, homolog of Drosophila E(spl)830.65187.22572.06187.05
Kctd10potassium channel tetramerisation domain containing 10950.07187.043585.62187.28
Dlx5distal-less homeobox 5635.75183.642170.24183.99
Clec14aC-type lectin domain family 14, member a1421.76183.142104.76183.21
Cox7bcytochrome c oxidase subunit VIIb13.33178.489.93170.27
CraddCASP2 and RIPK1 domain containing adaptor with death domain719.92176.35484.24176.15
Coq6coenzyme Q6 homolog (yeast)98.49175.5091.40175.27
Ppp2caprotein phosphatase 2 (formerly 2A), catalytic subunit, α isoform23.95169.9215.10162.69
Zfp455zinc finger protein 4551373.60168.352255.80168.43
Donsondownstream neighbor of SON322.96162.56196.84162.02
Tmietransmembrane inner ear618.70161.70818.52161.81
Gps1G protein pathway suppressor 114.83160.9013.31158.67
Akt2thymoma viral proto-oncogene 264.02158.9434.28155.29
Deaf1deformed epidermal autoregulatory factor 1 (Drosophila)193.09157.27316.93157.80
Eif2s1eukaryotic translation initiation factor 2, subunit 1 α472.78157.09307.81156.79
Dnajc14DnaJ (Hsp40) homolog, subfamily C, member 14332.10155.60410.38155.75
Arrdc1arrestin domain containing 1461.59155.23406.91155.15
NkapNFKB activating protein−23.74−146.421.120.04
Sall1sal-like 1 (Drosophila)515.64142.36977.24142.58
Rnmtl1RNA methyltransferase like 1286.57141.15314.45141.22
Ybx2Y box protein 2799.94139.73707.63139.70
Arl4aADP-ribosylation factor-like 4A203.58139.32117.37138.48
Tbc1d19TBC1 domain family, member 19 (Tbc1d19), mRNA [NM_144517]259.29136.98213.37136.79
Lrrc2leucine rich repeat containing 21145.66135.041901.54135.11
Gpr183G protein-coupled receptor 1831187.49131.432459.92131.53
Abcc3ATP-binding cassette, sub-family C (CFTR/MRP), member 3988.46131.051673.97131.14
Ccdc124coiled-coil domain containing 12416.46121.3513.70118.72
Rgl1ral guanine nucleotide dissociation stimulator-like 1113.72118.88182.93119.55
Mki67antigen identified by monoclonal antibody Ki 67903.22101.46977.16101.48
Hba-a1hemoglobin α, adult chain 13.6991.051.8254.31
Hba-a2hemoglobin α, adult chain 22.5490.321.7764.14
TankTRAF family member-associated Nf-kappa B activator381.5587.10407.5587.12
Slc10a1solute carrier family 10 (sodium/bile acid cotransporter family), member 11.8283.411.8283.18
H2-Gs10MHC class I like protein GS1010.3882.279.0280.23
Rnf168ring finger protein 16887.2280.35106.7580.63
Elk1ELK1, member of ETS oncogene family1.5965.371.7575.96
Exosc2exosome component 21.7462.551.3939.39
Cxcl11chemokine (C-X-C motif) ligand 112.8261.182.6859.47
Ttrtransthyretin−2.00−52.38−1.76−40.27
Kif5akinesin family member 5A1.6952.291.3026.92
Hccsholocytochrome c synthetase3.3750.043.1448.37
Zfp628zinc finger protein 628−1.34−47.431.3628.96
Atp6v0cATPase, H+ transporting, lysosomal V0 subunit C−1.40−43.10−1.18−11.73
Glagalactosidase, α2.0441.871.6131.24
Psapprosaposin−1.66−41.72−1.16−7.03
Dccdeleted in colorectal carcinoma2.3041.293.4052.26
Oxtoxytocin−16.11−39.141.020.00
Tubb2atubulin, β 2A−1.35−38.27−1.41−36.56
Ddndendrin−2.34−38.13−1.28−3.17
Ptmsparathymosin−1.44−37.911.198.51
Acot7acyl-CoA thioesterase 7−1.71−35.06−1.15−3.35
Ptgdsprostaglandin D2 synthase (brain)1.4534.821.2311.03
Thy1thymus cell antigen 1, theta−1.32−34.561.1811.17
Sparcl1SPARC-like 11.4233.771.010.13
Pkm2pyruvate kinase, muscle−1.57−33.58−1.13−3.69
Cox8acytochrome c oxidase, subunit VIIIa1.3132.711.1917.94
Table 3. The top 15 genes in each phenotype. The numbers are the gene commanding heights (GCH) in the indicated phenotypes. Top genes have the highest GCHs.
Table 3. The top 15 genes in each phenotype. The numbers are the gene commanding heights (GCH) in the indicated phenotypes. Top genes have the highest GCHs.
GeneDescriptionMRL/+MRL/lprFn14ko
Tfb2mtranscription factor B2, mitochondrial23439
Hiatl1hippocampus abundant transcript-like 18542
Mettl11amethyltransferase like 11A7314
Gipc1GIPC PDZ domain containing family, member 16821
Homezhomeodomain leucine zipper-encoding gene6022
Syt11synaptotagmin XI5823
Eef1geukaryotic translation elongation factor 1 γ5442
Brsk1BR serine/threonine kinase 15322
Gga2golgi associated, γ adaptin ear containing, ARF binding protein 25012
St3gal5ST3 β-galactoside α-2,3-sialyltransferase 54912
Dusp22dual specificity phosphatase 224621
Avl9AVL9 homolog4223
Foxq1forkhead box Q14211
Lmnalamin A4231
Anpepalanyl (membrane) aminopeptidase4111
Xrcc4X-ray repair complementing defective repair in Chinese hamster cells 44431
Tubb4tubulin, β 42432
Arhgap23Rho GTPase activating protein 232363
Ecdecdysoneless homolog5355
Cacfd1calcium channel flower domain containing 15292
NacadNAC α domain containing10262
Pnldc1poly(A)-specific ribonuclease (PARN)-like domain containing 13232
Mrpl40mitochondrial ribosomal protein L401223
Serf2small EDRK-rich factor 26222
Rpa2replication protein A27223
Unc79unc-79 homolog4205
Pogkpogo transposable element with KRAB domain8193
Gabpb1GA repeat binding protein, β 12185
Mad2l2MAD2 mitotic arrest deficient-like 231810
Ccdc101coiled-coil domain containing 1012182
Card19caspase recruitment domain family, member 192227
Ppfia1protein tyrosine phosphatase, receptor type, f polypeptide, interacting protein (liprin), α 16526
Map4k5mitogen-activated protein kinase kinase kinase kinase 52225
Ddx59DEAD (Asp-Glu-Ala-Asp) box polypeptide 592125
Fbxo25F-box protein 252121
Rg9mtd3RNA (guanine-9-) methyltransferase domain containing 33220
IscuIscU iron-sulfur cluster scaffold homolog, nuclear gene encoding mitochondrial protein3320
Gm1440760S ribosomal protein L27a pseudogene3219
Riok2RIO kinase 23519
Slc25a47solute carrier family 25, member 473219
Mapre3microtubule-associated protein, RP/EB family, member 39219
Pomcpro-opiomelanocortin-α3518
Polr1epolymerase (RNA) I polypeptide E6318
Olfr1347olfactory receptor 13473417
Tor2atorsin family 2, member A9116
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Iacobas, D.A.; Wen, J.; Iacobas, S.; Putterman, C.; Schwartz, N. TWEAKing the Hippocampus: The Effects of TWEAK on the Genomic Fabric of the Hippocampus in a Neuropsychiatric Lupus Mouse Model. Genes 2021, 12, 1172. https://doi.org/10.3390/genes12081172

AMA Style

Iacobas DA, Wen J, Iacobas S, Putterman C, Schwartz N. TWEAKing the Hippocampus: The Effects of TWEAK on the Genomic Fabric of the Hippocampus in a Neuropsychiatric Lupus Mouse Model. Genes. 2021; 12(8):1172. https://doi.org/10.3390/genes12081172

Chicago/Turabian Style

Iacobas, Dumitru A., Jing Wen, Sanda Iacobas, Chaim Putterman, and Noa Schwartz. 2021. "TWEAKing the Hippocampus: The Effects of TWEAK on the Genomic Fabric of the Hippocampus in a Neuropsychiatric Lupus Mouse Model" Genes 12, no. 8: 1172. https://doi.org/10.3390/genes12081172

APA Style

Iacobas, D. A., Wen, J., Iacobas, S., Putterman, C., & Schwartz, N. (2021). TWEAKing the Hippocampus: The Effects of TWEAK on the Genomic Fabric of the Hippocampus in a Neuropsychiatric Lupus Mouse Model. Genes, 12(8), 1172. https://doi.org/10.3390/genes12081172

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop