RNA Sequence Analyses throughout the Course of Mouse Cardiac Laminopathy Identify Differentially Expressed Genes for Cell Cycle Control and Mitochondrial Function

Lamin A/C (LMNA) gene mutations are a known cause of familial dilated cardiomyopathy, but the precise mechanisms triggering disease progression remain unknown. We hypothesize that analysis of differentially expressed genes (DEGs) throughout the course of Lmna knockout (Lmna−/−)-induced cardiomyopathy may reveal novel Lmna-mediated alterations of signaling pathways leading to dilated cardiomyopathy. Although Lmna was the only DEG down-regulated at 1 week of age, we identified 730 and 1004 DEGs in Lmna−/− mice at 2 weeks and 1 month of age, respectively. At 2 weeks, Lmna−/− mice demonstrated both down- and up-regulation of the key genes involving cell cycle control, mitochondrial dysfunction, and oxidative phosphorylation, as well as down-regulated genes governing DNA damage repair and up-regulated genes involved in oxidative stress response, cell survival, and cardiac hypertrophy. At 1 month, the down-regulated genes included those involved in oxidative phosphorylation, mitochondrial dysfunction, nutrient metabolism, cardiac β-adrenergic signaling, action potential generation, and cell survival. We also found 96 overlapping DEGs at both ages involved in oxidative phosphorylation, mitochondrial function, and calcium signaling. Impaired oxidative phosphorylation was observed at early disease stage, even before the appearance of disease phenotypes, and worsened with disease progression, suggesting its importance in the pathogenesis and progression of LMNA cardiomyopathy. Reduction of oxidative stress might therefore prevent or delay the development from Lmna mutation to LMNA cardiomyopathy.


