Witnessing Genome Evolution: Experimental Reconstruction of Endosymbiotic and Horizontal Gene Transfer

Taken to extremes: Loss of plastid rpl32 in Streptophyta and Cuscuta's unconventional solution for its replacement

The evolution of plant genomes is riddled with exchanges of genetic material within one plant (endosymbiotic gene transfer/EGT) and between unrelated plants (horizontal gene transfer/HGT). These exchanges have left their marks on plant genomes. Parasitic plants with their special evolutionary niche provide ample examples for these processes because they are under a reduced evolutionary pressure to maintain autotrophy and thus to conserve their plastid genomes. On the other hand, the close physical connections with different hosts enabled them to acquire genetic material from other plants. Based on an analysis of an extensive dataset including the parasite Cuscuta campestris and other parasitic plant species, we identified a unique evolutionary history of rpl32 genes coding for an essential plastid ribosomal subunit in Cuscuta. Our analysis suggests that the gene was most likely sequestered by HGT from a member of the Oxalidales order serving as host to an ancestor of the Cuscuta subgenus Grammica. Oxalidales had suffered an ancestral EGT of rpl32 predating the evolution of the genus Cuscuta. The HGT subsequently relieved the plastid rpl32 from its evolutionary constraint and led to its loss from the plastid genome. The HGT-based acquisition in Cuscuta is supported by a high sequence similarity of the mature L32 protein between species of the subgenus Grammica and representatives of the Oxalidales, and by a surprisingly conserved transit peptide, whose functionality in Cuscuta was experimentally verified. The findings are discussed in view of an overall pattern of EGT events for plastid ribosomal subunits in Streptophyta.

Assembly and analysis of the complete mitochondrial genome of Carya illinoinensis to provide insights into the conserved sequences of tRNA genes

Carya illinoinensis is an economically important nut tree, and its chloroplast (cp.) genome has been reported; however, its mitochondrial (mt) genome remains unknown. In the present study, we assembled the first mt genome of C. illinoinensis. The circular mt genome of C. illinoinensis is 495,205 bp long, with 37 protein-coding genes(PCGs), 24 tRNA genes, and 3 rRNA genes. All the tRNAs could be folded into typical cloverleaf secondary structures, with lengths of 58-88 bp. A conserved U-U-C-x-A-x2 consensus nucleotide sequence was discovered in the Ψ-loops of tRNA sequences. In addition, 447 dispersed repeats were detected, as well as found 482 RNA editing sites and 9,960 codons in the mt genome. Furthermore, a total of 27 DNA sequences with a length of 43,277 bp were transferred from the cp. to the mt genome, and eight integrated cp-derived genes (trnL-CAA, trnV-GAC, trnD-GUC, trnW-CCA, trnN-GUU, trnH-GUG, trnM-CAU, and rps7) were identified. We also obtained 1,086 hits, including 364.023 kp of nuclear genome sequences, that were transferred to the mt genome. To determine the evolutionary position of C. illinoinensis, we conducted a phylogenetic analysis of the mitogenomes of C. illinoinensis and 14 other taxa. The results strongly suggested that C. illinoinensis and Fagus sylvatica formed a single clade with 100% bootstrap support. This study sequenced comprehensive data on the C. illinoinensis mitochondrial genome and provided insights into the conserved sequences of tRNA genes, which could facilitate evolutionary research in other Carya trees in the future.

Structure, function, and assembly of PSI in thylakoid membranes of vascular plants

The photosynthetic apparatus is formed by thylakoid membrane-embedded multiprotein complexes that carry out linear electron transport in oxygenic photosynthesis. The machinery is largely conserved from cyanobacteria to land plants, and structure and function of the protein complexes involved are relatively well studied. By contrast, how the machinery is assembled in thylakoid membranes remains poorly understood. The complexes participating in photosynthetic electron transfer are composed of many proteins, pigments, and redox-active cofactors, whose temporally and spatially highly coordinated incorporation is essential to build functional mature complexes. Several proteins, jointly referred to as assembly factors, engage in the biogenesis of these complexes to bring the components together in a step-wise manner, in the right order and time. In this review, we focus on the biogenesis of the terminal protein supercomplex of the photosynthetic electron transport chain, PSI, in vascular plants. We summarize our current knowledge of the assembly process and the factors involved and describe the challenges associated with resolving the assembly pathway in molecular detail.

The plastomes of Lepismium cruciforme (Vell.) Miq and Schlumbergera truncata (Haw.) Moran reveal tribe-specific rearrangements and the first loss of the trnT-GGU gene in Cactaceae

BACKGROUND

Recent studies have revealed atypical features in the plastomes of the family Cactaceae, the largest lineage of succulent species adapted to arid and semi-arid regions. Most plastomes sequenced to date are from short-globose and cylindrical cacti, while little is known about plastomes of epiphytic cacti. Published cactus plastomes reveal reduction and complete loss of IRs, loss of genes, pseudogenization, and even degeneration of tRNA structures. Aiming to contribute with new insights into the plastid evolution of Cactaceae, particularly within the tribe Rhipsalideae, we de novo assembled and analyzed the plastomes of Lepismium cruciforme and Schlumbergera truncata, two South American epiphytic cacti.

METHODS AND RESULTS

Our data reveal many gene losses in both plastomes and the first loss of functionality of the trnT-GGU gene in Cactaceae. The trnT-GGU is a pseudogene in L. cruciforme plastome and appears to be degenerating in the tribe Rhipsalideae. Although the plastome structure is conserved among the species of the tribe Rhipsalideae, with tribe-specific rearrangements, we mapped around 200 simple sequence repeats and identified nine nucleotide polymorphism hotspots, useful to improve the phylogenetic resolutions of the Rhipsalideae. Furthermore, our analysis indicated high gene divergence and rapid evolution of RNA editing sites in plastid protein-coding genes in Cactaceae.

CONCLUSIONS

Our findings show that some characteristics of the Rhipsalideae tribe are conserved, such as plastome structure with IRs containing only the ycf2 and two tRNA genes, structural degeneration of the trnT-GGU gene and ndh complex, and lastly, pseudogenization of rpl33 and rpl23 genes, both plastid translation-related genes.

Evolution and maintenance of mtDNA gene content across eukaryotes

Across eukaryotes, most genes required for mitochondrial function have been transferred to, or otherwise acquired by, the nucleus. Encoding genes in the nucleus has many advantages. So why do mitochondria retain any genes at all? Why does the set of mtDNA genes vary so much across different species? And how do species maintain functionality in the mtDNA genes they do retain? In this review, we will discuss some possible answers to these questions, attempting a broad perspective across eukaryotes. We hope to cover some interesting features which may be less familiar from the perspective of particular species, including the ubiquity of recombination outside bilaterian animals, encrypted chainmail-like mtDNA, single genes split over multiple mtDNA chromosomes, triparental inheritance, gene transfer by grafting, gain of mtDNA recombination factors, social networks of mitochondria, and the role of mtDNA dysfunction in feeding the world. We will discuss a unifying picture where organismal ecology and gene-specific features together influence whether organism X retains mtDNA gene Y, and where ecology and development together determine which strategies, importantly including recombination, are used to maintain the mtDNA genes that are retained.

Functional conservation of GENOMES UNCOUPLED1 in plastid-to-nucleus retrograde signaling in tomato

