Genome-Wide Association Mapping Reveals That Specific and Pleiotropic Regulatory Mechanisms Fine-Tune Central Metabolism and Growth in Arabidopsis

Genome-wide association studies identify loci controlling specialized seed metabolites in Arabidopsis

Plants synthesize specialized metabolites to facilitate environmental and ecological interactions. During evolution, plants diversified in their potential to synthesize these metabolites. Quantitative differences in metabolite levels of natural Arabidopsis (Arabidopsis thaliana) accessions can be employed to unravel the genetic basis for metabolic traits using genome-wide association studies (GWAS). Here, we performed metabolic GWAS on seeds of a panel of 315 A. thaliana natural accessions, including the reference genotypes C24 and Col-0, for polar and semi-polar seed metabolites using untargeted ultra-performance liquid chromatography-mass spectrometry. As a complementary approach, we performed quantitative trait locus (QTL) mapping of near-isogenic introgression lines between C24 and Col-0 for specific seed specialized metabolites. Besides common QTL between seeds and leaves, GWAS revealed seed-specific QTL for specialized metabolites, indicating differences in the genetic architecture of seeds and leaves. In seeds, aliphatic methylsulfinylalkyl and methylthioalkyl glucosinolates associated with the ALKENYL HYDROXYALKYL PRODUCING loci (GS-ALK and GS-OHP) on chromosome 4 containing alkenyl hydroxyalkyl producing 2 (AOP2) and 3 (AOP3) or with the GS-ELONG locus on chromosome 5 containing methylthioalkyl malate synthase (MAM1) and MAM3. We detected two unknown sulfur-containing compounds that were also mapped to these loci. In GWAS, some of the annotated flavonoids (kaempferol 3-O-rhamnoside-7-O-rhamnoside, quercetin 3-O-rhamnoside-7-O-rhamnoside) were mapped to transparent testa 7 (AT5G07990), encoding a cytochrome P450 75B1 monooxygenase. Three additional mass signals corresponding to quercetin-containing flavonols were mapped to UGT78D2 (AT5G17050). The association of the loci and associating metabolic features were functionally verified in knockdown mutant lines. By performing GWAS and QTL mapping, we were able to leverage variation of natural populations and parental lines to study seed specialized metabolism. The GWAS data set generated here is a high-quality resource that can be investigated in further studies.

Natural variation in the SVP contributes to the pleiotropic adaption of Arabidopsis thaliana across contrasted habitats

Coordinated plant adaptation involves the interplay of multiple traits driven by habitat-specific selection pressures. Pleiotropic effects, wherein genetic variants of a single gene control multiple traits, can expedite such adaptations. Until present, only a limited number of genes have been reported to exhibit pleiotropy. Here, we create a recombinant inbred line (RIL) population derived from two Arabidopsis thaliana (A. thaliana) ecotypes originating from divergent habitats. Using this RIL population, we identify an allelic variation in a MADS-box transcription factor, SHORT VEGETATIVE PHASE (SVP), which exerts a pleiotropic effect on leaf size and drought-versus-humidity tolerance. Further investigation reveals that a natural null variant of the SVP protein disrupts its normal regulatory interactions with target genes, including GRF3, CYP707A1/3, and AtBG1, leading to increased leaf size, enhanced tolerance to humid conditions, and changes in flowering time of humid conditions in A. thaliana. Remarkably, polymorphic variations in this gene have been traced back to early A. thaliana populations, providing a genetic foundation and plasticity for subsequent colonization of diverse habitats by influencing multiple traits. These findings advance our understanding of how plants rapidly adapt to changing environments by virtue of the pleiotropic effects of individual genes on multiple trait alterations.

Hello darkness, my old friend: 3-KETOACYL-COENZYME A SYNTHASE4 is a branch point in the regulation of triacylglycerol synthesis in Arabidopsis thaliana

