Allopolyploid origin and diversification of the Hawaiian endemic mints.

Island systems provide important contexts for studying processes underlying lineage migration, species diversification, and organismal extinction. The Hawaiian endemic mints (Lamiaceae family) are the second largest plant radiation on the isolated Hawaiian Islands. We generated a chromosome-scale reference genome for one Hawaiian species, Stenogyne calaminthoides, and resequenced 45 relatives, representing 34 species, to uncover the continental origins of this group and their subsequent diversification. We further resequenced 109 individuals of two Stenogyne species, and their purported hybrids, found high on the Mauna Kea volcano on the island of Hawai'i. The three distinct Hawaiian genera, Haplostachys, Phyllostegia, and Stenogyne, are nested inside a fourth genus, Stachys. We uncovered four independent polyploidy events within Stachys, including one allopolyploidy event underlying the Hawaiian mints and their direct western North American ancestors. While the Hawaiian taxa may have principally diversified by parapatry and drift in small and fragmented populations, localized admixture may have played an important role early in lineage diversification. Our genomic analyses provide a view into how organisms may have radiated on isolated island chains, settings that provided one of the principal natural laboratories for Darwin's thinking about the evolutionary process.

The genome and population genomics of allopolyploid Coffea arabica reveal the diversification history of modern coffee cultivars.

Coffea arabica, an allotetraploid hybrid of Coffea eugenioides and Coffea canephora, is the source of approximately 60% of coffee products worldwide, and its cultivated accessions have undergone several population bottlenecks. We present chromosome-level assemblies of a di-haploid C. arabica accession and modern representatives of its diploid progenitors, C. eugenioides and C. canephora. The three species exhibit largely conserved genome structures between diploid parents and descendant subgenomes, with no obvious global subgenome dominance. We find evidence for a founding polyploidy event 350,000-610,000 years ago, followed by several pre-domestication bottlenecks, resulting in narrow genetic variation. A split between wild accessions and cultivar progenitors occurred ~30.5 thousand years ago, followed by a period of migration between the two populations. Analysis of modern varieties, including lines historically introgressed with C. canephora, highlights their breeding histories and loci that may contribute to pathogen resistance, laying the groundwork for future genomics-based breeding of C. arabica.

Complex Polyploids: Origins, Genomic Composition, and Role of Introgressed Alleles.

Introgression allows polyploid species to acquire new genomic content from diploid progenitors or from other unrelated diploid or polyploid lineages, contributing to genetic diversity and facilitating adaptive allele discovery. In some cases, high levels of introgression elicit the replacement of large numbers of alleles inherited from the polyploid's ancestral species, profoundly reshaping the polyploid's genomic composition. In such complex polyploids it is often difficult to determine which taxa were the progenitor species and which taxa provided additional introgressive blocks through subsequent hybridization. Here, we use population-level genomic data to reconstruct the phylogenetic history of Betula pubescens (downy birch), a tetraploid species often assumed to be of allopolyploid origin and which is known to hybridize with at least four other birch species. This was achieved by modeling of polyploidization and introgression events under the multispecies coalescent and then using an approximate Bayesian computation (ABC) rejection algorithm to evaluate and compare competing polyploidization models. We provide evidence that B. pubescens is the outcome of an autoploid genome doubling event in the common ancestor of B. pendula and its extant sister species, B. platyphylla, that took place approximately 178,000-188,000 generations ago. Extensive hybridization with B. pendula, B. nana, and B. humilis followed in the aftermath of autopolyploidization, with the relative contribution of each of these species to the B. pubescens genome varying markedly across the species' range. Functional analysis of B. pubescens loci containing alleles introgressed from B. nana identified multiple genes involved in climate adaptation, while loci containing alleles derived from B. humilis revealed several genes involved in the regulation of meiotic stability and pollen viability in plant species.

Tree architecture: A strigolactone-deficient mutant reveals a connection between branching order and auxin gradient along the tree stem.