Retrograde signaling between plastids and the nucleus is vital for chloroplast biogenesis and environmental responses. GENOMES UNCOUPLED1 (GUN1) was proposed to be a central integrator of multiple retrograde signaling pathways in the model plant Arabidopsis thaliana (Arabidopsis). However, the function of GUN1 orthologs in other plant species has not been well studied. Here, we found that many GUN1 orthologs from the Solanaceae family have a short N-terminus before the first pentatricopeptide repeat (PPR) motif which is predicted as intrinsically disordered regions (IDRs). Functional analyses of tomato (Solanum lycopersicum L.) GUN1 (SlGUN1), which does not contain N-terminal IDRs, show that it can complement the GUN phenotype of the Arabidopsis gun1 mutant (Atgun1). However, in contrast to the AtGUN1 protein, which does contain the N-terminal IDRs, the SlGUN1 protein is highly accumulated even after chloroplast biogenesis is completed, suggesting that the N-terminal IDRs may determine the stability of the GUN1 protein. Furthermore, we generated tomato Slgun1 genome-edited mutants via the CRISPR-Cas9 system. The Slgun1 mutants exhibited a typical GUN phenotype under lincomycin (Lin) or norflurazon (NF) treatment. Moreover, Slgun1 mutants are hypersensitive to low concentrations of Lin or NF. Taken together, our results suggest that, although lacking the N-terminal IDRs, SlGUN1 plays conserved roles in plastid retrograde signaling in tomato plants.

Horizontal gene transfer in eukaryotes: aligning theory with data

Horizontal gene transfer (HGT), or lateral gene transfer, is the non-sexual movement of genetic information between genomes. It has played a pronounced part in bacterial and archaeal evolution, but its role in eukaryotes is less clear. Behaviours unique to eukaryotic cells - phagocytosis and endosymbiosis - have been proposed to increase the frequency of HGT, but nuclear genomes encode fewer HGTs than bacteria and archaea. Here, I review the existing theory in the context of the growing body of data on HGT in eukaryotes, which suggests that any increased chance of acquiring new genes through phagocytosis and endosymbiosis is offset by a reduced need for these genes in eukaryotes, because selection in most eukaryotes operates on variation not readily generated by HGT.

Tandem gene duplication selected by activation of horizontally transferred gene in bacteria

Horizontal gene transfer occurs frequently in bacteria, but the mechanism driving activation and optimization of the expression of horizontally transferred genes (HTGs) in new recipient strains is not clear. Our previous study found that spontaneous tandem DNA duplication resulted in rapid activation of HTGs. Here, we took advantage of this finding to develop a novel technique for tandem gene duplication, named tandem gene duplication selected by activation of horizontally transferred gene in bacteria (TDAH), in which tandem duplication was selected by the activation of horizontally transferred selectable marker gene. TDAH construction does not contain any reported functional elements based on homologous or site-specific recombination and DNA amplification. TDAH only contains an essential selectable marker for copy number selection and 9-bp-microhomology border sequences for precise illegitimate recombination. One transformation and 3 days were enough to produce a high-copy strain, so its procedure is simple and fast. Without subsequent knockout of the endogenous recombination system, TDAH could also generate the relatively stable high-copy tandem duplication for plasmid-carried and genome-integrated DNA. TDAH also showed an excellent capacity for increase gene expression and worked well in different industrial bacteria. We also applied TDAH to select the optimal high copy number of ribA for vitamin B2 production in E. coli; the yield was improved by 3.5 times and remained stable even after 12 subcultures. TDAH is a useful tool for recombinant protein production and expression optimization of biosynthetic pathways. KEY POINTS: • We develop a novel and efficient technique (TDAH) for tandem gene duplication in bacterium. TDAH is based on the mechanism of HTG rapid activation. TDAH does not contain any reported functional elements based on homologous recombination and DNA amplification. TDAH only contains an essential selectable marker for copy number selection, so its construction and procedure are very simple and fast. • TDAH is the first reported selected and stable tandem-gene-duplication technique in which the selected high-copy plasmid-carried and genome-integrated DNA could remain stable without the subsequent knockout of recombination system. • TDAH showed an excellent capacity for regulating gene expression and worked well in different industrial bacteria, indicating it is a useful tool for recombinant protein production and expression optimization of biosynthetic pathways. • TDAH was applied to select the optimal high copy number of ribA for vitamin B2 production in E. coli; the yield was improved by 3.5-fold and remained stable even after 12 subcultures.

Nuclear genome of dinoflagellates: Size variation and insights into evolutionary mechanisms

Recent progress in high-throughput sequencing technologies has dramatically increased availability of genome data for prokaryotes and eukaryotes. Dinoflagellates have distinct chromosomes and a huge genome size, which make their genomic analysis complicated. Here, we reviewed the nuclear genomes of core dinoflagellates, focusing on the genome and cell size. Till now, the genome sizes of several dinoflagellates (more than 25) have been measured by certain methods (e.g., flow cytometry), showing a range of 3-250 pg of genomic DNA per cell. In contrast to their relatively small cell size, their genomes are huge (about 1-80 times the human haploid genome). In the present study, we collected the genome and cell size data of dinoflagellates and compared their relationships. We found that dinoflagellate genome size exhibits a positive correlation with cell size. On the other hand, we recognized that the genome size is not correlated with phylogenetic relatedness. These may be caused by genome duplication, increased gene copy number, repetitive non-coding DNA, transposon expansion, horizontal gene transfer, organelle-to-nucleus gene transfer, and/or mRNA reintegration into the genome. Ultimate verification of these factors as potential causative mechanisms would require sequencing of more dinoflagellate genomes in the future.

Cryo-EM structures of the plant plastid-encoded RNA polymerase

Chloroplasts are green plastids in the cytoplasm of eukaryotic algae and plants responsible for photosynthesis. The plastid-encoded RNA polymerase (PEP) plays an essential role during chloroplast biogenesis from proplastids and functions as the predominant RNA polymerase in mature chloroplasts. The PEP-centered transcription apparatus comprises a bacterial-origin PEP core and more than a dozen eukaryotic-origin PEP-associated proteins (PAPs) encoded in the nucleus. Here, we determined the cryo-EM structures of Nicotiana tabacum (tobacco) PEP-PAP apoenzyme and PEP-PAP transcription elongation complexes at near-atomic resolutions. Our data show the PEP core adopts a typical fold as bacterial RNAP. Fifteen PAPs bind at the periphery of the PEP core, facilitate assembling the PEP-PAP supercomplex, protect the complex from oxidation damage, and likely couple gene transcription with RNA processing. Our results report the high-resolution architecture of the chloroplast transcription apparatus and provide the structural basis for the mechanistic and functional study of transcription regulation in chloroplasts.

Characteristics and Comparative Analysis of the Special-Structure (Non-Single-Circle) Mitochondrial Genome of Capsicum pubescens Ruiz & Pav

Chilean peppers, cultivated from Capsicum pubescens, are globally renowned as popular vegetable and spice crops. C. pubescens belongs to the Capsicum L. (pepper) family and is one of the five pepper cultivars grown in China. In this study, we assembled and annotated the complete mt genome of C. pubescens. We investigated several aspects of its genome, including characteristics, codon usage, RNA editing sites, repeat sequences, selective pressure, gene clusters, and phylogenetic relationships. Furthermore, we compared it with other plant mt genomes. The data we obtained will provide valuable information for studying evolutionary processes in the Capsicum genus and will assist in the functional analysis of Capsicum mitogenomes.

Targeted knockout of a conserved plant mitochondrial gene by genome editing

Fusion proteins derived from transcription activator-like effectors (TALEs) have emerged as genome editing tools for mitochondria. TALE nucleases (TALENs) have been applied to delete chimaeric reading frames and duplicated (redundant) genes but produced complex genomic rearrangements due to the absence of non-homologous end-joining. Here we report the targeted deletion of a conserved mitochondrial gene, nad9, encoding a subunit of respiratory complex I. By generating a large number of TALEN-mediated mitochondrial deletion lines, we isolated, in addition to mutants with rearranged genomes, homochondriomic mutants harbouring clean nad9 deletions. Characterization of the knockout plants revealed impaired complex I biogenesis, male sterility and defects in leaf and flower development. We show that these defects can be restored by expressing a functional Nad9 protein from the nuclear genome, thus creating a synthetic cytoplasmic male sterility system. Our data (1) demonstrate the feasibility of using genome editing to study mitochondrial gene functions by reverse genetics, (2) highlight the role of complex I in plant development and (3) provide proof-of-concept for the construction of synthetic cytoplasmic male sterility systems for hybrid breeding by genome editing.

Plastid dsRNA transgenes trigger phased small RNA-based gene silencing of nuclear-encoded genes

