Genomes of Subaerial Zygnematophyceae Provide Insights into Land Plant Evolution.

The transition to a terrestrial environment, termed terrestrialization, is generally regarded as a pivotal event in the evolution and diversification of the land plant flora that changed the surface of our planet. Through phylogenomic studies, a group of streptophyte algae, the Zygnematophyceae, have recently been recognized as the likely sister group to land plants (embryophytes). Here, we report genome sequences and analyses of two early diverging Zygnematophyceae (Spirogloea muscicola gen. nov. and Mesotaenium endlicherianum) that share the same subaerial/terrestrial habitat with the earliest-diverging embryophytes, the bryophytes. We provide evidence that genes (i.e., GRAS and PYR/PYL/RCAR) that increase resistance to biotic and abiotic stresses in land plants, in particular desiccation, originated or expanded in the common ancestor of Zygnematophyceae and embryophytes, and were gained by horizontal gene transfer (HGT) from soil bacteria. These two Zygnematophyceae genomes represent a cornerstone for future studies to understand the underlying molecular mechanism and process of plant terrestrialization.

Harnessing landrace diversity empowers wheat breeding.

Harnessing genetic diversity in major staple crops through the development of new breeding capabilities is essential to ensure food security1. Here we examined the genetic and phenotypic diversity of the A. E. Watkins landrace collection2 of bread wheat (Triticum aestivum), a major global cereal, by whole-genome re-sequencing of 827 Watkins landraces and 208 modern cultivars and in-depth field evaluation spanning a decade. We found that modern cultivars are derived from two of the seven ancestral groups of wheat and maintain very long-range haplotype integrity. The remaining five groups represent untapped genetic sources, providing access to landrace-specific alleles and haplotypes for breeding. Linkage disequilibrium-based haplotypes and association genetics analyses link Watkins genomes to the thousands of identified high-resolution quantitative trait loci and significant marker-trait associations. Using these structured germplasm, genotyping and informatics resources, we revealed many Watkins-unique beneficial haplotypes that can confer superior traits in modern wheat. Furthermore, we assessed the phenotypic effects of 44,338 Watkins-unique haplotypes, introgressed from 143 prioritized quantitative trait loci in the context of modern cultivars, bridging the gap between landrace diversity and current breeding. This study establishes a framework for systematically utilizing genetic diversity in crop improvement to achieve sustainable food security.

Development of a next generation SNP genotyping array for wheat.

High-throughput genotyping arrays have provided a cost-effective, reliable and interoperable system for genotyping hexaploid wheat and its relatives. Existing, highly cited arrays including our 35K Wheat Breeder's array and the Illumina 90K array were designed based on a limited amount of varietal sequence diversity and with imperfect knowledge of SNP positions. Recent progress in wheat sequencing has given us access to a vast pool of SNP diversity, whilst technological improvements have allowed us to fit significantly more probes onto a 384-well format Axiom array than previously possible. Here we describe a novel Axiom genotyping array, the 'Triticum aestivum Next Generation' array (TaNG), largely derived from whole genome skim sequencing of 204 elite wheat lines and 111 wheat landraces taken from the Watkins 'Core Collection'. We used a novel haplotype optimization approach to select SNPs with the highest combined varietal discrimination and a design iteration step to test and replace SNPs which failed to convert to reliable markers. The final design with 43 372 SNPs contains a combination of haplotype-optimized novel SNPs and legacy cross-platform markers. We show that this design has an improved distribution of SNPs compared to previous arrays and can be used to generate genetic maps with a significantly higher number of distinct bins than our previous array. We also demonstrate the improved performance of TaNGv1.1 for Genome-wide association studies (GWAS) and its utility for Copy Number Variation (CNV) analysis. The array is commercially available with supporting marker annotations and initial genotyping results freely available.

Three-nucleotide periodicity of nucleotide diversity in a population enables the identification of open reading frames.

