A human mitofusin 2 mutation can cause mitophagic cardiomyopathy
Abstract
Cardiac muscle has the highest mitochondrial density of any human tissue, but mitochondrial dysfunction is not a recognized cause of isolated cardiomyopathy. Here, we determined that the rare mitofusin (MFN) 2 R400Q mutation is 15–20× over-represented in clinical cardiomyopathy, whereas this specific mutation is not reported as a cause of MFN2 mutant-induced peripheral neuropathy, Charcot–Marie–Tooth disease type 2A (CMT2A). Accordingly, we interrogated the enzymatic, biophysical, and functional characteristics of MFN2 Q400 versus wild-type and CMT2A-causing MFN2 mutants. All MFN2 mutants had impaired mitochondrial fusion, the canonical MFN2 function. Compared to MFN2 T105M that lacked catalytic GTPase activity and exhibited normal activation-induced changes in conformation, MFN2 R400Q and M376A had normal GTPase activity with impaired conformational shifting. MFN2 R400Q did not suppress mitochondrial motility, provoke mitochondrial depolarization, or dominantly suppress mitochondrial respiration like MFN2 T105M. By contrast to MFN2 T105M and M376A, MFN2 R400Q was uniquely defective in recruiting Parkin to mitochondria. CRISPR editing of the R400Q mutation into the mouse Mfn2 gene induced perinatal cardiomyopathy with no other organ involvement; knock-in of Mfn2 T105M or M376V did not affect the heart. RNA sequencing and metabolomics of cardiomyopathic Mfn2 Q/Q400 hearts revealed signature abnormalities recapitulating experimental mitophagic cardiomyopathy. Indeed, cultured cardiomyoblasts and in vivo cardiomyocytes expressing MFN2 Q400 had mitophagy defects with increased sensitivity to doxorubicin. MFN2 R400Q is the first known natural mitophagy-defective MFN2 mutant. Its unique profile of dysfunction evokes mitophagic cardiomyopathy, suggesting a mechanism for enrichment in clinical cardiomyopathy.
Editor's evaluation
The study presents a valuable finding on the link between mitofusin function (MFN2) and Parkin-mediated mitophagy, and that the combination of mitophagy and mitochondrial fusion defects is the basis of cardiomyopathy. The data based on the in vivo models are compelling and provide clinically relevant associations. This finding is important for understanding mitochondrial defects as a basis of heart pathologies.
https://doi.org/10.7554/eLife.84235.sa0eLife digest
Mitochondria are organelles with an essential role in providing energy to the cells of the body. If damaged, they are repaired by fusing and exchanging contents with sister mitochondria in a process that requires mitofusin proteins.
While mutations in the gene for mitofusin 2 have been linked to nerve damage, they do not appear to affect the heart – despite high concentrations of mitochondria in heart muscle cells. However, previous research showed that experimentally disrupting the programmed removal of mitochondria, a process also regulated by mitofusin 2, can cause heart muscle disease known as cardiomyopathy. This suggests that mutations affecting different mitofusin 2 roles might harm individual cell types in different ways.
To investigate, Franco et al. carried out a genetic screen of people with cardiomyopathy, identifying a rare mitofusin 2 mutation, called R400Q, that was more common in this group. Experiments showed that R400Q caused cardiomyopathy in mice and affected mitochondrial repair and replacement, but not movement. By contrast, a mutation linked to Charcot-Marie-Tooth disease type 2A – which causes nerve damage – affected mitochondrial movement but not clearance, leading to nerve cell damage but not cardiomyopathy. This led Franco et al. to suggest that mitochondrial movement is central to nerve cell health, whereas mitochondrial repair and replacement plays an important role in cardiac development.
Genetic cardiomyopathies affect around 1 in 500 people, but only half of the gene mutations responsible are known. These results suggest that mutations affecting mitochondrial quality control factors could be involved, highlighting a direction for future studies into modifiers of cardiomyopathy.
Introduction
Hearts rely upon mitochondria-derived ATP to fuel cardiac development, excitation–contraction–coupling, and myocardial repair after damage or senescence. For this reason, cardiac involvement is common in diseases like MELAS and MERRF caused by genetic defects of the mitochondrial genome (mtDNA) (Hsu et al., 2016). However, the vast majority of the ~1000 constituent mitochondrial proteins are not encoded by mtDNA, but instead by nuclear genes. Among nuclear-encoded mitochondrial proteins having central roles in mitochondrial health and homeostasis are mitofusins (MFN) 1 and 2 (Santel and Fuller, 2001). Gene ablation and mutant transgene expression studies targeting cardiac myocytes indicate that cardiac MFNs are critical mediators of reparative mitochondrial fusion (Chen et al., 2011; Kasahara et al., 2013), mitochondrial culling for quality control (Song et al., 2015; Song et al., 2017), and mitochondrial replacement for metabolic plasticity (Gong et al., 2015).
Clinical and experimental data relating MFN functioning to cardiac health appear somewhat discordant. On the one hand, the cited Mfn gene knockout studies unambiguously demonstrate that mitofusins are essential for heart development, contraction, and repair. On the other hand, cardiac involvement has not been reported in human conditions caused by natural loss-of-function MFN gene mutations. Indeed, dozens of different pathological MFN2 mutations cause neuronal rather than cardiac pathology, manifested as Charcot–Marie–Tooth disease type 2A (CMT2A) (Züchner et al., 2004) or optic degeneration (Züchner et al., 2006). We have proposed that mitochondrial motility defects induced by mutant MFN2 are responsible for clinical neuropathies because neurons have long axonal and dendritic processes that make them critically dependent upon mitochondrial transport (Dorn, 2019; Dorn, 2021). By contrast, impaired mitophagy seems to underlie developmental and senescent cardiomyopathies provoked by Mfn gene ablation or engineered MFN2 mutant expression (Bhandari et al., 2014; Gong et al., 2015; Dorn, 2016).
MFN2 multifunctionality reflects its ability to bind different protein partners as determined by specific phosphorylation events (Dorn, 2020; Li et al., 2022). Thus, MFN2 oligomerization with MFN1 or MFN2 mediates mitochondrial fusion (Chen et al., 2010; Dorn and Dang, 2022), MFN2 binding of Parkin mediates mitophagy (Chen and Dorn, 2013; Li et al., 2022), and MFN2 binding to Miro regulates mitochondrial transport (Misko et al., 2010). We reasoned that mutations affecting different MFN2 domains could evoke a phenotypic spectrum consisting of mitochondrial fragmentation from loss of fusion, mitochondrial dysmotility from impaired transport, and/or mitochondrial degeneration from interrupted mitophagy.
Here, as part of an ongoing effort to understand the relationship between MFN2 mutation site, protein dysfunction, and disease manifestation, we screened cardiomyopathy research cohorts for MFN2 mutations linked to cardiac hypertrophy or heart failure. A largely overlooked, very rare MFN2 mutation, MFN2 R400Q (Eschenbacher et al., 2012), was over-represented in cardiomyopathy among individuals of African descent. We found that fusogenicity of mutant MFN2 Q400 is depressed to the same extent as CMT2A mutants MFN2 T105M and M376 A/V. Uniquely, MFN2 Q400 is defective in recruiting Parkin to mitochondria, thereby impacting an important mitophagic pathway. Knockin (KI) mice engineered to carry Mfn2 Q400 develop a perinatal cardiomyopathy, whereas Mfn2 T105M KI mice develop a CMT2A-like neuropathy and Mfn2 M376V KI mice appear phenotypically normal. To our knowledge, MFN2 R400Q is the first example of a naturally occurring MFN2 mutation that primarily impacts mitophagy. These findings designate mitophagic dysfunction as a new mechanism for cardiomyopathy, supporting targeted investigations of other genetic mitophagy defects as cardiac risk modifiers. Finally, our data support the concept that different patterns of mitofusin dysfunction underlie organ-specific manifestations of clinical diseases caused by MFN2 mutations.
Results
The rare MFN2 R400Q mutation is over-represented in clinical cardiomyopathy
Rare mutations can be overlooked in standard genotype–phenotype association studies, requiring targeted analyses of exome/whole-genome sequence data (Auer and Lettre, 2015; Bomba et al., 2017). Here, to define the landscape of potential disease-causing MFN2 variants (MacArthur et al., 2012) we surveyed the gnomAD database (v2.1.1) of 141,456 individuals for rare (minor allele frequency [MAF] of 0.001–0.01) nonsynonymous MFN2 sequence variants. A total of 352 rare nonsynonymous MFN2 variants are reported (MAF 3.98E-06–6.87 E-03) of which 171 are singletons. Notably, 76 of the rare MFN2 variants are implicated in CMT2A (Stuppia et al., 2015; Figure 1—source data 1).
To detect any MFN2 variants associated with human heart conditions, MFN2 gene exons were resequenced in DNA samples from 398 subjects with idiopathic/non-ischemic dilated cardiomyopathy (DCM) from the Cincinnati Heart Study (Liggett et al., 2008; Matkovich et al., 2010; Cappola et al., 2011), 286 subjects with hypertrophic cardiomyopathy (HCM) from the Houston Hypertrophic Cardiomyopathy cohort (Osio et al., 2007; Li et al., 2017), and 861 non-affected controls from the Cincinnati cohort. Additionally, 424 individuals (281 probands) with familial cardiomyopathy (FCM) from the DCM Precision Medicine Study (Huggins et al., 2022; Trachtenberg et al., 2022) were screened. As summarized in Table 1, 11 non-synonymous MFN2 variants were detected, 6 of which encode putative CMT2A mutations. Also, 8 of the 11 MFN2 variants were singletons and therefore of minimal interest. MFN2 R400Q, which is not a known CMT2A mutation and is bioinformatically predicted to be damaging, was detected in three unrelated individuals: two with HCM and one with DCM (Table 1).
The MFN2 R400Q mutation (Figure 1A and B) is reported almost exclusively in African populations in the gnomAD database (Figure 1—source data 1). Indeed, all three subjects with the R400Q variant in our cardiomyopathy study populations identified as African American. The overall MFN2 R400Q MAF of 0.0015 in cardiomyopathy, when compared to MAF of 0.000074 from the gnomAD database, represents ~20-fold over-representation in the combined heart disease populations (p<0.00001). If only the gnomAD African-derived population is used to account for racial differences in R400Q prevalence (10/24,962), over-representation was ~15-fold (p-value = 0.0033). None of the non-affected controls (625 Caucasian and 236 African American subjects) had the MFN2 R400Q mutation. Because a previous description of this MFN2 mutant described impaired fusogenicity conferring retinal and myocardial abnormalities in transgenic Drosophila (Eschenbacher et al., 2012), we surmised that, like other loss-of-function MFN2 mutants (Flannick et al., 2014; Steinberg et al., 2015), MFN2 Q400 could have pathophysiological relevance. We therefore designated it a candidate cardiomyopathy modifier and proceeded with detailed investigations.
MFN2 Q400 is a dominant suppressor of mitochondrial fusion
Many CMT2A MFN2 mutations, including prototypical MFN2 T105M, are located within the GTPase domain (Figure 1A, green region) and lack catalytic GTPase activity (Franco et al., 2022). By contrast, MFN2 R400Q (Figure 1B) is located in the hydrophobic core within the so-called first heptad repeat domain (HR1; Figure 1A, red region). To better understand the relationship between MFN2 mutation location, protein dysfunction, and pathophysiological consequence, we compared fusogenic activities of mutant MFN2 Q400 to wild-type (WT) MFN2 and two positionally diverse human CMT2A mutants, MFN2 T105M in the GTPase domain and MFN2 M376V/A in the HR1 domain (Figure 1A; Franco et al., 2020; Zhou et al., 2021). An additional CMT2A GTPase domain mutation, MFN2 R94Q, and laboratory-engineered MFN2 K109A (Detmer and Chan, 2007), were included as complementary comparators (see Figure 1, Figure 1—figure supplement 1). Each human MFN2 was expressed at equal levels in murine mitofusin null (Mfn1/Mfn2 double knockout [DKO]) cells using adenoviral transduction (Figure 1C; Li et al., 2019). Importantly, using this system the transfected mitofusin proteins are properly inserted into mitochondria and fully functional (Franco et al., 2022).
Mitochondrial aspect ratio (length/width) is commonly used to gauge the relative activity of mitochondrial fusion vs. fission. The normal mitochondrial aspect ratio of ~5–6 in WT murine embryonic fibroblasts (MEFs) is reduced to ~2 in MEFs lacking both mitofusins (Mfn DKO) (compare white bars in Figure 1D and E). An increase in aspect ratio in Mfn DKO cells reflects intrinsic fusogenic activity of recombinant MFN2 proteins expressed therein (Franco et al., 2016). Consistent with an absolute requirement for GTP hydrolysis during mitochondrial fusion, all three GTPase domain MFN2 mutants exhibited diminished fusogenicity compared to WT MFN2 (Figure 1D, Figure 1—figure supplement 1; left panel, green bars). Intrinsic fusogenicity of MFN2 M376A and R400Q was similarly impaired (Figure 1D, red bars).
Most CMT2A-associated MFN2 mutants do not simply exhibit reduced intrinsic fusogenicity, but dominantly suppress mitochondrial fusion mediated by normal mitofusins 1 and 2. We assayed suppression of mitochondrial fusion by expressing our panel of MFN2 mutants in WT MEFs, which have normal endogenous Mfn1 and Mfn2 and, therefore, a characteristic baseline mitochondrial aspect ratio of 5–6 (Figure 1E, white bar). As expected, expression of WT MFN2 caused mitochondrial hyper-fusion, increasing mitochondrial aspect ratio to ~8 (Figure 1E, black bar). By comparison, all three GTPase domain MFN2 mutants and MFN2 R400Q suppressed normal mitochondrial fusion (i.e. lowered mitochondrial aspect ratio vs. adeno-βgal transduced) in WT cells. Uniquely among the expressed MFN2 mutants, MFN2 M376A did not alter baseline aspect ratio of WT MEF (Figure 1E, Figure 1—figure supplement 1, right panel). Thus, MFN2 R400Q is similar to GTPase domain MFN2 mutants in that it both lacks intrinsic fusogenic activity and dominantly suppresses normal mitochondrial fusion, whereas MFN2 M376A is a straightforward functional null protein without dominant inhibitory effects.
MFN2 HR1 mutants have normal GTPase activity and do not impair mitochondrial respiration
MFN1 and MFN2 are dynamin family GTPases; catalytic GTPase activity is required for mitochondrial fusion (Chan, 2006; Santel and Fuller, 2001; Franco et al., 2022). We asked how suppression of fusion by our panel of MFN2 mutants related to differences in GTPase activity. Adenovirally expressed WT MFN2 accounted for approximately two-thirds of total mitochondrial GTPase activity in transduced mitofusin DKO cells (Figure 2A, black bar). As expected, MFN2 T105M, R94Q, and K109A each lacked GTPase activity (Figure 2A, Figure 2—figure supplement 1, left panel, green bars). Remarkably, GTPase activity of MFN2 R400Q and M376A was the same as WT MFN2, that is, normal (Figure 2A, red bars). Thus, diminished GTPase activity cannot be the mechanism for impaired fusogenicity of these two MFN2 HR1 domain mutants.
Mitochondrial respiration and oxidative phosphorylation of ADP to ATP are driven by an electrochemical gradient across the inner mitochondrial membrane, the so-called ‘proton pump’ (Madeira, 2018). Functional electrochemical inner membrane integrity can be assessed using dyes that are concentrated within mitochondria in a voltage-dependent manner; loss of inner membrane electrical polarization abrogates mitochondrial staining. We used tetramethylrhodamine ethyl ester (TMRE) to assess mitochondrial polarization in WT MEFs expressing our panel of MFN2 mutants. Dissipation of the inner membrane electrochemical gradient (mitochondrial depolarization) was provoked by MFN2 T105M, R94Q, and K109A, but not by MFN2 M376A or R400Q (Figure 2B, Figure 2—figure supplement 1, right panel). Thus, for these MFN2 mutants, the ability to induce mitochondrial depolarization corresponds with impaired catalytic GTPase activity.
The functional consequences of MFN2 mutations on catalytic GTPase activity and electrochemical integrity were further assessed by measuring mitochondrial respiration in DKO MEFs transduced with adeno-MFN2 WT, T105M, M376A, or R400Q. As expected, introducing WT MFN2 into cells lacking mitofusins enhanced all phases of mitochondrial respiration (Figure 2C, black). By comparison, MFN2 T105M failed to enhance mitochondrial respiration (Figure 2C, green). Strikingly, MFN2 M376A and R400Q were intermediate in their effects, increasing mitochondrial respiration vs. the complete absence of mitofusins (i.e. vs. β-gal), but seemingly not to the same extent as WT MFN2 (Figure 2C, red). These data provide further evidence supporting the central relationship between mitochondrial GTPase activity, polarization, and respiration. In the context of our results presented in Figure 1, these results implicate a different, unidentified mechanism by which MFN2 HR1 domain mutants impair mitochondrial fusion and suppress mitochondrial respiration.
MFN2 HR1 domain mutations cause loss of conformational malleability
Mitofusin protein conformation changes upon activation, transitioning from a closed basal state to an open state that favors mitochondrial fusion (Franco et al., 2016). The closed conformation is maintained by peptide–peptide interactions involving MFN2 amino acids 367–384 in the HR1 domain (Rocha et al., 2018; Dorn, 2019). We considered that M376 mutations within the peptide–peptide interaction domain, and the proximity of the MFN2 R400 mutation site to the interacting domain, might alter the ability of these MFN2 mutant proteins to switch conformation. Consistent with this notion, GTPase-defective CMT2A mutant MFN2 T105M transitioned normally from open to closed state, as measured by antagonist and agonist peptide-induced changes in Forster resonance energy transfer (FRET; Dang et al., 2020; Figure 3A, green bars). By contrast, the MFN2 M376A and R400Q FRET signals did not change in the presence of MFN antagonist and agonist peptides that normally provoke the closed or open conformation, respectively (Figure 3A, red bars; Franco et al., 2016). Thus, these MFN2 HR1 domain mutants are defective in their ability to move between open fusion-permissive and closed fusion-limiting conformations, which can explain how mitochondrial fusion is impaired by these mutations despite normal catalytic GTPase activity.
MFN2 HR1 domain mutations do not impair mitochondrial motility
The appellative function of mitofusins is to mediate mitochondrial fusion, but MFN2 has additional non-canonical tasks including to promote mitochondrial transport (Misko et al., 2010; Rocha et al., 2018; Franco et al., 2020). Indeed, some CMT2A-linked MFN2 mutations dominantly suppress mitochondrial motility, thereby contributing to die-back of long peripheral nerves that require distal mitochondrial delivery for cell renewal and regeneration (Misko et al., 2010; Rocha et al., 2018; Franco et al., 2020). Directed mitochondrial transport can readily be measured in cultured neurons (Dang et al., 2022; Dorn and Dang, 2022). Accordingly, we expressed our representative MFN2 mutants in cultured dorsal root ganglion (DRG) neurons derived from wild-type mice. Consistent with previous observations that mitochondrial fusion and motility can be co-regulated by MFN2 (Franco et al., 2022), mitochondrial transport through neuronal processes was suppressed by GTPase-defective MFN2 mutants T105M and K109A (Figure 3B, Figure 3—figure supplement 1A, green bar), and an MFN2 mutant engineered to mimic phosphorylation at sites that suppress mitochondrial fusion, MFN2 EE (Chen and Dorn, 2013; Figure 3—figure supplement 1A, gray bar MFN2 E111, E442). Mitochondrial motility in DRG neurons was unaffected by HR1 domain mutants M376A and R400Q (Figure 3B, red bars) or by an MFN2 mutant engineered to prevent phosphorylation at the fusion-regulatory sites, MFN2 AA (Chen and Dorn, 2013; Figure 3—figure supplement 1A, A111A442). Thus, combined abnormalities of mitochondrial fusion and motility exhibited by GTPase-defective MFN2 mutants are dissociated in MFN2 HR1 domain mutants.
MFN2 Q400 is non-functional as a mitochondrial Parkin receptor
MFN2 can initiate mitophagy by switching from a fusion protein into a ‘Parkin receptor’ that recruits cytosolic Parkin to damaged mitochondria (Chen and Dorn, 2013; Xiong et al., 2019; Li et al., 2022). Mitophagic vs. fusogenic properties are determined by the phosphorylation status of MFN2 T111 and S442 (Li et al., 2022). Prior studies have suggested that MFN2 functioning as a Parkin receptor and mitophagy effector is important to heart development and/or contractile function (Chen and Dorn, 2013; Gong et al., 2015; Song et al., 2015). Mitochondrial Parkin recruitment by MFN2 WT, T105M, M376A, and R400Q was compared to that evoked by the MFN2 phosphorylation site mutants MFN2 EE that constitutively bind Parkin and MFN2 AA that is incapable of binding Parkin (Chen and Dorn, 2013). Parkin translocation from cytosol to mitochondria was stimulated by uncoupling mitochondrial respiration with carbonyl cyanide p-trifluoromethoxyphenylhydrazone (FCCP) and visualized as the appearance of mitochondrial Parkin aggregates (Figure 3—figure supplement 1B). Baseline cytosolic Parkin localization in the absence of FCCP was not affected by WT MFN2 or any of the naturally occurring MFN2 mutants (Figure 3C, Figure 3—figure supplement 1C, FCCP), whereas MFN2 EE promoted spontaneous Parkin aggregation as expected (Figure 3—figure supplement 1C; Chen and Dorn, 2013; Li et al., 2022). Parkin translocation in response to FCCP was unaffected by MFN2 WT, T105M, K109A, and M376A, but was markedly reduced in cells expressing MFN2 R400Q (Figure 3C, Figure 3—figure supplement 1C,+FCCP). Indeed, suppression of Parkin recruitment by MFN2 R400Q was similar in magnitude to that of MFN2AA (Chen and Dorn, 2013; Gong et al., 2015; Figure 3C, Figure 3—figure supplement 1C, +FCCP). Taken together, the results shown in 1—3 reveal a single characteristic of the MFN2 R400Q mutation that distinguishes it from the CMT2A-linked MFN2 mutants: it suppresses Parkin recruitment to depolarized mitochondria.
Mfn2 Q400 knock-in mice develop lethal perinatal cardiomyopathy
Accumulating evidence suggests that MFN-mediated mitochondrial fusion is essential to early embryonic heart development, and that MFN2-Parkin-mediated mitophagy is required for metabolic maturation of perinatal myocardium (Kasahara et al., 2013; Gong et al., 2015). Thus, MFN2 Q400-induced dysfunction of either fusion or mitophagy could underlie its apparent over-representation in human heart disease. The human MFN2 R400Q mutation is sufficiently rare that relating MFN2 Q400 dysfunction to a particular clinical manifestation would require linkage analysis in large families that were not identified in our FCM cohort. Alternately, we reasoned that manifestation of a cardiac phenotype unique to mice carrying the MFN2 R400Q mutation would provide unbiased evidence linking it to in vivo heart disease. To provide continuity with our in vitro studies comparing MFN2 R400Q to CMT2A-associated mutations at T105 and M376, we used CRISPR/Cas9 to engineer human mutations at all three MFN2 sites into the highly homologous mouse Mfn2 gene (Figure 4—figure supplement 1).
Heterozygous knock-in mice from each line were viable, and therefore interbred to expand colonies. Mutant allele frequency of resulting offspring, expected to be ~50%, was significantly lower in Mfn2 T105M mice, and trended lower in Mfn2 R400Q mice (Figure 4A). Genotyping of 7-day-old pups revealed no homozygous mutant Mfn2 T105M mice and an abnormally low percentage of homozygous mutant Mfn2 R400Q mice, whereas mutant allele distribution of Mfn2 M376V mice followed Mendelian norms (Figure 4B). Analysis of embryos from timed crosses revealed no homozygous Mfn2 T105M mice as early as E12.5, recapitulating early embryonic lethality described in germline Mfn1 and Mfn2 knockout mice (Chen et al., 2003); heterozygous Mfn2 T105M mice had normal hearts and a neurological phenotype to be reported separately. Heterozygous and homozygous Mfn2 M376V mice likewise had normal hearts. We interrogated the apparent paucity of homozygous Mfn2 R400Q KI mice by interbreeding surviving homozygotes (Mfn2 Q/Q400 mice). Litter size from homozygous crosses, reflecting survival of Q/Q400 mice, was normal at E18.5. However, approximately one-third of Q/Q400 fetuses died between E18.5 and birth (P0), and another approximately one-third of Q/Q400 pups succumbed by P7, after which there was no further abnormal mortality (Figure 4C).
Histological examination of two litters of P0 Q/Q400 mice by a blinded Washington University veterinary pathologist revealed ‘foci of cardiomyocyte vacuolation’ in the myocardia of over half of the hearts (Figure 4D). No other organs showed any histological or structural abnormalities. Ultrastructural examination revealed loss of myofibrillar integrity with mitochondrial scarcity, disorganization, fragmentation, and degeneration (Figure 4E, Figure 4—figure supplement 2) Thus, Mfn2 Q/Q400 mice develop a highly specific, incompletely penetrant late fetal/early perinatal cardiomyopathy.
Transcriptional and metabolomic profiling of MFN2 Q400 cardiomyopathy identifies signatures specific for cardiac pathology
Transcriptional profiling can provide insights into the mechanisms underlying cardiac phenotypes produced by genetic variation (Aronow et al., 2001). We performed RNA sequencing on E18.5 MFN2 Q/Q400 mouse hearts and their WT controls; E18.5 was selected because there was little mouse drop-out at this time point (see Figure 4C). Unsupervised clustering of dysregulated Q/Q400 cardiac mRNAs (p-value<0.02, FDR < 0.05, fold change >2 vs. wild-type) segregated hearts into those having near-normal mRNA levels (designated nominal; Q400n) and those manifesting an abnormal transcriptional profile (designated anomalous; Q400a) (Figure 5A). Interrogation of mRNA levels of genes central to mitochondrial function, including respiratory enzymes encoded by the mitochondrial genome and nuclear-encoded mitochondrial biogenesis and mitochondrial dynamics genes (Figure 5—figure supplements 1 and 2), confirmed segregation between normal and abnormal hearts. Gene ontology (GO) analysis of dysregulated Q400a transcripts categorized mRNAs with lower expression in Q400a hearts overwhelmingly to metabolic and mitochondrial disease pathways (Figure 5B). By contrast, genes with increased expression in the same hearts were linked to multiple injury responses (Figure 5B), likely reflecting a compensatory reaction to myocardial degeneration. The distinction between Q400a and Q400n hearts was noteworthy for SOD2 and major cardiac expressed GPX mRNAs, respectively encoding mitochondrial superoxide dismutase and glutathione peroxidases that are central to mitigating mitotoxicity from reactive oxygen species (ROS) (Figure 5C).
The abnormal metabolic transcriptional signature of Q400a mouse hearts prompted metabolomics studies of individual newborn mouse hearts. As with gene expression, metabolomics segregated MFN2 Q/Q400 mice into Q400a and Q400n (Figure 5D). Metabolic differentiation between Q400a and Q400n was largely driven by abnormally low levels of a preferred metabolic substrate for fetal myocardium, lactate (Lopaschuk et al., 1992; Piquereau and Ventura-Clapier, 2018) in Q400a hearts. Levels of C2 fatty acid acylcarnitines and several amino acids were also diminished (Figure 5—figure supplement 3). Specific interrogation of metabolic pathway gene expression in MFN2 Q400 fetal hears revealed broad downregulation of mRNAs encoding fatty acid entry, glycolysis, tricarboxylic cycle, and oxidative phosphorylation pathway proteins (Figure 5E–H). Thus, the development of cardiomyopathy in Mfn2 Q/Q400 fetuses and neonates corresponds with signature abnormalities of myocardial metabolism.
Cardiac myocyte mitophagy is suppressed by MFN2 Q400
Our in vitro studies revealed MFN2 Q400 to be defective in recruiting Parkin to mitochondria, suggesting that this mutation might adversely impact Parkin-mediated mitophagy. This notion gained support from our observation that the MFN2 Q400 fetal heart profile of abnormal fatty acid entry, glycolysis, tricarboxylic cycle, and oxidative phosphorylation pathway gene expression (see Figure 5) recapitulates transcriptional changes evoked by experimentally interrupting MFN2-Parkin-mediated mitophagy in perinatal hearts (Gong et al., 2015). Importantly, this mRNA profile has little similarity to the transcriptional signature provoked by suppressing mitochondrial fusion in embryonic myocardium (Kasahara et al., 2013). Given that MFN2 Q400 had only a modest effect on mitochondrial respiration (see Figure 2C), our findings suggested that the disturbances in myocardial metabolism induced by MFN2 Q400 might be a consequence of interrupting MFN2-Parkin-mediated mitophagy.
To better understand the impact of MFN2 R400Q on mitophagy in fetal heart cells, we expressed WT or mutant MFN2s in H9c2 rat cardiomyoblasts that retain features of embryonic cardiac myocytes (Hescheler et al., 1991; Kuznetsov et al., 2015; Figure 6A). First, we defined the time courses for FCCP-induced mitochondrial Parkin translocation (measured with mcParkin; Narendra et al., 2008), mitochondrial engulfment by autophagosomes (measured as mitochondria-LC3-GFP overlay; Jackson et al., 2005), and for mitochondrial delivery to lysosomes (measured using MitoQC; Williams et al., 2017) in H9c2 cells expressing WT MFN2. As expected, different stages of mitophagy progressed with different temporal patterns (Figure 6B). Remarkably, each of the stages of mitophagy was markedly delayed in H9C2 cells expressing MFN2 Q400 at an equivalent level (Figure 6C). Mitophagy slowly increased over longer time periods in MFN2 Q400 cells, likely reflecting MFN2-independent mitochondrial Parkin recruitment (Ordureau et al., 2014).
As in Mfn2 DKO cells (see Figure 3C), CMT2A-associated MFN2 T105M and M376A expressed in H9c2 cells promoted normal FCCP-stimulated Parkin recruitment to mitochondria, whereas Parkin aggregation was impaired in cardiomyoblasts expressing MFN2 R400Q (Figure 6D, Figure 6—figure supplement 1A); mitophagy measured as mitolysosome formation followed the same pattern (Figure 6E, Figure 6—figure supplement 1B). Mitophagy-defective MFN2 Q400 was less effective than WT MFN2 R400 in promoting H9c2 cardiomyoblast mitochondrial respiration (Figure 6—figure supplement 2), recapitulating its effects in MEFs (see Figure 2C). These findings further support the concept that the primary functional defect responsible for cardiac-specific abnormalities in MFN2 Q400 knock-in mice is impaired Parkin-mediated mitophagy.
Mfn2 R400Q suppresses compensatory mitophagy induced by doxorubicin and increases sensitivity to doxorubicin-induced cardiomyopathy
Taken together, the above results point to defective mitophagy as a central mechanism underlying murine perinatal cardiomyopathy induced by the Mfn2 R400Q mutation. Based on prior reports (Gong et al., 2015), developmentally regulated Mfn2-Parkin-mediated mitophagy is essential for myocardial metabolic maturation in the perinatal period, whereas in adult hearts the same mitophagy mechanism mediates mitochondrial quality control (Song et al., 2015; Song et al., 2017). We asked whether, despite having normal cardiac function under basal conditions, surviving Q/Q400 mice might have a residual mitophagy defect manifested only in the face of an extrinsic mitochondrial insult. Doxorubicin, an anti-cancer chemotherapeutic with well-characterized dose-dependent cardiotoxic effects mediated in large part by mitochondrial dysfunction (Wallace et al., 2020; Schirone et al., 2022), is such an insult. Therefore, we evaluated dose-dependent doxorubicin-induced cardiomyoblast toxicity in cultured H9c2 cells. Expression of MFN2 Q400 increased overall doxorubicin-stimulated H9c2 cell death and TUNEL apoptosis marker labeling (Figure 7A and B), while suppressing mitochondrial Parkin recruitment and mitolysosome formation (Figure 7C). Compared to WT MFN2 R400, expression of MFN2 Q400 increased baseline mitochondrial ROS to levels comparable to those provoked by doxorubicin in normal cells (Figure 7D).
Administration of doxorubicin to mice provokes mitochondrial injury that can lead to dilated cardiomyopathy (Dhingra et al., 2014). We compared the cardiac response to doxorubicin between young adult wild-type (R/R 400), heterozygous (R/Q 400), and surviving homozygous (Q/Q 400n) Mfn2 R400Q knock-in mice. Left ventricular (LV) chamber dimension at end-diastole (LVEDD) and calculated LV mass increased, and LV ejection performance (fractional shortening; FS) decreased to a greater extent in homozygous Mfn2 Q/Q 400 knock-in mice than WT or heterozygous mice (Figure 7E). Likewise, gravimetric heart weights were greater in doxorubicin-treated Mfn2 Q/Q400 mice than the other groups (Figure 7F). As in cultured H9c2 cardiomyoblasts, the exaggerated Mfn2 Q/Q400 mouse cardiotoxic response to doxorubicin was associated with diminished reactive mitophagy (Figure 7G).
Discussion
Our studies describe two major findings. First, we provide plausible evidence for a contribution of impaired MFN2-Parkin-mediated mitophagy to myocardial disease. Our conclusions are supported by a statistical genetic association in human cardiomyopathy, cardiac-specific phenotypes in knock-in mice expressing a candidate human mutation, and mechanistic studies defining the effects of the R400Q mutant MFN2 protein on mitophagy in cultured cardiomyoblasts and in vivo mouse hearts. Second, we describe a novel functional consequence of natural MFN2 mutation wherein an impaired ability to shift protein conformation leads to concomitant defects in both mitochondrial fusion and mitophagy, but a normal ability to promote mitochondrial transport. We surmise that the combination of mitophagy and fusion defects evokes cardiac disease, whereas combined mitochondrial motility and fusion defects described for many CMT2A-linked MFN2 mutations (including those studied herein) provoke peripheral neuropathy. This paradigm provides a mechanistic basis for tissue specificity of diseases evoked by MFN2 mutations.
Conventional hypertrophic and dilated cardiomyopathies are caused by mutations of sarcomeric proteins (Mazzarotto et al., 2020). The current data suggest that myocardial-specific clinical disease can also be induced by mutational disruption of mitophagy. As an initiator/effector of both mitochondrial fusion and Parkin-mediated mitophagy, MFN2 orchestrates cardiomyocyte fate in embryonic (Kasahara et al., 2013) and perinatal hearts (Gong et al., 2015). Under normal circumstances, mitofusin–mitofusin binding mediates mitochondrial fusion, and mitofusin–Parkin binding promotes mitophagy; these are mutually exclusive functions regulated by PINK1 kinase (Li et al., 2022). Thus, non-phosphorylated WT MFN2 promotes MFN-MFN-mediated fusion, but cannot bind Parkin. MFN2 phosphorylation on S378 extinguishes fusion (Rocha et al., 2018) and on T111 and S442 enables Parkin binding to initiate mitophagy (Chen and Dorn, 2013). MFN2 phosphorylation on these residues and its control of mitochondrial fate are mechanistically linked by altered MFN2 protein partnering as a consequence of phosphorylation-induced changes in MFN2 conformation (Dorn, 2020). Here, we considered that the location of the MFN2 R400Q and M376A mutations near PINK1-phosphorylatable MFN2 S378 and S442 might alter MFN2 conformational switching. Indeed, FRET studies revealed impaired transitioning of both MFN2 R400Q and M376A from the open/fusion-active conformation to the closed/fusion-inactive but mitophagy-permissive conformation. By comparison, GTPase-defective MFN2 mutants readily changed conformation. Our data do not indicate that these two HR2 domain mutations are permanently fixed in either conformational state, but that their transitioning from one conformation to the other is retarded.
Lack of conformational malleability did not distinguish between the cardiomyopathy- vs. the neuropathy-associated MFN2 mutants. Neither did impaired catalytic activity. While most CMT2A-associated MFN2 mutations occur in the GTPase domain, MFN2 M376A/V is in the HR1 domain and, like R400Q, has normal GTPase activity. Here, mitochondria expressing GTPase-defective MFN2 mutants were both fragmented (from impaired fusion) and depolarized (from degeneration of the mitochondrial respiratory chain). Consistent with respiratory impairment, Seahorse studies revealed severely depressed mitochondrial oxygen consumption in MEFs and H9c2 cells expressing GTPase-defective MFN2 mutants. By contrast, MFN2 HR1 domain mutants M376A and R400Q had normal GTPase activity; mitochondria expressing these mutants, while fragmented, were normally polarized and with only mildly impaired respiration compared to WT MFN2. These findings suggest that mitochondrial depolarization (rather than just fragmentation) may play an important pathological role in CMT2A, but not cardiomyopathy.
Our data provide a scientific foundation supporting consideration of human genetic mutations that suppress cardiomyocyte mitophagy as candidate risk modifiers for myocardial disease. For MFN2 R400Q, the genetic association between this mitophagy-modifying mutation and human cardiomyopathy is intriguing. However, the mutation is sufficiently rare that it is not possible from genetic data alone to define its contribution to clinical heart disease. While absence of MFN2 Q400 homozygosity in humans could reflect severe early cardiac pathology as manifested by our KI mice, it may simply be the consequence of its rarity. It is worth considering that the low prevalence of this mutation compared to non-synonymous MFN2 mutations in the same HR1 domain, and its over-representation in individuals of African descent, could also reflect negative selection throughout modern human evolution of an ancient and harmful gene variant. Consistent with this idea, the split phenotype of Mfn2 Q400 KI mice, that is, Q400a mice that succumbed to cardiomyopathy during the critical perinatal developmental period vs. Q400n mice that survived to adulthood with apparently normal hearts, belies the common underlying defect in cardiomyocyte mitophagy. We surmise that Q400 mice suffering from developmental mitophagic cardiomyopathy as late fetuses or early pups suffered from metabolic abnormalities and failed to upregulate ROS-managing enzymes. By comparison, Q400 mice that did not develop the spontaneous perinatal cardiomyopathy were better compensated against baseline ROS and metabolic derangements, yet remained hyper-sensitive to the extrinsic mitochondrial stress induced by doxorubicin because reactive mitophagy was impaired.
Finally, the Mfn2 Q400 KI mouse phenotype may provide fresh insights into the relative importance of mitochondrial fusion vs. mitophagy in hearts. We CRISPR engineered three human MFN2 mutations into the mouse genome. All three MFN2 mutations, T105M, M376A/V, and R400Q, exhibited impaired fusogenicity, but only MFN2 T105M provoked mitochondrial depolarization, abnormally low respiration, and impaired mitochondrial motility. This mutation induced early embryonic lethality when homozygous and caused a CMT2A-like neuropathy (to be described elsewhere) that spared the heart in heterozygous mice. Both M376A/V and R400Q were impaired in their ability to change conformation, explaining loss of function despite normal catalytic GTPase activity. However, only R400Q adversely affected Parkin localization to depolarized mitochondria, and therefore mitophagy, and of the three knock-in mouse lines only R400Q mice exhibited cardiac involvement.
In conclusion, the current findings point to MFN2 and potentially other mitophagy signaling pathways as candidate genetic modifiers of cardiomyopathy. In the context of MFN2 mutations, these data extend the paradigm that specific functional defects exhibited by mutant MFN2 proteins determine organ system disease involvement.
Materials and methods
Human study populations
Request a detailed protocolHeart failure (1117 Caucasian; 628 African American) and non-failing control subjects (625 Caucasian; 628 African American) from the Cincinnati Heart Study (Liggett et al., 2008; Matkovich et al., 2010; Cappola et al., 2011) and hypertrophic cardiomyopathy (240 Caucasian; 46 African American) had all MFN2 exons individually sequenced to identify DNA sequence variants. Additionally, 424 individuals (281 probands) with FCM from the DCM Precision Medicine Study (Huggins et al., 2022; Trachtenberg et al., 2022) were screened to detect any MFN2 mutations suitable for linkage analysis; none were identified. Overall and race-specific MFN2 allele frequencies were compared to non-failing controls from the Cincinnati Heart Study and to the gnomAD data set of 141,456 individuals from various racial backgrounds.
Cell lines
Request a detailed protocolMfn1/Mfn2 double null MEFs fibroblasts were purchased from American Type Culture Collection (ATCC Manassas, VA, CRL-2994). Normal MEFs were prepared by enzymatic dissociation from embryonic day E.13.5–14.5 C57BL/6J mice (The Jackson Laboratory, Cat# 000664). H9c2 rat cardiomyoblasts were purchased from American Type Culture Collection (ATCC, CRL-1446). MEFs and H9c2 cells were cultured at 37°C, 5% CO2-95% air in Dulbecco’s minimal essential medium (DMEM) containing glucose (4.5 g/l) (Thermo Fisher Scientific, 11965-084) with 10% (v/v) fetal bovine serum (FBS; Gibco, Gaithersburg, MD, Cat# 26140-079), 1× nonessential amino acids (Gibco, Cat# 11130051), 2 mM L-glutamine (Corning, NY, Cat# 34717007), and 1% (v/v) penicillin/streptomycin (Gibco, Cat# 15140-122). All cells were tested negative for mycoplasma. Adult mouse DRG neurons were prepared from 8- to 12-week-old C57BL/6J using enzymatic dissociation.
Viral vectors
Request a detailed protocolMFN2 mutants and FRET probes were generated using PCR mutagenesis, sequence verified, and sent to Vector BioLabs for custom adenovirus production. Ad-CMV-β-Gal was purchased from Vector BioLabs (Cat# 1080).
Protein structure modeling
Request a detailed protocolHypothetical structures of human MFN2 were generated using I-TASSER and Chimera UCSF. The putative closed conformation is based on structural homology with bacterial dynamin-like protein (PDB: 2J69). The Mfn2 open conformation was downloaded from the AlphaFold protein structure database (AF-095140-F1). Domain coloring is as follows: green GTPase (AA 94-265); red HR1 (AA 338-421); and blue HR-2 (AA 681-757).
Imaging
Request a detailed protocolLive cell studies of MEFs, H9c2 cells, or DRGs to assess mitochondrial morphology, respiratory function, mitophagy, and motility, and transmission electron microscopy studies were performed as previously described (Franco et al., 2016; Chen and Dorn, 2013; Misko et al., 2010 and Rocha et al., 2018). Briefly, confocal fluorescence microscopy was used to gauge mitochondrial aspect ratio (length/width) quantified based on MitoTracker Green staining using ImageJ; mitochondrial depolarization quantified as the number of green/total mitochondria in MitoTracker Green/TMRE double-stained images; mitochondrial Parkin aggregation quantified in cells expressing mcParkin with or without MitoTracker green; mitolysosome formation quantified as number of red dots per cell in Mito QC-stained images; mitochondrial engulfment by autophagosomes quantified as the ratio of yellow to total cells (green/red overlay) in cells expressing LC3-GFP and co-stained with MitoTracker Orange. Mitochondrial motility was assessed in processes of cultured mouse DRG neurons expressing mitochondrial targeted RFP using time-lapse confocal microscopy (1 frame/5 s).
Whole mouse embryos or adult mouse hearts were fixed with 4% formaldehyde solution in PBS, embedded in paraffin, and stained with hematoxylin-eosin (Sigma, GHS116 and HT110116) or TUNEL labeled using the manufacturer’s protocol (Promega, G3250).
Transmission electron microscopy of mouse heart tissues fixed in 4% paraformaldehyde, 2.5% glutaraldehyde in 0.1 M sodium cacodylate buffer, pH 7.4, used thin sections stained with osmium tetroxide/uranyl acetate. A JEOL electron microscope (JEM-1400, JEOL, Tokyo, Japan) at 1500×–5000× direct magnification was used for image acquisition. Mitochondrial size and surface density were quantified using ImageJ (NIH).
Western blotting
Request a detailed protocolSeventy-two hours after adenoviral transduction (MOI 50) with MFN2 mutants, Mfn1/Mfn2 DKO MEFs and H9c2 were pelleted and lysed in cell extraction buffer on ice (Thermo Fisher Scientific, Cat# FNN0011) with 1 mM PMSF, protease inhibitor (Roche, Cat# 05892970001) and phosphatase inhibitor (Roche, Cat# 04906837001). Primary antibodies used were anti-Mfn2 (1:500, Abcam ab56889) and anti-β-actin (1:3000, Proteintech, Cat# 66009-1). Horseradish peroxidase (HRP)-conjugated anti-mouse IgG (1:3000, cs7076) was from Cell Signaling Technology.
GTPase activity assay
Request a detailed protocolMfn1/Mfn2 double-null MEFs were transduced with ad-MFN2 mutants (MOI 50). Mitochondria were isolated as described (Dang et al., 2022; Dang et al., 2020; Detmer and Chan, 2007; Dhingra et al., 2014), prepared on ice, and used fresh. Then, 100 μg of mitochondrial protein in triplicate was incubated in GTPase Buffer (Promega #V7681; Madison, WI) with 10 μM GTP and 1 mM DTT. Also, 1 μM MiM 111, 1 μM Chimera or DMSO (vehicle) were added as indicated and the reactants incubated at room temperature for 90 min in 96-well plates. The Promega GTPase-Glo assay kit was used to measure GTPase activity following the manufacturer’s instructions. Luminescence was quantified on a Promega GloMax Luminometer.
MFN2 FRET assays
Request a detailed protocolMfn1/Mfn2 double-null MEFs were transduced with FRET ad-MFNs. Then, 72 hr later mitochondria were isolated and FRET measured in a 96-well format as described (Dang et al., 2020).
Mitochondrial respiration
Request a detailed protocolOxygen consumption rate (OCR) of mitofusin mutants expressed using adenoviral transduction of Mfn DKO MEFs or H9c2 cardiomyoblasts used a Seahorse XFe24 Extracellular Flux Analyzer (Seahorse Bioscience, Billerica, MA). Briefly, cells were plated on Seahorse XF24-well cell culture microplates for OCR measurements 72 hr after adenoviral transduction. Sensor cartridge hydration was performed overnight at 37°C. After basal respiration measurements, the following were autoinjected in sequence: 1 μM oligomycin to inhibit ATP synthase, 1 μM FCCP to uncouple oxidative phosphorylation, and 0.5 μM rotenone/antimycin A to abrogate electron transport (non-mitochondrial OCR). Each experiment averaged five or more replicate wells, and each experiment was repeated with a minimum of four biological replicates. ATP-linked respiration is basal OCR–OCR after oligomycin injection. Maximal mitochondrial respiration is OCR after FCCP–OCR after rotenone/antimycin A.
Flow cytometry of mitochondrial ROS production
Request a detailed protocolH9C2 cardiomyoblasts were transduced with adenoviral MFN2 constructs as described. Then, 72 hr later cells were co-stained with MitoSOX Red (5 μM for 10 min at 37°C) and MitoTracker Green FM 200 nM (Thermo Fisher, Cat# M36008 and #M7514, respectively). Fluorescence was analyzed on a ZE5 Flow Cytometer (Bio-Rad). Data were analyzed with FlowJo software version 10 and are presented as mean fluorescence intensity of five independent experiments.
Creation of Mfn2 Q400 knock-in mice using Crisp/Cas9
Request a detailed protocolThe knock-in strategy is schematically depicted in Figure 4—figure supplement 1A. CRISPR gRNAs for in vitro testing were identified using http://crispr.mit.edu/. Mfn2 sgRNA was cloned into BbsI digested plasmid pX330 (Addgene # 42230). sgRNA activity was validated in vitro by transfection of NIH3T3 cells using ROCHE Xtremegene HP. Cell pools were harvested 48 hr later for genomic DNA prep, followed by Sanger sequencing of products spanning the gRNA/Cas9 cleavage site and TIDE analysis (https://tide.nki.nl/) of sequence trace files. T7 sgRNA template was PCR-amplified, gel-purified, and in vitro-transcribed with the MEGAshortscript T7 kit (Life Technologies). T7 Cas9 template was PCR-amplified, gel-purified, and in vitro-transcribed with the T7 mMessage mMachine Ultra kit (Life Technologies). After transcription, RNA was purified with Megaclear kit (Life Technologies). A 200 nt ssODN donor DNA with the mutation centered within the oligo was ordered from IDT as an ultramer oligo. Injection concentrations were 50 ng/μl Cas9, 25 ng/μl gRNA, and 20 ng/μl ssODN. Founders were identified using QIAGEN pyrosequencer and Pyromark Q96 2.5.7 software. All animal experiments were approved by institutional IACUC protocols. B6/CBA F1 mice at 3–4 wk of age (JAX Laboratories, Bar Harbor, ME) were superovulated by intraperitoneal injection of 5 IU pregnant mare serum gonadotropin, followed 48 hr later by intraperitoneal injection of 5 IU human chorionic gonadotropin (PMS from Sigma, HGC from Millipore USA). Mouse zygotes were obtained by mating B6/CBA stud males with superovulated B6/CBA females at a 1:1 ratio. Single-cell fertilized embryos were injected into the pronucleus and cytoplasm of each zygote using standard techniques.
gRNA sequence: 5′ AAGAGCGGCAAGACCgGCTG 3
Antisense ssODN sequence: 5′cttcattctcacCTGCCTTTCCACTTCCTCCGTAATCTGCTTAATTCGCAGCTTGTAGTCTTGAGCCAGGAGCTCCAGCTGCTTGTCAATAAACCTCAGCtGGTCTTGCCGCTCTTCCCGCATTTCTAGGCAATAAACCctgagaggacaaaagcactgcttagaaacccgtgtccccacaagagcaagcagagaagcca 3′
Creation of Mfn2 T105M and M376V knock-in mice
Request a detailed protocolCRISPR engineering was also performed at Cyagen, Inc to insert two CMT2A-associated mutations into the mouse genome. The overall procedures are similar to that described above.
For Mfn2 T105M: gRNA (forward strand) sequence was 5′ TTTCATCCTAGGACGAGCAATGG 3′. The p.T105M (ACG to ATG) in the donor oligo (see Figure 4—figure supplement 1B) was introduced into exon 6 by homology-directed repair. Two synonymous mutations, p.S106 (AGC to TCA) and p.S110 (AGC to TCG), were introduced simultaneously to prevent the binding and recutting of the sequence by gRNA after homology-directed repair; the SalI restriction site was used for genotyping.
For Mfn2 M376V: gRNA (forward strand) sequence was 5′ CTCCTGAGCTGCGATGTGCAGGG 3′. The p.M376V (ATG to GTG) in the donor oligo (see Figure 4—figure supplement 1C) was introduced into exon 12 by homology-directed repair. One synonymous mutation p.S378 (TCC to TCA) was also introduced to prevent the binding and recutting of the sequence by gRNA after homology-directed repair.
Pups were genotyped and correct integration of the fragment carrying the point mutation (p.T105M) or (p.M376V) was confirmed by Sanger sequencing. The correctly integrated mutant F0 founder mice were back-crossed with WT C57BL/6J mice to generate F1s. F1 heterozygous animals, also confirmed by Sanger sequencing, were interbred to build mouse colonies. All subsequent generation knock-mouse genotypes for all three lines were determined by Sanger sequencing.
Omics studies
Request a detailed protocolTotal RNA from Mfn2Q/Q400 and C57BL/6J E.18.5 or P1 hearts was isolated using TRIzol Reagent (Cat# 15596-026, Ambion RNA Life Technologies) and treated with DNase I (Invitrogen, Cat# 18068015). mRNA library preparation and sequencing were performed at the Genome Technology Access Center (GTAC) of Washington University St. Louis, MO (USA) using the RiboErase protocol with indexing and pooled sequencing on a NovaSeq S4. Sequencing depth was 30M reads per sample. RNA-seq reads were aligned to the Mus musculus mm-10 assembly with STAR. Refseq was used as annotation model to quantify gene counts normalized to counts per million (CPM). Pair-wise comparisons between transcriptomic profiles from each experimental group were analyzed using GSA in Partek Flow and filtered a priori CPM ≥ 10 for at least one biological sample, FDR < 0.05, p-value<0.02, fold change >2, <-2. Unsupervised hierarchical clustering using Euclidean distance and average distance was performed using Partek Flow. KEGG Pathway analyses used Database for Annotation, Visualization and Integrated Discovery (DAVID v6.8). The GEO accession number is GSE214984.
Metabolomics analyses of acylcarnitines, organic acids, and amino acids in P0 hearts were performed as described (Gong et al., 2015).
Doxorubicin experiments
Request a detailed protocolIn vitro – H9c2 cells were transduced with adenovirus encoding MFN2 R400 (WT) or Q400 at 100 MOI 48 hr prior to addition of doxorubicin; immunoblotting confirmed equal levels of MFN2 expression. MFN2-expressing cells were cultured in the presence of increasing doxorubicin concentrations (0.1–50 uM), or in the concentration specified, for 1, 12, or 24 hr as indicated. Live-dead staining used the kit from Invitrogen (Cat# L3224). TUNEL labeling used the DeadEnd Fluorometric kit from Promega (Cat# G3250). Mito QC mitolysosome labeling and mc Parkin aggregation studies are described above.
In vivo – 8–10-week-old male WT (R/R400), heterozygous (R/Q400 knock-in), and homozygous (Q/Q400 knock-in) mice underwent baseline echocardiographic assessments prior to receiving a single dose of doxorubicin (20 ug/g via intraperitoneal injection; Dhingra et al., 2014). Four days later, each mouse was administered adenovirus-Mito QC (3.3 × 1011 particles in 100 ul). Echocardiography was performed 4 and 7 d post doxorubicin; mice were terminated after the 7-day echo study. Hearts were dissected, removed, placed in ice-cold physiological saline (0.9% NaCl), and weighed; heart weights were indexed to tibial length to account for individual variations in weight loss induced by doxorubicin. Mito QC staining was evaluated by confocal fluorescence imaging in 12 um frozen sections of OCT-embedded specimens.
Data presentation and statistical analyses
Request a detailed protocolData are reported as means ± SEM. Sample number (n) indicates the number of independent biological samples. Two-group comparisons used Student’s t-test; multiple group comparisons used one-way ANOVA with Tukey’s post hoc test for individual statistical comparisons. Survival was evaluated by Kaplan–Meier analysis and the log-rank (Mantel–Cox) test. Comparisons of population MAF used chi-square testing. p<0.05 was considered statistically significant.
Data availability
Sequencing Data have been deposited in GEO under accession code GSE214984. - The putative closed conformation is based on structural homology with bacterial dynamin-like protein (PDB: 2J69). The Mfn2 open conformation was downloaded from the AlphaFold protein structure database (AF-095140-F1).
-
NCBI Gene Expression OmnibusID GSE214984. A human mitofusin 2 mutation can cause mitophagic cardiomyopathy.
References
-
Mitochondrial fusion and fission in mammalsAnnual Review of Cell and Developmental Biology 22:79–99.https://doi.org/10.1146/annurev.cellbio.22.010305.104638
-
Mitofusins Mfn1 and Mfn2 coordinately regulate mitochondrial fusion and are essential for embryonic developmentThe Journal of Cell Biology 160:189–200.https://doi.org/10.1083/jcb.200211046
-
Mitochondrial fusion is essential for organelle function and cardiac homeostasisCirculation Research 109:1327–1331.https://doi.org/10.1161/CIRCRESAHA.111.258723
-
Complementation between mouse Mfn1 and Mfn2 protects mitochondrial fusion defects caused by CMT2A disease mutationsThe Journal of Cell Biology 176:405–414.https://doi.org/10.1083/jcb.200611080
-
Parkin-dependent mitophagy in the heartJournal of Molecular and Cellular Cardiology 95:42–49.https://doi.org/10.1016/j.yjmcc.2015.11.023
-
Evolving concepts of mitochondrial dynamicsAnnual Review of Physiology 81:1–17.https://doi.org/10.1146/annurev-physiol-020518-114358
-
Mitofusins as mitochondrial anchors and tethersJournal of Molecular and Cellular Cardiology 142:146–153.https://doi.org/10.1016/j.yjmcc.2020.04.016
-
Mitofusin activation enhances mitochondrial motility and promotes neuroregeneration in CMT2ANeural Regeneration Research 16:2201–2203.https://doi.org/10.4103/1673-5374.310684
-
Loss-of-function mutations in SLC30A8 protect against type 2 diabetesNature Genetics 46:357–363.https://doi.org/10.1038/ng.2915
-
Mitochondrial Dysfunction and Pharmacodynamics of Mitofusin Activation in Murine Charcot-Marie-Tooth Disease Type 2AThe Journal of Pharmacology and Experimental Therapeutics 383:137–148.https://doi.org/10.1124/jpet.122.001332
-
Reciprocal Regulation of Mitofusin 2-Mediated Mitophagy and Mitochondrial Fusion by Different PINK1 Phosphorylation EventsFrontiers in Cell and Developmental Biology 10:868465.https://doi.org/10.3389/fcell.2022.868465
-
Developmental changes in energy substrate use by the heartCardiovascular Research 26:1172–1180.https://doi.org/10.1093/cvr/26.12.1172
-
Overview of mitochondrial bioenergeticsMethods in Molecular Biology 1782:1–6.https://doi.org/10.1007/978-1-4939-7831-1_1
-
Association of an intronic, but not any exonic, FRMD4B sequence variant and heart failureClinical and Translational Science 3:134–139.https://doi.org/10.1111/j.1752-8062.2010.00220.x
-
Mitofusin 2 is necessary for transport of axonal mitochondria and interacts with the Miro/Milton complexThe Journal of Neuroscience 30:4232–4240.https://doi.org/10.1523/JNEUROSCI.6248-09.2010
-
Parkin is recruited selectively to impaired mitochondria and promotes their autophagyThe Journal of Cell Biology 183:795–803.https://doi.org/10.1083/jcb.200809125
-
Myozenin 2 is a novel gene for human hypertrophic cardiomyopathyCirculation Research 100:766–768.https://doi.org/10.1161/01.RES.0000263008.66799.aa
-
Maturation of cardiac energy metabolism during perinatal developmentFrontiers in Physiology 9:959.https://doi.org/10.3389/fphys.2018.00959
-
Control of mitochondrial morphology by a human mitofusinJournal of Cell Science 114:867–874.https://doi.org/10.1242/jcs.114.5.867
-
Loss-of-function variants in ABCA7 confer risk of Alzheimer’s diseaseNature Genetics 47:445–447.https://doi.org/10.1038/ng.3246
-
MFN2-related neuropathies: Clinical features, molecular pathogenesis and therapeutic perspectivesJournal of the Neurological Sciences 356:7–18.https://doi.org/10.1016/j.jns.2015.05.033
-
Mitochondrial determinants of doxorubicin-induced cardiomyopathyCirculation Research 126:926–941.https://doi.org/10.1161/CIRCRESAHA.119.314681
-
New methods for monitoring mitochondrial biogenesis and mitophagy in vitro and in vivoExperimental Biology and Medicine 242:781–787.https://doi.org/10.1177/1535370216688802
-
Restoring mitofusin balance prevents axonal degeneration in a Charcot-Marie-Tooth type 2A modelThe Journal of Clinical Investigation 131:e147307.https://doi.org/10.1172/JCI147307
-
Axonal neuropathy with optic atrophy is caused by mutations in mitofusin 2Annals of Neurology 59:276–281.https://doi.org/10.1002/ana.20797
Article and author information
Author details
Funding
National Institutes of Health (NIH R35 HL135736)
- Gerald W Dorn
National Institutes of Health (R42 NS113642)
- Gerald W Dorn
National Institutes of Health (R42 NS115184)
- Gerald W Dorn
National Institutes of Health (R01 HL128349)
- Daniel P Kelly
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
GWD conceived of and designed the research. GWD and AF wrote the manuscript. AF and JL performed mitochondrial function, GTPase and FRET assays, and knock-in mouse studies. RML and MS generated the CRISPR/Cas9 knock-in mouse. DPK and the Sanford-Burnham Medical Discovery Institute Metabolomics Core performed and analyzed the metabolomics studies. CdGS, MEM, ADS, AF JE, and XD analyzed RNAseq or human MFN2 genotype data. RH, AJM, and GWD provided human MFN2 sequencing from cardiomyopathy research cohorts. JE analyzed RNAseq in revised manuscript. All authors agree with their inclusion and place in the author list.
Supported by NIH R35 HL135736, R42 NS115184 and R42 NS113642 to GWD, R01 HL128349 to DPK, and by the Hope Center Transgenic Vectors Core at Washington University School of Medicine.
Ethics
All experimental procedures were approved by the Washington University in St. Louis, School of Medicine Animal Studies Committee; IACUC protocol number 22-0314 Exp:12/27/2025, .PI Gerald Dorn.
Copyright
© 2023, Franco et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 816
- views
-
- 142
- downloads
-
- 8
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Developmental Biology
Premature infants with bronchopulmonary dysplasia (BPD) have impaired alveolar gas exchange due to alveolar simplification and dysmorphic pulmonary vasculature. Advances in clinical care have improved survival for infants with BPD, but the overall incidence of BPD remains unchanged because we lack specific therapies to prevent this disease. Recent work has suggested a role for increased transforming growth factor-beta (TGFβ) signaling and myofibroblast populations in BPD pathogenesis, but the functional significance of each remains unclear. Here, we utilize multiple murine models of alveolar simplification and comparative single-cell RNA sequencing to identify shared mechanisms that could contribute to BPD pathogenesis. Single-cell RNA sequencing reveals a profound loss of myofibroblasts in two models of BPD and identifies gene expression signatures of increased TGFβ signaling, cell cycle arrest, and impaired proliferation in myofibroblasts. Using pharmacologic and genetic approaches, we find no evidence that increased TGFβ signaling in the lung mesenchyme contributes to alveolar simplification. In contrast, this is likely a failed compensatory response, since none of our approaches to inhibit TGFβ signaling protect mice from alveolar simplification due to hyperoxia while several make simplification worse. In contrast, we find that impaired myofibroblast proliferation is a central feature in several murine models of BPD, and we show that inhibiting myofibroblast proliferation is sufficient to cause pathologic alveolar simplification. Our results underscore the importance of impaired myofibroblast proliferation as a central feature of alveolar simplification and suggest that efforts to reverse this process could have therapeutic value in BPD.
-
- Developmental Biology
Cells called alveolar myofibroblasts, which have a central role in the development of the lung after birth, receive an orchestrated input from a range of different signaling pathways.