Impact Factor 4.151
2017 JCR, Clarivate Analytics 2018

Frontiers journals are at the top of citation and impact metrics

Original Research ARTICLE

Front. Genet., 07 January 2019 | https://doi.org/10.3389/fgene.2018.00684

Time- and Ventricular-Specific Expression Profiles of Genes Encoding Z-Disk Proteins in Pressure Overload Model of Left Ventricular Hypertrophy

Anastasia Knyazeva1*, Alexander Krutikov1, Alexey Golovkin1, Alexander Mishanin1, Georgii Pavlov2, Natalia Smolina1,3, Anastasia Hushkina1, Thomas Sejersen4, Gunnar Sjoberg4, Mikhail Galagudza1 and Anna Kostareva1,3
  • 1Almazov National Medical Research Centre, Saint Petersburg, Russia
  • 2Saint Petersburg State Academy of Veterinary Medicine, Saint Petersburg, Russia
  • 3Information Technologies and Programming Faculty, ITMO University, Saint Petersburg, Russia
  • 4Department of Women’s and Children’s Health, Center for Molecular Medicine, Karolinska Institute, Stockholm, Sweden

Mechanotransduction is an essential mechanism of transforming external mechanical stimulus to biochemical response. In cardiomyocytes mechanotransduction plays an important role in contraction, stretch sensing and homeostasis regulation. One of the major mechanosensitive area in cardiomyocytes, the Z-disk, consists of numbers of structural and signaling proteins, that may undergo conformational or gene expression changes under pathological stress conditions. In present study we examined a rat model of pressure overload cardiac hypertrophy validated by echocardiographic and histopathological examinations. We revealed, that during hypertrophy progression expression of several genes encoding Z-disk proteins (Actn2, Ldb3, Cmya5, Nebl) is different at early and late points of cardiac remodeling. Moreover, expression patterns of several genes are opposite in myocardium of overloaded left ventricle and hemodynamically unaffected right ventricle, and expression profiles in interventricular septum are more similar to right ventricle. Additionally, we revealed inconsistencies between mRNA and protein level changes of Actn2, one of the major structural Z-disk element. All these findings point out, that investigated Z-disk proteins participate in pathological stress adaptation through undergoing the gene expression changes, and suggest the novel important role of hypertrophic response modulation during different stages of cardiac remodeling.

Introduction

The mechanism of transforming external stretch signals to a biochemical response in cytoplasmic elements and nucleus is one of the important components of cell adaptation and homeostasis control. This process, called mechanotransduction, plays an important role in different type of mechanosensing cells, and its dysregulation may consequently lead to progression of various diseases (Jaalouk and Lammerding, 2009). In cardiac and skeletal muscles mechanotransduction is used to adapt to external stimuli by changes in contractile and metabolic properties of muscle cells. In cardiomyocytes mechanosensing and mechanotransducing structures are concentrated in specific cell areas, and vast number of proteins and protein complexes involved in mechanical signal transmission are localized at sarcolemma, intercalated disks and sarcomere (Epstein and Davis, 2003; Lyon et al., 2015). The latter consists of myofilaments and large amount of interacting proteins that provide together the myofilament stability and integrity. One of the major part of sarcomeric architecture is the Z-disk – a relatively small macromolecular structure, containing the huge network of proteins which is proposed to serve as one of the key mechanosensing areas in cardiomyocytes. These proteins form the complex framework structure, which determinates contractile properties in normal conditions and under strain (Pyle and Solaro, 2004; Gautel and Djinovic-Carugo, 2016; Tabish et al., 2017). Except for the structural functions, Z-disk proteins are involved in transmembrane signaling, and activation of several stress-signal pathways in cytoplasm via interactions with kinases, phosphatases and through intracellular calcium signaling (Hoshijima, 2006; Luther, 2009). Mutations in genes encoding Z-disk components may lead to hereditary cardiac and muscle disorders such as cardiomyopathies and myopathies (Knöll et al., 2011; Kostareva et al., 2016; Schänzer et al., 2018). As a consequence, one of proposed molecular mechanisms for cardiomyopathies suggests that the impaired sarcomere stretch response is due to mutations in Z-band components (Kimura, 2016; Dal Ferro et al., 2017). Additionally, several newly described by genome-wide association studies loci related to dilated cardiomyopathy (DCM), eccentric heart remodeling and heart failure include Z-line genes (Villard et al., 2011; Garnier et al., 2015). However, exact mechanism of gene expression regulation under hemodynamic load and the detailed role of Z-line-associated proteins in this process are still not completely understood. In the present study, we used an in vivo model of myocardial hypertrophy caused by pressure overload to analyze time-dependent dynamics of Z-line-associated gene expression in different myocardial chambers (left ventricle, right ventricle, and interventricular septum). We demonstrated that several genes, encoding Z-disk proteins, have distinct mRNA expression patterns between free left ventricle wall and interventricular septum with the latter having expression profiles more similar to hemodynamically uncompromised right ventricle. The observed changes did not directly lead to modulation of protein expression, pointing on protein degradation pathways as an important regulatory mechanism in hemodynamic overload conditions.