Accurate prediction of open reading frames (ORFs) is important for studying and using genome sequences. Ribosomes move along mRNA strands with a step of three nucleotides and datasets carrying this information can be used to predict ORFs. The ribosome-protected footprints (RPFs) feature a significant 3-nt periodicity on mRNAs and are powerful in predicting translating ORFs, including small ORFs (sORFs), but the application of RPFs is limited because they are too short to be accurately mapped in complex genomes. In this study, we found a significant 3-nt periodicity in the datasets of populational genomic variants in coding sequences, in which the nucleotide diversity increases every three nucleotides. We suggest that this feature can be used to predict ORFs and develop the Python package 'OrfPP', which recovers ~83% of the annotated ORFs in the tested genomes on average, independent of the population sizes and the complexity of the genomes. The novel ORFs, including sORFs, identified from single-nucleotide polymorphisms are supported by protein mass spectrometry evidence comparable to that of the annotated ORFs. The application of OrfPP to tetraploid cotton and hexaploid wheat genomes successfully identified 76.17% and 87.43% of the annotated ORFs in the genomes, respectively, as well as 4704 sORFs, including 1182 upstream and 2110 downstream ORFs in cotton and 5025 sORFs, including 232 upstream and 234 downstream ORFs in wheat. Overall, we propose an alternative and supplementary approach for ORF prediction that can extend the studies of sORFs to more complex genomes.

The Evolution of C4 Photosynthesis in Flaveria (Asteraceae): Insights from the Flaveria linearis Complex.

Flaveria is a leading model for C4 plant evolution due to the presence of a dozen C3-C4 intermediate species, many of which are associated with a phylogenetic complex centered around Flaveria linearis. To investigate C4 evolution in Flaveria, we updated the Flaveria phylogeny and evaluated gas exchange, starch δ13C, and activity of C4 cycle enzymes in 19 Flaveria species and 28 populations within the F. linearis complex. A principal component analysis identified six functional clusters: (1) C3, (2) sub-C2, (3) full C2, (4) enriched C2, (5) sub-C4, and (6) fully C4 species. The sub-C2 species lacked a functional C4 cycle, while a gradient was present in the C2 clusters from little to modest C4 cycle activity as indicated by δ13C and enzyme activities. Three Yucatan populations of F. linearis had photosynthetic CO2 compensation points equivalent to C4 plants but showed little evidence for an enhanced C4 cycle, indicating they have an optimized C2 pathway that recaptures all photorespired CO2 in the bundle sheath (BS) tissue. All C2 species had enhanced aspartate aminotransferase activity relative to C3 species and most had enhanced alanine aminotransferase activity. These aminotransferases form aspartate and alanine from glutamate and in doing so could help return photorespiratory nitrogen (N) from BS to mesophyll cells, preventing glutamate feedback onto photorespiratory N assimilation. Their use requires upregulation of parts of the C4 metabolic cycle to generate carbon skeletons to sustain N return to the mesophyll, and thus could facilitate the evolution of the full C4 photosynthetic pathway.

Wheat genetic resources have avoided disease pandemics, improved food security, and reduced environmental footprints: A review of historical impacts and future opportunities.

The use of plant genetic resources (PGR)-wild relatives, landraces, and isolated breeding gene pools-has had substantial impacts on wheat breeding for resistance to biotic and abiotic stresses, while increasing nutritional value, end-use quality, and grain yield. In the Global South, post-Green Revolution genetic yield gains are generally achieved with minimal additional inputs. As a result, production has increased, and millions of hectares of natural ecosystems have been spared. Without PGR-derived disease resistance, fungicide use would have easily doubled, massively increasing selection pressure for fungicide resistance. It is estimated that in wheat, a billion liters of fungicide application have been avoided just since 2000. This review presents examples of successful use of PGR including the relentless battle against wheat rust epidemics/pandemics, defending against diseases that jump species barriers like blast, biofortification giving nutrient-dense varieties and the use of novel genetic variation for improving polygenic traits like climate resilience. Crop breeding genepools urgently need to be diversified to increase yields across a range of environments (>200 Mha globally), under less predictable weather and biotic stress pressure, while increasing input use efficiency. Given that the ~0.8 m PGR in wheat collections worldwide are relatively untapped and massive impacts of the tiny fraction studied, larger scale screenings and introgression promise solutions to emerging challenges, facilitated by advanced phenomic and genomic tools. The first translocations in wheat to modify rhizosphere microbiome interaction (reducing biological nitrification, reducing greenhouse gases, and increasing nitrogen use efficiency) is a landmark proof of concept. Phenomics and next-generation sequencing have already elucidated exotic haplotypes associated with biotic and complex abiotic traits now mainstreamed in breeding. Big data from decades of global yield trials can elucidate the benefits of PGR across environments. This kind of impact cannot be achieved without widescale sharing of germplasm and other breeding technologies through networks and public-private partnerships in a pre-competitive space.