Results
Homozygous Lmna −/− mice exhibit severe growth retardation and cardiac dysfunction at 1 month of age. Homozygous mice harboring the targeted Lmna mutation, Lmna tm1Stw , were previously shown to have progressive muscular dystrophy and dilated cardiomyopathy [22][23][24] . We observed a similar disease progression from 1 week to 1 month. But unlike the former study, which showed a significant decrease in mouse weight in 2-week-old Lmna −/− mice 15 , we did not identify a significant difference in growth between wild type (WT) and Lmna −/− mice at 1 week and 2 weeks of age (Fig. 1a). However, by 1 month of age Lmna −/− mice exhibited severe growth retardation and cardiac dysfunction (Fig. 1a,b). This cardiac dysfunction was accompanied by increased myocardial fibrosis in Lmna −/− mice at 1 month of age (Fig. 1c). No appreciable differences in myocardial fibrosis were observed between Lmna −/− and WT mice at either 1 week or 2 weeks of age.
Early differentially expressed genes (DEGs) and related pathways in Lmna −/− mice. We first compared the gene expression profiles of WT and Lmna −/− mice at 1 week and 2 weeks of age. Although Lmna was the only gene with significantly reduced expression in 1-week-old Lmna −/− mice, we identified 730 DEGs (428 genes up-regulated, 293 genes down-regulated, 9 genes unmapped) in 2-week-old Lmna −/− mice. Based on these genes, 58 related canonical pathways (with p < 0.05) were identified through IPA pathway analysis. Figure 2a and Supplementary Table S1 list the six pathways with the highest p values. Table 1 lists representative key DEGs in the related pathways/functions among 2-week Lmna −/− mice. As shown in Table 1, Lmna −/− mice at 2 weeks of age demonstrated both down-and up-regulation of the key genes involving cell cycle regulation (e.g. chromosomal replication down-regulated, counteractive effects on G1/S and G2/M transition, M phase progression repressed), mitochondrial dysfunction, oxidative phosphorylation, and apoptosis. Meanwhile, Lmna −/− mice showed down-regulation in key genes in DNA damage response/repair pathways as well as up-regulation in genes involved in the oxidative stress response, acute phase response, cell survival/growth, cardiac hypertrophy, glycolysis, and triglyceride hydrolysis. We identified the up-regulation of FOXO3 (listed in Table 1 apoptosis pathway), which was also found in a previous study in 2-week Lmna −/− mice by Marian et al. 15 . However, instead of activation of the pathway in 2-week Lmna D300N mice shown by Marian et al. 16 , we identified overall down-regulation of DNA damage response pathway (z score = −1.265). But we did observe the overall up-regulation of TP53 pathway (z score = 0.302) due to up-regulation of TP53 effector (PERP) and other TP53 target genes, which was again similar to the findings from Marian et al. in 2-week Lmna D300N mice 16 .
Taken together, 2-week-old Lmna −/− mice demonstrated the activation of cell cycle arrest and apoptosis pathways and a failure of DNA repair in response to the targeted Lmna knockout mutation. Although the trend of cell cycle arrest and apoptosis seemed clear, some counteractive activation pathways were seen as well. For example, multiple cell survival factors were up-regulated to keep cells growing. In addition, while mitochondrial dysfunction (down-regulated genes for Complex III, IV, V) might compromise ATP generation and cause oxidative stress, other enzyme-coding genes for Complex I, III, IV, V and glycolysis were up-regulated to overcome ATP shortage, and the NRF2-mediated oxidative stress response pathway was dramatically activated to reduce oxidative damage. These counteractive effects might explain the lack of observed growth retardation at 2 weeks. In the meantime, up-regulation of the genes in cardiac hypertrophy, myofibroblast activation/fibrosis and acute phase response signaling might contribute to the early pathological changes before LMNA cardiomyopathy occurs.
Mitochondrial dysfunction and diminished nutrient metabolism associated with extensive late-stage DEGs in 1-month-old Lmna −/− mice. We identified 1004 DEGs between the gene expression profiles of WT and Lmna −/− mice at 1 month of age, with 699 genes had decreased expression levels while 290 genes had increased expression levels (15 genes unmapped) in Lmna −/− mice. Pathway analysis via IPA identified 66 canonical pathways (p < 0.05) based on these genes. Figure 2b and supplementary Table S2 list the top six pathways. Table 2 lists the examples of key DEGs in the related pathways/functions among 1-month-old Lmna −/− mice. Most impressively, these mice demonstrated extensive down-regulation of key genes associated with mitochondrial dysfunction, oxidative phosphorylation, a wide range of metabolic pathways (tricarboxylic acid cycle, glycolysis, glycogen synthesis, glycogenolysis, amino acid degradation, fatty acid β-oxidation, ketogenesis, and ketolysis), cardiac β-adrenergic signaling, G protein β γ signaling, caveolar-mediated endocytosis, and cardiac voltage-gated channels ( Table 2). Such alterations in pathways of metabolism and energy status in Lmna −/− mice may explain their remarkable growth retardation. Down-regulation of multiple voltage-gated channels that have been reported to be related to the generation of arrhythmias 25 -KCNA5, KCND2, KCNE1, SCN4A, SCN4B (listed in Table 2) -might be the underlying basis for the increased risk of arrhythmogenic events in LMNA cardiomyopathy. Meanwhile, we observed down-and up-regulated genes among these mice involved in cell survival signaling, cardiac hypertrophy, calcium signaling and autophagy (Table 2).
Several up-regulated genes in 1-month Lmna −/− mice are associated with sarcomere structure and ERK1/2 pathway. Compared with previous microarray studies by Worman et al. on mouse LMNA cardiomyopathy from Lmna H222P/H222P mutation 7,26 , we identified up-regulation of a novel gene associated with sarcomere structure, myomesin 2 (MYOM2). MYOM2 is a protein coding gene for M-band in sarcomere that plays an important role in maintaining sarcomere structure 27,28 . In addition, we identified up-regulation of dual specificity phosphatase 4 (DUSP4) and DUSP5 in 1-month Lmna −/− mice, while only DUSP4 was up-regulated in mice with Lmna H222P/ H222P -induced cardiomyopathy based on the study by Worman et al. 17 . DUSP4 and DUSP5 can be transcriptionally induced by ERK1/2, which suggests the possibility of ERK1/2 pathway activation in LMNA cardiomyopathy progression. In addition, we found some compensatory up-regulated genes in Lmna −/− mice, e.g. AMP-Activated Protein Kinase Subunit Gamma-2 (PRKAG2) was up-regulated by higher cellular AMP/ATP ratio to inhibit key enzymes of ATP consuming pathways and induces ATP generation, and cAMP-dependent protein kinase type I alpha regulatory subunit (PRKAR1A) was up-regulated by cAMP accumulation for ATP generation and cell survival.   www.nature.com/scientificreports www.nature.com/scientificreports/ down-regulated) at the two time points while 32 genes are regulated in opposite directions (30 genes up-regulated at 2 weeks but down-regulated at 1 month and 2 genes down-regulated at 2 weeks but up-regulated at 1 month). Pathway analysis identified 23 canonical pathways (p < 0.05) based on the overlapping DEGs. Figure 2c and supplementary Table S3 listed the top six canonical pathways. Table 3 lists key overlapping DEGs, which were involved in oxidative phosphorylation, mitochondrial function, calcium signaling, G protein β γ signaling, and caveolar-mediated endocytosis. Although expression changes of some listed key genes could cause functional impairment at both time points, compensatory expression changes of other genes could prevent the occurrence of actual functional impairment as in 2-week-old Lmna −/− mice. For example, the coding gene for Complex V enzyme MT-ATP6 was down-regulated in these mice (Table 3), which could affect ATP generation in mitochondria. However, the coding genes for Complex I/IV enzymes NDUFA13, NDUFS7, NDUFV1 and COX6A2 are simultaneously up-regulated (Table 3), producing more substrates for MT-ATP6 and protecting the mice from ATP shortage.
Other important overlapping DEGs include those involved in glycogen degradation (acid α-glucosidase, phosphofructokinase -PFKM), glycolysis (PFKM), and lactose degradation (prosaposin). All these were up-regulated in 2-week-old and down-regulated in 1-month-old mice, suggesting more glucose generation and usage in 2-week-old Lmna −/− mice and a dysfunction of glucose metabolism in 1-month-old mice. In addition, heat shock protein family B member 1 (HSPB1), nicotinamide riboside kinase (NMRK2) and dual specificity phosphatase 5 (DUSP5) were up-regulated at both time points. HSPB1 is important for stress resistance and actin organization while NMRK2 codes a key enzyme for NAD + synthesis from nicotinamide riboside (NAD + is one of the most important co-enzymes for redox reactions), and DUSP5, as a member of the dual specificity protein phosphatase subfamily, negatively regulates members of the mitogen-activated protein kinase (MAPK) superfamily (ERK1/2, stress-activated protein kinase/c-Jun N-terminal kinase -SAPK/JNK). DUSP1 and DUSP5 were up-regulated in To validate our RNA-sequencing results, we performed RT-qPCR and western blot for representative genes associated with the top 6 pathways (for 2 weeks, listed in Fig. 2a) as well as oxidation phosphorylation, mitochondrial function, fatty acid metabolism, sarcomere structure, and cardiac development for 1 month. RT-qPCR analysis confirmed decreased expression of Aurora Kinase A (AURKA), breast cancer type 1 susceptibility protein (BRCA1), cyclin dependent kinase 1 (CDK1), checkpoint kinase 1 (CHEK1), cyclin B1, minichromosome main-