Materials and Methods

Ethics Statement

All procedures were performed in accordance with the Guide for the Care and Use of Laboratory Animals published by the National Institute of Health and approved by the Local Ethics committees at Almazov National Medical Research Centre (Reference No. 16-3).

Experimental Procedure of Pressure-Overload Hypertrophy

Eight-weeks-old male specific-pathogen-free Wistar rats (RRID:RGD_13508588) underwent banding of the aortic arch to perform pressure overload (Cho et al., 2014). During surgical procedures anesthesia was maintained with 1.8–2.2% isoflurane (“Foran,” Abbott Laboratories). Surgical access was performed via a left anterolateral thoracotomy in 4th intercostal space with partial resection of the 3rd rib in the area from parasternal to the midclavicular line. For the access to the aortic arch a small hole in the left lobe of the thymus formed by means of the blunt tissue separation was used. Banding was performed with 6-0 Prolene (Ethicon) ligation around the aortic arch and a 22 G blunted needle. In sham-operated groups similar surgical procedures were performed, only without aortic ligation. Due to severe pressure overload conditions post-operative mortality was 3.4%. Animals were divided into groups according to model duration – 1 week (n = 16), 2 weeks (n = 12), 8 weeks (n = 14), 10 weeks (n = 7). Intact (n = 9) and sham-operated animals [1 week (n = 8), 2 weeks (n = 3), and 8 weeks (n = 4)] were examined as control. During all period animals received food and water ad libitum and care was performed according to Good Laboratory Practice (GLP) standards.

After corresponding model duration animals were sacrificed under 2.2% isoflurane by removing hearts. Hearts was dissected from atria and vessels, heart apex was mounted into Frozen Section Medium (Richard-Allan Scientific Neg-5, Thermo Fisher Scientific), frozen by submerging in liquid nitrogen and stored at -80°C. Remaining heart tissue was divided into left and right ventricles, interventricular septum, and the parts were immediately frozen separately at -80°C until use.

Echocardiography

During echocardiography rats were maintained with 1.5% isoflurane anesthesia. Ultrasound study was performed using Vevo-2100 (VisualSonics, Inc.) with scanning frequency 21 MHz and frame rate 100–120/s in 2D-mode. All measurements and calculations were performed according to recommendations from American Society of Echocardiography and the European Association of Cardiovascular Imaging (Lang et al., 2015). 2D- and M-images were obtained in longitudinal and short axes of the left ventricle (LV). M-mode cursor was positioned under control of 2D-imaging perpendicularly to left ventricle posterior wall (LVPW) and interventricular septum (IVS). Obtained images were analyzed offline. Left ventricle internal dimension (LVID), interventricular septal thickness at diastole (IVSd), left ventricle posterior wall thickness at diastole (LVPWd) were measured in two planes (short and long axes of left ventricle) in three cardiac cycles. Left ventricular mass (LVM) was calculated using mean value by Devereux and Reichek method (Devereux et al., 1986):

LVM(in mg)=0,8×1,04×[(IVSd+LVPWd+LVIDd)3LVIDd3]+0,6;

and left ventricle fractional shortening was calculated as

FS%=[(LVIDdLVIDd)/LVIDd]×100,

where d and s indicate diastole and systole.

Left ventricle ejection fraction was calculated as

LVEF(%)=[(EDVESV)/EDV]×100,

where EDV is end-diastolic volume, ESV is end-systolic volume.

Morphological Examination

12-μm thick cryostat sections of heart apices were prepared using cryotome The Tissue-Tek Cryo3 Flex (Sakura). Before staining sections were thawed and air dried for 30 min at room temperature. After blocking in 15% fetal calf serum (FCS) for 30 min at room temperature, samples were incubated overnight at +4°C with mouse monoclonal beta-sarcoglycan antibodies (Leica Biosystems, NCL-L-b-SARC). Signal was detected with DAB peroxidase substrate reaction using Novolink Polymer Detection System (Leica Biosystems) according to manufacturer’s recommendations. Sections were mounted with Richard-Allan Scientific Mounting Medium (Thermo Fisher Scientific). Imaging were acquired using light microscope Axio Observer D1 and ZEN software (Zeiss) with final magnification × 630. For analysis only cells with elliptical shape, distinct boundaries and approximately centrally located nucleus were chosen. Size of cardiomyocytes were quantified measuring minimal diameter of elliptical shaped cells in transverse section. Diameters were evaluated for 100 ± 20 cells from at least five microphotographs of an individual heart using NIH ImageJ software. Analysis was performed in a blind manner.

RNA Isolation and Real-Time PCR