Plastid transformation technology has been widely used to express traits of potential commercial importance, though the technology has been limited to traits that function while sequestered in the organelle. Prior research indicates that plastid contents can escape from the organelle, suggesting a possible mechanism for engineering plastid transgenes to function in other cellular locations. To test this hypothesis, we created tobacco (Nicotiana tabacum cv. Petit Havana) plastid transformants that express a fragment of the nuclear-encoded Phytoene desaturase (PDS) gene capable of catalyzing post-transcriptional gene silencing if RNA escapes into the cytoplasm. We found multiple lines of direct evidence that plastid-encoded PDS transgenes affect nuclear PDS gene silencing: knockdown of the nuclear-encoded PDS mRNA and/or its apparent translational inhibition, biogenesis of 21-nucleotide (nt) phased small interfering RNAs (phasiRNAs), and pigment-deficient plants. Furthermore, plastid-expressed dsRNA with no cognate nuclear-encoded pairing partner also produced abundant 21-nt phasiRNAs in the cytoplasm, demonstrating that a nuclear-encoded template is not required for siRNA biogenesis. Our results indicate that RNA escape from plastids to the cytoplasm occurs generally, with functional consequences that include entry into the gene silencing pathway. Furthermore, we uncover a method to produce plastid-encoded traits with functions outside of the organelle and open additional fields of study in plastid development, compartmentalization, and small RNA biogenesis.

Nuclear DNA segments homologous to mitochondrial DNA are obstacles for detecting heteroplasmy in sugar beet (Beta vulgaris L.)

Heteroplasmy, the coexistence of multiple mitochondrial DNA (mtDNA) sequences in a cell, is well documented in plants. Next-generation sequencing technology (NGS) has made it feasible to sequence entire genomes. Thus, NGS has the potential to detect heteroplasmy; however, the methods and pitfalls in heteroplasmy detection have not been fully investigated and identified. One obstacle for heteroplasmy detection is the sequence homology between mitochondrial-, plastid-, and nuclear DNA, of which the influence of nuclear DNA segments homologous to mtDNA (numt) need to be minimized. To detect heteroplasmy, we first excluded nuclear DNA sequences of sugar beet (Beta vulgaris) line EL10 from the sugar beet mtDNA sequence. NGS reads were obtained from single plants of sugar beet lines NK-195BRmm-O and NK-291BRmm-O and mapped to the unexcluded mtDNA regions. More than 1000 sites exhibited intra-individual polymorphism as detected by genome browsing analysis. We focused on a 309-bp region where 12 intra-individual polymorphic sites were closely linked to each other. Although the existence of DNA molecules having variant alleles at the 12 sites was confirmed by PCR amplification from NK-195BRmm-O and NK-291BRmm-O, these variants were not always called by six variant-calling programs, suggesting that these programs are inappropriate for intra-individual polymorphism detection. When we changed the nuclear DNA reference, a numt absent from EL10 was found to include the 309-bp region. Genetic segregation of an F2 population from NK-195BRmm-O x NK-291BRmm-O supported the numt origin of the variant alleles. Using four references, we found that numt detection exhibited reference dependency, and extreme polymorphism of numts exists among sugar beet lines. One of the identified numts absent from EL10 is also associated with another intra-individual polymorphic site in NK-195mm-O. Our data suggest that polymorphism among numts is unexpectedly high within sugar beets, leading to confusion about the true degree of heteroplasmy.

When and why are mitochondria paternally inherited?

In contrast with nuclear genes that are passed on through both parents, mitochondrial genes are maternally inherited in most species, most of the time. The genetic conflict stemming from this transmission asymmetry is well-documented, and there is an abundance of population-genetic theory associated with it. While occasional or aberrant paternal inheritance occurs, there are only a few cases where exclusive paternal inheritance of mitochondrial genomes is the evolved state. Why this is remains poorly understood. By examining commonalities between species with exclusive paternal inheritance, we discuss what they may tell us about the evolutionary forces influencing mitochondrial inheritance patterns. We end by discussing recent technological advances that make exploring the causes and consequences of paternal inheritance feasible.

The garden asparagus (Asparagus officinalis L.) mitochondrial genome revealed rich sequence variation throughout whole sequencing data

Garden asparagus (Asparagus officinalis L.) is a horticultural crop with high nutritional and medical value, considered an ideal plant for sex determination research among many dioecious plants, whose genomic information can support genetic analysis and breeding programs. In this research, the entire mitochondrial genome of A. officinalis was sequenced, annotated and assembled using a mixed Illumina and PacBio data. The garden asparagus circular mitochondrial genome measures 492,062 bp with a GC value of 45.9%. Thirty-six protein-coding genes, 17 tRNA and 6 rRNA genes were annotated, among which 8 protein-coding genes contained 16 introns. In addition, 254 SSRs with 10 complete tandem repeats and 293 non-tandem repeats were identified. It was found that the codons of edited sites located in the amino acids showed a leucine-formation trend, and RNA editing sites mainly caused the mutual transformation of amino acids with the same properties. Furthermore, 72 sequence fragments accounting for 20,240 bp, presentating 4.11% of the whole mitochondrial genome, were observed to migrate from chloroplast to mitochondrial genome of A. officinalis. The phylogenetic analysis showed that the closest genetic relationship between A. officinalis with onion (Allium cepa) inside the Liliaceae family. Our results demonstrated that high percentage of protein-coding genes had evolutionary conservative properties, with Ka/Ks values less than 1. Therefore, this study provides a high-quality garden asparagus mitochondrial genome, useful to promote better understanding of gene exchange between organelle genomes.

Chloroplast proteostasis: A story of birth, life, and death

Protein homeostasis (proteostasis) is a dynamic balance of protein synthesis and degradation. Because of the endosymbiotic origin of chloroplasts and the massive transfer of their genetic information to the nucleus of the host cell, many protein complexes in the chloroplasts are constituted from subunits encoded by both genomes. Hence, the proper function of chloroplasts relies on the coordinated expression of chloroplast- and nucleus-encoded genes. The biogenesis and maintenance of chloroplast proteostasis are dependent on synthesis of chloroplast-encoded proteins, import of nucleus-encoded chloroplast proteins from the cytosol, and clearance of damaged or otherwise undesired "old" proteins. This review focuses on the regulation of chloroplast proteostasis, its interaction with proteostasis of the cytosol, and its retrograde control over nuclear gene expression. We also discuss significant issues and perspectives for future studies and potential applications for improving the photosynthetic performance and stress tolerance of crops.

Evolutionary Overview of Aquaporin Superfamily

Aquaporins (AQPs) are present not only in three domains of life, bacteria, eukaryotes, and archaea, but also in viruses. With the accumulating arrays of AQP superfamily, the evolutional relationship has attracted much attention with multiple publications on "the genome-wide identification and phylogenetic analysis" of AQP superfamily. A pair of NPA boxes forming a pore is highly conserved throughout the evolution and renders key residues for the classification of AQP superfamily into four groups: AQP1-like, AQP3-like, AQP8-like, and AQP11-like. The complexity of AQP family has mostly been achieved in nematodes and subsequent evolution has been directed toward increasing the number of AQPs through whole-genome duplications (WGDs) to extend the tissue specific expression and regulation. The discovery of the intracellular AQP (iAQP: AQP8-like and AQP11-like) and substrate transports by the plasma membrane AQP (pAQP: AQP1-like and AQP3-like) have accelerated the AQP research much more toward the transport of substrates with complex profiles. This evolutionary overview based on a simple classification of AQPs into four subfamilies will provide putative structural, functional, and localization information and insights into the role of AQP as well as clues to understand the complex diversity of AQP superfamily.

Mitochondrial signal transduction