Plant lipids are important as alternative sources of carbon and energy when sugars or starch are limited. Here, we applied combined heat and darkness or extended darkness to a panel of ∼300 Arabidopsis (Arabidopsis thaliana) accessions to study lipid remodeling under carbon starvation. Natural allelic variation at 3-KETOACYL-COENZYME A SYNTHASE4 (KCS4), a gene encoding an enzyme involved in very long chain fatty acid (VLCFA) synthesis, underlies the differential accumulation of polyunsaturated triacylglycerols (puTAGs) under stress. Ectopic expression of KCS4 in yeast and plants proved that KCS4 is a functional enzyme localized in the endoplasmic reticulum with specificity for C22 and C24 saturated acyl-CoA. Allelic mutants and transient overexpression in planta revealed the differential role of KCS4 alleles in VLCFA synthesis and leaf wax coverage, puTAG accumulation, and biomass. Moreover, the region harboring KCS4 is under high selective pressure and allelic variation at KCS4 correlates with environmental parameters from the locales of Arabidopsis accessions. Our results provide evidence that KCS4 plays a decisive role in the subsequent fate of fatty acids released from chloroplast membrane lipids under carbon starvation. This work sheds light on both plant response mechanisms and the evolutionary events shaping the lipidome under carbon starvation.

Natural variation of respiration-related traits in plants

Plant respiration is one of the greatest global metabolic fluxes, but rates of respiration vary massively both within different cell types as well as between different individuals and different species. Whilst this is well known, few studies have detailed population-level variation of respiration until recently. The last 20 years have seen a renaissance in studies of natural variance. In this review, we describe how experimental breeding populations and collections of large populations of accessions can be used to determine the genetic architecture of plant traits. We further detail how these approaches have been used to study the rate of respiration per se as well as traits that are intimately associated with respiration. The review highlights specific breakthroughs in these areas but also concludes that the approach should be more widely adopted in the study of respiration per se as opposed to the more frequently studied respiration-related traits.

cis-prenyltransferase 3 and α/β-hydrolase are new determinants of dolichol accumulation in Arabidopsis

Dolichols (Dols), ubiquitous components of living organisms, are indispensable for cell survival. In plants, as well as other eukaryotes, Dols are crucial for post-translational protein glycosylation, aberration of which leads to fatal metabolic disorders in humans and male sterility in plants. Until now, the mechanisms underlying Dol accumulation remain elusive. In this study, we have analysed the natural variation of the accumulation of Dols and six other isoprenoids among more than 120 Arabidopsis thaliana accessions. Subsequently, by combining QTL and GWAS approaches, we have identified several candidate genes involved in the accumulation of Dols, polyprenols, plastoquinone and phytosterols. The role of two genes implicated in the accumulation of major Dols in Arabidopsis-the AT2G17570 gene encoding a long searched for cis-prenyltransferase (CPT3) and the AT1G52460 gene encoding an α/β-hydrolase-is experimentally confirmed. These data will help to generate Dol-enriched plants which might serve as a remedy for Dol-deficiency in humans.

Genome-wide association of the metabolic shifts underpinning dark-induced senescence in Arabidopsis

Dark-induced senescence provokes profound metabolic shifts to recycle nutrients and to guarantee plant survival. To date, research on these processes has largely focused on characterizing mutants deficient in individual pathways. Here, we adopted a time-resolved genome-wide association-based approach to characterize dark-induced senescence by evaluating the photochemical efficiency and content of primary and lipid metabolites at the beginning, or after 3 or 6 days in darkness. We discovered six patterns of metabolic shifts and identified 215 associations with 81 candidate genes being involved in this process. Among these associations, we validated the roles of four genes associated with glycine, galactinol, threonine, and ornithine levels. We also demonstrated the function of threonine and galactinol catabolism during dark-induced senescence. Intriguingly, we determined that the association between tyrosine contents and TYROSINE AMINOTRANSFERASE 1 influences enzyme activity of the encoded protein and transcriptional activity of the gene under normal and dark conditions, respectively. Moreover, the single-nucleotide polymorphisms affecting the expression of THREONINE ALDOLASE 1 and the amino acid transporter gene AVT1B, respectively, only underlie the variation in threonine and glycine levels in the dark. Taken together, these results allow us to present a very detailed model of the metabolic aspects of dark-induced senescence, as well as the process itself.

Genome-Wide Association Study Reveals Key Genes for Differential Lead Accumulation and Tolerance in Natural Arabidopsis thaliana Accessions

