YTHDC1 expression is highest in embryonic and early postnatal stages and gradually decreases with age
We assayed the temporal and spatial expression profiles of each transcriptome during mouse heart development from a published Gene Expression Omnibus database (www.ncbi.nlm.nih.gov/geo) with accession number GSE51483. Within this gene series, it was observed that the expression of Ythdc1 was highest during the embryonic stage. However, as age progresses from birth to adulthood, the expression gradually declined (Figure 1A/B), with a similar expression profile in both the left and right ventricles. To verify this expression pattern, mouse heart tissues were collected at various developmental stages, including the embryonic period. After Western-blot verification, the expression level of YTHDC1 was highest during embryonic development and the early postnatal phase, while its expression was significantly decreased in adulthood (Figure 1C/D), consistent with the data in the database. This expression profile suggests that YTHDC1 potentially plays a critical role in mouse cardiac development.
Cardiac-specific deletion of Ythdc1 causes abnormal heart development
To clarify the role YTHDC1 plays in heart development, we gathered human fetal hearts from the hospital aborted because of congenital heart disease. We also collected hearts without heart disease from aborted fetuses with similar gestational age as a control group for comparison. We observed a decrease in YTHDC1 expression in the heart tissue of a fetus with ventricular septal defect, aortic stenosis, and aortic arch dysplasia (Figure 2A). Additionally, we created heart specific Ythdc1 knockout mice (Ythdc1-CKO) by inserting loxp sites at both ends of Ythdc1 exons 5-7 and breeding them with Myh6-Cre+/- mice (Figure 2B). The results of the Western blot analysis have confirmed that Ythdc1 exhibited a specific knockdown of approximately 50% in the heart (Figure 2C). We collected liver, muscle, kidney, brain and lung tissues from Ythdc1-CKO and Ythdc1flox/flox mice, Western-blot analysis demonstrated that the YTHDC1 ablation had no significant impact on the expression level in these tissues except the heart (Figure 2D). We observed that the Ythdc1-CKO mice were born without any problems and followed the expected Mendelian inheritance patterns. However, 90% of Ythdc1-CKO mice died within 20 days after born (Figure 2E). To determine the cause of death, we dissected the mice and discovered that the hearts of the Ythdc1-CKO mice were enlarged (Figure 2F). We killed mice 15 days after birth, collected ventricular muscle tissue, and measured the body weight, heart weight, and heart weight/body weight ratio. Statistical analyses revealed a considerable decrease in body weight and a substantial increase in heart weight of the Ythdc1-CKO mice. Moreover, there was a significant increase in the ratio of heart weight to body weight (Figure 2G). The tissues were then fixed and examined using H&E staining and electron microscopy. H&E staining showed enlargement of the entire heart with thinning of the ventricular wall, similar to the characteristics of DCM (Figure 2H). Transmission electron microscopy reveals partial disorganization of muscle fibers, with certain Z-lines (Z) fractured and a minor absence of H-band (H). Additionally, certain mitochondria exhibit vacuolar degeneration (Figure 2I). To further assess cardiac function quantitatively, we performed echocardiography on days 1 and 15 after birth using a small animal ultrasound device (Vevo3100) to measure left ventricular diameter, thickness and other indicators. No significant differences between Ythdc1-CKO and Ythdc1flox/flox mice were observed in any echocardiographic or physiological parameters at 1 day of age. However, on day 15 after birth, a significant decrease in ejection fraction (EF) and fractional shortening (FS) was observed, along with a significant decrease in the end-diastolic diameter of the interventricular septum (IVS, d), while the end-diastolic diameter of the left ventricle (LVID, d) significantly increased (Figure 2J/K/L). In summary, our findings indicate that Ythdc1-CKO mice are born with normal hearts, then gradually develop enlarged hearts with severe heart failure similar to DCM, which directly results in the premature death of KO mice. These results demonstrate that YTHDC1 is of great importance in heart development, and specific deletion of YTHDC1 in the heart causes structural and functional changes in both human and mouse hearts.
Cardiac-specific deletion of Ythdc1 largely changes gene expression profile in the heart
To further explore the effect of heart-specific YTHDC1 deletion on the molecular mechanisms associated with cardiac development abnormalities, we conducted RNA-Seq analysis on heart tissue samples obtained from Ythdc1-CKO and Ythdc1flox/flox mice, aiming to examine the transcriptional profile comprehensively. Figure 3A illustrates that 437 genes were down-regulated and 599 genes were up-regulated in the Ythdc1-CKO heart (Figure 3A). GO analysis revealed that: The down-regulated genes were mainly associated with heart contraction, heart processes, regulation of ion transmembrane transport and regulation of ion transmembrane transport, regulation of blood circulation, calcium ion transport, cardiac muscle contraction, regulation of heart contraction, divalent metal ion transport, and divalent inorganic cation transport (Figure 3B), while the upregulated gene was associated with extracellular structure organization, extracellular matrix organization, positive regulation of the cell migration, positive regulation of cell motility, cartilage development, positive regulation of response to external stimulus, regulation of ossification, skeletal system development, and biomineral tissue development (Figure 3C). KEGG pathway enrichment analysis showed that pathways such as adrenergic signaling in cardiomyocytes, cardiac muscle contraction, neuroactive ligand-receptor interaction, arrhythmogenic right ventricular cardiomyopathy (ARVC), calcium signaling pathway, HCM, vascular smooth muscle contraction, the MAPK signaling pathway, DCM, and the serotonergic synapse were downregulated. (Figure 3D). The up-regulated genes were mainly related to the biosynthesis of amino acids, HIF-1 signaling pathway, ECM-receptor interaction, focal adhesion, PI3K-Akt signaling pathway, ErbB signaling pathway, TGF-β signaling pathway, HCM, and p53 signaling (Figure 3E). Previous studies have shown that Ythdc1 can regulate alternative splicing[18], so we used rMATS software to analyze alternative splicing events in RNA-seq data, and there were 250 significantly different alternative splicing events, among which: skipped exon (SE) accounted for 65.2%, retained intron (RI) accounted for 11.6%, mutually exclusive exon (MXE) accounted for 9.2%. alternative 3' splice site (A3SS) accounted for 8% and alternative 5' splice site (A5SS) for 6% (Figure 3F). GO analysis was performed on 250 differential alternative splicing events. The differentially spliced genes were mainly related to muscle contraction, muscle filament sliding, regulation of phagocytosis, sarcomere organization, chromatin organization, cardiac muscle contraction, ventricular cardiac muscle tissue morphogenesis, endocytosis, and actin cytoskeleton organization (Figure 3G). These data suggest that heart-specific YTHDC1 deletion dramatically alters gene expression profiles in the heart.
Cardiac-specific deletion of Ythdc1 significantly alters chromatin accessibility in the heart
Previous reports have indicated that several RBPs can participate in transcriptional regulation and influence chromatin through direct or indirect mechanisms[16, 28-30]. To further investigate whether YTHDC1 regulates chromatin accessibility and gene transcription in the heart, the transposase accessible chromatin sequencing (ATAC-seq) technique was used to map the open chromatin in the heart of Ythdc1-CKO mice. As expected, significant enrichments of open chromatin were observed near gene promoters and transcription start sites (TSS) within ATAC-seq peaks (Figure 4A/B). Analysis of the ATAC-seq data revealed that 2326 peak-related genes were up-regulated and 5566 peak-related genes were down-regulated (Figure 4C, Supplementary material online, Table S1). We identified several transcription factor (TF) binding motifs using the motif analysis software HOMER. The binding motifs of five transcription factors (NF1, Mef2a, Mef2b, Mef2c, and Mef2d) were switched off in Ythdc1-CKO mouse hearts (Figure 4D). These data suggest that YTHDC1 regulates chromatin accessibility in the heart. The following analysis of the ATAC-seq and RNA-seq data showed that 141 genes were down-regulated in both datasets (Figure 4E and Supplementary material online, Table S2). The KEGG analysis revealed the down-regulation of genes associated with adrenergic signaling in cardiomyocytes, MAPK signaling pathway, cardiac muscle contraction, calcium signaling pathway, and glutamatergic synapse (Figure 4F).