Due to their long lifespan, trees and bushes develop higher order of branches in a perennial manner. In contrast to a tall tree, with a clearly defined main stem and branching order, a bush is shorter and has a less apparent main stem and branching pattern. To address the developmental basis of these two forms, we studied several naturally occurring architectural variants in silver birch (Betula pendula). Using a candidate gene approach, we identified a bushy kanttarelli variant with a loss-of-function mutation in the BpMAX1 gene required for strigolactone (SL) biosynthesis. While kanttarelli is shorter than the wild type (WT), it has the same number of primary branches, whereas the number of secondary branches is increased, contributing to its bush-like phenotype. To confirm that the identified mutation was responsible for the phenotype, we phenocopied kanttarelli in transgenic BpMAX1::RNAi birch lines. SL profiling confirmed that both kanttarelli and the transgenic lines produced very limited amounts of SL. Interestingly, the auxin (IAA) distribution along the main stem differed between WT and BpMAX1::RNAi. In the WT, the auxin concentration formed a gradient, being higher in the uppermost internodes and decreasing toward the basal part of the stem, whereas in the transgenic line, this gradient was not observed. Through modeling, we showed that the different IAA distribution patterns may result from the difference in the number of higher-order branches and plant height. Future studies will determine whether the IAA gradient itself regulates aspects of plant architecture.

Phylogenetic Analysis of Allotetraploid Species Using Polarized Genomic Sequences.

Phylogenetic analysis of polyploid hybrid species has long posed a formidable challenge as it requires the ability to distinguish between alleles of different ancestral origins in order to disentangle their individual evolutionary history. This problem has been previously addressed by conceiving phylogenies as reticulate networks, using a two-step phasing strategy that first identifies and segregates homoeologous loci and then, during a second phasing step, assigns each gene copy to one of the subgenomes of an allopolyploid species. Here, we propose an alternative approach, one that preserves the core idea behind phasing-to produce separate nucleotide sequences that capture the reticulate evolutionary history of a polyploid-while vastly simplifying its implementation by reducing a complex multistage procedure to a single phasing step. While most current methods used for phylogenetic reconstruction of polyploid species require sequencing reads to be pre-phased using experimental or computational methods-usually an expensive, complex, and/or time-consuming endeavor-phasing executed using our algorithm is performed directly on the multiple-sequence alignment (MSA), a key change that allows for the simultaneous segregation and sorting of gene copies. We introduce the concept of genomic polarization that, when applied to an allopolyploid species, produces nucleotide sequences that capture the fraction of a polyploid genome that deviates from that of a reference sequence, usually one of the other species present in the MSA. We show that if the reference sequence is one of the parental species, the polarized polyploid sequence has a close resemblance (high pairwise sequence identity) to the second parental species. This knowledge is harnessed to build a new heuristic algorithm where, by replacing the allopolyploid genomic sequence in the MSA by its polarized version, it is possible to identify the phylogenetic position of the polyploid's ancestral parents in an iterative process. The proposed methodology can be used with long-read and short-read high-throughput sequencing data and requires only one representative individual for each species to be included in the phylogenetic analysis. In its current form, it can be used in the analysis of phylogenies containing tetraploid and diploid species. We test the newly developed method extensively using simulated data in order to evaluate its accuracy. We show empirically that the use of polarized genomic sequences allows for the correct identification of both parental species of an allotetraploid with up to 97% certainty in phylogenies with moderate levels of incomplete lineage sorting (ILS) and 87% in phylogenies containing high levels of ILS. We then apply the polarization protocol to reconstruct the reticulate histories of Arabidopsis kamchatica and Arabidopsis suecica, two allopolyploids whose ancestry has been well documented. [Allopolyploidy; Arabidopsis; genomic polarization; homoeologs; incomplete lineage sorting; phasing; polyploid phylogenetics; reticulate evolution.].

Genetic Diversity Analysis of Sapindus in China and Extraction of a Core Germplasm Collection Using EST-SSR Markers.

Sapindus is an important forest tree genus with utilization in biodiesel, biomedicine, and it harbors great potential for biochemical engineering applications. For advanced breeding of Sapindus, it is necessary to evaluate the genetic diversity and construct a rationally designed core germplasm collection. In this study, the genetic diversity and population structure of Sapindus were conducted with 18 expressed sequence tag-simple sequence repeat (EST-SSR) markers in order to establish a core germplasm collection from 161 Sapindus accessions. The population of Sapindus showed high genetic diversity and significant population structure. Interspecific genetic variation was significantly higher than intraspecific variation in the Sapindus mukorossi, Sapindus delavayi, and combined Sapindus rarak plus Sapindus rarak var. velutinus populations. S. mukorossi had abundant genetic variation and showed a specific pattern of geographical variation, whereas S. delavayi, S. rarak, and S. rarak var. velutinus showed less intraspecific variation. A core germplasm collection was created that contained 40% of genetic variation in the initial population, comprising 53 S. mukorossi and nine S. delavayi lineages, as well as single representatives of S. rarak and S. rarak var. velutinus. These results provide a germplasm basis and theoretical rationale for the efficient management, conservation, and utilization of Sapindus, as well as genetic resources for joint genomics research in the future.