Soil contamination by lead (Pb) has become one of the major ecological threats to the environment. Understanding the mechanisms of Pb transport and deposition in plants is of great importance to achieve a global Pb reduction. We exposed a collection of 360 Arabidopsis thaliana natural accessions to a Pb-polluted soil. Germination rates, growth, and leaf Pb concentrations showed extensive variation among accessions. These phenotypic data were subjected to genome wide association studies (GWAs) and we found a significant association on chromosome 1 for low leaf Pb accumulation. Genes associated with significant SNP markers were evaluated and we selected EXTENSIN18 (EXT18) and TLC (TRAM-LAG1-CLN8) as candidates for having a role in Pb homeostasis. Six Pb-tolerant accessions, three of them exhibiting low leaf Pb content, and three of them with high leaf Pb content; two Pb-sensitive accessions; two knockout T-DNA lines of GWAs candidate genes (ext18, tlc); and Col-0 were screened under control and high-Pb conditions. The relative expression of EXT18, TLC, and other genes described for being involved in Pb tolerance was also evaluated. Analysis of Darwinian fitness, root and leaf ionome, and TEM images revealed that Pb-tolerant accessions employ two opposing strategies: (1) low translocation of Pb and its accumulation into root cell walls and vacuoles, or (2) high translocation of Pb and its efflux to inactive organelles or intracellular spaces. Plants using the first strategy exhibited higher expression of EXT18 and HMA3, thicker root cell walls and Pb vacuolar sequestration, suggesting that these genes may contribute to the deposition of Pb in the roots. On the other hand, plants translocating high amounts of Pb showed upregulation of TLC and ABC transporters, indicating that these plants were able to properly efflux Pb in the aerial tissues. We conclude that EXT18 and TLC upregulation enhances Pb tolerance promoting its sequestration: EXT18 favors the thickening of the cell walls improving Pb accumulation in roots and decreasing its toxicity, while TLC facilitates the formation of dictyosome vesicles and the Pb encapsulation in leaves. These findings are relevant for the design of phytoremediation strategies and environment restoration.

Genome-wide association studies: assessing trait characteristics in model and crop plants

GWAS involves testing genetic variants across the genomes of many individuals of a population to identify genotype–phenotype association. It was initially developed and has proven highly successful in human disease genetics. In plants genome-wide association studies (GWAS) initially focused on single feature polymorphism and recombination and linkage disequilibrium but has now been embraced by a plethora of different disciplines with several thousand studies being published in model and crop species within the last decade or so. Here we will provide a comprehensive review of these studies providing cases studies on biotic resistance, abiotic tolerance, yield associated traits, and metabolic composition. We also detail current strategies of candidate gene validation as well as the functional study of haplotypes. Furthermore, we provide a critical evaluation of the GWAS strategy and its alternatives as well as future perspectives that are emerging with the emergence of pan-genomic datasets.

Improving Genomic Prediction Using High-Dimensional Secondary Phenotypes

In the past decades, genomic prediction has had a large impact on plant breeding. Given the current advances of high-throughput phenotyping and sequencing technologies, it is increasingly common to observe a large number of traits, in addition to the target trait of interest. This raises the important question whether these additional or “secondary” traits can be used to improve genomic prediction for the target trait. With only a small number of secondary traits, this is known to be the case, given sufficiently high heritabilities and genetic correlations. Here we focus on the more challenging situation with a large number of secondary traits, which is increasingly common since the arrival of high-throughput phenotyping. In this case, secondary traits are usually incorporated through additional relatedness matrices. This approach is however infeasible when secondary traits are not measured on the test set, and cannot distinguish between genetic and non-genetic correlations. An alternative direction is to extend the classical selection indices using penalized regression. So far, penalized selection indices have not been applied in a genomic prediction setting, and require plot-level data in order to reliably estimate genetic correlations. Here we aim to overcome these limitations, using two novel approaches. Our first approach relies on a dimension reduction of the secondary traits, using either penalized regression or random forests (LS-BLUP/RF-BLUP). We then compute the bivariate GBLUP with the dimension reduction as secondary trait. For simulated data (with available plot-level data), we also use bivariate GBLUP with the penalized selection index as secondary trait (SI-BLUP). In our second approach (GM-BLUP), we follow existing multi-kernel methods but replace secondary traits by their genomic predictions, with the advantage that genomic prediction is also possible when secondary traits are only measured on the training set. For most of our simulated data, SI-BLUP was most accurate, often closely followed by RF-BLUP or LS-BLUP. In real datasets, involving metabolites in Arabidopsis and transcriptomics in maize, no method could substantially improve over univariate prediction when secondary traits were only available on the training set. LS-BLUP and RF-BLUP were most accurate when secondary traits were available also for the test set.

Characterization of effects of genetic variants via genome-scale metabolic modelling