The analogy of mitochondria as powerhouses has expired. Mitochondria are living, dynamic, maternally inherited, energy-transforming, biosynthetic, and signaling organelles that actively transduce biological information. We argue that mitochondria are the processor of the cell, and together with the nucleus and other organelles they constitute the mitochondrial information processing system (MIPS). In a three-step process, mitochondria (1) sense and respond to both endogenous and environmental inputs through morphological and functional remodeling; (2) integrate information through dynamic, network-based physical interactions and diffusion mechanisms; and (3) produce output signals that tune the functions of other organelles and systemically regulate physiology. This input-to-output transformation allows mitochondria to transduce metabolic, biochemical, neuroendocrine, and other local or systemic signals that enhance organismal adaptation. An explicit focus on mitochondrial signal transduction emphasizes the role of communication in mitochondrial biology. This framework also opens new avenues to understand how mitochondria mediate inter-organ processes underlying human health.

Adaptation to an Intracellular Lifestyle by a Nitrogen-Fixing, Heterocyst-Forming Cyanobacterial Endosymbiont of a Diatom

The symbiosis between the diatom Hemiaulus hauckii and the heterocyst-forming cyanobacterium Richelia intracellularis makes an important contribution to new production in the world's oceans, but its study is limited by short-term survival in the laboratory. In this symbiosis, R. intracellularis fixes atmospheric dinitrogen in the heterocyst and provides H. hauckii with fixed nitrogen. Here, we conducted an electron microscopy study of H. hauckii and found that the filaments of the R. intracellularis symbiont, typically composed of one terminal heterocyst and three or four vegetative cells, are located in the diatom's cytoplasm not enclosed by a host membrane. A second prokaryotic cell was also detected in the cytoplasm of H. hauckii, but observations were infrequent. The heterocysts of R. intracellularis differ from those of free-living heterocyst-forming cyanobacteria in that the specific components of the heterocyst envelope seem to be located in the periplasmic space instead of outside the outer membrane. This specialized arrangement of the heterocyst envelope and a possible association of the cyanobacterium with oxygen-respiring mitochondria may be important for protection of the nitrogen-fixing enzyme, nitrogenase, from photosynthetically produced oxygen. The cell envelope of the vegetative cells of R. intracellularis contained numerous membrane vesicles that resemble the outer-inner membrane vesicles of Gram-negative bacteria. These vesicles can export cytoplasmic material from the bacterial cell and, therefore, may represent a vehicle for transfer of fixed nitrogen from R. intracellularis to the diatom's cytoplasm. The specific morphological features of R. intracellularis described here, together with its known streamlined genome, likely represent specific adaptations of this cyanobacterium to an intracellular lifestyle.

Born in the mitochondrion and raised in the nucleus: evolution of a novel tandem repeat family in Medicago polymorpha (Fabaceae)

Plant nuclear genomes harbor sequence elements derived from the organelles (mitochondrion and plastid) through intracellular gene transfer (IGT). Nuclear genomes also show a dramatic range of repeat content, suggesting that any sequence can be readily amplified. These two aspects of plant nuclear genomes are well recognized but have rarely been linked. Through investigation of 31 Medicago taxa we detected exceptionally high post-IGT amplification of mitochondrial (mt) DNA sequences containing rps10 in the nuclear genome of Medicago polymorpha and closely related species. The amplified sequences were characterized as tandem arrays of five distinct repeat motifs (2157, 1064, 987, 971, and 587 bp) that have diverged from the mt genome (mitogenome) in the M. polymorpha nuclear genome. The mt rps10-like arrays were identified in seven loci (six intergenic and one telomeric) of the nuclear chromosome assemblies and were the most abundant tandem repeat family, representing 1.6-3.0% of total genomic DNA, a value approximately three-fold greater than the entire mitogenome in M. polymorpha. Compared to a typical mt gene, the mt rps10-like sequence coverage level was 691.5-7198-fold higher in M. polymorpha and closely related species. In addition to the post-IGT amplification, our analysis identified the canonical telomeric repeat and the species-specific satellite arrays that are likely attributable to an ancestral chromosomal fusion in M. polymorpha. A possible relationship between chromosomal instability and the mt rps10-like tandem repeat family in the M. polymorpha clade is discussed.

Comprehensive analysis of plastid gene expression during fruit development and ripening of kiwifruit

Global survey of plastid gene expression during fruit ripening in kiwifruit provides cis-elements for the future engineering of the plastid genome of kiwifruit. A limitation in the application of plastid biotechnology for molecular farming is the low-level expression of transgenes in non-green plastids compared with photosynthetically active chloroplasts. Unlike other fruits, not all chloroplasts are transformed into chromoplasts during ripening of red-fleshed kiwifruit (Actinidia chinensis cv. Hongyang) fruits, which may make kiwifruit an ideal horticultural plant for recombinant protein production by plastid engineering. To identify cis-elements potentially triggering high-level transgene expression in edible tissues of the 'Hongyang' kiwifruit, here we report a comprehensive analysis of kiwifruit plastid gene transcription in green leaves and fruits at three different developmental stages. While transcripts of a few photosynthesis-related genes and most genetic system genes were substantially upregulated in green fruits compared with leaves, nearly all plastid genes were significantly downregulated at the RNA level during fruit development. Expression of a few genes remained unchanged, including psbA, the gene encoding the D1 polypeptide of photosystem II. However, PsbA protein accumulation decreased continuously during chloroplast-to-chromoplast differentiation. Analysis of post-transcriptional steps in mRNA maturation, including intron splicing and RNA editing, revealed that splicing and editing may contribute to regulation of plastid gene expression. Altogether, 40 RNA editing sites were verified, and 5 of them were newly discovered. Taken together, this study has generated a valuable resource for the analysis of plastid gene expression and provides cis-elements for future efforts to engineer the plastid genome of kiwifruit.

Phylogenomic analyses highlight innovation and introgression in the continental radiations of Fagaceae across the Northern Hemisphere

Northern Hemisphere forests changed drastically in the early Eocene with the diversification of the oak family (Fagaceae). Cooling climates over the next 20 million years fostered the spread of temperate biomes that became increasingly dominated by oaks and their chestnut relatives. Here we use phylogenomic analyses of nuclear and plastid genomes to investigate the timing and pattern of major macroevolutionary events and ancient genome-wide signatures of hybridization across Fagaceae. Innovation related to seed dispersal is implicated in triggering waves of continental radiations beginning with the rapid diversification of major lineages and resulting in unparalleled transformation of forest dynamics within 15 million years following the K-Pg extinction. We detect introgression at multiple time scales, including ancient events predating the origination of genus-level diversity. As oak lineages moved into newly available temperate habitats in the early Miocene, secondary contact between previously isolated species occurred. This resulted in adaptive introgression, which may have further amplified the diversification of white oaks across Eurasia.

Targeted introduction of heritable point mutations into the plant mitochondrial genome

The development of technologies for the genetic manipulation of mitochondrial genomes remains a major challenge. Here we report a method for the targeted introduction of mutations into plant mitochondrial DNA (mtDNA) that we refer to as transcription activator-like effector nuclease (TALEN) gene-drive mutagenesis (GDM), or TALEN-GDM. The method combines TALEN-induced site-specific cleavage of the mtDNA with selection for mutations that confer resistance to the TALEN cut. Applying TALEN-GDM to the tobacco mitochondrial nad9 gene, we isolated a large set of mutants carrying single amino acid substitutions in the Nad9 protein. The mutants could be purified to homochondriomy and stably inherited their edited mtDNA in the expected maternal fashion. TALEN-GDM induces both transitions and transversions, and can access most nucleotide positions within the TALEN binding site. Our work provides an efficient method for targeted mitochondrial genome editing that produces genetically stable, homochondriomic and fertile plants with specific point mutations in their mtDNA.

Ecological selection of bacterial taxa with larger genome sizes in response to polycyclic aromatic hydrocarbons stress