Fusarium graminearum effector FgEC1 targets wheat TaGF14b protein to suppress TaRBOHD-mediated ROS production and promote infection.

Fusarium head blight (FHB), caused by Fusarium graminearum, is a devastating disease of wheat globally. However, the molecular mechanisms underlying the interactions between F. graminearum and wheat remain unclear. Here, we identified a secreted effector protein, FgEC1, that is induced during wheat infection and is required for F. graminearum virulence. FgEC1 suppressed flg22- and chitin-induced callose deposition and reactive oxygen species (ROS) burst in Nicotiana benthamiana. FgEC1 directly interacts with TaGF14b, which is upregulated in wheat heads during F. graminearum infection. Overexpression of TaGF14b increases FHB resistance in wheat without compromising yield. TaGF14b interacts with NADPH oxidase respiratory burst oxidase homolog D (TaRBOHD) and protects it against degradation by the 26S proteasome. FgEC1 inhibited the interaction of TaGF14b with TaRBOHD and promoted TaRBOHD degradation, thereby reducing TaRBOHD-mediated ROS production. Our findings reveal a novel pathogenic mechanism in which a fungal pathogen acts via an effector to reduce TaRBOHD-mediated ROS production.

The genome of Eleocharis vivipara elucidates the genetics of C3-C4 photosynthetic plasticity and karyotype evolution in the Cyperaceae.

Eleocharis vivipara, an amphibious sedge in the Cyperaceae family, has several remarkable properties, most notably its alternate use of C3 photosynthesis underwater and C4 photosynthesis on land. However, the absence of genomic data has hindered its utility for evolutionary and genetic research. Here, we present a high-quality genome for E. vivipara, representing the first chromosome-level genome for the Eleocharis genus, with an approximate size of 965.22 Mb mainly distributed across 10 chromosomes. Its Hi-C pattern, chromosome clustering results, and one-to-one genome synteny across two subgroups indicates a tetraploid structure with chromosome count 2n = 4x = 20. Phylogenetic analysis suggests that E. vivipara diverged from Cyperus esculentus approximately 32.96 million years ago (Mya), and underwent a whole-genome duplication (WGD) about 3.5 Mya. Numerous fusion and fission events were identified between the chromosomes of E. vivipara and its close relatives. We demonstrate that E. vivipara has holocentromeres, a chromosomal feature which can maintain the stability of such chromosomal rearrangements. Experimental transplantation and cross-section studies showed its terrestrial culms developed C4 Kranz anatomy with increased number of chloroplasts in the bundle sheath (BS) cells. Gene expression and weighted gene co-expression network analysis (WGCNA) showed overall elevated expression of core genes associated with the C4 pathway, and significant enrichment of genes related to modified culm anatomy and photosynthesis efficiency. We found evidence of mixed nicotinamide adenine dinucleotide - malic enzyme and phosphoenolpyruvate carboxykinase type C4 photosynthesis in E. vivipara, and hypothesize that the evolution of C4 photosynthesis predates the WGD event. The mixed type is dominated by subgenome A and supplemented by subgenome B. Collectively, our findings not only shed light on the evolution of E. vivipara and karyotype within the Cyperaceae family, but also provide valuable insights into the transition between C3 and C4 photosynthesis, offering promising avenues for crop improvement and breeding.