Genome-scale metabolic networks for model plants and crops in combination with approaches from the constraint-based modelling framework have been used to predict metabolic traits and design metabolic engineering strategies for their manipulation. With the advances in technologies to generate large-scale genotyping data from natural diversity panels and other populations, genome-wide association and genomic selection have emerged as statistical approaches to determine genetic variants associated with and predictive of traits. Here, we review recent advances in constraint-based approaches that integrate genetic variants in genome-scale metabolic models to characterize their effects on reaction fluxes. Since some of these approaches have been applied in organisms other than plants, we provide a critical assessment of their applicability particularly in crops. In addition, we further dissect the inferred effects of genetic variants with respect to reaction rate constants, abundances of enzymes, and concentrations of metabolites, as main determinants of reaction fluxes and relate them with their combined effects on complex traits, like growth. Through this systematic review, we also provide a roadmap for future research to increase the predictive power of statistical approaches by coupling them with mechanistic models of metabolism.

Dissection of the domestication-shaped genetic architecture of lettuce primary metabolism

Lettuce (Lactuca sativa L.) is an important vegetable crop species worldwide. The primary metabolism of this species is essential for its growth, development and reproduction as well as providing a considerable direct source of energy and nutrition for humans. Here, through investigating 77 primary metabolites in 189 accessions including all major horticultural types and wild lettuce L. serriola we showed that the metabolites in L. serriola were different from those in cultivated lettuce. The findings were consistent with the demographic model of lettuce and supported a single domestication event for this species. Selection signals among these metabolic traits were detected. Specifically, galactinol, malate, quinate and threonate were significantly affected by the domestication process and cultivar differentiation of lettuce. Galactinol and raffinose might have been selected during stem lettuce cultivation as an adaption to the local environments in China. Furthermore, we identified 154 loci significantly associated with the level of 51 primary metabolites. Three genes (LG8749721, LG8763094 and LG5482522) responsible for the levels of galactinol, raffinose, quinate and chlorogenic acid were further dissected, which may have been the target of domestication and/or affected by local adaptation. Additionally, our findings strongly suggest that human selection resulted in reduced quinate and chlorogenic acid levels in cultivated lettuce. Our study thus provides beneficial genetic resources for lettuce quality improvement and sheds light on the domestication and evolution of this important leafy green.

Natural variation at FLM splicing has pleiotropic effects modulating ecological strategies in Arabidopsis thaliana

Investigating the evolution of complex phenotypes and the underlying molecular bases of their variation is critical to understand how organisms adapt to their environment. Applying classical quantitative genetics on a segregating population derived from a Can-0xCol-0 cross, we identify the MADS-box transcription factor FLOWERING LOCUS M (FLM) as a player of the phenotypic variation in plant growth and color. We show that allelic variation at FLM modulates plant growth strategy along the leaf economics spectrum, a trade-off between resource acquisition and resource conservation, observable across thousands of plant species. Functional differences at FLM rely on a single intronic substitution, disturbing transcript splicing and leading to the accumulation of non-functional FLM transcripts. Associations between this substitution and phenotypic and climatic data across Arabidopsis natural populations, show how noncoding genetic variation at a single gene might be adaptive through pleiotropic effects.

FLOWERING LOCUS M (FLM) is known as a repressor of Arabidopsis flowering. Here, the authors show that a single intronic substitution of FLM modulates leaf color and plant growth strategy along the leaf economics spectrum, as well as plays a role in plant adaptation.

The Gene scb-1 Underlies Variation in Caenorhabditis elegans Chemotherapeutic Responses

Pleiotropy, the concept that a single gene controls multiple distinct traits, is prevalent in most organisms and has broad implications for medicine and agriculture. The identification of the molecular mechanisms underlying pleiotropy has the power to reveal previously unknown biological connections between seemingly unrelated traits. Additionally, the discovery of pleiotropic genes increases our understanding of both genetic and phenotypic complexity by characterizing novel gene functions. Quantitative trait locus (QTL) mapping has been used to identify several pleiotropic regions in many organisms. However, gene knockout studies are needed to eliminate the possibility of tightly linked, non-pleiotropic loci. Here, we use a panel of 296 recombinant inbred advanced intercross lines of Caenorhabditis elegans and a high-throughput fitness assay to identify a single large-effect QTL on the center of chromosome V associated with variation in responses to eight chemotherapeutics. We validate this QTL with near-isogenic lines and pair genome-wide gene expression data with drug response traits to perform mediation analysis, leading to the identification of a pleiotropic candidate gene, scb-1, for some of the eight chemotherapeutics. Using deletion strains created by genome editing, we show that scb-1, which was previously implicated in response to bleomycin, also underlies responses to other double-strand DNA break-inducing chemotherapeutics. This finding provides new evidence for the role of scb-1 in the nematode drug response and highlights the power of mediation analysis to identify causal genes.