Polycyclic aromatic hydrocarbons (PAHs) are ubiquitous priority pollutants that cause great damage to the natural environment and health. Average genome size in a community is critical for shedding light on microbiome's functional response to pollution stress within an environment. Here, microcosms under different concentrations were performed to evaluate the selection of PAHs stress on the average genome size in a community. We found the distinct communities of significantly larger genome size with the increase of PAHs concentration gradients in soils, and consistent trends were discovered in soils at different latitudes. The abundance of Proteobacteria and Deinococcus-Thermus with relatively larger genomes increased along with PAHs stress and well adapted to polluted environments. In contrast, the abundance of Patescibacteria with a highly streamlined and smaller genome decreased, implying complex interactions between environmental selection and functional fitness resulted in bacteria with larger genomes becoming more abundant. Moreover, we confirmed the increased capacity for horizontal transfer of degrading genes between communities by showing an increased connection number per node positively related to the nidA gene along the concentration gradients in the co-occurrence network. Our findings suggest PAHs tend to select bacterial taxa with larger genome sizes, with significant consequences for community stability and potential biodegradation strategies.

Assembly and comparative analysis of the first complete mitochondrial genome of Acer truncatum Bunge: a woody oil-tree species producing nervonic acid

BACKGROUND

Acer truncatum (purpleblow maple) is a woody tree species that produces seeds with high levels of valuable fatty acids (especially nervonic acid). The species is admired as a landscape plant with high developmental prospects and scientific research value. The A. truncatum chloroplast genome has recently been reported; however, the mitochondrial genome (mitogenome) is still unexplored.

RESULTS

We characterized the A. truncatum mitogenome, which was assembled using reads from PacBio and Illumina sequencing platforms, performed a comparative analysis against different species of Acer. The circular mitogenome of A. truncatum has a length of 791,052 bp, with a base composition of 27.11% A, 27.21% T, 22.79% G, and 22.89% C. The A. truncatum mitogenome contains 62 genes, including 35 protein-coding genes, 23 tRNA genes and 4 rRNA genes. We also examined codon usage, sequence repeats, RNA editing and selective pressure in the A. truncatum mitogenome. To determine the evolutionary and taxonomic status of A. truncatum, we conducted a phylogenetic analysis based on the mitogenomes of A. truncatum and 25 other taxa. In addition, the gene migration from chloroplast and nuclear genomes to the mitogenome were analyzed. Finally, we developed a novel NAD1 intron indel marker for distinguishing several Acer species.

CONCLUSIONS

In this study, we assembled and annotated the mitogenome of A. truncatum, a woody oil-tree species producing nervonic acid. The results of our analyses provide comprehensive information on the A. truncatum mitogenome, which would facilitate evolutionary research and molecular barcoding in Acer.

Cryptic promoter activation occurs by at least two different mechanisms in the Arabidopsis genome

In gene-trap screening of plant genomes, promoterless reporter constructs are often expressed without trapping of annotated gene promoters. The molecular basis of this phenomenon, which has been interpreted as the trapping of cryptic promoters, is poorly understood. Here, we found that cryptic promoter activation occurs by at least two different mechanisms using Arabidopsis gene-trap lines in which a firefly luciferase (LUC) open reading frame (ORF) without an apparent promoter sequence was expressed from intergenic regions: one mechanism is 'cryptic promoter capturing', in which the LUC ORF captured pre-existing promoter-like chromatin marked by H3K4me3 and H2A.Z, and the other is 'promoter de novo origination', in which the promoter chromatin was newly formed near the 5' end of the inserted LUC ORF. The latter finding raises a question as to how the inserted LUC ORF sequence is involved in this phenomenon. To examine this, we performed a model experiment with chimeric LUC genes in transgenic plants. Using Arabidopsis psaH1 promoter-LUC constructs, we found that the functional core promoter region, where transcription start sites (TSSs) occur, cannot simply be determined by the upstream nor core promoter sequences; rather, its positioning proximal to the inserted LUC ORF sequence was more critical. This result suggests that the insertion of the coding sequence alters the local distribution of TSSs in the plant genome. The possible impact of the two types of cryptic promoter activation mechanisms on plant genome evolution and endosymbiotic gene transfer is discussed.

Engineering Metabolism in Nicotiana Species: A Promising Future

Molecular farming intends to use crop plants as biofactories for high value-added compounds following application of a wide range of biotechnological tools. In particular, the conversion of nonfood crops into efficient biofactories is expected to be a strong asset in the development of a sustainable bioeconomy. The 'nonfood' status combined with the high metabolic versatility and the capacity of high-yield cultivation highlight the plant genus Nicotiana as one of the most appropriate 'chassis' for molecular farming. Nicotiana species are a rich source of valuable industrial, active pharmaceutical ingredients and nutritional compounds, synthesized from highly complex biosynthetic networks. Here, we review and discuss approaches currently used to design enriched Nicotiana species for molecular farming using new plant breeding techniques (NPBTs).

Correction of frameshift mutations in the atpB gene by translational recoding in chloroplasts of Oenothera and tobacco

Translational recoding, also known as ribosomal frameshifting, is a process that causes ribosome slippage along the messenger RNA, thereby changing the amino acid sequence of the synthesized protein. Whether the chloroplast employs recoding is unknown. I-iota, a plastome mutant of Oenothera (evening primrose), carries a single adenine insertion in an oligoA stretch [11A] of the atpB coding region (encoding the β-subunit of the ATP synthase). The mutation is expected to cause synthesis of a truncated, nonfunctional protein. We report that a full-length AtpB protein is detectable in I-iota leaves, suggesting operation of a recoding mechanism. To characterize the phenomenon, we generated transplastomic tobacco lines in which the atpB reading frame was altered by insertions or deletions in the oligoA motif. We observed that insertion of two adenines was more efficiently corrected than insertion of a single adenine, or deletion of one or two adenines. We further show that homopolymeric composition of the oligoA stretch is essential for recoding, as an additional replacement of AAA lysine codon by AAG resulted in an albino phenotype. Our work provides evidence for the operation of translational recoding in chloroplasts. Recoding enables correction of frameshift mutations and can restore photoautotrophic growth in the presence of a mutation that otherwise would be lethal.

GUN control in retrograde signaling: How GENOMES UNCOUPLED proteins adjust nuclear gene expression to plastid biogenesis

Communication between cellular compartments is vital for development and environmental adaptation. Signals emanating from organelles, so-called retrograde signals, coordinate nuclear gene expression with the developmental stage and/or the functional status of the organelle. Plastids (best known in their green photosynthesizing differentiated form, the chloroplasts) are the primary energy-producing compartment of plant cells, and the site for the biosynthesis of many metabolites, including fatty acids, amino acids, nucleotides, isoprenoids, tetrapyrroles, vitamins, and phytohormone precursors. Signals derived from plastids regulate the accumulation of a large set of nucleus-encoded proteins, many of which localize to plastids. A set of mutants defective in retrograde signaling (genomes uncoupled, or gun) was isolated over 25 years ago. While most GUN genes act in tetrapyrrole biosynthesis, resolving the molecular function of GUN1, the proposed integrator of multiple retrograde signals, has turned out to be particularly challenging. Based on its amino acid sequence, GUN1 was initially predicted to be a plastid-localized nucleic acid-binding protein. Only recently, mechanistic information on the function of GUN1 has been obtained, pointing to a role in plastid protein homeostasis. This review article summarizes our current understanding of GUN-related retrograde signaling and provides a critical appraisal of the various proposed roles for GUNs and their respective pathways.

Molecular evolution and SSRs analysis based on the chloroplast genome of Callitropsis funebris

Chloroplast genome sequences have been used to understand evolutionary events and to infer efficiently phylogenetic relationships. Callitropsis funebris (Cupressaceae) is an endemic species in China. Its phylogenetic position is controversial due to morphological characters similar to those of Cupressus, Callitropsis, and Chamaecyparis. This study used next-generation sequencing technology to sequence the complete chloroplast genome of Ca. funebris and then constructed the phylogenetic relationship between Ca. funebris and its related species based on a variety of data sets and methods. Simple sequence repeats (SSRs) and adaptive evolution analysis were also conducted. Our results showed that the monophyletic branch consisting of Ca. funebris and Cupressus tonkinensis is a sister to Cupressus, while Callitropsis is not monophyletic; Ca. nootkatensis and Ca. vietnamensis are nested in turn at the base of the monophyletic group Hesperocyparis. The statistical results of SSRs supported the closest relationship between Ca. funebris and Cupressus. By performing adaptive evolution analysis under the phylogenetic background of Cupressales, the Branch model detected three genes and the Site model detected 10 genes under positive selection; and the Branch-Site model uncovered that rpoA has experienced positive selection in the Ca. funebries branch. Molecular analysis from the chloroplast genome highly supported that Ca. funebris is at the base of Cupressus. Of note, SSR features were found to be able to shed some light on phylogenetic relationships. In short, this chloroplast genomic study has provided new insights into the phylogeny of Ca. funebris and revealed multiple chloroplast genes possibly undergoing adaptive evolution.