Duplication and sub-functionalization of flavonoid biosynthesis genes plays important role in Leguminosae root nodule symbiosis evolution.

Gene innovation plays an essential role in trait evolution. Rhizobial symbioses, the most important N2-fixing agent in agricultural systems that exists mainly in Leguminosae, is one of the most attractive evolution events. However, the gene innovations underlying Leguminosae root nodule symbiosis (RNS) remain largely unknown. Here, we investigated the gene gain event in Leguminosae RNS evolution through comprehensive phylogenomic analyses. We revealed that Leguminosae-gain genes were acquired by gene duplication and underwent a strong purifying selection. Kyoto Encyclopedia of Genes and Genomes analyses showed that the innovated genes were enriched in flavonoid biosynthesis pathways, particular downstream of chalcone synthase (CHS). Among them, Leguminosae-gain type Ⅱ chalcone isomerase (CHI) could be further divided into CHI1A and CHI1B clades, which resulted from the products of tandem duplication. Furthermore, the duplicated CHI genes exhibited exon-intron structural divergences evolved through exon/intron gain/loss and insertion/deletion. Knocking down CHI1B significantly reduced nodulation in Glycine max (soybean) and Medicago truncatula; whereas, knocking down its duplication gene CHI1A had no effect on nodulation. Therefore, Leguminosae-gain type Ⅱ CHI participated in RNS and the duplicated CHI1A and CHI1B genes exhibited RNS functional divergence. This study provides functional insights into Leguminosae-gain genetic innovation and sub-functionalization after gene duplication that contribute to the evolution and adaptation of RNS in Leguminosae.

Dynamic Diversity of NLR Genes in Triticum and Mining of Promising NLR Alleles for Disease Resistance.

Bread wheat is an essential crop with the second-highest global production after maize. Currently, wheat diseases are a serious threat to wheat production. Therefore, efficient breeding for disease resistance is extremely urgent in modern wheat. Here, we identified 2012 NLR genes from hexaploid wheat, and Ks values of paired syntenic NLRs showed a significant peak at 3.1-6.3 MYA, which exactly coincided with the first hybridization event between A and B genome lineages at ~5.5 MYA. We provided a landscape of dynamic diversity of NLRs from Triticum and Aegilops and found that NLR genes have higher diversity in wild progenitors and relatives. Further, most NLRs had opposite diversity patterns between genic and 2 Kb-promoter regions, which might respectively link sub/neofunctionalization and loss of duplicated NLR genes. Additionally, we identified an alien introgression of chromosome 4A in tetraploid emmer wheat, which was similar to that in hexaploid wheat. Transcriptome data from four experiments of wheat disease resistance helped to profile the expression pattern of NLR genes and identified promising NLRs involved in broad-spectrum disease resistance. Our study provided insights into the diversity evolution of NLR genes and identified beneficial NLRs to deploy into modern wheat in future wheat disease-resistance breeding.

Genetic insights into the evolution of genera with the eastern Asia-eastern North America floristic disjunction: a transcriptomics analysis.

PREMISE: Large disjunctions in species distributions provide excellent opportunities to study processes that shape biogeographic patterns. One such disjunction is the eastern Asia-eastern North America (EA-ENA) floristic disjunction. For many genera with this disjunction, species richness is greater in EA than in ENA; this pattern has been attributed, in part, to higher rates of molecular evolution and speciation in EA. Longer branch lengths have been found in some EA clades, relative to their ENA sister clades, suggesting that the EA lineages have evolved at a higher rate, possibly due to environmental heterogeneity, potentially contributing to the species richness anomaly. METHODS: To evaluate whether rates of molecular evolution are elevated in EA relative to ENA, we used transcriptomes from species in 11 genera displaying this disjunction. Rates of molecular evolution were estimated for up to 385 orthologous nuclear loci per genus. RESULTS: No statistically significant differences were identified in pairwise comparisons between EA and ENA sister species, suggesting equal rates of molecular evolution for both species; the data also suggest similar selection pressures in both regions. For larger genera, evidence likewise argues against more species-rich clades having higher molecular evolutionary rates, regardless of region. Our results suggest that genes across multiple gene ontology categories are evolving at similar rates under purifying selection in species in both regions. CONCLUSIONS: Our data support the hypothesis that greater species richness in EA than ENA is due to factors other than an overall increase in rates of molecular evolution in EA.