The Mechanosensitive Ion Channel MSL10 Potentiates Responses to Cell Swelling in Arabidopsis Seedlings

The ability to respond to unanticipated increases in volume is a fundamental property of cells, essential for cellular integrity in the face of osmotic challenges. Plants must manage cell swelling during flooding, rehydration, and pathogen invasion-but little is known about the mechanisms by which this occurs. It has been proposed that plant cells could sense and respond to cell swelling through the action of mechanosensitive ion channels. Here, we characterize a new assay to study the effects of cell swelling on Arabidopsis thaliana seedlings and to test the contributions of the mechanosensitive ion channel MscS-like10 (MSL10). The assay incorporates both cell wall softening and hypo-osmotic treatment to induce cell swelling. We show that MSL10 is required for several previously demonstrated responses to hypo-osmotic shock, including a cytoplasmic calcium transient within the first few seconds, accumulation of ROS within the first 30 min, and increased transcript levels of mechano-inducible genes within 60 min. We also show that cell swelling induces programmed cell death within 3 h in a MSL10-dependent manner. Finally, we show that MSL10 is unable to potentiate cell swelling-induced death when phosphomimetic residues are introduced into its soluble N terminus. Thus, MSL10 functions as a phospho-regulated membrane-based sensor that connects the perception of cell swelling to a downstream signaling cascade and programmed cell death.

UDP-Sugar Producing Pyrophosphorylases: Distinct and Essential Enzymes With Overlapping Substrate Specificities, Providing de novo Precursors for Glycosylation Reactions

Nucleotide sugars are the key precursors for all glycosylation reactions and are required both for oligo- and polysaccharides synthesis and protein and lipid glycosylation. Among all nucleotide sugars, UDP-sugars are the most important precursors for biomass production in nature (e.g., synthesis of cellulose, hemicellulose, and pectins for cell wall production). Several recent studies have already suggested a potential role for UDP-Glc in plant growth and development, and UDP-Glc has also been suggested as a signaling molecule, in addition to its precursor function. In this review, we will cover primary mechanisms of formation of UDP-sugars, by focusing on UDP-sugar metabolizing pyrophosphorylases. The pyrophosphorylases can be divided into three families: UDP-Glc pyrophosphorylase (UGPase), UDP-sugar pyrophosphorylase (USPase), and UDP-N-acetyl glucosamine pyrophosphorylase (UAGPase), which can be distinguished both by their amino acid sequences and by differences in substrate specificity. Substrate specificities of these enzymes are discussed, along with structure-function relationships, based on their crystal structures and homology modeling. Earlier studies with transgenic plants have revealed that each of the pyrophosphorylases is essential for plant survival, and their loss or a decrease in activity results in reproductive impairment. This constitutes a problem when studying exact in vivo roles of the enzymes using classical reverse genetics approaches. Thus, strategies involving the use of specific inhibitors (reverse chemical genetics) are also discussed. Further characterization of the properties/roles of pyrophosphorylases should address fundamental questions dealing with mechanisms and control of carbohydrate synthesis and may allow to identify targets for manipulation of biomass production in plants.

High-resolution magic angle spinning NMR studies for metabolic characterization of Arabidopsis thaliana mutants with enhanced growth characteristics

Developing smart crops which yield more biomass to meet the increasing demand for plant biomass has been an active area of research in last few decades. We investigated metabolic alterations in two Arabidopsis thaliana mutants with enhanced growth characteristics that were previously obtained from a collection of plant lines expressing artificial transcription factors. The metabolic profiles were obtained directly from intact Arabidopsis leaves using high-resolution magic angle spinning (HR-MAS) NMR. Multivariate analysis showed significant alteration of metabolite levels between the mutants and the wild-type Col-0. Interestingly, most of the metabolites that were reduced in the faster-growing mutants are generally involved in the defence against stress. These results suggest a growth-defence trade-off in the phenotypically engineered mutants. Our results further corroborate the idea that plant growth can be enhanced by suppressing defence pathways.