Knockdown of the plastid-encoded acetyl-CoA carboxylase gene uncovers functions in metabolism and development

De novo fatty acid biosynthesis in plants relies on a prokaryotic-type acetyl-CoA carboxylase (ACCase) that resides in the plastid compartment. The enzyme is composed of four subunits, one of which is encoded in the plastid genome, whereas the other three subunits are encoded by nuclear genes. The plastid gene (accD) encodes the β-carboxyltransferase subunit of ACCase and is essential for cell viability. To facilitate the functional analysis of accD, we pursued a transplastomic knockdown strategy in tobacco (Nicotiana tabacum). By introducing point mutations into the translational start codon of accD, we obtained stable transplastomic lines with altered ACCase activity. Replacement of the standard initiator codon AUG with UUG strongly reduced AccD expression, whereas replacement with GUG had no detectable effects. AccD knockdown mutants displayed reduced ACCase activity, which resulted in changes in the levels of many but not all species of cellular lipids. Limiting fatty acid availability caused a wide range of macroscopic, microscopic, and biochemical phenotypes, including impaired chloroplast division, reduced seed set, and altered storage metabolism. Finally, while the mutants displayed reduced growth under photoautotrophic conditions, they showed exaggerated growth under heterotrophic conditions, thus uncovering an unexpected antagonistic role of AccD activity in autotrophic and heterotrophic growth.

A highly conserved mechanism for the detoxification and assimilation of the toxic phytoproduct L-azetidine-2-carboxylic acid in Aspergillus nidulans

Plants produce toxic secondary metabolites as defense mechanisms against phytopathogenic microorganisms and predators. L-azetidine-2-carboxylic acid (AZC), a toxic proline analogue produced by members of the Liliaceae and Agavaciae families, is part of such a mechanism. AZC causes a broad range of toxic, inflammatory and degenerative abnormalities in human and animal cells, while it is known that some microorganisms have evolved specialized strategies for AZC resistance. However, the mechanisms underlying these processes are poorly understood. Here, we identify a widespread mechanism for AZC resistance in fungi. We show that the filamentous ascomycete Aspergillus nidulans is able to not only resist AZC toxicity but also utilize it as a nitrogen source via GABA catabolism and the action of the AzhA hydrolase, a member of a large superfamily of detoxifying enzymes, the haloacid dehalogenase-like hydrolase (HAD) superfamily. This detoxification process is further assisted by the NgnA acetyltransferase, orthologue of Mpr1 of Saccharomyces cerevisiae. We additionally show that heterologous expression of AzhA protein can complement the AZC sensitivity of S. cerevisiae. Furthermore, a detailed phylogenetic analysis of AzhA homologues in Fungi, Archaea and Bacteria is provided. Overall, our results unravel a widespread mechanism for AZC resistance among microorganisms, including important human and plant pathogens.

The Genomic Impact of Mycoheterotrophy in Orchids

Mycoheterotrophic plants have lost the ability to photosynthesize and obtain essential mineral and organic nutrients from associated soil fungi. Despite involving radical changes in life history traits and ecological requirements, the transition from autotrophy to mycoheterotrophy has occurred independently in many major lineages of land plants, most frequently in Orchidaceae. Yet the molecular mechanisms underlying this shift are still poorly understood. A comparison of the transcriptomes of Epipogium aphyllum and Neottia nidus-avis, two completely mycoheterotrophic orchids, to other autotrophic and mycoheterotrophic orchids showed the unexpected retention of several genes associated with photosynthetic activities. In addition to these selected retentions, the analysis of their expression profiles showed that many orthologs had inverted underground/aboveground expression ratios compared to autotrophic species. Fatty acid and amino acid biosynthesis as well as primary cell wall metabolism were among the pathways most impacted by this expression reprogramming. Our study suggests that the shift in nutritional mode from autotrophy to mycoheterotrophy remodeled the architecture of the plant metabolism but was associated primarily with function losses rather than metabolic innovations.

Horizontal genome transfer by cell-to-cell travel of whole organelles

Recent work has revealed that both plants and animals transfer genomes between cells. In plants, horizontal transfer of entire plastid, mitochondrial, or nuclear genomes between species generates new combinations of nuclear and organellar genomes, or produces novel species that are allopolyploid. The mechanisms of genome transfer between cells are unknown. Here, we used grafting to identify the mechanisms involved in plastid genome transfer from plant to plant. We show that during proliferation of wound-induced callus, plastids dedifferentiate into small, highly motile, amoeboid organelles. Simultaneously, new intercellular connections emerge by localized cell wall disintegration, forming connective pores through which amoeboid plastids move into neighboring cells. Our work uncovers a pathway of organelle movement from cell to cell and provides a mechanistic framework for horizontal genome transfer.

Pas de Trois: An Overview of Penta-, Tetra-, and Octo-Tricopeptide Repeat Proteins From Chlamydomonas reinhardtii and Their Role in Chloroplast Gene Expression

Penta-, Tetra-, and Octo-tricopeptide repeat (PPR, TPR, and OPR) proteins are nucleus-encoded proteins composed of tandem repeats of 35, 34, and 38-40 amino acids, respectively. They form helix-turn-helix structures that interact with mRNA or other proteins and participate in RNA stabilization, processing, maturation, and act as translation enhancers of chloroplast and mitochondrial mRNAs. These helical repeat proteins are unevenly present in plants and algae. While PPR proteins are more abundant in plants than in algae, OPR proteins are more abundant in algae. In Arabidopsis, maize, and rice there have been 450, 661, and 477 PPR proteins identified, respectively, which contrasts with only 14 PPR proteins identified in Chlamydomonas reinhardtii. Likewise, more than 120 OPR proteins members have been predicted from the nuclear genome of C. reinhardtii and only one has been identified in Arabidopsis thaliana. Due to their abundance in land plants, PPR proteins have been largely characterized making it possible to elucidate their RNA-binding code. This has even allowed researchers to generate engineered PPR proteins with defined affinity to a particular target, which has served as the basis to develop tools for gene expression in biotechnological applications. However, fine elucidation of the helical repeat proteins code in Chlamydomonas is a pending task. In this review, we summarize the current knowledge on the role PPR, TPR, and OPR proteins play in chloroplast gene expression in the green algae C. reinhardtii, pointing to relevant similarities and differences with their counterparts in plants. We also recapitulate on how these proteins have been engineered and shown to serve as mRNA regulatory factors for biotechnological applications in plants and how this could be used as a starting point for applications in algae.

Regulation of Mitochondrial Respiratory Chain Complex Levels, Organization, and Function by Arginyltransferase 1

Arginyltransferase 1 (ATE1) is an evolutionary-conserved eukaryotic protein that localizes to the cytosol and nucleus. It is the only known enzyme in metazoans and fungi that catalyzes posttranslational arginylation. Lack of arginylation has been linked to an array of human disorders, including cancer, by altering the response to stress and the regulation of metabolism and apoptosis. Although mitochondria play relevant roles in these processes in health and disease, a causal relationship between ATE1 activity and mitochondrial biology has yet to be established. Here, we report a phylogenetic analysis that traces the roots of ATE1 to alpha-proteobacteria, the mitochondrion microbial ancestor. We then demonstrate that a small fraction of ATE1 localizes within mitochondria. Furthermore, the absence of ATE1 influences the levels, organization, and function of respiratory chain complexes in mouse cells. Specifically, ATE1-KO mouse embryonic fibroblasts have increased levels of respiratory supercomplexes I+III₂+IV_n. However, they have decreased mitochondrial respiration owing to severely lowered complex II levels, which leads to accumulation of succinate and downstream metabolic effects. Taken together, our findings establish a novel pathway for mitochondrial function regulation that might explain ATE1-dependent effects in various disease conditions, including cancer and aging, in which metabolic shifts are part of the pathogenic or deleterious underlying mechanism.