Enhancement of de novo sequencing, assembly and annotation of the Mongolian gerbil genome with transcriptome sequencing and assembly from several different tissues.

BACKGROUND: The Mongolian gerbil (Meriones unguiculatus) has historically been used as a model organism for the auditory and visual systems, stroke/ischemia, epilepsy and aging related research since 1935 when laboratory gerbils were separated from their wild counterparts. In this study we report genome sequencing, assembly, and annotation further supported by transcriptome sequencing and assembly from 27 different tissues samples. RESULTS: The genome was sequenced using Illumina HiSeq 2000 and after assembly resulted in a final genome size of 2.54 Gbp with contig and scaffold N50 values of 31.4 Kbp and 500.0 Kbp, respectively. Based on the k-mer estimated genome size of 2.48 Gbp, the assembly appears to be complete. The genome annotation was supported by transcriptome data that identified 31,769 (> 2000 bp) predicted protein-coding genes across 27 tissue samples. A BUSCO search of 3023 mammalian groups resulted in 86% of curated single copy orthologs present among predicted genes, indicating a high level of completeness of the genome. CONCLUSIONS: We report the first de novo assembly of the Mongolian gerbil genome enhanced by assembly of transcriptome data from several tissues. Sequencing of this genome and transcriptome increases the utility of the gerbil as a model organism, opening the availability of now widely used genetic tools.

African Orphan Crops Consortium (AOCC): status of developing genomic resources for African orphan crops.

MAIN CONCLUSION: The African Orphan Crops Consortium (AOCC) successfully initiated the ambitious genome sequencing project of 101 African orphan crops/trees with 6 genomes sequenced, 6 near completion, and 20 currently in progress. Addressing stunting, malnutrition, and hidden hunger through nutritious, economic, and resilient agri-food system is one of the major agricultural challenges of this century. As sub-Saharan Africa harbors a large portion of the severely malnourished population, the African Orphan Crops Consortium (AOCC) was established in 2011 with an aim to reduce stunting and malnutrition by providing nutritional security through improving locally adapted nutritious, but neglected, under-researched or orphan African food crops. Foods from these indigenous or naturalized crops and trees are rich in minerals, vitamins, and antioxidant, and are an integral part of the dietary portfolio and cultural, social, and economic milieu of African farmers. Through stakeholder consultations supported by the African Union, 101 African orphan and under-researched crop species were prioritized to mainstream into African agri-food systems. The AOCC, through a network of international-regional-public-private partnerships and collaborations, is generating genomic resources of three types, i.e., reference genome sequence, transcriptome sequence, and re-sequencing 100 accessions/species, using next-generation sequencing (NGS) technology. Furthermore, the University of California Davis African Plant Breeding Academy under the AOCC banner is training 150 lead African scientists to breed high yielding, nutritious, and climate-resilient (biotic and abiotic stress tolerant) crop varieties that meet African farmer and consumer needs. To date, one or more forms of sequence data have been produced for 60 crops. Reference genome sequences for six species have already been published, 6 are almost near completion, and 19 are in progress.

Is there foul play in the leaf pocket? The metagenome of floating fern Azolla reveals endophytes that do not fix N2 but may denitrify.