For total RNA isolation pieces of frozen samples were homogenized in TissueLyzer (QIAGEN) for 8 min at 50 Hz in Extract RNA reagent (Evrogen). Next purification steps were performed following instructions for Extract RNA reagent. Quality and purity of RNA were validated using NanoDrop 3300 SpectroPhotometer (Thermo Fisher Scientific) and electrophoresis in 1% agarose gel. RNA samples were treated with DNAseI (Thermo Fisher Scientific). cDNA template was synthesized using Random (dN)10-primer (Evrogen) and MMLV RT kit (Evrogen). Quantitative Real-Time PCR (qRT-PCR) was performed on cDNA template using commercial TaqMan Gene Expression Assays for Nppa, Fhl1, Fhl2, Actn2, Synpo2, Myoz2, Nebl, Cmya5, Ldb3, Csrp3, Ilk, and Actb genes (Applied Biosystems). Expression of Hprt1 was evaluated by qRT-PCR using oligonucleotide gene-specific primer pair. The list of TaqMan assays and primers is provided in Table 1. mRNA relative expression was quantified applying ΔΔCt-method using Actb and Hprt1 as housekeeping control (Livak and Schmittgen, 2001).

TABLE 1
www.frontiersin.org

Table 1. Oligonucleotide primers and Assay ID.

Immunoblotting