Mitochondrial pathways in human health and aging

Mitochondria are eukaryotic organelles known best for their roles in energy production and metabolism. While often thought of as simply the 'powerhouse of the cell,' these organelles participate in a variety of critical cellular processes including reactive oxygen species (ROS) production, regulation of programmed cell death, modulation of inter- and intracellular nutrient signaling pathways, and maintenance of cellular proteostasis. Disrupted mitochondrial function is a hallmark of eukaryotic aging, and mitochondrial dysfunction has been reported to play a role in many aging-related diseases. While mitochondria are major players in human diseases, significant questions remain regarding their precise mechanistic role. In this review, we detail mechanisms by which mitochondrial dysfunction participate in disease and aging based on findings from model organisms and human genetics studies.

Phylogenetic and evolutionary features of the plastome of Tropaeolum pentaphyllum Lam. (Tropaeolaceae)

Complete plastome sequence of Tropaeolum pentaphyllum revealed molecular markers, hotspots of nucleotide polymorphism, RNA editing sites and phylogenetic aspects Tropaeolaceae Juss. ex DC. comprises approximately 95 species across North and South Americas. Tropaeolum pentaphyllum Lam. is an unconventional and endangered species with occurrence in some countries of South America. Although this species presents nutritional, medicinal and ornamental uses, genetic studies involving natural populations or promising genotypes are practically non-existent. Here, we report the nucleotide sequence of T. pentaphyllum plastome. It represents the first complete plastome sequence of the family Tropaeolaceae to be fully sequenced and analyzed in detail. The sequencing data revealed that the T. pentaphyllum plastome is highly similar to the plastomes of other Brassicales. Notwithstanding, our analyses detected some specific features concerning events of IR expansion and structural changes in some genes such as matK, rpoA, and rpoC2. We also detected 251 SSR loci, nine hotspots of nucleotide polymorphism, and two specific RNA editing sites in the plastome of T. pentaphyllum. Moreover, plastid phylogenomic inference indicated a closed relationship between the families Tropaeolaceae and Akaniaceae, which formed a sister group to Moringaceae-Caricaceae. Finally, our data bring new molecular markers and evolutionary features to be applied in the natural population, germplasm collection, and genotype selection aiming conservation, genetic diversity evaluation, and exploitation of this endangered species.

The plastid NAD(P)H dehydrogenase-like complex: structure, function and evolutionary dynamics

The thylakoid NAD(P)H dehydrogenase-like (NDH) complex is a large protein complex that reduces plastoquinone and pumps protons into the lumen generating protonmotive force. In plants, the complex consists of both nuclear and chloroplast-encoded subunits. Despite its perceived importance for stress tolerance and ATP generation, chloroplast-encoded NDH subunits have been lost numerous times during evolution in species occupying seemingly unrelated environmental niches. We have generated a phylogenetic tree that reveals independent losses in multiple phylogenetic lineages, and we use this tree as a reference to discuss possible evolutionary contexts that may have relaxed selective pressure for retention of ndh genes. While we are still yet unable to pinpoint a singular specific lifestyle that negates the need for NDH, we are able to rule out several long-standing explanations. In light of this, we discuss the biochemical changes that would be required for the chloroplast to dispense with NDH functionality with regards to known and proposed NDH-related reactions.

phylostratr: a framework for phylostratigraphy

MOTIVATION

The goal of phylostratigraphy is to infer the evolutionary origin of each gene in an organism. This is done by searching for homologs within increasingly broad clades. The deepest clade that contains a homolog of the protein(s) encoded by a gene is that gene's phylostratum.

RESULTS

We have created a general R-based framework, phylostratr, to estimate the phylostratum of every gene in a species. The program fully automates analysis: selecting species for balanced representation, retrieving sequences, building databases, inferring phylostrata and returning diagnostics. Key diagnostics include: detection of genes with inferred homologs in old clades, but not intermediate ones; proteome quality assessments; false-positive diagnostics, and checks for missing organellar genomes. phylostratr allows extensive customization and systematic comparisons of the influence of analysis parameters or genomes on phylostrata inference. A user may: modify the automatically generated clade tree or use their own tree; provide custom sequences in place of those automatically retrieved from UniProt; replace BLAST with an alternative algorithm; or tailor the method and sensitivity of the homology inference classifier. We show the utility of phylostratr through case studies in Arabidopsis thaliana and Saccharomyces cerevisiae.

AVAILABILITY AND IMPLEMENTATION

Source code available at https://github.com/arendsee/phylostratr.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The draft mitochondrial genome of Magnolia biondii and mitochondrial phylogenomics of angiosperms

The mitochondrial genomes of flowering plants are well known for their large size, variable coding-gene set and fluid genome structure. The available mitochondrial genomes of the early angiosperms show extreme genetic diversity in genome size, structure, and sequences, such as rampant HGTs in Amborella mt genome, numerous repeated sequences in Nymphaea mt genome, and conserved gene evolution in Liriodendron mt genome. However, currently available early angiosperm mt genomes are still limited, hampering us from obtaining an overall picture of the mitogenomic evolution in angiosperms. Here we sequenced and assembled the draft mitochondrial genome of Magnolia biondii Pamp. from Magnoliaceae (magnoliids) using Oxford Nanopore sequencing technology. We recovered a single linear mitochondrial contig of 967,100 bp with an average read coverage of 122 × and a GC content of 46.6%. This draft mitochondrial genome contains a rich 64-gene set, similar to those of Liriodendron and Nymphaea, including 41 protein-coding genes, 20 tRNAs, and 3 rRNAs. Twenty cis-spliced and five trans-spliced introns break ten protein-coding genes in the Magnolia mt genome. Repeated sequences account for 27% of the draft genome, with 17 out of the 1,145 repeats showing recombination evidence. Although partially assembled, the approximately 1-Mb mt genome of Magnolia is still among the largest in angiosperms, which is possibly due to the expansion of repeated sequences, retention of ancestral mtDNAs, and the incorporation of nuclear genome sequences. Mitochondrial phylogenomic analysis of the concatenated datasets of 38 conserved protein-coding genes from 91 representatives of angiosperm species supports the sister relationship of magnoliids with monocots and eudicots, which is congruent with plastid evidence.

Mitonuclear genomics and aging

Our cells operate based on two distinct genomes that are enclosed in the nucleus and mitochondria. The mitochondrial genome presumably originates from endosymbiotic bacteria. With time, a large portion of the original genes in the bacterial genome is considered to have been lost or transferred to the nuclear genome, leaving a reduced 16.5 Kb circular mitochondrial DNA (mtDNA). Traditionally only 37 genes, including 13 proteins, were thought to be encoded within mtDNA, its genetic repertoire is expanding with the identification of mitochondrial-derived peptides (MDPs). The biology of aging has been largely unveiled to be regulated by genes that are encoded in the nuclear genome, whereas the mitochondrial genome remained more cryptic. However, recent studies position mitochondria and mtDNA as an important counterpart to the nuclear genome, whereby the two organelles constantly regulate each other. Thus, the genomic network that regulates lifespan and/or healthspan is likely constituted by two unique, yet co-evolved, genomes. Here, we will discuss aspects of mitochondrial biology, especially mitochondrial communication that may add substantial momentum to aging research by accounting for both mitonuclear genomes to more comprehensively and inclusively map the genetic and molecular networks that govern aging and age-related diseases.

The evolution of metabolism: How to test evolutionary hypotheses at the genomic level