Insights into bear evolution from a Pleistocene polar bear genome.

The polar bear (Ursus maritimus) has become a symbol of the threat to biodiversity from climate change. Understanding polar bear evolutionary history may provide insights into apex carnivore responses and prospects during periods of extreme environmental perturbations. In recent years, genomic studies have examined bear speciation and population history, including evidence for ancient admixture between polar bears and brown bears (Ursus arctos). Here, we extend our earlier studies of a 130,000- to 115,000-y-old polar bear from the Svalbard Archipelago using a 10× coverage genome sequence and 10 new genomes of polar and brown bears from contemporary zones of overlap in northern Alaska. We demonstrate a dramatic decline in effective population size for this ancient polar bear's lineage, followed by a modest increase just before its demise. A slightly higher genetic diversity in the ancient polar bear suggests a severe genetic erosion over a prolonged bottleneck in modern polar bears. Statistical fitting of data to alternative admixture graph scenarios favors at least one ancient introgression event from brown bears into the ancestor of polar bears, possibly dating back over 150,000 y. Gene flow was likely bidirectional, but allelic transfer from brown into polar bear is the strongest detected signal, which contrasts with other published work. These findings may have implications for our understanding of climate change impacts: Polar bears, a specialist Arctic lineage, may not only have undergone severe genetic bottlenecks but also been the recipient of generalist, boreal genetic variants from brown bears during critical phases of Northern Hemisphere glacial oscillations.

Chromosome-level genome assembly of the diploid blueberry Vaccinium darrowii provides insights into its subtropical adaptation and cuticle synthesis.

Vaccinium darrowii is a subtropical wild blueberry species that has been used to breed economically important southern highbush cultivars. The adaptive traits of V. darrowii to subtropical climates can provide valuable information for breeding blueberry and perhaps other plants, especially against the background of global warming. Here, we assembled the V. darrowii genome into 12 pseudochromosomes using Oxford Nanopore long reads complemented with Hi-C scaffolding technologies, and we predicted 41 815 genes using RNA-sequencing evidence. Syntenic analysis across three Vaccinium species revealed a highly conserved genome structure, with the highest collinearity between V. darrowii and Vaccinium corymbosum. This conserved genome structure may explain the high fertility observed during crossbreeding of V. darrowii with other blueberry cultivars. Analysis of gene expansion and tandem duplication indicated possible roles for defense- and flowering-associated genes in the adaptation of V. darrowii to the subtropics. Putative SOC1 genes in V. darrowii were identified based on phylogeny and expression analysis. Blueberries are covered in a thick cuticle layer and contain anthocyanins, which confer their powdery blue color. Using RNA sequencing, we delineated the cuticle biosynthesis pathways of Vaccinium species in V. darrowii. This result can serve as a reference for breeding berries whose colors are appealing to customers. The V. darrowii reference genome, together with the unique traits of this species, including its diploid genome, short vegetative phase, and high compatibility in hybridization with other blueberries, make V. darrowii a potential research model for blueberry species.

Two divergent haplotypes from a highly heterozygous lychee genome suggest independent domestication events for early and late-maturing cultivars.

Lychee is an exotic tropical fruit with a distinct flavor. The genome of cultivar 'Feizixiao' was assembled into 15 pseudochromosomes, totaling ~470 Mb. High heterozygosity (2.27%) resulted in two complete haplotypic assemblies. A total of 13,517 allelic genes (42.4%) were differentially expressed in diverse tissues. Analyses of 72 resequenced lychee accessions revealed two independent domestication events. The extremely early maturing cultivars preferentially aligned to one haplotype were domesticated from a wild population in Yunnan, whereas the late-maturing cultivars that mapped mostly to the second haplotype were domesticated independently from a wild population in Hainan. Early maturing cultivars were probably developed in Guangdong via hybridization between extremely early maturing cultivar and late-maturing cultivar individuals. Variable deletions of a 3.7 kb region encompassed by a pair of CONSTANS-like genes probably regulate fruit maturation differences among lychee cultivars. These genomic resources provide insights into the natural history of lychee domestication and will accelerate the improvement of lychee and related crops.

Genotype-Specific Expression and NLR Repertoire Contribute to Phenotypic Resistance Diversity in Plantago lanceolata.