Frozen heart tissue samples were homogenized in TissueLyzer for 15 min at 50 Hz in lysis buffer [50 mM Tris-HCl pH 8.0, 1% Triton X-100, 0.1% sodium dodecyl sulfate (SDS), 150 mM sodium chloride] supplemented with protease inhibitors cocktail (Sigma-Aldrich), and soluble fraction of lysates was served. Protein concentration was measured using DC Protein Assay (Bio-Rad). Lysates were run in SDS- polyacrylamide gel. After transfer membranes were primary stained with 0.1% Ponceau S (Sigma-Aldrich) and further incubated with α-actinin (Santa-Cruz, Cat#sc-166524) or p62/SQSTM1 (Sigma-Aldrich, Cat# P0067) primary antibodies overnight. After incubation with spice-specific secondary antibodies (Immun-Star Goat Anti-Mouse (GAM) or Goat Anti-Rabbit (GAR) HRP Conjugate, Bio-Rad) and development with SuperSignal West Femto substrate (Thermo Fisher Scientific) chemiluminescence was detected using Fusion Fix (Vilbert Lourmat). To avoid the bias linked to uneven increase in protein content under myocyte hypertrophy, the normalization was performed using not a single house-keeping gene product such as GAPDH or β-actin, but a total protein level detected by Ponceau S staining and normalized as relative optical densities. Analysis was performed using ImageJ software.

Statistical Analysis

Statistical analysis was performed in GraphPad Prism software. Normal distribution was tested using Shapiro–Wilk test. For paired comparison of normally distributed datasets Student’s t-test was used; for datasets that do not have normal distribution Mann–Whitney U-test was used. Relationship between LVM index and Nppa mRNA level in left ventricle was accessed by calculation of Pearson coefficient. Significant differences were considered if P < 0.05. All datasets are presented as dot-plot with mean ± standard deviation.

Results

Dynamics of Myocardial Hypertrophy in Aortic Constriction Model

A rat model of cardiac hypertrophy was established by controlled constriction of the aortic arch. Myocardial weighing revealed significant increase in heart weight-to-body weight ratio by 10 weeks after experimental procedure compared to intact animals (3.48 ± 0.78 and 2.34 ± 0.20 correspondingly) (Table 2). Increase in LVM index along with no significant LVIDd increase confirmed the development of concentric hypertrophy by week 10 of aortic banding. Progression of cardiac remodeling was supported by strong tendency of LVPW increase after 8 and 10 weeks even without reaching a level of statistical significance. Similar tendency was also noted for IVS thickness. As expected, ejection fraction and fractional shortening decreased with time, again without reaching a statistical significance. Myocardial hypertrophy was additionally confirmed by histological analysis of heart tissue. Counting the minimal diameter (Dmin) of cardiomyocytes revealed continuous increase in myocyte size with significant increase after 8 and 10 weeks of the aortic banding (Figures 1A,B). To validate cardiac hypertrophy on molecular level we measured mRNA expression Nppa, encoding natriuretic peptides A, using qRT-PCR. Earlier, this gene was shown to be a reliable marker of cardiomyocyte hypertrophic response and activation of embryonic gene program (Sergeeva and Christoffels, 2013). Increased expression of Nppa was shown in LV and IVS myocardium on all time points of aortic constriction compared to intact and sham operated animals (Figures 1C,D). In right ventricle qRT-PCR analysis revealed low level of Nppa expression regardless of the time point and duration of aortic banding (data not shown). Furthermore, we observed a positive correlation between mRNA level of Nppa in LV and LVM index calculated by echocardiography (r = 0.66, P < 0.01) (Figure 1E). Thus, morphological examination and gene expression signature documented the development of myocardial hypertrophy by weeks 8–10 in current experimental model.

TABLE 2
www.frontiersin.org

Table 2. Echocardiographic parameters obtained after different periods of aortic banding.

FIGURE 1
www.frontiersin.org

Figure 1. Morphological and molecular evidences of myocardial hypertrophy progression. (A) Immunohistochemical staining of beta- sarcoglycan allowing to detect outer membrane showing clear increase of cells size in transverse orientation by 10-weeks group compared to intact. (B) Increase in cell diameter (Dmin) illustrating progressive cardiomyocytes enlargement, and significant increase after 8 and 10 weeks of aortic banding performing. Nppa expression was upregulated in LV (C) and IVS (D) after 1, 2, 8, and 10 weeks of model duration compared to intact or sham-control groups. (E) Positive linear correlation was found for Nppa mRNA level and left ventricular mass (LVM) indexed to body weight. The analysis included experimental groups after 8 and 10 weeks of aortic constriction, 8 week’s sham-operated and intact animals, r indicates Pearson coefficient; for all P < 0.05.

Expression Patterns of Genes, Coding Z-Disk Proteins, in Different Parts of Myocardium

Considering the role of Z-disk as mediator of external mechanical stimulus, we chose several genes, encoding both structural and signaling Z-disk-associated proteins (Actn2, Ldb3, Cmya5, Fhl1, Fhl2, Csrp3, Nebl, Myoz2), to assess gene expression profiles under pressure overload conditions at different time points of LV hypertrophy development. The Ilk-encoding gene was added this list as it was previously bioinformatically shown to be tightly linked to Z-disk-associated genes (Kostareva et al., 2016). Firstly, we studied time-dependent changes of mRNA level in left ventricular free wall myocardium (LV). Figure 2 illustrates that a set of genes – Actn2 (encoding alpha-actinin-2), Cmya5 (myospryn), Ldb3 (LIM domain-binding protein 3) and Nebl (nebulette) – had similar expression profiles after aortic banding procedure. Particularly, expression of Actn2, Cmya5, and Nebl was significantly reduced after 1 week of aortic banding, normalized to control level by week 2 and subsequently decreased by weeks 8 and 10 (Figure 2A). Thus, we observed a significant decline in mRNA expression of some Z-disk associated genes at early stage after aortic banding and upon hypertrophy progression in LV free wall.

FIGURE 2
www.frontiersin.org

Figure 2. Altered mRNA expression of genes, encoding Z-disk proteins, in different parts of myocardium. (A) Actn2, Cmya5, and Nebl expression patterns after different period of aortic banding measured by qRT-PCR – downregulation after 1 week with subsequent return to normal level and further decrease after 10 weeks. Expression Ldb3 declined after 1, 2, and 10 weeks of aortic banding. (B) Altered expression of genes, encoding Z-disk proteins, in right ventricle (RV). Actn2, Cmya5, and Ldb3 expression was gradually increasing after 2 weeks up to 10 weeks of aortic constriction. Despite the gradual decrease of Nebl expression in RV, we showed no further Nebl reduction to the 10 week’s time point. Notably, expression pattern in RV is clearly distinct from pattern in LV. (C) Expression pattern of genes, encoding Z-disk proteins, in IVS are distinctly more close to RV than to LV profiles. Cmya and Nebl mRNA levels were elevated on later phases of hypertrophy progression (after 10 weeks). Opposite, Actn2 is upregulated only on the early phase (after 1 week). Ldb3 mRNA level was significantly increased after 1 week as well as after 10 weeks of aortic constriction. For all P < 0.05.

Analysis of pressure overload-induced gene expression changes in right ventricle revealed a pattern opposite to that observed in LV. A gradual increase of Actn2, Cmya5, and Ldb3 expression was observed starting from week 1 to week 10 of experiment (Figure 2B). Of note, the early changes became significant 1 week later compared to LV (by experimental week 2), and expression of genes that were significantly downregulated in LV after 10 weeks, was strongly upregulated in RV at similar time point.

Further, we evaluated mRNA expression of analyzed genes, encoding Z-disk proteins, in IVS myocardial samples. We revealed that mRNA expression pattern of Actn2, Cmya5, Nebl, and Ldb3 in IVS myocardium was more similar to the pattern in RV than LV free wall. In particular, Cmya5, Ldb3, and Nebl expression became elevated by 10 weeks of coarctation, while Actn2 and Ldb3 mRNA levels were significantly upregulated already after 1 week of aortic banding followed by reduction to normal level and increase at the later time points (Figure 2C).

Other genes selected for this study did not demonstrate defined change in expression pattern, and did not differ significantly between various time points and parts of myocardium (Supplementary Figure 1).

Relation Between mRNA Level and Protein Expression in Different Parts of Myocardium Under Hemodynamic Pressure

To check whether observed alterations in gene mRNA expression lead to changes also on protein level we investigated the amount of Actn2 in myocardial tissue by immunoblotting at 1 and 10 week’s time points which correspond to maximum change on gene expression level. Contrary to reduction of Actn2 mRNA in LV, Actn2 protein level demonstrated increase on 10 weeks after aortic banding compared to control animals (Figures 3A,B). In contrast, despite the increase of Actn2 mRNA in IVS and RV, Actn protein expression did not significantly change at any investigated time points. Therefore, alteration of mRNA expression of Z-disk associated genes on early and later time points of myocardial hypertrophy does not always result in corresponding change in protein level, possibly due to modulation on protein turnover pathways. To check the possible involvement of protein turnover pathways in modulation of protein level under pressure overload conditions, p62 expression was evaluated in LV, RV and IVS myocardium by immunoblotting. We observed increase of p62 level after 10 weeks in LV myocardium, while no alterations were detected for RV and IVS (Supplementary Figure 2). Thus, changes in mRNA expression can be modulated through protein degradation system upon development of LV hypertrophy.

FIGURE 3
www.frontiersin.org

Figure 3. Expression of Actn in left, right ventricle and interventricular septum myocardium. Representative images of Western blotting of Actn2 in LV, IVS, and RV myocardium and total protein staining (Ponceau S) as loading control (A). Densitometry analysis revealed significant increase of Actn2 expression only on 10 week’s time point in LV, and not in RV and IVS myocardium. Expression levels of Actn2 were normalized to Ponceau S staining and control group (B). For all P < 0.05.

Discussion

In cardiomyocytes mechanotransduction serves for transmission of the mechanical stimulus from membrane to the nucleus and involves several cell compartments – sarcolemmal protein complexes, intercalated disks, sarcomere and nuclear envelope (Lyon et al., 2015). Z-disks, located at the lateral borders of sarcomeres, are considered as one of the nodal points of signal transmission. A number of Z-associated proteins due to their domain organization may perform both structural and signaling properties and, thus, participate in regulation of gene expression (Frank et al., 2006). The potential signaling role of Z-disk is to perceive stress stimuli from sarcomere stretch and transmit it to the nucleus in order to adapt cardiomyocyte contractility, compliance and energetics to hemodynamic load (Hoshijima, 2006). Some of the proteins, initially located to Z-band region only, like muscle LIM protein (MLP) or myopodin (SYNPO2), can directly shuttle to the nucleus and enhance the action of transcription factors (Weins et al., 2001; Boateng et al., 2007).

In our study we identified time-dependent alterations of mRNA level of genes encoding Z-disk-associated proteins after induction of pressure-induced cardiac hypertrophy. We revealed that under LV pressure overload, expression of Actn2, Cmya5, Ldb3, and Nebl in LV myocardium was repressed at an early time point (1 week) with subsequent return to basal level, and thereafter again decrease at a late time point (10 weeks). Similar acute-phase regulation of cytoskeletal genes in LV was previously demonstrated by van den Bosch et al. (2006) in transverse aortic constriction (TAC) model. Together with our results, observed changes of transcript level of cytoskeletal and myofibrillar components in LV at the early time points may reflect an initial adaptive response for severe hemodynamic load.

We further detected the divergent gene expression profiles in left, right ventricles and interventricular septum during aortic constriction. As expected, after TAC gene expression patterns in RV and LV were different due to different hemodynamic stress. Surprisingly, IVS expression pattern of Z-disk associated genes was more similar to RV, rather than to LV. During embryonic heart development, LV and RV have a distinct cellular source: LV originates from primary heart field, whereas RV is developed from secondary heart field (Zaffran et al., 2004). Hemodynamic and molecular characteristics in RV stress response are less studied compared to LV (Voelkel et al., 2006). Transcriptome analysis in comparative investigations of pulmonary artery constriction (RV hypertrophy) and TAC (LV hypertrophy) revealed more abundant expression changes in RV than in LV and, moreover, underlined some ventricular-specific expression signatures of genes, responsible for cytoskeletal organization (Kreymborg et al., 2010; Friehs et al., 2013).

The development of IVS start from protrusions of both ventricles, and at a cellular level IVS has equal origin with contribution of cell growth from both LV and RV (Meilhac et al., 2004). Analysis of progenitor cardiac cells identified in the IVS demonstrated the presentation of this population as well in LV and notably not in RV, further supporting IVS similarity to LV origin (Stadtfeld et al., 2007). Increase of Nppa mRNA expression in IVS myocardium in our model confirms the activation of hypertrophic response program after TAC. However, expression pattern of genes, encoding for Z-disk proteins in IVS was more similar to RV rather than to LV. The latter fact suggests that the regulatory role of these genes differs in free LV wall and IVS. This observation might be important for unraveling molecular mechanisms of inherited cardiac hypertrophy where a single gene defect often leads to a marked difference in hypertrophic response in IVS and LV free wall. Thus, the obtained data may be of importance for understanding the molecular pathogenesis of hypertrophic cardiomyopathy (HCM) as well as other inherited cardiac pathologies. While the majority of HCM cases are caused by pathogenic variants in sarcomeric genes, the concomitant rare variants in structural and cytoskeletal genes and their expression profile may influence the phenotypic variability of cardiomyopathies (Lopes et al., 2015).

Changes of mRNA level can have impact not only on protein synthesis directly but also on regulatory processes. Previous investigation of transgenic mice with diastolic dysfunction revealed increase of desmin and α-actinin protein expression without altered transcriptional level (Sheng et al., 2016). In the present study, we also observed discrepancies in LV between mRNA and protein level of α-actinin at the late time point (10 weeks) of hemodynamic load. Observed divergence may partially be explained by concomitant increase of p62 protein level reflecting the repressed protein degradation (Klionsky et al., 2012) further supporting the possibility of regulatory rather than only structural role for Z-disk associated genes (Zhu et al., 2007; Gu et al., 2016). However, other protentional explanation may exist, for example, different life times of the molecules and different mechanisms of degradation.

In summary, we revealed that mRNA expression of genes encoding Z-disk proteins, is mostly reduced in LV myocardium on early and late phases of pressure overload hypertrophy. We demonstrated that expression levels of the same genes in RV have opposite to LV profile with IVS expression profile being more similar to RV rather than free LV. We observed that alteration in mRNA expression upon myocardial hypertrophy does not always correspond to protein level which could be modulated by protein turnover pathways and point on the fact that transcriptional changes do not entirely aim on adjusting protein level, but also on regulating the hypertrophic response during cardiac remodeling. Future integrated data using high throughput screening methods including RNA-sequencing and proteomic approach may provide further information about mechano-dependent gene expression profile in cardiovascular pathology.

Author Contributions

AKn, AKr, AG, AM, GP, and AH performed the experiments. AKn, AKr, GP, and AH performed the statistical data analysis. AKo, AG, NS, and MG designed the study. AKn drafted the manuscript. TS, GS, and AKo revised and approved the final manuscript. All authors agreed to be accountable for all aspects of the work.

Funding

This work was funded by Russian Science Foundation (Grant/Award No. 14-15-00745-II).

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We acknowledge Institute of Experimental Medicine in Almazov National Medical Research Centre Animal for performing all laboratory animal care.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2018.00684/full#supplementary-material

References

Boateng, S. Y., Belin, R. J., Geenen, D. L., Margulies, K. B., Martin, J. L., Hoshijima, M., et al. (2007). Cardiac dysfunction and heart failure are associated with abnormalities in the subcellular distribution and amounts of oligomeric muscle LIM protein. Am. J. Physiol. Heart Circ. Physiol. 7342, 259–269. doi: 10.1152/ajpheart.00766.2006

PubMed Abstract | CrossRef Full Text | Google Scholar

Cho, J. S., Cho, E. J., Lee, J., Choi, H. D., Park, K. C., Lee, K. H., et al. (2014). Myocardial mechanics in a rat model with banding and debanding of the ascending aorta. J. Cardiovasc. Ultrasound 22, 189–195. doi: 10.4250/jcu.2014.22.4.189

PubMed Abstract | CrossRef Full Text | Google Scholar

Dal Ferro, M., Stolfo, D., Altinier, A., Gigli, M., Perrieri, M., Ramani, F., et al. (2017). Association between mutation status and left ventricular reverse remodelling in dilated cardiomyopathy. Heart 103, 1704–1710. doi: 10.1136/heartjnl-2016-311017

PubMed Abstract | CrossRef Full Text | Google Scholar

Devereux, R. B., Alonso, D. R., Lutas, E. M., Gottlieb, G. J., Campo, E., Sachs, I., et al. (1986). Echocardiographic assessment of left ventricular hypertrophy: comparison to necropsy findings. Am. J. Cardiol. 57, 450–458. doi: 10.1016/0002-9149(86)90771-X

CrossRef Full Text | Google Scholar

Epstein, N. D., and Davis, J. S. (2003). Sensing stretch is fundamental. Cell 112, 147–150. doi: 10.1016/S0092-8674(03)00037-0

CrossRef Full Text | Google Scholar

Frank, D., Kuhn, C., Katus, H. A., and Frey, N. (2006). The sarcomeric Z-disc: a nodal point in signalling and disease. J. Mol. Med. 84, 446–468. doi: 10.1007/s00109-005-0033-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Friehs, I., Cowan, D. B., Choi, Y. H., Black, K. M., Barnett, R., Bhasin, M. K., et al. (2013). Pressure-overload hypertrophy of the developing heart reveals activation of divergent gene and protein pathways in the left and right ventricular myocardium. Am. J. Physiol. Heart Circ. Physiol. 304, H697–H708. doi: 10.1152/ajpheart.00802.2012

PubMed Abstract | CrossRef Full Text | Google Scholar

Garnier, S., Hengstenberg, C., Lamblin, N., Dubourg, O., De Groote, P., Fauchier, L., et al. (2015). Involvement of BAG3 and HSPB7 loci in various etiologies of systolic heart failure: results of a european collaboration assembling more than 2000 patients. Int. J. Cardiol. 189, 105–107. doi: 10.1016/j.ijcard.2015.04.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Gautel, M., and Djinovic-Carugo, K. (2016). The sarcomeric cytoskeleton: from molecules to motion. J. Exp. Biol. 219, 135–145. doi: 10.1242/jeb.124941

PubMed Abstract | CrossRef Full Text | Google Scholar

Gu, J., Hu, W., Song, Z. P., Chen, Y. G., Zhang, D. D., and Wang, C.-Q. (2016). Rapamycin inhibits cardiac hypertrophy by promoting autophagy via the MEK/ERK/Beclin-1 pathway. Front. Physiol. 7:104. doi: 10.3389/fphys.2016.00104

PubMed Abstract | CrossRef Full Text | Google Scholar

Hoshijima, M. (2006). Mechanical stress-strain sensors embedded in cardiac cytoskeleton: z disk, titin, and associated structures. Am. J. Physiol. Heart Circ. Physiol. 290, H1313–H1325. doi: 10.1152/ajpheart.00816.2005

PubMed Abstract | CrossRef Full Text | Google Scholar

Jaalouk, D. E., and Lammerding, J. (2009). Mechanotransduction gone awry. Nat. Rev. Mol. Cell. Biol. 10, 63–73. doi: 10.1038/nrm2597

PubMed Abstract | CrossRef Full Text | Google Scholar

Kimura, A. (2016). Molecular genetics and pathogenesis of cardiomyopathy. J. Hum. Genet. 61, 41–50. doi: 10.1038/jhg.2015.83

PubMed Abstract | CrossRef Full Text | Google Scholar

Klionsky, D. J., Abdalla, F. C., Abeliovich, H., Abraham, R. T., Acevedo-Arozena, A., Adeli, K., et al. (2012). Guidelines for the use and interpretation of assays for monitoring autophagy. Autophagy 8, 445–544. doi: 10.4161/auto.19496

CrossRef Full Text | Google Scholar

Knöll, R., Buyandelger, B., and Lab, M. (2011). The sarcomeric Z-disc and Z-discopathies. J. Biomed. Biotechnol. 2011, 569–628. doi: 10.1155/2011/569628

PubMed Abstract | CrossRef Full Text | Google Scholar

Kostareva, A., Kiselev, A., Gudkova, A., Frishman, G., Ruepp, A., Frishman, D., et al. (2016). Genetic spectrum of idiopathic restrictive cardiomyopathy uncovered by next-generation sequencing. PLoS One 11:e0163362. doi: 10.1371/journal.pone.0163362

PubMed Abstract | CrossRef Full Text | Google Scholar

Kreymborg, K., Uchida, S., Gellert, P., Schneider, A., Boettger, T., Voswinckel, R., et al. (2010). Identification of right heart-enriched genes in a murine model of chronic outflow tract obstruction. J. Mol. Cell. Cardiol. 49, 598–605. doi: 10.1016/j.yjmcc.2010.07.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Lang, R. M., Badano, L. P., Mor-Avi, V., Afilalo, J., Armstrong, A., Ernande, L., et al. (2015). Recommendations for cardiac chamber quantification by echocardiography in adults: an update from the American society of echocardiography and the European association of cardiovascular imaging. Eur. Heart J. Cardiovasc. Imaging 16, 233–271. doi: 10.1093/ehjci/jev014

PubMed Abstract | CrossRef Full Text | Google Scholar

Livak, K. J., and Schmittgen, T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCT method. Methods 25, 402–408. doi: 10.1006/meth.2001.1262

PubMed Abstract | CrossRef Full Text | Google Scholar

Lopes, L. R., Syrris, P., Guttmann, O. P., O’Mahony, C., Tang, H. C., Dalageorgou, C., et al. (2015). Novel genotype-phenotype associations demonstrated by high-throughput sequencing in patients with hypertrophic cardiomyopathy. Heart 101, 294–301. doi: 10.1136/heartjnl-2014-306387

PubMed Abstract | CrossRef Full Text | Google Scholar

Luther, P. K. (2009). The vertebrate muscle Z-disc: sarcomere anchor for structure and signalling. J. Muscle Res. Cell Motil. 30, 171–185. doi: 10.1007/s10974-009-9189-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Lyon, R. C., Zanella, F., Omens, J. H., and Sheikh, F. (2015). Mechanotransduction in cardiac hypertrophy and failure. Circ. Res. 116, 1462–1476. doi: 10.1161/CIRCRESAHA.116.304937

PubMed Abstract | CrossRef Full Text | Google Scholar

Meilhac, S. M., Esner, M., Kerszberg, M., Moss, J. E., and Buckingham, M. E. (2004). Oriented clonal cell growth in the developing mouse myocardium underlies cardiac morphogenesis. J. Cell Biol. 164, 97–109. doi: 10.1083/jcb.200309160

PubMed Abstract | CrossRef Full Text | Google Scholar

Pyle, W. G., and Solaro, R. J. (2004). At the crossroads of myocardial signaling: the Role of Z-discs in intracellular signaling and cardiac function. Circ. Res. 94, 296–305. doi: 10.1161/01.RES.0000116143.74830.A9

PubMed Abstract | CrossRef Full Text | Google Scholar

Schänzer, A., Rupp, S., Gräf, S., Zengeler, D., Jux, C., Akintürk, H., et al. (2018). Dysregulated autophagy in restrictive cardiomyopathy due to Pro209Leu mutation in BAG3. Mol. Genet. Metab. 123, 388–399. doi: 10.1016/j.ymgme.2018.01.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Sergeeva, I. A., and Christoffels, V. M. (2013). Regulation of expression of atrial and brain natriuretic peptide, biomarkers for heart development and disease. Biochim. Biophys. Acta 1832, 2403–2413. doi: 10.1016/j.bbadis.2013.07.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Sheng, J. J., Feng, H. Z., Pinto, J. R., Wei, H., and Jin, J.-P. (2016). Increases of desmin and α-actinin in mouse cardiac myofibrils as a response to diastolic dysfunction. J. Mol. Cell. Cardiol. 99, 218–229. doi: 10.1016/j.yjmcc.2015.10.035

PubMed Abstract | CrossRef Full Text | Google Scholar

Stadtfeld, M., Ye, M., and Graf, T. (2007). Identification of interventricular septum precursor cells in the mouse embryo. Dev. Biol. 302, 195–207. doi: 10.1016/j.ydbio.2006.09.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Tabish, A. M., Azzimato, V., Alexiadis, A., Buyandelger, B., and Knöll, R. (2017). Genetic epidemiology of titin-truncating variants in the etiology of dilated cardiomyopathy. Biophys. Rev. 9, 207–223. doi: 10.1016/j.ydbio.2006.09.025

PubMed Abstract | CrossRef Full Text | Google Scholar

van den Bosch, B. J. C., Lindsey, P. J., Burg, C. M. M., Vlies, S. A., Lips, D. J., van der Vusse, G. J., et al. (2006). Early and transient gene expression changes in pressure overload-induced cardiac hypertrophy in mice. Genomics 88, 480–488. doi: 10.1016/j.ygeno.2006.04.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Villard, E., Perret, C., Gary, F., Proust, C., Dilanian, G., Hengstenberg, C., et al. (2011). A genome-wide association study identifies two loci associated with heart failure due to dilated cardiomyopathy. Eur. Heart J. 32, 1065–1076. doi: 10.1093/eurheartj/ehr105

PubMed Abstract | CrossRef Full Text | Google Scholar

Voelkel, N. F., Quaife, R. A., Leinwand, L. A., Barst, R. J., McGoon, M. D., Meldrum, D. R., et al. (2006). Right ventricular function and failure: report of a national heart, lung, and blood institute working group on cellular and molecular mechanisms of right heart failure. Circulation 114, 1883–1891. doi: 10.1161/CIRCULATIONAHA.106.632208

PubMed Abstract | CrossRef Full Text | Google Scholar

Weins, A., Schwarz, K., Faul, C., Barisoni, L., Linke, W. A., and Mundel, P. (2001). Differentiation- and stress-dependent nuclear cytoplasmic redistribution of myopodin, a novel actin-bundling protein. J. Cell Biol. 155, 393–403. doi: 10.1083/jcb.200012039

PubMed Abstract | CrossRef Full Text | Google Scholar

Zaffran, S., Kelly, R. G., Meilhac, S. M., Buckingham, M. E., and Brown, N. A. (2004). Right ventricular myocardium derives from the anterior heart field. Circ. Res. 95, 261–268. doi: 10.1161/01.RES.0000136815.73623.BE

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhu, H., Tannous, P., Johnstone, J. L., Kong, Y., Shelton, J. M., Richardson, J. A., et al. (2007). Cardiac autophagy is a maladaptive response to hemodynamic stress. J. Clin. Invest. 117, 1782–1793. doi: 10.1172/JCI27523

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: Z-disk, cardiac hypertrophy, mechanotransduction, gene expression, pressure overload

Citation: Knyazeva A, Krutikov A, Golovkin A, Mishanin A, Pavlov G, Smolina N, Hushkina A, Sejersen T, Sjoberg G, Galagudza M and Kostareva A (2019) Time- and Ventricular-Specific Expression Profiles of Genes Encoding Z-Disk Proteins in Pressure Overload Model of Left Ventricular Hypertrophy. Front. Genet. 9:684. doi: 10.3389/fgene.2018.00684

Received: 03 October 2018; Accepted: 07 December 2018;
Published: 07 January 2019.

Edited by:

M. Z. A. Bhuiyan, Lausanne University Hospital (CHUV), Switzerland

Reviewed by:

Andreas Brodehl, Heart and Diabetes Center North Rhine-Westphalia, Germany
Boris Shenkman, Institute of Biomedical Problems (RAS), Russia

Copyright © 2019 Knyazeva, Krutikov, Golovkin, Mishanin, Pavlov, Smolina, Hushkina, Sejersen, Sjoberg, Galagudza and Kostareva. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Anastasia Knyazeva, knyazeva_aa@almazovcentre.ru