The origin of primordial metabolism and its expansion to form the metabolic networks extant today represent excellent systems to study the impact of natural selection and the potential adaptive role of novel compounds. Here we present the current hypotheses made on the origin of life and ancestral metabolism and present the theories and mechanisms by which the large chemical diversity of plants might have emerged along evolution. In particular, we provide a survey of statistical methods that can be used to detect signatures of selection at the gene and population level, and discuss potential and limits of these methods for investigating patterns of molecular adaptation in plant metabolism.

Mitochondrial genomes of the early land plant lineage liverworts (Marchantiophyta): conserved genome structure, and ongoing low frequency recombination

BACKGROUND

In contrast to the highly labile mitochondrial (mt) genomes of vascular plants, the architecture and composition of mt genomes within the main lineages of bryophytes appear stable and invariant. The available mt genomes of 18 liverwort accessions representing nine genera and five orders are syntenous except for Gymnomitrion concinnatum whose genome is characterized by two rearrangements. Here, we expanded the number of assembled liverwort mt genomes to 47, broadening the sampling to 31 genera and 10 orders spanning much of the phylogenetic breadth of liverworts to further test whether the evolution of the liverwort mitogenome is overall static.

RESULTS

Liverwort mt genomes range in size from 147 Kb in Jungermanniales (clade B) to 185 Kb in Marchantiopsida, mainly due to the size variation of intergenic spacers and number of introns. All newly assembled liverwort mt genomes hold a conserved set of genes, but vary considerably in their intron content. The loss of introns in liverwort mt genomes might be explained by localized retroprocessing events. Liverwort mt genomes are strictly syntenous in genome structure with no structural variant detected in our newly assembled mt genomes. However, by screening the paired-end reads, we do find rare cases of recombination, which means multiple concurrent genome structures may exist in the vegetative tissues of liverworts. Our phylogenetic analyses of the nuclear encoded double stand break repair protein families revealed liverwort-specific subfamilies expansions.

CONCLUSIONS

The low repeat recombination level, selection, along with the intensified nuclear surveillance, might together shape the structural evolution of liverwort mt genomes.

Intergenomic gene transfer in diploid and allopolyploid Gossypium

BACKGROUND

Intergenomic gene transfer (IGT) between nuclear and organellar genomes is a common phenomenon during plant evolution. Gossypium is a useful model to evaluate the genomic consequences of IGT for both diploid and polyploid species. Here, we explore IGT among nuclear, mitochondrial, and plastid genomes of four cotton species, including two allopolyploids and their model diploid progenitors (genome donors, G. arboreum: A₂ and G. raimondii: D₅).

RESULTS

Extensive IGT events exist for both diploid and allotetraploid cotton (Gossypium) species, with the nuclear genome being the predominant recipient of transferred DNA followed by the mitochondrial genome. The nuclear genome has integrated 100 times more foreign sequences than the mitochondrial genome has in total length. In the nucleus, the integrated length of chloroplast DNA (cpDNA) was between 1.87 times (in diploids) to nearly four times (in allopolyploids) greater than that of mitochondrial DNA (mtDNA). In the mitochondrion, the length of nuclear DNA (nuDNA) was typically three times than that of cpDNA. Gossypium mitochondrial genomes integrated three nuclear retrotransposons and eight chloroplast tRNA genes, and incorporated chloroplast DNA prior to divergence between the diploids and allopolyploid formation. For mitochondrial chloroplast-tRNA genes, there were 2-6 bp conserved microhomologies flanking their insertion sites across distantly related genera, which increased to 10 bp microhomologies for the four cotton species studied. For organellar DNA sequences, there are source hotspots, e.g., the atp6-trnW intergenic region in the mitochondrion and the inverted repeat region in the chloroplast. Organellar DNAs in the nucleus were rarely expressed, and at low levels. Surprisingly, there was asymmetry in the survivorship of ancestral insertions following allopolyploidy, with most numts (nuclear mitochondrial insertions) decaying or being lost whereas most nupts (nuclear plastidial insertions) were retained.

CONCLUSIONS

This study characterized and compared intracellular transfer among nuclear and organellar genomes within two cultivated allopolyploids and their ancestral diploid cotton species. A striking asymmetry in the fate of IGTs in allopolyploid cotton was discovered, with numts being preferentially lost relative to nupts. Our results connect intergenomic gene transfer with allotetraploidy and provide new insight into intracellular genome evolution.

The plastomes of Astrocaryum aculeatum G. Mey. and A. murumuru Mart. show a flip-flop recombination between two short inverted repeats

The plastomes of Astrocaryum murumuru and A. aculeatum revealed a lineage-specific structural feature originated by flip-flop recombination, non-synonymous substitutions in conserved genes and several molecular markers. Astrocaryum murumuru Mart. and A. aculeatum G.Mey. are two palm species of Amazon forest that are economically important as source of food, oil and raw material for several applications. Genetic studies aiming to establish strategies for conservation and domestication of both species are still in the beginning given that the exploitation is mostly by extractive activity. The identification and characterization of molecular markers are essential to assess the genetic diversity of natural populations of both species. Therefore, we sequenced and characterized in detail the plastome of both species. We compared both species and identified 32 polymorphic SSR loci, 150 SNPs, 46 indels and eight hotspots of nucleotide diversity. Additionally, we reported a specific RNA editing site found in the ccsA gene, which is exclusive to A. murumuru. Moreover, the structural analysis in the plastomes of both species revealed a 4.6-kb inversion encompassing a set of genes involved in chlororespiration and plastid translation. This 4.6-kb inversion is a lineage-specific structural feature of the genus Astrocaryum originated by flip-flop recombination between two short inverted repeats. Furthermore, our phylogenetic analysis using whole plastomes of 39 Arecaceae species placed the Astrocaryum species sister to Acrocomia within the tribe Cocoseae. Finally, our data indicated substantial changes in the plastome structure and sequence of both species of the genus Astrocaryum, bringing new molecular markers, several structural and evolving features, which can be applied in several areas such as genetic, evolution, breeding, phylogeny and conservation strategies for both species.

Living with Two Genomes: Grafting and Its Implications for Plant Genome-to-Genome Interactions, Phenotypic Variation, and Evolution

Plant genomes interact when genetically distinct individuals join, or are joined, together. Individuals can fuse in three contexts: artificial grafts, natural grafts, and host-parasite interactions. Artificial grafts have been studied for decades and are important platforms for studying the movement of RNA, DNA, and protein. Yet several mysteries about artificial grafts remain, including the factors that contribute to graft incompatibility, the prevalence of genetic and epigenetic modifications caused by exchanges between graft partners, and the long-term effects of these modifications on phenotype. Host-parasite interactions also lead to the exchange of materials, and RNA exchange actively contributes to an ongoing arms race between parasite virulence and host resistance. Little is known about natural grafts except that they can be frequent and may provide opportunities for evolutionary innovation through genome exchange. In this review, we survey our current understanding about these three mechanisms of contact, the genomic interactions that result, and the potential evolutionary implications.

Rapid functional activation of a horizontally transferred eukaryotic gene in a bacterial genome in the absence of selection

Horizontal gene transfer has occurred between organisms of all domains of life and contributed substantially to genome evolution in both prokaryotes and eukaryotes. Phylogenetic evidence suggests that eukaryotic genes horizontally transferred to bacteria provided useful new gene functions that improved metabolic plasticity and facilitated adaptation to new environments. How these eukaryotic genes evolved into functional bacterial genes is not known. Here, we have conducted a genetic screen to identify the mechanisms involved in functional activation of a eukaryotic gene after its transfer into a bacterial genome. We integrated a eukaryotic selectable marker gene cassette driven by expression elements from the red alga Porphyridium purpureum into the genome of Escherichia coli. Following growth under non-selective conditions, gene activation events were indentified by antibiotic selection. We show that gene activation in the bacterial recipient occurs at high frequency and involves two major types of spontaneous mutations: deletion and gene amplification. We further show that both mechanisms result in promoter capture and are frequently triggered by microhomology-mediated recombination. Our data suggest that horizontally transferred genes have a high probability of acquiring functionality, resulting in their maintenance if they confer a selective advantage.