High levels of phenotypic variation in resistance appears to be nearly ubiquitous across natural host populations. Molecular processes contributing to this variation in nature are still poorly known, although theory predicts resistance to evolve at specific loci driven by pathogen-imposed selection. Nucleotide-binding leucine-rich repeat (NLR) genes play an important role in pathogen recognition, downstream defense responses and defense signaling. Identifying the natural variation in NLRs has the potential to increase our understanding of how NLR diversity is generated and maintained, and how to manage disease resistance. Here, we sequenced the transcriptomes of five different Plantago lanceolata genotypes when inoculated by the same strain of obligate fungal pathogen Podosphaera plantaginis. A de novo transcriptome assembly of RNA-sequencing data yielded 24,332 gene models with N50 value of 1,329 base pairs and gene space completeness of 66.5%. The gene expression data showed highly varying responses where each plant genotype demonstrated a unique expression profile in response to the pathogen, regardless of the resistance phenotype. Analysis on the conserved NB-ARC domain demonstrated a diverse NLR repertoire in P. lanceolata consistent with the high phenotypic resistance diversity in this species. We find evidence of selection generating diversity at some of the NLR loci. Jointly, our results demonstrate that phenotypic resistance diversity results from a crosstalk between different defense mechanisms. In conclusion, characterizing the architecture of resistance in natural host populations may shed unprecedented light on the potential of evolution to generate variation.

Colonic Mucosal Microbiota and Association of Bacterial Taxa with the Expression of Host Antimicrobial Peptides in Pediatric Ulcerative Colitis.

Inflammatory bowel diseases (IBD), ulcerative colitis (UC) and Crohn's disease (CD), are chronic debilitating disorders of unknown etiology. Over 200 genetic risk loci are associated with IBD, highlighting a key role for immunological and epithelial barrier functions. Environmental factors account for the growing incidence of IBD, and microbiota are considered as an important contributor. Microbiota dysbiosis can lead to a loss of tolerogenic immune effects and initiate or exacerbate inflammation. We aimed to study colonic mucosal microbiota and the expression of selected host genes in pediatric UC. We used high-throughput 16S rDNA sequencing to profile microbiota in colonic biopsies of pediatric UC patients (n = 26) and non-IBD controls (n = 27). The expression of 13 genes, including five for antimicrobial peptides, in parallel biopsies was assessed with qRT-PCR. The composition of microbiota between UC and non-IBD differed significantly (PCoA, p = 0.001). UC children had a decrease in Bacteroidetes and an increase in several family-level taxa including Peptostreptococcaceae and Enterobacteriaceae, which correlated negatively with the expression of antimicrobial peptides REG3G and DEFB1, respectively. Enterobacteriaceae correlated positively with the expression siderophore binding protein LCN2 and Betaproteobacteria negatively with DEFB4A expression. The results indicate that reciprocal interaction of epithelial microbiota and defense mechanisms play a role in UC.

The photoreceptor UVR8 mediates the perception of both UV-B and UV-A wavelengths up to 350 nm of sunlight with responsivity moderated by cryptochromes.

The photoreceptors UV RESISTANCE LOCUS 8 (UVR8) and CRYPTOCHROMES 1 and 2 (CRYs) play major roles in the perception of UV-B (280-315 nm) and UV-A/blue radiation (315-500 nm), respectively. However, it is poorly understood how they function in sunlight. The roles of UVR8 and CRYs were assessed in a factorial experiment with Arabidopsis thaliana wild-type and photoreceptor mutants exposed to sunlight for 6 or 12 hr under five types of filters with cut-offs in UV and blue-light regions. Transcriptome-wide responses triggered by UV-B and UV-A wavelengths shorter than 350 nm (UV-Asw ) required UVR8 whereas those induced by blue and UV-A wavelengths longer than 350 nm (UV-Alw ) required CRYs. UVR8 modulated gene expression in response to blue light while lack of CRYs drastically enhanced gene expression in response to UV-B and UV-Asw . These results agree with our estimates of photons absorbed by these photoreceptors in sunlight and with in vitro monomerization of UVR8 by wavelengths up to 335 nm. Motif enrichment analysis predicted complex signaling downstream of UVR8 and CRYs. Our results highlight that it is important to use UV waveband definitions specific to plants' photomorphogenesis as is routinely done in the visible region.

ELIMÄKI Locus Is Required for Vertical Proprioceptive Response in Birch Trees.