Discussion
Our study is the first to utilize RNA sequencing and pathway analysis to set up detailed gene expression and pathway profiles over the course of Lmna −/− induced LMNA cardiomyopathy, which may lead to further mechanistic studies of the pathogenesis and progression of LMNA cardiomyopathy. Based on an average lifespan of 5 weeks for Lmna −/− mice, we selected three different time points that correspond to the early (1 week), middle (2 weeks), and late (1 month) stages of the disease to investigate the detailed molecular mechanisms over time in hearts isolated from WT and Lmna −/− mice. The significantly increased number of genes with expression changes between 1 and 2 weeks (from 1 gene to 730 genes) and between 2 weeks to 1 month of age (from 730 genes to 1004 genes) was associated with the phenotypic observation with significant cardiac dysfunction and fibrosis.
The major findings of this study are the gene expression and related pathophysiological changes at different stages of LMNA cardiomyopathy. We have further identified how these changes evolved over time with disease progression. More importantly, we have identified 96 overlapping DEGs in 2-week-old and 1-month-old Lmna −/− mice, which are mainly involved in oxidative phosphorylation, mitochondrial dysfunction, calcium signaling, and G protein beta gamma signaling. These findings have provided insights on the pathogenesis of LMNA cardiomyopathy due to loss of LMNA function.
We found the DEGs in 2-week-old Lmna −/− mice were mainly involved in cell cycle regulation, DNA damage response, mitochondrial function, oxidative phosphorylation, oxidative stress response, and counteractive  www.nature.com/scientificreports www.nature.com/scientificreports/ apoptosis/survival signaling. Although these mice seemed to fail in repairing their DNA damage, i.e. Lmna mutation (all related DEGs down-regulated) with down-regulated chromosomal replication and a trend of delayed cell cycle progression, we still observed some down-regulated genes meant to promote cell cycle progression. Similarly, up-regulated gene expression for mitochondrial complex I, III, IV and V enzymes might offset the consequences caused by down-regulated genes for complex III, IV, V enzymes for ATP synthesis from mitochondria. While oxidative stress may exist because of impaired oxidative phosphorylation, the oxidative stress response pathway was activated (all related DEGs up-regulated) to reduce oxidative injury. A similar pattern was observed with apoptosis and survival factors. There were up-regulated genes on both sides. These compensatory mechanisms might explain why the 2-week-old mice didn't show any signs of LMNA cardiomyopathy or growth retardation yet.
For 1-month-old Lmna −/− mice, we observed significant growth and cardiac phenotypes associated with extensive down-regulation of genes involving oxidation phosphorylation, mitochondrial function, almost all nutrient metabolism (e.g. TCA cycle, amino acid degradation, and fatty acid β-oxidation), cardiac β-adrenergic signaling, G protein beta gamma signaling. Overwhelming oxidative stress and ATP deficiency caused by seriously impaired oxidation phosphorylation and mitochondrial dysfunction, together with extensively reduced nutrient metabolism might explain the remarkable growth retardation of Lmna −/− mice at this stage. Up-regulated cardiac β-adrenergic signaling (z score = 1) and dysregulated G protein beta gamma signaling (z score =0) in the heart could also contribute to cardiac dysfunction we observed. Our study provided a unique model that could identify comprehensive and accurate gene expression levels using RNA sequence analysis to further understand signaling pathways involved in the progressive dilated cardiomyopathy process from Lmna mutation.
Our study is the first to report that mitochondrial dysfunction and impaired oxidative phosphorylation can occur in Lmna −/− mice as early as 2 weeks of age, even before an apparent phenotype of LMNA cardiomyopathy. At late disease stage (1 month of age), when the phenotypes are dramatic, these impairments are more prominent. 47.5% genes involving oxidative phosphorylation and 35.2% genes involving mitochondrial dysfunction are down-regulated. Mitochondrial dysfunction has been described in dilated or hypertrophic cardiomyopathies and heart failure but its exact mechanism is unclear [29][30][31][32] . Mitochondrial dysfunction further promotes increased reactive oxygen species 29,33 . Mitochondrial dysfunction and increased oxidative stress have been noted in progeria laminopathies [34][35][36][37] . A specific gene targeted therapy or drug therapy to reduce oxidative stress has shown improvement of progeria phenotype 34,35 . It was previously shown that there was increased oxidative stress in heart failure and efforts had been made to target the oxidative stress pathway to prevent heart failure progression 38,39 . Our study suggested the possibility that reduction of oxidative stress might prevent or delay the development from Lmna mutation to LMNA cardiomyopathy.
Increased reactive oxidative stress and mitochondrial dysfunction can also affect ERK1/2 as well as mTOR pathways 40,41 . The over-activation of ERK1/2 and mTOR pathways has been well described in LMNA cardiomyopathy 8,9,11,13 . Our RNA sequence analysis shows increased gene expression level of dual specificity phosphatase 1 -DUSP1 (dephosphorylates and inactivates ERK2) and DUSP5 (inactivates ERK1) at 2 weeks of age as well as up-regulated DUSP4 (inactivates ERK1, ERK2 and JNK), DUSP5 and down-regulated DUSP7, DUSP8 (inactivates SAPK/JNK and p38) and DUSP18, suggesting over-activation of ERK1/2 signaling pathway. We also find that mitogen-activated protein kinase kinase kinase 6 -MAP3K6 (upstream kinase of MAPK pathways) and mitogen-activated protein kinase-activated protein kinase 2 -MAPKAPK2 (regulated through direct phosphorylation by p38 MAPK) are up-regulated in 2-week-old Lmna −/− mice. Meanwhile, we observe multiple up-regulated genes involving insulin like growth factor 1 (IGF1) signaling including Akt2 in 2-week-old Lmna −/− mice, which may promote mTOR activation. In addition, we identify a novel coding gene for nicotinamide riboside kinase  www.nature.com/scientificreports www.nature.com/scientificreports/ (NMRK2) that is significantly up-regulated in Lmna −/− mice. NMRK2 has been shown to play a key role in myogenesis especially through its interactions with integrins thereby activating the integrin-mediated signaling pathways 42,43 . ERK1/2 pathway is a well-known downstream effector of integrin-mediated signaling pathway and our finding provides further insights of ERK1/2 signaling mechanisms in LMNA cardiomyopathy.
Besides identifying the importance of oxidative stress in the pathogenesis of LMNA cardiomyopathy, our study is the first to report the upregulation of counteractive factors and signaling pathways for both cell apoptosis and survival in 2-week Lmna −/− mice before the appearance of disease phenotypes. The upregulated pro-apoptotic factors include forkhead box O3 (FOXO3), BCL2 like 11 (BCL2L11, BIM), BCL2 interacting protein 3 (BNIP3), TP53 apoptosis effector (PERP), et al. (Table 1). Some of these factors are targets of FOXO3 and TP53, suggesting the upregulation of FOXO3 and TP53 pathways. Upregulated TP53 pathway (p = 0.00233, z score 0.302) is identified by IPA pathway analysis based on 2-week DEGs. The upregulated survival factors include protein kinase B beta (Akt2), pyruvate dehydrogenase kinase 2 (PDK2), BCL2 like 1 (BCL2L1, BCL-XL), insulin receptor substrate 2 (IRS2), insulin like growth factor binding protein 3/5/6 (IGFBP3/5/6), et al. (Table 1). We observed significantly upregulated IGF-1 pathway (p = 0.00545, z score 2.236), NRF2-mediated oxidative stress response (p = 1.59E-6, z score 3.207) and acute phase response signaling (p = 9.42E-4, z score 3.317), which may work together to promote cell survival and growth. It is reasonable to think these may be compensatory mechanisms under the background of DNA damage (deletion mutation), failure of DNA damage repair and counteractive factors on cell cycle progression in order to keep 2-week-old Lmna −/− mice in temporarily similar overall growth as that of their WT littermates. Interventions which inhibit cell apoptosis or enhance survival pathways may delay the appearance of disease phenotypes and improve the lifespan of Lmna −/− mice.  15 . They also found suppression of FOXO1/3 activation by shRNA down-regulated the target gene transcription, improved apoptosis and prolonged survival of 2-week mice by 2-fold 15 . Conformed to this study, our study also finds the up-regulation of FOXO3 and its target genes (BCL2 L11/ BIM, BCL6, BNIP3 for apoptosis/ autophagy, P21CIP1, PLK3, down-regulated PLK1/4 and cyclin B1/2 for cell cycle arrest). We don't detect the up-regulation of the target genes for DNA repair and reactive oxygen species detoxification. Although FOXO1 gene expression is not found up-regulated in our analysis, we do identify up-regulated expression of its target genes as well. No up-regulation of FOXO4 and its target genes are found in 2-week Lmna −/− mice. But we observe mostly down-regulated target genes for FOXO 1/3/4 in 1-month Lmna −/− mice although no gene expression level changes for these three FOXO factors.
FOXO3/1 proteins can be activated by DNA damage or oxidative stress 44 . But when the intensity of these stimuli increases beyond a certain threshold, they activate pro-apoptotic genes and genes for cell cycle arrest. Meanwhile, Akt is a major upstream regulator of FOXO factors [44][45][46] . Our study shows up-regulation of Akt2 (phosphorylates and inactivates FOXO) through IGF1 pathway, suggesting over-activation of FOXO3/1. In addition, over-activation of ERK1/2 (phosphorylates and inactivates FOXO3) in these mice as discussed above supports the existence of over-activation of FOXO3 47 . Combined with downstream and upstream changes, our study provides some potential mechanisms on FOXO3/1 activation in 2-week Lmna −/− mice. www.nature.com/scientificreports www.nature.com/scientificreports/ The activation of E2F transcription factor/DNA damage response/TP53 pathway has been found in a 2-week Lmna D300N mouse model 16 . In our study, E2F is predicted to be inhibited (p = 1.30E-8, z score −2.798) in 2-week Lmna −/− mice because all of its 10 target genes (among DEGs) are inhibited. For DNA damage response, all the related DEGs in 2-week Lmna −/− mice are down-regulated. However, our study shows TP53 is predicted to be activated (p = 5.88E-24, z score 3.977) with 46 target genes up-regulated in 2-week Lmna −/− mice. The data differences between our study and the previous study may be due to the difference of mouse models and different lifespan: Lmna D300N mice at 30 days vs. Lmna −/− mice at 42 days. At 2 weeks, these two kinds of mice may be at different disease stages, and their data may not be comparable.
In summary, our study identified failure of DNA repair, dysregulated cell cycle progression and mitochondrial dysfunction/impaired oxidative phosphorylation associated with loss of LMNA, thereby providing mechanistic insights to LMNA cardiomyopathy and heart failure. Compared with previous gene profiling using microarray in Lmna H222P/H222P -induced cardiomyopathy (at 10 weeks of age) 7,26 and RNA sequencing study in 2-week Lmna −/− mice 15 and Lmna D300N mice 16 , our study not only confirmed the previously known genes and pathways related to LMNA cardiomyopathy but also identified novel genes and pathways that could contribute to the pathogenesis and progression of LMNA cardiomyopathy. It will be interesting to further investigate whether we can apply the same analysis to peripheral samples for its better practicality and feasibility. This application will provide us with detailed information at genomic level which can guide us to better understand this devastating disease in human subjects and develop advanced treatments including gene therapy.