Dinitrogen fixation by Nostoc azollae residing in specialized leaf pockets supports prolific growth of the floating fern Azolla filiculoides. To evaluate contributions by further microorganisms, the A. filiculoides microbiome and nitrogen metabolism in bacteria persistently associated with Azolla ferns were characterized. A metagenomic approach was taken complemented by detection of N2 O released and nitrogen isotope determinations of fern biomass. Ribosomal RNA genes in sequenced DNA of natural ferns, their enriched leaf pockets and water filtrate from the surrounding ditch established that bacteria of A. filiculoides differed entirely from surrounding water and revealed species of the order Rhizobiales. Analyses of seven cultivated Azolla species confirmed persistent association with Rhizobiales. Two distinct nearly full-length Rhizobiales genomes were identified in leaf-pocket-enriched samples from ditch grown A. filiculoides. Their annotation revealed genes for denitrification but not N2 -fixation. 15 N2 incorporation was active in ferns with N. azollae but not in ferns without. N2 O was not detectably released from surface-sterilized ferns with the Rhizobiales. N2 -fixing N. azollae, we conclude, dominated the microbiome of Azolla ferns. The persistent but less abundant heterotrophic Rhizobiales bacteria possibly contributed to lowering O2 levels in leaf pockets but did not release detectable amounts of the strong greenhouse gas N2 O.

Redefining the structural motifs that determine RNA binding and RNA editing by pentatricopeptide repeat proteins in land plants.