Tree architecture has evolved to support a top-heavy above-ground biomass, but this integral feature poses a weight-induced challenge to trunk stability. Maintaining an upright stem is expected to require vertical proprioception through feedback between sensing stem weight and responding with radial growth. Despite its apparent importance, the principle by which plant stems respond to vertical loading forces remains largely unknown. Here, by manipulating the stem weight of downy birch (Betula pubescens) trees, we show that cambial development is modulated systemically along the stem. We carried out a genetic study on the underlying regulation by combining an accelerated birch flowering program with a recessive mutation at the ELIMÄKI locus (EKI), which causes a mechanically defective response to weight stimulus resulting in stem collapse after just 3 months. We observed delayed wood morphogenesis in eki compared with WT, along with a more mechanically elastic cambial zone and radial compression of xylem cell size, indicating that rapid tissue differentiation is critical for cambial growth under mechanical stress. Furthermore, the touch-induced mechanosensory pathway was transcriptionally misregulated in eki, indicating that the ELIMÄKI locus is required to integrate the weight-growth feedback regulation. By studying this birch mutant, we were able to dissect vertical proprioception from the gravitropic response associated with reaction wood formation. Our study provides evidence for both local and systemic responses to mechanical stimuli during secondary plant development.

PROTEIN PHOSPHATASE 2A-B'γ Controls Botrytis cinerea Resistance and Developmental Leaf Senescence.

Plants optimize their growth and survival through highly integrated regulatory networks that coordinate defensive measures and developmental transitions in response to environmental cues. Protein phosphatase 2A (PP2A) is a key signaling component that controls stress reactions and growth at different stages of plant development, and the PP2A regulatory subunit PP2A-B'γ is required for negative regulation of pathogenesis responses and for maintenance of cell homeostasis in short-day conditions. Here, we report molecular mechanisms by which PP2A-B'γ regulates Botrytis cinerea resistance and leaf senescence in Arabidopsis (Arabidopsis thaliana). We extend the molecular functionality of PP2A-B'γ to a protein kinase-phosphatase interaction with the defense-associated calcium-dependent protein kinase CPK1 and present indications this interaction may function to control CPK1 activity. In presenescent leaf tissues, PP2A-B'γ is also required to negatively control the expression of salicylic acid-related defense genes, which have recently proven vital in plant resistance to necrotrophic fungal pathogens. In addition, we find the premature leaf yellowing of pp2a-b'γ depends on salicylic acid biosynthesis via SALICYLIC ACID INDUCTION DEFICIENT2 and bears the hallmarks of developmental leaf senescence. We propose PP2A-B'γ age-dependently controls salicylic acid-related signaling in plant immunity and developmental leaf senescence.

Differences in growth and gas exchange between southern and northern provenances of silver birch (Betula pendula Roth) in northern Europe.

Due to its ubiquity across northern latitudes, silver birch (Betula pendula Roth) is an attractive model species for studying geographical trait variation and acclimation capacity. Six birch provenances from 60 to 67°N across Finland were grown in a common garden and studied for provenance and genotype variation. We looked for differences in height growth, photosynthetic gas exchange and chlorophyll content index (CCI) and compared the gas exchange of early and late leaves on short and long shoots, respectively. The provenances stratified into southern and northern groups. Northern provenances attained less height growth increment and had higher stomatal conductance (gs) and lower intrinsic water-use efficiency (WUE, Anet/gs) than southern provenances, whereas net photosynthesis (Anet) or CCI did not show clear grouping. Short shoot leaves had lower gs and higher WUE than long shoot leaves in all provenances, but there was no difference in Anet between shoot types. The separation of the provenances into two groups according to their physiological responses might reflect the evolutionary history of B. pendula. Latitudinal differences in gas exchange and water use traits can have plausible consequences for global carbon and water fluxes in a warming climate.

The avocado genome informs deep angiosperm phylogeny, highlights introgressive hybridization, and reveals pathogen-influenced gene space adaptation.