Methods
Animals. Animal use and the study protocol were approved by the Institutional Animal Care and Use Committee at the Cleveland Clinic in accordance with the United States Public Health Service Policy on the Humane Care and Use of Animals, and the NIH Guide for the Care and Use of Laboratory Animals. Heterogeneous C57BL/6 mice harboring the targeted Lmna mutation, Lmna tm1Stw , were purchased from the Jackson Laboratory (stock# 009125). Heterozygotes and littermate wild type (WT) C57BL/6 mice were inbred and LMNA gene mutation was confirmed by genotyping 22,23,48 . 8 WT and 8 Lmna −/− mice were used for each time point at 1 week, 2 weeks and 1 month of age. Hearts were isolated from WT and homozygous Lmna −/− mice at each time point. Hearts were rinsed with cold saline to remove blood, snap frozen in liquid nitrogen, and stored at −80 °C for future use or fixed in 10% formalin for histology examination. transthoracic echocardiography. Mice at 1 month of age were anesthetized with isoflurane inhalation and placed on a table. Echocardiography was performed using a GE Vivid 7 ultrasound with a 14 MHz transducer (GE Healthcare, Chicago, IL). Left ventricular shortening fraction was measured in 2D mode.
Histological examination. Formalin-fixed, paraffin-embedded left ventricular (LV) sections were prepared and stained with Picrosirus Red. Five to six regions were randomly selected from each section and fibrosis areas were analyzed by an independent researcher, viewing in the whole heart longitudinal and cross sections with Image Pro Plus v7.0 software (Media Cybernetics Inc, Rockville, MD). Fibrosis area was presented as ratio of fibrotic areas to the whole myocardium areas 49 .
RnA isolation and RnA sequence analysis. The hearts were harvested at three time points: 1 week, 2 weeks and 1 month of age, as postnatal growth of the homozygous mice is still not retarded at one week old of age but is severely retarded at 1 month of age. The entire heart tissue was used for RNA extraction to avoid variance caused by sampling different part of the heart tissue. Three hearts from WT and Lmna −/− at each time point (total 18 hearts) were used for this experiment. Total RNA was extracted and purified using RNeasy Fibrous Tissue Mini Kit (QIAGEN, Germantown, MD). RNA concentrations were measured using NanoDrop Spectrophotometers. Quality of the RNA samples was analyzed using the Agilent Bioanalyzer and all samples submitted for sequencing had a RNA Integrity number (RIN) > 8. Truseq Stranded Total RNA -RiboZero library preparation and RNA sequencing were performed by the Cleveland Clinic Lerner Research Institute Genomic Core. RNA sequencing was performed using the Illumina SBS v2 chemistry with 2 × 100 bp pair-end reads. The RNA-Sequencing data were aligned using STAR (v2.5.2b) program with default parameters 51 . Mapped read pairs were assigned to genes by collapsing all transcripts into a single gene model and then counting the number of reads that fully overlap the resulting exons. Reads that mapped to multiple locations were only counted once and those mapping to Scientific RepoRtS | (2020) 10:6632 | https://doi.org/10.1038/s41598-020-63563-x www.nature.com/scientificreports www.nature.com/scientificreports/ ambiguous regions were excluded. Reads uniquely mapped were considered for further analysis. DEseq. 2 (v3.8), a program using count-based matrices to identify differentially expressed genes (DEGs), was used to determine gene signatures 51 . To improve the reliability and accuracy of differential expression analysis, only genes with raw counts >5 in all individually sequenced samples were examined and compared between WT vs. Lmna −/− at 1 week of age, at 2 weeks of age and at 1 month of age. For differentially expressed genes, INGENUITY Pathway Analyses (IPA, Qiagen, Hilden Germany) were performed, including enrichment analysis for canonical pathway, disease and biological function, toxicology function.
Validation of the RnA-sequencing data. The most changed genes in the interesting signaling pathways were further confirmed using RT-qPCR and Western blot. RT-qPCR was performed using the TaqMan protocol in StepOnePlus Real-Time PCR System (Applied Biosystems, Foster City, CA). The TaqMan gene expression assays for the following genes were purchased from Applied Biosystems: AURKA (Mm01248177_m1), BRCA1 Statistical analysis. Data are expressed as the mean ± standard deviation (SD). Graphpad (Prism Software) was used for statistical analysis. Statistical significance was determined by 2-tailed Student's t-test with a value of p < 0.05 considered significant.

Data availability
The RNA sequencing data generated and analyzed during this study have been deposited in NCBI's Gene Expression Omnibus (GEO) and are accessible through GEO series accession number GSE133693 (https://www. ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE133693) and token number ahqvgqmypbexjgn.