The pentatricopeptide repeat (PPR) proteins form one of the largest protein families in land plants. They are characterised by tandem 30-40 amino acid motifs that form an extended binding surface capable of sequence-specific recognition of RNA strands. Almost all of them are post-translationally targeted to plastids and mitochondria, where they play important roles in post-transcriptional processes including splicing, RNA editing and the initiation of translation. A code describing how PPR proteins recognise their RNA targets promises to accelerate research on these proteins, but making use of this code requires accurate definition and annotation of all of the various nucleotide-binding motifs in each protein. We have used a structural modelling approach to define 10 different variants of the PPR motif found in plant proteins, in addition to the putative deaminase motif that is found at the C-terminus of many RNA-editing factors. We show that the super-helical RNA-binding surface of RNA-editing factors is potentially longer than previously recognised. We used the redefined motifs to develop accurate and consistent annotations of PPR sequences from 109 genomes. We report a high error rate in PPR gene models in many public plant proteomes, due to gene fusions and insertions of spurious introns. These consistently annotated datasets across a wide range of species are valuable resources for future comparative genomics studies, and an essential pre-requisite for accurate large-scale computational predictions of PPR targets. We have created a web portal (http://www.plantppr.com) that provides open access to these resources for the community.

Genome sequence and analysis of the tuber crop potato.

Potato (Solanum tuberosum L.) is the world's most important non-grain food crop and is central to global food security. It is clonally propagated, highly heterozygous, autotetraploid, and suffers acute inbreeding depression. Here we use a homozygous doubled-monoploid potato clone to sequence and assemble 86% of the 844-megabase genome. We predict 39,031 protein-coding genes and present evidence for at least two genome duplication events indicative of a palaeopolyploid origin. As the first genome sequence of an asterid, the potato genome reveals 2,642 genes specific to this large angiosperm clade. We also sequenced a heterozygous diploid clone and show that gene presence/absence variants and other potentially deleterious mutations occur frequently and are a likely cause of inbreeding depression. Gene family expansion, tissue-specific expression and recruitment of genes to new pathways contributed to the evolution of tuber development. The potato genome sequence provides a platform for genetic improvement of this vital crop.

Integrated genome sequence and linkage map of physic nut (Jatropha curcas L.), a biodiesel plant.

The family Euphorbiaceae includes some of the most efficient biomass accumulators. Whole genome sequencing and the development of genetic maps of these species are important components in molecular breeding and genetic improvement. Here we report the draft genome of physic nut (Jatropha curcas L.), a biodiesel plant. The assembled genome has a total length of 320.5 Mbp and contains 27,172 putative protein-coding genes. We established a linkage map containing 1208 markers and anchored the genome assembly (81.7%) to this map to produce 11 pseudochromosomes. After gene family clustering, 15,268 families were identified, of which 13,887 existed in the castor bean genome. Analysis of the genome highlighted specific expansion and contraction of a number of gene families during the evolution of this species, including the ribosome-inactivating proteins and oil biosynthesis pathway enzymes. The genomic sequence and linkage map provide a valuable resource not only for fundamental and applied research on physic nut but also for evolutionary and comparative genomics analysis, particularly in the Euphorbiaceae.

Transcriptomic, proteomic and metabolic changes in Arabidopsis thaliana leaves after the onset of illumination.

BACKGROUND: Light plays an important role in plant growth and development. In this study, the impact of light on physiology of 20-d-old Arabidopsis leaves was examined through transcriptomic, proteomic and metabolomic analysis. Since the energy-generating electron transport chains in chloroplasts and mitochondria are encoded by both nuclear and organellar genomes, sequencing total RNA after removal of ribosomal RNAs provides essential information on transcription of organellar genomes. The changes in the levels of ADP, ATP, NADP(+), NADPH and 41 metabolites upon illumination were also quantified. RESULTS: Upon illumination, while the transcription of the genes encoded by the plastid genome did not change significantly, the transcription of nuclear genes encoding different functional complexes in the photosystem are differentially regulated whereas members of the same complex are co-regulated with each other. The abundance of mRNAs and proteins encoded by all three genomes are, however, not always positively correlated. One such example is the negative correlation between mRNA and protein abundances of the photosystem components, which reflects the importance of post-transcriptional regulation in plant physiology. CONCLUSION: This study provides systems-wide datasets which allow plant researchers to examine the changes in leaf transcriptomes, proteomes and key metabolites upon illumination and to determine whether there are any correlations between changes in transcript and protein abundances of a particular gene or pathway upon illumination. The integration of data of the organelles and the photosystems, Calvin-Benson cycle, carbohydrate metabolism, glycolysis, the tricarboxylic acid cycle and respiratory chain, thereby provides a more complete picture to the changes in plant physiology upon illumination than has been attained to date.

Phosphorylation and Dephosphorylation of the Presequence of Precursor MULTIPLE ORGANELLAR RNA EDITING FACTOR3 during Import into Mitochondria from Arabidopsis.

The nucleus-encoded mitochondria-targeted proteins, multiple organellar RNA editing factors (MORF3, MORF5, and MORF6), interact with Arabidopsis (Arabidopsis thaliana) PURPLE ACID PHOSPHATASE2 (AtPAP2) located on the chloroplast and mitochondrial outer membranes in a presequence-dependent manner. Phosphorylation of the presequence of the precursor MORF3 (pMORF3) by endogenous kinases in wheat germ translation lysate, leaf extracts, or STY kinases, but not in rabbit reticulocyte translation lysate, resulted in the inhibition of protein import into mitochondria. This inhibition of import could be overcome by altering threonine/serine residues to alanine on the presequence, thus preventing phosphorylation. Phosphorylated pMORF3, but not the phosphorylation-deficient pMORF3, can form a complex with 14-3-3 proteins and HEAT SHOCK PROTEIN70. The phosphorylation-deficient mutant of pMORF3 also displayed faster rates of import when translated in wheat germ lysates. Mitochondria isolated from plants with altered amounts of AtPAP2 displayed altered protein import kinetics. The import rate of pMORF3 synthesized in wheat germ translation lysate into pap2 mitochondria was slower than that into wild-type mitochondria, and this rate disparity was not seen for pMORF3 synthesized in rabbit reticulocyte translation lysate, the latter translation lysate largely deficient in kinase activity. Taken together, these results support a role for the phosphorylation and dephosphorylation of pMORF3 during the import into plant mitochondria. These results suggest that kinases, possibly STY kinases, and AtPAP2 are involved in the import of protein into both mitochondria and chloroplasts and provide a mechanism by which the import of proteins into both organelles may be coordinated.