Modulation of ACD6 dependent hyperimmunity by natural alleles of an Arabidopsis thaliana NLR resistance gene

Plants defend themselves against pathogens by activating an array of immune responses. Unfortunately, immunity programs may also cause unintended collateral damage to the plant itself. The quantitative disease resistance gene ACCELERATED CELL DEATH 6 (ACD6) serves to balance growth and pathogen resistance in natural populations of Arabidopsis thaliana. An autoimmune allele, ACD6-Est, which strongly reduces growth under specific laboratory conditions, is found in over 10% of wild strains. There is, however, extensive variation in the strength of the autoimmune phenotype expressed by strains with an ACD6-Est allele, indicative of genetic modifiers. Quantitative genetic analysis suggests that ACD6 activity can be modulated in diverse ways, with different strains often carrying different large-effect modifiers. One modifier is SUPPRESSOR OF NPR1-1, CONSTITUTIVE 1 (SNC1), located in a highly polymorphic cluster of nucleotide-binding domain and leucine-rich repeat (NLR) immune receptor genes, which are prototypes for qualitative disease resistance genes. Allelic variation at SNC1 correlates with ACD6-Est activity in multiple accessions, and a common structural variant affecting the NL linker sequence can explain differences in SNC1 activity. Taken together, we find that an NLR gene can mask the activity of an ACD6 autoimmune allele in natural A. thaliana populations, thereby linking different arms of the plant immune system.

Plants defend themselves against pathogens by activating immune responses. Unfortunately, these can cause unintended collateral damage to the plant itself. Nevertheless, some wild plants have genetic variants that confer a low threshold for the activation of immunity. While these enable a plant to respond particularly quickly to pathogen attack, such variants might be potentially dangerous. We are investigating one such variant of the immune gene ACCELERATED CELL DEATH 6 (ACD6) in the plant Arabidopsis thaliana. We discovered that there are variants at other genetic loci that can mask the effects of an overly active ACD6 gene. One of these genes, SUPPRESSOR OF NPR1-1, CONSTITUTIVE 1 (SNC1), codes for a known immune receptor. The SNC1 variant that attenuates ACD6 activity is rather common in A. thaliana populations, suggesting that new combinations of the hyperactive ACD6 variant and this antagonistic SNC1 variant will often arise by natural crosses. Similarly, because the two genes are unlinked, outcrossing will often lead to the hyperactive ACD6 variants being unmasked again. We propose that allelic diversity at SNC1 contributes to the maintenance of the hyperactive ACD6 variant in natural A. thaliana populations.

Plastidial NAD-Dependent Malate Dehydrogenase: A Moonlighting Protein Involved in Early Chloroplast Development through Its Interaction with an FtsH12-FtsHi Protease Complex

Malate dehydrogenases (MDHs) convert malate to oxaloacetate using NAD(H) or NADP(H) as a cofactor. Arabidopsis thaliana mutants lacking plastidial NAD-dependent MDH (pdnad-mdh) are embryo-lethal, and constitutive silencing (miR-mdh-1) causes a pale, dwarfed phenotype. The reason for these severe phenotypes is unknown. Here, we rescued the embryo lethality of pdnad-mdh via embryo-specific expression of pdNAD-MDH. Rescued seedlings developed white leaves with aberrant chloroplasts and failed to reproduce. Inducible silencing of pdNAD-MDH at the rosette stage also resulted in white newly emerging leaves. These data suggest that pdNAD-MDH is important for early plastid development, which is consistent with the reductions in major plastidial galactolipid, carotenoid, and protochlorophyllide levels in miR-mdh-1 seedlings. Surprisingly, the targeting of other NAD-dependent MDH isoforms to the plastid did not complement the embryo lethality of pdnad-mdh, while expression of enzymatically inactive pdNAD-MDH did. These complemented plants grew indistinguishably from the wild type. Both active and inactive forms of pdNAD-MDH interact with a heteromeric AAA-ATPase complex at the inner membrane of the chloroplast envelope. Silencing the expression of FtsH12, a key member of this complex, resulted in a phenotype that strongly resembles miR-mdh-1. We propose that pdNAD-MDH is essential for chloroplast development due to its moonlighting role in stabilizing FtsH12, distinct from its enzymatic function.