The avocado, Persea americana, is a fruit crop of immense importance to Mexican agriculture with an increasing demand worldwide. Avocado lies in the anciently diverged magnoliid clade of angiosperms, which has a controversial phylogenetic position relative to eudicots and monocots. We sequenced the nuclear genomes of the Mexican avocado race, P. americana var. drymifolia, and the most commercially popular hybrid cultivar, Hass, and anchored the latter to chromosomes using a genetic map. Resequencing of Guatemalan and West Indian varieties revealed that ∼39% of the Hass genome represents Guatemalan source regions introgressed into a Mexican race background. Some introgressed blocks are extremely large, consistent with the recent origin of the cultivar. The avocado lineage experienced 2 lineage-specific polyploidy events during its evolutionary history. Although gene-tree/species-tree phylogenomic results are inconclusive, syntenic ortholog distances to other species place avocado as sister to the enormous monocot and eudicot lineages combined. Duplicate genes descending from polyploidy augmented the transcription factor diversity of avocado, while tandem duplicates enhanced the secondary metabolism of the species. Phenylpropanoid biosynthesis, known to be elicited by Colletotrichum (anthracnose) pathogen infection in avocado, is one enriched function among tandems. Furthermore, transcriptome data show that tandem duplicates are significantly up- and down-regulated in response to anthracnose infection, whereas polyploid duplicates are not, supporting the general view that collections of tandem duplicates contribute evolutionarily recent "tuning knobs" in the genome adaptive landscapes of given species.

Mechanistic insights into the evolution of DUF26-containing proteins in land plants.

Large protein families are a prominent feature of plant genomes and their size variation is a key element for adaptation. However, gene and genome duplications pose difficulties for functional characterization and translational research. Here we infer the evolutionary history of the DOMAIN OF UNKNOWN FUNCTION (DUF) 26-containing proteins. The DUF26 emerged in secreted proteins. Domain duplications and rearrangements led to the appearance of CYSTEINE-RICH RECEPTOR-LIKE PROTEIN KINASES (CRKs) and PLASMODESMATA-LOCALIZED PROTEINS (PDLPs). The DUF26 is land plant-specific but structural analyses of PDLP ectodomains revealed strong similarity to fungal lectins and thus may constitute a group of plant carbohydrate-binding proteins. CRKs expanded through tandem duplications and preferential retention of duplicates following whole genome duplications, whereas PDLPs evolved according to the dosage balance hypothesis. We propose that new gene families mainly expand through small-scale duplications, while fractionation and genetic drift after whole genome multiplications drive families towards dosage balance.

Tissue-specific study across the stem reveals the chemistry and transcriptome dynamics of birch bark.

Tree bark is a highly specialized array of tissues that plays important roles in plant protection and development. Bark tissues develop from two lateral meristems; the phellogen (cork cambium) produces the outermost stem-environment barrier called the periderm, while the vascular cambium contributes with phloem tissues. Although bark is diverse in terms of tissues, functions and species, it remains understudied at higher resolution. We dissected the stem of silver birch (Betula pendula) into eight major tissue types, and characterized these by a combined transcriptomics and metabolomics approach. We further analyzed the varying bark types within the Betulaceae family. The two meristems had a distinct contribution to the stem transcriptomic landscape. Furthermore, inter- and intraspecies analyses illustrated the unique molecular profile of the phellem. We identified multiple tissue-specific metabolic pathways, such as the mevalonate/betulin biosynthesis pathway, that displayed differential evolution within the Betulaceae. A detailed analysis of suberin and betulin biosynthesis pathways identified a set of underlying regulators and highlighted the important role of local, small-scale gene duplication events in the evolution of metabolic pathways. This work reveals the transcriptome and metabolic diversity among bark tissues and provides insights to its development and evolution, as well as its biotechnological applications.

Arabidopsis RCD1 coordinates chloroplast and mitochondrial functions through interaction with ANAC transcription factors.

Reactive oxygen species (ROS)-dependent signaling pathways from chloroplasts and mitochondria merge at the nuclear protein RADICAL-INDUCED CELL DEATH1 (RCD1). RCD1 interacts in vivo and suppresses the activity of the transcription factors ANAC013 and ANAC017, which mediate a ROS-related retrograde signal originating from mitochondrial complex III. Inactivation of RCD1 leads to increased expression of mitochondrial dysfunction stimulon (MDS) genes regulated by ANAC013 and ANAC017. Accumulating MDS gene products, including alternative oxidases (AOXs), affect redox status of the chloroplasts, leading to changes in chloroplast ROS processing and increased protection of photosynthetic apparatus. ROS alter the abundance, thiol redox state and oligomerization of the RCD1 protein in vivo, providing feedback control on its function. RCD1-dependent regulation is linked to chloroplast signaling by 3'-phosphoadenosine 5'-phosphate (PAP). Thus, RCD1 integrates organellar signaling from chloroplasts and mitochondria to establish transcriptional control over the metabolic processes in both organelles.